A True-Amplitude Imaging Method Based on Gaussian Beam Migration and Demigration

2020 ◽  
Vol 177 (10) ◽  
pp. 4707-4718
Author(s):  
Shaoyong Liu ◽  
Rushan Wu ◽  
Bo Feng ◽  
Huazhong Wang ◽  
Song Guo
2019 ◽  
Vol 16 (6) ◽  
pp. 1301-1319 ◽  
Author(s):  
Rui Zhang ◽  
Jian-Ping Huang ◽  
Su-Bin Zhuang ◽  
Zhen-Chun Li

Abstract For large-scale 3D seismic data, target-oriented reservoir imaging is more attractive than conventional full-volume migration, in terms of computation efficiency. Gaussian beam migration (GBM) is one of the most robust depth imaging method, which not only keeps the advantages of ray methods, such as high efficiency and flexibility, but also allows us to solve caustics and multipathing problems. But conventional Gaussian beam migration requires slant stack for prestack data, and ray tracing from beam center location to subsurface, which is not easy to be directly applied for target-oriented imaging. In this paper, we modify the conventional Gaussian beam migration scheme, by shooting rays from subsurface image points to receivers to implement wavefield back-propagation. This modification helps us to achieve a better subsurface illumination in complex structure and allows simple implementation for target reservoir imaging. Significantly, compared with the wavefield-based GBM, our method does not reconstruct the subsurface snapshots, which has higher efficiency. But the proposed method is not as efficient as the conventional Gaussian beam migration. Synthetic and field data examples demonstrate the validity and the target-oriented imaging capability of our method.


2020 ◽  
Author(s):  
Shaoyong Liu ◽  
Wenting Zhu ◽  
Ru-Shan Wu ◽  
Huazhong Wang ◽  
Song Guo

Geophysics ◽  
2015 ◽  
Vol 80 (1) ◽  
pp. S43-S54 ◽  
Author(s):  
Nizare El Yadari

Acoustic Gaussian beam migration is an attractive imaging method because it is flexible with input geometry, efficient, and accurate in imaging multipath arrivals. However, one of the hurdles that this method must overcome in production processing is its extension to use multimeasurement data, as recently allowed by novel acquisition technologies. This is inevitable when the compensation of the ghost effect is best corrected within a true-amplitude imaging process, a necessity for amplitude-variation-with-offset work. For this purpose, I introduced a novel formalism for vector-acoustic imaging, based on Green’s function theory, which can remove the ghost effect and produce amplitudes on reflectors that are proportional to the reflection coefficients. I established a theoretical framework with Gaussian beam representations of Green’s functions, including the weighted beam-stacking approach that reduced the cost of computation. I extended my formulas to use the steep-descent (i.e., stationary phase) approximation. Then, I explained the impact of this approximation on the illumination and the event continuity and sharpness. I also studied the special case of acoustic imaging corresponding to using single-measurement (i.e., pressure) data. I applied the derived formulations to realistic synthetic multisensor data (North Sea) using a research code of Gaussian beam migration. The numerical examples demonstrated that I can improve the illumination of the final images and obtain wide-bandwidth reflectivity maps.


Geophysics ◽  
2019 ◽  
Vol 84 (4) ◽  
pp. S329-S340 ◽  
Author(s):  
Yubo Yue ◽  
Paul Sava ◽  
Zhongping Qian ◽  
Jidong Yang ◽  
Zhen Zou

Gaussian beam migration (GBM) is an effective imaging method that has the ability to image multiple arrivals while preserving the advantages of ray-based methods. We have extended this method to linearized least-squares imaging for elastic waves in isotropic media. We have dynamically transformed the multicomponent data to the principal components of different wave modes using the polarization information available in the beam migration process, and then we use Gaussian beams as wavefield propagator to construct the forward modeling and adjoint migration operators. Based on the constructed operators, we formulate a least-squares migration scheme that is iteratively solved using a preconditioned conjugate gradient method. With this method, we can obtain crosstalk-attenuated multiwave images with better subsurface illumination and higher resolution than those of the conventional elastic Gaussian beam migration. This method also allows us to achieve a good balance between computational cost and imaging accuracy, which are both important requirements for iterative least-squares migrations. Numerical tests on two synthetic data sets demonstrate the validity and effectiveness of our proposed method.


