張波，管志川，勝亞楠，王慶，許傳斌

（中國(guó)石油大學(xué)（華東）石油工程學(xué)院）

深水油氣井井筒內(nèi)流體特性對(duì)密閉環(huán)空壓力的影響

張波，管志川，勝亞楠，王慶，許傳斌

（中國(guó)石油大學(xué)（華東）石油工程學(xué)院）

針對(duì)密閉環(huán)空壓力威脅深水油氣井安全生產(chǎn)的問(wèn)題，建立了密閉環(huán)空壓力計(jì)算模型，研究了井筒內(nèi)流體特性對(duì)環(huán)空壓力的影響規(guī)律，并分析了各可控因素的敏感性和工程可行性。建立了基于相容性原則的體積平衡方程組和基于井筒-地層耦合傳熱的井筒溫度計(jì)算模型，實(shí)現(xiàn)多層次含液密閉環(huán)空壓力計(jì)算。環(huán)空壓力隨著環(huán)空液體的膨脹壓縮比減小而降低；降低環(huán)空飽和度能從根本上消除環(huán)空壓力，提出了極限環(huán)空飽和度計(jì)算公式；產(chǎn)出液比熱容和流量的增加使環(huán)空壓力上升；環(huán)空壓力與產(chǎn)出液井底溫度呈線性關(guān)系，產(chǎn)出液井口溫度能夠反映環(huán)空壓力大??；含水率上升使環(huán)空壓力動(dòng)態(tài)增長(zhǎng)。環(huán)空飽和度的調(diào)控敏感性遠(yuǎn)高于其他因素，降低環(huán)空液體導(dǎo)熱系數(shù)具有較高的工程可行性。研發(fā)可釋放環(huán)空液體的水下井口設(shè)備、高可壓縮性材料和井下隔熱液體可有效控制環(huán)空壓力。圖12表1參28

深水油氣井；環(huán)空壓力；流體特性；環(huán)空流體；極限環(huán)空飽和度

0 引言

深水油氣井異常環(huán)空壓力可分為流體熱膨脹引起的密閉環(huán)空壓力［1］和氣竄導(dǎo)致的持續(xù)環(huán)空壓力［2］。受固井技術(shù)及地層信息不確定性的影響［3-4］，深水油氣井形成多層次的密閉含液環(huán)空，在投產(chǎn)以后產(chǎn)生密閉環(huán)空壓力。墨西哥灣、巴西、西非和南中國(guó)海海域的深水油氣井均存在不同程度的密閉環(huán)空壓力［5-7］。與陸上或近海平臺(tái)不同，深水油氣井所采用的水下井口無(wú)法釋放套管環(huán)空壓力。BP公司Marlin油田深水井生產(chǎn)套管破裂并最終導(dǎo)致油井廢棄，環(huán)空壓力就是原因之一［8］。墨西哥灣Pompano A-31井在鉆進(jìn)過(guò)程中，套管在環(huán)空壓力作用下變形致使卡鉆［9］。此外，井筒內(nèi)壓力的變化會(huì)破壞井筒密封完整性［10］，誘發(fā)持續(xù)環(huán)空壓力。

目前關(guān)于深水油氣井密閉環(huán)空壓力的研究集中于壓力風(fēng)險(xiǎn)監(jiān)測(cè)［11-12］和計(jì)算優(yōu)化［13-15］，并對(duì)相關(guān)的防治措施開(kāi)展了實(shí)驗(yàn)研究［16-17］及應(yīng)用分析［18-19］。而井筒內(nèi)的高溫地層產(chǎn)出流體是井筒溫度改變的熱量來(lái)源，環(huán)空液體則是壓力產(chǎn)生的載體，因此研究井筒內(nèi)流體特性對(duì)環(huán)空壓力的影響具有重要意義。本文建立基于相容性原則與多環(huán)空相互作用的體積平衡方程組和以井筒-地層耦合傳熱為基礎(chǔ)的井筒溫度場(chǎng)計(jì)算模型，引入比變異系數(shù)對(duì)井筒內(nèi)流體特性的影響規(guī)律和可控因素的敏感性及工程可行性進(jìn)行評(píng)價(jià)。

1 深水油氣井密閉環(huán)空壓力計(jì)算模型

1.1 多層次含液密閉環(huán)空壓力

