張芨強(qiáng) 李曉平 羅 誠(chéng) 袁 淋 杜知洋

1.中海石油(中國(guó))有限公司湛江分公司， 廣東 湛江 524057；2.西南石油大學(xué)油氣藏地質(zhì)及開發(fā)工程國(guó)家重點(diǎn)實(shí)驗(yàn)室， 四川 成都 610500；3.中國(guó)石油川慶鉆探工程公司地質(zhì)勘探開發(fā)研究院， 四川 成都 610056;4.中國(guó)石化西南油氣分公司川東北采氣廠， 四川 閬中 637402

?

基于橢圓流模型計(jì)算壓裂后氣井產(chǎn)能新思路

張芨強(qiáng)1,2李曉平2羅 誠(chéng)3袁 淋4杜知洋2

1.中海石油(中國(guó))有限公司湛江分公司， 廣東 湛江 524057；2.西南石油大學(xué)油氣藏地質(zhì)及開發(fā)工程國(guó)家重點(diǎn)實(shí)驗(yàn)室， 四川 成都 610500；3.中國(guó)石油川慶鉆探工程公司地質(zhì)勘探開發(fā)研究院， 四川 成都 610056;4.中國(guó)石化西南油氣分公司川東北采氣廠， 四川 閬中 637402

通常可將低滲透氣藏中垂直裂縫井的流動(dòng)劃分為裂縫區(qū)域外的橢圓徑向流和裂縫內(nèi)的線性流?；诜€(wěn)定滲流理論和橢圓流模型，引入橢圓坐標(biāo)系下的拉梅系數(shù)和擬壓力，通過對(duì)影響氣井產(chǎn)能多種因素的考慮，建立了垂直裂縫井在有限導(dǎo)流下的產(chǎn)能預(yù)測(cè)新模型。新模型計(jì)算結(jié)果表明：?jiǎn)?dòng)壓力梯度和應(yīng)力敏感使氣井產(chǎn)能減??；滑脫效應(yīng)使氣井產(chǎn)能增加；裂縫導(dǎo)流能力和裂縫長(zhǎng)度的增加使氣井產(chǎn)能增加，但一定的裂縫長(zhǎng)度對(duì)應(yīng)一個(gè)最佳的裂縫導(dǎo)流能力；隨著地層污染帶半徑增大和滲透率降低，氣井產(chǎn)能呈下降趨勢(shì)且下降幅度也在增大。利用新模型公式計(jì)算出的無(wú)阻流量與實(shí)際產(chǎn)能測(cè)試結(jié)果誤差較小，說明推導(dǎo)的產(chǎn)能預(yù)測(cè)方程具有一定的可靠性。

低滲透氣藏；產(chǎn)能方程；橢圓流；拉梅系數(shù)；啟動(dòng)壓力梯度；滑脫效應(yīng)；應(yīng)力敏感

0 前言

低滲透氣藏滲流實(shí)驗(yàn)結(jié)果表明，氣體在啟動(dòng)壓力梯度[1]、應(yīng)力敏感[2]以及滑脫效應(yīng)[3]等因素的影響下會(huì)出現(xiàn)非線性滲流特征[4-6]。對(duì)于低滲透氣藏中的氣井，為了獲取較高的產(chǎn)能，必須實(shí)施壓裂改造。近些年來(lái)，計(jì)算壓裂后形成垂直裂縫井產(chǎn)能的方法主要有利用保角變換法計(jì)算[7-8]，利用橢圓流模型和雙線性流模型計(jì)算[9-11]，將儲(chǔ)層中氣體的流動(dòng)考慮為徑向流和線性流組合模型計(jì)算[12-13]，以及其他垂直裂縫井的產(chǎn)能計(jì)算方法[14-17]。其中，利用橢圓流模型的計(jì)算主要采用了等價(jià)發(fā)展矩形的思想，相當(dāng)于把平面橢圓形的流動(dòng)區(qū)域簡(jiǎn)化為矩形，這樣會(huì)導(dǎo)致計(jì)算結(jié)果的偏差。本文從實(shí)際流動(dòng)區(qū)域入手，在橢圓流模型的基礎(chǔ)上引入了橢圓坐標(biāo)系下的拉梅系數(shù)和擬壓力，同時(shí)采用在橢圓坐標(biāo)系下對(duì)橢圓積分的方法，重新研究了低滲透氣藏垂直裂縫井產(chǎn)能計(jì)算方法。