Energies ◽  
2021 ◽  
Vol 14 (14) ◽  
pp. 4105
Author(s):  
Shaoyong Liu ◽  
Wenting Zhu ◽  
Zhe Yan ◽  
Peng Xu ◽  
Huazhong Wang

The estimation of the subsurface acoustic impedance (AI) model is an important step of seismic data processing for oil and gas exploration. The full waveform inversion (FWI) is a powerful way to invert the subsurface parameters with surface acquired seismic data. Nevertheless, the strong nonlinear relationship between the seismic data and the subsurface model will cause nonconvergence and unstable problems in practice. To divide the nonlinear inversion into some more linear steps, a 2D AI inversion imaging method is proposed to estimate the broadband AI model based on a broadband reflectivity. Firstly, a novel scheme based on Gaussian beam migration (GBM) is proposed to produce the point spread function (PSF) and conventional image of the subsurface. Then, the broadband reflectivity can be obtained by implementing deconvolution on the image with respect to the calculated PSF. Assuming that the low-wavenumber part of the AI model can be deduced by the background velocity, we implemented the AI inversion imaging scheme by merging the obtained broadband reflectivity as the high-wavenumber part of the AI model and produced a broadband AI result. The developed broadband migration based on GBM as the computational hotspot of the proposed 2D AI inversion imaging includes only two GBM and one Gaussian beam demigraton (Born modeling) processes. Hence, the developed broadband GBM is more efficient than the broadband imaging using the least-squares migrations (LSMs) that require multiple iterations (every iteration includes one Born modeling and one migration process) to minimize the objective function of data residuals. Numerical examples of both synthetic data and field data have demonstrated the validity and application potential of the proposed method.


Geophysics ◽  
1990 ◽  
Vol 55 (11) ◽  
pp. 1416-1428 ◽  
Author(s):  
N. Ross Hill

Just as synthetic seismic data can be created by expressing the wave field radiating from a seismic source as a set of Gaussian beams, recorded data can be downward continued by expressing the recorded wave field as a set of Gaussian beams emerging at the earth’s surface. In both cases, the Gaussian beam description of the seismic‐wave propagation can be advantageous when there are lateral variations in the seismic velocities. Gaussian‐beam downward continuation enables wave‐equation calculation of seismic propagation, while it retains the interpretive raypath description of this propagation. This paper describes a zero‐offset depth migration method that employs Gaussian beam downward continuation of the recorded wave field. The Gaussian‐beam migration method has advantages for imaging complex structures. Like finite‐difference migration, it is especially compatible with lateral variations in velocity, but Gaussian beam migration can image steeply dipping reflectors and will not produce unwanted reflections from structure in the velocity model. Unlike other raypath methods, Gaussian beam migration has guaranteed regular behavior at caustics and shadows. In addition, the method determines the beam spacing that ensures efficient, accurate calculations. The images produced by Gaussian beam migration are usually stable with respect to changes in beam parameters.


Geophysics ◽  
2018 ◽  
Vol 83 (5) ◽  
pp. T281-T289 ◽  
Author(s):  
Qianru Xu ◽  
Weijian Mao

We have developed a fast ray-tracing method for multiple layered inhomogeneous anisotropic media, based on the generalized Snell’s law. Realistic geologic structures continuously varying with embedded discontinuities are parameterized by adopting cubic B-splines with nonuniformly spaced nodes. Because the anisotropic characteristic is often closely related to the interface configuration, this model parameterization scheme containing the natural inclination of the corresponding layer is particularly suitable for tilted transverse isotropic models whose symmetry axis is generally perpendicular to the direction of the layers. With this model parameterization, the first- and second-order spatial derivatives of the velocity within the interfaces can be effectively obtained, which facilitates the amplitude computation in dynamic ray tracing. By using complex initial conditions for the dynamic ray system and taking the multipath effect into consideration, our method is applicable to Gaussian beam migration. Numerical experiments of our method have been used to verify its effectiveness, practicability, and efficiency in memory storage and computation.


Sign in / Sign up

Export Citation Format

Share Document