高壓含水氣井兩相流態(tài)節(jié)流判斷與控制方法

劉新福，劉春花，李清平，姚海元，郝忠獻，劉永飛

(1.青島理工大學 機械與汽車工程學院，山東 青島 266520；2.中國石油大學(華東) 機電工程學院，山東青島 266580；3.中海油研究總院有限責任公司，北京 100028；4.中國石油勘探開發(fā)研究院，北京 100083)

井筒氣水兩相管流節(jié)流前后動態(tài)特性分析與控制，對保障致密氣等含水氣井高效穩(wěn)定安全開采和解決因水合物淤積造成管阻問題具有重要意義[1-3]。由于氣井井下節(jié)流噴嘴的降壓作用，井口壓力得到有效調(diào)節(jié)，井口設備與管道承壓下降[4-6]；此外，節(jié)流后的含水氣流一般處于臨界流動狀態(tài)，井口設備的壓力波動不會對井底壓力產(chǎn)生影響，減少了壓力波動和井底出砂。

對于氣液兩相流態(tài)節(jié)流過程，目前主要通過研究節(jié)流中流體質(zhì)量、動量、能量與相態(tài)的變化，提出了經(jīng)驗和機理相結合的算法。Shao Huaishuang 等[7]在Perkins 和Al-safran 模型基礎上，提出一種節(jié)流閥多相流態(tài)和質(zhì)量流量預測方法，該方法考慮了流體經(jīng)噴嘴節(jié)流時的交替氣液絕熱條件和氣液兩相間的滑移現(xiàn)象，但模型中兩相流為均相流，且僅適用于臨界流與亞臨界流。C.M.Carstensen 等[8]在Gilbert 和Ros 模型基礎上，修正非平衡系數(shù)并解釋節(jié)流閥快速降壓的中間狀態(tài)，提出運用非平衡流動的概念預測凍結流動和絕熱膨脹的節(jié)流過程，但卻忽略了組分、溫壓、體積與流體溫度等因素的影響。王振嘉等[9]綜合產(chǎn)量對比分析法和油套環(huán)空壓差預測法，并結合井下節(jié)流噴嘴尺寸、下井深度等參數(shù)，提出一種適用于“低滲透率、低孔隙率、低豐度、強非均質(zhì)性”特征儲層井下節(jié)流氣井的積液量預測方法，其適用范圍受限。安永生等[10]建立了井下節(jié)流氣井滲流、噴嘴氣流數(shù)學模型，并依此推導數(shù)值求解源匯項方程，完成氣井滲流和噴嘴氣流的耦合處理，僅適用于遵守達西定律的單相氣體流動，且假設節(jié)流前后為等溫過程。曾焱等[11]建立了節(jié)流氣井儲層和井筒耦合生產(chǎn)動態(tài)模型，擬合試井解釋或生產(chǎn)動態(tài)數(shù)據(jù)，并依此預測節(jié)流氣井的流量和壓力變化情況，但僅適用于單相氣體流動工況。

S.Naseri 等[12]在GA-RBF 模型基礎上，提出一種亞音速氣流經(jīng)噴嘴和孔板式節(jié)流閥時人工智能預測流量系數(shù)方法，該方法采用徑向基函數(shù)神經(jīng)網(wǎng)絡結合遺傳算法，并考慮了雷諾數(shù)以及節(jié)流噴嘴直徑與管徑比值等因素。王熒光等[13]模擬井下節(jié)流過程壓力和溫度變化情況，依據(jù)NETOPT 網(wǎng)絡模型優(yōu)化節(jié)流噴嘴內(nèi)徑等結構參數(shù)。這些研究可歸納為網(wǎng)絡優(yōu)化模型，但由于算法的實際應用效果欠佳，目前工程上還沒有得到普遍應用。

目前，大部分研究成果不適用于高壓含水氣井或存在工程上應用受限等問題，筆者綜合等熵絕熱膨脹、焓與動能非平衡轉化等過程，提出一種適用于致密氣、濕氣井等高壓含水氣井開采工況下的井筒氣水兩相管流節(jié)流前后動態(tài)特性分析與控制方法，并通過數(shù)值模擬和在鄂爾多斯盆地東緣大寧-吉縣區(qū)塊井場驗證，分析高壓含水氣井節(jié)流過程隨噴嘴尺寸、含水率以及節(jié)流前壓力和溫度動態(tài)變化狀況等，以期達到保證高壓含水氣井高效穩(wěn)定開采和防治水合物的目的。

