薛軍 張茂林 宋惠馨 王賀 張藝鐘 彭軍 唐雯

摘 ?????要：氣井最小臨界攜液流量的準(zhǔn)確計(jì)算對(duì)于確定氣井合理配產(chǎn)、優(yōu)化氣田開發(fā)方案具有非常重要的意義。國內(nèi)外學(xué)者在Turner模型的基礎(chǔ)上，對(duì)臨界攜液系數(shù)進(jìn)行修正，從而推導(dǎo)出了不同的臨界攜液流量計(jì)算模型，但這些模型均將界面張力和天然氣偏差因子取為常數(shù)，忽略了溫度和壓力對(duì)它們的影響。因此，對(duì)液滴模型進(jìn)行了修正，通過經(jīng)驗(yàn)公式計(jì)算界面張力，運(yùn)用DAK方法計(jì)算天然氣偏差因子，提出了考慮實(shí)際界面張力和天然氣偏差因子的氣井臨界攜液流量模型。應(yīng)用修正前后的3種常規(guī)模型分別對(duì)某氣田的臨界攜液流量進(jìn)行計(jì)算對(duì)比，結(jié)果表明修正的臨界攜液模型能夠更加準(zhǔn)確預(yù)測(cè)氣井狀態(tài)，判斷井筒是否積液，較常規(guī)模型具有更高的準(zhǔn)確性和科學(xué)性，適用于氣田實(shí)際生產(chǎn)開發(fā)。

關(guān) ?鍵 ?詞：臨界攜液流量;界面張力;天然氣偏差因子

中圖分類號(hào)：TE122????????文獻(xiàn)標(biāo)識(shí)碼：?A ??????文章編號(hào)： 1671-0460（2019）11-2614-04

Correction of Calculation Method of Critical Fluid Carrying Flow in Gas Well

XUE Jun¹， ZHANG Mao-lin²， SONG Hui-xin¹， WANG He¹，

ZHANG Yi-zhong²， PENG Jun³， TANG Wen⁴

（1. School of Petroleum Engineering， Yangtze university， Hubei Wuhan 430100， China;

2. Collaborative Innovation Center for Unconventional Oil and Gas， Yangtze University， Hubei Wuhan 430100， China;

3. Jianghan Oilfield Branch Exploration and Development Research Institute， Hubei Wuhan 430073， China;

4. Southwest Oil and Gas Field Chongqing Gas Field， Chongqing 400021， China）

Abstract： ?The accurate calculation of the minimum critical fluid carrying flow of gas wells is of great significance for the determination of reasonable gas well production allocation and optimization of gas field development program. On the basis of Tuner model， scholars at home and abroad modified the critical fluid carrying coefficient， thus different critical fluid carrying flow calculation models were derived. However， these models take the interfacial tension and natural gas deviation factors as constants and ignore the influence of temperature and pressure on them. In this paper， the droplet model was modified， the interfacial tension?was?calculated by empirical formula， the deviation factor of natural gas was also calculated by DAK method，and the gas well critical fluid carrying flow model considering the actual interfacial tension and the deviation factor of natural gas was proposed. The critical flow rate of a gas field was?calculated by three general models before and after?modification. The results showed?that the modified critical fluid transport model was?more accurate than the conventional model in predicting gas well status and determining whether there was?fluid accumulation. So it is applicable to the actual production and development of gas field.

Key words： critical fluid carrying rate; interfacial tension; deviation factor of natural gas