深水油氣井在固井后形成多層次的密閉環(huán)空（見(jiàn)圖1），環(huán)空內(nèi)充滿殘留的鉆井液或完井液。投產(chǎn)后液體與套管之間的熱物性差異導(dǎo)致環(huán)空難以容納受熱膨脹的液體，進(jìn)而環(huán)空壓力上升以壓縮液體。根據(jù)體積相容性原則，溫度、壓力作用下的環(huán)空液體與密閉環(huán)空體積始終保持相等。據(jù)此可建立多環(huán)空條件下的體積平衡方程組：

圖1 深水油氣井結(jié)構(gòu)示意圖

通過(guò)該方程組可求得各環(huán)空的壓力值，其中環(huán)空體積變化是一個(gè)關(guān)于溫度、壓力的函數(shù)。考慮到相鄰環(huán)空的影響，溫度和壓力導(dǎo)致的環(huán)空徑向變化可分別表示為［20］：

整個(gè)環(huán)空的體積變化為：

1.2 投產(chǎn)后環(huán)空溫度變化

半穩(wěn)態(tài)方法在防蠟［21］、氣驅(qū)采油［22］和超臨界二氧化碳［23］等領(lǐng)域有著廣泛的應(yīng)用，計(jì)算簡(jiǎn)便且能保持較好的精度［19］。因此，基于半穩(wěn)態(tài)方法取如圖1所示的控制體，則進(jìn)出微元體的能量符合守恒定律：

（5）式中壓降由重力壓力降和摩阻壓力降組成：

（6）式中f為無(wú)因次水力摩阻系數(shù)［24］：

井筒內(nèi)以及井筒和地層之間的傳熱過(guò)程均符合徑向熱流守恒原理，則可得徑向熱流量為：

（8）式中，無(wú)因次地層溫度TD的計(jì)算公式主要有Ramey公式、Butler公式、Chiu公式和Hassan公式等。研究表明Hassan公式在整個(gè)生產(chǎn)時(shí)間內(nèi)的求解精度較高［25］，該公式在井筒溫度計(jì)算［26］和密閉環(huán)空壓力分析［27］中也取得了較好的應(yīng)用效果，因此采用Hassan公式計(jì)算無(wú)因次地層溫度：

聯(lián)立（5）式、（6）式和（8）式，可得環(huán)空溫度：

其中

則環(huán)空平均溫度變化值為：

2 環(huán)空液體的影響

西非海域某深水油井產(chǎn)出物為油水混合物，井身結(jié)構(gòu)如圖1所示，相關(guān)參數(shù)見(jiàn)表1。油套環(huán)空（環(huán)空1）能夠進(jìn)行放壓，因此用生產(chǎn)套管與技術(shù)套管間的環(huán)空（環(huán)空2）分析井筒內(nèi)流體對(duì)環(huán)空壓力的影響。

表1 計(jì)算參數(shù)

2.1 環(huán)空液體等壓膨脹系數(shù)和等溫壓縮系數(shù)

圖2顯示，環(huán)空壓力隨著液體等壓膨脹系數(shù)的降低和等溫壓縮系數(shù)的提高而減小。定義環(huán)空液體膨脹壓縮比：

由（12）式可知，EC值越小，環(huán)空壓力越小。而從工程角度來(lái)講，環(huán)空液體等溫壓縮系數(shù)是可控的。

圖2 液體膨脹性與可壓縮性對(duì)環(huán)空壓力的影響

圖3顯示，不同生產(chǎn)時(shí)間和工況下環(huán)空壓力均隨著液體等溫壓縮系數(shù)的增加而降低，但是不同范圍的敏感程度不同。以圖3中基礎(chǔ)數(shù)據(jù)曲線為例，環(huán)空液體等溫壓縮系數(shù)在（3～6）×10-4MPa-1變化時(shí)，調(diào)控效果明顯；當(dāng)?shù)葴貕嚎s系數(shù)超過(guò)8×10-4MPa-1以后，調(diào)控效果不佳。

圖3 環(huán)空壓力與液體等溫壓縮系數(shù)的關(guān)系曲線

2.2 環(huán)空液體導(dǎo)熱系數(shù)

