东曲矿8 号煤CO2 和CH4 竞争吸附特性分子模拟研究

邓小鹏 ，相建华

（1.山西能源学院，山西 晋中 030600；2.太原理工大学 矿业工程学院，山西 太原 030024）

CO2-ECBM 开采是一种相对高效的煤层气开采技术[1]。煤层中CO2和CH4竞争吸附的微观机制是制约CO2-ECBM 的主要机制[2-3]。煤是一种复杂的大分子化合物，因此，从分子水平研究煤层中CO2和CH4的竞争吸附机制至关重要。

与其他大分子有机化合物不同，煤大分子具有分子组成多样化和化学结构复杂的特征[4-5]。因此，合理地构建煤的大分子模型是在分子水平上研究煤的结构和煤中气体的竞争吸附机理的前提条件。目前，常用X 射线光电子能谱（XPS）、拉曼光谱（Raman spectroscopy）、高分辨率透射电子显微镜（HRTEM）、13C 核磁共振（13C-NMR）和傅里叶变换红外光谱（FTIR）等技术对煤的分子结构进行研究和分析[6-9]。相建华等[10]基于成庄矿煤的工业分析、元素分析、13C-NMR、XPS 等实验结果，构建了其大分子结构模型，并采用分子力学和分子动力学方法对所得煤结构模型进行结构优化，结果表明分子间芳香层片的定向排列是高煤级煤结构中短程有序的原因；王绍清等[11]基于HRTEM研究了煤的石墨化过程，指出石墨微晶形成过程还伴随着结构缺陷的演化；NIEKERK 等[12]利用HRTEM 测量了煤的芳香条纹的大小和分布，并根据13C-NMR 在芳香骨架上加入硫、氮、氧、脂肪链侧链和交联键，从而建立了2 种南非煤的分子模型；MOKONE 等[13]利用岩相分析、元素分析、氦密度、13C-NMR、HRTEM 等方法构建了含氧、氮、硫的多环芳烃分子结构，提供了更好、更准确的芳香烃测量方法。随着HRTEM 等技术在煤结构分析中的应用，学者们对煤大分子碳骨架结构的认知有了更进一步的提高，这些技术的进步都有利于更全面的表征煤大分子的结构特征，从而为气体在煤中的竞争吸附模拟奠定基础。

学者们针对CO2和CH4在煤中的竞争吸附特性展开了一系列研究。ZHENG 等[14]利用核磁共振手段研究了亚烟煤和无烟煤中CO2与CH4的竞争吸附特性，结果表明CO2的注入可使吸附CH4的解吸效率提高14%～26%；WEI 等[15]对中煤级煤中CO2和CH4气体的竞争吸附进行了实验研究，结果表明CH4浓度的下降速率随着CO2注入压力的增加表现为先增大后趋于稳定的趋势；DANG等[16]利用分子模拟方法研究了褐煤中CO2和CH4的吸附行为，结果表明CO2对CH4的吸附选择性随着压力的增加表现为先迅速下降后趋于稳定的趋势，这说明气体压力是影响竞争吸附行为的重要因素。此外，也有学者认为气体的分子动力学直径也是CH4和CO2在煤中竞争吸附行为的主导因素[17]。虽然有部分学者采取一系列方法研究了CO2和CH4在煤中的竞争吸附行为，但对于CO2和CH4在分子水平的竞争吸附机制仍需要进一步探讨。

合理的大分子模型是模拟CO2对CH4在煤中竞争吸附性能的基础。为此，基于东曲矿8 号煤的大分子结构模型开展了CO2和CH4的单组分与双组分竞争吸附研究，分析了不同吸附摩尔比条件下的吸附量和吸附选择性参数，并从相互作用能和势能分布角度揭示了竞争吸附机理，研究结果对于实现CO2置换CH4和提高煤层气采收率具有重要的意义。

1 模型结构及模拟方法

东曲矿8 号煤的工业分析与元素分析结果为：水分含量2.12%，灰分含量4.06%，挥发分含量13.24%，最大镜质组反射率1.82%，C 元素含量91.52%，H 元素含量4.04%，O 元素含量2.46%，N 元素含量1.43%；S 元素含量0.56%。

通过HRTEM、13C-NMR、XPS 和FTIR 构建了东曲矿8 号煤的大分子结构（C1 258H644O24N16S2），模型构建过程如图1。

图1 煤大分子结构模型的构建流程图Fig.1 Flowchart for the construction of structural model of coal macromolecules