氣井最小攜液臨界流量的精確計(jì)算對(duì)于氣田開發(fā)方案編制具有十分重要的意義^[1]。在有水氣藏開發(fā)中后期，氣藏壓力逐步降低、氣井產(chǎn)水量逐漸增大，出水導(dǎo)致井筒內(nèi)壓力損耗增加，達(dá)到同樣的產(chǎn)量需要更高的攜液舉升井底壓力，或同樣的井底壓力只能獲得較低的井口產(chǎn)量^[2]，使得氣流難以攜帶井底產(chǎn)出水到達(dá)地面，從而在井底形成積液^[3]。氣井積液會(huì)增大井底回壓，導(dǎo)致氣井產(chǎn)量下降，嚴(yán)重時(shí)甚至?xí)饩畨核?，從而造成停產(chǎn)^[4-7]。為了避免產(chǎn)生積液，氣井產(chǎn)量必須大于臨界攜液產(chǎn)量，通過運(yùn)用臨界攜液理論計(jì)算臨界攜液產(chǎn)量（最低產(chǎn)量），以此評(píng)價(jià)出水氣井的合理產(chǎn)量范圍，確定氣井的合理配產(chǎn)^[8，9]。

1996年Turner^[10]通過對(duì)垂直管流中液相的流動(dòng)方式的研究，提出了球形液滴模型來預(yù)測(cè)積液的形成，并推導(dǎo)出了臨界攜液流量計(jì)算公式。之后許多學(xué)者^[11-14]在液滴模型基礎(chǔ)上分別在模型系數(shù)、液相流動(dòng)方式、液滴形狀等方面作了大量工作。2001年李閔^[15]認(rèn)為液滴在高速氣流中運(yùn)動(dòng)時(shí)由于壓差作用會(huì)變形成橢球體，由此推導(dǎo)出氣井連續(xù)排液臨界攜液模型。2007年王毅忠^[16]根據(jù)高速氣流中運(yùn)動(dòng)的球帽形液滴，認(rèn)為模型中液滴應(yīng)為球帽狀，由此推導(dǎo)出氣井最小攜液臨界流量計(jì)算公式。但是液滴模型中氣水界面張力通常被認(rèn)為是常數(shù)60?mN/m，而實(shí)驗(yàn)表明其數(shù)值隨壓力與溫度的變化而變化^{[17，18]};天然氣偏差因子Z通常取常數(shù)0.88，但其數(shù)值明顯受到溫度和壓力的影響^[19]。因此，在求解氣井臨界攜液流量時(shí)，有必要考慮實(shí)際氣水界面張力和天然氣偏差因子的影響。

1 ?氣井臨界攜液流量計(jì)算公式的修正

液滴模型主要有Turner 模型^[10]、李閔模型^[15]、和王毅忠模型^[16]，各模型的攜液臨界流速不同，見表1。

3種液滴模型中，假設(shè)條件分別為球形、橢球形和球帽形3種狀態(tài)的液滴，曳力系數(shù)C_d是流體作用于顆粒上的曳力對(duì)顆粒在其運(yùn)動(dòng)方向上的投影面積與流體動(dòng)壓力乘積的比值^[20]，與液滴形態(tài)有關(guān)，故不同狀態(tài)的液滴相對(duì)應(yīng)取不同的曳力系數(shù)C_d^[16]。Turner模型中球形液滴曳力系數(shù)取為0.44，李閔模型中橢球形取為1.0，王毅忠模型中球帽形液滴根據(jù)實(shí)驗(yàn)數(shù)據(jù)將曳力系數(shù)取為1.17，將對(duì)應(yīng)的C_d代入各模型臨界流速公式，得到各模型的模型系數(shù)a，見表2。

各模型的攜液臨界流速均可用下公式^[21]表示：

?????（1）

（2）

（3）

式中：ν_g-攜液臨界流速，m/s;

q_c-攜液臨界流量，10⁴m³/d;

ρ_g-天然氣密度，kg/m³;

ρ_l-地層水密度，kg/m³;

σ-氣水界面張力，N/m;

a-模型方程系數(shù);

γ_g-天然氣相對(duì)密度;

p_wf-井底流壓，MPa;

T-井底溫度，K;

Z-天然氣偏差因子;

A-油管截面積，m²;

D-油管直徑，m（注：也可用井口壓力、井口溫度、及井

口壓力溫度下的Z）。

氣水表面張力可由以下經(jīng)驗(yàn)公式^[22-24]進(jìn)行計(jì)算：

?????（4）

s₁ =76exp（-0.036 257 5p）;s₂ =52.5-0.870 18p） ?（5）

