田云，王志彬，李穎川，白慧芳，李克智

（1.西南石油大學(xué)石油與天然氣工程學(xué)院，四川 成都 610500；2.西南石油大學(xué)油氣藏地質(zhì)及開(kāi)發(fā)工程國(guó)家重點(diǎn)實(shí)驗(yàn)室，四川 成都 610500；3.中國(guó)石化華北分公司工程技術(shù)研究院，河南 鄭州 450006）

速度管排水采氣井筒壓降模型的評(píng)價(jià)及優(yōu)選

田云1，王志彬2，李穎川2，白慧芳1，李克智3

（1.西南石油大學(xué)石油與天然氣工程學(xué)院，四川 成都 610500；2.西南石油大學(xué)油氣藏地質(zhì)及開(kāi)發(fā)工程國(guó)家重點(diǎn)實(shí)驗(yàn)室，四川 成都 610500；3.中國(guó)石化華北分公司工程技術(shù)研究院，河南 鄭州 450006）

準(zhǔn)確計(jì)算速度管中氣液兩相流壓力降，是速度管排水采氣工藝優(yōu)化設(shè)計(jì)、生產(chǎn)動(dòng)態(tài)及排液效果分析的基礎(chǔ)。文中以O(shè)rville Gaither 1963年以天然氣/水作為實(shí)驗(yàn)流動(dòng)介質(zhì)，在管徑25.4，31.75 mm速度管中測(cè)試的單相液流、氣液兩相流的實(shí)驗(yàn)數(shù)據(jù)為基礎(chǔ)，首先優(yōu)選了摩阻系數(shù)計(jì)算方法，其次對(duì)速度管中氣液兩相流的壓降模型進(jìn)行評(píng)價(jià)和優(yōu)選，最后利用大牛地3口水平井速度管流壓測(cè)試數(shù)據(jù)進(jìn)行了驗(yàn)證。利用單相液流測(cè)壓數(shù)據(jù)進(jìn)行摩阻系數(shù)優(yōu)選表明，AGA方法計(jì)算結(jié)果與實(shí)際情況最吻合，誤差最?。▔航灯骄^對(duì)誤差為22.37%）；8個(gè)常用氣液兩相管流壓降模型評(píng)價(jià)表明，Gray模型準(zhǔn)確性最好，其次為Ansari模型。這為速度管排水采氣井優(yōu)選出了可靠的兩相流壓降計(jì)算模型，有助于提高工藝設(shè)計(jì)及排液效果診斷水平。

速度管；排水采氣；井筒；壓降；評(píng)價(jià)與優(yōu)選

對(duì)于低壓、低滲、低產(chǎn)及有水氣藏，在生產(chǎn)中后期，部分井甚至在生產(chǎn)初期，氣井產(chǎn)量就開(kāi)始低于攜液臨界氣流量［1］，井筒逐漸積液，嚴(yán)重影響氣井的正常生產(chǎn)，部分氣井甚至?xí)霈F(xiàn)積液停產(chǎn)現(xiàn)象。為此，需通過(guò)降壓帶液、泡排、氣舉、優(yōu)選管柱（速度管）等排水采氣工藝措施，來(lái)維持氣井連續(xù)生產(chǎn)［2-3］。

速度管柱排水采氣工藝技術(shù)能較好地解決積液氣井存在的問(wèn)題，目前，其研究主要是依據(jù)連續(xù)攜液理論進(jìn)行速度管柱尺寸的選擇及工藝效果評(píng)價(jià),而對(duì)速度管內(nèi)的壓降和摩阻分布研究涉及較少。筆者利用Orville Gaither［4］1963年公開(kāi)發(fā)表的實(shí)驗(yàn)數(shù)據(jù)和大牛地3口典型井的測(cè)試數(shù)據(jù)，對(duì)速度管井筒壓降模型進(jìn)行了優(yōu)選和評(píng)價(jià)，這為速度管排水采氣井工藝優(yōu)化設(shè)計(jì)和速度管排液效果分析提供了必要的理論依據(jù)。

1 速度管排水采氣工藝概述

