張銳堯 李 軍 柳貢慧，2 楊宏偉

(1.中國石油大學(xué)(北京)石油工程學(xué)院 2.北京工業(yè)大學(xué))

0 引 言

隨著鉆井技術(shù)的進一步發(fā)展以及對能源需求量的增加，深水油氣勘探開發(fā)正成為世界石油工業(yè)的主要增長點和科技創(chuàng)新前沿[1-3]。盡管深水區(qū)域油氣資源豐富，但其地質(zhì)環(huán)境復(fù)雜，如深水地區(qū)海水低溫與地層高溫所形成的復(fù)雜溫度場環(huán)境，淺部地層存在弱膠結(jié)、易坍塌的特點以及低破裂壓力與高孔隙壓力所形成的安全密度窗口窄等[4-8]，使得深水鉆井過程中經(jīng)常會出現(xiàn)溢流與漏失等復(fù)雜工況。特別是當(dāng)多壓力系統(tǒng)存在時，即地層出現(xiàn)多個漏失層與溢流層的情況，則井筒內(nèi)部的鉆井液與地層流體會發(fā)生傳質(zhì)與傳熱，導(dǎo)致鉆井液的熱物性參數(shù)改變，進一步影響井筒壓力。

目前，國內(nèi)外的許多學(xué)者通過分析法或數(shù)值法對井筒溫度場進行了研究。H.J.RAMEY[9]將地層傳熱視為非穩(wěn)態(tài)而井筒傳熱視為穩(wěn)態(tài)進行溫度場計算。C.S.KABIR等[10]在不考慮鉆井液性能變化以及鉆井液流動產(chǎn)熱的條件下建立了溫度場計算模型。L.R.RAYMOND[11]所建立的溫度場模型也沒有考慮鉆井液循環(huán)過程中所產(chǎn)生的熱量，因此最終計算結(jié)果比實際值小。D.W.MARSHALL等[8]基于偏微分方程建立了從井筒到地層的井筒溫度預(yù)測模型。ZHANG Z.等[12-13]針對裸眼地層發(fā)生溢流以及關(guān)井階段，建立了井筒與地層的溫度場計算模型，通過環(huán)空溫度場的變化規(guī)律可以確定溢流發(fā)生的位置。LI M.B.等[14]建立了水平井瞬態(tài)傳熱計算模型，并利用有限體積法進行求解。盡管對井筒和地層溫度分布的相關(guān)研究已經(jīng)有豐碩的成果,但目前很少有學(xué)者對深水鉆井過程中多壓力系統(tǒng)條件下的井筒溫度場進行過研究。為此，本文對深水鉆井中多壓力系統(tǒng)在不同井深以及在不同漏失和溢流條件下的井筒溫度及鉆井液熱物性參數(shù)進行數(shù)值計算和敏感性分析，以期為深水鉆井提供一定的理論支撐，從而降低作業(yè)風(fēng)險與成本。

1 多壓力系統(tǒng)條件下井筒傳熱與傳質(zhì)模型建立

深水鉆井過程中所涉及到的傳熱過程包括：鉆柱內(nèi)鉆井液與鉆柱內(nèi)壁產(chǎn)生對流換熱；鉆柱內(nèi)壁與外壁產(chǎn)生熱傳導(dǎo)；環(huán)空中鉆井液與鉆柱外壁以及井壁產(chǎn)生對流換熱。為了建立一套完整的多壓力系統(tǒng)條件下的井筒傳熱與傳質(zhì)物理模型，需做如下假設(shè)：①鉆柱內(nèi)與環(huán)空中的流動過程考慮熱傳導(dǎo)與對流換熱，將套管、水泥環(huán)以及地層作為一個整體區(qū)域，僅考慮其在徑向上的熱傳導(dǎo)；②地層流體進入環(huán)空后與鉆井液均勻混合，溢流層與漏失層相互交替均勻分布,同類型的各子層中流體速率(漏失速率或溢流速率)相等；③忽略井筒流體在徑向上的溫度變化，僅考慮沿軸向的溫度變化；④海水與地層的熱物性參數(shù)為常數(shù)，不隨溫度而變化。根據(jù)以上假設(shè)建立物理模型，如圖1所示。

