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

?

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

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

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

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

井筒溫度場(chǎng)；高溫高壓氣體；參數(shù)計(jì)算

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

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

1 井筒溫度分布模型

1.1 流體溫度模型

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

圖1 井筒微元示意圖

將（1）式化簡(jiǎn)為：

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

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

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

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

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

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

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

流量公式：

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

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

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

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

在上述公式中：

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

1.3.1 定壓比熱采用文獻(xiàn)[6]給出的模型

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

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

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

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

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

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

2 氣井井筒溫度分布實(shí)例分析

2.1 實(shí)例計(jì)算

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

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

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

圖3 氣井井筒溫度分布

2.2 井筒溫度影響因素

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

2.2.2 生產(chǎn)時(shí)間

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

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

2.2.3 產(chǎn)量

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

3 結(jié)論

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

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

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

[2] 姜學(xué)海，王耀峰，竇益華．高溫高壓深井磨損套管應(yīng)力熱–結(jié)構(gòu)耦合場(chǎng)分析[J]．石油機(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ì)算氣井井筒溫度分布的新方法[J]．西南石油大學(xué)學(xué)報(bào)（自然科學(xué)版），1999，21（1）：56–58．

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

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

[8] 張波，管志川，張琦．深水油氣井開采過(guò)程環(huán)空壓力預(yù)測(cè)與分析[J]．石油學(xué)報(bào)，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ù)模型綜 合評(píng)價(jià)與選用[J]．石油鉆采工藝，2001，23（2）：42–46．

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

[12] 張弘，申瑞臣，梁奇敏，等．地下儲(chǔ)氣庫(kù)注采井溫度壓力耦合分析[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è)于陜西機(jī)械學(xué)院機(jī)械設(shè)計(jì)與制造專業(yè)，現(xiàn)從事石油機(jī)械研究及應(yīng)用工作。

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

編輯：趙川喜