1 氣水兩相井下節(jié)流熱力學模型

氣水兩相流經(jīng)節(jié)流噴嘴時，井下節(jié)流過程的能量變化基本方程[14-15]為：

氣水兩相流體流經(jīng)節(jié)流噴嘴的時間較短，節(jié)流前后氣水兩相流體幾乎不與井筒及其近井圍巖發(fā)生熱功和機械功的交換，則節(jié)流過程氣水兩相流體的能量組成為：

由此推導出噴嘴出口處兩相流的速度為：

井下節(jié)流前后，流經(jīng)噴嘴的氣水兩相流體的焓轉化為動能，節(jié)流后兩相流速得以提升，等熵絕熱過程中，節(jié)流后的兩相流體壓力下降，而比熱容增大；與此同時，節(jié)流后兩相流體的熵不變，但兩相流體的溫度下降。

2 氣水兩相井下節(jié)流場數(shù)學模型

不考慮井下節(jié)流過程的位能差、系統(tǒng)與外界的熱交換，且節(jié)流噴嘴內(nèi)部流體處于等熵絕熱狀態(tài)，依據(jù)節(jié)流過程能量守恒方程[16-18]，氣水兩相管流井下節(jié)流過程的能量變化(圖1)為：

圖1 高壓含水氣井筒流體能量隨井深變化情況Fig.1 Depth-varying fluid energy in wellbores of high-pressure gas wells producing water

單位質(zhì)量氣水兩相流體的能量組成為：

節(jié)流過程中系統(tǒng)熱量Q保持不變，則有：

基于熵的流體熱力學微分方程為：

由此推導出節(jié)流過程系統(tǒng)熱量微分方程為：

井下節(jié)流過程中，氣相等壓熱容(Cp)[19-20]與氣相等容熱容(CVg)之間的關系式為：

考慮等壓熱容、等容熱容等因素后的節(jié)流系統(tǒng)熱量為：

其中，

引入節(jié)流后出口與節(jié)流前入口的壓力比β，減小壓力比β，節(jié)流過程氣水兩相流體質(zhì)量流量有增大的趨勢，且節(jié)流后出口壓力降至一定程度后，流體速度存在音速界限并保持不變，即兩相流體達到音速的臨界流動狀態(tài)[21-24]，臨界流動狀態(tài)單位質(zhì)量兩相流體能量為：

井下節(jié)流過程中物質(zhì)的平衡方程[25-28]為：

井下節(jié)流過程流體的質(zhì)量流量為：

井下節(jié)流入口壓力與出口壓力之間的壓力比低于臨界壓力比[29-31]時，氣水兩相流體處于臨界流動狀態(tài)，質(zhì)量流量達到最大并保持不變，即：

3 高壓含水氣井節(jié)流特性實例分析

3.1 基本參數(shù)

選取鄂爾多斯盆地東緣大寧-吉縣區(qū)塊致密氣高壓含水氣井開展井筒氣水兩相流態(tài)井下節(jié)流前后動態(tài)特性分析。依據(jù)井下技術作業(yè)公司的檢測報告，大寧-吉縣區(qū)塊井場測試的井筒各部分物性參數(shù)與生產(chǎn)參數(shù)為：節(jié)流前壓力14、15、16、17、18 MPa，節(jié)流前溫度313、318、323、328、333 K，氣體密度0.68 kg/m3，含水率0、0.1、0.2、0.3、0.4 mol/mol，噴嘴內(nèi)徑3.0、3.5、4.0、4.5、5.0 mm，油管內(nèi)徑63 mm，油套內(nèi)徑139.7 mm。

3.2 節(jié)流前后兩相流動態(tài)特性

高壓含水氣井節(jié)流前后兩相流的動態(tài)特性會受到噴嘴尺寸等結構參數(shù)的影響。節(jié)流前壓力18 MPa、溫度333 K 和含水率0.1 mol/mol 時，氣水兩相流體節(jié)流過程質(zhì)量流量隨噴嘴內(nèi)徑和壓力比的動態(tài)變化關系如圖2 所示。減小壓力比β，節(jié)流過程兩相流體質(zhì)量流量先是呈指數(shù)級增大而后基本保持恒定，壓力比的閾值為0.55 左右。

