海洋双层管双梯度钻井井筒温度分布规律研究

王江帅 付盼 胡旭辉 宫臣兴 邓嵩 唐政 殷文

摘要：井筒温度对天然气水合物的存在状态有显著影响。为了准确掌握海洋双层管双梯度钻井井筒温度分布规律，基于传热学理论，充分考虑岩屑对内管混合流体热物性参数的影响，以及隔水管、套管、裸眼等井段环空内不同流体类型对传热的影响，引入综合导热系数，建立了适用于双层管双梯度钻井的井筒温度场预测模型，并探究了井筒温度分布规律。研究结果发现：对双层管双梯度钻井井筒而言，其在海底附近的温度达到最低值，在水合物固态流化开采时应特别注意该井段的井筒温度分布；注入流体温度、排量和钻井液比热容的增大，会导致内管流体温度剖面整体增大；注入流体温度主要影响井口附近的内管流体温度，而排量和钻井液比热容主要影响海底附近的内管流体温度，可以通过调节上述参数来改善不同井段井简温度分布，以减少水合物的生成。研究结果可为水合物的安全高效开采提供理论支撑。

关键词：天然气水合物；双层管双梯度钻井；井筒温度；综合导热系数；内管；外管环空

中图分类号：TE52 文献标识码：A DOI：10.16082／j．enki．issn．1001-4578.2022.12.008

Distribution Law of Wellbore Temperature in Offshore Dual-layer DEG

Wang Jiangshuai' Fu Pan Hu Xuhui Gong Chenxing5 Deng Song' Tang Zheng' Yin Wen'

（1.School of Petroleum and Natural Gas Engineering，Changzhou University; 2. CNPC Engineering Technology R&D Company Limited;3.KEMBL Petroleum Technology Derelopment Co.，Lid.;4.Oil and Gas Technology Research Institute，PetroChina Chan- gqing Oilfield Company; 5. State Engineering Laboratory of Exploration and Derelopmentof Lou-Permeability Oil & Gas Field）

Abstract：Wellbore temperature has a significant influence on the state of natural gas hydrate. In this paper， a wellbore temperature field prediction model for dual-layer dual-gradient drilling （DGD） was built according to the theory of heat transfer by incorporating the comprehensive thermal conductivity and considering the influence of cut- tings on the thermal physical parameters of the mixed fluid in the inner pipe and the influence of different fluids in the annulus of riser， casing and open hole on the heat transfer. With this model，the distribution law of wellbore temperature in offshore dual-layer DGD was investigated. It is found that， for the dual-layer DGD wellbore，the temperature near the seabed is the lowest，which should be carefully considered in the process of hydrate solid flu- idization production.The increase of injected fluid temperature，displacement and specific heat capacity of drilling fluid lifts the temperature profile of inner pipe fluid. Specifically， the temperature of the injectedI fluid mainly af- fects the temperature of the inner pipe fluid near the wellhead，while the displacement and the specific heat capaci- ty of drilling fluid mainly affect the temperature of the inner pipe fluid near the seabed.Therefore，it is possible to improve the wellbore temperature distribution in different well sections by adjusting the above parameters to mini-mize the formation of hydrate.'The research provide theoretical support for the safe and efficient exploitation hy-drate.

Keywords： natural gas hydrate;dual-layer dual-gradient drilling; wellbore temperature; comprehensivethermal conductivity ; inner pipe; outer pipe annulus

0引言

为了有效解决深水钻井的窄密度窗口[1—6]和浅层水合物开采面临的地层漏失严重、开采效率低[7]等难题，西南石油大学能源装备研究院创新性提出了一种基于双层管双梯度的深海油气及水合物开发技术方案[8]，并建立了该钻井方式下的井筒压力调控理论及策略。然而，该钻井方式下井筒流体物性参数变化、井筒压力计算[9—10]以及后续开采过程中水合物在双层管内是否发生分解、分解位置以及分解量等均受井筒温度的显著影响。因此，亟待开展海洋双层管双梯度钻井井筒温度场研究，以期为井筒压力的准确预测和保障水合物安全、高效的开采提供理论支撑。

