3. 氮气物性计算

连续油管氮气举升排液的设计计算，起点不是井筒多相流模型，而是氮气本身的物性。注气压力从几兆帕升到二十几兆帕、井筒温度从地面十几摄氏度变化到井底上百摄氏度，氮气的密度可以相差数倍，压缩因子偏离 1 的程度也随工况显著变化。若直接套用理想气体状态方程，在 15 MPa、360 K 条件下密度会被低估约 8%～12%，进而使环空静气柱压力、临界携液流量和注气量全部产生系统性偏差。因此本节的任务是：明确氮气密度、压缩因子 Z、黏度 μ 三个参数的计算路径，并建立一张可被后续章节直接调用的物性参数表。

一、三个物性参数的计算路径

压缩因子 Z 是核心中间量。工程上常用 Standing-Katz 图版对应的 Dranchuk-Abou-Kassem（DAK）相关式，以拟对比压力 p_pr=p/p_c、拟对比温度 T_pr=T/T_c 为自变量迭代求解。氮气的临界参数取 p_c=3.394 MPa、T_c=126.2 K，则 15 MPa、360 K 对应 p_pr≈4.42、T_pr≈2.85，属于中高压、中高温区，Z 约在 1.05～1.10 之间。低压高温时 Z→1，高压低温时 Z 明显大于 1，这是氮气在举升工况下的典型特征。

密度 ρ 由真实气体状态方程给出：ρ = pM/(ZRT)，其中 M=28.0134 g/mol，R=8.314 J/(mol·K)。可见密度对 p 近似线性、对 Z 与 T 成反比。计算顺序必须是先求 Z，再求 ρ，不能颠倒。

黏度 μ 采用“低压黏度 + 高压修正”的两步法：先按 Sutherland 或简化的低压氮气黏度关联式求 μ₀(T)，再用 Jossi-Stiel-Thodos 或 Lucas 方法以对比密度 ρ_r=ρ/ρ_c 修正，得到实际黏度。氮气黏度对压力不敏感、对温度较敏感，360 K 时约 19～20 μPa·s，高压下可增至 22～24 μPa·s。黏度主要用于雷诺数、摩阻系数与携液液滴沉降速度计算，精度要求低于密度。

二、物性参数表的建立方法

实际做法不是逐点手算，而是先确定工况网格，再批量计算。以某示例井为例（以下数据为教学假设，非实测）：井口注气压力 8 MPa、温度 293 K；井底压力 18 MPa、温度 373 K。取压力网格 8、10、12、14、16、18 MPa，温度网格 293、313、333、353、373 K，共 30 个节点。对每个节点依次执行：换算绝压与开氏温度 → 求 p_pr、T_pr → DAK 迭代求 Z → 求 ρ → 求 μ₀ 并修正得 μ。计算可用表格公式或脚本实现，核心是保证迭代收敛判据一致（如 |ΔZ|<1×10⁻⁵）。

下表给出该示例网格中部分节点的计算结果，用于说明数量级与变化趋势（示例值，供教学参考）：

从表中可读出三条工程判断：其一，压力由 8 MPa 升至 18 MPa，密度增大约 1.9 倍，静气柱压力贡献不可忽略；其二，Z 在 1.03～1.09 区间，用理想气体替代会带来 3%～9% 的密度误差，在深井中累积为可观的注气压力偏差；其三，黏度变化平缓，对摩阻的影响弱于密度。因此后续计算中密度必须用真实气体式，黏度可用表值插值。

三、计算中的常见偏差与校核

第一类偏差来自压力基准。现场压力表读数为表压，若直接代入状态方程，8 MPa 表压实际约 8.1 MPa 绝压，低压段误差约 1%，高压段虽相对误差减小，但 Z 的拟对比压力偏移会放大密度误差。第二类偏差来自温度基准，井筒温度剖面若用线性假设，在注气量较大时误差可达 10 K 以上，密度误差约 3%。第三类偏差来自迭代初值，DAK 在 p_pr>5 时收敛变慢，初值取 Z=1 可能不收敛，建议初值取 1.05 并限制迭代次数。

校核方法有三条：一是与 NIST 或专业物性软件在 3～5 个节点上对比，密度相对偏差应小于 1%；二是检查 Z 的单调性，等温线上 Z 随压力应单调上升（本工况区间内成立），若出现回落说明迭代发散；三是用理想气体密度做量级对照，偏差应在 3%～12% 的合理区间，若超过 20% 则需检查单位与临界参数。

完成本节后，读者应能独立完成“给定 p、T 网格 → 求 Z → 求 ρ → 求 μ → 形成参数表”的完整流程，并清楚每个参数的误差来源与校核手段。这张表将作为后续静气柱压力计算、临界携液流量判断和注气量设计的统一数据基础。