0 引 言
在我国能源体系中，煤炭是储量最大的不可再生能源，同时也是重要的基础能源和工业原料［1-2］。正确认识煤的分子结构特征对实现煤炭合理利用及高效转化具有重大意义。 许多现代分析测试技术可用于研究煤的结构特征，如傅里叶红外光谱（FTIR）、X 射线衍射（XRD）、X 射线光电子能谱（XPS）、拉曼光谱、核磁共振碳谱（13 C-NMR）和高分辨率透射镜（HRTEM）等。 近半个世纪以来，国内外煤炭研究学者们致力研究煤的结构与性质，发现煤中有大量的亲水基团和缩合芳环，且随着变质程度的升高，煤中芳香环数增加而脂肪结构不断减少［3-5］。 同时，学者们利用化学试验研究和现代测试技术对煤的分子结构进行研究，提出了130 多种煤的大分子结构模型［6］，例如，著名的Wiser 模型［7］、Given 模型［8］、Shinn 模型［9］等。 这些结构反映出不同变质阶段的煤分子结构，在一定程度上解释了煤的浮选、溶胀裂解以及气体吸附等行为［10-13］。 随着近几年计算机技术的快速发展，在研究复杂的分子结构及体系时，计算化学对试验的辅助及应用越来越成熟［14-15］。 构建煤的分子模型是对其进行分子模拟的基础，通过分子模拟计算获得煤的合理分子结构可为进一步选择浮选药剂提供理论基础。
以淮北矿业集团青东煤矿煤为研究对象，通过工业分析、元素分析、13C-NMR、XPS 等测试对其分子结构进行分析，建立青东煤的单个大分子结构模型及聚集态结构模型。 为寻找单个煤分子最低能量构型，利用分子模拟软件Materials studio（MS）7.0 对单个分子模型结构进行几何优化及退火动力学模拟研究，并构建出煤的聚集态结构模型，从分子尺度为研究青东煤的浮选药剂吸附机理提供了模型基础。
1 样品及试验方法
1.1 样品采集和制备
煤样采自于淮北矿业集团青东煤矿，采样方法遵循GB/T 482—2008《煤层煤样采取方式》。 在实验室内，按GB/T 478—2008《煤炭浮沉试验方法》对青东煤进行煤岩分离，选取出密度级为-1.30 g/cm3的煤样在70 ℃的条件下干燥至恒重。
1.2 工业分析和元素分析
青东煤样的工业分析和元素分析见表1。 采用德国Vario EL cube 型元素分析仪测定样品的C、H、S、N 元素含量，采用差减法计算O 元素含量。 并根据工业分析和元素分析结果计算青东煤样的氢碳原子比H/C、氧碳原子比O/C、氮碳原子比N/C 和硫碳原子比S/C（表1）。
表1 青东煤样的工业分析和元素分析
Table1 Proximate andelemental analysis of Qingdong coal sample
元素分析/%工业分析/%Mad Aad Vdaf H/C O/C N/C S/C w（C） w（H） w（O） w（N） w（S）0.53 1.88 25.68 90.44 4.96 2.91 1.43 0.26 0.66 0.026 0.016 0.001
1.3 核磁共振碳谱测试
采用Bruker AVANCE III600 型核磁共振仪对青东煤样进行13C-NMR 分析。 测试参数如下：外径为4 mm 的ZrO2探头，核磁共振频率75.47 MHz，转速7.0 kHz，脉宽4.2×10-6 s，采样时间0.05 s，碳氢交叉极化接触时间3 ms，谱宽3 000 Hz，数据采集累加次数6 000 次。
1.4 X 射线光电子能谱测试
采用美国热电Thermo Escalab 250XI 型光电子能谱仪对青东煤样进行XPS 测试。 测试参数如下：单色Al Ka （hv ＝1 486.6 eV），功率150 W，X 射线束斑500 μm，采用污染碳（284.8 eV）对电荷进行校正，通能窄扫20 eV，宽扫100 eV，真空度为10-8 mPa。
2 讨论与结果
2.1 核磁共振碳谱分析
煤的13C-NMR 谱图可划分为3 个区域［16-17］：①化学位移在0 ～80 的脂碳区；②化学位移在80 ～170 的芳碳区；③化学位移在200 左右的羰基碳区。由于煤中官能团种类丰富且结构复杂，13C-NMR 谱图上的峰可由多个子峰叠加形成，因此需要对谱图进行分峰处理，以确定各个官能团的峰位及相对含量［18-20］。 利用Origin8.0 软件对青东煤样的13 CNMR 谱图进行分峰拟合，结果如图1 所示。
图1 青东煤样的13C-NMR 分峰拟合
Fig.1 13C-NMR spectra of Qingdong coal
根据青东煤中各官能团的峰位及其相对百分含量，计算出煤样的12 个重要结构参数。 由于受核磁测试试验边带效应的影响，青东煤样的羰基碳区整体移至215 左右，作用强度比实际偏大［21］。 研究发现，煤化程度越高，羰基峰强度越偏离实际强度，将对煤炭芳香率与脂碳率的计算产生不可忽视的误差。 为降低该种误差，对青东煤样的芳香率f′αl 与脂碳率fαl 进行修正［22-23］，结果见表2。
根据表2 中的结构参数百分比，计算青东煤样分子结构中芳香化合物的桥碳与周碳比（XBP）：
＝0.35，该参数反映了煤的分子结构中芳香化合物缩聚程度的平均值，可用来计算芳香簇尺寸［24-25］。
表2 青东煤样的结构参数
Table 2 Structure parameters of Qingdong coal
fα fC α f′α fH α fN α fP α fS α fB α fαl f*αl fHαl fO αl 66.94 2.22 64.72 43.15 19.80 0.01 3.51 16.28 33.06 11.00 22.06 1.80
注：fα 为sp2 杂化碳；fαl 为sp3 杂化碳；fCα 为与羰基相连的芳香碳；f′α 为芳香碳；fHα 为质子化芳碳；fNα 为非质子化碳；fPα 为与氧原子相连的芳碳；fSα 为被烷基取代的芳碳；fBα 为桥接芳碳；f*αl 为脂肪烃中的甲基；fHαl 为亚甲基、季碳或次亚甲基；fOαl 为与氧原子相连的脂肪碳。
2.2 X 射线光电子能谱分析
煤样中碳、氧、氮、硫以多种官能团形式存在，其XPS 谱图经多个子峰叠加形成。 采用Case XPS 软件对4 个XPS 谱图进行分峰拟合，依据各结合能的归属，确定它们在样品中的存在形式［26-29］。 青东煤的碳原子在结构中以4 种形态（图2 和表3）存在：①284.63 eV 峰归属于芳香石墨化碳（C ═C）；②285.07 eV峰归属于脂肪族碳（C—H）；③285.46 eV 峰归属于羟基结构（C—O）和羰基结构（C ═O）。 计 算 子 峰 的 相 对 面 积， C ═C 结 构 碳 占65.15%，即青东煤样中的碳主要以芳香结构碳为主，与核磁共振碳谱图结果一致。 氧原子在青东煤表面结构中以2 种形态（图3 和表4）存在：①532.08 eV 峰为羟基氧（C—O）或醚氧（—O—），子峰相对面积45.39%；②533.20 eV 和533.94 eV 峰均为羰基氧（C ═O），子峰相对面积54.61%。 由此可知，青东煤样中的氧原子以羰基和酚羟基2 种形式存在。 氮原子在青东煤中的存在形式为吡啶和吡咯（图4 和表5）：①398.64 eV 峰为吡啶型氮（N—6），相对面积为59.48%；②400.49 eV 峰为吡咯型氮（N—5），相对面积40.52%。 硫原子在煤样中以亚砜型硫化物和噻吩型硫化物形态存在。 165.28 eV 峰为亚砜型硫化物，164.16 eV 峰为噻吩型硫化物（图5 和表6）。
图2 青东煤样的XPS C（1s）谱图
Fig.2 XPS C（1s） spectra of Qingdong coal
表3 青东煤样的XPS C（1s）数据
Table 3 XPS C（1s） data of Qingdong coal
284.630 C ═C 65.15 285.070 C—H 22.60 285.467 C—O， C ═O 12.25
图3 青东煤样的XPS O（1s）谱图
Fig.3 XPS O（1s） spectra of Qingdong coal
表4 青东煤样的XPS O（1s）数据
Table 4 XPS O（1s） data of Qingdong coal
峰位/eV 归属 相对面积/%533.94 C ═O， COO— 22.76 533.20 C ═O 31.85 532.08 C—O 45.39
图4 青东煤样的XPS N（1s）谱
Fig.4 XPS N（1s） spectra of Qingdong coal
图5 青东煤样的XPS S（2p）谱
Fig.5 XPS S（2p） spectra of Qingdong coal
表5 青东煤样的XPS N（1s）数据
Table 5 XPS N（1s） data of Qingdong coal
400.49 N-5 59.48 398.64 N-6 40.52
表6 青东煤样的XPS S（2p）数据
Table 6 XPS S（2p） data of Qingdong coal
峰位/eV 归属 相对面积/%165.28 亚砜型硫 29.35 164.16 噻吩型硫 70.65
3 青东煤分子结构及模型构建
3.1 芳香结构
文献［30-31］研究表明，当煤中碳质量分数为81%～91%时，芳香环平均环数为2 ～5 个。 青东煤样的碳质量分数为90.44%，XBP 为0.35，2 环芳香化合物萘的XBP为0.25，3 环的芳香化合物蒽的XBP 为0.4［32］，通过调整芳香环的类型和数量，使结构模型中的XBP 接近0.35。 当青东煤结构中萘的数量为2个，蒽为4 个时，结构中的XBP ＝0.357。 基于ACD/CNMR Predictor 软件只能用于不超过255 个原子（不包括H 原子数目）的分子，因此模型的芳香结构类型和数量见表7， 其中芳香碳原子数量为97 个［33］。
表7 青东煤样化学结构模型中芳香结构单元的类型
Table 7 Types of aromatic structure unit in chemical structural model of Qingdong coal
芳香结构 数量/个 芳香结构 数量/个images/BZ_253_1433_972_1522_1075.png2images/BZ_253_1911_968_2000_1078.png1images/BZ_253_1389_1171_1566_1274.png2images/BZ_253_1910_1163_2001_1282.png1images/BZ_253_1344_1370_1611_1473.png3images/BZ_253_1828_1354_2083_1488.png1
3.2 脂肪结构
亚甲基、次亚甲基、脂环结构和烷基侧链是煤中脂肪结构的主要存在形式［34］。 研究表明，烷基侧链的长度会随煤化程度的增加而减小［35］。 当碳质量分数为90.4%时，烷基侧链中原子的平均数量为1 ～2 个，说明在青东煤结构中烷基侧链不宜过长，以短链为主［36］。 由表2 可知，
，
，说明煤结构中亚甲基、季碳和次甲基的相对百分含量大于甲基的相对百分含量。 由煤样的芳香度为64.72%及模型中芳香原子数为97，可初步确定脂肪碳原子数为53。
3.3 杂原子结构
由上可初步确定青东煤样结构的碳原子总数为150，结合元素分析中各元素与碳元素原子比可推算出煤样模型中氧、氮、硫原子数量分别为3、2、0。 研究表明，煤中含氧官能团包括甲氧基、羧基、羟基和羰基等［37］。 随着煤化程度的增高，甲氧基含量迅速减少，在硬褐煤中基本消失；羧基主要存在于褐煤中，而羟基和羰基存在于整个烟煤阶段［38］。 结合13C-NMR 及XPS 分析结果可确定青东煤结构模型中3 个氧原子的存在形式为2 个羰基和1 个酚羟基。
吡啶型氮和吡咯型氮是煤结构中氮的主要存在形式［39］。 结合XPS 结果可知，青东模型中的2 个氮原子分别以吡啶型氮和吡咯型氮的形式存在。 XPS结果还表明，煤样中的有机硫包括亚砜型和噻吩型硫化物。 由于该煤样中硫质量分数仅为0.26%，计算得出的硫原子个数不足1 个，且模型中原子数量有限，因此构建模型时不再考虑硫原子。
3.4 青东煤分子结构模型
采用分子结构绘图软件Chem draw/Chem3D 构建青东煤平面结构模型。 应用ACD/CNMR Predictor软件计算出模型结构的化学位移，将数据导入gNMR 软件获得模型的计算13 C-NMR 谱图［40］。 利用Origin8.0 软件将模型的计算13C-NMR 谱图与试验测试谱图数据进行匹配，通过调整模型中的各个结构单元的连接方式，对结构模型进行修正［41］。 青东煤结构模型计算13 C-NMR 谱图和试验13 C-NMR谱图对比效果如图6 所示。
图6 青东煤的结构模型计算13C-NMR和试验13C-NMR 谱
Fig.6 Experimental and calculated13C-NMR spectrum of Qingdong coal
由图6 可知，模型计算谱图和试验谱图整体吻合较好，但由于在核磁共振碳谱测试试验过程中存在边带效应，使实验谱图中的碳氧区吸收峰强度偏大，因此在化学位移200～230 区域计算谱图吸收峰强度明显低于试验谱图。 青东煤分子结构的分子式为C142H128N2O3，其分子量为1 910.60。 碳、氮及氧元素质量分数分别为89.27%、1.47%、2.51%，与元素分析中各元素相对百分含量接近；而结构中的氢质量分数为6.75%，与元素分析中的4.96%偏差较大，这是因为在构建煤结构时未考虑硫元素造成的比例偏差。
表8 青东煤样结构模型的结构参数
Table 8 Structure parameters of macro-molecular structure model of Qingdong coal
分子式 分子量元素质量分数/%w（C） w（H） w（O） w（N） w（S）C142H128N2O3 1 910.60 89.27 6.75 2.51 1.47 0
综上所述，构建的青东煤分子平均结构模型结果与元素分析、XPS 分析及13 C-NMR 分析结果吻合度较好，可反映出青东煤的平均分子结构特征。
4 单个分子结构模型的分子模拟
4.1 单个分子结构模型
为确定单分子结构模型的最低能量构型，利用MS 软件对青东煤单个分子结构模型进行几何优化及退火动力学模拟。 将青东煤大分子结构（图7）导入MS 软件，使其结构饱和加氢。 几何优化采用Forcite 模块，优化方法为Smart Minimizer 方法。 参数设置如下：基于原子的总能量计算库仑能和范德华能，原子的净电荷采用电荷平衡法获得［42］，迭代步数为5 000 步，收敛标准采用Fine。 几何优化后，继续采用Forcite 模块的Anneal 进行高温弛豫。 模拟过程在NVT 系综下进行，选择Nose 控温法；初始温度设为300 K，最高温度设为800 K；升温速率为60 K/次，模拟时间为1 fs；每次循环结束后对输出结果再次进行分子力学优化，其设置参数如前所述。该力场计算定义［43］为
其中：Ebond为键伸缩能；Ecross为交叉项相互作用能；Etorsion为二面角扭转能；Eangle 为键角能；Eopp 离平面相互作用项；Elj 为范德华能；Eelect 为静电相互作用；EH为氢键作用能。 其中Ebond、Ecross、Etorsion、Eangle、Eopp为价电子能；Elj、Eelect、EH为非成键能。
选取模拟结果中能量最低的结构作为青东煤的最优几何构型，如图8 所示。 可以看出，经过分子力学和分子动力学优化后，单个分子结构为达到空间结构上官能团之间的最小斥力，部分桥键、脂肪键等化学键发生了明显的扭转［44］。 由于芳香环之间π-π相互作用，相邻芳香片层之间趋于平行排列，层间距增大，呈现出明显的立体感［45］。
表9 是青东煤结构模型优化前后的能量组成，优化前总能量为2 121.14 kJ/mol，优化后为1 255.85 kJ/mol，总能量大幅度降低。 优化后的模型中，键伸缩能及范德华能降低，其他项均升高，而总能量明显降低，说明该结构中的键伸缩能及范德华能占主导地位，其中键伸缩能属于价电子能，范德华能属于非成键能。 由于煤是由多个芳环组成的大分子结构，在模拟的过程中，青东煤分子模型由二维的平面结构转化为三维的立体结构，分子结构中部分化学键扭转，键角变化，造成键角能及扭转能增大，同时引起键伸缩势能的降低［46］。 优化后的立体结构中，芳香片层之间的平行排列产生较大的色散力，是范德华能下降的主要原因［47］。
图7 青东煤平面分子结构模型
Fig.7 Plane model molecular structure in Qingdong coal
图8 单个煤分子的最低能量结构模型
Fig.8 Energy-minimum conformation of a single coal model
表9 青东煤结构模型的能量组成 kJ/mol
Table 9 Potential energy for conformation of Qingdong coal
项目 总能量 Ebond Eangle Etorsion Elj Eelect初始 2121.14 931.43 73.92 517.04 634.14 -36.596最终 1255.85 208.40 268.86 632.37 75.25 47.26
4.2 聚集态结构模型
利用MS 软件中的Amorphous Cell（AC）模块将18 个优化后的青东煤大分子随机添加到尺寸为40Å×40Å×40Å（1Å ＝0.1 m）的晶胞中，并添加三维周期性边界条件［48］。 设定结构模型的密度为1.25 g/cm3。 首先利用分子力学进行几何优化，然后进行退火动力学模拟对该体系进行弛豫，温度设置为298～1 098 K，并使最终温度保持在298 K［49-50］。 选定能量最低的构型体系作为周期性边界条件下青东煤分子聚集态结构模型，如图9 所示。 从聚集态结构模型中可以看出，经过分子力学和分子动力学模拟优化后，分子结构发生明显弯曲、扭转。 由于煤是非晶体物质，无长程周期性，受周围分子的制约，原本近似平香碳结构发生扭曲形变，片层结构行排列的片状芳杂乱，使整个凝聚态结构模型更加紧凑［51］。 从图9 知青东煤大分子分布较为均匀，表面含氧官能团分布均匀，可反映出青东煤样表面的结构特征。
图9 青东煤聚集态结构模型
Fig.9 Aggregation-state for Qingdong coal molecules
5 结 论
1）根据青东煤的核磁共振碳谱得出表征煤结构的12 个结构参数，并进行修正；根据参数计算出煤分子结构中芳香桥碳与周碳比为0.35。13C-NMR及XPS 结果表明，青东煤的碳原子在结构中主要以芳香石墨化碳结构为主，其XPS 谱图子峰相对面积占65.15%；煤样中芳香结构单元包括2 个苯环，2个萘，4 个蒽；煤分子结构中的杂原子别以2 个羰基和1 个酚羟基、1 个吡啶和1 个吡咯存在。
2） 利 用Chem draw/Chem3D 及ACD/CNMR Predictor 构建了与青东煤的13 C-NMR 试验谱图拟合较好的结构模型；最终构建的青东煤单分子结构分子式为C142H128N2O3，分子量为1 910.60。
3）利用MS 软件对单个分子结构模型进行几何优化和分子动力学模拟后，桥键、脂肪键等化学键发生了扭转，分子内芳香片层之间的π-π 相互作用使相邻芳香片层之间趋于近似平行排列；优化前后分子总势能由2 121.14 kJ/mol 下降到1 255.85 kJ/mol，键伸缩能及范德华能占主导地位。
4）构建出聚集态结构模型，经过分子力学和分子动力学模拟优化后，聚集态结构模型中的大分子受周围分子的制约，原本近似平行排列的片状芳香碳结构发生扭曲形变，片层结构杂乱，整个凝聚态结构模型紧凑；构建的聚集态结构模型可为后续探索浮选药剂吸附模拟提供模型基础。
［1］ 许 莉.世界煤炭资源供需形势分析［J］.中国煤炭地质，2012， 24（6）：74-76.XU Li.World coal resource supply and demand analysis［J］.Coal Geology of China， 2012， 24（6）：74-76.
［2］ 许家林.煤矿绿色开采20 年研究及进展［J］.煤炭科学技术，2020， 48（9）：1-15.XU Jialin.Research and progress of coal mine green mining in 20 years［J］.Coal Science and Technology， 2020， 48（9）：1-15.
［3］ OKOLOG N， NEOMAGUS H W J P， EVERSON R C， et al.Chemical-structural properties of South African bituminous coals：insights from wide angle XRD-carbon fraction analysis，ATR-FT⁃IR， solid state 13C-NMR， and HRTEM techniques［J］.Fuel，2015， 158：779-792.
［4］ 张小东， 张 鹏.不同煤级煤分级萃取后的XRD 结构特征及其演化机理［J］.煤炭学报， 2014， 39（5）：941-946.ZHANG Xiaodong，ZHANG Peng.Characteristics of XRD parameter for different ranks of coals under fractional extraction and its evolu⁃tion mechanism［J］.Journal of China Coal Society，2014，39（5）：941-946.
［5］ 罗陨飞， 李文华.中低变质程度煤显微组分大分子结构的XRD研究［J］.煤炭学报， 2004， 29（3）：339-341.LUO Yunfei， LI Wenhua.X-ray diffraction analysis on the different macerals of several low-to-medium metamorpic grade coals［J］.Journal of China Coal Society， 2004， 29（3）：339-341.
［6］ MATHEWS J P，DUIN A C T V， CHAFFEE A L.The utility of coal molecular models［J］.Fuel Processing Technology， 2011， 92（4）：718-728.
［7］ WISER W H.Conversion of bituminous coal to liquids and gases：chemistry and representative processes.Magnetic Resonance［M］.Netherlands：Springer， 1984.
［8］ GIVEN P H.Structure of bituminous coals：evidence from distribu⁃tion of hydrogen［J］.Nature， 1959， 184：980-981.
［9］ SHINN J H.Visualization of complex hydrocarbon reaction systems［J］.Preprints of Papers American Chemical Society Division of Fuel Chemistry， 1996， 41（6）：418.
［10］ TAKANOHASHUI 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.
［11］ ZHENG M， LI X， LIU J，et al.Initial chemical reaction simu⁃lation of coal pyrolysis via ReaxFF molecular dynamics［J］.Energy ＆ Fuels， 2013， 27（6）：2942-2951.
［12］ ZHANG Z， YAN K.Molecular dynamics simulation of oxygen dif⁃fusion in dry and water-containing brown coal［J］.Molecular Physics， 2011， 109（19）：2367-2374.
［13］ ZHANG Z， WANG C， YAN K.Adsorption of collectors on model surface of Wiser bituminous coal：a molecular dynamics simulation study［J］.Minerals Engineering， 2015， 79：31-39.
［14］ 王英学，梁俊玺，胡勤政，等.计算化学与分子模拟方法及MS软件简介［J］.山东化工， 2019， 48（02）：131-132.WANG Yingxue， LIANG Junxi， HU Qinzheng et al.Review on the research of different composite coatings ［ J］. Shandong Chemical Industry， 2019， 48（02）：131-132.
［15］ 文玉华， 朱如曾， 周富信， 等.分子动力学模拟的主要技术［J］.力学进展， 2003， 33（1）：65-73.WENG Yuhua， ZHU Ruzeng， ZHOU Fuxin， et al.Main tech⁃niques of molecular dynamics simulation ［J］.Advances in Me⁃chanics， 2003， 33（1）：65-73.
［16］ 徐秀峰， 张蓬洲.高分辨固体 13C-NMR 和XPS 技术表征碳的骨架结构［J］.煤炭转化， 1995， 18（4）：57-62.XU Xiufeng， ZHANG Pengzhou.The study of carbon structure by solid 13C-NMR and XPS［J］.Coal Conversion， 1995， 18（4）：57-62.
［17］ 叶超辉， 李新安.煤的固体高分辨13C-NMR 谱［J］.科学通报，1985， 30（20）：1545-1547.YE Chaohui， LI Xinan.High resolution of solid state13 C-NMR spectrum in coal ［J］. Chinese Sci Bull， 1985， 30 （20）：1545-1547.
［18］ 彭立才， 韩德馨， 邵文斌， 等.柴达木盆地北缘侏罗系烃源岩干酪根13 C-核磁共振研究［J］.石油学报， 2002， 23（2）：34-37. PENG Licai， HAN Dexin， SHAO Wenbin， et al.13 C-NMR re⁃search on the kerogens of Jurassic hydrocabon source rock in the Northen Edge， Qaidam basin［J］.Acta Petrolei Sinica， 2002， 23（2）：34-37.
［19］ 郑昀辉， 戴中蜀.用NMR 研究低温热处理对低煤化度煤化学组成结构的影响［J］.煤炭转化， 1997， 20（4）：54-59.ZHENG Yunhui， DAI Zhangshu.Using NMR to research the in⁃fluence of low temperature heat treatment on the chemical compo⁃nent and structure of low rank coal［J］.Coal Conversion， 1997，20（4）：54-59.
［20］ TREWHELLA M J， POPLETT I J F， GRINT A.Structure of green river oil shalekerogen：determination using solid state 13CNMR.spectroscopy［J］.Fuel， 1986， 65（4）：541-546.
［21］ 贾建波.神东煤镜质组结构模型的构建及其热解甲烷生成机理的分子模拟［D］.太原：太原理工大学， 2010.
［22］ 贾建波， 曾凡桂， 孙蓓蕾.神东2-2煤镜质组大分子结构模型13 C-NMR 谱的构建与修正［J］.燃料化学学报， 2011， 39（9）：652-657.JIA Jianbo， ZENG Fanhui， SUN Beilei.Construction and modifi⁃cation of macromolecular structure model for vitrinite from Shendong 2-2 coal［J］.Journal of Fuel Chemistry and Tech⁃nology， 2011， 39（9）：652-657.
［23］ 孙蓓蕾.13C-NMR 谱的构建与修正［J］.燃料化学学报， 2011，39（9）：652-657.SUN Beilei.Construction and modification of 13C-NMR spectra［J］.Journal of Fuel Chemistry and Technology， 2011， 39（9）：652-657.
［24］ 王 丽， 张蓬洲.用固体核磁共振和电子能谱研究我国高硫煤的结构［J］.燃料化学学报， 1996， 24（6）：539-543.WANG Li， ZHANG Pengzhou. Study on structural characterization of three Chinese coals of high organic sulphur content using XPS and solid state NMR spectroscopy［J］.Journal of Fuel Chemistry and Technology， 1996， 24（6）：539-543.
［25］ 李鹏鹏.杜儿坪2 号煤结构模型构建及其分子模拟［D］.太原：太原理工大学， 2014.
［26］ KOZLOWSKI M.XPS study of reductively and non-reductively modified coals［J］.Fuel， 2004， 83（3）：259-265.
［27］ GRZYBEK T， PIETRZAK R， WACHOWSKA H.X-ray photoe⁃lectron spectroscopy study of oxidized coals with differentsulphur content［J］.Fuel Processing Technology， 2002， 77（1）：1-7.
［28］ GARDNER S D， SINGAMSETTY C S K， BOOTH G L， et al.Surface characterization of carbon fibers using angle-resolved XPS and ISS［J］.Carbon， 1995， 33（5）：587-595.
［29］ 徐秀峰， 张蓬洲.用XPS 表征氧、氮、硫元素的存在形态［J］.煤炭转化， 1996， 19（1）：72-77.XU Xiufeng， ZHANG Pengzhou.The XPS study of forms of oxy⁃gen， nitrogen and sulphur elements in gas coal［J］.Coal Conver⁃sion， 1996， 19（1）：72-77.
［30］ 谢克昌.煤的结构与反应性［M］.北京：科学出版社， 2002.
［31］ 相建华， 曾凡桂， 李 彬， 等.成庄无烟煤大分子结构模型及其分子模拟［J］.燃料化学学报， 2013， 41（4）：391-399.XIANG Jianhua， ZENG Fangui， LI Bin，et al.Construction of macromolecular structural model of anthracite from Chengzhuang coal mine and its molecular simulation［J］.Journal of Fuel Chem⁃istry and Technology， 2013， 41（4）：391-399.
［32］ 相建华， 曾凡桂， 梁虎珍， 等.不同变质程度煤的碳结构特征及其演化机制［J］.煤炭学报， 2016， 41（6）：1498-1506.XIANG Jianhua， ZENGFangui， LIANG Huzhen， et al.Carbon structure characteristics and evolution mechanism of different rank coals［J］.Journal of China Coal Society， 2016， 41（6）：1498-1506.
［33］ 相建华， 曾凡桂， 梁虎珍， 等.兖州煤大分子结构模型构建及其分子模拟［J］.燃料化学学报， 2011， 39（7）：481-488.XIANG Jianhua， ZENG Fangui， LIANG Huzhen， et al. Model construction of the macromolecular structure of Yanzhou coal and its molecular simulation［J］.Journal of Fuel Chemistry and Tech⁃nology， 2011， 39（7）：481-488.
［34］ 葛 涛， 张明旭， 马祥梅.新阳炼焦煤结构的FTIR 和XPS 谱学表征［J］.光谱学与光谱分析， 2017，37（8）：2406-2411.GE Tao， ZHANG Mingxu， MA Xiangmei.XPS and FTIR spectros⁃copy characterization about the structure of coking coal in Xinyang［J］.Spectroscopy and Spectral Analysis， 2017，37（8）：2406-2411.
［35］ 张 科， 姚素平， 胡文瑄， 等.煤红外光谱的精细解析及其煤化作用机制［J］.煤田地质与勘探， 2009， 37（6）：8-13.ZHANG Ke，YAO Suping， HU Wenxuan，et al.Analysis on infra⁃red spectra characteristic of coal and discussion of coalification mechanism［J］.Coal Geology ＆ Exploration， 2009， 37（6）：8-13.
［36］ 张双全.煤化学［M］.徐州：中国矿业大学出版社， 2009.
［37］ 朱学栋， 朱子彬， 韩崇家， 等.煤中含氧官能团的红外光谱定量分析［J］.燃料化学学报， 1999，27（4）：335-339.ZHU Xuedong， ZHUZibin， HAN Chongjia， et al.Quantitative a⁃nalysis of oxygen functional groups in coal by infrared spectroscopy［J］.Journal of Fuel Chemistry and Technology，1999（4）：335-339.
［38］ 张 莉， 曾凡桂， 相建华.内蒙五牧场矿区11 号煤层原煤大分子结构特征及其形成机制［J］.燃料化学学报， 2013， 41（11）：1294-1302.ZHANG Li， ZENG Fangui， XIANG Jianhua. Macromolecular structure and formation mechanism of raw coal in coal seam 11 of Wumuchang district， Inner Mongolia ［ J］. Journal of Fuel Chemistry and Technology， 2013， 41（11）：1294-1302.
［39］ 周 强.中国煤中硫氮的赋存状态研究［J］.洁净煤技术，2008， 53（1）：73-77.ZHOU Qiang.Study on occurrence mode of sulfur and nitrogen in coal in China［J］.Clean Coal Technology，2008，53（1）：73-77.
［40］ TAKANOHASHI T，KAWASHIMA H.Construction of amodel structure for Upper Freeport coal using13C-NMR chemical shift calculations［J］.Energy ＆ Fuels， 2002， 16（2）：379-387.
［41］ 姜永泼.屯兰2 号镜煤大分子聚集态结构模型的构建及分子模拟［D］.太原：太原理工大学， 2018.
［42］ 马延平， 相建华， 李美芬， 等.柳林3 号镜煤吡啶残煤大分子结构模型及分子模拟［J］.燃料化学学报， 2012， 40（11）：1300-1309.MA Yanping， XIANG Jianhua， LI Meifen， et al.Macromolecular structural model of the pyridine extracted residue of vitrain from No.3 coalbed，Liulin and molecular simulation［J］.Journal of Fuel Chemistry and Technology， 2012， 40（11）：1300-1309.
［43］ Zhang W， He M， Wei H， et al. Molecular dynamics simulations of interaction between sub-bituminous coal and water［J］.Molec⁃ular Simulation， 2018，44（9）：769-773.
［44］ 李 军，冯 杰， 李文英， 等.强弱还原煤聚集态对其可溶性影响的分子力学和分子动力学分析［J］.物理化学学报，2008， 24（12）：2297-2303.LI Jun， FENG Jie， LI Wenying，et al. Determining influence of the aggregative state of deoxidized coal on its extraction by molec⁃ular mechanics and molecular dynamics analysis［J］.Acta Physico-Chimica Sinica， 2008， 24（12）：2297-2303.
［45］ MARZEC A.Intermolecular interactions of aromatic hydrocarbons in carbonaceous materials：A molecular and quantum mechanics［J］.Carbon， 2000， 38（13）：1863-1871.
［46］ 马汝嘉， 张 帅， 侯丹丹，等.陕西凤县高煤级煤分子结构模型的 构 建 与 结 构 优 化［J］. 煤 炭 学 报， 2019， 44（6）：1827-1835.MA Rujia， ZHANG Shuai， HOU Dandanet al.Model construction and optimization of molecule structure of high-rank coal in Feng County，Shaanxi Province［J］.Journal of China Coal Society，2019， 44（6）：1827-1835.
［47］ CARLSON G A.Computer simulation of the molecular structure of bituminous coal［J］.Energy ＆ Fuels， 1992， 6（6）：771-778.
［48］ DONG T， MURATA S， WIURA M，et al. Computer-aided mo⁃lecular design study of coal model molecules.3.density simulation for model structures of bituminous Akabira coal［J］.Energy ＆ Fu⁃els， 1993， 7（6）：1123-1127.
［49］ 贾建波.神东煤镜质组结构模型的构建及其热解甲烷生成机理的分子模拟［D］.太原：太原理工大学， 2010.
［50］ 吴文忠.神东煤惰质组结构特征及其与CH4、CO2和H2O 相互作用的分子模拟［D］.太原：太原理工大学， 2010.
［51］ 马延平.柳林3 号煤的超分子构建及分子模拟［D］.太原：太原理工大学， 2012.
