绪论油井流入动态与多相流.ppt

《绪论油井流入动态与多相流.ppt》由会员分享，可在线阅读，更多相关《绪论油井流入动态与多相流.ppt（90页珍藏版）》请在金锄头文库上搜索。

1、采油工程原理与设计采油工程原理与设计第一章 油井流入动态与井筒多相流计算 油井流入动态 井筒气液两相流基本概念 计算气液两相垂直管流方法 油气从油藏流入井底和在井筒中的流动是油气开采油气从油藏流入井底和在井筒中的流动是油气开采的两个基本流动过程。油井流入动态和井筒多相流动规的两个基本流动过程。油井流入动态和井筒多相流动规律是油井各种举升方式设计和生产动态分析所需要的共律是油井各种举升方式设计和生产动态分析所需要的共同理论基础。同理论基础。第一节 油井流入动态（IPR曲线） 油井流入动态曲线（IPR曲线）： 表示产量与井底流压关系的曲线，简称表示产量与井底流压关系的曲线，简称IPRIPR曲线。曲

2、线。油井流入动态： 油井产量与井底流动压力的关系。油井产量与井底流动压力的关系。它反映了油藏向井的它反映了油藏向井的供油能力，反映了油藏压力、油层物性、流体物性、完井质供油能力，反映了油藏压力、油层物性、流体物性、完井质量等对油层渗流规律的影响，是采油工程与油藏工程的衔接量等对油层渗流规律的影响，是采油工程与油藏工程的衔接点。通过油井流入动态研究为油藏工程提供检验资料；为采点。通过油井流入动态研究为油藏工程提供检验资料；为采油工程的下一步工作提供依油工程的下一步工作提供依据；检查钻井、固井、完井和各据；检查钻井、固井、完井和各项工艺措施等技术水平的优劣。项工艺措施等技术水平的优劣。图图1-1

3、1-1 典型的流入动态曲线典型的流入动态曲线一、单相液体流入动态(基于达西定律)供给边缘压力不变圆形地层供给边缘压力不变圆形地层中心一口井的产量公式为：中心一口井的产量公式为：圆形封闭油藏圆形封闭油藏，拟稳态条件下的油井产量公式为：，拟稳态条件下的油井产量公式为：采油工程原理与设计采油工程原理与设计泄油面积形状与油井的位置系数泄油面积形状与油井的位置系数 对于非圆形封对于非圆形封闭泄油面积的油井闭泄油面积的油井产量公式，可根据产量公式，可根据泄油面积和油井位泄油面积和油井位置进行校正。置进行校正。采油采油( (液液) )指数指数： 单位生产压差下油井产油单位生产压差下油井产油( (液液) )量

4、，反映油层性质、厚量，反映油层性质、厚度、流体物性、完井条件及泄油面积等与产量之间关系的度、流体物性、完井条件及泄油面积等与产量之间关系的综合指标。综合指标。 单相流动时，油层物性及流体性质基本不随压力变化，单相流动时，油层物性及流体性质基本不随压力变化，产量公式可表示为产量公式可表示为: :采油指数采油指数J的获得：的获得： 油藏参数计算油藏参数计算 试井资料：试井资料：测得测得3 35 5个稳定工作制度下的产量及其流个稳定工作制度下的产量及其流压，便可绘制该井的实测压，便可绘制该井的实测IPRIPR曲线曲线 对于单相液体流动的直线型对于单相液体流动的直线型IPRIPR曲线，采油指数可定义为

5、曲线，采油指数可定义为产油量与生产压差之比，或者单位生产压差下的油井产油量；产油量与生产压差之比，或者单位生产压差下的油井产油量；也可定义为每增加单位生产压差时，油井产量的增加值，或也可定义为每增加单位生产压差时，油井产量的增加值，或油油井井IPRIPR曲线斜率的负倒数。曲线斜率的负倒数。注意事项：注意事项：采油工程原理与设计采油工程原理与设计图图1-1 1-1 典型的流入动态曲线典型的流入动态曲线 对于非直线型对于非直线型IPRIPR曲线，由于其斜率不是定值，按上述几曲线，由于其斜率不是定值，按上述几种定义所求得的采油指数则不同。所以，对于具有非直线型种定义所求得的采油指数则不同。所以，对于

6、具有非直线型IPRIPR曲线的油井，在使用采油指数时，应该说明相应的流动压曲线的油井，在使用采油指数时，应该说明相应的流动压力，不能简单地用某一流压下的采油指数来直接推算不同流压力，不能简单地用某一流压下的采油指数来直接推算不同流压下的产量。下的产量。当油井产量很高时，井底附近将出现当油井产量很高时，井底附近将出现非达西渗流非达西渗流：胶结地层的胶结地层的紊流速度系数紊流速度系数：非胶结地层非胶结地层紊流速度系数紊流速度系数：采油工程原理与设计采油工程原理与设计 如果在单相流动条件出现非达西渗滤，也可利用试井如果在单相流动条件出现非达西渗滤，也可利用试井所得的产量和压力资料求得所得的产量和压力

