基于遥感的内陆水体水质监测研究进展

　　摘要：从遥感数据、反演方法和水质参数三方面综述了水质遥感监测的研究进展，介绍了国内外常用遥感数据，对比了分析法、经验法、半经验法、机器学习和综合法5种反演方法的优缺点，总结了叶绿素a、悬浮物、有色可溶性有机物等光敏参数和非光敏参数的研究进展。目前在内陆水体水质遥感监测研究中还存在卫星传感器的针对性、反演算法的时空局限性、水质参数光谱特征的复杂性、大气校正的精确性和特殊类型水体的水质监测等问题;指出未来水质遥感监测应围绕新型遥感数据、通用反演模型、不同光谱特征、精确大气校正和特殊水体分类等方面开展。

　　关键词：水质监测;内陆水体;遥感数据;反演方法;水质参数

　　随着全球气候变化和人类活动的影响，许多内陆水体面临水质下降、水体富营养化和水生态系统破坏等一系列问题[1]。水质监测作为水生态环境管理的一个重要环节，可以有效保护水环境、控制水污染和维护水环境的健康[2]。传统水质监测主要采用现场水样采集和实验室测定的方法，虽然可监测的参数多，精度较高，但费时费力，经济成本高昂，并且单点尺度的水质难以代表整个水域的水质状况。

　　遥感技术具有范围广、速度快和低成本等优势，满足实时和大尺度的水质监测需求[3]，同时可以揭示常规方法难以发现的污染物和污染物的迁移和分布特征。水质遥感监测通过研究卫星遥感数据和实测水质数据之间的关系，构建水质参数的反演模型[4]，以此获取整个水域水质的空间分布和变化。

　　水质检测论文范例：农村安全供水水质监测研究

　　1自20世纪70年代以来，水质遥感监测逐渐从定性分析发展到定量反演，建立了分析、经验、半经验和机器学习等多种反演模型，大量新型卫星的发射为水质监测提供了更多遥感数据源，多种水质参数的反演取得了优异成果[5-6]。但是，与海洋遥感监测相比，水域光学复杂性、卫星传感器针对性、大气校正精确性以及其它一些尚未解决的问题给内陆水体水质遥感监测带来了诸多挑战，因此需要进一步促进对内陆水体水质遥感监测研究的理解和应用。本文从遥感数据、反演方法和水质参数三方面阐述内陆水体水质遥感监测研究的研究进展，提出当前水质遥感监测研究中存在的主要问题以及未来工作需要关注的重点，以期为内陆水体的水质监测提供参考。

　　1水质遥感监测常用数据源

　　遥感数据源分为多光谱数据源和高光谱数据源。在内陆水体的水质遥感监测中，常用的多光谱数据源包括美国陆地卫星系列(Landsat1-8)、中分辨率成像光谱仪(moderateresolutionimagingspectroradiometer，MODIS)和Sentinel-2。高光谱数据源可以在一定程度上提高水质反演的准确性与精度，但覆盖范围小，监测成本高[7]，因此对大面积的内陆水体进行全面观测存在缺陷。由于Landsat系列、MODIS和Sentinel-2均为免费的多光谱数据源，而国产卫星系列在近几年发展迅速，本文重点阐述这三种遥感数据源和部分国产卫星系列在内陆水体水质监测中的应用。

　　1.1Landsat系列

　　多光谱扫描仪(MultispectralScanner，MSS)是Landsat1-5和Landsat7携带的传感器。Carpenter[8]利用MSS对淡水湖的浊度和藻类色素浓度进行了建模和预测，证明了该数据源可用来反演水体中的水质参数。

　　但MSS影像波段较宽，空间分辨率较低，不适合用于监测叶绿素含量较低、悬浮物含量很高的内陆水体[9]。专题制图仪(ThematicMapper，TM)是Landsat4和Landsat5携带的传感器。自发射至今，TM工作状态良好，数据获取便利，虽然不是专门针对内陆水体水质遥感研究，但仍因其较高的空间、光谱和辐射分辨率在内陆水体水质监测中被广泛使用。Hafeez等[10]通过比较LandsatTM反射率数据和原位反射率数据评估多种机器学习模型的性能，结果表明神经网络模型可用于利用光学复杂的沿海水域的卫星数据改善水质监测和评估。增强型专题制图仪(EnhancedThematicMapperPlus，ETM+)是Landsat7携带的传感器。

