COMSOL热流固耦合在煤体吸附膨胀研究中的应用

1. COMSOL热流固耦合在煤体研究中的核心价值

煤体吸附膨胀变形研究一直是能源与地质工程领域的重点难点。传统实验方法成本高、周期长，且难以捕捉微观尺度的物理化学过程。COMSOL Multiphysics作为全球领先的多物理场仿真平台，其热流固耦合模块（Thermo-Hydro-Mechanical, THM）通过耦合温度场、渗流场和应力场，为这类复杂问题提供了突破性的解决方案。

在实际煤层气开采中，煤体孔隙结构会随气体吸附/解吸过程发生动态变化。以甲烷为例，每吸附1mmol/g的甲烷，煤体体积膨胀率可达0.1%-0.5%。这种微米级的变形会显著改变孔隙结构和渗透率，进而影响气体运移效率。COMSOL的独特优势在于：

• 内置的煤岩变形本构方程（如Langmuir吸附膨胀模型）
• 非线性孔隙度-渗透率动态关联算法
• 多场耦合的隐式求解器

关键提示：在设置吸附膨胀参数时，需特别注意Langmuir体积常数(V_L)和压力常数(P_L)的取值。不同煤阶（褐煤、烟煤、无烟煤）这两个参数差异可达2-3个数量级。

2. 煤体吸附膨胀建模的关键技术路线

2.1 物理模型构建要点

建立准确的煤体吸附膨胀模型需要分层次考虑多个物理过程。以西南石油大学的路广团队研究为例，其建模流程包含三个核心环节：

1. 基质变形控制方程：

   math复制\nabla \cdot \sigma + \rho_b g = 0
   

   其中吸附应变项ε_s需通过扩展Langmuir方程引入：

   math复制ε_s = \frac{ε_L P}{P + P_L} \cdot \frac{C}{C_{max}}
   

   （C为吸附气体浓度，C_max为最大吸附容量）

2. 气体运移方程：
   采用修正的达西定律考虑Klinkenberg效应：

   math复制v = -\frac{k_{abs}}{\mu} (1 + \frac{b}{P}) \nabla P
   

   其中滑脱因子b与孔径分布直接相关。

3. 热力学耦合：
   吸附热效应会导致局部温度变化，需联立能量守恒方程：

   math复制(\rho C_p)_{eff} \frac{\p