CN102879821A - 一种针对地震叠前道集的同相轴精细拉平处理方法 - Google Patents

Abstract

本发明涉及地球物理勘探中的地震叠前道集处理方法，尤其涉及一种针对地震叠前道集的同相轴精细拉平处理方法，包括对采集的地震数据进行前期处理；抽取想要进行同相轴精细拉平处理的叠前道集，用I行、J列的二维数组D来表示；设定同相轴精细拉平处理参数；以N_w为窗口大小，以N_m为窗口移动量，逐窗口判定是否在某个波形对比窗口的中心点设置拉平种子点，并计算拉平种子点移动量；计算拉伸采样坐标；根据拉伸采样坐标实施同相轴精细拉平。通过本发明实施例的方法，该方法对于地震资料的处理和解释都具有重要的意义，并具有稳定、高效的特点。

Description

一种针对地震叠前道集的同相轴精细拉平处理方法

技术领域

本发明涉及地球物理勘探中的地震叠前道集处理方法，尤其涉及一种针对地震叠前道集的同相轴精细拉平处理方法。

背景技术

地震勘探是利用地下介质的弹性差异，通过观测和分析大地对人工激发地震波的响应，推断地下岩层的性质和形态的地球物理勘探方法。地震勘探是钻探前勘测石油与天然气资源的重要手段，在煤田和工程地质勘查、区域地质研究和地壳研究等方面，也得到广泛应用。

地震勘探主要分为地震资料采集、处理与解释三个阶段。在地震资料采集阶段，多次覆盖已成为最基本的野外采集方法，所谓的多次覆盖是指，在地表适当地选取多个对地震波激发点和接收点，进而得到来自地下同一反射点的信号，将这些信号组合在一起就构成了一个共反射点道集（一种地震叠前道集）；在地震资料处理阶段，以同相轴拉平后的共反射点道集为基础，对其进行叠加求和（共反射点叠加）可以提高地震资料的信噪比；在地震资料解释阶段，以同相轴拉平后的共反射点道集为基础，通过分析振幅随入射角的变化（AVA分析）可以得到地下的各项同性弹性参数信息，通过分析振幅随方位角的变化（AVAz分析）可以得到地下的各项异性弹性参数信息。

无论是地震资料处理阶段的共反射点叠加，还是地震资料解释阶段的AVA和AVAz分析，都是以同相轴拉平后的共反射点道集为基础的，同相轴不平会影响后续工作的质量，甚至得到错误的结果。处理上，常根据共反射点道集的时距曲线方程进行同相轴拉平（动校正），虽然动校正可以保证共反射点道集同相轴基本拉平，但受各种因素（如由于地表高程起伏剧烈导致的静校正不准确、由于水平层状各项同性介质假设条件导致的常规时距曲线方程误差）的影响，动校正往往无法实现共反射点道集同相轴精细拉平，继而，降低了共反射点叠加的信噪比，并导致由AVA和AVAz分析得到的弹性参数不准确甚至错误。

发明内容

为了解决现有技术中同相轴不平的问题，从而造成地震勘探结果不准确的问题，提供了一种针对地震叠前道集的同相轴精细拉平处理方法。

本发明实施例提供的一种针对地震叠前道集的同相轴精细拉平处理方法，包括，

步骤101，对采集的地震数据进行前期处理；

步骤102，抽取想要进行同相轴精细拉平处理的叠前道集，用I行、J列的二维数组D来表示，其中，I为每个地震道的采样点数，J为叠前道集的地震道数，二维数组D，其第i行、第j列的值D(i,j)表示叠前道集的第j个地震道的第i个采样点的振幅值；

步骤103，设定同相轴精细拉平处理参数：N_w、N_m、N_s、Tol_g、Tol_i，其中，所述的N_w为波形对比窗口的大小，N_w为整数；所述的N_m为波形对比窗口的移动量，N_m为整数；所述的N_s为波形对比的搜索半径，N_s为整数；所述的Tol_g为群容差，是判定是否在某个波形对比窗口的中心点设置拉平种子点的依据；所述的Tol_i为个体容差，是优选某个地震道的相似地震道的依据；

步骤104，以N_w为窗口大小，以N_m为窗口移动量，逐窗口判定是否在某个波形对比窗口的中心点设置拉平种子点，并计算拉平种子点移动量；