7、资料求得C C和和D D值。值。 由试井资料绘制的由试井资料绘制的 直线的斜率为直线的斜率为D D，其截，其截距则为距则为C C。 二、油气两相渗流时的流入动态 o o、BoBo、KroKro都是压力的函数。用上述方法绘制都是压力的函数。用上述方法绘制IPRIPR曲线十分曲线十分繁琐。繁琐。通常结合生产资料来绘制通常结合生产资料来绘制IPRIPR曲线（曲线（VogelVogel方法）或简方法）或简化的化的FetkovichFetkovich方法。方法。平面径向流，直井平面径向流，直井油气两相渗流油气两相渗流时油井产量公式为：时油井产量公式为：(一)垂直井油气两相渗流时的流入动态1.1.Voge

8、l Vogel 方法方法(1968)(1968)假设条件假设条件：a.a.圆形封闭油藏，油井位于中心；溶解气驱油藏。圆形封闭油藏，油井位于中心；溶解气驱油藏。b.b.均质油层，含水饱和度恒定；均质油层，含水饱和度恒定；c.c.忽略重力影响；忽略重力影响；d.d.忽略岩石和水的压缩性；忽略岩石和水的压缩性； e.e.油、气组成及平衡不变；油、气组成及平衡不变；f.f.油、气两相的压力相同；油、气两相的压力相同；g.g.拟稳态下流动，在给定的某一瞬间，各点的脱气原油流量拟稳态下流动，在给定的某一瞬间，各点的脱气原油流量相同。相同。VogelVogel方程方程Vogel 曲线 a.a.计算计算c.c

9、.根据给定的流压及计算的相应产量绘制根据给定的流压及计算的相应产量绘制IPRIPR曲线。曲线。b.b.给定不同流压，计算相应的产量：给定不同流压，计算相应的产量： 已知地层压力和一个工作点：已知地层压力和一个工作点：利用利用VogelVogel方程绘制方程绘制IPRIPR曲线的步骤曲线的步骤 油藏压力未知，已知两个工作点油藏压力未知，已知两个工作点a.a.油藏平均压力的确定油藏平均压力的确定 b.b.计算计算d.d.根据给定的流压及计算的相应产量绘制根据给定的流压及计算的相应产量绘制IPRIPR曲线曲线c.c.给定不同流压，计算相应的产量给定不同流压，计算相应的产量Vogel曲线与数值模拟IP

10、R曲线的对比a.a.按按VogelVogel方程计算的方程计算的IPRIPR曲线，最大误差出现在用小生产压差下的测试资料曲线，最大误差出现在用小生产压差下的测试资料来预测最大产量。一般，误差低于来预测最大产量。一般，误差低于5 5。虽然，随着采出程度的增加，到开。虽然，随着采出程度的增加，到开采末期误差上升到采末期误差上升到2020，但其绝对值却很小。，但其绝对值却很小。b.b.如果用测试点的资料按直线外推时，最大误差可达如果用测试点的资料按直线外推时，最大误差可达 70 708080，只是在开，只是在开采末期约采末期约30%30%。C.C.采出程度采出程度N N对油井流入动态影响大，而对油井

11、流入动态影响大，而kh/kh/、B B0 0、k k、S S0 0等参数对其影响不等参数对其影响不大。大。图1-4 不同方法计算的油井IPR曲线1-用测试点按直线外推；2-计算机计算的；3-用Vogel方程计算的采油工程原理与设计采油工程原理与设计2.2.费特柯维奇方法费特柯维奇方法溶解气驱油藏：溶解气驱油藏：假设假设 与压力与压力 成直线关系，则：成直线关系，则：式中：式中：则：则：令：令：当当 时：时：所以：所以：3.3.非完善井非完善井VogelVogel方程的修正方程的修正油水井的非完善性：油水井的非完善性： 打开性质不完善；如射孔完成打开性质不完善；如射孔完成 打开程度不完善；如未全

12、部钻穿油层打开程度不完善；如未全部钻穿油层 打开程度和打开性质双重不完善打开程度和打开性质双重不完善 油层受到损害油层受到损害 酸化、压裂等措施酸化、压裂等措施 改变油井的完善性，从而增加或降低井底附改变油井的完善性，从而增加或降低井底附近的压力降，影响油井流入动态关系。近的压力降，影响油井流入动态关系。采油工程原理与设计采油工程原理与设计完善井和非完善井周围的压力分布示意图完善井和非完善井周围的压力分布示意图完善井完善井: :非完善井非完善井: :令：令：非完善井附加压力降非完善井附加压力降: :则：则：采油工程原理与设计采油工程原理与设计油井的流动效率（FE）：油井的理想生产压差与实际生产

