刘立民,刘定进,李 博

(中国石油化工股份有限公司石油物探技术研究院,江苏南京211103)

目前,国内外常用的逆时偏移成像算法包括有限差分法和有限元法。有限差分法因实现过程简单且计算高效而被广泛应用,但是,当山地地表起伏较大或地质构造复杂时,难以确保成像精度,且有限差分逆时偏移的适应性和灵活性较差。有限元法具有良好的边界适应性及模拟精度,可以适应山地起伏地表等复杂边界问题,但面对大规模复杂三维问题时,难以在波场模拟精度和运算效率方面取得平衡。国内外学者从地下介质特征、波场传播算子、吸收边界条件以及并行计算效率等方面对逆时偏移方法进行了大量研究,逐步提高成像精度及运行效率。冀国强等[1]研究基于粘声介质的逆时偏移方法,同时考虑因地层吸收而导致振幅衰减与相位频散,提出一种正则化形式的稳定传播粘声逆时偏移方法,提升深层构造的分辨率与可信度。曲英铭等[2]提出面向高陡构造的粘声棱柱波逆时偏移方法,改善了高陡构造的照明和成像效果。巩向博等[3]研究基于稀疏约束的最小二乘逆时偏移方法,提高小尺度散射体的成像分辨率。李青阳等[4]提出一种新的卷积完全匹配层(NCPML)吸收边界条件,用于逆时偏移后,计算效率和内存占用上较常规的完全匹配层(SPML)吸收边界条件更优。段心标[5]提出基于GPU的傅里叶有限差分逆时偏移方法,保护了低频信息的成像且提高了成像效率。……