- Home
- A-Z Publications
- Geophysical Prospecting
- Previous Issues
- Volume 41, Issue 2, 1993
Geophysical Prospecting - Volume 41, Issue 2, 1993
Volume 41, Issue 2, 1993
-
-
ASPECTS OF 1D SEISMIC MODELLING USING THE GOUPILLAUD PRINCIPLE1
Authors EVERT SLOB and ANTON ZIOLKOWSK1AbstractA reflection response function for a 1D discretized earth model can be obtained using ray‐theory and Z‐transforms with the Goupillaud model. This is usually done by taking the source function as a plane wave impinging normally on the layered earth. Two important problems have been tackled with this basic idea. The first, extraction of the source wavelet, and the second, a description of the free‐surface related problems.
In the Goupillaud model, the one‐way traveltime in each layer is taken to be the same time interval At, which is also the time unit for the Z‐transform. The two‐way traveltime in any layer is 2Δt, corresponding to a multiplication by Z2. The reflection impulse response therefore contains only even powers of Z. The convolution of the reflection response with the wavelet yields a seismogram whose Z‐transform contains both odd and even powers of Z. However, even though the seismogram contains more coefficients than unknowns, the wavelet cannot be extracted, because the coefficients are not independent: later coefficients are functions of earlier ones, which does not make sense physically. To overcome this physical problem for the reflection seismogram, the two‐way traveltime through the layer should be Δt. It is then impossible to extract the wavelet, as there are fewer coefficients in the seismogram than unknowns.
Szaraniec has proposed a modification to the Goupillaud model, known as the odd‐depth model, that includes the free surface and a top layer whose two‐way traveltime Δt is half the two‐way traveltime 2Δt of all the other layers. Using what Szaraniec calls the fundamental identity of the odd‐depth model, it is possible to extract the source wavelet from the seismogram. We show that this fundamental identity holds only if reflection coefficients of deeper interfaces are functions of the reflection coefficients of shallower interfaces; that is, for extremely improbable geologies.
Neither of these approaches offers a solution to the deconvolution problem. It is better to obtain the source signature from measurements in the field. Only Szaraniec's model offers the possibility of tackling the problem of the free surface but because of an inherent flaw in the model, it fails to address the problem.
-
-
-
VISCOELASTIC SEISMIC RESPONSES OF 2D RESERVOIR MODELS1
Authors I. B. KANG and GEORGE A. McMECHANAbstractThe use of relaxation mechanisms has recently made it possible to simulate viscoelastic (Q) effects accurately in time‐domain numerical computations of seismic responses. As a result, seismograms may now be synthesized for models with arbitrary spatial variations in compressional‐ and shear‐wave quality factors (Q9, and Qs, as well as in density (ρ) and compressional‐ and shear‐wave velocities (Vp, and Vs).
Reflections produced by Q contrasts alone may have amplitudes as large as those produced by velocity contrasts. Q effects, including their interaction with Vp, Vs and p, contribute significantly to the seismic response of reservoirs. For band‐limited data at typical seismic frequencies, the effects of Q on reflectivity and attenuation are more visible than those on dispersion.
Synthetic examples include practical applications to reservoir exploration, evaluation and monitoring. Q effects are clearly visible in both surface and offset vertical seismic profile data. Thus, AVO analyses that neglect Q may produce erroneous conclusions.
-
-
-
BOREHOLE STONELEY WAVE PROPAGATION ACROSS PERMEABLE STRUCTURES1
Authors X. M. TANG and C. H. CHENGAbstractThis study investigates the propagation of borehole Stoneley waves across permeable structures. By modelling the structure as a zone intersecting the borehole, a simple 1D theory is formulated to treat the interaction of the Stoneley wave with the structure. This is possible because the Stoneley wave is a guided wave, with no geometric spreading as it propagates along the borehole. The interaction occurs because the zone and the surrounding formation possess different Stoneley wavenumbers. Given appropriate representations of the wavenum‐ber, the theory can be applied to treat a variety of structures, including a fluid‐filled fracture. Of special interest are the cases of permeable porous zones and fracture zones. The results show that, while Stoneley wave reflections are generated, strong Stoneley wave attenuation is produced across a very permeable zone. This result is particularly important in explaining the observed strong Stoneley wave attenuation at major fractures where it has been difficult to explain the attenuation in terms of the single planar fracture theory. In addition, by using a simple and sufficiently accurate theory to model the effects of the permeable zone, a fast and efficient method is developed to characterize the fluid transport properties of a permeable fracture zone.
-
-
-
A COMPARISON OF LABORATORY AND FIELD MEASUREMENTS OF P‐WAVE ANISOTROPY1
Authors M. S. SAMS, M. H. WORTHINGTON, M. S. KING and M. SHAMS KHANSHIRAbstractP‐wave velocity anisotropy determined from a cross‐hole survey at the Imperial College borehole test site compares favourably with that measured in the laboratory on core from the holes. The holes penetrate a layered sequence of sandstones, shales and carbonates of the Namurian Upper Limestone Group. The Laboratory measurements of the vertical and horizontal velocities of core samples indicate that the shales exhibit P‐wave anisotropies of over 20% but that the sandstones and limestones are only slightly anisotropic. These discrete measurements have been used in combination with wireline data to produce a log of P‐wave anisotropy. Including the anisotropic information vastly improves the match between observed and synthetic traveltimes from the cross‐hole data set. This implies that there is little frequency dependence of intrinsic P‐wave anisotropy.
Inversion of the cross‐hole traveltimes highlights the need for good angular coverage in order to resolve the anisotropy parameters. The observed P‐wave anisotropy of the field data is due to the combined effect of sedimentary layering and the intrinsic anisotropy of the rocks. The intrinsic anisotropy is found to be the dominant factor.
-
-
-
EXAMPLES OF RESOLUTION IMPROVEMENT IN GEOELECTRICAL SOUNDINGS APPLIED TO GROUNDWATER INVESTIGATIONS1
More LessAbstractOne simulation and two field examples from New Jersey illustrate resolution improvement in geoelectrical soundings applied to groundwater exploration. Layered‐earth parameter resolution is derived from data obtained with the commonly used methods of resistivity, induced polarization (IP) and transient electromagnetic (TEM) soundings. Resolution improvement is achieved by simultaneous inversion of two or more data sets and by constraining parameters of the inverse problem.
A quantitative analysis showing the contribution of IP data to the resolution of geo‐electric sections is presented. Simultaneous inversion of simple IP data with conventional resistivity and resistivity‐TEM data sets resulted in improved parameter resolution. IP data improved resolution in three ways: (1) by decoupling correlated layered‐earth parameters, (2) by adding information to a geological interpretation about a second physical property, namely chargeability, and (3) by increasing the electrical information available.
-
-
-
RADON DEPTH MIGRATION1
Authors S. T. HILDEBRAND and R. J. CARROLLAbstractA depth migration method is presented that uses Radon‐transformed common‐source seismograms as input. It is shown that the Radon depth migration method can be extended to spatially varying velocity depth models by using asymptotic ray theory (ART) to construct wavefield continuation operators. These operators downward continue an incident receiver‐array plane wave and an assumed point‐source wavefield into the subsurface. The migration velocity model is constrained to have longer characteristic wavelengths than the dominant source wavelength such that the ART approximations for the continuation operators are valid.
This method is used successfully to migrate two synthetic data examples:
- 1 a point diffractor, and
- 2 a dipping layer and syncline interface model.
It is shown that the Radon migration method has a computational advantage over the standard Kirchhoff migration method in that fewer rays are computed in a main memory implementation.
-
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)