[發(fā)明專利]一種碳酸鹽巖層間巖溶儲層預測方法及裝置

有效

• 申請?zhí)枺?/td>: 202010250672.3
• 申請日: 2020-04-01
• 公開（公告）號: CN111427085B
• 公開（公告）日: 2023-02-10
• 發(fā)明（設計）人: 常少英;胡安平;姚根順;沈安江;魯慧麗;曹鵬;鄭劍鋒;曹曉初;邵冠銘;陳婭娜;張豪;丁振純
• 申請（專利權）人: 中國石油天然氣股份有限公司
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 北京三友知識產(chǎn)權代理有限公司 11127
• 代理人: 張德斌;閆加賀

1.一種碳酸鹽巖層間巖溶儲層預測方法，其特征在于，包括如下步驟：

S101、獲取地震純波數(shù)據(jù)和目的層的相對應時間構(gòu)造層位數(shù)據(jù)；

S102、通過對地震純波數(shù)據(jù)進行相對波阻抗反演，獲得相對波阻抗；

S103、通過對地震純波數(shù)據(jù)進行地層吸收因子計算，獲得地層吸收因子數(shù)據(jù)體；

S104、構(gòu)建相對波阻抗與地層吸收因子數(shù)據(jù)體的融合屬性計算公式；

S105、利用最小二乘法優(yōu)化融合屬性計算公式，獲得優(yōu)化后的融合屬性；

S106、利用優(yōu)化后的融合屬性對儲層進行預測；

步驟S104中所述相對波阻抗與地層吸收因子數(shù)據(jù)體的融合屬性計算公式具體為：

式II中，

F(t)表示融合屬性；

Ip(t)表示相對波阻抗；

Q表示吸收因子；

α表示低頻數(shù)值；

β表示頻率滑動視窗；

步驟S105中所述利用最小二乘法優(yōu)化融合屬性計算公式，獲得優(yōu)化后的融合屬性，具體包括：

統(tǒng)計實際鉆井儲層測井的相對波阻抗和地層吸收因子數(shù)據(jù)體，利用最小二乘法優(yōu)化融合屬性計算公式，獲得油氣層與水層分界門檻值以及優(yōu)化后的頻率滑動視窗；

根據(jù)優(yōu)化后的頻率滑動視窗，獲得優(yōu)化后的地層吸收因子數(shù)據(jù)體；

將優(yōu)化后的地層吸收因子數(shù)據(jù)體結(jié)合融合屬性計算公式，獲得優(yōu)化后的融合屬性。

2.根據(jù)權利要求1所述的碳酸鹽巖層間巖溶儲層預測方法，其特征在于，步驟S101還包括獲取其他數(shù)據(jù)，所述其他數(shù)據(jù)包括測井聲波、密度、電阻率、伽馬和孔隙度數(shù)據(jù)。

3.根據(jù)權利要求1所述的碳酸鹽巖層間巖溶儲層預測方法，其特征在于，步驟S102中所述通過對地震純波數(shù)據(jù)進行相對波阻抗反演，獲得相對波阻抗，具體為：

基于褶積模型，通過對地震純波數(shù)據(jù)利用道積分進行相對波阻抗反演，得到相對波阻抗的計算公式如式I，根據(jù)式I計算獲得相對波阻抗；

ln Ip(t)＝2∫S(t)dt I；

式I中，

S(t)表示地震信號；

Ip(t)表示相對波阻抗；

t表示時間。

4.根據(jù)權利要求1所述的碳酸鹽巖層間巖溶儲層預測方法，其特征在于，所述利用最小二乘法優(yōu)化融合屬性計算公式，獲得優(yōu)化后的融合屬性，具體通過求解下式III進行優(yōu)化：

式III中，

m表示實際鉆井儲層測井的數(shù)目；

Ip(t)表示相對波阻抗；

Q表示吸收因子；

α表示低頻數(shù)值；

β表示頻率滑動視窗；

ε表示油氣層與水層分界門檻值。

5.根據(jù)權利要求1所述的碳酸鹽巖層間巖溶儲層預測方法，其特征在于，步驟S106中所述利用優(yōu)化后的融合屬性對儲層進行預測，具體為將優(yōu)化后的融合屬性和油氣層與水層分界門檻值進行比較，得出預測結(jié)論。

6.根據(jù)權利要求5所述的碳酸鹽巖層間巖溶儲層預測方法，其特征在于，所述優(yōu)化后的融合屬性大于油氣層與水層分界門檻值的地方為含油氣地區(qū)。