现有的传热模型大多是用来预测正常循环[10—11]、漏失[12]、溢流[3]等不同工况下的井筒温度分布，这些模型针对的是“钻杆内流入、环空流出”这种井筒流动传热方式，适用于常规钻井方式。然而，双层管双梯度钻井过程中钻井液由外管环空注入，经钻头后与岩屑一起由内管返出。此外隔水管环空、套管环空、裸眼环空分别被海水、隔离液、钻井液 3 种不流动的流体充满[8]。由此可见，双层管双梯度钻井的井筒流动传热方式为“环空流入、钻杆内流出”，且不同井段环空流体类型差异较大，导致井筒—地层之间的传热非常复杂。现有的井筒温度场预测模型已不再适用。为此，笔者基于传热学理论，充分考虑岩屑对内管混合流体热物性参数的影响，以及隔水管、套管、裸眼等井段环空内不同流体类型对传热的影响，并将外管内壁至原始地温位置之间的区域，等效为一个传热单元，引入综合导热系数，建立了适用于双层管双梯度钻井的井筒温度场预测模型，并对井筒温度分布规律进行分析。

1物理模型

双层管双梯度钻井井筒物理模型如图 1 所示。由图 1 可知，钻井过程中钻井液由外管环空注入，经钻头喷嘴后与岩屑一起由内管返出，环空被海水、隔离液、钻井液3种流体充满，并处于静止状态。在此过程中，地层与外管环空钻井液进行热交换，外管环空钻井液与内管混合流体进行热交换。

2数学模型

2.1 基本假设

在建立双层管井筒温度场数学模型时，进行以下假设：①内管和外管环空内流体温度在径向上不发生变化；②井筒中心线向外一定距离处地层温度为原始地温（C.S.HOLMES and S.C.SWIFT 认为 不受扰动的地层半径值为3.05 m[14]，本文选取该值作为原始地温位置与井眼中心轴线的距离），同时将外管内壁至原始地温位置之间的区域等效为一个传热单元，引入综合导热系数；③热源项计算时，考虑由于井筒流体流动摩阻所产生的热量；④综合导热系数计算时，考虑海水、隔离液、钻井液等不同流体导热系数差异。

2.2 温度场模型

与常规钻井方式不同，双层管双梯度钻井流体从外管环空注入，经钻头后从内管返出。因此，基于热力学第一定律，分别针对外管环空和内管（见图2） 2种流动传热区域内的流体建立热传导控制方程。

2.2.1外管环空流体

其中：

式中：Q为流入外管环空控制体的流体热量，J；Qa，ou为流出外管环空控制体的流体热量， J； 9a.p为单位时间内外管环空流体向内管流体的导热量，J；9..，为单位时间内外管环空流体向地层的导热量，J；S为由于流体流动时压耗所产生的热量，J； Qachang为外管环空控制体流体热量的改变量，J；C.为流入外管环空控制体流体的比热容，J/（kg.℃）； qm为流入外管环空控制体的流体质量流量， kg/s； T.为流入外管环空控制体的流体温度，℃； T，.L+AL为流出外管环空控制体的流体温度，℃； D，为内管外径， m； A，为双层管的导热系数， W/（m.℃）； D，为原始地层温度位置处的圆柱外径，m；入l为外管内壁至原始地层温度位置处传热体的综合导热系数， W/（m.℃）； Twell为外管内壁至原始地层温度位置处的传热体温度，℃； A，为外管环空流体流动的面积， m2； p，为流入外管环空控制体的流体密度，kg/m'； T，，为控制体流入端对应的内管流体温度，℃；b，为内管的壁厚，m； Te为控制体流入端对应的原始地层温度，℃；b.l为外管内壁至原始地层温度位置处的距离，m； △L为控制体的长度，m； △t 为时间步长，s。

