邢旭峰，王忠林，黃妙芬，黃穎恩，莊 煬

石油污染水體油濃度及后向散射系數(shù)垂向變化特性分析

邢旭峰，王忠林，黃妙芬，黃穎恩，莊 煬

（廣東海洋大學(xué)數(shù)學(xué)與計算機學(xué)院，廣東 湛江 524088）

【目的】探討石油污染水體油濃度及后向散射系數(shù)垂向分布特性，為研究石油類污染水體水下光場變化奠定基礎(chǔ)。【方法】于2018年8月25－27日在大連港采用站點固定、時間隨潮汐變化的觀測方式，測量不同深度的油物質(zhì)和懸浮物濃度，以及含油水體的后向散射系數(shù)（backscattering coefficient,b）?！窘Y(jié)果】油濃度隨深度的變化比較復(fù)雜，未呈現(xiàn)明顯的逐漸下降或者逐漸升高的趨勢；不同波段的b隨著深度的變化趨勢是一致的；在量級上，b值最小對應(yīng)的波段是700 nm，其次是590 nm波段，最大是510、470、442和420 nm，且這4個波段的b值差別不大；水體后向散射系數(shù)b隨深度的變化可以分為陡增型、平緩增加型和單峰型等三種類型?！窘Y(jié)論】在含油水體中油濃度和后向散射系數(shù)的垂向變化受到潮汐的影響，處于不同漲落潮時不同來向的海水其組分含量不同；石油物質(zhì)通過吸附在懸浮顆粒物表面形成一個雙層結(jié)構(gòu)，共同影響著水體的后向散射系數(shù)。

石油類物質(zhì)濃度；懸浮物；后向散射系數(shù)；垂向變化；大連港

在水色遙感研究領(lǐng)域，掌握石油污染水體油濃度及后向散射系數(shù)（backscattering coefficient,b）垂向分布特性，是進行石油水體輻射傳輸特性研究的重要前提，同時考慮隨波長和深度變化的吸收和散射特性是利用輻射傳輸方程來研究石油類污染水體光學(xué)特性的關(guān)鍵環(huán)節(jié)[1]。

國內(nèi)外眾多學(xué)者已經(jīng)開展了針對水體不同組分的垂向特性研究。龐重光等[2]以1959至2006年共24個航次近5000組懸浮物濃度數(shù)據(jù)為基礎(chǔ)，對黃、東海海域的懸浮物濃度垂直分布特征進行研究并建立回歸模型；何劍波等[3]對千島湖水體葉綠素a濃度垂向分布特征進行研究，發(fā)現(xiàn)在春夏秋季，葉綠素a濃度在垂向上呈現(xiàn)先增加后降低的特征；Xiu等[4]研究了I類水體葉綠素濃度垂向非均勻性對遙感反射比的影響，得出了當表層葉綠素濃度相同時，在葉綠素濃度垂向均勻和非均勻的情況下，水體中的光學(xué)特性具有明顯差異的結(jié)論；Stramska等[5]利用輻射傳輸數(shù)值模式HydroLight研究了與葉綠素濃度垂向分布相關(guān)的水柱的固有光學(xué)特性；梁其椿等[6]利用HydroLight模擬生成了3200組藻類垂向高斯分布模型的數(shù)據(jù)，構(gòu)建了基于BP神經(jīng)網(wǎng)絡(luò)（Back-propagation Network）來估算巢湖藻類垂向分布結(jié)構(gòu)的算法；馬孟梟等[7]利用Hydrolight輻射傳輸模式研究了水體不同組分垂向分布對水體漫衰減系數(shù)d的影響；徐文龍等[8]利用海洋光學(xué)觀測的高垂向分辨率剖面數(shù)據(jù)，分析了葉綠素a濃度垂向剖面的時空分布特征及其與海洋動力環(huán)境要素的關(guān)系；Yin等[9]研究了城市富營養(yǎng)化水體含氮組分的垂向分布特征；Xing等[10]利用浮標數(shù)據(jù)對南中國海的顆粒物和溶解性光學(xué)性質(zhì)的垂向變化進行了分析。綜上分析，相關(guān)研究主要集中在各種水體組分或者光學(xué)參數(shù)方面，對油類水體的研究還鮮見報道。

