CN104390928B - 一种掺杂掺假煤炭的近红外光谱识别方法 - Google Patents
一种掺杂掺假煤炭的近红外光谱识别方法 Download PDFInfo
- Publication number
- CN104390928B CN104390928B CN201410579312.2A CN201410579312A CN104390928B CN 104390928 B CN104390928 B CN 104390928B CN 201410579312 A CN201410579312 A CN 201410579312A CN 104390928 B CN104390928 B CN 104390928B
- Authority
- CN
- China
- Prior art keywords
- coal
- sample
- spectrum
- detection model
- original
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 239000003245 coal Substances 0.000 title claims abstract description 146
- 238000000034 method Methods 0.000 title claims abstract description 71
- 238000002329 infrared spectrum Methods 0.000 title abstract description 14
- 238000001228 spectrum Methods 0.000 claims abstract description 114
- 238000001514 detection method Methods 0.000 claims abstract description 81
- 238000012360 testing method Methods 0.000 claims abstract description 11
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 28
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 15
- 229910052799 carbon Inorganic materials 0.000 claims description 15
- 238000007796 conventional method Methods 0.000 claims description 15
- NINIDFKCEFEMDL-UHFFFAOYSA-N Sulfur Chemical compound [S] NINIDFKCEFEMDL-UHFFFAOYSA-N 0.000 claims description 12
- 239000005864 Sulphur Substances 0.000 claims description 10
- 238000004611 spectroscopical analysis Methods 0.000 claims description 10
- 238000003860 storage Methods 0.000 claims description 9
- 238000011156 evaluation Methods 0.000 claims description 3
- 238000007621 cluster analysis Methods 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 claims 1
- 239000003610 charcoal Substances 0.000 claims 1
- 238000004458 analytical method Methods 0.000 abstract description 31
- 238000001320 near-infrared absorption spectroscopy Methods 0.000 abstract description 11
- 239000003153 chemical reaction reagent Substances 0.000 abstract description 8
- 230000000694 effects Effects 0.000 abstract description 2
- 239000002802 bituminous coal Substances 0.000 description 21
- RHZUVFJBSILHOK-UHFFFAOYSA-N anthracen-1-ylmethanolate Chemical compound C1=CC=C2C=C3C(C[O-])=CC=CC3=CC2=C1 RHZUVFJBSILHOK-UHFFFAOYSA-N 0.000 description 18
- 239000003830 anthracite Substances 0.000 description 18
- 238000012545 processing Methods 0.000 description 17
- 239000002956 ash Substances 0.000 description 10
- 238000012937 correction Methods 0.000 description 10
- 239000000126 substance Substances 0.000 description 10
- 230000008569 process Effects 0.000 description 8
- 230000003595 spectral effect Effects 0.000 description 8
- 238000003556 assay Methods 0.000 description 7
- 238000009826 distribution Methods 0.000 description 5
- 239000000203 mixture Substances 0.000 description 5
- 238000000513 principal component analysis Methods 0.000 description 4
- 238000011282 treatment Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000009841 combustion method Methods 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 239000003077 lignite Substances 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000010238 partial least squares regression Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000010521 absorption reaction Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 239000000571 coke Substances 0.000 description 2
- 238000002485 combustion reaction Methods 0.000 description 2
- 230000000052 comparative effect Effects 0.000 description 2
- 239000000470 constituent Substances 0.000 description 2
- 238000007405 data analysis Methods 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011049 filling Methods 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 239000005416 organic matter Substances 0.000 description 2
- 238000010561 standard procedure Methods 0.000 description 2
- 229910052717 sulfur Inorganic materials 0.000 description 2
- 239000011593 sulfur Substances 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- -1 volatile matter Substances 0.000 description 2
- 235000002918 Fraxinus excelsior Nutrition 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000001311 chemical methods and process Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 230000002301 combined effect Effects 0.000 description 1
- 238000013329 compounding Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000012850 discrimination method Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000009931 harmful effect Effects 0.000 description 1
- 238000001093 holography Methods 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 125000004435 hydrogen atom Chemical group [H]* 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000011017 operating method Methods 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 229910052698 phosphorus Inorganic materials 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000012847 principal component analysis method Methods 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 239000010453 quartz Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 1
- 238000011895 specific detection Methods 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 238000003756 stirring Methods 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种掺杂掺假煤炭的近红外光谱识别方法。是利用NIRS分析仪扫描收集不同煤种的标准煤炭样品的近红外原始光谱，分析待测煤炭样品的近红外光谱，将待测煤炭样品的近红外光谱与标准煤炭样品的近红外原始光谱进行比较，识别待测煤炭样品是否掺杂掺假。本发明简单易行地实现一次性同时分析得到7项关键指标的检测并结合光谱图和指标检测结果判断识别煤炭样品是否掺杂掺假。本发明提供的新方法检测得到煤炭样品中7项关键成分以及做出判断的谱图依据只需几分钟，而且无需称样、无需使用化学试剂或高温、高压、大电流等测试条件，不会产生化学、生物或电磁污染，不会对操作人员和环境造成不良影响。
Description
技术领域
本发明涉及煤炭分析检测技术领域，更具体地，涉及一种掺杂掺假煤炭的近红外光谱识别方法。
背景技术
中国是世界第一产煤大国，煤炭资源是我国的第一大能源，是工业发展所需的基础能源材料。我国煤炭资源非常丰富，占我国总能源的70％左右。
但由于我国煤炭资源本身灰分较高、发热量较低和洗煤能力跟不上煤炭产量的增长速度，以及一些中间商掺配劣质煤等原因，供给用户的商品煤质量始终满足不了生产的需求，目前缺乏有效解决煤炭混配、掺杂等难题的检测技术手段。
发明内容
本发明要解决的技术问题是针对现有煤炭掺假现象识别技术手段的不足，提供一种掺杂掺假煤炭的近红外光谱识别方法。
本发明的目的通过以下技术方案予以实现：
本发明研究总结不同煤种的煤炭样品在近红外DDS原始光谱上表现出不同的光谱特性，利用所述光谱特性实现在测定煤炭成分含量的同时鉴定样品是无烟煤或者是烟煤或者是贫瘦煤等煤炭品种，鉴定煤炭的品质，有效解决煤炭混配、掺杂等的识别难题。
具体地，所述一种掺杂掺假煤炭的近红外光谱识别方法包括以下步骤：
S1.用NIRS分析仪扫描收集不同煤种的标准煤炭样品的近红外原始光谱；
S2.将待测煤炭样品依次直接装满NIRS分析仪的进样器(不需称样)，启动扫描键，NIRS分析仪自动记录存储样品光谱，将所得待测煤炭样品的光谱与S1所得原始光谱进行比较，确定待测煤炭样品归属谱图类型，识别待测煤炭样品是否掺杂掺假。
优选地，S1所述不同煤种为无烟煤、烟煤、贫瘦煤、水汽煤、动力煤和/或褐煤。
进一步优选地，本发明还可以通过分析待测待测煤炭样品光谱的同时具体检测该样品中的7项指标的方法，更加准确地识别待测煤炭样品是否掺杂掺假。包括以下步骤：
S01.收集和制备若干个煤炭样品，常规方法分别测定每个样品的内水、灰分、挥发分、全硫、固定碳和可燃体的含量及发热量7项指标；
S02.用NIRS分析仪扫描收集所述煤炭样品的光谱数据和曲线；
S03.对S02所得样品的光谱数据进行处理，经回归计算获得7项指标的定标方程，修正和验证内水后建立检测模型；
S04.将待测煤炭样品依次直接装满NIR仪的进样器(不需称样)，启动扫描键，NIR仪自动记录存储样品光谱，将所得待测煤炭样品的光谱与S1所得原始光谱进行比较，确定待测煤炭样品归属谱图类型；选择相应的检测模型，获得检测结果。
优选地，所述7项指标为内水、挥发分、灰分、固定碳、可燃体、发热量和全硫。
本发明的有益效果如下：
煤炭现有检测标准的原理都是经典或现代的化学法和物理化学法：内水、挥发分、灰分、固定碳、可燃体都是经典重量分析方法，涉及天平、烘箱、高温炉等设备，碳化、灰化、恒重、称量、计算是经常需要进行的步骤，相当繁琐费时；发热量、全硫，目前最先进的是采用高温高压燃烧法，但是现代化的物理化学仪器的操作、维护、标定也很繁重。
长期以来，近红外光谱技术用来分析纯粹的有机物。因为近红外光谱的波数在4000cm-1以上(即2500nm以下)，因此，只有振动频率在2000cm-1以上的振动，才可能在近红外区内产生一级倍频，而能够在2000cm-1以上产生基频振动的主要是含氢官能团，如C-H、N-H、S-H和O-H的伸缩振动。几乎有机物中所有含氢基团的信息，都能在近红外光谱中得以反映。
煤炭是一种可燃性岩石。本发明首次将近红外光谱技术应用于煤炭这种由大部分有机物质和部分矿物质及水分组成的混合物质体系，并成功建立煤炭样品的检测分析方法，从而在一定程度上证明近红外光谱技术可以应用于分析无机物质，克服了现有技术偏见。
本发明利用煤炭样品的近红外光谱特性实现在测定煤炭成分含量的同时鉴定样品是无烟煤或者是烟煤或者是贫瘦煤等煤炭品种，更好地鉴定煤炭的品质，有效解决煤炭混配、掺杂等难题。
进一步地，本发明还可以通过分析待测待测煤炭样品光谱的同时具体检测该样品中的7项指标的方法，结合现有不同煤种关键指标的参考值，更加准确地识别待测煤炭样品是否掺杂掺假。
本发明与现行方法相比，具有显著的优越性。本发明仅需要采用一台近红外光谱分析仪，就可以代替现有技术多种、多台分析仪器，只需要磨样设备，不需要分析天平，通常只需一个人操作，并在几分钟内，通过采集一次被测样品的光谱，就可以同时完成煤炭样品掺假掺杂的鉴定以及煤炭样品7项关键指标的定量测定。在光谱采集过程中除消耗电能外，不需消耗任何试剂和标准物质，这样可以节省大量仪器设备的购置、操作、维修等费用，节省大量的时间和人力，极大地降低分析成本，显著提高检测工作的效率。
附图说明
图1无烟煤的NIRDDS原始光谱。
图2烟煤的NIRDDS原始光谱。
图3贫瘦煤的NIRDDS原始光谱。
图4水汽煤的NIRDDS原始光谱。
图5褐煤的NIRDDS原始光谱。
图6 509个无烟煤样品的近红外漫反射DDS原始光谱总图。
图7 134个烟煤样品的近红外漫反射DDS原始光谱总图。
图8 155个水汽煤样品的近红外漫反射DDS原始光谱总图。
图9 60个贫瘦煤样品的近红外漫反射DDS原始光谱总图。
图10 134个动力煤样品的近红外漫反射DDS原始光谱总图。
图11 179个未明确分组的样品的近红外漫反射DDS原始光谱总图。
图12 202个无烟煤样品的近红外漫反射DDS原始光谱总图(2008年样品)。
图13 235个未明确分组的样品的近红外漫反射DDS原始光谱总图(2008年样品)。
图14标样原始谱图(上为103f号无烟煤，下为101L号烟煤)。
图15标准煤样原始谱图。
图16 2008年无烟煤和标样无烟煤的原始谱图。
图17 2005-2008年烟煤和标样烟煤的原始谱图。
图18 2008年无烟煤和标样无烟煤原始光谱经(NONE+0011)处理图。
图19 2008年无烟煤和标样无烟煤原始光谱经(NONE+1441)处理图。
图20 2008年无烟煤和标样无烟煤原始光谱经(D+1441)处理图。
图21 2008年无烟煤和标样无烟煤原始光谱经(D+1441)处理图。
图22 2005-2008年烟煤和标样烟煤的原始光谱经(NONE+0011)处理图。
图23 2005-2008年烟煤和标样烟煤的原始光谱经(D+1441)处理图。
图24 2005-2008年烟煤和标样烟煤的原始光谱经(D+0011)处理图。
图25 509个无烟煤样品原始光谱经(NONE+1441)处理图。
图26 509个无烟煤样品原始光谱经(NONE+0011)处理图。
图27 509个无烟煤样品原始光谱经(D+1441)处理图。
图28 Y型(374个)原始光谱图。
图29 W型(1367个)原始光谱图。
图30 P型(436个)原始谱图。
图31 X型(52个)原始光谱。
图32 Y型原始谱图的GH值分布图。
图33 W型原始谱图的GH值分布图。
图34 P型原始谱图的GH值分布图。
图35 X型原始谱图的GH值分布图。
图36经交互验证的Y型原始光谱内水检测模型(MPLS+D+1441)。
图37经交互验证的W型原始光谱检内水测模型(PLS+D+0011)。
图38经交互验证的P型原始光谱内水检测模型(MPLS+D+0011)。
图39经交互验证的X型原始光谱检内水测模型(MPLS+NONE+2441)。
具体实施方式
下面结合附图和具体实施例进一步说明本发明。本发明实施例采用的设备可参照《高温管式炉燃烧法分析煤和焦碳分析样品中硫含量的标准试样方法》(ASTMD-4239-2010e1)、《煤的工业分析方法》(GB/T212-2008)、《煤的发热量测定方法》(GB/T213-2008)中所列设备和试剂。其他除非特别说明，本发明实施例采用的试剂和设备为本领域常规使用的试剂和设备。
实施例1
S1.用NIRS分析仪扫描收集不同煤种的标准煤炭样品的近红外原始光谱；本实施例以无烟煤、烟煤、贫瘦煤、水汽煤、动力煤和/或褐煤为例，收集的原始光谱见附图1至附图5所示。
S2.将待测煤炭样品依次直接装满NIRS分析仪的进样器(不需称样)，启动扫描键，NIRS分析仪自动记录存储样品光谱，将所得待测煤炭样品的光谱与S1所得原始光谱进行比较，确定待测煤炭样品归属谱图类型，识别待测煤炭样品是否掺杂掺假。
实施例2
本实施例提供一种同时检测煤炭样品中7项指标并识别待测煤炭样品是否掺杂掺假的方法，包括以下步骤：
S01.收集和制备若干个煤炭样品，常规方法分别测定每个样品的内水、灰分、挥发分、全硫、固定碳和可燃体的含量及发热量7项指标；
S02.用NIRS分析仪扫描收集所述煤炭样品的光谱数据和曲线；
S03.对S02所得样品的光谱数据进行处理，经回归计算获得7项指标的定标方程，修正和验证内水后建立检测模型；
将待测煤炭样品依次直接装满NIR仪的进样器(不需称样)，启动扫描键，NIR仪自动记录存储样品光谱，将所得待测煤炭样品的光谱与S1所得原始光谱进行比较，确定待测煤炭样品归属谱图类型；选择相应的检测模型，获得检测结果。
其中，S01所述收集和制备煤炭样品的方法是参照《商品煤样采取方法》(GB475-2008)和《煤样的制备方法》(GB474-2008)中的要求进行。所述常规方法测定样品的内水含量的方法是按《煤的工业分析方法》(GB/T212-2008)进行，所述常规方法测定样品发热量的方法是按《煤的发热量测定方法》(GB/T213-2008)中的要求进行，所述常规方法测定样品中的固定碳的方法是按《煤的工业分析方法》(GB/T212-2008)中的要求进行，所述常规方法测定样品的灰分含量的方法是按《煤的工业分析方法》(GB/T212-2008)进行，所述常规方法测定样品的挥发分含量的方法是按《煤的工业分析方法》(GB/T212-2008)进行，所述常规方法测定样品中的可燃体的方法是可燃体参照《煤的工业分析方法》(GB/T212-2008)中的要求进行，所述常规方法测定样品的全硫含量的方法是按《高温管式炉燃烧法分析煤和焦碳分析样品中硫含量的标准试样方法》(ASTMD-4239-2010e1)进行。
S02所述用NIRS分析仪扫描收集所述样品的光谱数据是将样品(不需称样)，依次直接装满NIR仪的进样器，采用数字光栅系统进行扫描，通过NIRS仪自动记录和存储样品光谱。
本实施例采用丹麦FOSS ANALYTICAL A/S公司生产的SY-3650－2型近红外分析仪，单色器：单光束全息数字光栅，波长范围1100～2500nm；光源：钨丝灯；检测器：PBS，自动移动输送多点定位检测，检测点大于32个；操作温度：15～32℃；分析时间：每分钟用连续光谱扫描32个子样品以上并完成其光谱分析；系统噪音：噪音信号值小于2×10-5AU；结果报告：在终端屏幕上显示结果和显示“已超越校正范围”的样品，同时由打印机报告结果，也可通过RS-232C接口与外部电子计算机相联。
操作方法：仔细搅拌煤样，装取适量的试样于石英制的样品槽中，装填过程中务必使样品厚度均匀，盖上小纸板，轻轻按压，使样品紧密分布。装填完毕后，观察样品槽表面的样品情况，如果发现有缝隙或样品有松动的现象，应重新装填试样。将装填好的样品槽置于进样轨道上，点击“扫描”，通过上下往复运动的进样方式，完成扫描检测和据的输出。整个检测操作的过程中，无需称样，无需添加任何化学试剂。
首先，采用数字光栅系统(Digital Dispersive System，DDS)进行扫描采集样品光谱。本发明科学地进行分类，将煤炭样品分别按照无烟煤、烟煤、水汽煤、贫瘦煤、动力煤加以分组，对没有具体分类的样品则单独列为一组。煤样不需要称样，依次直接装满NIR仪的进样器，然后把满出进样器的样品铲去，近红外分析仪自动记录存储样品的DDS光谱，得原始谱图，见附图1至附图8所示。由附图1至附图8所示原始光谱图可知，只有附图4显示的贫瘦煤样品的近红外漫反射DDS原始光谱的变化趋势比较接近或一致，其他煤种如无烟煤、烟煤、水汽煤、动力煤的原始谱图都出现比较明显差异的变化趋势。由于附图6至图11为近三个月集中采集的谱图，综合考虑分析样品贮藏时间和环境对其原始光谱存在的影响，本发明收集2008年煤样及时采集谱图，得图12至图13，可以看到，图12与图6基本一致，图13各种煤炭样品的近红外漫反射DDS原始光谱基本涵盖了前述无烟煤、烟煤、水汽煤、动力煤的原始谱图的变化趋势。这说明我们的样品贮存条件或光谱采集时间对原始光谱影响不明显，可以满足本发明研究需要。为更好考察煤炭原始光谱，我们对煤炭标准样品采用同样方式采集原始光谱，得图14至图15。并将所得标准样品原始光谱与相应煤种的原始谱图合在一起进行观察，得图16至图17。由图16至图17可以看到光谱变化趋势与前述相应煤种基本一致。
接着，本发明对图16至图17分别进行光谱预处理，得图18至图24。光谱预处理方式分别为：NONE+0011、NONE+1441、D+1441、D+1441、NONE+0011、D+1441、D+0011(其中，NONE表示未进行光谱预处理，即为原始谱图，D表示趋势变换法，1441表示以每4个光谱点为间隔作一阶导数的数学处理，0011表示不作任何导数处理。例如，D+1441为对原始光谱图进行趋势变换法预处理后，再以每4个光谱点为间隔作一阶导数进行数学处理)。由图18至图24可以看到经预处理后的谱图变化趋势更加明显、清晰，但同一煤种的各种不同处理方式所得效果均普遍存在明显不同的峰型和变化。总体看来，光谱预处理并未能获得比较一致的谱图。对前述采集到的各种煤原始光谱进行预处理，所得情况基本相同，以509个无烟煤样品原始光谱经(NONE+1441)处理图为例，处理结果请见附图25至附图27所示，其他样品预处理图在此不一一赘述(图略)。
S2是将收集到的样品原始光谱按其不同变化趋势分类，按照峰型和变化趋势相同或接近的图谱合在一起，进行归类，分别得到Y型原始光谱、W型原始光谱、P型原始光谱和X型原始光谱。本发明将前述采集到的原始光谱按其不同变化趋势分组，再将峰型和变化趋势相同或比较接近的图谱合在一起，进行归类，得谱图28至图31。图28的谱图标记为Y型原始光谱，代表大多数的水汽煤和部分烟煤、部分动力煤等煤种，共有374个样品。图29被标记为W型原始光谱，代表大多数的无烟煤、部分烟煤、部分动力煤等煤种，共有1367个样品。图30被标记为P型原始光谱，代表大多数的贫瘦煤、部分烟煤、部分动力煤等煤种，共有436个样品。图31被标记为X型原始光谱，代表非常见煤种，共有195个样品。
S03所述对S02所得样品的光谱数据进行数学处理是采用WinISI软件进行光谱分析和建立检测模型，将S02所得样品的光谱数据导入NIR仪，确定检测模型，打印检测结果。光谱预处理分别采用趋势变换法、标准正态变量转换法、多元离散校正、反相多元离散校正等方法的一种或多种，结合煤炭样品特性分析，最终确定最佳处理方法；回归校正方法采用逐步回归分析法(SMLR)、主成分分析法(PCA)和最小偏差分析法(PLS)，将数据降维，以消除众多信息共存中相互重叠的信息部分并最终作到对光谱的量化。
利用检测模型检测一组未知待测成分含量的样品中内水含量，再将NIR法所得的检测值与常规物理化学法检测值进行比较和评价。两种方法的比较结果用预测标准偏差(SEP,Standard Error of Prediction)和相应的决定系数(RSQp)或相关系数Rp进行衡量。
所述的检测模型的建立方法包括以下步骤：
S031.GH值分析：GH值即马氏距离，是光谱数据回归校正的得分三维图中，每个样品距离中心样品点的距离。GH值在近红外光谱数据分析中通常被设定为3.0，含义是标准变异单位的3倍，即就是近似于标准差(SD)的2.84倍，这意味着将有大约10％的样品的GH值会大于3.0。如果哪个样品的GH值大于3.0，该样品需剔除，另作分析。本发明采用主成分分析法(PCR)对Y型原始光谱、W型原始光谱、P型原始光谱和X型原始光谱分别进行聚类分析；结果见附图32至附图35所示，由附图32至35分析可知，四种类型原始谱图的GH值绝大多数小于3.0，Y型大于3.0的样品有8个，W型大于3.0的样品有10个，P型大于3.0的样品有13个，X大于3.0的样品共有9个，本发明将所述GH值大于3.0的样品剔除后，用GH值小于3.0的样品集分别建立相应类型谱图的定标(检测)模型，其中Y型366个，W型1357个，P型423个，X型186个。
S032.通过计算S031所述定标(检测)模型的SEC值和RSQ值；
利用WinISI软件中光谱处理和回归校正方法对上述四类原始光谱进行光谱处理和数据分析。光谱预处理分别采用趋势变换法、标准正态变量转换法、多元离散校正、反相多元离散校正等方法的一种或多种，结合煤炭样品特性分析，最终确定最佳处理方法；数学处理采用以每4个光谱点的间隔(Gap)作一阶导数(1441)或二阶导数(2441)方法，(0011)表示不做任何导数处理。回归校正方法采用逐步回归分析法(Stepwise Mutiple LinearRegression，SMLR)、主成分分析法(Principal Component Analysis，PCR)、最小偏差分析法(Partial Least Squares Regression,PLS)和偏最小偏差分析法(Modified PartialLeast Squares Regression，MPLS)计算所述定标(检测)模型的SEC值和RSQ值。
SEC是定标标准偏差(Standard Error of Calibration，SEC)，指通过建立的定标模型对定标样品集进行预测所得到的近红外分析值和传统化学法分析值的标准偏差，是回归读数和实际读数吻合程度的标志。SEC越低，说明近红外分析结果与传统分析结果越吻合，可信度越高。RSQ(R squared)为定标决定系数，是相关系数(Rc)的平方，指定标模型对定标样品集变异所能描述出的百分率，表示近红外分析值与常规方法分析值线性关系密切的程度。一个好的检测模型要求有低的SEC和高的RSQ(Rc)。根据本领域现有公知常识，理论上，具备最低SEC值和最高RSQ值的模型就应该是最佳的检测模型。但是本发明经过大量长期的实验总结发现，并非最低SEC值和最高RSQ值的模型就是最佳的检测模型，但最终能被定为最佳检测模型的一定具有较低的SEC值和较高的RSQ值。为此本发明对建立的全部检测模型逐一进行评价试验。以Y型原始谱图为例，其他可参照Y型进行实验总结，不一一赘述。煤炭样品的内水光谱处理和回归校正结果分别见表1所示(Y型原始谱图)：
表1 Y型原始谱图光谱处理和回归校正效果
表2 Y型检测模型的交互验证结果
S033.对检测模型进行评价实验确定最佳检测模型。由表1和表2可以看出，不论哪种类型的煤炭光谱，每一个煤质指标的检测模型，其交互验证结果的特征值--SECV值、(1-VR)值与相应模型的SEC值、RSQ值都有所不同，SECV值都略大于SEC值，(1-VR)值都略小于RSQ值。SECV值、(1-VR)值比SEC值、RSQ值更能反映检测模型将来用于预测未知样品的准确度，这是因为SEC值、RSQ值反映的是检测模型对定标样品集进行预测所得到的近红外分析值和传统化学法分析值的标准偏差、近红外分析值与常规方法分析值线性关系密切的程度；而SECV值、(1-VR)值反映的是在交互验证计算时，检测模型对没有参加定标建模的样品进行预测所得到的近红外分析值和传统化学法分析值的标准偏差、没有参加定标建模的样品近红外分析值与常规方法分析值线性关系密切的程度。由此可见，一个好的检测模型不但要求有低的SEC和高的RSQ(Rc)，更要求有低的SECV值和高的(1-VR)或Rv值。
不论哪种类型的煤炭光谱，每一个煤质指标的检测模型，一个拥有最低SECV值和最高(1-VR)或Rv值的检测模型，一定具有较低的SEC值和较高的RSQ值。本发明基于更多种除煤炭样品以外的样品检测研究中总结的结与本发明研究结果一致，本发明四种光谱类型的每个煤质指标具有最低SECV值和最高(1-VR)值的检测模型判定为最佳检测模型，并将主要信息总结分析于表3。
表3最佳检测模型的主要特征参数及其交互验证结果
由表3可见，利用WinISI软件对煤炭的近红外漫反射原始光谱进行多种光谱处理和回归校正，结果以采用标准正态变量转换(SNV)、趋势变换(D)、无散射(NONE)处理和一次光滑化(0011)、以4个光谱点的间隔(Gap)作一阶导数(1441)、二阶导数(2441)处理以及MPLS、PLS回归校正的组合效果最佳。所得最佳检测模型的主要特征值和相应的回归方法、光谱处理方式详见表3所列。各个煤质指标的最佳NIR检测模型有692项参数，其中包含一个常数项，本实施例以内水为例，检测模型的参数和方程曲线详见表4至表7，方程曲线(检测模型)如图36至图39所示，其它检测模型的参数和方程曲线图略。
表4 Y型原始光谱检测模型参数(内水％)
表5 W型原始光谱检测模型参数(内水％)
表6 P型原始光谱检测模型参数(内水％)
表7 X型原始光谱检测模型参数(内水％)
本实施例利用实施例确立的检测模型检测一组未知待测成分含量的煤炭样品中成分含量，再将NIR法所得的检测值与常规物理化学法检测值进行比较和评价。两种方法的比较结果用预测标准偏差(SEP，Standard Error of Prediction)和相应的决定系数(RSQp)或相关系数Rp来衡量。
随机选用一些煤炭样品，不需称样，依次直接装满NIR仪的进样器，启动扫描键，NIR仪自动记录存储样品光谱。在判断样品归属谱图类型后，点击相应的检测模型，就可以打印检测结果，快速完成煤炭样品的内水、灰分、挥发分、全硫、固定碳和可燃体的含量及发热量7项指标的检测工作。本应用试验Y、W、P、X四种类型采用的样品数、常规法得到的实验值(均值和标准偏差)、近红外法得到的扫描值(均值和标准偏差)，SEP和相应的RSQp、Rp值分别见表8所示。
表8各种煤炭采用本发明方法分析检测结果
由表8可见，本发明煤炭中内水、灰分、挥发分、全硫、固定碳和可燃体的含量及发热量7项指标检测方法的预测标准偏差、决定系数和相关系数等指标均与相应的定标和验证结果相近。具体常规法物理化学实验值、NIR法扫描值及两者之差值、现行标准方法的再现性限差见表9至表12，每类型每项目只列出10个样品的结果，以内水、发热量、固定碳为例，其他大量实验结果不能在此一一赘述。
表9 Y型检测模型检测结果—内水(％)
表10 W型检测模型初步应用结果—内水(％)
表11 P型检测模型初步应用结果--内水(％)
表12 X型检测模型初步应用结果—内水(％)
表13 Y型检测模型初步应用结果—发热量(MJ/Kg)
表14 W型检测模型初步应用结果—发热量(MJ/Kg)
表15 P型检测模型初步应用结果—发热量(MJ/Kg)
表16 X型检测模型初步应用结果—发热量(MJ/Kg)
表17 Y型检测模型初步应用结果—固定碳(％)
表18 W型检测模型初步应用结果—固定碳(％)
表19 P型检测模型初步应用结果—固定碳(％)
表20 X型检测模型初步应用结果—固定碳(％)
由上述代表性的表9至表20可见，本发明建立的煤炭中指标检测方法的准确度符合现行标准方法再现性限差的要求。操作非常简便快速，50秒左右就可以完成一个样品的扫描和内水、灰分、挥发分、全硫、固定碳和可燃体的含量及发热量7项指标检测，包括数据的输出。整个检测操作的过程中，无需称样，无需添加任何化学试剂，具有快速、便捷、无污染的特点。
本发明简单易行地实现一次性同时分析得到7项关键指标的检测并结合光谱图和指标检测结果判断识别煤炭样品是否掺杂掺假。本发明提供的新方法检测得到煤炭样品中7项关键成分以及做出判断的谱图依据只需几分钟，而且无需称样、无需使用化学试剂或高温、高压、大电流等测试条件，不会产生化学、生物或电磁污染，不会对操作人员和环境造成不良影响。就每一项单个的指标比较来说：水分的检测，与现行烘箱法相比：Y型原始谱图定标标准偏差为0.14％，定标相关系数为0.9969；交互验证标准差0.15％，交互验证相关系数0.9965；初步应用的标准偏差为0.17％，相关系数为0.995。W型原始谱图定标标准偏差为0.25％，定标相关系数为0.9021；交互验证标准差0.28％，交互验证相关系数0.8756；初步应用的标准偏差为0.23％，相关系数为0.875。P型原始谱图定标标准偏差为0.41％，定标相关系数为0.9814；交互验证标准差0.47％，交互验证相关系数0.9746；初步应用的标准偏差为0.83％，相关系数为0.944。X型原始谱图定标标准偏差为0.12％，定标相关系数为0.9981；交互验证标准差0.31％，交互验证相关系数0.9874；初步应用的标准偏差为0.13％，相关系数为0.997。
灰分的检测，与现行高温炉法相比：Y型原始谱图定标标准偏差为0.34％，定标相关系数为0.9955；交互验证标准差0.35％，交互验证相关系数0.9954；初步应用的标准偏差为0.77％，相关系数为0.980。W型原始谱图定标标准偏差为0.51％，定标相关系数为0.9851；交互验证标准差0.74％，交互验证相关系数0.9687；初步应用的标准偏差为0.57％，相关系数为0.982。P型原始谱图定标标准偏差为1.29％，定标相关系数为0.9785；交互验证标准差1.62％，交互验证相关系数0.9659；初步应用的标准偏差为1.32％，相关系数为0.970。X型原始谱图定标标准偏差为1.48％，定标相关系数为0.9879；交互验证标准差2.05％，交互验证相关系数0.9763；初步应用的标准偏差为3.88％，相关系数为0.926。
挥发分检测，与现行高温炉法相比：Y型原始谱图定标标准偏差为0.27％，定标相关系数为0.9922；交互验证标准差0.29％，交互验证相关系数0.9911；初步应用的标准偏差为0.37％，相关系数为0.986。W型原始谱图定标标准偏差为0.22％,定标相关系数为0.9618；交互验证标准差0.24％,交互验证相关系数0.9521；初步应用的标准偏差为0.25％，相关系数为0.928。P型原始谱图定标标准偏差为0.93％,定标相关系数为0.9889；交互验证标准差1.05％，交互验证相关系数0.9858；初步应用的标准偏差为0.89％，相关系数为0.975。X型原始谱图定标标准偏差为0.59％,定标相关系数为0.9991；交互验证标准差0.65％，交互验证相关系数0.9989；初步应用的标准偏差为1.61％，相关系数为0.994。
全硫的检测，与高温燃烧红外吸收法相比：Y型原始谱图定标标准偏差为0.06％，定标相关系数为0.9838；交互验证标准差0.07％，交互验证相关系数0.9822；初步应用的标准偏差为0.08％，相关系数为0.976。W型原始谱图定标标准偏差为0.04％，定标相关系数为0.9496；交互验证标准差0.05％，交互验证相关系数0.9379；初步应用的标准偏差为0.05％，相关系数为0.917。P型原始谱图定标标准偏差为0.17％，定标相关系数为0.9572；交互验证标准差0.21％，交互验证相关系数0.9330；初步应用的标准偏差为0.22％，相关系数为0.927。X型原始谱图定标标准偏差为0.11％，定标相关系数为0.8988；交互验证标准差0.16％，交互验证相关系数0.7412；初步应用的标准偏差为0.12％，相关系数为0.849。
发热量的检测，与现行氧弹燃烧法相比：Y型原始谱图定标标准偏差为0.12MJ/Kg，定标相关系数为0.9960；交互验证标准差0.14MJ/Kg，交互验证相关系数0.9951；初步应用的标准偏差为0.26MJ/Kg，相关系数为0.983。W型原始谱图定标标准偏差为0.22MJ/Kg，定标相关系数为0.9815；交互验证标准差0.31MJ/Kg，交互验证相关系数0.9610；初步应用的标准偏差为0.27MJ/Kg，相关系数为0.972。P型原始谱图定标标准偏差为0.48MJ/Kg，定标相关系数为0.9734；交互验证标准差0.60MJ/Kg，交互验证相关系数0.9594；初步应用的标准偏差为0.88MJ/Kg，相关系数为0.971。X型原始谱图定标标准偏差为0.45％，定标相关系数为0.9830；交互验证标准差0.72MJ/Kg，交互验证相关系数0.9582；初步应用的标准偏差为1.20MJ/Kg，相关系数为0.901。
固定碳的检测，与现有烘箱、高温炉法相比：Y型原始谱图定标标准偏差为0.29％，定标相关系数为0.9948；交互验证标准差0.30％，交互验证相关系数0.9931；初步应用的标准偏差为0.43％，相关系数为0.986。W型原始谱图定标标准偏差为2.20％，定标相关系数为0.8499；交互验证标准差2.49％，交互验证相关系数0.8036；初步应用的标准偏差为2.11％，相关系数为0.835。P型原始谱图定标标准偏差为1.34％，定标相关系数为0.9731；交互验证标准差1.65％，交互验证相关系数0.9587；初步应用的标准偏差为3.10％，相关系数为0.941。X型原始谱图定标标准偏差为1.66％，定标相关系数为0.9883；交互验证标准差2.16％,交互验证相关系数0.9802，初步应用的标准偏差为2.64％，相关系数为0.972。
可燃体的检测，与现有烘箱、高温炉法相比：Y型原始谱图定标标准偏差为0.34％,定标相关系数为0.9954；交互验证标准差0.37％，交互验证相关系数0.9948；初步应用的标准偏差为0.60％，相关系数为0.988。W型原始谱图定标标准偏差为2.35％，定标相关系数为0.8414；交互验证标准差2.62％，交互验证相关系数0.7981；初步应用的标准偏差为2.19％，相关系数为0.819。P型原始谱图定标标准偏差为1.32％，定标相关系数为0.9685；交互验证标准差1.51％，交互验证相关系数0.9586；初步应用的标准偏差为4.67％，相关系数为0.867。X型原始谱图定标标准偏差为1.34％，定标相关系数为0.9874；交互验证标准差1.92％，交互验证相关系数0.9738；初步应用的标准偏差为3.32％，相关系数为0.930。
Claims (1)
1.一种同时检测煤炭中7项指标并识别煤炭是否掺杂掺假的方法，其特征在于，包括以下步骤：
S01.收集和制备若干个煤炭样品，常规方法分别测定每个样品的内水、灰分、挥发分、全硫、固定碳和可燃体的含量及发热量7项指标；
S02.用数字光栅近红外漫反射仪扫描收集所述煤炭样品的光谱数据和曲线；将收集到的样品原始光谱按照将峰型和变化趋势相同或接近的图谱合在一起的方式，进行归类，分别得到Y型原始光谱、W型原始光谱、P型原始光谱和X型原始光谱；
S03.对步骤S02所得样品的光谱数据进行处理，经回归计算获得7项指标的定标方程，修正和验证内水后建立检测模型；
S04.将待测煤炭样品依次直接装满数字光栅近红外漫反射仪的进样器，不需称样，启动扫描键，采用数字光栅系统进行扫描，数字光栅近红外漫反射仪自动记录存储样品光谱，将所得待测煤炭样品的光谱与用数字光栅近红外漫反射仪扫描收集的不同煤种的标准煤炭样品的近红外原始光谱进行比较，确定待测煤炭样品归属谱图类型；选择相应的检测模型，获得检测结果；
所述检测模型的建立方法包括以下步骤：
S031.GH值分析：采用主成分分析法对Y型原始光谱、W型原始光谱、P型原始光谱和X型原始光谱分别进行聚类分析，将所述GH值大于3.0的样品剔除后，用GH值小于3.0的样品集分别建立相应类型谱图的检测模型；
S032.计算步骤S031所述检测模型的SEC值和RSQ值；
S033.对检测模型进行评价实验确定最佳检测模型；其中，具有最低SECV值和最高1-VR值的检测模型判定为最佳检测模型。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410579312.2A CN104390928B (zh) | 2014-10-24 | 2014-10-24 | 一种掺杂掺假煤炭的近红外光谱识别方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410579312.2A CN104390928B (zh) | 2014-10-24 | 2014-10-24 | 一种掺杂掺假煤炭的近红外光谱识别方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104390928A CN104390928A (zh) | 2015-03-04 |
| CN104390928B true CN104390928B (zh) | 2018-03-20 |
Family
ID=52608858
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410579312.2A Expired - Fee Related CN104390928B (zh) | 2014-10-24 | 2014-10-24 | 一种掺杂掺假煤炭的近红外光谱识别方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104390928B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106990066B (zh) * | 2017-02-22 | 2020-03-13 | 中国科学院遥感与数字地球研究所 | 一种鉴别煤炭种类的方法及装置 |
| CN108469425A (zh) * | 2018-02-05 | 2018-08-31 | 管维新 | 一种近红外光谱原理分析煤质的方法 |
| CN109211831A (zh) * | 2018-09-26 | 2019-01-15 | 天津出入境检验检疫局化矿金属材料检测中心 | 一种无烟煤、烟煤和褐煤的鉴别方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN201207037Y (zh) * | 2008-05-16 | 2009-03-11 | 辽宁东信盛大科技有限公司 | 近红外煤质在线检测分析仪 |
| CN202285003U (zh) * | 2011-08-23 | 2012-06-27 | 中煤西安设计工程有限责任公司 | 近红外煤粉水分在线测量装置 |
| CN103728268A (zh) * | 2013-12-17 | 2014-04-16 | 沈阳工程学院 | 近红外煤质在线测量装置 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001116689A (ja) * | 1999-10-15 | 2001-04-27 | Japan Science & Technology Corp | 石炭の水分値推定方法 |
| JP4416896B2 (ja) * | 2000-02-02 | 2010-02-17 | 東北電力株式会社 | 石炭成分分析装置センサーの標準板及び該標準板を用いた器差補正方法 |
| EP2965053A1 (en) * | 2013-03-07 | 2016-01-13 | Siemens Aktiengesellschaft | Systems and methods for boosting coal quality measurement statement of related cases |
-
2014
- 2014-10-24 CN CN201410579312.2A patent/CN104390928B/zh not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN201207037Y (zh) * | 2008-05-16 | 2009-03-11 | 辽宁东信盛大科技有限公司 | 近红外煤质在线检测分析仪 |
| CN202285003U (zh) * | 2011-08-23 | 2012-06-27 | 中煤西安设计工程有限责任公司 | 近红外煤粉水分在线测量装置 |
| CN103728268A (zh) * | 2013-12-17 | 2014-04-16 | 沈阳工程学院 | 近红外煤质在线测量装置 |
Non-Patent Citations (2)
| Title |
|---|
| 基于近红外光谱分析的电厂煤质分类快速检测方法;丁仁杰 等;《电站系统工程》;20070131;第23卷(第1期);32-34 * |
| 应用NIRS分析技术快速检测煤炭质量;苏彩珠 等;《检验检疫科学》;20071231;第17卷(第6期);34-35 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104390928A (zh) | 2015-03-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104390927B (zh) | 煤炭样品中灰分的快速检测方法 | |
| CN101893561B (zh) | 一种润滑油新油质量近红外光谱快速测定方法 | |
| CN102879353B (zh) | 近红外检测花生中蛋白质组分含量的方法 | |
| CN102590129B (zh) | 近红外检测花生中氨基酸含量的方法 | |
| CN101769866B (zh) | 一种牛奶成分的检测装置和方法 | |
| CN104749132A (zh) | 一种测定面粉中偶氮甲酰胺的含量的方法 | |
| CN109540837B (zh) | 近红外快速检测苎麻叶片木质纤维素含量的方法 | |
| CN112179871A (zh) | 一种酱类食品中己内酰胺含量无损检测的方法 | |
| Fodor et al. | Estimation of middle distillate fuel properties by FT-IR | |
| CN108760647A (zh) | 一种基于可见/近红外光谱技术的小麦带菌量线检测方法 | |
| CN104390928B (zh) | 一种掺杂掺假煤炭的近红外光谱识别方法 | |
| CN103134763B (zh) | 由红外光谱预测原油密度的方法 | |
| CN110376154A (zh) | 基于光谱校正的水果在线检测方法及系统 | |
| CN104297201A (zh) | 一种快速、准确定量调和油中各种油份比例的检测方法 | |
| CN204255854U (zh) | 一种同时检测煤炭样品中7项指标的系统 | |
| CN104568811B (zh) | 煤炭样品中发热量的新型快速检测方法 | |
| CN110231306A (zh) | 一种无损、快速测定奇亚籽蛋白质含量的方法 | |
| CN102954946B (zh) | 由红外光谱测定原油硫含量的方法 | |
| CN104849234A (zh) | 基于近红外光谱分析吡虫啉原药主成分含量的测定方法 | |
| CN105092526A (zh) | 一种基于近红外光谱技术的二元掺伪芝麻油含量的快速测定方法 | |
| CN104535525B (zh) | 煤炭样品中的固定碳的快速检测方法 | |
| CN104568821B (zh) | 一种同时检测煤炭样品中7项指标的方法 | |
| CN104390936B (zh) | 煤炭样品中挥发分的快速检测方法 | |
| CN104374735B (zh) | 煤炭样品中全硫的新型快速检测方法 | |
| US7420170B2 (en) | Fourier transform infrared (FTIR) chemometric method to determine cetane number of diesel fuels containing fatty acid alkyl ester additives |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Su Caizhu Inventor after: Li Peishen Inventor after: Zheng Jianguo Inventor after: Zhang Haifeng Inventor after: Tang Yong Inventor after: Ruan Xiongjie Inventor after: Qiu Minmin Inventor after: Zhang Weimin Inventor after: Xie Wenxi Inventor after: Ma Lina Inventor after: Chen Ruiyang Inventor before: Su Caizhu Inventor before: Li Guowei Inventor before: Liu Erlong Inventor before: Zheng Jianguo Inventor before: Zheng Shuyun Inventor before: Qiu Minmin Inventor before: Cai Yingjun Inventor before: Yao Baihui |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180320 Termination date: 20181024 |