Showing papers in "Geophysical Prospecting in 1994"

The relationships between the velocities, attenuations and petrophysical properties of reservoir sedimentary rocks1

Abstract: The authors have measured the velocities and attenuations of compressional and shear waves in 29 water-saturated samples of sandstones and shales at a confining pressure of 60 MPa and at frequencies of about 0.85 MHz. The measurements were made using a pulse echo method in which the samples (diameter 5 cm, length 1.5 cm to 2.5 cm) were placed between perspex buffer rods inside a high-pressure cell. The velocity of each seismic wave was determined from the travel time difference of equivalent phase points (corrected for diffraction effects) of the signals reflected from the top and from the base of each samples. Attenuation was determined in a similar way by comparison of the diffraction corrected amplitudes of the signals. The attenuation data are presented as quality factors'': Q[sub p] and Q[sub s] for compressional and shear waves respectively and shear waves respectively. The results show that Qs is strongly correlated with V[sub s], that Q[sub p] is weakly correlated with V[sub p], and that Q[sub p] is strongly correlated with Q[sub s] [center dot] Q[sub p] is strongly dependent on the volume percentage of the assemblage of intra-pore minerals, whether they are clays or carbonates. It is concluded that themore » attenuation mechanism is due to the local fluid flow arising from the differential dilation of the solid rock frame and the intra-pore mineral assemblage, which is a results of their very different elastic moduli.« less

2D finite‐difference elastic wave modelling including surface topography1

Abstract: A 2D numerical finite-difference algorithm accounting for surface topography is presented. Higher-order, dispersion-bounded, cost-optimized finite-difference operators are used in the interior of the numerical grid, while non-reflecting absorbing boundary conditions are used along the edges. Transformation from a curved to a rectangular grid achieves the modelling of the surface topography. We use free-surface boundary conditions along the surface. In order to obtain complete modelling of the effects of wave propagation, it is important to account for the surface topography, otherwise near-surface effects, such as scattering, are not modelled adequately. Even if other properties of the medium, for instance randomization, can improve numerical simulations, inclusion of the surface topography makes them more realistic.

A study of the geoelectrical properties of peatlands and their influence on ground‐penetrating radar surveying1

Abstract: Geophysical surveys and chemical analyses on cores were carried out in three Ontario peatlands, from which we have gained a better understanding of the peat properties that control the geophysical responses. The electrical conductivity depends linearly on the concentration of total dissolved solids in the peat pore waters and the pore waters in turn bear the ionic signatures of the underlying mineral sediments. The ionic concentration, and thus the electrical conductivity, increase linearly from the surface to basement. The average bulk electrical conductivity of peatlands at Ellice Marsh, near Stratford, and at Wally Creek Area Forest Drainage Project, near Cochrane, are of the order of 25 mS/m. The Mer Bleue peatland, near Ottawa, has extremely high electrical conductivity, reaching levels of up to 380 mS/m near the base of the peat. The Mer Bleue peatland water has correspondingly high values of total dissolved solids, which originate from the underlying Champlain Sea glaciomarine clays. The dielectric permittivity in peats is largely controlled by the bulk water content. Ground penetrating radar can detect changes in water content greater than 3%, occurring within a depth interval less than 15 cm. The principal peatland interfaces detected are the near-surface aerobic to anaerobic transition and the peat to mineral basement contact. The potential for the successful detection of the basement contact using the radar can be predicted using the radar instrument specifications, estimates of the peatland depth, and either the bulk peat or the peat pore water electrical conductivities. Predicted depths of penetration of up to 10 m for Ellice Marsh and Wally Creek exceed the observed depths of 1 to 2 m. At Mer Bleue, on the other hand, we observe that, as predicted, a 100 MHz signal will penetrate to the base of a 2 m thick peat but a 200 MHz signal will not.

Elastic modelling on a grid with vertically varying spacing1

