楊俊杰， 黃妙芬， 駱蔚健， 王忠林， 邢旭峰

1. 廣東海洋大學(xué)電子與信息工程學(xué)院， 廣東 湛江 524088

2. 廣東海洋大學(xué)數(shù)學(xué)與計算機(jī)學(xué)院， 廣東 湛江 524088

3. 廣東海洋大學(xué)海洋與氣象學(xué)院， 廣東 湛江 524088

引 言

目前， 關(guān)于海洋油污染的研究主要集中在油膜的遙感探測和識別上， 主要的技術(shù)手段包括熱紅外[1]、 雷達(dá)[2]、 光學(xué)遙感[3]， 而關(guān)于水中油對太陽光在水體輻射傳輸過程中的影響研究甚少。 油屬于光學(xué)活性物質(zhì)[4]， 水中油能夠改變太陽光在水體中的輻射傳輸過程。 這不僅會引起上行輻射信息的變化， 也會影響到水上傳感器的入瞳信號強(qiáng)度， 還會造成水體光場的輻亮度信息變化。 光場指的是光譜輻亮度在水體中各個位置隨角度分布的完整描述， 水下光場的研究是解決水下能見度、 水下觀測系統(tǒng)、 水下成像、 水下目標(biāo)的目視檢測及水下初級生產(chǎn)力等問題的前提條件[5-6]。 輻亮度是描述光的一種輻射測量學(xué)術(shù)用語， 也是海洋光學(xué)關(guān)注的基本輻射物理量之一， 所有其他的輻射物理量(如輻照度、 漫衰減系數(shù)、 光合有效輻射等)都可以通過輻亮度推導(dǎo)得出[7]。 這些參數(shù)可以反映不同水體可獲得太陽能量的大小或光子數(shù)量的多少， 進(jìn)而影響水下生物的生長過程及海氣界面的碳交換過程。 輻亮度可用于全面描述水體光場結(jié)構(gòu)隨空間、 時間、 方向及波長的變化[7-8]。 水體的水下光場及輻亮度分布則是由海水的固有光學(xué)特性、 其垂向結(jié)構(gòu)及光照條件所決定的[5]。 輻亮度與水體組分固有光學(xué)特性之間的關(guān)系可由輻射傳輸方程進(jìn)行表達(dá)和計算[9]。 油污染事件發(fā)生后， 在海洋物理過程[10]的作用下或在人為污染應(yīng)急處理[11]之后， 部分油類污染物質(zhì)與海水混合形成油水乳化物， 以水中油的形式存在于海水水體中， 水中油的出現(xiàn)及其在海水中的垂向分布會改變水體的光學(xué)特性。 因此， 含油水體的水下光場研究是描述太陽輻射在含油水體中輻射傳輸過程的基礎(chǔ)工作。 其中， 輻亮度隨海水深度、 波長及太陽入射天頂角的變化特征是研究含油水體水下光場及進(jìn)一步開展含油水體對海洋生地化循環(huán)過程影響的核心內(nèi)容； 同時， 上行輻亮度還是水色傳感器的重要信號源， 也是利用水色遙感技術(shù)探測水體組分構(gòu)成的關(guān)鍵物理量。

水下光場及水體中輻亮度隨空間、 時間、 波長和方向變化的研究方法主要有兩種， 分別是理論模擬和野外實測兩種方法。 同野外實測方法相比， 理論模擬方法具有能夠描述更多細(xì)節(jié)特征、 計算結(jié)果精確度更高、 代價更低等優(yōu)點， 并且， 理論方法還能用于模擬計算野外實測無法開展測量的情況。 其中， 比較常見的理論模擬方法有Hydrolight和Monte Carlo。 和Monte Carlo方法相比， Hydrolight采用不變嵌入法[7, 12]計算水體輻射傳輸方程， 具有計算速度快、 精度高等優(yōu)點而被廣泛應(yīng)用于海洋光學(xué)及水色遙感領(lǐng)域。

