1、多相管流理论与计算课程总学时：课程总学时：3232学时学时课程主体和基础内容：课程主体和基础内容：李宾飞李宾飞 张黎明张黎明 课程前沿及延伸部分：课程前沿及延伸部分：李明忠李明忠多相流在油气田开发中的应用（多相流在油气田开发中的应用（2-42-4学时）学时）李兆敏李兆敏泡沫流体流动规律研究及应用（泡沫流体流动规律研究及应用（2 2学时）学时）王卫阳王卫阳多相流测试技术（多相流测试技术（4 4学时）学时）课程安排课程安排 第一章 概 论第一节 引言多相流理论多相流理论多相流体力学理论多相流体力学理论相相 相是体系中具有相同化学组成和物理性质的一部相是体系中具有相同化学组成和物理性质的一部 分，分
2、，与体系的其它均匀部分由界面隔开与体系的其它均匀部分由界面隔开多相流动体系多相流动体系 两相或两相以上的流体一同参与流动的体系两相或两相以上的流体一同参与流动的体系多相管流多相管流多相多相流体在流体在管道管道中的流动中的流动Multiphase FlowMultiphase Flow多相流动多相流动 两相或两相以上的流体一同参与的流动两相或两相以上的流体一同参与的流动油水混合物：油水混合物：从相的角度从相的角度 区分为区分为 油相、水相油相、水相 从物质形态看从物质形态看同属液态同属液态 相相物质形态物质形态(气态、液态、固态气态、液态、固态)盐水：盐水：从相的角度划分只有从相的角度划分只有
3、液相液相 从组分划分可分为：从组分划分可分为：水和水和NaCl 相相组分组分(指混合物中的各个成分指混合物中的各个成分)原油？原油？第一节 引言 水水-冰冰 液固两相流液固两相流 水水-水蒸汽水蒸汽 气液两相流气液两相流 泥浆泥浆 液固两相流液固两相流 油油-气气 气液两相流气液两相流 油油-气气-水水 气液液三相流气液液三相流 单工质单工质（相同化学成分）（相同化学成分）多工质多工质（不同化学成分）（不同化学成分）多相流体流动遵循的规律与单相流体并不相同，需引入新多相流体流动遵循的规律与单相流体并不相同，需引入新的理论来反映多相流体流动规律。的理论来反映多相流体流动规律。多相体系：多相体系：
4、水水-冰、水冰、水水蒸汽、泥浆、油气、油气水等水蒸汽、泥浆、油气、油气水等多相体系的类型多相体系的类型第一节 引言 热能工程：热能工程：锅炉系统、制冷系统、热管锅炉系统、制冷系统、热管 航天技术：航天技术：平衡温差、保护设备平衡温差、保护设备 核工业：核工业：汽液两相流动汽液两相流动 石油工业：石油工业：两相渗流计算、井控设计、采油工艺设计、两相渗流计算、井控设计、采油工艺设计、油气集输油气集输 化工行业：化工行业：工艺设计工艺设计 其它行业：其它行业：水利、粉状物管线输送水利、粉状物管线输送多相流应用领域多相流应用领域第一节 引言油气是深埋于地下的油气是深埋于地下的混合流体混合流体矿藏，因此
5、，油气藏的开发与矿藏，因此，油气藏的开发与开采离不开流体力学理论及其分枝开采离不开流体力学理论及其分枝多相流理论。多相流理论。举例来说，举例来说，油气井压力控制、油气管流计算、举升参数设计、油气井压力控制、油气管流计算、举升参数设计、工况分析、集输设计等，工况分析、集输设计等，都离不开多相流的理论与计算方法。都离不开多相流的理论与计算方法。多相流理论是贯穿于石油开采全过程的基本理论一、多相流理论在石油工业中的地位和作用一、多相流理论在石油工业中的地位和作用石油工程石油工程(油气井工程和油气田开发工程油气井工程和油气田开发工程)以及以及油气储运工油气储运工程程都与多相流理论有着极为密切的联系。都
6、与多相流理论有着极为密切的联系。钻井工程钻井工程：油气井压力控制油气井压力控制 (含气泥浆的压降计算含气泥浆的压降计算)采油工程采油工程：采油方式优化设计、采油设备的工况分析采油方式优化设计、采油设备的工况分析 （油气混合物在井筒中流动的压降和温度计算）（油气混合物在井筒中流动的压降和温度计算）储运工程储运工程：油气集输管线的设计油气集输管线的设计 （油气混合物在管线中压降和温度的计算）（油气混合物在管线中压降和温度的计算）许多工程设计都将计算多相流体在管道中流动的许多工程设计都将计算多相流体在管道中流动的压降和温度压降和温度。一、多相流理论在石油工业中的地位和作用一、多相流理论在石油工业中的
7、地位和作用自喷采油示意图(1)油井能否自喷？油井能否自喷？(2)一定产量的油气混合物到一定产量的油气混合物到达井口时的剩余压力？达井口时的剩余压力？tpwftppptpwfptpp井口油压井口油压井底流压井底流压两相流压降两相流压降htpp htpp 自喷生产自喷生产机械采油机械采油(人工举升人工举升)自喷采油设计中的两相流计算自喷采油设计中的两相流计算 气举采油系统示意图气举采油系统示意图 依靠从地面注依靠从地面注入井内的高压气体入井内的高压气体与油层产出流体在与油层产出流体在井筒中混合，井筒中混合，利用利用气体的膨胀使井筒气体的膨胀使井筒中的混合液密度降中的混合液密度降低，低，将流到井内的
8、将流到井内的原油举升到地面。原油举升到地面。气举采油设计中的两相流计算气举采油设计中的两相流计算气举采油及井筒压力分布设计的原则设计的原则：最大限度地发挥油藏的潜最大限度地发挥油藏的潜力和地面设备的能力，获力和地面设备的能力，获得最高的产油量。得最高的产油量。tpwftppp气举采油设计中的两相流计算气举采油设计中的两相流计算常规有杆泵生产地面驱动螺杆泵生产有杆泵采油设计中的两相流计有杆泵采油设计中的两相流计算算 近年来，随着陆地油气资源的日趋减少，石油和天然气勘近年来，随着陆地油气资源的日趋减少，石油和天然气勘探开发已转向广阔的海洋。早期开发的海上油田都位于浅海水探开发已转向广阔的海洋。早期
9、开发的海上油田都位于浅海水域，采用的是固定平台采油系统。域，采用的是固定平台采油系统。1965 1965 年前，钻探活动一般年前，钻探活动一般小于小于90m 90m 水深。水下采油技术的发展不仅使深海石油开采成为水深。水下采油技术的发展不仅使深海石油开采成为可能，同时大大降低了深海采油成本，目前深水可能，同时大大降低了深海采油成本，目前深水(4001500(4001500米米)和超深水和超深水(超过超过15001500米米)海域的油气开发，已成为美国、英国、海域的油气开发，已成为美国、英国、挪威、巴西等国竞相开展的热点。挪威、巴西等国竞相开展的热点。海洋及深水采油中的多相流计算海洋及深水采油中
10、的多相流计算典型的海洋采油系统固定平台张力腿平台半潜式平台海底管线海底管线浮式生产储油（浮式生产储油（FPSOFPSO）船）船 单点系泊单点系泊单点系泊单点系泊 但当水深达到但当水深达到2000ft(约约600m)时，其温度将下降至时，其温度将下降至40(4.4)左右。油气左右。油气混合物的热量大量散失，温度会迅速降低，原油和油水乳状液的粘度会明显混合物的热量大量散失，温度会迅速降低，原油和油水乳状液的粘度会明显增加，天然气水合物、蜡、沥青质等固态物质也可能在管道中析出并沉积。增加，天然气水合物、蜡、沥青质等固态物质也可能在管道中析出并沉积。流动保障流动保障 Flow Assurance海洋及
11、深水采油中的多相流计算海洋及深水采油中的多相流计算“流动保障流动保障”确保油气的无阻塞流动并使系统的运行费用达确保油气的无阻塞流动并使系统的运行费用达到最低。到最低。Pipe-in-Pipe保温材料海洋及深水采油中的多相流计算海洋及深水采油中的多相流计算管线管束(flowline bundles)海洋及深水采油中的多相流计算海洋及深水采油中的多相流计算渤海平均水深渤海平均水深18m,最深最深83m 黄海平均水深黄海平均水深44m,最深最深140m 东海平均水东海平均水深深370m,最深最深2719m 南海平均水深南海平均水深1212m,最深最深5377m石油工业中的多相流技术研究始于石油工业中
12、的多相流技术研究始于1950年左右年左右,早期研究者早期研究者大多采用实验研究的方法，所使用的数据主要来自室内实验大多采用实验研究的方法，所使用的数据主要来自室内实验和油田实际生产和油田实际生产 80年代初，计算机的引入极大地促进了多相流的发展年代初，计算机的引入极大地促进了多相流的发展二、多相流理论的研究简史控制方程（机理模型）物理现象流动机理70年代，石油工业开始采用已在其它工业领域中使用的一些年代，石油工业开始采用已在其它工业领域中使用的一些物理机理来预测多相流的流型物理机理来预测多相流的流型80年代中期应用高新技术及仪器进行多相流的模拟试验，年代中期应用高新技术及仪器进行多相流的模拟试
13、验，期望深入认识多相流动现象及流动机理，从而改进模型，提期望深入认识多相流动现象及流动机理，从而改进模型，提高精度。高精度。核密度计、超声波传感器、电导和光导探针、电容传感器、核密度计、超声波传感器、电导和光导探针、电容传感器、激光多普勒测速仪、高速摄像机等。激光多普勒测速仪、高速摄像机等。西安交通大学西安交通大学 动力工程多相流国家重点实验室动力工程多相流国家重点实验室 目前，目前，双流体瞬态模拟双流体瞬态模拟方法和精确描述物理现象的方法和精确描述物理现象的稳态机稳态机理模型理模型是多相管流研究的主要方法是多相管流研究的主要方法二、多相流理论的研究简史1、多相流问题未得到解析解、多相流问题未
14、得到解析解2、油气水三相流的研究不过深入、油气水三相流的研究不过深入3、水平井段变质量流动研究较少、水平井段变质量流动研究较少5、缺乏专用研究仪器、缺乏专用研究仪器4、缺乏向下流动的综合机理模型、缺乏向下流动的综合机理模型三、目前存在的问题 水平井示意图水平井示意图前次课回顾相的概念相的概念 相是体系中具有相同化学组成和物理性质的一相是体系中具有相同化学组成和物理性质的一部分部分 ，与体系的其它均匀部分有界面隔开，与体系的其它均匀部分有界面隔开相相物质形态物质形态(气态、液态、固态气态、液态、固态)例如：水-冰 液固两相流 水-水蒸气 气液两相流 泥浆 液固两相流 油-气 气液两相流 油-气-
15、水 气液液三相流 单工质单工质多工质多工质相相组分组分(指混合物中的各个成分指混合物中的各个成分 )气液两相流基本概念气液两相流基本概念 气液两相流的模型气液两相流的模型 油藏流体高压物性的计算油藏流体高压物性的计算 多相流体的温度分布计算多相流体的温度分布计算 垂直气液两相管流压降计算垂直气液两相管流压降计算 水平气液两相管流水平气液两相管流压降计算压降计算 倾斜气液两相管流的压降计算倾斜气液两相管流的压降计算 专题讲座专题讲座课程主要内容课程主要内容 第二节 气液两相管流的基本特征与研究方法细分散体系细分散体系 细小的液滴或气泡均匀分散在连续相中细小的液滴或气泡均匀分散在连续相中 粗分散体
16、系粗分散体系 较大的气泡或液滴分散在连续相中较大的气泡或液滴分散在连续相中 混合流动型混合流动型 两相均非连续相两相均非连续相 分层流动分层流动 两相均为连续相两相均为连续相 一、基本特征1 1气液两相流的分类气液两相流的分类 Gas-liquid two phase flowGas-liquid two phase flow体系中存在相界面体系中存在相界面 两相之间也存在力的作用，出现质两相之间也存在力的作用，出现质量和能量的交换时伴随着机械能的损失量和能量的交换时伴随着机械能的损失 两相的分布情况多种多样两相的分布情况多种多样 两相流动中两相介质的分布两相流动中两相介质的分布状况称为流型状
17、况称为流型 两相流动中存在滑脱现象两相流动中存在滑脱现象 相间速度的差异称为滑脱相间速度的差异称为滑脱，滑脱将产生附加的能量损失滑脱将产生附加的能量损失 沿程流体体积流量有很大变化，质量流量不变沿程流体体积流量有很大变化，质量流量不变2 气液两相流的基本特征气液两相流的基本特征一、基本特征垂直气液两相流流型 水平气液两相流流型多相流体的主要流型多相流体的主要流型一、基本特征不能简单地用层流或紊流来描述气液两相流不能简单地用层流或紊流来描述气液两相流水平或倾斜流动不是轴对称的水平或倾斜流动不是轴对称的由于相界面的存在增加了研究的复杂性由于相界面的存在增加了研究的复杂性总能量方程中应考虑与表面形成
18、的能量问题总能量方程中应考虑与表面形成的能量问题多相流动中各相的温度、组分的浓度都不是均匀的，多相流动中各相的温度、组分的浓度都不是均匀的，相之间有传热和传质相之间有传热和传质各相流速不同，出现滑脱问题，是多相流研究的核心各相流速不同，出现滑脱问题，是多相流研究的核心与重点与重点3.研究气液两相流应注意的几个问题研究气液两相流应注意的几个问题一、基本特征1.经验方法经验方法 从气液两相流动的物理概念出发，或者使用从气液两相流动的物理概念出发，或者使用因次分析法，或者根据流动的基本微分方程式，得到反映因次分析法，或者根据流动的基本微分方程式，得到反映某一特定的两相流动过程的一些无因次参数，然后依
19、据实某一特定的两相流动过程的一些无因次参数，然后依据实验数据整理出描述这一流动过程的经验关系式。验数据整理出描述这一流动过程的经验关系式。优点：优点：使用方便，在一定使用方便，在一定条件下能取得好的结果条件下能取得好的结果缺点：缺点：使用有局限性，且使用有局限性，且很难从其中得出更深层次很难从其中得出更深层次的关系的关系穆迪（Moody）图 二、研究方法2 半经验方法半经验方法 根据所研究的气液两相流动过程的特点，采用根据所研究的气液两相流动过程的特点，采用适当的假设和简化，再从两相流动的基本方程式出发，求得描适当的假设和简化，再从两相流动的基本方程式出发，求得描述这一流动过程的函数关系式，最
20、后用实验方法确定出函数关述这一流动过程的函数关系式，最后用实验方法确定出函数关系式中的经验系数。系式中的经验系数。优点：优点：有一定的理论基础，应用广泛有一定的理论基础，应用广泛缺点：缺点：存在简化和假设，具有不准确性存在简化和假设，具有不准确性二、研究方法3 理论分析方法理论分析方法 针对各种流动过程的特点，应用流体力学针对各种流动过程的特点，应用流体力学方法对其流动特性进行分析，进而建立起描述这一流动过程方法对其流动特性进行分析，进而建立起描述这一流动过程的解析关系式。的解析关系式。优点：优点：以理论分析为基础，可以得到解析关系式以理论分析为基础，可以得到解析关系式缺点：缺点：建立关系式困
21、难，求解复杂建立关系式困难，求解复杂二、研究方法1.1.定义：定义：两相流动中两相介质的分布状况称为流型或两相流动两相流动中两相介质的分布状况称为流型或两相流动结构结构。流型是两相流的一个重要的研究方面流型是两相流的一个重要的研究方面 各相介质的体积比例各相介质的体积比例 介质的流速介质的流速 各相的物理及化学性质各相的物理及化学性质(密度、粘度、界面张力等密度、粘度、界面张力等)流道的几何形状流道的几何形状 壁面特性壁面特性 管道的安装方式管道的安装方式三、流动型态2.影响流型的因素第一类划分方法：根据两相介质分布的外形划分第一类划分方法：根据两相介质分布的外形划分3.流型分类垂直气液两相流
22、流型 水平气液两相流流型三、流动型态第二类划分方法：按流动的数学模型或流体的分散程度划分第二类划分方法：按流动的数学模型或流体的分散程度划分 分散流、间歇流和分离流三种分散流、间歇流和分离流三种 第一类划分方法较为直观第一类划分方法较为直观第二类划分方法便于进行数学处理第二类划分方法便于进行数学处理第一类划分方法第一类划分方法泡状流泡状流弹状流或团状流弹状流或团状流层状流层状流波状流波状流段塞流或冲击流段塞流或冲击流环状流环状流雾状流雾状流第二类划分方法第二类划分方法分散流分散流间歇流间歇流分离流分离流分离流分离流间歇流间歇流分离流分离流分散流分散流三、流动型态4.流型图 描述描述(表示表示)
23、流型变化及其界限的图。流型变化及其界限的图。把流型变换的实验数据把流型变换的实验数据加以总结归纳后，按照两个或多个主要的流动参数绘成曲线，加以总结归纳后，按照两个或多个主要的流动参数绘成曲线，便可以得到流型图。便可以得到流型图。贝克的流型分布图贝克的流型分布图阿齐兹阿齐兹-戈威尔戈威尔-福格拉锡流型分布图福格拉锡流型分布图三、流动型态第三节 气液两相流动的特性参数一、流量质量流量质量流量:单位时间内流过过流断面的流体质量，单位时间内流过过流断面的流体质量，kg/slgGGG体积流量体积流量:单位时间内流过过流断面的流体体积，单位时间内流过过流断面的流体体积，m3/slgQQQ二、流速实际速度实
24、际速度(velocity):单位相面积所通过的该相体积流量单位相面积所通过的该相体积流量,m/s 气相实际速度气相实际速度gggAQv 液相实际速度液相实际速度lllAQv lgAAAgggQG lllQG gAlA折算速度折算速度(superficial velocity):又称表观流速，又称表观流速，假定管道全被一相占据时的流动速度假定管道全被一相占据时的流动速度,m/s气相折算速度气相折算速度AQvgsg液相折算速度液相折算速度AQvlsl两相混合物两相混合物(平均平均)速度速度sgsllgmvvAQQv混合物的质量速度混合物的质量速度AGgggAvQ lllAvQ AvQsggAvQs
25、ll二、流速滑差滑差(滑脱速度滑脱速度)Slip velocity Slippage velocity气液两相实际速度之差气液两相实际速度之差。lgsvvvv滑动滑动(滑移滑移)比比 Slip ratio气相实际速度与液相实际速度之比。气相实际速度与液相实际速度之比。lgvvS 三、滑差和滑动比三、滑差和滑动比质量含气率：质量含气率：单位时间内流过过流断面的混合物总质量单位时间内流过过流断面的混合物总质量G中气相质量所占的份额。中气相质量所占的份额。glggGGGGGxglllGGGGGx1四、含气率与含液率四、含气率与含液率质量含液率：质量含液率：单位时间内流过过流断面的混合物总质量中单位时
26、间内流过过流断面的混合物总质量中液相质量所占的份额。液相质量所占的份额。体积含气率：体积含气率：单位时间流过过流断面两相流体单位时间流过过流断面两相流体(混合物混合物)总总体积体积Q中气相所占的份额。中气相所占的份额。lgggQQQQQ体积含液率：体积含液率：单位时间流过过流断面两相流体单位时间流过过流断面两相流体(混合物混合物)总总体积体积Q中液相所占的份额。中液相所占的份额。lgllQQQQQ1四、含气率与含液率四、含气率与含液率与x间的关系lglggxxxQQQ1gllggGGGx1lllQG gggQG 四、含气率与含液率四、含气率与含液率真实含气率真实含气率(void fractio
27、n):又称又称截面含气率截面含气率或或空隙率空隙率，为任一流动截面内气相面积占总面积的份额为任一流动截面内气相面积占总面积的份额(气相面气相面积与管道总面积之比积与管道总面积之比)。AAg真实含液率真实含液率(liquid holdup):又称又称截面含液率截面含液率或或持液率持液率，为任一流动截面内液相面积占总面积的份额为任一流动截面内液相面积占总面积的份额(液相面液相面积与管道总面积之比积与管道总面积之比)AAHlllgAAAgAlA四、含气率与含液率四、含气率与含液率 ,时当glvvlgglglggllgggglggQQQvvQQQvQvQvQAAA ,glvv一般地，SvvQQQQQQ
28、vvQQQlglgglglgg)1(与与的比较的比较 四、含气率与含液率四、含气率与含液率lgglgglggLLLLALALAVVVlgglgglggAAALALALAVVV快关阀法测量真实含气率：易于实现，只能得到平均值，快关阀法测量真实含气率：易于实现，只能得到平均值，且不能在线测量。且不能在线测量。gVlVgLlL快关阀示意图快关阀示意图真实含气率的简单测量真实含气率的简单测量 Ag Al L 四、含气率与含液率四、含气率与含液率流动密度流动密度:单位时间内流过过流断面的混合物质量与体单位时间内流过过流断面的混合物质量与体 积之比。积之比。QG/lgggllglQQQQGG)1(五、流动
29、密度与真实密度真实密度：真实密度：在流道上取微段，微段内两相流体的质量在流道上取微段，微段内两相流体的质量 与容积之比与容积之比 lgllggllggAAALALALA)1(Ag Al L ,)(则时无滑脱当glvvglvv 一般地，lg)1(lg)1(lg)1(glvv 流动密度与真实密度的比较lg)1(五、流动密度与真实密度第一章小结相的定义相的定义气液两相流流型的两类划分方法气液两相流流型的两类划分方法气液两相流的基本特征气液两相流的基本特征两相流的主要研究方法两相流的主要研究方法流型、流型图流型、流型图流量（质量、体积）流量（质量、体积）流速（实际速度、折算速度、混合物速度）流速（实际
30、速度、折算速度、混合物速度）滑脱速度、滑动比滑脱速度、滑动比含气率与含液率含气率与含液率(质量、体积、真实质量、体积、真实)及其相互关系及其相互关系密度（流动密度与真实密度）密度（流动密度与真实密度）前次课回顾相的概念相的概念 相是体系中具有相同化学组成和物理性质的一相是体系中具有相同化学组成和物理性质的一部分部分 ，与体系的其它均匀部分有界面隔开，与体系的其它均匀部分有界面隔开相相物质形态物质形态(气态、液态、固态气态、液态、固态)例如：水-冰 液固两相流 水-水蒸气 气液两相流 泥浆 液固两相流 油-气 气液两相流 油-气-水 气液液三相流 单工质单工质多工质多工质相相组分组分(指混合物中
31、的各个成分指混合物中的各个成分 )体系中存在相界面体系中存在相界面 两相之间也存在力的作用，出两相之间也存在力的作用，出现质量和能量的交换时伴随着机械能的损失现质量和能量的交换时伴随着机械能的损失 两相的分布情况多种多样两相的分布情况多种多样 两相流动中两相介质的分两相流动中两相介质的分布状况称为流型布状况称为流型 两相流动中存在滑脱现象两相流动中存在滑脱现象 相间速度的差异称为滑脱相间速度的差异称为滑脱，滑脱将产生附加的能量损失滑脱将产生附加的能量损失 沿程流体体积流量有很大变化，质量流量不变沿程流体体积流量有很大变化，质量流量不变气液两相流的基本特征气液两相流的基本特征前次课回顾垂直气液两
32、相流流型 水平气液两相流流型多相流体的主要流型多相流体的主要流型前次课回顾体系中存在相界面体系中存在相界面 两相之间也存在力的作用，出现质两相之间也存在力的作用，出现质量和能量的交换时伴随着机械能的损失量和能量的交换时伴随着机械能的损失 两相的分布情况多种多样两相的分布情况多种多样 两相流动中两相介质的分布两相流动中两相介质的分布状况称为流型状况称为流型 两相流动中存在滑脱现象两相流动中存在滑脱现象 相间速度的差异称为滑脱，相间速度的差异称为滑脱，滑脱将产生附加的能量损失滑脱将产生附加的能量损失 沿程流体体积流量有很大变化，质量流量不变沿程流体体积流量有很大变化，质量流量不变气液两相流的基本特
33、征气液两相流的基本特征课程回顾一、流量质量流量质量流量:单位时间内流过过流断面的流体质量，单位时间内流过过流断面的流体质量，kg/slgGGG体积流量体积流量:单位时间内流过过流断面的流体体积，单位时间内流过过流断面的流体体积，m3/slgQQQ二、流速实际速度实际速度(velocity):单位相面积所通过的该相体积流量单位相面积所通过的该相体积流量,m/s 气相实际速度气相实际速度gggAQv 液相实际速度液相实际速度lllAQv lgAAAgggQG lllQG gAlA折算速度折算速度(superficial velocity):又称表观流速，又称表观流速，假定管道全被一相占据时的流动速
34、度假定管道全被一相占据时的流动速度,m/s气相折算速度气相折算速度AQvgsg液相折算速度液相折算速度AQvlsl两相混合物两相混合物(平均平均)速度速度sgsllgmvvAQQv混合物的质量速度混合物的质量速度AGgggAvQ lllAvQ AvQsggAvQsll滑差滑差(滑脱速度滑脱速度)Slip velocity Slippage velocity气液两相实际速度之差气液两相实际速度之差。lgsvvvv滑动滑动(滑移滑移)比比 Slip ratio气相实际速度与液相实际速度之比。气相实际速度与液相实际速度之比。lgvvS 质量含气率：质量含气率：单位时间内流过过流断面的混合物总质量单位
35、时间内流过过流断面的混合物总质量G中气相质量所占的份额。中气相质量所占的份额。glggGGGGGxglllGGGGGx1质量含液率：质量含液率：单位时间内流过过流断面的混合物总质量中单位时间内流过过流断面的混合物总质量中液相质量所占的份额。液相质量所占的份额。体积含气率：体积含气率：单位时间流过过流断面两相流体单位时间流过过流断面两相流体(混合物混合物)总总体积体积Q中气相所占的份额。中气相所占的份额。lgggQQQQQ体积含液率：体积含液率：单位时间流过过流断面两相流体单位时间流过过流断面两相流体(混合物混合物)总总体积体积Q中液相所占的份额。中液相所占的份额。lgllQQQQQ1与x间的关
36、系lglggxxxQQQ1gllggGGGx1lllQG gggQG 真实含气率真实含气率(void fraction):又称又称截面含气率截面含气率或或空隙率空隙率，为任一流动截面内气相面积占总面积的份额为任一流动截面内气相面积占总面积的份额(气相面气相面积与管道总面积之比积与管道总面积之比)。AAg真实含液率真实含液率(liquid holdup):又称又称截面含液率截面含液率或或持液率持液率，为任一流动截面内液相面积占总面积的份额为任一流动截面内液相面积占总面积的份额(液相面液相面积与管道总面积之比积与管道总面积之比)AAHlllgAAAgAlA ,时当glvvlgglglggllggg
37、glggQQQvvQQQvQvQvQAAA ,glvv一般地，SvvQQQQQQvvQQQlglgglglgg)1(与与的比较的比较 流动密度流动密度:单位时间内流过过流断面的混合物质量与体单位时间内流过过流断面的混合物质量与体 积之比。积之比。QG/lgggllglQQQQGG)1(真实密度：真实密度：在流道上取微段，微段内两相流体的质量在流道上取微段，微段内两相流体的质量 与容积之比与容积之比 lgllggllggAAALALALA)1(,)(则时无滑脱当glvvglvv 一般地，lg)1(lg)1(lg)1(glvv 流动密度与真实密度的比较lg)1(第二章 气液两相流的模型 气液两相流
38、动的规律较单相流复杂得多气液两相流动的规律较单相流复杂得多,常采用简常采用简化的流动模型进行处理，以便探讨其流动规律。常用化的流动模型进行处理，以便探讨其流动规律。常用的模型有的模型有均相流动模型均相流动模型、分相流动模型分相流动模型和和漂移流动模漂移流动模型、流动型态模型型、流动型态模型等。等。第一节 均相流动模型（1）气相和液相的实际速度相等，即）气相和液相的实际速度相等，即（2）两相介质已达到热力学平衡状态，压力、）两相介质已达到热力学平衡状态，压力、密度互为单值函数。密度互为单值函数。vvvgl对于泡状流和雾状流，具有较高的精确性对于泡状流和雾状流，具有较高的精确性对于弹状流和段塞流，
39、需要进行时间平均修正对于弹状流和段塞流，需要进行时间平均修正对于层状流、波状流和环状流，则误差较大对于层状流、波状流和环状流，则误差较大两个假定两个假定特点特点定义：把气液两相混合物看成均匀介质，其物定义：把气液两相混合物看成均匀介质，其物 性参数取两相的均值而建立的模型性参数取两相的均值而建立的模型两相流基本方程式以单相流基本方程式为基础。单两相流基本方程式以单相流基本方程式为基础。单相流基本方程式理论上总结为三个基本方程式：相流基本方程式理论上总结为三个基本方程式：连续性方程连续性方程动量方程动量方程能量方程能量方程质量守恒质量守恒动量守恒动量守恒能量守恒能量守恒一、均流模型的基本方程式一
40、、均流模型的基本方程式1连续方程式连续方程式根据质量守恒定律根据质量守恒定律常数vAG 2动量方程式动量方程式取一维流段来研究，取一维流段来研究，根据动量根据动量定理，可得动量方程式：定理，可得动量方程式：稳定的一维均相流动GdvgAdzdFAdpsin2-23能量方程式能量方程式根据机械能守恒定律，有根据机械能守恒定律，有02sin2dEvpgzd密度密度可用两相混合物的比容可用两相混合物的比容 表示为表示为vv/1所以所以vppdpvvpdvpd)(则能量方程式为则能量方程式为0)2(sin2dEvddpvvpdgdzdE 单位质量两相流体单位质量两相流体的机械能损失的机械能损失2-32-
41、5压差的表达式压差的表达式dEvpdvdgdzdpv)2(sin2二、均流模型的压力梯度微分方程式在动量方程式中在动量方程式中GdvgAdzdFAdpsinDdzdFww流体与管壁的剪切应力流体与管壁的剪切应力2-72-6其中fvw221代入代入2-7得得DdzvfDdzdFw22 f 范宁摩阻系数范宁摩阻系数4212vw 摩擦阻力系数摩擦阻力系数穆迪（穆迪（Moody）图）图 Re16 Re64 f层流区层流区4f221v单位体积流体的动能 1lglllggvvxvGGvGvGQvvAGAGAGAQvvvlg1所以22lglvvxvdAGvdAGvAGGdGdvGdvgAdzdFAdpsin
42、Gdv2-132-11由于假定两相流动已达到热力学平衡状态由于假定两相流动已达到热力学平衡状态dpdpdvxdpdpdvxdxvvvvxvdvdlglglgl1因因0dpdvldpdpdvxdxvvvdglg则则所以所以dpdpdvxdxvvAGGdvglg2)(pfv 2-15将将dF、Gdv的表达式代入动量方程式，的表达式代入动量方程式，得得dpdpdvxdxvvAGvgAdzDdzvfAdpglg22sin2两边同除两边同除 得得AdzdzdpdpdvxdzdxvvAGvgDvfdzdpglg222sin42整理可得整理可得dpdvxAGdzdxvvAGvvxvgvvxvAGDfdzdp
43、glglgllgl2221sin2因许多参数因许多参数f 、沿程变化，无法用解析法积分，沿程变化，无法用解析法积分，得用差分法分段计算得用差分法分段计算压力梯度微分方程式lgvvdpdvg三、三、气液两相流能量平衡方程建立气液两相流能量平衡方程建立1能量平衡方程推导能量平衡方程推导根据能量守衡定律写出两个根据能量守衡定律写出两个流动断面间的能量平衡关系流动断面间的能量平衡关系量从断面流出的流体能的能量在断面和之间耗失体额外做的功在断面和之间对流进入断面的流体能量倾斜管流能量平衡关系示意图2222221121112sin2sinVpmvmgZUqVpmvmgZU此式除内能U外，其它参数可测得动能
44、22mv势能mgh膨胀能pV内能 U倾斜管流能量平衡关系示意图于是得能量平衡方程式差分形式于是得能量平衡方程式差分形式0)()2()sin(2qpVmvmgZU将所得能量平衡方程式写成微分形式：将所得能量平衡方程式写成微分形式：0)(sindqpVddZmgmvdvdU2222221121112sin2sinVpmvmgZUqVpmvmgZU内能内能U虽然不能直接测量和计算其绝对值，但可求得两种状虽然不能直接测量和计算其绝对值，但可求得两种状态下的相对变化。根据热力学第一定律，对于可逆过程：态下的相对变化。根据热力学第一定律，对于可逆过程：pdVdUdqpdVdqdU或或而对于不可逆过程来讲：
45、而对于不可逆过程来讲：pdVdUdqdqr以以 表示摩擦消耗的功，则表示摩擦消耗的功，则:wdlwrdldq wdlpdVdqdU 摩擦产生的热量 摩擦消耗的功rdqwdl0sinwdldZmgmvdvVdp代入代入 并整理得并整理得dU0)(sindqpVddZmgmvdvdU0sindZIdgdZdvvdZdpw两边同除以两边同除以VdZsindZIdgdZdvvdZdpw总压力降总压力降动能损失动能损失重位损失重位损失摩擦损失摩擦损失sin)(gdZdp重位dZdIdZdpw摩擦)(dZdvvdZdp加速度)(加速度摩擦重位)()()(dZdpdZdpdZdpdZdp2)(2vddzdp
46、摩擦dZdp2sin2vddZdvvgdZdp由并取为正值0,0重位dZdp22vddxdvvdxdp 9022vddhdvvgdhdp对于水平管流对于垂直管流以h表示高度x表示流向坐标讨论讨论1sin2sin2mmmmmmmvddZdvvgdZdp对于多相混合物通用表达式，求解的关键是确定 、及 mmvm2部分相关参数的计算部分相关参数的计算（1 1）两相介质的平均密度）两相介质的平均密度在均流模型中，两相介质的密在均流模型中，两相介质的密度取气液两相密度的平均值，度取气液两相密度的平均值，而求其平均值的主要以下两种而求其平均值的主要以下两种lg)1(lg)1(按体积含气率按体积含气率计算计
47、算按空隙率计算按空隙率计算流动密度流动密度(无滑脱密度无滑脱密度)真实密度真实密度(有滑脱密度有滑脱密度)滑脱损失实际上属于重力损失（2 2）两相介质的平均粘度）两相介质的平均粘度在均流模型中，两相介质的粘度在均流模型中，两相介质的粘度是气液两相粘度的平均值，而求是气液两相粘度的平均值，而求其平均值的方法很多，常用的有其平均值的方法很多，常用的有以下几类以下几类按质量含气率按质量含气率计算计算西克奇蒂公式西克奇蒂公式麦克亚当斯公麦克亚当斯公式式杜克勒公式杜克勒公式戴维森公式戴维森公式lgmxx1lgmxx11llggmxvvxvv111lglmvvx各有特色和适用范围按空隙率计算按空隙率计算l
48、gm1按体积含气率按体积含气率计算计算lgm11lgm（3 3）均流模型摩擦阻力折算系数）均流模型摩擦阻力折算系数按均流模型进行气液两相流动摩阻压差计算时，常把按均流模型进行气液两相流动摩阻压差计算时，常把两相流动两相流动摩擦阻力的计算与单相流动摩擦阻力的计算关联起来摩擦阻力的计算与单相流动摩擦阻力的计算关联起来，即常使，即常使用用全液相折算系数全液相折算系数、分液相折算系数分液相折算系数或或分气相折算系数分气相折算系数。全液相折算系数全液相折算系数水平管道内的两相流动水平管道内的两相流动,均匀流动，均匀流动，管管径为径为D，截面积为，截面积为A，流段长度为，流段长度为dz。速度速度v沿流程不
49、变，质量流量为沿流程不变，质量流量为G，此，此时，没有重位压差与加速度。时，没有重位压差与加速度。DdzvfDdzdFw22且两相流动dpAdFDdzvfAdpdFl2 20000200220020020 2 2vfvfDdzvfDdzvfdpdpdFdFll定义全液相折算全液相折算系数系数020dpdp单相流的摩擦力单相流的摩擦力全液相流动全液相流动设管道的设管道的D、A和和dz不变，通过管道的质量流不变，通过管道的质量流量仍为量仍为G，但流体为单一的液体，没有气相，但流体为单一的液体，没有气相，此时的流体密度为此时的流体密度为 ，速度为，速度为 l0v分液相折算系数管道的管道的D、A和和d
50、z不变，通过管道的流体也不变，通过管道的流体也为单一的液体，但其质量流量等于两相流为单一的液体，但其质量流量等于两相流动中液相的质量流量。此时的流体密度动中液相的质量流量。此时的流体密度为为 ，速度为，速度为 。lslvDdzvfAdpdFsllslslsl2 222222 2 2sllslsllslslsllvfvfDdzvfDdzvfdpdpdFdF分液相折算分液相折算系数系数slldpdp2单相液流的摩擦力单相液流的摩擦力分液相流动分液相流动分气相折算系数DdzvfAdpdFsglsgsgsg2 222222 2 2sggsgsggsgsgsggvfvfDdzvfDdzvfdpdpdFd