步骤105，计算拉伸采样坐标；

步骤106，根据拉伸采样坐标实施同相轴精细拉平。

根据所述一种针对地震叠前道集的同相轴精细拉平处理方法的一个进一步的方面，对采集的地震数据进行前期处理中包括，对采集的地震数据进行去噪、静校正、动校正的前期处理。

根据所述一种针对地震叠前道集的同相轴精细拉平处理方法的再一个进一步的方面，所述叠前道集包括共反射点道集和方位角道集。

根据所述一种针对地震叠前道集的同相轴精细拉平处理方法的另一个进一步的方面，所述逐窗口判定是否在某个波形对比窗口的中心点设置拉平种子点，并计算拉平种子点移动量具体包括，

步骤401，对于第k个波形对比窗口，计算地震道与地震道之间最优移动量的二维数组S，其中二维数组S中的第j₁行、第j₂列的值S(j₁,j₂)表示叠前道集的第j₁个地震道和第j₂个地震道的最优移动量，并计算地震道与地震道之间最大相关系数的二维数组C，其中二维数组C中的第j₁行、第j₂列的值C(j₁,j₂)表示叠前道集的第j₁个地震道和第j₂个地震道的最大相关系数；所述的地震道与地震道之间最优移动量S(j₁,j₂)和地震道与地震道之间最大相关系数C(j₁,j₂)的计算方法为：设i_s＝(k-1)N_m，i_e＝i_s+N_w，定义集合

定义长度为（N_w+1）的向量f，f的第i个元素f(i)＝D(i_s+i-1,j₁)，

定义2N_s+1个长度为（N_w+1）的向量g_l，

g_l的第i个元素g_l(i)＝D(i_s+l+i-1,j₂)，则最大相关系数C(j₁,j₂)由以下公式给出：

$C(j_1,j_2)=\underset{l\∈\hat{L}}{\max }\frac{f\·g_l}{\sqrt{f\·f}\sqrt{g_l\·g_l}}$

最优移动量S(j₁,j₂)，即取得C(j₁,j₂)时的l值，由以下公式给出：

$S(j_1,j_2)=\arg \underset{l\∈\hat{L}}{\max }\frac{f\·g_l}{\sqrt{f\·f}\sqrt{g_l\·g_l}}$

步骤402，对于第k个波形对比窗口，判定是否在该波形对比窗口的中心点设置拉平种子点，其中，所述的判定是否设置拉平种子点的方法为：求取二维数组C的所有元素的平均值c_m：

$c_m=\frac{1}{J^2}\underset{j_1=1}{\overset{J}{\mathrm{\Σ}}}\underset{j_2=1}{\overset{J}{\mathrm{\Σ}}}C(j_1,j_2)$

若c_m≥Tol_g，则需要在该波形对比窗口的中心点设置拉平种子点，记录该波形对比窗口的中心点 $i_c=\mathrm{round}\left(\frac{i_s+i_e}{2}\right)=\mathrm{round}((k-1)N_m+\frac{N_w}{2}),$ 其中函数round(x)表示对x进行四舍五入取整，转步骤103，若c_m＜Tol_g，则不需要在该波形对比窗口的中心点设置拉平种子点，记k＝k+1，转步骤101继续；

步骤403，对于第k个波形对比窗口，计算拉平种子点移动量，定义长度为J的向量m用来存储每个地震道的拉平种子点移动量，其中，所述的计算拉平种子点移动量的方法为：

选取基准地震道，定义集合

对于任意一个定值j₁，

第j₁道的相似道集合为 ${\hat{J}}_{j_1}=\left\{j_2\right|C(j_1,j_2)\≥{\mathrm{Tol}}_i,j_2\∈\hat{J}\},$ 则基准地震道为 $j_b=\arg \underset{j_1\∈\hat{J}}{\max }\mathrm{card}\left({\hat{J}}_{j_1}\right),$ 其中函数表示集合的元素的个数；

优选基准地震道的相似地震道，相似地震道集合为

计算每个地震道的拉平种子点移动量，

令 $j_{\min }=\underset{j\∈{\hat{J}}_{j_b}}{\min }j,$ $j_{\max }=\underset{j\∈{\hat{J}}_{j_b}}{\max }j,$ 则