本文利用在大連港海域現(xiàn)場采集的石油類物質(zhì)固有光學(xué)特性及水中油濃度數(shù)據(jù)， 并基于Hydrolight進(jìn)行含水中油水體水下光場數(shù)值模擬， 旨在揭示水中油的敏感光譜特性， 并分析水中油對純海水水體上行輻亮度光譜的影響， 從水色遙感機(jī)理角度分析探測水中油污染的可行性， 為進(jìn)一步利用遙感技術(shù)探測水中油污染信息奠定基礎(chǔ)。

1 實驗部分

1.1 研究區(qū)域及觀測站位

研究區(qū)位于中國遼寧大連港海域， 如圖1所示。 大連港坐落在遼東半島的南端， 屬于港闊水深、 冬季不凍的萬噸貨輪港， 其中新港港區(qū)、 北良港區(qū)和長興島港區(qū)是主要的油品碼頭。 原油卸裝及中轉(zhuǎn)過程發(fā)生的原油泄漏、 大型船舶運輸過程的排放及石油污染事故均對大連港海域的水環(huán)境產(chǎn)生危害。 圖1中的A， B和C為本次研究的觀測站位， A站位為歷史石油管道爆炸事故點， B站位位于大連港航道上， C站位遠(yuǎn)離事故發(fā)生地點和船運航道。

圖1 研究區(qū)位及觀測點位

1.2 數(shù)據(jù)測量與光場模擬

2018年8月25日—27日， 分別依每個整點時刻從7:00—17:00在圖1的三個站點進(jìn)行測量， 測量的參數(shù)包括水中油濃度、 吸收系數(shù)、 后向散射系數(shù)等。 水中油濃度、 吸收系數(shù)及后向散射系數(shù)的測量方法和分離算法已被詳細(xì)闡述[13-15]， 不再贅述。

本文將現(xiàn)場測量及處理后的數(shù)據(jù)作為Hydrolight的輸入， 模擬水深為50 m(間隔為0.6 m)的含水中油(濃度分別為0.1， 0.2， 0.3， 0.5， 1.0， 1.5， 2.0， 3.0， 5.0和10.0 mg·m-3)水體、 在不同太陽天頂角(0°， 15°， 30°， 45°， 60°， 75°和90°)條件下的水下光場， 波長區(qū)間為300～1 000 nm(分辨率為5 nm)。 文獻(xiàn)[15]中列舉了更詳細(xì)的單一組分含水中油水體光場模擬的Hydrolight參數(shù)設(shè)置[15]。

2 結(jié)果與討論

2.1 水中油上行輻亮度光譜特性

圖2為含水中油水體不同水深(剛好在水面之下0-， 10.2， 20.4， 30.0， 40.2和50.0 m)處的上行輻亮度光譜(300～1 000 nm)， 圖2(a—l)分別表示含不同濃度(0.1， 1， 5和10 mg·m-3)的水中油水體在不同天頂角(0°， 45°和90°)入射條件下的上行輻亮度光譜。 首先， 在300～1 000 nm光譜區(qū)間， 不同濃度的水中油對水體上行輻亮度光譜的影響主要集中在可見光波段(380～760 nm)， 而在紫外和紅外波段， 由于純水的高吸收特性， 水中油對水體的上行輻亮度光譜影響微弱， 甚至可以忽略(尤其當(dāng)水中油濃度較低時)。 其次， 當(dāng)水中油濃度較低時， 含水中油水體主要體現(xiàn)出純海水的上行輻亮度光譜特性[圖2(a—c)]， 峰值出現(xiàn)在417.5 nm附近。 其主要原因在于純水的吸收系數(shù)在可見光波段大致呈遞增的趨勢， 藍(lán)綠波段純水的吸收能力弱， 在紅光波段純水的吸收系數(shù)快速增加。 隨著水中油濃度的增加， 上行輻亮度光譜峰值逐漸向長波方向(442.5～472.5 nm)移動， 且在藍(lán)光波段范圍內(nèi)， 水中油上行輻亮度光譜均出現(xiàn)一個明顯的吸收谷(432.5 nm)。 此外， 隨著水中油濃度的增加， 水中油的吸收效應(yīng)逐漸凸顯， 上行輻亮度光譜在藍(lán)光波段的量值逐漸降低。