代入控制方程，推导后得出外管环空流体温度计算方程：

其中：

2.2.2 内管流体

其中：

（5）

式中： Qpin为流入内管控制体的流体热量，J；Qp.ou为流出内管控制体的流体热量， J； Qp.change为内管控制体流体热量的改变量， J； C，为流入内管控制体流体的比热容， J/（kg.℃）； qm为流入内管控制体的流体质量流量，kg/s； Tp.+s为流入内管控制体的流体温度，℃；A，为内管流体流动的面积， m2； p，为流入内管控制体的流体密度， kg/m2。

代入控制方程，推导后得出内管流体温度计算方程：

其中：

2.3 辅助方程

2. 3. 1 热物性参数

由于钻井液从井底往上返出时有岩屑的进入，因此改变了内管混合流体的热物性参数。

混合流体密度：

混合流体的质量流量：

混合流体的比热容：

式中：p.为岩屑的密度， kg/m'； Q，为钻井液排量，m2/s；Q. 为岩屑的产生速率，m2/s； C. 为岩屑的比热容，J/（kg·℃）；m，为混合流体质量流量， kg/s； m，为返出流体质量流量， kg/s； m，为岩屑的质量流量，kg/so

2. 3. 2 综合导热系数

由于建立数学模型时，将外管内壁与原始地温之间的区域视为一个传热单元，因此针对该传热单元引入综合导热系数入ll，计算公式如下：

式中：L为计算点的深度， m； Dwater为水深， m；Da为套管鞋深度， m; bai为外管壁厚， m; Asea 为海水导热系数，W/（m·℃）；b2为隔水管段环空间隙， m； Arr为隔水管导热系数， W/（m.℃）; bri为隔水管壁厚，m; 入为隔离液导热 系数， W/（m·℃）； Aeaing为套管导热系数，W/（m·℃）; beaing为套管壁厚， m; An为水泥 环导热系数，W/（m·℃）；beeme为水泥环壁厚，m; 入;为地层导热系数， W/（m.℃）;A。为钻井 液导热系数，W/（m·℃）。

2.3.3 热源项计算

热源项[1?]由下式计算：

式中：△P，为单位长度的摩阻压降，Pa/m； A，为钻头喷嘴的截面积，m2； E，为钻头破岩效率；M为轴向上的钻柱扭矩， N.m； w为转速， r/s

3模型求解

通过迭代方式对井筒温度场模型进行求解，具体求解步骤如下：①假定内管流体温度为原始地温7.=T.L，利用公式（3）计算出外管环空流体温度分布T； ②将外管环空流体温度T，代入公式（7）计算内管流体温度分布T；③将计算得到的内管流体温度T，，再次代入公式（3）重新计算外管环空流体温度分布T；④将外管环空流体温度71代入公式（7）中重新计算内管流体温度分布 7p.。重复步骤③和④，当内管流体温度分布满足以下判断依据时，认为计算结果稳定。

式中： D为井深， m； i为时间节点；j为空间节点；E 为误差控制参数，无量纲。

4实例分析

使用上述建立的数学模型，开展了双层管双梯度钻井井筒温度场数值模拟研究。模拟井为一口深水直井，基本参数包括：井深 2000 m，水深 1500m，套管下深1 800 m，隔水管外径0.533 m，隔水管内径 0. 508 m，套管外径 0.468 m，套管内径0.404 m，外管外径 0.308 m，外管内径 0.248 m，内管外径 0. 127 m，内管内径 0. 109 m，钻井液密度1 100 kg/m2，排量 30 1/s，钻井液比热容3 180J/（kg·℃），岩石密度 2 650 kg/m3，岩石比热容920J/（kg·℃），管柱导热系数43 W/（m·℃）， 岩石导热系数 2.25 W/（m·℃），水泥环导热系数1.7 W/（m .℃），钻井液导热系数1 W/（m·℃），海水导热系数 0.8 W/（m·℃），隔离液密度1 440 kg/m2，隔离液导热系数 0.5W/（m·℃），隔离液比热容 1 840J/（kg·℃），机械钻速 2 m/h，注入流体温度 30 ℃，海平面温度30℃，地温梯度0.028 ℃/m。

