水平井二维稳定渗流场分析及应用

玄建

中国石化胜利油田分公司勘探开发研究院，山东 东营 257000

引 言

随着钻井技术的不断进步，水平井在油气田开发中得到广泛应用，尤其是在低渗透、页岩油气藏开发中发挥了较大作用[1-4]。很多学者研究水平井时，重点研究水平井产量问题。Joshi通过垂直及水平面内的二维渗流问题，推导出均质油藏水平井稳态产能方程[5]。也有学者对水平井二维渗流进行了深入研究：李晓军等通过保角变换研究断块油藏水平井与直井组合井网渗流场[6]，张强等按照椭球体模型分析稳定渗流场[7]，袁淋等重点研究水平井及其酸化措施后的产能问题[8-10]。总体来看，目前，国内外的研究成果公式极其复杂，难以推广应用。在薄层油藏开发中，为提高开发效果，应用水平井开发是一种常用方法[11-12]，这种情况可简化为水平井二维渗流，因此，研究水平井二维稳定渗流具有重要意义。本次研究通过儒可夫斯基变换，把水平井二维稳定渗流转换为直井径向二维稳定渗流，开展水平井二维稳定渗流场中等势线、流线以及渗流速度研究，并据此开展水平段长度优化、水平段方向优化、水平井垂直注水井一注一采井网优化，并给出油气田开发过程中部署水平井的具体应用方法。

1 水平井二维稳定渗流场分析

水平井二维稳定渗流场是渗流力学分析中的重要内容，研究水平井二维稳定渗流场具有重要意义。由于水平井渗流场较复杂，因此应用保角变换分析。保角变换是一种化难为易的方法，将边界形状比较复杂的渗流场转换为较为简单的渗流场。根据儒可夫斯基公式变换

式中：w-保角变换后变量，m；

c-椭圆簇半焦距，m；

z-保角变换前变量，m。

z平面内的同心圆簇变换为w平面内以（–c，0）、（c，0）为焦点的椭圆簇，特别地，当z平面内的同心圆半径逼近c时，w平面内的椭圆逼近从（–c，0）到（c，0）的一条割线，由于水平井可视为无限导流管，因此这条割线视为二维渗流场中的水平井，那么，渗流场中的等势线为共焦点的椭圆簇

式中：x-横坐标，m；

y-纵坐标，m；

a-椭圆半长轴，m；

b-椭圆半短轴，m。

相应地，z平面内通过圆心的射线簇变换为w平面内以（–c，0）、（c，0）为焦点的双曲线簇，即水平井渗流场中的流线簇

式中：a′-双曲线半实轴，m；

b′-双曲线半虚轴，m。

根据保角变换，求得水平井产量公式

式中：Q-水平井产量，m3/s；

pe-油藏边界压力，Pa；

pc-井底压力，Pa；

µ-流体黏度，mPa·s；

K-渗透率，D；

h-油藏厚度，m；

ae-油藏边界椭圆半长轴，m；

be-油藏边界椭圆半短轴，m。

在z平面内，任一点渗流速度为

式中：vz-z平面内渗流速度，m/s；

θ-流线夹角，rad；

r-半径，m；

dlz-z平面内等势线上两点reiθ、rei(θ+dθ)的距离，m。

在w平面内，任一点渗流速度为

式中：vw-w平面内渗流速度，m/s；

dlw-dlz在w平面内的变换，m。

由渗流速度公式可以看出，当r等于1时，代表水平井井壁上的渗流速度，即流体流入水平井时，在水平井中点（θ=π/2 或θ=3π/2）渗流速度最小，愈靠近水平井两端（θ=0或θ=π）渗流速度愈大（图1）。

图1 水平井井壁不同位置上的渗流速度图Fig.1 The map of seepage velocity on the borehole wall of horizontal well

2 应用算例

2.1 水平段长度优化