1 壓裂井滲流物理模型

圖1 氣井壓裂后滲流物理模型

氣井在壓裂后改變了地層中流體的滲流方式，形成了橢圓流動(dòng)，見圖1。對(duì)壓裂后的氣井做出如下假設(shè)：

1)壓裂后形成沿井眼對(duì)稱分布的有限導(dǎo)流垂直裂縫。

2)裂縫剖面為與氣層等高的矩形。

3)地層中存在裂縫外的橢圓流動(dòng)區(qū)域和裂縫內(nèi)的線性流動(dòng)區(qū)域。

4)流體為單相、均質(zhì)流體，滲流過程中等溫且無(wú)任何特殊的物理化學(xué)現(xiàn)象發(fā)生。

5)地層中存在污染。

2 壓裂井產(chǎn)能預(yù)測(cè)模型

2.1 地層橢圓流

直井壓裂后形成的垂直裂縫井會(huì)改變氣體在地層中的滲流方式，誘發(fā)以裂縫端點(diǎn)為焦點(diǎn)的平面二維橢圓流[18]。由幾何學(xué)知識(shí)可知直角坐標(biāo)和橢圓坐標(biāo)的關(guān)系為：

x=afcos(η)y=bfsin(η)

(1)

af=Lfch(ξ)bf=Lfsh(ξ)

(2)

式中：x、y分別為直角坐標(biāo)；η、ξ分別為橢圓坐標(biāo)；af、bf分別為橢圓的長(zhǎng)半軸、短半軸；Lf為裂縫半長(zhǎng)，m；ch(ξ)、sh(ξ)分別為ξ的雙曲余弦、雙曲正弦函數(shù)。

根據(jù)氣體的穩(wěn)定滲流理論，在橢圓坐標(biāo)系下考慮啟動(dòng)壓力梯度、應(yīng)力敏感以及滑脫效應(yīng)的氣體運(yùn)動(dòng)方程[5-6，19]為：

(3)

引入橢圓坐標(biāo)系下的拉梅系數(shù)(hξ、hη)和漢密爾頓算子p，帶入式(3)后對(duì)η在[0，2π]積分可化簡(jiǎn)為：

在橢圓坐標(biāo)系下，對(duì)產(chǎn)量公式兩邊同時(shí)積分，有：

(5)

式中：qsc為氣井壓裂后地面產(chǎn)量，m3/d；Bg為氣體體積系數(shù)，Bg=ZTpsc/(pTsc)；A為滲流面積，m2；h為氣層厚度，m。

聯(lián)立式(4)和(5)，化簡(jiǎn)可得：

式中：Z為氣體偏差因子；psc為標(biāo)況下氣體的壓力，0.101 325MPa；T為地層溫度，K;Tsc為標(biāo)況下氣體溫度，293K。

其中：S=(Ki/Ks-1)(ξs-ξf)，此時(shí)ξf對(duì)應(yīng)長(zhǎng)軸為L(zhǎng)f的橢圓，由式(2)可知ξf=0，而ξe可通過等壓橢圓族方程進(jìn)行求解[20-22]，可得ξe=ln(2re/Lf)。

式中：ξf為橢圓坐標(biāo)系下裂縫處坐標(biāo)；ξe為橢圓坐標(biāo)系下氣藏外邊界坐標(biāo)；ξs為橢圓坐標(biāo)系下污染帶坐標(biāo)；re為供給半徑，m；pf為裂縫兩端的壓力，MPa；S為污染帶的表皮因子；Ks為污染帶滲透率，10-3μm2。