根据所述一种针对地震叠前道集的同相轴精细拉平处理方法的另一个进一步的方面，在计算拉伸采样坐标中，定义I行、J列的二维数组X_new用来存储拉伸采样坐标，所述的计算拉伸采样坐标的方法为：假设需要在Q个波形对比窗口的中心点设置拉平种子点，设其中第q个中心点以及它所对应的拉平种子点移动量分别为i_q和

q＝1,2,…,Q，定义集合

则对于任意一个定值j，

根据所述一种针对地震叠前道集的同相轴精细拉平处理方法的另一个进一步的方面，根据拉伸采样坐标实施同相轴精细拉平进一步包括，定义I行、J列的二维数组X用来存储原始采样坐标，定义I行、J列的二维数组D^new用来存储同相轴精细拉平处理后的叠前道集；

设置原始采样坐标，具体地，X(i,j)＝i,i＝1,2,…,I;j＝1,2,…,J；

对于任意一个定值j，以X的第j列X_j为自变量、以D的第j列D_j为函数值，构建三次样条插值函数

则 $D^{\mathrm{new}}(i,j)=f_{X_j,D_j}\left(X^{\mathrm{new}}(i,j)\right).$

根据所述一种针对地震叠前道集的同相轴精细拉平处理方法的另一个进一步的方面，在根据拉伸采样坐标实施同相轴精细拉平之后还包括步骤107，用D^new代替步骤102中的D，并在步骤103中设定新的同相轴精细拉平处理参数，再次执行步骤104-步骤106。

根据所述一种针对地震叠前道集的同相轴精细拉平处理方法的另一个进一步的方面，重复执行所述步骤107两至三次。

根据所述一种针对地震叠前道集的同相轴精细拉平处理方法的另一个进一步的方面，所述在步骤103中设定新的同相轴精细拉平处理参数进一步包括，按照参数设定的先后次序N_w、N_m、N_s逐渐减小，Tol_g、Tol_i保持不变或增大。

通过本发明实施例的方法，该方法对于地震资料的处理和解释都具有重要的意义，并具有稳定、高效的特点。

附图说明

此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不构成对本发明的限定。在附图中：

图1是本发明实施例提供的一种针对地震叠前道集的同相轴精细拉平处理方法流程图；

图2是本发明实施例提供的一个方位角道集；

图3是本发明实施例提供的移动波形对比窗口的示意图；

图4是本发明实施例提供的用来存储地震道与地震道之间最大相关系数的二维数组C的示意图；

图5是本发明实施例提供的用来存储地震道与地震道之间最优移动量的二维数组S的示意图；

图6是本发明实施例提供的欲设置拉平种子点的波形对比窗口的中心点，及其对应的拉平种子点移动量的示意图；

图7是本发明实施例提供的求取拉伸采样坐标的示意图一；

图8是本发明实施例提供的求取拉伸采样坐标的示意图二；

图9a是本发明实施例中拉平种子点在处理前的叠前道集上的位置示意图；

图9b是本发明实施例中经过一次同相轴精细拉平处理后的示意图；

图10a是本发明实施例提供的同相轴精细拉平处理前的效果图；

图10b是本发明实施例提供的经过一次同相轴精细拉平处理后的效果图；

图10c是本发明实施例提供的经过二次同相轴精细拉平处理后的效果图；

图10d是本发明实施例提供的经过三次同相轴精细拉平处理后的效果图。

具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合实施例和附图，对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。

如图1为本发明实施例提供的一种针对地震叠前道集的同相轴精细拉平处理方法流程图。

包括步骤101，对采集的地震数据进行去噪、静校正、动校正的前期处理操作。

通过去噪来消除面波、多次波等干扰波；通过静校正消除由于地表不平坦、激发点深度不同等问题对于时距曲线形态的影响；通过动校正使时距曲线与地下界面的形态一致。

步骤102，抽取想要进行同相轴精细拉平处理的叠前道集。这种叠前道集一般是共反射点道集（道与道的入射角不同），也可以是方位角道集（道与道的方位角不同），它们都可以用I行、J列的二维数组D来表示。如图2所示，其给出了一个方位角道集。

