[發(fā)明專利]一種湖相碳酸鹽巖孔隙度的計算方法

有效

• 申請?zhí)枺?/td>: 202110042589.1
• 申請日: 2021-01-13
• 公開（公告）號: CN112904441B
• 公開（公告）日: 2022-04-15
• 發(fā)明（設(shè)計）人: 史瑞其;張世鑫;賈懷存;杜向東;韓文明;閆青華;李丹;孫林潔
• 申請（專利權(quán)）人: 中國海洋石油集團有限公司;中海油研究總院有限責(zé)任公司
• 主分類號: G01V9/00
• 分類號: G01V9/00
• 代理公司: 北京紀(jì)凱知識產(chǎn)權(quán)代理有限公司 11245
• 代理人: 謝斌

本發(fā)明公開了一種湖相碳酸鹽巖孔隙度的計算方法，首先，獲取湖相碳酸鹽巖的地震反演模型；其次，通過多道空間約束地震反演方法提高湖相碳酸鹽巖的地震反演模型的空間表征能力，獲取湖相碳酸鹽巖的三維波阻抗參數(shù)；最后，利用BP神經(jīng)網(wǎng)絡(luò)算法實現(xiàn)湖相碳酸鹽巖三維波阻抗參數(shù)與孔隙度參數(shù)的轉(zhuǎn)化，獲取湖相碳酸鹽巖的孔隙度參數(shù)。本發(fā)明公開的湖相碳酸鹽巖孔隙度的計算方法，通過多道空間約束反演方法、BP神經(jīng)網(wǎng)絡(luò)算法等計算手段成功實現(xiàn)了將湖相碳酸鹽巖的地震資料數(shù)據(jù)中三維波阻抗參數(shù)轉(zhuǎn)化為孔隙度參數(shù)，提高了縱波阻抗的反演精度，避免了多元擬合計算孔隙度造成的誤差，能用于預(yù)測湖相碳酸鹽巖孔隙度。

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

本發(fā)明涉及油氣勘探技術(shù)領(lǐng)域，具體涉及一種湖相碳酸鹽巖孔隙度的計算方法。

背景技術(shù)

湖相碳酸鹽巖是分布最為廣泛的一類陸相碳酸鹽巖，巴西大坎波斯盆地深水鹽下裂谷期油氣資源豐富，目的層就為湖相碳酸鹽巖儲層。由于深水油氣勘探鉆井成本高，具有高孔隙度物性的湖相碳酸鹽巖儲層是實現(xiàn)高效勘探的前提，因此孔隙度預(yù)測的精度至關(guān)重要。目前，在油氣勘探階段進行孔隙度計算的方法有兩種，一是經(jīng)驗公式法，即基于孔隙度與縱波阻抗等彈性參數(shù)的擬合關(guān)系，利用多元擬合算法實現(xiàn)孔隙度的計算；二是地震反演方法，即建立包含孔隙度的反演目標(biāo)方程，通過地震反演實現(xiàn)孔隙度的計算。

由于巴西大坎波斯盆地裂谷期火山活動劇烈，湖相碳酸鹽巖儲層與火成巖交互發(fā)育，且湖相碳酸鹽巖沉積相帶橫向變化快，空間非均質(zhì)性強。巖性復(fù)雜會導(dǎo)致孔隙度與彈性參數(shù)的擬合相關(guān)度差，因此采用經(jīng)驗公式法、利用常規(guī)多元擬合計算會產(chǎn)生較大誤差。常規(guī)地震反演的方法主要是針對“層狀”碎屑巖儲層，對于具有強非均質(zhì)性的“塊狀”湖相碳酸鹽巖儲層的空間表征能力不足，直接反演精度較低，因此利用常規(guī)地震反演計算孔隙度的方法不適用于巴西大坎波斯盆地湖相碳酸鹽巖儲層。

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

本發(fā)明的目的在于提供一種湖相碳酸鹽巖孔隙度的計算方法，用以解決現(xiàn)有的經(jīng)驗公式法和地震反演方法都不能很好地預(yù)測湖相碳酸鹽巖孔隙度的問題。

本發(fā)明公開了一種湖相碳酸鹽巖孔隙度的計算方法，包括以下步驟：

步驟S1：選定湖相碳酸鹽巖儲層，在地層框架約束下構(gòu)建湖相碳酸鹽巖的波阻抗模型，獲取湖相碳酸鹽巖的地震反演模型；

步驟S2：通過多道空間約束地震反演方法和湖相碳酸鹽巖的地震反演模型提高湖相碳酸鹽巖地震反演波阻抗的空間表征能力，獲取湖相碳酸鹽巖的三維波阻抗參數(shù)；

步驟S3：利用BP神經(jīng)網(wǎng)絡(luò)算法實現(xiàn)湖相碳酸鹽巖三維波阻抗參數(shù)與孔隙度參數(shù)的轉(zhuǎn)化，獲取湖相碳酸鹽巖的孔隙度參數(shù)。

優(yōu)選地，所述步驟S1包括以下步驟：

步驟S11：選定湖相碳酸鹽巖儲層，獲取湖相碳酸鹽巖的全疊加地震數(shù)據(jù)、測井?dāng)?shù)據(jù)、地層速度模型和層位數(shù)據(jù)，其中測井?dāng)?shù)據(jù)包括波阻抗參數(shù)；

步驟S12：根據(jù)湖相碳酸鹽巖的地下層位的層序關(guān)系并利用層位數(shù)據(jù)構(gòu)建湖相碳酸鹽巖的地層層位框架；

步驟S13：根據(jù)測井?dāng)?shù)據(jù)開展地層速度與波阻抗參數(shù)的交會統(tǒng)計分析，獲得湖相碳酸鹽巖的地層速度與波阻抗參數(shù)的巖石物理擬合關(guān)系，將湖相碳酸鹽巖的地層速度模型轉(zhuǎn)化為初始波阻抗模型；

步驟S14：對湖相碳酸鹽巖的井信息進行地質(zhì)統(tǒng)計學(xué)分析，并據(jù)此建立湖相碳酸鹽巖的變差函數(shù)；

步驟S15：在湖相碳酸鹽巖的地層層位框架下，以所述初始波阻抗模型為協(xié)克里金約束量，結(jié)合湖相碳酸鹽巖的變差函數(shù)，利用協(xié)克里金算法更新波阻抗模型獲得湖相碳酸鹽巖的更新波阻抗模型；

步驟S16：對湖相碳酸鹽巖的更新波阻抗模型數(shù)據(jù)進行0-10Hz低通濾波處理，獲得湖相碳酸鹽巖的低頻地震反演模型，其矩陣形式記為M。

優(yōu)選地，所述步驟S2包括以下步驟：