2.2 裂縫區(qū)域線性流

氣體在裂縫中流動(dòng)的物理模型見圖2。氣體從裂縫中流動(dòng)到井底時(shí)，存在高速非達(dá)西效應(yīng)，F(xiàn)orchheimerPH[23]通過實(shí)驗(yàn)提出了以下二次方程描述高速非達(dá)西流動(dòng)：

(9)

(10)

(11)

式中：Kf為裂縫滲透率，10-3μm2；v1為流體在裂縫中的流速，m/d;Wf為裂縫寬度，m；βg為氣體紊流系數(shù)，m-1；Mair為氣體相對(duì)分子質(zhì)量；γg為氣體相對(duì)密度；R為摩爾氣體常數(shù)，0.008 471MPa·m3/(kmol·K)；ρg為氣體密度，kg/m3。

圖2 氣體在裂縫中流動(dòng)的物理模型

將式(10)和(11)帶入式(9)，對(duì)x從rw到Lf積分，p從pwf到pf積分，且μg和Z的值按平均壓力處的值計(jì)算，并轉(zhuǎn)化為礦場(chǎng)實(shí)用單位制得

最后聯(lián)立式(8)、(12)即得低滲透氣藏有限導(dǎo)流垂直裂縫井的產(chǎn)能預(yù)測(cè)模型：

(13)

對(duì)于方程組(13)，首先輸入相關(guān)的計(jì)算參數(shù)pe、pwf等，接下來(lái)給出一組pf值，將pf值分別帶入到式(13)的兩個(gè)方程中，利用牛頓下山法和二次方程求根公式法分別求解出pf和qsc對(duì)應(yīng)的兩條關(guān)系曲線，曲線交點(diǎn)上的值即為井底流壓pwf對(duì)應(yīng)的產(chǎn)量qsc，改變pwf的值，便可求得無(wú)阻流量qAOF等參數(shù)。

3 實(shí)例分析

3.1 產(chǎn)能計(jì)算

3.2 產(chǎn)能影響因素分析

根據(jù)3.1節(jié)的壓裂氣井的基本參數(shù)，分析影響氣井無(wú)阻流量的敏感性因素。

3.2.1 非達(dá)西效應(yīng)影響因素

啟動(dòng)壓力梯度、應(yīng)力敏感和滑脫效應(yīng)對(duì)氣井無(wú)阻流量的影響見圖3～5。從圖中可知，啟動(dòng)壓力梯度和應(yīng)力敏感與氣井無(wú)阻流量都成負(fù)相關(guān)，而滑脫效應(yīng)與氣井無(wú)阻流量成正相關(guān)。說明啟動(dòng)壓力梯度和應(yīng)力敏感會(huì)阻礙氣體的流動(dòng)，在氣體的滲流過程中產(chǎn)生一種“附加阻力”；滑脫效應(yīng)會(huì)促進(jìn)氣體的流動(dòng)，在氣體的滲流過程中產(chǎn)生一種“附加動(dòng)力”。

圖3 啟動(dòng)壓力梯度對(duì)氣井無(wú)阻流量的影響

圖4 應(yīng)力敏感對(duì)氣井無(wú)阻流量的影響

圖5 滑脫效應(yīng)對(duì)氣井無(wú)阻流量的影響

3.2.2 裂縫參數(shù)

裂縫參數(shù)對(duì)氣井無(wú)阻流量的影響見圖6。從圖6可知，裂縫長(zhǎng)度和裂縫導(dǎo)流能力對(duì)氣井無(wú)阻流量的影響相似，它們的增加都會(huì)使無(wú)阻流量增大；但在一定裂縫長(zhǎng)度下，氣井無(wú)阻流量隨著裂縫導(dǎo)流能力的增加會(huì)很快達(dá)到一個(gè)上限值。說明裂縫長(zhǎng)度和裂縫導(dǎo)流能力的增加增大了氣體流向井筒的通道；一定的裂縫長(zhǎng)度對(duì)應(yīng)一個(gè)最佳的裂縫導(dǎo)流能力。

