天然地震数值模拟与实例对比
|
摘要
为了解决长时间正演模拟的不稳定性问题,实现天然地震波场长时间数值模拟.在高阶有限差分数值模拟的基础上,给出了多轴向完全匹配吸收边界(M-PML)二维划分方式,讨论了M-PML吸收边界转换系数P的取值对其吸收能力与截断误差对数值模拟有效信号的影响.实现了二维弹性波高阶有限差分的长时间数值模拟.采用2013年12月16日湖北省巴东县地震及余震信息,对该地震数据进行了数值模拟,并与武汉、秭归两台站数据进行了分析和讨论,验证该正演模拟方法的正确性和有效性.为天然地震波场传播规律、成像和震相识别及震源定位等研究提供了基础.
In order to resolve the long-time forward modeling instability and realize the long-time numerical simulation of natural seismic wave field,based on staggered-grid high-order difference method,the 2-D model partitioning way of M-PML is provided,the value of the M-PML absorbing boundary conversion coefficient Pand its influence on the absorption capacity and the influence on effective signal made by truncation error are discussed,the 2Delastic wave long time numerical simulation using higher-order finite difference is realized.In addition,combining the data of the 5.5 magnitude earthquake which occurred in Hubei Badong on December 16,2013 and its aftershock recorded by Wuhan and Zigui Seismograph Stations,the practical significance of the forward modeling method is verified.The results can provide a basis for the research of the propagation law,imaging,phase identification and hypocentral location of natural seismic wave field.
In order to resolve the long-time forward modeling instability and realize the long-time numerical simulation of natural seismic wave field,based on staggered-grid high-order difference method,the 2-D model partitioning way of M-PML is provided,the value of the M-PML absorbing boundary conversion coefficient Pand its influence on the absorption capacity and the influence on effective signal made by truncation error are discussed,the 2Delastic wave long time numerical simulation using higher-order finite difference is realized.In addition,combining the data of the 5.5 magnitude earthquake which occurred in Hubei Badong on December 16,2013 and its aftershock recorded by Wuhan and Zigui Seismograph Stations,the practical significance of the forward modeling method is verified.The results can provide a basis for the research of the propagation law,imaging,phase identification and hypocentral location of natural seismic wave field.
引文
Berenger,J.P.,1994.A Perfectly Match Layer for the Absorption of Electromagnetic Waves.Journal of Computation Physics,114(2):185-200.doi:10.1006/jcph.1994.1159.
Collino,F.,Tsogka,C.,2001.Application of the PerfectlyMatched Absorbing Layer Model to the Linear Elastodynamic Problem in Anisotropic Heterogeneous Media.Geophysics,66(1):294-307.doi:10.1190/1.1444908
Dong,L.G.,Ma,Z.T.,Cao,J.Z.,2000,A Study on Stability ofthe Staggered-Grid High-Order Difference Method ofFirst-Order Elastic Wave Equation.Chinese Journal ofGeophys.,43(6):856-864(in Chinese with English abstract).
Meza-Fajardo,K.C.,Papageorgiou,A.S.,2008.A Nonconvolutional,Aplit-Field,Perfectly Matched Layer for WavePropagation in Isotropic and Anisotropic Elastic Media:Stability Analysis.Bulletin of the Seismological Society of America,98(4):1811-1836.doi:10.1785/0120070223
Mittet,R.,2002.Free-Surface Boundary Conditions for Elastic Staggered-Grid Modeling Schemes.Geophysics,67(5):1616-1623.doi:10.1190/1.1512752
Pei,Z.L.,2005.Numerical Simulation of Elastic Wave Equation in 3-D Isotropic Media with Staggered-Grid HighOrder Difference Method.GPP,44(4):308-315(inChinese).
Xu,Y.X.,Xia,J.H.,Miller,R.D.,2007.Numerical Investigation of Implementation of Air-Earth Boundary by Acoustic-Elastic Boundary Approach.Geophysics,72(5):SM147-SM153.doi:10.1190/1.2753831
Zhao,X.,Li,Q.,Cai,J.A.,2007.On Minimum 1-D VelocityModel Applied in Three Gorges Reservoir Area.Journal of Geodesy and Geodynamics,27(Spec.):1-7(inChinese with English abstract).
董良国,马在田,曹景忠,2000.一阶弹性波方程交错网格高阶差分解法稳定性研究.地球物理学报,43(6):856-864.
裴正林,2005.三维各向同性介质弹性波方程交错网格高阶有限差分法模拟.石油物探,44(4):308-315.
赵旭,李强,蔡晋安,2007.三峡库首区最小一维速度模型研究.大地测量与地球动力学,27(专刊):1-7.
Collino,F.,Tsogka,C.,2001.Application of the PerfectlyMatched Absorbing Layer Model to the Linear Elastodynamic Problem in Anisotropic Heterogeneous Media.Geophysics,66(1):294-307.doi:10.1190/1.1444908
Dong,L.G.,Ma,Z.T.,Cao,J.Z.,2000,A Study on Stability ofthe Staggered-Grid High-Order Difference Method ofFirst-Order Elastic Wave Equation.Chinese Journal ofGeophys.,43(6):856-864(in Chinese with English abstract).
Meza-Fajardo,K.C.,Papageorgiou,A.S.,2008.A Nonconvolutional,Aplit-Field,Perfectly Matched Layer for WavePropagation in Isotropic and Anisotropic Elastic Media:Stability Analysis.Bulletin of the Seismological Society of America,98(4):1811-1836.doi:10.1785/0120070223
Mittet,R.,2002.Free-Surface Boundary Conditions for Elastic Staggered-Grid Modeling Schemes.Geophysics,67(5):1616-1623.doi:10.1190/1.1512752
Pei,Z.L.,2005.Numerical Simulation of Elastic Wave Equation in 3-D Isotropic Media with Staggered-Grid HighOrder Difference Method.GPP,44(4):308-315(inChinese).
Xu,Y.X.,Xia,J.H.,Miller,R.D.,2007.Numerical Investigation of Implementation of Air-Earth Boundary by Acoustic-Elastic Boundary Approach.Geophysics,72(5):SM147-SM153.doi:10.1190/1.2753831
Zhao,X.,Li,Q.,Cai,J.A.,2007.On Minimum 1-D VelocityModel Applied in Three Gorges Reservoir Area.Journal of Geodesy and Geodynamics,27(Spec.):1-7(inChinese with English abstract).
董良国,马在田,曹景忠,2000.一阶弹性波方程交错网格高阶差分解法稳定性研究.地球物理学报,43(6):856-864.
裴正林,2005.三维各向同性介质弹性波方程交错网格高阶有限差分法模拟.石油物探,44(4):308-315.
赵旭,李强,蔡晋安,2007.三峡库首区最小一维速度模型研究.大地测量与地球动力学,27(专刊):1-7.
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心