圖2 高壓含水氣井節(jié)流過程隨噴嘴尺寸動態(tài)變化狀況Fig.2 Dynamic variations of the throttling process in high-pressure gas wells producing water under different inner diameters and setting depths of chokeling nozzles

噴嘴內(nèi)徑等結構參數(shù)對高壓含水氣井兩相流節(jié)流過程質(zhì)量流量的影響相對較大，且給定壓力比β時，節(jié)流過程質(zhì)量流量隨噴嘴內(nèi)徑的不斷增大而顯著提升，這與安永生等[10]所得出的“井下節(jié)流噴嘴直徑越大，氣井產(chǎn)量越高”的結論一致。圖2a 中，噴嘴內(nèi)徑由3.0 mm擴至5.0 mm 時，氣水兩相流體的臨界質(zhì)量流量由1.45×104kg/d 增至4.05×104kg/d，臨界質(zhì)量流量提升幅度達到179.3%。圖2b 給出不同壓力比β工況節(jié)流后溫度變化情況，當β=1 時，節(jié)流后溫度等于節(jié)流前溫度；而后隨壓力比不斷減小，節(jié)流后溫度呈現(xiàn)類“拋物線”式下降的趨勢。其主要原因在于節(jié)流前入口壓力一定的情況下，節(jié)流后出口壓力的下降會造成壓力比β的減小，節(jié)流過程氣水兩相流體質(zhì)量流量有增大的趨勢，綜合等壓熱容、等容熱容等因素作用，節(jié)流之后出現(xiàn)兩相流體溫度降低的現(xiàn)象。此外，節(jié)流前入口壓力一定時，噴嘴結構參數(shù)設計不合理易使得節(jié)流后出口壓力發(fā)生“低壓”工況(圖2b 中壓力比小于0.1)，由此會造成節(jié)流之后的兩相流體溫度降至零度以下；而當節(jié)流后出口壓力一定時，不斷增加典型“高壓”含水氣井的節(jié)流前入口壓力(14～18 MPa)，也易出現(xiàn)節(jié)流后兩相流體溫度降至零度以下的現(xiàn)象。

圖2c 中，給定壓力比β時，噴嘴下入深度越深，節(jié)流后溫度越高，這與曾焱等[11]所得出的“井下節(jié)流溫度隨噴嘴深度的不斷增大而逐漸提升”的結論一致。此外，圖2c 中氣水兩相流體節(jié)流后溫度隨噴嘴下入深度和壓力比的不斷減小而呈現(xiàn)類“拋物線”式逐漸下降的趨勢；且噴嘴下入深度較淺且壓力較小時，節(jié)流之后的兩相流體溫度可降至-293 K～-333 K 而出現(xiàn)低溫工況，其主要原因在于高壓含水氣井節(jié)流前后兩相流動態(tài)特性分析充分考慮了等熵絕熱、等壓熱容、等容熱容、不同井深節(jié)流能量等因素的綜合作用，而且低溫高壓工況的出現(xiàn)易導致水相態(tài)發(fā)生變化并伴隨水合物生成的情況，從而不利于水合物生成的防治；圖2c 中噴嘴下入深度400 m 且壓力比0.01 的工況下，氣水兩相流體的節(jié)流后溫度最低可降至-325.7 K。為此，需合理設計高壓含水氣井噴嘴的下入深度，避免出現(xiàn)井深超1.6 km 而節(jié)流作業(yè)中噴嘴下入深度不足1.0 km 的情況；并且需合理制定節(jié)流前后的壓力等生產(chǎn)制度，防止發(fā)生因壓力比過小而造成節(jié)流后溫度過低等問題；同時，充分利用地層溫度對兩相流體進行升溫，以保證氣水兩相流體節(jié)流過程中的溫度和節(jié)流后的溫度高于水合物的臨界生成溫度。

