Impedance‐type approximations of the P–P elastic reflection coefficient: Modeling and AVO inversion

Geophysics ◽  
2004 ◽  
Vol 69 (2) ◽  
pp. 592-598 ◽  
Author(s):  
Lúcio T. Santos ◽  
Martin Tygel

The normal‐incidence elastic compressional reflection coefficient admits an exact, simple expression in terms of the acoustic impedance, namely the product of the P‐wave velocity and density, at both sides of an interface. With slight modifications a similar expression can, also exactly, express the oblique‐incidence acoustic reflection coefficient. A severe limitation on the use of these two reflection coefficients in analyzing seismic reflection data is that they provide no information on shear‐wave velocities that refer to the interface. We address the natural question of whether a suitable impedance concept can be introduced for which arbitrary P–P reflection coefficients can be expressed in a form analogous to their acoustic counterparts. Although no closed‐form exact solution exists, our analysis provides a general framework for which, under suitable restrictions of the medium parameters, possible impedance functions can be derived. In particular, the well‐established concept of elastic impedance and the recently introduced concept of reflection impedance can be better understood. Concerning these two impedances, we examine their potential for modeling and for estimating the AVO indicators of intercept and gradient. For typical synthetic examples, we show that the reflection impedance formulation provides consistently better results than those obtained using the elastic impedance.

Geophysics ◽  
2006 ◽  
Vol 71 (5) ◽  
pp. E49-E55 ◽  
Author(s):  
Jonathan E. Downton ◽  
Charles Ursenbach

Contrary to popular belief, a linearized approximation of the Zoeppritz equations may be used to estimate the reflection coefficient for angles of incidence up to and beyond the critical angle. These supercritical reflection coefficients are complex, implying a phase variation with offset in addition to amplitude variation with offset (AVO). This linearized approximation is then used as the basis for an AVO waveform inversion. By incorporating this new approximation, wider offset and angle data may be incorporated in the AVO inversion, helping to stabilize the problem and leading to more accurate estimates of reflectivity, including density reflectivity.


2016 ◽  
Vol 4 (4) ◽  
pp. T613-T625 ◽  
Author(s):  
Qizhen Du ◽  
Bo Zhang ◽  
Xianjun Meng ◽  
Chengfeng Guo ◽  
Gang Chen ◽  
...  

Three-term amplitude-variation with offset (AVO) inversion generally suffers from instability when there is limited prior geologic or petrophysical constraints. Two-term AVO inversion shows higher instability compared with three-term AVO inversion. However, density, which is important in the fluid-type estimation, cannot be recovered from two-term AVO inversion. To reliably predict the P- and S-waves and density, we have developed a robust two-step joint PP- and PS-wave three-term AVO-inversion method. Our inversion workflow consists of two steps. The first step is to estimate the P- and S-wave reflectivities using Stewart’s joint two-term PP- and PS-AVO inversion. The second step is to treat the P-wave reflectivity obtained from the first step as the prior constraint to remove the P-wave velocity related-term from the three-term Aki-Richards PP-wave approximated reflection coefficient equation, and then the reduced PP-wave reflection coefficient equation is combined with the PS-wave reflection coefficient equation to estimate the S-wave and density reflectivities. We determined the effectiveness of our method by first applying it to synthetic models and then to field data. We also analyzed the condition number of the coefficient matrix to illustrate the stability of the proposed method. The estimated results using proposed method are superior to those obtained from three-term AVO inversion.


Geophysics ◽  
2011 ◽  
Vol 76 (1) ◽  
pp. SA19-SA34 ◽  
Author(s):  
Shohei Minato ◽  
Toshifumi Matsuoka ◽  
Takeshi Tsuji ◽  
Deyan Draganov ◽  
Jürg Hunziker ◽  
...  

Crosswell reflection method is a high-resolution seismic imaging method that uses recordings between boreholes. The need for downhole sources is a restrictive factor in its application, for example, to time-lapse surveys. An alternative is to use surface sources in combination with seismic interferometry. Seismic interferometry (SI) could retrieve the reflection response at one of the boreholes as if from a source inside the other borehole. We investigate the applicability of SI for the retrieval of the reflection response between two boreholes using numerically modeled field data. We compare two SI approaches — crosscorrelation (CC) and multidimensional deconvolution (MDD). SI by MDD is less sensitive to underillumination from the source distribution, but requires inversion of the recordings at one of the receiver arrays from all the available sources. We find that the inversion problem is ill-posed, and propose to stabilize it using singular-value decomposition. The results show that the reflections from deep boundaries are retrieved very well using both the CC and MDD methods. Furthermore, the MDD results exhibit more realistic amplitudes than those from the CC method for downgoing reflections from shallow boundaries. We find that the results retrieved from the application of both methods to field data agree well with crosswell seismic-reflection data using borehole sources and with the logged P-wave velocity.