速度管排水采氣工藝，是將優(yōu)選的較小管徑油管（即速度管）作為生產(chǎn)管柱，安裝于井筒，進(jìn)而降低井筒臨界攜液流速，增大井筒中氣體的流速，減少氣液的滑脫損失，提高氣井帶液能力，以達(dá)到排除井底積液、維持氣井平穩(wěn)生產(chǎn)的目的［5］。速度管排水采氣技術(shù)，工藝可靠，可不壓井作業(yè)，減輕了對(duì)地層的傷害；一次性投入低，依靠氣井自身能量排水，節(jié)省了人力、物力。該工藝技術(shù)在四川、大牛地、長(zhǎng)慶等氣田得到了廣泛的試驗(yàn)和應(yīng)用，取得了較好的排水采氣效果，為改善低產(chǎn)低效井的開(kāi)發(fā)效果提供了新的技術(shù)手段。隨著油氣田的不斷開(kāi)發(fā)，低產(chǎn)氣井逐年增多，該技術(shù)應(yīng)用前景廣闊［6-7］。

速度管排水采氣工藝設(shè)計(jì)的理論基礎(chǔ)主要是連續(xù)攜液理論［8］，以及速度管氣液兩相管流壓降計(jì)算模型。對(duì)于產(chǎn)氣量較高、氣液比較高、井筒流型為環(huán)狀流的氣井，可對(duì)比不同管柱尺寸的連續(xù)攜液氣流量與氣井產(chǎn)氣量的關(guān)系，優(yōu)化管柱尺寸及下深；但對(duì)于產(chǎn)氣量較低、氣液比較低、井筒流型為段塞流的氣井，若安裝速度管后仍不能連續(xù)攜液，則可對(duì)比分析井筒的滑脫損失、摩擦壓降及總壓降，以確定合理的速度管尺寸：這就需要正確選取氣液兩相管流壓降模型，以提高速度管設(shè)計(jì)的準(zhǔn)確性和可靠性。同時(shí)，速度管氣液兩相壓降計(jì)算模型，也是速度管排水采氣井動(dòng)態(tài)分析、井筒積液量計(jì)算及排液效果診斷的理論基礎(chǔ)［9-10］。

2 摩阻系數(shù)計(jì)算方法優(yōu)選

速度管內(nèi)氣液混合物流速較高，摩阻壓降較大。準(zhǔn)確計(jì)算摩阻壓降是計(jì)算速度管總壓降的重要方面［11］，而摩阻系數(shù)是計(jì)算摩阻壓降的重要參數(shù)，其計(jì)算的準(zhǔn)確性直接決定了摩阻壓降計(jì)算的準(zhǔn)確性；因此，有必要對(duì)多相管流的摩阻系數(shù)進(jìn)行分析。

研究表明，氣液兩相管流摩阻壓降可采用單相氣流或單相液流的摩阻系數(shù)進(jìn)行近似計(jì)算［12］。Weymouth 1912年從生產(chǎn)實(shí)踐中歸納出計(jì)算井筒摩阻壓降的純經(jīng)驗(yàn)公式。Moody 1944年發(fā)表了各種自然粗糙管道的摩阻系數(shù)圖版，被公認(rèn)為單相流體沿程摩阻系數(shù)的基準(zhǔn)圖版，被廣泛應(yīng)用［13］。Panhandle認(rèn)為，摩阻系數(shù)λ是雷諾數(shù)Re的函數(shù)，理論上只適用于水力光滑區(qū)（這是由于美國(guó)的管子很光滑（取粗糙度k=0.020 mm））。之后，Panhandle提出了可用于雷諾數(shù)較大的阻力平方區(qū)，即 Panhandle B式。美國(guó)燃?xì)鈪f(xié)會(huì)（American Gas Association，簡(jiǎn)稱AGA）認(rèn)為：在紊流流速較低（不完全紊流）時(shí)，摩阻系數(shù)λ只與雷諾數(shù)Re有關(guān)；在紊流流速較高（完全紊流）時(shí)，摩阻系數(shù)只是相對(duì)粗糙度的函數(shù)。

