上海航天空間技術(shù)有限公司

上海水環(huán)境模擬與水生態(tài)修復(fù)工程技術(shù)研究中心 鄧成晨 鄭順安 杭 君 王 婧

淀山湖又稱薛淀湖，屬吞吐性淺水湖泊，位于上海市青浦區(qū)與江蘇省蘇州市昆山市交界處，緯度范圍為31°04′N 至31°12′N，經(jīng)度范圍為120°53′E 至121°01′E，總面積約62km2。其中，屬于上海境內(nèi)的淀山湖湖區(qū)面積約46.7km2（約占75.3%），是上海最大的淡水湖泊。淀山湖平均水深2.1m，最大水深3.6m，湖泊容積1.3×108m3，年均水溫20.6℃，主要承泄太湖來水。太湖水由西北向東南經(jīng)急水港、大朱厙等河港進(jìn)入湖體，然后經(jīng)攔路港、淀浦河等河流泄入黃浦江，停留時(shí)間約29 天，約占黃浦江水量的17%。淀山湖是黃浦江上游重要的水源保護(hù)地，其水質(zhì)狀況對上海市飲用水源地的飲用水安全起著重要的作用。

近年來，由于環(huán)湖農(nóng)藥化肥及鄉(xiāng)鎮(zhèn)工業(yè)廢水等排入，湖水局部呈現(xiàn)富營養(yǎng)化現(xiàn)象。2017年急水港橋的數(shù)據(jù)表明，淀山湖水質(zhì)處于劣V 類水平，主要污染物是總氮?？偟怯绊懺孱惙敝车闹饕獱I養(yǎng)鹽之一，總氮含量上升，導(dǎo)致湖泊富營養(yǎng)化，進(jìn)而使藻類等浮游生物過度繁殖，使得近些年淀山湖藍(lán)藻水華暴發(fā)現(xiàn)象頻發(fā)，水質(zhì)惡化，對居民生活造成嚴(yán)重影響。因此，快速、全面掌握藍(lán)藻水華暴發(fā)情況，對于藍(lán)藻水華治理和預(yù)警有著非常重要的作用。

隨著衛(wèi)星遙感技術(shù)的快速發(fā)展，通過衛(wèi)星遙感技術(shù)實(shí)現(xiàn)對藍(lán)藻水華的大范圍監(jiān)測成為可能，而衛(wèi)星遙感技術(shù)具有宏觀性、綜合性和實(shí)效性等優(yōu)勢，已被廣泛應(yīng)用于水體藍(lán)藻水華監(jiān)測中。

一、數(shù)據(jù)來源

據(jù)上海水務(wù)局報(bào)道稱，2018年7月14日中午，螢火蟲環(huán)境保育志愿者小組在金澤鎮(zhèn)巡視期間路過環(huán)淀山湖生態(tài)帶（蔡浜村段）湖濱時(shí)，記錄到近岸水域長約1.7km，寬約5～20m 的明顯藍(lán)藻水華污染帶（如圖1所示）。這與青浦區(qū)水務(wù)局2017年8月15日印發(fā)的《青浦區(qū)淀山湖藍(lán)藻水華應(yīng)急處置預(yù)案》（以下簡稱“預(yù)案”）中提到的淀山湖水華暴發(fā)時(shí)間和區(qū)域一致。預(yù)案指出，“淀山湖水華的發(fā)生期集中在6月?10月，其中暴發(fā)期集中在7月?9月。藍(lán)藻受風(fēng)向作用向下風(fēng)口移動而堆積。根據(jù)淀山湖岸線地形和夏季風(fēng)向特點(diǎn)，一般在東南風(fēng)作用下，藍(lán)藻會在金澤鎮(zhèn)蔡浜村淀山湖防洪大堤處形成堆積。”且該污染帶位于上海市引用水源二級保護(hù)區(qū)內(nèi)。

圖1 志愿者拍攝的環(huán)淀山湖生態(tài)帶（蔡浜村段）湖濱藍(lán)藻水華暴發(fā)情況

通過查詢事發(fā)時(shí)段的高分衛(wèi)星數(shù)據(jù)發(fā)現(xiàn)，2018年7月15日，“高分”五號衛(wèi)星（GF-5）的可見光短波紅外高光譜相機(jī)（AHSI）有效載荷對上海市西部進(jìn)行了拍攝，影像覆蓋整個(gè)淀山湖區(qū)域范圍（如圖2所示）；7月16日，“高分”一號衛(wèi)星（GF-1）的寬幅相機(jī)（WFV4）對淀山湖區(qū)域進(jìn)行拍攝；7月19日、20日、21日分別有“高分”衛(wèi)星對淀山湖區(qū)域進(jìn)行了拍攝，但因天氣問題，大部分湖面被云層覆蓋，數(shù)據(jù)質(zhì)量欠佳，故筆者選取7月15日的GF-5 衛(wèi)星數(shù)據(jù)和7月16日的GF-1 衛(wèi)星數(shù)據(jù)進(jìn)行淀山湖區(qū)域藍(lán)藻水華反演。