水色遙感傳感器的光譜分辨率高， 光譜范圍一般覆蓋400～900 nm區(qū)間。 圖2表明不同濃度的水中油均可引起水體上行輻亮度光譜(可見光波段)的明顯變化， 其波譜響應(yīng)區(qū)間與水色傳感器的波譜設(shè)置相對應(yīng)， 屬于水色遙感的研究光譜范疇， 為利用水色遙感技術(shù)探測含水中油水體及反演水中油濃度提供了重要的光譜依據(jù)。 同時， 利用水色傳感器獲取含水中油水體的上行輻射的光譜信息， 可以為建立反演水中油濃度的反演模型提供大量的數(shù)據(jù)源。

圖2 上行輻亮度光譜

2.2 水中油上行輻亮度的深度響應(yīng)機(jī)制

水體組分對不同波長的太陽輻射具有不同的吸收和散射能力， 因此， 不同波長的太陽輻射具有相異的穿透能力。 由于石油污染事件多發(fā)生在近海海域， 本文設(shè)置水體的最大水深為50 m， 分析三個水色遙感常用波段(442.5， 557.5和677.5 nm)的上行輻亮度隨水體深度(剛好在水面以下0-～50 m)的變化特征(圖3)。 圖3(a—f)分別表示兩種濃度(1和10 mg·m-3)水中油水體在三種天頂角(0°， 45°和90°)入射條件下的上行輻亮度隨深度的變化曲線。 上行輻亮度隨著水深(剛好在水面以下0-～50 m)的增加呈現(xiàn)出先遞減、 而后在接近水體下界面之前遞增的趨勢。 此外， 含高濃度(10 mg·m-3)水中油水體的上行輻亮度隨水深的增加(呈遞增趨勢之前的深度)而快速減少， 而含低濃度水中油水體的上行輻亮度隨水深增加而降低的速率則相對較慢(尤其在藍(lán)綠波段)， 該現(xiàn)象表明含高濃度水中油水體對光的衰減作用強(qiáng)于含低濃度水中油水體。

圖3 不同波長(442.5, 557.5和677.5 nm)的上行輻亮度隨水體深度的變化曲線

采用最大水深50 m的下界面， 下界面將部分入射至此處的太陽輻射反射回水體， 這是造成越接近水體的下界面， 上行輻亮度不減反增的根本原因。 本文含水中油水體屬于近岸或光學(xué)淺水水體(optically shallow waters)， 含不同濃度的水中油水體剛好在水面之下(0-)上行輻亮度之間的差異(圖2和圖3)表明， 水體中不同濃度的水中油對剛好在水面之下(0-)上行輻亮度的貢獻(xiàn)是不同的； 同時， 也不可忽視水體下界面對剛好在水面之下(0-)上行輻亮度的貢獻(xiàn)(圖3)。 水體輻射傳輸利用水體中不同深度內(nèi)的水體組分與光的相互作用， 對光的下行和上行傳輸過程進(jìn)行描述， 是水色遙感的核心機(jī)理。 不同濃度的水中油對光線在水體中輻射傳輸過程的影響與水色三要素一致， 均會引起光在傳輸過程中不同程度的衰減， 符合利用水色遙感技術(shù)探測水體組分的輻射傳輸機(jī)理。 這表明可以將水中油視作一種新的水體組分， 基于水體輻射傳輸機(jī)理開展含水中油水體的水色遙感反演研究。

2.3 剛好位于水面之上和水面之下上行輻亮度的太陽天頂角響應(yīng)機(jī)理

水體上行輻亮度需經(jīng)過水汽界面的透射后再經(jīng)大氣傳輸進(jìn)入傳感器， 傳感器除了能接收到水體的上行輻射信息之外， 還受到水汽界面對太陽輻射反射的影響。 圖4是剛好在水面之上(帶點實線)和剛好在水面之下(虛線)的上行輻亮度隨天頂角(0°， 15°, 30°， 45°， 60°， 75°及90°)的變化曲線圖。 隨著天頂角的增大， 輻亮度光譜的量值逐漸變小， 表明太陽天頂角對輻亮度光譜的量值具有影響作用。 此外， 當(dāng)天頂角為0°時， 剛好在水面之上(0+)上行輻亮度往往超過剛好在水面之下(0-)上行輻亮度的1～3個數(shù)量級。

圖4 剛好位于水面之上和水面之下的上行輻亮度隨天頂角變化曲線圖

