地震波交错网格高阶差分数值模拟研究

周学明; 李庆春; 李志华; 马 婷

PDF(1164 KB)

铁道工程学报 ›› 2011, Vol. 28 ›› Issue (8) : 1-6.

路网研究

地震波交错网格高阶差分数值模拟研究

周学明1，李庆春2，李志华1，马婷2

作者信息 +

Research on Staggered － grid High － order Finite － difference Numerical Simulation of Seismic Wave

ZHOU Xue － ming1，LI Qing － chun2，LI Zhi － hua1，MA Ting2

Author information +

文章历史 +

摘要

研究目的: 地震波数值模拟技术是勘探地球物理学中的重要组成部分，研究通过弹性波一阶速度－应力方程，采用交错网格高阶有限差分法实现了地震波在各向同性介质中的高精度的数值模拟，并采用完全匹配层( PML) 吸收边界来消除边界反射，可取得较好的效果。

研究结论: 通过模型的正演计算和复杂模型的处理结果表明，交错网格高阶有限差分法数值模拟是一种快速有效的地震波数值模拟方法。可用于为地震勘探选择采集参数，优化观测系统提供依据，以提高采集数据的质量，有效地解决铁路勘察过程中的地质问题。

Abstract

Research purposes: The technology of seismic wave numerical simulation is an important part in the geophysical prospecting． In this paper，based on the elastic － velocity － stress equation，the numerical simulation of seismic wave propagation in isotropic media is carried out with staggered － grid high － order finite － difference method， and the perfectly matched layer( PML) is used to absorpt boundary to eliminate the boundary reflection． The result of simulation is rather perfect．

Research conclusions: The results of forward computation and the seismic process of a complex model show the staggered － grid high － order numerical simulation is a rapid and efficient simulation method for seismic wave propagation． This method can be used for the selection of acquisition parameters in seismic exploration and for the optimization of the observation system to improve the quality of the acquisition data and solve the geological problems in railway exploration．

导出引用

周学明, 李庆春, 李志华, 等. 地震波交错网格高阶差分数值模拟研究[J]. 铁道工程学报, 2011, 28(8): 1-6

ZHOU Xue － ming, LI Qing － chun, LI Zhi － hua, et al. Research on Staggered － grid High － order Finite － difference Numerical Simulation of Seismic Wave[J]. Journal of Railway Engineering Society, 2011, 28(8): 1-6

中图分类号： P631

参考文献

［1］ Virieux，J． SH － wave Propagation in Heterogenerous Media: Velocity － stress Finite － difference Method［J］． Geophysics，1984( 11) : 1933 － 1957．

［2］ Virieux，J． P － SV Wave Propagation in Heterogeneous Media: Velocity － stress Finite － difference Method［J］． Geophysic，1986( 4) : 889 － 901．

［3］ Lavender，A． R． Fourth － order Finite － difference P － SV Seismograms［J］． Geophysics，1988 ( 11 ) : 1425 － 1436．

［4］ Dablain，M． A． The Application of High － order Differnceing to the Scalar Wave Equation ［J ］． Geophysics，1986( 1) : 54 － 66．

［5］ Collino，F．，Tsogka，C． Application of the Perfectly Matched Absorbing Layer Model to the Linear Elastodynamic Problem in Anisotropic Heterogeneous Media［J］． Geophysics，2001( 1) : 294 － 307．

［6］董良国，马在田，曹景忠，等．一阶弹性波方程交错网格高阶差分解法［J］．地球物理学报，2000( 3) : 411 －419．

Dong Liangguo，Ma Zaitian，Cao Jingzhong，etc． A Staggered － grid High － order Difference Method of One － order Elastic Wave Equation［J］． Chinese Journal of Geophysics( in Chinese) ，2000( 3) : 411 － 419．

［7］董良国，马在田，曹景忠．一阶弹性波方程交错网格高阶差分解法稳定性研究［J］．地球物理学报，2000( 6) : 856 － 864．

Dong Liangguo，Ma Zaitian，Cao Jingzhong． A Study on Stablity of the Staggered － grid High － order Difference Method of First － order Elastic Wave Equation ［J］． Chinese Journal of Geophysics ( in Chinese) ，2000( 6) : 856 － 864．

［8］周海滨，李志华．困难地形条件下地震反射在煤矿采空区的技术研究［J］．铁道工程学报，2010( 2) : 18 －21．

Zhou Haibin，Li Zhihua． Technical Research on the Seismic Reflection in Coal Seam Gob under Difficult Terrain Condition［J］． Journal of Railway Engineering Society，2010( 2) : 18 － 21．

［9］李志华，朱旭东．复杂地质隧道的综合物探技术［J］．铁道工程学报，2008( 8) : 58 －61．

Li Zhihua，Zhu Xudong． The Application of Integrated Geophysical Exploration Technology in the Tunnel Under the Complicated Geological Conditions ［J］． Journal of Railway Engineering Society，2008( 8) : 58 － 61．

［10］裴正林，牟永光．非均匀介质地震波传播交错网格高阶有限差分法模拟［J］．石油大学学报( 自然科学版) ，2003 ( 6) : 17 －21

Pei Zhenglin，Mou Yongguang． Staggered － grid High － order Finit － dfference Modeling for Seismic Propagation in HeterogEneous Media［J］． Journal of the Unversity of Petroleum，China，2003( 6) : 17 － 21．

［11］王守东，声波方程完全匹层吸收边界［J］．石油地球物理勘探，2003( 1) : 31 －34

Wang Shoudong． Perfectly Matched Layer Absorbing Condition for Acoustic Equation［J］． Oil Geophysical Prospecting，2003( 1) : 31 － 34