13、压差之比。油井的理想生产压差与实际生产压差之比。油层受污染的或不完善井，油层受污染的或不完善井，完善井完善井, ,增产措施后的超完善井，增产措施后的超完善井，利用流动效率计算非完善直井流入动态的方法利用流动效率计算非完善直井流入动态的方法图图1-6 1-6 Standing Standing 无因次无因次IPRIPR曲线曲线 StandingStanding方法方法( (FE=0.5FE=0.5 1.5)1.5)StandingStanding方法计算不完善井方法计算不完善井IPRIPR曲线的步骤：曲线的步骤：a.a.根据已知根据已知P Pr r和和P Pwfwf计算在计算在FE=1FE=1时

14、最大产量时最大产量b.b.预测不同流压下的产量预测不同流压下的产量c.c.根据计算结果绘制根据计算结果绘制IPRIPR曲线曲线 HarrisonHarrison方法方法图图1-7 Harrison无因次无因次IPR曲线曲线(FE1)b.b.求求FEFE对应的最大产量，即对应的最大产量，即P Pwfwf=0=0时的产量时的产量HarrisonHarrison方法计算超完善井方法计算超完善井IPRIPR曲线的步骤：曲线的步骤：a.a.计算计算FE=1FE=1时的时的q qomax(FE=1)omax(FE=1) 先求先求P Pwfwf/P/Pr r，然后查图然后查图1-71-7中对应的中对应的FE

15、FE曲线上的相应值曲线上的相应值q qo o/q/qomax(FE=1)omax(FE=1)。则则c.c.计算不同流压下的产量计算不同流压下的产量d.d.根据计算结果绘制根据计算结果绘制IPRIPR曲线曲线采油工程原理与设计采油工程原理与设计非完善井的IPR曲线绘制采油工程原理与设计采油工程原理与设计采油工程原理与设计采油工程原理与设计 ( (二二) )斜井和水平井的斜井和水平井的IPRIPR曲线曲线 ChengCheng对溶解气驱油藏中斜井和水平井进行了对溶解气驱油藏中斜井和水平井进行了数值模拟，并用回归的方法得到了类似数值模拟，并用回归的方法得到了类似VogelVogel方程方程的不同井斜

16、角井的的不同井斜角井的IPRIPR回归方程：回归方程：P P=P=Pwfwf/P/Pr r； q q=q=qo o/q/qomaxomax ；A A、B B、C C为取决于为取决于井斜角的系数。井斜角的系数。BendakhliaBendakhlia等用两种三维三相黑油模拟器研究了等用两种三维三相黑油模拟器研究了多种情况下溶解气驱油藏中水平井的流入动态关多种情况下溶解气驱油藏中水平井的流入动态关系。得到了不同条件下系。得到了不同条件下IPRIPR曲线。曲线。BendakhliaBendakhlia用公式来拟合用公式来拟合IPRIPR曲线图版，发现吻合曲线图版，发现吻合很好。很好。图1-8 拟合的

17、IPR曲线与实际曲线的对比 _拟合的IPR曲线,实际曲线 图1-9 参数v、n与采收率系数之间的关系采油工程原理与设计采油工程原理与设计其它水平井产能计算模型：其它水平井产能计算模型：Borisov模型：模型：Giger模型：模型：Joshi模型：模型：Renard & Dupuy模型：模型：（1 1）基本公式）基本公式 当油藏压力高于饱和压力，而流动压力低于饱当油藏压力高于饱和压力，而流动压力低于饱和压力时，油藏中将同时存在单相和两相流动，拟稳和压力时，油藏中将同时存在单相和两相流动，拟稳态条件下产量的一般表达式为态条件下产量的一般表达式为: :三、三、 时的流入动态时的流入动态采油工程原理

18、与设计采油工程原理与设计井周围压力和气体饱和分布井周围压力和气体饱和分布采油工程原理与设计采油工程原理与设计 组合型组合型IPRIPR曲线曲线（2 2）实用计算方法）实用计算方法( (组合型组合型IPRIPR方法方法) )流压等于饱和压力时的产量为：流压等于饱和压力时的产量为： 当当 时，由于油藏中全部为时，由于油藏中全部为单相单相液体流动。液体流动。 当当 后，油藏中出现后，油藏中出现两相两相流动。流动。流入动态公式为：流入动态公式为：直线段直线段采油指数采油指数四、油气水三相四、油气水三相IPR IPR 曲线曲线PetrobrasPetrobras提出了计算三相流动提出了计算三相流动IPR

