典型煤的近红外反射光谱特征的方向性研究

杨 恩，王世博

1.江苏建筑职业技术学院智能制造学院，江苏 徐州 221116 2.中国矿业大学机电工程学院，江苏 徐州 221116

引 言

煤炭是我国经济和社会发展的主要能源和重要原材料，利用先进的科技手段高效、绿色、安全地探测和开发煤炭具有重要的意义。高光谱遥感作为一种可见-近红外波段光谱分辨率较高、地物探测更为精确的遥感技术，近年来已在煤矿区的煤炭资源调查、环境监测等领域得到了一定的研究和应用[1-3]。由于地物的反射光谱特征具有方向性，即地物的反射光谱特征随探测器探测角度的变化而变化，因此获得地物的方向反射光谱特征是实现地物高光谱遥感精确识别结果的重要基础[4-7]。获得煤、矸石、植被、水体等地物的方向反射光谱特征对于煤矿区高光谱遥感更为精确的探测结果具有重要作用。

目前，为校准机载和卫星传感器对地物的高光谱遥感数据，同时为提供地面尺度的地物反射光谱数据，已有多位学者测定了多种地物的方向反射光谱特征，并采用了不同的测量装置。Knowles-Middleton等[4]最早采用了一个封闭的半圆形盒子，里面设有链传动轨道，通过光电倍增管输出电压测定了雪面定向反射特征。Schopfer等[8]采用双视角测角支架和两台ASD FieldSpec-3光纤光谱仪测定了人造地物的方向反射光谱特征。Cloutis等[5]采用ASD FieldSpec Pro HR光纤光谱仪、准直光源、测角仪等装置测定了碳质球粒陨石的方向反射光谱特征。赵子傑等[6]采用ASD FieldSpec-3光纤光谱仪、NENULGS角度测量系统等测定了沙地表面的方向反射光谱特征。路鹏等[7]为研究月球表面主要矿物的方向反射特征，采用ASD FieldSpec-3 Hi-Res光谱仪测定了斜长石和单斜辉石的方向反射光谱特征。以上文献报道了对于岩石、矿物、植被、水体、人造物等地物的方向反射光谱特征的研究，然而煤作为露天煤矿等煤矿区的主要地物，鲜有文献报道相关煤的方向反射光谱特征。研究典型煤的方向反射光谱特征将为矿区煤炭高光谱遥感最优探测几何的选择、煤炭资源精准探测等提供参考。

选取了我国无烟煤、烟煤、褐煤三大类煤中的无烟煤一号、贫煤、气煤、褐煤二号4种典型煤样，利用方向反射测量球坐标实验装置测定了每种煤半球空间各个反射方向的近红外波段(1 000～2 500 nm)反射光谱曲线，获得了每种煤特征波长处反射率的空间分布特征，分析了每种煤前向反射(180°探测方位角)和后向反射(0°探测方位角)特征波长反射率的变化规律。

1 实验部分

1.1 煤样处理

从宁夏汝其沟煤矿、山东东风煤矿、山东兴隆庄煤矿、吉林舒兰煤矿分别收集了无烟煤一号、贫煤、气煤、褐煤二号4种典型煤样各1块。将每种煤样切割成10 cm×10 cm×5 cm的方体试样，并对其中一侧10 cm×10 cm表面打磨平整，采集此表面反射光谱后，将此表面煤样切割下来碾成粉末，取其中一部分粉末依次进行工业分析、煤灰X射线荧光(XRF)分析，其中XRF分析所用仪器为德国Bruker Spectrometer S8 Tiger。按煤阶由高到低顺序排列，4种煤的空气干燥基固定碳、挥发分、灰分、水分工业分析结果及煤灰XRF分析结果如表1所示。

表1 4种典型煤样的主要成分Table 1 Main compositions of the four typical coal samples

1.2 煤近红外方向反射光谱测定

如图1所示，搭建方向反射测量球坐标实验装置，该装置主要由球坐标支架、100 W卤钨灯光源、准直镜、Y型石英光纤、650 nm激光器、PTFE参照白板、荷兰Avantes AvaSpec近红外光纤光谱仪、采集计算机组成，其中近红外光纤光谱仪波长范围为1 000～2 500 nm，Y型光纤合并端连接准直镜，一个分支端连接光谱仪用于光谱数据采集，另一个分支端连接激光器用于指示探测位置，准直镜固定于球坐标支架反射角滑轨用于收集反射光。

