Poststack impedance inversion considering the diffractive component of the wavefield

Geophysics ◽  
2020 ◽  
Vol 85 (1) ◽  
pp. R11-R28 ◽  
Author(s):  
Kun Xiang ◽  
Evgeny Landa

Seismic diffraction waveform energy contains important information about small-scale subsurface elements, and it is complementary to specular reflection information about subsurface properties. Diffraction imaging has been used for fault, pinchout, and fracture detection. Very little research, however, has been carried out taking diffraction into account in the impedance inversion. Usually, in the standard inversion scheme, the input is the migrated data and the assumption is taken that the diffraction energy is optimally focused. This assumption is true only for a perfectly known velocity model and accurate true amplitude migration algorithm, which are rare in practice. We have developed a new approach for impedance inversion, which takes into account diffractive components of the total wavefield and uses the unmigrated input data. Forward modeling, designed for impedance inversion, includes the classical specular reflection plus asymptotic diffraction modeling schemes. The output model is composed of impedance perturbation and the low-frequency model. The impedance perturbation is estimated using the Bayesian approach and remapped to the migrated domain by the kinematic ray tracing. Our method is demonstrated using synthetic and field data in comparison with the standard inversion. Results indicate that inversion with taking into account diffraction can improve the acoustic impedance prediction in the vicinity of local reflector discontinuities.

Geophysics ◽  
2009 ◽  
Vol 74 (5) ◽  
pp. R59-R67 ◽  
Author(s):  
Igor B. Morozov ◽  
Jinfeng Ma

The seismic-impedance inversion problem is underconstrained inherently and does not allow the use of rigorous joint inversion. In the absence of a true inverse, a reliable solution free from subjective parameters can be obtained by defining a set of physical constraints that should be satisfied by the resulting images. A method for constructing synthetic logs is proposed that explicitly and accurately satisfies (1) the convolutional equation, (2) time-depth constraints of the seismic data, (3) a background low-frequency model from logs or seismic/geologic interpretation, and (4) spectral amplitudes and geostatistical information from spatially interpolated well logs. The resulting synthetic log sections or volumes are interpretable in standard ways. Unlike broadly used joint-inversion algorithms, the method contains no subjectively selected user parameters, utilizes the log data more completely, and assesses intermediate results. The procedure is simple and tolerant to noise, and it leads to higher-resolution images. Separating the seismic and subseismic frequency bands also simplifies data processing for acoustic-impedance (AI) inversion. For example, zero-phase deconvolution and true-amplitude processing of seismic data are not required and are included automatically in this method. The approach is applicable to 2D and 3D data sets and to multiple pre- and poststack seismic attributes. It has been tested on inversions for AI and true-amplitude reflectivity using 2D synthetic and real-data examples.


2016 ◽  
Vol 5 (2) ◽  
pp. 199
Author(s):  
Yuanyin Zhang ◽  
Zandong Sun ◽  
Zhijun Jin ◽  
Ning Dong ◽  
Yequan Chen ◽  
...  

For the modeling of complex reservoirs with strong heterogeneity, for instance the deeply buried paleokarst reservoirs in the Tarim Basin, the traditional method by lateral interpolation and extrapolation of measured logs between well locations with the guiding of interpreted seismic horizons is driven by distance and often leads to non-geologic solutions, while the past improvements via adding seismic velocities or attributes information are still not accurate due to the resolution limitation or AVO (amplitude versus offset) effects contamination. In this paper, we present an amplitude-based modeling method by utilizing the heterogeneous information from seismic data to guide the geological model construction, based on the inverted pure P-wave data which have removed the AVO effects. The proposed method is applied in the impedance inversion of the paleokarst carbonate reservoirs in the Tarim Basin, where the reservoirs are characterized by substantial heterogeneity. Both the constructed Low frequency model (LFM) and the inverted impedance results of proposed method are more correlative with drilling data than that of traditional method. This method is more beneficial for strong heterogeneous reservoirs description especially in well insufficient or absent areas, suggested by the comparisons with traditional methods in the ZG8 area.


Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. S187-S200 ◽  
Author(s):  
Dmitrii Merzlikin ◽  
Sergey Fomel ◽  
Mrinal K. Sen

Diffraction imaging aims to emphasize small-scale subsurface heterogeneities, such as faults, pinch-outs, fracture swarms, channels, etc. and can help seismic reservoir characterization. The key step in diffraction imaging workflows is based on the separation procedure suppressing higher energy reflections and emphasizing diffractions, after which diffractions can be imaged independently. Separation results often contain crosstalk between reflections and diffractions and are prone to noise. We have developed an inversion scheme to reduce the crosstalk and denoise diffractions. The scheme decomposes an input full wavefield into three components: reflections, diffractions, and noise. We construct the inverted forward modeling operator as the chain of three operators: Kirchhoff modeling, plane-wave destruction, and path-summation integral filter. Reflections and diffractions have the same modeling operator. Separation of the components is done by shaping regularization. We impose sparsity constraints to extract diffractions, enforce smoothing along dominant local event slopes to restore reflections, and suppress the crosstalk between the components by local signal-and-noise orthogonalization. Synthetic- and field-data examples confirm the effectiveness of the proposed method.


