基于液膜臨界速度的井筒積液預測新方法*

豆寧輝 柯 珂 鮑洪志 付瑋琪

(1. 中國石油大學(華東)石油工程學院 山東青島 266580； 2. 中國石化石油工程技術研究院 北京 100101；3. 中國礦業大學 煤炭資源與安全開采國家重點實驗室 江蘇徐州 221116)

在天然氣開發過程中，隨著氣田生產時間的增加，天然氣產量隨之遞減，而氣田的含水率逐漸上升。當地層產出水無法被天然氣連續帶出井筒時，地層水將在井底形成積液[1]。井筒積液量的增加會增加井底壓力，減少地層與井筒間的壓差，從而降低天然氣產量，嚴重時井底積液甚至會將氣井壓死[2]。在海上氣田中要增加開采效率，氣井多以定向井和水平井的形式存在，一旦出現積液造成氣井停產，危害極大，修井成本極高[3]。因此，預測氣井積液的發生可以幫助現場工程師準確判斷氣井生產狀況，及時采取攜液措施可以有效避免積液發生。

1969年，Turner等[4]建立了基于氣液受力模型的井底積液預測模型，該模型認為液滴在垂直井筒中下落速度受拖曳力和重力的影響，并提出了最小氣體流速的概念，當井筒內氣體流速大于最小氣體流速，井筒內無積液，反之，則存在積液。通過對比現場數據發現，在現場應用中所建立模型預測的結果需要上浮20%。Coleman等[5]運用Exxon公司的現場數據校驗了Turner模型，通過對比發現，井深直徑和井口壓力會嚴重影響液滴的臨界速度，但氣體密度、井筒溫度和氣液間界面張力對臨界速度的影響不大。Guo等[6]在Turner模型的基礎上，提出了液滴動能理論，認為使液滴上升的最小氣體流速所對應的是氣相流體的最小動能。Guo等研究認為即使turner模型所預測的結果上浮20%,該模型結果仍然低估了井底積液的狀況。Zhou等[7]基于turner模型建立了液滴濃度模型，臨界氣體流速會影響井筒內的液滴濃度和持液率，該模型運用井筒持液率的概念對井底積液情況進行判定。Veeken等[8]指出turner模型的預測結果比現場數據要低40%，同時實驗發現，井底積液不能全歸咎于井筒內液滴的下落，在環霧流中液膜的逆向流動更容易造成井底積液的發生。目前有關井筒積液研究多基于turner模型，通過判定垂直井筒中液滴的下落速度從而決定井底的積液狀況，但是針對定向井和基于液膜逆流理論的研究還不多。因此，筆者在前人實驗數據的基礎上，通過數值模擬和理論分析的方法，對定向井井底積液預測理論進行深入研究，對比液膜逆流法和最小壓降法的預測精度，以期為深水氣井積液預測提供新思路。

1 多相流動模型

本文基于井筒多相流動模型分別對垂直井筒、傾角75°和60°井筒進行積液預測。Zhang等[9-10]基于段塞流動力學建立了全角度全流型的氣液兩相流模型，統稱Unified模型。在Unified模型中當液膜長度無限大時，段塞流即會轉變為分層流或環霧流，當段塞長度無限大時，段塞流即會轉變為泡狀流，因此該模型可以用于預測環霧流中的氣液流動狀態。

段塞流由多個單位段塞組成，單位段塞由液膜部分和段塞部分構成，其中液膜部分由液膜和氣核部分組成，如圖1所示。在段塞流中，段塞向前部方向運動時，段塞部分液體流向下一段液膜，而損失的液體則由前一部分液膜所補充。同理，液膜部分損失的液體由前一部分的段塞所補充。因此，可以建立段塞流的連續性方程，液相、氣相連續性方程如式(1)、(2)。

圖1 段塞流示意圖[9]

lUvSL=lSHLSvS+lF(HLFvF+HLCvC)

(1)

lUvSG=lS(1-HLS)vS+lF(1-HLF-HLC)vC

(2)