19、IPR曲线的方法。曲线的方法。综合综合IPRIPR曲线的实质曲线的实质: : 按按含水率含水率取纯油取纯油IPRIPR曲线和水曲线和水IPRIPR曲线的曲线的加权加权平均值平均值。当已知测试点。当已知测试点计算采液指数时，是按计算采液指数时，是按产量加权平均；当预测产量加权平均；当预测产量或流压时是按流压产量或流压时是按流压加权平均。加权平均。 图图1-12 油气水三相油气水三相IPR 曲线曲线( (一一) )采液指数计算采液指数计算已知一个测试点已知一个测试点( ( 、 ) )和饱和压力和饱和压力 及油藏压力及油藏压力 当当 时：时：思考题：推导思考题：推导 时的采液指数计算式。时的采液指数

20、计算式。 当当 时：时： 当当 时：时：其中：其中：直线段直线段采液指数采液指数 ，则：，则： ，则按流压加权平均进行推导：，则按流压加权平均进行推导： ( (二二) )某一产量下的流压计算某一产量下的流压计算所以：所以：因为：因为：所以：所以：若若 ，则综合，则综合IPRIPR曲线的斜率可近似曲线的斜率可近似为常数。为常数。思考题：试推导思考题：试推导五、五、多层油藏油井流入动态多层油藏油井流入动态（1 1）多油层油井流入动态）多油层油井流入动态迭加型迭加型IPRIPR图图1-13 1-13 多层油藏油井流入动态多层油藏油井流入动态（2 2）含水油井流入动态）含水油井流入动态图图1-14 1

21、-14 含水油井流入动态与含水变化含水油井流入动态与含水变化 ( ) ( )图图1-15 1-15 含水油井流入动态曲线含水油井流入动态曲线 ( ) ( )采油工程原理与设计采油工程原理与设计采油工程原理与设计采油工程原理与设计小小 结结(1) (1) 介绍的方法阐明了油井流入动态的物理意义，也是目介绍的方法阐明了油井流入动态的物理意义，也是目前现场最常用的计算方法。前现场最常用的计算方法。(2) (2) 油井流入动态研究主要有三种途径：油井流入动态研究主要有三种途径： 基于基于VogelVogel、FetkovichFetkovich、PetrobrasPetrobras方法的完善。方法的完

22、善。 建立不同类型油藏和井底条件的渗流模型。建立不同类型油藏和井底条件的渗流模型。 利用单井流入动态的油藏数值模拟技术。利用单井流入动态的油藏数值模拟技术。(3) (3) 油井流入动态是采油工程各项技术措施设计、分析与油井流入动态是采油工程各项技术措施设计、分析与评价的依据。评价的依据。第二节第二节 井筒气液两相流基本概念井筒气液两相流基本概念井筒多相流理论：井筒多相流理论：研究各种举升方式油井生产规律基本理论研究各种举升方式油井生产规律基本理论研究特点研究特点：流动复杂性、无严格数学解流动复杂性、无严格数学解研究途径研究途径：基本流动方程基本流动方程 实验资料相关因次分析实验资料相关因次分析

23、 近似关系近似关系一、一、井筒气液两相流动的特性井筒气液两相流动的特性( (一一) )气液两相流动与单相液流的比较气液两相流动与单相液流的比较流动型态（流动结构、流型）：流动型态（流动结构、流型）： 流动过程中油、气的分布状态。流动过程中油、气的分布状态。( (二二) )气液混合物在垂直管中的流动结构变化气液混合物在垂直管中的流动结构变化影影响响流流型型的的因因素素各相介质的体积比例各相介质的体积比例介质的流速介质的流速各相的物理及化学性质各相的物理及化学性质(密度、粘度界面张力等密度、粘度界面张力等)流道的几何形状流道的几何形状壁面特性壁面特性管道的安装方式管道的安装方式采油工程原理与设计采

24、油工程原理与设计流动型态的划分方法：两类流动型态的划分方法：两类第一类划分方法：第一类划分方法：根据两相介质分布的外形划分，根据两相介质分布的外形划分，包括泡包括泡状流、弹状流或团状流、（层状流、波状流）、段塞流或状流、弹状流或团状流、（层状流、波状流）、段塞流或冲击流、环状流、雾状流冲击流、环状流、雾状流 垂直气液两相流流型垂直气液两相流流型 水平气液两相流流型水平气液两相流流型采油工程原理与设计采油工程原理与设计第二类划分方法：第二类划分方法：按流动的数学模型或流体的分散程度划分，按流动的数学模型或流体的分散程度划分，包括包括分散流、间歇流、分离流分散流、间歇流、分离流 两种分类方法比较两

