ALBERT

Notes: Total magnetic intensity contour maps for the study region (between 2°E to 10°E and 56°N to 60°N) were digitized and converted to a regular grid of 285 × 285 points. The study area measures approximately 444 km × 444 km and the grid spacing is thus 1. 56 km. The International Geomagnetic Reference Field for 1975 was gridded for the above-used net, and from the two data sets a further grid of the ▵T field was generated. A large number of profiles were constructed which were suitable for depth determinations. The regular grid ▵T data is also convenient for the computation of the second vertical derivative. Using the method of vertical prisms of Vacquier et al. (1963), a large suite of curvature-depth indices was measured to complement the depths obtained from the intensity slopes and from boreholes which reach the crystalline basement. The depth to the magnetic basement has been contoured, and the resulting map is shown to be in good agreement with what is known about the deeper geology of the study area.The work reported here is part of a research project supported by Amoco Norway, BP Petroleum Development Ltd, Elf Aquitaine, Esso Exploration and Production, Norwegian Gulf, Norsk Hydro, Mobil Exploration Norway, Norwegian Petroleum Directorate, Royal Norwegian Council for Scientific and Industrial Research (NTNF), Norske Shell, and Statoil.

Notes: Continuous, single-channel reflection profiling has been carried out in PVC-lined boreholes, primarily with the aim of ascertaining the position of an old subsurface gas storage tunnel on a proposed dam site.Tube wave reflection patterns thus generated have been interpreted in terms of sediment rigidity and shear wave velocity, and these results could be compared with some independent data. It is interesting to note that, within the well section penetrating Tertiary clays, the velocity of the hydraulic transients apparently was not affected by the PVC casing, which might be explained by a tight coupling between casing and clay wall. In such situations, tube waves turn out a straightforward tool for the determination of shear wave velocity and the derivation of dynamic elastic moduli of unconsolidated sediments.Further applications of the study of the distribution of seismic velocities on the dam site dealt with the consolidation history of the clays. A level of abnormally low P-wave velocities has been detected and interpreted as a gas-charged horizon which, by its coincidence with the base level of clay diapirs, might be considered to have contributed to clay flowage in past geological times. Data about maximum past burial depth, derived from shear wave velocities, turned out to be in agreement with results from consolidation testing.

Notes: The most common source of seismic energy is an explosion at some depth in a borehole. The radiated waves are reflected not only at the subsurface layers but also at the free surface. The earth's surface acts as a generator of both P- and S-waves.If the source depth is much less than the dominant wavelength the reflected waves resemble closely the waves generated by a single force. Theoretical seismograms were computed with different methods to look for the relevance of the surface-reflected waves. The numerical experiments show reflected shear waves even for small shotpoint—receiver distances. Due to their polarization these waves can be detected most easily on in-line horizontal geophones. The existence of these waves was examined during a conventional survey in Northern Germany. Conventional data analysis shows a large variability in the νp/νs ratio. The method used here produced a shear-wave section with a rather good signal-to-noise ratio down to 4 s S-wave reflection time.

Notes: An inverse problem is one in which the parameters of a model are determined from measured seismic data. Important to the solution of inverse problems is the issue of whether or not a solution exists. In this paper we show, in a constructive manner, that a solution does exist to the specific inverse problem of determining the parameters of a horizontally stratified, lossless, isotropic and homogeneous layered system that is excited by a non-normal incidence (NNI) plane wave. Mode conversion between P- and S-waves is included.We develop a seven-step layer-recursive procedure for determining all of the parameters for layer j. These parameters are P-wave and S-wave velocities and angles of incidence, density, thickness, traveltimes, and reflection- and transmission-coefficient matrices. Downward continuation of data from the top of one layer to the top of the next lower layer is an important step in our procedure, just as it is in normal incidence (NI) inversion. We show that, in order to compute all parameters of layer j, we need to (and can) compute some parameters for layer j+ 1. This is a non-causal phenomenon that seems to be necessary in NNI inversion but is not present in NI inversion.