Geophysics ◽  
2016 ◽  
Vol 81 (3) ◽  
pp. R57-R74 ◽  
Author(s):  
Santi Kumar Ghosh ◽  
Animesh Mandal

Because seismic reflection data are band limited, acoustic impedance profiles derived from them are nonunique. The conventional inversion methods counter the nonuniqueness either by stabilizing the answer with respect to an initial model or by imposing mathematical constraints such as sparsity of the reflection coefficients. By making a nominal assumption of an earth model locally consisting of a stack of homogeneous and horizontal layers, we have formulated a set of linear equations in which the reflection coefficients are the unknowns and the recursively integrated seismic trace constitute the data. Drawing only on first principles, the Zoeppritz equation in this case, the approach makes a frontal assault on the problem of reconstructing reflection coefficients from band-limited data. The local layer-cake assumption and the strategy of seeking a singular value decomposition solution of the linear equations counter the nonuniqueness, provided that the objective is to reconstruct a smooth version of the impedance profile that includes only its crude structures. Tests on synthetic data generated from elementary models and from measured logs of acoustic impedance demonstrated the efficacy of the method, even when a significant amount of noise was added to the data. The emergence of consistent estimates of impedance, approximating the original impedance, from synthetic data generated for several frequency bands has inspired our confidence in the method. The other attractive outputs of the method are as follows: (1) an accurate estimate of the impedance mean, (2) an accurate reconstruction of the direct-current (DC) frequency of the reflectivity, and (3) an acceptable reconstruction of the broad outline of the original impedance profile. These outputs can serve as constraints for either more refined inversions or geologic interpretations. Beginning from the restriction of band-limited data, we have devised a method that neither requires a starting input model nor imposes mathematical constraints on the earth reflectivity and still yielded significant and relevant geologic information.


2020 ◽  
Author(s):  
Irene DeFelipe ◽  
Juan Alcalde ◽  
Monika Ivandic ◽  
David Martí ◽  
Mario Ruiz ◽  
...  

Abstract. Seismic reflection data (normal incidence and wide-angle) are unique assets for Solid Earth Science as they provide critical information about the physical properties and structure of the lithosphere, as well as about the shallow subsurface for exploration purposes. The resolution of these seismic data is highly appreciated, however they are logistically complex and expensive to acquire and their geographical coverage is limited. Therefore, it is essential to make the most of the data that has already been acquired. The collation and dissemination of seismic open access data is then key to promote accurate and innovative research and to enhance new interpretations of legacy data. This work presents the Seismic DAta REpository (SeisDARE), which is, to our knowledge, one of the first comprehensive open access online databases that stores seismic data registered with a permanent identifier (DOI). The datasets included here are openly accessible online and guarantee the FAIR (Findable, Accessible, Interoperable, Reusable) principles of data management, granting the inclusion of each dataset into a statistics referencing database so its impact can be measured. SeisDARE includes seismic data acquired in the last four decades in the Iberian Peninsula and Morocco. These areas have attracted the attention of international researchers in the fields of geology and geophysics due to the exceptional outcrops of the Variscan and Alpine orogens and wide foreland basins; the crustal structure of the offshore margins that resulted from a complex plate kinematic evolution; and the vast quantities of natural resources contained within. This database has been built thanks to a network of national and international institutions, promoting a multidisciplinary research, and is open for international data exchange and collaborations. As part of this international collaboration, and as a model for inclusion of other global seismic datasets, SeisDARE also hosts seismic data acquired in Hardeman County, Texas (USA), within the COCORP project (Consortium for Continental Reflection Profiling). SeisDARE aims to make easily accessible old and recently acquired seismic data and to establish a framework for future seismic data management plans. The SeisDARE is freely available at https://digital.csic.es/handle/10261/101879, bringing endless research and teaching opportunities to the scientific, industrial and educational communities.


Geophysics ◽  
2020 ◽  
Vol 85 (3) ◽  
pp. R251-R262 ◽  
Author(s):  
Ligia Elena Jaimes-Osorio ◽  
Alison Malcolm ◽  
Ali Gholami

