4. 多相流基本概念与流型判别
油井井筒中，原油、天然气和水同时向上流动，这种气液（或气液固）共存并具有相间界面的流动称为多相流。井筒温压剖面计算之所以比单相流复杂，根本原因就在于多相流中气液两相的体积占比随压力、温度变化而剧烈改变，且两相之间存在速度差。本节围绕相分数、表观速度、滑脱效应三个基本参数展开，并给出垂直井筒四种常见流型的判别思路。
一、相分数：描述“各相占多少”
相分数是描述多相流组成最基础的参数，分为体积相分数和质量相分数两类。工程上最常用的是持液率（liquid holdup，符号 HL）和持气率（HG），定义为某一段井筒控制体内液相或气相所占的体积份额：
H_L = V_L / (V_L + V_G) H_G = 1 - H_L与之对应的是入口体积含液率（无滑脱持液率）λL = qL/(qL+qG)，它只由地面产量折算到井筒条件下的体积流量决定，不考虑两相速度差。真实持液率 HL 与 λL 的差异，正是滑脱效应存在的直接证据。持液率直接影响混合流体密度 ρm = HLρL + HGρG，进而决定重力压降项，因此它是井筒压力计算中最关键的中间变量之一。
二、表观速度：描述“各相流得多快”
表观速度（superficial velocity）是指假定某一相单独充满整个过流断面时的等效流速，用该相在井筒条件下的体积流量除以管子总截面积：
v_sL = q_L / A v_sG = q_G / A
v_m = v_sL + v_sG （混合表观速度）表观速度不是真实相速度。真实液相速度 vL = vsL/HL，真实气相速度 vG = vsG/HG。由于气液两相在井筒中向上流动时气相通常更快，vG > vL，二者之差即为滑脱速度 vs = vG − vL。表观速度的实用价值在于：它把复杂的多相流动折算成两个可加的“单相等效流量”，是流型判别图版（如 Beggs-Brill、Duns-Ros、Taitel-Dukler）最常用的横纵坐标输入。
三、滑脱效应：两相速度差的来源与后果
滑脱效应指气相因密度小、浮升作用强而比液相上升更快，造成就地持液率高于无滑脱持液率的现象。其后果有三：一是混合密度增大，重力压降上升；二是气相实际停留时间缩短，井筒内气量分布被重新分配；三是摩阻计算需用真实相速度而非表观速度。滑脱程度通常用滑脱速度或持液率比 HL/λL 衡量，在泡状流和段塞流中滑脱最明显，在环状流高速气芯中滑脱反而相对减弱。
四、流型判别方法与实操思路
垂直井筒流型判别通常采用“图版法+机理判据”相结合的方式。以常用的 Beggs-Brill 流型图为例，横坐标为入口体积含液率 λL，纵坐标为混合表观速度 vm 与弗劳德数 Fr 的组合，划分出分层、过渡、间歇、分散四区；映射到垂直井筒，则对应泡状、段塞、搅动、环状四类。实际计算时可按以下顺序操作：
- 取控制体：沿井筒自下而上分段（如每 50 m 一段），逐段取当地压力 p、温度 T。
- 折算体积流量：用气体状态方程把地面气量折算到当地 p、T，得到 qG；油水体积用体积系数 Bo、Bw 折算。
- 算表观速度：由 qL、qG 和油管截面积 A 求 vsL、vsG、vm。
- 算判别参数：求 λL、Fr = vm/√(gD)、以及气相体积分数。
- 查图或代入判据：确定该段流型，再选用对应的持液率与摩阻相关式。
需要强调的是，同一口井从井底到井口往往经历多种流型：井底压力高、气体被压缩，常呈泡状流或段塞流；向上压力降低、气体膨胀，表观气速迅速增大，逐渐过渡为搅动流乃至环状流。因此流型判别必须逐段进行，不能用一个全井统一流型。这也是井筒温压剖面计算需要分段迭代的根本原因。
五、示例：一段井筒的流型判别
设某垂直井段油管内径 62 mm（A ≈ 3.02×10⁻³ m²），当地条件下 qL = 0.0012 m³/s，qG = 0.018 m³/s。则 vsL ≈ 0.40 m/s，vsG ≈ 5.96 m/s，vm ≈ 6.36 m/s，λL ≈ 0.0625。取 g = 9.81、D = 0.062 m，Fr ≈ 6.36/√0.608 ≈ 8.2。该点落在高气速、低含液率区域，判为环状流，应选用环状流的液膜摩阻与持液率关系式计算压降。若同一井段压力升高使 qG 降为 0.003 m³/s，则 vsG ≈ 1.0 m/s，vm ≈ 1.4 m/s，Fr ≈ 1.8，流型转为段塞流，持液率显著上升，重力压降项随之增大。上述数值为教学示例假设，实际应用需以现场 PVT 与产量数据为准。
掌握相分数、表观速度与滑脱效应这三个基本量，并坚持逐段判别流型，才能为后续持液率计算、摩阻压降计算和温压耦合迭代打下可靠基础。