4.1 井筒温度分布规律

图3展示了双层管双梯度钻井时外管环空流体温度和内管流体温度分布情况。由图 3 可以看出：循环时钻井液以一定温度从外管环空注入，向下流动过程中受到海水温度场的影响，外管环空流体温度逐渐降低；进入泥线以下地层时，受地层温度场的影响，外管环空流体温度有所升高。总体来说，在井筒上部（即图3中①区域），热传导方向为外管环空流体→地层；在井筒下部（即图3中②区域），热传导方向为地层→外管环空流体。此外，外管环空流体经钻头后沿内管向上返出，该过程中内管流体受周围环境温度的影响在流动方向上呈现出先减小后增大的变化趋势。

另外，结合双层管双梯度钻井进行水合物固态流化开采工艺，对井筒温度分布进行深入分析。由于受到海水温度场的影响，水合物固态流化开采时，在井下被破碎并分离后，沿内管向上流动过程中在海底泥线附近区域（低温高压环境）易形成结晶。一旦有水合物晶体形成，将会影响井筒流动，严重时会给井筒流动保障带来安全隐患。因此，建议优化相关可控参数来增大海底泥线附近井筒温度，避免因井筒温度过低而形成大量的水合物晶体。

4.2 不同因素对井筒温度分布的影响

4. 2.1 注入流体温度

图 4 展示了 3 种不同注入流体温度条件下的内管流体温度分布情况。由图 4 可以看出，随着注入流体温度的升高，内管流体温度剖面整体增大。内管流体温度的升高有利于抑制水合物的生成，从而缩小水合物在井筒内的生成区域，保障井筒流动。

由图 4 还可以看出，井口附近内管流体温度受注入流体温度的影响较大，而海底、井底附近内管流体温度受注入流体温度的影响相对较小。具体来看，注入流体温度从 20 ℃ 升高到 40 ℃时，井口（0 m）温度升高了18.25 ℃，而海底（1 500 m）温度仅升高了 1.48 ℃，井底（2000 m）温度仅仅增加了 1.02 ℃。这是因为，对于井口流体来说，内管流体温度受外管环空注入流体温度直接影响，然而对于井底来说，经过长时间循环后，井筒换热充分，因此注入流体温度对井底温度的影响非常有限。

4. 2.2 排量

图5展示了3 种不同排量条件下的内管流体温度分布情况。由图 5 可以看出，随着排量的增加，内管流体温度剖面整体增大，有利于抑制水合物的生成，从而缩小水合物在井筒内的生成区域，保障井筒流动。与常规陆地钻井不同的是，对于深水浅部地层钻井，排量的增大并没有导致井底温度减小，反而使其增大。这是因为，排量越大井筒内流体与地层的换热时间越少，受深水海底和浅部地层低温的影响更小，故内管流体温度增加。

由图 5 还可以看出，海底附近内管流体温度受排量的影响较大，而井口、井底附近内管流体温度受排量的影响相对较小。具体来看，排量从 20 L/s增加到40 I/s 时，井深1 500 m 处的海底温度升高了3.64 ℃，井底（2000 m）温度升高了1.18 ℃，而井口（0 m）温度仅仅变化了 0.27 ℃。也就是说，相比于井口、井底，排量的增加可以更有效地改善海底附近井筒内的低温环境，从而有效抑制水合物在海底附近的井筒内形成。

4. 2. 3 钻井液比热容

图6展示了3种不同比热容钻井液在内管流体温度分布情况。由图 6 可以看出，随着钻井液比热容的增加，内管流体温度剖面整体增大，同样有利于抑制水合物的生成，从而缩小水合物在井筒内的生成区域，保障井筒流动性。这是因为，钻井液比热容越大，井筒内流体与流体之间的热对流增强越明显，相比于海底低温，注入流体温度较高，有利于内管流体温度增加。

