秦彥斌，朱禮濤，鄭 杰，竇益華

?

高溫高壓氣井生產(chǎn)過程井筒溫度場分析

秦彥斌，朱禮濤，鄭 杰，竇益華

（西安石油大學(xué)機械工程學(xué)院，陜西西安 710065）

高溫高壓氣井在生產(chǎn)過程中受到地層高溫流體的影響，井筒溫度原有的平衡被打破，井筒溫度重新分布會引起環(huán)空壓力增高，威脅井筒安全服役和井筒的完整性。為了準(zhǔn)確預(yù)測井筒溫度，基于質(zhì)量、動量、能量守恒、傳熱學(xué)、井筒傳熱理論，再考慮氣體焦耳-湯姆遜效應(yīng)、氣體溫度、壓力、密度及物性參數(shù)的影響，建立井筒溫度預(yù)測模型；將流體物性參數(shù)根據(jù)不同的溫度壓力分段計算，可提高模型計算的精確性。最后，通過實例計算分析了環(huán)境溫度的影響因素。

井筒溫度場；高溫高壓氣體；參數(shù)計算

高溫高壓（HTHP）井是以普通橡膠密封性能來界定的，指井底溫度高于150 ℃、壓力高于70 MPa的井。由高溫高壓井的定義可知，高溫高壓氣井生產(chǎn)過程中溫度、壓力較高，井筒附近的溫度分布影響流體的物性參數(shù)，在不同溫度、壓力條件下氣體的熱物性參數(shù)變化較大，所對應(yīng)的流體流動性差別加大。另外，由于流體和井筒周圍地層之間存在著溫差，必然會向周圍地層導(dǎo)熱；在徑向?qū)岬倪^程中，環(huán)空溫度穩(wěn)態(tài)被打破，引起環(huán)空溫度上升及環(huán)空帶壓現(xiàn)象，環(huán)空壓力隨著溫差的加大而急劇增加。環(huán)空溫度升高，套管軸向變形增加，有上抬井口的風(fēng)險[1-2]。因此，建立精確的井筒溫度預(yù)測模型，對管柱安全服役及井筒完整性至關(guān)重要。國外學(xué)者Ramey[3]建立了單相理想氣體在井筒流動時溫度分布模型，Hasan 和 Kabir[4]提出了氣舉井溫度分布的半解析解，國內(nèi)學(xué)者毛偉[5]、朱德武[6]、楊進[7]、張波[8]等人建立了井筒瞬態(tài)溫度計算模型。以上學(xué)者在建井溫度預(yù)測模型時，沒有考慮氣體在流動過程中壓力、溫度與氣體的密度、流速、物性參數(shù)之間的聯(lián)系；實際情況下氣體流動時，流速和密度與溫度、壓力相互影響，溫度和壓力與物性參數(shù)相互影響。聯(lián)系流速、密度、溫度及壓力的關(guān)系，結(jié)合溫度和壓力與流體的物性參數(shù)的關(guān)系，以少量含硫天然氣單向流為研究對象，筆者建立流體在油管中的瞬態(tài)溫度模型，對于高含硫井，可改變流體參數(shù)進行計算。

結(jié)合實際井身結(jié)構(gòu)，根據(jù)能量、質(zhì)量、動量守恒以及傳熱學(xué)建立井筒溫度預(yù)測模型，通過實例計算分析井筒溫度影響因素。

1 井筒溫度分布模型

1.1 流體溫度模型

高溫高壓氣井在生產(chǎn)時，井筒內(nèi)流體溫度高于地層溫度，井筒內(nèi)的流體與地層之間會發(fā)生熱量傳遞，從而引起井筒溫度的改變。為了簡化計算模型，做出如下假設(shè)：①井筒中為穩(wěn)態(tài)傳熱，地層為非穩(wěn)態(tài)傳熱；②不考慮軸向傳熱，只考慮徑向傳熱；③地層溫度呈線性分布；④流體在井筒內(nèi)穩(wěn)定流動。

圖1 井筒微元示意圖

將（1）式化簡為：

將（2）式以焓變的形式表達，假設(shè)流體穩(wěn)態(tài)流動，可得微分式：

微元段井筒熱量損失就等于流體傳向井筒外壁的熱量，可得：

