Compared with the transverse isotropic(TI)medium,the orthorhombic anisotropic medium has both horizontal and vertical symmetry axes and it can be approximated as a set of vertical fissures developed in a group of hori...Compared with the transverse isotropic(TI)medium,the orthorhombic anisotropic medium has both horizontal and vertical symmetry axes and it can be approximated as a set of vertical fissures developed in a group of horizontal strata.Although the full-elastic orthorhombic anisotropic wave equation can accurately simulate seismic wave propagation in the underground media,a huge computational cost is required in seismic modeling,migration,and inversion.The conventional coupled pseudo-acoustic wave equations based on acoustic approximation can be used to significantly reduce the cost of calculation.However,these equations usually suffer from unwanted shear wave artifacts during wave propagation,and the presence of these artifacts can significantly degrade the imaging quality.To solve these problems,we derived a new pure P-wave equation for orthorhombic media that eliminates shear wave artifacts while compromising computational efficiency and accuracy.In addition,the derived equation involves pseudo-differential operators and it must be solved by 3D FFT algorithms.In order to reduce the number of 3D FFT,we utilized the finite difference and pseudo-spectral methods to conduct 3D forward modeling.Furthermore,we simplified the equation by using elliptic approximation and implemented 3D reverse-time migration(RTM).Forward modeling tests on several homogeneous and heterogeneous models confirm that the accuracy of the new equation is better than that of conventional methods.3D RTM imaging tests on three-layer and SEG/EAGE 3D salt models confirm that the ORT media have better imaging quality.展开更多
Least-squares reverse-time migration(LSRTM) formulates reverse-time migration(RTM) in the leastsquares inversion framework to obtain the optimal reflectivity image. It can generate images with more accurate amplitudes...Least-squares reverse-time migration(LSRTM) formulates reverse-time migration(RTM) in the leastsquares inversion framework to obtain the optimal reflectivity image. It can generate images with more accurate amplitudes, higher resolution, and fewer artifacts than RTM. However, three problems still exist:(1) inversion can be dominated by strong events in the residual;(2) low-wavenumber artifacts in the gradient affect convergence speed and imaging results;(3) high-wavenumber noise is also amplified as iteration increases. To solve these three problems, we have improved LSRTM: firstly, we use Hubernorm as the objective function to emphasize the weak reflectors during the inversion;secondly, we adapt the de-primary imaging condition to remove the low-wavenumber artifacts above strong reflectors as well as the false high-wavenumber reflectors in the gradient;thirdly, we apply the L1-norm sparse constraint in the curvelet-domain as the regularization term to suppress the high-wavenumber migration noise. As the new inversion objective function contains the non-smooth L1-norm, we use a modified iterative soft thresholding(IST) method to update along the Polak-Ribie re conjugate-gradient direction by using a preconditioned non-linear conjugate-gradient(PNCG) method. The numerical examples,especially the Sigsbee2 A model, demonstrate that the Huber inversion-based RTM can generate highquality images by mitigating migration artifacts and improving the contribution of weak reflection events.展开更多
Research on seismic anisotropy and attenuation plays a significant role in exploration geophysics. To enhance the imaging quality for complicated structures, we develop several effective improvements for anisotropic a...Research on seismic anisotropy and attenuation plays a significant role in exploration geophysics. To enhance the imaging quality for complicated structures, we develop several effective improvements for anisotropic attenuation effects in reverse-time migration (Q-RTM) on surface and vertical seismic profiling (VSP) acquisition geometries. First, to suppress pseudo-shear wave artifact and numerical instability of the commonly used anisotropic pseudo-acoustic wave equations, an optimized pure P-wave dispersion relation is derived and the corresponding pure-mode wave equation is solved by combining the finite-difference and Possion methods. Second, a simplified anisotropic pure-mode visco-acoustic wave equation (PVAWE) based on standard linear solid model is established. Third, a time-dispersion correlation strategy is applied to improve the modeling accuracy. Fourth, we extend a target-oriented scheme to anisotropic attenuated modeling and imaging. Instead of the conventional wavefield modeling and RTM, the proposed approach can extract available wavefield information near the target regions and produce high imaging resolution for target structures. Last, both anisotropic surface and VSP Q-RTMs are executed by combining optimized PVAWE, time-dispersion correlation and target-oriented algorithm. Modeling examples demonstrate the advantages of our schemes. Moreover, our modified Q-compensated imaging workflow can be regarded as a supplement to the classical anisotropic RTM.展开更多
Conventional seismic exploration method based on post-stack data usually fails to identify the distribution of fractured and caved carbonate reservoirs in the Tarim Basin,so the rich pre-stack information should be ap...Conventional seismic exploration method based on post-stack data usually fails to identify the distribution of fractured and caved carbonate reservoirs in the Tarim Basin,so the rich pre-stack information should be applied to the prediction of carbonate reservoirs.Amplitude-preserved seismic data processing is the foundation.In this paper,according to the feature of desert seismic data （including weak reflection,fast attenuation of high frequency components,strong coherent noises,low S/N and resolution）,a set of amplitude-preserved processing techniques is applied and a reasonable processing flow is formed to obtain the high quality data.After implementing a set of pre-stack amplitude-preserved processing,we test and define the kernel parameters of amplitude-preserved Kirchhoff PSTM （pre-stack time migration） and subsequent gathers processing,in order to obtain the amplitude-preserved gathers used to the isotropic pre-stack inversion for the identification of caved reservoirs.The AVO characteristics of obtained gathers fit well with the synthetic gathers from logging data,and it proves that the processing above is amplitudepreserved.The azimuthal processing techniques,including azimuth division and binning enlargement,are implemented for amplitude-preserved azimuthal gathers with the uniform fold.They can be used in the anisotropic inversion to detect effective fractures.The processing techniques and flows are applied to the field seismic data,and are proved available for providing the amplitude-preserved gathers for carbonate reservoir prediction in the Tarim Basin.展开更多
High-precision seismic imaging is the core task of seismic exploration,guaranteeing the accuracy of geophysical and geological interpretation.With the development of seismic exploration,the targets become more and mor...High-precision seismic imaging is the core task of seismic exploration,guaranteeing the accuracy of geophysical and geological interpretation.With the development of seismic exploration,the targets become more and more complex.Imaging on complex media such as subsalt,small-scale,steeply dipping and surface topography structures brings a great challenge to imaging techniques.Therefore,the seismic imaging methods range from stacking-to migration-to inversion-based imaging,and the imaging accuracy is becoming increasingly high.This review paper includes:summarizing the development of the seismic imaging;overviewing the principles of three typical imaging methods,including common reflection surface(CRS)stack,migration-based Gaussian-beam migration(GBM)and reverse-time migration(RTM),and inversion-based least-squares reverse-time migration(LSRTM);analyzing the imaging capability of GBM,RTM and LSRTM to the special structures on three typical models and a land data set;outlooking the future perspectives of imaging methods.The main challenge of seismic imaging is to produce high-precision images for low-quality data,extremely deep reservoirs,and dual-complex structures.展开更多
Reverse Time Migration(RTM)Surface Ofset Gathers(SOGs)are demonstrated to deliver more superior residual dip information than ray-based approaches.It appears more powerful in complex geological settings,such as salt a...Reverse Time Migration(RTM)Surface Ofset Gathers(SOGs)are demonstrated to deliver more superior residual dip information than ray-based approaches.It appears more powerful in complex geological settings,such as salt areas.Still,the computational cost of constructing RTM SOGs is a big challenge in applying it to 3D feld data.To tackle this challenge,we propose a novel method using dips of local events as a guide for RTM gather interpolation.The residual-dip information of the SOGs is created by connecting local events from depth-domain to time-domain via ray tracing.The proposed method is validated by a synthetic experiment and a feld example.It mitigates the computational cost by an order of magnitude while producing comparable results as fully computed RTM SOGs.展开更多
基金supported by the Marine S&T Fund of Shandong Province for Pilot National Laboratory for Marine Science and Technology(No.2021QNLM020001)the Major Scientific and Technological Projects of Shandong Energy Group(No.SNKJ2022A06-R23)+2 种基金the Funds of Creative Research Groups of China(No.41821002)National Natural Science Foundation of China Outstanding Youth Science Fund Project(Overseas)(No.ZX20230152)the Major Scientific and Technological Projects of CNPC(No.ZD2019-183-003)。
文摘Compared with the transverse isotropic(TI)medium,the orthorhombic anisotropic medium has both horizontal and vertical symmetry axes and it can be approximated as a set of vertical fissures developed in a group of horizontal strata.Although the full-elastic orthorhombic anisotropic wave equation can accurately simulate seismic wave propagation in the underground media,a huge computational cost is required in seismic modeling,migration,and inversion.The conventional coupled pseudo-acoustic wave equations based on acoustic approximation can be used to significantly reduce the cost of calculation.However,these equations usually suffer from unwanted shear wave artifacts during wave propagation,and the presence of these artifacts can significantly degrade the imaging quality.To solve these problems,we derived a new pure P-wave equation for orthorhombic media that eliminates shear wave artifacts while compromising computational efficiency and accuracy.In addition,the derived equation involves pseudo-differential operators and it must be solved by 3D FFT algorithms.In order to reduce the number of 3D FFT,we utilized the finite difference and pseudo-spectral methods to conduct 3D forward modeling.Furthermore,we simplified the equation by using elliptic approximation and implemented 3D reverse-time migration(RTM).Forward modeling tests on several homogeneous and heterogeneous models confirm that the accuracy of the new equation is better than that of conventional methods.3D RTM imaging tests on three-layer and SEG/EAGE 3D salt models confirm that the ORT media have better imaging quality.
基金supported by National Key R&D Program of China (No. 2018YFA0702502)NSFC (Grant No. 41974142, 42074129, and 41674114)+1 种基金Science Foundation of China University of Petroleum (Beijing) (Grant No. 2462020YXZZ005)State Key Laboratory of Petroleum Resources and Prospecting (Grant No. PRP/indep-42012)。
文摘Least-squares reverse-time migration(LSRTM) formulates reverse-time migration(RTM) in the leastsquares inversion framework to obtain the optimal reflectivity image. It can generate images with more accurate amplitudes, higher resolution, and fewer artifacts than RTM. However, three problems still exist:(1) inversion can be dominated by strong events in the residual;(2) low-wavenumber artifacts in the gradient affect convergence speed and imaging results;(3) high-wavenumber noise is also amplified as iteration increases. To solve these three problems, we have improved LSRTM: firstly, we use Hubernorm as the objective function to emphasize the weak reflectors during the inversion;secondly, we adapt the de-primary imaging condition to remove the low-wavenumber artifacts above strong reflectors as well as the false high-wavenumber reflectors in the gradient;thirdly, we apply the L1-norm sparse constraint in the curvelet-domain as the regularization term to suppress the high-wavenumber migration noise. As the new inversion objective function contains the non-smooth L1-norm, we use a modified iterative soft thresholding(IST) method to update along the Polak-Ribie re conjugate-gradient direction by using a preconditioned non-linear conjugate-gradient(PNCG) method. The numerical examples,especially the Sigsbee2 A model, demonstrate that the Huber inversion-based RTM can generate highquality images by mitigating migration artifacts and improving the contribution of weak reflection events.
基金supported by the National Key R&D Program of China(2021YFA0716902)National Natural Science Foundation of China(42004119,42174156)+1 种基金the Fundamental Research Funds for the Central Universities,CHD(300102261306)the National Engineering Research Center of Offshore Oil and Gas Exploration,No.6 Courtyard,Taiyanggong South Street,Chaoyang District,Beijing,100028.
文摘Research on seismic anisotropy and attenuation plays a significant role in exploration geophysics. To enhance the imaging quality for complicated structures, we develop several effective improvements for anisotropic attenuation effects in reverse-time migration (Q-RTM) on surface and vertical seismic profiling (VSP) acquisition geometries. First, to suppress pseudo-shear wave artifact and numerical instability of the commonly used anisotropic pseudo-acoustic wave equations, an optimized pure P-wave dispersion relation is derived and the corresponding pure-mode wave equation is solved by combining the finite-difference and Possion methods. Second, a simplified anisotropic pure-mode visco-acoustic wave equation (PVAWE) based on standard linear solid model is established. Third, a time-dispersion correlation strategy is applied to improve the modeling accuracy. Fourth, we extend a target-oriented scheme to anisotropic attenuated modeling and imaging. Instead of the conventional wavefield modeling and RTM, the proposed approach can extract available wavefield information near the target regions and produce high imaging resolution for target structures. Last, both anisotropic surface and VSP Q-RTMs are executed by combining optimized PVAWE, time-dispersion correlation and target-oriented algorithm. Modeling examples demonstrate the advantages of our schemes. Moreover, our modified Q-compensated imaging workflow can be regarded as a supplement to the classical anisotropic RTM.
基金financially supported by National Basic Research Program of China(No.2011CB201100)
文摘Conventional seismic exploration method based on post-stack data usually fails to identify the distribution of fractured and caved carbonate reservoirs in the Tarim Basin,so the rich pre-stack information should be applied to the prediction of carbonate reservoirs.Amplitude-preserved seismic data processing is the foundation.In this paper,according to the feature of desert seismic data （including weak reflection,fast attenuation of high frequency components,strong coherent noises,low S/N and resolution）,a set of amplitude-preserved processing techniques is applied and a reasonable processing flow is formed to obtain the high quality data.After implementing a set of pre-stack amplitude-preserved processing,we test and define the kernel parameters of amplitude-preserved Kirchhoff PSTM （pre-stack time migration） and subsequent gathers processing,in order to obtain the amplitude-preserved gathers used to the isotropic pre-stack inversion for the identification of caved reservoirs.The AVO characteristics of obtained gathers fit well with the synthetic gathers from logging data,and it proves that the processing above is amplitudepreserved.The azimuthal processing techniques,including azimuth division and binning enlargement,are implemented for amplitude-preserved azimuthal gathers with the uniform fold.They can be used in the anisotropic inversion to detect effective fractures.The processing techniques and flows are applied to the field seismic data,and are proved available for providing the amplitude-preserved gathers for carbonate reservoir prediction in the Tarim Basin.
基金supported by seismic wave propagation and imaging(SWPI)group of China University of Petroleum(East China)supported by National Natural Science Foundation of China(42174138,41904101,42074133)+3 种基金Natural Science Foundation of Shandong Province(ZR2019QD004)Funds for the Central Universities(19CX02010A)the Major Scientific and Technological Projects of CNPC(ZD 2019183-003)Talent introduction fund of China University of Petroleum(East China)(20180041)。
文摘High-precision seismic imaging is the core task of seismic exploration,guaranteeing the accuracy of geophysical and geological interpretation.With the development of seismic exploration,the targets become more and more complex.Imaging on complex media such as subsalt,small-scale,steeply dipping and surface topography structures brings a great challenge to imaging techniques.Therefore,the seismic imaging methods range from stacking-to migration-to inversion-based imaging,and the imaging accuracy is becoming increasingly high.This review paper includes:summarizing the development of the seismic imaging;overviewing the principles of three typical imaging methods,including common reflection surface(CRS)stack,migration-based Gaussian-beam migration(GBM)and reverse-time migration(RTM),and inversion-based least-squares reverse-time migration(LSRTM);analyzing the imaging capability of GBM,RTM and LSRTM to the special structures on three typical models and a land data set;outlooking the future perspectives of imaging methods.The main challenge of seismic imaging is to produce high-precision images for low-quality data,extremely deep reservoirs,and dual-complex structures.
基金This study is jointly supported by the National Key R&D Program of China(2017YFC1500303 and 2020YFA0710604)the Science Foundation of China University of Petroleum,Beijing(2462019YJRC007 and 2462020YXZZ047)the Strategic Cooperation Technology Projects of CNPC and CUPB(ZLZX2020-05).
文摘Reverse Time Migration(RTM)Surface Ofset Gathers(SOGs)are demonstrated to deliver more superior residual dip information than ray-based approaches.It appears more powerful in complex geological settings,such as salt areas.Still,the computational cost of constructing RTM SOGs is a big challenge in applying it to 3D feld data.To tackle this challenge,we propose a novel method using dips of local events as a guide for RTM gather interpolation.The residual-dip information of the SOGs is created by connecting local events from depth-domain to time-domain via ray tracing.The proposed method is validated by a synthetic experiment and a feld example.It mitigates the computational cost by an order of magnitude while producing comparable results as fully computed RTM SOGs.