图1 煤方向反射实验台及光谱采集示意图(a)：方向反射测量球坐标实验装置；(b)：煤方向反射光谱采集示意图Fig.1 Directional reflection platform and spectra acquisition diagram for coal(a)：Spherical coordinate device for directional reflection measurment; (b):Directional reflectance spectra acquisition diagram for coal

对于每种煤块样，方向反射测量球坐标实验装置采用卤钨聚光灯模拟太阳照射试样被测表面，卤钨灯光束入射方向固定，其方位角为0°，入射角为45°，准直镜探头模拟试样反射光谱探测接收器，其接收方位角(φ)设定0°，45°，90°，135°和180° 5个位置，每个方位角下设定0°，15°，30°，45°，60°和75° 6个接收反射角(θo)，光源照射距离和准直镜探测距离均为0.5 m，即对于每种试样共采集30个探测方位的反射光谱数据。在任一探测方向，测定试样被测表面内不同区域反射光谱数据，取其均值作为试样在该探测方向的反射光谱数据。按对称取试样在180°～360°方位角相应探测方位的光谱数据[6-7]，以此获得每种煤样近似半球空间各反射方向的近红外波段反射光谱。

2 结果与讨论

2.1 煤近红外方向反射光谱曲线特征

实验所得无烟煤一号、贫煤、气煤、褐煤二号4种典型煤样表面各探测方位角(φ)下各反射角(θo)方向的近红外波段(1 000～2 500 nm)光谱反射率曲线分别如图2、图3、图4、图5所示(两端边缘波段1 000～1 049和2 450～2 500 nm光谱曲线受光谱仪暗电流影响，噪声较大，图中已去除)。对于同一种煤，不同探测方向的近红外光谱曲线表现出了基本相似的波形，但整体反射率大小和局部波形特征有一定的差异性。受反射能量越多，接收信噪比越高的影响，所得反射光谱曲线具有随整体反射率升高，光谱吸收谷特征更加明显的规律。同时，4种煤样按煤阶由高到低顺序，各探测方向近红外反射光谱波形整体上符合吸收谷特征增多、反射率增大的趋势[9-10]。

图3 贫煤各位角(φ)下的近红外反射光谱曲线(a)：φ=0°；(b)：φ=45°；(c)：φ=90°；(d)：φ=135°；(e)：φ=180°Fig.3 Near-infrared spectral reflectance curves of meager coal under different azimuth angles (φ)(a)：φ=0°；(b)：φ=45°；(c)：φ=90°；(d)：φ=135°；(e)：φ=180°

如图2(e)、图3(e)、图4(e)、图5(e)所示，在180°探测方位角，即前向反射方向时，4种煤样均表现出随探测反射角变化，反射光谱曲线变化较为明显、强度差异较大的特征，整体规律为随反射角增大，反射光谱曲线上升。而在0°探测方位角，即后向反射方向，如图2(a)、图3(a)、图4(a)、图5(a)所示，4种煤样均表现出随探测反射角变化，反射光谱曲线变化相对不明显、强度差异较小的特征，无烟煤一号后向反射光谱曲线整体最大反射率出现在15°反射角附近，贫煤、气煤、褐煤二号后向反射光谱曲线整体最大反射率均出现在45°～60°反射角附近。

如图2(b—d)、图3(b—d)、图4(b—d)、图5(b—d)所示，45°，90°，135° 3个方位角下，4种煤样反射光谱曲线强度无明显的随反射角变化规律。

图4 气煤各方位角(φ)下的近红外反射光谱曲线(a)：φ=0°；(b)：φ=45°；(c)：φ=90°；(d)：φ=135°；(e)：φ=180°Fig.4 Near-infrared spectral reflectance curves of gas coal under different azimuth angles (φ)(a)：φ=0°；(b)：φ=45°；(c)：φ=90°；(d)：φ=135°；(e)：φ=180°

图5 褐煤二号各方位角(φ)下的近红外反射光谱曲线(a)：φ=0°；(b)：φ=45°；(c)：φ=90°；(d)：φ=135°；(e)：φ=180°Fig.5 Near-infrared spectral reflectance curves of No.2 lignite under different azimuth angles (φ)(a)：φ=0°；(b)：φ=45°；(c)：φ=90°；(d)：φ=135°；(e)：φ=180°

2.2 煤近红外特征波长反射率空间分布特征

