水是维持城市生态环境平衡的重要因素,遥感监测手段已经成为评价城市水环境的重要手段。基于2020年4月珠海一号影像和2020年7月哨兵二号影像,建立水质参数反演模型,对哈尔滨市松北区水体的叶绿素a浓度、悬浮物浓度、透明度和综合营养化指数进行季度监测。研究结果表明:松北区银水湾为污染水体,叶绿素a浓度、悬浮物浓度和综合营养指数值都较高,与其相邻的松花江水质情况较好,水质参数均在正常水体值的范围内。

随着我国城镇化的快速发展,工业生产、农业灌溉和人民日常用水量日趋增加。人们生产、生活产生的工业废水和生活污水被排放到城市河道中,造成城市水体水质日益恶化,严重破坏城市河流生态系统,对城市景观和城市发展造成较大影响,对区域的生态环境、人居环境和人们的身体健康造成了严重的危害。监测和治理水体污染是目前亟待解决的城市环境问题。

传统水质监测方法主要采用地面水站监测,监测范围有限,难以实现长时间连续的大范围监测。而遥感技术可以从全空间领域,多时间序列快速、连续地体现整个水域的水质状况。水体中的叶绿素a浓度、悬浮物浓度和透明度是水质监测的重要指标参数。对此,国内外学者开展了大量的相关研究,王爽等人以Landsat-8/OLI、FY-3A/MERSI和HJ-1B/CCD遥感影像为数据源,结合实测的叶绿素a浓度和水体光谱特征,得到胶州湾的叶绿素a浓度的空间分布情况。宋庆君等人选取太湖为试验区,对比分析TASSAN模型、比值模型、近红外估算模型反演出水体悬浮物浓度,最终得出结论利用近红外(812 nm)波峰高度来估算太湖悬浮物的效果最好。Keiner等人利用Landsat-TM数据,对水体叶绿素、悬浮物和可溶性有机物3个水质参数进行反演。Dall等人利用MODIS影像,构建3个特征波段反演模型,反演水体的叶绿素a浓度。

总结目前的研究成果可知,利用遥感技术进行水质监测的研究已经较为成熟,但是多数研究主要是针对水体的单因子进行反演,如叶绿素a浓度。本文基于哨兵二号(Sentinel-2)卫星影像和珠海一号(ZH-1)卫星影像,对哈尔滨市松北区银水湾等河段的水体的绿素a浓度、悬浮物浓度、透明度和综合富营养化指数进行反演,形成每季度水体监测,并对研究区域内水质进判别和评价,辅助有关部门进行定量污染源解析,服务水体治理措施,可有效地节约成本,提高治理效果。

松北区内河流纵横,水量丰枯不均,主要江河有松花江及其支流金水河和呼兰河,区内河水平缓而宽阔,区内陆河、渠、湖泊较为丰富。本研究以松北区银水湾作为研究区,如图1所示。

银水湾位于哈尔滨松北区,银水湾支渠具有排涝、生态环境、景观、旅游等综合功能,局部放大图如图2所示。支渠总长度为6.494 km,平均水面宽度60 m,水面面积0.37 km2,常年河水呈黄色,污染严重,被称为“金水湾”,现场照片如图3所示。

本文采用4月10日的ZH-1卫星影像和7月24日的Sentinel-2卫星影像,对松北区银水湾,采用多种遥感影像数据源,进行4月和7月的每季度水质监测。

Sentinel-2卫星携带一个多光谱成像仪,可覆盖13个光谱波段,幅宽为290km。13个波段,光谱覆盖范围从可见光到短红外,空间分辨率分别为10m、20 m和60m。两颗互补,重访周期为5d。在光学数据中,Sentinel-2是首颗包含3个“红边”波段的光学卫星,可以提供有关植被状态的关键信息,本文采用2020年7月24日的Sentinel-2影像,对波段2、3、4、8进行合成,空间分辨率为10m。