井筒微元段水泥環(huán)外壁向地層傳遞的熱量：

根據(jù)熱力學(xué)定律可知：

由氣體比定壓熱熔以及焦耳–湯姆遜系數(shù)定義可將（6）式改寫為：

單相流的壓力梯度公式為：

將（4）、（5）、（7）、（8）式帶入（3）式得：

流量公式：

將（12）式帶入（11）式，則其通解為：

1.2 環(huán)空溫度計算

圖2 高溫高壓井結(jié)構(gòu)簡圖及傳熱示意圖

由（4）與（5）式可得：

在上述公式中：

1.3 物性參數(shù)的計算

1.3.1 定壓比熱采用文獻[6]給出的模型

1.3.2 天然氣的壓縮因子采用分段計算模型

在計算天然氣壓縮因子時需要先得到對比壓力和對比溫度。下面先介紹計算天然氣的對比壓力和對比溫度的算法。

在實際生產(chǎn)過程中，天然氣并不是單一物質(zhì)，大多含有非烴類物質(zhì)，為控制計算模型誤差，本文對計算天然氣擬臨界壓力和擬臨界溫度采用組分析的方法：

對比壓力與對比溫度就是壓力與溫度和其對應(yīng)的臨界參數(shù)之比：

國內(nèi)學(xué)者李相方將石油工業(yè)中幾種主流的天然氣壓縮因子計算模型，分別將幾種模型從低溫、低壓到高溫、高壓的條件下進行數(shù)值計算，對計算結(jié)果進行系統(tǒng)的比較并分析計算精度，最后給出了不同壓力、溫度下采用不同的計算模型[10]：

焦耳通過大量的實驗對氣體的內(nèi)能進行研究并提出了焦耳定律。焦耳–湯姆遜系數(shù)的定義是在等焓過程中降低單位壓力所產(chǎn)生的溫度變化或者是溫度隨時間的變化率。在文獻調(diào)研的過程中，許多文獻在計算井筒溫度時都用到焦耳–湯姆遜系數(shù)，基本上根據(jù)VDW、RK、SRK和PR四個立方狀態(tài)方程推得。本文由計算精度較高的SRK立方狀態(tài)方程推得計算焦耳–湯姆遜系數(shù)公式：

2 氣井井筒溫度分布實例分析

2.1 實例計算

根據(jù)典型的高溫高壓氣井井身結(jié)構(gòu)，推導(dǎo)出井筒環(huán)空溫度函數(shù)，結(jié)合井筒相關(guān)參數(shù)分析環(huán)空溫度影響因素。本文天然氣的相對氣密度取值0.55；某高溫高壓氣井井身結(jié)構(gòu)如圖2所示，結(jié)構(gòu)參數(shù)見表1。

圖3為產(chǎn)量110×104m3/d井筒溫度以及地層溫度分布。由圖可知流體溫度沿井筒上升過程中逐漸降低，環(huán)空溫度升高，因此，環(huán)空溫度的準(zhǔn)確預(yù)測對井筒完整性至關(guān)重要。

表1 典型高溫高壓井井身結(jié)構(gòu)

圖3 氣井井筒溫度分布

2.2 井筒溫度影響因素

2.2.1 焦耳–湯姆遜效應(yīng)

2.2.2 生產(chǎn)時間

生產(chǎn)初期井筒溫度主要受地層溫度的影響。圖5為井筒某深度隨生產(chǎn)時間溫度變化的曲線，0點生產(chǎn)時期為初期流體熱量未傳遞到地層的時刻。由圖可見，在生產(chǎn)時間小于30 d的時候，井筒溫度迅速上升，之后隨著生產(chǎn)時間的增加，溫度上升緩慢。井筒溫度上升主要在生產(chǎn)初期，為了保證井筒的完整性，在生產(chǎn)初期針對環(huán)空溫度要采取必要的管控措施。

圖5 井筒溫度隨生產(chǎn)時間的變化趨勢

2.2.3 產(chǎn)量

圖6 不同產(chǎn)量下井筒溫度分布

3 結(jié)論