圖1 深水鉆井多壓力系統(tǒng)條件下的井筒傳熱與傳質(zhì)物理模型

該物理模型將主要的傳熱區(qū)域分為3個部分：一是鉆柱內(nèi)鉆井液所在傳熱區(qū)域(區(qū)域1)，二是環(huán)空內(nèi)鉆井液所在傳熱區(qū)域(區(qū)域2)，三是套管、水泥環(huán)與地層等只存在徑向?qū)岬膮^(qū)域(區(qū)域3)。根據(jù)漏失層與溢流層位置，將多壓力系統(tǒng)在軸向上劃分為多個子層，漏失層的漏失速率為vl(l=2,4,6，……，n)，溢流層的溢流速率為vk(k=1,3,5，……，m)，其中，n、m分別為漏失層和溢流層的層數(shù)(見圖1c)。vl與vk單位均為m3/s。多壓力系統(tǒng)位于中部地層時，其位置如圖1b所示，多壓力系統(tǒng)位于井底往上連續(xù)地層時，其位置位于井筒最下部。

1.1 數(shù)學(xué)模型

1.1.1 鉆柱內(nèi)的傳熱數(shù)學(xué)模型

鉆井液在鉆柱內(nèi)部流動過程中會與鉆柱內(nèi)壁產(chǎn)生對流換熱，同時鉆柱內(nèi)壁會與鉆柱外壁之間產(chǎn)生熱傳導(dǎo)。根據(jù)熱力學(xué)第一定律，建立如式(1)所示的數(shù)學(xué)模型。式(1)左邊表示Δt時間內(nèi)單元體內(nèi)能量的變化，式(1)右邊第1項表示Δt時間內(nèi)鉆井液與鉆柱內(nèi)壁之間的對流換熱量，右邊第2項表示Δt時間內(nèi)流入單元體的熱量，右邊第3項表示Δt時間內(nèi)鉆井液因流動摩擦所產(chǎn)生的熱量。

πUapdpi[Ta(y,t)-Tp(y,t)]+

(1)

式中：Uap為環(huán)空內(nèi)鉆井液與鉆柱內(nèi)鉆井液的綜合對流換熱系數(shù),W/(m2·K)；dpi為鉆柱內(nèi)徑,mm；Tp為鉆柱溫度,℃；Ta為環(huán)空溫度,℃；ρp為鉆柱內(nèi)鉆井液密度,kg/m3；Qcp為鉆柱內(nèi)的鉆井液微元體流動產(chǎn)生的熱量，W/m3；cp為鉆柱內(nèi)鉆井液比熱容,J/(kg·K)；Qm為鉆井液排量,m3/s；t為時間，s；y為鉆井液微元體沿井眼方向的長度，m。

1.1.2 環(huán)空中的傳熱數(shù)學(xué)模型

當(dāng)多壓力系統(tǒng)位于中部地層時(見圖1b)，鉆井液從鉆柱內(nèi)注入，經(jīng)鉆頭后再由環(huán)空返回地面。此時，下部環(huán)空位于多壓力系統(tǒng)的上游段，因此下部環(huán)空的鉆井液熱物性參數(shù)不受多壓力系統(tǒng)的影響，其傳熱過程與正常循環(huán)時相同。下部環(huán)空中的鉆井液繼續(xù)上返，當(dāng)流過多壓力系統(tǒng)段的漏失層位時，部分鉆井液會進入地層中，而在溢流層位則有地層流體進入環(huán)空中，并與鉆井液均勻混合后上返至地面。在整個循環(huán)過程中，上部環(huán)空中的鉆井液熱物性參數(shù)會受到一定的影響，而下部環(huán)空中的鉆井液熱物性參數(shù)則不受多壓力系統(tǒng)的影響。根據(jù)上、下環(huán)空中的傳熱情況不同，本文分別建立不同的傳熱方程。

