井筒多相流课件

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,0,井筒多相流理论,：,研究各种举升方式油井生产规律基本理论,研究特点,：,流动复杂性、无严格数学解,研究途径,：,基本流动方程,实验资料相关因次分析,近似关系,第二节,井筒气液两相流基本概念,第1页/共42页,井筒多相流理论：研究特点：流动复杂性、无严格数学解研究途径：,1,一、,井筒气液两相流动的特性,(一)气液两相流动与单相液流的比较,Comparison between single phase and two phase flow,第2页/共42页,一、井筒气液两相流动的特性(一)气液两相流动与单相液流的比较,2,流动型态（流动结构、流型）：,流动过程中油、气的分布状态。,(二)气液混合物在垂直管中的流动结构变化,Flow Regime,纯液流,Liquid flow,当井筒压力大于饱和压力时，天然气溶解在原油中，产液呈单相液流。,影响流型的因素：,气液体积比、流速、气液界面性质等。,第3页/共42页,流动型态（流动结构、流型）：(二)气液混合物在垂直管中的流动,3,泡流,Bubble Flow,井筒压力稍低于饱和压力时，溶解气开始从油中分离出来，气体都以小气泡分散在液相中。,滑脱现象：,Slippage,混合流体流动过程中，由于流体间的密度差异，引起的小密度流体流速大于大密度流体流速的现象。,如：油气滑脱、气液滑脱、油水滑脱等。,特点,：,气体是分散相，液体是连续相；,气体主要影响混合物密度，对摩擦阻力影响不大；,滑脱现象比较严重。,第4页/共42页,泡流 Bubble Flow滑脱现象：Slippage 特,4,段塞流,Slug Flow,当混合物继续向上流动，压力逐渐降低，气体不断膨胀，小气泡将合并成大气泡，直到能够占据整个油管断面时，井筒内将形成一段液一段气的结构。,特点：,气体呈分散相，液体呈连续相；,一段气一段液交替出现；,气体膨胀能得到较好的利用；,滑脱损失变小；,摩擦损失变大。,第5页/共42页,段塞流 Slug Flow特点：气体呈分散相，液体呈连续相,5,环流,Circular Flow,油管中心是连续的气流而管壁为油环的流动结构。,特点：,气液两相都是连续相；,气体举油作用主要是靠摩擦携带；,滑脱损失变小；,摩擦损失变大。,第6页/共42页,环流 Circular Flow特点：气液两相都是连续相；,6,雾流,Mist Flow,气体的体积流量增加到足够大时，油管中内流动的气流芯子将变得很粗，沿管壁流动的油环变得很薄，绝大部分油以小油滴分散在气流中。,特点,：,气体是连续相，液体是分散相；,气体以很高的速度携带液滴喷出井口；,气、液之间的相对运动速度很小；,气相是整个流动的控制因素。,第7页/共42页,雾流 Mist Flow特点：气体是连续相，液体是分散相；,7,总结：,油井生产中可能出现的流型自下而上依次为：纯油(液)流、泡流、段塞流、环流和雾流。,实际上，在同一口井内，一般不会出现完整的流型变化。,图1-17 油气沿井筒喷出时的流型变化示意图,纯油流；,泡流；,段塞流；,环流；,雾流,第8页/共42页,总结：图1-17 油气沿井筒喷出时的流型变化示意图第8页/,8,实际计算：直接求存在滑脱混合物密度或包括滑脱在内的摩擦阻力系数。,(三)滑脱损失概念,因滑脱而产生的附加压力损失称为滑脱损失,。,Slippage pressure drop,单位管长上滑脱损失为：,图1-18 气液两相流流动断面简图,滑脱损失的实质：,液相的流动断面增大引起混合物密度的增加。,第9页/共42页,实际计算：直接求存在滑脱混合物密度或包括滑脱在内的摩擦阻力系,9,二、井筒气液两相流能量平衡方程,及压力分布计算步骤,Wellbore Pressure Profile Calculation,两个流动断面间的能量平衡关系：,(一)能量平衡方程推导,图2-19 倾斜管流能量平衡关系示意图,具有能量：内能、位能、动能、膨胀能,第10页/共42页,二、井筒气液两相流能量平衡方程 两个流动断面间的能量平衡关系,10,图2-19 倾斜管流能量平衡关系示意图,倾斜多相管流断面1和断面2的流体的能量平衡关系为：,第11页/共42页,图2-19 倾斜管流能量平衡关系示意图倾斜多相管流断面1和,11,适合于各种管流的通用压力梯度方程：,则：,令：,第12页/共42页,适合于各种管流的通用压力梯度方程：则：令：第12页/共42页,12,井筒多相垂直管流压力分布,图1-17 油气沿井筒喷出时的流型变化示意图,纯油流；,泡流；,段塞流；,环流；,雾流,压力,计算过程复杂,压力计算与流体物性参数有关,流体物性参数等是压力的函数,工程上采用迭代方法编程计算,按深度增量和压力增量迭代,学习关键是掌握计算的原理,未知数：,密度、流速、摩擦阻力系数,第13页/共42页,井筒多相垂直管流压力分布图1-17 油气沿井筒喷出时的流型,13,以计算段下端压力为起点，重复步，计算下一段的深,度和压力，直到各段的累加深度等于管长为止。,(2)多相垂直管流压力分布计算步骤,重复的计算，直至 。,1)按深度增量迭代的步骤,已知任一点(井口或井底)的压力作为起点，任选一个合适,的压力降作为计算的压力间隔,p,。,估计一个对应的深度增量,h,。,计算该管段的平均温度及平均压力，并确定流体性质参数。,判断流型，并计算该段的压力梯度,dp/dh。,计算对应于的该段管长(深度差),h,。,计算该段下端对应的深度及压力。,第14页/共42页,以计算段下端压力为起点，重复步，计算下一段的深(2),14,2)按压力增量迭代的步骤（略）,思考题：,根据上述步骤整理出计算压力分布的程序流程框图。,说明：,a.,计算压力分布过程中，温度和压力是相关的；,b.,流体物性参数计算至关重要，但目前方法精度差；,c.,不同的多相流计算方法差别较大，因此在实际应用中有必要根据油井的实际情况筛选精度相对高的方法。,第15页/共42页,2)按压力增量迭代的步骤（略）思考题：根据上述步骤整理出计算,15,第三节,Orkiszewski,方法,综合了,Griffith&Wallis,和,Duns&Ros,等,方法,处理过渡性流型时，采用,Ros,方法(内插法),针对每种流动型态提出存容比及摩擦损失的计算方法,提出了四种流型，即泡流、段塞流、过渡流及环雾流,把,Griffith,段塞流相关式改进后推广到了高流速区,于,1967年提出，适用于垂直管流计算,第16页/共42页,第三节 Orkiszewski方法 综合了Griffith,16,图1-24,Orkiszewski,方法计算流程框图,第17页/共42页,图1-24 Orkiszewski方法计算流程框图第17页/,17,出现雾流时，气体体积流量远大于液体体积流量。根据气体定律，动能变化可表示为：,一、,压力降公式及流动型态划分界限,由垂直管流能量方程可知，压力降是摩擦能量损失、势能变化和动能变化之和：,所以压降计算式为：,未知数,：,第18页/共42页,出现雾流时，气体体积流量远大于液体体积流量。根据气体定律，动,18,表1-3,Orkiszewski,方法流型划分界限,不同流动型态下 和 的计算方法不同。,第19页/共42页,表1-3 Orkiszewski方法流型划分界限不同流动型态,19,二、平均密度及摩擦损失梯度的计算,气相存容比(,含气率,),H,g,：,管段中气相体积与管段容积之比值。,液相存容比(,持液率,),H,L,：,管段中液相体积与管段容积之比值。,(1)泡流,平均密度:,求混合物密度需要先计算含气率或持液率,泡流滑脱现象严重，若无滑脱：,第20页/共42页,二、平均密度及摩擦损失梯度的计算气相存容比(含气率)Hg：,20,滑脱速度,：,气相流速与液相流速之差,。,（由实验测定）,则：,泡流摩擦损失梯度,按液相进行计算：,真实流速与表观流速的关系：,第21页/共42页,滑脱速度：气相流速与液相流速之差。（由实验测定）则：泡流摩擦,21,图1-21 摩擦阻力系数曲线,（教材,p37,）,图1-21,第22页/共42页,图1-21 摩擦阻力系数曲线（教材p37）图1-21第22,22,(2)段塞流,平均密度,:,段塞流的摩擦梯度：,段塞流计算中，关键是滑脱速度,v,s,的计算。,目前，,v,s,的计算方法有两种：查图迭代法和经验公式法。,第23页/共42页,(2)段塞流平均密度:段塞流的摩擦梯度：段塞流计算中，关键是,23,泡流雷诺数：,图1-22,C,1,N,b,曲线,雷诺数：,图1-23,C,2,N,Re,曲线,滑脱速度的计算,迭代法,滑脱速度：,第24页/共42页,泡流雷诺数：图1-22 C1Nb曲线雷诺数：图1-23,24,滑脱速度的计算,经验公式计算法,详见教材,p3893,公式（1-54）（1-58）。,(3)过渡流,过渡流的混合物平均密度及摩擦梯度是先按段塞流和雾流分别进行计算，然后用内插方法来确定相应的数值。,第25页/共42页,滑脱速度的计算经验公式计算法详见教材p3893公式（1,25,雾流混合物,平均密度,计算公式与泡流相同：,由于雾流的气液无相对运动速度，即滑脱速度接近于零，基本上没有滑脱。,雾流摩擦系数可根据气体雷诺数和液膜相对粗糙度查图得。,摩擦梯度:,(4)雾流,所以：,第26页/共42页,雾流混合物平均密度计算公式与泡流相同：由于雾流的气液无相对运,26,图1-24,Orkiszewski,方法计算流程框图,第27页/共42页,图1-24 Orkiszewski方法计算流程框图第27页/,27,第四节,Beggs&Brill,方法,水和空气、聚丙烯管实验基础上总结的,方法,建立流型分布图，将七种流型归为三类，增加了过渡流,计算时先按水平管流计算，然后采用倾斜校正系数校,正成相应的倾斜管流,倾斜度-90,+90，分上坡和下坡流动,1973年提出，适用于,水平、垂直和任意倾斜,管流计算,第28页/共42页,第四节 Beggs&Brill 方法 水和空气、聚丙烯管,28,Beggs&Brill,两相水平管流型,分离流,分层流,波状流,环状流,间歇流,团状流,段塞流,分散流,泡 流,雾 流,第29页/共42页,Beggs&Brill 两相水平管流型分离流分层流波状流,29,一、基本方程,单位质量气液混合物稳定流动的机械能量守恒方程为：,(1)位差压力梯度,：消耗于混合物静水压头的压力梯度。,(2)摩擦压力梯度,：克服管壁流动阻力消耗的压力梯度。,假设条件,:气液混合物既未对外作功，也未受外界功。,第30页/共42页,一、基本方程 单位质量气液混合物稳定流动的机械能量守恒方程为,30,(3)加速度压力梯度,：由于动能变化而消耗的压力梯度。,忽略液体压缩性、考虑到气体质量流速变化远远小于气体密度变化，则：,(4)总压力梯度（,Beggs-Brill,方法的基本方程,）,第31页/共42页,(3)加速度压力梯度：由于动能变化而消耗的压力梯度。忽略液体,31,图1-26,Beggs-Brill,流型分布图（教材,p,45）,二、,Beggs&Brill,方法的流型分布图及流型判别式,第32页/共42页,图1-26 Beggs-Brill流型分布图（教材p45,32,表2-4,Beggs-Brill,法流型判别条件,第33页/共42页,表2-4 Beggs-Brill法流型判别条件第33页/共,33,三、持液率及混合物密度确定,(1)持液率,Beggs&Brill,方法计算倾斜管流时首先按水平管计算,然后进行倾斜校正。,表1-6,a、b、c,常数表,第34页/共42页,三、持液率及混合物密度确定(1)持液率表1-6 a、b、c,34,实验结果表明，倾斜校正系数与倾斜角、无滑脱持液率、弗洛德数及液体速度数有关。,图1-27 不同,E,L,下的倾斜校正系数,第35页/共42页,实验结果表明，倾斜校正系数与倾斜角、无滑脱持液率、