在本步骤中，所述I为每个地震道的采样点数，例如I＝401；所述J为叠前道集的地震道数，例如J＝6。所述的二维数组D，其第i行、第j列的值D(i,j)表示叠前道集的第j个地震道的第i个采样点的振幅值。

步骤103，设定同相轴精细拉平处理参数：N_w、N_m、N_s、Tol_g、Tol_i。

其中，所述的N_w为波形对比窗口的大小，N_w为整数，例如N_w＝50个采样点数，图2纵向单位是采样点数。

所述的N_m为波形对比窗口的移动量，N_m为整数，例如设定N_m＝17个采样点数；

所述的N_s为波形对比的搜索半径，N_s为整数，例如设定N_s＝6个采样点数；

所述的Tol_g为群容差，它是判定是否在某个波形对比窗口的中心点设置拉平种子点的依据，当波形对比窗口内地震道相关系数的平均值大于Tol_g时，则需要在该波形对比窗口的中心点设置拉平种子点，反之亦反，例如设定Tol_g＝0.75；

所述的Tol_i为个体容差，它是优选某个地震道的相似地震道的依据，当波形对比窗口内某个地震道与基准地震道的相似系数大于Tol_i时，则该地震道为基准地震道的相似地震道，例如设定Tol_i＝0.85。

步骤104，以N_w为窗口大小，以N_m为窗口移动量，逐窗口判定是否在某个波形对比窗口的中心点设置拉平种子点，并计算拉平种子点移动量。

图3给出了移动波形对比窗口的示意图，图中用梯形框标出了最初的5个波形对比窗口的位置。

具体地，对于第k个波形对比窗口（从道集的第1个采样点算起），可以通过以下步骤4.1到步骤4.3操作：

步骤4.1，对于第k个波形对比窗口，计算地震道与地震道之间的最优移动量及最大相关系数。

定义J行、J列的二维数组S用来存储地震道与地震道之间的最优移动量。定义J行、J列的二维数组C用来存储地震道与地震道之间的最大相关系数。图4给出了当k＝3时用来存储地震道与地震道之间最大相关系数的二维数组C的示意图，所述的二维数组C，其第j₁行、第j₂列的值C(j₁,j₂)表示叠前道集的第j₁个地震道和第j₂个地震道的最大相关系数。图5给出了当k＝3时用来存储地震道与地震道之间最优移动量的二维数组S的示意图，所述的二维数组S，其第j₁行、第j₂列的值S(j₁,j₂)表示叠前道集的第j₁个地震道和第j₂个地震道的最优移动量。

所述的地震道与地震道之间最优移动量和最大相关系数，其计算方法为：以计算S(j₁,j₂)和C(j₁,j₂)为例，设i_s＝(k-1)N_m，i_e＝i_s+N_w，定义集合

定义长度为（N_w+1）的向量f，f的第i个元素f(i)＝D(i_s+i-1,j₁)，

定义2N_s+1个长度为（N_w+1）的向量g_l，

g_l的第i个元素g_l(i)＝D(i_s+l+i-1,j₂)，

则最大相关系数C(j₁,j₂)由以下公式给出：

$C(j_1,j_2)=\underset{l\∈\hat{L}}{\max }\frac{f\·g_l}{\sqrt{f\·f}\sqrt{g_l\·g_l}}$

最优移动量S(j₁,j₂)，即取得C(j₁,j₂)时的l值，由以下公式给出：

$S(j_1,j_2)=\arg \underset{l\∈\hat{L}}{\max }\frac{f\·g_l}{\sqrt{f\·f}\sqrt{g_l\·g_l}}$

上式中的向量运算是指，a·b代表向量a与b的内积运算，即将向量a与b的对应元素相乘后再求和，其中a、b仅仅是一个变量代号，旨在说明向量运算。

步骤4.2，对于第k个波形对比窗口，判定是否在该波形对比窗口的中心点设置拉平种子点。

所述的判定是否设置拉平种子点的方法为：求取二维数组C的所有元素的平均值c_m：

$c_m=\frac{1}{J^2}\underset{j_1=1}{\overset{J}{\mathrm{\Σ}}}\underset{j_2=1}{\overset{J}{\mathrm{\Σ}}}C(j_1,j_2)$

