单程波动方程叠前深度偏移并行算法
|
摘要
根据单程波动方程波场外推原理,文中总结了共炮道集数据混合域波场外推与成像算法流程,该流程主要包括地震数据预处理、波动方程速度分析和波场外推成像3部分;鉴于单程波动方程波场外推算法和数据流具有非常适合粗粒度的任务级并行计算,处理器间的通讯较少,各个计算节点上可能有多个处理任务,且各节点的处理任务不等,负载极不平衡等特点,本文采用了主从模式动态负载平衡并行算法,该算法可以充分利用多个处理节点的空闲资源,提高作业的计算效率;将混合域单程波动方程波场外推算法与主从模式动态负载平衡并行算法相结合,实现了软件级断点恢复功能。应用国产"Stseis叠前深度成像处理软件系统"对胜利探区的实际资料进行了生产性应用处理,取得了较好的叠前成像结果。
According to wavefield extrapolation principle of one-way wave equation,the paper summarized the flow of wavefield extrapolation and imaging algorithm of common shot gather in mixed domain,which mainly included three parts: seismic data preprocessing,wave equation velocity analysis and wavefield extrapolation imaging.In view of the characters of wavefield extrapolation algorithm of one-way wave equation and data flow: very suitable for parallel computation of gross-grain task grade,less communication between processor,multiple processing tasks,unequal processing task and extremely unbalanced load at each computing node,the paper adopted master/slaver dynamic load balanced parallel algorithm that can fully use idle resources of multiple processing nodes,improving computing efficiency of job;combining the wavefield extrapolation algorithm of one-way wave equation in mixed domain with master/slaver dynamic load balanced parallel algorithm realized the restoration function of break point of software grade.Using Chinese-built "Stseis prestack depth imaging processing software system" to carry out productive application processing for practical data in Shengli exploration area achieved better prestack imaging results.
According to wavefield extrapolation principle of one-way wave equation,the paper summarized the flow of wavefield extrapolation and imaging algorithm of common shot gather in mixed domain,which mainly included three parts: seismic data preprocessing,wave equation velocity analysis and wavefield extrapolation imaging.In view of the characters of wavefield extrapolation algorithm of one-way wave equation and data flow: very suitable for parallel computation of gross-grain task grade,less communication between processor,multiple processing tasks,unequal processing task and extremely unbalanced load at each computing node,the paper adopted master/slaver dynamic load balanced parallel algorithm that can fully use idle resources of multiple processing nodes,improving computing efficiency of job;combining the wavefield extrapolation algorithm of one-way wave equation in mixed domain with master/slaver dynamic load balanced parallel algorithm realized the restoration function of break point of software grade.Using Chinese-built "Stseis prestack depth imaging processing software system" to carry out productive application processing for practical data in Shengli exploration area achieved better prestack imaging results.
引文
[1]Bevc D.Imaging complex structures with semirecur-sive Kirchhoff migration.Geophysics,1997,62(2):577~588
[2]Biondi B and Fomel S.Azimuth moveout for 3-D pres-tack imaging.Geophysics,1998,63(2):574~588
[3]Biondi B and Palacharla G.3-D prestack migrationof common-azimuth data.Geophysics,1996,61(6):1822~1832
[4]徐兆涛等.三维波动方程叠前深度偏移在资料处理中的应用.石油物探,2003,42(增刊):31~34
[5]程玖兵等.频率—空间域有限差分法叠前深度偏移.地球物理学报,2001,44(3):389~395
[6]马在田.地震成像技术——有限差分法偏移.北京:石油工业出版社,1989
[7]辛可锋等.叠前地震数据的规则化.石油地球物理勘探,2002,37(4):311~317
[8]马在田等.计算地球物理学概论.上海:同济大学出版社,1997
[9]李振春等.波动方程法共成像点道集偏移速度建模.地球物理学报,2003,46(1):86~93
[10]杨淑卿等.基于P道集的波动方程偏移速度分析技术及其应用.石油物探,2003,42(增刊):31~34
[11]Fomel S.Time-migration velocity analysis by velocitycontinuation.Geophysics,2003,68(5):1662~1672
[12]Zhenyue Liu and Bleistein N.Migration velocity anal-ysis:Theory and an iterative algorithm.Geophysics,1995,60(1):142~153
[13]马在田.高阶方程偏移的分裂算法.地球物理学报,1983,26(4):377~388
[14]王华忠.三维地震波场最优外推算子及其在深度偏移成像中的应用(同济大学博士论文),1997
[15]张关泉.波场分裂、平方根算子与偏移.反射地震学论文集.上海:同济大学出版社,2000,6~20
[16]李幼铭等.地震反问题评述.地球物理学进展,1992,7(3):23~31
[17]Claerbout J F.Imaging the earth’s interior.Black-well Scientific Publications,1985
[18]Huang L J et al.Extended local Born Fourier migra-tion.Geophysics,1999a,64(5):1524~1534
[19]Jin S et al.3-D wave equation GSP prestack depth mi-gration with application to a deep-water dataset fromGulf of Mexico.Expanded Abstracts of72nd SEGMtg,2002,1316~1319
[20]陈国良.并行计算——结构、算法、编程.北京:高等教育出版社,1999
[21]都志辉.高性能计算并行编程技术——MPI并行程序设计.北京:清华大学出版社,2001
[22]徐兆涛等.复杂地质体高精度成像方法研究与应用.胜利油气地球物理技术论文集.北京:石油工业出版社,2004,3~15
[23]Gazdag J et al.Migration of seismic data by phaseshift plus interpolation.Geophysics,1984,49(2):124~131
[24]徐兆涛.波动方程叠前深度成像理论与应用研究(同济大学博士论文),2006
[2]Biondi B and Fomel S.Azimuth moveout for 3-D pres-tack imaging.Geophysics,1998,63(2):574~588
[3]Biondi B and Palacharla G.3-D prestack migrationof common-azimuth data.Geophysics,1996,61(6):1822~1832
[4]徐兆涛等.三维波动方程叠前深度偏移在资料处理中的应用.石油物探,2003,42(增刊):31~34
[5]程玖兵等.频率—空间域有限差分法叠前深度偏移.地球物理学报,2001,44(3):389~395
[6]马在田.地震成像技术——有限差分法偏移.北京:石油工业出版社,1989
[7]辛可锋等.叠前地震数据的规则化.石油地球物理勘探,2002,37(4):311~317
[8]马在田等.计算地球物理学概论.上海:同济大学出版社,1997
[9]李振春等.波动方程法共成像点道集偏移速度建模.地球物理学报,2003,46(1):86~93
[10]杨淑卿等.基于P道集的波动方程偏移速度分析技术及其应用.石油物探,2003,42(增刊):31~34
[11]Fomel S.Time-migration velocity analysis by velocitycontinuation.Geophysics,2003,68(5):1662~1672
[12]Zhenyue Liu and Bleistein N.Migration velocity anal-ysis:Theory and an iterative algorithm.Geophysics,1995,60(1):142~153
[13]马在田.高阶方程偏移的分裂算法.地球物理学报,1983,26(4):377~388
[14]王华忠.三维地震波场最优外推算子及其在深度偏移成像中的应用(同济大学博士论文),1997
[15]张关泉.波场分裂、平方根算子与偏移.反射地震学论文集.上海:同济大学出版社,2000,6~20
[16]李幼铭等.地震反问题评述.地球物理学进展,1992,7(3):23~31
[17]Claerbout J F.Imaging the earth’s interior.Black-well Scientific Publications,1985
[18]Huang L J et al.Extended local Born Fourier migra-tion.Geophysics,1999a,64(5):1524~1534
[19]Jin S et al.3-D wave equation GSP prestack depth mi-gration with application to a deep-water dataset fromGulf of Mexico.Expanded Abstracts of72nd SEGMtg,2002,1316~1319
[20]陈国良.并行计算——结构、算法、编程.北京:高等教育出版社,1999
[21]都志辉.高性能计算并行编程技术——MPI并行程序设计.北京:清华大学出版社,2001
[22]徐兆涛等.复杂地质体高精度成像方法研究与应用.胜利油气地球物理技术论文集.北京:石油工业出版社,2004,3~15
[23]Gazdag J et al.Migration of seismic data by phaseshift plus interpolation.Geophysics,1984,49(2):124~131
[24]徐兆涛.波动方程叠前深度成像理论与应用研究(同济大学博士论文),2006
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心