基于频率域高阶有限差分法的正演模拟及并行算法

基于频率域高阶有限差分法的正演模拟及并行算法

殷文^{1, 2}

1.东营职业学院，山东东营 257091;2.中石化胜利油田有限公司物探研究院，山东东营 257022

收稿日期:2007-02-01 修回日期:1900-01-01 出版日期:2008-01-26 发布日期:2008-01-26
通讯作者: 殷文

Forward Modeling and Parallel Algorithm Based on High-Order Finite-Difference Method in Frequency Domain

YIN Wen^{1, 2}

1.Vocation College of Dongying City, Dongying, Shandong 257091, China;2.Geophysical Research Institute, Shengli Oilfield Company Limited SINOPEC, Dongying, Shandong 257022, China

Received:2007-02-01 Revised:1900-01-01 Online:2008-01-26 Published:2008-01-26
Contact: YIN Wen

摘要/Abstract

摘要： 在弹性波频率空间域有限差分数值模拟方面，差分网格及边界条件是影响弹性波模拟成功与否的关键，为了压制数值模拟中的网格频散，采用25点有限差分算子，建立了有限差分矩阵方程，且借鉴匹配层衰减边界条件思想，设计了弹性波频率空间域有限差分数值模拟算法。由于采用高阶有限差分法来提高差分格式的精度，将会导致计算量显著增加，为此，对频率空间域有限差分弹性波数值模拟方法，采用流水线技术与分治策略进行了并行算法研究，提高了计算效率，使得在合理的计算时间内更精确地模拟弹性波在弹性介质中的传播过程。

关键词: 正演模拟, 有限差分算子, 完全匹配层, 并行算法

Abstract: For finite difference numerical simulation in elastic-wave frequencyspace domain, difference grids and boundary condition are the key for elasticwave modeling. 25-point finite-different operators were adopted to suppress grid dispersion in numerical simulation and finite-difference matrix equations were constructed. Making use of the idea of match-layer attenuation boundary condition, the numerical simulation algorithm of elastic-wave finite difference in frequency-space domain was designed. Besides, the computing loads were increased observably by using high-order finite-difference method to improve accuracy. So for elastic-wave numerical simulation of finite difference in frequency-space domain, parallel algorithm was also studied in combination with divided-conquer and pipelining technique to enhance computational efficiency, so that elastic wave propagation in elastic medium was simulated much more accurately in reasonable computation time.

Key words: forward modeling, finite-difference operators, perfectly matched layer, parallel algorithm

中图分类号:

P631.4

殷文. 基于频率域高阶有限差分法的正演模拟及并行算法[J]. J4, 2008, 38(1): 144-0151.

YIN Wen. Forward Modeling and Parallel Algorithm Based on High-Order Finite-Difference Method in Frequency Domain[J]. J4, 2008, 38(1): 144-0151.

参考文献

相关文章 14

[1]	王泰涵, 黄大年, 马国庆, 李野, 林松. 基于并行预处理算法的三维重力快速反演[J]. 吉林大学学报(地球科学版), 2018, 48(2): 384-393.
[2]	安振芳, 张进, 张建中. 海洋三维VC观测系统优化设计[J]. 吉林大学学报(地球科学版), 2018, 48(1): 271-284.
[3]	刘志强, 孙建国, 孙辉, 刘明忱, 高正辉, 石秀林. 曲线坐标系下的完全匹配层吸收边界条件[J]. 吉林大学学报(地球科学版), 2017, 47(6): 1875-1884.
[4]	刘海飞, 柳杰, 高寒, 郭荣文, 童孝忠, 麻昌英. 五极纵轴激电测深三维有限元正演模拟[J]. 吉林大学学报(地球科学版), 2016, 46(3): 884-892.
[5]	李庆洋, 李振春, 黄建平, 李娜, 苏在荣. 基于贴体全交错网格的起伏地表正演模拟影响因素[J]. 吉林大学学报(地球科学版), 2016, 46(3): 920-929.
[6]	朱超, 夏志远, 王传武, 宋光永, 魏学斌, 王鹏, 王海峰, 王波. 致密油储层甜点地震预测[J]. 吉林大学学报(地球科学版), 2015, 45(2): 602-610.
[7]	赵建国，史瑞其，陈竞一，赵维俊，王宏斌，潘建国. 辅助微分方程完全匹配层在声波方程数值模拟中的应用[J]. 吉林大学学报(地球科学版), 2014, 44(2): 675-682.
[8]	汤文武，柳建新，童孝忠. 电导率连续变化的线源FCSEM有限元正演模拟[J]. 吉林大学学报(地球科学版), 2013, 43(5): 1646-1654.
[9]	刘斌, 李术才, 李树忱, 聂利超, 钟世航, 宋杰, 刘征宇. 隧道含水构造电阻率法超前探测正演模拟与应用[J]. J4, 2012, 42(1): 246-253.
[10]	韩利, 韩立国, 李翔, 王德利, 崔杰. 二阶声波方程频域PML边界条件及频域变网格步长并行计算[J]. J4, 2011, 41(4): 1226-1232.
[11]	赵迎月, 顾汉明, 李宗杰, 杨子川, 蔡成国. 塔中地区奥陶系典型地质体地震识别模式[J]. J4, 2010, 40(6): 1262-1270.
[12]	董航，刘四新，王春晖，曾昭发，鹿琪，王者江，易兵. 探地雷达测量近地表含水量的研究[J]. J4, 2009, 39(1): 163-0167.
[13]	翁爱华,王雪秋. 高精度地面核磁共振数据正演计算[J]. J4, 2007, 37(3): 620-0623.
[14]	刘四新，曾昭发，徐波. 地质雷达数值模拟中有损耗介质吸收边界条件的实现[J]. J4, 2005, 35(03): 378-0381.

Metrics

Viewed

Full text

Abstract

Cited

Shared

Discussed