國際上對水色三要素（黃色物質(zhì)、葉綠素、懸浮物）的研究體系已經(jīng)比較成熟，將“未形成明顯油膜的水中油”作為一個新的水色因子引入該體系，針對石油類污染水體的“表觀光學(xué)量、固有光學(xué)量（吸收系數(shù)和散射系數(shù)）、熒光特性”等方面的研究工作已經(jīng)展開[11-15]。黃妙芬等[16]以雙臺子和繞陽河為例進行了石油類污染的水體后向散射系數(shù)b特性研究，發(fā)現(xiàn)在含有石油類污染的水體中b和懸浮物濃度接近線性關(guān)系；宋慶君等[17]進行了石油類濃度對b光譜影響的研究，發(fā)現(xiàn)bb隨著石油類濃度的增加而減?。稽S妙芬等[18]以石英砂作為懸浮物，進行油砂混合配比試驗，發(fā)現(xiàn)石油類吸附在懸浮物表面，增大了顆粒物粒徑，導(dǎo)致后向散射系數(shù)增大。目前關(guān)于含油水體吸收系數(shù)和散射系數(shù)的研究僅停留在隨波長的變化，未考慮隨深度的變化，而應(yīng)用輻射傳輸模式Hydrolight來研究石油類污染水體輻射傳輸特性，作為輸入量的吸收系數(shù)和散射系數(shù)要求是隨波長和深度變化的。本研究以石油類污染水體后向散射系數(shù)b為研究對象，探討其油濃度和b的垂向變化特性，以期為進一步進行石油污染水體水下光場研究奠定基礎(chǔ)。

1 材料與方法

1.1 研究區(qū)域與試驗時間

試驗區(qū)域位于遼東半島南端的大連港，港區(qū)內(nèi)有較大規(guī)模的石油化工基地，生產(chǎn)活動及油輪運輸不可避免地給周邊海域帶來油污染，是含油水體光學(xué)特性及油含量變化特征研究的理想試驗場[19]。本研究在該區(qū)域內(nèi)選擇了三個試驗點，分別位于曾經(jīng)發(fā)生過石油管道爆炸事故點（站點A）、航道（站點B）及島嶼東側(cè)水域（站點C），三個站點的深度分別為21.5 m、24.8 m和33.0 m。三種觀測點分屬不同類型：第一種是周邊有石油企業(yè)，生產(chǎn)活動會產(chǎn)生水體油污染；第二種是在航道上，過往的輪船以運輸原油的貨船為主，船舶排污水會導(dǎo)致水體油污染；第三種是選擇島的東邊，相對來說隨著潮汐變化受油污染的時間會短一些。

試驗時間為2018年8月25－27日，每天觀測時段為7:00－17:00，逢整點時間進行觀測。

1.2 現(xiàn)場數(shù)據(jù)采集與方法描述

1.2.1 采樣方法 水樣采集使用的是國家海洋技術(shù)中心研發(fā)的有機玻璃材質(zhì)采水器。在采水器的繩子上標有刻度，以便按要求采集不同深度（0、3、5、10、15 m）的水樣，所采集的水樣用于油濃度、懸浮物濃度、吸收系數(shù)的測量。為了保證數(shù)據(jù)的質(zhì)量，在8:00、10:00、12:00、14:00和16:00等觀測時刻同時采集平行樣以進行測量比對。

所采集的用于油濃度測量的水樣在船上直接用儀器進行測量；用于懸浮物濃度和吸收系數(shù)測量的水樣在船上先進行過濾，過濾后的樣本分別保存在冰柜和液氮罐中，待試驗結(jié)束后帶到岸上實驗室進行測量。流速和流向、后向散射系數(shù)的測量采用現(xiàn)場布放儀器直接獲取的方式。

1.2.2 油濃度測量 油濃度測量使用美國特納TD-500D便攜式紫外熒光測油儀進行。該儀器測定原理與分子熒光光度法（《SL 366-2006水質(zhì)石油類的測定分子熒光光度法》）相同。

1.2.3 懸浮物濃度測量 懸浮物濃度測量按照中華人民共和國國家標準《GB 17378.4-2007 海洋監(jiān)測規(guī)范第4部分：海水分析》規(guī)范進行。