1963年Orville Gaither利用天然氣作為流動(dòng)介質(zhì)，對(duì)單相氣體在速度管中的流動(dòng)規(guī)律進(jìn)行了研究。實(shí)驗(yàn)測(cè)試的數(shù)據(jù)如下：氣流量1.37×104～5.31×104m3/d；天然氣相對(duì)密度0.62；平均溫度15～23℃；油管尺寸25.40，31.75 mm；井口壓力0.59～4.54 MPa；井底壓力1.52～5.55 MPa；壓降0.29～4.30 MPa；實(shí)驗(yàn)組數(shù)為25（管徑25.40 mm）+27（管徑31.75 mm）。

粗糙度取0.032 mm，利用Moody，Panhandle A，Panhandle B，AGA，Weymouth摩阻系數(shù)計(jì)算方法計(jì)算了單相氣體在油管和環(huán)空中的流動(dòng)壓力降，并與測(cè)試值進(jìn)行了對(duì)比（見(jiàn)圖1）。

圖1 不同方法下的壓降計(jì)算值與單相氣流測(cè)試壓降對(duì)比

不同摩阻系數(shù)計(jì)算方法計(jì)算的井筒壓降與測(cè)試壓降的誤差統(tǒng)計(jì)如表1所示。從表1可知，利用AGA方法計(jì)算的壓降誤差相對(duì)較小。

表1 單相氣體壓降計(jì)算誤差統(tǒng)計(jì) %

3 兩相流井筒壓降計(jì)算

Orville Gaither1963年利用天然氣和地層水作為流動(dòng)介質(zhì)，對(duì)氣液兩相流體在速度管中的流動(dòng)規(guī)律進(jìn)行了研究。

實(shí)驗(yàn)測(cè)試的數(shù)據(jù)參數(shù)范圍如表2所示。

表2 氣液兩相流動(dòng)實(shí)驗(yàn)參數(shù)范圍

利用文獻(xiàn)中的實(shí)驗(yàn)數(shù)據(jù)，將粗糙度取為0.032 mm，選取AGA方法計(jì)算摩阻系數(shù)，各兩相管流模型計(jì)算的壓降與測(cè)試值對(duì)比如圖2所示。

圖2 兩相管流模型計(jì)算壓力降與測(cè)試值對(duì)比

表3為各模型預(yù)測(cè)壓降與實(shí)測(cè)壓降對(duì)比得到的壓降平均誤差、壓降平均絕對(duì)誤差和壓降標(biāo)準(zhǔn)差。其中：壓降平均誤差表示模型的整體偏差，平均絕對(duì)誤差表示平均誤差的大小，標(biāo)準(zhǔn)差表示模型計(jì)算結(jié)果的離散程度［14-15］。從表3可知，Gray模型的各項(xiàng)誤差最小，在各模型中性能最佳；從平均相對(duì)誤差和標(biāo)準(zhǔn)差可知，Gray模型的對(duì)中性和離散程度均最好。計(jì)算壓降平均誤差相對(duì)較小的兩相管流模型依次是 Gray模型、Ansari模型、No-Slip模型。

表3 油管內(nèi)流動(dòng)壓降誤差統(tǒng)計(jì) %

4 實(shí)例分析與驗(yàn)證

以大牛地幾口典型的水平井為例。在3口水平井安裝速度管后進(jìn)行了流壓測(cè)試，其基礎(chǔ)數(shù)據(jù)見(jiàn)表4。

表4 水平井速度管流壓測(cè)試基礎(chǔ)數(shù)據(jù)

將壓力計(jì)下放點(diǎn)置于造斜點(diǎn)附近，完成水平井垂直段的流壓測(cè)試。利用流壓測(cè)試數(shù)據(jù)對(duì)兩相管流模型再次進(jìn)行了評(píng)價(jià)，測(cè)試壓降與計(jì)算值對(duì)比見(jiàn)表5。

表5 水平井速度管測(cè)試壓降與模型計(jì)算壓降對(duì)比 MPa

壓降誤差統(tǒng)計(jì)見(jiàn)表6。

表6 水平井速度管流動(dòng)壓降誤差統(tǒng)計(jì) %

本文主要以壓降平均絕對(duì)誤差為依據(jù)，并參考其他誤差，比較各壓降模型的準(zhǔn)確性。

