计算油藏相渗曲线的新方法及应用_黄祥峰

石 油 地 质 与 工 程

2013年1月           PETROLEUM GEOLOGY ANDENGINEERING        第27卷 第1期 

()文章编号:1673-8217201301-0053-03

计算油藏相渗曲线的新方法及应用

黄祥峰1,张光明1,郭俊磊2,王 鹤3,覃兆辉4

(长江大学石油工程学院,湖北荆州4渤海石油装备(天津)中成机械制造公司;1.34023;2.

中国石化胜利油田分公司黄河钻井四公司;中国石化华北分公司)3.4.

摘要:借助于多元线性拟合,根据油田生产动态数据,推导出油水相对渗透率与含水饱和度的定量关系式。实例分利用该计算方法得到的相渗曲线能够更好地反映油藏的整体特性,该计算方法可对油藏的含水上升规律析表明,

做出较准确预测。

关键词:相渗曲线;动态数据;含水饱和度;相对渗透率中图分类号:TE313        文献标识码:A

  油水两相相对渗透率资料是研究油水两相渗流

是油田开发参数计算、动态分析以及油藏数的基础,

值模拟研究不可或缺的重要资料。目前,油水两相相对渗透率曲线主要通过取心岩样在室内实验室测

1]

。但是,一方面由于储层非均质性的影响,岩心定[

       

Krobs-we

=aeKrw

)、()由(可得:12

()3

室内实验获取的相渗曲线无法准确地反映油藏的整另一方面,钻井取心带来的污染也给实验室体特性,

最终使得实验获取的相渗曲线的测定带来了偏差,

相对渗透率资料无法真实代表油藏的平均相对渗透率。目前研究相对渗透率比较常用的方法主要有利用测井资料估算、毛管压力曲线计算、水驱曲线计算、非稳态法计算以及相关经验公式等方法,但这些方法都有一定的局限性。本次研究将利用油藏生产物性参数以及实验得到的相渗曲线等资动态数据、

料,通过分流方程以及油水两相相对渗透率与含水饱和度的指数关系式计算相渗曲线。

Boo))=lnnna-bSw  l+l-1

fwsμw

3]

:又有含水饱和度与采出程度的关系[       R=

()4

wwi

1-Swi

)、()由(推导出:45

()5

oo))lnlnb(1-SwR+lna-bSw+-1=-i)i

fwwsμ

()6oo,令Y=ln)ln-1)X=R,A=+  

fwwsμ

,式中Y、1-SwB=lna-bSwX为复合参-b(i)i,

)数,即可将(式转变为:A、B为常数;6()X+B7       Y=A

对于一个特定的油藏,故Y、A、B为常数,X呈线性关系,根据生产动态数据对其进行线性拟合,可确定出A、最终可求解出a、将a、B的值,b;b代入()式可以确定与不同含水饱和度下油水两相相对3

渗透率比值。

1.2 相对渗透率曲线的确定

由文献可知油水两相相对渗透率的表达

]4-5

:式[

1 相渗曲线的公式推导

1.1 确定油水相对渗透率比值与含水饱和度的关系

根据油水两相渗流的达西定律,在不考虑重力和毛管力的情况下,地下水油比公式为:      M=

QwKrwo

=×QoKroμw

()1

因Bw≈1,分流方程关系式为:

w/w

= fws=/Qw/Bw+QoBo

()2

owro

1+×

Krowμ

绝大多数沉积岩油水两相相对渗透率比值与流

[2]

Cwi

ww)S    Krw=Krw(or)

1-Swi-Sor

()8

收稿日期:2012-09-20

作者简介:黄祥峰,在读研究生,主要从事油藏开发1985年生,方面的学习和研究。

体含水饱和度之间的函数关系如下:

·54·石 油 地 质 与 工 程           2013年 第1期

orwC

)oSw    Kro=Kro(i)

1-Sor-Swi

式中CCw为常数。o、

)、()由(可得:89

()9

适应性。由A、有a=3B推导出a、b的值,7357.91,)将a、式得到油水相对渗透率b=14.9425,b代入(3

比值和含水饱和度的关系,求出不同含水饱和度下油水两相相对渗透率的比值,然后利用(式进行11)得到A1=0.多元线性拟合,827,A2=-0.965,A3求出参数C将C=0.890;0.827,Cw=0.965;o=o、

)和(式中求出不同含水饱和度下的油Cw代入(89)水两相相对渗透率。

将理论计算的油水两相相对渗透率曲线与实际的油水两相相对渗透率曲线进行对比,见图1。总的来讲,理论计算和实验的水相相渗曲线吻合得很好,但理论的油相渗透率值比实验值大,说明实验的相渗曲线不能完全反映实际地层的油水两相相渗规律。显然,结合油田生产动态数据计算的两相相渗更考虑了曲线不仅反映了实际岩石和流体的性质,实际油藏的非均质性,通过油田实际生成动态数据得到的油水两相相渗曲线将更具有代表性

r1oorw

)=C)l   l-ggoKr1-Swor-Swi