25、种分类方法比较 第一类划分方法较为直观第一类划分方法较为直观第二类划分方法便于进行数学处理第二类划分方法便于进行数学处理第第一一类类划划分分方方法法泡状流泡状流弹状流或团状流弹状流或团状流层状流层状流波状流波状流段塞流或冲击流段塞流或冲击流环状流环状流雾状流雾状流第第二二类类划划分分方方法法分散流分散流间歇流间歇流分离流分离流分离流分离流间歇流间歇流分离流分离流分散流分散流两类划分结果两类划分结果的对应关系的对应关系采油工程原理与设计采油工程原理与设计 纯液流纯液流 当井筒压力大于饱和压力时，天然气溶解在原油中，当井筒压力大于饱和压力时，天然气溶解在原油中，产液呈单相液流。产液呈单相液流。泡流

26、泡流 井筒压力稍低于饱和压力时，溶解气开始从井筒压力稍低于饱和压力时，溶解气开始从油中分离出来，气体都以小气泡分散在液相中。油中分离出来，气体都以小气泡分散在液相中。滑脱现象：滑脱现象： 混合流体流动过程中，由于流体间的密度混合流体流动过程中，由于流体间的密度差异，引起的小密度流体流速大于大密度流体差异，引起的小密度流体流速大于大密度流体流速的现象。流速的现象。 如：油气滑脱、气液滑脱、油水滑脱等。如：油气滑脱、气液滑脱、油水滑脱等。特点特点：气体是分散相，液体是连续相；气体是分散相，液体是连续相； 气体主要影响混合物密度，对摩擦阻力影响不大；气体主要影响混合物密度，对摩擦阻力影响不大； 滑脱

27、现象比较严重。滑脱现象比较严重。段塞流段塞流 当混合物继续向上流动，压力逐渐降低，当混合物继续向上流动，压力逐渐降低，气体不断膨胀，小气泡将合并成大气泡，直到气体不断膨胀，小气泡将合并成大气泡，直到能够占据整个油管断面时，井筒内将形成一段能够占据整个油管断面时，井筒内将形成一段液一段气的结构。液一段气的结构。特点：特点：气体呈分散相，液体呈连续相；气体呈分散相，液体呈连续相； 一段气一段液交替出现；一段气一段液交替出现； 气体膨胀能得到较好的利用；气体膨胀能得到较好的利用； 滑脱损失变小；滑脱损失变小； 摩擦损失变大。摩擦损失变大。环流 油管中心是连续的气流而管壁为油环的流动结构。特点：气液两

28、相都是连续相； 气体举油作用主要是靠摩擦携带； 滑脱损失变小； 摩擦损失变大。雾流 气体的体积流量增加到足够大时，油管中内流动的气流芯子将变得很粗，沿管壁流动的油环变得很薄，绝大部分油以小油滴分散在气流中。特点：气体是连续相，液体是分散相； 气体以很高的速度携带液滴喷出井口； 气、液之间的相对运动速度很小； 气相是整个流动的控制因素。总结： 油井生产中可能出现的流型自下而上依次为：纯油(液)流、泡流、段塞流、环流和雾流。 实际上，在同一口井内，一般不会出现完整的流型变化。图1-17 油气沿井筒喷出时的流型变化示意图纯油流；泡流；段塞流；环流；雾流实际计算：直接求存在滑脱混合物密度或包括滑脱在内

29、的摩实际计算：直接求存在滑脱混合物密度或包括滑脱在内的摩擦阻力系数。擦阻力系数。( (三三) )滑脱损失概念滑脱损失概念因滑脱而产生的附加压力损失称为滑脱损失因滑脱而产生的附加压力损失称为滑脱损失。单位管长上滑脱损失为：单位管长上滑脱损失为：图图1-18 1-18 气液两相流流动断面简图气液两相流流动断面简图滑脱损失的实质滑脱损失的实质： 液相的流动断面增大引起液相的流动断面增大引起混合物密度的增加。混合物密度的增加。二、井筒气液两相流能量平衡方程二、井筒气液两相流能量平衡方程及压力分布计算步骤及压力分布计算步骤 两个流动断面间的能量平衡关系：两个流动断面间的能量平衡关系：( (一一) )能量

30、平衡方程推导能量平衡方程推导图图2-19 2-19 倾斜管流能量平衡关系示意图倾斜管流能量平衡关系示意图图图2-19 2-19 倾斜管流能量平衡关系示意图倾斜管流能量平衡关系示意图倾斜多相管流断面倾斜多相管流断面1 1和断面和断面2 2的流体的能量平衡关系为：的流体的能量平衡关系为：适合于各种管流的通用压力梯度方程：适合于各种管流的通用压力梯度方程：则：则：令：令：采油工程原理与设计采油工程原理与设计以计算段下端压力为起点，重复以计算段下端压力为起点，重复步，计算下一段的深步，计算下一段的深度和压力，直到各段的累加深度等于管长为止。度和压力，直到各段的累加深度等于管长为止。(2)(2)多相垂直

