Construction and optimization of lignite molecular structure model
-
摘要: 针对煤分子建模方法较多且不统一的现状,本文从微观角度出发,选取对分子结构影响较小的物理方法,探究了褐煤大分子结构并对其进行优化,提出了一种较为简便的煤大分子建模方法。基于元素分析、核磁共振碳谱测试及X射线光电子能谱测试的实验方式,利用分子动力学模拟软件对褐煤进行了分析研究并搭建分子结构模型。研究表明,褐煤大分子中芳香结构多为并五苯,脂肪碳结构多为亚甲基及次甲基结构,烷烃链多为环烷烃结构;氧原子以醚键氧形式为主,其次为羧基氧和羰基氧;氮原子以吡啶氮结构形式存在。优化后煤分子模型结构更加紧凑,能量明显变小,范德华势能作为非键势能的主要组成部分及保持煤分子结构稳定的主要因素,变化最明显。为了构建煤分子晶胞,加入周期性边界条件,根据能量变化得到晶胞密度为1.2 g/cm3,总能量为1 140.624 kJ/mol,与实际情况基本相符,验证了该建模方法的实效性。该研究为直观认识煤大分子结构提供了方法,并对煤与瓦斯突出、煤自燃等灾害事故的机理研究及防治具有重要意义。Abstract: In view of the current situation that there are many but no uniform coal molecular modeling methods, this paper explores the structuve of lignite macromolecules and optimizes them from a microscopic perspective by selecting physical methods that have less impact on molecular structure.A relatively simple coal macromolecule modeling method was put forward.Based on the experimental methods of elemental analysis, 13C-NMR and XPS, molecular dynamics simulation software method was used to analyze and study lignite and build molecular structure models.The results show that the aromatic structure of lignite macromolecules is mostly pentaphenyl.The fat carbon structure is mostly methylene and methine, and the alkane chain is mostly cycloalkanes.Oxygen atoms are mostly ether bond oxygen, followed by carboxyl oxygen and carbonyl oxygen.Nitrogen atom exist in the form of Pyridinic nitrogen N-6.After the optimization, the structure of the coal molecular model is more compact and the energy is obviously reduced.As the main component of the non-bonding potential energy and the main factor to maintain the stability of the coal molecular structure, the van der Waals potential energy has the most obvious change.In order to construct coal molecular cell, periodic boundary conditions were added.According to the energy variation, the cell density was 1.2 g/cm3 and the total energy was 1 140.624 kJ/mol, which was basically consistent with the actual situation and verified the effectiveness of the modeling method.This study provides a method for the direct understanding of the macromolecular structure of coal, and is of great significance for the mechanism research and prevention of coal and gas outburst, coal spontaneous combustion and other disasters.
-
Key words:
- lignite /
- macromolecular modeling /
- molecular dynamics /
- structure optimization
-
表 1 13C-NMR谱图化学位移的结构归属
Table 1. Structural attribution of chemical shifts in 13C-NMR spectra
碳结构名称 主要归属 典型化学位移/10-6 碳结构名称 主要归属 典型化学位移/10-6 脂甲基碳 14~16 质子化芳碳 100~129 芳环甲基碳 16~22 桥头芳香碳 129~137 亚甲基碳、次甲基碳 22~36 侧枝芳香碳 137~148 次甲基碳、季碳 36~50 氧取代芳香碳 148~165 氧接甲基及氧接亚甲基碳 50~56 羧基碳 165~188 氧接次甲基碳 60~70 醛、醌、苯基酮的羰基碳 188~205 氧接季碳 75~90 链烷酮、环烷酮的羰基碳 205~220 表 2 煤样元素构成
Table 2. Element composition of coal sample
质量/mg C/% H/% O/% N/% S/% 1.653 58.88 4.36 25.38 0.65 0.21 表 3 结构归属及相对面积
Table 3. Structure and relative area
结构归属 结构比例/% 质子化芳碳 22.161 桥头芳香碳 19.698 亚甲基碳、次甲基碳 19.255 次甲基碳、季碳 11.987 侧枝芳香碳 10.166 氧取代芳香碳 7.151 氧接甲基及氧接亚甲基碳 4.697 芳环甲基碳 1.455 链烷酮、环烷酮的羰基碳 1.446 氧接次甲基碳 1.130 脂甲基碳 0.771 醛、醌、苯基酮的羰基碳 0.082 表 4 褐煤13C-NMR结构参数
Table 4. 13C-NMR structural parameters in lignite
fa faC f′a faN faH faP 60.74 1.53 59.18 37.02 22.16 7.15 faS faB fal fal* falH falO 10.17 19.7 39.29 14.21 19.26 5.83 表 5 优化后褐煤大分子结构能量
Table 5. Energy of macromolecular structure of lignite after optimization
kJ/mol 总能量 非成键能 价电子能 氢键能 范德华能 库仑能 键伸缩能 键角能 扭转能 反转能 优化前 45 036.044 0.000 31 114.270 -53.144 12 462.340 364.826 1 126.894 20.862 优化后 3 680.514 0.000 1 939.305 -141.607 487.150 522.143 852.746 20.779 表 6 密度为1.2 g/cm3褐煤大分子结构能量
Table 6. Energy of macromolecular structure for lignite after optimization under 1.2 g/cm3
密度/(g·cm-3) 总能量/(kJ·mol-1) 非成键能/(kJ·mol-1) 价电子能/(kJ·mol-1) 氢键能 范德华能 氢键能 范德华能 氢键能 范德华能 氢键能 1.2 1 140.624 0 2 090.887 -3 353.010 443.876 833.981 1 058.200 66.693 -
[1] 中国煤炭工业协会. 中华人民共和国推荐性国家标准: 中国煤炭分类GB/T 5751—2009[S]. 北京: 中国标准出版社, 2010. [2] 林雄超, 王彩红, 田斌, 等. 脱灰对两种烟煤半焦碳结构及CO2气化反应性的影响[J]. 中国矿业大学学报, 2013, 42(6): 1040-1046. doi: 10.3969/j.issn.1000-1964.2013.06.023Lin Xiongchao, Wang Caihong, Tian Bin, et al. Effects of de-ashing on the micro-structural transformation and CO2 reactivity of two Chinese bituminous coal chars[J]. Journal of China University of Mining & Technology, 2013, 42(6): 1040-1046. doi: 10.3969/j.issn.1000-1964.2013.06.023 [3] 孙成功, 李保庆, Colin E·Snape. 煤中有机硫形态结构和热解过程硫变迁特性的研究[J]. 燃料化学学报, 1997, 25(4): 71-75. https://www.cnki.com.cn/Article/CJFDTOTAL-RLHX704.013.htmSun Chenggong, Li Baoqing, Colin E·Snape. Characteri-zation of organic sulfur forms in some Chinese coals by high pressure tpr and sulphur transfer during hydropyrolysis[J]. Journal of Fuel Chemistry and Technology, 1997, 25(4): 71-75. https://www.cnki.com.cn/Article/CJFDTOTAL-RLHX704.013.htm [4] 戴中蜀, 郑昀辉. 低煤化度煤经低温热解后各基团变化的研究[J]. 煤炭转化, 1997, 20(1): 54-58. https://www.cnki.com.cn/Article/CJFDTOTAL-MTZH199701009.htmDai Zhongshu, Zheng Yunhui. The study of groups variation of low rank coal during low temperature pyrolysis by FTIR spectrography[J]. Coal Conversion, 1997, 20(1): 54-58. https://www.cnki.com.cn/Article/CJFDTOTAL-MTZH199701009.htm [5] 吴国光, 王祖讷. 低阶煤的热重-傅里叶变换红外光谱的研究[J]. 中国矿业大学学报, 1998, 27(2): 72-75. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGKD802.016.htmWu Guoguang, Wang Zuna. Research on TG-FTIR of low rank coals[J]. Journal of China University of Mining & Technology, 1998, 27(2): 72-75. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGKD802.016.htm [6] 何满潮, 韩宗芳, 杨华. 不同温度下高岭石变形及破坏机理的分子动力学模拟[J]. 矿业科学学报, 2019, 4(1): 8-16. http://kykxxb.cumtb.edu.cn/article/id/191He Manchao, Han Zongfang, Yang Hua. Molecular dynamics simulation of deformation and failure mechanism of kaolinite at different temperatures[J]. Journal of Mining Science and Technology, 2019, 4(1): 8-16. http://kykxxb.cumtb.edu.cn/article/id/191 [7] Chermin H A G, Van K D W. Chemical structure and properties of coal. XVⅡ-A mathematical model of coal pyrolysis[J]. Fuel, 1957, 36(1): 85-104. http://www.researchgate.net/publication/284060154_Chemical_structure_and_properties_of_coal-XVII_A_mathematical_model_of_coal_pyrolysis [8] Given P H. The distribution of hydrogen in coals and its relation to coal structure[J]. Fuel, 1960, 39(2): 147. [9] Wiser W H. Reported in division of fuel chemistry[J]. Preprints, 1971, 20(1): 122. [10] Shinn J H. From coal to single-stage and two-stage products: a reactive model of coal structure[J]. Fuel, 1984, 63(9): 1187-1196. doi: 10.1016/0016-2361(84)90422-8 [11] Lian Lulu, Qin Zhihong, Li Chunsheng, et al. Molecular model construction of the dense medium component scaffold in coal for molecular aggregate simulation[J]. ACS Omega, 2020, 5(22): 13375-13383. doi: 10.1021/acsomega.0c01575 [12] Wu Dun, Zhang Hui, Hu Guangqing, et al. Fine characterization of the macromolecular structure of huainan coal using XRD, FTIR, 13C-CP/MAS NMR, SEM, and AFM techniques[J]. Molecules, 2020, 25(11): 2661. doi: 10.3390/molecules25112661 [13] 洪迪昆. 准东煤热解及富氧燃烧的反应分子动力学研究[D]. 武汉: 华中科技大学, 2018. [14] 魏帅, 严国超, 张志强, 等. 晋城无烟煤的分子结构特征分析[J]. 煤炭学报, 2018, 43(2): 555-562. https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB201802031.htmWei Shuai, Yan Guochao, Zhang Zhiqiang, et al. Molecular structure analysis of Jincheng anthracite coal[J]. Journal of China Coal Society, 2018, 43(2): 555-562. https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB201802031.htm [15] Wu Siyuan, Jin Zhixin, Deng Cunbao. Molecular simulation of coal-fired plant flue gas competitive adsorption and diffusion on coal[J]. Fuel, 2019, 239: 97-96. doi: 10.1016/j.fuel.2018.10.146 [16] 聂尧, 赵越超. 煤中多组分混合气体竞争吸附研究现状及工程应用[J]. 矿业科学学报, 2020, 5(1): 45-57. http://kykxxb.cumtb.edu.cn/article/id/264Nie Yao, Zhao Yuechao. Research status and engineering application of competitive adsorption of multicomponent mixed gases in coal[J]. Journal of Mining Science and Technology, 2020, 5(1): 45-57. http://kykxxb.cumtb.edu.cn/article/id/264 [17] 王擎, 王智超, 贾春霞, 等. 基于固体13C核磁共振技术对油砂沥青质结构的研究[J]. 化工进展, 2014, 33(6): 1392-1396. https://www.cnki.com.cn/Article/CJFDTOTAL-HGJZ201406007.htmWang Qing, Wang Zhichao, Jia Chunxia, et al. Study on structural features of oil sands with solid state 13C NMR[J]. Chemical Industry and Engineering Progress, 2014, 33(6): 1392-1396. https://www.cnki.com.cn/Article/CJFDTOTAL-HGJZ201406007.htm [18] Solum M S, Pugmire R J, Grant D M, et al. N-15 CPMAS NMR of the argonne premium coals[J]. Energy & Fuels, 1997, 11(2): 491-494. http://www.sciencedirect.com/science/article/pii/S0140670197844219 [19] 郑昀辉, 戴中蜀. 用NMR研究低温热处理对低煤化度煤化学组成结构的影响[J]. 煤炭转化, 1997, 20(4): 54-59. https://www.cnki.com.cn/Article/CJFDTOTAL-MTZH199704009.htmZheng Yunhui, Dai Zhongshu. Using NMR to research the influence of low temperature plyolysis on the chemical component and structure of low rank coal[J]. Coal Conversion, 1997, 20(4): 54-59. https://www.cnki.com.cn/Article/CJFDTOTAL-MTZH199704009.htm [20] Song Y, Jiang B, Li J H. Retraction note to: simulations and experimental investigations of the competitive adsorption of CH4 and CO2 on low-rank coal vitrinite[J]. Journal of Molecular Modeling, 2019, 25(6): 178. doi: 10.1007/s00894-019-4074-8 [21] Song Yu, Jiang Bo, Li Wu. Molecular simulation of CH4/CO2/H2O competitive adsorption on low rank coal vitrinite[J]. Physical Chemistry Chemical Physics, 2017, 19(27): 17773-17788. doi: 10.1039/C7CP02993D [22] Xiang J H, Zeng F G, Liang H Z, et al. Molecular simulation of the CH4/CO2/H2O adsorption onto the molecular structure of coal[J]. Science China Earth Sciences, 2014, 57(8): 1749-1759. doi: 10.1007/s11430-014-4849-9 [23] 张洋, 吴兵, 李滕滕, 等. 二次氧化煤燃烧和灭火过程的实验研究[J]. 矿业科学学报, 2020, 5(3): 272-277. http://kykxxb.cumtb.edu.cn/article/id/289Zhang Yang, Wu Bing, Li Tengteng, et al. Experimental study on combustion and extinguishing of secondary oxidation coal[J]. Journal of Mining Science and Technology, 2020, 5(3): 272-277. http://kykxxb.cumtb.edu.cn/article/id/289 [24] Takanohashi T, Nakamura K, Terao Y, et al. Computer simulation of solvent swelling of coal molecules: Effect of different solvents[J]. Energy & Fuels, 2000, 14(2): 393-399. doi: 10.1021/ef990147f [25] Carlson G A. Computer simulation of the molecular structure of bituminous coal[J]. Energy & Fuels, 1992, 6(6): 771-778. doi: 10.1021/ef00036a012