[發明專利]一種非均質邊水氣藏氣井見水時間預測方法及裝置

有效

• 申請號: 201911167256.0
• 申請日: 2019-11-25
• 公開（公告）號: CN112836438B
• 公開（公告）日: 2022-08-30
• 發明（設計）人: 彭鑫嶺;李繼強;付德奎;戚志林;梁梅生;黃小亮;陳建;蔣光跡
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司中原油田普光分公司
• 主分類號: G06F30/28
• 分類號: G06F30/28;G06F119/14
• 代理公司: 鄭州睿信知識產權代理有限公司 41119
• 代理人: 吳敏

1.一種非均質邊水氣藏氣井見水時間預測方法，其特征在于，該預測方法包括以下步驟：

1)獲取研究區內非均質邊水氣藏地質資料和各單井測試資料，確定氣藏儲層各層的滲透率和厚度物性參數；

2)根據距邊水氣藏氣水界面距離最近的井組的試井及動態生產測試資料，確定所述井組中各井的儲層物性參數和生產動態參數；

3)根據確定的井組中各井的儲層物性參數和生產動態參數、以及氣、水兩相滲流理論建立氣井見水時間與邊水距離的關系，確定近似均質邊水氣藏見水時間預測模型；

4)根據研究區內各井氣藏儲層中各層的滲透率和厚度物性參數確定各井儲層平均滲透率，利用最大滲透率和平均滲透率確定各井儲層滲透率突進系數；

5)根據非均質氣藏滲透率最大的小層最早見水的特性，建立滲透率最大的小層見水時間預測模型，根據近似均質邊水氣藏見水時間和非均質氣藏滲透率最大的小層的見水時間與儲層突進系數的關系，建立非均質邊水氣藏氣井見水時間預測模型，以該非均質邊水氣藏氣井見水時間預測模型進行預測。

2.根據權利要求1所述的非均質邊水氣藏氣井見水時間預測方法，其特征在于，所述步驟2)中近似均質邊水氣藏見水時間預測模型的建立過程如下：

A.根據井組各井的儲層物性參數和生產動態參數確定各井水驅前緣處氣相的壓降梯度和水相的壓降梯度；

B.根據水驅前緣處氣相和水相的壓力梯度相等的原則，建立水相滲流速度與氣相滲流速度的關系；

C.建立氣井見水時間與水驅前緣距離氣井距離在單位時間內的變化量之間的關系；

D.根據步驟C中的關系和步驟B中的關系建立氣井見水時間與邊水距離的積分關系，得到確定近似均質下邊水氣藏見水時間預測模型。

3.根據權利要求2所述的非均質邊水氣藏氣井見水時間預測方法，其特征在于，所述步驟2)中確定的近似均質邊水氣藏見水時間預測模型為：

其中t_bt為近似均質邊水氣藏見水時間；φ為儲層孔隙度；S_wc為束縛水飽和度；S_gr為殘余氣飽和度；h為儲層有效厚度；q_sc為氣井產量；B_g為氣體體積系數；μ_g為氣相粘度；μ_w為水相粘度；k為儲層絕對滲透率；k_rg(S_wc)為水驅前緣處氣相的相對滲透率；k_rw(S_gr)為水驅前緣處水相的相對滲透率；r_w為氣井井筒半徑；r_fi為氣井距水驅前緣的初始距離。

4.根據權利要求1所述的非均質邊水氣藏氣井見水時間預測方法，其特征在于，所述步驟5)中建立的非均質邊水氣藏氣井見水時間預測模型為：

其中t_btmax為非均質邊水氣藏見水時間；k_ma為儲層滲透率突進系數；φ為儲層孔隙度；S_wc為束縛水飽和度；S_gr為殘余氣飽和度；h為儲層有效厚度；q_sc為氣井產量；B_g為氣體體積系數；μ_g為氣相粘度；μ_w為水相粘度；k為儲層絕對滲透率；k_rg(S_wc)為水驅前緣處氣相的相對滲透率；k_rw(S_gr)為水驅前緣處水相的相對滲透率；r_w為氣井井筒半徑；r_fi為氣井距水驅前緣的初始距離。

5.根據權利要求1所述的非均質邊水氣藏氣井見水時間預測方法，其特征在于，所述儲層滲透率突進系數采用的計算公式為：

其中k_ma為儲層滲透率突進系數；為儲層平均滲透率；為儲層最大滲透率。