发明名称 华南沿海赤潮灾害应急监测与预警方法
摘要 本发明涉及一种华南沿海赤潮灾害应急监测与预警方法,包括:基于相应水域光谱数据得到反射率数据,利用得到的反射率数据计算荧光高度并采取6波段算法提取水域颗粒物后向散射系数,基于得到的水域颗粒物后向散射系数再次采用6波段算法估算出叶绿素浓度[chl]’,最后根据所得的荧光高度值、颗粒物后向散射系数和叶绿素浓度值判读赤潮爆发几率。本发明所述的方法对华南沿海各地的赤潮爆发几率的判别能力较好,并可以为应对赤潮事件提供决策支持。
申请公布号 CN105445233A 申请公布日期 2016.03.30
申请号 CN201510761278.5 申请日期 2015.11.09
申请人 国家海洋局南海预报中心 发明人 杨超宇;白毅平;罗军
分类号 G01N21/53(2006.01)I;G01N21/55(2014.01)I 主分类号 G01N21/53(2006.01)I
代理机构 北京兆君联合知识产权代理事务所(普通合伙) 11333 代理人 刘俊玲
主权项 华南沿海赤潮灾害应急监测与预警方法,其特征在于,包括以下步骤:A:按照所需时间、区域测量得到相应水体光谱数据后计算得到相应水体的反射率数据;B:利用步骤A所得反射率数据计算荧光高度(FLH),其中FLH的计算采用667nm、678nm、748nm三个波段进行比值计算;计算FLH的具体方法如下式(1)所示:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>F</mi><mi>L</mi><mi>H</mi><mo>=</mo><msub><mi>nL</mi><mi>w</mi></msub><mrow><mo>(</mo><mn>678</mn><mo>)</mo></mrow><mo>-</mo><msub><mi>nL</mi><mi>w</mi></msub><mrow><mo>(</mo><mn>667</mn><mo>)</mo></mrow><mo>-</mo><mo>&lsqb;</mo><msub><mi>nL</mi><mi>w</mi></msub><mrow><mo>(</mo><mn>748</mn><mo>)</mo></mrow><mo>-</mo><msub><mi>nL</mi><mi>w</mi></msub><mrow><mo>(</mo><mn>667</mn><mo>)</mo></mrow><mo>&rsqb;</mo><mfrac><mrow><mn>678</mn><mo>-</mo><mn>667</mn></mrow><mrow><mn>748</mn><mo>-</mo><mn>667</mn></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000842528210000011.GIF" wi="1653" he="127" /></maths>C:利用步骤A所得反射率数据采取6波段算法提取海水颗粒物后向散射系数,具体步骤如下:步骤A得到的已经过大气校正的反射率数据r<sub>rs</sub>表达为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>r</mi><mrow><mi>r</mi><mi>s</mi></mrow></msub><mo>=</mo><msub><mi>g</mi><mn>1</mn></msub><mfrac><msub><mi>b</mi><mi>b</mi></msub><mrow><mi>a</mi><mo>+</mo><msub><mi>b</mi><mi>b</mi></msub></mrow></mfrac><mo>+</mo><msub><mi>g</mi><mn>2</mn></msub><msup><mrow><mo>(</mo><mfrac><msub><mi>b</mi><mi>b</mi></msub><mrow><mi>a</mi><mo>+</mo><msub><mi>b</mi><mi>b</mi></msub></mrow></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000842528210000012.GIF" wi="1254" he="150" /></maths>其中g<sub>1</sub>、g<sub>2</sub>分别为0.0949和0.0797;a是海水总吸收系数,b<sub>b</sub>是海水总后向散射系数;其中a包含纯海水吸收、CDM吸收和浮游植物吸收的总和;其中浮游植物吸收可以表达为:a<sub>ph</sub>(λ)=a<sup>*</sup>[chl]         (3)其中,a<sup>*</sup>是叶绿素比吸收系数(m<sup>2</sup>/mg);[chl]是叶绿素浓度(mg/m<sup>3</sup>);(这里,a<sup>*</sup>在412,443,488,531,555,667nm),参数选取参见表1表1:模型参数<img file="FDA0000842528210000013.GIF" wi="1422" he="533" />CDOM吸收与波长存在一个指数关系,表达为:a<sub>CDM</sub>(λ)=a<sub>CDM</sub>(443)exp[‑S(λ‑443)]        (4)S是a<sub>CDM</sub>光谱斜率,可以通过反射率比值计算得出:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>s</mi><mo>=</mo><mn>0.015</mn><mo>+</mo><mfrac><mn>0.002</mn><mrow><mn>0.6</mn><mo>+</mo><mfrac><mrow><msub><mi>r</mi><mrow><mi>r</mi><mi>s</mi></mrow></msub><mrow><mo>(</mo><mn>443</mn><mo>)</mo></mrow></mrow><mrow><msub><mi>r</mi><mrow><mi>r</mi><mi>s</mi></mrow></msub><mrow><mo>(</mo><mn>555</mn><mo>)</mo></mrow></mrow></mfrac></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000842528210000014.GIF" wi="1236" he="175" /></maths>所述的海水总后向散射系数由纯水后向散射与颗粒物后向散射系数计算得出,其中颗粒物后向散射系数通过下式计算得出:b<sub>bp</sub>(λ)=b<sub>bp</sub>(λ<sub>0</sub>)(λ/λ<sub>0</sub>)<sup>‑Y</sup>          (6)Y是后向散射光谱指数,在这里取常数1(λ<sub>0</sub>=443nm)固定参数S与Y,将式(3)‑(6)代入式(2),解非线性方程,估算出叶绿素浓度[chl]、b<sub>bp</sub>(443)和a<sub>CDM</sub>(443);D:将步骤C中提取得到的b<sub>bp</sub>(443)作为输入参数,同时固定S,将所述的式(3)‑(6)再次代入式(2),再次解非线性方程,估算出叶绿素浓度[chl]’、a<sub>CDM</sub>(443)’和Y’;E:根据步骤B所得的荧光高度值和步骤D所得的叶绿素浓度值判读赤潮爆发几率,具体包括:E1.根据步骤B得到的结果做第一级判读:FLH大于等于0.05mW cm‑<sup>2</sup>μm<sup>‑1</sup>sr<sup>‑1</sup>的区域判别为高度疑似,反之判别为低度疑似;E2.根据步骤D得到的叶绿素浓度[chl]’值对步骤F1判别为高度疑似的区域进一步做第二级判读:[chl]’值高于5mg/L的区域判别为赤潮爆发几率高,反之判别为赤潮爆发几率低。
地址 510310 广东省广州市海珠区新港中路353号