[发明专利]一种针对地震叠前道集的同相轴精细拉平处理方法

权利要求书

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(j1,j2)=maxl∈L^f·glf·fgl·gl]]>

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

S(j1,j2)=argmaxl∈L^f·glf·fgl·gl]]>

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

cm=1J2Σj1=1JΣj2=1JC(j1,j2)]]>

若c_m≥Tol_g，则需要在该波形对比窗口的中心点设置拉平种子点，记录该波形对比窗口的中心点ic=round(is+ie2)=round((k-1)Nm+Nw2),]]>其中函数round(x)表示对x进行四舍五入取整，转步骤103，若c_m＜Tol_g，则不需要在该波形对比窗口的中心点设置拉平种子点，记k＝k+1，转步骤101继续；

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

选取基准地震道，定义集合对于任意一个定值j₁，第j₁道的相似道集合为J^j1={j2|C(j1,j2)≥Toli,j2∈J^},]]>则基准地震道为jb=argmaxj1∈J^card(J^j1),]]>其中函数表示集合的元素的个数；

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

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

令jmin=minj∈J^jbj,]]>jmax=maxj∈J^jbj,]]>则