富营养化水库表层水体藻类分布Landsat遥感监测方法
【专利摘要】本发明涉及一种富营养化水库表层水体藻类分布Landsat遥感监测方法，基于Landsat遥感数据，结合浮游藻类指数与实测藻类生物量的拟合模型，构建藻类水华分布的决策树，识别富营养化水库表层水体藻类生物量以及时空分布特征。本发明步骤简单，可操作性强，能够大范围的监测水体表层藻类的分布情况，且藻类信息提取较快，不受地面条件限制。为富营养化水库管理提供技术基础，以期为水库中藻类水华的防治提供科学依据。
【专利说明】
富营养化水库表层水体藻类分布Landsat遥感监测方法
技术领域
[0001] 本发明涉及富营养化水库表层水体藻类分布遥感监测的开发技术领域，基于 Landsat遥感图像分析，开发适应于富营养化水库藻类快速监测的监测方法。
【背景技术】
[0002] 自20世纪60年代以来，世界各地区的湖泊、水库、河流和海湾相继出现了水体富营 养化问题，引起了科学家的广泛关注和研究。0ECD对18个国家的湖泊和水库的调查结果显 示，富营养型所占比例高达65%，仅18%为贫营养型。其中，在美国就有70%的湖泊和水库 为富营养型。根据欧洲环境署公布的信息，自1970年以来监测调查的湖泊和水库中富营养 化和重富营养化所占比例均在40%以上;对401个河流监测站1986~1996年的数据进行分 析，发现33%的站点处于重富营养状态。除了上述地区，澳大利亚、日本、墨西哥、巴西、阿根 廷、印度、伊朗、南非、埃及等众多国家也都发生过严重的水体藻华事件乃至藻毒素中毒事 件。可以说，水体富营养化问题已成为21世纪全球环境可持续发展面临的巨大挑战之一。
[0003] 随着水域营养化问题的加剧，湖泊水库中藻类水华频繁暴发，生态灾害事件频繁 发生，严重影响了水域生态系统功能的发挥，制约了区域经济的可持续发展。为了能够控制 水体中藻类水华大规模爆发，减少区域经济损失，实现快速探测到水体中藻类的浓度以及 时空分布特征具有重要意义。传统意义上的藻类监测依赖于水体实地采样，固定，显微镜观 测，但是此方法检测时间长，且无法识别藻类的时空分布情况。因此为了解水体中藻类密度 和时空分布特征，控制藻类水华大规模爆发，建立快速、有效的氨基酸测试方法显得尤为重 要。
【发明内容】
[0004] 本发明的目的在于提供一种富营养化水库表层水体藻类分布Landsat遥感监测方 法，利用遥感技术精确获取水库表层水体中藻类的浓度表征指标的分布特征，继而准确估 算藻类的浓度范围以及分布区域，分析藻类水华的发生状况及趋势，评估水库的富营养化 状况，为水库的管理提供科学依据。
[0005] 为了实现本发明的技术目的，采用如下技术方案：
[0006] -种富营养化水库表层水体藻类分布Landsat遥感监测方法，包括:通过野外调查 采样分析获得表层水体藻类生物量及其空间分布情况;基于水库水体和藻类的光谱特性， 运用浮游藻类指数(FAI)提取水库表层水体的藻类分布情况，结合实测的水库表层水体中 藻类生物量，构建水库表层水体中藻类分布的遥感监测方法。基于该方法，可以准确获取富 营养化水库表层水体藻类生物量的月际变化规律及其空间分布特征。
[0007] 作为进一步的实施方式，具体的方法包括如下步骤：
[0008] (1 )、确定富营养化水库表层水体的藻类生物量
[0009] 现场使用HQM-1有机玻璃采水器在距离水体表面0.5处取500mL水样，放入聚乙烯 塑料瓶中，带回实验室储存于4°C的冰箱中，在48h内运用热乙醇法测定水体表层的叶绿素a (Chl-a)用于表征水库表层水体中藻类生物量。具体的测定方法为:取250ml水样过滤，滤膜 向内对折放入5ml离心管；取250ml玻璃三角瓶装适量90%乙醇在控温水浴锅预热，预热温 度为85°C;取出样品，立即加入4ml热乙醇，水域2min后，将样品避光6h;，用25mm玻璃纤维滤 膜过滤萃取液定容到l〇ml;在分光光度计上用90 %乙醇作参比液进行比色，先在665nm波长 测消光率E665，再在750nm波长测消光率E75Q，然后在样品比色皿中加1滴lmo 1 /L盐酸进行酸 化，加盖摇匀，lmin后重新在665nm波长测消光率A665 ,再在750nm波长测消光率A75Q;计算 Chi-a浓度。
[0010] Chl_a = 0 · 446*[ (E665-E750)-(A665-A750)]
[0011] (2)、富营养化水库表层水体藻类分布Landsat遥感监测方法构建
[0012] 为了消除遥感图像在成像过程中产生的扭曲，变形等问题，对Landsat遥感图像进 行条带噪声去除，波段融合，辐射定标和大气校正处理;计算水库表层水体的FAI指数值;结 合实测的水库表层Chl-a浓度值，进行多元回归分析，构建水库表层水体藻类监测模型。
[0013] 浮游藻类指数FAI
[0015] 式中Rrc;,B3，Rrc;,B4，R rc;,B5分别为LandsatETM+遥感图像经大气校正后的波段3，波段 4和波段5 的反射率值，其中 λΒ3 = 660nm; λΒ4=830nm; λΒ5 = 1650nm 〇
[0016] (3)、基于水库表层水体藻类的遥感监测模型构建藻类的时空分布特征
[0017] 根据以往研究，划定藻类水华等级为Type 1-藻类低浓度区（Oyg/L〈Chl-a〈2〇yg/ 〇，了7口6 2-藻类中浓度区（2(^8凡〈〇11-3〈5(^8凡），17口6 3-藻类高浓度区（5(^8凡〈〇11-&〈 155yg/L)，Type 4-藻类水华重灾区（Chl-a>155yg/L)。结合上文中得到的藻类监测模型，获 得FAI指数阈值。基于上述步骤和方法，建立富营养化水库表层水体藻类分布的决策树，识 别藻类水华的空间分布。
[0018] 本发明与传统富营养化水库表层水体藻类生物量的检测方法相比，能够实现对富 营养化水库藻类进行大范围监测，快速获得表层水体中藻类信息，且不受地面条件限制，步 骤简单，可操作性强。此外，本发明针对中富营养化程度水库藻类生物量及分布特征的监测 效果明显。
【附图说明】
[0019] 图1为实例2中青草沙水库表层水体Chl-a浓度时空分布图
[0020] 图2为实例2中青草沙水库叶绿素a与FAI逐步回归模型 [0021]图3为实施例2中青草沙水库2014年-2015年藻类水华分布图
[0022] 图中标号：水体-蓝色，Typel-绿色，Type2_黄色，Type3_红色，Type4_紫色。
【具体实施方式】
[0023] 本发明涉及一种富营养化水库表层水体藻类分布Landsat遥感监测方法，基于 Landsat遥感数据，结合浮游藻类指数与实测藻类生物量的拟合模型，构建藻类水华分布的 决策树，识别富营养化水库表层水体藻类生物量以及时空分布特征。本发明步骤简单，可操 作性强，能够大范围的监测水体表层藻类的分布情况，且藻类信息提取较快，不受地面条件 限制。为富营养化水库管理提供技术基础，以期为水库中藻类水华的防治提供科学依据。下 面通过实施例进一步说明本发明，但本发明的保护范围不限于所述内容。
[0024] 实施例1(理论依据）
[0025] -种富营养化水库表层水体藻类分布Landsat遥感监测方法，包括:通过野外调 查采样分析获得表层水体藻类生物量及其空间分布情况;基于水库水体和藻类的光谱特 性，运用浮游藻类指数(FAI)提取水库表层水体的藻类分布情况，结合实测的水库表层水体 中藻类生物量，构建水库表层水体中藻类分布的遥感监测方法。基于该方法，可以准确获取 富营养化水库表层水体藻类生物量的月际变化规律及其空间分布特征。
[0026]具体的方法包括如下步骤：
[0027] (1)、确定富营养化水库表层水体的藻类生物量
[0028] 现场使用HQM-1有机玻璃采水器在距离水体表面0.5处取500mL水样，放入聚乙烯 塑料瓶中，带回实验室储存于4°C的冰箱中，在48h内运用热乙醇法测定水体表层的叶绿素a (Chl-a)用于表征水库表层水体中藻类生物量。具体的测定方法为:取250ml水样过滤，滤膜 向内对折放入5ml离心管；取250ml玻璃三角瓶装适量90%乙醇在控温水浴锅预热，预热温 度为85°C;取出样品，立即加入4ml热乙醇，水域2min后，将样品避光6h;，用25mm玻璃纤维滤 膜过滤萃取液定容到l〇ml;在分光光度计上用90 %乙醇作参比液进行比色，先在665nm波长 测消光率E665，再在750nm波长测消光率E75Q，然后在样品比色皿中加1滴lmo 1 /L盐酸进行酸 化，加盖摇匀，lmin后重新在665nm波长测消光率A665 ,再在750nm波长测消光率A75Q;计算 Chi-a浓度。
[0029] Chl_a = 0 · 446*[ (E665-E75Q)-(A665-A750)]
[0030] (2)、富营养化水库表层水体藻类分布Landsat遥感监测方法构建
[0031] 为了消除遥感图像在成像过程中产生的扭曲，变形等问题，对Landsat遥感图像进 行条带噪声去除，波段融合，辐射定标和大气校正处理;计算水库表层水体的FAI指数值;结 合实测的水库表层Chl-a浓度值，进行多元回归分析，构建水库表层水体藻类监测模型。 [0032]浮游藻类指数FAI
[0034] 式中Rrc,B3，Rrc,B4，Rrc,B5分别为LandsatETM+遥感图像经大气校正后的波段3，波段 4和波段5 的反射率值，其中 λΒ3 = 660nm; λΒ4=830nm; λΒ5 = 1650nm 〇 [0035] (3)、基于水库表层水体藻类的遥感监测模型构建藻类的时空分布特征
[0036]根据以往研究，划定藻类水华等级为Type 1-藻类低浓度区（Oyg/L〈Chl-a〈2〇yg/ 〇，了7口6 2-藻类中浓度区（2(^8凡〈〇11-3〈5(^8凡），17口6 3-藻类高浓度区（5(^8凡〈〇11-&〈 155yg/L)，Type 4-藻类水华重灾区（Chl-a>155yg/L)。结合上文中得到的藻类监测模型，获 得FAI指数阈值。基于上述步骤和方法，建立富营养化水库表层水体藻类分布的决策树，识 别藻类水华的空间分布。
[0037]实施例1(实例及验证）
[0038] 选取青草沙水库（31° 2Y N-31。25' N，121° 35' E-120。43' E)为研究对象，针对本发 明进行研究和验证。青草沙水库位于长江口中心位置，总面积约67.2平方公里，河段长约 180km，平均水深8m，总库容为6.44亿m3。在2014年1月-2015年11月期间，用麦哲伦GPS在青 草沙水库布设S1-S14采样点（表1)，使用HQM-1有机玻璃采水器在距离水体表面0.5处取 500mL水样，用热乙醇法测定各采样点的Ch 1 -a浓度。
[0039] 对研究期间有效的Landsat遥感图像进行条带噪声去除，波段融合，辐射定标和大 气校正，其中大气校正选取FLAASH大气校正模型，在ENVI软件FLAASH模块中，采用中炜度夏 季标准大气模型（mid-latitude summer atmosphere model)和城市气溶胶模型（burden aerosol model)，传感器高度选用649.093km，地面高程选用3.4m进行相应的图像处理。
[0040] 经过遥感图像预处理后，计算研究区域的FAI指数值，并识别Landsat遥感图像上 S1-S14的FAI值，对实测Chl-a浓度值与FAI值进行多元回归分析，构建藻类分布Landsat遥 感监测模型，y = 0.126x-0.5,其中y为ln(Chl-a);x为FAI值。
[0041] 表1青草沙水库采样点经炜度分布表
[0044] 基于该模型以及FAI指数阈值，分别为了7?61(-0.4384〈？六1〈-0.1225);了 7?62(- 0 · 1225〈FAI〈-0.0071); Type3(-0.0071〈FAI〈0.1355); Type4(FAI>0.1355)。基于上述步骤 和方法，建立青草沙水库表层水体藻类分布的决策树，识别青草沙水库藻类水华的空间分 布。
【主权项】
1. 一种富营养化水库表层水体藻类分布Landsat遥感监测方法，其特征在于，具体步骤 如下： (1 )、确定富营养化水库表层水体的藻类生物量 通过野外实地监测获得水库表层水体藻类生物量的浓度(Chl-a)，W及藻类的空间分 布特征； (2) 、富营养化水库表层水体藻类分布Landsat遥感监测方法构建 为了消除遥感图像在成像过程中产生的扭曲，变形等问题，对Landsat遥感图像进行条 带噪声去除，波段融合，福射定标和大气校正处理;计算水库表层水体的FA巧旨数值;结合实 测的水库表层化1-a浓度值，进行多元回归分析，构建水库表层水体藻类监测模型； 浮游藻类指数FAI式中Rrc,B3，Rrc,B4，Rrc,化分别为LandsatETM+遥感图像经大气校正后的波段3，波段4和 波段5 的反射率值，其中 λΒ3 = 660nm; λΒ4=830nm; λβ日=1650nm; (3) 、基于水库表层水体藻类的遥感监测模型构建藻类的时空分布特征 根据划定藻类水华等级Type 1-藻类低浓度区（0yg/L<^l-a<20yg/L)，Type 2-藻类中 浓度区（2化邑/1<〇^-曰<504邑/1)，了7口6 3-藻类高浓度区（5化邑/1<化1-曰<1554邑/1)，了7口6 4- 藻类水华重灾区（化l-a〉155yg/L)，并结合所述的藻类监测模型，获得FA巧旨数阔值，建立富 营养化水库表层水体藻类分布的决策树，识别藻类水华的空间分布。
【文档编号】G06K9/00GK106092929SQ201610398656
【公开日】2016年11月9日
【申请日】2016年6月7日
【发明人】李灵慧, 龚庆碗, 吴婷, 金毅麟, 李建华
【申请人】同济大学