ABSTRACT

When well kicking occurs in horizontal section, it will be difficult to discover at early stage of well kick. However, the volume of overflows will suddenly increase, when the overflows passes into vertical section. Considering the characteristic of wellbore temperature field, drift model and the coupling between wellhole and formation, combining basic equation of gas-water two-phase flow, and then build a horizontal well kick theoretical model of high temperature and pressure well. This paper simulates different moment, unbalanced formation and length of barefoot interval of high temperature and pressure of horizontal wells on sea. The simulation results show that annulus pressure and void fraction vary quickly in the early stage of overflow; with the increase of formation pressure coefficient, the pit gain, annulus pressure and void fraction vary more quickly; the pit gain, annulus pressure and void fraction vary slightly in short barefoot interval; while the pit gain, annulus pressure and void fraction vary obviously in short barefoot interval.

Keywords:High Temperature and Pressure, Horizontal Wells, Multiphase Flow, Overflow Simulation

高温高压水平井溢流模拟分析

蒋凯¹，丁前²，朱海峰¹，夏强¹，孙小辉²

¹中国海洋石油国际有限公司，北京

²中国石油大学(华东)，山东 青岛

收稿日期：2016年12月2日；录用日期：2016年12月17日；发布日期：2016年12月23日

摘 要

水平井溢流后，气体在水平井段膨胀量小，溢流不容易发现，进入垂直井筒后，气体上升膨胀导致井涌快速发展，给井控造成困难。针对海上高温高压水平井的具体工况，考虑井筒温度场、漂移模型的特点以及水平井段与地层的耦合，结合气液两相流基本方程，建立了高温高压水平井井涌理论模型，并对不同时刻、不同地层欠压实、不同裸眼长度的海上高温高压水平井井涌进行了模拟。模拟结果表明：溢流早期，环空压力和截面含气率变化比较快；地层压力系数越高，泥浆池增量、环空压力和截面含气率变化越迅速；裸眼长度较小时，泥浆池增量、环空压力和截面含气率变化不明显，裸眼长度较大时，泥浆池增量、环空压力和截面含气率变化很快。

关键词 :高温高压，水平井，多相流，溢流模拟

1. 引言

溢流是水平井钻井作业过程中的常见问题，特别是高温高压水平井，当水平段发生溢流初期，气体膨胀量很小，压力保持基本恒定，关井套压等于或接近于关井立压，泥浆池增量很小，井涌监测困难，突发性强 [1] ；加上地层高温高压，若发生溢流，井控难度极大 [2] 。为了对高温高压水平井溢流进行有效控制，必须对溢流进行模拟分析。国内外学者对水平井的稳态温度和压力 [3] [4] ，以及高温高压井控问题做了大量工作，形成了比较完善了计算模型，但是对高温高压水平井还缺少相应的研究。由于高温高压水平井钻井过程中复杂温度场以及气侵发生后水平段多相流压力场的复杂性，致使无法准确地对高温高压水平井进行模拟分析。本文通过建立高温高压水平井溢流模拟的温度场和压力场计算理论模型，借助计算机算法，研究高温高压水平井钻井过程中井涌动态规律，从而对高温高压水平井溢流进行准确地模拟。

2. 高温高压水平井井涌理论模型

2.1. 气液两相流压力计算通用控制方程

压力是进行溢流模拟的关键参数，对高温高压水平井溢流进行模拟就必须对其压力场作准确计算。溢流一旦发生，整个井筒内就变成气、液两相环境，为了准确地对压力场进行计算，必须考虑各相流体的影响，得到气液两相流动的控制方程组，然后对其求解得到井筒中的压力分布 [5] 。计算气液两相流压力时，控制方程组主要由连续性方程、动量守恒方程组成 [6] ：

由连续性方程和动量守恒方程可以看出，其中求解的关键参数有气相速度和液相速度，漂移模型是计算气相速度和液相速度的常用模型，该最早由N. Zuber和J. A. Findlay提出 [7] ，后来学者在此基础上进行了修正，提出了水平井筒气–液两相流漂移模型 [8] ；而水平井由垂直段、造斜段及水平段三部分组成，气相趋于向上井壁聚集，其滑脱速度受井眼轨迹的限制，使上述漂移模型在计算全井段时结果出现较大偏差。在调研大量漂移模型的基础上，发现H. Shi和J. A. Holmes提出的气液两相流模型在井斜范围为0˚~88˚内能够很好的与实验结果相拟合 [9] ，基本满足了工程需要。

