最小二乘Kirchhoff射线束叠前时间偏移

刘俊成， 陈祥忠 ， 郭井学， 王瑞星， 惠梦琳， 黄申硕

1 吉林大学地球探测科学与技术学院， 长春 130026 2 北京桔灯地球物理勘探股份有限公司， 北京 102200 3 中国极地研究中心， 上海 200129

0 引言

地震偏移成像是地震数据处理过程中的重要环节(Schneider, 1978; Stolt et al., 1987; Baysal et al., 1983; Claerbout, 1992; Gray et al., 2000; 王华忠等，2009；孙小东等，2020)，其不但可以用于恢复地下的构造信息，还可以估算地下的岩性和物性参数用于后续的地震反演和储层预测.按照偏移速度的不同定义方式，地震偏移方法可以分为两大类：时间偏移和深度偏移.深度偏移可以适应地下介质速度的剧烈横向变化，是复杂介质成像的首选技术手段，但其实际应用效果严重依赖于速度建模的精度.时间偏移虽然无法对复杂构造进行准确的成像，但当速度横向变化相对平缓时，其依然可以满足大部分工区的成像处理需要.此外，相对于深度偏移，时间偏移的计算效率更高，所需偏移速度也更容易估算.因此，即便在计算机硬件高度发达的今天，时间偏移依然是业界广泛使用的地震成像方法.

众所周知，常规的偏移方法只是地震波场正演算子的共轭转置(Tarantola, 1984; Claerbout, 1985)，对于有限观测系统所采集的带限地震数据，其只能生成模糊的地下构造成像结果.此外，复杂的地下介质构造和非规则的地震数据空间采样，还会导致偏移假象和非规则的成像照明，严重影响成像振幅信息的可靠性.基于线性化地震反演理论的最小二乘偏移(Least-Squares Migration, LSM)是解决上述问题的有效手段(Nemeth et al., 1999; Duquet et al., 2000; Dai and Schuster, 2013; Hu et al., 2016; Yue et al., 2019；周东红等，2020；刘国章等，2020；……