煤在近红外波段的反射光谱吸收谷特征主要受其所含的水分、矿物、有机成分的影响。煤中的水分在1 400和1 900 nm波长附近有吸收谷特征[10-13]。煤中的主要矿物成分最为常见的有黏土矿物、石英，黏土矿物中—OH、石英结合H2O也对煤1 400和1 900 nm波长附近的吸收谷特征产生影响，同时黏土矿物中的Al—OH对煤2 200 nm波长附近的吸收谷特征产生主要影响[10-12]。煤有机成分中脂肪侧链CH3、CH2、芳香结构C—H等弯曲、伸缩振动的合频、倍频，所产生的近红外波段反射光谱吸收谷特征主要出现在1 700和2 300 nm波长附近[11,14]。基于以上主要成分对煤近红外波段反射光谱吸收谷特征的影响，选取每种煤样在半球空间各方向近红外反射光谱的1 400，1 700，1 900，2 200和2 300 nm 5个特征波长处的反射率进行分析，得到无烟煤一号、贫煤、气煤、褐煤二号特征波长处的反射率空间分布极坐标云图，如图6、图7、图8、图9所示。

图6、图7、图8、图9中各特征波长处反射率空间分布极坐标云图表明4种典型煤在近红外波段各特征波长处均表现出了一定的双向反射特征。180°探测方位角前向反射时，在各特征波长处4种典型煤均表现出了较明显的反射热点特征，即在某个方向反射能量较大，在反射率空间分布极坐标云图中，这个方向的颜色相对于周围其他方向的颜色更趋向于暖色(红色)，故此，称其为热点。整体上反射角越大，反射热点特征越明显。而0°探测方位角后向反射时，在各特征波长处4种典型煤均表现出了相对较弱的反射热点特征，整体表现为无烟煤一号后向反射热点相对于贫煤、气煤、褐煤二号更为明显，贫煤、气煤、褐煤二号按煤阶降低的顺序，后向反射热点特征整体上逐渐变明显。

如图6、图7、图8所示，无烟煤一号、贫煤、气煤3种煤中，同一种煤在各特征波长间表现出了相似的反射率分布特征，其中无烟煤一号在后向反射方向15°反射角附近出现局部热点；贫煤在后向反射方向45°反射角附近出现局部热点，前向反射方向40°反射角附近出现局部热点；气煤在前向反射方向50°反射角附近出现局部热点。如图9所示，褐煤二号在1 400，1 700和1 900 nm 3个特征波长处表现出相对较明显的后向反射特征，后向反射方向45°～60°反射角附近出现局部热点，2 200和2 300 nm 2个特征波长处的后向反射特征相对不明显。此外，褐煤二号还在1 400，1 700和1 900 nm 3个特征波长处，90°和270°方位角的45°反射角方向附近出现局部热点。

图6 无烟煤一号特征波长(λ)反射率极坐标云图(a)：λ=1 400 nm；(b)：λ=1 700 nm；(c)：λ=1 900 nm；(d)：λ=2 200 nm；(e)：λ=2 300 nmFig.6 Polar nephograms of reflectance at characteristic wavelengths (λ)of No.1 anthracite(a)：λ=1 400 nm；(b)：λ=1 700 nm；(c)：λ=1 900 nm；(d)：λ=2 200 nm；(e)：λ=2 300 nm

图7 贫煤特征波长(λ)反射率极坐标云图(a)：λ=1 400 nm；(b)：λ=1 700 nm；(c)：λ=1 900 nm；(d)：λ=2 200 nm；(e)：λ=2 300 nmFig.7 Polar nephograms of reflectance at characteristic wavelengths (λ)of meager coal(a)：λ=1 400 nm；(b)：λ=1 700 nm；(c)：λ=1 900 nm；(d)：λ=2 200 nm；(e)：λ=2 300 nm

图8 气煤特征波长(λ)反射率极坐标云图(a)：λ=1 400 nm；(b)：λ=1 700 nm；(c)：λ=1 900 nm；(d)：λ=2 200 nm；(e)：λ=2 300 nmFig.8 Polar nephograms of reflectance at characteristic wavelengths (λ)of gas coal(a)：λ=1 400 nm；(b)：λ=1 700 nm；(c)：λ=1 900 nm；(d)：λ=2 200 nm；(e)：λ=2 300 nm

图9 褐煤二号特征波长(λ)反射率极坐标云图(a)：λ=1 400 nm；(b)：λ=1 700 nm；(c)：λ=1 900 nm；(d)：λ=2 200 nm；(e)：λ=2 300 nmFig.9 Polar nephograms of reflectance at characteristic wavelengths (λ)of No.2 lignite(a)：λ=1 400 nm；(b)：λ=1 700 nm；(c)：λ=1 900 nm；(d)：λ=2 200 nm；(e)：λ=2 300 nm