1.2.4 吸收系數(shù)測量 吸收光譜測定采用日本日立UV-3900可見光分光光度計，測定時波長設(shè)置為250 ~ 800 nm，在樣品的制備、測量和分析過程都遵循《海洋光學(xué)調(diào)查技術(shù)規(guī)程》[20]。吸收系數(shù)的測量主要是為后向散射系數(shù)b進行sigma校正時提供參數(shù)。

1.2.5 后向散射系數(shù)測量 水體后向散射系數(shù)bb采用美國Hobilabs公司生產(chǎn)6通道后向散射儀（HydroScat-6 Sprctral Backscattering Sensor，HS6，140°）進行，6個通道對應(yīng)的波長分別為420、442、470、510、590、700 nm。后向散射系數(shù)的測量方法參照《海洋光學(xué)調(diào)查技術(shù)規(guī)程》，采用沿深度布放的方式。

現(xiàn)場測量的后向散射系數(shù)需要進行sigma校正[21]，校正公式如下：

式中，b為校正后的后向散射系數(shù)（m-1），為校正系數(shù)（無單位），bu為未校正的后向散射系數(shù)（m-1）。其中，有兩種計算方法，公式如下：

式（2）和（3）中的0、1、2和exp是與HS-6有關(guān)的定標參數(shù)(無量綱)，可從定標文件中得到，而bb為水體中的非水成分所造成的后向散射系數(shù)信號衰減，其計算公式如下：

根據(jù)HS-6的設(shè)計和光學(xué)傳輸幾何原理，信號傳輸過程中衰減量包括100%的吸收（）和40%的散射（）[22]。式（4）中和分別為吸收系數(shù)和散射系數(shù)，單位為m-1，其中由分光光度計測得，是顆粒物后向散射系數(shù)與后向散射概率之比。

1.2.6 流速與流向測量 流速與流向測量所使用的設(shè)備是國產(chǎn)的直讀式海流計，型號為SLC9-2。按照標準《GB/T 12763.2-2007 海洋調(diào)查規(guī)范第2部分：海洋水文觀測》進行，采用船只錨定測流法。

2 結(jié)果與討論

2.1 油濃度隨深度變化特征分析

對于油濃度采用插值法將其深度分辨率與處理后的后向散射系數(shù)統(tǒng)一起來，深度間隔取0.1 m，得到大連港三個觀測站點的油濃度隨深度的變化情況，如圖1所示，圖中橫坐標P為油濃度（mg/L），縱坐標為深度（m），圖例為曲線對應(yīng)的觀測時刻。由圖1可見，站點A、B和C油濃度分別在0.7 ~ 4.1 mg/L、1.0 ~ 9.3 mg/L和0.2 ~ 7.8 mg/L之間，B站點全天都處于高值狀態(tài)；整體上看，三個站點油濃度隨深度變化較復(fù)雜，無明顯增加或者減小趨勢。

表1為觀測期間海水流向及漲落潮詳情表。分析表1可見，7:00－10:00所有站點都是處于漲潮，10:00－16:00時都處于落潮，雖然同樣是漲落潮，但流向差別很大，說明海水的來向不同，其組分濃度不同，因而油含量也不同，加上海水的運動增加了油含量隨深度變化的復(fù)雜性。

對于油濃度隨深度的變化，在3 m深度附近大部分情況下有一個升高的現(xiàn)象。石油類物質(zhì)主要以漂浮油、分解油、分散油和乳化油的形式存在于水體中，漂浮油以油膜形式存在，分解油和分散油主要以黃色物質(zhì)的形式存在，而乳化油主要以顆粒物的形式附著于水體的懸浮顆粒物上。3 m深度附近處于海水混合層，受到風(fēng)浪的影響大，漲落潮的影響海水交換也比較頻繁，這可能是導(dǎo)致3 m處附近部分情況下油濃度增加的原因之一。

（a）站點A；（b）站點B；（c）站點C

表1 不同站點流向變化詳情

注：↗表示漲潮，↘表示落潮

2.2 含油水體后向散射系數(shù)bb隨深度變化特征分析

