水平氣井積液診斷方法及攜液模型研究進展

黃 斌，張庭瑋，傅 程，付思強，賀世博，黃立凱

(1.東北石油大學 石油工程學院，黑龍江 大慶 163318； 2.大慶油田博士后科研工作站，黑龍江 大慶 163413；3.大慶油田 第二采油廠，黑龍江 大慶 163414)

在天然氣田的生產開發過程中，天然氣井的積液現象十分普遍，這種現象使得氣井產量下降甚至造成氣井不能正常生產。對于已經投產的氣藏，大概有80%的有水氣藏，中石油四大氣區出水井約有6 000口，產量占總產量的60%。因此，研究氣井積液的診斷方法和臨界攜液模型的梳理與發展趨勢意義重大。

要解決氣井積液問題，需要對氣井積液進行準確判斷，如今水平氣井積液診斷方法主要包括生產數據分析法、生產測試法、理論分析法，本文對這些方法進行整理分析，同時針對不同的井段建立了不同的氣井攜液模型，解決了不同井段的臨界攜液流量計算問題，增強了模型計算的準確性。

1 氣井積液機理

在天然氣井的生產早期，井筒內的各種流體以環狀流形式流動，這時井筒內的氣體都以較高的流速運行[1]。隨著石油資源開發年限的逐步延長，氣體攜液量逐漸下降，當氣體的流速低于臨界攜液的流速時，液體就開始快速地聚集，產量也隨之開始大幅度地降低，井筒中的流態由環狀流轉變成了段塞流。這時由于氣體的產量繼續減少，氣流不再把液體直接攜帶到井口，最終形成了井底的積液[2]。氣井積液過程如圖1所示。

圖1 氣井積液過程Fig.1 Gas well liquid accumulation process

2 氣井積液診斷方法

2.1 生產數據分析法

該方法是將產液量、產氣量等氣田常規數據與生產數據進行比較，若出現以下幾種情況：產量急劇減小、氣井出液量急劇下降、井口油壓或套壓急劇下降、油套壓差增加、井底壓力或其壓力梯度快速增加、氣井間噴等，則可判斷為積液。

油套壓差法[3]是沒有節流器的生產井進行井筒積液判斷最常用的一種方法，井筒是否積液可以通過油套壓差大小直觀判斷，并能初步估算氣井積液量。

油管內液面的深度：

(1)

套管環形空間的液面深度：

(2)

式中，Hlt為油管液面深度；Hlc為套管液面深度；ρ為液體密度；Pc1、Pc2分別為井筒積液前后套壓；Pt1、Pt2分別為井筒積液前后油壓。

對于正常氣井，當油套壓差增大，說明油管中的流動損失很大，攜液能力不足，舉升不正常，積液較多，液體不能完全被氣體攜帶出來，或者短時間內油壓和套壓急劇降低也是井筒積液的表現特征[4]。

但是該方法對于有封隔器的氣井并不適用，而且油套壓差下降的原因不止氣井積液一種，可能造成對于氣井積液情況的誤判。

2.2 生產測試法

在生產或關井狀態下，在井筒內進行壓力剖面試驗或采用其他儀器對氣液界面進行探測，以判斷井內是否有積液[5]。這種方法的優勢是能精確診斷，缺點是無法長期持續監測，不能對即將積液的氣井起到預警作用，發現時氣井已經積液。另外，該方法需要進行作業和配備有關儀器，增加了開發費用。

實測壓力梯度曲線法[6]是生產測試法的主要方法，目前現場使用的壓力測試過程比較簡單，是在關井以及在生產過程中，對深度不同位置的壓力進行測試。壓力梯度曲線隨著井身以及密度等因素會發生顯著的變化，以此來判斷氣井中的積液情況。氣體的重率計算公式如下：

(3)

式中，γ為氣體在井底狀況下的重率；γSC為氣體在地面標準狀況下的重率；P為井底壓力；PSC為地面標況下的壓力；Z為地面標況下的壓縮因子；ZSC為在井底條件下的壓縮因子；T為地面標況下的溫度；TSC為井底溫度。

式(3)對井筒內任何位置的氣體重率都可以進行計算。氣液混合的壓力越大，表明井中含水的量越大，反映到壓力梯度圖中時，梯度曲線越陡。

