CN108279217A - 一种基于太赫兹时域光谱的煤岩判别方法 - Google Patents
一种基于太赫兹时域光谱的煤岩判别方法 Download PDFInfo
- Publication number
- CN108279217A CN108279217A CN201810404583.2A CN201810404583A CN108279217A CN 108279217 A CN108279217 A CN 108279217A CN 201810404583 A CN201810404583 A CN 201810404583A CN 108279217 A CN108279217 A CN 108279217A
- Authority
- CN
- China
- Prior art keywords
- coal
- sample
- terahertz
- ash content
- domain spectroscopy
- 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.)
- Granted
Links
- 239000003245 coal Substances 0.000 title claims abstract description 213
- 238000001328 terahertz time-domain spectroscopy Methods 0.000 title claims abstract description 63
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000001228 spectrum Methods 0.000 claims abstract description 71
- 238000000862 absorption spectrum Methods 0.000 claims abstract description 70
- 238000005065 mining Methods 0.000 claims abstract description 15
- 230000009467 reduction Effects 0.000 claims abstract description 12
- 238000005516 engineering process Methods 0.000 claims abstract description 8
- 230000003595 spectral effect Effects 0.000 claims abstract description 8
- 239000002956 ash Substances 0.000 claims description 89
- 239000011159 matrix material Substances 0.000 claims description 30
- 235000002918 Fraxinus excelsior Nutrition 0.000 claims description 25
- 238000004611 spectroscopical analysis Methods 0.000 claims description 25
- 239000011435 rock Substances 0.000 claims description 20
- 238000000605 extraction Methods 0.000 claims description 13
- 238000012360 testing method Methods 0.000 claims description 12
- 239000000203 mixture Substances 0.000 claims description 8
- 239000000047 product Substances 0.000 claims description 8
- 230000005540 biological transmission Effects 0.000 claims description 7
- 238000009499 grossing Methods 0.000 claims description 7
- 238000010521 absorption reaction Methods 0.000 claims description 6
- 239000006227 byproduct Substances 0.000 claims description 6
- 239000000284 extract Substances 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 5
- RHZUVFJBSILHOK-UHFFFAOYSA-N anthracen-1-ylmethanolate Chemical class C1=CC=C2C=C3C(C[O-])=CC=CC3=CC2=C1 RHZUVFJBSILHOK-UHFFFAOYSA-N 0.000 claims description 4
- 230000008033 biological extinction Effects 0.000 claims description 4
- 239000002802 bituminous coal Substances 0.000 claims description 4
- 230000014509 gene expression Effects 0.000 claims description 4
- 239000003077 lignite Substances 0.000 claims description 4
- 230000003287 optical effect Effects 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 4
- 238000000513 principal component analysis Methods 0.000 claims description 3
- 238000003860 storage Methods 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 2
- 238000007689 inspection Methods 0.000 claims 1
- 238000001514 detection method Methods 0.000 description 11
- 238000004458 analytical method Methods 0.000 description 5
- 239000000126 substance Substances 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 239000000463 material Substances 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 239000010883 coal ash Substances 0.000 description 3
- -1 geological condition Substances 0.000 description 3
- 229920000573 polyethylene Polymers 0.000 description 3
- 239000004698 Polyethylene Substances 0.000 description 2
- 239000003830 anthracite Substances 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 239000002817 coal dust Substances 0.000 description 2
- 239000000470 constituent Substances 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000009776 industrial production Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 229910001218 Gallium arsenide Inorganic materials 0.000 description 1
- 241000258241 Mantis Species 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000003818 cinder Substances 0.000 description 1
- 238000009841 combustion method Methods 0.000 description 1
- 238000000151 deposition Methods 0.000 description 1
- 239000003989 dielectric material Substances 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000010903 husk Substances 0.000 description 1
- 239000000843 powder Substances 0.000 description 1
- 229910052594 sapphire Inorganic materials 0.000 description 1
- 239000010980 sapphire Substances 0.000 description 1
- 229910052710 silicon Inorganic materials 0.000 description 1
- 239000010703 silicon Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/3581—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using far infrared light; using Terahertz radiation
- G01N21/3586—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using far infrared light; using Terahertz radiation by Terahertz time domain spectroscopy [THz-TDS]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N2021/3595—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using FTIR
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Health & Medical Sciences (AREA)
- Toxicology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种基于太赫兹时域光谱的煤炭灰分判别方法。尤其适用于采煤机自主调高和智能化开采技术领域使用。制作纯煤样品和不同比例煤岩混合样品若干，采用太赫兹时域光谱仪采集参考信号和不同灰分煤的太赫兹时域光谱；利用快速傅里叶变换将煤太赫兹时域光谱转换成频域光谱；从太赫兹频域光谱中提取出煤的吸收光谱；利用PCA对吸收光谱数据进行降维，随后利用非线性拟合对处理后的煤光谱进行建模，确立煤炭种类、煤中灰分和其太赫兹吸收光谱的关系曲线，构建基于太赫兹光谱的煤炭灰分模型，从而实现基于太赫兹光谱技术的煤炭灰分判别。其步骤简单，识别效率高，识别准确，实现采煤机的截割状态的准确识别，自动控制采煤机滚筒调高。
Description
技术领域
本发明涉及一种基于太赫兹时域光谱的煤炭灰分判别方法，尤其适用于采煤机自主调高和智能化开采技术领域使用。
背景技术
煤的灰分含量是煤工业分析中的一项重要指标，煤中灰分的含量直接影响了煤的发热量和煤质，也影响了商品煤的售价。传统的灰分分析方法为燃烧法即判别燃烧后煤中的残余部分，这是一种离线的检测方法，无法满足日常生产和当前的工业生产需求。此外，由于煤中的灰分也来自与矸石和岩层，因此准确的识别出煤中的灰分，也可以对煤岩界面识别和采煤机滚筒调高起到重要的指导意义。
为了解决上述问题，提出了一种基于太赫兹谱的煤炭灰分判别方法。由于太赫兹是一个比较特殊的频段，其太赫兹光子能级和许多电介质以及材料的振动能级相当，因此当太赫兹波照射煤岩后其光子能量会被吸收，其光强将会产生变化从而形成太赫兹光谱，通过研究太赫兹光谱将能获取检测物质的物理信息(如吸收光谱)，理论实验表明煤中灰分含量的变化会导致其太赫兹吸收谱产生变化，且灰分含量是导致其变化的主要因素，因此可以用太赫兹吸收光谱表征煤炭中的灰分含量。
相对于其它的光谱技术，太赫兹光谱可以直接得到电场的相位，而不需要利用Kramers-Kronig方程求解相位角，极大的减少了计算量和复杂度。因此该方法是一种较新颖的研究方法，且展现出了巨大的潜力。
然而与其它的标准物质不同，煤岩介质在太赫兹频段没有吸收峰，这也表明煤和岩的组成成分复杂且没有主导吸收的物质，直接使用太赫兹吸收光谱进行煤炭灰分研究得到的结果并不理想。此外，研究表明矿井中煤和岩的赋存条件、地质情况、成分变化等因素都会导致太赫兹光谱产生变化。
发明内容
针对上述技术问题，提供一种步骤简单，使用方便，效率高，判别效果好的基于太赫兹时域光谱的煤岩识别方法。
为实现上述技术目的，本发明的基于太赫兹时域光谱的煤炭灰分判别方法，其特征在于包括如下步骤：
a.挑选不同品种的煤，利用挑选出的煤加入纯岩样制作不同灰分煤样品，利用太赫兹时域光谱仪为对空气直接检测得到无样本时的信号作为参考信号，再利用太赫兹时域光谱仪采集所有煤样品的太赫兹时域光谱数据并分类存储起来，无样本时的信号为对空气直接检测得到的信号，所述的太赫兹时域光谱数据有200维，并将所有煤样品按照中国煤炭分类标准分成：无烟煤、烟煤、褐煤几类；
b.利用快速傅里叶变换将所有煤样品的太赫兹时域光谱转换为太赫兹频域光谱；从所有转换后的太赫兹频域光谱中提取吸收光谱，并对吸收光谱数据进行平滑和加窗预处理，选取数据噪声超过1THz以上的煤样品吸收光谱，将选取的煤样品吸收光谱中0.4-1.0THz频段光谱数据作为优选的吸收光谱；
c.将所有煤样品，包括不同灰分煤样品和纯煤样品的0.4-1.0THz频段吸收光谱进行标注分类从而生成不同灰分煤样品和纯煤样品的煤光谱数据库，利用主成分分析/PCA技术对煤光谱数据库中的所有数据进行降维和特征提取，使获得的太赫兹时域光谱数据从200维降到了20维；
d.利用非线性拟合法对所有降维后的0.4-1.0THz频段吸收光谱进行建模，根据所有煤样品的三类信息：煤炭种类、煤中灰分和其太赫兹吸收光谱确立三类信息的关系曲线，从而构建基于太赫兹光谱的煤炭灰分模型；
e.根据实际需要采集被测区域的煤样品后利用太赫兹时域光谱仪获取煤样品的太赫兹时域光谱，重复步骤b、c、d，将被测区域煤样品的光谱数据带入太赫兹光谱的煤炭灰分模型从而对煤样品中的灰分进行判别。
太赫兹时域光谱采集装置为天津大学组建的系统设备。
每次采集煤样品的太赫兹时域光谱时选取3-5个不同的点位进行采集，每个点位重复采集3-4次光谱，所述煤样品的为粉末状或块状。
在检测过程中，采用太赫兹时域光谱系统的透射扫描模块获取以干燥空气为背景的参考信号时域光谱和不同灰分含量煤样品的太赫兹时域光谱，干燥空气的空气湿度低于5％。
所述利用快速傅里叶变换将所有煤样品的太赫兹时域光谱转换为太赫兹频域光谱，根据太赫兹频域光谱提取折射率光谱n(ω)和吸收光谱α(ω)采用THz光学参数提取模型：
T(ω)＝Esam(ω)/Eref(ω)
其中，T(ω)为煤样品的透射系数光谱，Eref(ω)为参考信号的太赫兹频谱，Esam(ω)为煤样品太赫兹频谱，ω为角频率，d为煤样品的厚度，c为光速，n(ω)为煤样品的折射率，α(ω)为煤样品的吸收系数，ρ(ω)项为煤样品信号与参考信号幅度的比值，项为测试样品信号与参考信号的相位差，k(ω)为测试样品复折射率的虚部，也称为消光系数；测试样品信号为不同灰分的煤炭信号，直接测出来的是时域信号，通过傅里叶变换到频域，然后利用上述公式，可以得到煤的吸收光谱；经傅里叶变换后的太赫兹频谱范围为0.1～3.5THz。
所述主成分分析方法为：
1)设不同灰分煤样品的吸收光谱为矩阵X，矩阵X包括n个不同灰分含量的煤样品和p个吸收光谱数据，其中n个不同灰分含量的煤样品为列向量，p个吸收光谱数据向量为行向量，表示同一样本的光谱数据维度：
式中，ei＝(xi1,xi2,L,xip)T∈Rp,(i＝1,2,L,n)为第i个样本，xj＝(x1j,x2j,L,xpj)T∈Rp,(j＝1,2,L,p)为第j个变量；原始矩阵X可以由e向量矩阵和x向量矩阵内积表示；代表原始矩阵X由e向量矩阵和x向量矩阵内积表示；
2)将矩阵X变换改写成：
F＝AXT
改写后的F矩阵的列向量表示不同灰分含量煤样品的个数，行向量表示同一煤样品的光谱数据维度；将p个吸收光谱数据向量x1,x2,…,xp变成新向量F1,F2,…,Fp，将矩阵X以2个矩阵向量内积的形式表示，即：第i个样本ei＝(xi1,xi2,L,xip)T∈Rp,(i＝1,2,L,n)变换后的向量为：F(i)＝(Fi1,Fi2,L,Fip)T∈Rp,(i＝1,2,L,n)，其中，式中R为自然数，表示取值范围，a为A矩阵中的行向量，T为矩阵转置；
3)利用公式：确定A矩阵中的行向量aj，i表示第i个样本，Fij表示第i个样本对应的Fj，使与aj对应的Fj的方差最大，式中，Var(Fj)为煤样品变换后的值与均值的均方差，是样本变换后Fj的均值，V是X的协方差矩阵，V是实对称矩阵，假设V有p个非负特征值λ1,λ2,Lλp，且按降序排列，特征值所对应的特征向量分别为a1,a2,Lap，且各个特征向量之间正交，则公式可变成：
4)设F1＝Xa1为第一主成份，以此类推计算出所有主成份，将主成份按降序排列，因此取前几个对原始数据贡献大的20个主成份，使样本维数降低。
所述加窗预优选后的吸收光谱数据，提取信噪比较高的0.4-1THz频段光谱数据；所述平滑处理利用Savitzky-Golay平滑方法处理光谱曲线，消除吸收光谱中的系统和环境噪声。
将太赫兹时域光谱仪设置在采煤机滚筒上实时采集采煤机滚筒附近物体的灰分，当检测到的灰分若灰分值超过设定阀值时，说明当前采煤机滚筒切割到岩层，需要及时调节采煤机滚筒的高度。
有益效果：
本发明采用太赫兹时域光谱仪采集参考信号和不同灰分含量煤的太赫兹时域光谱；利用快速傅里叶变换将煤太赫兹时域光谱转换成频域光谱；从太赫兹频域光谱中提取出煤的吸收光谱。采用加窗、平滑和PCA等技术对光谱数据进行预处理，随后利用非线性拟合对处理后的煤光谱进行建模与判别，从而实现基于太赫兹光谱技术的煤炭灰分判别。
本发明专利研究利用太赫兹光谱技术和模式识别方法监测煤中的灰分含量，提出了基于太赫兹时域光谱的煤炭灰分含量在线监测方法，相对于传统的离线检测方法而言，这是一种高效、在线、无损的检测方法，可以满足日常生产和当前的工业生产需求。此外，由于煤中的灰分也来自与矸石和岩层，因此准确的识别出煤中的灰分，也可以对煤岩界面识别和采煤机滚筒调高起到重要的指导意义。其检测效果好，检测效率高，具有广泛的实用性。
附图说明
图1为本发明的不同灰分含量的煤炭样品示意图；
图2为本发明的太赫兹时域光谱采集系统结构示意图；
图3为本发明的不同灰分时的煤炭吸收光谱示意图；
图4为本发明的煤中灰分含量与其吸收光谱的关系示意图；
图5为本发明的煤炭灰分含量预测图。
图中：1、发射天线，2、接收天线，3、飞秒激光器，4、样本检测仓，5、延时线，6、分束镜。
具体实施方式
下面通过具体实施例对本发明作进一步详述，以下实施例只是描述性的，不是限定性的，不能以此限定本发明的保护范围
本发明基于太赫兹时域光谱的煤岩识别方法，使用的太赫兹时域光谱采集装置为天津大学组建的系统设备：如图2所示，太赫兹时域光谱仪包括用于产生光源的飞秒激光器3、产生太赫兹脉冲的发射天线1、样本检测仓4、收太赫兹脉冲的接收天线2以及调整飞秒激光和太赫兹脉冲之间时间差的延时线5五个部分，并通过分束镜6将光束相互传输；所述的样本检测仓位4于发射天线1和接收天线2之间，发射天线太赫兹光垂直作用到样本检测仓4的煤岩样品上；飞秒激光器为美国相干公司的Mantis飞秒激光器，中心波长800nm，带宽大于70nm，重复频率80MHz；发射天线材质为高阻砷化镓，接收天线材质为蓝宝石上硅,系统频率范围：0.1-3.5THz。
具体步骤如下：
a.如图1所示，挑选不同品种的煤，利用挑选出的煤加入纯岩样制作不同灰分煤样品，利用太赫兹时域光谱仪为对空气直接检测得到无样本时的信号作为参考信号，再利用太赫兹时域光谱仪采集所有煤样品的太赫兹时域光谱数据并分类存储起来，每次采集煤样品的太赫兹时域光谱时选取3-5个不同的点位进行采集，每个点位重复采集3-4次光谱，所述煤样品的为粉末状或块状；无样本时的信号为对空气直接检测得到的信号，所述的太赫兹时域光谱数据有200维，并将所有煤样品按照中国煤炭分类标准分成：无烟煤、烟煤、褐煤几类；采用太赫兹时域光谱系统的透射扫描模块获取以干燥空气为背景的参考信号时域光谱和不同灰分含量煤样品的太赫兹时域光谱，干燥空气的空气湿度低于5％；
b.利用快速傅里叶变换将所有煤样品的太赫兹时域光谱转换为太赫兹频域光谱；从所有转换后的太赫兹频域光谱中提取吸收光谱，并对吸收光谱数据进行平滑和加窗预处理，选取数据噪声超过1THz以上的煤样品吸收光谱，将选取的煤样品吸收光谱中0.4-1.0THz频段光谱数据作为优选的吸收光谱；所述加窗预优选后的吸收光谱数据，提取信噪比较高的0.4-1THz频段光谱数据；所述平滑处理利用Savitzky-Golay平滑方法处理光谱曲线，消除吸收光谱中的系统和环境噪声；
利用快速傅里叶变换将所有煤样品的太赫兹时域光谱转换为太赫兹频域光谱，根据太赫兹频域光谱提取折射率光谱n(ω)和吸收光谱α(ω)采用THz光学参数提取模型：
T(ω)＝Esam(ω)/Eref(ω)
其中，T(ω)为煤样品的透射系数光谱，Eref(ω)为参考信号的太赫兹频谱，Esam(ω)为煤样品太赫兹频谱，ω为角频率，d为煤样品的厚度，c为光速，n(ω)为煤样品的折射率，α(ω)为煤样品的吸收系数，ρ(ω)项为煤样品信号与参考信号幅度的比值，项为测试样品信号与参考信号的相位差，k(ω)为测试样品复折射率的虚部，也称为消光系数；测试样品信号为不同灰分的煤炭信号，直接测出来的是时域信号，通过傅里叶变换到频域，然后利用上述公式，可以得到煤的吸收光谱；经傅里叶变换后的太赫兹频谱范围为0.1～3.5THz；
c.将所有煤样品，包括不同灰分煤样品和纯煤样品的0.4-1.0THz频段吸收光谱进行标注分类从而生成不同灰分煤样品和纯煤样品的煤光谱数据库，利用主成分分析/PCA技术对煤光谱数据库中的所有数据进行降维和特征提取，使获得的太赫兹时域光谱数据从200维降到了20维；
d.利用非线性拟合法对所有降维后的0.4-1.0THz频段吸收光谱进行建模，根据所有煤样品的三类信息：煤炭种类、煤中灰分和其太赫兹吸收光谱确立三类信息的关系曲线，从而构建基于太赫兹光谱的煤炭灰分模型；
e.根据实际需要采集被测区域的煤样品后利用太赫兹时域光谱仪获取煤样品的太赫兹时域光谱，重复步骤b、c、d，将被测区域煤样品的光谱数据带入太赫兹光谱的煤炭灰分模型从而对煤样品中的灰分进行判别。
将太赫兹时域光谱仪设置在采煤机滚筒上实时采集采煤机滚筒附近物体的灰分，当检测到的灰分若灰分值超过设定阀值时，说明当前采煤机滚筒切割到岩层，需要及时调节采煤机滚筒的高度。
实施例一：
步骤(1)分别从不同品种的煤岩中挑选并制作不同灰分含量的煤样品，采用太赫兹时域光谱仪采集参考信号和煤岩样品的太赫兹时域光谱数据，所述的数据包括：无烟煤、烟煤、褐煤、砂岩和页岩；
在步骤1中，所采用的煤岩样品来自全国各地的矿区，将大块的煤块和岩块经过自然干燥、破碎、通过80目筛后，将岩粉混入煤粉中，制作不同灰分含量的煤样品；
在步骤1中，实施本发明方法所采用的太赫兹光谱采集装置包括了飞秒激光器、发射天线、样本检测仓、接收天线以及延时线五个部分；
在步骤1中，样品仓采用长度为50mm、宽度50mm、高度为30mm、厚度0.5mm的聚乙烯薄片制成。选用聚乙烯材料是因为在太赫兹波段聚乙烯可以看成透明介质；
在步骤1中，样品仓可通过定时旋转的方式倒出原有的煤样品，然后复位装入新的煤样品，样品仓位于太赫兹发射光路上，使太赫兹光垂直作用到煤样品上；
步骤(2)利用快速傅里叶变换将太赫兹时域光谱转换为太赫兹频域光谱；从太赫兹频域光谱中提取吸收光谱，煤中的灰分含量会影响到其吸收光谱的特性曲线，如图3所示。对吸收光谱进行平滑和加窗预处理，经过比较后选取0.4-1.0THz频段光谱数据；
作为优选，在步骤(2)时，通过快速傅里叶变换将参考信号和样品信号转到频域，经傅里叶变换后的太赫兹频谱范围为0.1～3.5THz。基于太赫兹频域光谱提取透射光谱T(ω)、折射率光谱n(ω)和吸收光谱α(ω)采用THz光学参数提取模型，计算公式如下：
T(ω)＝Esam(ω)/Eref(ω)
其中，透射系数为T(ω)，Eref(ω)和Esam(ω)分别为参考信号太赫兹频谱和样品太赫兹频谱，样品的厚度为d，c为光速，n(ω)为样品的折射率，α(ω)为样品的吸收系数，ρ(ω)项为测试样品信号与参考信号幅度的比值，项为测试样品信号与参考信号的相位差，k(ω)为测试样品复折射率的虚部，也称为消光系数。
作为优选，在步骤(2)时，通过加窗技术对吸收光谱进行优选，提取信噪比较高的0.4-1THz频段光谱数据，并利用Savitzky-Golay平滑方法处理光谱曲线，消除系统和环境噪声。
步骤(3)将所有的吸收光谱数据X标注为不同灰分含量的光谱数据，煤炭中不同的灰分含量与吸收光谱存在一定的对应关系，然而与其它的标准物质不同，煤岩介质在太赫兹频段没有吸收峰，这也表明煤和岩的组成成分复杂且没有主导吸收的物质，直接使用太赫兹吸收光谱进行煤炭灰分研究得到的结果并不理想。此外，研究表明矿井中煤和岩的赋存条件、地质情况、成分变化等因素都会导致太赫兹光谱产生变化。因此，采用PCA技术对光谱数据进行降维和特征提取。
在步骤3中，将不同灰分含量的煤吸收光谱数据X通过PCA进行降维和特征提取。
在步骤3中，主成分分析方法为：
1)设不同灰分煤样品的吸收光谱为矩阵X，矩阵X包括n个不同灰分含量的煤样品和p个吸收光谱数据，其中n个不同灰分含量的煤样品为列向量，p个吸收光谱数据向量为行向量，表示同一样本的光谱数据维度：
式中，ei＝(xi1,xi2,L,xip)T∈Rp,(i＝1,2,L,n)为第i个样本，xj＝(x1j,x2j,L,xpj)T∈Rp,(j＝1,2,L,p)为第j个变量；原始矩阵X可以由e向量矩阵和x向量矩阵内积表示；代表原始矩阵X由e向量矩阵和x向量矩阵内积表示；
2)将矩阵X变换改写成：
F＝AXT
改写后的F矩阵的列向量表示不同灰分含量煤样品的个数，行向量表示同一煤样品的光谱数据维度；将p个吸收光谱数据向量x1,x2,…,xp变成新向量F1,F2,…,Fp，将矩阵X以2个矩阵向量内积的形式表示，即：第i个样本ei＝(xi1,xi2,L,xip)T∈Rp,(i＝1,2,L,n)变换后的向量为：F(i)＝(Fi1,Fi2,L,Fip)T∈Rp,(i＝1,2,L,n)，其中，式中R为自然数，表示取值范围，a为A矩阵中的行向量，T为矩阵转置；
3)利用公式：确定A矩阵中的行向量aj，i表示第i个样本，Fij表示第i个样本对应的Fj，使与aj对应的Fj的方差最大，式中，Var(Fj)为煤样品变换后的值与均值的均方差，是样本变换后Fj的均值，V是X的协方差矩阵，V是实对称矩阵，假设V有p个非负特征值λ1,λ2,Lλp，且按降序排列，特征值所对应的特征向量分别为a1,a2,Lap，且各个特征向量之间正交，则公式可变成：
在步骤3中，将不同灰分含量的煤样品吸收光谱数据X作为输入，输入前的X光谱数据维度为200维，利用上述主成分分析方法进行降维和特征提取，输出光谱X1维度为200维,经过计算和验证发现可以使用20维的输出光谱数据X1来代替原来的吸收光谱数据X。
步骤(4)利用煤样品的PCA输出光谱数据X1作为新的输入采用非线性拟合方法进行建模。确立煤炭种类、煤中灰分和其太赫兹吸收光谱X1的关系曲线，构建基于太赫兹光谱的煤炭灰分判别模型，如图4所示；
在步骤4中，非线性拟合方法可以利用matlab拟合工具箱实现.
在步骤4中，可以利用标准差RMSE(Root Mean Square Error)和相关系数R-square来评估拟合结果。其中R-square＝1-SSE/SST，SSE(Sum of Squares for Error)为残差平方和，SST(Sum of Squares for Total)为总平方和。在本例中，利用该方法拟合出的曲线所得到的煤炭灰分预测值与实际值的标准差为RMSE＝0.989，R-square＝0.98，如图5所示。
步骤(5)采集如步骤(1)中不同品种的煤粉，获取煤炭样本的太赫兹时域光谱，经步骤(2)处理后带入步骤(3)和(4)所建模型中，进行煤炭灰分在线判别。
当该方法用于采煤机调高时，通过步骤1-5实时判别当前灰分，当灰分超过设定阈值时，说明当前切割到岩层，需要立刻调整采煤机滚筒高度。
Claims (8)
1.一种基于太赫兹时域光谱的煤炭灰分判别方法，其特征在于包括如下步骤：
a.挑选不同品种的煤，利用挑选出的煤加入纯岩样制作不同灰分煤样品，利用太赫兹时域光谱仪为对空气直接检测得到无样本时的信号作为参考信号，再利用太赫兹时域光谱仪采集所有煤样品的太赫兹时域光谱数据并分类存储起来，无样本时的信号为对空气直接检测得到的信号，所述的太赫兹时域光谱数据有200维，并将所有煤样品按照中国煤炭分类标准分成：无烟煤、烟煤、褐煤几类；
b.利用快速傅里叶变换将所有煤样品的太赫兹时域光谱转换为太赫兹频域光谱；从所有转换后的太赫兹频域光谱中提取吸收光谱，并对吸收光谱数据进行平滑和加窗预处理，选取数据噪声超过1THz以上的煤样品吸收光谱，将选取的煤样品吸收光谱中0.4-1.0THz频段光谱数据作为优选的吸收光谱；
c.将所有煤样品，包括不同灰分煤样品和纯煤样品的0.4-1.0THz频段吸收光谱进行标注分类从而生成不同灰分煤样品和纯煤样品的煤光谱数据库，利用主成分分析/PCA技术对煤光谱数据库中的所有数据进行降维和特征提取，使获得的太赫兹时域光谱数据从200维降到了20维；
d.利用非线性拟合法对所有降维后的0.4-1.0THz频段吸收光谱进行建模，根据所有煤样品的三类信息：煤炭种类、煤中灰分和其太赫兹吸收光谱确立三类信息的关系曲线，从而构建基于太赫兹光谱的煤炭灰分模型；
e.根据实际需要采集被测区域的煤样品后利用太赫兹时域光谱仪获取煤样品的太赫兹时域光谱，重复步骤b、c、d，将被测区域煤样品的光谱数据带入太赫兹光谱的煤炭灰分模型从而对煤样品中的灰分进行判别。
2.根据权利要求1所述基于太赫兹时域光谱的煤炭灰分判别方法，其特征在于：太赫兹时域光谱采集装置为天津大学组建的系统设备。
3.根据权利要求1所述基于太赫兹时域光谱的煤炭灰分判别方法，其特征在于：每次采集煤样品的太赫兹时域光谱时选取3-5个不同的点位进行采集，每个点位重复采集3-4次光谱，所述煤样品的为粉末状或块状。
4.根据权利要求1所述基于太赫兹时域光谱的煤炭灰分判别方法，其特征在于：在检测过程中，采用太赫兹时域光谱系统的透射扫描模块获取以干燥空气为背景的参考信号时域光谱和不同灰分含量煤样品的太赫兹时域光谱，干燥空气的空气湿度低于5％。
5.根据权利要求1所述基于太赫兹时域光谱的煤炭灰分判别方法，其特征在于：所述利用快速傅里叶变换将所有煤样品的太赫兹时域光谱转换为太赫兹频域光谱，根据太赫兹频域光谱提取折射率光谱n(ω)和吸收光谱α(ω)采用THz光学参数提取模型：
T(ω)＝Esam(ω)/Eref(ω)
其中，T(ω)为煤样品的透射系数光谱，Eref(ω)为参考信号的太赫兹频谱，Esam(ω)为煤样品太赫兹频谱，ω为角频率，d为煤样品的厚度，c为光速，n(ω)为煤样品的折射率，α(ω)为煤样品的吸收系数，ρ(ω)项为煤样品信号与参考信号幅度的比值，项为测试样品信号与参考信号的相位差，k(ω)为测试样品复折射率的虚部，也称为消光系数；测试样品信号为不同灰分的煤炭信号，直接测出来的是时域信号，通过傅里叶变换到频域，然后利用上述公式，可以得到煤的吸收光谱；经傅里叶变换后的太赫兹频谱范围为0.1～3.5THz。
6.根据权利要求1所述基于太赫兹时域光谱的煤炭灰分判别方法，其特征在于所述主成分分析方法为：
1)设不同灰分煤样品的吸收光谱为矩阵X，矩阵X包括n个不同灰分含量的煤样品和p个吸收光谱数据，其中n个不同灰分含量的煤样品为列向量，p个吸收光谱数据向量为行向量，表示同一样本的光谱数据维度：
式中，ei＝(xi1,xi2,L,xip)T∈Rp,(i＝1,2,L,n)为第i个样本，xj＝(x1j,x2j,L,xpj)T∈Rp,(j＝1,2,L,p)为第j个变量；原始矩阵X可以由e向量矩阵和x向量矩阵内积表示；代表原始矩阵X由e向量矩阵和x向量矩阵内积表示；
2)将矩阵X变换改写成：
F＝AXT
改写后的F矩阵的列向量表示不同灰分含量煤样品的个数，行向量表示同一煤样品的光谱数据维度；将p个吸收光谱数据向量x1,x2,…,xp变成新向量F1,F2,…,Fp，将矩阵X以2个矩阵向量内积的形式表示，即：第i个样本ei＝(xi1,xi2,L,xip)T∈Rp,(i＝1,2,L,n)变换后的向量为：F(i)＝(Fi1,Fi2,L,Fip)T∈Rp,(i＝1,2,L,n)，其中，式中R为自然数，表示取值范围，a为A矩阵中的行向量，T为矩阵转置；
3)利用公式：确定A矩阵中的行向量aj，i表示第i个样本，Fij表示第i个样本对应的Fj，使与aj对应的Fj的方差最大，式中，Var(Fj)为煤样品变换后的值与均值的均方差，是样本变换后Fj的均值，V是X的协方差矩阵，V是实对称矩阵，假设V有p个非负特征值λ1,λ2,Lλp，且按降序排列，特征值所对应的特征向量分别为a1,a2,L ap，且各个特征向量之间正交，则公式可变成：
4)设F1＝Xa1为第一主成份，以此类推计算出所有主成份，将主成份按降序排列，因此取前几个对原始数据贡献大的20个主成份，使样本维数降低。
7.根据权利要求1所述基于太赫兹时域光谱的煤炭灰分判别方法，其特征在于：所述加窗预优选后的吸收光谱数据，提取信噪比较高的0.4-1THz频段光谱数据；所述平滑处理利用Savitzky-Golay平滑方法处理光谱曲线，消除吸收光谱中的系统和环境噪声。
8.根据上述任意权利要求所述基于太赫兹时域光谱的煤炭灰分判别方法，其特征在于：将太赫兹时域光谱仪设置在采煤机滚筒上实时采集采煤机滚筒附近物体的灰分，当检测到的灰分若灰分值超过设定阀值时，说明当前采煤机滚筒切割到岩层，需要及时调节采煤机滚筒的高度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810404583.2A CN108279217B (zh) | 2018-04-28 | 2018-04-28 | 一种基于太赫兹时域光谱的煤岩判别方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810404583.2A CN108279217B (zh) | 2018-04-28 | 2018-04-28 | 一种基于太赫兹时域光谱的煤岩判别方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108279217A true CN108279217A (zh) | 2018-07-13 |
| CN108279217B CN108279217B (zh) | 2020-12-22 |
Family
ID=62812066
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810404583.2A Active CN108279217B (zh) | 2018-04-28 | 2018-04-28 | 一种基于太赫兹时域光谱的煤岩判别方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108279217B (zh) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109544602A (zh) * | 2018-11-27 | 2019-03-29 | 济南大学 | 一种基于快速傅里叶变换的燃烧动态监测方案 |
| CN109709063A (zh) * | 2019-02-14 | 2019-05-03 | 重庆市计量质量检测研究院 | 基于太赫兹时域光谱技术的翡翠鉴定方法 |
| CN110080766A (zh) * | 2019-04-30 | 2019-08-02 | 中国矿业大学 | 综采工作面煤岩识别装置及方法 |
| CN110455742A (zh) * | 2019-08-16 | 2019-11-15 | 中央民族大学 | 利用太赫兹时域光谱技术鉴别黄曲霉毒素b1和b2的方法 |
| CN111307745A (zh) * | 2020-04-13 | 2020-06-19 | 中国科学院上海技术物理研究所杭州大江东空间信息技术研究院 | 一种煤粉灰分含量测量方法 |
| CN111912816A (zh) * | 2020-08-11 | 2020-11-10 | 玉林师范学院 | 太赫兹光谱与均值漂移算法结合的转基因作物检测方法 |
| CN111965121A (zh) * | 2020-08-25 | 2020-11-20 | 吉林大学 | 一种基于光谱分析技术的煤炭焦渣的快速分类检测方法 |
| CN113504200A (zh) * | 2021-08-04 | 2021-10-15 | 桂林电子科技大学 | 一种基于太赫兹时域光谱的香豆素类食品添加剂鉴别方法 |
| CN114264627A (zh) * | 2021-12-20 | 2022-04-01 | 河南工业大学 | 太赫兹传感器及其使用方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130087709A1 (en) * | 2011-10-07 | 2013-04-11 | Heidy Hodex Visbal Mendoza | Mercury gas sensing using terahertz time-domain spectroscopy |
| CN105279379A (zh) * | 2015-10-28 | 2016-01-27 | 昆明理工大学 | 基于凸组合核函数主成分分析的太赫兹光谱特征提取方法 |
| CN106525759A (zh) * | 2016-10-11 | 2017-03-22 | 中国农业大学 | 一种基于衰减全反射太赫兹介电谱鉴别蜂蜜品种的方法 |
| CN106645014A (zh) * | 2016-09-23 | 2017-05-10 | 上海理工大学 | 基于太赫兹光谱的物质识别方法 |
| CN107561033A (zh) * | 2017-09-21 | 2018-01-09 | 上海理工大学 | 基于太赫兹光谱的混合物中关键物质定性和定量测定方法 |
-
2018
- 2018-04-28 CN CN201810404583.2A patent/CN108279217B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130087709A1 (en) * | 2011-10-07 | 2013-04-11 | Heidy Hodex Visbal Mendoza | Mercury gas sensing using terahertz time-domain spectroscopy |
| CN105279379A (zh) * | 2015-10-28 | 2016-01-27 | 昆明理工大学 | 基于凸组合核函数主成分分析的太赫兹光谱特征提取方法 |
| CN106645014A (zh) * | 2016-09-23 | 2017-05-10 | 上海理工大学 | 基于太赫兹光谱的物质识别方法 |
| CN106525759A (zh) * | 2016-10-11 | 2017-03-22 | 中国农业大学 | 一种基于衰减全反射太赫兹介电谱鉴别蜂蜜品种的方法 |
| CN107561033A (zh) * | 2017-09-21 | 2018-01-09 | 上海理工大学 | 基于太赫兹光谱的混合物中关键物质定性和定量测定方法 |
Non-Patent Citations (1)
| Title |
|---|
| 王昕 等: "基于太赫兹光谱技术的煤岩识别方法", 《煤矿开采》 * |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109544602B (zh) * | 2018-11-27 | 2023-02-21 | 济南大学 | 一种基于快速傅里叶变换的燃烧动态监测方法 |
| CN109544602A (zh) * | 2018-11-27 | 2019-03-29 | 济南大学 | 一种基于快速傅里叶变换的燃烧动态监测方案 |
| CN109709063A (zh) * | 2019-02-14 | 2019-05-03 | 重庆市计量质量检测研究院 | 基于太赫兹时域光谱技术的翡翠鉴定方法 |
| CN110080766A (zh) * | 2019-04-30 | 2019-08-02 | 中国矿业大学 | 综采工作面煤岩识别装置及方法 |
| CN110080766B (zh) * | 2019-04-30 | 2024-03-01 | 中国矿业大学 | 综采工作面煤岩界面识别装置及方法 |
| CN110455742A (zh) * | 2019-08-16 | 2019-11-15 | 中央民族大学 | 利用太赫兹时域光谱技术鉴别黄曲霉毒素b1和b2的方法 |
| CN110455742B (zh) * | 2019-08-16 | 2022-03-01 | 中央民族大学 | 利用太赫兹时域光谱技术鉴别黄曲霉毒素b1和b2的方法 |
| CN111307745A (zh) * | 2020-04-13 | 2020-06-19 | 中国科学院上海技术物理研究所杭州大江东空间信息技术研究院 | 一种煤粉灰分含量测量方法 |
| CN111912816A (zh) * | 2020-08-11 | 2020-11-10 | 玉林师范学院 | 太赫兹光谱与均值漂移算法结合的转基因作物检测方法 |
| CN111965121A (zh) * | 2020-08-25 | 2020-11-20 | 吉林大学 | 一种基于光谱分析技术的煤炭焦渣的快速分类检测方法 |
| CN111965121B (zh) * | 2020-08-25 | 2023-03-21 | 吉林大学 | 一种基于光谱分析技术的煤炭焦渣的快速分类检测方法 |
| CN113504200A (zh) * | 2021-08-04 | 2021-10-15 | 桂林电子科技大学 | 一种基于太赫兹时域光谱的香豆素类食品添加剂鉴别方法 |
| CN114264627A (zh) * | 2021-12-20 | 2022-04-01 | 河南工业大学 | 太赫兹传感器及其使用方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108279217B (zh) | 2020-12-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108279217A (zh) | 一种基于太赫兹时域光谱的煤岩判别方法 | |
| CN108458989A (zh) | 一种基于太赫兹多参数谱的煤岩识别方法 | |
| CN104897607A (zh) | 便携式近红外光谱食品建模与快速检测一体化方法和系统 | |
| CN102004088B (zh) | 一种基于神经网络的煤质特性在线测量方法 | |
| CN105008898A (zh) | 用于组合的libs和ir吸收光谱法研究的系统和方法 | |
| CN101210873A (zh) | 一种利用太赫兹时域光谱快速检测植物油纯度的方法及设备 | |
| CN102564989A (zh) | 一种基于太赫兹光谱的煤炭无损快速检测方法 | |
| CN107703097A (zh) | 利用衰减全反射探头与近红外光谱仪构建快速预测原油性质的模型的方法及其应用 | |
| CN107860743A (zh) | 利用反射式近红外光纤探头构建快速预测原油性质的模型的方法及其应用 | |
| CN105699319A (zh) | 一种基于高斯过程的近红外光谱煤全水分快速检测方法 | |
| CN105044024A (zh) | 一种基于近红外光谱技术对葡萄果实进行无损检测的方法 | |
| CN105092436A (zh) | 一种沉积物粒度光谱分析方法及装置 | |
| CN104596975A (zh) | 近红外漫反射光谱技术测定造纸法再造烟叶木质素的方法 | |
| Yu et al. | A novel method of on-line coal-rock interface characterization using THz-TDs | |
| CN105486661A (zh) | 一种基于高斯过程的近红外光谱煤灰分快速检测方法 | |
| Wang et al. | Wood species identification using terahertz time-domain spectroscopy | |
| Begum et al. | Reflectance spectroscopy based rapid determination of coal quality parameters | |
| AU2008328624B2 (en) | Automatic waste treatment process | |
| Hadiwijaya et al. | Water content prediction of ‘crystal’guava using visible-near infrared spectroscopy and chemometrics approach | |
| Olale et al. | Limitations to use of infrared spectroscopy for rapid determination of carbon-nitrogen and wood density for tropical species | |
| CN105067556A (zh) | 一种快速检测生物液态样品的光学装置及光学方法 | |
| Liu et al. | Classification of sand grains by terahertz time-domain spectroscopy and chemometrics | |
| Liu et al. | Automatic Coal-Rock Recognition by Laser-Induced Breakdown Spectroscopy Combined with an Artificial Neural Network | |
| Rahman et al. | A comprehensive review on wood chip moisture content assessment and prediction | |
| Duan et al. | Prediction of Wood Moisture Content Based on THz Time-Domain Spectroscopy. |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |