[发明专利]一种基于谱反演的砂体厚度预测方法及系统

权利要求书

1.一种基于谱反演的砂体厚度预测方法，包括：

从原始地震数据中获取地震子波和地层反射系数；

建立目标函数，并根据所述目标函数、所述地震子波、奇偶分解后的地层反射系数进行谱反演，以获取偶分量反射系数序列；

根据所述地震子波、所述偶分量反射系数序列和所述目标函数合成新地震数据；

根据所述新地震数据预测砂体厚度。

2.根据权利要求1所述的预测方法，其特征在于，获取所述偶分量反射系数序列进一步包括：

建立包括地震记录、地震子波、奇偶分解后的地层反射系数的目标函数；

在指定频率范围内，根据所述目标函数对每道地震记录进行运算，以确定对应地层反射系数的偶分量反射系数和奇分量反射系数的位置点和幅值大小；

反复迭代修改偶分量反射系数和奇分量反射系数的位置点和幅值直到满足预设要求，以确定偶分量反射系数序列。

3.根据权利要求2所述的预测方法，其特征在于，所述目标函数为：

O(re,ro,T)=Σf=fLfH[αe|Re[S(t,f)/W(t,f)]-Σi=1N/22re(i,N-i+1)cos(πfTi)|]+Σf=fLfH[αo|Im[S(t,f)/W(t,f)]-Σi=1N/22ro(i,N-i+1)sin(πfTi)|]]]>

其中，S(t,f)为地震记录；W(t,f)为地震子波；r_e为偶分量反射系数序列；r_o为奇分量反射系数序列；f_L为低截频；f_H为高截频；N为时窗内包含的地层反射系数的个数；Ti为第i层与第N-i+1层之间的间隔；α_e和α_o为加权系数，t为采样时间，f为主频，Re取实部，Im取虚部，T为地震道计算时间，r_e(I,N-i+1)和r_o(I,N-i+1)分别表示第i层与第N-i+1层之间的偶分量反射系数序列和奇分量反射系数序列。