一节油井流入动态IPR曲线课件

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,*,*,第一节 油井流入动态（IPR曲线）,教学目的：,掌握油井流入动态、采油指数等相关定义；并掌握单相流体流动、油气两相渗流、单相与油气两相渗流同时存在、油气水三相以及多油层情况下油井流入动态的绘制方法。,教学重点、难点：,教学重点,1、油井流入动态的定义以及计算方法,2、不同条件下油井流入动态的计算,第一节 油井流入动态（IPR曲线）教学目的：教学重点、难,1,教学难点,教法说明：,教学内容：,1、单相与两相渗流同时存在时油井流入动态的计算,2、油气水三相流动时油井流入动态的计算,课堂讲授并辅助以多媒体课件展示相关的数据和图表。,单相液体的流入动态,油气两相渗流时的流入动态,p,r,p,b,p,wf,时的流入动态,油气水三相流入动态,多层油藏油井流入动态,教学难点教法说明：教学内容：1、单相与两相渗流同时存在时油井,2,油井流入动态,：,油井产量(,q,o,),与井底流动压力(,p,wf,),的关系，反映了油藏向该井供油的能力。,基本概念,油井流入动态曲线,：,表示产量与流压关系的曲线，简称IPR曲线。,I,nflow,P,erformance,R,elationship,Curve,油井流入动态：基本概念油井流入动态曲线：,3,图1-1 典型的流入动态曲线,IPR曲线基本形状与油藏,驱动类型,有关。即使在同一驱动方式下，还将取决于,油藏压力,、,油层厚度,、,渗透率,及,流体物理性质,等。,p,r,q,o,max,图1-1 典型的流入动态曲线IPR曲线基本形状与油藏驱动类型,4,油井生产系统组成,油层到井底的流动,(地层渗流),井底到井口的流动,(井筒多相管流),井口到分离器,(地面水平或倾斜管流),油井生产的三个基本流动过程,气液两相流基本理论,油井流入动态,油井生产系统组成油层到井底的流动井底到井口的流动井口到分离器,5,一、,单相液体流入动态,供给边缘压力不变、圆形地层中心一口井,的产量公式为：,（1-1）,圆形封闭油藏、拟稳态条件下,产量公式为：,（1-2）,一、单相液体流入动态供给边缘压力不变、圆形地层中心一口井,6,图1-2 泄油面积形状与油井的位置系数,对于非圆形封闭泄油面积的油井产量公式，可根据泄油面积和油井位置进行校正。,图1-2 泄油面积形状与油井的位置系数对于非圆形封闭泄油面,7,单相流动时，油层物性及流体性质基本不随压力,变化。,采油指数,可,定义为:,单位生产压差下的油井产油量，是反映油层性质、厚度、流体参数、完井条件及泄油面积等与产量之间的关系的综合指标。,生产压差,直线型,单相流动时，油层物性及流体性质基本不随压力采油指数可定义为:,8,采油指数,J,的获得：,试井资料：,测得35个稳定工作制度下的产量及其流压，便可绘制该井的实测IPR曲线，取,其斜率的负倒数,油藏参数计算,对于单相液体流动的直线型IPR曲线，采油指数可定义为,产油量与生产压差之比,，也可定义为,每增加单位生产压差时，油井产量的增加值,，或,IPR曲线斜率的负倒数,。,注意事项：,因此，对于具有非直线型IPR曲线的油井，在使用采油指数时，应该说明相应的流动压力，不能简单地用某一流压下的采油指数来直接推算不同流压下的产量。,采油指数J的获得：试井资料：测得35个稳定工作制度下的产量,9,当油井产量很高时，井底附近将出现,非达西渗流,：,胶结地层的,紊流速度系数,：,非胶结地层,紊流速度系数,：,C,、,D,值也可用试井资料获取,当油井产量很高时，井底附近将出现非达西渗流：胶结地层的紊流速,10,二、,油气两相渗流时的流入动态,o,、B,o,、K,ro,都是压力的函数。用上述方法绘制IPR曲线十分繁琐。,通常结合生产资料来绘制,IPR,曲线。,平面径向流，直井,油气两相渗流,时油井产量公式为：,(一)垂直井油气两相渗流时的流入动态,二、油气两相渗流时的流入动态o、Bo、Kro都是压力的,11,1.Vogel 方法,(1968),假设条件,：,a.圆形封闭油藏，油井位于中心；,b.均质油层，含水饱和度恒定；,c.忽略重力影响；,d.忽略岩石和水的压缩性；,e.油、气组成及平衡不变；,f.油、气两相的压力相同；,g.拟稳态下流动，在给定的某一瞬间，各点的脱气原油流量相同。,数值模拟结果的总结,1.Vogel 方法(1968)假设条件：数值模拟结果的总,12,归一化曲线,归一化曲线,13,Vogel方程,经典方程,Vogel方程经典方程,14,a.计算,c.根据给定的流压及计算的相应产量绘制IPR曲线,b.给定不同流压，计算相应的产量：,、已知地层压力和一个工作点(,q,o,(test),p,wf,(test),),利用Vogel方程绘制IPR曲线的步骤,a.计算c.根据给定的流压及计算的相应产量绘制IPR曲线b,15,、已知两个工作点，油藏压力未知,a.油藏平均压力的确定,：已知或利用两组,q,o,p,wf,测试计算，即,b.计算,c.由流入动态关系式计算相关参数,、已知两个工作点，油藏压力未知a.油藏平均压力的确定：已,16,图2-4 计算的溶解气驱油藏油井IPR曲线,1-用测试点按直线外推；2-计算机计算值；3-用Vogel方程计算值,Vogel曲线与数值模拟IPR曲线的对比,图2-4 计算的溶解气驱油藏油井IPR曲线1-用测试点按,17,对比结果：,按Vogel方程计算的IPR曲线，最大误差出现在用小生产压差下的测试资料来预测最大产量时，但一般误差低于5。,虽然随着采出程度的增加到开采末期误差上升到20左右，但其绝对值却很小,。,如果用测试点的资料按直线外推，最大误差可达 7080，只是在开采末期约30%。,采出程度,N,p,对油井流入动态影响大，而,kh/、B,o,、k、S,o,等对其影响不大。,对比结果：按Vogel方程计算的IPR曲线，最大误差出现在用,18,2.费特柯维奇方法,溶解气驱油藏,假设(,k,ro,/,o,B,o,)与压力,p,成线性关系，则,其中，,2.费特柯维奇方法溶解气驱油藏假设(kro/oBo)与压力,19,所以：,当 时：,令：,费特柯维奇,基本方程,所以：当 时：令：费特柯维奇,20,3.不完善井Vogel方程的修正,油水井的不完善性：,射孔完成,打开性质不完善,；,未全部钻穿油层,打开程度不完善,；,打开程度和打开性质,双重不完善,；,在钻井或修井过程中油层受到损害或进行酸化、压裂等措施，从而改变油井的完善性。,3.不完善井Vogel方程的修正油水井的不完善性：射孔完成,21,图1-5 完善井和非完善井周围,的压力分布示意图,图1-5 完善井和非完善井周围,22,油井的,流动效率,FE,：,油井的理想生产压差与实际生产压差之比,为“正”称“正”表皮，油井不完善；,为“负”称“负”表皮，油井超完善。,油井的流动效率FE：油井的理想生产压差与实际生产压差之比为“,23,完善井,非完善井,令：,非完善井表皮附加压力降,于是,完善井非完善井令：非完善井表皮附加压力降于是,24,表皮系数或井壁阻力系数,S,油层受污染的或不完善井，,完善井,增产措施后的超完善井，,表皮系数,S,通常由试井方法获得,表皮系数或井壁阻力系数S油层受污染的或不完善井，完善井,增产,25,利用流动效率计算直井流入动态的方法,Standing方法(,1970),(,FE=0.5,1.5,),图1-6 FE 1时的无因次IPR曲线(standing IPR曲线),利用流动效率计算直井流入动态的方法Standing方法(,26,a.根据已知,p,r,和,p,wf,计算在FE=1时最大产量,standing方法计算不完善井IPR曲线的步骤：,b.预测不同流压下的产量,c.根据计算结果绘制IPR曲线,a.根据已知pr和pwf计算在FE=1时最大产量standi,27,Harrison方法,(,FE=1,2.5,),图1-7 Harrison无因次IPR曲线(FE1),Harrison方法(FE=1 2.5)图1-7,28,Harrison方法可用来计算高流动效率井的IPR曲线和预测低流压下的产量。其计算步骤如下：,a.计算,FE,=1时的,q,o,max,先求,p,wf,/,p,r,，然后查图1-7中对应的FE曲线上的相应值,q,o,/,q,omax,(FE=1)，,则,b.计算不同流压下的产量,c.根据计算结果绘制IPR曲线,d.求FE对应的最大产量，即,p,wf,=0时的产量,Harrison方法可用来计算高流动效率井的IPR曲,29,(二)斜井和水平井的IPR曲线,1990年，Cheng对溶解气驱油藏中斜井和水平井进行了数值模拟，并用回归的方法得到了类似Vogel方程的不同井斜角井的IPR回归方程：,p,=,p,wf,/,p,r,；q=,q,o,/,q,o,max,；A、B、C为取决于井斜角的系数,优点：,使用简单，仅需一组测点，便可得IPR曲线,缺点：,方程没有归一化，,(二)斜井和水平井的IPR曲线1990年，Cheng对溶解,30,1989年，Bendakhlia等用两种三维三相黑油模拟器研究了多种情况下溶解气驱油藏中水平井的流入动态关系。得到了不同条件下IPR曲线。,图1-8 拟合的IPR曲线与实际曲线的对比,_拟合的IPR曲线,实际曲线,1989年，Bendakhlia等用两种三维三相黑油模拟器研,31,曲线表明：,早期的IPR曲线近似于直线,随着采出程度增加,曲度增加,接近衰竭时曲度稍有减小。,Bendakhlia,建议用以下公式来拟合IPR曲线图：,曲线表明：早期的IPR曲线近似于直线,随着采出程度增加,曲度,32,图1-9 参数,v、n,与采出程度之间的关系,图1-9 参数v、n与采出程度之间的关系,33,IPR曲线的应用,油井流入动态反映了油藏向该井供油的能力。,根据测试资料确定IPR曲线。,根据IPR曲线确定流压和产量的对应关系。,p,r,q,o,max,IPR曲线的应用油井流入动态反映了油藏向该井供油的能力。根据,34,三、,p,r,p,b,p,wf,时的流入动态,（1）基本公式,当油藏压力高于饱和压力，而流动压力低于饱和压力时，油藏中将同时存在单相和两相流动，拟稳态条件下产量的一般计算表达式为:,需要分段积分,三、prpbpwf时的流入动态（1）基本公式需要分段积分,35,图1-11 组合型IPR曲线,（2）实用计算方法,图1-11 组合型IPR曲线（2）实用计算方法,36,当,p,r,p,b,时，由于油藏中全部为,单相,液体流动,流入动态公式为：,流压等于饱和压力时的产量为：,当,p,r,pb时，由于油藏中全部为单相液体流动流压等于饱和,37,A-油相IPR曲线,B-水相IPR曲线,C-油气水三相综合IPR曲线,四、油气水三相IPR 曲线,Petrobras,提出了计算三相流动IPR曲线的方法,综合IPR曲线的实质,:,是按,含水率,取纯油IPR曲线和水IPR曲线的,加权平均值,。当已知测试点计算曲线，可按产量加权平均；当预测产量或流压，可按流压加权平均。,图1-12 油气水三相IPR 曲线,A-油相IPR曲线四、油气水三相IPR 曲线Petrobr,38,(一)采液指数计算,(,由测试点确定曲线,),已知,p,r,、,p,b,和一个测试点,p,wf,(test),、,q,t,(test),(1),(2),(一)采液指数计算(由测试点确定曲线)已知 pr、pb和一,39,(二)某一产量,q,t,下的流压,p,wf,计算,(1),(二)某一产量qt下的流压pwf计算(1),40,(2),(2),41,因为：,所以：,综合IPR曲线的斜率可近似为常数,(3),图1-12 油气水三相IPR 曲线,因为：所以：综合IPR曲线的斜率可近似为常数(3)图1-1,42,五、,多层油藏油井流入动态,（1）多油层油井流入动态,图1-13 多层油藏油井流入动态,流压低于,14MPa,后，只有,第三个层,工作；,流压降低到,1