（1）高溫、高壓氣井在生產(chǎn)過程中，井筒流體、密閉環(huán)空溫度受地層高溫流體的影響，呈現(xiàn)規(guī)律性變化。環(huán)空溫度隨產(chǎn)量增加急劇上升，到達一定產(chǎn)量時可看作環(huán)空溫度不再受產(chǎn)量的影響。生產(chǎn)時間對環(huán)空溫度的影響主要集中在前期，當(dāng)生產(chǎn)時間超過30 d后，環(huán)空溫度基本穩(wěn)定，因此，在管柱設(shè)計時期提前做好防護措施。

（2）環(huán)空溫度對高溫、高壓氣井井筒完整性的影響體現(xiàn)在環(huán)空溫度的改變導(dǎo)致環(huán)空附壓增加，為了更好地保證井筒長期有效地運行，需要對環(huán)空壓力做進一步研究。

[1] 竇益華，許愛榮，張福祥，等．高溫高壓深井試油完井問題綜述[J]．石油機械，2008，36（9）：140–142．

[2] 姜學(xué)海，王耀峰，竇益華．高溫高壓深井磨損套管應(yīng)力熱–結(jié)構(gòu)耦合場分析[J]．石油機械，2009，37（6）：19–23．

[3] RAMEY H J J. Wellbore heat transmission[J]．Journal of Petroleum Technology，1962，14（4）：427–435．

[4] Hasan A R，Kabir C S．Heat transfer during two–phase flow in wellbores [C]．SPE22 948，l991：695–708．

[5] 毛偉，梁政．計算氣井井筒溫度分布的新方法[J]．西南石油大學(xué)學(xué)報（自然科學(xué)版），1999，21（1）：56–58．

[6] 朱德武，何漢平．凝析氣井井筒溫度分布計算[J]．天然氣工業(yè)，1998，18（1）：60–62．

[7] 楊進，唐海雄，劉正禮，等．深水油氣井套管環(huán)空壓力預(yù)測模型[J]．石油勘探與開發(fā)，2013，40（5）：616–619．

[8] 張波，管志川，張琦．深水油氣井開采過程環(huán)空壓力預(yù)測與分析[J]．石油學(xué)報，2015，36（8）：1 012–1 017．

[9] HASAN A R，IZGEC B，KABIR C S．Ensuring sustained production by managing annular–pressure buildup[C]． Society of Petroleum Engineers，2009．

[10] 李相方，莊湘琦，剛濤，等．天然氣偏差系數(shù)模型綜 合評價與選用[J]．石油鉆采工藝，2001，23（2）：42–46．

[11] 毛偉，張立德．焦耳–湯姆遜系數(shù)計算方法研究[J]．特種油氣藏，2002，9（5）：44–46．

[12] 張弘，申瑞臣，梁奇敏，等．地下儲氣庫注采井溫度壓力耦合分析[J]．科學(xué)技術(shù)與工程， 2017，17（31）：66–73．

Analysis of wellbore temperature field in production process of HTHP gas wells

QIN Yanbin, ZHU Litao, ZHENG Jie, DOU Yihua

(School of Mechanical Engineering, Xi'an Petroleum University, Xi'an, Shannxi 710065, China)

Affected by the high temperature fluid in the formation during the production process of high temperature and pressure gas (HTHP) gas wells, the original balance of the wellbore temperature will be broken and the redistribution of the wellbore temperature will cause the annulus pressure to increase, which will threaten the safe service of the wellbore and the integrity of the wellbore. In order to accurately predict the temperature of the wellbore, based on mass, momentum, energy conservation, heat transfer theory and wellbore heat transfer theory, by considering the influence of gas Joule-Thomson effect, gas temperature, pressure, density and physical parameters, a wellbore temperature prediction model will be established. The fluid property parameters will be calculated according to different temperature pressures so as to improve the accuracy of the model calculation.

wellbore temperature field;high temperature and pressure gas; parameter calculation

1673–8217（2019）02–0095–06

TE827

A

2018–11–06

秦彥斌，碩士，副教授，1970年生，1992年畢業(yè)于陜西機械學(xué)院機械設(shè)計與制造專業(yè)，現(xiàn)從事石油機械研究及應(yīng)用工作。

國家自然科學(xué)基金項目“頁巖氣水平井壓裂與生產(chǎn)套管變形機理及其控制機制研究”（NSFC51674199）。

編輯：趙川喜