Notes: Previously ignored characteristics of the seismic recording instrument are presently experienced as limitations as more sophisticated interpretive methods using wider frequency ranges are developed to extract stratigraphic information from seismic land data for hydrocarbon and mineral exploration. Most of these limitations arise from inadequate characteristics of the first element of the seismic instrument: the geophone. A geophone does not faithfully follow the motion of the earth for higher frequencies due to poor geophone-earth coupling. This filtering effect brings about time shifts that are dependent on the frequency and the soil type. A geophone can also produce spurious outputs, brought about by the motion of the suspended part of the geophone, with a magnitude comparable to that of the desired output. The suspension is made very compliant to obtain the required sensitivity. A compliant suspension, however, gives a large sag. The geophone can therefore only be used in one position, tolerating little tilt. A compliant suspension also widens the traveling range of the movable part. Minor sensitivity changes with travel are then noticeable as nonlinearity, since the surface wave is large with respect to the reflected wave. A compliant suspension is usually realized in the form of thin, spirally shaped spring-spiders. Such suspensions exhibit transverse or rotational resonances that are in or close to the seismic frequency band. Excited by ground roll, they can produce considerable undesirable output.The novel geophone we describe is a light-weight (17 g) acceleration-sensitive transducer which gives good ground coupling and partial correction for the increasing damping in the earth with increasing frequencies. It employs internal hybrid electronics for a magnetodynamic velocity-nulling feedback system. Velocity nulling makes the movable part of the geophone virtually rigid with respect to the housing. This makes the geophone characteristics independent of the suspension. The springs used are stiff in a transverse and rotational direction so that the suspension resonances are well outside the useful frequency band. This suspension also allows the geophone to be used in any orientation while being only sensitive to the vibration component along the main axis. The feedback system makes the sensitivity flat within 1 dB from 2 Hz to 500 Hz, with a phase tolerance smaller than 5°. The geophone is robust, has no moving internal wires, employs a current output [sensitivity 1 mA/(m s−2)] and internal gain so that the signal-to-cable-noise ratio is improved. This type of output allows parallel connection without any interaction between the geophones.

Notes: If the conductivity of any one of the layers of a horizontally stratified earth varies exponentially with depth with or without a discontinuity at the interface, the corresponding expressions for apparent resistivity for Wenner- and Schlumberger-sounding arrays can be formulated. The general case has been broadly divided into three categories for mathematical simplicity. All previous discussions of this problem can be regarded as particular cases of the present study.

Notes: In migration procedures, the velocity profile of the subsurface is the most important input information. Since, in general, this information is only approximately known, errors in the migration output due to errors in the velocity input occur in all practical applications. In migration, velocity errors and depth errors can be interchanged. This interchange property is perfect in the paraxial approximation. From this result it follows that migration with incorrect velocities may still yield correctly migrated data if the imaging principle is modified. This attractive property can be used in the stripping version of migration (recursive migration).

Notes: The presence of the water layer in marine seismic prospecting provides an effective waveguide for acoustic energy trapped between the sea-bed and the sea-surface. This energy persists to large ranges and can be the dominant early feature on far-offset traces. On airgun records, there is commonly a lower frequency set of arrivals following the water-trapped waves. These arrivals are not as obvious with higher frequency watergun sources. By using a combination of intercept-time/slowness (τ—p) mapping on observational data and theoretical modelling, we are able to identify the origin of the events. If a very rapid increase in a seismic wavespeed occurs beneath the sea-bed sediments, a new waveguide is formed bounded by the sea surface and this transition zone. The low frequency waves are principally guided within this thicker waveguide. Numerical filtering in the τ—p domain followed by trace reconstruction is very effective in removing the low frequency noise.

Notes: A method for interval velocity analysis is formulated on the basis of wavefield extrapolation, i.e., on the basis of wave-equation migration. When this scheme is applied to multioffset seismic sections or to an ensemble of CMP gathers, it allows for the proper treatment of dipping events. The underlying assumptions are that local velocities should be derived from data associated with events within the interval under consideration. To minimize the effect of the region above the layer of interest, the data are first extrapolated to the top of the analysis interval. Subsequent analysis of these data then pertains to the events within this interval. Velocity estimation consists of repeated wavefield extrapolations through the analysis interval using a set of trial velocities. The optimal velocity is chosen on the basis of coherency measures designed to express the collective phase agreement among a set of offset Fourier modes. The reliability of this approach to interval velocity estimation is demonstrated on synthetic multi-offset data.