賈永祿 孫高飛 聶仁仕 黃小云 劉 浩.“油氣藏地質(zhì)及開發(fā)工程”國(guó)家重點(diǎn)實(shí)驗(yàn)室·西南石油大學(xué), 四川 成都 60500;.中國(guó)石油大慶油田有限責(zé)任公司試油試采分公司, 黑龍江 大慶 634

?

基于應(yīng)力敏感性的視均質(zhì)火山巖氣藏試井模型

賈永祿1孫高飛1聶仁仕1黃小云2劉 浩1
1.“油氣藏地質(zhì)及開發(fā)工程”國(guó)家重點(diǎn)實(shí)驗(yàn)室·西南石油大學(xué), 四川 成都 610500;2.中國(guó)石油大慶油田有限責(zé)任公司試油試采分公司, 黑龍江 大慶 163412

中國(guó)含油氣盆地中火山巖氣藏資源豐富,勘探潛力大?；鹕綆r儲(chǔ)層分布范圍廣、基質(zhì)物性差、滲透率低,后期改造作用強(qiáng)。由于氣體流動(dòng)過程中存在應(yīng)力敏感性,因此常規(guī)試井模型誤差較大。針對(duì)這一問題,考慮到模型存在很強(qiáng)的非線性,建立了基于應(yīng)力敏感性視均質(zhì)火山巖氣藏不穩(wěn)定試井模型,并重新定義部分參數(shù),通過攝動(dòng)變換技術(shù)、拉普拉斯變換、Stehfest數(shù)值反演等方法對(duì)模型求解。推導(dǎo)了無窮大邊界、封閉邊界下模型的半解析解,獲得試井圖版,在此基礎(chǔ)上進(jìn)一步研究了變形介質(zhì)滲透率模量等參數(shù)對(duì)壓力動(dòng)態(tài)曲線的影響。研究表明,基于應(yīng)力敏感性的試井曲線中壓力和壓力導(dǎo)數(shù)曲線均上翹,壓力導(dǎo)數(shù)曲線在0.5線之上;封閉邊界和應(yīng)力敏感性都會(huì)引起曲線抬升;氣藏儲(chǔ)層應(yīng)力敏感的存在使氣藏壓力降落不穩(wěn),產(chǎn)量遞減趨快,單井控制程度有限,因此在氣藏開采過程中采用水平井壓裂改造十分必要,同時(shí)應(yīng)保持氣藏內(nèi)部原始?jí)毫σ岳跉饩€(wěn)產(chǎn)和延長(zhǎng)開采時(shí)間。

火山巖氣藏;應(yīng)力敏感;視均質(zhì);試井模型;特征

0 前言

1 物理模型和數(shù)學(xué)模型

1.1 物理模型

應(yīng)力敏感性是指巖石骨架所受有效應(yīng)力的改變而引起巖石滲透率的變化程度[3]。FattI[7]、戈?duì)柌贾Z夫[8]指出巖石的孔隙度受壓力變化的影響很小,其變化可忽略,而滲透率受到的影響不可忽略。本文研究的火山巖氣藏中裂縫被碎屑物質(zhì)或者化學(xué)物質(zhì)充填并壓實(shí),這種裂縫對(duì)應(yīng)力的敏感性與基質(zhì)基本一致,實(shí)際上,被碎屑物充填了的裂縫已經(jīng)不是真正意義上的裂縫,而是具有裂縫形態(tài)的(微型)基質(zhì)巖石層[9]。

視均質(zhì)火山巖儲(chǔ)層既是儲(chǔ)集空間又是滲流通道,其物理模型的假設(shè)條件為:1)儲(chǔ)層為等厚無窮大或圓形封閉地層,忽略應(yīng)力敏感性引起儲(chǔ)層厚度的變化;2)視均質(zhì)介質(zhì)中的流體為單相可壓縮氣體,流動(dòng)服從等溫達(dá)西滲流;3)忽略重力和毛細(xì)管力的影響;4)氣井貫穿整個(gè)儲(chǔ)層并以恒定產(chǎn)量q生產(chǎn);5)考慮巖石滲透率應(yīng)力敏感性的影響。

1.2 數(shù)學(xué)模型

巖石變形后的滲透率k與初始滲透率ki滿足:

k(p)=kie-γ*(ψi-ψ)

(1)

引入擬體積系數(shù):

(2)

火山巖氣藏氣體擴(kuò)散方程為:

(3)

