[發(fā)明專利]一種全波形速度場(chǎng)約束下的鹽下碳酸鹽巖儲(chǔ)層反演方法

有效

• 申請(qǐng)?zhí)枺?/td>: 202110483109.5
• 申請(qǐng)日: 2021-04-30
• 公開（公告）號(hào): CN113109875B
• 公開（公告）日: 2022-12-20
• 發(fā)明（設(shè)計(jì)）人: 董洪超;王宗俊;范廷恩;高云峰;蔡文濤;聶妍;劉向南;樊鵬軍;田楠;周建楠
• 申請(qǐng)（專利權(quán)）人: 中海石油（中國(guó)）有限公司;中海石油（中國(guó)）有限公司北京研究中心
• 主分類號(hào): G01V1/50
• 分類號(hào): G01V1/50
• 代理公司: 北京紀(jì)凱知識(shí)產(chǎn)權(quán)代理有限公司 11245
• 代理人: 王灝增

本發(fā)明涉及一種全波形速度場(chǎng)約束下的鹽下碳酸鹽巖儲(chǔ)層反演方法，其包括以下步驟：1)目標(biāo)工區(qū)巖石物理分析及井震標(biāo)定，2)全波形速度場(chǎng)約束下低頻模型的構(gòu)建，3)基于部分疊加地震數(shù)據(jù)進(jìn)行疊前同時(shí)反演，4)反演結(jié)果檢驗(yàn)及分析。本發(fā)明在求取低頻模型時(shí)，采用全波形速度場(chǎng)約束下的多井插值法進(jìn)行低頻模型構(gòu)建，準(zhǔn)確補(bǔ)充地震所缺低頻成分；同時(shí)基于巖石物理分析，確定IP?IS參數(shù)(縱波阻抗減橫波阻抗)為反演敏感彈性參數(shù)；通過疊前反演，得到IP?IS彈性參數(shù)體，進(jìn)一步提升鹽下碳酸鹽巖儲(chǔ)層預(yù)測(cè)精度。

技術(shù)領(lǐng)域

本發(fā)明涉及地質(zhì)領(lǐng)域，具體的說(shuō)，涉及一種全波形速度場(chǎng)約束下的鹽下碳酸鹽巖儲(chǔ)層反演方法。

背景技術(shù)

目前，反演已廣泛應(yīng)用于勘探階段的儲(chǔ)層定性或定量預(yù)測(cè)和開發(fā)階段的儲(chǔ)量計(jì)算、井網(wǎng)部署、油藏動(dòng)態(tài)監(jiān)測(cè)等方面，已成為地下儲(chǔ)層刻畫過程中必不可少的方法之一。反演過程中，通過低頻模型補(bǔ)償?shù)卣鹳Y料中缺失的低頻信息是開展巖性識(shí)別、儲(chǔ)層定量解釋的重要步驟。一般情況下，低頻模型是以地震解釋層位和沉積規(guī)律為約束，將測(cè)井信息在整個(gè)數(shù)據(jù)體范圍內(nèi)進(jìn)行內(nèi)插和外推得到的，一定程度上補(bǔ)償了地震資料中缺失的低頻信息。但是，在井間儲(chǔ)層橫向變化劇烈、井點(diǎn)位置地震反射能量差異較大的區(qū)域內(nèi)，補(bǔ)充的低頻信息反而會(huì)給反演結(jié)果定量解釋帶來(lái)誤差。

碳酸鹽巖儲(chǔ)層在全球油氣分布中占有重要地位，碳酸鹽巖儲(chǔ)集層構(gòu)成的油氣田儲(chǔ)量大、單井產(chǎn)量高，容易形成大型油氣田。和碎屑巖儲(chǔ)層相比，碳酸鹽巖儲(chǔ)層更為復(fù)雜，空間分布隨機(jī)性強(qiáng)，突變性大，儲(chǔ)層縱向巨厚，多期發(fā)育，橫向沉積相帶變化快，非均質(zhì)性強(qiáng)，儲(chǔ)層與圍巖阻抗差異不明顯。同時(shí)，鹽丘和圍巖間存在巨大的速度差異，造成鹽下地震波場(chǎng)十分復(fù)雜，地震資料的分辨率和信噪比低，導(dǎo)致鹽下碳酸鹽巖儲(chǔ)層儲(chǔ)層識(shí)別和預(yù)測(cè)難度大。

