AVO modeling and the locally converted shear wave

Geophysics ◽  
1994 ◽  
Vol 59 (8) ◽  
pp. 1237-1248 ◽  
Author(s):  
James L. Simmons ◽  
Milo M. Backus

The locally converted shear wave is often neglected in ray‐trace modeling when reproduction of the AVO response of potential hydrocarbon reservoirs is attempted. Primaries‐only ray‐trace modeling in which the Zoeppritz equations describe the reflection amplitudes is most common. The locally converted shear waves, however, often have a first‐order effect on the seismic response. This fact does not appear to be widely recognized, or else the implications are not well understood. Primaries‐only Zoeppritz modeling can be very misleading. Interference between the converted waves and the primary reflections from the base of the layers becomes increasingly important as layer thicknesses decrease. This interference often produces a seismogram that is very different from one produced under the primaries‐only Zoeppritz assumption. For primaries‐only modeling of thin layers, synthetic seismograms obtained by use of a linearized approximation to the Zoeppritz equations to describe the reflection coefficients are more accurate than those obtained by use of the exact Zoeppritz reflection coefficients. A real‐data example consisting of an assemblage of very thin layers has recently been discussed in the literature. Inferences as to the true earth properties based on the predicted amplitude variation with offset are in error because the primaries‐only assumption is invalid. For one of the models, primaries‐only modeling predicts an amplitude increase of approximately a factor of three from the near trace to the far trace. Reflectivity modeling predicts an amplitude decrease with offset. The O’Doherty‐Anstey effect suggests that transmission loss for primary reflections should not be included in normal‐incidence synthetic seismograms if the short‐period reverberations are not also included. The same principle holds for prestack modeling. Similarly, the Zoeppritz equations should not be used for synthetic seismograms without including the locally converted shear wave.

Geophysics ◽  
1996 ◽  
Vol 61 (6) ◽  
pp. 1575-1588 ◽  
Author(s):  
James L. Simmons ◽  
Milo M. Backus

A practical approach to linear prestack seismic inversion in the context of a locally 1-D earth is employed to use amplitude variation with offset (AVO) information for the direct detection in hydrocarbons. The inversion is based on the three‐term linearized approximation to the Zoeppritz equations. The normal‐incidence compressional‐wave reflection coefficient [Formula: see text] models the background reflectivity in the absence of hydrocarbons and incorporates the mudrock curve and Gardner’s equation. Prediction‐error parameters, [Formula: see text] and [Formula: see text], represent perturbations in the normal‐incidence shear‐wave reflection coefficient and the density contribution to the normal incidence reflectivity, respectively, from that predicted by the mudrock curve and Gardner’s equation. This prediction‐error approach can detect hydrocarbons in the absence of an overall increase in AVO, and in the absence of bright spots, as expected in theory. Linear inversion is applied to a portion of a young, Tertiary, shallow‐marine data set that contains known hydrocarbon accumulations. Prestack data are in the form of angle stack, or constant offset‐to‐depth ratio, gathers. Prestack synthetic seismograms are obtained by primaries‐only ray tracing using the linearized approximation to the Zoeppritz equations to model the reflection amplitudes. Where the a priori assumptions hold, the data are reproduced with a single parameter [Formula: see text]. Hydrocarbons are detected as low impedance relative to the surrounding shales and the downdip brine‐filled reservoir on [Formula: see text], also as positive perturbations (opposite polarity relative to [Formula: see text]) on [Formula: see text] and [Formula: see text]. The maximum perturbation in [Formula: see text] from the normal‐incidence shear‐wave reflection coefficient predicted by the a priori assumptions is 0.08. Hydrocarbon detection is achieved, although the overall seismic response of a gas‐filled thin layer shows a decrease in amplitude with offset (angle). The angle‐stack data (70 prestack ensembles, 0.504–1.936 s time range) are reproduced with a data residual that is 7 dB down. Reflectivity‐based prestack seismograms properly model a gas/water contact as a strong increase in AVO and a gas‐filled thin layer as a decrease in AVO.