31、管流压力分布计算步骤多相垂直管流压力分布计算步骤重复重复的计算，直至的计算，直至 。1)1)按深度增量迭代的步骤按深度增量迭代的步骤已知任一点已知任一点( (井口或井底井口或井底) )的压力作为起点，任选一个合适的压力作为起点，任选一个合适的压力降作为计算的压力间隔的压力降作为计算的压力间隔 p p。估计一个对应的深度增量估计一个对应的深度增量 h h 。计算该管段的平均温度及平均压力，并确定流体性质参数。计算该管段的平均温度及平均压力，并确定流体性质参数。判断流型，并计算该段的压力梯度判断流型，并计算该段的压力梯度dp/dhdp/dh。计算对应于的该段管长计算对应于的该段管长( (深度差深度

32、差) ) h h。计算该段下端对应的深度及压力。计算该段下端对应的深度及压力。采油工程原理与设计采油工程原理与设计2)2)按压力增量迭代的步骤（略）按压力增量迭代的步骤（略）思考题：思考题：根据上述步骤整理出计算压力分根据上述步骤整理出计算压力分布的程序流程框图。布的程序流程框图。说明：说明：a. 计算压力分布过程中，温度和压力是相关的；计算压力分布过程中，温度和压力是相关的；b. 流体物性参数计算至关重要，但目前方法精度差；流体物性参数计算至关重要，但目前方法精度差；c. 不同的多相流计算方法差别较大，因此在实际应用中不同的多相流计算方法差别较大，因此在实际应用中有必要根据油井的实际情况筛选

33、精度相对高的方法。有必要根据油井的实际情况筛选精度相对高的方法。第三节第三节 OrkiszewskiOrkiszewski方法方法 综合了综合了Griffith & Wallis Griffith & Wallis 和和 Duns & Ros Duns & Ros 等等方法方法 处理过渡性流型时，采用处理过渡性流型时，采用RosRos方法方法( (内插法内插法) ) 针对每种流动型态提出存容比及摩擦损失的计算方法针对每种流动型态提出存容比及摩擦损失的计算方法 提出了四种流型，即泡流、段塞流、过渡流及环雾流提出了四种流型，即泡流、段塞流、过渡流及环雾流 把把GriffithGriffith段塞流

34、相关式改进后推广到了高流速区段塞流相关式改进后推广到了高流速区 19671967年提出，适用于垂直管流计算年提出，适用于垂直管流计算采油工程原理与设计采油工程原理与设计出现雾流时，气体体积流量远大于液体体积流量。根据气出现雾流时，气体体积流量远大于液体体积流量。根据气体定律，动能变化可表示为：体定律，动能变化可表示为：一、一、压力降公式及流动型态划分界限压力降公式及流动型态划分界限由垂直管流能量方程可知，压力降是摩擦能量损失、势能变由垂直管流能量方程可知，压力降是摩擦能量损失、势能变化和动能变化之和：化和动能变化之和：所以压降计算式为：所以压降计算式为：在雾流情况下才有明显意义在雾流情况下才有

35、明显意义表1-3 Orkiszewski方法流型划分界限不同流动型态下 和 的计算方法不同。二、平均密度及摩擦损失梯度的计算二、平均密度及摩擦损失梯度的计算气相存容比气相存容比( (含气率含气率) )H Hg g ：管段中气相体积与管段容管段中气相体积与管段容积之比值。积之比值。液相存容比液相存容比( (持液率持液率) )H HL L ：管段中液相体积与管段容管段中液相体积与管段容积之比值。积之比值。(1)(1)泡流泡流平均密度平均密度: :？滑脱速度滑脱速度：气相流速与液相流速之差气相流速与液相流速之差。则：则：泡流摩擦损失梯度泡流摩擦损失梯度按液相进行计算：按液相进行计算：Griffith

36、Griffith实验测得：实验测得：0.244m/s0.244m/s图图1-21 1-21 摩擦阻力系数曲线（教材摩擦阻力系数曲线（教材p37p37）图1-21(2)(2)段塞流段塞流平均密度平均密度: :段塞流的摩擦梯度：段塞流的摩擦梯度：段塞流计算关键是滑脱速度段塞流计算关键是滑脱速度vs和液体分布系数和液体分布系数的计算。的计算。液体分布系数液体分布系数的计算见课本的计算见课本p39表表1-4和公式（和公式（1-58）。）。滑脱速度滑脱速度vs的计算方法有两种：查图迭代法和经验公式法。的计算方法有两种：查图迭代法和经验公式法。采油工程原理与设计采油工程原理与设计泡流雷诺数：泡流雷诺数：图

