胡勇， 潘冬明*， 吴如山， 韩立国， 段超然

1 中国矿业大学， 资源与地球科学学院， 徐州 221008 2 Modeling and Imaging Laboratory, University of California, Santa Cruz 95060, USA 3 吉林大学， 地球探测科学与技术学院， 长春 130026 4 常州工学院， 土木建筑工程学院， 常州 213000

0 引言

随着油气勘探技术的不断发展，同时对地下构造探测精度要求也逐渐提高，如今正从构造勘探阶段逐步转向岩性勘探阶段，其中地震数据偏移成像方法在油气勘探领域占有重要的地位.近年来Kirchhoff 偏移、高斯束偏移、单程波偏移和逆时偏移等地震成像方法得到了快速发展(秦宁等，2015；You et al.，2019；豆辉和徐逸鹤，2019；Sun et al., 2020).在复杂介质中，逆时偏移(RTM)方法具有成像精度高的优点，能够适用于横向速度变化剧烈的区域，获取地下构造高精度反射系数(Sun and Zhang, 2009；马方正等，2016；You et al.，2017).但是逆时偏移的伴随算子是正演算子的共轭转置，而并非其逆算子(Claerbout，1991)，因此无法获得真振幅成像结果.

为了进一步提高逆时偏移的成像精度，发展了最小二乘逆时偏移方法(LSRTM)，其在反演的理论框架下，利用反偏移数据来不断地与观测数据相匹配，并通过优化算法进行迭代，最终获取地下构造的高精度成像结果.关于LSRTM的发展，首先由Nemeth等(1999)在最小二乘目标函数的框架下实现了Kirchhoff偏移.随后，Dai等(2012)、Dai和Schuster(2013)利用多源地震数据进行测试，很大程度上提高了LSRTM的计算效率.任浩然等(2013)将Hessian算子应用于LSRTM中，获得了相对保真的成像结果，改善了地震偏移的成像精度.Tan和Huang(2014)在原始LSRTM的基础上，提出了波场分离成像条件，并通过更新震源波场实现陡倾角断层成像.此外，为了克服……