高阶交错网格有限差分法的地震波场数值模拟
|
摘要
基于地震波波动方程和有限差分法,针对多种不同的复杂地质模型,进行高阶交错网格地震波场模拟仿真研究.首先对地震波波动方程原始形式进行推演,得到一阶速度应力弹性波方程,然后用时空二元高阶交错差分网格的方法对其进行离散化,最后对三层背斜模型、断层模型和含油模型进行了数值模拟实验,获取相应的波场快照和地震道图像进行对比分析,并讨论了差分阶次对于仿真精度的影响.仿真结果表明:算法计算精度准确,计算效率高,适合于模拟非均匀介质、复杂构造和复杂地质体的地震波场.
Based on seismic wave equation and finite-differential method,seismic wave field numerical modeling with high-order staggered mesh is studied for different complicated geographical mode.First,seismic wave equation is derived to obtain one-order stress-velocity equation,then,the equation is discretized in time and space by staggered mesh finite-differential method.This method is used to simulate three-level strata sloping downward model,fault model and oil layer model and the differential order of the simulation accuracy is discussed.The results indicated this method is suit for simulating non-homogeneous media and complicated geographical mode with high accuracy and efficiency.
Based on seismic wave equation and finite-differential method,seismic wave field numerical modeling with high-order staggered mesh is studied for different complicated geographical mode.First,seismic wave equation is derived to obtain one-order stress-velocity equation,then,the equation is discretized in time and space by staggered mesh finite-differential method.This method is used to simulate three-level strata sloping downward model,fault model and oil layer model and the differential order of the simulation accuracy is discussed.The results indicated this method is suit for simulating non-homogeneous media and complicated geographical mode with high accuracy and efficiency.
引文
[1]Sen M K,Stoffa P L.Global optimization methods in geo-physical inversion[M].Amsterdam:Elsevier,1995.
[2]Jervis M,Sen M K,Stoffa P L.Prestack migration veloci-ty estimation using nonlinear methods[J].Geophysics,1996,61(1):138-150.
[3]Liu Z P,Liu J Q.Seismic-controlled nonlinear extrapola-tion of well parameters using neural networks[J].Geo-physics,1998,63(6):2035-2041.
[4]崔炯成,吴媚,李维新,等.高分辨率非线性地震波阻抗反演方法和应用[J].石油物探,2007,46(1):1-12.CUI Jiongcheng,WU Mei,LI Weixin,et al.Nonlinearhigh-resolution impedance inversion method and its applica-tion[J].Geophysical Prospecting for Petroleum,2007,46(1):1-12.(in Chinese)
[5]韩小俊,赵宪生,王建民.波阻抗反演结果的影响因素分析[J].大庆地质石油与开发,2005,24(4):92-93.HAN Xiaojun,ZHAO Xiansheng,WANG Jianmin.Anal-ysis of influential factors for wave impedance inversion re-sults[J].Petroleum Geology&Oilfield Development inDaqing[J].2005,24(4):92-93.(in Chinese)
[6]刘诚,郭科.基于模式识别的非线性地震资料解释方法[J].物探化探计算技术,2008,30(1):27-29.LIU Cheng,GUO Ke.Non-linear method of quake dataexplanation based on pattern recognition[J].ComputingTechniques for Geophysical and Geochemical Exploration,2008,30(1):27-29.(in Chinese)
[7]Hicks G J.Arbitrary source and receiver positioning in fi-nite-difference schemes using Kaiser windowed sinc func-tions[J].Geophysics,2002,67(1):156-166.
[8]刘喜武,刘洪.地震盲反褶积综述[J].地球物理学进展,2003,18(2):203-209.LIU Xiwu,LIU Hong.Survey on seismic blind deconvolu-tion[J].Progress In Geophysics,2003,18(2):203-209.(in Chinese)
[9]Barnes A E.Weighted average seismic attributes[J].Geo-physics,2000,65(1):275-285.
[10]Smith D N,Ferguson J F.Constrained inversion of seis-mic refraction data using the controlled random search[J].Geophysics,2000,65(5):1622-1630.
[11]左莹.基于高阶交错网格的有限差分地震波场数值模拟[D].西安:长安大学,2009.ZUO Ying.Seismic wave field numerical modeling basedon high-order finite-differentical method[D].Xi’an:Chang’an University,2009.(in Chinese)
[12]Collino F,Tsogka C.Application of the perfectly matchedabsorbing layer model to the linear elastodynamic problemin anisotropic heterogeneous media[J].Geophysics,2001,66(1):294-307.
[13]Carcione J,Herman G C,Kroode A P E.Seismic model-ing[J].Geophysics,2002,67(4):1304-1325.
[14]贾跃玮.基于MATLAB的地震正演模型实现[J].石油工业计算机应用,2009,17(2):10-13.JIA Yuewei.The seismic forward model based on Matlab[J].Computer Applications of Petroleum,2009,17(2):10-13.(in Chinese)
[2]Jervis M,Sen M K,Stoffa P L.Prestack migration veloci-ty estimation using nonlinear methods[J].Geophysics,1996,61(1):138-150.
[3]Liu Z P,Liu J Q.Seismic-controlled nonlinear extrapola-tion of well parameters using neural networks[J].Geo-physics,1998,63(6):2035-2041.
[4]崔炯成,吴媚,李维新,等.高分辨率非线性地震波阻抗反演方法和应用[J].石油物探,2007,46(1):1-12.CUI Jiongcheng,WU Mei,LI Weixin,et al.Nonlinearhigh-resolution impedance inversion method and its applica-tion[J].Geophysical Prospecting for Petroleum,2007,46(1):1-12.(in Chinese)
[5]韩小俊,赵宪生,王建民.波阻抗反演结果的影响因素分析[J].大庆地质石油与开发,2005,24(4):92-93.HAN Xiaojun,ZHAO Xiansheng,WANG Jianmin.Anal-ysis of influential factors for wave impedance inversion re-sults[J].Petroleum Geology&Oilfield Development inDaqing[J].2005,24(4):92-93.(in Chinese)
[6]刘诚,郭科.基于模式识别的非线性地震资料解释方法[J].物探化探计算技术,2008,30(1):27-29.LIU Cheng,GUO Ke.Non-linear method of quake dataexplanation based on pattern recognition[J].ComputingTechniques for Geophysical and Geochemical Exploration,2008,30(1):27-29.(in Chinese)
[7]Hicks G J.Arbitrary source and receiver positioning in fi-nite-difference schemes using Kaiser windowed sinc func-tions[J].Geophysics,2002,67(1):156-166.
[8]刘喜武,刘洪.地震盲反褶积综述[J].地球物理学进展,2003,18(2):203-209.LIU Xiwu,LIU Hong.Survey on seismic blind deconvolu-tion[J].Progress In Geophysics,2003,18(2):203-209.(in Chinese)
[9]Barnes A E.Weighted average seismic attributes[J].Geo-physics,2000,65(1):275-285.
[10]Smith D N,Ferguson J F.Constrained inversion of seis-mic refraction data using the controlled random search[J].Geophysics,2000,65(5):1622-1630.
[11]左莹.基于高阶交错网格的有限差分地震波场数值模拟[D].西安:长安大学,2009.ZUO Ying.Seismic wave field numerical modeling basedon high-order finite-differentical method[D].Xi’an:Chang’an University,2009.(in Chinese)
[12]Collino F,Tsogka C.Application of the perfectly matchedabsorbing layer model to the linear elastodynamic problemin anisotropic heterogeneous media[J].Geophysics,2001,66(1):294-307.
[13]Carcione J,Herman G C,Kroode A P E.Seismic model-ing[J].Geophysics,2002,67(4):1304-1325.
[14]贾跃玮.基于MATLAB的地震正演模型实现[J].石油工业计算机应用,2009,17(2):10-13.JIA Yuewei.The seismic forward model based on Matlab[J].Computer Applications of Petroleum,2009,17(2):10-13.(in Chinese)
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心