某常规薄层圆形油藏面积约3.1 km2（半径1 km），具有恒定边界压力20 MPa，厚度3.0 m，单储系数15.0×104t/（km2·m），渗透率0.3 D，地下原油黏度5.0 mPa·s。在油藏中心部署一口水平井，生产压差2.0 MPa，水平井产量根据水平井二维稳定渗流公式计算得到。水平井日产油量和日产油量增长率随水平段长度变化情况见图2。从图2可见，当水平段生产井段（如：射孔完井时的射孔井段）小于100 m时，水平井日产油量增长率较大，当水平段生产井段（射孔井段）超过100 m后，水平井日产油量增长率小于0.5%。当油藏半径不同（2，4，8 km），油藏其他参数均相同时，水平井日产油量增长率与半径为1 km的油藏一致，即当水平段生产井段（射孔井段）超过100 m后，水平井日产油量增长率小于0.5%，因此该类油藏部署水平井时，建议水平段生产井段（射孔井段）不宜超过100 m。

图2 水平井日产油量和日产油量增长率随水平段长度变化图Fig.2 The change of daily output and its increasing rate with the length of horizontal section

2.2 水平段方向优化

油气田开发过程中部署水平井时通常根据数值模拟及油藏工程方法优化水平段方向[13-14]，而根据水平井二维稳定渗流场优化水平段方向也同样是一种重要的方法。根据水平井二维稳定渗流速度公式计算得到水平井二维稳定渗流场中的等势线和流线分布（图3）。渗流场中任意r处，在水平段的延伸方向（θ=0、θ=π）渗流速度最大，渗流速度可简化为

在水平段的中垂线方向（θ=π/2、θ=3π/2）渗流速度最小，渗流速度可简化为

因此，水平井二维稳定渗流场中，当r相等时，水平段延伸方向的渗流速度始终大于中垂线上的渗流速度，因此，与一口直井的二维径向稳定渗流不同，在水平井二维稳定渗流场中存在主流线，且主流线为水平段延伸方向（θ=0、θ=π）的直线，因此在部署水平井时需要考虑水平段的方向。

图3 水平井二维稳定渗流场中等势线及流线分布图Fig.3 The map of the equipotential lines and the streamlines in the two-dimensional steady seepage field of horizontal well

垦71块东营组油藏是具有气顶和较强边水能量的构造油藏，地质储量642.0×104t，含油面积4.8km2。开发初期，气窜和边水锥进现象较为普遍，为减缓开发过程中的气窜和边水锥进，计划在该油藏腰部部署水平井开发，那么气窜方向和边水锥进方向必须避开水平井渗流场中的主流线。由于构造油藏中的油气界面、油水界面平行于构造线，因此部署水平井时，水平段方向应与构造线平行，才能实现该类油藏水平井开发的最佳效果。2016年，垦71块东营组油藏在腰部部署6口水平井，方案实施后，综合含水由96.2%降至81.7%，开发效果较好。

2.3 水平井垂直注水井一注一采井网优化

水平井采油、垂直注水井注水是常见的一注一采井网，在开发方案设计中常遇到井位优化问题，很多学者的研究较为复杂[15-20]，使用不便。笔者在水平井二维稳定渗流场分析基础上，通过平移、旋转及儒可夫斯基变换，得到垂直注水井水平井形成的一注一采二维稳定渗流场流线分布图（图4）。

图4 水平井垂直注水井不同方位组合一注一采流线分布图Fig.4 The streamlines of different patterns with one horizontal well and one vertical well

图4中，垂直注水井的位置不同，则流线分布不同，但主流线均在以（–c，0）、（c，0）为焦点的双曲线上，对判断水平井来水方向和见水位置具有一定指导意义。特别地，当垂直注水井在水平井中垂线上时（θ=π/2 或θ=3π/2），主流线即为中垂线；当垂直注水井在水平井延伸方向上时（θ=0或θ=π），主流线与水平井方向一致，当垂直注水井在任一位置时（θ=θ0），主流线在双曲线上。根据推导的水平井渗流场渗流速度公式，结合水平井垂直注水井不同方位组合形成的一注一采流线分布图，得出水平井中垂线上是部署垂直注水井的最优位置。当垂直注水井部署在中垂线上时，与其它位置的垂直注水井相比，水驱波及面积相对较大，主流线渗流速度相对较小，水驱前缘推进相对较慢，能够延长水平井无水采油时间，提高经济效益。