發(fā)明內(nèi)容

針對(duì)上述問題，本發(fā)明的目的是提供一種全波形速度場(chǎng)約束下的鹽下碳酸鹽巖儲(chǔ)層反演方法，基于全波形反演得到的速度場(chǎng)約束建立低頻模型，補(bǔ)充地震所缺低頻成分；同時(shí)基于巖石物理分析，確定IP-IS參數(shù)(縱波阻抗-減橫波阻抗參數(shù))為反演敏感彈性參數(shù)；通過疊前反演，得到IP-IS彈性參數(shù)體，實(shí)現(xiàn)鹽下碳酸鹽巖儲(chǔ)層精細(xì)預(yù)測(cè)。

為實(shí)現(xiàn)上述目的，本發(fā)明提供一種全波形速度場(chǎng)約束下的鹽下碳酸鹽巖儲(chǔ)層反演方法，對(duì)于目標(biāo)工區(qū)基于全波形反演和巖石物理分析的結(jié)果進(jìn)行疊前反演得到得到IP-IS彈性參數(shù)體(即縱波阻抗-減橫波阻抗參數(shù)體)，實(shí)現(xiàn)鹽下碳酸鹽巖儲(chǔ)層精細(xì)預(yù)測(cè)。

所述的全波形速度場(chǎng)約束下的鹽下碳酸鹽巖儲(chǔ)層反演方法，包括以下步驟：

1)目標(biāo)工區(qū)巖石物理分析及井震標(biāo)定；

2)全波形速度場(chǎng)約束下低頻模型的構(gòu)建，其中，全波形速度場(chǎng)是針對(duì)目標(biāo)工區(qū)的特有的低頻模型建立方法；

3)基于部分疊加地震數(shù)據(jù)對(duì)目標(biāo)工區(qū)進(jìn)行疊前反演；

4)將步驟3)得到的反演結(jié)果進(jìn)一步處理，得到IP-IS彈性參數(shù)體，作為反演的最終結(jié)果，對(duì)儲(chǔ)層進(jìn)行精細(xì)預(yù)測(cè)。

所述步驟1)中還包括，在測(cè)井資料一致性處理及橫波預(yù)測(cè)的基礎(chǔ)上，進(jìn)行巖石物理分析，明確能夠區(qū)分儲(chǔ)層及非儲(chǔ)層的反演敏感彈性參數(shù)；根據(jù)測(cè)井?dāng)?shù)據(jù)、部分疊加地震數(shù)據(jù)進(jìn)行井震標(biāo)定，提取多井子波。

所述步驟2)中采用全波形速度場(chǎng)約束下的多井插值法進(jìn)行低頻模型構(gòu)建，具體過程包括：

21)首先，速度場(chǎng)應(yīng)為全波形反演得到的速度場(chǎng)，并對(duì)其進(jìn)行去野值、平滑、插值處理；

22)提取目的層段測(cè)井縱波阻抗曲線，對(duì)其低通濾波，保留10hz以下有效成分，分別統(tǒng)計(jì)目的層段全波形速度以及低通濾波后測(cè)井縱波阻抗的均值和方差，利用正態(tài)分布參數(shù)估計(jì)公式，將全波形速度場(chǎng)轉(zhuǎn)化為測(cè)井縱波阻抗數(shù)據(jù)體；

23)應(yīng)用簡(jiǎn)單克里格法消除井點(diǎn)位置的縱波阻抗預(yù)測(cè)誤差，得到縱波阻抗低頻模型；