由图6还可以看出，海底附近内管流体温度受钻井液比热容的影响较大，而井口、井底附近内管流体温度受钻井液比热容的影响相对较小。具体来看，钻井液比热容从 2 180J/（kg·℃）增加到4 180J/（kg·℃）时，井深 1 500 m 处的海底温 度升高了3.45 ℃，井底（2 000 m）温度升高了1.13 ℃，而井口（0m）温度仅仅变化了0.25 ℃。也就是说，相比于井口、井底，钻井液比热容的增加可以明显地改善海底附近井筒内的低温环境，对有效抑制水合物在海底附近的井筒内形成有一定的作用。

综上分析，可以得出以下 3 点认识：①对双层管双梯度钻井井筒而言，其在海底附近的温度达到最低值，在水合物固态流化开采时应注意增大海底附近井筒温度，避免因温度过低而形成大量的水合物晶体，从而阻碍井筒流动。②注入流体温度、排量和钻井液比热容的增大，均会使得内管流体温度剖面整体增大，是通过调节井筒温度剖面来抑制水合物生成的重要手段。③注入流体温度主要影响的是井口附近的内管流体温度，而排量和钻井液比热容主要影响的是海底附近的内管流体温度，可以通过调节上述参数来改善不同井段井筒温度分布，以减少水合物的生成。

5结论

（1）针对双层管双梯度钻井方式，充分考虑岩屑对内管混合流体热物性参数的影响以及隔水管、套管、裸眼等井段环空内不同流体类型对传热的影响，并将外管内壁至原始地温位置之间的区域等效为一个传热单元，引入综合导热系数，建立了该钻井方式下的井筒温度场预测模型。

（2）对双层管双梯度钻井井筒而言，其在海底附近的温度达到最低值，在水合物固态流化开采时应特别注意该段的井筒温度分布。

（3）注入流体温度、排量和钻井液比热容的增大，均会使得内管流体温度剖面整体增大；其中，注入流体温度主要影响的是井口附近的内管流体温度，而排量和钻井液比热容主要影响的是海底附近的内管流体温度，可以通过调节上述参数来改善不同井段井筒温度分布，以减少水合物的生成。

参 考 文 献

[1]孙宝江，张振楠.南海深水钻井完井主要挑战与对策[J]. 石油钻探技术， 2015， 43 （4）： 1—7.

SUN B J， ZHANG Z N. Challenges and countermeas- ures for the drilling and completion of deepwater wells in the South China Sea [J]. Petroleum Drilling Tech- niques， 2015， 43 （4）： 1-7.

[2]赵涛，黄元元，贾向锋，等.我国海洋油气钻井装备技术现状及展望[J].石油机械， 2022， 50 （4）：56—62.

ZHAO T， HUANG Y Y， JIA X F， et al. Offshore oil/gas drilling equipment in China： review and pros- pect [J]. China Petroleum Machinery， 2022， 50 （4）： 56-62.

[3] 王江帅，李军，柳贡慧，等.基于井下分离的深水双梯度钻井参数优化[J].石油勘探与开发，201946 （4）： 776-781.

WANG JS， LIJ， LIU GH， et al. Parameters optimi- zation in deepwater dual-gradient drilling based on downhole separation [J]. Petroleum Exploration and Development， 2019， 46 （4）： 776-781.

[4]张锐尧，李军，柳贡慧，等.深水钻井多压力系统条件下的井筒温度场研究[J].石油机械， 2021，49 （7）： 77-85.

ZHANG RY， LIJ， LIU GH， et al. Research on the wellbore temperature field under the multiple pressure system during deep water drilling [J]. China Petroleum Machinery， 2021， 49 （7）： 77-85.

