COMSOL模拟煤层封存二氧化碳。 该案列分别从水平井和垂直井两个角度分析煤层封存二氧化碳，分析裂隙和基质压力以及煤层渗透率，附有详细的SP讲解。
地下六百米深处的煤层正默默上演着神奇的气体收纳术。以山西某煤矿为例，我们尝试用COMSOL还原二氧化碳注入的真实场景——当高压气体顺着井筒钻进煤层缝隙，岩层像吸饱水的海绵般发生着肉眼不可见的形变。
先看水平井模型。建立几何结构时特别要注意煤岩的裂隙网络，这里有个实用技巧：直接在材料属性里用概率函数生成随机裂隙分布。代码片段中这个for循环控制着裂隙密度随深度变化的梯度：
for (int i=0; i<layerCount; i++) { double depthFactor = 1 - (currentDepth/maxDepth); fractureDensity[i] = baseDensity * Math.pow(depthFactor, 2.3); }相当于在每层岩体中，越靠近地表裂隙越稀疏。配合达西定律模块设置时，初始渗透率别直接填常数，而是关联裂隙开度的二次方函数k=∫(a^3/12μ)ds，这样才能真实反映裂隙对导流能力的指数级影响。
垂直井的模型更考验网格划分技巧。在井筒周围需要做三层边界层网格，最内层单元尺寸控制在0.1米级。这里有个骚操作——在求解器设置里勾选"非线性渐变"，配合自动时间步进，能有效避免气体突窜导致的发散问题。监测点数据导出时记得用探针功能，实时捕获基质孔隙压力波动，像这样截取瞬态数据：
probe1 = mphprobe(model,'pore_pressure','selection',[23,45,67]); timeData = probe1.getTime(); pressureData = probe1.getData();对比两种井型的结果很有意思。水平井注气第30天，距井筒50米处的裂隙压力已达8MPa，而基质仅5.2MPa，这1.8MPa压差驱动着气体持续渗入微孔隙。渗透率演化曲线显示，前两周k值暴涨40%，之后趋于平稳——对应着裂隙撑开、煤体膨胀的动态平衡过程。
垂直井的表现则大相径庭。在800米注气段，气体优先沿原生裂隙纵向突进，形成倒锥形压力分布。有个反直觉现象：注气速率提高20%时，封存量反而下降15%。究其根源，过快的注入导致裂隙过度扩展，反而连通了外围低渗区形成泄漏通道。
当把吸附膨胀效应纳入计算时，模型开始展现煤层的"记忆效应"。在材料设置里勾选Swelling Strain选项后，渗透率衰减呈现先快后慢的特征。这里有个防坑指南：膨胀系数别直接照搬文献值，应该用现场岩芯的应力-应变曲线反演参数，否则可能偏差一个数量级。
最终对比显示，水平井在单位时间封存量上占优，但垂直井的波及范围更均匀。实际工程中得玩"排列组合"，比如采用从林式井网，让水平井做主封存通道，垂直井作压力调节器。下次如果再优化模型，可以试试在注气后期切换氮气驱替，说不定能多压进去15%的二氧化碳。