如图4所示，当k＝3时（第3个波形对比窗口），c_m＝0.78。若c_m≥Tol_g，则需要在该波形对比窗口的中心点设置拉平种子点，记录该波形对比窗口的中心点其中函数round(x)表示对x进行四舍五入取整，转步骤4.3继续；若c_m＜Tol_g，则不需要在该波形对比窗口的中心点设置拉平种子点，记k＝k+1，转步骤4.1继续。

步骤4.3，对于第k个波形对比窗口，计算拉平种子点移动量。定义长度为J的向量m用来存储每个地震道的拉平种子点移动量。

所述的计算拉平种子点移动量的方法为：

步骤4.3.1，首先选取基准地震道。具体地，定义集合

对于任意一个定值j₁，

第j₁道的相似道集合为 ${\hat{J}}_{j_1}=\left\{j_2\right|C(j_1,j_2)\≥{\mathrm{Tol}}_i,j_2\∈\hat{J}\},$ 则基准地震道为 $j_b=\arg \underset{j_1\∈\hat{J}}{\max }\mathrm{card}\left({\hat{J}}_{j_1}\right),$ 其中函数

表示集合

的元素的个数。

步骤4.3.2，优选基准地震道的相似地震道，具体地，相似地震道集合为

${\hat{J}}_{j_b}=\left\{j_2\right|C(j_b,j_2)\≥{\mathrm{Tol}}_i,j_2\∈\hat{J}\};.$

步骤4.3.3，计算每个地震道的拉平种子点移动量。

令 $j_{\min }=\underset{j\∈{\hat{J}}_{j_b}}{\min }j,$ $j_{\max }=\underset{j\∈{\hat{J}}_{j_b}}{\max }j,$ 则

如图4和图5所示，当k＝3时，

每一行元素值大于Tol_i的元素个数分别为

由于在第1、2、5、6行有最大值4，因此可以取其中之一为基准地震道，这里取j_b＝1，这时

j_min＝1，j_max＝6，m(1)＝0，m(2)＝0，m(5)＝3，m(6)＝0，m(3)和m(4)则需要由前面四个点线性内插求出，因此有m(3)＝1，m(4)＝2。

步骤105，计算拉伸采样坐标，定义I行、J列的二维数组X_new用来存储拉伸采样坐标。

所述的计算拉伸采样坐标的方法为：假设在完成步骤104所述的工作后，需要在Q个波形对比窗口的中心点设置拉平种子点，设其中第q个中心点以及它所对应的拉平种子点移动量分别为i_q和

q＝1,2,…,Q，定义集合

则对于任意一个定值j，

图6给出了欲设置拉平种子点的波形对比窗口的中心点，及其对应的拉平种子点移动量。图7给出了求取拉伸采样坐标的示意图一，其中从第二行到倒数第二行对应于上述公式的第一种情况，第一行对应于于上述公式的第二种情况，最后一行对应于上述公式的第三种情况，图8给出了求取拉伸采样坐标的示意图二，以图7所示的每一列点为基础，通过线性内插，就可以求出剩余的拉伸采样坐标，最终所有点的拉伸采样坐标如图8所示，其中，左图横坐标为地震道号，左图纵坐标为原始采样坐标，由于数据点很多，这里采用灰度图方法表示拉伸采样坐标，右图的颜色条标示了不同灰度代表的拉伸采样坐标值。图9a给出了拉平种子点在处理前的叠前道集上的位置。

步骤106，根据拉伸采样坐标实施同相轴精细拉平。定义I行、J列的二维数组X用来存储原始采样坐标，定义I行、J列的二维数组D^new用来存储同相轴精细拉平处理后的叠前道集。图9b给出了经过一次同相轴精细拉平处理后的叠前道集。对比图9a和图9b可以看出处理后同相轴被拉平了。

所述的根据拉伸采样坐标实施同相轴精细拉平的方法为：步骤6.1，首先设置原始采样坐标，具体地，X(i,j)＝i,i＝1,2,…,I;j＝1,2,…,J；步骤6.2，对于任意一个定值j，

以X的第j列X_j为自变量、以D的第j列D_j为函数值，构建三次样条插值函数

则 $D^{\mathrm{new}}(i,j)=f_{X_j,D_j}\left(X^{\mathrm{new}}(i,j)\right).$

