一种表层海水中多环芳烃的遥感监测方法与流程
本发明属于海洋环境遥感领域，具体涉及海洋污染物、海洋生态环境监测、海洋水色遥感、水质遥感等领域。该方法将海洋水色参数有色溶解有机物(coloreddissolvedorganicmaterial，cdom，或称黄色物质yellowsubstance)的卫星遥感反演技术拓展到海洋表层水体中多环芳烃
(polycyclicaromatichydrocarbons，简称pahs)的监测系统中，使用遥感手段实现海洋环境中持久性有机污染物的大范围、周期性监测。
背景技术：
多环芳烃(polycyclicaromatichydrocarbons，pahs)是指分子中含有两个或以上苯环的芳香族碳氢化合物，是一类普遍存在于环境中的持久性有机污染物(persistentorganicpollutants，pops)中的一种，具有强烈的致癌、致畸、致突变作用，由于广泛存在、种类繁多以及强烈致癌性而备受关注。目前已发现200多种pahs同系物，是环境中致癌化学物质中最大的一类，在环境中虽然含量较少，但广泛存在于各种环境介质中，主要来源于煤、石油、木材的不完全燃烧和石油产品以及在石油开采、使用过程中产生的泄露和排放，严重危害生态环境和人类健康。20世纪80年代，美国环境保护署(usepa)选定毒性显著的16种列入优先控制污染物名单。
目前对于海洋环境中pahs的监测方式主要采用现场采样和实验室定量分析相结合的方式进行，传统的测量方法需要派遣大量工作人员出海逐个站点进行样品采集，并对样品进行一系列的有机分析预处理工作，再对pahs的含量进行检测，从样品采集、样品处理到样品测定，由于海水中pahs浓度极低，对测量技术要求较高，整个流程耗费相当长的时间，需要投入大量的人力、财力、物力。海洋面积广阔，自然条件恶劣，加之测量费用过高，传统的监测方法效率低，成本高，只能做到点、线监测，无法实现对海洋环境中的pahs大面积监测。同时海洋环境中的pahs随着运动的水体具有明显的动态变化特征，传统的监测方法无法做到短周期、高频率地监测。无论从采样广度还是采样周期，传统的监测方法都无法满足大范围、短周期的监测需求。
遥感具有大尺度、周期性、快速同步获取水体信息、成本相对低廉等优点，随着技术的发展，运用遥感手段可以高效地监测海洋水体中各种组分含量的分布和动态变化，从而克服传统现场观测方法的不足，在水环境检测中的应用越来越广泛，海洋是时刻处于运动变化中的水体，对传感器在时间和空间分辨率上具有更高的要求。地球静止海洋水色成像仪goci(geostationaryoceancolorimager)是世界上第一个为海洋水色观测而设计的静止轨道卫星传感器，搭载在首颗地球静止气象卫星coms
(communication，ocean，meteorologicalsatellite)上，主要任务是观测以朝鲜半岛为中心的一定范围内的海洋环境的变化，对该区域的海洋生态系统进行长期的和短期的监测。每天8景成像的时间分辨率为水色遥感动态监测提供了良好的契机。goci共设置了8个波段，中心波长分别为412nm、443nm、490nm、555nm、660nm、680nm、745nm、865nm。goci的精度非常高，其辐射校正误差小于3.8％，空间分辨率为500m，
时间分辨率为1h，可以较好地捕捉水体的空间、光谱特征。
使用遥感手段可以监测到海洋环境中叶绿素、悬浮物、溶解性有机物等水体水色重要参数，目前对于水体悬浮颗粒物tsm、有色溶解有机质cdom、色素、叶绿素等遥感监测研究较多，技术已相对成熟。而pahs在光谱上表现较差，很难通过遥感手段直接获取或者进行反演，大量基础研究表明，pahs的环境行为与水体中的可遥感组分(如tsm、cdom等)密切相关，为利用遥感技术进行监测提供了理论可能，cdom是水色遥感中三大水质参数之一，以cdom为中介，使用遥感手段间接地实现研究海域pahs大范围、高频次监测，从而可以节约大量的人力、经济成本。
综上所述，为了弥补当前海洋环境中pahs的监测能力和技术手段的不足，本发明以河口海湾为研究区，以表层海水中pahs的观测分析为基础，通过分析cdom浓度与表层海水中pahs(溶解态pahs)的定量关系，开发以cdom浓度对溶解态pahs进行定量表征的关键技术，根据cdom提取模型研发基于goci数据的海洋表层水体中pahs遥感监测方法，从而实现对河口海湾水体中溶解态pahs污染分布状况进行大范围、周期性的重复监测，为今后河口海湾pahs浓度分布的遥感监测提供预研成果。
技术实现要素：
本发明的目的是：针对现有传统技术的不足以及海洋生态环境监测和科学研究的需要，本发明提出的一种使用遥感手段进行河口海湾海洋污染物的监测方法，尤其适用于河口海湾表层水体中溶解态pahs的浓度监测，具有大面积、高频率的特征。本发明通过建立河口海湾水色参数cdom浓度定量表征海洋溶解态pahs的含量算法模型，并应用到世界上首颗地球静止气象卫星传感器goci影像上，获取河口海湾表层海水中溶解态pahs的遥感反演浓度大面积分布情况。
本发明解决其关键技术问题所采用的技术方案包括以下几个步骤：
步骤(1)、对待监测河口海湾区域提取水色参数cdom浓度的实测数据、表层海水中pahs浓度的实测数据；同时获取表层水体遥感反射率rrs的实测数据；
所述的cdom浓度无法直接测得，用cdom吸收系数表示，通过现场采集水样并过滤，采用国际上通用的“光透射测量方法”，使用紫外-可见分光光度计uv-3600测定cdom吸光度d(λ)，根据样品的吸光度，得到cdom的吸收系数acdom(λ),单位为m-1：
所述的pahs浓度采用液相色谱法测量，为现有成熟技术，即通过现场采集水样并过滤，将过滤后的水样按照《水质多环芳烃的测定液液萃取和固相萃取高效液相色谱法》(hj478-2009)完成萃取、净化、浓缩等预处理工作，然后采用液相色谱法测量pahs浓度。
所述的pahs浓度为usepa选定的多环芳烃中16种优先控制物质总浓度，16种多环芳烃分别为萘(nap)、苊烯(acy)、二氢苊(ace)、芴(flu)、菲(phe)、蒽(ant)、荧蒽(flua)、芘(pyr)、苯并[a]蒽(baa)、
所述的表层水体遥感反射率rrs是采用asd手持式光谱仪水面以上测量法获取，通过asd手持式光谱仪分别测量水体、天空、灰板的光谱曲线，遥感反射率由公式(2)计算所得；
其中，lsw为海水表面辐亮度，仪器测量水体所得；lsky为天空漫散射光，仪器测量天空获得；r为菲涅尔反射率，即气-水界面对天空光的反射率，受到太阳位置、观测几何、风速风向以及海面粗糙度等因素的影响，一般取值0.028；lp为标准灰板反射信号，由仪器对准灰板测量所得；ρp为标准板的反射率，由实验室校正所得。
步骤(2)、根据现有成熟的地球静止海洋水色成像仪goci，获取goci同步l1b遥感数据；使用envi软件中的flassh模块对gocil1b数据的进行flaash大气校正。
所述的goci遥感数据与cdom浓度、pahs浓度、rrs为同步数据，为官网上下载的l1b一级影像数据；
所述的goci同步遥感数据的大气校正，具体如下：
所述的使用envi软件中的flaash模块实现gocil1b数据的flaash大气校正，具体如下：
第一步，在envi中输入bil或者是bip格式的gocil1b数据；
第二步，打开flaash模块，输入goci传感器各个波段的中心波长和半峰全宽等参数；
第三步，继续输入大气参数，包括大气能见度和气溶胶类型；
第四步，根据goci的光谱响应函数制作光谱响应函数光谱库文件；
第五步，输入光谱响应函数光谱库文件，运行flaash模块，制作生成遥感反射率数据；
步骤(3)、根据实测pahs浓度含量和cdom浓度定量关系进行定量分析，构建cdom定量表征表层水体中pahs浓度的算法模型，根据goci波段设置以及经验，选用443nm处吸收系数表示cdom浓度：
pahs的浓度和cdom在浓度具有一定的正相关性(r2＝0.87，p<0.01)，随着acdom(443nm)的增加，pahs浓度呈指数增长趋势，海水中的cdom对pahs产生一定的增容作用，能够改变pahs在海洋环境中的分布情况。
步骤(4)、根据rrs的实测数据，以及cdom浓度的实测数据，建立cdom浓度估算模型，，见公式(4)：
443nm处cdom吸收系数acdom(443)＝a(443)-aw(443)-ap(443)(4)；①水体总吸收系数
其中悬浮颗粒后向散射系数
参数
水面以下遥感反射率rrs＝rrs/(0.52+1.7rrs)(8)
水面以上遥感反射率rrs通过公式(2)所得；
555nm处悬浮颗粒后向散射系数
其中aw(555)表示555nm纯水吸收系数；
参数
443、555nm处纯水后向散射系数bbw(443)、bbw(555)为固定值，根据ioccg的网站查询，分别取值为0.0024,0.0009。
443nm处参数
55nm处参数
②443nm处纯水吸收系数aw(443)，为固定值，根据ioccg的网站查询，取值为0.0072。
③443nm处浮游植物和非藻类颗粒物的吸收系数之和
ap(443)＝j1bbp(555)j2j1＝0.63，j2＝0.88(14)
步骤(5)、根据步骤(4)cdom浓度估算模型使用matlab实现cdom浓度估算，结合步骤(2)大气校正得到的rrs数据，从而获取goci的cdom遥感产品；
步骤(6)、根据步骤(3)cdom定量表征pahs浓度的算法模型，结合步骤(5)goci的cdom遥感产品，从而获得goci的pahs浓度分布情况。
本发明的有益效果如下：
本发明能够实现使用遥感手段监测近海水域pahs的污染状况，发明主要是为了大面积、高频次地监测海洋污染物的分布特征，提出基于实测数据分析的河口海湾表层水体中pahs的分布状况遥感监测模型，并应用到goci卫星数据，借助海洋水体中可遥组分cdom实现pahs浓度的遥感监测。
本发明有利于提高包括pahs在内的海洋环境污染物监测技术水平，通过技术推广、实际应用示范，将卫星遥感等空间信息技术应用服务更好地纳入到浙江设环境保护监测的决策链中。对于减少、消除和预防海洋持久性有机污染物带来的健康及环境风险，维系人类健康繁衍和维护生态环境安全，促进可持续发展，建立和谐社会具有重要理论意义和实际应用价值。该发明在海洋污染物、海洋生态监测等行业具有较高的推广价值和应用前景。
附图说明
图1为本发明流程图。
图2为2014年5月7日泛杭州湾海域cdom浓度卫星遥感反演结果；其中编号00，01，02，03，05，06与goci数据编号一致，分别为对应goci传感器成像时间8:30，9:30，10:30，11:30，13:30，14:30的影像结果。
图3为2014年5月7日泛杭州湾海域表层水体溶解态pahs浓度分布；其中编号00，01，02，03，05，06与goci数据编号一致，分别为对应goci传感器成像时间8:30，9:30，10:30，11:30，13:30，14:30的影像结果。
具体实施方式
下面结合附图对本发明作进一步的说明。
参照图1，一种海洋表层水体溶解态pahs的遥感监测方法，所述遥感监测pahs方法包括以下步骤：
步骤(1)、对待监测河口海湾区域提取水色参数cdom浓度的实测数据、表层海水中pahs浓度的实测数据；同时获取表层水体遥感反射率rrs的实测数据；
所述的cdom浓度无法直接测得，用cdom吸收系数表示，通过现场采集水样并过滤，采用国际上通用的“光透射测量方法”，使用紫外-可见分光光度计uv-3600测定cdom吸光度d(λ)，根据样品的吸光度，得到cdom的吸收系数acdom(λ)(单位m-1)：
所述的pahs浓度采用液相色谱法测量，为现有成熟技术，即通过现场采集水样并过滤，将过滤后的水样按照《水质多环芳烃的测定液液萃取和固相萃取高效液相色谱法》(hj478-2009)完成萃取、净化、浓缩等预处理工作，然后采用液相色谱法测量pahs浓度。
所述的pahs浓度为usepa选定的多环芳烃中16种优先控制物质总浓度，16种多环芳烃分别为萘(nap)、苊烯(acy)、二氢苊(ace)、芴(flu)、菲(phe)、蒽(ant)、荧蒽(flua)、芘(pyr)、苯并[a]蒽(baa)、
所述的表层水体遥感反射率rrs是采用asd手持式光谱仪水面以上测量法获取，通过asd手持式光谱仪分别测量水体、天空、灰板的光谱曲线，遥感反射率由公式(2)计算所得；
其中，lsw为海水表面辐亮度，仪器测量水体所得；lsky为天空漫散射光，仪器测量天空获得；r为菲涅尔反射率，即气-水界面对天空光的反射率，受到太阳位置、观测几何、风速风向以及海面粗糙度等因素的影响，取值0.028；lp为标准灰板反射信号，由仪器对准灰板测量所得；ρp为标准板的反射率，由实验室校正所得。
步骤(2)、根据现有成熟的地球静止海洋水色成像仪goci，获取goci同步l1b遥感数据；使用envi软件中的flassh模块对gocil1b数据的进行flaash大气校正。
所述的goci遥感数据与cdom浓度、pahs浓度、rrs为同步数据，为官网上下载的l1b一级影像数据；
所述的goci同步遥感数据的大气校正，具体如下：
所述的使用envi软件中的flaash模块实现gocil1b数据的flaash大气校正，具体如下：
第一步，在envi中输入bil或者是bip格式的gocil1b数据；
第二步，打开flaash模块，输入goci传感器各个波段的中心波长和半峰全宽等参数；
第三步，继续输入大气参数，包括大气能见度和气溶胶类型；
第四步，根据goci的光谱响应函数制作光谱响应函数光谱库文件；
第五步，输入光谱响应函数光谱库文件，运行flaash模块，制作生成遥感反射率数据；
步骤(3)、根据实测pahs浓度含量和cdom浓度定量关系进行定量分析，构建cdom定量表征表层水体中pahs浓度的算法模型，根据goci波段设置以及经验，选用443nm处吸收系数表示cdom浓度：
pahs的浓度和cdom在浓度具有一定的正相关性(r2＝0.87，p<0.01)，随着acdom(443nm)的增加，pahs浓度呈指数增长趋势，海水中的cdom对pahs产生一定的增容作用，能够改变pahs在海洋环境中的分布情况。
步骤(4)、根据rrs的实测数据，以及cdom浓度的实测数据，建立cdom浓度估算模型，，见公式(4)：
443nm处cdom吸收系数acdom(443)＝a(443)-aw(443)-ap(443)(4)；②水体总吸收系数
其中悬浮颗粒后向散射系数
参数
水面以下遥感反射率rrs＝rrs/(0.52+1.7rrs)(8)
水面以上遥感反射率rrs通过公式(2)所得；
555nm处悬浮颗粒后向散射系数
555nm处水体总吸收系数a(555)，公式如下：
其中aw(555)表示555nm纯水吸收系数；
参数
443、555nm处纯水后向散射系数bbw(443)、bbw(555)为固定值，根据ioccg的网站查询，分别取值为0.0024,0.0009。
443nm处参数
55nm处参数
②443nm处纯水吸收系数aw(443)，为固定值，根据ioccg的网站查询，取值为0.0072。
③443nm处浮游植物和非藻类颗粒物的吸收系数之和
ap(443)＝j1bbp(555)j2j1＝0.63，j2＝0.88(14)
步骤(5)、根据步骤(4)cdom浓度估算模型使用matlab实现cdom浓度估算，结合步骤(2)大气校正得到的rrs数据，从而获取goci的cdom遥感产品；
步骤(6)、根据步骤(3)cdom定量表征pahs浓度的算法模型，结合步骤(5)goci的cdom遥感产品，从而获得goci的pahs浓度分布情况见图3。泛杭州湾海域表层水体溶解态pahs遥感监测结果与实测pahs基本一致，反演效果良好，实现了泛杭州湾海域pahs大面积、高频率的遥感监测。空间分布上，长江口以及舟山东南部海域表层海水中溶解态pahs的含量明显高于其他海域，一方面来自由江河径流携带而来的近岸地表径流的陆源输入。对于远离陆地的外海区，受海水的混合稀释作用影响，江河携带的有机物质大大减少，随着cdom浓度的减少pahs的含量总体偏低。日变化上，影像时间8:30，在长江口以及舟山东南部海域pahs呈较高浓度分布，随着时间推移，较高的pahs浓度呈现南移趋势，至10:30，长江口附近pahs浓度减小，舟山东南部海域浓度明显增加。较高的pahs浓度继续南移，至下午泛杭州湾海域pahs浓度含量均减小。
- 还没有人留言评论。精彩留言会获得点赞！