下部環(huán)空鉆井液傳熱方程如式(2)所示。式(2)左邊第1項表示Δt時間內(nèi)環(huán)空單元體內(nèi)能的變化，式(2)右邊第1項和第2項分別表示環(huán)空內(nèi)鉆井液與鉆柱外壁以及井壁之間對流換熱的熱量，右邊第3項表示流入單元體內(nèi)的熱量，第4項表示環(huán)空內(nèi)鉆井液流動過程中摩擦所產(chǎn)生的熱量。

πUapdpo[Ta(y,t)-Tp(y,t)]+

πUafdw[Ta(y,t)-Tf(y,t)]+

(2)

式中：dw為井眼內(nèi)徑，mm；dpo為鉆柱外徑，mm；ca為環(huán)空內(nèi)鉆井液比熱容，J/(kg·K)；Uaf為環(huán)空內(nèi)鉆井液與井壁的綜合對流換熱系數(shù),W/(m2·K)；Tf為地層溫度,℃；ρa為環(huán)空鉆井液密度,kg/m3；Qca為正常循環(huán)時環(huán)空中鉆井液微元體摩擦產(chǎn)生的熱量,W/m3。

當(dāng)多壓力系統(tǒng)位于中部地層時，因為漏失與溢流同時存在，所以在溢流層位會有地層流體進入環(huán)空，而漏失層位會有部分鉆井液進入地層中，從而使多壓力系統(tǒng)所在井深處的環(huán)空內(nèi)鉆井液熱物性參數(shù)發(fā)生變化，并且隨著鉆井液上返，整個上部環(huán)空中的傳熱過程都會受到影響。因此建立如下傳熱方程：

(3)

當(dāng)多壓力系統(tǒng)位于井底往上連續(xù)地層時，在井底附近環(huán)空中，部分鉆井液會因為漏失進入地層中，同時在溢流層位中會有地層流體進入環(huán)空中；然后鉆井液與地層流體在井底附近的環(huán)空區(qū)域混合均勻后，上返至地面。因此，從井底開始出現(xiàn)溢流、漏失至上返到地面的過程中，整個環(huán)空中的鉆井液熱物性參數(shù)都會受到多壓力系統(tǒng)的影響，于是得到如下傳熱方程：

(4)

1.1.3 地層中的傳熱數(shù)學(xué)模型

地層中存在流體與否對于傳熱過程有較大的影響。此次研究只考慮存在流體的情況，地層中流體的流動可以看作是單向不可壓縮平面徑向穩(wěn)定滲流，其傳熱數(shù)學(xué)模型見文獻[15]。

1.1.4 密度模型

在鉆井循環(huán)過程中，由于鉆井液的熱物性參數(shù)與溫度、壓力之間相互影響，所以本文利用多元非線性回歸分析方法對W.C.MCMORDIE等[16]的水基鉆井液試驗數(shù)據(jù)進行處理，得到鉆井液的密度ρ與溫度T、壓力p的關(guān)系式：

(5)

1.2 初始條件與邊界條件

1.2.1 初始條件

海水的溫度分布受季節(jié)、水流和深度等多種因素的綜合影響，而海水的溫度變化又會影響井筒的傳熱過程。因此，本文將V.DOKHANI等[17]的海水溫度分布數(shù)學(xué)模型作為初始條件。

1.2.2 邊界條件

井口溫度可以看作是鉆柱內(nèi)部的初始溫度，其邊界條件可以表示為:

T0(y=0,t=0)=Tin

(6)

根據(jù)假設(shè)條件③，不考慮井筒溫度在徑向上的變化，則在井底處鉆柱內(nèi)鉆井液、鉆頭處以及環(huán)空中的鉆井液位于同一井深位置，三者對應(yīng)的溫度相等，其邊界條件可以表示為:

T0(y=H,t)=T1(y=H,t)=T2(y=H,t)

(7)

遠(yuǎn)離井筒的環(huán)境溫度分布邊界條件為:

T(x→∞,y,t)=Tsurf+Gh

(8)

式中：T0為鉆柱內(nèi)鉆井液在井底的溫度,℃；Tin為鉆柱的入口溫度,℃。T1為鉆柱壁在井底的溫度,℃；T2為環(huán)空內(nèi)鉆井液在井底的溫度,℃；G為地溫梯度,℃/m；H為井深，m；h為泥線到目標(biāo)地層的深度，m。