　　ETM+因增加了一个全色波段，在红外波段的分辨率更高，数据信息量更丰富。Lin等[11]测试由LandsatETM+和机器学习算法确定湖泊藻类生物量的准确性，结果表明该模型在分析藻类状况的时间变化方面很有效。但Landsat7的扫描行校正器(ScanLinesCorrector，SLC)曾发生故障，致使部分图像出现重叠和丢失[12]，异常数据需要通过模型校正，这限制了其在水质遥感监测中的应用。

　　陆地成像仪(OperationalLandImager，OLI)-热红外传感器(ThermalInfraredSensor，TIRS)是Landsat8携带的传感器。与ETM+相比，OLI增加了一个短波近红外波段，主要用于卷云检测和数据质量评价等，TIRS主要用于地表温度反演和土壤湿度评价等。Landsat8波段组合多，数据应用范围广。Cao等[13]基于Landsat8影像，利用极限梯度增强树模型开发了一种估计浑浊的湖水中叶绿素a的算法，绘制了2013-2018年中国东部数百个大于1km2的湖泊中叶绿素a的时空变化。

　　1.2MODIS

　　MODIS是搭载在美国EOS卫星上最重要的光学传感器，可用于对陆表、生物圈、固态地球、大气和海洋进行长期全球观测[14]。MODIS具有较高的时间分辨率，其中Terra和Aqua卫星分别在一天中的上午和下午过境，观测整个地球表面。MODIS数据接收相对简单，全球免费获取，相比其他卫星遥感数据提供了有助于辐射校正的大气廓线数据，在未来的水质监测中具有巨大的发展潜力。

　　Chang等[15]利用MODIS影像和遗传规划模型研究了美国坦帕湾总磷的时空变化规律，阐明了由自然灾害和人为扰动引起的耦合动态影响。Wang等[16]利用2000—2017年MODIS数据揭示了中国水质清晰度的长期时空变化趋势，探讨了不同地区的湖泊清晰度和经济发展以及气候变化的关系。

　　2水质遥感监测方法

　　20世纪70年代初遥感监测技术开始应用于水体研究，从单一的水域识别逐渐发展到水质参数的定量反演。20世纪80年代后，多光谱数据和高光谱数据的涌现推动了遥感技术的快速进步，随着水质参数的光谱特性和算法研究的不断深入，水质遥感监测方法经历了分析法、经验法、半经验法、机器学习和综合法的发展历程。

　　2.1分析法

　　分析法是以生物-光学模型为基础的水质反演方法。该方法通过辐射传输模型模拟光在大气和水体中的传播过程，利用遥感测得的水体反射率计算水质参数的特征吸收系数和后向散射系数，建立反射光谱与水质参数之间的关系[22]。分析法具有严格的物理意义，模型稳定性好，反演精度较高，适用性强，不需要大量的地面实测水质数据作支撑。

　　但这一方法基于已知水体中水质参数的光谱特性，在初次建立模型时需要测量固有光学量和表观光学量，如Li等[23]在浅水内陆有色可溶性有机物的时空变化研究中测量了CDOM吸收系数和不同波长的遥感反射率等。这些参数需要通过实验手段或野外观测得出，其中部分参数随监测水域的不同而变化[24]，因此分析法在水质监测实际应用中的效果并不理想。

　　2.2经验法

　　经验法是以多光谱遥感数据为基础的水质反演方法。该方法依据遥感影像数据和地面实测数据之间的关系，选取相关系数最高的波段或波段组合用于统计回归分析，获得水质参数的最佳反演模型。Majid等[25]利用经验算法和神经网络算法估算香港沿海地区的叶绿素a和悬浮物浓度，其中用到的波段组合方法包括相加、相减、相乘、相除、平均和取对数。

　　经验法利用简单易用的回归模型构建相对复杂的关系，在一定程度上能提高水质反演的精度，但也存在一些缺陷[26]。首先，经验法容易受到地域和时间的限制，通用性较差;其次，经验法需要以大量实测数据为基础才能达到理想的水质反演精度;此外，经验法在实测数据的影响下只能反演一定浓度范围内的水质参数，超出这个范围后，反演结果的误差明显增大;最后，经验法因缺乏物理依据往往导致遥感影像和水质数据的关系得不到有效保证，模型的准确性难以保障。

　　2.3半经验法

　　半经验法是以高光谱遥感数据为基础的水质反演方法。该方法通过水质参数的光谱特征和实测数据的关系获得最佳波段或波段组合，选取合适的数学统计分析模型(线性回归、指数/对数函数回归、多项式回归、主成分分析和聚类分析等)反演水质参数。林剑远[27]利用半经验法定量反演城市河网水质，建立了四种水质参数的线性、二次多项式、指数和幂模型。半经验法充分利用了水质参数的光谱特征，兼顾了遥感影像数据和实测水质数据之间的统计关系[28]，是目前最常用的水质反演方法。但半经验法仍有较强的时空局限性，针对不同类型、不同地区和不同季节的水体需要寻找合适的数学统计模型，提高模型的通用性。

