An accurate and wide-angle one-way propagator for wavefield extrapolation is an important topic for research on wave-equation prestack depth migration in the presence of large and rapid velocity variations. Based on t...An accurate and wide-angle one-way propagator for wavefield extrapolation is an important topic for research on wave-equation prestack depth migration in the presence of large and rapid velocity variations. Based on the optimal separable approximation presented in this paper, the mixed domain algorithm with forward and inverse Fourier transforms is used to construct the 3D one-way wavefield extrapolation operator. This operator separates variables in the wavenumber and spatial domains. The phase shift operation is implemented in the wavenumber domain while the time delay for lateral velocity variation is corrected in the spatial domain. The impulse responses of the one-way wave operator show that the numeric computation is consistent with the theoretical value for each velocity, revealing that the operator constructed with the optimal separable approximation can be applied to lateral velocity variations for the case of small steps. Imaging results of the SEG/EAGE model and field data indicate that the new method can be used to image complex structure.展开更多
Multiple prediction and subtraction techniques based on wavefield extrapolation are effective for suppressing multiple related to water layers. In the conventional wavefield extrapolation method,the multiples of the s...Multiple prediction and subtraction techniques based on wavefield extrapolation are effective for suppressing multiple related to water layers. In the conventional wavefield extrapolation method,the multiples of the seismic data are predicted from the known total wave field by the Green function convoluted with each point of the bottom. However,only the energy near the stationary phase point has an effect on the summation result when the convolutional gathers are added. The research proposed a stationary phase point extraction method based on high-resolution radon transform. In the radon domain,the energy near the stationary phase point is directly added along the convolutional gathers curve,which is a valid solution to the problem of the unstable phase of the events of multiple. The Curvelet matching subtraction technique is used to remove the multiple,which improved the accuracy of the multiple predicted by the wavefield extrapolation and the artifacts appearing around the events of multiple are well eliminated. The validity and feasibility of the proposed method are verified by the theoretical and practical data example.展开更多
Seismic anisotropy has been extensively acknowledged as a crucial element that influences the wave propagation characteristic during wavefield simulation,inversion and imaging.Transversely isotropy(TI)and orthorhombic...Seismic anisotropy has been extensively acknowledged as a crucial element that influences the wave propagation characteristic during wavefield simulation,inversion and imaging.Transversely isotropy(TI)and orthorhombic anisotropy(OA)are two typical categories of anisotropic media in exploration geophysics.In comparison of the elastic wave equations in both TI and OA media,pseudo-acoustic wave equations(PWEs)based on the acoustic assumption can markedly reduce computational cost and complexity.However,the presently available PWEs may experience SV-wave contamination and instability when anisotropic parameters cannot satisfy the approximated condition.Exploiting pure-mode wave equations can effectively resolve the above-mentioned issues and generate pure P-wave events without any artifacts.To further improve the computational accuracy and efficiency,we develop two novel pure qP-wave equations(PPEs)and illustrate the corresponding numerical solutions in the timespace domain for 3D tilted TI(TTI)and tilted OA(TOA)media.First,the rational polynomials are adopted to estimate the exact pure qP-wave dispersion relations,which contain complicated pseudo-differential operators with irrational forms.The polynomial coefficients are produced by applying a linear optimization algorithm to minimize the objective function difference between the expansion formula and the exact one.Then,the developed optimized PPEs are efficiently implemented using the finite-difference(FD)method in the time-space domain by introducing a scalar operator,which can help avoid the problem of spectral-based algorithms and other calculation burdens.Structures of the new equations are concise and corresponding implementation processes are straightforward.Phase velocity analyses indicate that our proposed optimized equations can lead to reliable approximation results.3D synthetic examples demonstrate that our proposed FD-based PPEs can produce accurate and stable P-wave responses,and effectively describe the wavefield features in complicated TTI and TOA media.展开更多
Imaging sea-bed sediment layers from echo data, which are collected by a system composed of a seismic profiler and a hydrophone streamer towed behind the profiler, is a way to reconstruct the structure of sedimeat lay...Imaging sea-bed sediment layers from echo data, which are collected by a system composed of a seismic profiler and a hydrophone streamer towed behind the profiler, is a way to reconstruct the structure of sedimeat layers with acoustic wav equation. The equation which describes the wave propagation is used for backward extrapolation of echo data observed at sea surface. When the medium is homogeneous or horizontally layered, time imaging approach is valid. However, in the case where a considerable lateral variation in velocity exists, the image section processed with the time approach does not represent the real structure, because of distortions caused by thin-lens effect similar as in optics. In this case, depth imaging approach must be used for both the time-shift correction of refraction terms and the convergence of diffractions simultaneously as wavefields are downward continued. As a result, the good image can be derived to determine the structure of sea-bed sediment layers.展开更多
基金This research is sponsored by China National Natural Science Foundation （N0. 40474047）.
文摘An accurate and wide-angle one-way propagator for wavefield extrapolation is an important topic for research on wave-equation prestack depth migration in the presence of large and rapid velocity variations. Based on the optimal separable approximation presented in this paper, the mixed domain algorithm with forward and inverse Fourier transforms is used to construct the 3D one-way wavefield extrapolation operator. This operator separates variables in the wavenumber and spatial domains. The phase shift operation is implemented in the wavenumber domain while the time delay for lateral velocity variation is corrected in the spatial domain. The impulse responses of the one-way wave operator show that the numeric computation is consistent with the theoretical value for each velocity, revealing that the operator constructed with the optimal separable approximation can be applied to lateral velocity variations for the case of small steps. Imaging results of the SEG/EAGE model and field data indicate that the new method can be used to image complex structure.
基金Supported by the National Science and Technology Major Project(No.2016ZX05026-002-003)the National Natural Science Foundation of China(No.41374108)
文摘Multiple prediction and subtraction techniques based on wavefield extrapolation are effective for suppressing multiple related to water layers. In the conventional wavefield extrapolation method,the multiples of the seismic data are predicted from the known total wave field by the Green function convoluted with each point of the bottom. However,only the energy near the stationary phase point has an effect on the summation result when the convolutional gathers are added. The research proposed a stationary phase point extraction method based on high-resolution radon transform. In the radon domain,the energy near the stationary phase point is directly added along the convolutional gathers curve,which is a valid solution to the problem of the unstable phase of the events of multiple. The Curvelet matching subtraction technique is used to remove the multiple,which improved the accuracy of the multiple predicted by the wavefield extrapolation and the artifacts appearing around the events of multiple are well eliminated. The validity and feasibility of the proposed method are verified by the theoretical and practical data example.
基金supported by the National Key R&D Program of China(2021YFA0716902)National Natural Science Foundation of China(NSFC)under contract number 42374149 and 42004119National Science and Technology Major Project(2024ZD1002907)。
文摘Seismic anisotropy has been extensively acknowledged as a crucial element that influences the wave propagation characteristic during wavefield simulation,inversion and imaging.Transversely isotropy(TI)and orthorhombic anisotropy(OA)are two typical categories of anisotropic media in exploration geophysics.In comparison of the elastic wave equations in both TI and OA media,pseudo-acoustic wave equations(PWEs)based on the acoustic assumption can markedly reduce computational cost and complexity.However,the presently available PWEs may experience SV-wave contamination and instability when anisotropic parameters cannot satisfy the approximated condition.Exploiting pure-mode wave equations can effectively resolve the above-mentioned issues and generate pure P-wave events without any artifacts.To further improve the computational accuracy and efficiency,we develop two novel pure qP-wave equations(PPEs)and illustrate the corresponding numerical solutions in the timespace domain for 3D tilted TI(TTI)and tilted OA(TOA)media.First,the rational polynomials are adopted to estimate the exact pure qP-wave dispersion relations,which contain complicated pseudo-differential operators with irrational forms.The polynomial coefficients are produced by applying a linear optimization algorithm to minimize the objective function difference between the expansion formula and the exact one.Then,the developed optimized PPEs are efficiently implemented using the finite-difference(FD)method in the time-space domain by introducing a scalar operator,which can help avoid the problem of spectral-based algorithms and other calculation burdens.Structures of the new equations are concise and corresponding implementation processes are straightforward.Phase velocity analyses indicate that our proposed optimized equations can lead to reliable approximation results.3D synthetic examples demonstrate that our proposed FD-based PPEs can produce accurate and stable P-wave responses,and effectively describe the wavefield features in complicated TTI and TOA media.
文摘Imaging sea-bed sediment layers from echo data, which are collected by a system composed of a seismic profiler and a hydrophone streamer towed behind the profiler, is a way to reconstruct the structure of sedimeat layers with acoustic wav equation. The equation which describes the wave propagation is used for backward extrapolation of echo data observed at sea surface. When the medium is homogeneous or horizontally layered, time imaging approach is valid. However, in the case where a considerable lateral variation in velocity exists, the image section processed with the time approach does not represent the real structure, because of distortions caused by thin-lens effect similar as in optics. In this case, depth imaging approach must be used for both the time-shift correction of refraction terms and the convergence of diffractions simultaneously as wavefields are downward continued. As a result, the good image can be derived to determine the structure of sea-bed sediment layers.
