3. Langmuir吸附模型理论推导
Langmuir吸附模型是煤层气领域描述煤岩吸附等温线最常用的理论工具。它把复杂的多孔介质吸附过程抽象为一个可解析的物理图像,从而用两个参数——最大吸附量和吸附平衡常数——概括整条等温线的形状。理解这一模型的推导逻辑,比记住公式本身更重要,因为只有弄清假设从何而来,才能判断在什么条件下拟合结果可信、在什么条件下会系统性偏离。
单分子层吸附假设的物理图像
Langmuir模型建立在四个核心假设之上。第一,固体表面存在数量有限、能量均一的吸附位点,每个位点最多容纳一个吸附质分子,因此吸附只能形成单分子层,不存在多层堆叠。第二,吸附质分子之间没有横向相互作用,已吸附分子不影响邻近空位的吸附能力。第三,吸附与脱附是动态可逆过程,当吸附速率等于脱附速率时达到吸附平衡。第四,表面均匀,所有吸附位点的吸附亲和力相同。
把这套图像映射到煤岩上:煤的微孔表面可视为大量吸附位点的集合,甲烷分子被范德华力束缚在孔隙壁上,形成近似单层的吸附相。需要强调的是,煤岩实际表面存在非均质性和微孔填充效应,因此这些假设只是近似,这正是后文讨论局限性的出发点。
Langmuir方程的推导
设表面总吸附位点数为 N,已占据位点数为 Na,覆盖度 θ = Na/N。气相压力为 p,吸附速率正比于气相分子与空位的碰撞频率,即 ka·p·(1−θ);脱附速率正比于已占据位点数,即 kd·θ。平衡时两者相等:
k_a · p · (1 − θ) = k_d · θ
整理得 θ = (ka/kd)·p / [1 + (ka/kd)·p]。令 b = ka/kd,并注意到吸附量 V = VL·θ(VL 为饱和吸附量),即得 Langmuir 方程的标准形式:
V = V_L · b · p / (1 + b · p)
其中 V 为单位质量煤在压力 p 下的吸附量(常用 m³/t),VL 为 Langmuir 体积(最大单层吸附量),b 为 Langmuir 压力常数(MPa⁻¹)。为便于线性拟合,常改写为倒数形式:
p / V = p / V_L + 1 / (V_L · b)
以 p/V 对 p 作图,斜率为 1/VL,截距为 1/(VL·b),由斜率和截距即可反算 VL 与 b。这是实验室处理等温线数据最常用的线性回归方式。
吸附平衡常数 b 的物理意义
b 的量纲是压力的倒数,其倒数 1/b 具有明确的物理含义:当 p = 1/b 时,覆盖度 θ = 0.5,即吸附量达到饱和吸附量的一半。因此 1/b 可理解为“半饱和压力”,它衡量吸附位点被占据的难易程度。b 值越大,1/b 越小,说明在更低压力下即可达到半饱和,煤对甲烷的亲和力越强。
从热力学角度,b 与吸附平衡的自由能变化相关,可写成 b = b0·exp(−ΔH/(RT)) 的形式,其中 ΔH 为吸附焓。这意味着 b 随温度升高而减小,温度越高吸附越弱,这与“升温脱附”的工程经验一致。在煤层气开采中,b 的大小直接影响等温线在原始储层压力附近的斜率,进而影响降压解吸的难易判断。
适用条件与局限性
Langmuir 模型适用于以微孔为主、吸附以单层物理吸附为主导的体系,在煤岩甲烷吸附的中低压段通常拟合良好。但在实际应用中需注意以下局限:
- 高压段偏差:煤岩存在微孔填充和少量多层吸附,高压下实测吸附量可能继续缓慢上升,而 Langmuir 曲线已趋于平台,导致 VL 被低估或 b 被高估。
- 表面非均质性:煤表面吸附位点能量不均一,强吸附位点优先占据,实际等温线在低压段可能比 Langmuir 曲线更陡,此时可考虑双位点 Langmuir 或 Freundlich 类模型。
- 多组分竞争:Langmuir 方程描述单组分吸附,注 CO₂ 驱替或含水分条件下需引入扩展 Langmuir 方程处理竞争吸附。
- 温度单一性:VL 与 b 均随温度变化,一组参数只对应一个温度,跨温度使用必须重新拟合。
综上,Langmuir 模型以清晰的单分子层物理图像和简洁的两参数形式,为煤岩吸附等温线提供了可解释、可外推的定量框架。掌握其推导逻辑与 b 的物理含义,并在拟合中主动识别高压段偏差与非均质性影响,才能让这一经典模型在煤层气储层评价中发挥应有的作用。