一种二类水体固有光学参数优化模拟方法

1.本发明涉及一种二类水体固有光学参数优化模拟方法，面向我国典型的内陆二类湖泊太湖，基于野外实地采样数据以及2seacolor模型优化模拟得到的研究区固有光学参数，实现水色要素高精度同步反演。


背景技术：

2.水资源短缺、水生态损害、水环境污染等新问题日益突出，基于“三水统筹”的水生态安全管理面临新的挑战。水生态监测是水生态环境管理的“顶梁柱”，是生态文明建设的重要基础支撑之一。生态环境部发布的《关于开展生态环境遥感监测试点工作的通知》(监测函〔2019〕6号)要求试点省份加快构建生态环境遥感调查、监测与评估体系。《重点流域水生态环境保护规划(2021-2025年) (征求意见稿)》明确指出，要基于流域特色要素，制定统筹水资源、水生态、水环境的指标体系，以高水平保护引导推动高质量发展。
3.地表水体分为一类水体和二类水体。大洋开阔水体主要为一类水体，近岸河口水体主要为二类水体。由于二类水体更接近人类生产生活区，所以受人类活动影响更大，污染严重。近年来，我国多个湖泊出现富营养化、有机污染、水域面积缩小、盐渍化、生态系统失衡等问题(程琦，2012)。水体富营养化不但会破坏湖泊的整体功能，破坏湖泊生态系统的多样性，而且还会危害浮游植物，导致水华暴发。藻类产生的毒素通过饮用水厂影响人类健康，威胁人类生存，水面上的藻类散发出异味，严重影响河湖景观。
4.长期以来，近岸二类水体水色三要素的反演一直是水色光学遥感中的一个难题。水色三要素主要包括叶绿素(chlorophyll，以chlorophyll-a为主chl-a)、非藻类固体悬浮物(suspended particulate matter,spm)和有色可溶性有机物 (chromophoric dissolved organic matter or coloreddissolved organic material, cdom)(唐军武等，2003)。目前水色三要素的遥感监测方法主要分为经验方法、半经验/半分析方法和分析方法三种，从辐射传输角度进行反演的模型很少，往往以单组分水色要素反演模型为主，多组分水色要素同步反演模型少。
5.鉴于目前二类水体水色三要素遥感监测方法存在的不足，本发明将研究焦点集中在实测的叶绿素a浓度、悬浮物浓度、黄色物质在440nm处的吸收系数和水面光谱的关系上，通过水色三要素与遥感反射率建立数学关系，提出一种二类水体固有光学参数优化模拟方法，进而有效地弥补传统监测方法在时间和空间尺度上的不连续性，提供连续、准确的水质参数变化数据。


技术实现要素：