圖4是環(huán)空壓力隨油套環(huán)空（環(huán)空A）液體導(dǎo)熱系數(shù)的變化曲線，可以看出：①環(huán)空壓力隨環(huán)空液體導(dǎo)熱系數(shù)降低而減小。這是由于環(huán)空液體導(dǎo)熱系數(shù)的降低增加了井筒徑向傳熱熱阻，減緩了環(huán)空溫度上升速度。②當(dāng)環(huán)空液體導(dǎo)熱系數(shù)大于0.60 W/（m·K）時(shí)，不同情況下曲線差異較大，最高差值可達(dá)40.19 MPa，其中產(chǎn)液量的影響尤為顯著；當(dāng)導(dǎo)熱系數(shù)在0.20 W/（m·K）以下時(shí)，差異顯著減小，在0.05 W/（m·K）處基本處于同一數(shù)值范圍。較小的環(huán)空液體導(dǎo)熱系數(shù)對(duì)環(huán)空壓力具有較好的調(diào)控效果。

圖4 環(huán)空壓力與環(huán)空液體導(dǎo)熱系數(shù)的關(guān)系曲線

2.3 環(huán)空飽和度

環(huán)空飽和度即環(huán)空液體與環(huán)空的體積比。圖5顯示，環(huán)空壓力在環(huán)空飽和度小于一定值時(shí)為零，之后呈線性增長(zhǎng)關(guān)系。以圖5中基礎(chǔ)數(shù)據(jù)曲線為例，零點(diǎn)位置為0.972，之后環(huán)空壓力從零增長(zhǎng)至57.32 MPa。這是因?yàn)檩^低的環(huán)空飽和度能夠緩解環(huán)空液體與環(huán)空的體積矛盾，容納發(fā)生熱膨脹的液體。圖5還顯示，零點(diǎn)位置隨著生產(chǎn)時(shí)間和產(chǎn)液量的增大而左移。

由圖5可知，調(diào)整環(huán)空飽和度能夠有效控制甚至消除環(huán)空壓力。定義任何工況下環(huán)空壓力始終保持為零的環(huán)空飽和度為極限環(huán)空飽和度，此時(shí)環(huán)空體積變化只與溫度有關(guān)。由于井筒的最高溫度等于油藏溫度，極限環(huán)空飽和度可由下式求取：

圖5 環(huán)空壓力與環(huán)空飽和度的關(guān)系曲線

3 地層產(chǎn)出液的影響

3.1 產(chǎn)出液熱容流率

熱容流率為產(chǎn)出液比熱容與質(zhì)量流速的乘積。圖6是環(huán)空壓力隨產(chǎn)出液比熱容和日產(chǎn)液量的變化趨勢(shì)，可知環(huán)空壓力隨熱容流率的增大而增大。這是因?yàn)閺较騻鳠崴俣入S著熱容流率增大而加快。

圖6 產(chǎn)出液熱容流率對(duì)環(huán)空壓力的影響

產(chǎn)出液比熱容取決于油藏液體的組分及性質(zhì)，而日產(chǎn)液量是可控的。圖7顯示了日產(chǎn)液量對(duì)環(huán)空壓力的影響，可以看出環(huán)空壓力對(duì)不同范圍的日產(chǎn)液量敏感度不同。以圖7中基本數(shù)據(jù)曲線為例，產(chǎn)液量由50 t/d增加到200 t/d，環(huán)空壓力增加29 MPa；而當(dāng)產(chǎn)液量從500 t/d增長(zhǎng)到800 t/d時(shí)，環(huán)空壓力僅增加3.64 MPa。這表明，較大的日產(chǎn)液量削弱了比熱容對(duì)環(huán)空壓力的影響。

圖7 環(huán)空壓力與日產(chǎn)液量的關(guān)系曲線

3.2 產(chǎn)出液溫度

產(chǎn)出液是井筒溫度場(chǎng)再分布的熱量來(lái)源。由圖8可知，環(huán)空壓力隨產(chǎn)出液井底溫度呈線性增加，不同條件下的曲線僅斜率發(fā)生輕微變化。因此，深水高溫油氣井要重視密閉環(huán)空壓力帶來(lái)的風(fēng)險(xiǎn)。

圖8 環(huán)空壓力與產(chǎn)出液井底溫度的關(guān)系曲線

生產(chǎn)實(shí)踐證明環(huán)空壓力與產(chǎn)出液井口溫度之間存在相關(guān)性［19］。產(chǎn)出液井口溫度隨著生產(chǎn)的進(jìn)行不斷變化，圖9所對(duì)應(yīng)的生產(chǎn)時(shí)間為10～800 d，可以看出環(huán)空壓力與產(chǎn)出液井口溫度近似呈線性關(guān)系。因此，可以通過(guò)監(jiān)測(cè)井口產(chǎn)出液溫度來(lái)估算環(huán)空壓力的大小。

