多场耦合下煤层中CO2封存与甲烷驱替研究：涉及流固耦合、二元气体竞争吸附及多场动态变化
comsol流固耦合二氧化碳封存-煤层注气驱替甲烷 研究煤体吸附膨胀变形 涉及 渗流场 温度场 变形场 二元气体竞争吸附 孔隙率 渗透率动态变化 甲烷、二氧化碳封存量 参考The coupling mechanism of the thermal-hydraulic-mechanical fields in CH4-bearing coal and its application in the CO2-enhanced coalbed methane recovery 论文 文件中包括：讲解视频 模型（自定义物理场输入函数）参考论文 不
煤层注CO₂驱替甲烷这事儿听着像给煤做SPA，不过里头藏着的物理场纠缠可比按摩复杂多了。咱们得同时盯着气体怎么挤进煤缝（渗流场）、温度变化怎么捣乱（温度场）、煤块自己还膨胀收缩（变形场）。COMSOL里玩这个得先搞明白自定义物理场怎么喂数据，就拿煤吸附CO₂膨胀的模型来说，得自己写个函数表达煤颗粒变胖的幅度。
function d_swell = adsorption_swell(P_CO2, T) % 引用论文中的吸附应变经验公式 Q_CO2 = 32e3; % 吸附热J/mol R = 8.314; epsilon_max = 0.015; % 最大膨胀应变 d_swell = epsilon_max * (P_CO2.*exp(Q_CO2/(R*T))) ./ (1 + P_CO2.*exp(Q_CO2/(R*T))); end这代码其实在说：温度越高，煤吸CO₂越来劲，但吸饱了就开始摆烂（符合Langmuir型曲线）。注意指数项里藏着温度对吸附能力的调控，25℃和60℃工况下煤的膨胀量能差出两倍多。
渗透率动态变化才是真·戏精。煤一膨胀，孔隙就被挤压，咱们得用应力-渗透率耦合模型：
// COMSOL中自定义渗透率函数 if (isLinearized) { permeability = k0 * (1 + alpha_phi*(solid.dispXx + solid.dispYy + solid.dispZz)); } else { // Palmer-Mansoori修正模型 permeability = k0 * Math.exp(-3*effective_stress/(E_coal*(1-nu_coal))); }这里有个骚操作——当煤体受压变形超过临界值，渗透率断崖式下跌。模拟时得盯着应力状态自动切换计算公式，否则计算结果会在突变区域疯狂震荡。
二元气体竞争吸附得像处理情敌关系似的。CO₂和CH4在煤表面抢位置，得用扩展Langmuir方程：
q_CO2 = (b_CO2*P_CO2)/(1 + b_CO2*P_CO2 + b_CH4*P_CH4) q_CH4 = (b_CH4*P_CH4)/(1 + b_CO2*P_CO2 + b_CH4*P_CH4)但实际参数得考虑温度影响，比如60℃时CO₂的吸附系数b值会比常温下降40%，这直接导致高温工况下驱替效率下降——所以注CO₂前得先评估煤层温度场是否适合。
comsol流固耦合二氧化碳封存-煤层注气驱替甲烷 研究煤体吸附膨胀变形 涉及 渗流场 温度场 变形场 二元气体竞争吸附 孔隙率 渗透率动态变化 甲烷、二氧化碳封存量 参考The coupling mechanism of the thermal-hydraulic-mechanical fields in CH4-bearing coal and its application in the CO2-enhanced coalbed methane recovery 论文 文件中包括：讲解视频 模型（自定义物理场输入函数）参考论文 不
搞完这些物理场，最后算封存量时得注意煤基质与裂隙的双重存储机制。裂隙里的自由气体用达西定律算，吸附态气体得结合煤比表面积参数：
adsorbed_CO2 = volume_coal * rho_coal * q_CO2 * M_CO2 / Vm free_CO2 = porosity * P_CO2 * M_CO2 / (Z * R * T)有个坑点：当渗透率剧烈变化时，自由相占比会突然跃升，这时候时间步长必须调小到秒级，否则质量不守恒误差能吃掉30%的结果精度。
模型验证时，拿现场气水比数据和模拟结果对比。发现有个月的数据对不上，折腾一周才发现是没考虑注气导致的温度下降——煤吸附CO₂是放热过程，局部温升能到5-8℃，这反过来又影响吸附平衡。加上热力学耦合后，产气曲线终于和实测吻合了。
玩这个模型最大的感悟：别相信默认求解器设置！流固耦合的非线性得用手动调整牛顿迭代步长，还得开着误差估计随时救场。有一次求解器在渗透率突变处卡住，改成间断伽辽金法才搞定，算完感觉头发都白了几根。