3. 固定床压降基本方程

在催化剂装填与床层分布分析中，压降是最直接、最敏感的工程指标之一。装填不均、颗粒破碎、粉尘夹带或局部架桥，往往首先表现为床层压降异常升高或分布不均。要判断这些现象、估算风机或压缩机负荷，就必须掌握固定床压降的基本计算工具——Ergun方程。本节围绕它的结构、适用条件、各参数物理意义以及工程上常用的简化形式展开，帮助读者从"会套公式"走向"会判断公式是否适用、结果是否可信"。

一、Ergun方程的结构：黏性项与惯性项的叠加

Ergun方程描述流体通过颗粒填充床时的总压降，其经典形式为：

ΔP/L = 150·μ·(1-ε)²·u / (ε³·d_p²) + 1.75·ρ·(1-ε)·u² / (ε³·d_p)

式中 ΔP 为床层压降（Pa），L 为床层高度（m），μ 为流体动力黏度（Pa·s），ρ 为流体密度（kg/m³），u 为按空床截面积计算的表观流速（m/s），ε 为床层空隙率（无量纲），d_p 为颗粒当量直径（m）。

方程由两项叠加而成，这一结构是理解它的关键。第一项与流速 u 成正比，来源于流体与颗粒表面的黏性剪切摩擦，称为黏性项（层流项）；第二项与 u² 成正比，来源于流体在颗粒间反复收缩、扩张、绕流产生的惯性损失，称为惯性项（湍流项）。两项共用 (1-ε)/ε³ 这一空隙率修正因子，反映了"固体越多、通道越窄、弯曲越剧烈，阻力越大"的几何本质。工程上常把两项之比视为判断流型的依据：黏性项占主导时床层接近层流，惯性项占主导时接近湍流，二者相当则处于过渡区——而工业固定床反应器绝大多数正落在过渡区，这正是必须用完整 Ergun 方程而非单一达西定律的原因。

二、适用条件：什么时候可以放心用

Ergun方程并非万能，其推导建立在若干假设之上，工程使用前应逐条核对：

• 颗粒床层、非空管：适用于由离散颗粒堆积而成的固定床，不适用于蜂窝整体式催化剂、规整填料或空管。整体式载体应改用通道型压降关联式。
• 颗粒尺寸相对均匀、球形度可表征：d_p 通常取等体积当量直径 d_p = 6V_p/S_p，非球形颗粒需引入球形度 φ_s 修正，实际当量直径取 φ_s·d_p。粒径分布过宽时，细粉填充大颗粒间隙，实测压降常高于按平均粒径计算值。
• 床层空隙率均匀：Ergun方程假设 ε 沿床层均一。若存在贴壁空隙率偏大、局部压实或分层装填，需分段计算或引入壁效应修正（如 Mehta-Schwartz 修正）。
• 雷诺数范围：以颗粒雷诺数 Re_p = ρ·u·d_p/μ 衡量，一般 Re_p 在约 1～1000 范围内适用性较好；Re_p 过小接近纯层流可用 Blake-Kozeny 式，过大进入完全湍流则惯性项系数可能偏离 1.75。
• 不可压缩或低马赫数流动：气相压降较大（如 ΔP/P 超过约 10%）时，气体密度沿床层明显变化，需按可压缩流分段积分，不能直接用入口密度一次算到底。

三、参数物理意义：每个量都在"说话"

理解参数意义，才能在结果异常时反推原因。μ 与 ρ 由流体物性和操作温度压力决定，气相压降对温度尤其敏感——温度升高 μ 增大、ρ 减小，黏性项上升而惯性项下降，净效应需实算。u 是空床表观流速，即体积流量除以床层截面积，不是颗粒间隙内的真实流速；间隙真实流速约为 u/ε，这也是压降随 ε 减小而急剧上升的原因。ε 是最敏感的参数，方程中它以 (1-ε)/ε³ 出现，ε 从 0.40 降到 0.35，该因子增大约 40%，压降随之显著上升，这正是装填密实度、颗粒破碎和粉尘堵塞影响压降的数学根源。d_p 出现在分母，粒径减半压降约增至 4 倍（黏性项）甚至更高，因此细粉夹带对压降的破坏力极大。

四、工程简化形式：抓住主导项

工程估算中常根据流型舍去次要项。低流速、液相或细颗粒体系，Re_p 较小，惯性项可忽略，方程简化为 Blake-Kozeny 形式：

ΔP/L ≈ 150·μ·(1-ε)²·u / (ε³·d_p²)   （黏性项主导）

高流速、大颗粒、气相体系，惯性项占优，可简化为 Burke-Plummer 形式：

ΔP/L ≈ 1.75·ρ·(1-ε)·u² / (ε³·d_p)   （惯性项主导）

需要强调的是，简化只用于快速估算或趋势判断，正式设计核算、压降异常诊断仍应使用完整方程。一个实用的判断做法是先算 Re_p，再算两项比值，若某项占比超过约 80% 方可安全简化。

五、示例演算：从参数到判断

以下为一组教学假设参数，用于演示计算流程，非实测数据。设气相通过催化剂床层：μ = 2.0×10⁻⁵ Pa·s，ρ = 1.2 kg/m³，u = 0.5 m/s，ε = 0.40，d_p = 3 mm = 0.003 m，L = 2 m。

先算几何因子：(1-ε)/ε³ = 0.60/0.064 = 9.375。黏性项 = 150×2.0×10⁻⁵×9.375×0.5/(0.003²) ≈ 1562 Pa/m；惯性项 = 1.75×1.2×9.375×0.25/0.003 ≈ 1641 Pa/m。两项相当，处于过渡区，必须用完整方程。总压降 ≈ (1562+1641)×2 ≈ 6406 Pa。此时 Re_p = 1.2×0.5×0.003/2.0×10⁻⁵ = 90，落在过渡区，与两项比值判断一致。

若运行中发现实测压降明显高于此估算值，应优先排查：ε 是否因装填过密或粉尘堵塞而下降、d_p 是否因破碎产生细粉、u 是否因流量偏大而升高。按前述敏感性，ε 下降 0.05 或细粉使有效 d_p 减半，都足以让压降成倍上升。把 Ergun方程当作"参数敏感性分析器"，而不只是算一个数，才是它在装填与分布分析中的真正价值。