Geophysics ◽  
1981 ◽  
Vol 46 (7) ◽  
pp. 972-983 ◽  
Author(s):  
Håvar Gjøystdal ◽  
Bjørn Ursin

When reflection data are available from a grid of crossing seismic lines, it is possible to construct normal incidence time maps from interpreted stacked sections and then apply three‐dimensional (3-D) ray‐tracing techniques following the normal‐incidence raypaths down to the various reflectors. The main disadvantage of this well‐known “time map migration” procedure is that interval velocities must be known a priori, and they must be estimated in advance by some approximate method. A technique is presented here which combines the above procedure with an inversion algorithm, providing direct calculations of interval velocities from the additional use of nonzero offset traveltime observations. A generalized linear inversion scheme is used, making possible a complete calculation of interval velocities and reflection interfaces, the latter represented by bicubic spline functions. To test the method in practice, we have applied it to (1) synthetic data generated from a constructed model, and (2) real data obtained from marine seismic sections. In the latter case, velocities and reflector depths obtained were compared to those obtained directly from a well log in the area. These results show a reasonably good resolution for layers that are not too deep relative to the shot/receiver offsets used. For deep and/or thin layers, the results are not satisfactory. This indicates the general limitation of seismic reflection data to resolve interval velocity, even in the presence of horizontally layered structure.


Geophysics ◽  
2021 ◽  
pp. 1-50
Author(s):  
German Garabito ◽  
José Silas dos Santos Silva ◽  
Williams Lima

In land seismic data processing, the prestack time migration (PSTM) image remains the standard imaging output, but a reliable migrated image of the subsurface depends on the accuracy of the migration velocity model. We have adopted two new algorithms for time-domain migration velocity analysis based on wavefield attributes of the common-reflection-surface (CRS) stack method. These attributes, extracted from multicoverage data, were successfully applied to build the velocity model in the depth domain through tomographic inversion of the normal-incidence-point (NIP) wave. However, there is no practical and reliable method for determining an accurate and geologically consistent time-migration velocity model from these CRS attributes. We introduce an interactive method to determine the migration velocity model in the time domain based on the application of NIP wave attributes and the CRS stacking operator for diffractions, to generate synthetic diffractions on the reflection events of the zero-offset (ZO) CRS stacked section. In the ZO data with diffractions, the poststack time migration (post-STM) is applied with a set of constant velocities, and the migration velocities are then selected through a focusing analysis of the simulated diffractions. We also introduce an algorithm to automatically calculate the migration velocity model from the CRS attributes picked for the main reflection events in the ZO data. We determine the precision of our diffraction focusing velocity analysis and the automatic velocity calculation algorithms using two synthetic models. We also applied them to real 2D land data with low quality and low fold to estimate the time-domain migration velocity model. The velocity models obtained through our methods were validated by applying them in the Kirchhoff PSTM of real data, in which the velocity model from the diffraction focusing analysis provided significant improvements in the quality of the migrated image compared to the legacy image and to the migrated image obtained using the automatically calculated velocity model.


Geophysics ◽  
1978 ◽  
Vol 43 (6) ◽  
pp. 1099-1110 ◽  
Author(s):  
Albert C. Reynolds

Many finite difference models in use for generating synthetic seismograms produce unwanted reflections from the edges of the model due to the use of Dirichlet or Neumann boundary conditions. In this paper we develop boundary conditions which greatly reduce this edge reflection. A reflection coefficient analysis is given which indicates that, for the specified boundary conditions, smaller reflection coefficients than those obtained for Dirichlet or Neumann boundary conditions are obtained. Numerical calculations support this conclusion.


2021 ◽  
Author(s):  
Yi Luan ◽  
Rui Ding ◽  
Wenshen Gu ◽  
Xiaofan Zhang ◽  
Xinliang Chen ◽  
...  

