不同厚度海上油膜高光譜遙感波段敏感性研究

劉丙新，李 穎，張至達，周 新

(大連海事大學航海學院，遼寧 大連 116026)

不同厚度海上油膜高光譜遙感波段敏感性研究

劉丙新，李 穎，張至達，周 新

(大連海事大學航海學院，遼寧 大連 116026)

利用地物反射率光譜結合星載高光譜遙感環境噪聲估計法，研究了星載高光譜遙感數據對較薄、中等、較厚原油油膜的響應特性，并提出了敏感波段.結果表明：對于較薄油膜(10 μm和50 μm)，Hyperion傳感器的第8—52波段能夠有效識別；對300 μm和1 000 μm的中等厚度油膜，Hyperion傳感器的探測通道為第26—53波段；對于1 500 μm和2 000 μm的較厚油膜， Hyperion的第15—23波段以及第27—49波段數據探測性能較好.

高光譜遙感；油膜厚度；敏感性

海上溢油是海洋污染的主要形式之一，溢油發生后，其位置、種類、面積和相對厚度等相關信息[1]是公眾和媒體十分關注的.在應用遙感進行溢油檢測方面，國內外對多光譜、熱紅外、雷達等均有廣泛的研究與應用[2-3]，但由于海洋環境復雜、海面大氣影響及水體對電磁波的散射與吸收作用，海面油膜遙感信息較弱，導致海面油膜信息提取時存在“同物異譜，同譜異物”的現象.隨著高光譜遙感技術的出現和發展，針對海面油膜信息的高光譜遙感探測技術研究也得到不斷發展，該技術可獲取地物目標近乎連續的反射光譜，從而根據光譜特征差異來區分海水與溢油目標.對于海上油膜的光譜特征，國內已進行了很多研究，并取得了一系列的成果[4-6].此外，Lammonglia等[7]對原油、油膜以及乳化后的油的光譜特征進行了實驗室測量與分析，指出三種狀態下均可表現出可見光近紅外的光譜特征；趙冬至等[5]對潤滑油、柴油和原油在不同厚度下的反射率進行了分析，并指出柴油隨厚度增加反射率增大，而潤滑油和原油則反之；張永寧等[8]根據油膜光譜特征選擇利用衛星通道以達到區分油水的目的.但是利用星載高光譜遙感數據對油膜和海水進行識別，尤其是不同油膜厚度敏感性的研究還鮮見報道.本文借鑒M.Wettle等[9]的方法，以Hyerion高光譜遙感數據為例，通過光譜響應函數對實測光譜進行重采樣，并計算Hyperion圖像數據的環境噪聲等效反射比，將地面實測數據進行歸一化處理，以此為基礎評價了該數據對不同厚度油膜和海水光譜差異的敏感性.遙感數據對油膜厚度的敏感性是傳感器進行油膜探測潛力的標志之一，是利用光學遙感進行溢油監測與識別的重要前期研究，對選擇適當的遙感數據、有效識別海上油膜信息具有一定的參考意義.

1 數據與方法

1.1 油膜光譜反射率測量

本文所采用的設備是ASD FieldSpec?3地物光譜儀，其光譜范圍為350～2 500 nm，采樣間隔為1.4 nm (350～1 000 nm)和2 nm (1 000～2 500 nm).實驗水槽采用70 cm高的圓桶，內涂黑色啞光漆、外用遮光布遮嚴，以模擬深水環境.測量地點選擇在周圍無高出建筑物、無雜散光的樓頂平臺.測量時天空晴朗無云，風力微弱，以保證水面基本處于相對靜止狀態.

利用100 mL注射器向水面緩慢滴入定量油，根據滴入油的體積和油膜所形成的面積，計算油膜厚度.數據采集前，儀器自動進行暗電流校正，測量過程中每隔3～5 min進行一次優化；測量時探頭距觀測表面50 cm，并保持垂直向下，視場角3°；參考板采用與ASD FieldSpec?3相配套的漫反射標準參考板；每次測量重復10次，剔除異常曲線后取目標反射亮度的平均值(Do)；在其前后準同步測量參考板反射太陽輻射的平均值(Dr).通過比值法計算測量目標的反射率[10]

式中：Ro，Rr分別為測量目標和參考板的反射率；Do，Dr分別為測量目標和參考板反射亮度的測量平均值；λ為波譜采樣點的中心波長.

1.2 Hyperion數據獲取與處理

研究中所需的影像數據采用了Hyperion于2006年5月6日獲取的影像數據，該數據包括了遼東灣區域.Hyperion是地球觀測一號星(EO-1)上搭載的推掃式高光譜成像光譜儀，其光譜范圍為356～2 577 nm，共有242個波段，光譜分辨率約為10 nm.經過輻射定標處理的波段數為198個，實際可用波段為196個[11].考慮到水體在大于1 000 nm后對光的強吸收作用，本文選擇波長小于1 000 nm的Hyperion波段進行分析.為獲取Hyperion數據上地物的反射率數據，本文利用ENVI 4.5的FLAASH軟件對原始數據進行大氣校正.

