煤炭加工与利用
煤作为一种组成极其复杂的大分子物质体系,其结构决定其反应性及综合利用的途径,因此,煤化学结构模型的研究有助于在分子水平上研究煤的反应路径和机理[1]。碳是煤结构单元的骨架,是煤中有机质的主要组成元素,煤中碳质量分数一般在70%以上。氧是煤中最丰富的杂原子,对煤的亲水性、表面电性、表面作用具有很大影响[2]。碳、氧结构是煤化学和煤化工的重要研究内容。李霞等分别利用红外光谱和拉曼光谱对低中煤级煤结构演化进行了表征,获知煤结构演化特征与第1,2次煤化作用发生的跃变密切相关[3-4]。刘钦甫等对中高煤阶煤结构的红外光谱学特征进行了研究,认为在此煤化阶段煤中芳香烃类物质增多[5]。相建华、石开仪、马汝嘉等先后构建了山西、内蒙古、辽宁、陕西等地的无烟煤及长焰煤大分子结构模型,并探讨了其形成机制[6-8]。
炼焦煤是冶金、铸造及化工等部门的重要原料,煤的焦化是煤炭综合利用的重要途径。中国焦炭产量居世界第1位,因此,炼焦煤是最受关注的煤炭资源之一。炼焦煤占中国查明煤炭资源储量的27%左右,其中约有1/2 的肥煤、瘦煤为高硫煤,1/3 的焦煤为高硫、高灰煤[9]。但是,对高硫炼焦煤中碳氧结构的认知还不够深入,大分子结构模型鲜有报道。
13C 交叉极化/魔角旋转-核磁共振(13C CP/MAS NMR)、X射线光电子能谱(XPS)和傅里叶红外光谱(FTIR)等光谱学仪器对煤表面化学特性均具有高度识别能力,能够实现在对煤样不造成破坏的条件下,直接获取煤中碳结构参数,被广泛应用于煤结构的研究中[10-14]。笔者针对山西高阳炼焦煤,通过13C CP/MAS NMR、XPS和FTIR的联合表征和比较研究,获取煤中碳结构和氧结构精准参数,为构建炼焦煤大分子结构模型提供依据。
采集山西高阳炼焦煤,按照GB/T474—2008制样,取200目样品密封保存。利用HD-GF500全自动工业分析仪测定煤中水分、灰分、挥发分。通过Multi EA 4000元素分析仪和SDS 601定硫仪分别测定煤中C,H,N和S质量分数,根据元素分析数据计算碳氢、碳氧原子比。煤质分析结果见表1。
表1 煤质分析
Table 1 Coal quality analysis
工业分析/%MadAadVdaf元素分析/%CdafHdafOdafNdafSdaf原子比H/CO/C1.6012.1334.6784.785.305.521.293.110.750.05
FTIR测试设备是IR-Tracer100傅里叶变换红外光谱仪,KBr为载体,煤样和载体以1∶100的比例进行混合研磨,在压片机上压制成0.1~1.0 mm的薄片,测试范围为4 000~400 cm-1,分辨率为16 cm-1,扫描次数累加到32次。
利用ESCALAB 250Xi型X-光电子能谱仪进行XPS测试,Mg/Al Kα双阳极非单色化X射线源,离子枪能量为100~4 000 eV,X射线单色器选择分析区域为30~400 μm,步长5 μm。能量分辨率优于0.8 eV。
13C CP/MAS-NMR测试在JNM-ECZ600R核磁共振谱仪上完成,CP/MAS固体双共振探头,3.2 mm ZrO2转子,MAS转速为(12±0.03) kHz,扫描累加1 100次,脉宽为0.1 μs,脉冲延迟时间为3 s。
选择Peakfit V4软件对FTIR 和13C CP/MAS-NMR图谱进行拟合,XPS图谱拟合软件是XPS Peak4.1。
高阳炼焦煤的13C CP/MAS-NMR拟合谱如图1所示,拟合优度R2为0.998,拟合效果良好。主要由化学位移为0~90×10-6和100×10-6~160×10-6的两个峰群组成,分别表征煤中的脂肪碳和芳香碳结构[15]。根据拟合特征峰的化学位移和面积获取的碳结构及其质量分数见表2。
图1 高阳炼焦煤13C CP/MAS-NMR拟合图谱
Fig.113C CP/MAS-NMR fitting spectrum of Gaogang coking coal
表2 煤中碳结构及其质量分数
Table 2 Content of carbon structure in coal
化学位移/10-6官能团结构质量分数15R—CH30.0421Ar—CH30.0233,43—CH2—,—CH,C0.0950O—CH3,O—CH2—0.0162,71O—CH0.0181R—O—R0.01105,115,121,125Ar—H0.26129,133,139Bridgehead(C—C),Ar—R0.48150Ar—O0.04167,181—COOH0.01194—CO0.02
根据煤中不同碳结构及其质量分数,利用式(1)~(6)计算芳碳率等煤结构参数[16-17],结果见表3。
fa+fal=1
(1)
(2)
(3)
(4)
(5)
(6)
式中,分别表示芳碳率、脂碳率、甲基碳、亚甲基碳和次甲基碳、氧接脂碳、羧基碳和羰基碳、芳碳、质子化芳碳、非质子化芳碳、烷基取代芳碳、芳香桥碳、氧接芳碳;Xb为芳核平均结构尺寸,为芳香化合物的桥碳与周碳的比值(桥碳比)。
为为0.48,分别占脂肪碳结构和芳香碳结构的60%以上。表明高阳炼焦煤中脂肪碳结构以亚甲基和次甲基为主,具有较多的脂肪侧链;而芳香碳结构主要以非质子化芳碳形式存在。芳核平均结构尺寸为0.43,因此,煤中芳香结构以三环和四环缩合为主。
表3 高阳煤13C CP/MAS-NMR 结构参数
Table 3 13C CP/MAS-NMR structural parameters of Gaoyang coal
fafalf∗alfHalfOalfCaf′afHafNafSafBafPaXb0.770.230.060.140.030.030.740.260.480.120.320.040.43
2.2.1 芳香烃结构
700~900 cm-1是煤中芳香烃结构的FTIR吸收振动波数范围,苯环的取代方式决定了芳香烃结构的拟合峰归属最多为6个[18]。高阳煤芳香烃结构拟合谱如图2(a)所示,R2为0.998,根据特征峰波数和面积分别解析峰结构归属及相对质量分数,结果见表4。6个特征峰分别归属于苯环二取代、三取代、四取代和五取代4种结构,其中,苯环五取代结构相对质量分数达到41.42%,是煤中最主要的芳香结构,苯环四取代和三取代结构质量分数分别为30.65%,19.82%,苯环二取代结构不足10%。
图2 FTIR拟合谱
Fig.2 FTIR fitting spectrum
表4 煤中芳香结构FTIR解析
Table 4 FTIR analysis for aromatic hydrocarbon structure in coal
峰波数/cm-1峰面积相对质量分数/%化学结构8742.9135.19苯环五取代8610.516.23苯环五取代8170.536.42苯环四取代8002.0024.23苯环四取代7521.6419.82苯环三取代7400.678.11苯环二取代
2.2.2 脂肪烃结构
2 800~3 000 cm-1是煤中脂肪烃结构的波数区间。脂肪烃结构拟合谱图如图2(b)所示,R2达到0.999,拟合度很高,拟合峰归属及其相对质量分数见表5。根据煤中脂氢的类型,脂肪结构的FTIR比较合理的拟合峰个数为6~9[19]。高阳煤脂肪烃的9个拟合峰分别归属于甲基、亚甲基和次甲基,其中,亚甲基是最主要的脂肪烃结构,占比为41.85%,甲基、次甲基占比分别为29.86%,28.29%。可见,含有较多的烷基侧链和环状脂肪烃。
表5 煤中脂肪烃结构FTIR解析
Table 5 FTIR analysis for aliphatic hydrocarbon structure in coal
峰波数/cm-1峰面积相对质量分数/%化学结构2 9990.140.87RCH32 9620.724.50RCH32 9481.7510.81RCH32 9292.1913.68RCH32 9172.0612.87R2CH22 9042.2814.24R3CH2 8852.2514.05R3CH2 8611.8311.43R2CH22 8442.8117.55R2CH2
2.2.3 煤FTIR结构参数
在利用FTIR解析计算煤结构参数时,近似认为煤中只有芳香氢和脂肪氢两类氢原子存在。根据煤中芳香烃和脂肪烃结构的FTIR 解析结果及煤质分析数据计算高阳煤结构的芳氢率和芳碳率。
(7)
(8)
其中,n(Har)为煤芳香结构的氢原子数;n(H),n(C)分别为煤结构中的总氢原子数和总碳原子数;I为波数区间内吸收峰面积;far-F为FTIR中煤的芳碳率;n(Cal),n(Hal)分别表示煤中的脂肪碳原子数、脂肪结构中的氢原子数;n(Hal)/n(Cal)一般取经验值1.8。计算获取高阳炼焦煤中芳氢率和芳碳率分别为0.34和0.73。
2.2.4 羟基基团
羟基在 FTIR中的吸收振动波数为3 700~3 200 cm-1。FTIR拟合谱图如图2(c)所示,R2为0.996,拟合度略低于芳香烃结构和脂肪烃结构,主要是因为羟基的实验谱中出现了较多的振动峰。峰结构及相对质量分数表征结果见表6。煤中羟基主要有6种,在断裂、交联键时具有很强的活化效应,主要存在于端基和侧链中,因此,羟基的FTIR拟合一般有4~7个峰[20]。图2(c)共有6个拟合峰,分别归属于游离羟基、羟基 π 氢键、羟基自缔合氢键、醚氧键与羟基形成的氢键四种结构。与芳环上π电子形成的羟基 π 氢键是高阳炼焦煤最主要的羟基结构,占比为73.20%。煤中游离羟基质量分数很少,仅占2.42%,表明煤中脂链的环化与官能团有强烈的缩合作用,减弱了游离羟基的伸缩振动。高阳煤中羟基自缔合氢键、醚氧键与羟基形成的氢键质量分数较高,占羟基总量的24.38%,氢键是煤缔合模型的重要标志,因此,煤中缔合结构以多聚体形式为主。
表6 煤中羟基基团FTIR解析
Table 6 FTIR analysis for hydroxyl groups in coal
峰波数/cm-1峰面积相对质量分数/%化学结构3 6860.350.33游离羟基3 6490.810.77游离羟基3 6231.381.32游离羟基3 49747.0045.00OH…π3 41529.4528.20OH…π3 35014.0213.42OH…OH3 28111.4410.96OH…O
2.2.5 含氧官能团
1 800~1 000 cm-1是煤中除羟基以外(如羧基、羰基、醚氧键等)含氧官能团的红外波数范围。由于该范围还包括CS,Si—O—Si,Si—O—C等杂原子以及—CH2,—CH3、芳核CC等脂肪烃和芳香烃的吸收峰[5]。因此含氧官能团是FTIR最复杂的拟合区间,一般有13~18个拟合峰[21]。而高阳煤在1 030~1 020 cm-1 区域出现了较强的硅伸缩振动吸收峰,为了减少干扰将拟合区间调整为 1 800~1 060 cm-1。拟合谱图如图2(d)所示,共有15个特征峰,解析结果见表7。R2为0.996,拟合度与煤中羟基相当,拟合效果理想。含氧官能团主要以共轭羰基和酚羟基的形式存在,芳基醚和羧基含量较少。
表7 煤中含氧官能团 FTIR解析
Table 7 Oxygen containing functional groups in coal FTIR analytical
峰波数/cm-1峰面积相对质量分数/%化学结构1 7121.871.37羧基1 6735.363.92羧基1 64715.2911.18共轭羰基1 61325.218.45共轭羰基1 58112.829.38芳核CC1 5488.346.10芳核CC1 5083.082.25芳核CC1 4677.325.35—CH3,—CH21 43414.6210.70—CH3,—CH21 4072.641.93—CH3,—CH21 38215.8411.59—CH31 3362.882.11酚羟基1 3090.700.51酚羟基1 28714.0110.25酚羟基1 0776.714.91芳基醚
2.3.1 碳结构
煤中碳谱的XPS拟合,一般会出现芳碳、脂碳、酚碳/醚基和羧基4类拟合峰和表征结构[22]。高阳煤XPS碳谱的拟合谱如图3(a)所示,共有3个拟合峰,电子结合能、化学结构及各种结构相对质量分数见表8。
表8 煤中碳结构XPS解析
Table 8 XPS analysis for carbon structure in coal
电子结合能/eV 化学结构峰面积相对质量分数/%284.73C—C65 733.4876.70285.30C—H15 016.0317.52286.00C—O4 951.485.78
图3 碳结构和氧结构的XPS拟合图谱
Fig.3 XPS fitting spectrum of carbon and oxygen structure
XPS碳谱分析结果显示,芳构碳占比为76.70%,与13C CP/MAS-NMR芳碳率为0.77,FTIR芳碳率为0.73的分析结果非常接近。醚碳/酚碳质量分数为5.78%,没有出现羧基的拟合特征峰,醚基和酚基等含氧官能团质量分数为5.78%,表明煤样被氧化程度不高。
2.3.2 氧结构
XPS不能区分含氧结构中的醚基和羟基,碳氧单键可以认为是两种氧结构的总量。GRZYBEK T等[23]认为电子结合能为531.3和 532.8 eV分别是煤中碳氧双键和碳氧单键的吸收峰,周刚等[24]则认为(531.3±0.2),(532.8±0.3)与(534.1±0.4) eV 等3处吸收峰分别归属于羧基、碳氧单键和碳氧双键。
高阳煤XPS氧谱的拟合谱如图3(b)所示,532.90,531.80,534.10 eV分别归属于碳氧单键、羰基和羧基的特征峰比较合理,表征结果见表9。高阳煤中主要以碳氧单键结构存在,530.00 eV 附近没有出现吸收峰,说明无机氧质量分数极少,煤样的氧化程度较低,与碳谱的解析结果一致。羰基和羧基XPS和FTIR分析结果中的比例分别为4.5∶1和5.6∶1,因此,煤中羰基是羧基个数的5倍左右。
2.4.1 芳香碳结构
高阳炼焦煤中碳质量分数为84.78%,芳香结构单元以3~5个缩合环数为主,芘、菲和蒽是中等煤化度烟煤基本单元核的主要结构环[25]。根据高阳煤的Xb设计分子模型中芳香结构单元及个数见表10,分子模型中芳香烃碳原子个数为118。f′a和分别为0.74和0.03,因此,羰基碳和羧基碳之和为5个。
表9 煤中氧结构XPS解析
Table 9 XPS analysis for oxygen structure in coal
电子结合能/eV 化学结构峰面积相对质量分数/%531.80CO13 650.9535.25532.90C—O22 028.8056.90534.10OC—O3 040.677.85
表10 高阳炼焦煤分子模型中芳香结构单元
Table 10 Aromatic structural unit in coal molecular structure model
2.4.2 脂肪碳结构
煤中脂肪结构以烷基侧链、环烷烃和氢化芳烃的形式存在,侧链长度随煤化程度的增加而迅速减小。根据芳碳率的计算结果,分子模型中脂肪烃碳原子个数为22~27。其中,甲基碳、亚甲基碳、次甲基碳的个数分别为7~8,9~11,6~8。高阳煤氢碳原子个数比为0.75,因此,分子模型中的氢原子数为109~113。
2.4.3 氧原子结构
根据煤质分析结果及分子模型碳原子数计算氧原子个数为7,高阳炼焦煤中硫质量分数达到3.11%,属于高硫煤。煤中硫主要有无机硫、硫醇、硫醚、噻吩、砜、亚砜等结构,烟煤中有机硫质量分数较高,考虑煤中砜或者亚砜可能会占用氧原子,因此分子模型中羰基和酚羟基个数为6、醚氧键和有机硫共同占用1个氧原子比较合理。
(1)高阳炼焦煤13C CP/MAS-NMR,FTIR,XPS芳碳率的解析结果非常接近,分别为0.77,0.73,0.77。煤中和Xb分别为0.06,0.14,0.03,0.03,0.26,0.12,0.32,0.04和0.43。
(2)高阳炼焦煤中的芳香烃结构以苯环五取代、苯环四取代和苯环三取代为主,其中苯环五取代结构相对质量分数达到41.42%。亚甲基是最主要的脂肪烃结构,占比为41.85%,煤中含有较多的烷基侧链和环状脂肪烃。含氧结构主要以羰基和酚羟基的形式存在,芳基醚和羧基质量分数较少。
(3)芳香烃碳原子个数为118的煤分子模型中,脂肪烃碳原子个数为27~32,其中,甲基碳、亚甲基碳、次甲基碳、羰基和羧基的个数分别为7~8,9~11,6~8,5。氧原子个数为7,其中,羰基和酚羟基为6、醚氧键和有机硫中的砜或者亚砜结构共用1个氧原子。
(4)高阳炼焦煤精准分子模型的构建需要有更多的碳氧原子组成,从而准确地获取羰基、羟基、羧基、醚氧键等化学结构信息,同时需要对硫、氮等杂原子结构做进一步的研究。
[1] 王凤,李光跃,李莹莹,等.煤化学结构模型研究进展及应用[J].洁净煤技术,2016,22(1):26-32.
WANG Feng,LI Guangyue,LI Yingying,et al.Application and research progress of coal chemical structure model[J].Clean Coal Technology,2016,22(1):26-32.
[2] QIN Zhihong.New advances in coal structure model[J].International Journal of Mining Science and Technology,2018,(28):541-559.
[3] 李霞,曾凡桂,王威,等.低中煤级煤结构演化的FTIR表征[J].煤炭学报,2015,40(12):2900-2908.
LI Xia,ZENG Fangui,WANG Wei,et al.FTIR characterization of structural evolution in low-middle rank coals[J].Journal of China Coal Society,2015,40(12):2900-2908.
[4] 李霞,曾凡桂,王威,等.低中煤级煤结构演化的拉曼光谱表征[J].煤炭学报,2016,41(9):2298-2304.
LI Xia,ZENG Fangui,WANG Wei,et al.Raman characterization of structural evolution in the low-middle rank coals[J].Journal of China Coal Society,2016,41(9):2298-2304.
[5] 刘琬玥,刘钦甫,刘霖松,等.沁水盆地北部中高煤阶煤结构的FTIR特征研究[J].煤炭科学技术,2019,47(2):181-187.
LIU Wanyue,LIU Qinfu,LIU Linsong,et al.Study on FTIR features of middle and high rank coal structure in north part of Qinshui Basin[J].Coal Science and Technology,2019,47(2):181-187.
[6] 相建华,曾凡桂,李彬,等.成庄无烟煤大分子结构模型及其分子模拟[J].燃料化学学报,2013,41(4):391-399.
XIANG Jianhua,ZENG Fangui,LI Bin,et al.Construction of macromolecular structural modelof anthracite from Chengzhuang coal mine and its molecular simulation[J].Journal of Fuel Chemistry and Technology,2013,41(4):391-399.
[7] 石开仪,陶秀祥,李志,等.利用红外光谱构建抚顺煤大分子结构模型[J].高分子通报,2013,(3):61-66.
SHI Kaiyi,TAO Xiuxiang,LI Zhi,et al.Study of construction of Fushun coal macromolecule structural model by infrared spectroscopy[J].Polymer Notification,2013,(3):61-66.
[8] 马汝嘉,张帅,侯丹丹,等.陕西凤县高煤级煤分子结构模型的构建与结构优化[J].煤炭学报,2019,44(6):1827-1835.
MA Rujia,ZHANG Shuai,HOU Dandan,et 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.
[9] 唐跃刚,贺鑫,程爱国,等.中国煤中硫含量分布特征及其沉积控制[J].煤炭学报,2015,40(9):1977-1988.
TANG Yuegang,HE Xin,CHENG Aiguo,et al.Occurrence and sedimentary control of sulfur in coals of China[J].Journal of China Coal Society,2015,40(9):1977-1988.
[10] DING Dianshi,LIU Guijian,FU Biao.Influence of carbon type on carbon isotopic composition of coal from the perspective of solid-state13C NMR[J].Fuel,2019,245(6):174-180.
[11] JING Zhenhua,RODRIGUES Sandra,STROUNINA Ekaterina,et al.Use of FTIR,XPS,NMR to characterize oxidative effects of NaClO on coal molecular structures[J].International Journal of Coal Geology,2019,201(1):1-13.
[12] LEVI Gianluca,SENNECA Osvalda,CAUSA Mauro,et al.Prob-ing the chemical nature of surface oxides during coal char oxidation by high-resolution XPS[J].Carbon,2015,90(8):181-196.
[13] 赵云刚,李美芬,曾凡桂,等.伊敏褐煤不同化学组分结构特征的红外光谱研究[J].煤炭学报,2018,43(2):546-553.
ZHAO Yungang,LI Meifen,ZENG Fangui,et al.FTIR study of structural characteristics of different chemical components from Yimin Lignite[J].Journal of China Coal Society,2018,43(2):546-553.
[14] JIANG Jingyu,YANG Weihua,CHENG Yuanping,et al.Molecular structure characterization of middle-high rank coal via XRD,Raman and FTIR spectroscopy:Implications for coalification[J].Fuel,2019,239(3):559-572.
[15] 陈丽诗,王岚岚,潘铁英,等.固体核磁碳结构参数的修正及其在煤结构分析中的应用[J].燃料化学学报,2017,45(10):1153-1163.
CHEN Lishi,WANG Lanlan,PAN Tieying,et al.Calibration of solid state NMR carbon structural parameters and application in coal structure analysis[J].Journal of Fuel Chemistry and Technology,2017,45(10):1153-1163.
[16] 刘鹏,王岚岚,张德祥.先锋褐煤在水热处理过程中的结构演绎[J].燃料化学学报,2016,44(2):129-137.
LIU Peng,WANG Lanlan,ZHANG Dexiang.Structural evolution of Xianfeng lignite during hydrothermal treatment[J].Journal of Fuel Chemistry and Technology,2016,44(2):129-137.
[17] 潘铁英,张玉兰,苏克曼.波谱解析法[M].上海:华东理工大学出版社,2009.
[18] 葛涛,马祥梅,张明旭.新阳炼焦煤结构的FTIR和XPS谱学表征[J].光谱学与光谱分析,2017,37(8):2146-2152.
GE Tao,MA Xiangmei,ZHANG Mingxu.XPS and FTIR spectroscopy characterization about the structure of coking coal in Xinyang[J].Spectroscopy and Spectral Analysis,2017,37(8):2146-2152.
[19] 安文博,王来贵,刘向峰,等.基于FTIR和XRD法分析阜新长焰煤结构特征[J].高分子通报,2018,(3):67-74.
AN Wenbo,WANG Laigui,LIU Xiangfeng,et al.Analysis the structural characteristics of fuxin long flame coal based on FTIR and XRD Experiments[J].Polymer Notification,2018,(3):67-74.
[20] ZHANG Zheng,QIN Yong,ZHUANG Xinguo,et al.Poroperm characteristics of high-rank coals from Southern Qinshui Basin by mercury intrusion,SEM-EDS,nuclear magnetic resonance and relative permeability analysis[J].Journal of Natural Gas Science and Engineering,2018,51:116-128.
[21] KRUSZEWSKI Lukasz,NARUSIS Fabianska,CIESIELCZUK Justyna,et al.First multi-tool exploration of a gas-condensate-pyrolysate system from the environment of burning coal mine heaps:An in situ FTIR and laboratory GC and XRD study based on Upper Silesian materials[J].Science of The Total Environment,2018,640-641(11):1044-1071.
[22] LIU Fenrong,LI Wen,GUO Huiqing,et al.XPS study on the change of carbon-containing groups and sulfur transformation on coal surface[J].Journal of Fuel Chemistry and Technology,2011,39(2):81-84.
[23] GRZYBEK Teresa,PIETRZAK Robert,WACHOWSKA Helena.X-ray photoelectron spectroscopy study of oxidized coals with different sulphur content[J].Fuel Process Technol,2002,(77):1-7.
[24] 周刚,徐翠翠,邱晗.基于核磁-能谱实验的中变质烟煤煤尘低润湿性分析[J].化工进展,2016,35(10):3441-3446.
ZHOU Gang,XU Cuicui,QIU Han.Analysis of the low wettability about the bituminous coal dust with medium metamorphic grade based on NMR and XPS experiment[J].Chemcial Industry and Engineering Progress,2016,35(10):3441-3446.
[25] 魏帅,严国超,张志强,等.晋城无烟煤的分子结构特征分析[J].煤炭学报2018,43(2):555-562.
WEI 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.
GE Tao,WANG Meng,LI Fen,et al.Study and spectroscopy characterization of carbon & oxygen structures of Gaoyang coking coal[J].Journal of China Coal Society,2021,46(3):1024-1031.