三维弹性波数值模拟中改进的NPML研究

doi:10.6038/pg20160655

地球物理学进展 ›› 2016, Vol. 31 ›› Issue (6): 2762-2766.doi: 10.6038/pg20160655

三维弹性波数值模拟中改进的NPML研究

谢波^1,2, 熊章强^1,3, 张大洲^1,3

1. 中南大学地球科学与信息物理学院, 长沙 410083;
2. 四川大学锦城学院, 成都 611731;
3. 中南大学有色金属成矿预测与地质环境监测教育部重点实验室, 长沙 410083

收稿日期:2016-02-12 修回日期:2016-06-13 出版日期:2016-12-20 发布日期:2016-12-20
作者简介:谢波,男,1979年生,博士研究生,主要从事地震学及数值模拟相关工作.(E-mail:40252126@qq.com)
基金资助:
国家自然科学基金项目（41274123）资助.

Reformulated NPML in 3D elastic-wave numerical simulation

XIE Bo^1,2, XIONG Zhang-qiang^1,3, ZHANG Da-zhou^1,3

1. School of Geosciences and Info-Physics, Central South University, Changsha 410083, China;
2. JinCheng College of SiChuan University, ChengDu 611731, China;
3. Key Laboratory of Metallogenic Prediction of Nonferrous Metals and Geological Environment Monitoring, Ministry of Education, Changsha 410083, China

Received:2016-02-12 Revised:2016-06-13 Online:2016-12-20 Published:2016-12-20

摘要/Abstract

摘要：

在地震波传播数值模拟的过程中，需要使用吸收边界条件达到衰减人为边界反射的目的，其中完全匹配层（PML）吸收边界条件被认为是目前最理想的吸收边界条件.但在三维地震波传播数值模拟中直接应用常规的PML吸收边界条件存在编程实现复杂、计算资源需求高等不足之处.本文基于传统的NPML吸收边界条件，提出了一种改进的方法.该方法通过引入伸展函数和辅助变量，把PML介质中的三维弹性波动方程分解成正常计算项和衰减项，其中衰减项的计算只在PML区域内部进行，其值只与坐标轴的方向相关，可通过离散数值迭代求解.分析和试验结果表明，该方法避免了卷积运算，减少了辅助变量的个数，具有理论方程简洁、编程实现简单、计算资源需求少等优点，能更有效地解决三维地震波传播数值模拟中的吸收边界问题.

Abstract:

The absorbing boundary conditions to avoid the artificial boundary reflections must be required in the numerical simulation of seismic wave propagation. The perfectly matched layer has been generally proved to be the best available absorbing boundary condition. However, the classical PML absorbing boundary condition in 3D Elastic-Wave numerical simulation has some disadvantages:it is difficult to implement and needs more computing resources. In this paper, we present a reformulated method based on the traditional NPML absorbing boundary condition. we divided the 3D elastic-wave equation in PML media into the normal calculation and decay by introducing the extension function and auxiliary variables. The decay needs to be calculated only when the calculation area is located in the PML area. It's value relates to the axis direction, and it can be calculated through discrete numerical iteration computing. Analysis and experimental results show that the reformulated NPML given in this paper, since it avoids the convolution and decreases the auxiliary variables, has some advantages:the theoretical equation is more simple, easier to implement and needs less computing resources, so it is more effective than the traditional PML to resolve the absorbing boundary condition in 3D elastic-wave numerical simulation.

中图分类号:

P631.4

谢波, 熊章强, 张大洲. 2016. 三维弹性波数值模拟中改进的NPML研究. 地球物理学进展, 31(6): 2762-2766. doi: 10.6038/pg20160655.

XIE Bo, XIONG Zhang-qiang, ZHANG Da-zhou. 2016. Reformulated NPML in 3D elastic-wave numerical simulation. Progress in Geophysics. 31(6): 2762-2766. doi: 10.6038/pg20160655.

参考文献