37、图1-22 1-22 C C1 1N Nb b曲线曲线雷诺数：雷诺数：图图1-23 1-23 C C2 2N NReRe曲线曲线滑脱速度的计算滑脱速度的计算查图迭代法查图迭代法滑脱速度：滑脱速度：采油工程原理与设计采油工程原理与设计滑脱速度的计算滑脱速度的计算经验公式计算法经验公式计算法详见教材详见教材p3893公式（公式（1-54）（1-57）。）。(3)(3)过渡流过渡流 过渡流的混合物平均密度及摩擦梯度是先按段塞流和雾流过渡流的混合物平均密度及摩擦梯度是先按段塞流和雾流分别进行计算，然后用内插方法来确定相应的数值。分别进行计算，然后用内插方法来确定相应的数值。雾流混合物雾流混合物平均密度

38、平均密度计算公式与泡流相同：计算公式与泡流相同：由于雾流的气液无相对运动速度，即滑脱速度接近于零，由于雾流的气液无相对运动速度，即滑脱速度接近于零，基本上没有滑脱。基本上没有滑脱。雾流摩擦系数可根据气体雷诺数和液膜相对粗糙度查图得。雾流摩擦系数可根据气体雷诺数和液膜相对粗糙度查图得。摩擦梯度摩擦梯度: :(4)(4)雾流雾流所以：所以：图图1-24 1-24 OrkiszewskiOrkiszewski方法计算流程框图方法计算流程框图第四节第四节 Beggs & Brill Beggs & Brill 方法方法 水和空气、聚丙烯管实验基础上总结的水和空气、聚丙烯管实验基础上总结的方法方法 建立

39、流型分布图，将七种流型归为三类，增加了过渡流建立流型分布图，将七种流型归为三类，增加了过渡流 计算时先按水平管流计算，然后采用倾斜校正系数校计算时先按水平管流计算，然后采用倾斜校正系数校正成相应的倾斜管流正成相应的倾斜管流 倾斜度倾斜度 -90 -90 +90+90，分上坡和下坡流动，分上坡和下坡流动 19731973年提出，适用于年提出，适用于水平、垂直和任意倾斜水平、垂直和任意倾斜管流计算管流计算Beggs & Brill Beggs & Brill 两相水平管流型两相水平管流型分离流分离流分层流分层流波状流波状流环状流环状流间歇流间歇流团状流团状流段塞流段塞流分散流分散流泡泡 流流雾雾

40、流流一、基本方程一、基本方程 单位质量气液混合物稳定流动的机械能量守恒方程为：单位质量气液混合物稳定流动的机械能量守恒方程为：(1)(1)位差压力梯度位差压力梯度：消耗于混合物静水压头的压力梯度。：消耗于混合物静水压头的压力梯度。(2)(2)摩擦压力梯度摩擦压力梯度：克服管壁流动阻力消耗的压力梯度。：克服管壁流动阻力消耗的压力梯度。假设条件假设条件: :气液混合物既未对外作功，也未受外界功。气液混合物既未对外作功，也未受外界功。(3)(3)加速度压力梯度加速度压力梯度：由于动能变化而消耗的压力梯度。：由于动能变化而消耗的压力梯度。 忽略液体压缩性、考虑到气体质量流速变化远远小于气体忽略液体压缩

41、性、考虑到气体质量流速变化远远小于气体密度变化，则：密度变化，则：(4)(4)总压力梯度（总压力梯度（Beggs-BrillBeggs-Brill方法的基本方程方法的基本方程）图图1-26 1-26 Beggs-BrillBeggs-Brill流型分布图（教材流型分布图（教材p p4545）二、二、 Beggs & BrillBeggs & Brill方法的流型分布图及流型判别式方法的流型分布图及流型判别式分离流分离流间歇流间歇流过渡流过渡流分散流分散流表表2-4 2-4 Beggs-BrillBeggs-Brill法流型判别条件法流型判别条件三、持液率及混合物密度确定三、持液率及混合物密度确

42、定(1)(1)持液率持液率Beggs & BrillBeggs & Brill方法计算倾斜管流时首先按水平管计方法计算倾斜管流时首先按水平管计算算, ,然后进行倾斜校正。然后进行倾斜校正。表表1-6 1-6 a a、b b、c c常数表常数表 实验结果表明，倾斜校正系数与倾斜角、无滑实验结果表明，倾斜校正系数与倾斜角、无滑脱持液率、弗洛德数及液体速度数有关。脱持液率、弗洛德数及液体速度数有关。图图1-27 1-27 不同不同E EL L下的倾斜校正系数下的倾斜校正系数根据实验结果回归的倾斜校正系数的相关式为根据实验结果回归的倾斜校正系数的相关式为: :对于垂直管对于垂直管: :系数系数C C与