2 模型驗證與敏感性分析

通過有限差分方法對模型進行離散，然后結(jié)合南海某區(qū)塊的鉆井?dāng)?shù)據(jù)，通過高斯-賽德爾迭代的方法對模型進行求解。將求解結(jié)果與現(xiàn)有經(jīng)典模型以及現(xiàn)場測量數(shù)據(jù)進行對比分析，實現(xiàn)對本文模型的驗證。對不同影響因素下的環(huán)空溫度、井底溫度以及鉆井液密度等進行數(shù)值計算和敏感性分析。

2.1 模型驗證

選擇南海某區(qū)塊的鉆井?dāng)?shù)據(jù)進行計算，其關(guān)鍵參數(shù)為[18-19]：井深5 094 m，水深1 521 m，入口溫度15 ℃，地表溫度20 ℃，地溫梯度0.024 ℃/m，鉆井液密度1 200 kg/m3，地層巖石密度2 640 kg/m3，海水密度1 020 kg/m3，鉆井液比熱容3 935 J/(kg·K)，鉆井液與巖石的導(dǎo)熱系數(shù)分別為1.73和2.10 W/(m·K)，海水的導(dǎo)熱系數(shù)為0.60 W/(m·K)，孔隙流體密度為1 400 kg/m3，孔隙流體導(dǎo)熱系數(shù)為0.8 W/(m·K)，所用鉆頭直徑為215.9 mm。

將以上參數(shù)代入本文模型以及其他幾個模型進行計算，得到如圖2a所示結(jié)果。通過模型對比，本文模型與張正模型的計算結(jié)果很接近。因為Marshall模型沒有考慮井身結(jié)構(gòu)的影響，所以下部環(huán)空鉆井液從地層中獲得的熱量多，上部環(huán)空鉆井液向海水傳遞的熱量多。因此，下部環(huán)空溫度高于本文計算結(jié)果，上部環(huán)空溫度低于本模型計算結(jié)果。

楊謀模型主要適用于傳統(tǒng)的陸地鉆井，并且模型假設(shè)條件差異導(dǎo)致計算結(jié)果與本文模型有較大誤差。張正模型對漏失工況下的環(huán)空溫度進行了計算，并與其他經(jīng)典模型以及現(xiàn)場數(shù)據(jù)進行了對比[20-22]，驗證了該模型的合理性與精確度。本文模型不僅考慮了單獨漏失與單獨溢流情況，同時還考慮了存在多壓力系統(tǒng)時的工況，并且本文模型與張正模型計算的環(huán)空溫度誤差最大僅為3.31 ℃，滿足誤差小于5%的要求，驗證了本文模型的正確性，說明本文模型適用范圍更廣。圖2b為循環(huán)穩(wěn)定后某段時間的井口測量溫度與本文模型計算的井口溫度對比，因為模型是在假設(shè)條件下計算，所以二者之間存在一定誤差，計算值與測量值誤差在5%以內(nèi)，同樣滿足精度要求。綜上所述，本文模型預(yù)測結(jié)果具有準(zhǔn)確性與合理性。

圖2 多壓力系統(tǒng)條件下井筒溫度場模型驗證

2.2 敏感性分析

當(dāng)多壓力系統(tǒng)分別位于中部地層(2 500～3 000 m)或井底往上連續(xù)地層(井底往上500 m)時，結(jié)合上述鉆井?dāng)?shù)據(jù)對不同溢流和漏失速率條件下的環(huán)空溫度、環(huán)空與鉆柱內(nèi)的溫度差、井底溫度以及環(huán)空內(nèi)鉆井液密度的變化進行了計算和敏感性分析。

2.2.1 多壓力系統(tǒng)對環(huán)空溫度的影響