3.3 產(chǎn)出液含水率

實(shí)際生產(chǎn)中，油井產(chǎn)出液含水率不斷上升，最終到達(dá)相對(duì)平穩(wěn)狀態(tài)。含水率的上升意味著產(chǎn)出液比熱容和日產(chǎn)液量的增加。如圖10所示，某海上油田投產(chǎn)以后含水率由零上升至96.23%，平均單井產(chǎn)液量由初期的500 m3/d上升至1 000 m3/d以上。結(jié)合3.1中的分析可知，含水率的上升最終會(huì)導(dǎo)致環(huán)空壓力的增加。這說(shuō)明，深水油氣井的環(huán)空壓力上升過(guò)程是一個(gè)動(dòng)態(tài)過(guò)程，在實(shí)際的預(yù)測(cè)中要考慮產(chǎn)出液含水率的變化，提高分析評(píng)價(jià)的可靠性。

圖9 環(huán)空壓力與產(chǎn)出液井口溫度的關(guān)系曲線

圖10 某海上油田生產(chǎn)情況［28］

4 環(huán)空壓力調(diào)控敏感性與可行性

環(huán)空壓力調(diào)控敏感性為某一影響因素變化一定數(shù)值時(shí)所對(duì)應(yīng)的環(huán)空壓力在一定范圍內(nèi)的波動(dòng)程度。標(biāo)準(zhǔn)差系數(shù)能夠反映一組數(shù)據(jù)的變動(dòng)程度，同時(shí)消除單位差異所帶來(lái)的影響。因此引入影響因素與環(huán)空壓力的標(biāo)準(zhǔn)差系數(shù)比值來(lái)表征某一因素敏感性，稱為比變異系數(shù)：

按照敏感性定義和比變異系數(shù)計(jì)算方法，首先對(duì)各個(gè)因素取值進(jìn)行歸零化：

環(huán)空液體膨脹性、產(chǎn)出液含水率、溫度及比熱容均為非可控因素，因此不再對(duì)這些因素進(jìn)行評(píng)價(jià)。其余可控因素的敏感性如圖11所示，敏感性越大意味著調(diào)控效果越好。由圖11可知，環(huán)空飽和度、環(huán)空液體導(dǎo)熱系數(shù)、環(huán)空液體等溫壓縮系數(shù)和日產(chǎn)液量的敏感性依次降低，其中環(huán)空飽和度的敏感性遠(yuǎn)大于后3個(gè)因素，后3個(gè)因素在不同范圍的敏感性差異性較大。

圖11 各影響因素敏感性

圖12定性地描述了各調(diào)控措施的工程可行性。具體分析如下。

①安裝在套管外側(cè)的可壓縮物質(zhì)可在一定程度上降低環(huán)空飽和度，但只能保護(hù)單一環(huán)空且不能從根本上消除環(huán)空壓力。

②環(huán)空中注入多元醇、可溶性鹽等能有效降低液體導(dǎo)熱系數(shù)。同理，安裝隔熱油管或套管也是可行的。

③在環(huán)空液體中混入氮?dú)猓?8］或可壓縮的玻璃微球能夠調(diào)節(jié)其壓縮性。但深水油氣井的水深一般超過(guò)500 m，超深水可達(dá)1 000 m以上，氮?dú)饣虿Ａ⑶蛟谧⑷脒^(guò)程中首先受到隔水管內(nèi)液柱壓力的作用，造成氮?dú)獾目蓧嚎s性大幅降低或玻璃微球提前破裂，且調(diào)控效果隨注入量的增加而降低。

④深水油氣田的成本回收依賴于高產(chǎn)量，因此無(wú)論是從經(jīng)濟(jì)效益還是調(diào)控效果角度，降低日產(chǎn)液量均不具有可行性。

圖12 調(diào)控措施工程可行性與調(diào)控效果

5 結(jié)論

環(huán)空壓力隨著環(huán)空液體膨脹壓縮比的減小而降低，但變化趨勢(shì)逐漸放緩。降低環(huán)空液體的導(dǎo)熱系數(shù)能夠有效調(diào)控環(huán)空壓力，調(diào)控效果隨著導(dǎo)熱系數(shù)的減小而增強(qiáng)。降低環(huán)空飽和度能從根本上消除環(huán)空壓力。產(chǎn)出液熱容流率的增加會(huì)導(dǎo)致環(huán)空壓力上升，效果隨日產(chǎn)液量增加而減弱。環(huán)空壓力與產(chǎn)出液溫度呈線性關(guān)系，利用產(chǎn)出液井口溫度可對(duì)環(huán)空壓力進(jìn)行估算。含水率升高致使環(huán)空壓力動(dòng)態(tài)上升。