3 结 语

（1）在水平井二维稳定渗流场分析基础上，得到水平井二维稳定渗流过程中渗流速度参数公式。特别地，在水平井井壁上，水平井中点处渗流速度最小，愈靠近水平井两端渗流速度愈大。

（2）常规油藏水平井生产井段（射孔井段）不宜超过100 m，实际应用中，应当根据油藏的具体参数优化得到最佳的水平段长度；水平段方向应于构造线平行。

（3）水平井与垂直注水井形成的一注一采二维稳定渗流场中的流线分布图，主流线在共焦点的双曲线上，对判断水平井来水方向和见水位置具有一定指导意义。分析得出垂直注水井在水平井的中垂线上时，水驱波及面积相对较大，开发效果较好，是最佳部署位置。

[1]王欢，廖新维，赵晓亮，等.超低渗透油藏分段多簇压裂水平井产能影响因素与渗流规律 以鄂尔多斯盆地长8超低渗透油藏为例[J].油气地质与采收率，2014，21（6）：107-110.doi：10.3969/j.issn.1009-9603.2014.06.027 WANG Huan,LIAO Xinwei,ZHAO Xiaoliang,et al.A study on productivity and flow regimes of segmented multi-cluster fractured horizontal well in ultra-low perme-ability reservoir:A case of Chang 8 ultra-low permeability reservoir in Ordos Basin[J].Petroleum Geology and Recovery Efficiency,2014,21(6):107–110.doi:10.3969/j.-issn.1009-9603.2014.06.027

[2]徐庆岩，于靖之，蒋文文，等.特低渗透油藏压裂水平井与直井联合开采井网产能预测新方法[J].科学技术与工程，2014，14（1）：30-34.doi：10.3969/j.issn.1671-1815.2014.18.007 XU Qingyan, YU Jingzhi, JIANG Wenwen, et al. A new method of predicting productivity of combined fractured horizontal well and vertical well for ultra-low permeability reservoir[J]. Science Technology and Engineering, 2014,14(1): 30–34. doi: 10.3969/j.issn.1671-1815.2014.18.007

[3]卞晓冰，张景臣，张士诚.低/特低渗透油藏压裂水平井布缝模式研究[J].科学技术与工程，2015，15（11）：58-62.doi：10.3969/j.issn.1671-1815.2015.11.010 BIAN Xiaobing,ZHANG Jingchen,ZHANG Shicheng.Research of fracture distribution pattern for fractured horizontal wells in low&ultra-low permeability reservoirs[J].Science Technology and Engineering,2015,15(11):58–62.doi:10.3969/j.issn.1671-1815.2015.11.-010

[4]郭小哲，周长沙.页岩气藏压裂水平井渗流数值模型的建立[J].西南石油大学学报（自然科学版），2014，36（5）：90-96.doi：10.11885/j.issn.1674-5086.2013.09.-11.07 GUO Xiaozhe,ZHOU Changsha.Seepage numerical model for fractured horizontal well in shale gas reservoir[J].Journal of Southwest Petroleum University(Science&Technology Edition),2014,36(5):90–96.doi:10.-11885/j.issn.1674-5086.2013.09.11.07

[5]邓森，王怒涛，饶庆锋，等.考虑摩阻损失的水平井长度优化[J].大庆石油地质与开发，2017，36（1）：61-64.doi：10.3969/j.issn.1000-3754.2017.01.010 DENG Sen,WANG Nutao,RAO Qingfeng,et al.Optimization of the horizontal well length considering the frictionloss[J].Petroleum Geology and Oilfield Development in Daqing,2017,36(1):61–64.doi:10.3969/j.issn.1000-3754.2017.01.010