通過對所獲取的b值進行分析，發(fā)現(xiàn)含油水體b隨深度的變化大致可以歸為三種類型：（1）陡增型，其特點是隨著深度增加，b值變化不大，但在11.5 m處附近突然劇增，如圖2（a）所示。（2）平緩增加型：其特點是隨著深度增加，b值緩慢增加，如圖2（b）所示。（3）單峰型：該類型又可分為三種情況，峰值表層型，即峰值出現(xiàn)在深度0.5 ~ 2.5 m之間，如圖2（c）所示；峰值中層型，即峰值出現(xiàn)在深度7 ~ 10 m之間，如圖2（d）所示；峰值下層型，即峰值出現(xiàn)在深度10 ~ 12.5 m之間，如圖2（e）所示。

另外，從圖2還可以看到，（1）不同波段所對應(yīng)的b隨著深度的變化趨勢是一致的；（2）從量級上看，波長為700 nm所對應(yīng)的b最小，其次是590 nm，510、470、442和420 nm這4個波段的b值在量級上差別不大。

2.3 油物質(zhì)和懸浮顆粒物對bb的影響

考慮到油濃度和懸浮物濃度測定時所采樣的最大深度為15 m，而HS-6儀器布放的最大深度超過了15 m，為了統(tǒng)一起來，將b數(shù)據(jù)的最大深度取至15 m。由于b是6波段，為了分析b的影響，擬選擇其中一個波段的數(shù)據(jù)進行討論。前人的研究建立了b關(guān)于波長的參數(shù)化模型[23]，模型如下：

式（5）中，λ0為參考波長（nm），λ0取值為555 nm或HS-6的6個波段中與之差值最小的波段；n為散射波長變化指數(shù)，取值會因水體組分不同而有所差異。由于其他波段的bb可以通過該參數(shù)化模型進行計算，因而用555 nm對應(yīng)bb進行討論具有一定的代表性，但由于后向散射測量儀器HS-6沒有555 nm波段，故本文選取了與555 nm最接近的590 nm波段進行分析。

2.3.1 時間變化影響特征分析 將各個站點相同時間的所有深度的油濃度p、懸浮物濃度s和b觀測值分別求平均，得到每個站點相應(yīng)參數(shù)均值隨時間變化的數(shù)據(jù)，如圖3所示。

根據(jù)表1可知，三個站點在8:00－10:00點為漲潮，10:00－16:00為落潮。分析圖3可以看到，站點A的p在漲潮時降低，落潮時增加，而s和b則是漲潮時增加，落潮時降低。站點B的三個參數(shù)變化趨勢基本是一致的，都是漲潮時降低，落潮時增加，但14:00之后s和p突然降低，而b保持增加。C站點的p在12:00之后突然降低，s和b是先增加后降低。這是由于漲落潮海水來向不同，其組分不同，加上C站點位于島的東側(cè)，局部環(huán)流使得海水來向變得更加復(fù)雜所致。

分析圖3還可以看到，s和b的變化趨勢基本是一致的，但這種一致性的趨勢會受到p的影響而產(chǎn)生波動。這主要是由于油物質(zhì)對b的影響是通過懸浮物顆粒進行的，兩者的共同作用影響著含油水體的b，而且在含油水體中懸浮物濃度與b呈現(xiàn)出線性關(guān)系所致。

2.3.2 空間變化影響特征分析 圖4為三個站點在五個觀測時刻懸浮物濃度s隨深度的變化曲線，圖中s為懸浮物濃度（mg/L），為深度（m）。分析圖4可以看到，s在多數(shù)情況下隨深度有一定的增加趨勢，說明懸浮物的沉降速度大于水動力對其向上運輸?shù)乃俣?；總體上看，A站點的s高于B站點和C站點，B站點和C站點的s差別不大。

在石油類污染水體中，油和懸浮物混合后，油會附著在顆粒物的表面，形成一個雙層結(jié)構(gòu)，從而導(dǎo)致水體后向散射特性的改變，因而利用HS-6測量得到的b是油和懸浮物共同作用的結(jié)果。圖5為三個站點的b隨深度變化曲線，對比分析圖4和圖5可見，從整體上看，對應(yīng)站點的s與b隨深度的變化趨勢基本是一致的，這是由于油物質(zhì)的附著改變了顆粒物的折射系數(shù)和粒徑，從而改變后向散射系數(shù)所導(dǎo)致，也體現(xiàn)出了在含油水體中，懸浮物濃度與b的呈現(xiàn)出線性關(guān)系。