環(huán)空液體壓縮性及導(dǎo)熱系數(shù)、環(huán)空飽和度和日產(chǎn)液量作為可控因素，其敏感性和工程可行性具有差異。環(huán)空飽和度的敏感性遠(yuǎn)高于其他因素。環(huán)空液體導(dǎo)熱系數(shù)、等溫壓縮系數(shù)和日產(chǎn)液量的敏感性依次降低，但處于同一數(shù)量級(jí)。從工程角度來(lái)看，環(huán)空液體的導(dǎo)熱系數(shù)具有較高的調(diào)控可行性，環(huán)空液體可壓縮性次之。降低環(huán)空飽和度的工程可行性較低，降低日產(chǎn)液量基本不具有工程可行性。

為降低環(huán)空壓力所帶來(lái)的風(fēng)險(xiǎn)，應(yīng)研發(fā)具有釋放套管環(huán)空液體能力的水下井口設(shè)備，實(shí)現(xiàn)對(duì)環(huán)空飽和度的調(diào)節(jié)。同時(shí)研制高可壓縮性材料和新型隔熱材料，進(jìn)而合理配置環(huán)空液體的可壓縮性和導(dǎo)熱性。

符號(hào)注釋：

a——管柱內(nèi)徑，m；b——管柱外徑，m；C——由井深結(jié)構(gòu)確定的待定系數(shù)；Cf——流體比熱容，J/（kg·K）；Cr——比變異系數(shù)，f；dto——油管直徑，m；E——套管的彈性模量，MPa；EC——環(huán)空液體膨脹壓縮比，MPa/K；f——無(wú)因次水力摩阻系數(shù)；g——重力加速度，取9.8 m/s2；gf——地溫梯度，K/m；h——井深，m；i——環(huán)空編號(hào)；kT——液體等溫壓縮系數(shù)，MPa-1；n——因素取值個(gè)數(shù)；p——壓力，MPa；pan——環(huán)空壓力，MPa；pi——管柱內(nèi)壓，MPa；po——管柱外壓，MPa；pzi——進(jìn)入微元體的流體壓力，Pa；pzo——流出微元體的流體壓力，Pa；pana——環(huán)空壓力的平均值，MPa；Ra——油管內(nèi)壁的粗糙度，m；Ran——環(huán)空熱阻，m·K/W；Re——雷諾數(shù)；Rto——徑向傳熱總熱阻，m·K/W；r——計(jì)算點(diǎn)處的半徑，m；rw——井筒半徑，m；San——環(huán)空極限飽和度，f；T——溫度，K；T0——泥線溫度，K；TD——無(wú)因次地層溫度；Te——地層溫度，K；Tf——油管內(nèi)流體溫度，K；t——生產(chǎn)時(shí)間，s；tD——無(wú)因次生產(chǎn)時(shí)間，tD=tαe/rw2；Tan——環(huán)空平均溫度，K；Tan0——環(huán)空初始平均溫度，K；Van——環(huán)空體積，m3；Vf——環(huán)空液體初始體積，m3；vf——油管內(nèi)流體流速，m/s；Wf——流量，kg/s；xifa——因素變化數(shù)值的平均值；xj——因素x的第j個(gè)取值；xj′——因素x第j個(gè)取值的歸零化結(jié)果；z——計(jì)算點(diǎn)處坐標(biāo)，m；zi——環(huán)空底部位置，m；zo——環(huán)空頂部位置，m；α——液體等壓膨脹系數(shù)，K-1；αe——地層熱擴(kuò)散系數(shù)，m2/s；αs——鋼材的線性膨脹系數(shù)，K-1；Δpan——環(huán)空壓力變化，MPa；Δrp——壓力引起的管徑變化，m；Δrpi——壓力引起的內(nèi)徑變化，m；Δrpo——壓力引起的內(nèi)徑變化，m；Δrt——溫度引起的管徑變化，m；Δrti——溫度引起的內(nèi)徑變化，m；Δrto——溫度引起的內(nèi)徑變化，m；ΔT——環(huán)空溫度變化，K；ΔTd——計(jì)算點(diǎn)的溫度變化，K；ΔVan——環(huán)空體積變化，m3；θ——井斜角，（°）；λe——地層導(dǎo)熱系數(shù)，W/（m·K）；μ——套管泊松比；ρf——流體密度，kg/m3；σp——環(huán)空壓力的標(biāo)準(zhǔn)差，MPa；σx——因素變化數(shù)值的標(biāo)準(zhǔn)差；Φki——進(jìn)入微元體的動(dòng)能，J/s；Φko——流出微元體的動(dòng)能，J/s；Φr——徑向熱流量，J/s。

