资源描述
单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,采油工程原理与设计,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,采油工程原理与设计,(讲课稿),讲课教师：刘洪,绪 论,采油工程是什么?,采油工程专业方向的水平如何?,采油工程师的要求?,?采油工程?课程介绍什么？,据评估，我国石油天然气可开采资源量为,110,160,10,8,t,，按,150,10,8,t,估算，人均为,12t/,人，而世界人均量为,68t/,人。,我国探明石油可采储量为,10,8,t,，人均为人，相当于世界人均量的,12%,。,中国油气资源状况：少、分散,油区：20多个,陆上如大庆、胜利、,辽河、新疆等,海上如渤海、东海、,南海东部、南海西部等,2021年产量：,石油：亿吨,天然气：830亿m3,与2021年相比增长7.7%。,进口石油：2.04104t,原油对外依存度首次超过50%。,效益：,油价：美元/桶,本钱：25美元/桶左右,中国油气生产状况：,效益良好、供不应求、形势严峻,预计2021年全国石油年产量为1.8108t；,预计2021年石油消耗量力求控制在4.3108t。,世界,50,大石油公司最新排名中，中国石油、中国石化、中国海油分列第,5,、第,25,和第,48,位。,海外油气产量首超亿吨，权益产量达,5500,万吨,海上石油产量稳步增长，占全国石油总产量的比重从,2000,年的,6.7%,升至,11%,。,长庆油田生产油气量首次突破,3000,万吨，达,3006,万吨，跃居成为我国第二大油气田，成为近,10,年我国油气产量增长最快的油田。,1.,采油工程的定义,采油工程：,油田开采过程中根据开发目标通过生产井和注入井所采取的工程技术措施的总称。,油气高效入井,经济有效举升,油气田开发效益,1经济有效举升,实质：,实现流体向高处流动,保证生产系统运行效率和油气资源开发经济效益。,实现入井油气等流体向地面的流动。,2油气高效入井,技术措施：,通过各种物理、化学作用处理油气层。,酸处理技术,压裂技术,波场处理技术,微生物技术,热力技术,采油工程系统组成,油藏：具有一定储存和流动特性的孔隙或裂缝介质系统，,人工建造系统：井底、井筒、井口装置、采油设备、注水设备以及地面集输、别离和储存设备等。,经济有效地提高油井,产量和原油采收率,采油工,程目标,采油工程特点：,涉及的技术面广、综合性强而又复杂,与油藏工程、地面工程和钻井工程等紧密联系,工作对象是条件随油藏动态不断变化的采、注井,难度大,涉及油田开发的重要决策和经济效益,2.,采油工程技术水平情况,采油工程技术是决定油田产量上下、采油速度快慢、最终采收率大小、经济效益的优劣等重要问题的关键技术。,中国采油工程技术大致经历了5个阶段：,15060年代初的探索、试验阶段,26070年代分层开采工艺配套技术开展阶段,37080年代多种油藏类型采油工艺技术开展阶段,48090年代采油工程新技术重点突破开展阶段,590年代中期以来采油系统工程形成和开展阶段,采油工程技术面临难题：,高含水油田开采,低渗透油田开发,西部沙漠超深油田开采,海洋和滩海油田开发建设,稠油和超稠油油田开采,要求：,采油工程技术到达更高的水平，以创造更高的效益。,3.,采油工程的任务,当前和未来采油工程的两大任务：,如何采用先进适用的技术提高已开发油气田采收率、,提高油气井产量以及降低油气生产本钱？,如何依靠技术进步和技术创新将未投入开发的油田,投入开发？,学习目标：,掌握采油工程师的工作技能，接受严格的采油工程师根本训练，具有分析和解决采油工程实际问题的能力，以胜任工程设计、生产管理、科技开发和应用研究等方面的工作。,只有这样，才能适应和效劳于不断开展的石油工业，才能运用自己的知识去迎接石油工业新技术的挑战。,采油工艺,4.,课程内容简介,采油方法,增产维,护措施,自喷采油,气举采油,有杆泵采油,无杆泵采油,注水,酸化、酸压,砂、蜡、水、稠、低,工作原理,计算方法,工程设计,分 析,水力压裂,采油工程课程内容体系,第一章 油井流入动态与井筒多相流动计算,第二章 自喷与气举采油,第三章 有杆泵采油,第四章 无杆泵采油,第五章 注水,第六章 水力压裂技术,第七章 酸处理技术,第八章 复杂条件下的开采技术,第九章 完井方案设计与试油,第十章 采油工程方案设计概要,第一章 油井流入动态与井筒多相流计算,油井流入动态,井筒气液两相流根本概念,计算气液两相垂直管流方法,油气从油藏流入井底和在井筒中的流动是油气开采的两个根本流动过程。油井流入动态和井筒多相流动规律是油井各种举升方式设计和生产动态分析所需要的共同理论根底。,第一节 油井流入动态IPR曲线,油井流入动态曲线IPR曲线：,表示产量与井底流压关系的曲线，简称IPR曲线。,油井流入动态,：,油井产量与井底流动压力的关系。,它反映了油藏向井的供油能力，反映了油藏压力、油层物性、流体物性、完井质量等对油层渗流规律的影响，是采油工程与油藏工程的衔接点。通过油井流入动态研究为油藏工程提供检验资料；为采油工程的下一步工作提供依,据；检查钻井、固井、完井和各项工艺措施等技术水平的优劣。,图1-1 典型的流入动态曲线,一、,单相液体流入动态,(基于达西定律),供给边缘压力不变圆形地层中心一口井的产量公式为：,圆形封闭油藏,，拟稳态条件下的油井产量公式为：,泄油面积形状与油井的位置系数,对于非圆形封闭泄油面积的油井产量公式，可根据泄油面积和油井位置进行校正。,采油(液)指数,：,单位生产压差下油井产油(液)量，反映油层性质、厚度、流体物性、完井条件及泄油面积等与产量之间关系的综合指标。,单相流动时，油层物性及流体性质根本不随压力变化，产量公式可表示为:,采油指数,J,的获得：,油藏参数计算,试井资料：,测得35个稳定工作制度下的产量及其流压，便可绘制该井的实测,IPR,曲线,对于单相液体流动的直线型,IPR,曲线，采油指数可定义为产油量与生产压差之比，或者单位生产压差下的油井产油量；也可定义为每增加单位生产压差时，油井产量的增加值，或,油井,IPR,曲线斜率的负倒数。,本卷须知：,图1-1 典型的流入动态曲线,对于非直线型IPR曲线，由于其斜率不是定值，按上述几种定义所求得的采油指数那么不同。所以，对于具有非直线型IPR曲线的油井，在使用采油指数时，应该说明相应的流动压力，不能简单地用某一流压下的采油指数来直接推算不同流压下的产量。,当油井产量很高时，井底附近将出现,非达西渗流,：,胶结地层的,紊流速度系数,：,非胶结地层,紊流速度系数,：,如果在单相流动条件出现非达西渗滤，也可利用试井所得的产量和压力资料求得,C,和,D,值。,由试井资料绘制的 直线的斜率为,D,，其截距则为,C,。,二、,油气两相渗流时的流入动态,o、Bo、Kro都是压力的函数。用上述方法绘制IPR曲线十分繁琐。通常结合生产资料来绘制IPR曲线Vogel方法或简化的Fetkovich方法。,平面径向流，直井,油气两相渗流,时油井产量公式为：,(一)垂直井油气两相渗流时的流入动态,1.,Vogel,方法(1968),假设条件,：,a.,圆形封闭油藏，油井位于中心；溶解气驱油藏。,b.,均质油层，含水饱和度恒定；,c.,忽略重力影响；,d.,忽略岩石和水的压缩性；,e.,油、气组成及平衡不变；,f.,油、气两相的压力相同；,g.,拟稳态下流动，在给定的某一瞬间，各点的脱气原油流量相同。,Vogel,方程,Vogel,曲线,a.,计算,c.,根据给定的流压及计算的相应产量绘制,IPR,曲线。,b.,给定不同流压，计算相应的产量：,地层压力和一个工作点：,利用,Vogel,方程绘制,IPR,曲线的步骤,油藏压力未知，两个工作点,a.油藏平均压力确实定,b.,计算,d.,根据给定的流压及计算的相应产量绘制,IPR,曲线,c.,给定不同流压，计算相应的产量,Vogel曲线与数值模拟IPR曲线的比照,a.,按,Vogel,方程计算的,IPR,曲线，最大误差出现在用小生产压差下的测试资料来预测最大产量。一般，误差低于5。虽然，随着采出程度的增加，到开采末期误差上升到20，但其绝对值却很小。,b.,如果用测试点的资料按直线外推时，最大误差可达 7080，只是在开采末期约30%。,C.,采出程度N对油井流入动态影响大，而kh/、B,0,、k、S,0,等参数对其影响不大。,图1-4 不同方法计算的油井,IPR,曲线1-用测试点按直线外推；2-计算机计算的；3-用,Vogel,方程计算的,溶解气驱油藏：,假设 与压力,成直线关系，则：,式中：,则：,令：,当 时：,所以：,Vogel,方程的修正,油水井的非完善性：,翻开性质不完善；如射孔完成,翻开程度不完善；如未全部钻穿油层,翻开程度和翻开性质双重不完善,油层受到损害,酸化、压裂等措施,改变油井的完善性，从而增加或降低井底附近的压力降，影响油井流入动态关系。,完善井和非完善井周围的压力分布示意图,完善井:,非完善井:,令：,非完善井附加压力降:,则：,油井的流动效率FE：,油井的理想生产压差与实际生产压差之比。,油层受污染的或不完善井，,完善井,增产措施后的超完善井，,利用流动效率计算非完善直井流入动态的方法,图1-6,Standing,无因次,IPR,曲线,Standing,方法(,1.5),Standing,方法计算不完善井,IPR,曲线的步骤：,a.,根据已知,P,r,和,P,wf,计算在,FE=1,时最大产量,b.,预测不同流压下的产量,c.,根据计算结果绘制,IPR,曲线,Harrison,方法,图1-7,Harrison,无因次,IPR,曲线(,FE1),b.,求,FE,对应的最大产量，即,P,wf,=0,时的产量,Harrison,方法计算超完善井,IPR,曲线的步骤：,a.,计算,FE=1,时的,q,omax(FE=1),先求,P,wf,/P,r,，,然后查图1-7中对应的,FE,曲线上的相应值,q,o,/q,omax(FE=1)。,则,c.,计算不同流压下的产量,d.,根据计算结果绘制,IPR,曲线,(二)斜井和水平井的,IPR,曲线,Cheng,对溶解气驱油藏中斜井和水平井进行了数值模拟，并用回归的方法得到了类似,Vogel,方程的不同井斜角井的,IPR,回归方程：,P,=P,wf,/P,r,； q,=q,o,/q,omax,；A、B、C,为取决于井斜角的系数。,Bendakhlia,等用两种三维三相黑油模拟器研究了多种情况下溶解气驱油藏中水平井的流入动态关系。得到了不同条件下,IPR,曲线。,Bendakhlia,用公式来拟合,IPR,曲线图版，发现吻合很好。,图1-8 拟合的IPR曲线与实际曲线的比照 _拟合的IPR曲线,实际曲线,图1-9 参数v、n与采收率系数之间的关系,其它水平井产能计算模型：,Borisov,模型：,Giger,模型：,Joshi,模型：,Renard & Dupuy,模型：,1根本公式,当油藏压力高于饱和压力，而流动压力低于饱和压力时，油藏中将同时存在单相和两相流动，拟稳态条件下产量的一般表达式为:,三、,时的流入动态,井周围压力和气体饱和分布,组合型,IPR,曲线,2实用计算方法(组合型IPR方法),流压等于饱和压力时的产量为：, 当 时，由于油藏中全部为,单相,液体流动。, 当 后，油藏中出现,两相,流动。,流入动态公式为：,直线段,采油指数,四、油气水三相,IPR,曲线,Petrobras,提出了计算三相流动,IPR,曲线的方法。,综合IPR曲线的实质:,按含水率取纯油IPR曲线和水IPR曲线的加权平均值。当测试点计算采液指数时，是按产量加权平均；当预测产量或流压时是按流压加权平均。,图1-12 油气水三相,IPR,曲线,(一)采液指数计算,已知一个测试点( 、 )和饱和压力 及油藏压力, 当 时：,思考题：推导 时的采液指数计算式。, 当 时：, 当 时：,其中：,直线段,采液指数, ，则：, ，则按流压加权平均进行推导：,(二)某一产量下的流压计算,所以：,因为：,所以：,若 ，则综合,IPR,曲线的斜率可近似为常数。,思考题：试推导,五、,多层油藏油井流入动态,1多油层油井流入动态迭加型IPR,图1-13 多层油藏油井流入动态,2含水油井流入动态,图1-14 含水油井流入动态与含水变化,( ),图1-15 含水油井流入动态曲线,( ),小 结,(1) 介绍的方法说明了油井流入动态的物理意义，也是目前现场最常用的计算方法。,(2) 油井流入动态研究主要有三种途径：,基于,Vogel,、Fetkovich、Petrobras,方法的完善。,建立不同类型油藏和井底条件的渗流模型。,利用单井流入动态的油藏数值模拟技术。,(3) 油井流入动态是采油工程各项技术措施设计、分析与评价的依据。,第二节 井筒气液两相流根本概念,井筒多相流理论：,研究各种举升方式油井生产规律根本理论,研究特点,：,流动复杂性、无严格数学解,研究途径：根本流动方程,实验资料相关因次分析,近似关系,一、,井筒气液两相流动的特性,(一)气液两相流动与单相液流的比较,流动型态流动结构、流型：,流动过程中油、气的分布状态。,(二)气液混合物在垂直管中的流动结构变化,影响流型的因素,各相介质的体积比例,介质的流速,各相的物理及化学性质(密度、粘度界面张力等),流道的几何形状,壁面特性,管道的安装方式,流动型态的划分方法：两类,第一类划分方法：根据两相介质分布的外形划分，包括泡状流、弹状流或团状流、层状流、波状流、段塞流或冲击流、环状流、雾状流,垂直气液两相流流型,水平气液两相流流型,第二类划分方法：按流动的数学模型或流体的分散程度划分，包括分散流、间歇流、别离流,两种分类方法比较,第一类划分方法较为直观,第二类划分方法便于进行数学处理,第一类划分方法,泡状流,弹状流或团状流,层状流,波状流,段塞流或冲击流,环状流,雾状流,第二类划分方法,分散流,间歇流,别离流,别离流,间歇流,别离流,分散流,两类划分结果,的对应关系, 纯液流,当井筒压力大于饱和压力时，天然气溶解在原油中，产液呈单相液流。,泡流,井筒压力稍低于饱和压力时，溶解气开始从油中别离出来，气体都以小气泡分散在液相中。,滑脱现象：,混合流体流动过程中，由于流体间的密度差异，引起的小密度流体流速大于大密度流体流速的现象。,如：油气滑脱、气液滑脱、油水滑脱等。,特点：气体是分散相，液体是连续相；,气体主要影响混合物密度，对摩擦阻力影响不大；,滑脱现象比较严重。,段塞流,当混合物继续向上流动，压力逐渐降低，气体不断膨胀，小气泡将合并成大气泡，直到能够占据整个油管断面时，井筒内将形成一段液一段气的结构。,特点：,气体呈分散相，液体呈连续相；,一段气一段液交替出现；,气体膨胀能得到较好的利用；,滑脱损失变小；,摩擦损失变大。,环流,油管中心是连续的气流而管壁为油环的流动结构。,特点：,气液两相都是连续相；,气体举油作用主要是靠摩擦携带；,滑脱损失变小；,摩擦损失变大。,雾流,气体的体积流量增加到足够大时，油管中内流动的气流芯子将变得很粗，沿管壁流动的油环变得很薄，绝大局部油以小油滴分散在气流中。,特点,：,气体是连续相，液体是分散相；,气体以很高的速度携带液滴喷出井口；,气、液之间的相对运动速度很小；,气相是整个流动的控制因素。,总结：,油井生产中可能出现的流型自下而上依次为：纯油(液)流、泡流、段塞流、环流和雾流。,实际上，在同一口井内，一般不会出现完整的流型变化。,图1-17 油气沿井筒喷出时的流型变化示意图,纯油流；,泡流；,段塞流；,环流；,雾流,实际计算：直接求存在滑脱混合物密度或包括滑脱在内的摩擦阻力系数。,(三)滑脱损失概念,因滑脱而产生的附加压力损失称为滑脱损失,。,单位管长上滑脱损失为：,图1-18 气液两相流流动断面简图,滑脱损失的实质,：,液相的流动断面增大引起混合物密度的增加。,二、井筒气液两相流能量平衡方程,及压力分布计算步骤,两个流动断面间的能量平衡关系：,(一)能量平衡方程推导,图2-19 倾斜管流能量平衡关系示意图,图2-19 倾斜管流能量平衡关系示意图,倾斜多相管流断面1和断面2的流体的能量平衡关系为：,适合于各种管流的通用压力梯度方程：,则：,令：,以计算段下端压力为起点，重复步，计算下一段的深,度和压力，直到各段的累加深度等于管长为止。,(2)多相垂直管流压力分布计算步骤,重复的计算，直至 。,1)按深度增量迭代的步骤,任一点(井口或井底)的压力作为起点，任选一个适宜,的压力降作为计算的压力间隔p。,估计一个对应的深度增量,h,。,计算该管段的平均温度及平均压力，并确定流体性质参数。,判断流型，并计算该段的压力梯度,dp/dh。,计算对应于的该段管长(深度差),h,。,计算该段下端对应的深度及压力。,2)按压力增量迭代的步骤略,思考题：,根据上述步骤整理出计算压力分布的程序流程框图。,说明：,a. 计算压力分布过程中，温度和压力是相关的；,b. 流体物性参数计算至关重要，但目前方法精度差；,c. 不同的多相流计算方法差异较大，因此在实际应用中有必要根据油井的实际情况筛选精度相对高的方法。,第三节,Orkiszewski,方法,综合了,Griffith & Wallis,和,Duns & Ros,等,方法,处理过渡性流型时，采用,Ros,方法(内插法),针对每种流动型态提出存容比及摩擦损失的计算方法,提出了四种流型，即泡流、段塞流、过渡流及环雾流,把Griffith段塞流相关式改进后推广到了高流速区,1967年提出，适用于垂直管流计算,出现雾流时，气体体积流量远大于液体体积流量。根据气体定律，动能变化可表示为：,一、,压力降公式及流动型态划分界限,由垂直管流能量方程可知，压力降是摩擦能量损失、势能变化和动能变化之和：,所以压降计算式为：,在雾流情况下才有明显意义,表1-3,Orkiszewski,方法流型划分界限,不同流动型态下 和 的计算方法不同。,二、平均密度及摩擦损失梯度的计算,气相存容比(含气率,),H,g,：,管段中气相体积与管段容积之比值。,液相存容比(持液率),H,L,：,管段中液相体积与管段容积之比值。,(1)泡流,平均密度:,？,滑脱速度,：,气相流速与液相流速之差,。,则：,泡流摩擦损失梯度,按液相进行计算：,Griffith,实验测得：,0.244m/s,图1-21 摩擦阻力系数曲线（教材,p37,）,图1-21,(2)段塞流,平均密度,:,段塞流的摩擦梯度：,段塞流计算关键是滑脱速度vs和液体分布系数的计算。,液体分布系数的计算见课本p39表1-4和公式1-58。,滑脱速度vs的计算方法有两种：查图迭代法和经验公式法。,泡流雷诺数：,图1-22,C,1,N,b,曲线,雷诺数：,图1-23,C,2,N,Re,曲线,滑脱速度的计算,查图迭代法,滑脱速度：,滑脱速度的计算,经验公式计算法,详见教材p3893公式1-541-57。,(3)过渡流,过渡流的混合物平均密度及摩擦梯度是先按段塞流和雾流分别进行计算，然后用内插方法来确定相应的数值。,雾流混合物,平均密度,计算公式与泡流相同：,由于雾流的气液无相对运动速度，即滑脱速度接近于零，根本上没有滑脱。,雾流摩擦系数可根据气体雷诺数和液膜相对粗糙度查图得。,摩擦梯度:,(4)雾流,所以：,图1-24,Orkiszewski,方法计算流程框图,第四节,Beggs & Brill,方法,水和空气、聚丙烯管实验根底上总结的方法,建立流型分布图，将七种流型归为三类，增加了过渡流,计算时先按水平管流计算，然后采用倾斜校正系数校,正成相应的倾斜管流,倾斜度 -90,+90，分上坡和下坡流动,1973年提出，适用于,水平、垂直和任意倾斜,管流计算,Beggs & Brill,两相水平管流型,分离流,分层流,波状流,环状流,间歇流,团状流,段塞流,分散流,泡 流,雾 流,一、根本方程,单位质量气液混合物稳定流动的机械能量守恒方程为：,(1)位差压力梯度,：消耗于混合物静水压头的压力梯度。,(2)摩擦压力梯度,：克服管壁流动阻力消耗的压力梯度。,假设条件,:气液混合物既未对外作功，也未受外界功。,(3)加速度压力梯度,：由于动能变化而消耗的压力梯度。,忽略液体压缩性、考虑到气体质量流速变化远远小于气体密度变化，那么：,(4)总压力梯度Beggs-Brill方法的根本方程,图1-26,Beggs-Brill,流型分布图（教材,p,45）,二、,Beggs & Brill,方法的流型分布图及流型判别式,别离流,间歇流,过渡流,分散流,表2-4,Beggs-Brill,法流型判别条件,三、持液率及混合物密度确定,(1)持液率,Beggs & Brill,方法计算倾斜管流时首先按水平管计算,然后进行倾斜校正。,表1-6,a、b、c,常数表,实验结果说明，倾斜校正系数与倾斜角、无滑脱持液率、弗洛德数及液体速度数有关。,图1-27 不同,E,L,下的倾斜校正系数,根据实验结果回归的倾斜校正系数的相关式为:,对于垂直管:,系数,C,与无滑脱持液率 、弗洛德数和液相速度数有关。,表1-6 系数,d、e、f、g,其中：,对于过渡流型，先分别用别离流和间歇流计算，之后采用内插法确定其持液率。,利用持液率计算流动条件下混合物实际密度：,四、阻力系数,气液两相流阻力系数与无滑脱气液两相流阻力系数的比值与持液率和无滑脱持液率(入口体积含液率)之间的关系：,当1,y1.2,时,其中：,两相流动的雷诺数：,也可用,Moody,图上的光滑管曲线来确定:,气液两相流阻力系数:,Beggs & Brill,方法计算流程框图(,p49),小 结,(1) 模拟计算多相管流规律的数学相关式及图版研究很多。纵观这许多数学相关式，其基本通式一般都是从基本能量守恒方程出发建立的：,(2) 对Poettmann-Carpenter方法、Fanch-Brown相关式、Baxendell-Thomas 相关式、Hagedron-Brown 关系式、Duns-Ros相关式、Orkiszewski相关式、Beggs-Brill 相关式、Dukler相关式、Mukherjee-Brill相关式、 Aziz相关式、Eaton相关式、Ansari相关式等十二种方法进行了比照分析，不同的方法有其适用条件和精度，可根据具体油田实际选用。,(3) 数学相关式大体分为三种类型：,在计算井筒流体混合密度时不考虑液体滞留量的影响，而液体滞留量与管壁摩阻损失用一个经验摩阻系数来表达，不区分多相流体的流态分布情况。,在计算多相流体混合物密度时考虑液体滞留量的影响，而摩阻系数要依据液体和气体的组成特征来确定,考虑液体滞留量的影响，摩阻系数取决于多相流体中连续相的特征。同时考虑了不同流态的影响，其中包括泡状流、段塞流、环流及雾状流。,END,
展开阅读全文