步骤107，根据同相轴精细拉平处理的需要，用步骤106中的D^new代替步骤102中的D，并在步骤103中设定新的同相轴精细拉平处理参数，再次执行步骤104-步骤106。

所述的根据同相轴精细拉平处理的需要，其目的在于提高同相轴精细拉平处理的程度，一般地，步骤107需要重复执行2-3次。本实例在第一次同相轴精细拉平处理的基础上，定新的同相轴精细拉平处理参数，又重复执行步骤107两次，图10a-图10d给出了整个拉平处理过程的效果图，可以看到，经过三次同相轴精细拉平处理，叠前道集的同相轴得到明显拉平，这对于后续的地震资料处理和解释具有重要的意义。

在本步骤107中，所述的在步骤103中设定新的同相轴精细拉平处理参数，其一般原则是：按照参数设定的先后次序N_w、N_m、N_s应该逐渐减小，Tol_g、Tol_i可以保持不变或适当增大；例如可以在本实例第二次同相轴精细拉平处理的参数为N_w＝30，N_m＝10，N_s＝3，Tol_g＝0.75，Tol_i＝0.85；本实例第三次同相轴精细拉平处理的参数为N_w＝18，N_m＝6，N_s＝2，Tol_g＝0.85，Tol_i＝0.88。

本发明实施例提供了一种针对地震叠前道集的同相轴精细拉平处理方法，核心是提出了一种稳定的道集同相轴精细拉平种子点选取方法，该方法对于地震资料的处理和解释都具有重要的意义，并具有稳定、高效的特点。

以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (9)

1.一种针对地震叠前道集的同相轴精细拉平处理方法，其特征在于包括，

步骤101，对采集的地震数据进行前期处理；

步骤102，抽取想要进行同相轴精细拉平处理的叠前道集，用I行、J列的二维数组D来表示，其中，I为每个地震道的采样点数，J为叠前道集的地震道数，二维数组D，其第i行、第j列的值D(i,j)表示叠前道集的第j个地震道的第i个采样点的振幅值；

步骤103，设定同相轴精细拉平处理参数：N_w、N_m、N_s、Tol_g、Tol_i，其中，所述的N_w为波形对比窗口的大小，N_w为整数；所述的N_m为波形对比窗口的移动量，N_m为整数；所述的N_s为波形对比的搜索半径，N_s为整数；所述的Tol_g为群容差，是判定是否在某个波形对比窗口的中心点设置拉平种子点的依据；所述的Tol_i为个体容差，是优选某个地震道的相似地震道的依据；

步骤104，以N_w为窗口大小，以N_m为窗口移动量，逐窗口判定是否在某个波形对比窗口的中心点设置拉平种子点，并计算拉平种子点移动量；

步骤105，计算拉伸采样坐标；

步骤106，根据拉伸采样坐标实施同相轴精细拉平。

2.根据权利要求1所述的方法，其特征在于，对采集的地震数据进行前期处理中包括，对采集的地震数据进行去噪、静校正、动校正的前期处理。

3.根据权利要求1所述的方法，其特征在于，所述叠前道集包括共反射点道集和方位角道集。

4.根据权利要求1所述的方法，其特征在于，所述逐窗口判定是否在某个波形对比窗口的中心点设置拉平种子点，并计算拉平种子点移动量具体包括，

步骤401，对于第k个波形对比窗口，计算地震道与地震道之间最优移动量的二维数组S，其中二维数组S中的第j₁行、第j₂列的值S(j₁,j₂)表示叠前道集的第j₁个地震道和第j₂个地震道的最优移动量，并计算地震道与地震道之间最大相关系数的二维数组C，其中二维数组C中的第j₁行、第j₂列的值C(j₁,j₂)表示叠前道集的第j₁个地震道和第j₂个地震道的最大相关系数；所述的地震道与地震道之间最优移动量S(j₁,j₂)和地震道与地震道之间最大相关系数C(j₁,j₂)的计算方法为：设i_s＝(k-1)N_m，i_e＝i_s+N_w，定义集合

定义长度为（N_w+1）的向量f，f的第i个元素f(i)＝D(i_s+i-1,j₁)，

定义2N_s+1个长度为（N_w+1）的向量g_l，

g_l的第i个元素g_l(i)＝D(i_s+l+i-1,j₂)，