Abstract: We present a discrete modelling scheme which solves the elastic wave equation on a grid with vertically varying grid spacings. Spatial derivatives are computed by finite-difference operators on a staggered grid. The time integration is performed by the rapid expansion method. The use of variable grid spacings adds flexibility and improves the efficiency since different spatial sampling intervals can be used in regions with different material properties. In the case of large velocity contrasts, the use of a non-uniform grid avoids spatial oversampling in regions with high velocities. The modelling scheme allows accurate modelling up to a spatial sampling rate of approximately 2.5 gridpoints per shortest wavelength. However, due to the staggering of the material parameters, a smoothing of the material parameters has to be applied at internal interfaces aligned with the numerical grid to avoid amplitude errors and timing inaccuracies. The best results are obtained by smoothing based on slowness averaging. To reduce errors in the implementation of the free-surface boundary condition introduced by the staggering of the stress components, we reduce the grid spacing in the vertical direction in the vicinity of the free surface to approximately 10 gridpoints per shortest wavelength. Using this technique we obtain accurate results for surface waves in transversely isotropic media.

The potential of a noise‐reducing antenna for surface NMR groundwater surveys in the earth's magnetic field1

Abstract: A method of non-invasive NMR in the earth's field has been developed and is now used for groundwater surveys to depths of investigation of 100 m or more. A circular wire loop of diameter 100 m, laid out on the ground, is employed to excite and receive the NMR signal in the earth's field. However, in areas with high electromagnetic noise, the NMR measurements may be inaccurate. To overcome this problem, a noise-reducing figure-of-eight-shaped antenna, consisting of two touching coils each of diameter 50 m, has been utilized. Using this antenna, the NMR signal has been calculated for different depths of water-saturated layers with various inclinations of the geomagnetic field. The model calculations and experimental data have been compared and found to be mutually consistent. The two-coil antenna is shown to be suitable for studies at depths of up to 30–40 m, which is of practical importance for engineering geology.

Inversion of reflection seismograms by differential semblance analysis: algorithm structure and synthetic examples1

Abstract: Seismograms predicted from acoustic or elastic earth models depend very non-linearly on the long wavelength components of velocity. This sensitive dependence demands the use of special variational principles in waveform-based inversion algorithms. The differential semblance variational principle is well-suited to velocity inversion by gradient methods, since its objective function is smooth and convex over a large range of velocity models. An extension of the adjoint state technique yields an accurate estimate of the differential semblance gradient. Non-linear conjugate gradient iteration is quite successful in locating the global differential semblance minimum, which is near the ordinary least-squares global minimum when coherent data noise is small. Several examples, based on the 2D primaries-only acoustic model, illustrate features of the method and its performance.

Using fractal crustal magnetization models in magnetic interpretation1

Abstract: Evidence from borehole susceptibility logs and the spectral analysis of aeromagnetic data suggests that the three-dimensional distribution of magnetization within the crust can be described as fractal. This property can be exploited in magnetic interpretation methods which explicitly require statistical information on the spatial variation of magnetization. Specifically, we address the problem of magnetic source depth estimation through downward continuation and gridding aeromagnetic survey data using the method of kriging. When magnetic data are continued downwards the depth at which the power spectrum flattens out (the ‘white’ depth) can be taken to be an estimate of the top of the source distribution. This procedure assumes that individual sources are uncorrelated with each other. Taking into account the correlation of the magnetization using a fractal description leads to a reduction in this depth estimate. Gridding of randomly distributed magnetic measurements using kriging requires an estimate of the covariance of the data. Compared with the assumption of a white (uncorrelated) magnetization distribution, using fractal covariances for kriging produces gridded estimates which more closely reflect the statistics of the underlying magnetization process and produce maps with a justifiable degree of smoothness.

Pore pressure profiles in fractured and compliant rocks1

Abstract: Fluid permeability in fractured rocks is sensitive to pore-pressure changes. This dependence can have large effects on the flow of fluids through rocks. The authors define the permeability compliance [gamma] = 1/k([partial derivative]k/[partial derivative]p[sub p])[sub pc], which is the sensitivity of the permeability k to the pore pressure p[sub p] at a constant confining pressure p[sub c], and solve the specific problems of constant pressure at the boundary of a half-space, a cylindrical cavity and a spherical cavity. The results show that when the magnitude of permeability compliance is large relative to other compliances, diffusion is masked by a piston-like pressure profile. The authors expect this phenomenon to occur in highly fractured and compliant rock systems where [gamma] may be large. The pressure profile moves rapidly when fluids are pumped into the rock and very slowly when fluids are pumped out. Consequently, fluid pressure, its history and distribution around injection and production wells may be significantly different from pressures predicted by the linear diffusion equation. The propagation speed of the pressure profile, marked by the point where [partial derivative]p[sub p]/[partial derivative]x is a maximum, decreases with time approximately as [radical]t, and the amplitude of the profile also dissipates with timemore » (or distance). The effect of permeability compliance can be important for fluid injection into and withdrawal from reservoirs. For example, excessive drawdown could cause near-wellbore flow suffocation. Also, estimates of the storage capacity of reservoirs may be greatly modified when [gamma] is large. The large near-wellbore pressure gradients caused during withdrawal by large [gamma] can cause sanding and wellbore collapse due to excessive production rates.« less

