解决方案
【MS应用实例】西安科技大学李树刚/林海飞成果展示：含吸附甲烷的煤注入N2和CO2的分子模拟
摘要
本文采用蒙特卡罗（GCMC）模拟方法研究N2和CO2替换甲烷的动力学机制，确定了含吸附甲烷的最低能量煤模型。研究了不同温度下，注入N2和CO2后CH4的脱附行为。结果表明，CO2和N2主要通过占据吸附位置方式来驱离甲烷气体。CO2注入使CH4的扩散活度降低了20.53%，CO2对甲烷的脱附作用优于N2。
前言
全球能源作为一个重要的基础性问题，已被列入人类生产生活的日程。因此，我们应该关注如何最好地利用不可再生资源，包括天然气和石油煤。煤层气(CBM)主要由甲烷(CH4)组成，CO2、N2、O2、水蒸气等气体含量较低。因此，在现有的清洁能源中，它的利用被认为是一种节能、性价比高的替代能源。从分子水平确定煤与瓦斯的相互作用机理一直是煤层气研究开发中的一个基础性问题。近几年， KongL等人在单轴应力和分层预压缩成型条件下进行了CO2 /N2注入煤层气排放的模拟实验，推导出了煤层气对CO2 /N2耦合数学模型。通过计算机分子模拟可以从分子水平研究吸附剂与吸附质之间的微观吸附机理。Matranga和Liu X等人采用蒙特卡罗方法模拟了活性炭上的CH2, CO2和N2的吸附行为。结果表明，CO2的吸附量和吸附热均高于CH2、N2。因此，需要进一步分析N2和CO2注入对含吸附甲烷煤模型中CH4解吸动力学的影响，才得以更深入地了解CH4解吸机理。
理论方法：
西安科技大学李树刚/林海飞教授团队将优化方法设置为medium，所选力场为COMPASS(Materials studio），采用Ewald和基于原子的方法计算了范德华力和氢键，搭建的模型如图1所示。
各模型的吸附量、吸附热随压力的关系如图2所示。图2a为含吸附甲烷煤模型在不同温度下的吸附量与压力的关系。并给出了甲烷在不同温度下的吸附构型。图2b为不同温度下吸附热与压力的关系。结果表明，在不同温度下，模型中甲烷浓度的分布是一致的。
结果与讨论：
如图3所示，CH4-N2和CH4-CO2体系初始模型和解吸模型总能均随温度升高而增加，且与温度呈线性正相关。解吸后的CH4-N2和CH4-CO2构型下，体系总能量增加，表明脱附为吸热反应。
如图4所示，模型中，随着CO2和N2的注入，真空层中CH4的相对浓度增大。随着温度的升高，真空层中CO2、N2和CH4相对浓度均增加，总体关系为CH4>CO2>N2。这种变化说明部分CH4扩散到狭缝区域，CO2和N2的相对浓度分布仍集中在煤分子层中。
在图5中，CH4-N2和CH4-CO2模型中各气体的MSD随温度的升高而逐渐升高，这说明体系中的分子活性随温度升高而增大。
Fig. 5. The MSD curves and desorption configurations of different models at different temperatures.
由图6可知，CH4- N2和CH4- CO2模型中气体的MSD均随时间增加;在同一时间内，气体的MSD随温度的升高而升高。如图6所示，曲线拟合度R2均大于0.98，相对较高，从而保证了计算得到的扩散系数d的可靠性。随着温度的升高，气体的扩散系数逐渐增大。温度越高，分子的热力学能就越大。根据能量守恒原理，内能将转化为动能，使分子作剧烈运动。
Fig. 6. The fitting curves of diffusion coefficients of CH4, CO2, and N2 in the coal model at different temperatures.
结论
西安科技大学李树刚/林海飞教授团队表明：
(1) CO2和N2主要通过占据吸附位置驱赶CH4。CO2注入模型比N2注入模型显示出更多的游离甲烷分子。CH4-N2和CH4-CO2体系的初始模型和解吸模型总能量均随温度升高而增加。从微观角度看，CH4- CO2模型比CH4- N2模型更稳定。
(2) 随着温度的升高，真空层中CH4的平均相对浓度和运动速度增大。真空层中三种气体的平均相对浓度与平均速度分布的关系为CH4>CO2>N2。
(3) 三种气体的扩散活化能关系为CH4>N2>CO2。CH4在CH4- N2和CH4- CO2模型中的扩散活化能分别为5.26和4.18 kcal/mol。CO2注入模型中，CH4的扩散活化能降低了20.53%。
(4) CO2注入对甲烷脱附效果优于N2。本文可为煤层瓦斯喷射技术的研究和发展提供一定的依据。
(5) 该研究可为煤层气注气增产技术的发展提供理论依据。后续研究将结合具体的煤岩参数对研究结果进行细化。
文章详情请见：https://doi.org/10.1016/j.energy.2020.119686
本文是通过Materials Studio顺利完成，其中CO2/N2占据CH4位置等动力学模拟在FORCITE模块下实现。北京泰科是Materials studio官方代理商，更多软件详情或者技术支持，请咨询北京泰科。
电话：010-64951848
邮箱：sales@tech-box.com.cn
作者简介
李树刚 西安科技大学教授
博士，教授，博导。百千万人才工程国家级人选，国家有突出贡献中青年专家，享受国务院政府特殊津贴，陕西省“高层次人才特殊支持计划”教学名师。长期从事煤与瓦斯安全共采、采场矿山压力与岩层控制的教学与科研工作。
林海飞 西安科技大学教授
博士，教授，博导。获评陕西省青年科技新星、中国职业安全健康协会青年科技奖、陕西省煤炭工业优秀科技工作者。现任陕西省煤炭工业协会、陕西省煤炭学会安全工程专业委员会副秘书长。主要从事煤矿瓦斯灾害防治理论和技术的教学与科研工作，研发出固气耦合相似材料模拟系统及相应的实验材料，提出矩形梯台带抽采卸压瓦斯技术。
公司简介
北京泰科博思科技有限公司(Beijing Tech-Box S&T Co. Ltd.）成立于2007年，是国内领先的分子模拟及虚拟仿真综合解决方案提供商。
北京泰科博思科技有限公司与国际领先的模拟软件厂商、开发团队深入合作，为高校、科研院所和企业在材料、化工、药物、生命科学、环境、人工智能及数据挖掘、虚拟仿真教学等领域提供专业的整体解决方案。用户根据需要在我们的平台上高效的进行各种模拟实验，指导实际的生产设计。
北京泰科博思科技有限公司拥有一支一流的技术服务团队和资深的专家咨询团队，以客户真正需求出发，服务客户，为客户创造价值。我们秉承“职业、敬业、担当、拼搏、合作”的企业精神，致力于用国际领先的软件产品和专业全面的技术支持服务，成为客户可信赖的合作伙伴。
©北京泰科博思科技有限公司 / 京ICP备09107432号-1