帅庆伟，徐云霞，文鹏飞，沙志彬，万晓明

1.自然资源部海底矿产资源重点实验室，广州 510075

2.中国地质调查局广州海洋地质调查局，广州 510760

天然气水合物是一种存在于低温、高压环境下的甲烷聚合物。海洋天然气水合物主要分布于水深300～4 000 m的深水盆地[1]。琼东南盆地发育于南海北部西北端大陆边缘，位于南海北部陆坡深水区，具有天然气水合物形成需要的高压低温稳定域环境。盆地内具有生物成因和热成因天然气的巨大气源供给条件，是重要的天然气水合物远景区。盆地内具有非常丰富的泥底辟体、气烟囱及伴生的底辟断层、裂隙等理想的气体运移通道，导致地震反射层发生畸变，呈现出无反射、模糊反射、杂乱反射等特征。而琼东南丰富的含气特征，导致地震波能量被大量吸收而发生严重衰减，会在成像剖面上形成大量模糊、空白反射或同相轴下拉现象[2-6]。为了解决含气导致的成像模糊、时间偏移归位不准等问题，恢复真实的地层构造信息，本文对琼东南水合物资料进行叠前深度偏移处理。

速度层析既是深度偏移的过程也是深度偏移的结果，在深度偏移中占有非常重要的作用。目前常用的速度更新主要采用线性更新，在每次速度迭代中均需要进行一次深度偏移，重复性的偏移和拾取剩余曲率（RMO）信息会严重降低项目运行效率，导致计算成本较高。另一种速度更新方法是非线性反演方法，在速度更新过程中不需要进行多次迭代偏移，能极大地提高速度更新效率，V Farra于1988年提出基于拾取的剩余曲率进行反偏移获得运动学不变量进行非线性反演[7]，能有效提升速度反演效率；Clapp等将构造信息加入到反演中，通过引入模型中间变量，使得反演结果具有更精细的构造特征[8]；Emma Evans等提出在非线性层析过程中加入层位约束，可以通过硬约束和软约束来提升断层等构造的速度更新精度[9]；Guillaum等提出在层析反演中加入倾角约束[10]，能极大提高速度反演精度。

考虑提升水合物成像精度、恢复真实构造形态的同时考虑提升速度反演运算效率，本文尝试利用非线性速度层析方法对琼东南水合物数据进行速度建模，并进行深度偏移处理。成像结果有效地消除了含气导致的同相轴下拉现象，断层等构造特征成像清晰，有效地获得了真实的构造成像结果，取得了良好的应用效果。

1 非线性层析

线性层析每更新一次速度都要进行一次深度偏移，要获得最终理想的速度模型需要进行多次深度偏移，使得速度更新的效率相对较低。非线性层析在速度更新过程中仅需进行一次深度偏移，在深度偏移道集拾取剩余曲率，通过对拾取的剩余量进行反偏移从而获得各反射点的双程走时和射线参数，然后利用这些信息在叠前域获得运动学不变量。将运动学不变量与速度模型进行多次非线性反演，进行模型更新，直到获得理想的速度[11-13]。非线性层析流程如图1所示。

本文非线性层析实现过程主要分为4步：①剩余曲率拾取，②运动学反偏移，③混合模型建立，④非线性反演。

1.1 剩余曲率拾取

在共偏移距域进行剩余曲率（RMO）拾取，其描述从零偏移距到最大偏移距深度的误差，拾取的精度决定了后续速度层析的精度[14-16]。为保证提高拾取的精度，降低由噪声引起的拾取误差，首先要对共偏移距道集进行噪声压制，提高资料的信噪比。同时需要对拾取后的RMO进行质量控制，分析拾取的剩余曲率是否与道集同相轴一致，并且要将远偏移距处畸变的拾取剔除掉，如图2a所示。另外，在拾取RMO的过程中，也会同时生成剩余曲率（gamma）属性体与质量因子（semblance）属性体，可对拾取结果的质量进行控制。gamma属性体描述整体剩余曲率情况，蓝色表示进行深度偏移时速度偏小，红色表示速度偏高；semblance为构造约束属性，值越高说明拾取的可信度越高。通过两属性体的分析，可以直观地分析出哪些部分速度偏小，哪些地方速度偏高，哪些地方拾取更为可靠，如图2b与图2c所示。

图1 非线性层析流程Fig.1 Flow chart of non-linear tomography technology

1.2 运动学反偏移

运动学反偏移即将拾取的RMO深度误差以及gamma与semblance等属性参数进行反偏移，涉及到的概念为图偏移，图偏移能够将时间域的层位转换到深度域，在时间域里拾取的层位能够提供2种信息：例如在一个偏移剖面中，有一法向入射的射线如图3a所示，通过射线能够计算双程旅行时；与法向入射有一定夹角的射线如图3b所示，与法向入射相比有一定的时间坡度。通过这个拾取，就能找到沿产生这种射线的入射点XNR，也因此能够通过图偏移获得实际的速度模型。在模型更新过程中，能够消除速度模型对原来的深度域观测模型的依赖，从而提升模型更新的效率。

图2 剩余量拾取a.道集中 RMO 拾取量，b.剩余曲率属性，c.质量因子属性。Fig.2 The remaining quantity picking up a.RMO quantity is picked up from gathers，b.gamma attribute，c.semblance attribute.

图3 运动学反偏移原理Fig.3 Principle of kinematic reverse migration

1.3 混合模型建立