　　2.4机器学习

　　机器学习是指通过某些算法指导计算机利用已知数据得出适当模型，并利用此模型对新数据进行分析或者预测的过程[29]。近年来，随着人工智能技术的发展，越来越多的研究把机器学习理论融入到水质遥感监测中。水体中复杂的光谱特征导致水质遥感监测在本质上是一个非线性的反演过程。机器学习具有较强的适应性、组织性和容错性，可以通过持续不断的学习、校正和误差反馈提升遥感反演模型的精度和泛化能力，正适合模拟遥感影像和水质参数间错综复杂的关系。

　　目前，应用于水质反演的机器学习模型包括随机森林[30]、支持向量机[31]、神经网络[32]等。机器学习在许多水质遥感监测研究中表现优异，但构建遥感反演模型需要大量的训练样本，同时校正参数会增加模型训练的时间成本，因此如何平衡模型复杂程度和计算效率至关重要[33]。此外，机器学习的理论基础有待完善，模型结果的可解释性有待提高，模型的普遍适用性有待增强。

　　2.5综合法

　　单一的水质遥感监测方法均存在各自的局限性，如叶绿素a的经验模型通常基于与其光吸收和发射特性相关的物理原理，一般依赖光谱波段中的绿色、蓝色、红色和近红外波段[34]。但不同水体的物理特征、组成成分和藻类种群具有很大差异，这些因素都会影响经验模型的适用性。综合法是指通过比较或结合几种甚至更多的水质遥感监测方法，发挥每种水质监测方法的优势，在充分利用水质参数光谱特征的基础上提高水质反演精度，增强模型的通用性。

　　Hansen[35]使用经验法和多元线性逐步回归方法建立叶绿素a分季节的遥感反演模型，结果表明与传统的经验模型相比，数据驱动的机器学习方法在模型拟合上有良好的改进，利用其它波段所含信息有助于解释湖泊特定的光学特性，如悬浮物、水色组分和季节特定的藻类种群。Pahlevan等[36]通过高光谱数据反演不同区域的叶绿素a，选择传统的固有光学量和叶绿素a算法与混合密度网络算法做对比，显著改进了叶绿素a的反演结果，证明了混合密度网络模型产生内陆和沿海水域高精度水生遥感产品的潜力。

　　3水质遥感监测参数

　　根据水中物质对光辐射的吸收和散射性质，水质遥感监测参数分为光敏参数和非光敏参数。光敏参数是指有明显光学特性的水质参数，主要包括叶绿素a、悬浮物和有色可溶性有机物等;非光敏参数是指无明显光学特性的水质参数，主要包括化学需氧量(COD)、生化需氧量(BOD)、总磷(TP)和总氮(TN)等。

　　3.1叶绿素a

　　叶绿素a是衡量水体富营养化程度的重要指标之一，不仅体现浮游植物的生物量和初级生产力[37]，也影响水体的反射光谱特征。国内外学者针对不同研究区域建立了许多叶绿素a的遥感反演模型[38-39]，但叶绿素a具有明显的季节变化特征，同时悬浮物和有色可溶性有机物的存在会使水体光谱特征发生显著变化，掩盖叶绿素a的光谱吸收特性[40]，因此这些模型仍存在精度不够高和适用性不强的问题。

　　4存在的问题

　　当前，我国地表水环境质量总体保持持续改善的势头，但从水生态环境保护的整体性来看，不平衡不协调的问题依然突出。水质遥感监测在水生态环境保护中表现出广阔的应用前景，但内陆水体水质遥感监测的基础研究仍然不够充分，理论和方法也不够成熟，尚有许多问题存在：a.卫星传感器的分辨率问题。不同遥感数据源在时间分辨率、空间分辨率、光谱分辨率和辐射分辨率上各有优势和缺陷，这对水质参数提取的精度有着重要影响。内陆水体面积较小，组成成分复杂，受人为活动的影响比较严重，需要长期、连续和高精度的动态监测。

　　b.水质反演算法的时空局限性。内陆水体光学特性具有较强的区域性和季节性特点，即使同一水域也有光学深水区和光学浅水区之分。虽然许多遥感反演算法大幅提升了内陆水体水质反演的精度，但因缺乏理论基础导致可解释性较差，仍然有较强的时空局限性。c.水质参数的反演种类和反演机理问题。光敏参数的研究已经日趋成熟，但水体中叶绿素a、悬浮物和CDOM间相互的光谱特征干扰仍受到广泛关注，此外非光敏参数的反演种类依然较少，反演机理尚不明晰，反演模型仍需发展，模型精度也有待提高。

