《油井流入动态幻灯片.ppt》由会员分享,可在线阅读,更多相关《油井流入动态幻灯片.ppt(56页珍藏版)》请在得力文库 - 分享文档赚钱的网站上搜索。
1、油井流入动态第1页,共56页,编辑于2022年,星期日油井生产系统组成油层到井底的流动(地层渗流)井底到井口的流动(井筒多相管流)井口到分离器(地面水平或倾斜管流)油井生产的三个基本流动过程气液两相流基本理论油井流入动态第2页,共56页,编辑于2022年,星期日第一节 油井流入动态(IPR曲线)油井流入动态:油井产量(qo)与井底流动压力(pwf)的关系,反映了油藏向该井供油的能力。基本概念油井流入动态曲线:表示产量与流压关系的曲线,简称IPR曲线。Inflow Performance Relationship Curve第3页,共56页,编辑于2022年,星期日4 反映了油藏向井的供油能力,
2、反映了油藏压力、油层物性、流体反映了油藏向井的供油能力,反映了油藏压力、油层物性、流体物性、完井质量等对油层渗流规律的影响,是采油工程与油藏工程的物性、完井质量等对油层渗流规律的影响,是采油工程与油藏工程的衔接点。衔接点。作用:作用:为油藏工程提供检验资料;为油藏工程提供检验资料;为采油工程的下一步工作提供依据;为采油工程的下一步工作提供依据;检查钻井、完井和各项工艺措施等技术水平的优劣。检查钻井、完井和各项工艺措施等技术水平的优劣。物理意义:物理意义:第一节第一节 油井流入动态油井流入动态油井流入动态:油井流入动态:油井产量与井底流动压力的关系。油井产量与井底流动压力的关系。第4页,共56页
3、,编辑于2022年,星期日图1-1 典型的流入动态曲线IPR曲线基本形状与油藏驱动类型驱动类型有关。即使在同一驱动方式下,还将取决于油藏压力、油层厚度、渗透率及流体物理性质等。prqomax第5页,共56页,编辑于2022年,星期日一、单相液体流入动态供给边缘压力不变、圆形地层中心一口井的产量公式为:(1-1)圆形封闭油藏、拟稳态条件下产量公式为:(1-2)第6页,共56页,编辑于2022年,星期日图1-2 泄油面积形状与油井的位置系数对于非圆形封闭泄油面积的油井产量公式,可根据泄油面积和油井位置进行校正。第7页,共56页,编辑于2022年,星期日单相流动时,油层物性及流体性质基本不随压力单相
4、流动时,油层物性及流体性质基本不随压力变化。变化。采油指数采油指数可定义为:单位生产压差下的油井产油量,是反映油层性质、厚度、流体参数、完井条件及泄油面积等与产量之间的关系的综合指标。生产压差生产压差直线型直线型第8页,共56页,编辑于2022年,星期日采油指数J的获得:试井资料:试井资料:测得测得3 35 5个稳定工作制度下的产量及其流压,便可个稳定工作制度下的产量及其流压,便可绘制该井的实测绘制该井的实测IPRIPR曲线,取曲线,取其斜率的负倒数其斜率的负倒数油藏参数计算油藏参数计算对对于于单单相相液液体体流流动动的的直直线线型型IPRIPR曲曲线线,采采油油指指数数可可定定义义为为产产油
5、油量量与与生生产产压压差差之之比比,也也可可定定义义为为每每增增加加单单位位生生产产压压差差时时,油油井井产产量量的的增增加加值值,或或IPRIPR曲线斜率的负倒数曲线斜率的负倒数。注意事项:因此,对于具有非直线型因此,对于具有非直线型IPRIPR曲线的油井,在使用采油指数时,应该曲线的油井,在使用采油指数时,应该说明相应的流动压力,不能简单地用某一流压下的采油指数来直接推说明相应的流动压力,不能简单地用某一流压下的采油指数来直接推算不同流压下的产量。算不同流压下的产量。第9页,共56页,编辑于2022年,星期日当油井产量很高时,井底附近将出现非达西渗流非达西渗流:胶结地层的紊流速度系数:非胶
6、结地层紊流速度系数:C、D值也可用试井资料获取值也可用试井资料获取第10页,共56页,编辑于2022年,星期日二、油气两相渗流时的流入动态o、Bo、Kro都是压力的函数。用上述方法绘制IPR曲线十分繁琐。通常结合生产资料来绘制通常结合生产资料来绘制 IPR曲线。曲线。平面径向流,直井油气两相渗流时油井产量公式为:(一)垂直井油气两相渗流时的流入动态第11页,共56页,编辑于2022年,星期日1.1.Vogel Vogel 方法方法(1968)(1968)假设条件:a.圆形封闭油藏,油井位于中心;b.均质油层,含水饱和度恒定;c.忽略重力影响;d.忽略岩石和水的压缩性;e.油、气组成及平衡不变;
7、f.油、气两相的压力相同;g.拟稳态下流动,在给定的某一瞬间,各点的脱气原油流量相同。数值模拟结果的数值模拟结果的总结总结第12页,共56页,编辑于2022年,星期日归一化归一化曲线曲线第13页,共56页,编辑于2022年,星期日Vogel方程经典方程经典方程第14页,共56页,编辑于2022年,星期日 a.计算c.根据给定的流压及计算的相应产量绘制IPR曲线b.给定不同流压,计算相应的产量:、已知地层压力和一个工作点、已知地层压力和一个工作点(qo(test),pwf(test)利用Vogel方程绘制IPR曲线的步骤第15页,共56页,编辑于2022年,星期日、已知两个工作点,油藏压力未知、
8、已知两个工作点,油藏压力未知a.油藏平均压力的确定:已知或利用两组qopwf 测试计算,即 b.计算c.由流入动态关系式计算相关参数第16页,共56页,编辑于2022年,星期日图2-4 计算的溶解气驱油藏油井IPR曲线1-用测试点按直线外推;2-计算机计算值;3-用Vogel方程计算值Vogel曲线与数值模拟IPR曲线的对比第17页,共56页,编辑于2022年,星期日对比结果:对比结果:按按Vogel方方程程计计算算的的IPR曲曲线线,最最大大误误差差出出现现在在用用小小生生产产压压差差下下的的测测试试资资料料来来预预测测最最大大产产量量时时,但但一一般般误误差差低低于于5。虽虽然然随随着着采
9、采出出程程度度的的增增加加到到开开采采末末期期误误差差上上升升到到20左左右右,但其绝对值却很小但其绝对值却很小。如如果果用用测测试试点点的的资资料料按按直直线线外外推推,最最大大误误差差可可达达 7080,只是在开采末期约,只是在开采末期约30%。采采出出程程度度 Np 对对油油井井流流入入动动态态影影响响大大,而而kh/、Bo、k、So等对其影响不大。等对其影响不大。第18页,共56页,编辑于2022年,星期日2.费特柯维奇方法溶解气驱油藏溶解气驱油藏假设假设(k kroro/o oB Bo o)与压力与压力p p 成线性关系,则成线性关系,则其中,其中,第19页,共56页,编辑于2022
10、年,星期日所以:当 时:令:令:费特柯维奇费特柯维奇基本方程基本方程第20页,共56页,编辑于2022年,星期日3.3.不完善井不完善井VogelVogel方程的修正方程的修正油水井的不完善性:射孔完成射孔完成打开性质不完善打开性质不完善;未全部钻穿油层未全部钻穿油层打开程度不完善打开程度不完善;打开程度和打开性质打开程度和打开性质双重不完善双重不完善;在钻井或修井过程中油层受到损害或进行酸化、压裂等在钻井或修井过程中油层受到损害或进行酸化、压裂等措施,从而改变油井的完善性。措施,从而改变油井的完善性。第21页,共56页,编辑于2022年,星期日 图1-5 完善井和非完善井周围 的压力分布示意
11、图第22页,共56页,编辑于2022年,星期日油井的流动效率FE:油井的理想生产压差与实际生产压差之比油井的理想生产压差与实际生产压差之比为为“正正”称称“正正”表皮,油井不完善;表皮,油井不完善;为为“负负”称称“负负”表皮,油井超完善。表皮,油井超完善。第23页,共56页,编辑于2022年,星期日完善井完善井非完善井非完善井令:非完善井表皮附加压力降于是于是第24页,共56页,编辑于2022年,星期日表皮系数或井壁阻力系数表皮系数或井壁阻力系数S S油层受污染的或不完善井,油层受污染的或不完善井,完善井完善井,增产措施后的超完善井,增产措施后的超完善井,表皮系数表皮系数S S通常由试井方法
12、获得通常由试井方法获得第25页,共56页,编辑于2022年,星期日利用流动效率计算直井流入动态的方法Standing方法(1970)(FE=0.5 1.5)图1-6 FE 1时的无因次IPR曲线(standing IPR曲线)第26页,共56页,编辑于2022年,星期日a.根据已知根据已知pr和和pwf计算在计算在FE=1时最大产量时最大产量standing方法计算不完善井IPR曲线的步骤:b.b.预测不同流压下的产量预测不同流压下的产量c.c.根据计算结果绘制根据计算结果绘制IPRIPR曲线曲线第27页,共56页,编辑于2022年,星期日Harrison方法(FE=1 2.5)图1-7 Ha
13、rrison无因次IPR曲线(FE1)第28页,共56页,编辑于2022年,星期日 HarrisonHarrison方法可用来计算高流动效率井的方法可用来计算高流动效率井的IPRIPR曲线和曲线和预测低流压下的产量。其计算步骤如下:预测低流压下的产量。其计算步骤如下:a.计算FE=1时的qomax 先求先求pwf/pr,然后查图然后查图1-7中对应的中对应的FE曲线上的相应值曲线上的相应值 qo/qomax(FE=1),则则b.计算不同流压下的产量c.根据计算结果绘制IPR曲线d.求FE对应的最大产量,即pwf=0时的产量第29页,共56页,编辑于2022年,星期日(二)斜井和水平井的IPR曲
14、线1990年,Cheng对溶解气驱油藏中斜井和水平井进行了数值模拟,并用回归的方法得到了类似Vogel方程的不同井斜角井的IPR回归方程:p=pwf/pr;q=qo/qomax ;A、B、C为取决于井斜角的系数优点:优点:使用简单,仅需一组测点,便可得使用简单,仅需一组测点,便可得IPR曲线曲线缺点:缺点:方程没有归一化,方程没有归一化,第30页,共56页,编辑于2022年,星期日19891989年,年,BendakhliaBendakhlia等用两种三维三相黑油模拟器研究了多等用两种三维三相黑油模拟器研究了多种情况下溶解气驱油藏中水平井的流入动态关系。得到了不种情况下溶解气驱油藏中水平井的流
15、入动态关系。得到了不同条件下同条件下IPRIPR曲线。曲线。图图1-8 1-8 拟合的拟合的IPRIPR曲线与实际曲线的曲线与实际曲线的对比 _拟合的IPR曲线,实际曲线第31页,共56页,编辑于2022年,星期日曲线表明:早期的IPR曲线近似于直线,随着采出程度增加,曲度增加,接近衰竭时曲度稍有减小。Bendakhlia建议用以下公式来拟合IPR曲线图:第32页,共56页,编辑于2022年,星期日 图图1-9 参数参数v、n与采出程度之间的关系第33页,共56页,编辑于2022年,星期日IPR曲线的应用曲线的应用油井流入动态反映了油藏向该井供油的能力。1.根据测试资料确定根据测试资料确定IP
16、R曲线。曲线。2.根据根据IPR曲线确定曲线确定流压和产量的对应流压和产量的对应关系。关系。prqomax第34页,共56页,编辑于2022年,星期日三、prpbpwf时的流入动态(1)基本公式 当油藏压力高于饱和压力,而流动压力低于饱和压力时,油藏中将同时存在单相和两相流动,拟稳态条件下产量的一般计算表达式为:需要分段需要分段积分积分第35页,共56页,编辑于2022年,星期日图1-11 组合型IPR曲线(2)实用计算方法第36页,共56页,编辑于2022年,星期日 当prpb时,由于油藏中全部为单相液体流动流入动态公式为:流压等于饱和压力时的产量为:当prpb后,油藏中出现两相流动流入动态
17、公式为:采油指数:第37页,共56页,编辑于2022年,星期日A-油相IPR曲线B-水相IPR曲线 C-油气水三相综合IPR曲线四、油气水三相IPR 曲线Petrobras提出了计算三相流动IPR曲线的方法综合IPR曲线的实质:是按含水率取纯油IPR曲线和水IPR曲线的加权平均值。当已知测试点计算曲线,可按产量加权平均;当预测产量或流压,可按流压加权平均。图1-12 油气水三相IPR 曲线第38页,共56页,编辑于2022年,星期日(一)采液指数计算(由测试点确定曲线)已知 pr、pb和一个测试点pwf(test)、qt(test)(1)(2)图1-12 油气水三相IPR 曲线第39页,共56
18、页,编辑于2022年,星期日 (二)某一产量qt下的流压pwf计算(1)第40页,共56页,编辑于2022年,星期日(2)图1-12 油气水三相IPR 曲线第41页,共56页,编辑于2022年,星期日因为:所以:综合IPR曲线的斜率可近似为常数(3)图1-12 油气水三相IPR 曲线第42页,共56页,编辑于2022年,星期日五、多层油藏油井流入动态(1)多油层油井流入动态图1-13 多层油藏油井流入动态流压低于流压低于1414MPaMPa后,只后,只有有第三个层第三个层工作;工作;流压降低到流压降低到1212MPaMPa和和1010MPaMPa后,则后,则I I层和层和IIII层层陆续出油。
19、陆续出油。总的总的IPRIPR曲线曲线则是则是分层的迭加分层的迭加。其特。其特点是:随着流压的降低,点是:随着流压的降低,由于参加工作的小层数由于参加工作的小层数增多,产量将大幅度增增多,产量将大幅度增加,采油指数也随之增加,采油指数也随之增大。大。第43页,共56页,编辑于2022年,星期日(2)含水油井流入动态图1-14 含水油井流入动态与含水变化 ()图1-15 含水油井流入动态曲线 ()第44页,共56页,编辑于2022年,星期日小 结(1)上述介绍的方法阐明了油井流入动态的物理意义,也是目前现场最常用的计算方法。(2)油井流入动态研究主要有三种途径:基于Vogel、Fetkovich
20、、Petrobras方法的完善。建立不同类型油藏和井底条件的渗流模型。利用单井流入动态的油藏数值模拟技术。(3)油井流入动态是采油工程各项技术措施设计、分析与评价的依据。第45页,共56页,编辑于2022年,星期日第二节 井筒气液两相流基本概念相的概念相的概念 相是体系中具有相同化学组成和物理性质的一部 分,与体系的其它均匀部分有界面隔开例如:例如:水水-冰系统、泥浆、油冰系统、泥浆、油-气气-水等均是多相体系水等均是多相体系油气是深埋于地下的油气是深埋于地下的流体流体矿藏矿藏采油设备的优化设计和工况分析、油气集输设计等采油设备的优化设计和工况分析、油气集输设计等都离不开都离不开气液气液两相流
21、两相流的理论与计算方法的理论与计算方法随压力的降低,溶解气将不断从原油中逸出,因此,井筒中将随压力的降低,溶解气将不断从原油中逸出,因此,井筒中将不可避免地出现不可避免地出现气液两相流动气液两相流动。第46页,共56页,编辑于2022年,星期日一、井筒气液两相流动的特性(一)气液两相流动与单相液流的比较第47页,共56页,编辑于2022年,星期日流动型态(流动结构、流型):流动过程中油、气的分布状态。流动过程中油、气的分布状态。(二)气液混合物在垂直管中的流动结构变化 纯液流 当井筒压力大于饱和压力时,天然气溶解在原油中,产液呈单相液流。影响流型的因素:气液体积比、流速、气液界面性质等。气液体
22、积比、流速、气液界面性质等。第48页,共56页,编辑于2022年,星期日泡流 井筒压力稍低于饱和压力时,溶解气开始从油中井筒压力稍低于饱和压力时,溶解气开始从油中分离出来,气体都以小气泡分散在液相中。分离出来,气体都以小气泡分散在液相中。滑脱现象:混合流体流动过程中,由于流体间的密度差异,引起的小密度流体流速大于大密度流体流速的现象。如:油气滑脱、气液滑脱、油水滑脱等。特点:气体是分散相,液体是连续相;气体主要影响混合物密度,对摩擦阻力影响不大;滑脱现象比较严重。第49页,共56页,编辑于2022年,星期日段塞流 当混合物继续向上流动,压力逐渐降低,气当混合物继续向上流动,压力逐渐降低,气体不
23、断膨胀,小气泡将合并成大气泡,直到能够体不断膨胀,小气泡将合并成大气泡,直到能够占据整个油管断面时,井筒内将形成一段液一段占据整个油管断面时,井筒内将形成一段液一段气的结构。气的结构。特点:气体呈分散相,液体呈连续相;气体呈分散相,液体呈连续相;一段气一段液交替出现;一段气一段液交替出现;气体膨胀能得到较好的利用;气体膨胀能得到较好的利用;滑脱损失变小;滑脱损失变小;摩擦损失变大。摩擦损失变大。第50页,共56页,编辑于2022年,星期日环流 油管中心是连续的气流而管壁为油管中心是连续的气流而管壁为油环的流动结构。油环的流动结构。特点:气液两相都是连续相;气液两相都是连续相;气体举油作用主要是
24、靠摩擦携带;气体举油作用主要是靠摩擦携带;摩擦损失变大。摩擦损失变大。第51页,共56页,编辑于2022年,星期日雾流 气体的体积流量增加到足够大时,油管气体的体积流量增加到足够大时,油管中内流动的气流芯子将变得很粗,沿管壁中内流动的气流芯子将变得很粗,沿管壁流动的油环变得很薄,绝大部分油以小油流动的油环变得很薄,绝大部分油以小油滴分散在气流中。滴分散在气流中。特点:气体是连续相,液体是分散相;气体是连续相,液体是分散相;气体以很高的速度携带液滴喷出井口;气体以很高的速度携带液滴喷出井口;气、液之间的相对运动速度很小;气、液之间的相对运动速度很小;气相是整个流动的控制因素。气相是整个流动的控制
25、因素。第52页,共56页,编辑于2022年,星期日总结:油井生产中可能出现的流型自油井生产中可能出现的流型自下而上依次为:纯油下而上依次为:纯油(液液)流、泡流、流、泡流、段塞流、环流和雾流。段塞流、环流和雾流。实际上,在同一口井内,一般不实际上,在同一口井内,一般不会出现完整的流型变化。会出现完整的流型变化。图1-17 油气沿井筒喷出时的流型变化示意图纯油流;泡流;段塞流;环流;雾流第53页,共56页,编辑于2022年,星期日(三)滑脱损失概念因滑脱而产生的附加压力损失称为滑脱损失。图1-18 气液两相流流动断面简图第54页,共56页,编辑于2022年,星期日无滑脱无滑脱实际实际由于有滑脱时,气体流速大,液体流速小,为了保持体积流量不变,气气体体过过流流断断面面将将减小,而液体的过流断面将增加减小,而液体的过流断面将增加。由于滑脱存在:第55页,共56页,编辑于2022年,星期日单位管长上滑脱损失为:单位管长上滑脱损失为:滑脱损失的实质:液相的流动断面增大引起混合物密度的增加。图1-18 气液两相流流动断面简图第56页,共56页,编辑于2022年,星期日
限制150内