3. 井筒压力场控制方程与流体力学基础

井筒压力场是多相流计算的“骨架”：温度剖面决定流体物性，而压力剖面决定相态、密度、溶解气量乃至举升效率。要可靠地算出沿井深的压力分布，必须先把流体力学的基本守恒关系落到井筒这一具体对象上，明确压降由哪几部分组成、每一项在什么条件下占主导，再据此组织建模思路。本节从守恒方程出发，逐项拆解压降组成，并给出可操作的建模流程。

一、坐标系与基本假设

先约定参照系，全文与图示保持一致：取井口为起点、沿井眼轴线向下为 s 正方向，垂深为 z（向下为正），井斜角为 θ（与垂线夹角）。于是沿程微元 ds 对应的垂深增量为 dz = ds·cosθ。压力 p 沿流动方向的变化率写作 dp/ds，其值为负表示压力随下行而升高（即压降为正）。工程上常把“压降梯度”定义为 −dp/ds，本节沿用这一约定。

建模通常采用一维稳态两流体或均相假设：流体沿轴向流动，同一截面上各相压力相等、温度相同；忽略径向变化与瞬态项；井筒与地层间存在热交换，但压力方程本身只通过物性（密度、黏度、界面张力）与温度耦合。这些假设在常规生产井的稳定生产段是合理的，但在开井/关井瞬态、气侵等工况下需要引入瞬态项，属于后续扩展内容。

二、质量守恒与动量守恒

对井筒微元控制体，质量守恒要求流入与流出的质量流量之差等于控制体内质量累积。稳态下累积项为零，对气相、液相分别写出：

d(ρ_g · v_sg · A)/ds = 0      （气相，无相间传质时）
d(ρ_l · v_sl · A)/ds = 0      （液相）
Σ(ρ_k · v_sk) = G_m = const   （总质量流量守恒）

其中 v_sg、v_sl 为表观速度，A 为流通截面积。若存在溶解气析出或凝析，则相间传质项不为零，需在方程右端加入源项，这也是黑油模型与组分模型的分界点。

动量守恒描述压力梯度与各项力之间的平衡。对微元段沿轴向列动量方程，整理后得到工程上最常用的压力梯度表达式：

dp/ds = −ρ_m · g · cosθ − (f · ρ_m · v_m²) / (2 · D) − ρ_m · v_m · (dv_m/ds)

式中 ρ_m 为混合密度，v_m 为混合流速，f 为范宁摩阻系数，D 为管径。三项依次对应重力、摩擦与加速度，下面逐项讨论其物理含义与量级特征。

三、三项压降的物理含义与量级判断

重力项是举升流体所需克服的静液柱压力，与混合密度和井斜角余弦成正比。垂直井中它通常占总压降的绝大部分；水平井段 cosθ≈0，重力项趋近于零，压力主要由摩擦控制。这也解释了为什么水平井的井底压力对流量变化更敏感。

摩擦项与流速平方成正比、与管径成反比。多相流中摩阻系数 f 不再是单相 Moody 图上的常数，而依赖流型（泡状、段塞、环状等）与持液率。工程上常先用 Beggs-Brill、Hagedorn-Brown 等经验相关式估算持液率与摩阻倍数，再回代计算。若忽略流型修正，高气液比井的摩擦压降可能被低估数倍。

加速度项来源于流体沿程膨胀（气相为主）导致的动能变化。在井口附近压力低、气体膨胀剧烈，dv_m/ds 较大，该项可达总压降的 5%~15%；而在井底高压段通常很小。对高产气井，忽略加速度项会系统性高估井底压力。

四、井筒压力场建模思路

压力场求解本质上是沿程积分一个常微分方程，但方程右端依赖当地物性与持液率，而物性又依赖当地压力与温度，因此必须迭代。可操作的流程如下：

1. 离散井筒：沿井眼轴线划分 N 个微元段，段长 10~50 m；井口、井底、射孔段、变径处、气举阀位置应设节点。
2. 确定边界条件：通常已知井口压力与温度（或井底压力），另一端的温度由井筒传热模型给出，形成 p-T 耦合的初值。
3. 物性计算：由当地 p、T 计算气相偏差因子 Z、黏度、液相体积系数与溶解气油比，得到 ρ_g、ρ_l 与界面张力。
4. 持液率与流型：用所选多相流相关式计算持液率 H_L，进而得到 ρ_m = H_Lρ_l + (1−H_L)ρ_g 与 v_m。
5. 三项压降求和：按上式算出该段 Δp，更新段末压力，判断收敛（|Δp_新−Δp_旧| < ε，ε 取 0.01~0.1 kPa）。
6. 逐段推进：以段末状态作为下一段入口，重复直至井底，得到完整 p(s) 剖面。

以一个教学示例说明量级（以下为合理假设的示例参数，非实测）：某垂直油井，井深 2500 m，油管内径 62 mm，产液量 80 m³/d，气液比 120 m³/m³，井口压力 1.2 MPa。取中段某 50 m 微元，估算 ρ_m≈620 kg/m³、v_m≈1.1 m/s、f≈0.028，则重力项约 0.30 MPa、摩擦项约 0.012 MPa、加速度项约 0.002 MPa，重力占绝对主导。若把该井改为水平段，cosθ≈0，同样流量下摩擦项成为主要贡献，压降梯度明显下降但沿程温降与相态变化更复杂。

把上述关系串起来看，井筒压力场建模的核心不是某一个公式，而是“守恒方程定框架、三项压降定贡献、物性与持液率定系数、分段迭代定结果”的闭环。掌握这一闭环，后续引入温度耦合、瞬态项与组分模型时，只需在相应位置扩展，而不必推翻整体结构。