弹性介质中二维波动方程数值模拟

赵英杰+王智亮+孙雄伟+田继东+何丹

摘要：本文对于弹性介质中二维波动方程，介绍了求解地震波方程的ONADM方法，然后用该方法对二维波动方程在均匀介质和水平层状速度模型进行数值模拟，结果表明，在均匀介质中得到的了清晰、无频散的波场快照，水平层状速度模型的地震记录清晰，直达波和反射波区分明显。可见，ONADM方法是一种有效的地震波模拟算法，对提高地震勘探精度，提高地震防震减灾具有一定的指导意义。

关键词：波动方程；ONADM；波场快照

4. 结论

本文在采用NAD类算法，在粗网格的条件下，对均匀介质和水平层状介质进行模拟，获得了清晰的波场快照，表明该方法能够较好地压制数值频散，得到清晰的地震波场快照和地震记录。该方法对我们认识地下介质构造，提高地震勘探精度获得更加准确的地震成像结果具有重要意义。同时，运用该方法进行地震波场模拟对防震减灾也具有一定的指导意义。

参考文献：

[1] 宋国杰.三维弹性波方程的改进近似解析离散化方法及波场模拟.硕士毕业论文.

[2] 杨桂通，张善元.弹性动力学.中国铁道出版社， 1988.

[3] Yang D H， Teng J M， Zhang Z J， et al. A nearly analytic discrete method for acoustic and elastic wave equations in anisotropic media. Bull.Seis.Soc.Am.， 2003，93（2）：882—890.

[4] 牛濱华. 半空间介质与地震波传播[M]. 石油工业出版社， 2002.

[5] 宋亚辉， 孙琪皓， 聂振邦，等. 裂隙介质下天然气水合物CO2-CH4相平衡实验研究[J]. 西部资源， 2016（1）：25-26.

[6] 孔然， 唐喜， 谢梦谦. 地震勘探中的低信噪比处理分析[J]. 西部资源， 2016（5）.

致谢：

1.大北区块构造数值模拟研究，中石油集团股份有限公司科技公关重大专项，2011zx05046—3.

2.41674141，分数阶粘性地震波场的高精度有限差分算法研究，国家自然科学基金面上项目.

3.西南石油大学地震波传播与成像（2015CXTD07）青年科技创新团队.

4.X151514KLX0项目名称：各向异性弹性介质中的拟声波控制方程逼近研究.