海表面對太陽直射輻射的反射會淹沒水體組分的弱信號， 在含水中油水體中也具有該效果(圖4和表1)。 表1以圖4(d)(10 mg·m-3)為例， 計算在不同波長及天頂角條件下剛好在水面之上(0+)上行輻亮度占剛好在水面之下(0-)上行輻亮度的比率。 當(dāng)太陽天頂角為0°時， 太陽直射入射反射參與貢獻(xiàn)的剛好在水面之上(0+)上行輻亮度是剛好在水面之下(0-)上行輻亮度的幾十甚至幾百倍。 太陽直射入射反射信號屬于強(qiáng)信號， 不包含水體組分的組成信息。 水色傳感器具有一定的側(cè)擺能力， 可以避免太陽直射反射信號進(jìn)入傳感器， 保證傳感器接收到的有用信號主要由水體組分(本文中的水分子和水中油)貢獻(xiàn)。

水體向上輻射傳輸?shù)男盘柲芰枯^低， 通常清潔水體在影像上呈現(xiàn)暗像元的特征， 因而， 水色傳感器均具有較高的信噪比特點。 如2.1所述， 藍(lán)光波段含最強(qiáng)的水中油光譜信號，然而， 隨太陽天頂角及水中油濃度的增加， 藍(lán)光波段上行輻亮度量值逐漸降低， 勢必造成傳感器入瞳信號減弱。 含水中油水體上行輻亮度量值隨水中油濃度及太陽天頂角增加而減少的現(xiàn)象表明水色衛(wèi)星上午10點至下午14點的過境特點能保證有足夠多含水中油水體的上行輻亮度進(jìn)入水色傳感器， 且水色傳感器的高信噪比特征也適合屬于暗目標(biāo)體的含水中油水體探測。 所以， 水色衛(wèi)星提供的遙感數(shù)據(jù)將是探測水中油污染和反演水中油濃度的重要數(shù)據(jù)源之一。

表1 剛好在水面之上(0+)上行輻亮度占剛好在水面之下(0-)上行輻亮度的比率(%)

3 結(jié) 論

結(jié)合大連港現(xiàn)場實測數(shù)據(jù)及輻射傳輸模型Hydrolight對含水中油水體的上行輻亮度進(jìn)行模擬， 分析水中油對上行輻亮度光譜的影響及上行輻亮度隨深度及天頂角的變化特征， 結(jié)果表明：

(1)含水中油水體的上行輻亮度光譜響應(yīng)區(qū)間與水色傳感器的波譜設(shè)置相對應(yīng)， 主要在可見光波段， 且其峰值變化集中在藍(lán)光波段(442.5～472.5 nm)， 為利用水色遙感技術(shù)探測含水中油水體及反演水中油濃度提供了重要的光譜依據(jù)， 以便發(fā)揮水色遙感技術(shù)的范圍廣、 快速、 成本低且便于長期動態(tài)監(jiān)測的優(yōu)點。

(2)與海表油膜通過油類污染物質(zhì)改變海表反射率的遙感探測機(jī)理不同， 含水中油水體上行輻亮度隨深度變化特征表明各個水深處的含水中油水層及下界面均對上行輻射產(chǎn)生不同程度的貢獻(xiàn)。 在水媒介中， 水中油對光在水體中輻射傳輸過程的影響機(jī)理與水色三要素是一致的， 不僅會在下行傳輸過程中造成不同程度的衰減， 還會對上行輻射產(chǎn)生影響。 因此， 開展水色遙感探測和反演水中油具備了機(jī)理基礎(chǔ)， 可將水中油視作一種新的水體組分， 基于輻射傳輸機(jī)理研究水中油對太陽輻射上行及下行傳輸過程的影響。

(3)水色衛(wèi)星搭載的傳感器具有高信噪比的特點和一定的側(cè)擺能力， 且其過境時間一般在上午10點和下午14點之間， 使得水色傳感器既能保證避開太陽直射反射信號以接收到由水體組分貢獻(xiàn)的有用信號， 還能最大程度感應(yīng)到由水體組分變化引起的暗像元變化。 因此， 水色遙感數(shù)據(jù)是建立反演水中油濃度算法的可靠數(shù)據(jù)源之一。