3D self-imaging condition for finite aperture objects

1996 ◽  
Vol 132 (1-2) ◽  
pp. 41-47 ◽  
Author(s):  
R Moignard ◽  
J.L de Bougrenet de la Tocnaye
Geophysics ◽  
1984 ◽  
Vol 49 (11) ◽  
pp. 1896-1907 ◽  
Author(s):  
P. Temme

A comparison of common‐midpoint (CMP), single‐shot, and plane‐wave migration was made for simple two‐dimensional structures such as a syncline and a horizontal reflector with a laterally variable reflection coefficient by using synthetic seismograms. The seismograms were calculated employing the finite‐difference technique. CMP sections were simulated by 18-fold stacking and plane‐wave sections by slant stacking. By applying a finite‐difference scheme, the synthetic wave field was continued downward. The usual imaging condition of CMP migration was extended in order to carry out migration of single‐shot and plane‐wave sections. The reflection coefficient was reconstructed by comparing the migrated wave field with the incident wave field at the reflector. The results are: (1) all three migration techniques succeeded in reconstructing the reflector position; (2) as a consequence of the finite aperture of the geophone spread, only segments of the reflector could be reconstructed by single‐shot and plane‐wave migration; (3) for single‐shot and plane‐wave migration the reflection coefficient could be obtained; and (4) CMP migration may lead to incorrect conclusions regarding the reflection coefficient.


Geophysics ◽  
2008 ◽  
Vol 73 (2) ◽  
pp. S47-S61 ◽  
Author(s):  
Paul Sava ◽  
Oleg Poliannikov

The fidelity of depth seismic imaging depends on the accuracy of the velocity models used for wavefield reconstruction. Models can be decomposed in two components, corresponding to large-scale and small-scale variations. In practice, the large-scale velocity model component can be estimated with high accuracy using repeated migration/tomography cycles, but the small-scale component cannot. When the earth has significant small-scale velocity components, wavefield reconstruction does not completely describe the recorded data, and migrated images are perturbed by artifacts. There are two possible ways to address this problem: (1) improve wavefield reconstruction by estimating more accurate velocity models and image using conventional techniques (e.g., wavefield crosscorrelation) or (2) reconstruct wavefields with conventional methods using the known background velocity model but improve the imaging condition to alleviate the artifacts caused by the imprecise reconstruction. Wedescribe the unknown component of the velocity model as a random function with local spatial correlations. Imaging data perturbed by such random variations is characterized by statistical instability, i.e., various wavefield components image at wrong locations that depend on the actual realization of the random model. Statistical stability can be achieved by preprocessing the reconstructed wavefields prior to the imaging condition. We use Wigner distribution functions to attenuate the random noise present in the reconstructed wavefields, parameterized as a function of image coordinates. Wavefield filtering using Wigner distribution functions and conventional imaging can be lumped together into a new form of imaging condition that we call an interferometric imaging condition because of its similarity to concepts from recent work on interferometry. The interferometric imaging condition can be formulated both for zero-offset and for multioffset data, leading to robust, efficient imaging procedures that effectively attenuate imaging artifacts caused by unknown velocity models.


Geophysics ◽  
2017 ◽  
Vol 82 (2) ◽  
pp. S111-S127 ◽  
Author(s):  
Qizhen Du ◽  
ChengFeng Guo ◽  
Qiang Zhao ◽  
Xufei Gong ◽  
Chengxiang Wang ◽  
...  

The scalar images (PP, PS, SP, and SS) of elastic reverse time migration (ERTM) can be generated by applying an imaging condition as crosscorrelation of pure wave modes. In conventional ERTM, Helmholtz decomposition is commonly applied in wavefield separation, which leads to a polarity reversal problem in converted-wave images because of the opposite polarity distributions of the S-wavefields. Polarity reversal of the converted-wave image will cause destructive interference when stacking over multiple shots. Besides, in the 3D case, the curl calculation generates a vector S-wave, which makes it impossible to produce scalar PS, SP, and SS images with the crosscorrelation imaging condition. We evaluate a vector-based ERTM (VB-ERTM) method to address these problems. In VB-ERTM, an amplitude-preserved wavefield separation method based on decoupled elastic wave equation is exploited to obtain the pure wave modes. The output separated wavefields are both vectorial. To obtain the scalar images, the scalar imaging condition in which the scalar product of two vector wavefields with source-normalized illumination is exploited to produce scalar images instead of correlating Cartesian components or magnitude of the vector P- and S-wave modes. Compared with alternative methods for correcting the polarity reversal of PS and SP images, our ERTM solution is more stable and simple. Besides these four scalar images, the VB-ERTM method generates another PP-mode image by using the auxiliary stress wavefields. Several 2D and 3D numerical examples are evaluated to demonstrate the potential of our ERTM method.


2021 ◽  
Author(s):  
Matthew Bailey ◽  
Adrian Hilton ◽  
Jean-Yves Guillemaut

Geophysics ◽  
2018 ◽  
Vol 83 (6) ◽  
pp. S569-S577 ◽  
Author(s):  
Yang Zhao ◽  
Houzhu Zhang ◽  
Jidong Yang ◽  
Tong Fei

