[發明專利]一種基于多尺度耦合的潛山裂縫氣藏產量預測方法

有效

• 申請號: 202210063731.5
• 申請日: 2022-01-20
• 公開（公告）號: CN114528729B
• 公開（公告）日: 2023-07-18
• 發明（設計）人: 梁豪;郭書生
• 申請（專利權）人: 中海石油（中國）有限公司海南分公司
• 主分類號: G06F30/23
• 分類號: G06F30/23;G06F30/28;G06T5/00;G06T7/10;G06T7/55;G06F113/08;G06F119/14
• 代理公司: 成都智涌知識產權代理事務所(普通合伙) 51313
• 代理人: 周正輝

1.一種基于多尺度耦合的潛山裂縫氣藏產量預測方法，包括：

(1)對于低滲透率儲集介質，通過巖心尺度進行建立數字巖心模型，獲取巖心尺度的孔隙度等效滲透率k_m和形狀因子α，并將孔隙度等效滲透率k_m和形狀因子α作為低滲透率儲集介質的基礎參數；

(2)對于高滲透率流動介質，通過測井尺度進行測井解釋，通過測井解釋資料獲取測井尺度的裂縫信息，獲取高滲流動介質的等效滲透率K_f和孔隙度并將等效滲透率K_f和孔隙度作為高滲透率流動介質的基礎參數；

(3)基于步驟(1)中獲取的低滲透率儲集介質基礎參數以及步驟(2)中獲取的高滲透率流動介質基礎參數，采用氣體竄流量q關聯低滲透率儲集介質和高滲透率流動介質的氣體，并隨機分布大尺度裂縫，采用有限元方法建立氣井生產尺度的雙重介質耦合模型，進行潛山裂縫氣藏產量預測；

所述氣體竄流量q的獲取方法為：

式中，α為低滲透率儲集介質的形狀因子；k_m為低滲透率儲集介質的滲透率，p_m為低滲透率儲集介質的壓力，p_f為高滲透率流動介質的壓力，μ為氣體粘度，ρ為氣體密度。

2.如權利要求1所述的多尺度耦合的潛山裂縫氣藏產量預測方法，所述步驟(1)中采用CT掃描方法獲取低滲透率儲集介質的基礎參數。

3.如權利要求2所述的多尺度耦合的潛山裂縫氣藏產量預測方法，所述步驟采用CT掃描方法獲取低滲透率儲集介質的基礎參數，還包括：

通過CT成像裝置完成CT圖像重構，然后依次進行圖像濾波、圖像分割處理、裂縫提取，獲取數字巖心三維圖像，通過所述數字巖心三維圖像獲取巖心尺度的孔隙度等效滲透率k_m和形狀因子α_m。

4.如權利要求3所述的多尺度耦合的潛山裂縫氣藏產量預測方法，所述濾波方法為非局部均勻濾波法。

5.如權利要求4所述的多尺度耦合的潛山裂縫氣藏產量預測方法，所屬步驟采用非局部均勻濾波法對CT圖像進行濾波，還包括：

1)以點A為中心，取其周圍邊長為s的正方形作為搜索框，再取邊長為d的正方形作為相似框，ds；

2)找出搜索框內所有邊長為d的相似框組合，記錄每個相似框的中心位置，選取其中某一個相似框的中心為點B，將點B與點A周圍的相似框做減法，同時加入由高斯核計算得到的加權值，從而獲取由各點的差值與權重的積組成的二維數組；

3)將所述二維數組求和后平均，得到相似框的中心點B相對于目標點A的擬權重值，采用指數減函數獲取得到真實的權重值；

4)結合各個點的灰度值，計算出濾波后A點的灰度；

5)對所有點進行步驟1)-4)的處理，以此類推，得到所有點經過濾波處理后的灰度值，實現整個圖像的非局部均勻濾波。

6.如權利要求1所述的多尺度耦合的潛山裂縫氣藏產量預測方法，采用裂縫系統復雜系數對裂縫系統復雜性進行表征，裂縫系統復雜系數為裂縫整體效率與裂縫平均集聚系數的幾何平均值。

7.如權利要求2所述的多尺度耦合的潛山裂縫氣藏產量預測方法，所述步驟(2)中還包括通過測井尺度REV的方法來獲取高滲透率流動介質的等效滲透率K_f。

8.如權利要求1所述的多尺度耦合的潛山裂縫氣藏產量預測方法，所述步驟(2)中高滲透率流動介質的孔隙度為測井尺度下的裂縫孔隙度減去巖心尺度下的低滲透率儲集介質的孔隙度

9.如權利要求1所述的多尺度耦合的潛山裂縫氣藏產量預測方法，所述測井尺度REV的獲取方法為：

REV邊長＝裂縫密度*裂縫長度*sinθ。