1.3 遙感影像環境噪聲等效反射比計算

對傳感器-大氣-目標系統噪聲的準確估計，能夠提高環境信息提取的準確性.為評估遙感系統提取環境變量的準確性和精確性，需要對遙感圖像的環境噪聲等效輻射率和環境噪聲等效反射比進行計算.環境噪聲當量由傳感器信噪比SNR和特定觀測圖像的特點構成.本文利用環境噪聲等效反射比來對其進行描述，該值與儀器的信噪比有關，同時受到大氣變量、汽-水界面、太陽光和天空光的散射等影響[9].公式為

RE=σ(R).

式中：RE是環境噪聲等效反射比值；σ(R)是圖像上一個覆蓋內部盡量均勻、光學深度很大水域的窗口內各波段反射率的標準差，通過調整窗口大小(3×3，5×5，7×7等)，使σ(R)達到收斂[12].在窗口位置選擇上，采用Wettle等人提出的局部收斂自動定位算法(automated local convergence locator，ALCL)[13].

1.4 油膜厚度敏感性評價

通過光譜測量獲得不同厚度情況下油膜的連續光譜，并以此計算其反射率，記為Rm.通過Hyperion數據光譜響應函數，對Rm進行濾波，得到Hyperion數據對不同厚度油膜的響應，結果光譜反射率記為Rr，并將Rr歸一化至Hyperion的環境噪聲等效反射比值光譜.以此將反射光譜拉伸到遙感數據可區分的量級，記為S.結果光譜揭示了Hyperion數據在區分油膜厚度方面的理論極限值，見圖1.

圖1 數據處理流程圖

2 結果與討論

2.1 反射率光譜特征

圖2是海水和不同厚度油膜的反射率光譜特征，油膜的厚度分別為10，50，300，1 000，1 500和2 000 μm.通過光譜曲線可知，在整個波段范圍內，油膜越薄，其測得的反射率值越大.同時可知，較薄的油膜(10，50 μm)的反射率總是高于海水的反射率，尤其在427～583 nm的可見光近紅外波段，油膜和海水的反射率光譜曲線差異較大，而在波長498～579 nm之間，較厚的油膜(1 500，2 000 μm)的反射率比海水的低.對于所有厚度油膜，在771 nm處具有明顯的峰值.

圖2 海水和不同厚度油膜反射率光譜曲線

2.2 實測反射率光譜經響應函數濾波后的光譜特征

圖3是不同厚度油膜的實測反射率光譜經Hyperion的光譜響應函數濾波所形成的光譜曲線，其中包含測量時背景海水的反射率光譜曲線.圖4是油膜與海水的實測反射率經傳感器響應函數處理并利用環境噪聲等效反射比值歸一化處理后的光譜曲線，圖5是經環境噪聲等效反射比值歸一化處理后不同厚度油膜的環境噪聲等效反射比值分別減去海水的環境噪聲等效反射比值所得到的光譜曲線.

圖3 傳感器波譜響應函數處理后反射率光譜

根據Wettle等[8]的研究，以環境噪聲等效反射比值為單位進行度量，若兩個地物在某波長位置的環境噪聲等效反射比值度量值絕對值大于1，在傳感器所獲取的圖像數據上，二者在該波長范圍內能夠有比較明顯的區別，則傳感器在該波長位置所獲取的數據能夠區別這兩個地物.由圖5可發現，較薄的油膜(10，50 μm)和海水的環境噪聲等效反射比值差異在427～875 nm間均超過1，即在此波段范圍內Hyperion數據能夠區分薄油膜和海水的差異，從而識別海上的薄油膜.在498～579 nm和620～854 nm范圍內，較厚的油膜(1 500，2 000 μm)和海水的環境噪聲等效反射比值的差異也均超過1，因此在該波長范圍內利用Hyperion數據進行厚油膜的識別更加有效.對于中間厚度(300，1 000 μm)的油膜，在427～590 nm范圍內與海水的環境噪聲等效反射比值的差異并不是很大，但是差異先為正，后為負，并且在610～885 nm范圍內與海水的環境噪聲等效反射比值差異均大于1.理論上通過Hyperion傳感器能夠比較好地識別海上原油油膜，薄油膜和厚油膜的環境噪聲等效反射比值差異主要體現在508～610 nm范圍內，并且能夠區分薄油膜和厚油膜.

圖4 歸一化至環境噪聲等效反射比的反射率曲線

圖5 油膜與海水歸一化反射率差異

3 結論

本文利用Hyperion高光譜數據對海上溢油油膜厚度探測的能力進行了評估，利用遙感數據對傳感器監測油膜厚度能力的分析，得到了進行油膜厚度監測的最有效的波段范圍，節約了經濟和時間成本.