圖3和圖4分別表示多壓力系統(tǒng)位于中部地層以及井底往上連續(xù)地層時的環(huán)空溫度分布。從圖3和圖4可以看出：當(dāng)多壓力系統(tǒng)位于中部地層時，環(huán)空鉆井液溫度受到影響最大的井段為2 500～3 500 m(圖3區(qū)域1)；當(dāng)多壓力系統(tǒng)位于井底往上連續(xù)地層時(井底往上500 m)，則從井底至1 500 m的環(huán)空內(nèi)鉆井液溫度都受到多壓力系統(tǒng)的影響(圖4區(qū)域2)。無論多壓力系統(tǒng)位于中部地層或井底往上連續(xù)地層，在不同溢、漏工況(溢流與漏失速率不同)時，只要溢流速率大于漏失速率，則環(huán)空溫度都要高于正常循環(huán)時的溫度，反之環(huán)空溫度低于正常循環(huán)時的溫度。圖5和圖6分別表示多壓力系統(tǒng)位于中部地層與井底往上連續(xù)地層時，溢流同存條件與單獨溢流或者漏失條件對環(huán)空溫度的影響。從圖5和圖6可以看出：當(dāng)多壓力系統(tǒng)位于中部地層時，環(huán)空溫度分布受影響區(qū)域主要為多壓力系統(tǒng)往上的環(huán)空區(qū)域(圖5區(qū)域1)；當(dāng)多壓力系統(tǒng)位于井底往上連續(xù)地層時，從井底直到海底泥線處環(huán)空溫度都受到顯著影響，所以后者受到的影響范圍更廣；當(dāng)多壓力系統(tǒng)位置一定時，在相同井深處，出現(xiàn)單獨溢流時環(huán)空溫度最高，而單獨漏失時環(huán)空溫度最低，多壓力系統(tǒng)(溢漏同存)對環(huán)空溫度分布的影響在上述兩者之間。

圖3 位于中部地層時的環(huán)空溫度

圖4 位于井底往上連續(xù)地層時的環(huán)空溫度

圖5 位于中部地層且不同溢漏條件下的環(huán)空溫度

圖6 在井底往上連續(xù)地層且不同溢漏條件下的環(huán)空溫度

2.2.2 多壓力系統(tǒng)對環(huán)空與鉆柱內(nèi)溫度差的影響

圖7和圖8分別表示多壓力系統(tǒng)位于中部地層或井底往上連續(xù)地層時，環(huán)空與鉆柱內(nèi)溫度差隨井深的分布。

圖7 位于中部地層時環(huán)空與鉆柱內(nèi)的溫度差

圖8 位于井底往上連續(xù)地層時環(huán)空與鉆柱內(nèi)的溫度差

由圖7可知：當(dāng)多壓力系統(tǒng)位于中部地層時，溫度差的變化最顯著的井段位于多壓力系統(tǒng)區(qū)域附近(中部井段)；當(dāng)多壓力系統(tǒng)位于井底往上連續(xù)地層時，其對溫度差的影響從井底區(qū)域一直到泥線附近都比較顯著。無論多壓力系統(tǒng)位于中部地層還是井底往上連續(xù)地層，只要出現(xiàn)溢流，則溫度差都明顯增加，而出現(xiàn)漏失時，溫度差明顯減小，多壓力系統(tǒng)條件下(溢漏同存)的溫度差隨井深的分布介于兩者之間。在多壓力系統(tǒng)條件下，當(dāng)溢流速率大于漏失速率時，在相同井深位置，其溫度差分布要大于溢流速率小于漏失速率條件下的分布，而正常循環(huán)時的溫度差分布在兩者之間。

2.2.3 多壓力系統(tǒng)對井底溫度的影響

圖9和圖10分別為多壓力系統(tǒng)位于中部地層與井底往上連續(xù)地層時的井底溫度。

圖9 位于中部地層時的井底溫度

圖10 位于井底往上連續(xù)地層時的井底溫度

