[发明专利]一种基于范数方程提高分辨率的反褶积方法

权利要求书

1.一种基于范数方程提高分辨率的反褶积方法，其特征在于，包括以下步骤：

步骤1、建立符合地震记录的普通模型，

步骤2、设定所述地震记录普通模型中的噪音为零，得到地震记录标准模型，

步骤3、通过所述标准模型，得到频率模型，并通过所述频率模型得到实际振幅模型和实际振幅谱，

步骤4、利用L1/L2范数方程，设定标准振幅谱，并求取所述实际振幅谱与所述标准振幅谱的误差能量，

步骤5、求取误差能量最小时的实际振幅谱，对所述实际振幅谱进行傅里叶反变换得到对应的时间函数，

步骤6、将所述时间函数带入所述标准模型，得到地表一致性反褶积因子，将所述反褶积因子与地震道褶积完成反褶积处理。

2.根据权利要求1所述的反褶积方法，其特征在于，步骤1中所述的地震记录普通模型为：

x_ijk(t)＝ξ_ijk(t)*s_i(t)*g_j(t)*m_(i+j)/2(t)*o_k(t)+n_ijk(t)

其中：s_i(t)为第i炮位置的子波，g_j(t)为第j个检波点位置的子波，m_(i+j)/2(t)为共中心点的子波，o_k(t)为第k个偏移距对应的子波，x_ijk(t)为对应的地震记录、ξ_ijk(t)为反射系数、n_ijk(t)为噪音。

3.根据权利要求2所述的反褶积方法，其特征在于，步骤2中所述的地震记录标准模型为：

x_ijk(t)＝ξ_ijk(t)*s_i(t)*g_j(t)*m_(i+j)/2(t)*o_k(t)

其中：s_i(t)为第i炮位置的子波，g_j(t)为第j个检波点位置的子波，m_(i+j)/2(t)为共中心点的子波，o_k(t)为第k个偏移距对应的子波，x_ijk(t)为对应的地震记录、ξ_ijk(t)为反射系数。

4.根据权利要求3所述的反褶积方法，其特征在于，步骤3具体包括以下步骤：

步骤3.1、对所述标准模型进行傅里叶变换，转换到频率域，得到频率模型，

步骤3.2、设定该模型对应的子波为最小相位，将所述频率模型分解到振幅谱，得到实际振幅谱，

步骤3.3、对所述实际振幅谱取对数得到实际振幅模型。

5.根据权利要求4所述的反褶积方法，其特征在于，步骤3.3中所述的实际振幅模型为：

lnA_ijk＝lnA_s+lnA_g+lnA_m+lnA_o

其中：A_ijk为实际振幅谱，A_s为该炮点位置的子波的振幅谱，A_g为对应的检波点位置子波的振幅谱，A_m为共中心点子波的振幅谱，A_o为偏移距对应的子波的振幅谱。

6.根据权利要求5所述的反褶积方法，其特征在于，步骤4中所述的误差能量其中，A_ijk为实际振幅谱，为模型振幅谱。

7.根据权利要求6所述的反褶积方法，其特征在于，步骤5具体包括以下步骤：

步骤5.1、令求得各分量对数谱，其中A_S为中心点子波的实际振幅谱，A_g为检波点子波实际振幅，A_m为共中心点子波的实际振幅，A_o为偏移距对应子波的实际振幅，

步骤5.2、当误差能量最小时，采用雅克比迭代法求取中心点子波、检波点子波、共中心点的子波、偏移距对应的子波的振幅谱，

步骤5.3、对所述中心点子波、检波点子波、共中心点的子波、偏移距对应的子波的振幅谱进行傅里叶反变换得到对应的时间函数s_i(t)、g_j(t)、m_(i+j)/2(t)、o_k(t)，其中：s_i(t)为第i炮位置的子波，g_j(t)为第j个检波点位置的子波，m_(i+j)/2(t)为共中心点的子波，o_k(t)为第k个偏移距对应的子波。