第一作者: 徐 奥(1992-),女,硕士研究生,主要从事环境遥感方面的研究。Email: xuao_2010@163.com。
目前我国累计堆存的铁尾矿量高达约50亿t,所产生的尾矿粉尘污染极其严重。为此,使用辽宁鞍山矿区铁尾矿粉尘进行实景观测实验,采用人工模拟降尘与光谱测量手段,研究了降尘量对植物叶片光谱的影响规律; 并利用降尘量与植物叶片光谱相关性最好的优势波段与铁元素的光谱吸收特征,分别建立了叶面降尘量的基于植物叶片的优势波段反射率和基于铁元素的独特光谱吸收指数的2种定量反演模型。研究表明: 随着铁尾矿粉尘降尘量逐渐增大,植物叶片反射光谱曲线与粉尘的光谱曲线差异逐渐减小; 2种反演方式中降尘量与植物叶片光谱变量均呈极显著相关,但基于铁元素光谱吸收特征的定量反演模型精度更高。研究结果可为高光谱遥感应用于矿区降尘量定量监测提供基础模型与技术依据。
In China, iron tailings dumps have been accumulated up to about 5 billion tons. The tailings have led to extremely serious dust pollution. Therefore, dust effects on leaf spectra were studied on the basis of the observation of real experiments with Anshan mine tailings dust and by means of artificial simulated dust and spectral measurements. The dust samples of iron tailings were taken from the Anshan mining area. The quantitative inversion of foliar dustfall was realized by using the band of the best correlation between the dustfall and the vegetation leaf spectrum and the characteristics of absorption spectra of iron respectively. The results show that, when the dustfall of iron tailings on leaf increased, the differences of spectral curve between leaf and dust decreased. In the two inversion methods, dustfall and vegetation leaf spectral variables were significantly related to each other. Furthermore, the precision of the inversion modeling according to spectral characteristics of iron is higher than that of the one according to best correlation band. The results could provide basic model and technical basis for quantifying the amount of mining dust monitoring with hyperspectral remote sensing.
粉尘是影响我国空气质量的主要污染物之一[1]。在一些依靠固体矿床开采发展起来的矿区城市, 粉尘污染现象普遍比较严重, 尾矿库表面产生的粉尘可传播到很远距离[2], 会对周边的生态环境造成不良影响, 导致土壤污染、土地退化、植被破坏, 甚至直接威胁到人畜生存安全[3]。降尘量是衡量一个区域生态环境质量的重要因子之一, 对其进行快速、准确监测可为地区生态环境质量评估提供基础数据和评价依据。传统的降尘量测定工序过于复杂、工作量大、效率低, 且得到的数据不够准确[4, 5]。遥感方法具有探测范围大、客观性强, 并能迅速反映动态变化的优势, 比一般探测手段的可靠性更强[6]。目前, 已有学者针对降尘量遥感监测进行了研究。Baddock等[7]以澳大利亚艾尔湖的4次典型沙尘暴为研究对象, 评价了运用MODIS数据进行粉尘源观测的可行性; Gangale等[8]研究了大气红外探测器获取的光谱数据, 验证了火山灰在红外区(800~1 200 nm)存在的特殊波段, 推断出火山灰粒度分布、红外透明度和矿物成分含量; Markowicz等[9]根据波兰观测站的火山灰遥感观测数据, 研究了火山灰光学厚度的光谱分布特征, 并与雷达观测数据进行了比较。
植物在陆表广泛覆盖, 既是卫星遥感观测和记录的第一表层, 也是遥感图像表征的直接信息[10]。因此, 以植物为背景对象研究降尘量的遥感监测方法具有普遍意义。光谱特征是利用遥感技术对相关目标进行定性与定量探测的重要依据。目前, 已有学者针对降尘覆盖下的植物叶片的光谱特征进行研究。谭德军等[11]对重庆市万盛矿区粉尘污染的典型农作物进行了高光谱反射率测定, 分析了粉尘污染导致的植物总体光谱反射率和在1 000 nm及1 200 nm附近的光谱小吸收谷异常, 揭示了受粉尘污染的最严重区域和植物种类; 王涛等[12]对叶面除尘前、后的反射光谱和一阶导数光谱的差异进行了分析比较, 表明叶面降尘可提高可见光波段的反射率、降低近红外波段的反射率, 且对可见光波段的影响要大于近红外波段; 万余庆等[13]对长寿草和芋头叶进行了煤粉尘覆盖条件下的光谱测试研究, 表明煤粉尘量越大、植物自身的光谱特征越弱、粉尘的光谱特性越强。但是, 专门针对铁尾矿降尘的研究还很少。我国铁矿资源丰富, 但铁矿资源品位低, 铁尾矿排放量大, 每生产1 t铁精矿平均排出2.5 t铁尾矿[14]。据统计, 在1949―2009年间, 全国铁尾矿排放量约为62亿t, 特别是2003年以后尾矿排放量呈大幅度上升趋势, 2009年排出的尾矿量达5.5亿t, 2010年排放量超过6亿t[3]。当前, 对排出的大部分铁尾矿均采取了堆存处置, 导致铁尾矿的堆存量巨大。截至2010年, 我国累计堆存的铁尾矿量高达50亿t左右[15], 且随着铁矿开采和铁尾矿排放的不断继续, 铁尾矿堆存量还将日益增大[16]。辽宁鞍山地区是我国最重要的钢铁工业基地之一[17], 拥有丰富的铁矿资源[18]。但是, 鞍山市东部和南部的铁矿尾矿库在主导风向作用下, 常使周边地区降尘量远超过8 t/(月· km2)的标准, 污染十分严重[17]。本文以鞍山齐大山铁矿尾矿砂为粉尘样本, 通过实地监测实验, 采用不同级差控制植物叶面铁尾矿粉尘降尘量, 测量不同降尘量情况下的叶片光谱; 进而统计分析降尘量反演的优势波段, 建立植物粉尘降尘量反演模型; 同时, 基于铁尾矿粉尘中铁元素独特的光谱反射特性建立降尘量反演模型, 并对2种反演模型的精度进行对比。本文旨在通过实地观测实验为高光谱遥感应用于矿区降尘量定量监测提供基础模型与技术依据。
本文实验采集的植物样品为绿萝叶片。摘取绿萝叶片后, 立即将其固定于实验测试的水平平台上, 以确保叶面平整。粉尘样本采集自鞍钢矿业公司齐大山铁矿尾矿库, 主要成分及其百分比含量如下: SiO2, 82.28%; TFe(全铁), 9.90%; FeO, 1.62%; MgO, 0.85%; Al2O3, 0.73%; CaO, 0.66%。粉尘的粒度特征为: 粒径≤ 19.348 μ m占5.31%, (19.348, 283.704) μ m占89.38%, ≥ 283.704 μ m占5.31%。
根据前期实验观测分析, 降尘量达240 g/m2(未超过实际矿区周边降尘量)时, 植物叶片光谱特征基本表现为粉尘光谱特征。辽宁省环保厅规定降尘量标准为8 t/(月· km2)[17], 本次实验设置的粉尘浓度级差为辽宁省标准的一半, 即4 g/m2, 共62组(0~244 g/ m2)。实验时, 将绿萝叶片放置于其外接矩形大小的盒内, 按4 g/m2级差逐级将铁尾矿粉尘均匀覆盖到矩形区域上。图1示出不同降尘量下的植物叶片。
本实验使用 SVC HR-1024光谱仪(350~2 500 nm)采集绿萝叶片的反射光谱信息。仪器参数见表1。
![]() | 表1 SVC HR-1024基本参数 Tab.1 Basic parameters of SVC HR-1024 |
实验光源为卤素灯, 高度角为60° , 距目标50 cm; 光谱仪测量时间为2 s; 使用4° 镜头, 镜头距叶片高度设为55 cm, 以确保叶片充满视场(图2)。
每添加1次粉尘覆盖, 测量1次光谱曲线, 并保持测量位置和角度不变, 共测得到62条光谱曲线。
用Excel的随机数发生器将62条实测光谱曲线数据随机分成2组, 每组各31条, 第1组用于建模, 第2组用于验证。首先, 用第1组光谱测试曲线进行光谱-叶面降尘量相关性分析, 并选择相关系数绝对值最大的波段作为优势波段, 对降尘量进行线性回归分析, 建立降尘量的统计反演模型; 同时, 根据铁尾矿砂的吸收特征, 选择900 nm附近的光谱吸收指数(spectral absorption index, SAI)作为光谱变量, 建立降尘量的特征反演模型; 然后, 利用第2组数据对上述反演模型进行验证, 评估2种模型的预测能力和稳定性。
本文选择均方根误差(root mean square error, RMSE)、平均相对误差(mean relative error, MRE)和决定系数(R2)作为精度指标, 即
式中: N为样本数目; xi为反演得到的降尘量; yi为实测的降尘量。RMSE越小, 表明拟合结果越好; MRE越小, 表示分析结果越接近真实值。R2表示回归分析趋势线的预测值与对应的实测数据之间的拟合程度, 取值范围为0~1; 当趋势线的R2等于或接近1时, 其可靠性最高。
图3示出铁尾矿光谱曲线(粗黑线)、叶片光谱曲线(粗绿线)和不同降尘量下的叶片反射光谱曲线(细彩线, 越偏离粗黑线、降尘量越少)。
从图3可以看出, 在整体上随着铁尾矿粉尘降尘量的逐渐增大, 植物叶片的反射光谱曲线与铁尾矿粉尘光谱曲线的差异逐渐减小; 当降尘量达到设定的最大值时, 植物叶片与粉尘的光谱曲线之间差异最小。具体表现为: 叶片在350~690 nm范围内, 光谱反射率逐渐升高, 绿峰逐渐消失, 光谱曲线趋于平缓; 在708 nm附近出现交点, 红边斜率逐渐减小; 在750~1 250 nm范围内, 光谱反射率逐渐降低, 红外反射坪下降; 在1 384.3 nm附近出现交点; 在1 395.5~1 576.5 nm以及1 763~2 500 nm范围内, 光谱反射率逐渐升高; 在1 580~1 860 nm范围内, 反射峰逐渐降低; 在1 400 nm和1 900 nm附近, 水分吸收谷逐渐消失, 趋于平缓。
为选取合适的反演波段, 对350~2 500 nm范围内所有波段的叶片反射率与铁尾矿粉尘降尘量进行了相关性分析(图4)。
因实测光谱曲线样本数目为31个, 故当|R|> 0.456时, 可满足99%的置信度。如图4所示, 没有被粉色栅栏遮盖的波段均满足99%的置信度。可见, 植物叶片的粉尘降尘量与其光谱反射率在350~711.1 nm(含绿峰), 1 395.5~1 576.5 nm(含水分吸收谷)及1 763~2 500 nm(含水分吸收谷)范围内呈比较显著的正相关; 在711.1~1 395.5nm(含红外反射坪), 1 591~1 738.3 nm范围内呈比较显著的负相关; 而在正、负相关关系转换的过渡波段范围内(如1 576.5~1 591 nm和1 727.6~1 770.1 nm, 植物粉尘降尘量与其光谱反射率的相关性较低, 未能通过显著性检验。
根据植物粉尘降尘量与光谱反射率之间的相关性, 提取出相关系数绝对值最大的波段, 称之为“ 优势波段” 。本文提取的优势波段的波长为2 239.2 nm, 该波段中植物粉尘降尘量与光谱反射率相关系数达到0.982, 呈极显著相关(图4)。据此, 对植物粉尘降尘量与光谱反射率进行建模, 即
式中: Dus为粉尘降尘量, g/m2; Ref为光谱反射率, %。二者的相关系数为0.963, 呈极显著相关。
将第2组曲线相应波段的光谱反射率输入上述反演模型, 计算得到粉尘降尘量; 将其与实测的叶片降尘量进行对比, 得到RMSE=13.288 g/m2, MRE=18.7%, R2=0.960(图5)。
通过非吸收基线方程和比值处理, 剔除非吸收物质的光谱贡献, 王晋年等[19]建立了光谱吸收指数(spectral absorption index, SAI), 并据此测定某一特定波长的相对光谱吸收深度, 给出了综合描述光谱吸收特征的新“ 度量” 。在鞍山铁矿山尾矿中, Fe2O3的含量为14.37%, 仅次于Si
其计算公式为
式中: λ 1, λ 2和λ m分别为2个肩部S1, S2和谷点m对应的波长; ρ 1, ρ 2和ρ m分别为2个肩部S1, S2和谷点m对应的反射率; ρ 为2个肩部连线上的谷点对应的反射率。
通过计算光谱曲线的一阶导数, 提取出900 nm附近的肩部波段759.8 nm和980.3 nm, 谷点波段886.2 nm; 进而计算植物在不同铁尾矿粉尘降尘量下的SA
式中: Dus为粉尘降尘量, g/m2; SA
将第2组光谱反射率数据输入上述模型进行验证, 得到反演的植物粉尘降尘量; 并与实测的粉尘降尘量进行对比, 得到RMSE=12.041 g/m2, MRE=18.2%, R2=0.975(图7)。
2种方法的R2均较高, 说明反演模型具有较高的可靠性。与基于植物粉尘降尘量优势波段反射率的反演方法相比, 基于铁粉尘光谱吸收指数的反演精度相对更高, RMSE降低9.4%, MRE降低0.5%。其原因在于: 优势波段虽然具有普适性, 但缺少对铁尾矿粉尘中特征元素Fe3+的特殊关注; 而运用光谱吸收指数SA
通过实验测定并分析了铁尾矿粉尘在不同降尘量条件下植物叶片的光谱特征变化, 据此建立了铁尾矿降尘量的反演模型, 并进行了验证对比。得出如下结论:
1)植物叶片在铁尾矿粉尘降尘量逐渐增大的情况下, 植物光谱特征(如绿峰、红边、红外反射坪及水分吸收谷)逐渐消失。
2)分别利用2 239.2 nm优势波段反射率和铁光谱吸收指数SA
需要指出的是, 大范围的航空与航天遥感降尘监测是基于植被的冠层尺度, 而本文只是基于地面实验的叶片尺度研究了铁尾矿粉尘降尘量的高光谱反演方法, 与大范围、多尺度的遥感应用还有一定差距。因此, 后续研究将重点关注植物叶片与冠层尺度之间反演结果的转换, 为实现降尘量的大范围航空、航天遥感监测提供技术支持。
The authors have declared that no competing interests exist.
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|
[14] |
|
[15] |
|
[16] |
|
[17] |
|
[18] |
|
[19] |
|
[20] |
|
[21] |
|
[22] |
|