计算井筒压力的另一个关键因素是溢流量，水平段各位置处溢流速率会影响水平井筒内的压力分布，而井筒内压力分布反过来影响从气藏的径向流入量的大小及分布，溢流气体的侵入过程和井筒内的多相流动过程相互影响，需要对水平井段与地层进行耦合以便对溢流过程进行准确描述。

由于钻井过程中循环压耗的存在，从跟端(靠近垂直井筒段的一端)到指端(远离垂直井筒段的一端)，环空压力不断减小，地层溢流量不断增大，使水平井段的溢流量产生不均匀变化，因此应对水平段进行积分求出溢流量 [10] 。

(1)

2.2. 水平井温度场方程

传统溢流模拟过程中通常忽略温度场的影响，主要关注压力的变化。但是，溢流发展过程中，气体的物理性质，如膨胀性，与环空温度分布密切相关。考虑深水环境及水平井钻进工况的特点，需要对溢流期间的井筒温度场进行准确描述。

传统上，对于石油工程井筒内流动的温度场分布，已经进行了广泛的研究，主要集中于单相流状态下的计算，如正常钻进工况 [11] [12] [13] ，生产工况 [14] [15] [16] ，气井测试工况 [3] ，固井工况 [17] 。但溢流过程中，井筒内存在气液两相流动，而温度分布对于流体的物性及流动特征影响很大，需要进行准确描述。

考虑酸性气体的焓变，考虑复杂的环境温度场，结合井筒内多相流动特点，推导了针对于海上高温高压水平井工况下的井筒温度场计算模型(如图1所示)。

考虑溢流气体在井筒环空流动过程中的相态变化，利用比焓进行水平井井筒温度场模型推导。基本假设为：考虑钻杆和环空内的流动为一维流动；不考虑钻头摩擦生热，忽略纵向热量传递；不考虑岩屑的影响；考虑钻杆内为液态单相流，环空内为气液两相流。

钻杆内钻井液以及环空内溢流气体和钻井液混合物稳态流动下，流体的能量守恒方程为：

(2)

式中，Q为单位质量的流体向外传递的热量，J/kg；h为单位质量流体的焓，J/kg；为流体流速，m/s；g为重力加速度，m/s²；为井斜角，rad；L为井筒任意深度，m。

由钻杆流入环空的热量为

(3)

其中

式中，和分别为钻杆和环空内流体的温度，K；为钻杆内流体的比热容，；为钻杆内钻井液的质量流速，kg/s。

钻杆内为钻井液单相，环空内为钻井液和溢流气体的气液两相，则

(4)

(5)

式中，为截面含气率；为截面持液率；为定压比热容比热容，；为焦耳–汤普逊系数，K/MPa；为环空内的压力，MPa。

图1. 溢流期间井筒传热示意图

由环空流入钻杆和地层的热量和为

(6)

其中

式中，为外界环境温度，K；为环空内流体的比热容，；为环空内钻井液的质量流速，Kg/s；为环空内溢流气体的质量流速，Kg/s；为无因次温度；为地层导热系数，；为海水温度，K；为水平井跟端地层温度，K；为地温梯度，℃/100 m；H为井深，m；为海水深度，m；为水平井跟端井深，m。

将式(2)、(3)、(4)联立，得到环空内的流体温度场方程为

(7)

将式(2)、(5)、(6)联立，得到钻杆内的流体温度场方程为

(8)

其中

式中，和分别为钻杆和环空流体的传热系数，；和分别为钻杆内径和套管内径，m；为钻杆内钻井液流速，m/s；为环空内溢流气体的表观流速，m/s；为环空内钻井液的表观流速，m/s，为环空内流体的混合平均流速，m/s。

2.3. 理论模型的离散及求解

