張立剛， 譚甲興， 陶澤俊， 王倩穎， 李 浩， 婁絹平， 孫道坤， 薛東陽

( 1. 東北石油大學(xué) 提高油氣采收率教育部重點實驗室，黑龍江 大慶 163318； 2. 大慶鉆探工程公司 井下作業(yè)工程公司，吉林 松原 138000； 3. 大慶油田有限責(zé)任公司 第二采油廠，黑龍江 大慶 163414 )

0 引言

油田開發(fā)過程中，原油結(jié)蠟問題普遍存在，當(dāng)油井井筒溫度處于析蠟點和凝點溫度區(qū)域時，部分蠟晶析出并沉積在管壁上，形成較長井段的結(jié)蠟層[1-3]，析出的石蠟可能發(fā)生在原油流經(jīng)的任何場所，如生產(chǎn)油管、潛油泵、輸油管、分離器、油罐及其他地面設(shè)備中，造成油流通道面積減小，阻力增加，流量下降，嚴重時甚至停產(chǎn)[4-5]?，F(xiàn)有的清蠟技術(shù)主要有熱處理法、化學(xué)反應(yīng)法、機械清除法、磁流體控制技術(shù)、微生物代謝產(chǎn)物清蠟法和超聲波處理技術(shù)等，其中反循環(huán)熱洗是目前主要的有桿泵井熱處理法清蠟方式[6-8]。反循環(huán)熱洗是在井口向套管內(nèi)注入熱水，熱水沿油管外壁流下，在油管和抽油桿環(huán)空中將蠟逐漸全部熔化并及時排出，由于以接近抽油機理論排量的方式注入熱水并及時排出，抽油井無恢復(fù)期或恢復(fù)期很短，對儲層的污染可降低至最低[9-10]。采液過程及熱洗過程中井筒溫度場的分布，是結(jié)蠟位置分析和熱洗方案制定的依據(jù)，如果熱洗過程中結(jié)蠟段的井筒溫度低于蠟熔化溫度，將無法實現(xiàn)有效清蠟，為保障反循環(huán)熱洗清蠟效果，有必要分析不同工況下采液過程和熱洗過程中的井筒溫度場分布。

采出液與洗井液在井筒中的流動過程及與周圍環(huán)境介質(zhì)的熱交換過程是典型的熱流耦合過程，不同的產(chǎn)液參數(shù)和熱洗參數(shù)組合對井筒溫度場影響很大。根據(jù)傳熱學(xué)原理，王梟建立井筒內(nèi)循環(huán)介質(zhì)對流傳導(dǎo)方程，編制熱洗過程中油套環(huán)空內(nèi)的熱水溫度計算軟件，求解井身結(jié)構(gòu)上的溫度曲線[11]。張克舫等以環(huán)空中的注入水和油管中的返回水為研究對象，假設(shè)套管內(nèi)和油管內(nèi)垂直于流動方向的任一截面內(nèi)的流體溫度是均勻的，推導(dǎo)熱洗水在注入過程中及返回過程中的耦合流動和傳熱數(shù)學(xué)模型，采用數(shù)學(xué)分析法求解熱洗水在注入過程和采出過程中的溫度分布的分析解[12]。廖凱[13]、金峰等[14]采用井溫測井方法，測試洗井前、洗井中及洗井后的井筒溫度場，評價抽油井熱洗過程中油套環(huán)空中溫度場的變化。宋奇等[15]等將油管壁、油套環(huán)空、套管壁和水泥環(huán)產(chǎn)生的熱阻等效為串聯(lián)總熱阻，推導(dǎo)地層—井筒傳熱模型，計算產(chǎn)液過程中的井筒溫度場分布，并預(yù)測結(jié)蠟位置。這些研究多關(guān)注地層和井筒的傳熱過程，對套管—套管、油管環(huán)空流體—油管—油管、抽油桿環(huán)空流體—抽油桿之間軸向和徑向上的熱傳導(dǎo)和熱對流進行簡化，制約對熱洗過程中套管和油管環(huán)空、油管和抽油桿環(huán)空溫度場分布的認識。

筆者以有桿泵井為研究對象，考慮地層—水泥環(huán)—套管—套管、油管環(huán)空流體—油管—油管和抽油桿環(huán)空流體—抽油桿之間的耦合傳熱過程，模擬分析采液過程和反循環(huán)熱洗過程井筒溫度場的分布，為結(jié)蠟位置分析、反循環(huán)熱洗溫度、壓力和排量等方案的制定提供依據(jù)。

1 熱流耦合模型

基于有桿泵井采出液與洗井液在井筒中的流動過程及與周圍環(huán)境介質(zhì)的熱交換過程，做基本假設(shè)：(1)固井狀況良好，套管、水泥環(huán)和地層間各界面無脫開及相互滑移；(2)反循環(huán)熱洗過程中地層無吸液；(3)井筒沒有采取加熱措施，不存在內(nèi)熱源；(4)忽略結(jié)蠟層對管壁導(dǎo)熱系數(shù)的影響。