[6]李晓军，张琪，路智勇.断块油藏水平井与直井组合井网渗流解析解及其应用[J].中国石油大学学报（自然科学版），2010，34（3）：85-88.doi：10.3969/j.issn.1673-5005.2010.03.018 LI Xiaojun,ZHANG Qi,LU Zhiyong.Analytic solution of horizontal-vertical composed well pattern seepage in fault-block reservoir and its application[J].Journal of China University of Petroleum(Edition of Natural Science),2010,34(3):85–88.doi:10.3969/j.issn.1673-5005.2010.-03.018

[7]张强，刘生国，张俊松，等.考虑流体压缩性的应力敏感油藏水平井稳定渗流模型[J].岩性油气藏，2012，24（4）：121-124.doi：10.3969/j.issn.1673-8926.-2012.04.024 ZHANG Qiang, LIU Shengguo, ZHANG Junsong, et al.A steady seepage flow model of horizontal wells based on fluid compressibility and stress sensitivity[J]. Lithologic Reservoirs, 2012, 24(4): 121–124. doi: 10.3969/j.issn.-1673-8926.2012.04.024

[8]袁淋，王朝明，李晓平，等.致密砂岩气藏气水同产水平井产能公式推导及应用[J].岩性油气藏，2016，28（3）：121-126.doi：10.3969/j.issn.1673-8926.2016.03.016 YUAN Lin,WANG Chaoming,LI Xiaoping,et al.Deduction and application of productivity formula for horizontal well with producing water in tight sandstone gas reservoir[J].Lithologic Reservoirs,2016,28(3):121–126.doi:10.3969/j.issn.1673-8926.2016.03.016

[9]袁淋，李晓平.水平井酸化后产能公式推导及应用[J].岩性油气藏，2014，26（1）：127-130.doi：10.3969/j.-issn.1673-8926.2014.01.022 YUAN Lin,LI Xiaoping.Deduction and application of productivity formula of horizontal wells after acidizing[J].Lithologic Reservoirs,2014,26(1):127–130.doi:10.-3969/j.issn.1673-8926.2014.01.022

[10]袁淋，李晓平，张璐，等.水平井稳态产能公式对比与分析[J].岩性油气藏，2013，25（6）：127-132.doi：10.-3969/j.issn.1673-8926.2013.06.024 YUAN Lin,LI Xiaoping,ZHANG Lu,et al.Analysis and correlation of steady-state productivity formulas for horizontal wells[J].Lithologic Reservoirs,2013,25(6):127–132.doi:10.3969/j.issn.1673-8926.2013.06.024

[11]李英明，任海洋.水平井技术在洼70块超稠油薄层油藏的应用[J].断块油气田，2012，19（S1）：40-43.doi：10.6056/dkyqt2012z1011 LI Yingming,REN Haiyang.Application of horizontal well in thin reservoir with super-heavy oil in Block Wa70[J].Fault-Block Oil&Gas Field,2012,19(S1):40–43.doi:10.6056/dkyqt2012z1011

[12]聂彬，刘月田.特低丰度薄层油藏水平井开发效果评价体系研究[J].油气井测试，2013，22（5）：14-17.doi：10.3969/j.issn.1004-4388.2013.05.004 NIE Bin,LIU Yuetian.Research on horizontal well development effect assessment in super-low abundance thin reservoir[J].Well Testing,2013,22(5):14–17.doi:10.-3969/j.issn.1004-4388.2013.05.004

[13]祝效华，谌柯宇，毛宇，等.水平井段脉冲携岩数值模拟研究[J].西南石油大学学报（自然科学版），2016，38（6）：164-171.doi：10.11885/j.issn.1674-5086.2015.-01.14.01 ZHU Xiaohua,CHEN Keyu,MAO Yu,et al.Numerical simulation of cuttings pulse carrying in a horizontal well[J].Journal of Southwest Petroleum University(Science&Technology Edition),2016,38(6):164–171.doi:10.11885/j.issn.1674-5086.2015.01.14.01