Using the two-way elastic-wave equation, elastic reverse time migration (ERTM) is superior to acoustic RTM because ERTM can handle mode conversions and S-wave propagations in complex realistic subsurface. However, ERTM results may not only contain classical backscattering noises, but they may also suffer from false images associated with primary P- and S-wave reflections along their nonphysical paths. These false images are produced by specific wave paths in migration velocity models in the presence of sharp interfaces or strong velocity contrasts. We have addressed these issues explicitly by introducing a primary noise removal strategy into ERTM, in which the up- and downgoing waves are efficiently separated from the pure-mode vector P- and S-wavefields during source- and receiver-side wavefield extrapolation. Specifically, we investigate a new method of vector wavefield decomposition, which allows us to produce the same phases and amplitudes for the separated P- and S-wavefields as those of the input elastic wavefields. A complex function involved with the Hilbert transform is used in up- and downgoing wavefield decomposition. Our approach is cost effective and avoids the large storage of wavefield snapshots that is required by the conventional wavefield separation technique. A modified dot-product imaging condition is proposed to produce multicomponent PP-, PS-, SP-, and SS-images. We apply our imaging condition to two synthetic models, and we demonstrate the improvement on the image quality of ERTM.


Geophysics ◽  
2019 ◽  
Vol 84 (2) ◽  
pp. S95-S111 ◽  
Author(s):  
Wei Zhang ◽  
Ying Shi

Elastic reverse time migration (RTM) has the ability to retrieve accurately migrated images of complex subsurface structures by imaging the multicomponent seismic data. However, the imaging condition applied in elastic RTM significantly influences the quality of the migrated images. We evaluated three kinds of imaging conditions in elastic RTM. The first kind of imaging condition involves the crosscorrelation between the Cartesian components of the particle-velocity wavefields to yield migrated images of subsurface structures. An alternative crosscorrelation imaging condition between the separated pure wave modes obtained by a Helmholtz-like decomposition method could produce reflectivity images with explicit physical meaning and fewer crosstalk artifacts. A drawback of this approach, though, was that the polarity reversal of the separated S-wave could cause destructive interference in the converted-wave image after stacking over multiple shots. Unlike the conventional decomposition method, the elastic wavefields can also be decomposed in the vector domain using the decoupled elastic wave equation, which preserves the amplitude and phase information of the original elastic wavefields. We have developed an inner-product imaging condition to match the vector-separated P- and S-wave modes to obtain scalar reflectivity images of the subsurface. Moreover, an auxiliary P-wave stress image can supplement the elastic imaging. Using synthetic examples with a layered model, the Marmousi 2 model, and a fault model, we determined that the inner-product imaging condition has prominent advantages over the other two imaging conditions and generates images with preserved amplitude and phase attributes.


Geophysics ◽  
2019 ◽  
Vol 84 (2) ◽  
pp. KS13-KS27 ◽  
Author(s):  
Daniel Rocha ◽  
Paul Sava ◽  
Jeffrey Shragge ◽  
Ben Witten

In passive seismic monitoring of microseismicity, full-wavefield imaging offers a robust approach for the estimation of source location and mechanism. With multicomponent data and the full 3D anisotropic elastic wave equation, the coexistence of P- and S-modes at the source location in time-reversal wavefield extrapolation allows the development of imaging conditions that identify the source position and radiation pattern. We have developed an imaging condition for passive wavefield imaging that is based on energy conservation and is related to the source mechanism. Similar to the correlation between the decomposed P- and S-wavefields — the most common imaging condition used in passive elastic wavefield imaging — our proposed imaging condition compares the different modes present in the displacement field producing a strong and focused correlation at the source location without costly wave-mode decomposition at each time step. Numerical experiments demonstrate the advantages of the proposed imaging condition (compared to PS correlation with decomposed wave modes), its sensitivity with respect to velocity inaccuracy, and its quality and efficacy in estimating the source location.


2018 ◽  
Vol 125 (6) ◽  
pp. 1720-1730 ◽  
Author(s):  
Eric T. Geier ◽  
Kent Kubo ◽  
Rebecca J. Theilmann ◽  
Gordon K. Prisk ◽  
Rui C. Sá

The location of lung regions with compromised ventilation (often called ventilation defects) during a bronchoconstriction event may be influenced by posture. We aimed to determine the effect of prone versus supine posture on the spatial pattern of methacholine-induced bronchoconstriction in six healthy adults (ages 21–41, 3 women) using specific ventilation imaging. Three postural conditions were chosen to assign the effect of posture to the drug administration and/or imaging phase of the experiment: supine methacholine administration followed by supine imaging, prone methacholine administration followed by supine imaging, and prone methacholine administration followed by prone imaging. The two conditions in which imaging was performed supine had similar spatial patterns of bronchoconstriction despite a change in posture during methacholine administration; the odds ratio for recurrent constriction was mean (SD) = 7.4 (3.9). Conversely, dissimilar spatial patterns of bronchoconstriction emerged when posture during imaging was changed; the odds ratio for recurrent constriction between the prone methacholine/supine imaging condition and the prone methacholine/prone imaging condition was 1.2 (0.9). Logistic regression showed that height above the dependent lung border was a significant negative predictor of constriction in the two supine imaging conditions ( P < 0.001 for each) but not in the prone imaging condition ( P = 0.20). These results show that the spatial pattern of methacholine bronchoconstriction is recurrent in the supine posture, regardless of whether methacholine is given prone or supine but that prone posture during imaging eliminates that recurrent pattern and reduces its dependence on gravitational height. NEW & NOTEWORTHY The spatial pattern of methacholine bronchoconstriction in the supine posture is recurrent and skewed toward the dependent lung, regardless of whether inhaled methacholine is administered while supine or while prone. However, both the recurrent pattern and the gravitational skew are eliminated if imaging is performed prone. These results suggest that gravitational influence on regional lung inflation and airway topography at the time of measurement play a role in determining regional bronchoconstriction in the healthy lung.


Sign in / Sign up

Export Citation Format

Share Document