考虑非达西渗流的底水锥进临界产量计算模型

doi:10.7623/syxb201201014

石油学报 ›› 2012, Vol. 33 ›› Issue (1): 106-111.DOI: 10.7623/syxb201201014

考虑非达西渗流的底水锥进临界产量计算模型

宋兆杰 1,2　李相方 1　李治平 2　李元生 1　尹邦堂 1

1　中国石油大学石油工程教育部重点实验室　北京　102249；2　中国地质大学能源学院　北京　100083

收稿日期:2011-05-22 修回日期:2011-09-27 出版日期:2012-01-25 发布日期:2012-03-22
通讯作者: 宋兆杰
作者简介:宋兆杰，男，1985年11月生，2010年获中国石油大学（北京）油气田开发工程专业硕士学位，现为中国地质大学（北京）在读博士研究生，主要从事石油天然气开发研究工作。
基金资助:
国家自然科学基金项目（No.50971428）、国家重大科技专项（2011ZX05030-005-03）和中央高校基本科研业务费专项资金（2011PY0206）资助。

A model for calculating critical production rates of water coning with consideration of non-Darcy flow

SONG Zhaojie 1, 2　LI Xiangfang 1　LI Zhiping 2　LI Yuansheng 1　YIN Bangtang 1

Received:2011-05-22 Revised:2011-09-27 Online:2012-01-25 Published:2012-03-22

摘要/Abstract

摘要：

底水锥进是底水气藏开发的重要问题。有关临界产量计算模型的物理模型一般假设地层流体为达西渗流，忽略非达西渗流压降。对于高产气井，近井地带流体渗流速度很大，非达西渗流阻力不可忽略。笔者在临界产量计算模型中考虑高速非达西渗流对底水锥进的影响，引入拟表皮系数，建立了考虑非达西渗流的底水锥进临界产量计算模型,计算中避免了非达西渗流半径的复杂计算过程。为量化非达西渗流的影响，应用数值模拟方法，通过建立高渗、中渗、低渗底水气藏数值模型，模拟产能试井过程并建立产能方程，计算了考虑非达西渗流的底水锥进临界产量。结果表明：低渗气藏非达西渗流对临界产量计算结果的影响可以忽略；当渗透率为30 mD、射孔厚度为208.3m、避水高度为31 m时，未考虑非达西渗流的临界产量误差为14.46%，非达西渗流影响不可忽略。因此当渗透率大于30 mD时，应选用考虑非达西渗流的底水锥进临界产量计算模型，更符合实际储层流体渗流规律。

关键词: 底水气藏, 临界产量, 底水锥进, 非达西渗流, 拟表皮系数, 数值模拟, 产能方程

Abstract:

Water coning is an important issue in the development of bottom-water-drive gas reservoirs. Physical models for calculating critical production rates generally assume that formation fluids are Darcy flow in reservoirs, which ignores the pressure drop caused by non-Darcy flow. However, the fluid flow velocity near wellbore areas in high production gas wells is so great that the pressure drop caused by non-Darcy flow should not be ignored. We investigated the influence of non-Darcy flow on water coning in models for calculating critical production rates and derived a new calculation model with full consideration of non-Darcy flow by introducing a pseudo-skin factor. This model avoids complex calculations of the non-Darcy flow radius, adopts numerical simulation methods to quantify the influence of non-Darcy flow, and calculates critical production rates with full consideration of the non-Darcy flow bottom-water coning through establishing numerical models for low-, middle-and high-permeability gas reservoirs driven by bottom water, respectively, simulating deliverability testing processes and upbuilding deliverability equations. The calculation results indicated that the influence of non-Darcy flow on critical production rates could be ignored in a low-permeability gas reservoir, while an error of the critical production rate without consideration of non-Darcy flow reached 14.46% when the permeability was 30 mD, the perforation thickness was 208.3 m and the water avoidance thickness was 31 m. Therefore, a model for calculating critical production rates of water coning with consideration of non-Darcy flow should be applied so as to be more accordant with the actual flow law of reservoir fluids when the permeability is higher than 30 mD.