在判斷氣井積液的幾種方法中，實測壓力梯度曲線法是最直觀的方法。每一氣井沒有發生積液時的壓力梯度不同，不能單獨看曲線中的梯度值，也需要關心曲線的波動和拐點[7]。

2.3 理論分析計算法

理論分析計算法包括臨界流量計算法、凝析水量計算法。臨界流量計算法將在下一節詳述。凝析水量計算法適合生產初期無自由水產出的井。根據井口與井底的條件，分別計算凝析水量，再與理論計算凝析水量進行比較就可以判斷氣井是否積液[8]。

凝析水量計算公式[9]是井底條件下的天然氣飽和含水量和井口條件下天然氣飽和含水量的差值。

Qw=Qg×W

(4)

W=1.601 9×10-4A[0.32(0.056 25T+1)]BC

A=3.4+418.027 8/P

B=3.214 7+3.853 7×10-2P-4.775 2×

10-4×P2

C=1-0.489 3S-1.757S2

式中，Qw為氣井凝析水量；Qg為氣井產氣量，104m3/d；W為水氣比；P為壓力；S為氯化鈉含量；C為礦化度校正系數。

這種判斷方法的假設是地層無自由水產出且地層自然條件下天然氣完全被水蒸氣飽和，不適用于有邊/底水產出的氣井，因此該方法的適用性較差。氣井積液判斷方法分析見表1。

表1 氣井積液判斷方法分析Tab.1 Analysis of judgment method of gas well liquid-loading

3 水平井氣井攜液模型

臨界攜液模型是預測氣井積液的重要手段，因此臨界攜液模型的建立對于天然氣井的生產具有重要意義。本文對于水平井三段的臨界攜液模型分別進行了梳理與總結。

3.1 直井段

(1)液滴模型。Turner[10]最先提出液滴和液膜兩種解釋氣井積液的理論，然后比較了兩種物理模型，結果表明，液滴模型能更加準確的反應液體在井筒內的運動情況，并上調20%的安全系數。Turner模型認為井筒內隨氣流運動的液滴是圓球型的，是后來許多學者研究氣井攜液理論的基礎，已成為預測臨界攜液流量的經典算法。李閩等[11]認為井筒內氣流攜帶的液滴是橢球型的，相比于Turner模型，更加符合現場數據。根據以上研究可以發現，很多學者認為井筒內運動的液滴形狀不變，曳力系數和關系式系數都是定值，這就導致計算模型具有一定的局限性。

目前，有學者認為氣流攜帶的液滴會發生變形，曳力系數也不是一個定值，液滴變形也是攜液模型需要考慮的因素[12]。細化了臨界攜液模型的影響因素，彌補了現有的直井液滴模型曳力系數為定值的缺陷，適用條件更廣。

(2)液膜模型。也有部分學者認為直井內液相還是以液膜形式存在，潘杰等[13]基于液膜理論，認為液膜模型中也會夾帶液滴，液滴夾帶判據為臨界液膜流量和臨界氣相流速。李金潮等[14]根據液膜在不同氣速范圍內的速度分布規律，將液膜與管壁剪切應力為零時對應的氣速作為氣井積液臨界氣速，并考慮管徑、液相流速、氣芯中液滴夾帶等因素的影響，構建了適用于垂直氣井積液預測的零剪切應力模型。

(3)積液位置研究。Sutton等[15]認為臨界攜液流速的計算條件是位于井口或井底。一般情況下，如果井口壓力大于1 000 psia，則應使用井口條件來計算臨界速度。低于100 psia的低壓作業應根據井底條件來計算臨界速度。王毅忠等[16]提出最小攜液臨界流量的變化與產液氣井的類型相關，而且產水天然氣井的計算條件一般都是以井底的情況為準，而產凝析油的氣井則應以井口最小臨界攜液流量情況為準。積液位置的研究為臨界攜液模型的計算提供了確切的計算條件。

(4)可視化實驗。魏納[17]對液滴模型、液膜模型進行可視化研究。實驗結果表明，液滴與液膜在管流中是不斷交替出現的，液膜模型雖然在垂直管流中不占主導地位，但也對氣井積液影響巨大。可視化實驗能夠了解氣井積液過程中氣液兩相流的變化過程，以及液滴在管流中的形狀變化情況。