高壓含水氣井節(jié)流前后的動態(tài)特性也受到壓力等生產(chǎn)參數(shù)的影響，噴嘴內(nèi)徑5.0 mm、節(jié)流前溫度333 K和含水率0.1 mol/mol 工況下，氣水兩相流體節(jié)流過程質(zhì)量流量以及節(jié)流后溫度隨節(jié)流前壓力和壓力比的動態(tài)變化關系如圖3 所示。

圖3 高壓含水氣井節(jié)流過程隨壓力動態(tài)變化狀況Fig.3 Dynamic variations of the chokeling process in high-pressure gas wells producing water under different pressures

高壓含水氣井氣水兩相流體節(jié)流過程質(zhì)量流量隨節(jié)流前壓力的增大而緩慢提升，且較噴嘴內(nèi)徑的影響較低，圖3a 中，節(jié)流前壓力由14 MPa 增至18 MPa 時，氣水兩相流體的臨界質(zhì)量流量由3.17×104kg/d 升至4.05×104kg/d，臨界質(zhì)量流量提升的幅度僅為27.8%。圖3b中，節(jié)流過程給定壓力比β時，節(jié)流前壓力越高，節(jié)流后溫度越低，且其影響程度隨壓力比下降而顯著增大。

高壓含水氣井節(jié)流前后的動態(tài)特性還受到溫度等生產(chǎn)參數(shù)的影響，噴嘴內(nèi)徑5.0 mm、節(jié)流前壓力18 MPa和含水率0.1 mol/mol 工況下，氣水兩相流體節(jié)流過程質(zhì)量流量以及節(jié)流溫差隨節(jié)流前溫度和壓力比的動態(tài)變化關系如圖4 所示。

圖4 高壓含水氣井節(jié)流過程隨溫度動態(tài)變化狀況Fig.4 Dynamic variations of the chokeling process in highpressure gas wells producing water under different temperatures

節(jié)流前溫度等生產(chǎn)參數(shù)對高壓含水氣井氣水兩相流體節(jié)流過程質(zhì)量流量以及節(jié)流溫差的影響均較小，節(jié)流過程質(zhì)量流量隨節(jié)流前溫度的升高而略有降低，圖4a中，節(jié)流前溫度由313 K 升至333 K 時，氣水兩相流體的臨界質(zhì)量流量由4.27×104kg/d 降至4.05×104kg/d，下降幅度僅為5.15%。圖4b 中，節(jié)流過程給定壓力比β時，節(jié)流前溫度越高，節(jié)流溫差ΔT=Tout-Tin越低，且其影響程度隨壓力比下降而顯著增大。

高壓含水氣井節(jié)流前后的動態(tài)特性同時會受到含水率等物性參數(shù)的影響，噴嘴內(nèi)徑5.0 mm、節(jié)流前壓力18 MPa 和溫度333 K 工況下，氣水兩相流體節(jié)流過程質(zhì)量流量以及節(jié)流后溫度隨含水率和壓力比的動態(tài)變化關系如圖5 所示。

圖5 高壓含水氣井節(jié)流過程隨含水率動態(tài)變化狀況Fig.5 Dynamic variations of the chokeling process in highpressure gas wells producing water under different water contents

含水率等物性參數(shù)對高壓含水氣井氣水兩相流體節(jié)流過程質(zhì)量流量的影響相對較小，與節(jié)流前壓力的影響接近，圖5a 中，含水率由0 增至0.4 mol/mol 時，氣水兩相流體的臨界質(zhì)量流量由4.04×104kg/d 升至5.01×104kg/d，增幅為23.76%。圖5b 中，節(jié)流過程給定壓力比時，含水率越高，節(jié)流后溫度也越高。

3.3 節(jié)流判斷結果對比

氣水兩相流體井筒運移節(jié)流過程臨界質(zhì)量流量受噴嘴內(nèi)徑的影響最大，這與安永生等[10]所提出的“定氣嘴尺寸生產(chǎn)”的概念一致，含水率和節(jié)流前壓力的影響次之，而節(jié)流前溫度和噴嘴下入深度的影響最小。