從圖9和圖10可以看出，當(dāng)其他條件不變時，隨著循環(huán)時間的延長，井底溫度逐漸降低，然后趨于穩(wěn)定。這是因為鉆井液在循環(huán)過程中，井筒與地層之間不斷進行熱量交換，下部環(huán)空中的鉆井液會吸收地層中的熱量，然后傳遞給上部環(huán)空，但是隨著循環(huán)時間的延長，井筒與地層經(jīng)過充分的熱量交換，最終地層溫度保持相對恒定，井底溫度也趨于穩(wěn)定。當(dāng)循環(huán)相同時間，如果只出現(xiàn)溢流時，則井底溫度最高；只出現(xiàn)漏失時，井底溫度最低。當(dāng)多壓力系統(tǒng)位于井底時，如果溢流速率大于漏失速率，則其井底溫度低于單獨溢流時的井底溫度，而高于正常循環(huán)時的井底溫度。如果溢流速率小于漏失速率，則其井底溫度高于單獨漏失時的井底溫度，低于正常循環(huán)時的井底溫度。與多壓力系統(tǒng)位于中部地層時相比，當(dāng)其位于井底時，如果都出現(xiàn)溢流速率小于漏失速率的情況，則后者對井底溫度的影響更大。

2.2.4 多壓力系統(tǒng)對環(huán)空內(nèi)鉆井液密度的影響

圖11和圖12分別表示多壓力系統(tǒng)位于中部地層和井底往上連續(xù)地層時，不同溢流或漏失工況對鉆井液密度的影響。

圖11 位于中部地層時的鉆井液密度

圖12 位于井底地層時的鉆井液密度

從圖11和圖12可見：在多壓力系統(tǒng)條件下，當(dāng)溢流速率大于漏失速率時，相同井深處鉆井液密度最大，與正常循環(huán)時相比顯著增加；當(dāng)溢流速率小于漏失速率時，鉆井液密度仍然大于正常循環(huán)時的密度，但小于溢流速率大于漏失速率工況下的密度；僅出現(xiàn)溢流時，鉆井液密度在溢流速率大于漏失速率與溢流速率小于漏失速率兩種工況的密度之間，但是明顯大于只有漏失或正常循環(huán)時的鉆井液密度。僅出現(xiàn)漏失時，由于鉆井液的排量減小，但是其熱物性參數(shù)沒有改變，所以與正常循環(huán)時相比，密度基本保持不變。無論多壓力系統(tǒng)位于中部地層還是井底往上連續(xù)地層，只要出現(xiàn)溢流工況，則環(huán)空中鉆井液的密度都顯著增加，并且溢流速率大于漏失速率時，環(huán)空鉆井液密度最大，而漏失工況對環(huán)空鉆井液密度基本沒有影響。

3 結(jié) 論

根據(jù)多壓力系統(tǒng)位于不同位置時井筒內(nèi)傳熱與傳質(zhì)的特點，并依據(jù)熱力學(xué)第一定律，推導(dǎo)了多壓力系統(tǒng)條件下的井筒溫度場數(shù)學(xué)模型。分別采用有限差分法和高斯塞德爾迭代法對模型進行離散和求解。將本文模型與現(xiàn)有經(jīng)典模型以及實測數(shù)據(jù)進行對比，驗證了模型的可靠性。最后對影響井筒溫度與熱物性參數(shù)的因素進行了數(shù)值計算和敏感性分析，得出以下結(jié)論。

(1)多壓力系統(tǒng)位于中部地層時，從該區(qū)域到泥線處的環(huán)空溫度及鉆井液熱物性參數(shù)會受到顯著影響；當(dāng)其位于井底往上連續(xù)地層時，則從井底直到泥線處整個環(huán)空溫度分布以及鉆井液熱物性參數(shù)都會受到顯著影響。

(2)單獨溢流或漏失對環(huán)空溫度和井底溫度的影響程度要大于兩者同存時的情況，并且相同井深位置，只存在溢流時環(huán)空溫度最高，只存在漏失時環(huán)空溫度最低，多壓力系統(tǒng)條件下(溢漏同存)介于兩者之間。

(3)在多壓力系統(tǒng)條件下，只要出現(xiàn)溢流工況，則環(huán)空中鉆井液的密度都顯著增加，并且在溢流速率大于漏失速率時，環(huán)空鉆井液密度最大，而漏失工況對環(huán)空鉆井液密度基本沒有影響。