ZH-1高光谱卫星采用推扫成像方式,单次成像范围为150km×400(km/min),空间分辨率为10m,光谱分辨率为2.5nm,波长范围为400~1 000nm,32个波段。目前,单颗高光谱卫星重访周期为6d,8颗高光谱卫星综合重访周期缩约为1d。本文采用2020年4月10日的ZH-1影像数据,并对数据进行预处理,用于后续水质参数反演工作。

本文主要以高光谱ZH-1卫星影像和多光谱Sentinel-2卫星影像为主要数据源,并对遥感影像进行辐射校正和正射校正等数据预处理。利用自动化提取水体,得到研究区水体分布范围,并对研究区水体的叶绿素a浓度、悬浮物浓度和透明度等水质参数进行反演,在此基础上得到研究区水质的综合污染情况的空间分布图和污染水体的专题报告,具体研究路线 技术路线

为了清晰表示出水体参数的变化特征,需提取出研究区的水体分布情况。本文利用易康(eCognition)解译软件,采用面向对象的方法,利用遥感影像的绿波段和近红外波段构建归一化水体指数(Normalized Difference Water Index,NDWI)对影像水体进行自动化提取,如图5所示。

水体中水质参数的成分和浓度不同,对应遥感影像波段的光谱信息也不同,因此可以根据遥感影像的光谱信息与水体水质参数关系建立模型,得到水体水质参数浓度的空间分布关系。本文利用模型分析法反演出研究区内水体的叶绿素a浓度、悬浮物浓度、透明度以及水体综合营养指数等水质参数,从而得到整个研究区内水质空间分布情况。

水体叶绿素a浓度是评定水体富营养化和水体污染程度的重要指标,水体的叶绿素a浓度不同,对应遥感影像的光谱信息也不同。一般情况下,含有叶绿素a的水体的光谱特性为在近红外波段值最高,在红光波段反射率呈现最低状态。因此,本文采用影像的红光波段和近红外波段,反演出水体叶绿素a浓度,构建的反演模型为:

式中,TLIchi-a为基于叶绿素a浓度的营养状态指数;TLIsd为基于透明度的营养状态指数。

分别对两种营养指数进行归一化操作,然后选取两个营养指数的最大值作为新的综合营养指数,叠加到遥感影像上进行最佳阈值选取,得到污染水体的分布情况。

由上述可知,研究区内水体的叶绿素a浓度主要在30~45mg·m-3之间,银水湾流域的叶绿素a浓度高于相邻松花江流域的叶绿素a浓度,银水湾的平均叶绿素a浓度为42.38mg·m-3,松花江的平均叶绿素a浓度为33.17mg·m-3。

由上述结果可知,水体透明度与水体悬浮物浓度呈负相关,银水湾流域的水体透明度低于松花江的水体透明度,银水湾平均透明度为19.36cm,松花江流域水体透明度为36.47cm。

水体综合营养指数主要综合利用叶绿素a浓度、透明度,分别对其进行归一化处理并取最大值,然后进行最佳阈值选取,分别得到水体的4月和7月每季度监测的污染状况,如图9和图10所示。其中,珠海一号卫星高光谱卫星,光谱信息更为精细化,可将其综合营养指数根据选取最佳阈值,污染程度进行3级显示,4月份珠海一号影像和7月份哨兵二号卫星影像的综合营养指数统计值统计情况,分别见表4和表5。

水体的水质参数对于遥感影像的红外和近红外波段敏感,在此波段内污染水体与非污染水体光谱特征明显不同,可辅助识别水体的污染状况。本文主要是基于现有水质参数反演模型对研究区的水质参数进行反演,缺少实地采集的水质数据,今后可同步获取影像与实地采集数据,更精确和客观地得到水体污染状况分布。且本文主要采用的是10m空间分辨率的遥感影像,对于一些水面较窄的水体很难进行水质情况的判别,尤其城市中的重点排污渠等,在今后的相关工作中应结合更高分辨率的遥感影像。

发表回复

您的电子邮箱地址不会被公开。 必填项已用*标注