对多相流控制方程以及温度场方程直接求解析解是很困难的，而且深水钻井多温度梯度的环境要求对方程进行离散求解，因此上述方程需采用数值算法求其离散解具体求解过程如图2。

3. 高温高压水平井井涌动态实例分析

针对某高温高压气田，依据本文提出的求解方法对海上高温高压水平井井涌模拟分析模型进行求解。该井参数为：该气田属于高温高压体系，最大油藏压力57.3 MPa，最高油藏温度151.7℃。目标井深4000 m，造斜点位置2800 m，在3700 m进入水平段；地层破裂压力梯度0.0215 MPa/m，井控设备承压能力13.6 MPa，海水密度1025 Kg/m³，钻井液密度1797 Kg/m³，塑性粘度0.03 Pa.s，打开油层厚度500 m，气相渗透率30 md，天然气相对密度0.917，钻头直径149.5 mm，钻井液排量0.016 m³/s，水深55.5 m，隔水管内径508 mm，套管内径164 m，地热梯度为3.5℃/100 m；该井井深结构如图3所示。应用基础数据，根据前面的求解方法对高温高压水平井溢流模型进行求解，基于计算结果进行了如下讨论。

图2. 求解过程

图3. 计算井井身结构

3.1. 不同时刻井涌模拟分析

图4井涌10 min时刻，不同时刻环空压力随井深变化的对比。从图4可以看出，井深小于2000 m时，不同时刻环空压力基本不变化；井深大于2000 m时，随井深增加，不同时刻环空压力差增大。对比0 min、5 min、10 min时刻环空压力可以看出，溢流早期环空压力变化比较快。由于溢流早期溢流量小，气体以气泡形态分布于钻井液内，因而其运移速度小，气体在水平段聚集，故环空压力变化比较大；随着溢流时间增加，气泡开始聚并，在水平井环空上端面形成连续相，由于气体滑脱速度的存在，气相运移速度增大，一部分气体运移出水平段，故环空压力变化相对减弱。因此水平井钻井过程中若监测到环空压力波动，而后波动很小，不应放松警惕，以免造成井涌甚至井喷。

图5为井涌10 min时刻，不同时刻截面含气率随井深变化的对比。从图5可以看出，井深小于2000 m时，由于溢流还未向上运移到该井段，不同时刻截面含气率相同；井深大于2000 m时，随井深增加，不同时刻截面含气率先增大后减小，在井深3600 m左右达到最大值，即在水平井跟端达到最大值，因此当溢流气体运移出水平段，进入倾斜垂直井段后，由于环空压力的降低，气体会不断膨胀，加上地层高压，会导致井涌突然失控造成井喷。对比0 min、5 min、10 min时刻截面含气率可以看出，溢流早期截面含气率变化比较快。由于溢流早期溢流量小，气体在钻井液内以气泡形态存在，因而其运移速度小，

图4. 不同时刻环空压力随井深变化曲线

图5. 不同时刻截面含气率随井深变化曲线

气体在水平段聚集，故环空截面含气率比较大；随着溢流时间增加，气泡开始聚并，在水平井环空上端面形成连续相，由于气体滑脱速度的存在，气相运移速度增大，一部分气体运移出水平段，造成环空含气量相对减小，故环空截面含气率相对降低。

3.2. 不同地层压力条件井涌模拟分析

图6为井涌发生后，不同地层压力条件泥浆池增量随时间变化的对比。从图6可以看出，泥浆池增量随溢流时间的增加而增大；对比地层压力系数分别为2.0、2.1、2.2时泥浆池增量可以看出，地层压力系数越大，泥浆池增量越大；地层压力越高，进入井筒的气体越多，导致环空压力越低，进而又促进气体进入井筒，因此，高压地层易造成泥浆池增量的迅速增大，引发严重的井喷事故。

图7为井涌10 min时刻，不同地层压力条件环空压力随井深变化的对比。从图7可以看出，井深小于3600 m时，不同地层压力条件环空压力基本不变化；井深大于3600 m时，随井深增加，不同地层压力条件环空压力差增大。对比地层压力系数分别为2.0、2.1、2.2时环空压力可以看出，地层压力系数越

图6. 不同地层压力下泥浆池增量随时间变化曲线