Key words: gas reservoir driven by bottom water, critical production, water coning, non-Darcy flow, pseudo-skin factor, numerical simulation, deliverability equation

宋兆杰　李相方　李治平　李元生　尹邦堂. 考虑非达西渗流的底水锥进临界产量计算模型[J]. 石油学报, 2012, 33(1): 106-111.

SONG Zhaojie　LI Xiangfang　LI Zhiping　LI Yuansheng　YIN Bangtang. A model for calculating critical production rates of water coning with consideration of non-Darcy flow[J]. Editorial office of ACTA PETROLEI SINICA, 2012, 33(1): 106-111.

参考文献

[1]张勇，杜志敏，郭肖，等．底水凝析气藏薄油环开采数值模拟方法研究[J]．石油学报，2009，30(1)：88-91．

Zhang Yong,Du Zhimin,Guo Xiao,et al.Numerical simulation method for thin oil ring in condensate gas reservoir with bottom water[J].Acta Petrolei Sinica,2009,30(1):88-91.

[2]王敬，刘慧卿，刘松原，等．非均质底水油藏水平井水淹规律研究[J]．石油学报，2010，31(6)：970-973．

Wang Jing,Liu Huiqing,Liu Songyuan,et al.A flooding law in horizontal wells of heterogeneous reservoirs with bottom water[J].Acta Petrolei Sinica,2010,31(6):970-973.

[3]Meyer H I,Garder A O.Mechanics of two immiscible fluids in porous media[J].Journal of Applied Physics,1954,25(11):1400-1406.

[4]Hagoort J.Fundamental of gas reservoir engineering[M].New York:Elsevier Science Publishers B V,1988:202-205.

[5]谢兴礼，罗凯，宋文杰．凝析气新的产能方程研究[J]．石油学报，2001，22(3)：36-42．

Xie Xingli,Luo Kai,Song Wenjie.A novel equation for modeling gas condensate well deliverability[J].Acta Petrolei Sinica,2001,22(3):36-42.

[6]郑祥克，陶永建，涂彬，等．利用IPR方法确定启动压力[J]．石油学报，2003，24(2)：80-82．

Zheng Xiangke,Tao Yongjian,Tu Bin,et al.Determination of threshold pressure with IPR method[J].Acta Petrolei Sinica,2003,24(2):80-82.

[7]汪永利，蒋廷学，曾斌．气井压裂后稳态产能的计算[J]．石油学报，2003，24(4)：65-68．

Wang Yongli,Jiang Tingxue,Zeng Bin.Productivity performances of hydraulically fractured gas well[J].Acta Petrolei Sinica,2003,24(4):65-68.

[8]Verruijt A.Theory of groundwater flow[M].New York:Macmillan Press Ltd,1970:70-76.

[9]李士伦．天然气工程[M]．北京：石油工业出版社，2000：117-132．

Li Shilun.Natural gas engineering[M].Beijing:Petroleum Industry Press,2000:117-132.

[10]陈元千，董宁宇．确定气井湍流系数和湍流表皮系数的新方法[J]．断块油气田，2001，8(1)：20-23．

Chen Yuanqian,Dong Ningyu.New method of determining turbulence factor and turbulence skin factor[J].Fault-Block Oil & Gas Field,2001,8 (1):20-23.

[11]Noman R,Shrimanker N,Archer J S.Estimation of the coefficient of inertial resistance in high rate gas wells[R].SPE 14207,1982.

[12]Tessem R.High velocity coefficient dependence of rock properties[M].Norway:Thesis Pet.Inst,NTH,1980:26-30.

[13]Firoozabodi A,Katz D.An analysis of high velocity gas flow through porous media[J].JPT,1979(2):211-216.

[14]Jones S C.Using the inertial coefficient β to characterize heterogeneity in reservoir rock[R].SPE 16949,1987.

[15]Katz D L,Cornell D.Handbook of natural gas engineering[M].London:McGraw-Hill Book Company Inc.,1959:49-50.

[16]李治平，邬云龙，青永固．气藏动态分析与预测方法[M]．北京:石油工业出版社，2002：20-29．

