📚 课程目录第3章多相流管流与压力温度分布

3. 多相流管流与压力温度分布

电潜泵（ESP）从井下吸入口到地面出油管线，整段举升通道内流动的几乎都是油气水多相混合物。泵的吸入压力、排出压力、气体体积分数、黏度和流体密度，都取决于这段管流中的持液率、压降与温度分布。因此，多相流管流计算不是选泵之后的校核工作，而是确定泵挂深度、泵型、级数和运行频率的前置输入。本节按“持液率—压降—温度剖面—对泵工况的影响”这条物理链条展开，并给出可操作的示例计算思路。

一、持液率：多相管流的核心状态量

持液率 H_L 定义为管道某一截面上液相所占面积与总流通面积之比，即 H_L = A_L/A。与之对应，截面含气率 α_g = 1 − H_L。需要区分两个容易混淆的量：体积含液率（入口液相体积流量与总流量之比）和截面持液率。由于气液两相存在滑脱——气体密度小、上升快，液相被“滞留”——在垂直上升管中截面持液率通常明显大于体积含液率。滑脱越强，持液率越高，管内混合密度越大，静压梯度也越大。

影响持液率的因素包括：气液比、液相流量、管径、流体黏度、管斜角以及压力温度（通过改变气液物性）。工程上常用 Beggs-Brill、Hagedorn-Brown、Orkiszewski 等经验相关式计算。以 Beggs-Brill 为例，其思路是先按无滑脱计算，再用倾角修正和滑脱修正因子修正持液率，因此对水平、倾斜和垂直段可统一处理，适合井筒—地面管线一体化计算。

二、压降计算：摩擦、重力与加速度三项

多相管流的总压降由三部分组成：重力压降（举升液柱）、摩擦压降和加速度压降。对垂直井筒，重力项往往占主导；对长距离水平出油管线，摩擦项占主导。总压梯度可写为：

dp/dz = ρ_m·g·sinθ + f_m·ρ_m·v_m²/(2D) + ρ_m·v_m·dv_m/dz
其中：ρ_m = ρ_L·H_L + ρ_g·(1−H_L)   （混合密度）
      v_m = (q_L + q_g)/A             （混合流速）
      f_m 由混合雷诺数与持液率相关式确定

由于 ρ_m、H_L、f_m 都随压力和温度变化，压降方程必须分段迭代求解：把井筒沿深度分成若干段（如每 50～100 m 一段），从已知的井口或井底压力出发，逐段假设压降、更新物性、重算持液率，直到收敛。这正是多相流管流计算比单相水力计算复杂得多的原因。

示例（合理假设，非实测）：某井产液 120 m³/d、气液比 80 m³/m³、井口压力 1.2 MPa、井深 1800 m、油管 62 mm。分段迭代后若算得井底流压约 9～10 MPa 量级，则泵吸入口压力需在此基础上扣除沉没段压降，用于后续气体影响判断。实际数值必须以本井 PVT 与管流软件结果为准。

三、温度剖面模型

井筒温度并非简单地等于地层温度。流体从地层流向井筒时，在井底附近温度接近原始地层温度；沿井筒上升过程中，流体向周围地层散热，同时受到地温梯度控制。常用模型把井筒传热视为一维稳态径向传热：以地温梯度给出远处地层温度 T_e(z)，用总传热系数 U 描述油管—环空—套管—水泥环—地层的热阻，得到：

dT_f/dz = (T_e(z) − T_f) / L_R  +  焦耳-汤姆逊项 + 摩擦生热项
L_R 为松弛距离，取决于质量流量、比热和总传热系数 U

对 ESP 井，还需特别关注两点：一是电机与电缆发热会抬高泵段附近流体温度，电机效率损失最终转化为热量进入井液；二是焦耳-汤姆逊效应在气体膨胀或节流处引起温度变化。温度剖面直接影响液相黏度、气体溶解度和泵内流体物性，因此温度模型必须与压降模型耦合迭代，而不是单独计算。

四、这些分布如何落到泵工况上

把上述结果代入泵的选型与运行，可以归纳为四条直接影响。

第一，吸入压力与沉没度。泵吸入口压力等于该点流压，由井底流压减去沉没段的重力与摩擦压降得到。若持液率被低估，混合密度偏小，算出的吸入压力偏高，实际沉没度不足，泵容易因供液不足而气蚀。工程上通常要求沉没度留有余量，具体数值按区块经验与气体含量确定。

第二，游离气与气锁。泵吸入口处的游离气体积分数由该点压力、温度和溶解气量决定。压力越低、温度越高，游离气越多。当游离气体积分数超过泵型允许范围（不同泵型差异较大，常见在 10%～30% 量级），需加装气体分离器或加深泵挂。这一步判断完全依赖管流计算给出的吸入点压力与持液率。

第三，混合密度与所需扬程。泵需要提供的总扬程要克服井筒混合液柱、地面管线压降和井口回压。混合密度 ρ_m 由持液率决定，持液率算错，扬程需求就错，级数随之偏差。实际做法是把压降剖面按深度积分，得到泵入口到井口的所需总压差，再换算为扬程。

第四，黏度与泵效。温度剖面决定液相黏度。温度偏低时黏度上升，泵的容积效率和扬程都会下降，同时摩擦压降增大。对高含蜡或高黏油田，温度剖面还是判断是否需要伴热或化学降黏的依据。

五、可执行的计算组织方式

实际工作中，建议按以下顺序组织一次完整的管流—泵工况联算，而不是把各模型割裂使用：

准备输入：井身结构、油管/套管尺寸、产液量、气液比、含水率、PVT 数据、地温梯度、井口压力与温度。
沿井深分段（建议 50～100 m），从井口向下或从井底向上逐段迭代，每段更新压力、温度与流体物性。
每段用多相流相关式计算持液率与压降，同时用温度模型更新该段温度，二者交替迭代至收敛。
提取泵吸入口处的压力、温度、游离气体积分数和混合密度，作为泵选型与气体处理的判据。
把压降剖面积分得到所需总扬程，结合泵型特性曲线确定级数与运行频率，并复核电机负载与电缆压降。

常见误区：用单相液柱公式估算吸入压力，或把井口产量直接当作泵排量而不考虑泵内气体和滑脱。前者会高估沉没度，后者会低估泵的实际负载。多相流条件下，任何简化都必须说明适用条件，并在结果上留出安全余量。

综上，持液率决定混合密度，混合密度与流速决定压降，压降与传热共同决定温度剖面，而温度又反过来改变物性和持液率。这条耦合链的终点，就是泵的吸入条件、扬程需求和运行稳定性。把这条链算清楚，ESP 的泵挂深度、泵型、级数和频率才有可靠依据。