则最大相关系数C(j₁,j₂)由以下公式给出：

$C(j_1,j_2)=\underset{l\∈\hat{L}}{\max }\frac{f\·g_l}{\sqrt{f\·f}\sqrt{g_l\·g_l}}$

最优移动量S(j₁,j₂)，即取得C(j₁,j₂)时的l值，由以下公式给出：

$S(j_1,j_2)=\arg \underset{l\∈\hat{L}}{\max }\frac{f\·g_l}{\sqrt{f\·f}\sqrt{g_l\·g_l}}$

步骤402，对于第k个波形对比窗口，判定是否在该波形对比窗口的中心点设置拉平种子点，其中，所述的判定是否设置拉平种子点的方法为：求取二维数组C的所有元素的平均值c_m：

$c_m=\frac{1}{J^2}\underset{j_1=1}{\overset{J}{\mathrm{\Σ}}}\underset{j_2=1}{\overset{J}{\mathrm{\Σ}}}C(j_1,j_2)$

若c_m≥Tol_g，则需要在该波形对比窗口的中心点设置拉平种子点，记录该波形对比窗口的中心点 $i_c=\mathrm{round}\left(\frac{i_s+i_e}{2}\right)=\mathrm{round}((k-1)N_m+\frac{N_w}{2}),$ 其中函数round(x)表示对x进行四舍五入取整，转步骤103，若c_m＜Tol_g，则不需要在该波形对比窗口的中心点设置拉平种子点，记k＝k+1，转步骤101继续；

步骤403，对于第k个波形对比窗口，计算拉平种子点移动量，定义长度为J的向量m用来存储每个地震道的拉平种子点移动量，其中，所述的计算拉平种子点移动量的方法为：

选取基准地震道，定义集合

对于任意一个定值j₁，

第j₁道的相似道集合为 ${\hat{J}}_{j_1}=\left\{j_2\right|C(j_1,j_2)\≥{\mathrm{Tol}}_i,j_2\∈\hat{J}\},$ 则基准地震道为 $j_b=\arg \underset{j_1\∈\hat{J}}{\max }\mathrm{card}\left({\hat{J}}_{j_1}\right),$ 其中函数

表示集合

的元素的个数；

优选基准地震道的相似地震道，相似地震道集合为

计算每个地震道的拉平种子点移动量，

令 $j_{\min }=\underset{j\∈{\hat{J}}_{j_b}}{\min }j,$ $j_{\max }=\underset{j\∈{\hat{J}}_{j_b}}{\max }j,$ 则

5.根据权利要求4所述的方法，其特征在于，在计算拉伸采样坐标中，定义I行、J列的二维数组X_new用来存储拉伸采样坐标，所述的计算拉伸采样坐标的方法为：假设需要在Q个波形对比窗口的中心点设置拉平种子点，设其中第q个中心点以及它所对应的拉平种子点移动量分别为i_q和

q＝1,2,…,Q，定义集合

则对于任意一个定值j，

6.根据权利要求4所述的方法，其特征在于，根据拉伸采样坐标实施同相轴精细拉平进一步包括，定义I行、J列的二维数组X用来存储原始采样坐标，定义I行、J列的二维数组D^new用来存储同相轴精细拉平处理后的叠前道集；

设置原始采样坐标，具体地，X(i,j)＝i,i＝1,2,…,I;j＝1,2,…,J；

对于任意一个定值j，

以X的第j列X_j为自变量、以D的第j列D_j为函数值，构建三次样条插值函数

则 $D^{\mathrm{new}}(i,j)=f_{X_j,D_j}\left(X^{\mathrm{new}}(i,j)\right).$

7.根据权利要求6所述的方法，其特征在于，在根据拉伸采样坐标实施同相轴精细拉平之后还包括步骤107，用D^new代替步骤102中的D，并在步骤103中设定新的同相轴精细拉平处理参数，再次执行步骤104-步骤106。

8.根据权利要求7所述的方法，其特征在于，重复执行所述步骤107两至三次。

9.根据权利要求7所述的方法，其特征在于，所述在步骤103中设定新的同相轴精细拉平处理参数进一步包括，按照参数设定的先后次序N_w、N_m、N_s逐渐减小，Tol_g、Tol_i保持不变或增大。

Images