图1 中，HRTEM 用于确定煤大分子的碳骨架结构，13C-NMR、XPS 和FTIR 用于确定脂肪结构和官能团结构的分布。通过Materials Studio 软件对煤大分子结构进行分子力学与分子动力学优化，力场设置为COMPASS Ⅱ，最终得到东曲矿8 号煤的最低能量构型（图1d）。由于煤中存在一些难以清除的小分子[17]，其模型的密度（1.40 g/cm3）略低于实测密度（1.42 g/cm3）。在Atom Volume &Surfaces 模块中运用半径为1.65 nm 的圆球探针对孔隙空间进行探测（图1e），所得煤大分子中的孔隙体积为6 493.21 Å3（1 Å=10-10m）。

采用Materials Studio 软件中的吸附模块进行了吸附模拟。模拟温度和最大吸附压力分别设置为298.15 K 和30 MPa，根据吸附质气体压力确定气体摩尔比。为了研究CH4和CO2的竞争吸附特性，采用了单组分与不同摩尔比的双组分气体吸附方案。气体吸附模拟参数见表1。

表1 气体吸附模拟参数Table 1 Gas adsorption simulation parameters

基于蒙特卡洛方法对煤大分子模型进行CH4和CO2等温吸附模拟，计算固定压力下的等温吸附量来表征气体在特定压力下的吸附性能。模拟过程采用COMPASS Ⅱ力场，平衡步数为100 000步，生产步数为1 000 000 步，并计算平均吸附量。通过Langmuir 公式获得了气体的吸附等温线[18]。朗格缪尔拟合公式为：

式中：V为气体吸附量，mmol/g；a为压力趋于无穷大时的最大吸附量，mmol/g；b为吸附常数，1/MPa；p为吸附压力，MPa。

2 模拟结果

2.1 单组分气体吸附

东曲矿8 号煤单组分气体的等温吸附曲线如图2，单组分等温吸附曲线的朗格缪尔拟合参数见表2。

表2 单组分等温吸附曲线的朗格缪尔拟合参数Table 2 Langmuir fitting parameters of single component isothermal adsorption curves

图2 单组分气体的等温吸附曲线Fig.2 Isothermal adsorption curves for single component gas

单组分气体吸附条件下CH4和CO2的气体吸附量均随着压力的增大而呈现先快速增加后稳定的趋势，符合I 型等温吸附特征。单组分吸附条件下CO2的吸附量显著大于CH4的吸附量，为煤中CH4的置换提供了理论依据。CO2的吸附常数b的值也大于CH4的吸附常数，这说明CO2比CH4更容易达到最大吸附量。主要原因是CO2与煤大分子之间的作用力较强，且CO2分子之间具有较强的相互作用力[19]。

2.2 双组分竞争吸附

2.2.1 气体吸附量

双组分竞争吸附下的等温吸附曲线如图3，双组分竞争吸附等温吸附曲线的朗格缪尔拟合参数见表3。

表3 双组分竞争吸附等温吸附曲线的朗格缪尔拟合参数Table 3 Langmuir fitting parameters of dual component competitive adsorption isothermal adsorption curves

图3 双组分竞争吸附的等温吸附曲线Fig.3 Isothermal adsorption curves for dual component competitive adsorption

双组分竞争吸附下的等温吸附规律同样符合Ⅰ型等温线特征。在CO2/CH4竞争吸附体系中，总吸附量随着CO2的摩尔分数的增大而增大；当CO2的摩尔分数为60% 时总吸附量达到最大值（4.21 mmol/g），当CO2的摩尔分数为40%时总吸附量最小（4.06 mmol/g）。对比3 种不同摩尔比吸附条件下的吸附常数可以发现，CO2的摩尔分数越大，竞争吸附越容易达到最大吸附量。双组分气体的吸附特性与竞争吸附过程中的吸附选择性有关。

2.2.2 吸附选择性

在竞争吸附体系中，吸附选择性系数是评价气体间竞争吸附的重要指标之一。气体竞争吸附体系中CO2对CH4的吸附选择性系数SCO2/CH4的定义如下[20]：

式中：SCO2/CH4为CO2对CH4的吸附选择性系数；xCO2、xCH4为CO2和CH4在吸附相中的摩尔分数；yCO2、yCH4为体积相中CO2和CH4的摩尔分数。

当SCO2/CH4＞1 时，表明CO2容易吸附于煤大分子中，选择性系数越高，吸附效果就越好。CO2和CH4的在煤大分子中的双组分竞争吸附的选择性系数如图4。

图4 不同吸附体系中CO2 对CH4 的吸附选择性系数Fig.4 Adsorption selectivity coefficient of CO2 to CH4 in different adsorption systems

对比不同摩尔分数条件下的吸附选择性系数随气体压力的变化规律可以发现：不同条件下CO2对CH4的选择性吸附系数始终大于1；当CO2的摩尔分数为60%和50%时，吸附选择性系数随着气体压力的增大而增大；当CO2的摩尔分数为40%时，吸附选择性系数随着气体压力的增大而呈减小趋势；不同摩尔比的CO2和CH4竞争吸附体系中，吸附选择性系数SCO2/CH4随着CO2摩尔分数的减小而增大。因此，较低摩尔分数的CO2具有较强的竞争吸附性能，有利于采收率的提高。