3.2 斜井段

(1)液膜模型。Belfroid等[18]考慮了傾角對于氣井積液預測的影響，將傳統Turner模型和Fiedler形狀函數[19]相結合，提出了適用于帶有傾角的井筒的臨界攜液模型，該模型在預測水平井斜井段積液方面優于Turner模型。Chen Dechun等[20]將井筒內的流體劃分為筒心的氣芯和管壁附近的液膜，建立了斜井井筒的液膜模型，分析了修正項對井斜度、油管直徑等不同井參數的敏感性。相比Belfroid模型，陳德春模型考慮了管壁對液滴的影響。利用Belfroid的數據對模型進行驗證，在與數據的匹配方面甚至優于Belfroid模型。

(2)液滴模型。Shi Juntai等[21]通過實驗研究發現，井筒內的較大液滴呈半漢堡形，在斜井井筒中，半漢堡液滴的形態用半漢堡液滴的平面與垂直于流動方向的橫截面之間的夾角來描述，它隨液滴的傾斜角度而變化。這一研究將臨界攜液流速與液滴寬度聯系到一起，適用于多裂縫水平井的積液判斷。

(3)可視化實驗。劉永輝等[22]建立了一套可視化氣水兩相流動模擬實驗裝置，利用高速攝像機觀察液膜反轉現象，與之對應的流速便是臨界攜液流速。將實驗與理論模型相結合，能夠更準確地預測氣井積液的發生，更加貼合現場實際情況。

3.3 水平段

對于水平段來說，管內氣液兩相流型一般為分層流[23]，此時的管底液膜厚度遠大于管頂。通過實驗，肖高棉等[24]發現水平段對于水平井攜液影響并不大，但在之后的研究中又發現，由于水平段流型是分層流，因此管底液膜厚度遠大于管頂液膜厚度，這會導致水平段積液，此時，界面不穩定波動在水平段是主要的影響因素。水平段對應的攜液模型有攜帶沉降模型和K-H波動模型[25]等。水平管分層流模型如圖2所示。

圖2 水平管分層流模型示意Fig.2 Diagram of horizontal pipe layered flow model

Lin等[25]利用線性穩定理論，發現K-H波動理論在非黏條件下能正確預測分層流的穩定流動，建立K-H波動不穩定發生的水平管連續攜液模型。Andritsos等[26]研究發現，液滴霧化大概是在氣速形成K-H不穩定波動的兩倍時，建立的連續攜液模型是基于K-H不穩定波動理論。

攜帶沉降機理[27]是井筒內的高速氣流將液體以液滴的形式夾帶到管頂，從而引起液膜厚度的變化，而管頂液膜又由于重力沿管壁回流至管底。管中液滴的形成和攜帶與氣液界面的波動密切相關，氣液界面波動將氣泡卷入到液膜中，而浮力會使氣泡上升并向周圍的液膜聚集。此時的液膜由于受到重力作用、氣液界面的剪切力與氣泡破裂而變得非常薄，最終導致液泡破裂，形成一些小液滴。其中一部分液滴被氣流直接帶走，從而引起液滴攜帶；另一部分液滴沉降在管壁四周形成液膜，如圖3所示。

圖3 攜帶沉降機理示意Fig.3 Carrying sedimentation mechanism diagram

4 結論與展望

(1)幾種氣井積液的診斷方法都有各自的適用性，只使用一種方法并不能完全準確地預測氣井積液，而且每種方法都有一定的局限性。

(2)攜液模型目前存在的問題是對于液滴尺寸無法準確測量，液膜模型存在的問題則是無法準確測量氣液界面的剪切力及液膜和管壁面的剪切力共同作用下的液膜所受的實際剪切力。

(3)許多對于臨界攜液模型的改進，都是以Turner模型為基礎進行改進，他們之間的差別需要進一步分析。盡管臨界攜液模型有很多，但是沒有能夠完全使用于大部分氣井的模型，模型的合理性和具體適用條件還需要進一步確定。需要將理論模型與實驗結合起來，通過實驗進一步了解攜液機理，提出適用范圍更廣的攜液模型。