2015 ◽  
Vol 3 (1) ◽  
pp. SF69-SF79 ◽  
Author(s):  
Ioan Sturzu ◽  
Alexander Mihai Popovici ◽  
Tijmen Jan Moser ◽  
Sudha Sudhakar

Diffraction imaging is recognized as a new approach to image small-scale fractures in shale and carbonate reservoirs. By identifying the areas with increased natural fracture density, reservoir engineers can design an optimal well placement program that targets the sweet spots (areas with increased production), and minimizes the total number of wells used for a prospective area. High-resolution imaging of the small-scale fractures in shale reservoirs such as Eagle Ford, Bakken, Utica, and Woodbine in the US, and Horn River, Montney, and Utica in Canada improves the prospect characterization and predrill assessment of the geologic conditions, improves the production and recovery efficiency, reduces field development cost, and decreases the environmental impact of developing the field by using fewer wells to optimally produce the reservoir. We evaluated several field data examples using a method of obtaining images of diffractors using specularity filtering that could be performed in depth and time migration. Provided that a good migration velocity was available, we used the deviation of ray scattering from Snell’s law to attenuate reflection energy in the migrated image. The resulting diffraction images reveal much of the structural detail that was previously obscured by reflection energy.


Geophysics ◽  
2016 ◽  
Vol 81 (6) ◽  
pp. S399-S408 ◽  
Author(s):  
Jingtao Zhao ◽  
Suping Peng ◽  
Wenfeng Du ◽  
Xiaoting Li

Clarifying and locating small-scale discontinuities or inhomogeneities in the subsurface, such as faults and collapsed columns, plays a vital role in safe coal mining because these discontinuities or inhomogeneities may destroy the continuity of layers and result in dangerous mining accidents. Diffractions carry key information from these objects and therefore can be used for high-resolution imaging. However, diffracted/scattered waves are much weaker than reflected waves and consequently require separation before being imaged. We have developed a Mahalanobis-based diffraction imaging method by modifying the classic Kirchhoff formula with an exponential function to account for the dynamic differences between reflections and diffractions in the shot domain. The imaging method can automatically account for destroying of reflected waves, constructive stacking of diffracted waves, and strengthening of scattered waves. The method can overcome the difficulties in handling Fresnel apertures, and it is suitable for high-resolution imaging because of the consistency of the waveforms in the shot domain. Although the proposed method in principle requires a good migration velocity model for calculating elementary diffraction traveltimes, it is robust to an inaccurate migration velocity model. Two numerical experiments demonstrate the feasibility of the proposed method in removing reflections and highlighting diffractions, and one field application further confirms its efficiency in resolving masked faults and collapsed columns.


2019 ◽  
Vol 220 (3) ◽  
pp. 1569-1584
Author(s):  
Zhengwei Li ◽  
Jianfeng Zhang

SUMMARY Accurate identification of the locations and orientations of small-scale faults plays an important role in seismic interpretation. We have developed a 3-D migration scheme that can image small-scale faults using diffractions in time. This provides a resolution beyond the classical Rayleigh limit of half a wavelength in detecting faults. The scheme images weak diffractions by building a modified dip-angle gather, which is obtained by replacing the two dip angles dimensions of the conventional 2-D dip-angle gather with tangents of the dip angles. We build the modified 2-D dip-angle gathers by calculating the tangents of dip angles following 3-D prestack time migration (PSTM). In the resulting modified 2-D dip-angle gathers, the Fresnel zone related to the specular reflection exhibits an ellipse. Comparing with the conventional 2-D dip-angle gather, diffraction event related a fault exhibits a straight cylinder shape with phase-reversal across a line related the orientation of the fault. As a result, we can not only mute the Fresnel zones related to reflections, correct phase for edge diffractions and obtain the image of faults, but also detect the orientations of 3-D faults using the modified dip-angle gathers. Like the conventional dip-angle gathers, the modified dip-angle gathers can also be used to image diffractions resulting from other sources. 3-D Field data tests demonstrate the validity of the proposed diffraction imaging scheme.


Geophysics ◽  
2018 ◽  
Vol 83 (1) ◽  
pp. S47-S55 ◽  
Author(s):  
Parsa Bakhtiari Rad ◽  
Benjamin Schwarz ◽  
Dirk Gajewski ◽  
Claudia Vanelle

Diffraction imaging can lead to high-resolution characterization of small-scale subsurface structures. A key step of diffraction imaging and tomography is diffraction separation and enhancement, especially in the full prestack data volume. We have considered point diffractors and developed a robust and fully data-driven workflow for prestack diffraction separation based on wavefront attributes, which are determined using the common-reflection-surface (CRS) method. In the first of two steps, we apply a zero-offset-based extrapolation operator for prestack diffraction separation, which combines the robustness and stability of the zero-offset CRS processing with enhanced resolution and improved illumination of the finite-offset CRS processing. In the second step, when the finite-offset diffracted events are separated, we apply a diffraction-based time migration velocity model building that provides high-quality diffraction velocity spectra. Applications of the new workflow to 2D/3D complex synthetic data confirm the superiority of prestack diffraction separation over the poststack method as well as the high potential of diffractions for improved time imaging.


Sign in / Sign up

Export Citation Format

Share Document