圖6 裂縫參數(shù)對(duì)氣井無(wú)阻流量的影響

3.2.3 污染帶參數(shù)

污染帶參數(shù)對(duì)氣井無(wú)阻流量的影響見圖7。從圖7可知，隨著地層污染帶半徑增大和污染帶滲透率降低，氣井無(wú)阻流量呈下降趨勢(shì)且下降幅度增大。說明污染帶半徑的增大和污染帶滲透率的降低對(duì)氣體的滲流會(huì)產(chǎn)生更大的阻力；污染帶半徑越大，污染帶滲透率對(duì)無(wú)阻流量的影響就越大；污染帶滲透率越大，污染帶半徑對(duì)無(wú)阻流量的影響就越小。

圖7 污染帶參數(shù)對(duì)氣井無(wú)阻流量的影響

4 結(jié)論

1)建立了垂直裂縫井在有限導(dǎo)流下的產(chǎn)能預(yù)測(cè)模型，經(jīng)實(shí)例計(jì)算驗(yàn)證其可靠性較高，為求取壓裂氣井產(chǎn)能提供了一種新思路。

2)啟動(dòng)壓力梯度和應(yīng)力敏感對(duì)氣體滲流產(chǎn)生一種“附加阻力”，使氣井產(chǎn)能降低；而滑脫效應(yīng)產(chǎn)生一種“附加動(dòng)力”，使氣井產(chǎn)能增加。

3)裂縫導(dǎo)流能力和裂縫長(zhǎng)度的增加擴(kuò)大了氣體的流動(dòng)通道，使氣井產(chǎn)能增加；一定的裂縫長(zhǎng)度對(duì)應(yīng)一個(gè)最佳的裂縫導(dǎo)流能力。污染帶滲透率的減小和污染帶半徑的增大會(huì)阻礙氣體的流動(dòng)，使氣井產(chǎn)能降低；污染帶半徑越大，污染帶滲透率對(duì)氣井產(chǎn)能的影響越大；污染帶滲透率越小，污染帶半徑對(duì)氣井產(chǎn)能的影響越大。

[1] 王曉琴,吳 聚,冉 艷,等. 非線性滲流對(duì)異常高壓氣藏產(chǎn)能的影響[J].巖性油氣藏,2012,24(4):125-127.WangXiaoqin,WuJu,RanYan,etal.InfluenceofNon-linerFlowonProductivityofAbnormalHighPressureGasReservoir[J].LithologicReservoirs, 2012, 24 (4): 125-127.

[2] 樊懷才,李曉平,竇天財(cái)，等.應(yīng)力敏感效應(yīng)的氣井流量動(dòng)態(tài)特征研究[J].巖性油氣藏,2010,22(4):130-133.FanHuaicai,LiXiaoping,DouTiancai,etal.StudyofStressSensitivityEffectonFlowDynamicFeaturesofGasWells[J].LithologicReservoirs, 2010, 22 (4): 130-133.

[3] 羅瑞蘭,程林松,朱華銀,等.研究低滲氣藏氣體滑脫效應(yīng)需注意的問題[J].天然氣工業(yè),2007,27(4):92-94.LuoRuilan,ChengLinsong,ZhuHuayin,etal.ProblemsontheStudyofSlippageEffectinLow-permeabilityGasReservoirs[J].NaturalGasIndustry, 2007, 27 (4): 92-94.

[4] 張芨強(qiáng),李曉平,袁 淋,等.非達(dá)西滲流對(duì)低滲透氣藏氣水同產(chǎn)水平井產(chǎn)能的影響[J].巖性油氣藏,2014,26(6):120-125.ZhangJiqiang,LiXiaoping,YuanLin,etal.InfluenceofNon-DarcyFlowonDeliverabilityofGas-waterProducingHorizontalWellinLowPermeabilityGasReservoirs[J].LithologicReservoirs, 2014, 26 (6): 120-125.