圖2 GF-5衛(wèi)星2018年7月15日拍攝淀山湖及志愿者發(fā)現(xiàn)藍(lán)藻水華位置情況

GF-1 衛(wèi)星是我國第一顆自主研制的對地觀測高分辨率衛(wèi)星，于2013年4月26日發(fā)射升空，是中國“高分”系列的第一顆衛(wèi)星，該星配有2 臺2m 全色/8m 多光譜相機(jī)和4 臺16m 多光譜寬幅相機(jī)（WFV），設(shè)計(jì)壽命為5～8年。寬幅相機(jī)幅寬為800km，可實(shí)現(xiàn)大視野觀測。

GF-5 衛(wèi)星是我國重要的民用高光譜觀測衛(wèi)星，于2018年5月9日發(fā)射升空并于2019年4月10日因故失效，“高分”五號02 星已于2021年9月7日發(fā)射升空。GF-5 衛(wèi)星搭載的可見光短波紅外高光譜相機(jī)是世界上首臺兼顧寬覆蓋與寬譜段的高光譜相機(jī)，相機(jī)空間分辨率為30m，幅寬為60km，波譜范圍從可見光到短波紅外（400～2500nm），波譜通道數(shù)為330 個(gè)，可見光譜段的光譜分辨率為5nm，可對內(nèi)陸水體、生態(tài)環(huán)境、巖礦種類等進(jìn)行有效探測，為環(huán)境監(jiān)測、農(nóng)業(yè)資源、防災(zāi)減災(zāi)等行業(yè)提供高質(zhì)量、高可靠性的高光譜數(shù)據(jù)。

二、數(shù)據(jù)預(yù)處理和分析

湖泊水體中藍(lán)藻水華的生長受多種環(huán)境因子的影響和制約，而葉綠素a 的濃度則在很大程度上反應(yīng)藍(lán)藻水華的生長情況。因此，筆者基于GF-1/WFV4 和GF-5/AHSI 衛(wèi)星載荷數(shù)據(jù)，通過反演葉綠素a 濃度來反映藍(lán)藻水華生長及空間分布情況，采用的反演研究技術(shù)路線如圖3所示。

圖3 葉綠素a反演研究技術(shù)路線

首先對獲取的兩景高分影像進(jìn)行數(shù)據(jù)預(yù)處理，主要包括輻射定標(biāo)、大氣校正和幾何校正等。

輻射定標(biāo)是將遙感傳感器的數(shù)字量化值（DN）轉(zhuǎn)換為大氣外層表面發(fā)射率（或稱為表觀輻射亮度值）以消除傳感器誤差的過程，是定量化的基礎(chǔ)，其計(jì)算公式如式（1）所示。

式中，DNλ是指不同波段中地物的灰度值，Gain和Offset分別為定標(biāo)斜率（Wm-2sr-1mm-1）和絕對定標(biāo)系數(shù)偏移量（Wm-2sr-1μm-1），Lλ為表觀輻射亮度。

大氣校正的目的是消除由吸收和散射作用造成的輻射量誤差，反演出地物真實(shí)的反射率。筆者使用了2 種大氣校正方法，GF-1/WFV4 傳感器數(shù)據(jù)使用了基于輻射傳輸模型的FLAASH 大氣校正；GF-5/AHSI 傳感器數(shù)據(jù)利用了6S 輻射傳輸模型根據(jù)輸入的參數(shù)動態(tài)構(gòu)建大氣校正查找表的方式來進(jìn)行校正。

使用衛(wèi)星遙感技術(shù)監(jiān)測水體主要是利用水的光譜反射輻射特征（如反射率）與各種水體參數(shù)（如葉綠素a 濃度）的關(guān)系，建立一系列相關(guān)模型來提取或反演水體參數(shù)。

葉綠素a 對可見光和近紅外光波段具有特殊的“陡坡效應(yīng)”。在420～500nm 的波段范圍，葉綠素a 有一個(gè)比較平坦的吸收峰，水體的光譜反射率較低，在540～580nm 的波段范圍存在一個(gè)反射峰，在670nm 波長附近存在又一個(gè)吸收峰，當(dāng)葉綠素a 濃度較高時(shí)，水體的光譜反射率在該處出現(xiàn)谷值，在700nm 波長附近存在又一個(gè)反射峰，該反射峰是含藻類水體最顯著特征，用于判定水體中是否含有藻類葉綠素a。

目前，多光譜衛(wèi)星數(shù)據(jù)反演葉綠素a 濃度主要有經(jīng)驗(yàn)法、半經(jīng)驗(yàn)/半分析法以及分析法。經(jīng)驗(yàn)法主要以葉綠素a 濃度和遙感參數(shù)之間的統(tǒng)計(jì)關(guān)系為基礎(chǔ)來實(shí)現(xiàn)對水體葉綠素a 濃度的遙感反演。主要包括單波段模型、波段比值模型、多波段組合模型等。半經(jīng)驗(yàn)/半分析法是結(jié)合經(jīng)驗(yàn)法和輻射傳輸模型，采用理論分析與經(jīng)驗(yàn)統(tǒng)計(jì)相結(jié)合的方法描述模型的過程。主要包括三波段模型、四波段模型等。分析法是從水色遙感機(jī)理出發(fā)，構(gòu)建水體光譜模型，用以表示水體參數(shù)與光譜反射率之間的相應(yīng)關(guān)系，通過輻射傳輸模型模擬光在大氣中的傳輸，建立水體固有光學(xué)量和表觀光學(xué)量的關(guān)系方程，計(jì)算得到水體參數(shù)。