6.针对当前水色三要素遥感监测缺乏多组分水色要素同步反演方法、存在信息多散、汇集复杂、分析获取不及时等问题，本发明提出了一种二类水体固有光学参数优化模拟方法，为水生态保护提供可靠有效准确的信息源。
7.本发明的目的通过以下技术步骤实现：
8.步骤1)选取研究区并设计试验方案进行采样点处叶绿素a浓度、悬浮物浓度、黄色物质在440nm处的吸收系数以及水面光谱测量；
9.步骤2)计算各采样点的吸收系数和后向散射系数；
10.步骤3)利用2seacolor辐射传输模型模拟各采样点遥感反射率；
11.步骤4)计算模拟出的遥感反射率和实测采样点的遥感反射率均方根差并确定研究区固有光学参数值。
12.所述步骤1)的具体方法为：
13.a)选取研究区域；b)基于研究区设计湖面试验方案；c)水质参数实地测量。
14.所述步骤2)的具体方法为：
15.a)计算研究区各采样点的吸收系数：
16.a(λ)＝aw(λ)+a
chla
(λ)+as(λ)+a
cdom
(λ)
17.式中，a(λ)指波长λ处的水体总吸收系数；aw(λ)为水分子吸收系数、a
chla
(λ)为叶绿素a吸收系数、as(λ)为悬浮物的吸收系数、a
cdom
(λ)为黄色物质所产生的吸收系数，单位都是每米(m-1
)。
18.b)计算研究区各采样点的后向散射系数：
19.bb(λ)＝bw(λ)+b
chla
(λ)+bs(λ)
20.式中，bb(λ)指波长λ处的水体总散射系数；bw(λ)水分子散射系数；b
chla
(λ)为叶绿素a散射系数；bs(λ)为悬浮物的散射系数。
21.所述步骤3)的具体方法为：
22.a)将研究区的固有光学参数叶绿素a浓度的指数项x、悬浮物散射光谱斜率ns、和悬浮物在550nm处的散射系数bs(550)三个要素作为未知量；b)将三个未知量设定合理的初值、范围和步长，利用2seacolor模型不断循环得到对应于每一组x、ns、和bs(550)的遥感反射率。
23.其中，2seacolor模型的理论公式为(salama等，2015)：
[0024][0025][0026]
μw＝cosθ
′s[0027][0028][0029]
式中，是半无限介质的定向半球面反射率；
[0030]
x是水体总后向散射系数bb(单位m-1
)和总吸收系数a的比值(单位m-1
)；
[0031]
μw指平坦水面下太阳天顶角的余弦值；若水面上太阳天顶角为θs，则θ
′s＝arcsin(sinθs/nw)，其中nw是水的折射系数，取值1.33；
[0032]
r(0-)是水平面下辐亮度反射率；
[0033]
q是水平面向上的辐亮度反射率与水平面向下的辐亮度反射率的比值，取值 3.25；
[0034]rrs
为离水反射率。
[0035]
所述步骤4)的具体方法为：
[0036]
a)计算每一个采样点模拟得到的遥感反射率与同步实测的遥感反射率的均方根差(rmse)；b)找出最小的rmse，确定研究区的固有光学参数值。
附图说明
[0037]
图1为本发明实施的方法流程示意图；
[0038]
图2为试验采样点分布图；
[0039]
图3为模拟出的各采样点的遥感反射率图。
具体实施方式
[0040]
下面结合附图对本发明“一种二类水体固有光学参数优化模拟方法”作进一步阐述说明。
[0041]
(一)试验设计与水体样本采集
[0042]
首先，选取研究区域(以近岸二类水体太湖为例)；其次，在研究区范围内设计覆盖整个湖面的试验采样点，如图2所示。在太湖北部梅梁湾、竺山湾等水污染较严重区域加密采样点，在太湖南部较清洁水体区域适当减少采样点；最后，采用水样瓶进行水体采样用于实验室测量叶绿素a浓度、悬浮物浓度和黄色物质在440nm处的吸收系数以及采用fieldspec handheld2手持便携式光谱分析仪进行光谱数据采集。
[0043]
(二)吸收系数和散射系数的计算
[0044]
水体中的吸收系数考虑水体中的四种组分，即水分子、叶绿素a、悬浮物和黄色物质，其中纯水分子的吸收系数是波长λ的函数。基于研究区采样点实测的叶绿素a浓度、悬浮物浓度以及440nm处黄色物质的吸收系数，可以求得研究区总吸收系数。
[0045]
其中，叶绿素a的吸收系数的计算公式为：
[0046][0047]achla
(λ)＝(a0(λ)+a1(λ)ln a
chla
(440))
×achla
(440)
[0048]
式中，c
chla
为叶绿素a浓度，a0(λ)、a1(λ)为经验参数，a
chl
(440)为太湖叶绿素a波长440nm处的光谱吸收系数。
[0049]
悬浮物的吸收系数的计算公式为：
[0050][0051]
式中，cs为悬浮物浓度，as(440)为太湖悬浮物波长440nm处的光谱吸收系数，为悬浮物的归一化单位吸收系数。
[0052]
黄色物质的吸收系数的计算公式为：
[0053]acdom
(λ)＝a
cdom
(440)exp[-s
cdom
(λ-440)]
[0054]
式中，s
cdom
为黄色物质吸收光谱斜率，本研究根据实测数据确定的平均值为
0.015nm-1
。
[0055]
由于黄色物质对光的后向散射作用非常微弱，因此后向散射系数的计算只考虑水体中的三种组分，即水分子、叶绿素和悬浮物，其中纯水分子的散射系数是波长λ的函数。
[0056]
其中，叶绿素a的散射系数的计算公式为：
[0057][0058]
式中，为叶绿素a的归一化单位吸收系数，x为叶绿素a浓度的指数项，其优化范围为(0.01，2.00)，步长为0.01。
[0059]
悬浮物的散射系数的计算公式为：
[0060][0061]
式中，bs(550)为悬浮物550nm处的散射系数，设置其优化范围为(0.01，2.00)，步长为0.01；ns为悬浮物散射光谱斜率，设置其优化范围为(-1.00，1.00)，步长为0.01。
[0062]
将x、ns、和bs(550)三个研究区固有光学参数作为未知数，通过设置的初始值、范围和步长的情况下不断进行循环计算，得到与x、ns、和bs(550)相对应的各个采样点总散射系数的值。
[0063]
(三)遥感反射率模拟结果与固有光学参数确定
[0064]
基于步骤(二)中计算的吸收系数、散射系数以及本发明基于2seacolor 模型模拟得到的400-700nm波长范围的遥感反射率与实测遥感反射率进行均方根差计算。在均方根误差达到最小时，确定对应的x为0.41、ns为-1.15、bs(550) 为0.03即为适合研究区水色三要素反演的固有光学参数值。以波长为横轴，模拟得到的遥感反射率为纵轴作遥感反射率随波长变化曲线图，选取代表性采样点数据成图，如图3所示；
[0065]
其中，均方根差的计算公式为：
[0066][0067]
式中，i为遥感反射率波段数，为模拟出的遥感反射率，为实测的遥感反射率。