[5]KlinkenbergLJ.ThePermeabilityofPorousMediatoLiquidsandGases[J].APIDrillingandProductionPractice, 1941: 200-213.

[6]GeorgeDV,StanleyCJ.ApplicationofStress-DependentRockPropertiesinReservoirStudies[C]//Paper86979PresentedattheSPEInternationalThermalOperationsandHeavyOilSymposiumandWesternRegionalMeeting,16-18March2004，Bakersfield,California,USA.NewYork:SPE,2004.

[7] 汪永利,蔣廷學(xué),曾 斌.氣井壓裂后穩(wěn)態(tài)產(chǎn)能的計(jì)算[J].石油學(xué)報(bào),2003,24(4):65- 68.WangYongli,JiangTingxue,ZengBin.ProductivityPerformancesofHydraulicallyFracturedGasWell[J].ActaPetroleiSinica, 2003, 24 (4): 65-68.

[8] 蔣廷學(xué),單文文,楊艷麗.垂直裂縫井穩(wěn)態(tài)產(chǎn)能的計(jì)算[J].石油勘探與開發(fā),2001,28(2):61- 63.JiangTingxue,ShanWenwen,YangYanli.TheCalculationofStableProductionCapabilityofVerticallyFracturedWell[J].PetroleumExplorationandDevelopment, 2001, 28 (2): 61-63.

[9] 李 生,李 霞,曾志林,等.低滲透油藏垂直裂縫井產(chǎn)能評(píng)價(jià)[J].大慶石油地質(zhì)與開發(fā),2005,24(1):54-56.LiSheng,LiXia,ZengZhilin,etal.WellProductionEvaluationofVerticalFractureinLowPermeableReservoir[J].PetroleumGeology&OilfieldDevelopmentinDaqing, 2005, 24 (1): 54-56.

[10] 尹洪軍,劉 宇,付春權(quán).低滲透油藏壓裂井產(chǎn)能分析[J].特種油氣藏,2005,12(2):55- 56.YinHongjun,LiuYu,FuChunquan.ProductivityAnalysisforFracturedWellinLowPermeabilityReservoir[J].SpecialOilandGasReservoirs, 2005, 12 (2): 55-56.

[11] 熊 健,于路均,郭 平.非線性滲流低滲氣藏壓裂井的產(chǎn)能方程[J].天然氣與石油,2012,30(6):42-45.XiongJian,YuLujun,GuoPing.AnalysisonProductivityEquationofFracturedWellinLowPermeabilityGasReservoirwithNon-linearSeepage[J].NaturalGasandOil, 2012, 30 (6): 42-45.

[12] 楊正明,張 松,張訓(xùn)華,等.氣井壓后穩(wěn)態(tài)產(chǎn)能公式和壓裂數(shù)值模擬研究[J].天然氣工業(yè),2003,23(4):74-76.YangZhengming,ZhangSong,ZhangXunhua,etal.TheSteady-stateProductivityFormulaAfterFracturingforGasWellsandFracturingNumericalSimulation[J].NaturalGasIndustry, 2003, 23 (4): 74-76.

[13] 熊 健,曾 山,王紹平. 低滲透油藏變導(dǎo)流垂直裂縫井產(chǎn)能模型[J].巖性油氣藏,2013,25(6):122-126.XiongJian,ZengShan,WangShaoping.AProductivityModelofVerticallyFracturedWellwithVaryingConductivityforLowPermeabilityReservoirs[J].LithologicReservoirs, 2013, 25 (6): 122-126.

[14] 熊 健,劉向君,陳 朕. 低滲氣藏壓裂井動(dòng)態(tài)產(chǎn)能預(yù)測(cè)模型研究[J].巖性油氣藏,2013,25(2):82-85.XiongJian,LiuXiangjun,ChenZhen.DynamicProductivityPredictionModelforFracturingWellsinLowPermeabilityGasReservoir[J].LithologicReservoirs, 2013, 25 (2): 82-85.