［1］楊進(jìn)，唐海雄，劉正禮，等. 深水油氣井套管環(huán)空壓力預(yù)測(cè)模型［J］. 石油勘探與開(kāi)發(fā)，2013，40（5）: 616-619. YANG Jin，TANG Haixiong，LIU Zhengli，et al. Prediction model of casing annulus pressure for deepwater well drilling and completion operation［J］. Petroleum Exploration and Development，2013，40（5）: 616-619.

［2］張波，管志川，張琦，等. 高壓氣井環(huán)空壓力預(yù)測(cè)與控制措施［J］.石油勘探與開(kāi)發(fā)，2015，42（4）: 518-522. ZHANG Bo，GUAN Zhichuan，ZHANG Qi，et al. Prediction of sustained annular pressure and the pressure control measures for high pressure gas wells［J］. Petroleum Exploration and Development，2015，42（4）: 518-522.

［3］GUPTA N，BOGAERTS M，BELLABARBA M，et al. Overcoming ultradeepwater cementing challenges in the Caribbean［R］. SPE 169245，2014.

［4］許玉強(qiáng)，管志川，張洪寧，等. 基于含可信度地層壓力剖面的精細(xì)井身結(jié)構(gòu)設(shè)計(jì)方法［J］. 中國(guó)石油大學(xué)學(xué)報(bào)（自然科學(xué)版），2014，38（6）: 72-77. XU Yuqiang，GUAN Zhichuan，ZHANG Hongning，et al. Methods for precise casing design based formation profiles with different credibilities［J］. Journal of China University of Petroleum （Edition of Natural Science），2014，38（6）: 72-77.

［5］BELLARBY J，KOFOED S S，MARKETZ F. Annular pressure build-up analysis and methodology with examples from multifrachorizontal wells and HPHT reservoirs［R］. SPE 163557，2013.

［6］FERREIRA M V D，DOS SANTOS A R，VANZAN V. Thermally insulated tubing application to prevent annular pressure build in Brazil offshore fields［R］. SPE 151044，2012.

［7］黃小龍，嚴(yán)德，田瑞瑞，等. 深水套管環(huán)空圈閉壓力計(jì)算及控制技術(shù)分析［J］. 中國(guó)海上油氣，2014，26（6）: 61-65. HUANG Xiaolong，YAN De，TIAN Ruirui，et al. Research of calculation of casing annulus trapped pressure and its control techniques for deepwater well［J］. China Offshore Oil and Gas，2014，26（6）: 61-65.

［8］BRADFORD D W，F(xiàn)RITCHIE D G，GIBSON D H，et al. Marlin failure analysis and redesign: Part 1: Description of failure［J］. SPE Drilling & Completion，2004，19（2）: 101-104.

［9］PATTILLO P D，COCALES B W，MOREY S C. Analysis of an annular pressure buildup failure during drill ahead［J］. SPE Drilling & Completion，2006，21（4）: 242-247.

［10］初緯，沈吉云，楊云飛，等. 連續(xù)變化內(nèi)壓下套管-水泥環(huán)-圍巖組合體微環(huán)隙計(jì)算［J］. 石油勘探與開(kāi)發(fā)，2015，42（3）: 379-385. CHU Wei，SHEN Jiyun，YANG Yunfei，et al. Calculation of micro-annulus size in casing-cement sheath-formation system under continuous internal casing pressure change［J］. Petroleum Exploration and Development，2015，42（3）: 379-385.

［11］SULTAN N W，F(xiàn)AGET J B P，F(xiàn)JELDHEIM M，et al. Real-time casing annulus pressure monitoring in a subsea HP/HT exploration well［R］. OTC 19286，2008.

［12］JANDHYALA S R K，CHINEY A. Finite element approach to predict the effect of annular pressure buildup on wellbore materials［R］. OTC 24837，2014.

