吉林大学学报(地球科学版)

• 地球探测与信息技术 • 上一篇    下一篇

辅助微分方程完全匹配层在声波方程数值模拟中的应用

赵建国1,史瑞其2,陈竞一3,赵维俊4,王宏斌5,潘建国5   

  1. 1.中国石油大学(北京)地球物理与信息工程学院/油气资源与探测国家重点实验室/CNPC物探重点实验室, 北京102249;
    2.中国海洋石油研究总院, 北京100027;
    3.塔尔萨大学地球科学学院, 美国 塔尔萨74104;
    4.沈阳地质矿产研究所, 沈阳110034;
    5.中石油勘探开发研究院西北分院,兰州730000
  • 收稿日期:2013-05-06 出版日期:2014-03-26 发布日期:2014-03-26
  • 作者简介:赵建国(1976-),男,副教授,博士,主要从事电磁勘探、岩石物理、地震模拟研究,E-mail:jgzhao761215@yahoo.com.cn
  • 基金资助:

    国家自然科学基金项目(41274138);中国石油大学(北京)基础研究基金项目(KYJJ2012-05-25)

Application of Auxiliary Differential Equation Perfectly Matched Layer for Numerical Modeling of Acoustic Wave Equations

Zhao Jianguo1,Shi Ruiqi2,Chen Jingyi3,Zhao Weijun4,Wang Hongbin5,Pan Jianguo5   

  1. 1.College of Geophysics and Information Engineering, China University Petroleum, Beijing/State Key Laboratory of Oil Resources and Prospecting/ CNPC Key Lab of Geophysical Prospecting, Beijing102249, China;
    2.Research Institute, China National Offshore Oil Corp., Beijing100027, China;
    3.Department of Geoscience, University of Tulsa,Oklahoma74104, USA;
    4.Shenyang Institute of Geology and Mineral Resources, Shenyang110034, China;
    5.Northwest Institute of Oil Exploration and Development, CNPC, Lanzhou730000, China
  • Received:2013-05-06 Online:2014-03-26 Published:2014-03-26

摘要:

在地震波数值模拟中,为提高算法精度,需要使用高阶时间更新格式,而普通的非分裂完全匹配层(PML)吸收边界局限于低阶时间格式。辅助微分方程完全匹配层(ADEPML)是一种可以适应任意阶时间格式的非分裂完全匹配层技术,且可以直接应用复频移拉伸算子以提高PML在高角度入射时的效果。作者将ADEPML应用于声波方程四阶RungeKutta时间格式的数值模拟中,对其吸收效能进行了检验。数值模拟表明,复频移ADEPML在高角度入射时表现优于非复频移ADEPML。另外,不同辅助变量更新格式的吸收效果存在微小差异,显格式下计算结果与解析解吻合较好。长时间能量衰减计算表明ADEPML可以稳定至2 × 105时间步。

关键词: 声波方程, 完全匹配层, 辅助微分方程, 任意阶时间格式, 数值模拟

Abstract:

In seismic wave modeling, a high-order time marching scheme is required in order to increase the accuracy of the algorithm. However, ordinary unsplit perfectly matched layer(PML) absorbing boundaries are constrained to low-order time schemes. The auxiliary differential equation PML (ADE-PML) is a kind of unsplit PMLs that can be applied to arbitrary-order time schemes, and also can implement complex-frequency-shifted (CFS) coordinate tensor to enhance the performance of PML at high-angle incidence.We introduce the ADE-PML in the numerical modeling in fourth-order Runge-order time scheme for acoustic wave equations, and validate the efficiency. The numerical experiments demonstrate that CFS ADE-PML preforms better than non-CFS ADE-PML at high-angle incidence. Besides, there are minor distinctions among the absorbing efficiency of different update schemes for auxiliary variable, and the solution calculated by explicit scheme better agree with the analytical one. The long-time computation of energy decay shows that ADE-PML is stable up to 2 × 105 steps.

Key words: acoustic wave equation, perfectly matched layer, auxiliary differential equations, arbitrary-order time scheme, numerical modeling

中图分类号: 

  • P631.4