圖4 懸浮物濃度隨深度變化曲線

圖5 三個站點后向散射系數(shù)隨深度變化曲線

3 結(jié)論

本研究發(fā)現(xiàn)，在含油水體中油濃度和后向散射系數(shù)垂向變化受到潮汐的影響；石油物質(zhì)通過吸附在懸浮顆粒物表面形成一個雙層結(jié)構(gòu)，共同影響著水體的后向散射系數(shù)，且懸浮物濃度與后向散射系數(shù)呈現(xiàn)出線性變化關(guān)系。

本研究所得結(jié)論針對大連港海域，今后將拓展到更多的海域展開研究，以提高研究的精度和結(jié)論的適用性。

[1] 黃妙芬, 駱蔚健, 劉楊, 等. 應(yīng)用Hydrolight模式研究含油水體輻射傳輸特性的關(guān)鍵技術(shù)[J]. 海洋技術(shù)學(xué)報, 2019, 38(3): 7-14.

[2] 龐重光, 韓丹岫, 趙恩寶. 黃東海懸浮泥沙濃度的垂向分布特征及其回歸模型[J]. 泥沙研究, 2009(5): 69-75.

[3] 何劍波. 千島湖水溫、溶解氧及葉綠素a垂向特征研究[D]. 杭州: 浙江工業(yè)大學(xué), 2014.

[4] XIU P, LIU Y, TANG J. Variations of ocean colour parameters with nonuniform vertical profiles of chlorophyll concentration[J]. International Journal of Remote Sensing, 2008, 29(3): 831-849.

[5] STRAMSKA M, STRAMSKI D. Effects of a nonuniform vertical profile of chlorophyll concentration on remote-sensing reflectance of the ocean[J]. Applied Optics, 2005, 44(9): 1735-1747.

[6] 梁其椿, 張玉超, 薛坤, 等. 巢湖藻類高斯垂向分布結(jié)構(gòu)參數(shù)的遙感估算[J]. 湖泊科學(xué), 2017, 29(3): 546-557.

[7] 馬孟梟, 張玉超, 錢新, 等. 巢湖水體組分垂向分布特征及其對水下光場的影響[J]. 環(huán)境科學(xué), 2014, 35(5): 1698-1707.

[8] 徐文龍, 王桂芬, 周雯, 等. 南海東北部夏季葉綠素a濃度垂向變化特征及其對水動力過程的響應(yīng)[J]. 熱帶海洋學(xué)報, 2018, 37(5): 62-73.

[9] YIN X J, CHEN L J, TANG D Q, et al. Seasonal and vertical variations in the characteristics of the nitrogen-related functional genes in sediments from urban eutrophic lakes[J]. Applied Soil Ecology, 2019, 143: 80-88.

[10] XING X G, QIU G Q, BOSS E, et al. Temporal and vertical variations of particulate and dissolved optical properties in the South China sea[J]. Journal of Geophysical Research, 2019, 124: 3779-3795.

[11] HUANG M F, XING X F, SONG Q J, et al. A new algorithm of retrieving a petroleum substances absorption coefficient in sea water based on a remote sensing image[J]. Acta Oceanologica Sinica, 2016, 35(11): 97-104.

[12] HUANG M F, SONG Q J, XING X F, et al. A biooptical model of retrieving petroleum concentration in sea water [J]. Acta Oceanologica Sinica, 2014, 33(5): 81-85.

[13] 黃妙芬, 宋慶君, 陳利搏, 等. 基于歸一化遙感反射比反演水中石油含量模式研究[J]. 海洋技術(shù)學(xué)報, 2015, 34(1): 1-9.

[14] 黃妙芬, 王忠林, 劉遠, 等. 遼東灣海域有色可溶性有機物熒光特征分析[J]. 大連海洋大學(xué)學(xué)報, 2018, 33(2): 137-144.