　　d.大气校正模型问题。国内外学者提出了很多大气校正模型，但利用标准大气剖面数据直接进行大气校正的误差较大，且不同地区的大气剖面状态不同。因此，内陆水体的大气校正需要考虑气溶胶光学特性的复杂性和水面反射光校正的不确定性。e.特殊类型水体的水质监测问题。内陆水体的水环境问题依旧存在，针对黑臭水体、富营养化水体和浑浊度较高水体的水质反演仍是水质遥感监测研究的重点和难点。这些水体光谱特性复杂，水质影响因素众多，水质参数的遥感反演进展缓慢，遥感监测的精度较低。

　　5研究展望

　　近年来，水生态健康受到高度重视。生态环境部在2019年启动了重点流域“十四五”规划编制工作，印发了《重点流域水生态环境保护“十四五”规划编制技术大纲》，将重点流域规划名称由“水污染防治”调整为“水生态环境保护”，体现了新时期流域生态环境保护工作的新要求。利用遥感技术监测内陆水体的水质状况为水生态环境保护提供了新思路和新方法，不仅节省了大量人力、物力和财力，同时也帮助了环境保护部门高效掌握和监管水生态环境的实时状况。

　　目前，基于遥感的内陆水体水质监测研究仍有需要补充和加强的地方，未来应重点关注以下几个方面：

　　a.深入采用更多的新型遥感数据，融合多种遥感源以实现不同时空尺度下的水质反演。未来应充分结合微波遥感与可见光或红外数据，加强高光谱技术在内陆水体水质遥感监测中的应用，发展专门针对内陆水体(如湖泊、水库等)水质遥感监测的传感器。

　　b.生物-光学模型是建立在光学传输物理过程之上的通用模型，未来应深入研究水质参数的内在光学特性与表面反射率或离水辐亮度之间的理论关系，将生物-光学模型和经验法、半经验法、机器学习模型等相结合，发展没有时间和空间特殊性的反演模型。

　　c.深入研究水体中不同组分的光谱特征及其差异，发现各水质参数的光谱响应曲线特征，了解不同水质参数之间光谱相互影响的规律，同时扩大水质参数的监测种类，增加非光敏参数的可行性分析和定量遥感监测，建立不同水质参数的光谱特征数据库。

　　d.内陆水体大气校正算法应区分大气和水体物质对传感器总信息贡献的解耦方法，考虑水域上空不同类型气溶胶潜在的复杂混合，同时减少或消除水面反射光的干扰。未来应加强对水体光谱特性和机理模型的研究，发展针对内陆水体水质遥感的精确大气校正模型。

　　e.内陆水体的水质遥感监测应先通过水体分类构建不同时空尺度下的反演模型，再逐步扩大研究区域和研究对象，获得精度高和应用广的统一模型，最终形成一套完整的内陆水体遥感监测体系，为未来的水质遥感监测奠定坚实的基础。

　　参考文献：

　　[1]季永兴,韩非非,施震余,等.长三角一体化示范区水生态环境治理思考[J].水资源保护,2021,37(1):103-109.(JIYongxing,HANFeifei,SHIZhenyu,etal.Thinkingonwatereco-environmentmanagementinregionalintegrationdemonstrationareaofYangtzeRiverDelta[J].WaterResourcesProtection,2021,37(1):103-109.(inChinese))

　　[2]蓝志峰,肖桂荣.基于流域要素空间关系的水污染溯源研究[J].水资源保护,2019,35(1):56-62.(LANZhifeng,XIAOGuirong.Waterpollutiontraceabilitybasedonspatialrelationshipofbasinelements[J].WaterResourcesProtection,2019,35(1):56-62.(inChinese))

　　[3]LIJ,PEIY,ZHAOS,etal.AreviewofremotesensingforenvironmentalmonitoringinChina[J].RemoteSensing,2020,12(7):1130.

　　[4]XIONGY,RANY,ZHAOS,etal.RemotelyassessingandmonitoringcoastalandinlandwaterqualityinChina:Progress,challengesandoutlook[J].CriticalReviewsinEnvironmentalScienceandTechnology,2019(8):1-37.

　　作者：王波1，黄津辉1，郭宏伟1，许旺2，曾清怀2，麦有全2，祝晓瞳1，田上1