式中：σ-溫度為t℃時(shí)水的表面張力，mN/m;

σ₁-溫度為23.33?℃時(shí)水的表面張力，mN/m;

σ₂-溫度為137.78?℃時(shí)水的表面張力，mN/m;

t-某一時(shí)刻水的溫度，℃;

p-某一時(shí)刻壓力值，MPa。

1975年Dranchuk和Abou-Kassem在Starling和Carnahan修正的BWR狀態(tài)方程基礎(chǔ)上，對(duì)1500個(gè)原始Standing-Katz偏差因子數(shù)據(jù)進(jìn)行非線性回歸分析，提出了11個(gè)常系數(shù)的SC-BWR狀態(tài)方程計(jì)算天然氣偏差因子的DAK方法^{[25，26]}。其表達(dá)式為：

（6）

?????（7）

（8）

（9）

（10）

式中：ρ_r—天然氣擬對(duì)比密度;

T_r—天然氣擬對(duì)比溫度;

p_r—擬對(duì)比壓力;

Z—天然氣偏差因子;

a₁-a₁₁—系數(shù)（在DAK模型中各值分別為：a₁=0.326?5，

a₂=-1.070?0，a₃=-0.533?9，a₄=0.015?69，a₅=-0.051?65，

a₆=0.547?5，a₇=-0.736?1，a₈=0.184?4，a₉=0.105?6，

a₁₀=0.613?4，a₁₁=0.721?0）。

該方法為非線性隱式方程，使用牛頓-拉弗森迭代法求解。

將代入式（1），將代入式（2）、式（3），則得到考慮實(shí)際氣水界面張力和天然氣偏差因子的氣井臨界攜液流量公式。

???????（11）

????????????????（12）

（13）

修正后的公式將、考慮在內(nèi)，可以在一定程度上提高氣井臨界攜液流量計(jì)算結(jié)果的準(zhǔn)確性。

2 ?實(shí)例計(jì)算分析

表3基礎(chǔ)數(shù)據(jù)來自某氣田，天然氣相對(duì)密度0.584 5，臨界壓力4.68 MPa，臨界溫度193.6 K;氣水界面張力由實(shí)測(cè)溫度壓力數(shù)據(jù)計(jì)算所得，地層水密度為1 000 kg/m³。

由表3可以看出，對(duì)于氣井，是否考慮實(shí)際界面張力和天然氣偏差因子對(duì)于臨界攜液流量計(jì)算結(jié)果影響較大。利用Turner模型計(jì)算該氣田的臨界攜液流量存在的誤差較大，在判斷氣井是否積液時(shí)與實(shí)際結(jié)果偏差較大，修正后的Turner模型在一定程度上減小了誤差;利用李閔模型計(jì)算該氣田的臨界攜液流量，存在氣井實(shí)際產(chǎn)量在低于李閔模型計(jì)算的臨界攜液流量時(shí)，氣井仍然可以正常生產(chǎn)，未出現(xiàn)積液的現(xiàn)象，修正后的李閔模型計(jì)算的臨界攜液流量與實(shí)際生產(chǎn)資料更加符合;利用王毅忠模型計(jì)算該氣田的臨界攜液流量結(jié)果誤差較小，均可以較準(zhǔn)確地預(yù)測(cè)氣井狀態(tài)，但是修正后的王毅忠模型進(jìn)一步提高了計(jì)算結(jié)果的準(zhǔn)確性和可靠性，與實(shí)際情況更加吻合。

將實(shí)際氣井產(chǎn)量與各公式計(jì)算的最小攜液產(chǎn)量進(jìn)行比較。若實(shí)際產(chǎn)氣量等于最小攜液產(chǎn)量，數(shù)據(jù)點(diǎn)就會(huì)落在圖中的對(duì)角線上。如果公式適用，那么，接近積液氣井的數(shù)據(jù)應(yīng)該接近這條對(duì)角線，未積液氣井的數(shù)據(jù)應(yīng)在這條對(duì)角線的上面，而積液氣井的數(shù)據(jù)落在對(duì)角線的下面。由圖1可以看出，修正后的Turner模型計(jì)算結(jié)果偏大，不能準(zhǔn)確預(yù)測(cè)氣井積液狀態(tài);修正后的李閔模型和王毅忠模型計(jì)算結(jié)果跟生產(chǎn)實(shí)際值較為吻合，可以較準(zhǔn)確預(yù)測(cè)氣井積液狀態(tài)^[27，28]。

