双孔介质瓦斯流动模型与COMSOL模拟实践

1. 双孔介质瓦斯流动模型原理剖析

煤层瓦斯流动模拟本质上是一个多物理场耦合问题。与常规单孔介质模型不同，双孔介质模型将煤层结构划分为两个相互作用的系统：基质系统和裂隙系统。这种划分并非人为假设，而是基于煤体显微CT扫描显示的天然二元结构特征。

基质系统由煤的微孔隙构成，孔径通常在纳米级（2-50nm），瓦斯主要以吸附态存在于孔隙表面。这部分瓦斯释放遵循扩散机制，可以用修正的菲克第二定律描述：

code复制∂C_m/∂t = D_m∇²C_m - Q

其中D_m是等效扩散系数，包含了孔隙曲折度的影响因子。实测数据显示，典型烟煤的D_m值在10^-12~10^-10 m²/s量级，且随压力呈指数变化。

裂隙系统则相当于煤层的"高速公路网络"，其渗透率可比基质高出3-6个数量级。这里的流动符合达西定律，但需要考虑Klinkenberg效应——气体滑脱现象导致的有效渗透率升高：

code复制K_eff = K∞(1 + b/P)

b是滑脱因子，对于甲烷在煤裂隙中的流动，b值通常在0.1-1MPa范围。这个效应在低压抽采阶段尤为显著，可能使计算流量偏差达30%以上。

2. COMSOL模型搭建关键步骤

2.1 几何建模与网格划分

建议采用"外圆内方"的几何策略：外部计算域设为圆形（半径至少为预计抽采半径的3倍），中心布置抽采孔。这种设置既能保证边界效应最小化，又便于后处理分析。

网格划分需要特别注意：

• 基质区域采用较密的三角形网格（最大单元尺寸<0.1m）
• 裂隙网络周边进行边界层加密（3-5层，厚度系数0.3）
• 抽采孔附近设置局部细化（网格尺寸梯度比不超过1.5）

经验提示：先进行网格独立性验证，逐步加密直到关键参数（如累计抽采量）变化<2%

2.2 物理场耦合设置

在COMSOL中需建立两个PDE接口分别对应基质和裂隙系统。耦合项的实现有两种方式：

1. 弱形式PDE直接添加源项
2. 使用系数型边界偏微分方程

推荐第二种方法，其稳定性更好。具体参数设置：

code复制耦合系数 α = σ·k_m/D_m
σ为形状因子（典型值0.1-1mm^-2）
k_m为基质渗透率（10^-21~10^-19m²）

2.3 初始条件与边界条件

初始瓦斯压力场建议采用指数衰减分布：

code复制