KrSwio(i)SwSw)l   Cw+lgg

1-SwKrSi-Sorwor()10

Kr1SoorSw

),),令Y′=lX1=lX2ggKr1-Swor-Swi

ww),A1=Co,A2=-Cw,A3==lg

1-Swi-Sorrwo(i),式中Yl′、X1、X2为复合参数,A1、A2、gKrSw(or)]6-7

):则可将(式简化成[A3为常数,10

()′=A1X1+A2X2+A311     Y

在确定了油水两相相对渗透率的比值与含水饱利用(式进行多元线性拟合和度关系的前提下,11)最终求出C将C得出A1、A2、A3的值,Cw的值,o、o)、()和Cw代入(式,计算出不同含水饱和度下的89油水两相相对渗透率值。

根据以上分析,作如下总结:

)利用收集到的油藏动态数据和物性参数,(1

)对式(进行线性拟合,得到A、利用等式推7B的值,)将a、得到油水两相导出系数a、b的值,b代入式(3相对渗透率比值和含水饱和度的关系;

()求出不同含水饱和度下的油水相对渗透率2

利用(式进行多元线性拟合,得到A1、比值后,11)即可求出C将CA2、A3的值,Cw的值,Cw代入o、o、

())和(式中求出不同含水饱和度下的油水两相相89对渗透率。

图1 理论和实验相渗曲线对比

然后通  通过油藏动态数据回归得到相渗曲线,

)、()过(分别计算出不同含水饱和度下对应的含25水率和采出程度值,即可以得到不同采出程度时含根据实验所测的相渗曲线也可水上升规律。同时,

以预测油田含水上升规律,见图2。与通过实验得到的含水上升规律相比,理论计算的含水上升规律与油田实际的含水上升规律吻合的更好,故通过理论计算的含水上升规律曲线能对油田将来的含水上升规律进行更符合实况的预测。

2 实例应用

2.1 基本数据

某油藏的基本参数如下:地面原油密度0.854/原油体积系数1.含油面积5.cm,272,67km, g

平均有效厚度1平均有效孔隙度0.原油5.6m,120, 

·s,,粘度为2地层水粘度0.地层mPa32mPa·s  水的体积系数1.束缚水饱和度0.残余油饱和01,4,度0.实验得到的束缚水饱和度下的油相相对渗28,。透率1.原油地质储量40,20.42×10t

4 

3 认识和结论

()推导出了一种新的渗透率计算方法,其适用1

的范围更加广泛,对于油田的开发有重要的意义。()利用生产动态数据确定油藏相对渗透率曲2线,是把整个油藏作为一个系统,有效克服了油田取心污染对相渗透率测定带来的偏差,不仅反映了流体和岩石的性质,也同时反映了油藏的非均质性,其

2.2 相对渗透率曲线的处理和计算

)利用油藏动态数据和物性参数,对式(进行线7性拟合,得到系数A=-8.拟合9655,B=4.5513,相关系数0.表明此方法在该油田具有很好的9703,

                黄祥峰等.计算油藏相渗曲线的新方法及应用·55·

计算的油藏相渗曲线可以作为整个油藏相渗曲线的代表

——相关系数,数,无因次量;由相对渗透率曲线的特征a、b—决定,与岩石孔隙大小、分布及界面张力和润湿性等参数有关。

参考文献

[]]石油钻1J. 张玄奇.油水相对渗透率曲线的实验测定[

():采工艺,1994,16587-90.

[]北京:石油工业出版社,2M].1994: 何更生.油层物理[

247.

[]北京:石油工业出3M]. 杨胜之,魏俊之.油藏物理学[

版社,2004:239-249.

图2 预测含水率上升规律对比

[]等.用生产数据计算油藏相对渗4 王怒涛,陈浩,王陶,

]:透率曲线[西南石油学院学报,J.2005,27(5)26-27.

[]王世虎,姜泽菊.温度、界面张力对油水相对5 王玉斗,

]:油气田地面工程,渗透率的影响[J.2004,23(4)11,31.

[]等.利用动态数据计算相渗曲6 吕新东,冯文光,杨宇,

]():线的新方法[特种油气藏,J.2009,16565-66.[]等.利用相对渗透率曲线计算油田7 熊健,梁蕊,蔡洪,

],重庆科技学院学报(自然科学版)采收率[J.2011,13():384-85.

编辑:李金华  

()通过该方法计算出的相渗曲线可以对油藏3

的含水上升规律做出科学的预测。

符号注释

——地下产水量和地下产油量,/;——地Qw、Qo—maM—

——水相相对渗透率和油相相对渗下水油比,小数;KrKrw、o———水、·s;——地面脱气原油透率;油的粘度,mPao—o—μw、μρ/———残余油饱和度下水相相对渗透密度,cm-3;KrSgw(or)———束缚水含水饱和度下油相相对渗透率;率;KrSwo(i)——地上含水率,小数;——含水饱和度;小数;Sw—s—fw

——束缚水饱和度,——残余油饱和水,小数;小数;SwSi—or———原油体积系数,——地层水的体积系无因次量;Bo—Bw—

檳檳檳檳檳檳檳檳檳檳檳檳檳檳檳檳檳檳檳檳檳檳檳檳檳檳檳檳檳檳檳檳檳檳檳檳檳檳檳檳檳檳檳檳檳檳檳檳檳(上接第41页)

,onstrainedbrockhsics[J].TheLeadinEde  c   ypygg  

():2009,28121486-1491.

[]htaseohsi3avidC,SimonS,DavidC,etal.Ti D      -gggpy

:AVO]calinversionforreservoircharacterization[J.    ,():CSEGRecorder2010,35529-35. 

[]4ermanD M,CarlosT V.Detectionandsatialdeline G     -p

ationofthin-sandsedimentarseuenceswithoint     yqj stocasticrestackinversionofwelllosand3Dseismic        pg[],dataJ.RervoirEvaluation&Enineerinamlitude   pgg2010:246-264.

[]5inletonS.Theeffectsofseismicdataconditioninon S      gg 

restacksimultaneousimedanceinversion[J].The   pp,():LeadinEde2009,281772-781.gg 

[][]6analesLL.RandomnoisereductionJ.SEGExand C     -p

,():edAbstracts1984,3525-527. 

编辑:刘洪树  

[]7ulunaN,MaesanM,RoendeH H,Gatherflatten G   -yg 

[]():in.TheLeadinEdeJ.2007,26121538-1543. ggg []8wanH W,Velocitiesfromamlitudevariationswith S     p

[],():offsetJ.Geohsics2001,6661735-1743.py[]9anninA.Someasectsofamlituderecoverractical C     -gppp 

]beforeAVOandinversion[J.SEGExandedAb     -yp 

,stracts2010:393-396.[][]10ossCP,BealePL.SeismicoffsetbalancinJ.Ge R      -g

,():ohsics1994,59193-101.py

[],D11endrelJebeeH,PedersenR,etal.Estimation P    y

andinterretationofPandSimedancevolumesfrom        ppinversionofP-waveoffsetseismicdatasimultaneous      [],J.SEGExandedAbstracts2000:146-149.  p

Jan.2013             PETROLEUM GEOLOGYANDENGINEERING  

:;H;KewordsChaiwouareaonanchiformationLu  -pgyy ;;hcaoouformationreservoircharacteristicsdrocarboni   -gydentification

atherrocessinres39EffectsanalsisofCRPotimalon       -gpgpyp tackinversion 

(,BaoXiieetalExlorationandDevelomentInstitute      jpp,,,,HDainOilfieldLtdCo.PetroChinaDaineilonian qgqggjg )163712

arameterrestackElasticinversionfromseismicamli      -ppptudedatahasanextensivealicationinthelitholoandh         -ppgyy rediction.Inversiondrocarbonalorithmsaswellasthe      pgualitofinutdataoftenexhibitdirectimactsonitsaccu         -qypp 

,rac.Asthebasicdataofseismicinversiontherestack       ypualitofCRPatherdatahashueeffectsoninversionre         -qygg 

,sults.Basedontheofenhancininversionaccuracurose      ppgy theketechnoloiesofCRPatherdataconditioninare      yggg  valueofthisdataconditioninmethodisverresented.The       -gp ifiedbmatchinotimizedCRPdatawithAVOresonseof       ygpp  ,m,wellloeanwhiletheimrovementofinversionresults     gp,arewelldemonstratedbanalzinfromasectsofwavelets      yygp  /,SNandresidualerror.Exectedlandadesirableinver      -py,sionresulthavebeenachievedafterthissetofCRPdata         wasrestacksuccessfullaliedinAVAsimultaneousin      -pypp versioninTannanareaofTamuchaebasin.      g

:p;CrestackathersKewordsinversionRP;gotimal  py ;AVO,rocessinresonsereservoirrediction  pgpp

46Reservesevaluationtechnoloofcomlexfaulted-blocks    gyp inBohai 

(,CLiu WeionetalTianinBranchComanNOOC,   ygjpy ,)TanuTianin300452 ggj

HdrocarbonreservoirsinBohaiaremainlcomlex     yyp ,therocessofreservesevaluationfaulted-blocks.Durin    gp ,reatroblemshavebeenoccurredsuchasunitclassifica       -gp,tioninterfaceselectionandeffectivethicknessdetermina     -tion.Throuhseismicinterretationtechnoloandreser    -gpgy ,voirdescritiontechnolothesacialdistributionofmain      pgyp,sandbodieshasbeenmadeclearwhichlasoodfoundation        ygforeffectivethicknessdeterminationandoil-bearinarea     g ,thesametimeincombinationwithnewloinselection.At       gggressuretechnoloandcableformationtestininterretation    gygpp  ,,,datathelowerlimitofeffectivethicknessreservoirtesfluid      ypinterfacendeservoironnectivitaveeenonfirmed. a r c b cy h,Throuhsecificresearchasetofscientificreservesevaluation       gp,wmethodshavebeensummeduhichcanrovidereatuid        -ppgganceforreservesresearchofsimilaroilfields.      

:;;Kewordscomlexfaultedblockreservoirevaluation   py 

·Ⅲ·

;;oil-beariareaeffectivethicknessreservesaccurac  gy ermeabilit53Newmeasurementmethodsforrelativecurve     py itsalicationand  pp

(HuanXianfenetalColleeofPetroleumEnineer    -ggggg  ,Y,,H)inantzeUniversitJinzhouubei434023  ggygermeabilitOil-waterrelativedataisindisensable    pyp 

,warametersforcalculatindevelomenthichmakesd    -pgpy ,namicanalsisandnumericalreservoirstimulationsoitis       ycriticalermeabilittoobtainrelativedataeffectivel.The     pyy ermeabilitrelativecurvesobtainedblaboratorexeri   -pyyyp   roerlmentsinoilindustrcannotreflectreservoirinte       -ppyyg  resentritattributesandmeasureddeviationbecauseofthe       yp offormationheteroeneitandcoredurinollutionimact      gygpp ,corinThereforeanew methodshouldbeconsideredtore      -g.,solveit.Combinedwithfielddatatheresearcherformance      pbasedonfielderformancedataandlinearfittincanacuire        pgq therelationshibetweenoil-waterrelativeermeabilitand    ppy  ,watersaturation.Afterstudinfieldcasethisaerets     ygppg ustifiablerelativecurvebtheoreticalcalculaermeabilit    -jpyy  

,ttionwhichcanreflectreservoirinteritattributeshen     gy ,wbasedonitaterholducanberedicted.     pp 

:;Kewordsrelativecurvedaermeabiliterformance  -pypy  ;w;ermeabilittaatersaturationrelative  py

56Waterfloodindevelomenteffectevaluationofsandcon    -gpy  lomeratebottomwateroilreservoirofLei64faulted-block       g

(JianQinGraduateSchoolofNortheastPetroleum Uni    -g ,,H)versitDaineilonian163318yqggjg 

Lei64faulted-blockisahuemassiveblockofsand         gybottom waterthinoilreservoir.Bottom waterconlomerate    ginectionroducdevelomentwasadotedandahihlevel        -jpppgtionunderlow watercontenthasalwasbeenmaintained.       y,,,theastfewearshoweverduetoravitdrainaeFor      pygyg ,theenerdeclinedandthenthebottom waterinectioncan      gyj ,tnoteffectivelrelenishtheuerformationenerhe     ypppgy blockressureresentedroedndharleclined  d a p sppppy d,wseed.Throuhreservoircharacteristicsanalsisaterin    -pgyectiondnamiccharacteristicsanalsisandwaterinection      jyyj,effectevaluationareasonabledevelomentmannerhasbeen      provided.p

:;;KewordsLei64faulted-blockLiaohebasinwater  -y ;floodindevelomenteffectevaluation gp 

59Practiceandunderstandinaboutwaterinectionenhance     -gj mentofultra-lowermeabilitoilreservoirinYaointai      pyyg 

Gan Wenun(ChanlinOilProductionPlantofNorth    -jgg ,,,eastOilandGasBranchComanSinoecSonuanJilin     pypgy


© 2024 实用范文网 | 联系我们: webmaster# 6400.net.cn