[15] 黃妙芬, 宋慶君, 邢旭峰, 等. 石油污染水體熒光圖譜特征分析[J]. 光譜學(xué)與光譜分析, 2014, 34(9): 2466-2471.

[16] 黃妙芬, 宋慶君, 唐軍武, 等. 石油類污染水體后向散射特性分析: 以遼寧省盤錦市雙臺子河和繞陽河為例[J]. 海洋學(xué)報(中文版), 2009, 31(3): 12-20.

[17] 宋慶君, 黃妙芬, 唐軍武, 等. 水中石油類含量對后向散射系數(shù)光譜的影響[J]. 光譜學(xué)與光譜分析, 2010, 30(9): 2438-2442.

[18] 黃妙芬, 邢旭峰, 宋慶君, 等. 油砂混合水體后向散射系數(shù)光譜貢獻分離算法Ⅰ: 實驗理論[J].光譜學(xué)與光譜分析, 2017, 37(1):205-211.

[19] 黃妙芬, 劉楊, 邢旭峰, 等.大連港海域油含量時空變化特征分析[J]. 大連海洋大學(xué)學(xué)報, 2020, 35(2): 273-279.

[20] 國家海洋局908專項辦公室. 海洋光學(xué)調(diào)查技術(shù)規(guī)程[M]. 北京: 海洋出版社, 2006.

[21] 宋慶君, 唐軍武, 馬榮華. 水體后向散射系數(shù)校正方法研究[J]. 海洋技術(shù), 2008, 27(1): 48-52.

[22] BOSS E. Particulate backscattering ratio at LEO 15 and its use to study particle composition and distribution[J]. Journal of Geophysical Research Atmospheres, 2004, 109(C1): C01014.

[23] 宋慶君, 唐軍武. 黃海、東海海區(qū)水體散射特性研究[J]. 海洋學(xué)報, 2006, 28(4): 56-63.

Vertical Variation Characteristics of Petroleum Concentration and Backscattering Coefficient in Petroleum-polluted Water Body

XING Xu-feng, WANG Zhong-lin, HUANG Miao-fen, HUANG Yin-en, ZHUANG Yang

（,,524088,）

【Objective】The study lays a foundation for the study of underwater light field in petroleum-polluted water by grasping the vertical distribution characteristics of petroleum concentration and backscattering coefficient in petroleum-polluted water. 【Method】August 25th to 27th, 2018, in-situ data including petroleum concentration, suspended matter concentration and backscattering coefficient was adopted. The data were measured by using the mode that the station is fixed and the time changes with tide at Dalian port. 【Result】The regime of vertical variation of petroleum concentration is complicated, and the concentration does not decrease or increase gradually with depth increasing. The minimum and second minimum of backscattering coefficient of each station are at 700nm and 590nm respectively. By the way, magnitude of backscattering coefficient is similar at 420nm, 442nm, 470nm and 510nm, and backscattering coefficients of these stations vary with depth consistently. There are three kinds of vertical variations of total backscattering coefficient in this study area, they are abrupt growth, steady growth and single hump. 【Conclusion】The vertical variation of petroleum concentration and backscattering coefficient is mainly affected by tidal action due to the fact that ebb and flow of different time and stations could influence current in different way.In addition, the two-layer structure formed by the adsorption of petroleum substances on the surface of suspended particles affects the backscattering coefficient of water.

petroleum concentration; suspended particles; backscattering coefficient; vertical variation; Dalian port

P76

A

1673-9159（2020）04-0075-07

10.3969/j.issn.1673-9159.2020.04.011

2020-02-25

國家自然科學(xué)基金項目（41771384）；廣東海洋大學(xué)2017年“創(chuàng)新強校工程”自主創(chuàng)新能力提升項目（GDOU2017052501）；廣東海洋大學(xué)科研啟動經(jīng)費資助項目（E16187）

邢旭峰(1972－)，男，副教授，碩士，主要研究方向為海洋信息采集技術(shù)。E-mail: xingxf@gdou.edu.cn

邢旭峰，王忠林，黃妙芬，等. 石油污染水體油濃度及后向散射系數(shù)垂向變化特性分析[J]. 廣東海洋大學(xué)學(xué)報，2020，40（4）：75-81.