海洋雙層管雙梯度鉆井井筒溫度分布規律研究

王江帥 付盼 胡旭輝 宮臣興 鄧嵩 唐政 殷文

摘要：井筒溫度對天然氣水合物的存在狀態有顯著影響。為了準確掌握海洋雙層管雙梯度鉆井井筒溫度分布規律，基于傳熱學理論，充分考慮巖屑對內管混合流體熱物性參數的影響，以及隔水管、套管、裸眼等井段環空內不同流體類型對傳熱的影響，引入綜合導熱系數，建立了適用于雙層管雙梯度鉆井的井筒溫度場預測模型，并探究了井筒溫度分布規律。研究結果發現：對雙層管雙梯度鉆井井筒而言，其在海底附近的溫度達到最低值，在水合物固態流化開采時應特別注意該井段的井筒溫度分布；注入流體溫度、排量和鉆井液比熱容的增大，會導致內管流體溫度剖面整體增大；注入流體溫度主要影響井口附近的內管流體溫度，而排量和鉆井液比熱容主要影響海底附近的內管流體溫度，可以通過調節上述參數來改善不同井段井簡溫度分布，以減少水合物的生成。研究結果可為水合物的安全高效開采提供理論支撐。

關鍵詞：天然氣水合物；雙層管雙梯度鉆井；井筒溫度；綜合導熱系數；內管；外管環空

中圖分類號：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

（本文編輯 宋治國）