2.3 相互作用能

相互作用能指的是吸附前后系统总能量的变化量，CH4和CO2在煤大分子中的吸附主要由范德华能和静电能引起[21]。因此，进一步分析了单组分与双组分竞争吸附条件下的相互作用能的变化。不同吸附体系的相互作用能如图5。

图5 不同吸附体系的相互作用能Fig.5 Interaction energy of different adsorption systems

由图5（a）可知，单组分气体吸附条件下，相互作用能随着吸附量的增加而增加，其中，CH4吸附体系的范德华能介于-70.71 kcal/mol 和-271.25 kcal/mol 之间（1 kcal=4.186 kJ），CO2吸附体系的范德华能介于-68.87 kcal/mol 和-254.98 kcal/mol之间，2 种单组分气体吸附体系的范德华能并无明显差异；CH4吸附体系的静电能介于-0.44 kcal/mol和-2.01 kcal/mol 之间，而CO2吸附体系的静电能介于-34.07 kcal/mol 和-140.29 kcal/mol 之间。静电能在CO2与煤大分子的吸附过程中发挥了重要的作用。气体在煤中的等量吸附热可用于判断物理吸附与化学吸附[22]，与CH4吸附往较低的吸附热不同，CO2吸附往往具有较大的吸附热，这与静电作用引起的化学吸附具有密切的关系。

由图5（b）可知，双组分竞争吸附条件下，双组分竞争吸附体系的范德华能介于-15.34 kcal/mol和-295.03 kcal/mol 之 间，静 电 能 介 于-14.56 kcal/mol 和-90.86 kcal/mol 之间；不同摩尔比体系的范德华能与静电能均随着CO2摩尔分数的增加而显著增加；如总吸附量为1 mmol/g 的条件下，当CO2的摩尔分数由40%增加至60%，吸附体系的 范 德 华 能 由- 15.63 kcal/mol 增 加 至-72.66 kcal/mol，静电能则由-14.56 kcal/mol 增加至30.23 kcal/mol；且吸附量越大，由CO2摩尔比增加导致的范德华能和静电能的增加幅度也越大。这是由于体系中总吸附量随着CO2的摩尔分数的增大而增大，从而导致相互作用能的增加。

2.4 势能分布

气体分子的势能分布有利于判别气体在煤大分子中的竞争吸附位点[23]。吸附势的绝对值越大，吸附作用就越强，吸附位点也更具优势。单组分与双组分竞争吸附条件的势能分布如图6。

图6 不同吸附体系的势能分布Fig.6 Potential energy distribution of different adsorption systems

由图6 可知，势能分布形式基本表现为正态分布；单组分CO2的吸附势于-7.05 kcal/mol 时达到优势吸附位点，高于CH4的吸附势（-3.95 kcal/mol）。煤中含有一定数量的羧基和羟基的等含氧官能团，而CO2比CH4优先吸附这些极性官能团[10]，与煤大分子的相互作用更强，因此，CO2的吸附位点更具优势。

对比3 种不同摩尔比条件下的势能分布可以发现：CO2的摩尔比可以改变CH4吸附势的分布，当CO2的摩尔比由60%降低至50%和40%，CH4的吸附势由-4.15 kcal/mol 降低至-4.05 kcal/mol 和-3.85 kcal/mol。CO2的摩尔比的降低抑制了CH4的吸附。CH4吸附势的降低有利于CO2对CH4的竞争吸附，因此，较低的CO2摩尔比有利于采收率的提高，这与吸附选择性分析中得到的结果是一致的。

3 结 语

1）单组分吸附模拟中CO2的吸附量和吸附常数显著大于CH4的吸附量。双组分竞争吸附模拟中的总吸附量随着CO2的摩尔分数的增大而增大，吸附常数也表现出相同的规律。

2）不同摩尔比条件下的双组分吸附模拟中CO2对CH4的选择性吸附系数始终大于1。CO2摩尔分数越大，选择性吸附系数越小，较低摩尔分数的CO2具有较强的竞争吸附性能。

3）相互作用能随着吸附量的增大而显著增大，CO2吸附体系中较大的静电能促进了煤大分子对CO2的吸附，因此，不同摩尔比体系的相互作用能随着CO2摩尔分数的增加而显著增加。

4）单组分条件下CO2的吸附势大于CH4的吸附势。双组分竞争吸附条件下CH4的吸附势随着CO2摩尔比的降低而降低，这有利于CO2对CH4的竞争吸附，验证了吸附选择性分析的结果。