[1] 夏英杰, 陈华斌, 陈健, 姚明宇, 杨海. 不同施工参数下深层页岩水力压裂数值模拟[J]. 吉林大学学报(地球科学版), 2026, 56(3): 924-936.
[2] 涂君汕, 陈辉, 邓居智, 陈康, 李焱, 余辉. 水库堤坝对高密度电法的影响及校正技术[J]. 吉林大学学报(地球科学版), 2026, 56(3): 1013-1025.
[3] 束龙仓, 位书静, 澈丽木格, 温中琦, 刘波. 西辽河平原生态输水的地下水响应及不确定性量化[J]. 吉林大学学报(地球科学版), 2026, 56(2): 647-660.
[4] 黄兴国, 翁央央, 韩丽. 基于广义递归卷积的孔隙黏弹地震波正演模拟[J]. 吉林大学学报(地球科学版), 2026, 56(1): 377-385.
[5] 周柳湘, 余思琴, 陈俊华, 刘城, 陈义, 张鑫鑫, . 潜孔锤与冲击钻机联合钻进工艺成孔碎岩过程数值模拟[J]. 吉林大学学报(地球科学版), 2025, 55(5): 1608-1618.
[6] 苑成旺, 张敏, 张少龙, 秦磊, 郭海洋, 墨海滢, 信旸. 长春地铁基坑支护结构变形与地表沉降特征[J]. 吉林大学学报(地球科学版), 2025, 55(3): 879-892.
[7] 于子望, 卢帅屹, 白林, 郑天琪. CO2地质封存岩石力学问题研究进展[J]. 吉林大学学报(地球科学版), 2025, 55(3): 930-942.
[8] 程瑶, 陆丹丹, 赵龙飞.

微波加热油页岩储层的热响应 [J]. 吉林大学学报(地球科学版), 2025, 55(2): 387-400.

[9] 刁国君, 何昕.

地铁明挖车站下穿越既有建筑物的局部暗扩挖施工关键技术 [J]. 吉林大学学报(地球科学版), 2025, 55(2): 536-549.

[10] 蔺学旻, 肖红琳. 缝洞型储层声波远探测测井响应模拟及应用——以塔河油田碳酸盐岩缝洞型储层为例[J]. 吉林大学学报(地球科学版), 2025, 55(1): 312-327.
[11] 魏虹羽, 李世超, 王伟安. 地球动力学数值模拟算法的应用现状与展望[J]. 吉林大学学报(地球科学版), 2025, 55(1): 98-124.
[12] 吴泽坤, 何来胜, 白晓宇, 麻栋栋, 牛永昌, 赵广, 桑松魁, 闫楠, 张明义. 层状土中静压桩连续贯入现场试验与数值模拟[J]. 吉林大学学报(地球科学版), 2024, 54(4): 1291-1304.
[13] 庞志超, 肖华, 毛晨飞, 陈国军, 梁琬坤, 高明, 张啸. 准噶尔盆地南缘地区含膏质地层岩性特征及测井识别方法[J]. 吉林大学学报(地球科学版), 2024, 54(4): 1419-1431.
[14] 刘怀湜, 赵雪娇, 刘雨新, 房庭瑞, 张静, 嵇艳鞠. 基于改进粒子群优化的超顺磁效应多参数提取[J]. 吉林大学学报(地球科学版), 2024, 54(3): 993-1002.
[15] 韩长玉, 赵浩汀.

路基填筑短期内热棒高导热性影响 [J]. 吉林大学学报(地球科学版), 2024, 54(2): 570-580.

Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 尤敏鑫,刘建民. 同位素地球化学在峨眉山大火成岩省研究中的应用现状与进展[J]. 吉林大学学报(地球科学版), 2014, 44(4): 1231 -1243 .
[2] 杨春梅, 李洪奇,陆大卫,张方礼,高 原,邵英超. 不同驱替方式下岩石电阻率与饱和度的关系[J]. J4, 2005, 35(05): 667 -671 .
[3] 祝洪臣,张炯飞,权 恒. 大兴安岭中生代两期成岩成矿作用的元素、同位素特征及其形成环境[J]. J4, 2005, 35(04): 436 -0442 .
[4] 朱建伟, 赵刚, 刘博, 郭巍, 成俊. 油页岩测井识别技术及应用[J]. J4, 2012, 42(2): 289 -295 .
[5] 陈力,梁海安,张文娟,荣帆. 模糊数学方法在城市工程地质环境区划中的应用--以抚顺市城区为例[J]. J4, 2008, 38(5): 837 -0840 .
[6] 高桂梅,苏 克,王文颖,甘树才,刘招君. 吉林省桦甸油页岩中稀土元素和微量元素的研究[J]. J4, 2006, 36(6): 974 -0979 .
[7] 吴孔运,蒋忠诚,叶 晔. 不同植物群落对灰岩试块溶蚀速率的影响[J]. J4, 2007, 37(5): 967 -0971 .
[8] 周彦章,迟宝明,刘中培. 山东夏甸金矿床充水机理构造控制模式[J]. J4, 2008, 38(2): 255 -0260 .
[9] 张渊,刘连登,孙景贵,陈国华,张洪喜,闫复传,杨开春. 胶东西北部黄埠岭金矿床两期次叠加成矿[J]. J4, 2008, 38(1): 21 -0026 .
[10] 鲁程鹏, 束龙仓, 苑利波, 张蓉蓉, 黄币娟, 王彬彬. 基于示踪试验求解岩溶含水层水文地质参数[J]. J4, 2009, 39(4): 717 -721 .