基于保辛算法的声波叠前逆时偏移
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
叠前逆时偏移是目前成像精度最高的地震偏移方法之一,其实现过程中的一个重要步骤是数值求解全波方程,所以快速有效求解全波方程的数值算法对逆时偏移至关重要.四阶近似解析辛可分Runge-Kutta(NSPRK)方法是近年发展的一种具有高效率、高精度的数值求解波动方程的保辛差分方法,能在粗网格条件下有效压制数值频散,从而提高计算效率,节省计算机内存需求量.本文利用四阶NSPRK方法构造的基本思想,发展了具有六阶空间精度的NSPRK方法,并对新的六阶NSPRK方法进行了详细的稳定性和数值频散分析,以及计算效率比较和波场模拟.同时将该方法用于声波叠前逆时偏移中,得到一种时间上保辛、空间具有六阶精度、低数值频散、可应用大步长进行波场延拓并能长时计算的叠前逆时偏移方法,对Sigsbee2B模型进行了偏移成像,并和四阶NSPRK方法、传统的六阶差分方法、四阶Lax-Wendroff correction(LWC)方法进行了对比.数值结果表明,基于六阶NSPRK方法的叠前逆时偏移能得到更好的成像结果,是一种优于四阶NSPRK方法、传统的六阶差分方法、四阶LWC叠前逆时偏移的方法,尤其是在粗网格情况下具有更明显的优越性.
Prestack reverse time migration(RTM)is currently one of the most accurate seismic migration imaging methods.An important step of its implementation is to solve the full-wave equation,so the numerical algorithm of solving the full wave equation quickly and effectively is crucial for reverse-time migration.The nearly analytic symplectic partitioned Runge-Kutta(NSPRK)method is a high-efficiency and high-precision numerical symplectic difference method developed in recent years,which can effectively suppress the numerical dispersion when a coarse grid is used.So it can increase computation efficiency and save computer memory demand.In this article,we develop a sixth-order NSPRK method based on the fourth-order NSPRK method,analyze the stability and numerical dispersion of the new sixth-order NSPRK method,compare its calculation efficiency and conduct wave field simulation.Then we apply it to the acoustic prestack reverse time migration,obtaining a sixth-order symplectic prestack reverse time migration method with low numerical dispersion,which can be applied to wavefield extrapolation with large steps and compute in long time.We get migration images of the Sigsbee2 B model,and comparethem to the fourth-order NSPRK method,the sixth-order conventional finite difference method,and the fourth-order Lax-Wendroff correction method(LWC).Numerical results show that the method based on the NSPRK prestack reverse time migration can get better imaging results than the fourth-order NSPRK method,the sixth-order conventional finite difference method,and the fourth-order LWC prestack reverse time migration method,especially in the case of coarse meshes.
引文
Baysal E,Kosloff D D,Sherwood J W C.1983.Reverse timemigration.Geophysics,48(11):1514-1524,doi:10.1190/1.1441434.
    Chang W F,McMechan G A.1986.Reverse-time migration ofoffset vertical seismic profiling data using the excitation timeimaging condition.Geophysics,51(1):67-84,doi:10.1190/1.1442041.
    Chang W F,McMechan G A.1987.Elastic reverse-time migration.Geophysics,52(10):1365-1375,doi:10.1190/1.1442249.
    Chang W F,McMechan G A.1990.3Dacoustic prestack reversetime migration.Geophysical Prospecting,38(7):737-755,doi:10.1111/j.1365-2478.1990.tb01872.x.
    Chang W F,McMeehan G A.1994.3-D elastic prestack,reversetime depth migration.Geophysics,59(4):597-609,doi:10.1190/1.1443620.
    Chattopadhyay S,McMeehan G A.2008.Imaging conditions forprestack reverse-time migration.Geophysics,73(3):S81-S89,doi:10.1190/1.2903822.
    Dablain M A.1986.The application of high-order differencing to thescalar wave equation.Geophysics,51(1):54-66,doi:10.1190/1.1442040.
    Di Q Y,Wang M Y.1997.The study of finite element reverse-timemigration for elastic wave.Chinese J.Geophys.(in Chinese),40(4):570-579.
    Du Q Z,Qin T.2009.Multicomponent prestack reverse-timemigration of elastic waves in transverse isotropic medium.Chinese J.Geophys.(in Chinese),52(3):801-807.
    Duveneck E,Bakker P M.2011.Stable P-wave modeling forreverse-time migration in tilted TI media.Geophysics,76(2):S65-S75,doi:10.1190/1.3533964.
    Feng K,Qin M Z.2003.Symplectic Geometric Algorithms forHamiltonian Systems(in Chinese).Hangzhou:Zhejiang Science&Technology Press,185-206,271-278.
    Komatitsch D,Tromp J.2003.A perfectly matched layer absorbingboundary condition for the second-order seismic wave equation.Geophysical Journal International,154(1):146-153,doi:10.1046/j.1365-246X.2003.01950.x.
    Kondoh Y,Hosaka Y,Ishii K.1994.Kernel optimum nearlyanalytical discretization(KOND)algorithm applied to parabolicand hyperbolic equations.Computers&Mathematics withApplications,27(3):59-90,doi:10.1016/0898-1221(94)90047-7.
    Li X F,Li Y Q,Zhang M G,et al.2011.Scalar seismic-waveequation modeling by a multisymplectic discrete singularconvolution differentiator method.Bulletin of the SeismologicalSociety of America,101(4):1710-1718,doi:10.1785/0120100266.
    Li X F,Wang W S,Lu M W,et al.2012.Structure-preservingmodelling of elastic waves:a symplectic discrete singularconvolution differentiator method.Geophysical Journal International,188(3):1382-1392,doi:10.1111/j.1365-246X.2011.05344.x.
    Li Y Q,Li X F,Zhu T.2011.The seismic scalar wave fieldmodeling by symplectic scheme and singular kernel convolutionaldifferentiator.Chinese J.Geophys.(in Chinese),54(7):1827-1834,doi:10.3969/j.issn.0001-5733.2011.07.016.
    Liu H W,Li B,Liu H,et al.2010.The algorithm of high orderfinite difference pre-stack reverse time migration and GPUimplementation.Chinese J.Geophys.(in Chinese),53(7):1725-1733,doi:10.3969/j.issn.0001-5733.2010.07.024.
    Liu Y.2013.Globally optimal finite-difference schemes based onleast squares.Geophysics,78(4):T113-T132,doi:10.1190/GEO2012-0480.1.
    Liu Y K,Chang X,Jin D G,et al.2011.Reverse time migration ofmultiples for subsalt imaging.Geophysics,76(5):WB209-WB216,doi:10.1190/geo2010-0312.1.
    Ma X,Yang D H,Zhang J H.2010.Symplectic partitioned RungeKutta method for solving the acoustic wave equation.ChineseJ.Geophys.(in Chinese),53(8):1993-2003,doi:10.3969/j.issn.0001-5733.2010.08.026.
    McMechan G A.1983.Migration by extrapolation of timedependent boundary values.Geophysical Prospecting,31(3):413-420,doi:10.1111/j.1365-2478.1983.tb01060.x.
    Moczo P,Kristek J,Halada L.2000.3Dfourth-order staggeredgrid finite-difference schemes:Stability and grid dispersion.Bulletin of the Seismological Society of America,90(3):587-603,doi:10.1785/0119990119.
    Sun G.1997.A class of explicitly symplectic schemes for waveequations.Mathematica Numerica Sinica(in Chinese),(1):1-10.
    Wang L,Yang D H,Deng X Y.2009.A WNAD method forseismic stress-field modeling in heterogeneous media.ChineseJ.Geophys.(in Chinese),52(6):1126-1535,doi:10.3969/j.issn.0001-5733.2009.06.014.
    Whitmore N D.1983.Iterative depth migration by backward timepropagation.//The 53rd SEG Annual meeting ExpandedAbstracts.SEG-1983-0382.Las Vegas,Nevada.
    Wu W J,Lines L R,Lu H X.1996.Analysis of higher-order,finite-difference schemes in 3-D reverse-time migration.Geophysics,61(3):845-856,doi:10.1190/1.1444009.
    Xu K,Zhou B,McMechan G A.2010.Implementation of prestackreverse time migration using frequency-domain extrapolation.Geophysics,75(2):S61-S72,doi:10.1190/1.3339386.
    Yan H Y,Liu Y.2013.Acoustic prestack reverse time migrationusing the adaptive high-order finite-difference method in timespace domain.Chinese J.Geophys.(in Chinese),56(3):971-984,doi:10.6038/cjg20130325.
    Yang D H,Teng J W,Zhang Z J,et al.2003.A nearly analyticdiscrete method for acoustic and elastic wave equations inanisotropic media.Bulletin of the Seismological Society ofAmerica,93(2):882-890,doi:10.1785/0120020125.
    Yang D H,Peng J M,Lu M,et al.2006.Optimal nearly analyticdiscrete approximation to the scalar wave equation.Bulletin ofthe Seismological Society of America,96(3):1114-1130,doi:10.1785/0120050080.
    Yang D H,Song G J,Lu M.2007a.Optimally accurate nearlyanalytic discrete scheme for wave-field simulation in 3Danisotropic media.Bulletin of the Seismological Society ofAmerica,97(5):1557-1569,doi:10.1785/0120060209.
    Yang D H,Song G J,Chen S,et al.2007b.An improved nearlyanalytical discrete method:an efficient tool to simulate theseismic response of 2-D porous structures.Journal of Geophysics andEngineering,4(1):40-52,doi:10.1088/1742-2132/4/1/006.
    Yang D H,Wang N,Chen S,et al.2009.An explicit method basedon the implicit Runge-Kutta algorithm for solving waveequations.Bulletin of the Seismological Society of America,99(6):3340-3354,doi:10.1785/0120080346.
    Yang D H,Tong P,Deng X Y.2012.A central difference methodwith low numerical dispersion for solving the scalar waveequation.Geophysical Prospecting,60(5):885-905,doi:10.1111/j.1365-2478.2011.01033.x.
    Yoon K,Shin C,Suh S,et al.2003.3Dreverse time migrationusing the acoustic wave equation:An experience with the SEG/EAGE data set.The Leading Edge,22(1):38-41,doi:10.1190/1.1542754.
    Zhang J H,Yao Z X.2013.Optimized finite-difference operator forbroadband seismic wave modeling.Geophysics,78(1):A13-A18,doi:10.1190/GEO2012-0277.1.
    Zhang Y,Sun J.2009.Practical issues in reverse time migration:true amplitude gathers,noise removal and harmonic sourceencoding.First Break,26:29-35.
    Zhang Y,Zhang H Z,Zhang G Q.2011.A stable TTI reverse timemigration and its implementation.Geophysics,76(3):WA3-WA11,doi:10.1190/1.3533964.
    底青云,王妙月.1997.弹性波有限元逆时偏移技术研究.地球物理学报,40(4):570-579.
    杜启振,秦童.2009.横向各向同性介质弹性波多分量叠前逆时偏移.地球物理学报,52(3):801-807.
    冯康,秦孟兆.2003.哈密尔顿系统的辛几何算法.杭州:浙江科学技术出版社,185-206,271-278.
    李一琼,李小凡,朱童.2011.基于辛格式奇异核褶积微分算子的地震标量波场模拟.地球物理学报,54(7):1827-1834,doi:10.3969/j.issn.0001-5733.2011.07.016.
    刘红伟,李博,刘洪等.2010.地震叠前逆时偏移高阶有限差分算法及GPU实现.地球物理学报,53(7):1725-1733,doi:10.3969/j.issn.0001-5733.2010.07.024.
    马啸,杨顶辉,张锦华.2010.求解声波方程的辛可分Runge-Kutta方法.地球物理学报,53(8):1993-2003,doi:10.3969/j.issn.0001-5733.2010.08.026.
    孙耿.1997.波动方程的一类显式辛格式.计算数学,(1):1-10.
    王磊,杨顶辉,邓小英.2009.非均匀介质中地震波应力场的WNAD方法及其数值模拟.地球物理学报,52(6):1126-1535,doi:10.3969/j.issn.0001-5733.2009.06.014.
    严红勇,刘洋.2013.基于时空域自适应高阶有限差分的声波叠前逆时偏移.地球物理学报,56(3):971-984,doi:10.6038/cjg20130325.

版权所有:© 2023 中国地质图书馆 中国地质调查局地学文献中心