式(1)、(2)中：lU為單位段塞長度，m；vSL為表觀液體流速，m/s；vSG為表觀氣體流速，m/s；HLS為段塞部分持液率，無量綱；HLF為液膜部分持液率，無量綱；HLC為氣核持液率，無量綱；vT為段塞過渡速度，m/s；vS為段塞速度，m/s；vF為液膜速度，m/s；vC為氣核速度，m/s；lS為段塞長度，m；lF為液膜長度，m。

當液膜長度lF為正無窮(即段塞長度lS為0)，同時氣核與管壁無接觸，段塞流則發展為環霧流，環霧流的動量守恒方程如式(3)。

(ρL-ρC)gsinθ=0

(3)

式(3)中：τF為液膜的剪切應力，Pa；SF為液膜周長，m；A為氣液間表面積，m2；τI為氣液界面的剪切應力，Pa；SI為氣液界面周長，m；ρL為液體密度，kg/m3；ρC為氣核密度，kg/m3；g為重力，m/s2；θ為管道傾角,(°)。結合式(1)、(2)，可求解式(3)，輔助方程見參考文獻[9]。

2 液膜速度計算及井底積液預測方法

利用計算機Visual Basic.Net語言編寫程序進行求解時需要提供的基本參數主要有：氣相流速、液相流速、井筒直徑、井壁粗糙度、井筒井斜角、液體密度、氣體密度、液體黏度、氣體黏度、氣液界面張力、管內壓力、管內溫度等。Unified多相流動方程的求解過程是：首先根據井筒內氣液兩相流動狀態，判定流型，根據不同流型采用不同算法。針對環霧流，首先假設一個相應的液膜速度，在依據連續性方程和動量方程計算相應參數，利用牛頓迭代法計算液膜速度，當假設和計算的液膜速度的誤差小于千分之一時，該液膜速度即為相應流動條件下的解。利用計算所得的液膜速度與臨界液膜速度進行對比，如果液膜速度小于臨界液膜速度，井底存在積液；當液膜速度大于臨界液膜速度時，井底則不存在積液。Unified模型預測井底積液問題計算流程如圖2所示。

圖2 Unified模型預測井底積液計算流程

同時，Yuan[11]開展環霧流實驗，實驗觀察了在垂直井筒、傾角75°和60°井筒中井底積液情況。本文利用Yuan[11]的實驗數據與多相流動模型的計算結果相對比，分析液膜速度和井底積液之間的關系。

2.1 數值模擬

如圖3所示，流體在垂直井筒、傾角75°和60°井筒內的流動壓降存在最小值，隨著氣相流速的增加，井筒內壓降先降低，達到最低壓降值后呈緩慢增長。相同氣相流速條件下，液相流速越大，井筒內壓井越大，而在液相流速和氣相流速都相同時，井筒的傾角越接近垂直角，井筒內的壓降越大。Unified模型預測井筒內的流型為間歇流。Veeken[8]開展環霧流中井底積液實驗，指出當井筒內流動壓降達到最低值時，井底開始產生積液，該預測方法可稱為最低壓降預測法。同時，本文中當模型計算壓降接近最低值時，對比Yuan的實驗數據發現，在間歇流中井筒內也出現了積液現象。

圖3 不同傾斜角、表觀液速條件下壓降和表觀氣速的關系

不同氣液流速情況下垂直井筒、傾斜75°井筒和60°井筒內的積液情況如圖4所示，可以看出，表觀液速越大，井底積液發生的可能性越大。在垂直管中，當表觀液速等于0.1 m/s時，盡管表觀氣速達到27.5 m/s，管內仍然出現積液情況。同時，對比不同傾斜角對積液的影響，井筒傾斜角越大，管內越容易出現積液情況，如在相同的表觀液速(0.1 m/s)情況下，表觀氣速從12.5 m/s增加至30 m/s，60°井筒內的無積液情況要明顯多于垂直井筒，表觀氣速越大，越不容易出現井底積液。

圖4 不同氣液流速和井筒傾斜角條件下的井底積液情況預測(Yuan[11])

