技术领域
[0001] 本
发明属于
地震偏移成像领域，具体涉及一种基于直达波场的最小二乘逆时偏移梯度预条件方法。
背景技术
[0002] 精确的地下构造偏移成像可为后续的地震资料解释提供精确的构造信息，一直以来都是
地震勘探领域的研究热点。最小二乘逆时偏移作为一种高
精度成像方法，将常规逆时偏移成像问题视为最小二乘
框架下的反演问题，采用反演的方法来求取地下的反射系数模型，能够获得高精度、高保幅、高
分辨率的成像结果。近年来，最小二乘逆时偏移技术的研究与发展受到了业界的广泛关注。然而，在常规最小二乘逆时偏移中，由于存在几何扩散效应以及地下照明不均衡的问题，浅层的
地震波能量要强于中深层，使得
迭代过程中反射系数模型的浅层更新一直占据主导地位，导致了中深层反演收敛较慢且精度较低。因此，采用合适的梯度预条件方法以较好地均衡浅、中、深层能量就显得尤为重要。
[0003] 当前梯度预条件方法主要有深度梯度加权法、Hessian矩阵类方法及地震波能量加权法三大类。深度梯度加权法采用深度加权公式来弥补深层能量弱的
缺陷，即对浅层梯度乘以较小的加权系数，而对深层梯度乘以较大的加权系数。但由于深度加权公式中的加权系数选取缺少严格的理论
支撑，此种方法只能粗略地估计梯度与深度的关系，中深层能量补偿并不准确。Hessian矩阵类方法利用Hessian矩阵或近似Hessian矩阵来修正振幅能量。此类方法通常需要显式计算Hessian矩阵或近似Hessian矩阵，其会不可避免地带来巨大的计算及内存消耗。这使得该类方法难以广泛应用于实际地震
数据处理中。为避免Hessian矩阵或近似Hessian矩阵的计算和存储，Tan和Huang于2014年在地震勘探最小二乘逆时偏移领域中采用了一种基于地震波能量的梯度预条件方法(即为地震波能量加权法)，其通过计算地下各点的正时、逆时波场能量并将其作为几何扩散效应的校正因子，然后对梯度进行加权处理，近似消除了几何扩散与地下照明不均衡的影响，大幅度提高了中深层的成像效果，其实现流程如图1所示。然而，常规的地震波能量中除了所需的直达波场能量以外，还包含有反射波、多次波等其它波场能量，这使得几何扩散效应的校正因子仍然不够精确，进而使最终的成像结果并未达到最佳。
发明内容
[0004] 本发明要解决的技术问题在于提供一种基于直达波场的最小二乘逆时偏移梯度预条件方法。本发明在最小二乘逆时偏移的每次迭代过程中，首先基于背景速度模型与反射系数模型，采用
声波方程有限差分线性化正演模拟技术得到地震记录，然后采用无反射声波方程有限差分正演模拟技术，模拟得到以炮点
震源子波为扰动的正时“近似直达”波场能量，和以上述线性化正演模拟得到的地震记录为扰动的逆时“近似直达”波场能量，最后应用此正时“近似直达”波场能量和逆时“近似直达”波场能量对梯度进行预条件处理。本方法可有效提高梯度预条件的精度，显著提升最小二乘逆时偏移的中深层成像效果，可为后续地震勘探的解释工作奠定良好的
基础。
[0005] 本发明采取以下技术方案：
[0006] 一种基于直达波场的最小二乘逆时偏移梯度预条件方法，其特征在于所述方法具体包括以下步骤：
[0007] (1)基于
水平层状介质速度模型，采用全排列的观测方式进行地震
数据采集，炮点与检波点均在地表，然后基于速度模型与采集的地震数据进行最小二乘逆时偏移成像；
[0008] (2)在最小二乘逆时偏移的第i次迭代计算中，i>＝1，输入上次迭代更新后的反射系数模型Mi–1，其中第1次迭代时输入给定的反射系数模型M0，基于背景速度模型v0、反射系数模型Mi–1，采用声波方程有限差分线性化正演模拟技术，以震源子波w为扰动，模拟得到地震记录R；
[0009] (3)基于背景速度模型v0，采用无反射声波方程有限差分正演模拟技术，以震源子波w为扰动，模拟得到正时“近似直达”波场Us；无反射声波方程的具体表达式为：
[0010]
[0011] 式(1)中，P为应
力，x、z分别为空间坐标，v为声
波速度，t为时间，f为震源扰动项；
[0012] (4)计算正时“近似直达”波场能量Es，其具体表达式如下：
[0013]
[0014] 式(2)中，Es为正时“近似直达”波场能量，Us为正时“近似直达”波场，t为时间，T为正演模拟时长；
[0015] (5)基于背景速度模型v0，采用无反射声波方程有限差分正演模拟技术，无反射声波方程的具体表达式如式(1)所示，然后以第(2)步中得到的地震记录R为扰动，模拟得到逆时“近似直达”波场Ur；
[0016] (6)计算逆时“近似直达”波场能量Er，其具体表达式如下：
[0017]
[0018] 式(3)中，Er为逆时“近似直达”波场能量，Ur为逆时“近似直达”波场，t为时间，T为逆时延拓时长，正演模拟时长与逆时延拓时长一致；
[0019] (7)输入本次迭代计算得到的梯度值gi，对其进行预条件处理，其表达式为：
[0020]
[0021] 式(4)中，gi为本次迭代计算得到的梯度值，Es为正时“近似直达”波场能量，Er为逆时“近似直达”波场能量，Gi为基于直达波场梯度预条件处理后的梯度值；
[0022] 最小二乘逆时偏移每次迭代中的梯度预条件处理均需要重复(2)～(7)6个步骤。
[0024] 相对于常规方法而言，基于直达波场的最小二乘逆时偏移梯度预条件方法显著降低了地震波能量中反射波、多次波能量对于梯度预条件处理的影响，使得地下照明补偿因子与几何扩散效应的校正因子更为精确。本发明即在最小二乘逆时偏移的每次迭代过程中，首先基于背景速度模型与反射系数模型，采用声波方程有限差分线性化正演模拟技术得到地震记录，然后采用无反射声波方程有限差分正演模拟技术，模拟得到以炮点震源子波为扰动的正时“近似直达”波场能量，和以上述线性化正演模拟得到的地震记录为扰动的逆时“近似直达”波场能量，最后应用此正时“近似直达”波场能量和逆时“近似直达”波场能量对梯度进行预条件处理。该方法可有效克服常规基于地震波能量梯度预条件方法中波场能量不准确的缺陷，实现地下介质的高精度成像，从而为后续的地震资料解释提供精确的构造信息。
附图说明
[0025] 图1为基于常规地震波能量的最小二乘逆时偏移梯度预条件方法
流程图；
[0026] 图2为基于直达波场的最小二乘逆时偏移梯度预条件方法流程图；
[0027] 图3为水平层状介质真实速度模型；
[0028] 图4为水平层状介质背景速度模型；
[0029] 图5为水平层状介质理论反射系数模型；
[0030] 图6为基于常规地震波能量梯度预条件的最小二乘逆时偏移成像结果；
[0031] 图7为基于直达波场梯度预条件的最小二乘逆时偏移成像结果；
[0032] 图8为两种梯度预条件方法的偏移结果与理论反射系数模型抽道对比图(X＝1500m)。
[0033] 图9为两种梯度预条件方法的偏移结果与理论反射系数模型抽道对比图(中深层放大显示，X＝1500m)。
具体实施方式
[0034] 下面通过实例结合附图来对本发明的技术方案作进一步解释，但本发明的保护范围不受任何形式上的限制。
[0036] 本发明提出的一种基于直达波场的最小二乘逆时偏移梯度预条件方法，其实现流程如图2所示。本发明以水平层状介质模型(水平层状介质真实速度模型如图3所示)来模拟野外真实的地下介质情况，水平层状介质模型共四层，第一层速度为2500m/s，第二层速度为3000m/s，第三层速度为3500m/s，第四层速度为3800m/s，模型横向长度为3000m，纵向长度为2000m。对水平层状介质速度模型进行高斯平滑，得到背景模型(水平层状介质背景速度模型如图4所示)，由真实速度模型与背景速度模型即可得到理论反射系数模型(水平层状介质理论反射系数模型如图5所示)。
[0037] 下面详细阐述本发明具体实施方式：
[0038] (1)基于上述水平层状介质速度模型，采用全排列的观测方式进行地震数据采集，共61炮，每炮401道接收，炮间隔为50m，道间隔为5m，炮点与检波点均在地表，深度为0m。基于速度模型与采集的地震数据进行最小二乘逆时偏移成像。
[0039] (2)在最小二乘逆时偏移的第i(i>＝1)次迭代计算中，输入上次迭代更新后的反射系数模型Mi-1(第1次迭代时输入给定的反射系数模型M0)，基于背景速度模型v0、反射系数模型Mi-1(第1次迭代时为给定的反射系数模型M0)，采用声波方程有限差分线性化正演模拟技术，以震源子波w为扰动，模拟得到地震记录R；
[0040] (3)基于背景速度模型v0，采用无反射声波方程有限差分正演模拟技术，以震源子波w为扰动，模拟得到正时“近似直达”波场Us；无反射声波方程的具体表达式为：
[0041]
[0042] 式(1)中，P为
应力，x、z分别为空间坐标，v为声波速度，t为时间，f为震源扰动项。
[0043] (4)计算正时“近似直达”波场能量Es，其具体表达式如下：
[0044]
[0045] 式(2)中，Es为正时“近似直达”波场能量，Us为正时“近似直达”波场，t为时间，T为正演模拟时长；
[0046] (5)基于背景速度模型v0，采用无反射声波方程有限差分正演模拟技术，无反射声波方程的具体表达式如式(1)所示，然后以第(2)步中得到的地震记录R为扰动，模拟得到逆时“近似直达”波场Ur；
[0047] (6)计算逆时“近似直达”波场能量Er，其具体表达式如下：
[0048]
[0049] 式(3)中，Er为逆时“近似直达”波场能量，Ur为逆时“近似直达”波场，t为时间，T为逆时延拓时长，正演模拟时长与逆时延拓时长一致；
[0050] (7)输入本次迭代计算得到的梯度值gi，对其进行预条件处理，其表达式为：
[0051]
[0052] 式(4)中，gi为本次迭代计算得到的梯度值，Es为正时“近似直达”波场能量，Er为逆时“近似直达”波场能量，Gi为基于直达波场梯度预条件处理后的梯度值；
[0053] 最小二乘逆时偏移每次迭代中的梯度预条件处理均需要重复(2)～(7)6个步骤。
[0054] 图6为基于常规地震波能量梯度预条件的最小二乘逆时偏移成像结果，图7为基于直达波场梯度预条件的最小二乘逆时偏移成像结果。对比图6、图7可知，基于直达波场梯度预条件的最小二乘逆时偏移可获得更优的成像效果，尤其在中深层，其横向能量更为均衡，分辨率更高。为了更清晰地对比两种梯度预条件方法的偏移结果，将两个偏移结果与理论反射系数模型进行抽道对比，如图8、图9所示。由图8、图9可知，基于直达波场梯度预条件的最小二乘逆时偏移结果与理论反射系数模型更为接近，其保幅性更好。因此，采用本发明的梯度预条件方法来进行最小二乘逆时偏移可更为精确地求取地下反射系数模型，能够获得高精度、高保幅、高分辨率的成像结果，从而为后续的实际地震资料解释提供精确的构造信息。