43、无滑脱持液率与无滑脱持液率 、弗洛德数和液相速度数有关。、弗洛德数和液相速度数有关。表表1-6 1-6 系数系数d d、e e、f f、g g其中：其中：对于过渡流型对于过渡流型，先分别用分离流和间歇流计算，之后采用先分别用分离流和间歇流计算，之后采用内插法确定其持液率。内插法确定其持液率。利用持液率计算流动条件下混合物实际密度：利用持液率计算流动条件下混合物实际密度：四、阻力系数四、阻力系数 气液两相流阻力系数与无滑脱气液两相流阻力系数的比气液两相流阻力系数与无滑脱气液两相流阻力系数的比值与持液率和无滑脱持液率值与持液率和无滑脱持液率( (入口体积含液率入口体积含液率) )之间的关系：之间的

44、关系：当当11y1.2y1.2时时其中：其中：两相流动的雷诺数：两相流动的雷诺数： 也可用也可用MoodyMoody图上的光滑管曲线来确定图上的光滑管曲线来确定: :气液两相流阻力系数气液两相流阻力系数: :Beggs & BrillBeggs & Brill方法计算流程框图方法计算流程框图( (p49)p49)小小 结结(1) (1) 模拟计算多相管流规律的数学相关式及图版研究很多。纵模拟计算多相管流规律的数学相关式及图版研究很多。纵观这许多数学相关式，其基本通式一般都是从基本能量守恒方观这许多数学相关式，其基本通式一般都是从基本能量守恒方程出发建立的：程出发建立的：(2) (2) 对对Po

45、ettmann-CarpenterPoettmann-Carpenter方法、方法、Fanch-BrownFanch-Brown相关式、相关式、Baxendell-Thomas Baxendell-Thomas 相关式、相关式、Hagedron-Brown Hagedron-Brown 关系式、关系式、Duns-Duns-RosRos相关式、相关式、OrkiszewskiOrkiszewski相关式、相关式、Beggs-Brill Beggs-Brill 相关式、相关式、DuklerDukler相关式、相关式、Mukherjee-BrillMukherjee-Brill相关式、相关式、 Azi

46、zAziz相关式、相关式、EatonEaton相关式、相关式、AnsariAnsari相关式等十二种方法进行了对比分析，不同的相关式等十二种方法进行了对比分析，不同的方法有其适用条件和精度，可根据具体油田实际选用。方法有其适用条件和精度，可根据具体油田实际选用。采油工程原理与设计采油工程原理与设计(3) (3) 数学相关式大体分为三种类型：数学相关式大体分为三种类型：在在计计算算井井筒筒流流体体混混合合密密度度时时不不考考虑虑液液体体滞滞留留量量的的影影响响，而而液液体体滞滞留留量量与与管管壁壁摩摩阻阻损损失失用用一一个个经经验验摩摩阻阻系系数数来来表表达达，不不区分多相流体的流态分布情况。区

47、分多相流体的流态分布情况。在在计计算算多多相相流流体体混混合合物物密密度度时时考考虑虑液液体体滞滞留留量量的的影影响响，而而摩阻系数要依据液体和气体的组成特征来确定摩阻系数要依据液体和气体的组成特征来确定考考虑虑液液体体滞滞留留量量的的影影响响，摩摩阻阻系系数数取取决决于于多多相相流流体体中中连连续续相相的的特特征征。同同时时考考虑虑了了不不同同流流态态的的影影响响，其其中中包包括括泡泡状状流流、段塞流、环流及雾状流。段塞流、环流及雾状流。END采油工程原理与设计采油工程原理与设计作业1.1.IPR曲线的概念；不同流态的曲线形状，画图说明。2.IPR曲线的制作方法。3.采油指数的概念，制作方法；指数如何计算？达西与非达西的计算式区别在哪里？4.Vogel方法的实质是什么？实际应用的计算过程，通过实例说明其曲线绘制方法。5.泄油面积？区别“圆形封闭边界”和 “稳压边界”及“无限大地层”的压力含义。流动效率概念。6.滑脱损失？其实质？产生的原因？存在哪些流态中？7.井筒多相流的流态包括几种模式？自制表格对比四种流态特征？8.“Vogel等计算方法”的“文献查询”（关键词：无因次IPR曲线，多相流分析方法，采等油指数的计算方法(改变教学模式，建立“自主学习”的教学方式). 9.思考题：推导时的采液指数计算式。思考题：推导时的采液指数计算式。