由此可見，根據(jù)實(shí)際生產(chǎn)資料，修正前后的Turner模型均不適用于該氣田;修正前后的李閔模型雖然在判斷氣井是否積液時(shí)結(jié)論比較一致，但由計(jì)算結(jié)果可以看出，修正后的李閔模型更加符合實(shí)際;修正前后的王毅忠模型計(jì)算結(jié)果均適用于該氣田，但修正后的王毅忠模型準(zhǔn)確性和可靠性更高。

對(duì)于氣田，單井合理配產(chǎn)也是一個(gè)重要的開發(fā)指標(biāo)，尤其是對(duì)于出水氣井，在實(shí)際生產(chǎn)中合理配產(chǎn)的確定需充分考慮氣井最小臨界攜液流量，避免由于配產(chǎn)過低造成井筒積液，對(duì)氣層造成損害，產(chǎn)生不可估量的經(jīng)濟(jì)損失。因此，在確定合理配產(chǎn)時(shí)，應(yīng)該綜合分析修正前后模型的臨界攜液流量結(jié)果，避免氣井水淹，以保證氣田正常、有序開發(fā)。

3 ?結(jié)論

（1）在結(jié)合生產(chǎn)實(shí)踐的基礎(chǔ)上，對(duì)計(jì)算臨界攜液流量的常用模型：Turner模型、李閔模型和王毅忠模型在理論上進(jìn)行了修正，提出了考慮實(shí)際界面張力和天然氣偏差因子的臨界攜液流量計(jì)算公式。

（2）根據(jù)某氣田實(shí)際生產(chǎn)資料分析，是否考慮實(shí)際界面張力和天然氣偏差因子對(duì)臨界攜液流量的計(jì)算結(jié)果有較大影響?？紤]實(shí)際界面張力和天然氣偏差因子的李閔模型和王毅忠模型計(jì)算結(jié)果與氣井實(shí)際生產(chǎn)資料更加吻合，具有更高的準(zhǔn)確性、可靠性和科學(xué)性，可作為氣井積液的判斷標(biāo)準(zhǔn)。

（3）對(duì)于氣井，尤其是出水氣井，應(yīng)當(dāng)采用多種模型相結(jié)合以及修正前后模型相結(jié)合的方法，以此分析氣井?dāng)y液能力，進(jìn)行合理配產(chǎn)。

參考文獻(xiàn)：

[1]李元生， 滕賽男，楊志興，等. 考慮界面張力和液滴變形影響的攜液臨界流量模型 [J] . 石油鉆采工藝，2017， 39 （2）：218-223.

[2]秦濤， 柴小穎， 毛鳳華， 謝梅， 李興彥. 澀北一號(hào)氣田井筒積液診斷模式及判斷[J]. 青海石油， 2012， 30 （02）： 51-55.

[3]魏納，李穎川，李悅欽，等. 氣井積液可視化實(shí)驗(yàn)[J]. 鉆采工藝，2007，30 （3）：43-45.

[4]劉雙全，吳曉東，吳革生，等. 氣井井筒攜液臨界流速和流量的動(dòng)態(tài)分布研究[J]. 天然氣工業(yè)，2007， 27 （2）：104-106.

[5]譚曉華，李曉平. 考慮氣體連續(xù)攜液及液滴直徑影響的氣井新模型[J] . 水動(dòng)力學(xué)研究與進(jìn)展，2013， 28 （1）：41-47.

[6]楊蕾， 劉寧， 劉輝， 王華軍. 氣井井筒攜液模型研究與應(yīng)用[J]. 天然氣與石油， 2011， 29 （02）： 48-49+61+2.

