详细讲解课程：“comsol煤层瓦斯抽采及瓦斯流动理论、钻孔抽采案例、煤与瓦斯耦合
案例分析、注气驱替影响分析”
凌晨三点的实验室键盘声格外清脆，我正在调整煤层渗透率参数。屏幕上突然弹出瓦斯压力云图异
常波动的警告——这已经是本周第三次遇到多物理场耦合震荡了。点燃第七根香烟，我突然意识到，搞煤层
瓦斯模拟就像在煤矿里玩电子游戏，既要懂流体力学规则，又要会写作弊代码。
在COMSOL中建立煤层几何模型时，新手常犯的典型错误是忽略煤岩各向异性。记得去年有个研究生
把渗透率张量设成了对称矩阵，结果钻孔周围的瓦斯扩散轨迹画出了抽象派油画。正确的材料属性设置应
该是这样的：
```matlab
% 渗透率张量设置（m）
K_xx = 1e-15 * cos(theta)^2 + 1e-17 * sin(theta)^2;
K_yy = 1e-17 * cos(theta)^2 + 1e-15 * sin(theta)^2;
K_xy = (1e-15 - 1e-17)*sin(theta)*cos(theta);
```
这里的theta是煤层主渗透方向与坐标轴的夹角，这个参数不实测真的会死人。去年某矿的注气驱
替项目就因为这个角度偏差15度，导致实际抽采量比预测少了37%。
钻孔抽采案例最刺激的部分永远是边界条件设置。当你在达西定律接口里把抽采压力设为-0.95MP
a时，系统突然弹窗警告负压超限的样子，像极了ATM机吞卡时的冷漠表情。这时候需要祭出分段函数大法
：
```java
if (t < 3600*24, -0.8e6, -1.2e6)
```
这个时间相关的压力边界，完美模拟了矿上先缓后急的骚操作。记得把时间步长控制在0.1秒以内，
否则瓦斯浓度波阵面会变成锯齿状心电图。
说到煤与瓦斯耦合，固体力学模块里的塑性变形方程和达西流的联姻简直虐心。有次我把煤体膨胀
系数设反了，结果模型里的煤层在注气后不是膨胀而是塌缩，活生生把瓦斯抽采模拟成了黑洞吞噬。正确
的耦合方式应该是：
```matlab
epsilon_p = beta*(C - C0); // 孔隙应变与瓦斯浓度关系
perm = perm0 * exp(alpha*epsilon_p); // 渗透率动态变化
```
这个非线性关系让求解器经常闹脾气，需要准备三种以上迭代方法备胎。建议把牛顿迭代法的最大
步数调到50，再备好龙格库塔法当Plan B。