單波段模型法是反演葉綠素a 濃度最簡單的經(jīng)驗(yàn)?zāi)Ｐ退惴?，通過分析水體及葉綠素a 的光譜信息，尋找葉綠素a 反射的敏感波段。但由于在單波段經(jīng)驗(yàn)?zāi)Ｐ椭校髯V段的光譜反射率與葉綠素a 濃度相關(guān)性都不大，因此，難以構(gòu)建理想的葉綠素a 濃度單波段反演模型。

波段比值模型法是將水體參數(shù)和光譜反射率在兩個(gè)波段上的比值之間建立相應(yīng)關(guān)系。對于葉綠素a 濃度而言，近紅外波段與紅波段比值與之具有一定的相關(guān)性，其表達(dá)方式如式（2）所示。

式中，λ1和λ2分別為高分衛(wèi)星遙感影像的近紅外波段和紅波段，波段比值法在一類水體中進(jìn)行葉綠素a 濃度反演的精度尚可。

針對葉綠素a 濃度對比了多種單波段和波段比值算法，對相關(guān)波段進(jìn)行組合分析，得出高分衛(wèi)星數(shù)據(jù)紅波段反射率與綠波段和近紅外反射率之和的比值與淀山湖水體葉綠素a 濃度具有較好的相關(guān)性，構(gòu)建基于高分衛(wèi)星數(shù)據(jù)的淀山湖葉綠素a 濃度反演模型如式（3）所示。

式中，ρChla 為葉綠素a 濃度，b4 為高分衛(wèi)星的近紅外波段；b3 為高分衛(wèi)星的紅波段；b2 為高分衛(wèi)星的綠波段；a、b、c 分別為相關(guān)系數(shù)。經(jīng)過模型修正后的反演結(jié)果如圖4所示。

圖4 2018年7月15日GF-5和16日GF-1數(shù)據(jù)反演葉綠素a濃度成果圖

分析結(jié)果顯示，在志愿者巡視的金澤鎮(zhèn)蔡浜村淀山湖防洪大堤處附近的水域葉綠素a 濃度約為30mg/m3，與志愿者拍攝照片的藍(lán)藻水華情況基本吻合。

三、結(jié)論

湖泊富營養(yǎng)化是現(xiàn)階段我國主要的水環(huán)境問題，湖泊藍(lán)藻水華的發(fā)生規(guī)律、內(nèi)在機(jī)制以及預(yù)警研究已經(jīng)成為重要的研究課題。近年來，淀山湖藍(lán)藻水華頻發(fā)，開展相關(guān)的基礎(chǔ)研究，可以為富營養(yǎng)化防治及藍(lán)藻水華預(yù)警提供科學(xué)數(shù)據(jù)及決策依據(jù)。

筆者基于GF-5 高光譜數(shù)據(jù)和GF-1 的寬幅相機(jī)多光譜數(shù)據(jù)的葉綠素a 濃度的建模分析，得到了2018年7月15日和16日淀山湖葉綠素a 濃度分布，從建模得到的結(jié)果來看，2 天葉綠素a 的空間分布及濃度變化不大，同時(shí)印證了志愿者的巡視結(jié)果。

今后可從以下4 個(gè)方面深入研究。一是本文所使用的算法模型尚不成熟，特別是在高濃度葉綠素a 水體中，精度誤差較大，需要多方面的優(yōu)化和提高。二是高分衛(wèi)星因各種原因，目前尚無法做到對同一地域進(jìn)行高頻次監(jiān)測，同時(shí)，上海屬亞熱帶季風(fēng)氣候，在藍(lán)藻水華暴發(fā)的5月?9月，上海陰雨天較多，若要進(jìn)行長時(shí)間跨度密集監(jiān)測，研究淀山湖藍(lán)藻水華暴發(fā)整個(gè)流程，最終對防治和預(yù)警提供決策依據(jù)，需要獲取多源衛(wèi)星和載荷數(shù)據(jù)以及更多的樣本。三是藍(lán)藻水華在不同的氣候環(huán)境下，如光照、氮、溶解氧、水溫等，藍(lán)藻水華的上浮率不同，因此，在葉綠素a 濃度建模反演時(shí)，需要考慮以上因素。四是衛(wèi)星遙感反演水體中的葉綠素a 濃度主要是通過獲取含葉綠素a 水體的反射輻射特性進(jìn)行，因此，應(yīng)該考慮可能影響葉綠素a 濃度反演的因素，主要包括水體中的其他綠色植物、氣象環(huán)境、水面油污及其他漂浮物、總懸浮物、水體透明度、濁度，甚至浪花等。