[14]曹丽娜，李晓平，罗诚，等.裂缝性低渗气藏水平井不稳定产量递减探讨[J].西南石油大学学报（自然科学版），2017，39（3）：103-110.doi：10.11885/j.issn.1674-5086.2016.03.21.03 CAO Lina,LI Xiaoping,LUO Cheng,et al.Investigation of unstable production decline in fractured lowpermeability gas reservoir horizontal wells[J].Journal of Southwest Petroleum University(Science&Technology Edition),2017,39(3):103–11.doi:10.11885/j.issn.1674-5086.2016.03.21.03

[15]沈瑞，李兆国，段宝江，等.王窑油田西部低渗透油藏正方形反九点井网加密调整物理模拟[J].油气地质与采收率，2014，21（3）：89-91.doi：10.13673/j.cnki.cn37-1359/te.2014.03.024 SHEN Rui, LI Zhaoguo, DUAN Baojiang, et al. Physical simulation of square inverted nine-spot pattern infilling and adjustment for low permeability reservoirs in western Wangyao Oilfield[J]. Petroleum Geology and Recovery Efficiency, 2014, 21(3): 89–91. doi: 10.13673/j.cnki.-cn37-1359/te.2014.03.024

[16]杜殿发，郑洋，赵艳武，等.反九点联合井网水平井蒸汽驱产能公式[J].特种油气藏，2016，23（4）：90-92.doi：10.3969/j.issn.1006-6535.2016.04.020 DU Dianfa,ZHENG Yang,ZHAO Yanwu,et al.Horizontal well productivity equation of steam flooding with inverted nine-spot combination well pattern[J].Special Oil&Gas Reservoirs,2016,23(4):90–92.doi:10.3969/-j.issn.1006-6535.2016.04.020

[17]赵春森，肖丹凤，宋文玲，等.水平井与直井交错井网优化方法[J].石油勘探与开发，2005，32（1）：119-122.doi：10.3969/j.issn.1008-4800.2016.07.098 ZHAO Chunsen,XIAO Danfeng,SONG Wenling,et al.Optimization of the staggered well pattern of mixed horizontal and vertical wells[J].Petroleum Exploration and Development,2005,32(1):119–122.doi:10.3969/j.issn.-1008-4800.2016.07.098

[18]蒋廷学，卞晓冰，王海涛，等.深层页岩气水平井体积压裂技术[J].天然气工业，2017，37（1）：90-96.doi：10.3787/j.issn.1000-0976.2017.01.011 JIANG Tingxue,BIAN Xiaobing,WANG Haitao,et al.Volume fracturing of deep shale gas horizontal wells[J].Natural Gas Industry,2017,37(1):90–96.doi:10.3787/j.-issn.1000-0976.2017.01.011

[19]赵金洲，陈曦宇，李勇明，等.水平井分段多簇压裂模拟分析及射孔优化[J].石油勘探与开发，2017，44（1）：117-124.doi：10.11698/ped.2017.01.14 ZHAO Jinzhou,CHEN Xiyu,LI Yongming,et al.Numerical simulation of multi-stage fracturing and optimization of perforation in a horizontal well[J].Petroleum Exploration and Development,2017,44(1):117–124.doi:10.-11698/ped.2017.01.14

[20]田鸿照，孙野.直井与水平井组合吞吐转汽驱操作参数优选[J].岩性油气藏，2013，25（3）：127-130.doi：10.3969/j.issn.1673-8926.2013.03.023 TIAN Hongzhao,SUN Ye.Operation parameters optimization of steam flooding with the combination of vertical and horizontal wells[J].Lithologic Reservoirs,2013,25(3):127–130.doi:10.3969/j.issn.1673-8926.2013.03.-023