[7]李品， 江厚順， 王萌， 隋海慶， 吳岸. 氣井井筒積液預(yù)測(cè)模型優(yōu)選[J]. 當(dāng)代化工， 2016， 45 （08）： 1820-1823.

[8]羅睿. ?JHS低滲透氣藏氣藏工程研究[D]. 西南石油大學(xué)， 2014.

[9]潘杰， 王武杰， 魏耀奇， 陳軍斌， 王亮亮. 考慮液滴形狀影響的氣井臨界攜液流速計(jì)算模型[J]. 天然氣工業(yè)， 2018， 38 （01）： 67-73.

[10]TURNER R G，HUBBARD M G，DUKLER A E.Analysis and prediction of minimum flow rate for the continuous removal of liquids from gas wells[J]. Journal of Petroleum Technology， 1969， 21（11）： 1475-1482.

[11]COLEMAN S B， CLAY H B， MCCURDY D G， LILI L N.A new look at predicting gas-well load-up[J]. Journalof Petroleum Technology， 1991， 43（3）： 329-333.

[12]Nosseir，M.A. A New Approach for AccuratePrediction of Loading in Gas Well Under Different FlowingConditions[C].SPE 37408，1997.

[13]彭朝陽. 氣井?dāng)y液臨界流量研究[J]. 新疆石油地質(zhì)，2010，31（1）：72-74.

[14]熊健，李凌峰，張濤. 氣井?dāng)y液臨界流量的計(jì)算方法[J]. 天然氣與石油，2011，29（4）：54-56.

[15]李閩，郭平，譚天光. 氣井?dāng)y液新觀點(diǎn)[J]. 石油勘探與開發(fā)，2001， 28（5）：105-106.

[16]王毅忠，劉慶文. 計(jì)算氣井最小攜液臨界流量的新方法[J]. 大慶石油地質(zhì)與開發(fā)，2007， 26（6）：82-85.

[17]SUTTON R P.?An improved model for water-hydrocarbon surface tension at reservoir conditions[R].?SPE 124968，?2009.

[18]劉剛. 氣井?dāng)y液臨界流量計(jì)算新方法[J]. 斷塊油氣田，2014， 21（3）：339-343.

[19]郭布民，敬季昀，王杏尊. 氣井臨界攜液流量計(jì)算方法的改進(jìn)[J]. 斷塊油氣田，2018，25（4）：484-487.

[20]鄭思佳. 氣固兩相流中采氣管線關(guān)鍵易損件沖蝕特性研究[D].西南石油大學(xué)， 2016.

[21]紀(jì)成. N氣藏動(dòng)態(tài)分析及產(chǎn)能研究[D]. 西南石油大學(xué)， 2015.

[22]楊繼盛.采氣工藝基礎(chǔ)[M].北京：石油工業(yè)出版社，1994：330-334.

[23]陳元千. 油氣藏工程實(shí)用方法[M]. 北京： 石油工業(yè)出版社，1990：442-452.

[24]楊繼盛，劉建儀.采氣實(shí)用計(jì)算[M].北京：石油工業(yè)出版社，2004：42-46.

[25]Dranchuk P M，Abou-Kassem J H.Calculations of z factors for natural gases using equations of state[J].Journal Canadian Petroleum，?1975，?14?（3）：?34-36.

[26]胡建國，郭分喬，許進(jìn)進(jìn). 計(jì)算天然氣偏差因子的DAK方法的修正[J]. 石油與天然氣地質(zhì)，2013， 34（1）：120-123.

[27]鄭軍， 閆長(zhǎng)輝， 張文洪， 鄭峰， 徐衛(wèi)峰. 大牛地氣田氣井最小攜液產(chǎn)量研究[J]. 油氣地質(zhì)與采收率， 2011， 18 （01）： 70-73+116.

[28]李閩， 郭平， 張茂林， 梅海燕， 劉武， 李士倫. 氣井連續(xù)攜液模型比較研究[J]. 西南石油學(xué)院學(xué)報(bào)， 2002 （04）： 30-32+3.