從表5和表6可知：Gray模型壓降平均絕對(duì)誤差為13.70%，標(biāo)準(zhǔn)差為4.80%，是參與比較的模型中最佳的；其次為Hagedorn-Brown模型和Ansari模型。而常規(guī)管柱氣水兩相流壓降模型的評(píng)價(jià)結(jié)果表明，Hagedorn-Brown模型計(jì)算的誤差最小。這是由于該模型是唯一根據(jù)現(xiàn)場(chǎng)試驗(yàn)井?dāng)?shù)據(jù)推導(dǎo)而來(lái)，對(duì)于產(chǎn)水氣井條件，它明顯優(yōu)于其他基于實(shí)驗(yàn)數(shù)據(jù)的經(jīng)驗(yàn)方法［16］。

由此可見(jiàn)，速度管兩相流壓降模型評(píng)價(jià)結(jié)果與常規(guī)管柱模型評(píng)價(jià)結(jié)果存在一定的差異。出現(xiàn)這種差異的原因在于：在相同流量的情況下，速度管中氣液混合物流速比常規(guī)尺寸管柱流速高，流型也有一定的差異，井筒的滑脫損失程度有一定差異；加之小管徑流動(dòng)時(shí)的實(shí)驗(yàn)數(shù)據(jù)量相對(duì)較少，以上幾種經(jīng)驗(yàn)?zāi)Ｐ褪腔诖蠊苤鶜庖簝上嗔鲃?dòng)的實(shí)驗(yàn)數(shù)據(jù)建立起來(lái)的，未充分考慮管徑的影響，因而導(dǎo)致對(duì)速度管氣液兩相流壓降計(jì)算不準(zhǔn)，誤差較大。

為此，綜合考慮Orville Gaither的實(shí)驗(yàn)數(shù)據(jù)和大牛地速度管的流壓測(cè)試數(shù)據(jù)，以及各氣液兩相管流模型的評(píng)價(jià)結(jié)果，應(yīng)首選Gray模型對(duì)大牛地氣田速度管井筒壓力分布進(jìn)行計(jì)算，其次選用Ansari模型。

5 結(jié)論

1）速度管中氣液混合物流速較高，摩阻壓降較大，準(zhǔn)確選取摩阻系數(shù)計(jì)算方法是進(jìn)行速度管壓降計(jì)算的基礎(chǔ)。利用單相液流測(cè)壓數(shù)據(jù)進(jìn)行摩阻系數(shù)優(yōu)選表明，AGA方法計(jì)算結(jié)果與實(shí)際情況最吻合，誤差最小。

2）根據(jù)1963年Orville Gaither的實(shí)驗(yàn)數(shù)據(jù)和大牛地速度管的流壓測(cè)試數(shù)據(jù)，對(duì)8個(gè)常用氣液兩相管流壓降模型評(píng)價(jià)。結(jié)果表明，速度管排水采氣時(shí)，宜首選Gray模型，其次為Ansari模型進(jìn)行流壓計(jì)算。

［1］周瑞立，周艦，羅懿，等.低滲產(chǎn)水氣藏?cái)y液模型研究與應(yīng)用［J］.巖性油氣藏，2013，25（4）：123-128.

［2］田相雷，蔣海巖，岳艷如，等.氣舉排水采氣生產(chǎn)規(guī)律研究［J］.斷塊油氣田，2013，19（6）：760-763.

［3］白曉弘，田偉，田樹(shù)寶，等.低產(chǎn)積液氣井氣舉排水井筒流動(dòng)參數(shù)優(yōu)化［J］.斷塊油氣田，2014，21（1）：125-128.

［4］Gaither O D，Winkler H W，Kirkpatrick C V.Single-and two-phase fluid flow in small vertical conduits including annular configurations［J］.Journal of Petroleum Technology，1963，15（3）：309-320.

［5］周艦，羅懿，劉岳龍，等.大牛地水平井速度管排液工藝研究與應(yīng)用［J］.石油機(jī)械，2013，41（12）：102-105.

［6］趙彬彬，白曉弘，陳德見(jiàn)，等.速度管柱排水采氣效果評(píng)價(jià)及應(yīng)用新領(lǐng)域［J］.石油機(jī)械，2012，40（11）：62-65.

［7］林英松，蔣金寶，劉兆年，等.連續(xù)油管壓裂新技術(shù)［J］.斷塊油氣田，2008，15（2）：118-121.

