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

有效

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

1.一種湖相碳酸鹽巖孔隙度的計算方法，其特征在于，包括以下步驟：

步驟S1：選定湖相碳酸鹽巖儲層，在地層框架約束下構建湖相碳酸鹽巖的波阻抗模型，獲取湖相碳酸鹽巖的地震反演模型；具體包括以下步驟：

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

步驟S12：根據湖相碳酸鹽巖的地下層位的層序關系并利用層位數據構建湖相碳酸鹽巖的地層層位框架；

步驟S13：根據測井數據開展地層速度與波阻抗參數的交會統計分析，獲得湖相碳酸鹽巖的地層速度與波阻抗參數的巖石物理擬合關系，將湖相碳酸鹽巖的地層速度模型轉化為初始波阻抗模型；

步驟S14：對湖相碳酸鹽巖的井信息進行地質統計學分析，并據此建立湖相碳酸鹽巖的變差函數；

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

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

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

步驟S3：利用BP神經網絡算法實現湖相碳酸鹽巖三維波阻抗參數與孔隙度參數的轉化，獲取湖相碳酸鹽巖的孔隙度參數。

2.如權利要求1所述的湖相碳酸鹽巖孔隙度的計算方法，其特征在于，所述步驟S2包括以下步驟：

步驟S21：根據全疊加地震數據和測井數據提取地震子波，獲得子波矩陣W，并構建地震反演的褶積算子P，其中褶積算子的公式如下：

P＝kron(I,W) (式1)

式中，kron(I,W)表示I和W的Kroneker積，I表示單位矩陣；

步驟22：通過多道空間約束的地震反演目標方程建立三維地震數據與三維波阻抗參數的映射關系，其中多道空間約束的地震反演目標方程的表達式為：

式中，S表示全疊加地震數據的矩陣；

m為待反演的地下地層界面的反射系數向量；

λ₁為低頻地震反演模型的約束系數；

M為低頻地震反演模型的矩陣形式；

K＝kron(I,C)，其中t為終止時間，t₀為起始時間，τ為時間采樣點；

λ₂為多道空間約束系數；

||m||_2,1為多道空間約束項；

步驟23：對多道空間約束的地震反演目標方程進行最小化求解，獲得m的反演結果，并計算湖相碳酸鹽巖的波阻抗數據I_P，其中波阻抗數據I_P的表達式為：

式中，t為終止時間，t₀為起始時間，τ為時間采樣點。

3.如權利要求2所述的湖相碳酸鹽巖孔隙度的計算方法，其特征在于，所述步驟S3包括以下步驟：

步驟S31：對所述三維波阻抗參數進行條件化處理，包括按照設計的采樣率對空間進行插值和利用擴散濾波處理方法除去隨機噪音；

步驟S32：對湖相碳酸鹽巖的測井孔隙度數據進行條件化處理，包括0-90Hz帶通濾波和空間采樣粗化處理，使其空間采樣率與對三維波阻抗參數進行空間插值所采用的采樣率保持一致；

步驟S33：利用BP神經網絡的算法進行波阻抗參數與孔隙度參數的非線性映射計算，將湖相碳酸鹽巖的波阻抗參數轉化為湖相碳酸鹽巖的孔隙度參數。