根据前述4种典型煤样近红外反射光谱曲线随整体反射率升高，光谱吸收谷特征更加明显的规律，优先选择煤特征波长反射率空间分布极坐标云图局部热点所处的探测方向作为高光谱遥感传感器对矿区煤炭的最佳探测方向。

2.3 煤近红外特征波长反射率-反射角相关性

由煤近红外特征波长反射率空间分布特征可知，在180°探测方位角前向反射方向和0°探测方位角后向反射方向，随反射角变化煤有较复杂的反射强度分布，然而这两个探测方向同时又是地物反射光谱探测中主要使用的探测方向[5-7]。据此，进一步研究4种典型煤前、后向反射光谱探测中特征波长处反射率与反射角的关系。在前向反射方向和后向反射方向，对每种煤在1 400，1 700，1 900，2 200和2 300 nm 5个特征波长处进行反射率与反射角的相关性分析，如图10、图11、图12、图13所示。

从图10、图11、图12、图13中4种煤特征波长反射率-反射角分布点可知，在前向反射时，4种煤均呈现随反射角增大，特征波长反射率整体上增大的趋势；在后向反射时，4种煤均呈现随反射角增大，特征波长反射率整体上先增大后减小的趋势。在前向反射时，如图10(a)、图11(a)、图12(a)、图13(a)所示，采用线性函数拟合每种煤特征波长反射率与反射角关系；在后向反射时，如图10(b)、图11(b)、图12(b)、图13(b)所示，采用高斯函数拟合每种煤特征波长反射率与反射角关系。4种煤前向反射和后向反射所用两种拟合方法的拟合方程决定系数均大于0.75。

图10 无烟煤一号特征波长反射率-反射角(θo)回归曲线(a)：前向反射；(b)：后向反射Fig.10 Regression curves of reflectance-reflection angle at characteristic wavelengths of No.1 anthracite(a)：Forward direction;(b)：Backward direction

图11 贫煤特征波长反射率-反射角(θo)回归曲线(a)：前向反射；(b)：后向反射Fig.11 Regression curves of reflectance-reflection angle at characteristic wavelengths of meager coal(a)：Forward direction;(b)：Backward direction

图12 气煤特征波长反射率-反射角(θo)回归曲线(a)：前向反射；(b)：后向反射Fig.12 Regression curves of reflectance-reflection angle at characteristic wavelengths of gas coal(a)：Forward direction;(b)：Backward direction

图13 褐煤二号特征波长反射率-反射角(θo)回归曲线(a)：前向反射；(b)：后向反射Fig.13 Regression curves of reflectance-reflection angle at characteristic wavelengths of No.2 lignite(a)：Forward direction;(b)：Backward direction

此外，4种煤后向反射高斯函数拟合曲线波峰所对应反射角与特征波长反射率空间分布云图中后向反射热点所在反射角范围基本一致，并且高斯函数曲线波峰呈现出随煤阶降低向较大角度移动的规律。

3 结 论

针对露天煤矿、堆煤场等煤矿区地表的高光谱遥感精准煤炭资源调查、环境监测等问题，获得了典型煤的近红外方向反射光谱特征，所得结论如下：

(1)无烟煤一号、贫煤、气煤、褐煤二号4种典型煤，每种煤各方向近红外波段反射光谱曲线波形相似，但整体反射率大小和局部波形特征具有差异性，表现为反射率越大，吸收谷越明显。随反射角增大，4种煤在前向反射方向各波段反射率均增大，在后向反射方向各波段反射率变化较小。

(2)在近红外波段1 400，1 700，1 900，2 200和2 300 nm 5个特征波长处反射率空间分布中，无烟煤一号、贫煤、气煤、褐煤二号均具有双向反射特征，在前向反射方向，均表现出明显的反射热点；在后向反射方向，反射热点均较弱，其中无烟煤一号相对最明显，贫煤、气煤、褐煤二号按煤阶降低顺序，逐渐变明显。

(3)无烟煤一号、贫煤、气煤、褐煤二号在近红外波段1 400，1 700，1 900，2 200和2 300 nm 5个特征波长处的方向反射率，前向反射时，均随反射角增大而增大，近似为反射角的线性函数；后向反射时，均随反射角增大而先增大后减小，近似为反射角的高斯函数。