Conventional amplitude variation with offset (AVO) inversion analysis uses the Zoeppritz equations, which are based on a plane-wave approximation. However, because real seismic data are created by point sources, wave reflections are better modeled by spherical waves than by plane waves. Indeed, spherical reflection coefficients deviate from planar reflection coefficients near the critical and postcritical angles, which implies that the Zoeppritz equations are not applicable for angles close to critical reflection in AVO analysis. Elastic finite-difference simulations provide a solution to the limitations of the Zoeppritz approximation because they can handle near- and postcritical reflections. We have used a coupled acoustic-elastic local solver that approximates the wavefield with high accuracy within a locally perturbed elastic subdomain of the acoustic full domain. Using this acoustic-elastic local solver, the local wavefield generation and inversion are much faster than performing a full-domain elastic inversion. We use this technique to model wavefields and to demonstrate that the amplitude from within the local domain can be used as a constraint in the inversion to recover elastic material properties. Then, we focus on understanding how much the amplitude and phase contribute to the reconstruction accuracy of the elastic material parameters ([Formula: see text], [Formula: see text], and [Formula: see text]). Our results suggest that the combination of amplitude and phase in the inversion helps with the convergence. Finally, we analyze elastic parameter trade-offs in AVO inversion, from which we find that to recover accurate P-wave velocities we should invert for [Formula: see text] and [Formula: see text] simultaneously with fixed density.


Geophysics ◽  
2011 ◽  
Vol 76 (2) ◽  
pp. B55-B70 ◽  
Author(s):  
E. M. Takam Takougang ◽  
A. J. Calvert

To obtain a higher resolution quantitative P-wave velocity model, 2D waveform tomography was applied to seismic reflection data from the Queen Charlotte sedimentary basin off the west coast of Canada. The forward modeling and inversion were implemented in the frequency domain using the visco-acoustic wave equation. Field data preconditioning consisted of f-k filtering, 2D amplitude scaling, shot-to-shot amplitude balancing, and time windowing. The field data were inverted between 7 and 13.66 Hz, with attenuation introduced for frequencies ≥ 10.5 Hz to improve the final velocity model; two different approaches to sampling the frequencies were evaluated. The limited maximum offset of the marine data (3770 m) and the relatively high starting frequency (7 Hz) were the main challenges encountered during the inversion. An inversion strategy that successively recovered shallow-to-deep structures was designed to mitigate these issues. The inclusion of later arrivals in the waveform tomography resulted in a velocity model that extends to a depth of approximately 1200 m, twice the maximum depth of ray coverage in the ray-based tomography. Overall, there is a good agreement between the velocity model and a sonic log from a well on the seismic line, as well as between modeled shot gathers and field data. Anomalous zones of low velocity in the model correspond to previously identified faults or their upward continuation into the shallow Pliocene section where they are not readily identifiable in the conventional migration.


Geophysics ◽  
1995 ◽  
Vol 60 (4) ◽  
pp. 939-954 ◽  
Author(s):  
Subhashis Mallick

I cast the inversion of amplitude‐variation‐with‐offset (AVO) data into the framework of Bayesian statistics. Under such a framework, the model parameters and the physics of the forward problem are used to generate synthetic data. These synthetic data are then matched with the observed data to obtain an a‐posteriori probability density (PPD) function in the model space. The genetic algorithm (GA) uses a directed random search technique to estimate the shape of the PPD. Unlike the classical inversion methods, GA does not depend upon the choice of an initial model and is well suited for the AVO inversion. For the single‐layer AVO inversion where the amplitudes from a single reflection event are inverted, GA estimates the normal incidence reflection coefficient [Formula: see text] and the contrast of the Poisson’s ratio (Δσ) to reasonable accuracy, even when the signal‐to‐noise ratio is poor. Comparisons of single‐layer amplitude inversion using synthetic data demonstrate that GA inversion obtains more accurate results than does the least‐squares fit to the approximate reflection coefficients as is usually practiced in the industry. In the multilayer AVO waveform inversion, all or a part of the prestack data are inverted. Inversion of this type is nonunique for the estimation of the absolute values of velocities, Poisson’s ratios, and densities. However, by applying simplified approximations to the P‐wave reflection coefficient, I verify that [Formula: see text], the contrast in the acoustic impedance (ΔA), and the gradient in the reflection coefficient (G), can be estimated from such an inversion. From the GA estimated values of [Formula: see text], ΔA, and G, and from reliable estimates of velocity and Poisson’s ratio at the start time of the input data, an inverted model can be generated. I apply this procedure to marine data and demonstrate that the the synthetics computed from such an inverted model match the input data to reasonable accuracy. Comparison of the log data from a nearby well shows that the GA inversion obtains both the low‐ and the high‐frequency trends (within the bandwidth of seismic resolution) of the P‐wave acoustic impedance. In addition to P‐wave acoustic impedance, GA also obtains an estimate of the Poisson’s ratio, an extremely important parameter for the direct detection of hydrocarbons from seismic data.


Sign in / Sign up

Export Citation Format

Share Document