式中:k為滲透率,10-3μm2;p為氣藏壓力,MPa;Bp為擬體積系數(shù),無因次;Bg為氣體體積系數(shù),無因次;μ為流體黏度,mPa·s;Z為氣體偏差因子,無因次;ψ為真實(shí)氣體擬壓力,MPa2/(mPa·s);γ*為擬滲透率模量,(mPa·S)/MPa2;φi為孔隙度,無因次;Ct為儲(chǔ)層綜合壓縮系數(shù),1/MPa;下標(biāo)i表示原始狀態(tài)。

將式(3)無因次化,再結(jié)合初始條件和內(nèi)外邊界條件,組成應(yīng)力敏感性火山巖視均質(zhì)氣藏定產(chǎn)量生產(chǎn)下的不穩(wěn)定試井?dāng)?shù)學(xué)模型:

(4)

2 模型求解

擴(kuò)散方程(3)非線性,采用Pedrosa變換化簡(jiǎn)方程[10]:

(5)

應(yīng)用攝動(dòng)技術(shù)(PerturbationTechnique)變換式,取零階攝動(dòng)解:

(6)

對(duì)方程組(6)作tD→u的拉普拉斯變換,根據(jù)虛宗量的貝塞爾函數(shù)的性質(zhì)可得:

無窮大外邊界條件下的拉式空間解:

(7)

圓形封閉外邊界條件下的拉式空間解:

(8)

式中:I0、I1分別為0階與1階第一類虛宗量貝塞爾函數(shù);K0、K1分別為0階與1階第二類虛宗量貝塞爾函數(shù)。

對(duì)式(7)～(8)通過Stehfest數(shù)值反演[11],并代入Pedrosa變換函數(shù)得井底無因次擬壓力:

(9)

式中:S-1為Stehfest數(shù)值逆變換。

3 試井曲線特征分析

3.1 應(yīng)力敏感對(duì)無窮大邊界影響

圖1 無窮大外邊界下應(yīng)力敏感氣藏試井曲線

3.2 應(yīng)力敏感對(duì)圓形封閉邊界影響

值得注意的是,現(xiàn)場(chǎng)試井解釋中,用封閉邊界模型去擬合應(yīng)力敏感性引起的上翹曲線,會(huì)出現(xiàn)錯(cuò)誤,如果再用錯(cuò)誤封閉邊界半徑等數(shù)據(jù)去計(jì)算單井控制儲(chǔ)量等指標(biāo),結(jié)果誤差很大。

圖2 圓形封閉外邊界下應(yīng)力敏感氣藏試井曲線

3.3 封閉邊界對(duì)應(yīng)力敏感影響

火山巖氣藏儲(chǔ)層滲透率應(yīng)力敏感作用的存在使得氣藏壓力降落不穩(wěn),產(chǎn)量遞減加快,因此在氣藏開發(fā)中采用水平井壓裂改造十分必要[12]。

圖3 圓形封閉外邊界對(duì)應(yīng)力敏感性氣藏試井曲線的影響

4 結(jié)論

1)火山巖氣藏氣體與巖石變形之間的耦合關(guān)系模型能很好解釋應(yīng)力敏感性氣藏的滲流理論。通過基于應(yīng)力敏感性試井模型對(duì)原來試井模型的修正,可對(duì)火山巖應(yīng)力敏感性氣藏進(jìn)行較準(zhǔn)確的試井解釋,若不考慮這種應(yīng)力敏感特性，結(jié)果會(huì)產(chǎn)生較大偏差。

2)許多火山巖氣藏中的裂縫被碎屑物質(zhì)或化學(xué)物質(zhì)充填并壓實(shí),這種裂縫已經(jīng)不是真正意義上的裂縫,其試井曲線反映出均質(zhì)曲線特征,這種視均質(zhì)火山巖儲(chǔ)層介質(zhì)既是儲(chǔ)集空間又是滲流通道。

3)在現(xiàn)場(chǎng)試井解釋中,壓力導(dǎo)數(shù)曲線上翹,可以是應(yīng)力敏感性引起的,用封閉邊界模型擬合得到的邊界距離會(huì)偏小。

4)由于火山巖氣藏儲(chǔ)層滲透率應(yīng)力敏感性的作用使得氣藏壓力下降不穩(wěn),產(chǎn)量遞減趨快,并且單井控制程度有限,因此在氣藏開發(fā)中采用水平井壓裂改造十分必要,同時(shí)應(yīng)盡量保持氣藏內(nèi)部原始?jí)毫σ岳跉饩€(wěn)產(chǎn)和延長(zhǎng)開采時(shí)間。