Abstract Since the end of 2019, the COVID-19 epidemic has swept the world. With the widespread spread of the COVID-19 and the continuous emergence of mutated strains, the situation for the prevention and control of the COVID-19 epidemic remains severe. On May 21, 2021, Guangzhou City, Guangdong Province, notified the discovery of a new locally confirmed case. Guangzhou became the first city in mainland China to compete with the delta mutant strain. As a local hospital with strong nucleic acid detection capabilities, Sun Yat-sen University Sun Yat-sen Memorial Hospital took the lead in launching the construction and deployment of the Mobile Shelter Laboratories and large-scale screening work in Foshan and Zhanjiang, Guangdong Province. Through summarizing "practical" experience, observation and comparison data analysis, we use real data to verify a feasible solution for rapid expansion of detection capabilities in a short period of time. We hope that these experiences will have certain reference value for other countries or regions, especially the underdeveloped areas of medical and health care.


Energies ◽  
2020 ◽  
Vol 13 (12) ◽  
pp. 3074 ◽  
Author(s):  
Shulin Pan ◽  
Ke Yan ◽  
Haiqiang Lan ◽  
José Badal ◽  
Ziyu Qin

Conventional sparse spike deconvolution algorithms that are based on the iterative shrinkage-thresholding algorithm (ISTA) are widely used. The aim of this type of algorithm is to obtain accurate seismic wavelets. When this is not fulfilled, the processing stops being optimum. Using a recurrent neural network (RNN) as deep learning method and applying backpropagation to ISTA, we have developed an RNN-like ISTA as an alternative sparse spike deconvolution algorithm. The algorithm is tested with both synthetic and real seismic data. The algorithm first builds a training dataset from existing well-logs seismic data and then extracts wavelets from those seismic data for further processing. Based on the extracted wavelets, the new method uses ISTA to calculate the reflection coefficients. Next, inspired by the backpropagation through time (BPTT) algorithm, backward error correction is performed on the wavelets while using the errors between the calculated reflection coefficients and the reflection coefficients corresponding to the training dataset. Finally, after performing backward correction over multiple iterations, a set of acceptable seismic wavelets is obtained, which is then used to deduce the sequence of reflection coefficients of the real data. The new algorithm improves the accuracy of the deconvolution results by reducing the effect of wrong seismic wavelets that are given by conventional ISTA. In this study, we account for the mechanism and the derivation of the proposed algorithm, and verify its effectiveness through experimentation using theoretical and real data.


2016 ◽  
Vol 41 (1) ◽  
pp. 43-53 ◽  
Author(s):  
Ying-Chun Chang ◽  
Ho-Chih Cheng ◽  
Min-Chie Chiu ◽  
Yuan-Hung Chien

Abstract Research on plenums partitioned with multiple baffles in the industrial field has been exhaustive. Most researchers have explored noise reduction effects based on the transfer matrix method and the boundary element method. However, maximum noise reduction of a plenum within a constrained space, which frequently occurs in engineering problems, has been neglected. Therefore, the optimum design of multi-chamber plenums becomes essential. In this paper, two kinds of multi-chamber plenums (Case I: a two-chamber plenum that is partitioned with a centre-opening baffle; Case II: a three-chamber plenum that is partitioned with two centre-opening baffles) within a fixed space are assessed. In order to speed up the assessment of optimal plenums hybridized with multiple partitioned baffles, a simplified objective function (OBJ) is established by linking the boundary element model (BEM, developed using SYSNOISE) with a polynomial neural network fit with a series of real data – input design data (baffle dimensions) and output data approximated by BEM data in advance. To assess optimal plenums, a genetic algorithm (GA) is applied. The results reveal that the maximum value of the transmission loss (TL) can be improved at the desired frequencies. Consequently, the algorithm proposed in this study can provide an efficient way to develop optimal multi-chamber plenums for industry.


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.


Sign in / Sign up

Export Citation Format

Share Document