4. 井筒多相流基础与关联式选择
油井从井底到井口的流动，几乎不可能保持单相。当井底流压低于泡点压力时，溶解气开始析出，井筒内同时存在油、气，甚至水，形成气液两相或油气水三相流动。多相流与单相流最本质的差别在于：气液两相之间存在滑脱效应，气相密度小、流速快，液相被"拖在后面"，导致相同质量流量下混合物的实际密度和摩阻都随流型剧烈变化。因此，井筒压力分布计算的核心，是判断当前处于什么流型，再选择能描述该流型的压降关联式。
一、多相流的流型划分
在垂直井筒中，随井深自下而上、压力逐渐降低、气体不断膨胀，流型通常按以下顺序演变：
- 泡状流（Bubble Flow）：气相以分散小气泡形式存在于连续液相中，气泡近似均匀上浮。多见于高流压、低气液比井段，滑脱小，混合密度接近液相密度。
- 段塞流（Slug Flow）：气泡聚并成与管径同量级的泰勒气泡，气液段交替上升。这是自喷井筒中最常见的流型，压力波动明显，滑脱显著。
- 搅动流（Churn Flow）：气速进一步增大，段塞结构被破坏，气液剧烈掺混、上下振荡，是段塞流向环状流过渡的不稳定区。
- 环状流（Annular Flow）：气相在管中心形成高速气芯，液相以液膜形式沿管壁上升，并夹带液滴。多见于高气液比、井口附近井段。
流型不是固定属性，而是随井深连续变化的。同一口井从井底到井口可能依次经历泡状—段塞—搅动—环状流，因此严格来说压降计算应分段进行，每段按当地流型选取关联式。工程上常按井段长度或压力步长离散，逐段迭代。
二、常用压降关联式的特点
井筒压力梯度由重力项、摩阻项和加速度项组成，多相流关联式的差异主要体现在如何计算持液率和摩阻系数。工程上最常用的两类是 Hagedorn-Brown 与 Beggs-Brill。
Hagedorn-Brown 关联式是垂直管多相流的经典经验式，基于大量现场与实验数据拟合，通过四个无量纲数（液相速度数、气相速度数、管径数、液相黏度数）计算持液率，再结合摩阻压降。它对垂直井、中低气液比、泡状流至段塞流段表现稳定，是许多商业软件垂直井段的默认选项之一。其局限在于对高气液比环状流和水平/大斜度井的适应性较差。
Beggs-Brill 关联式的突出优点是引入管斜角修正，先按无滑脱条件判别流型（分层、间歇、分散），再用持液率与摩阻系数随倾角变化的经验关系修正，因此可统一处理垂直、倾斜和水平管段。对于定向井、水平井或井斜变化大的井，Beggs-Brill 通常比纯垂直关联式更合适。但在接近垂直、高气液比的环状流条件下，其精度可能不如专门的环状流模型。
三、关联式选择依据
选择关联式不能只看"哪个名气大"，而应结合井型、流型、气液比和可用数据综合判断。可参考以下决策逻辑：
| 条件 | 推荐关联式 | 理由 |
|---|---|---|
| 直井、中低气液比、泡状—段塞流 | Hagedorn-Brown | 垂直管经验数据充分，持液率计算稳定 |
| 定向井/水平井、井斜变化大 | Beggs-Brill | 含倾角修正，可统一处理各井段 |
| 高气液比、井口环状流段 | Gray 或 Duns-Ros 对比校核 | 专门针对高气速环状流，减少持液率低估 |
| 凝析气井、低液量 | Gray 关联式 | 针对气井少量凝析液工况拟合 |
| 无实测数据、初步估算 | 两种关联式并行计算对比 | 用结果离散度评估不确定性 |
四、一个分段计算的示例
假设某直井井底流压 28 MPa，井口回压 2 MPa，产液量 60 m³/d，气液比 350 m³/m³，井深 3000 m。按 0.5 MPa 压力步长自井底向井口迭代：
- 以井底压力为起点，计算当地气液物性（溶解气油比、体积系数、黏度、界面张力）。
- 用 Beggs-Brill 流型判别图判断该段流型；若为泡状或段塞流，改用 Hagedorn-Brown 计算持液率与摩阻。
- 计算该段压力梯度，累加得到下一节点压力，重复直至井口。
- 将计算井口压力与回压 2 MPa 对比，若偏差超过 5%，调整产量或关联式重新迭代。
示例中若井底段判为泡状流、中段为段塞流、井口段为环状流，则应在不同井段切换关联式，而不是全井统一使用一个公式。这种"分段判流型、分段选式"的做法，正是节点分析中井筒压力计算的基本功。
掌握流型演变规律与关联式适用边界后，井筒压力计算才能为后续节点分析提供可靠输入。实际工作中，建议把"流型判别—关联式选择—分段迭代—实测校核"作为固定检查链条，避免因关联式误选导致产能预测出现系统性偏差。