[發(fā)明專(zhuān)利]一種基于譜反演的砂體厚度預(yù)測(cè)方法及系統(tǒng)

在審

• 申請(qǐng)?zhí)枺?/td>: 201610143480.6
• 申請(qǐng)日: 2016-03-14
• 公開(kāi)（公告）號(hào): CN107193046A
• 公開(kāi)（公告）日: 2017-09-22
• 發(fā)明（設(shè)計(jì)）人: 鄔興威;司朝年;韓東;劉坤巖;解麗慧;王萍
• 申請(qǐng)（專(zhuān)利權(quán)）人: 中國(guó)石油化工股份有限公司;中國(guó)石油化工股份有限公司石油勘探開(kāi)發(fā)研究院
• 主分類(lèi)號(hào): G01V1/30
• 分類(lèi)號(hào): G01V1/30
• 代理公司: 北京聿宏知識(shí)產(chǎn)權(quán)代理有限公司11372
• 代理人: 張文娟,朱繪

1.一種基于譜反演的砂體厚度預(yù)測(cè)方法，包括：

從原始地震數(shù)據(jù)中獲取地震子波和地層反射系數(shù)；

建立目標(biāo)函數(shù)，并根據(jù)所述目標(biāo)函數(shù)、所述地震子波、奇偶分解后的地層反射系數(shù)進(jìn)行譜反演，以獲取偶分量反射系數(shù)序列；

根據(jù)所述地震子波、所述偶分量反射系數(shù)序列和所述目標(biāo)函數(shù)合成新地震數(shù)據(jù)；

根據(jù)所述新地震數(shù)據(jù)預(yù)測(cè)砂體厚度。

2.根據(jù)權(quán)利要求1所述的預(yù)測(cè)方法，其特征在于，獲取所述偶分量反射系數(shù)序列進(jìn)一步包括：

建立包括地震記錄、地震子波、奇偶分解后的地層反射系數(shù)的目標(biāo)函數(shù)；

在指定頻率范圍內(nèi)，根據(jù)所述目標(biāo)函數(shù)對(duì)每道地震記錄進(jìn)行運(yùn)算，以確定對(duì)應(yīng)地層反射系數(shù)的偶分量反射系數(shù)和奇分量反射系數(shù)的位置點(diǎn)和幅值大小；

反復(fù)迭代修改偶分量反射系數(shù)和奇分量反射系數(shù)的位置點(diǎn)和幅值直到滿(mǎn)足預(yù)設(shè)要求，以確定偶分量反射系數(shù)序列。

3.根據(jù)權(quán)利要求2所述的預(yù)測(cè)方法，其特征在于，所述目標(biāo)函數(shù)為：

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為偶分量反射系數(shù)序列；r_o為奇分量反射系數(shù)序列；f_L為低截頻；f_H為高截頻；N為時(shí)窗內(nèi)包含的地層反射系數(shù)的個(gè)數(shù)；Ti為第i層與第N-i+1層之間的間隔；α_e和α_o為加權(quán)系數(shù)，t為采樣時(shí)間，f為主頻，Re取實(shí)部，Im取虛部，T為地震道計(jì)算時(shí)間，r_e(I,N-i+1)和r_o(I,N-i+1)分別表示第i層與第N-i+1層之間的偶分量反射系數(shù)序列和奇分量反射系數(shù)序列。