通過對實測水體和不同油膜厚度數據重采樣并利用環境噪聲等效反射比值進行歸一化處理，再分別與海水的環境噪聲等效反射比值相減可得到光譜通道的敏感程度.對于較薄油膜(10，50 μm)，Hyperion傳感器的第8—52波段能夠有效識別；對300，1 000 μm的中等厚度油膜，Hyperion傳感器的探測通道為第26—53波段；對于較厚油膜(1 500，2 000 μm)， Hyperion的第15—23波段以及第27—49波段數據探測性能較好.

本文對Hyperion高光譜傳感器對海上原油油膜的探測性能進行了初步研究，獲取的星載高光譜遙感影像與實測光譜曲線完全同步，研究結果將更有實用性與針對性；另外研究精度還需要考慮具體海況、天氣條件及觀測角度等因素.

[1] FINGAS M，BROWN C E. Oil spill remote sensing：a review. Oil spill science and technology [M]. Boston:Gulf Professional Publishing，2011:111-169.

[2] BREKKE C，SOLBERG A H S. Oil spill detection by satellite remote sensing[J]. Remote Sensing of Environment，2005，95(2005)：1-13.

[3] DE BEUKELAER S M，MACDONALD I R，GUINNASSO N L，et al. Distinct side-scan sonar，RADARSAT SAR，and acoustic profiler signatures of gas and oil seeps on the Gulf of Mexico slope[J]. Geo-Marine Letters，2003，23(3)：177-186.

[4] 陸應誠，田慶久，王晶晶，等. 海面油膜光譜響應實驗研究[J]. 科學通報,2008，53(9)：1085-1088.

[5] 趙冬至，叢丕福. 海面溢油的可見光波段地物光譜特征研究[J]. 遙感技術與應用，2000，15(3)：160-164.

[6] 趙云升，吳太夏，羅楊潔，等. 水面溢油的多角度偏振與二向性反射定量關系研究[J]. 遙感學報,2006，10(3)：294-298.

[7] LAMMOGLIA T，FILHO C R D S. Spectroscopic characterization of oils yielded from Brazilian offshore basins：potential applications of remote sensing[J]. Remote Sensing of Environment，2011，115(10)：2525-2535.

[8] 張永寧，丁倩，高超，等. 油膜波譜特征分析與遙感監測溢油[J]. 海洋環境科學,2000，19(3)：5-10.

[9] WETTLE M，DANIEL P J，LOGAN G A，et al. Assessing the effect of hydrocarbon oil type and thickness on a remote sensing signal a sensitivity study based on the optical properties of two different oil types and the HYMAP and Quickbird sensors[J]. Remote Sensing of Environment,2009，113(2009)：2000-2010.

[10] 李穎，劉丙新，蘭國新，等. 有冰海區油膜光譜特征研究[J]. 光譜學與光譜分析，2010(4)：1018-1021.

[11] 陸應誠，田慶久，宋鵬飛，等. 海面油膜高光譜遙感信息提取[J]. 遙感學報,2009，13(4)：691-695.

[12] DEKKER A G，PETERS S. The use of the thematic mapper for the analysis of eutrophic lakes：a case study in the Netherlands[J]. International Journal of Remote Sensing，1993，14(5)：799-821.

[13] WETTLE M，BRANDO V E，DEKKER A G. A methodology for retrieval of environmental noise equivalent spectra applied to four Hyperion scenes of the same tropical coral reef[J]. Remote Sensing of Environment,2004，93(1-2)：188-197.

(責任編輯：方 林)

Study on the sensitivity of hyper-spectral imagery to detect oil film with different thickness

LIU Bing-xin，LI Ying，ZHANG Zhi-da，ZHOU Xin

(College of Navigation，Dalian Maritime University，Dalian 116026，China)

The reflectance of water and different crude oil film thickness was measured in this study. The spectral response functions of Hyperion sensor was used for filtering the measured reflectance to the theoretic spectral response. Then，these spectral response spectra were normalized to environmental noise equivalent reflectance. Final，subtract the environmental noise equivalent reflectance between crude oil film and water，according to which，the sensitivity of the sensor in oil film detecting could be evaluated. For thin oil film(10 μm and 50 μm)，the range for Hyperion sensor to identify the film is within the Band 8 to Band 52. Band 26—53 of Hyperion sensor are the most efficient ones to monitor the medium oi film(300 μm and 1 000 μm). For thick oil film(1 500 μm and 2 000 μm)，the range for Hyperion sensor to identify the film is within the wavelength from Band 15—23 and Band 27—49.

hyper-spectral remote sensing；oil film thickness；sensitivity

1000-1832(2015)04-0156-05

10.16163/j.cnki.22-1123/n.2015.04.032

2014-08-09

國家863計劃項目(2013AA7013037)；國家自然科學基金資助項目(51509030)；中央高校基本科研業務費資助項目(3132015006)；國家海洋公益項目(201205012，201305002)；遼寧省自然科學基金資助項目(2015020081).

劉丙新(1984—)，男，博士，副教授，主要從事遙感技術應用研究；通訊作者：李穎(1968—)，女，教授，主要從事海洋溢油遙感監測研究.

TP 722.4 [學科代碼] 420·2040

A