Berenger J P. 1994. A perfectly matched layer for the absorption of electromagnetic waves[J]. Journal of Computational Physics, 114(2):185-200.
Berenger J P. 2002. Application of the CFS PML to the absorption of evanescent waves in waveguides[J]. IEEE Microwave and Wireless Components Letters, 12(6):218-220.
Chew W C, Weedon W H. 1994. A 3D perfectly matched medium from modified Maxwell's equations with stretched coordinates[J]. Microwave and Optical Technology Letters, 7(13):599-604.
Clayton R, Engquist B. 1977. Absorbing boundary conditions for acoustic and elastic wave equations[J]. Bulletin of the Seismological Society of America, 67(6):1529-1540.
Drossaert F H, Giannopoulos A. 2007. A nonsplit complex frequency-shifted PML based on recursive integration for FDTD modeling of elastic waves[J]. Geophysics, 72(2):T9-T17.
Higdon R L. 1987. Numerical absorbing boundary conditions for the wave equation[J]. Mathematics of Computation, 49(179):65-90.
Kuzuogulu M, Mittra R. 1996. Frequency dependence of the constitutive parameters of causal perfectly matched anisotropic absorbers[J]. IEEE Microwave and Guided Wave Letters, 6(12):447-449.
Liao Z P, Wong H L, Yang B, et al. 1984. A transmitting boundary for transient wave analyses[J]. Scientia Sinica Series A-Mathematical Physical Astronomical & Technical Sciences, 27(10):1063-1076.
Ma S, Liu P C. 2006. Modeling of the perfectly matched layer absorbing boundaries and intrinsic attenuation in explicit finite-element methods[J]. Bulletin of the Seismological Society of America, 96(5):1779-1794.
Marfurt K J. 1984. Accuracy of finite-difference and finite-element modeling of the scalar and elastic wave equations[J]. Geophysics, 49(5):533-549.
Qin Z, Lu M H, Zheng X D, et al. 2009. The implementation of an improved NPML absorbing boundary condition in elastic wave modeling[J]. Applied Geophysics, 6(2):113-121.
Wang S D. 2003. Absorbing boundary condition for acoustic wave equation by perfectly matched layer[J]. Oil Geophysical Prospecting (in Chinese), 38(1):31-34.
Wrenger J P. 2002. Numerical reflection from FDTD-PMLs:A comparison of the split PML with the unsplit and CFS PMLs[J]. IEEE Transactions on Antennas and Propagation, 50(3):258-265.
Xiong Z Q, Mao C Y. 2011. Improved unsplit PML boundary conditions in acoustic wave numerical simulation[J]. Oil Geophysical Prospecting (in Chinese), 46(1):35-39.
Xiong Z Q, Zhang D Z, Qin Z, et al. 2008. Boundary conditions and case analysis of numerical modeling of Rayleigh wave[J]. Journal of Central South University (Science and Technology) (in Chinese), 39(4):824-830.
Xu Y, Zhang J F. 2008. An irregular-grid perfectly matched layer absorbing boundary for seismic wave modeling[J]. Chinese Journal of Geophysics (in Chinese), 51(5):1520-1526, doi:10.3321/j.issn:0001-5733.2008.05.026.
Zhao H B, Wang X M, Wang D, et al. 2007. Applications of the boundary absorption using a perfectly matched layer for elastic wave simulation in poroelastic media[J]. Chinese Journal of Geophysics (in Chinese), 50(2):581-591, doi:10.3321/j.issn:0001-5733.2007.02.032.
王守东. 2003. 声波方程完全匹配层吸收边界[J]. 石油地球物理勘探, 38(1):31-34.
熊章强, 毛承英. 2011. 声波数值模拟中改进的非分裂式PML边界条件[J]. 石油地球物理勘探, 46(1):35-39.
熊章强, 张大洲, 秦臻,等. 2008. 瑞雷波数值模拟中的边界条件及模拟实例分析[J]. 中南大学学报(自然科学版), 39(4):824-830.
徐义, 张剑锋. 2008. 地震波数值模拟的非规则网格PML吸收边界[J]. 地球物理学报, 51(5):1520-1526, doi:10.3321/j.issn:0001-5733.2008.05.026.
赵海波, 王秀明, 王东,等. 2007. 完全匹配层吸收边界在孔隙介质弹性波模拟中的应用[J]. 地球物理学报, 50(2):581-591, doi:10.3321/j.issn:0001-5733.2007.02.032.

[1]	陈楠,徐雷鸣,宋志翔. HTI介质方位AVO（AVAZ）正演模拟软件研制与应用[J]. 地球物理学进展, 2017, 32(3): 1331-1338.
[2]	陈中山,张昭,孟红星,牛小军,张灯亮. 复杂地表条件下的三维观测系统优化设计[J]. 地球物理学进展, 2017, 32(2): 591-595.
[3]	凌云,杜向东,曹思远. 基于Zener线性体的黏弹孔隙介质衰减频散特征分析[J]. 地球物理学进展, 2017, 32(1): 205-209.
[4]	沈洪涛,郭乃川,秦童,于茜,王伟. 地质统计学反演技术在超薄储层预测中的应用[J]. 地球物理学进展, 2017, 32(1): 248-253.
[5]	张瑞,文晓涛,李世凯,王文化,刘敏兵. 分频蚂蚁追踪在识别深层小断层中的应用[J]. 地球物理学进展, 2017, 32(1): 350-356.
[6]	张晓普,林君,杨泓渊,赵玉江,张帅帅. 陆上地震数据采集系统通讯技术现状及展望[J]. 地球物理学进展, 2016, 31(3): 1390-1398.
[7]	韩东,胡向阳,邬兴威,刘坤岩,司朝年,王鹏,解丽慧,韩科龙. 基于地质统计学反演的缝洞储集体物性定量评价[J]. 地球物理学进展, 2016, 31(2): 655-661.
[8]	马见青,李庆春,王卫东,王美丁. 多分量地震的时频域瞬时极化分析和滤波[J]. 地球物理学进展, 2015, 30(6): 2723-2729.
[9]	吴海波,董守华,黄亚平,程彦. 基于地震属性的煤层火成岩侵入预测[J]. 地球物理学进展, 2015, 30(3): 1376-1381.
[10]	吴海波,董守华,黄亚平,祁雪梅. 煤层火成岩侵入的反射波特征研究与应用[J]. 地球物理学进展, 2014, 29(6): 2779-2784.
[11]	王玉凤,何兵寿,张省. 基于纵横波解耦的双相介质弹性波方程正演模拟及结果分析[J]. 地球物理学进展, 2013, 28(6): 3030-3039.
[12]	李向云,曹思远,杨燕,邱林林. 改进的正交多项式变换在地震资料去噪中的应用[J]. 地球物理学进展, 2013, 28(4): 2085-2093.
[13]	李振春,李娜,黄建平,田坤,李庆洋. 裂缝介质横波分裂时差影响因素定量研究[J]. 地球物理学进展, 2013, 28(1): 240-249.
[14]	赖晓玲,李松林,张先康,成双喜. 玛沁—兰州—靖边剖面壳幔复杂性的研究[J]. 地球物理学进展, 2001, 16(2): 65-72.
[15]	刘财,冯晅,何敏,杨宝俊,赵建国. 非线性波动理论在地震勘探中的应用综述[J]. 地球物理学进展, 2000, 15(4): 68-78.

三维弹性波数值模拟中改进的NPML研究

Reformulated NPML in 3D elastic-wave numerical simulation

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

多维度评价

本文评价