非均质各向异性油藏水平井流入动态

大庆石油学院学报第36卷第3期2012年6月JOURNAL OF DAQING PET ROLEU M INSTIT UT E Vol.36No.3Jun.2012收稿日期:20120117;编辑:任志平基金项目:国家科技重大专项(2008ZX05031-02-03)作者简介:王海静(1983-),女,博士研究生,主要从事油气田地下工程力学方面的研究.非均质各向异性油藏水平井流入动态王海静1,薛世峰2,高存法1,仝兴华3( 1.南京航空航天大学航空宇航学院,江苏南京 210016; 2.中国石油大学储运与建筑工程学院,山东青岛266555; 3.山东大学威海分校机电与信息工程学院,山东威海 264209)摘 要:为分析渗透率非均质性对水平井流入动态的影响,根据质量和动量守恒定律,建立全耦合的非均质各向异性盒式油藏渗流与水平井筒变质量流的耦合模型;通过坐标变换建立各向异性油藏水平井非均质表皮计算模型,将非均质渗透率场转化为等效均质渗透率场和沿井筒变化的非均质表皮;应用源函数法和势的叠加原理得到油藏渗流模型的解,给出基于模型离散的顺序迭代解耦方法.结果表明:渗透率非均质性对水平井流量剖面和采油指数影响较大.水平井流量分布规律与近井渗透率变化规律基本相同;高渗区流量大,低渗区流量小;在等效渗透率相同的情况下,随着近井非均质性的增强,水平井平均采油指数逐渐降低,平均流量分布不均衡因数逐渐增大.该研究为非均质各向异性油藏水平井产能预测和完井参数优化提供了依据.关 键 词:非均质性;各向异性;油藏;水平井;耦合模型;渗透率场;源函数法;势的叠加原理;采油指数;不均衡因数;流量分布;等效渗透率中图分类号:TE355 文献标识码:A 文章编号:10001891(2012)030079070 引言随着世界范围内油气资源需求量的不断增长以及钻完井技术的日益完善,水平井成为高效开发油田的重要手段,在油气勘探开发中得到越来越广泛的应用.储层非均质性是影响水平井开发效果的重要因素,由此引发的储量动用不均衡、局部边底水过早突破是水平井开采过程中遇到的主要问题,越来越受到人们的关注.结合非均质油藏地质条件开展水平井流入动态研究,对于水平井产能预测以及完井方案设计具有重要意义.油藏内的流体经过水平井筒采出的过程按其流动规律可以划分为油藏渗流与井筒变质量管流两部分.1989年,Dikken B J [1]首次将井筒紊流压降引入水平井产能计算,分析水平井产能随长度的变化规律,指出在水平井油藏模拟和设计中应该考虑井筒压力损失的影响.随后,人们在水平井筒管流与油藏渗流的耦合方面开展大量研究.Ouyang L [2]和Penmatcha V R [3]等基于源函数法,建立非稳态的均质盒式油藏水平井筒)油藏流动耦合模型,.刘想平[4]、苏玉亮[5]等根据镜像反映原理和势的叠加原理,分别提出均质无限大油藏裸眼水平井和射孔水平井的稳态井筒)油藏耦合模型及求解方法,分析井筒压降对水平井稳态流量剖面的影响.李松泉等[6]应用Gr een 函数,建立非稳态的均质盒式气藏水平井筒)油藏流动耦合模型,分析不同阶段井筒压降对水平井流量剖面的影响.目前关于水平井流入动态的研究大多针对均质油藏,而对非均质油藏的研究鲜有报道.在建立非均质油藏)井筒耦合模型的基础上,笔者分析渗透率非均质性对水平井流入动态的影响,首先建立适用于非均质各向异性盒式油藏单相微可压缩流体不稳定流的油藏)井筒全耦合微分模型;然后提出渗透率非均质性定量评价方法,将非均质渗透率场转化为等效均质渗透率场和沿井筒变化的非均质表皮;采用源函数法求解油藏渗流方程,给出耦合模型的离散形式和解耦方法;通过算例验证文中模型及方法的准确性,分析渗透率非均质性对水平井压力/流量剖面以及产能的影响.网络出版时间：2012-07-23 09:24网络出版地址：/kcms/detail/23.1297.TE.20120723.0924.014.html1 油藏)井筒耦合模型图1 盒式油藏水平井系统盒式油藏中有1口水平井(见图1),井筒与油藏边界平行,长度为L,半径为r w .假设储层孔隙度<为常数,渗透率各向异性.流体单相微可压,作等温不稳定流动.以井筒趾端O 为原点,建立空间直角坐标系Oxyz,使Ox 轴与井筒轴线重合,Oy,Oz 平行于油藏边界.1.1 油藏渗流方程根据质量守恒定律和Darcy 渗流定律,建立油藏渗流控制方程:k x 52p 5x 2+k y 52p 5y 2+k z 52p 5z2=<L c t 5p 5t ,(1)式中:p 为流体压力;k x ,k y ,k z 分别为x ,y,z 方向的渗透率;<为孔隙度,无因次;L 为流体黏度;c t 为综合压缩系数;t 为时间.1.2 井筒变质量流方程定义/比流量0q 为单位长度井筒的径向流入量,单位为m 3/(s #m),则水平井筒变质量流的质量守恒方程和动量守恒方程为5Q 5t=-5(Q v)5x +Q q A ,(2)5(Q v)5t +55x A Q v 22=-5p 5x -Q f v 22d,(3)式中:Q 为流体密度;v 为井筒截面平均流速;A 为井筒横截面积;A 为动量修正因数[7];f 为考虑井壁入流影响的摩擦因数[8];d 为井筒内径.1.3 初边值条件假设整个系统初始压力均匀分布,即p |t=0=p 0,(4)式中:p 0为油藏初始压力.油藏外边界为定压或封闭,即p |#1=p 0,KL#àp #n |#2r =0.(5) 井筒趾端边界封闭,即v |x=0=0.(6) 井筒跟端边界定压或流量.当水平井以定压力生产时,有p |x=L =p wf ,(7)式中:L 为井筒长度;p wf 为井筒跟端压力.当水平井以定产量生产时,有Av |x=L =BQ,(8)式中:B 为体积因数;Q 为油井产量.1.4 耦合模型方程(1-3)构成非均质各向异性油藏中水平井的油藏)井筒耦合模型,含有3个未知量:压力p,比流量q 和井筒截面平均流速v.水平井油藏)井筒耦合问题就是求方程(1-3)满足初边值条件(4-8)的解.定解问题(1-8)由于具有较强的非线性,很难采用解析法求解,尤其在非均质油藏中.根据实际情况对模型进行适当简化处理后,可以采用半解析方法求解.大 庆 石 油 学 院 学 报 第36卷 2012年2非均质性定量评价图2 近井非均质影响区示意储层渗透率属区域化变量,其分布具有结构性和随机性双重特点.Dur lofsky L J 研究指出,渗透率非均质性对油井产能的影响可以分解为全局影响和近井影响[9].其中,全局影响体现为一个大尺度的全局等效均质渗透率K ={k x ,k y ,k z },近井影响体现为一个以井筒为轴心、半轴长为各方向渗透率相关长度一半的椭圆形影响区(见图2),并给出近井非均质表皮计算模型.然而,在推导非均质表皮时,Dur lofsky L J 直接应用H awkins 的各向同性油藏中表皮的计算方法,忽略渗透率各向异性对近井流动形态和压降的影响,严重影响该模型在各向异性油藏的精度.为建立各向异性条件下的非均质表皮计算模型,假设近井非均质影响椭圆与等压线重合,并作坐标变换将各向异性影响区转化为等效各向同性影响区:½y =y/I ani ,yz ,½z =z I ani ,yz ,(9)式中:I ani ,y z 为yz 平面的渗透率各向异性因数,I ani ,y z =K y /K z .取等效各向同性平面内椭圆形井筒的等效半径为两半轴长的平均值,根据H awkins 方法,建立各向异性油藏非均质表皮计算模型:S h =K y K h ,y -1l n1I ani ,y z +1r h ,y r w+r h ,yr w2+I 2ani ,y z -1,(10)式中:K h ,y 为近井非均质影响区y 方向的等效渗透率;r h ,y 为非均质影响椭圆的y 半轴;r w 为井筒半径.3 耦合模型解法3.1 模型简化源函数法具有计算简便、灵活性强的特点,与势的叠加原理相结合可以求解任意边界条件下平行六面体油藏的非稳态渗流问题[10-11],在均质油藏压力响应计算和试井分析中得到广泛应用[12-13].通过非均质性定量评价方法,将非均质渗透率场简化为一个常数等效渗透率和沿井筒变化的非均质表皮系数,从而采用源函数法求解油藏渗流模型.对于渗透率各向异性油藏,在应用源函数法之前,需通过坐标变换将其转化为等价的渗透率为K c =3K x K y K z 的各向同性油藏:x c =K y K zK x K y K zx ,y c =K x K zK x K y K zy ,z c =K x K yK x K y K zz ,(11)式中:上标/c0为坐标变换后参数.基于势的叠加原理,得到用源函数表示的定解问题(1)、(4-5)的井筒压力解:p w =p 0-L <cL c Q L 0QtqG(M c w f ,M c w ,t-S )d S d x -L q 2P K y K zS t ,(12)式中:G(M c w f ,M c w ,t -S )为等效各向同性油藏中井筒轴线上一点M c w 在井壁M c wf 处的瞬时源函数[11],随着油藏边界条件的不同而不同;S t 为由储层非均质性、污染、完井引起的总表皮系数,S t =S h +S d +S c ,无因次.其中,污染表皮系数S d 和完井表皮系数S c 的计算见文献[14-16].当井筒流量分布达到稳定时,式(12)写成p w =p 0-L<cL cQ Lq QtG(M c wf ,M c w ,t -S )d S d x -L q 2P K y K zS t .(13)此时,井筒内流动为稳定流,可以忽略质量守恒方程(2)和动量守恒方程(3)中的时间导数项,偏微分方程简化为常微分方程.若同时假设流体不可压缩,则式(2)和式(3)简化为第3期 王海静等:非均质各向异性油藏水平井流入动态d v d x=q,(14)d p w d x =-Q f v 22d -Q d d x A v 22.(15)3.2 模型离散简化后的油藏)井筒耦合模型(13-15)形式更简洁,各未知量之间的关系更清楚.由于模型方程具有非线性特征,仍然无法直接求解.因此,采用离散方法,将对连续变量的求解转化为对有限个离散变量的求解.将井筒离散为n 个微元L i (i =1,2,,,n),从趾端开始依次编号(见图1).将各井筒微元看作是流量均匀的线汇,取井壁中点M wf ,i (i =1,2,,,n)作为其压力计算点.根据式(13),建立各井筒微元的压力计算公式:p w ,i =p 0-6nj =1q j $L j <c t $L c jQtG(M c w f ,i ,L c j ,t-S )d S -L q i2P K y K zS t ,i (i =1,2,,,n),(16)式中:p w ,i 为第i 个井筒微元的压力;q j 为第j 个井筒微元的比流量;$L j 为第j 个井筒微元的长度;G(M c w f ,i ,L c j ,t -S )为等效各向同性油藏中井筒微元线汇L c j 在第i 个微元井壁中点M c w f ,i 处的瞬时源函数[11],随油藏边界条件的不同而不同;S t ,i 为第i 个井筒微元的总表皮.同样,根据式(6)、式(14-15),得到井筒变质量流的离散方程:v i =6nj =1q j $L j -0.5q i $L i /A (i =1,2,,,n),(17)p w,i -p w,j+1=-Qf 2d A 3q i +1(v 3i,j+1-v 3i+1)+A 3q i(v 3i -v 3i,j+1)-Q 2(A i v 2i -A i+1v 2i+1)(i =1,2,,,n -1),(18)式中:v i 为第i 个井筒微元中点处的截面平均流速;v i,i +1为第i 和i +1个井筒微元交点处的截面平均流速.v i,j+1=6nj =i+1q i $L j /A (i =1,2,,,n -1).(19) 内边界条件式(7-8)的离散形式:p w,l =p wf ,(20)6n i=1q i $L i =BQ.(21)控制方程(16-18)与边界条件式(20)或式(21)构成离散形式的油藏)井筒耦合模型.3.3 解耦方法油藏)井筒耦合模型(16-18)与式(20)或式(21),共3n 个方程,3n 个未知量,分别是各井筒微元的压力p i 、比流量q i 和截面速度v i (i =1,2,,,n),方程组有唯一解.鉴于模型的弱耦合特点,采用顺序迭代法求解.(1)在定压生产条件下,给压力赋初值p wf ,根据式(16-18)和式(20)循环计算比流量、截面速度和压力,直到满足计算精度为止;(2)在定产量生产条件下,给比流量赋初值BQ/L,根据式(16-18)和式(21)循环计算速度、压差、比流量,直到满足计算精度为止.4 算例分析一非均质垂向各向异性顶底封闭油藏,厚度为20m,初始压力为10MPa,孔隙度为0.3;渗透率场呈对数正态分布,等效水平渗透率为1L m 2,相关长度为200m,各向异性因数为10.地层原油黏度为10mPa #s,密度为850kg/m 3,体积因数为1.25,综合压缩系数为1.0@10-3MPa -1.水平井位于油藏中心,以800m 3/d 定产量生产,生产段长度为400m,井筒半径为0.1m,内径为0.121m,相对粗糙度为0.01.大 庆 石 油 学 院 学 报 第36卷 2012年定义近井渗透率非均质性定量评价指标/变异因数0C v 为标准差与平均值之比,即C v =1n 6ni =1(K i -ÇK )2/ÇK @100%,(22)式中:K i 为第i 个井筒微元的近井渗透率;ÇK 为近井平均渗透率;n 为井筒微元数目.4组对数正态分布,平均值为1L m 2,变异因数C v 分别为0,0.2,0.4和0.6的近井水平渗透率见表1.表1 近井水平渗透率分布距跟端距离/m 渗透率/10-3L m 2第1组第2组第3组第4组C v =0C v =0.2C v =0.4C v =0.6距跟端距离/m 渗透率/10-3L m 2第1组第2组第3组第4组C v =0C v =0.2C v =0.4C v =0.60~201000765594855200~22010001035133467620~401000695992644220~2401000675103258940~6010001037751597240~2601000123577390160~80100012007821470260~28010009011831212580~10010008279991826280~300100013117091618100~120100013039701954300~32010009805581018120~14010001300840693320~3401000918778455140~1601000800904688340~360100012441170908160~180100011512401941360~380100010335521129180~2001000881783463380~40010007091247450采用文中模型和方法,计算100d 时水平井筒的压力及流量分布,流量剖面已经稳定.井筒压力分布见图3.受管壁摩擦和径向入流影响,井筒内存在压力损失,趾端压力高,跟端压力低.由于油井定产量生产,4种渗透率情况下的跟、趾端压差基本相同,约为0.02MPa.4种情况下的井筒压力有所不同.第3组渗透率时井筒压力最高,第2组时井筒压力最低,其他组的井筒压力介于二者之间.井筒压力与变异因数之间无明显对应关系.井筒流量剖面见图4.均质情况下(第1组渗透率),井筒两端流量较高,中间流量较低,流量剖面呈现较明显的端部效应.由于井筒压力损失与生产压差相比较小,跟、趾端流量差异较小,流量剖面的均衡性较好.非均质情况下(第2、3、4组渗透率),水平井流量剖面出现明显波动,均衡性较差.高渗区流量大,低渗区流量小,流量分布规律与渗透率基本相同.渗透率变化幅度越大,则流量剖面波动越明显.对文中模型计算结果与有限元计算结果进行比较(见图4),2种方法得到的流量剖面变化规律基本一致.与有限元计算结果相比,文中模型计算结果略有不同,但最大误差不超过8%,可以满足油藏工程设计的需要.图3 井筒压力剖面图4 井筒流量剖面为定量评价近井非均质性对油井产能和井筒流量剖面均衡性的影响,定义采油指数PI 和不均衡因数C u :PI =Q/(p 0-p wf ),(23)第3期 王海静等:非均质各向异性油藏水平井流入动态C u =1n6n i =1(q i -q *i )2/½q @100%,(24)式中:q *i 为均质无限导流情况下第i 个井筒微元的比流量;½q 为各井筒微元的平均比流量.近井非均质性对产能的影响见图5,其中每个点代表同一变异因数下100组渗透率实现的采油指数平均值,误差线代表标准误.由图5可以看出,随着变异因数的增大,平均采油指数逐渐降低.当变异因数从0增加到2.0时,采油指数大约降低9%.此外,随着变异因数的增大,采油指数标准误逐渐增大,即产能的离散程度逐渐增强.近井非均质性对流量剖面均衡性的影响见图6,其中每个点代表同一变异因数下100组渗透率实现的不均衡因数平均值,误差线代表标准误.由图6可以看出,当变异因数为0时,即均质情况下,不均衡因数为0.012,仅由井筒压力损失引起,流量剖面的均衡性较好.随着变异因数的增大,标准误和平均流量不均衡因数逐渐增大,即随着非均质性的增强,平均流量不均衡性逐渐增强.图5 非均质性对产能的影响图6 非均质性对流量剖面均衡性的影响5 结论(1)建立非均质油藏渗流与井筒变质量流的耦合模型,提出基于非均质表皮和源函数的解耦方法.该模型考虑储层非均质性、渗透率各向异性、井筒压力损失、近井地层损害和完井工艺对产能的影响,适用于非均质各向异性盒式油藏水平井的产能评价与动态预测.(2)渗透率非均质性对水平井流量剖面和采油指数影响较大.水平井流量分布规律与近井渗透率变化规律基本相同.高渗区流量大,低渗区流量小.在等效渗透率相同的情况下,随着近井非均质性的增强,水平井平均采油指数逐渐降低,.(3)文中提出的非均质性定量评价方法为非均质各向异性油藏水平井产能预测提供了一个简单、快捷的途径;与数值模拟方法相比,该简化方法能够满足一般油藏工程设计的精度要求,但有待进一步改进.参考文献:[1] Dikken B J.Pressure d rop in horizontal wells and its effect on production perform ance [J].SPE J ou rnal of Petroleum T echnology,1990,42(11):426-1433.[2] Ouyan g L,Aziz K A.Sim plified approach to cou ple wellbore flow an d res ervoir in flow for arbitrary well configuration s[R].SPE48936,1998.[3] Penmatcha V R,Aziz K.A com preh ensive res ervoir/wellbore model for horizontal wells [R].SPE 39521,1998.[4] 刘想平,郭呈柱,蒋志祥,等.油层中渗流与水平井筒内流动的耦合模型[J].石油学报,1999,20(3):82-86.[5] 苏玉亮,张东,李明忠.油藏中渗流与水平井筒内流动的耦合数学模型[J].中国矿业大学学报,2007,36(6):752-758.[6] 李松泉,廉培庆,李秀生.水平井井筒和气藏耦合的非稳态模型[J ].西南石油大学学报:自然科学版,2009,31(1):53-57.[7] White F M.Flu id m ech anics [M].New York:McGraw-Hill Book Company,1986.[8] Ouyang L,Arbabi S,Aziz K.General wellbore flow model for horizontal,vertical,and slanted well com pletions [J].SPE Journal,1998,3(2):124-133.大 庆 石 油 学 院 学 报 第36卷 2012年第3期王海静等:非均质各向异性油藏水平井流入动态[9]Duilofsky L J.An approximate model for well productivity in heterog eneous porous media[J].M athematical Geology,2000,32(4):421-438.[10]Gringarten A C,Ramey J R,Henry J.Application of the P function to h eat conduction and fluid flow problem s[R].SPE3818,1971.[11]Gringarten A C,Ramey J R,H enry J.T he use of source and Green's fun ctions in s olving un steady-flow problems in res ervoirs[J].Society of Petroleum Engineers J ou rnal,1973,13(5):285-296.[12]曾昭英,徐魁生,冯凤萍,等.矩形油藏垂直裂缝试井数学模型及应用方法[J].大庆石油学院学报,2004,28(5):13-15.[13]王海涛,张烈辉.基于源函数法的油层部分射开顶底封闭各向异性的三重介质油藏不稳定渗流[J].大庆石油学院学报,2008,32(3):29-33.[14]Furu i K,Zhu D,H ill A D.A comprehen sive model of horizontal well completion performan ce[J].SPE Produ ction&Operations,2005,20(3):207-220.[15]Furu i K,Zhu D,H ill A D.A new skin-factor model for perforated horizontal wells[J].SPE Drilling&Completion,2008,23(3):205-215.[16]Yildiz T.Assessm ent of total skin factor in perforated wells[J].SPE Reservoir E valuation&Engineering,2006,9(1):61-76.(上接第45页)参考文献:[1]付广,孟庆芬,徐琴.乌尔逊凹陷南二段油气成藏与分布主控因素及有利区预测[J].吉林大学学报:地球科学版,2004,34(3):377-382,387.[2]陈守田,刘招君,崔凤林,等.海拉尔盆地含油气系统[J].吉林大学学报:地球科学版,2002,32(2):151-154.[3]张吉光,彭苏萍,张宝玺,等.乌尔逊)贝尔断陷油气藏类型与勘探方法探讨[J].石油勘探与开发,2002,29(3):48-50.[4]侯启军,冯子辉,霍秋立.海拉尔盆地乌尔逊凹陷石油运移模式与成藏期[J].地球科学:中国地质大学学报,2004,29(4):397-403.[5]霍秋立,汪振英,李敏,等.海拉尔盆地贝尔凹陷油源及油气运移研究[J].吉林大学学报:地球科学版,2006,36(3):377-383.[6]吴河勇,李子顺,冯子辉,等.海拉尔盆地乌尔逊)贝尔凹陷构造特征与油气成藏过程分析[J].石油学报,2006,27(S1):1- 6.[7]伍英,陈均亮,张莹.海拉尔)塔木察格盆地构造带与油气关系[J].大庆石油学院学报,2009,33(3):31-35.[8]刘志宏,万传彪,任延广,等.海拉尔盆地乌尔逊-贝尔凹陷的地质特征及油气成藏规律[J].吉林大学学报:地球科学版,2006,36(4):527-534.[9]董焕忠.海拉尔盆地乌尔逊凹陷北部晚期构造变形与油气聚集[J].大庆石油学院学报,2010,34(6):32-36.[10]张云峰,王磊,李晶.海拉尔盆地苏德尔特古潜山裂缝发育特征[J].大庆石油学院学报,2006,30(5):4- 6.[11]付晓飞,胡春明,李景伟.贝尔凹陷的形成机制及其油气分布规律[J].石油学报,2008,29(3):356-362.[12]姜茵华.海拉尔盆地乌尔逊凹陷南屯组沉积相分析[J].大庆石油学院学报,2011,35(4):12-16.[13]文全.海拉尔盆地下白垩统扇三角洲沉积特征[J].大庆石油学院学报,2011,35(3):17-21.[14]马中振,庞雄奇,王洪武,等.海拉尔盆地乌尔逊)贝尔凹陷断层控藏作用[J].西南石油大学学报:自然科学版,2009,31(6):27-31,205.[15]付广,刘云燕.贝尔凹陷/源断势0对油成藏与分布的控制作用[J].地质科学,2011,46(4):1068-1078.。