Wavefield decomposition in arbitrary direction and an imaging condition based on stratigraphic dip

Xuebao Guo,Ying Shi,Weihong Wang,Hongliang Jing,Zhen Zhang

DOI: https://doi.org/10.1190/geo2019-0617.1

IF: 3.264

2020-09-01

Geophysics

Abstract:In reverse time migration (RTM), wavefield decomposition can play an important role in addressing the issue of migration noise, especially low-frequency noise. The complete wavefield decomposition based on the Hilbert transform is a commonly used method in RTM, but it is accompanied by extra wavefield simulation and wavefield storage. We have developed three distinct methods. The first is a convenient method for wavefield decomposition, which is based on Poynting vectors. Only the unit vector in one direction is needed to realize the wavefield decomposition in an arbitrary direction by this method. It breaks through the limitation that the Hilbert transform-based method is applicable only to the up- and downgoing wave or left- and right-going wave decomposition, and the calculation cost is negligible compared with RTM. The second is a method based on the instantaneous wavenumber, which we developed for calculating the wave propagation direction. On the basis of wavefield decomposition, the imaging angle gather from the new method performs better than that of the Poynting vector method. Meanwhile, it also is used for generating the incident angle gather and dip angle gather. The latter expresses the dip angle of underground strata. More importantly, the above methods allow us to control the wavefield decomposition direction and three angles at any position underground. The third adopts a stratigraphic imaging condition method, and we briefly analyze the relationship between the new method and the inverse-scattering imaging condition. The stratigraphic imaging condition maps the results to the dip angle of the stratum through a spatial gradient wavefield, which can enhance the effective imaging information. The above three kinds of angle gathers also can be constructed by the stratigraphic imaging condition. Numerical experiments demonstrate that the imaging results and the angle gathers obtained by our proposed method have higher accuracy and resolution.

geochemistry & geophysics

What problem does this paper attempt to address?

This paper attempts to solve the following key problems:

1. **Suppression of low - frequency noise**: In reverse - time migration (RTM), low - frequency noise is a major problem. Traditional wave - field decomposition methods (such as methods based on the Hilbert transform) can effectively separate up - going and down - going waves or left - going and right - going waves, but they have a high computational cost and can only be applied to wave - field decomposition in specific directions.

Formula examples:
\[
𝒅𝑝𝑜𝑦𝑛 = -𝑃∇𝑃
\]
\[
cos𝜃𝑝𝑜𝑦𝑛 = \frac{𝒅𝑠𝑝𝑜𝑦𝑛 ∙ 𝒅𝑟𝑝𝑜𝑦𝑛}{|𝑑𝑠𝑝𝑜𝑦𝑛||𝑑𝑟𝑝𝑜𝑦𝑛|}
\]

2. **Wave - field decomposition in arbitrary directions**: The paper proposes a new method based on the Poynting vector, which can achieve wave - field decomposition in arbitrary directions. This method can achieve wave - field decomposition with only one unit vector, breaks through the limitations of traditional methods, and has a low computational cost.

3. **Imaging angle estimation**: The paper introduces the instantaneous wavenumber method to estimate the propagation angle of seismic waves, and then obtains the imaging angle, incident angle and formation dip. This method can generate more focused angle gathers, although it cannot distinguish between seismic waves propagating in parallel.

Formula examples:
\[
𝑃𝑧𝐻(𝑥,𝑧) = \frac{1}{𝜋} \int_{𝑅} \frac{𝑃(𝑥,𝑧―𝜏)}{𝜏} 𝑑𝜏
\]
\[
𝛼𝑖𝑛𝑠𝑡 = \tan^{-1}\left(\frac{𝑤1}{𝑤2}\right)
\]

4. **Improved imaging conditions**: The paper proposes imaging conditions based on the formation dip, which maps the imaging information to the formation dip, thereby enhancing the effective imaging information and suppressing non - dip noise. This imaging condition can be regarded as a special form of the inverse - scattering imaging condition.

Formula examples:
\[
𝐼stru(𝑥,𝑧) = \int [𝑛∙∇𝑃𝑠(𝑥,𝑧,𝑡)][𝑛∙∇𝑃𝑟(𝑥,𝑧,𝑡)] 𝑑𝑡
\]

5. **Numerical experiment verification**: Through numerical experiments using the Sigsbee 2A model and the Marmousi model, the effectiveness and accuracy of the proposed method are verified. The results show that the new method has higher performance in terms of imaging quality and resolution.

In summary, the main objective of this paper is to develop an efficient and flexible method for wave - field decomposition and imaging conditions to solve the low - frequency noise problem in RTM and improve imaging accuracy.