目前主流速度模型为网格模型，网格单元的大小控制了每一层内速度的变化，以及各地层之间的接触关系。一个精细的网格，能够精确地定义各层之间的位置，但需要非常大的数据空间。而混合模型包含了地质单元与许多地震属性，每一层都被定义为一个独立的属性体。由于属性体中网格单元的大小与层位位置无关，单元网格的大小可以定义得大一些，因此，混合模型比传统的网格模型更加节省资源；较之用网格来定义每一层的属性，混合模型采用样条插值思想，在保证模型精度的前提下，能够大大缩减模型大小[9,17-18]，如图4所示。假设一维速度模型如图4a所示，均匀分割的网格模型如要真实地还原速度模型，则需要较小的网格（图4b），较为稀疏的网格则不能真实接近速度模型如图4c所示，但采用样条插值思想，却能用较为稀疏的网格点真实地还原速度模型（图4d）。

所包含的属性信息包括骨架信息与倾角信息等，倾角信息是在初始深度偏移叠加剖面上做最小二乘搜索离散的局部地层斜率数据体，产生两个三维离散地层斜率数据体，分别为投影到主测线和联络测线方向的地层斜率[19-21]。图5为主测线方向的地层倾角剖面显示，蓝色表示地层倾向沿北东方向，红色表示沿北西向，白色表示地层较平，由图可见，计算的倾角与实际构造吻合较好，能够为非线性反演提供良好的约束。

图4 样条插值思想a.一维速度模型，b.小网格还原速度情况，c.大网格还原速度情况，d.样条插值还原速度情况。Fig.4 Spline interpolation a.1D velocity model，b.small grid reduction velocity model，c.big grid reduction velocity model，d.spline interpolation reduction velocity model.

图5 主测线方向的地层斜率剖面Fig.5 Stratigraphic inclining profile in the direction of inline

1.4 非线性反演

非线性层析，采用数据驱动的全三维层析反演更新速度模型，每一次模型更新均是一个正、反演迭代的过程，主要包括两个部分：正演模拟，以产生速度模型的RMO曲线；反演模拟，即扰动和建模RMO之间的线性关系，从而产生模型更新[22-23]。正演和反演是反复迭代进行的，直到最后的RMO量足够低。

图6a与图6b分别为初始速度模型和最终速度模型所计算的RMO量统计直方图，横坐标为gamma值，纵坐标为gamma值对应的拾取个数，gamma为1说明偏移道集很平；gamma小于1则速度偏小，道集上翘；gamma大于1则速度过大，道集下拉。最终速度模型相对初始速度模型而言，gamma值更趋近于1，说明速度更加准确。

2 应用效果分析

在前期道集精细处理的基础上，结合琼东南海域天然气水合物相对的地球物理特征以及地质构造特征，应用非线性层析方法对该区速度模型进行多次非线性速度层析，提升速度模型精度和深度偏移成果质量。整个层析过程中仅进行一次深度偏移，大大提升了速度层析的速度。

图7是水合物目标区的一条主测线的叠前时间偏移剖面，BSR的强振幅特征刻画清晰。在水合物稳定域之下，由于含有丰富的气体，导致地震反射同相轴出现下拉现象，不利于进行BSR之下的地质构造解释，不利于建立精细的水合物相关地质模型。图8为利用非线性层析技术获得速度模型后，获得的深度偏移成果，水合物分布位置约在1 900 m，水深约为 1 700 m，可知该区域水合物分布较浅，为海底以下200 m左右；水合物之下的含气区域，由于速度模型准确，经过深度偏移，由含气导致的同相轴下拉现象消除，恢复了真实的地层特征；且基底的陡倾角特征更加真实可靠。

图9为非线性层析前后速度对比图，图9a为cdp29350处的速度曲线对比图，初始深度域层速度模型为线性趋势，不能较好地区分速度增大和减小区域；更新后的曲线能很明显地分析出速度的变化趋势，更利于进行含水合物和含气区域的分析；图9b与图9c分别为进行深度偏移的初始模型与经过非线性层析后的最终模型在深度2 200 m处的切片显示。切片显示中速度值域范围为 1 508～2 000 m/s，蓝色为低速度值，红色为高速度值。速度层析以后的速度切片中有更多的细节特征，对由含气导致的低速异常范围刻画更加清晰，低速异常值范围约为1 650～1 750 m/s，通过更新后的速度，能更加直观地对水合物及游离气的空间分布范围进行刻画。

图6 RMO 量统计直方图a.初始速度模型计算的RMO量统计直方图，b.最终速度模型计算的RMO量统计直方图。Fig.6 RMO quantity statistical histogram a.statistical histogram of RMO quantity calculated by initial velocity model，b.statistical histogram of RMO quantity calculated by final velocity model.

图7 主测线叠前时间偏移剖面Fig.7 Prestack time migration profile along inline

图8 主测线叠前深度偏移剖面Fig.8 Prestack depth migration profile along inline

图9 非线性层析前后速度切片对比a.层析前后速度曲线对比，b.初始速度切片，c.非线性层析后速度切片。Fig.9 Comparison of velocity slices before and after nonlinear tomography a.the comparison of velocity curves before and after tomography，b.initial velocity slice，c.Velocity slice after nonlinear tomography.

3 结束语

非线性层析充分运用运动学不变量信息，能有效避免速度层析过程中的多次迭代深度偏移，能够极大提高速度层析的效率；且非线性层析过程中采用混合模型方式，有效地融入倾角、地层骨架等地震属性信息，采用样条插值的模式进行模型建立，能有效地减小模型空间，提升模型层析精度。

在实际天然气水合物资料应用中，非线性层析获得的速度模型能有效地对速度异常体进行精细刻画；基于此进行的深度偏移，其成果能有效地消除含气导致的地层同相轴下拉现象，恢复含气区域的地层真实形态；能有效地对倾角区域进行清晰成像，使得地下地质构造特征更加真实可靠。