［8］李閩，郭平，劉武，等.氣井連續(xù)攜液模型比較研究［J］.斷塊油氣田，2002，24（6）：39-41，91.

［9］趙界，李穎川，劉通，等.大牛地地區(qū)致密氣田氣井積液判斷新方法［J］.巖性油氣藏，2013，25（1）：122-125.

［10］代志華.氣井井筒積液診斷及排液效果分析［J］.化學(xué)工程與裝備，2011，23（3）：85-87.

［11］李士倫.天然氣工程［M］.北京：石油工業(yè)出版社，2008：109-117.

［12］喻西崇，馮叔初.多相管流中沿程摩阻系數(shù)分析［J］.油氣田地面工程，2001，20（2）：3-4.

［13］Moody L F.Friction factors for pipe flow［J］.Trans.Asme，1944，66（8）：671-684.

［14］廖開(kāi)貴，李穎川，楊志，等.產(chǎn)水氣藏氣液兩相管流動(dòng)態(tài)規(guī)律研究［J］.石油學(xué)報(bào)，2009，30（4）：607-612.

［15］劉通，任桂蓉，趙容懷.非環(huán)狀流氣井兩相流機(jī)理研究［J］.巖性油氣藏，2013，25（6）：103-106，132.

［16］李穎川，朱家富，秦勇.排水采氣井油管和環(huán)空兩相流壓降優(yōu)化模型［J］.石油學(xué)報(bào)，1999，20（2）：87-92.

（編輯 李宗華）

Evaluation and optimization of wellbore pressure drop model for drainage and gas recovery by velocity string

Tian Yun1,Wang Zhibin2,Li Yingchuan2,Bai Huifang1,Li Kezhi3
(1.School of Petroleum and Gas Engineering,Southwest Petroleum University,Chengdu 610500,China; 2.State key Laboratory of Oil and Gas Reservoir Geology and Exploitation,Southwest Petroleum University,Chengdu 610500, China;3.Research Institute of Engineering and Technology,Huabei Company,SINOPEC,Zhengzhou 450006,China)

Accurately calculating the pressure drop of two-phase flow in velocity string is the basis of optimizing design for drainage and gas recovery and analyzing production performance and drainage effect.In this paper,experiment data are from Orville Gather (1963)who carried out a single-phase and two-phase experiment in 25.4 and 31.75 mm velocity strings using natural gas/water as experimental medium.Firstly,the calculation method of friction factor is optimized.Secondly,the pressure drop of two-phase flow is evaluated and optimized.Finally,the optimum pressure drop of two-phase flow is verified by production data of three typical horizontal wells in Daniudi Gas Field after installing velocity string.Using the measured pressure data of pure liquid flow to optimize friction coefficient shows that the calculation of AGA method is the most consistent with actual results,whose error is the least（the average absolute error of pressure drop is 22.37%）.The evaluation of eight frequently-used two-phase flow pressure models shows that the priority selection to calculating wellbore pressure drop is Gray model and the second is Ansari model.This study can optimize a dependable calculation model of two-phase flow pressure drop for drainage and gas recovery by velocity strings and help to improve the technology design and the diagnosis level of drainage effect.

velocity string;drainage and gas recovery;wellbore;pressure drop;evaluation and optimization

國(guó)家科技重大專項(xiàng)課題“大牛地致密低滲氣田特殊結(jié)構(gòu)井滲流機(jī)理與采氣工藝研究”（2011ZX05045-05-01）

TE375

：A

10.6056/dkyqt201501029

2014-09-15；改回日期：2014-11-21。

田云，女，1989年生，在讀碩士研究生，研究方向?yàn)椴捎筒蓺饫碚撆c工程技術(shù)。E-mai：swputy@163.com。

田云，王志彬，李穎川，等.速度管排水采氣井筒壓降模型的評(píng)價(jià)及優(yōu)選［J］.斷塊油氣田，2015，22（1）：130-133.

Tian Yun，Wang Zhibin，Li Yingchuan，et al.Evaluation and optimization of wellbore pressure drop model for drainage and gas recovery by velocity string［J］.Fault-Block Oil&Gas Field，2015，22（1）：130-133.