Li Zhiping,Wu Yunlong,Qing Yonggu.The analysis and prediction methods of gas reservoir performance[M].Beijing:Petroleum Industry Press,2002:20-29.

[17]李传亮．修正DUPUIT临界产量公式[J]．石油勘探与开发，1993，20(4)：91-95．

Li Chuanliang.A modified DUPUIT formula for critical production rate[J].Petroleum Exploration and Development,1993,20(4):91-95.

考虑非达西渗流的底水锥进临界产量计算模型

A model for calculating critical production rates of water coning with consideration of non-Darcy flow

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

Metrics

本文评价

[1]	许林, 刘书杰, 许明标, 冯桓榰, 邢希金, 邓佳佳. 压差激活密封剂的微缺陷自适应修复行为及机理[J]. 石油学报, 2021, 42(5): 686-694.
[2]	邹玮, 余一欣, 刘金水, 蒋一鸣, 唐贤君, 陈石, 余浪. 东海盆地西湖凹陷中央反转构造带发育主控因素及宁波背斜形成过程[J]. 石油学报, 2021, 42(2): 176-185.
[3]	祝效华, 罗云旭, 刘伟吉, 高锐, 贾玉丹, 刘呈君. 等离子体电脉冲钻井破岩机理的电击穿实验与数值模拟方法[J]. 石油学报, 2020, 41(9): 1146-1162.
[4]	李文婧, 姜源, 单保东, 禹晓珊, 王超. 盐穴储气库注采运行时温效应对腔体稳定性的影响[J]. 石油学报, 2020, 41(6): 762-776.
[5]	窦晓峰, 宁伏龙, 李彦龙, 刘昌岭, 孙嘉鑫, 李杨, 李晓东, 赵颖杰, 张凌, 刘乐乐. 基于连续-离散介质耦合的水合物储层出砂数值模拟[J]. 石油学报, 2020, 41(5): 629-642.
[6]	计秉玉. 对油气藏工程研究方法发展趋势的几点认识[J]. 石油学报, 2020, 41(12): 1774-1778.
[7]	郭建春, 任冀川, 王世彬, 苟波, 赵俊生, 伍林. 裂缝性致密碳酸盐岩储层酸压多场耦合数值模拟与应用[J]. 石油学报, 2020, 41(10): 1219-1228.
[8]	韩强, 屈展, 叶正寅. 基于多尺度均匀化理论的页岩强度表征[J]. 石油学报, 2019, 40(7): 858-865.
[9]	苏皓, 雷征东, 李俊超, 龚斌, 赵文琪, 鞠斌山, 张荻萩, 秦鹤. 储集层多尺度裂缝高效数值模拟模型[J]. 石油学报, 2019, 40(5): 587-593,634.
[10]	解宁, 李文婧. 地下盐穴储气库盐岩热损伤机理[J]. 石油学报, 2019, 40(3): 357-369,382.
[11]	吴俊晨, 范宜仁, 曹军涛, 范卓颖, 王小龙, 吴飞. 基于大尺寸地层模型的砂岩储层油基钻井液侵入模拟[J]. 石油学报, 2019, 40(11): 1407-1414.
[12]	杨震, 马清明, 杨宁宁, 崔海波, 李运升. 基于正交天线的随钻方位电磁波电阻率成像响应特征模拟[J]. 石油学报, 2018, 39(9): 1063-1069.
[13]	徐保蕊, 蒋明虎, 赵立新, 王月文, 张晓光, 邓鑫. 螺旋分离器水流动特性的CFD模拟与PIV试验[J]. 石油学报, 2018, 39(2): 223-231.
[14]	张文彪, 段太忠, 刘彦锋, 徐睿, 杨志成, 张德民. 综合沉积正演与多点地质统计模拟碳酸盐岩台地——以巴西Jupiter油田为例[J]. 石油学报, 2017, 38(8): 925-934.
[15]	宋先知, 吕泽昊, 崔柳, 李根生, 刘昱, 胡晓东, 纪国栋. 高温射流反应腔结构设计与反应规律[J]. 石油学报, 2017, 38(5): 587-596.