- Home
- A-Z Publications
- Geophysical Prospecting
- Previous Issues
- Volume 44, Issue 3, 1996
Geophysical Prospecting - Volume 44, Issue 3, 1996
Volume 44, Issue 3, 1996
-
-
The relationship between the first Fresnel zone and the normalized geometrical spreading factor1
By Jianguo SunAbstractConventionally, the Fresnel zone and the geometrical spreading factor are investigated separately, because they belong to different theories of wave propagation. However, if the paraxial ray method is used for establishing the Fresnel–Kirchhoff diffraction formula for a laterally inhomogeneous multilayered medium, it can be shown that the normalized geometrical spreading factor is inversely proportional to the area of the first Fresnel zone associated with the reflection point. Therefore, if no diffracting edge cuts the first Fresnel zone, the geometrical optics approximation represents the principal part of the wavefield obtained by Fresnel–Kirchhoff diffraction theory. Otherwise, the geometrical optics approximation has to be corrected by adding edge diffractions. It is also shown that Kirchhoff‐type migration and geometrical spreading factor correction both reduce the first Fresnel zone to a zone with unit area.
-
-
-
Multiple suppression by 2D filtering in the parabolic τ–p domain: a wave‐equation‐based method1
Authors Binzhong Zhou and Stewart GreenhalghAbstractThe filter for wave‐equation‐based water‐layer multiple suppression, developed by the authors in the x‐t, the linear τ‐p, and the f‐k domains, is extended to the parabolic τ‐2 domain. The multiple reject areas are determined automatically by comparing the energy on traces of the multiple model (which are generated by a wave‐extrapolation method from the original data) and the original input data (multiples + primaries) in τ‐p space. The advantage of applying the data‐adaptive 2D demultiple filter in the parabolic τ‐p domain is that the waves are well separated in this domain. The numerical examples demonstrate the effectiveness of such a dereverberation procedure. Filtering of multiples in the parabolic τ‐p domain works on both the far‐offset and the near‐offset traces, while the filtering of multiples in the f‐k domain is effective only for the far‐offset traces.
Tests on a synthetic common‐shot‐point (CSP) gather show that the demultiple filter is relatively immune to slight errors in the water velocity and water depth which cause arrival time errors of the multiples in the multiple model traces of less than the time dimension (about one quarter of the wavelet length) of the energy summation window of the filter. The multiples in the predicted multiple model traces do not have to be exact replicas of the multiples in the input data, in both a wavelet‐shape and traveltime sense. The demultiple filter also works reasonably well for input data contaminated by up to 25% of random noise. A shallow water CSP seismic gather, acquired on the North West Shelf of Australia, demonstrates the effectiveness of the technique on real data.
-
-
-
Guided seismic waves in layered poroviscoelastic media for continuity logging applications: model studies1
More LessAbstractThe gas industry is continuing to concentrate its research and development efforts on new and advanced technology to improve reservoir descriptions through the producing life and development history of heterogeneous gas reservoirs. A very important aspect of this need is the ability to reduce the uncertainty of estimating probable reserves and to lower the operating costs to recover incremental reserves in producing and depleted gasfields. Established methods for reducing uncertainty in heterogeneous reservoir compartments, such as VSP and cross‐well techniques may enhance resolution, but they are currently not economically justifiable in on‐shore gasfields. Continuity logging using guided waves is an alternative approach to analysing inter‐well seismic data to confirm the continuity of heterogeneous gas reservoir compartments; in particular, the continuity of sand and shale stratigraphy in gas reservoirs.
The solution of a coupled system of differential equations based on Biot and homogenization theories is adapted to calculate guided seismic waves trapped in low‐velocity layers. The general solution is for a 3D source in a horizontally layered poroviscoelastic medium having isotropic and laterally homogeneous material properties. A unified representation of the medium that includes fluid‐solid interactions and viscoelastic losses is incorporated into the solution. The guided‐wave part of the vector wave field and fluid‐pressure of the complete wave motion in layered dissipative media is verified and used to simulate dispersion and attenuation of guided seismic waves for continuity logging applications. The results of this work suggest that the multimode wave solution is appropriate to simulate guided seismic wave signatures to indicate continuity of layered earth structures in poroviscoelastic reservoirs. In particular, the normal mode information can be used for planning continuity logging surveys and for interpreting the corresponding seismic data. Further, fluid‐pressure waveforms show that maximum amplitude normal modes can be detected at layer interfaces in fluid‐filled porous media, and the corresponding Airy phase wave groups may carry information on the formation permeability.
-
-
-
Linearized elastic parameter sections1
Authors Bjørn Ursin, Bjørn Olav Ekren and Egil TjålandAbstractContrasts in elastic parameters can be estimated directly from seismic data using offset‐dependent information in the PP reflection coefficient. A linear approximation to the PP reflection coefficient including three coefficients is fitted to the data, and relative contrasts in various elastic parameters are obtained from linear combinations of the estimated coefficients. Linearized elastic parameter sections for the contrasts in P‐wave impedance, P‐wave velocity, density, plane‐wave modulus, and the change in bulk modulus and shear modulus normalized with the plane‐wave modulus are estimated. If the average P‐ to S‐wave velocity ratio is known, linearized parameter sections including the contrast in the average P‐ to S‐wave velocity ratio and a fluid factor section can be computed. Applied to synthetic data, visual comparison of the estimated and true elastic parameter sections agree qualitatively, and the results are confirmed by an analysis of the standard deviation of the estimated parameters. The parameter sections obtained by inversion of a shallow seismic anomaly in the Barents Sea are promising, but the reliability is uncertain because neither well data nor regional trends are available.
-
-
-
3D Modelling of the electromagnetic response of geophysical targets using the FDTD method1
More LessAbstractA publicly available and maintained electromagnetic finite‐difference time domain (FDTD) code has been applied to the forward modelling of the response of 1D, 2D and 3D geophysical targets to a vertical magnetic dipole excitation. The FDTD method is used to analyse target responses in the 1 MHz to 100MHz range, where either conduction or displacement currents may have the controlling role. The response of the geophysical target to the excitation is presented as changes in the magnetic field ellipticity.
The results of the FDTD code compare favourably with previously published integral equation solutions of the response of 1D targets, and FDTD models calculated with different finite‐difference cell sizes are compared to find the effect of model discretization on the solution. The discretization errors, calculated as absolute error in ellipticity, are presented for the different ground geometry models considered, and are, for the most part, below 10% of the integral equation solutions.
Finally, the FDTD code is used to calculate the magnetic ellipticity response of a 2D survey and a 3D sounding of complicated geophysical targets. The response of these 2D and 3D targets are too complicated to be verified with integral equation solutions, but show the proper low‐ and high‐frequency responses.
-
-
-
EM target location1
By L.‐J. GeliusAbstractWe consider the problem of computing the most probable location of a target based on radar measurements of the subsurface. Our algorithm makes use of the maximum likelihood estimator (MLE), which represents a correlation between the measured data and synthetic data generated for the object of interest at different locations. Previous studies assume a plane‐wave acquisition geometry and target object(s) embedded in a uniform background. In this paper, a generalization of the MLE method is presented which is valid for discrete point sources (and receivers) and a 2D model (i.e. a 2.5D acquisition geometry). Within this formulation the treatment of a non‐uniform background model is also possible. We concentrate on geotechnical ground investigations and assume that the characteristic dimensions of the target object are in the range 1–2λ, (λ being the wavelength). The potential of the method is demonstrated employing cross‐hole radar data acquired in a controlled field experiment. The MLE result is also compared with the image obtained employing a full reconstruction method such as diffraction tomography.
-
-
-
Practical techniques for 3D resistivity surveys and data inversion1
Authors M.H. Loke and R.D. BarkerAbstractTechniques to reduce the time needed to carry out 3D resistivity surveys with a moderate number (25 to 100) of electrodes and the computing time required to interpret the data have been developed. The electrodes in a 3D survey are normally arranged in a square grid and the pole‐pole array is used to make the potential measurements. The number of measurements required can be reduced to about one‐third of the maximum possible number without seriously degrading the resolution of the resulting inversion model by making measurements along the horizontal, vertical and 45° diagonal rows of electrodes passing through the current electrode. The smoothness‐constrained least‐squares inversion method is used for the data interpretation. The computing time required by this technique can be greatly reduced by using a homogeneous half‐space as the starting model so that the Jacobian matrix of partial derivatives can be calculated analytically. A quasi‐Newton updating method is then used to estimate the partial derivatives for subsequent iterations. This inversion technique has been tested on synthetic and field data where a satisfactory model is obtained using a modest amount of computer time. On an 80486DX2/66 microcomputer, it takes about 20 minutes to invert the data from a 7 by 7 electrode survey grid. using the techniques described below, 3D resistivity surveys and data inversion can be carried out using commercially available field equipment and an inexpensive microcomputer.
-
-
-
Velocity sensitivity in transversely isotropic media1
Authors C.H. Chapman and D.E. MillerAbstractWe consider the problem of determining and predicting how the wave speeds in particular directions for a transversely isotropic (TI) medium depend on particular combinations of the density‐normalized moduli Aij. The expressions for the qP and qSV velocities are known to depend on four moduli. Normally, we can only determine three independent parameters from qP data, or two from qSZ data, as the others have much lower sensitivity. The resolvable parameters are conveniently described by axial and off‐axis parameters: for qP rays, P0°= A11, P90°= A33 and P45°=(A11+ A33)/4 + (A13+2A55)/2; and for qSV rays, S0°= S90°=A55 and S 45°= (A11+ A 33)/4‐ A13/2. These parameters control the magnitude of the squared‐velocities on the axes and at approximately 45°. For an arbitrary TI medium, if the medium is perturbed in a way that preserves a particular parameter, then slowness points in the associated direction and mode witl be approximately preserved in the new medium. we refer to these parameters as ‘push‐pins’, i.e. if a parameter is fixed, the associated part of the slowness surface is pinned in place.
Because, these five push‐pins only contain four independent moduli, we can only fix at most three push‐pins. Perturbing one of the other parameters inevitably perturbs the other. Numerical results illustrating the linkage between two push‐pins, when three are fixed, are presented.
So‐called anomalous TI media occur when the roles of the qP and qSV waves are reversed: in some directions the faster ray has transverse polarization. That, in turn, requires anomalous velocities at the push‐pins, i.e. S0° > P0°, S45° > P45° and/or S90° > P90° (equivalent to the usual anomalous conditions A11 < A55, < 0 and/or A33 < A55). In the Appendix, we confirm that anomalous sensitivities of the velocities at the five push‐pins only occur in such media, although the push‐pins still apply if interpreted appropriately. Truly anomalous sensitivities, in which push‐pins play no role, only occur in media near the boundary between normal and anomalous.
-
-
-
On the downward continuation of electromagnetic fields1
More LessAbstractA concise derivation is given for the downward continuation of the tangential fields on the surface to yield expressions for the fields at a specified depth. A homogeneous slab region is assumed for the analysis.
-
Volumes & issues
-
Volume 72 (2023 - 2024)
-
Volume 71 (2022 - 2023)
-
Volume 70 (2021 - 2022)
-
Volume 69 (2021)
-
Volume 68 (2020)
-
Volume 67 (2019)
-
Volume 66 (2018)
-
Volume 65 (2017)
-
Volume 64 (2015 - 2016)
-
Volume 63 (2015)
-
Volume 62 (2014)
-
Volume 61 (2013)
-
Volume 60 (2012)
-
Volume 59 (2011)
-
Volume 58 (2010)
-
Volume 57 (2009)
-
Volume 56 (2008)
-
Volume 55 (2007)
-
Volume 54 (2006)
-
Volume 53 (2005)
-
Volume 52 (2004)
-
Volume 51 (2003)
-
Volume 50 (2002)
-
Volume 49 (2001)
-
Volume 48 (2000)
-
Volume 47 (1999)
-
Volume 46 (1998)
-
Volume 45 (1997)
-
Volume 44 (1996)
-
Volume 43 (1995)
-
Volume 42 (1994)
-
Volume 41 (1993)
-
Volume 40 (1992)
-
Volume 39 (1991)
-
Volume 38 (1990)
-
Volume 37 (1989)
-
Volume 36 (1988)
-
Volume 35 (1987)
-
Volume 34 (1986)
-
Volume 33 (1985)
-
Volume 32 (1984)
-
Volume 31 (1983)
-
Volume 30 (1982)
-
Volume 29 (1981)
-
Volume 28 (1980)
-
Volume 27 (1979)
-
Volume 26 (1978)
-
Volume 25 (1977)
-
Volume 24 (1976)
-
Volume 23 (1975)
-
Volume 22 (1974)
-
Volume 21 (1973)
-
Volume 20 (1972)
-
Volume 19 (1971)
-
Volume 18 (1970)
-
Volume 17 (1969)
-
Volume 16 (1968)
-
Volume 15 (1967)
-
Volume 14 (1966)
-
Volume 13 (1965)
-
Volume 12 (1964)
-
Volume 11 (1963)
-
Volume 10 (1962)
-
Volume 9 (1961)
-
Volume 8 (1960)
-
Volume 7 (1959)
-
Volume 6 (1958)
-
Volume 5 (1957)
-
Volume 4 (1956)
-
Volume 3 (1955)
-
Volume 2 (1954)
-
Volume 1 (1953)