[1] 易曉忠,張 河,胥 青,等.克拉美麗氣田火山巖氣藏試井分析及應(yīng)用[J].新疆石油天然氣,2010,6(3):37-42.YiXiaozhong,ZhangHe,XuQing,etal.AnalysesofKelameiliGasFieldLavaReservoirWellTestandItsApplication[J].XinjiangOil&Gas,2010,6(3):37-42.

[2] 李?？?黃小云,蔣凱軍.慶深氣田火山巖氣藏試井曲線特征分析[J].油氣井測(cè)試,2012,21(1):16-18.LiHaike,HuangXiaoyun,JiangKaijun.TheCharacteristicAnalysisofWellTestCurveforQingshenVolcanicGasReservoir.[J].WellTesting,2012,21(1):16-18.

[3] 劉之的,苗福全,羅曉芳,等.火山巖裂縫型儲(chǔ)層應(yīng)力敏感性實(shí)驗(yàn)研究[J].天然氣地球科學(xué),2012,23(2):208-212.LiuZhidi,MiaoFuquan,LuoXiaofang,etal.ExperimentofStressSensitivityforIgneousFracturedReservoir[J].NaturalGasGeoscience,2012,23(2):208-212.

[4] 霍凌婧.火山巖氣藏復(fù)雜滲流機(jī)理研究[D].北京:中國(guó)科學(xué)院研究生院(滲流流體力學(xué)研究所),2009.HuoLingjing.MechanismResearchforComplexVolcanicRockGasReservoirPercolation[D].Beijing:GraduateSchoolofChineseAcademyofSciences(SeepageFluidMechanicsResearchInstitute),2009.

[5] 薛定諤.多孔介質(zhì)中的滲流物理[M].王鴻勛,譯.北京:石油工業(yè)出版社,1982：26-33.SchrodingerAE.TranslatedPorousMediaFlowPhysics[M].WangHongxun,trans.Beijing:PetroleumIndustryPress,1982：26-33.

[6]GladkovSO.Gas-kineticModelofHeatConductionofHeterogeneousSubstances[J].TechnicalPhysics,2008,53(7):828-832.

[7]FattI.PoreVolumeCompressibilitiesofSandstoneReservoirsRocks[J].AAPG,1958,42(8):1924-1957.

[8] 戈?duì)柌贾Z夫.異常油田開發(fā)[M].張樹寶,譯.北京:石油工業(yè)出版社,1987：6-14.GorbunovAT.AbnormalOilfieldDevelopment[M].ZhangShubao,trans.Beijing:PetroleumIndustryPress,1987：6-14.

[9] 李傳亮.裂縫性油藏的應(yīng)力敏感性及產(chǎn)能特征[J].新疆石油地質(zhì),2008,29(1):72-75.LiChuanliang.StressSensitivityandProductivityofFracturedReservoirs[J].XinjiangPetroleumGeology,2008,29(1):72-75.

[10]PedrosaOA.PressureTransientResponseinStress-sensitiveFormations[C]//Paper15115PresentedattheSPECaliforniaRegionalMeeting,2-4April1986,Oakland,California,USA.NewYork:SPE,1986.

[11] 賈永祿,趙必榮.拉普拉斯變換及數(shù)值反演在試井分析中的應(yīng)用[J].天然氣工業(yè),1992,12(1):60-64.JiaYonglu,ZhaoBirong.ApplicationofLaplaceTransformandNumericalInversiontoWellTest-analysis[J].NaturalGasIndustry,1992,12(1):60-64.

[12] 郭 洋,楊勝來.我國(guó)火山巖油氣藏壓裂技術(shù)研究進(jìn)展[J].天然氣與石油,2012,30(2):49-51.GuoYang,YangShenglai.CurrentStatusandDevelopmentTrendofVolcanicReservoirFracturingTechnologyinChina[J].NaturalGasandOil,2012,30(2):49-51.

2015-03-28

國(guó)家自然科學(xué)青年基金資助項(xiàng)目(51304164);高等學(xué)校博士學(xué)科點(diǎn)專項(xiàng)科研基金新教師類資助課題(20135121120001)

賈永祿(1948-),男,四川南充人,教授，博士生導(dǎo)師,學(xué)士，從事油氣滲流理論和油氣藏工程科研工作。

10.3969/j.issn.1006-5539.2015.04.008