[15] 羅天雨,趙金洲,郭建春.求取壓裂后氣井產(chǎn)能的橢圓流方法[J].天然氣工業(yè),2005,25(10):94-96.LuoTianyu,ZhaoJinzhou,GuoJianchun.EllipticalFlowMethodtoCalculateProductivityofGasWellsafterFracturing[J].NaturalGasIndustry, 2005, 25 (10): 94-96.

[16] 岳建偉,段永剛,青紹學(xué),等.含多條垂直裂縫的水平壓裂氣井產(chǎn)能研究[J].天然氣工業(yè),2004,24(10):102-104.YueJianwei,DuanYonggang,QingShaoxue,etal.StudyonProductionPerformanceofFracturedHorizontalGasWellswithSeveralVerticalFractures[J].NaturalGasIndustry, 2004, 24 (10): 102-104.

[17] 宋軍政,郭建春.計(jì)算氣井壓后產(chǎn)能新方法的應(yīng)用研究[J].鉆采工藝,2005,28(5):47- 49.SongJunzheng,GuoJianchun.TheAppliedResearchonaNewMethodofCalculatingDeliverabilityofFracturedGasWell[J].Drilling&ProductionTechnology, 2005, 28 (5): 47-49.

[18] 劉慈群.單相和兩相流體多維滲流問題[A].第十屆全國(guó)水動(dòng)力學(xué)學(xué)術(shù)會(huì)議論文集[C].北京:海洋出版社,1996,439-445.LiuCiqun.MultidimentionalSeepageProblemofSingle-phaseandTwo-phaseFlow[A].ProceedingsoftheTenthNationalCongressonFluidMechanics[C].Beijing:ChinaOceanPress, 1996, 439-445.

[19] 李曉平.地下油氣滲流力學(xué)[M].北京:石油工業(yè)出版社,2007:27-28.LiXiaoping.UndergroundOilandGasSeepageMechanics[M].Beijing:PetroleumIndustryPress, 2007: 27-28.

[20] 周 琦,姜漢橋,李志剛.低滲透氣藏垂直裂縫井產(chǎn)能計(jì)算[J].油氣井測(cè)試,2009,18(2):8-11.ZhouQi,JiangHanqiao,LiZhigang.DeliverabilityCalculationforGasWellswithVerticalFractureinLowPermeabilityReservoir[J].WellTesting, 2009, 18 (2): 8-11.

[21] 宋付權(quán),劉慈群,張盛宗.低滲透油藏中水平井的產(chǎn)能公式分析[J].大慶石油地質(zhì)與開發(fā),1999,18(3):33-35.SongFuquan,LiuCiqun,ZhangShengzong.AnalysisonProductivityFormulaofHorizontalWellsinLowPermeabilityReservoir[J].PetroleumGeology&OilfieldDevelopmentinDaqing, 1999, 18 (3): 33-35.

[22] 劉英憲,劉宗賓,周海燕,等.基于橢圓流模型的變形介質(zhì)油藏水平井產(chǎn)能研究[J].復(fù)雜油氣藏,2012,5(1):55-57.LiuYingxian,LiuZongbin,ZhouHaiyan,etal.StudyontheProductivityofHorizontalWellsinDeformedMediaReservoirBasedonEllipticalFlowTheory[J].ComplexHydrocarbonReservoirs, 2012, 5 (1): 55-57.

[23]ForchheimerPH.WasserbewegundurchBoden[J].Zeitsch-riftdesVereinesDeutscherIngenieure, 1901, 49: 1736-1749 & 50: 1781-1788.

2015-02-15

國(guó)家杰出青年科學(xué)基金項(xiàng)目“油氣滲流力學(xué)”(51125019)

張芨強(qiáng)(1990-)，男，四川南充人，碩士研究生，從事油藏工程及油田開發(fā)研究。

10.3969/j.issn.1006-5539.2015.06.011