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

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

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

赵建国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]. 吉林大学学报(地球科学版), 2022, 52(1): 162-.
[2] 韩复兴, 王若雯, 孙章庆, 高正辉. 地震声波数值模拟中人工边界条件的差别与组合[J]. 吉林大学学报(地球科学版), 2022, 52(1): 261-.
[3] 师文豪, 杨天鸿. 渗流应力耦合作用下顺倾向层状边坡各向异性渗流特征数值模拟[J]. 吉林大学学报(地球科学版), 2021, 51(6): 1783-1788.
[4] 余莉, 张钰, 王维玉, 韩子豪, 赵拓. 基坑装配式可回收支护和桩锚支护结构的受力与变形分析[J]. 吉林大学学报(地球科学版), 2021, 51(6): 1789-1800.
[5] 魏家斌, 王卫东, 吴江斌. 免共振沉桩过程对地表振动影响的FLAC3D数值模拟[J]. 吉林大学学报(地球科学版), 2021, 51(5): 1514-1522.
[6] 李立云, 王子英, 王晓静, 杜修力. 近铁路基坑通风井段变形特征及其机制分析[J]. 吉林大学学报(地球科学版), 2021, 51(5): 1441-1451.
[7] 蔡晓光, 徐洪路, 李思汉, 张少秋. 地震作用下返包式加筋土挡墙数值模拟[J]. 吉林大学学报(地球科学版), 2021, 51(5): 1416-1426.
[8] 李一赫, 王殿举, 于法浩, 刘志强. 下刚果盆地白垩系盐构造的形成演化[J]. 吉林大学学报(地球科学版), 2020, 50(6): 1628-1638.
[9] 吕雅馨, 骆祖江, 徐成华. 南京汤山地区地热水资源评价[J]. 吉林大学学报(地球科学版), 2020, 50(6): 1844-1853.
[10] 盛冲, 许鹤华, 张云帆, 张文涛, 任自强. 钙质砂水理性质及对岛礁淡水透镜体形成的影响[J]. 吉林大学学报(地球科学版), 2020, 50(4): 1127-1138.
[11] 段云星, 杨浩. 增强型地热系统采热性能影响因素分析[J]. 吉林大学学报(地球科学版), 2020, 50(4): 1161-1172.
[12] 孙可明, 张宇. 缝网间距对高温岩体储留层温度影响规律模拟[J]. 吉林大学学报(地球科学版), 2019, 49(6): 1723-1731.
[13] 孙超, 许成杰. 基坑开挖对周边环境的影响[J]. 吉林大学学报(地球科学版), 2019, 49(6): 1698-1705.
[14] 王常明, 李桐, 田书文, 李硕. 基于LAHARZ的泥石流堆积范围预测模型的建立及应用[J]. 吉林大学学报(地球科学版), 2019, 49(6): 1672-1679.
[15] 杨新乐, 秘旭晴, 张永利, 李惟慷, 戴文智, 王亚鹏, 苏畅. 注热联合井群开采煤层气运移采出规律数值模拟[J]. 吉林大学学报(地球科学版), 2019, 49(4): 1100-1108.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 张渊,刘连登,孙景贵,陈国华,张洪喜,闫复传,杨开春. 胶东西北部黄埠岭金矿床两期次叠加成矿[J]. J4, 2008, 38(1): 21 -0026 .
[2] 鲁程鹏, 束龙仓, 苑利波, 张蓉蓉, 黄币娟, 王彬彬. 基于示踪试验求解岩溶含水层水文地质参数[J]. J4, 2009, 39(4): 717 -721 .
[3] 路兴昌,张艳红. 基于三维激光扫描的空间地物建模[J]. J4, 2008, 38(1): 167 -0171 .
[4] 崔建军, 刘晓春, 胡娟, 曲玮. 桐柏杂岩中印支期变质岩包体的变质作用[J]. J4, 2009, 39(4): 618 -629 .
[5] 尤敏鑫,刘建民. 同位素地球化学在峨眉山大火成岩省研究中的应用现状与进展[J]. 吉林大学学报(地球科学版), 2014, 44(4): 1231 -1243 .
[6] 杨春梅, 李洪奇,陆大卫,张方礼,高 原,邵英超. 不同驱替方式下岩石电阻率与饱和度的关系[J]. J4, 2005, 35(05): 667 -671 .
[7] 张彦龙,葛文春,柳小明,张吉衡. 大兴安岭新林镇岩体的同位素特征及其地质意义[J]. J4, 2008, 38(2): 177 -0186 .
[8] 张玉玲,姚军,赵晓波,曹春英,郑松志. 复合型微生物絮凝剂产生菌YL3的优化条件[J]. J4, 2008, 38(5): 864 -0868 .
[9] 于 强,任战利. 鄂尔多斯盆地黄陵、东胜地区地温场对比[J]. J4, 2008, 38(6): 933 -0936 .
[10] 陈力,梁海安,张文娟,荣帆. 模糊数学方法在城市工程地质环境区划中的应用--以抚顺市城区为例[J]. J4, 2008, 38(5): 837 -0840 .