數(shù)值計算中描述熱流耦合的控制方程[16-18]包括連續(xù)性方程、動量方程和能量守恒方程：

·(Aρu)=0,

(1)

(2)

(3)

式中：A為流體過流面積；ρ為流體密度；u為流體流動速度；p為壓力；F為體積力(重力)；Cp為常壓比熱容；T為溫度；k為導(dǎo)熱系數(shù)；Qwall流體域和周圍環(huán)境熱交換量；fD為摩擦因數(shù)；dh為過流當(dāng)量直徑。

運用Churchill方程，可得層流、湍流和過度流區(qū)域的摩擦因數(shù)fD：

(4)

(5)

(6)

(7)

式中：Re為雷諾數(shù)；C、D為中間變量；e為管壁絕對粗糙度；d為水力直徑；μ為洗井液動力黏度。

2 產(chǎn)液過程

2.1 物理模型

根據(jù)井筒結(jié)構(gòu)，建立以井眼軸線為中心的抽油桿—油管—套管—水泥環(huán)—地層的軸對稱有限元計算模型(見圖1(a))，井筒長度為1 000.0 m，抽油桿直徑為22.0 mm，油管內(nèi)徑為62.0 mm，外徑為73.0 mm，套管內(nèi)徑為124.0 mm，外徑為139.7 mm，井眼直徑為215.9 mm。介質(zhì)材料參數(shù)：水和原油的比熱容分別為4 200、2 200 J/(kg·K)，導(dǎo)熱系數(shù)分別為0.50、0.32 W/(m·K)，套管和油管的比熱容和導(dǎo)熱系數(shù)分別為460 J/(kg·K)和18.50 W/(m·K)，水泥環(huán)和地層的導(dǎo)熱系數(shù)為5.00 W/(m·K)。邊界條件：遠場地層的溫度邊界為地層溫度，流場邊界條件施加于油管底部，按產(chǎn)液量施加流速邊界，溫度按油層溫度設(shè)定為48.5 ℃。對物理模型進行三角形網(wǎng)格剖分和迭代求解，可得產(chǎn)液過程井筒溫度場分布云圖(見圖1(b))(產(chǎn)液量為10 t/d，純原油)。由圖1(b)可知，采出液沿油管—抽油桿環(huán)空舉升過程中溫度不斷降低，油溫降到析蠟點后，原油中的高碳數(shù)蠟將結(jié)晶、析出并粘附在管壁上。

2.2 油管和抽油桿環(huán)空溫度場

模擬計算不同含水率和產(chǎn)液量下油管和抽油桿環(huán)空溫度場，含水率為90%，不同產(chǎn)液量下溫度場分布規(guī)律見圖2。由圖2(a)可知，隨產(chǎn)出液從井底向井口運移，溫度逐漸下降，不同井段溫降速率不同，井深700 m以下時，溫降速率較小，溫度變化平緩；井深700 m至井口時，溫降速率逐漸增大。由圖2(b)可知，同一截面位置，隨產(chǎn)液量增加，溫度升高，在產(chǎn)液量為0～85 t/d時，溫升速率較快；當(dāng)產(chǎn)液量超過85 t/d時，溫度趨于平穩(wěn)。

產(chǎn)液量為110 t/d時，不同含水率下溫度場分布規(guī)律見圖3。由圖3可知，隨含水率增加，井下各點溫度逐漸升高，不同深度位置溫升速率不同，井下越深位置溫度變化越平緩。結(jié)合析蠟溫度，可判斷不同工況下的結(jié)蠟位置，如果一口井的產(chǎn)液量越高，則含水率越高，出現(xiàn)的結(jié)蠟段位置越淺。

圖1 有桿泵井產(chǎn)液過程物理模型及井筒溫度場云圖Fig.1 Physical model and wellbore temperature field of fluid production process in rod pump well

圖2 不同產(chǎn)液量油管和抽油桿環(huán)空溫度場分布Fig.2 Temperature field in annular of tubing and sucker rod with different liquid yield

圖3 不同含水率油管和抽油桿環(huán)空溫度場分布Fig.3 Temperature field in annular of tubing and sucker rod with different moisture content

3 反循環(huán)熱洗過程

3.1 物理模型

反循環(huán)洗井過程中，熱水在油管和套管環(huán)空中向下流動，在油管和抽油桿環(huán)空返出，物理模型見圖4(a)，介質(zhì)材料屬性與2.1相同。邊界條件：遠場地層的溫度邊界為地層溫度，套管和油管環(huán)空頂部按洗井排量和溫度施加流速邊界和溫度邊界，油管和抽油桿環(huán)空頂部施加壓力邊界。徑向上洗井液與油管內(nèi)流體、套管、水泥環(huán)及地層換熱，軸向上存在熱對流，形成穩(wěn)定的溫度場。當(dāng)洗井排量為16 m3/d、洗井溫度為70 ℃時，井筒溫度場分布規(guī)律見圖4(b)。由圖4(b)可知，隨洗井液向下流動，油管和套管環(huán)空溫度逐漸降低，隨洗井液上返，油管內(nèi)溫度逐漸升高。