不同氣液流速情況下垂直井筒、傾斜75°井筒和60°井筒內的液膜速度分布規律如圖5。圖5根據有積液和無積液兩種情況對液膜速度的分布進行的梳理，可以發現，在有積液情況下，垂直井筒、傾角75°井筒和60°井筒內大部分的液膜流速分別低于-0.387 m/s、-0.297 m/s和-0.178 m/s。據此可以得出，當井底開始出現積液液膜開始逆流，同時井筒的傾斜角越小，相同氣液流速情況下，液膜逆流速度越小，井底出現積液的概率也越小。

圖5 不同氣液流速和積液情況下的液膜速度分布規律

2.2 模型對比及驗證

根據上文所述，井底積液的預測可以分別通過預測井筒最小壓降和液膜速度的計算所達到。基于Yuan實驗數據[11]，對比了在不同井筒傾斜角、不同表觀液速和不同表觀氣速條件下，最低壓降法和液膜逆流法的預測結果和實驗誤差，結果見表1。由表1數據對比發現，液膜逆流法的預測誤差要普遍低于最小壓降法。最小壓降法的平均誤差為38.22%，而液膜逆流法的平均誤差為11.56%。因此，液膜逆流法的預測結果更為準確。文中液膜速度為表觀流速。

表1 最小壓降法和液膜逆流法預測井底積液誤差

2.3 影響因素分析

當液膜速度為0時，井筒底部應該出現積液。但是，通過與Yuan實驗數據的對比發現，實際當液膜速度為-0.30 m/s，井底開始出現積液。在Yuan的實驗中，所觀察到的液膜速度為表觀流速，液膜速度是通過觀察液膜中氣泡的運行速度所決定的，而液相和氣泡間的滑脫速度被忽略不計了。但是在實驗中，氣泡和液膜間的滑脫速度實際存在，氣泡速度應等于液膜速度加滑脫速度，即模型所預測的液膜速度-0.30 m/s包括了氣泡的滑脫速度。Harmathy[12]建立了氣泡在液體上升過程中的滑脫速度計算模型，如式(4)，計算參數如表2所示。通過計算發現，在實驗條件下，氣泡與液膜的滑脫速度為0.25 m/s。由此得出當實驗中所觀測的氣泡速度為靜止時，實際液膜速度應為-0.25 m/s，而模型所預測的液膜速度-0.30 m/s，所對應實驗中的實際液膜速度為-0.05 m/s。在前人所做的多相流動實驗中，通常液膜首先發生振蕩現象，然后再開始逆流，液膜在管壁上的振蕩可以避免井底積液的發生，因此，當液膜完全逆流時，即當實際液膜速度為-0.05 m/s時，井底積液才會產生。

表2 滑脫速度計算模型基本參數及氣體滑脫速度

(4)

式(4)中：σL為液體界面張力，N/m。

圖6展示了不同井筒角度和不同液相流速下的液膜速度與氣相速度之間的關系。由圖6可知，在所有條件下隨著氣相速度的降低，液膜速度也隨之降低，該現象表明當氣相速度過低時井底會出現積液。同時，在井筒傾斜角為60°和液相速度為0.100 m/s時，相同的氣相速度條件下，液相速度越大，液膜的速度越小。因此，井筒的井斜角越大，液相速度越大，井筒內發生積液的可能性越大。

圖6 不同井筒傾斜角和表觀液速下的液膜速度和表觀氣速之間的關系曲線

3 結論

本文借助Visual Basic語言編制計算機程序，運用Unified多相流動模型針對傾角分別為垂直井筒、傾角75°和60°井筒開展了井底積液數值模擬研究，分別通過計算井筒內最低壓降和液膜速度來判定井底積液的產生。對比實驗數據發現，最小流動壓降法預測結果的平均誤差為38.22%，而液膜逆流法預測結果的平均誤差為11.56%，最大誤差僅為19%。因此，液膜逆流法的預測誤差更低，當液膜速度為-0.3 m/s時井底發生積液。同時，數值模擬結果表明，當井筒傾角越大和液體表觀速度越大時井筒內產生積液的可能性越大。本文研究結果為液膜流速預測深水氣井井筒內的積液情況提供了新思路。