A method to estimate the total magnetization direction from a distortion analysis of magnetic anomalies1

Abstract: Knowledge of the declination and inclination of the total and induced magnetization vectors is normally required for the interpretation and analysis of magnetic anomalies A new method of estimating the direction of the total magnetization vector of magnetized rocks from magnetic anomalies is proposed The unknown declination and inclination ( D * T and I * T ) can be found by applying a reduction-to-the-pole operator to the measured anomalies for different couples of total magnetization direction parameters ( D T and I T ) and by observing the variation of the anomaly minimum as a function of both D T and I * T and D * T are estimated using the maximum of this function Comparing our method to previous methods, one advantage is that our estimates are not zero-level dependent; furthermore, the method allows inclinations to be well estimated, with the same accuracy as declinations; finally declinations are not underestimated Our method is applied to a real case and meaningful results are obtained; it is shown that the feasibility of the method is improved by removing the low-frequency components

Seismic velocities in fractured rocks: an experimental verification of Hudson's theory1

Abstract: Flow of fluids in many hydrocarbon reservoirs aquifers is enhanced by the presence of cracks and fractures. These cracks could be detected by their effects on propagation of compressional and shear waves through the reservoir: several theories, including Hudson`s, claim to predict the seismic effects of cracks. Although Hudson`s theory has already been used to calculate crack densities from seismic survey`s, the predictions of the theory have not yet been tested experimentally on rocks containing a known crack distribution. This paper describes an experimental verification of the theory. The rock used, Carrara marble, was chosen for its uniformity and low porosity, so that the effect of cracks would not be obscured by other influences. Cracks were induced by loading of laboratory specimens. Velocities of compressional and shear waves were measured by ultrasound at 0.85 MHz in dry and water-saturated specimens at high and low effective pressures.The cracks were then counted in polished sections of the specimens. In ``dry`` specimens with both dry and saturated cracks, Hudson`s theory overpredicted observed crack densities by a constant amount that is attributed to the observed value being systematically underestimated. The theory made poor predictions for fully saturated specimens. Shear-wave splitting, caused by anisotropy duemore » to both crystal and crack alignment, was observed. Cracks were seen to follow grain boundaries rather than the direction of maximum compression due to loading. The results demonstrate that Hudson`s theory may be used in some cases to determine crack and fracture densities from compressional- and shear-wave velocity data.« less

Sign‐change correction for prestack migration of P‐S converted wave reflections1

Abstract: The polarization direction or 'sign’ of reflected converted P–S waves depends upon the angle of incidence of the incident P-wave. Sign reversal due to reversal of the angle of incidence is often encountered and is an impediment to P–S wave processing and imaging, because when P–S events or P-S migrated images with mixed signs are stacked, destructive interference occurs. We have created and demonstrated a means of correcting for this reversal. To do this, a P-wave angle of incidence is calculated for every point in the image space. This is done by calculating a P–S reflected waveform for every point, by extrapolating the reflected S-wavefield backwards from the receiver line, and then cross-correlating this waveform with the S-wave reflections observed at the receiver line. A multiplier, (sgn α) is assigned to each point in the image space, where α is the angle of incidence of the P-wave. The multiplier was applied to a set of prestack reverse time migration images derived from a cross-borehole physical elastic model data set. The improvement in the stacked image when the sign correction is applied is spectacular. The P-S image quality is comparable to, or better than, stacked migrated P-P images. The method appears to be applicable to all reflection modes and to all recording geometries.

Inversion for seismic anisotropy using genetic algorithms1

Abstract: A general inversion scheme based on a genetic algorithm is developed to invert seismic observations for anisotropic parameters. The technique is applied to the inversion of shear-wave observations from two azimuthal VSP data sets from the Conoco test site in Oklahoma. Horizontal polarizations and time-delays are inverted for hexagonal and orthorhombic symmetries. The model solutions are consistent with previous studies using trial and error matching of full waveform synthetics. The shear-wave splitting observations suggest the presence of a shear-wave line singularity and are consistent with a dipping fracture system which is known to exist at the test site. Application of the inversion scheme prior to full waveform modelling demonstrates that a considerable saving in time is possible whilst retaining the same degree of accuracy.

Using the pseudospectral technique on curved grids for 2D acoustic forward modelling1

Abstract: The Fourier pseudospectral method has been widely accepted for seismic forward modelling because of its high accuracy compared to other numerical techniques. Conventionally, the modelling is performed on Cartesian grids. This means that curved interfaces are represented in a ‘staircase fashion‘causing spurious diffractions. It is the aim of this work to eliminate these non-physical diffractions by using curved grids that generally follow the interfaces. A further advantage of using curved grids is that the local grid density can be adjusted according to the velocity of the individual layers, i.e. the overall grid density is not restricted by the lowest velocity in the subsurface. This means that considerable savings in computer storage can be obtained and thus larger computational models can be handled. One of the major problems in using the curved grid approach has been the generation of a suitable grid that fits all the interfaces. However, as a new approach, we adopt techniques originally developed for computational fluid dynamics (CFD) applications. This allows us to put the curved grid technique into a general framework, enabling the grid to follow all interfaces. In principle, a separate grid is generated for each geological layer, patching the grid lines across the interfaces to obtain a globally continuous grid (the so-called multiblock strategy). The curved grid is taken to constitute a generalised curvilinear coordinate system, where each grid line corresponds to a constant value of one of the curvilinear coordinates. That means that the forward modelling equations have to be written in curvilinear coordinates, resulting in additional terms in the equations. However, the subsurface geometry is much simpler in the curvilinear space. The advantages of the curved grid technique are demonstrated for the 2D acoustic wave equation. This includes a verification of the method against an analytic reference solution for wedge diffraction and a comparison with the pseudospectral method on Cartesian grids. The results demonstrate that high accuracies are obtained with few grid points and without extra computational costs as compared with Cartesian methods.

On reducing ambiguity in the interpretation of transient electromagnetic sounding data1

Abstract: An attempt to resolve non-uniqueness in the interpretation of transient electromagnetic (TEM) sounding data using measured data alone is made. It is shown in the various examples studied that sufficiently early time measurements can be the determining factor in reducing the ambiguity caused by model equivalence. The early delay times thus play a dual role in transient soundings: they are responsible for resolving shallow structures and they may eliminate the ambiguity in the interpretation of geoelectric parameters of deeper targets. This is illustrated by the results of a follow-up TEM survey at the Dead Sea coast of Israel where the use of supplementary early time measurements allowed non-uniqueness in the determination of the depth to fresh/saline groundwater interface to be resolved.

An analytic signal approach to the interpretation of total field magnetic anomalies1

Abstract: The analytic signal (AS) is defined as the square root of the sum of the squares of the vertical and the two horizontal derivatives of the total magnetic field ΔT. This paper verifies theoretically that peaks of the AS correlate directly with their magnetic causative bodies and are positioned symmetrically over them, i.e. the main feature of the AS is that it is independent of the inclination of the magnetic field. This avoids the difficulties that are often faced in the conventional process of reduction to pole for ΔT, when the direction of magnetization of the causative bodies is not known. In addition, the AS has characteristics similar to the derivative features of the magnetic field, so that it is very sensitive to edge effects of the causative magnetic bodies. The theoretical derivations are tested by comparison with calculations on models, and, in a field example from Hunan Province, China, the AS is applied successfully to the interpretation of ΔT, whereas the conventional process of reduction to pole fails, due to the reverse magnetization of the causative body.

Radar experiments in isotropic and anisotropic geological formations (granite and schists)1

Abstract: In order to understand various aspects of radar wave propagation, a survey of electromagnetic wave behaviour relative to the geological characteristics of the formations prospected was undertaken The sites chosen for the tests were a granite quarry and an underground schist working By investigating an electrically resistive isotropic site and a conductive anisotropic site, it was demonstrated that non-conventional use of a radar system (antennae raised, various orientations of the transmitter/receiver, etc) could improve data quality, and could allow information other than reflector depth to be collected (volume scattering intensity, isotropy, etc) By studying wave propagation velocities, we underlined the difficulties encountered in establishing a velocity versus depth law, despite recourse to seismic data processing, such as NMO corrections The results of field experiments, complemented by laboratory measurements of dielectric permittivities, clearly showed anisotropy effects: in the case of a path that is perpendicular to the schistosity plane, an electromagnetic wave propagates more slowly and is more attenuated than a wave parallel to the schistosity plane

An airborne tensor VLF system. From concept to realization1

Abstract: Radio signals from very low frequency (VLF) transmitters distributed world-wide have been used for several decades to study the lateral variations of the electrical conductivity in the upper few hundred metres of the earth's crust. Traditionally, in airborne applications, the total magnetic fields from one or two transmitters are measured to form the basis for construction of maps that primarily show those conductive structures that are parallel or subparallel to the direction to the transmitters. The tensor VLF technique described in this paper makes use of all signals available in a predefined frequency band to construct transfer functions relating the vertical magnetic field and the two horizontal magnetic field components. These transfer functions are uniquely determined for a particular measuring site and contain information about the lateral conductivity variations in all directions. First experiences with real field data, acquired during a test survey in Sweden, show that maps of the so-called peaker , the spatial divergence of the transfer functions, give an image of the conducting structures. Most of the structures can be correlated to small valleys filled with conducting sediments or valleys underlain by conductive fracture zones in the crystalline rocks.

Upward continuation of potential fields from a polyhedral surface1

Abstract: An equivalent source procedure is derived for upward continuation of unevenly spaced gravity and magnetic data. The dipole layer is placed on a topographic relief approximated by a polyhedral surface, the stations being the vertices of the triangular faces. The dipoles have linear magnitudes, being directed along the normal vector over each triangle. The unknown values of the dipole magnitudes at each station are obtained by a suitable modification of the usual integral equation considering the discontinuity of the normal vector at each vertex of the dipole surface. Profile data processing is also studied. A numerical test outlines the accuracy and the limitations of the model for the case of a magnetic field significantly perturbed by a rough topographic relief.

Gravity interpretation of sedimentary basins with hyperbolic density contrast

Abstract: The variation in the density of sediments with depth in a sedimentary basin can be represented by a hyperbolic function. Gravity anomaly expressions for a 2D vertical prism and an asymmetric trapezium with a hyperbolic density distribution are derived in a closed form. These are used in inverting the gravity anomaly of a sedimentary basin with variable density. Firstly, the basin is viewed as a series of prisms juxtaposed with each other. The initial thickness of each prism is obtained from the gravity anomaly at its centre, based on the gravity anomaly of an infinite slab with a hyperbolic density contrast. These thicknesses are improved, based on the differences between the observed and the calculated anomalies. For an improved rate of convergence of the solution, these thicknesses may alternatively be refined using the well-known ridge regression technique. Secondly, the basin is approximated by an asymmetric trapezium and its anomalies are inverted for the parameters of the trapezium using the ridge regression. Since this approximation serves to oversimplify the floor of the basin, it must be used only when the sediment-basement interface has minor undulations. The results of a hypothetical case and two field cases (the San Jacinto Graben, California and the Godavari Graben, southern India) are presented. In both field cases, the interpreted depths are comparable with the real ones, proving the validity of the assumption of a hyperbolic density distribution of the sediments in the two basins considered.

Normal moveout in focus1

Abstract: Improving the accuracy of NMO corrections and of the corresponding interval velocities entails implementing a better approximation than the formula used since the beginning of seismic processing. The exact equations are not practical as they include many unknowns. The approximate expression has only two unknowns, the reflection time and the rms velocity, but becomes inaccurate for large apertures of the recording system and heterogeneous vertical velocities. Several methods of improving the accuracy have been considered, but the gains do not compensate for the dramatic increase in computing time. Two alternative equations are proposed: the first containing two parameters, the reflection time and the focusing time, is not valid for apertures much greater than is the standard formula, but has a much faster computing time and does not stretch the far traces; the other, containing three parameters, the reflection time, like focusing time and the tuning velocity, retains high frequencies for apertures about twice those allowed by the standard equation. Its computing time can be kept within the same limits. NMO equations, old and new, are designed strictly for horizontal layering, but remain reliable as long as the rays travel through the same layers in both the down and up directions.more » An equation, similar to Dix's formula, is given to compute the interval velocities. The entire scheme can be automated to produce interval-velocity sections without manual picking.« less

A physical approach to computing magnetic fields1

Abstract: The superposition integral expressing the field due to a magnetic source body is relatively simple to evaluate in the case of a homogeneous magnetization. In practice this generally requires that any remnant component is uniform and the susceptibility of the body is sufficiently low to permit the assumption of a uniform induced magnetization. Under these conditions the anomalous magnetic field due to a polyhedral body can be represented in an intuitive and physically appealing manner. It is demonstrated that the components of the magnetic field H can be expressed as a simple combination of the potentials due to two elementary source distributions. These are, firstly, a uniform double layer (normally directed dipole moment density) located on the planar polygonal faces of the body and, secondly, a uniform line source located along its edges. In practice both of these potentials (and thus the required magnetic field components) are easily computed. The technique is applicable to polyhedra with arbitrarily shaped faces and the relevant expressions for the magnetic field components are suitable for numerical evaluation everywhere except along the edges of the body where they display a logarithmic singularity.

The fast Hankel transform1

Abstract: An efficient algorithm is presented to compute the Hankel transform. The algorithm yields simultaneously all the required weights for a given order of the Bessel function using the fast Fourier transform. An additional shift is introduced to the filter abscissa besides Koefoed's shift to give a better filter performance.

Definitions of apparent resistivity for the presentation of magnetotelluric sounding data1

Abstract: A new definition of apparent resistivity for the presentation of magnetotelluric sounding data is proposed. The new definition is based on the frequency-normalized impedance function. Both the existing and proposed definitions of apparent resistivity are analysed theoretically and are compared using model curves computed for a 1D earth model. Apparent resistivity curves computed using the proposed definition are a better approximation to the true resistivity values of the subsurface layers. In addition, the layers are more noticeable on the apparent resistivity curves, which is an advantage, especially for the ascending and descending type of apparent resistivity curve.

Statistical analysis of non‐stationary voltage recordings in geoelectrical prospecting1

Abstract: We study the geoelectrical problem of picking out the useful signal from voltage time series, monitored under conditions of a low signal-to-noise ratio and non-stationary noise. Statistical tests performed at different sites show that geoelectrical noise often belongs to the class of non-stationary phenomena with non-Gaussian probability distributions. In such cases, the application of conventional methods of geoelectrical useful signal extraction, based on the stationary white-noise assumption, gives biased estimates. For the on-line processing of geoelectrical recordings, we recommend the use of the periodogram technique combined with the Kolmogorov–Smirnov test, a suitable algorithm of which is described in detail. The suggested procedure allows data acquisition to stop as soon as the useful signal power is estimated with a relative error smaller than a predetermined value. Finally, we compare the suggested procedure with the autoregressive approach. The previously used and simpler periodogram method, applied to the solution of problems of this kind, appears to give better performances than the autoregressive analysis.

Transversely isotropic media equivalent to thin isotropic layers1

Abstract: Any set of isotropic layers is equivalent, in the long wavelength limit, to a unique transversely isotropic (TI) layer; to find the elastic moduli of that layer is a solved problem. The converse problem is to find a set of isotropic layers equivalent to a given TI media. Here, explicit necessary and sufficient conditions on the TI stiffness moduli for the existence of an equivalent set of isotropic layers are found by construction of a minimal decomposition consisting of either two or three isotropic constituent layers. When only two constituents are required, their elastic properties are uniquely determined. When three constituents are required, two have the same Poisson's ratio and the same thickness fraction, and even then there is a one-parameter family of satisfactory minimal decompositions. The linear slip model for fractured rock (aligned fractures in an isotropic background) yields a restricted range of transverse isotropy dependent on only four independent parameters. If the ratio of the normal to tangential fracture compliance is small enough, the medium is equivalent to thin isotropic layering and in general its minimal decomposition consists of three constituents.

A gravimetric 3D global inversion for cavity detection

Abstract: A gravimetric survey, covering a site 200 m square, was carried out in order to locate karstic cavities. After eliminating the regional trend using a polynomial fit, the residual is modelled by least-squares prediction. Correlated signals for several wavelengths are detected. The inversion of these anomalies is performed by a global 3D adjustment using spherical bodies as models. The adjustment is repeated in order to obtain a stable configuration. The results show the probable presence of a system of cavities and galleries. Data collected from boreholes and the subsequent appearance of sink-holes are consistent with the results.

Location of secondary faults from cross‐correlation of the second vertical derivative of gravity anomalies1

Abstract: A method is described to locate secondary faults, which can be difficult to identify on the Bouguer gravity map. The method is based on cross‐correlation between the theoretical anomaly due to a vertical step and the second vertical derivative of the Bouguer anomaly. Faults are located from the closed maxima and minima on the cross‐correlation contour map calculated for two perpendicular directions. One‐dimensional model computations show that the magnitude of the extremum of the cross‐correlation is related to the depth to the top of the hanging wall and the throw of the fault. Application of the method to the Bouguer gravity map of the former mouth of the Yellow River in the Shengli Oilfield area near the Bo Hai Sea shows the effectiveness of the method.

Cost-effective 3D one-pass depth migration1

Abstract: A crucial point in the processing of 3D seismic data is the migration step, both because of its 3D nature and the computational cost involved. The efficiency and accuracy of 3D migration are determined by the wavefield extrapolation technique employed. Wavefield extrapolation based on second-order differential operators of variable-length is very efficient and accurate at the same time. Compared to migration based on the McClellan transform and operator splitting, the use of variable-length second-order differential operators offers significant advantages. The 3D migration operator has an almost perfect circular symmetry. No positioning errors in the 45° azimuth between the in-line and cross-line directions are evident. The method is, in practice, only limited by spatial aliasing and does not require expensive interpolation of data to reduce numerical artifacts. This reduces the computational cost of 3D one-pass depth migration by a large factor.

Resistivity measurements on the sea bottom to map fracture zones in the bedrock underneath sediments1

Abstract: Refraction seismics with the shotpoints and the hydrophone cable on the sea-bottom, have become the standard geophysical method for investigating rock quality before constructing offshore tunnels in Norway. In connection with the construction of a sub-sea tunnel by the Norwegian Public Roads Administration, research work was carried out to compare two low-velocity zones, indicated by refraction seismics with other methods. A special resistivity cable for pole-dipole measurements on the sea-floor, with 10 m between the electrodes, was constructed. A 200 m long profile, crossing the two low-velocity zones, was measured with all combinations of electrode distances. The two zones were detected as low-resistivity zones. A special data processing technique to enhance the anomalies is described. Resistivity soundings in a seawater environment to detect fracture zones in the bedrock underneath the bottom sediments, are discussed. It is concluded that severely fractured zones, which may cause difficulties for the tunnel construction, can be detected both with sea-surface and sea-floor arrays using long electrode spacings.

Relation of in situ resistivity to water content in salt rocks1

Abstract: The investigation of water in salt-rock formations is of particular relevance to underground nuclear waste repositories. In the Asse salt-mine (Germany) a study into the relationship of in situ resistivity to water content has been made. Measurements were carried out in older rock-salt using an electrode array in boreholes, an electrode profile in a drift and small resistivity sensors in and around a drift seal. Further measurements were make on moist zones in a contact area of younger rock-salt and carnallitite and also in older rock-salt with anhydrite bands using electrode profiles in the drifts. The resistivities range from 10[sup 2] [Omega]m to 10[sup 6] [Omega]m. Corresponding probes have water contents from 0.01% to 1.3%. A definite relationship between resistivity and water content is revealed which can be described by Archie's law using a cementation factor m of 1.9. Porosities are between 0.08% and 1.4% and the saturations vary considerably. An explicitly influence of saturation on resistivity cannot be discovered using the present data. The results enable us to estimate the in situ water content and the order of the in situ porosity using resistivity surveys at different scales. This increases significantly the safety of a nuclear repository site.