图7. 不同地层压力下环空压力随井深变化曲线

高，环空压力变化越明显，故高压水平井溢流时，环空压力会出现比较明显的变化，若控制不及时，超过设备的承压极限，会造成井涌失控引发井喷。

图8为井涌10 min时刻，不同地层压力条件截面含气率随井深变化的对比。从图8可以看出，井深小于2000 m时，由于溢流还未向上运移到该井段，不同地层压力条件截面含气率相同；井深大于2000 m时，随井深增加，不同地层压力条件截面含气率先增大后减小，在井深3600 m左右达到最大值，即水平井跟端达到最大值。对比地层压力系数分别为2.0、2.1、2.2时截面含气率可以看出，地层压力越高，截面含气率变化越迅速，井控反应时间越短。

3.3. 不同裸眼段长度井涌模拟分析

图9为井涌发生后，不同裸眼段长度泥浆池增量随时间变化的对比。从图9可以看出，溢流早期泥

图8. 不同地层压力下截面含气率随井深变化曲线

图9. 不同裸眼段长度泥浆池增量随时间变化曲线

浆池增量很小，溢流后期泥浆池增量增加比较明显。对比裸眼段长度分别为100 m、300 m、500 m、700 m时泥浆池增量可以看出，裸眼段较小时，泥浆池增量变化不显著；裸眼段较大时，泥浆池增量显著增加。因此，裸眼段较小时，溢流很难监测；裸眼段较大时，溢流后期泥浆池增量能够明显的监测到溢流。

图10为井涌10 min时刻，不同裸眼段长度环空压力随井深变化的对比。从图10可以看出，井深小于3000 m时，不同裸眼段长度环空压力基本不变化；井深大于3000 m时，随井深增加，不同裸眼段长度环空压力差增大。对比裸眼段长度分别为100 m、300 m、500 m、700 m时环空压力可以看出，随着裸眼长度的增加，环空压力变化越明显。因此在长水平段钻井过程中，较长的裸眼长度容易造成较多的溢流涌入井筒，需要严密监测环空压力变化，及早做出预防措施，确保安全钻进。

图11为井涌10 min时刻，不同裸眼段长度截面含气率随井深变化的对比。从图11可以看出，井深

图10. 不同裸眼段长度环空压力随井深变化曲线

图11. 不同裸眼段长度截面含气率随井深变化曲线

小于2000 m时，截面含气率相同；井深大于2000 m时，随井深增加，不同裸眼段长度截面含气率先增大后减小。对比裸眼段长度分别为100 m、300 m、500 m、700 m时截面含气率可以看出，截面含气率随着裸眼长度的增加而增加，裸眼长度越长，截面含气率峰值对应的井深越浅，即随着裸眼长度的增加，溢流峰值能够更快地沿井眼环空向上运移。因此，在长水平段钻井过程中，截面含气率较高，溢流峰值运移速度较快，井涌突发性强，井控反应时间短，需要做好周密的预防措施。

4. 结论

(1) 溢流早期环空压力变化相对较快，截面含气率变化也相对较快。

(2) 高压水平井钻井过程中，井涌突发性强，井控反应时间短，钻井作业过程中应严密监测溢流。

(3) 在长水平段钻井过程中需要严密监测泥浆池增量、环空压力以及截面含气率的变化，及早做出预防措施，确保安全钻进。

(4) 溢流进入垂直管段后会因井筒压力降低而膨胀，地层高温会加速井涌发生的突发性；高温高压水平井由于地层高温、储层高压、水平钻井，一直是钻井作业的难题，本文通过对海上某高温高压气井进行溢流模拟，对高温高压水平井的安全顺利钻进提供了一定技术借鉴。

基金项目

国家重大专项(2016ZX05020006)，国家自然基金(51374232)资助。

文章引用

蒋凯,丁前,朱海峰,夏强,孙小辉. 高温高压水平井溢流模拟分析
Overflow Simulation Analysis of Horizontal Well with High Temperature and Pressure[J]. 流体动力学, 2016, 04(04): 82-92. http://dx.doi.org/10.12677/IJFD.2016.44010

参考文献 (References)