德兴铜矿矿山重金属污染修复效果高光谱遥感评价
Evaluating the remediation effect of heavy metal pollution in the Dexing copper mine based on hyperspectral remote sensing
通讯作者: 张登荣(1966-),男,博士,教授,研究方向为地质遥感、GIS开发应用。Email:13805747261@126.com。
责任编辑: 李瑜
收稿日期: 2022-05-16 修回日期: 2022-08-31
基金资助: |
|
Received: 2022-05-16 Revised: 2022-08-31
作者简介 About authors
王嘉芃(1991-),女,硕士,工程师,研究方向为遥感生态监测、环境调查评估。Email:
合理快速评价矿山重金属污染修复的效果对于矿山生态恢复与重建治理工作具有重要意义。以德兴铜矿为例,根据野外实测植被光谱,分析矿区内主要植被的典型光谱特征; 根据实验室化验的植被叶片内重金属含量,分析其重金属含量与光谱特征参数红边位置的关系; 利用 2003年和2009年2景Hyperion高光谱卫星数据计算矿区植被的红边位置,推断矿区植被富集重金属的情况,进而评价矿山重金属污染修复的效果。研究结果表明,在典型复垦区1号、2号尾矿库四周重金属污染修复取得了较好的效果; 与2003年相比,2009年重金属污染修复整体上取得了一定的成效,大部分区域被修复,但仍有部分新增污染区,需采取修复措施。该方法可快速、合理、大范围地评价矿区重金属污染修复的效果。
关键词:
Evaluating the remediation effect of heavy metal pollution in mines properly and rapidly holds considerable significance for ecological restoration and rehabilitation of mines. Based on the field-measured vegetation spectra, this study analyzed the typical spectral features of the main vegetation in the Dexing copper mining area. According to the heavy metal content in the leaves of vegetation tested in the laboratory, this study analyzed the relationship between heavy metal content and red edge position-a spectral characteristic parameter. This study calculated the red edge position of the vegetation in 2003 and 2009 using 2-scene Hyperion hyperspectral data, inferring the heavy metal enrichment in the vegetation of the mining area. Furthermore, this study evaluated the remediation effect of heavy metal pollution in the mining area. The results show that satisfactory results have been achieved from the remediation of heavy metal pollution around mine tailings nos. 1 and 2 in typical reclamation areas. Compared with 2003, 2009 witnessed generally satisfactory remediation effects of heavy metal pollution, with most areas being remedied and some newly polluted areas requiring remediation. The method proposed in this study can achieve a quick and reasonable evaluation of the remediation effect of large-scale heavy metal pollution in mining areas.
Keywords:
本文引用格式
王嘉芃, 徐建国, 沈家晓, 张登荣.
WANG Jiapeng, XU Jianguo, SHEN Jiaxiao, ZHANG Dengrong.
0 引言
德兴铜矿自20世纪80年代初就开始生态修复的试验研究,2014年正式晋升为国家级“绿色矿山”,矿区生态修复取得了很大成效。为了探究利用高光谱遥感技术评价矿山重金属污染修复效果的方法,以德兴铜矿为研究区,根据野外实测植被光谱和实验室化验的植被叶片内重金属含量,分析其重金属含量与光谱特征参数红边位置的关系; 以此为基础,利用Hyperion高光谱卫星数据计算矿区植被的红边位置,推断矿区植被富集重金属的情况,进而评价矿区重金属污染的生态修复效果。
1 研究区概况及数据源
1.1 研究区概况
图1
1.2 数据源及其预处理
1.2.1 野外实测植被光谱数据
野外实测光谱数据为2012年4月15日实测的德兴铜矿矿区的茅草、盐肤木的光谱数据。使用ASD Field Spec便携式地物光谱仪对研究区内的主要植物(茅草、盐肤木)叶片进行光谱采集,测量波段范围350~1 350 nm,光谱分辨率优于8 nm,视场角为25°,实际测量时间为当天9: 00—15: 00。测量植物叶片光谱前,先将光谱仪预热并进行白板校正; 测量时,天空晴朗基本无云,标准光纤探头距离植物叶片1.2 m左右,垂直叶片测量,采样范围约1 m2,每个采样点重复采样10次,以其平均值为该采样点植物叶片的光谱反射率。另外,在远离矿区的无重金属污染区域采集正常植被光谱数据作为背景值。
1.2.2 野外采集植物叶片重金属化验数据
重金属化验数据为2012年4月15日采集的茅草叶片重金属含量数据。野外采集了研究区内典型修复区1号、2号尾矿库及远离矿区的茅草叶片样本,分别用自来水和去离子水冲洗干净,沥干水后在105 ℃杀青30 min,在60 ℃的电热干燥箱中烘干至恒重; 研磨后,过100目尼龙筛网,混匀后准确称取4.0 g样品用干灰化法,盐酸溶解提取; 其中Cu,Zn和Mn的含量采用火焰原子吸收分光光度法测定,Cd和Pb的含量采用石墨炉法测定,得到茅草叶片样本中各重金属含量化验数据(表1)。
表1 不同尾矿库茅草叶片的重金属含量
Tab.1
采样点 | Pb | Cu | Cd | Mn | Zn |
---|---|---|---|---|---|
1号尾矿库 | 4.4 | 19.0 | 0.15 | 46 | 21 |
2号尾矿库 | 3.5 | 12.0 | 0.10 | 100 | 36 |
背景值 | 3.1 | 6.9 | 0.08 | 232 | 25 |
1.2.3 Hyperion高光谱遥感数据
高光谱卫星遥感数据为2003年1月12日和2009年4月18日的2景Hyperion数据。本次研究采用的Hyperion数据为L1R级产品,空间分辨率为30 m,242个波段,光谱区间在400~2 500 nm,光谱分辨率优于10 nm。在进行信息提取之前需对Hyperion高光谱数据进行预处理,主要包括辐射定标、条纹去除、几何纠正及大气校正。辐射定标采用EO-1使用手册指定的定标系数,将影像的DN值转化为绝对辐射亮度值; 检查每个波段数据质量时发现第9—12波段条纹严重,采用低通滤波方法去除条纹,抑制噪声; 将2景Hyperion数据与已校正的Landsat ETM数据进行配准,保证几何纠正的精度在半个像元以内; 利用USGS发布的Hyperion高光谱图像大气纠正算法(FLAASH)进行大气校正。经过预处理后2009年的Hyperion数据如图2(a)所示,图中十字位置的大气校正后植被波谱曲线如图2(b)所示。
图2
2 研究方法
依据前人相关研究[9],选定红边位置作为表征植被光谱特征变异的参量,利用红边位置变化分析矿区植被重金属含量,进而反映矿区重金属污染修复的效果。首先利用矿区植被野外实测的光谱数据分析其光谱特征,计算实测点红边位置,分析红边位置与叶片重金属含量的关系,佐证可以利用红边位置反映矿山重金属污染修复效果; 然后利用2景不同时相的Hyperion高光谱遥感数据计算研究区植被红边位置,分析不同红边位置计算方法的适用性,并确定表征正常情况的红边位置阈值; 最后比较2景影像间红边位置变化,依据重金属含量与红边位置变化关系分别推断2景影像中矿区重金属污染修复的效果。总体技术路线如图3所示。植被红边位置计算方法对结果有较大影响[27]。目前红边位置的计算方法主要分为2类: 一类是基于导数光谱,如最大一阶导数法、拉格朗日法和线性外推法等; 另一类则基于曲线拟合,如四点插值法、倒高斯模型法和多项式拟合法等[27-28]。薛利红等[27]认为最大一阶导数法与叶绿素含量的关系存在不连续性,拉格朗日法和倒高斯模型法在一定程度上存在对高叶绿素含量的饱和现象; 张永贺等[28]的研究表明五次多项式拟合法提取红边位置估算叶绿素含量的精度最高,但是算法过于复杂,四点内插法估算精度次之,但计算简便。Pu等[20]研究利用拉格朗日法、四点插值法、多项式拟合和倒高斯模型4种方法从Hyperion 影像中提取红谷位置和红边位置,结果表明,四点插值法是提取Hyperion 影像中这2个参数的最佳方法。因此,本文采用四点插值法[20,28],选用中心波长值为 671.02 nm,701.55 nm,742.25 nm 以及782.95 nm的4个波段,提取Hyperion数据中植被的红边位置。
图3
具体过程如下: 首先利用归一化植被指数(normalized difference vegetation index,NDVI)NDVI705>0剔除矿区的水体,公式为:
3 结果与分析
3.1 野外实测植被光谱数据特征分析
图4中1号和2号尾矿库内茅草及盐肤木2种植物的地面光谱均出现不同程度的“蓝移”现象,波谱变异特征明显: 叶绿素强反射峰反射增强,红谷吸收减弱; 700 nm附近的反射率上升的斜率减小,红边出现4~19 nm的“蓝移”现象,该特征在图4茅草的表现中尤为显著,这是植物重金属含量过高的典型变异特征。利用四点插值法可知,1号和2号尾矿库茅草的红边位置分别是712 nm和709 nm,背景值的红边位置为721 nm,1号尾矿库茅草红边位置蓝移9 nm,2号尾矿库茅草红边位置蓝移12 nm; 同理,1号和2号尾矿库盐肤木的红边位置分别是716 nm和714 nm,背景值的红边位置为720 nm,1号尾矿库盐肤木红边位置蓝移4 nm,2号尾矿库盐肤木红边位置蓝移6 nm。由表1可知,1号和2号尾矿库中茅草叶片的铜重金属含量远远超过了作为背景值叶片的含量,同时茅草光谱的红边位置明显小于作为背景值正常生长植被光谱的红边位置。结合前人研究[9],植被光谱红边位置变化与植被重金属含量密切相关,一般情况下红边位置越接近正常值,重金属含量越低,反之,红边位置越小,重金属含量越高,这为进一步利用Hyperion高光谱卫星遥感数据提取红边位置,评价研究区重金属污染修复效果奠定了实验基础。
图4
3.2 Hyperion影像时相对红边位置计算结果的影响分析
2003年和2009年Hyperion影像的时相相差3个月,为了分析时相差异对本文结果的影响,选取了研究区内相对均匀分布的2 585个正常植被像元样本数据,利用四点插值法计算正常植被像元的红边位置,分析2幅影像中正常植被的红边位置差异。结果表明,2003年正常植被红边位置均值为722 nm,2009年正常植被红边位置均值为724 nm,相比受重金属污染产生的红边位置蓝移程度而言,2幅影像因时相不同导致的红边位置变化差异不大。这可能与研究区处于亚热带湿润季风区,且以常绿树种为主要植被类型有关。马东辉等[29]指出冬季常绿植被具有相似的光谱特征,叶绿素变化稳定的常绿植被,其光谱随季节变化较为稳定; 吴继友等[30]指出红边特征随季节的长期大趋势变化与叶绿素含量随季节的大趋势变化相关,从春季4月底—6月中旬,红边位置呈微弱上升趋势。由于2003年影像时间为1月12日,2009年影像时间为4月18日,1月初—4月初同属植物生长的萌芽期,叶绿素含量变化较小,因此对研究区内常绿植被的光谱特征影响也较小。
3.3 矿区重金属污染修复效果高光谱遥感评价
图5
图5
2003年和2009年德兴铜矿植被红边位置
Fig.5
Red edge position of vegetation in Dexing copper mine in 2003 and 2009
3.3.1 典型复垦区
这一结果与2个尾矿库开展生态修复与土地复垦工作密不可分。1号尾矿库1986年停止使用,开始进行生态修复与土地复垦研究。早在1999年,张海星等[31]就在尾矿库边坡3 hm2试验地进行生态恢复模式试验,百喜草、无叶节节草和弯叶画眉草等7 种先锋草本植物已在尾砂上定居生长; 2001年,杨修等[32]也在1号尾矿库库面纯尾砂立地进行植被恢复与重建试验,研究了矿区自然植被的次生演替序列和13种试验先锋草种的优劣势,也选出了水蜡烛和假俭草等12种优良的试验先锋草种。截止到2010年,1号尾矿库已复垦面积达到0.775 3 km2,约占总面积的93%; 2号尾矿库仍处于工作状态,已复垦面积为0.579 3 km2,约占2号尾矿库总面积的48%。从典型复垦区来看,矿区重金属污染修复取得了较好的成效。
3.3.2 红边位置变异程度
图6
图6
2003年和2009年相同红边位置区间内像元个数对比
Fig.6
Comparisons of the number of pixels in the same red edge position intervals in 2003 and 2009
3.3.3 红边位置变异区域变化
对2003年和2009年红边位置变异(690~718 nm)区进行遥感解译,运用空间分析功能,获得了2009年相比2003年重金属污染已修复区域和新增污染待修复区(图7)。从图7可看出,矿区重金属污染已修复区多集中于1号尾矿库四周、2号尾矿库四周、西源废石场和铜厂矿区露天采场东南侧,这与德兴铜矿已开展生态修复工作的主要区域相吻合。而新增的重金属污染区则分布于铜厂矿区露天采场北侧山间谷地和富家坞矿区露天采场西南侧山间谷底,这多是由于2个采场内酸性水长期冲刷矿石,导致大量重金属元素溢出,堆积于地势低洼的山间谷地; 且从影像上可以分辨出这些山间谷底植被并不茂盛,对重金属的吸附作用并不明显。这应当引起矿区生态环境部门的警觉,需采取一定的措施对这些区域进行针对性的生态修复,以减轻重金属污染、保护矿区生态环境。
图7
图7
相比2003年,2009年已修复区和新增污染区分布
Fig.7
Compared with 2003, the distribution of repaired and newly polluted areas in 2009
4 结论与讨论
1)通过系统分析了矿区典型植物茅草和盐肤木的野外实测波谱特征与重金属含量的相关关系,以红边位置作为表征光谱特征变异的参量,利用不同时相的Hyperion高光谱卫星数据,获得了矿区植被红边位置信息; 根据矿区植被红边位置的变异程度,推断矿区植被富集重金属的情况,进而评价矿区重金属污染的生态修复效果。该方法具有理论基础,流程简单可行,可移植性强。
2)以德兴铜矿为例,分析讨论了矿区重金属污染在典型复垦区、红边位置变异程度和红边位置变异区域变化上的修复效果: 在典型复垦区1号和2号尾矿库四周重金属污染修复取得了较好的效果; 与2003年相比,2009年重金属污染修复整体上取得了一定的成效,大部分区域被修复,但仍有部分新增污染区。
但是,本次研究的方法仍然存在一定的缺陷: 一是,本次实验选用的Hyperion高光谱卫星遥感数据空间分辨率为30 m,导致计算红边位置时视所有植物为单一物种,忽略了不同物种之间的光谱特征差异性,这可能对红边位置变异阈值的选择有所影响; 二是,植物光谱的红边位置变化除了受矿区重金属污染的影响外,还可能受地表渗漏的酸、碱废水及大气中的氟化氢(HF)、二氧化硫(SO2)等有害气体的影响,虽然本次研究的植物同在一个研究区,但也无法保证除重金属含量外的其他自然条件完全相同,这可能导致在评价重金属污染修复效果时存在一定偏差; 三是,虽然从提取的红边位置信息中可清晰看出红边位置较正常位置是否有蓝移,也能在一定程度上反映该区的重金属含量高低,说明重金属污染修复的效果,但始终无法定量地反映出红边位置蓝移区域的重金属含量。因此,本文评价方法仍需经过更严谨的改进,进一步开展更多针对性研究工作,以求更精确地评价矿区重金属污染修复的效果,服务于矿山生态环境修复与治理工作。
参考文献
Assessment of rehabilitation projects results of a gold mine area using landscape function analysis
[J].DOI:10.1016/j.apgeog.2019.05.005 URL [本文引用: 1]
绿色发展视阈下矿山地质修复模式的升级与重塑
[J].
Upgrading and reshaping of mine geological restoration mode from the perspective of green development
[J].
海南岛2018年矿山地质环境遥感评价和生态修复对策
[J].
Remote sensing evaluation of mine geological environment of Hainan Island in 2018 and ecological restoration countermeasures
[J].
RESTOQUARRY:Indicators for self-evaluation of ecological restoration in open-pit mines
[J].DOI:10.1016/j.ecolind.2019.03.001 URL [本文引用: 1]
Interaction mechanism between floristic quality and environmental factors during ecological restoration in a mine area based on structural equation modeling
[J].DOI:10.1016/j.ecoleng.2018.09.021 URL [本文引用: 1]
基于Hyperion影像植被光谱的土壤重金属含量空间分布反演——以青海省玉树县为例
[J].
DOI:10.13287/j.1001-9332.201606.030
[本文引用: 1]
本研究主要探讨了利用Hyperion影像植被光谱估算土壤重金属含量的可行性.以野外采集的三江源区玉树县48个表层土壤样品As、Pb、Zn、Cd实验室测定含量值,以及从两景Hyperion影像提取的48个土壤样本点相应的176个植被光谱反射率波段及构建的5种植被指数为数据源,利用偏最小二乘回归方法(PLSR)建立土壤各重金属含量与上述两套Hyperion影像上提取的变量之间的估算模型.模型分别为176个植被光谱反射率波段与土壤各重金属含量间的估算模型(植被光谱反射率模型),和以5种植被指数作为自变量,与土壤各重金属含量建立的估算模型(综合植被指数模型).运用验证样本的4种重金属元素实测含量值的标准差与均方根误差的比值(RPD)作为检验标准,As、Pb两种模型RPD均小于1.4,不具备粗略估算能力;Zn、Cd两种模型RPD分别为1.53、1.46与1.46、1.42,均具备粗略估算能力.根据上述结果将Zn的光谱反射率估算模型与Hyperion影像相结合反演得到土壤重金属Zn含量的空间分布,Zn含量在214国道、308省道和乡镇附近偏高,主要受到较强的人类活动影响.表明运用Hyperion高光谱影像植被光谱反射率可以间接估算土壤Zn、Cd元素含量.
The estimating heavy metal concentrations in topsoil from vegetation reflectance spectra of Hyperion images:A case study of Yushu County,Qinghai,China
[J].
航天高光谱遥感应用研究进展(特邀)
[J].
Advances in application of space hyperspectral remote sensing(invited)
[J].DOI:10.3788/IRLA URL [本文引用: 1]
铜胁迫下植被光谱变化规律研究
[J].
Study on spectral change of vegetation under Cu stress
[J].
重金属污染胁迫下盐肤木的生化效应及波谱特征
[J].
Biogeochemistry responses and spectral characteristics of rhus chinensis mill under heavy metal contamination stress
[J].
重金属污染水稻的冠层反射光谱特征研究
[J].
Study on canopy spectral characteristics of paddy polluted by heavy metals
[J].
洎水河流域重金属污染区五节芒叶片光谱特征响应研究
[J].
The spectral response of typical vegetation leaves to heavy metal pollution in Jishui River basin
[J].
Wavelet-based detection of crop zinc stress assessment using hyperspectral reflectance
[J].DOI:10.1016/j.cageo.2010.11.019 URL [本文引用: 1]
黑龙江多金属矿区植物胁迫光谱及其与金属元素含量关系研究
[J].
Vegetation stress spectra and their relations with the contents of metal elements within the plant leaves in metal mines in Heilongjiang
[J].
Assessing leaf spectral properties of Phragmites australis impacted by acid mine drainage
[J].
重金属铜污染植被光谱响应特征研究
[J].
Spectral response of wheat and lettuce to copper pollution
[J].
Estimate of heavy metals in soil and streams using combined geochemistry and field spectroscopy in Wan-sheng mining area,Chongqing,China
[J].DOI:10.1016/j.jag.2014.06.013 URL [本文引用: 1]
Estimation of chlorophyll content in potato leaves based on spectral red edge position
[J].
WREP:A wavelet-based technique for extracting the red edge position from reflectance spectra for estimating leaf and canopy chlorophyll contents of cereal crops
[J].
A modified linear extrapolation approach towards red edge position detection and stress monitoring of wheat crop using hyperspectral data
[J].DOI:10.1080/01431161.2013.877616 URL [本文引用: 1]
Extraction of red edge optical parameters from Hyperion data for estimation of forest leaf area index
[J].
基于红边参数的植被叶绿素含量高光谱估算模型
[J].
Hyperspectral models for estimating vegetation chlorophyll content based on red edge parameter
[J].
基于红边位置变化特征的油菜种植区域提取
[J].
Extracting oilseed rape growing regions based on variation characteristics of red edge position
[J].
Monitoring plant response to phenanthrene using the red edge of canopy hyperspectral reflectance
[J].
利用不同红边位置算法估测玉米叶绿素含量
[J].
Estimation of corn chlorophyll content using different red edge position algorithms
[J].
德兴铜矿矿山污染高光谱遥感直接识别研究
[J].
Identification of mining pollution using Hyperion data at Dexing copper mine in Jiangxi Province,China
[J].
典型矿区土壤重金属污染对植被影响遥感监测分析——以江西省德兴铜矿为例
[J].
Monitoring and analysis of the impacts of soil heavy metal pollution on vegetation in typical mining areas using remote sensing imageries:A case study of Jiangxi Dexing copper mine
[J].
采用不同红边位置提取技术估测蔬菜叶绿素含量的比较研究
[J].
Comparative study on estimation of chlorophyll content in spinach leaves using various red edge position extraction techniques
[J].
基于红边位置的木荷叶片叶绿素含量估测模型研究
[J].
Estimation model of schima superba leaf chlorophyll content based on red edge position
[J].
南京冬季典型植被光谱特征分析
[J].
Research on spectral characteristics of winter typical vegetation in Nanjing
[J].
山东招远金矿区赤松针叶反射光谱红边的季节特征
[J].
Seasonal characteristics of spectral reflectance of korean pine leaves in the gold mine area of Zhaoyuan City in Shandong Province
[J].
德兴铜矿1号尾矿库废弃土地生态恢复试验研究
[J].
Study on the eco-recover test of waste land of 1# tailings bank in Dexing copper mine
[J].
/
〈 |
|
〉 |