圖4 有桿泵井反循環(huán)熱洗物理模型及井筒溫度場云圖Fig.4 Physical model and wellbore temperature field of indirect thermal washing process in rod pump well

3.2 油管和套管環(huán)空溫度場

洗井排量為16 m3/h時，不同洗井溫度下油管和套管環(huán)空溫度分布見圖5。由圖5(a)可知，當(dāng)洗井排量一定時，油管和套管環(huán)空溫度隨井深增加而逐漸降低，溫降速率逐漸變小而趨于平緩。由圖5(b)可知，隨洗井井口溫度升高，油管和套管環(huán)空各點溫度呈指數(shù)函數(shù)增加。

圖5 不同洗井溫度下油管和套管環(huán)空溫度分布Fig.5 Temperature field in annular of casing and tubing with different washing temperature

洗井溫度為80 ℃時，不同洗井排量下油管和套管環(huán)空溫度分布見圖6。由圖6可知，當(dāng)洗井溫度一定時，隨洗井排量增加，油管和套管環(huán)空溫降速率變低，各點溫度增高。

3.3 油管和抽油桿環(huán)空溫度場

洗井排量為16 m3/h時，不同洗井溫度下油管和抽油桿環(huán)空溫度分布見圖7。由圖7(a)可知，隨洗井液在油管和抽油桿環(huán)空向上返排，油管和抽油桿環(huán)空溫度逐漸升高，各點溫升速率不同，越往井口溫升速率越快。由圖7(b)可知，當(dāng)洗井排量一定時，隨洗井井口溫度增加，油管和抽油桿環(huán)空各點溫度呈指數(shù)函數(shù)增加。

圖6 不同洗井排量下油管和套管環(huán)空溫度分布Fig.6 Temperature field in annular of tubing and casing with different flow rate

圖7 不同洗井溫度下油管和抽油桿環(huán)空溫度分布Fig.7 Temperature field in annular of tubing and sucker rod with different washing temperature

洗井溫度為80 ℃時，不同洗井排量下油管和抽油桿環(huán)空溫度分布見圖8。由圖8可知，當(dāng)洗井溫度一定時，隨洗井排量增加，油管和抽油桿環(huán)空各點溫度升高。

4 數(shù)值模擬驗證

為了驗證數(shù)值模擬結(jié)果的可靠性，對L10-2516、L10-3102井進行洗井溫度測試。L10-2516和L10-3102井原油含蠟量較高，井下400 m以上全部結(jié)蠟。2口井采用反循環(huán)熱洗清蠟，熱洗排量為16 m3/h，熱洗溫度分別為81.5、91.0 ℃。在熱洗過程中，將5個井溫測試儀下入油管和套管環(huán)形空間。先記錄從井口120～920 m(泵掛深度)每200 m各點的溫度。待溫度穩(wěn)定后，將測試儀上提100 m測試，完成10個點的數(shù)據(jù)錄取，監(jiān)測熱洗液的溫度變化。將2口井洗井過程中的油管和套管環(huán)空溫度測試結(jié)果和數(shù)值模擬結(jié)果進行對比(見表1)。由表1可知，2口井反循環(huán)熱洗過程中的油管和套管環(huán)空溫度測試和數(shù)值模擬結(jié)果的符合率在96.59%～99.74%之間，根據(jù)數(shù)值模擬結(jié)果可以指導(dǎo)產(chǎn)液過程的結(jié)蠟位置判斷和洗井參數(shù)優(yōu)化。

圖8 不同洗井排量下油管和抽油桿環(huán)空溫度分布Fig.8 Temperature field in annular of tubing and sucker rod with different flow rate

表1 L10-2516和L10-3102井溫度測試與數(shù)值模擬結(jié)果

5 結(jié)論

(1)運用熱流耦合模型，通過有限元數(shù)值模擬分析油井生產(chǎn)過程中不同采液量、含水率等工況下井筒溫度分布規(guī)律。隨產(chǎn)出液從井底向井口運移，溫度逐漸降低，產(chǎn)液量越大，溫降速率度越低；含水率越高，溫降速率越慢。

(2)分析反循環(huán)熱洗過程中不同排量和洗井溫度等工況下油管內(nèi)及油管和套管環(huán)空溫度分布規(guī)律。自上而下，油管和套管環(huán)空溫度逐漸降低，油管內(nèi)溫度逐漸升高；隨洗井溫度和洗井排量增加，井筒內(nèi)各點溫度呈指數(shù)函數(shù)增加。

(3)2口井的數(shù)值模擬結(jié)果和測試結(jié)果符合率在96.59%～99.74%之間，表明數(shù)值模擬結(jié)果可以指導(dǎo)產(chǎn)液過程的結(jié)蠟位置判斷和洗井參數(shù)優(yōu)化。