［13］張波，管志川，張琦. 深水油氣井開(kāi)采過(guò)程環(huán)空壓力預(yù)測(cè)與分析［J］. 石油學(xué)報(bào)，2015，33（8）: 1012-1017. ZHANG Bo，GUAN Zhichuan，ZHANG Qi. Prediction and analysis on annular pressure of deepwater well in the production stage［J］. Acta Petrolei Sinica，2015，33（8）: 1012-1017.

［14］YIN Feng，GAO Deli. Improved calculation of multiple annuli pressure buildup in subsea HPHT wells［R］. SPE 170553，2014.

［15］張百靈，楊進(jìn)，黃小龍，等. 深水井筒環(huán)空壓力計(jì)算模型適應(yīng)性評(píng)價(jià)［J］. 石油鉆采工藝，2015，37（1）: 56-59. ZHANG Bailing，YANG Jin，HUANG Xiaolong，et al. Adaptability evaluation of calculation model of annular pressure of deepwater wellhole［J］. Oil Drilling & Production Technology，2015，37（1）: 56-59.

［16］OSGOUEI R E，MISKA S Z，OZBAYOGLU M E，et al. Annular pressure build up （APB） analysis: Optimization of fluid rheology［R］. SPE 170262，2014.

［17］周波，楊進(jìn)，劉正禮，等. 深水油氣井筒環(huán)空注氮控壓機(jī)理［J］. 石油勘探與開(kāi)發(fā)，2015，42（3）: 386-389. ZHOU Bo，YANG Jin，LIU Zhengli，et al. Mechanism of pressure management by injecting nitrogen in casing annulus of deepwater wells［J］. Petroleum Exploration and Development，2015，42（3）: 386-389.

［18］艾爽，程林松，劉紅君，等. 深層高溫高壓氣井臨界產(chǎn)量計(jì)算模型［J］. 計(jì)算物理，2015，32（3）: 327-333. AI Shuang，CHENG Linsong，LIU Hongjun，et al. A critical production rate model for deep high temperature high gas wells［J］. Chinese Journal of Computational Physics，2015，32（3）: 327-333.

［19］HASAN R，IZGEC B，KABIR S. Sustaining production by managing annular-pressure buildup［J］. SPE Operation & Production，2010，25（2）: 195-203.

［20］楊海波，曹建國(guó)，李洪波. 彈性與塑性力學(xué)簡(jiǎn)明教程［M］. 北京:清華大學(xué)出版社，2011. YANG Haibo，CAO Jianguo，LI Hongbo. Introduction to elasticity and plasticity［M］. Beijing: Tsinghua University Press，2011.

［21］宋洵成，管志川，韋龍貴，等. 保溫油管海洋采油井筒溫度壓力計(jì)算耦合模型［J］. 石油學(xué)報(bào)，2012，33（6）: 1064-1067. SONG Xuncheng，GUAN Zhichuan，WEI Longgui，et al. A coupled model of temperature-pressure calculation for offshore production wellbores with insulated tubing［J］. Acta Petrolei Sinica，2012，33（6）: 1064-1067.

［22］王瑞和，倪紅堅(jiān). 二氧化碳連續(xù)管井筒流動(dòng)傳熱規(guī)律研究［J］. 中國(guó)石油大學(xué)學(xué)報(bào)（自然科學(xué)版），2013，37（5）: 65-70. WANG Ruihe，NI Hongjian. Wellbore heat transfer law of carbon dioxide coiled tubing drilling［J］. Journal of China University of Petroleum （Edition of Natural Science），2013，37（5）: 65-70.

［23］李兆敏，張丁涌，衣懷峰，等. 多元熱流體在井筒中的流動(dòng)與傳熱規(guī)律［J］. 中國(guó)石油大學(xué)學(xué)報(bào)（自然科學(xué)版），2012，36（6）: 79-83，88. LI Zhaomin，ZHANG Dingyong，YI Huaifeng，et al. Flow and heat transfer regulation of muti-thermal fluids injection in wellbore［J］. Journal of China University of Petroleum （Edition of Natural Science），2012，36（6）: 79-83，88.

［24］CHEN A H. An explicit equation for （calculating the） friction factor in a pipe［J］. Industrial and Engineering Chemistry Research Fundamentals，1979，18（3）: 296-297.