[5]罗鸣，吴江，陈浩东，等.南海西部窄安全密度窗口超高温高压钻井技术[J].石油钻探技术，201947 （1）： 8-12.

LUO M， WU J， CHEN H D， et al. Ultra-High tem- perature high pressure drilling technology for narrow safety density window strata in the westen South China [J]. Petroleum Drilling Techniques， 2019， 47 （1）： 8-12.

[6]王江帅，李军，何岩峰，等.变梯度控压钻井井控过程中井口回压变化规律[J]. 石油学报，2021，42 （11）： 1499-1505.

WANG J S， LI J， HE Y F， et al. Variation law of wellhead back pressure under well control during varia-ble gradient managed pressure drilling [J]. Acta Petro- lei Sinica， 2021， 42 （11）： 1499-1505.

[7]付强，王国荣，周守为，等.海洋天然气水合物开采技术与装备发展研究[J]. 中国工程科学，2020，22 （6）： 32-39.

FU Q， WANG G R， ZHOU S W， et al. Development of marine natural gas hydrate mining technology and e- quipment [J]. Strategie Study of CAE， 2020， 22 （6）： 32-39.

[8]王国荣，钟林，刘清友，等.基于双层管双梯度深海油气及水合物开发技术研究[J].海洋工程装备与技术，2019，6（增刊 1）：225—233.

WANG G R， ZHONG L， LIU Q Y， et al. Research on marine petroleum and hydrate development technology based on dual gradient drilling of double-layer pipe [J]. Ocean Engineering Equipment and Technology， 2019，6 （S1）： 225-233.

[9]何淼，柳贡慧，李军，等.多相流全瞬态温度压力场耦合模型求解及分析[J].石油钻探技术，2015，43 （2）： 25-32.

HE M， LIU GH， LIJ， et al. Solution and analysis of fully transient temperature and pressure coupling model for multiphase flow [J]. Petroleum Drilling Tech- niques， 2015， 43 （2）： 25-32.

[10]王江帅，李军，柳贡慧，等.循环钻进过程中井筒温度场新模型[J].断块油气田， 2018， 25 （2）： 240—243.

WANG J S， LI J， LIU G H， et al. New model of wellbore temperature field during drilling process [J]. Fault-Block Oil & Gas Field， 2018， 25 （2）： 240- 243.

[11]WANG JS， LIJ， LIU G H， et al. Development and application of wellbore heat transfer model considering variable mass flow [J]. Underground Space， 2021， 6 （3）： 316-328.

[12]WANGJS， LIJ， LIU GH， et al. Development of a wellbore heat transfer model considering circulation loss [J]. Arabian Journal of Geosciences， 2020， 13（2）： 85.

[13]孙宝江，王雪瑞，孙小辉，等.井筒四相流动理论在深水钻完井工程与测试领域的应用与展望[J].天然气工业， 2020， 40 （12）： 95—105.

SUNBJ， WANG XR， SUNXH， et al. Application and prospect of the wellbore four-phase flow theory in the field of deepwater drilling and completion engineer- ing and testing [J]. Natural Gas Industry， 2020， 40 （12）： 95-105.

[14]HOLMES C S， SWIFT S C. Calculation of circulating mud temperatures [J]. Jourmal of Petroleum Technol- ogy， 1970， 22 （6）： 670-674.

[15]王江帅，李军，柳贡慧，等.变压力梯度下钻井环空压力预测[J].石油学报，2020，41 （4）：497—504.

WANG JS， LIJ， LIU G H， et al. Prediction of an- nulus pressure in variable pressure gradients drilling [J]. Acta Petrolei Sinica， 2020， 41 （4）： 497- 504.

第一作者简介：王江帅，生于1993年， 2021年获中国石油大学（北京）油气井工程博士学位，现主要从事深水钻井、控压钻井、井筒多相流与传热、石油工程人工智能方面研究工作。地址：（213100）江苏省常州市。E—mail：wjs125126@ 163. com。

收稿日期：2022—07—31

（本文编辑 宋治国）