相較于安永生、曾焱等前人研究[9-11]，一方面筆者團隊綜合等熵絕熱膨脹、焓與動能非平衡轉化等過程，提出了一種適用于致密氣等高壓含水氣井開采工況和氣水兩相流體節(jié)流前后動態(tài)特性分析方法，而不是僅適用于單相氣體流動工況、特定特征儲層節(jié)流積液量預測或假設節(jié)流前后為等溫過程；另一方面，本文不僅量化分析了噴嘴內(nèi)徑等結構參數(shù)對節(jié)流過程動態(tài)特性的影響，而且揭示出含水率和噴嘴深度等物性參數(shù)以及壓力和溫度等生產(chǎn)參數(shù)對高壓含水氣井節(jié)流特性的具體影響程度，初步解決了現(xiàn)有算法實際應用效果欠佳和工程上應用受限的問題，并為優(yōu)化井下節(jié)流器及其節(jié)流噴嘴的結構參數(shù)和保障氣井流動安全提供了理論依據(jù)。

4 結論

a.綜合等熵絕熱、等壓熱容、等容熱容、不同井深節(jié)流能量等因素，推導出氣水兩相流體經(jīng)噴嘴的能量、動能和溫度動態(tài)變化等熱力學微分方程組，提出了一種高壓氣水兩相節(jié)流動態(tài)特性分析與節(jié)流控制方法，為優(yōu)化井下節(jié)流器及其節(jié)流噴嘴的結構參數(shù)和保障氣井流動安全提供了理論依據(jù)，并可避免致密氣等高壓含水氣井因水合物淤積而引起的管阻問題。

b.增大噴嘴內(nèi)徑及其下入深度和節(jié)流前壓力同時降低節(jié)流前溫度，有利于提高兩相流體節(jié)流過程中臨界質(zhì)量流量，并提升高壓含水氣井的產(chǎn)氣量；而增加含水率雖然會小幅增加臨界質(zhì)量流量，但產(chǎn)氣量也隨之下降；且節(jié)流臨界質(zhì)量流量受噴嘴內(nèi)徑的影響最大，含水率和節(jié)流前壓力次之，而節(jié)流前溫度和噴嘴深度的影響最小。

c.在高壓含水氣井的實際生產(chǎn)中，應合理設計噴嘴下入深度和節(jié)流前后壓力，并控制含水率和節(jié)流溫度，同時充分利用地層溫度對兩相流體進行升溫，達到保障高壓含水氣井高效穩(wěn)定開采和防治水合物的目的。

d.致密氣、濕氣井等高壓含水氣井開采工況受多因素影響，且高壓氣水兩相流體節(jié)流前后流場及運移規(guī)律復雜，本次研究成果僅用于高壓含水氣井節(jié)流前后兩相流動態(tài)特性分析，后續(xù)可繼續(xù)開展高壓含水氣井不同開采工況噴嘴動態(tài)節(jié)流場仿真分析和節(jié)流噴嘴優(yōu)化設計工作。

符號注釋：

A為連接節(jié)流器的油管流道橫截面積，m2；Ain、Aout分別為節(jié)流前噴嘴入口處和節(jié)流后噴嘴出口處的節(jié)流器流道橫截面積，m2；CV為流體的等容熱容，J/(mol·K)；Cp、CVg、CVl分別為氣相等壓熱容、等容熱容和水相等容熱容，J/(mol·K)；fg、fl分別為氣相與水相的滯留率；fλ為滯留率；Hin、Hout分別為噴嘴入口和出口處流體焓，J/kg；J為熱功當量，J/kg；Mg為氣相的相對分子質(zhì)量，J/(kg·K)；n為節(jié)流系統(tǒng)熱量系數(shù)；p為流體壓力，MPa；pin、pout分別為節(jié)流前后流體壓力，MPa；Q為總交換熱，J/kg；R為氣體通用常數(shù)，J/(mol·K)；T為流體溫度，K；Tin、Tout分別為節(jié)流前后流體溫度，K；v為節(jié)流過程中流體速度，m/s；vin、vout分別為噴嘴入口和出口處流體的速度，m/s；V為流體的摩爾體積，L/mol；VM為井筒壁面流體的摩爾體積，L/mol；Vin、Vout分別為節(jié)流前后流體的摩爾體積，L/mol；w為節(jié)流過程流體的質(zhì)量流量，kg/d；W為總交換功，J/kg；Z為壓縮因子；ρl為水相密度，kg/m3。