［25］王彌康. 注蒸汽井井筒熱傳遞的定量計(jì)算［J］. 石油大學(xué)學(xué)報(bào)（自然科學(xué)版），1994，18（4）: 77-82. WANG Mikang. Quantitative calculation of wellbore heat transmission for steam injection wells［J］. Journal of the University of Petroleum （Edition of Natural Science），1994，18（4）: 77-82.

［26］HASAN R，KABIR C S. Heat transfer during two-phase flow in wellbores: Part I: Formation temperature［R］. SPE 22866，1991.

［27］GUAN Zhichuan，ZHANG Bo，WANG Qing，et al. Design of thermalinsulated pipes applied in deepwater well to mitigate annular pressure build-up［J］. Applied Thermal Engineering，2016，9: 129-136.

［28］鄒信波，羅東紅，許慶華，等. 海上特高含水老油田挖潛策略與措施: 以珠江口盆地陸豐凹陷LFD13-1油田為例［J］. 中國(guó)海上油氣，2012，24（6）: 28-33，37. ZOU Xinbo，LUO Donghong，XU Qinghua，et al. Some strategical measures to tap the potential in offshore mature oilfields with extra-high water cut: A case of LFD13-1 field in Lufeng depression，Pearl River Mouth basin［J］. China Offshore Oil and Gas，2012，24（6）: 28-33，37.

（編輯 胡葦瑋）

Impact of wellbore fluid properties on trapped annular pressure in deepwater wells

ZHANG Bo，GUAN Zhichuan，SHENG Yanan，WANG Qing，XU Chuanbin
（College of Petroleum Engineering，China University of Petroleum，Qingdao 266580，China）

To reduce the threat of trapped annular pressure on deepwater wells safe production，this study established a model of calculating trapped annular pressure，examined the influence of wellbore fluid properties on trapped annular pressure，and analyzed the sensitivities and engineering feasibilities of controllable factors. To realize the calculation of trapped annular pressure under multiple annuli with liquid，a volume balance matrix was built according to compatibility principle and a wellbore temperature computing model was built based on wellbore-formation coupled heat transfer. Annular pressure decreases as the expansion-compression ratio of annular fluid reduces. Decreasing annular saturation can eliminate annular pressure radically and then a formula was proposed to give extreme annular saturation. The increase of production fluid specific heat capacity and flow rate leads to enhancement of annular pressure. Annular pressure keeps a linear relation to production fluid hole bottom temperature and the wellhead temperature can reflect the value of annular pressure. The water ratio increase of production fluid causes dynamic increase of annular pressure. The sensitivity of annular saturation is much higher than other factors. Decreasing annular liquid thermal conductivity has relatively higher engineering feasibility. The annular pressure can be controlled effectively by developing subsea wellheads with the ability to release annular fluid，highly compressible materials and downhole thermal-insulated fluids.

deepwater wells； annular pressure； fluid properties； annular fluid； extreme annular saturation

國(guó)家自然科學(xué)基金“隔水管內(nèi)氣液兩相流聲傳播特性與深水鉆井氣侵檢測(cè)方法”（51574275）；教育部長(zhǎng)江學(xué)者創(chuàng)新團(tuán)隊(duì)項(xiàng)目“海洋油氣鉆完井理論與工程”（IRT1086）；中國(guó)石油大學(xué)研究生自主創(chuàng)新工程（YCXJ2016011）

TE21

A

1000-0747（2016）05-0799-07

10.11698/PED.2016.05.17

張波（1990-），男，山東濟(jì)南人，中國(guó)石油大學(xué)（華東）在讀博士研究生，主要從事油氣井流體力學(xué)和深水油氣井方面的研究工作。地址：山東省青島市經(jīng)濟(jì)技術(shù)開(kāi)發(fā)區(qū)長(zhǎng)江西路66號(hào)，中國(guó)石油大學(xué)（華東）石油工程學(xué)院，郵政編碼：266580。E-mail: zhangboupc@126.com

聯(lián)系作者：管志川（1959-），男，山東單縣人，博士，中國(guó)石油大學(xué)（華東）教授，主要從事油氣井力學(xué)、井下測(cè)控技術(shù)、深井超深井鉆井和深水鉆井等方面的教學(xué)與研究工作。地址：山東省青島市經(jīng)濟(jì)技術(shù)開(kāi)發(fā)區(qū)長(zhǎng)江西路66號(hào)，中國(guó)石油大學(xué)（華東）石油工程學(xué)院，郵政編碼：266580。E-mail: guanzhch@upc.edu.cn

2015-11-30

2016-04-21