樊彥國, 孟志河, 張 磊, 劉復(fù)生, 朱 浩

(1. 中國石油大學(xué)(華東) 地球資源與信息學(xué)院, 山東 青島 266555; 2. 中國石油天然氣華東勘察設(shè)計研究院巖土工程處, 山東 青島 266071)

基于遙感光譜反射率反演黃河口海域Ⅱ類水體懸浮泥沙濃度模型

樊彥國1, 孟志河2, 張 磊1, 劉復(fù)生1, 朱 浩1

(1. 中國石油大學(xué)(華東) 地球資源與信息學(xué)院, 山東 青島 266555; 2. 中國石油天然氣華東勘察設(shè)計研究院巖土工程處, 山東 青島 266071)

依據(jù)黃河口海域含沙水體野外實(shí)測光譜反射率為參考實(shí)驗(yàn)數(shù)據(jù), 結(jié)合水體表面光譜觀測數(shù)據(jù)建立了經(jīng)驗(yàn)回歸模型, 從而定量估測黃河口海域水體表層懸浮泥沙質(zhì)量濃度, 最終得出黃河口海域不同年份水體的懸浮泥沙質(zhì)量濃度圖, 這對分析河口海岸帶沖淤變化過程具有重要的意義。

黃河口; 懸浮泥沙濃度(SSC); 遙感; 光譜反射率

在海洋遙感研究中, 海洋水體按照光學(xué)性質(zhì)劃分為Ⅰ類水體(case 1 waters)和Ⅱ類水體 (case 2 waters)。其中, Ⅱ類水體的光學(xué)特性主要由懸浮物、黃色物質(zhì)(又稱有色可溶性有機(jī)物)決定; 近岸水體直接受到來自海岸帶和河口徑流攜帶的陸源物質(zhì)影響,屬于Ⅱ類水體。遙感監(jiān)測Ⅱ類水體的懸浮泥沙濃度(suspended sediment concentration, SSC)是海洋遙感領(lǐng)域的一個重要研究方向。掌握近海水體懸浮泥沙濃度和運(yùn)移特征是分析河口形態(tài)和演變規(guī)律的重要方面。Ⅱ類水體懸浮泥沙的濃度直接影響到水質(zhì)的透明度和水色等光學(xué)性質(zhì), 同時也影響到水生生態(tài)環(huán)境和河口海岸帶沖淤變化過程。借助遙感手段, 可以快速地獲取水體泥沙的信息, 對分析研究海岸帶沖淤變化具有重要的意義[1]。

1 建模依據(jù)

1.1 黃河口海域光譜特性

黃河口海域不同懸浮泥沙濃度水體光譜反射特性[2～4]：

1)黃河口含懸浮泥沙水體光譜反射率均小于0.07 sr-1。

2)不同懸浮泥沙濃度水體, 其光譜反射率曲線不同。

3)從波段來看, 在400～550 nm波段范圍內(nèi), 光譜反射率對懸浮泥沙質(zhì)量濃度的變化反映不敏感。在550～700 nm, 大體趨勢是懸浮泥沙質(zhì)量濃度升高, 水體的反射率也隨之提高, 但也存在懸浮泥沙質(zhì)量濃度為 403.2(14.9) mg/L水體的光譜反射率高于懸浮泥沙質(zhì)量濃度為1 207.3(72.9) mg/L和1 713.1(183.4) mg/L水體的光譜反射率的現(xiàn)象(其中, 括號外的數(shù)值為懸浮泥沙濃度, 括號內(nèi)的數(shù)值為有機(jī)物濃度)。

4)在750～900 nm波段, 光譜反射率隨水體中懸浮泥沙質(zhì)量濃度的增長呈現(xiàn)明顯的增長趨勢, 特別是位于810 nm左右的次主峰值, 變化更為顯著。當(dāng)懸浮泥沙質(zhì)量濃度小于 25 mg/L時, 次主峰幾乎不存在。這可能是由于該波段處于水體的強(qiáng)吸收區(qū), 當(dāng)水體懸浮泥沙質(zhì)量濃度過低時, 懸浮泥沙后向散射小, 水體表面的光譜反射率也很低, 甚至趨近于零。

1.2 理論依據(jù)

按照 Dekker等[5]的方法, 利用海水表面光譜測量數(shù)據(jù)可模擬Landsat TM/ETM+可見光至近紅外譜段的光譜響應(yīng), 結(jié)果如圖1所示。從圖1可知, 對于懸浮泥沙質(zhì)量濃度低并且有機(jī)物質(zhì)量濃度低的水體,在TM2波段范圍內(nèi), 光譜反射率隨懸浮泥沙質(zhì)量濃度的增加而顯著增大, 故而適合這類水體的遙感反演。當(dāng)水體懸浮泥沙質(zhì)量濃度較低(＜10 mg/L)時,TM4波段對懸浮泥沙質(zhì)量濃度的響應(yīng)不敏感, 而對懸浮泥沙質(zhì)量濃度較高的水體, 其光譜反射率隨懸浮泥沙質(zhì)量濃度的增長明顯, 因此TM4波段只適合高懸浮泥沙質(zhì)量濃度水體的反演。當(dāng)水體有機(jī)顆粒物質(zhì)質(zhì)量濃度較低(＜20 mg/L)時, TM3波段的光譜反射率隨懸浮泥沙質(zhì)量濃度的增加而增加, 但水體中有機(jī)顆粒物質(zhì)質(zhì)量濃度較高時, TM3光譜反射率與懸浮泥沙濃度的相關(guān)關(guān)系并不明顯。因此, 對于懸浮泥沙質(zhì)量濃度變異較大的黃河口海域, 很難基于單一波段準(zhǔn)確地反演水體懸浮泥沙質(zhì)量濃度。盡管研究表明黃河含懸浮泥沙水體光譜反射率與水體懸浮泥沙質(zhì)量濃度在MSS6和MSS7波段有較為滿意的對數(shù)關(guān)系, 但相關(guān)系數(shù)R2均小于 0.85。其他河口也有利用陸地衛(wèi)星MSS5或TM3單一波段數(shù)據(jù)成功反演懸浮泥沙質(zhì)量濃度的例子, 但都缺乏對特定水體進(jìn)行光譜特性的分析, 或水體中懸浮泥沙質(zhì)量濃度較低。為克服單一波段反演的不足, 許多學(xué)者嘗試?yán)枚鄠€波段數(shù)據(jù)的運(yùn)算結(jié)果來反演水體中的懸浮泥沙質(zhì)量濃度, 并取得了較好的效果[6]。

圖1 黃河口海域不同懸浮泥沙質(zhì)量濃度水體在TM各波段光譜特性Fig. 1 Simulated spectral responses of TM wave bands to the water with different suspended sediment contents at the Yellow River Estuary

2 懸浮泥沙質(zhì)量濃度信息反演模型

本文研究了懸浮泥沙質(zhì)量濃度與模擬所得的Landsat TM波段光譜反射率及其組合的回歸關(guān)系。

考慮到不同波段光譜反射率之比可部分消除懸沙顆粒折射系數(shù)和后向散射有效因子對光譜反射率的影響, 且應(yīng)用于遙感影像時, 可消除大氣影響的乘積效應(yīng), 故選擇建立RTM3/RTM2與懸浮泥沙質(zhì)量濃度之間的關(guān)系來估測水體表層的懸浮泥沙質(zhì)量濃度。

利用 SPSS統(tǒng)計軟件的回歸分析功能建立RTM3/RTM2與懸浮泥沙質(zhì)量濃度之間的關(guān)系。所用光譜反射率資料是1989年8月在黃河三角洲鄰近海域?qū)崪y而得。所用儀器為RS-2型四通道野外光譜儀, 其4個通道的光譜范圍與陸地衛(wèi)星中多光譜掃描儀 4個波段的范圍完全相同。光譜視場角15o, 在觀測時基本上消除了局部波浪起伏造成的誤差和照度變化造成的誤差。在實(shí)測中取兩次結(jié)果的平均值。每個光譜測量站均取水樣, 室內(nèi)分析其懸浮泥沙質(zhì)量濃度, 每個站位的光譜反射率及懸浮泥沙質(zhì)量濃度見表1。

表1 實(shí)測光譜反射率與水體懸浮泥沙質(zhì)量濃度的對數(shù)值Tab.1 Measured spectral reflectance and logarithm of suspended sediment concentrations

圖2 不同時間的黃河口海域懸浮泥沙質(zhì)量濃度Fig. 2 Suspended sediment concentrations of the Yellow River Estuary at different times

利用 SPSS軟件進(jìn)行統(tǒng)計分析得到懸浮泥沙質(zhì)量濃度M的數(shù)理統(tǒng)計模型如下：

應(yīng)用模型進(jìn)行圖像解譯, 得到不同時間的黃河口海域懸浮泥沙質(zhì)量濃度(圖2)。

3 結(jié)語

1)論文對黃河口海域含沙水體在其可見光至近紅外 4個波段的光譜特性進(jìn)行了模擬分析, 并結(jié)合表觀光譜觀測數(shù)據(jù)建立了經(jīng)驗(yàn)回歸模型, 從而可估測研究黃河口海域水體表層懸沙的質(zhì)量濃度, 最終得出該海域不同年份水體的懸浮泥沙質(zhì)濃度圖。

2)由于遙感影像獲取的是表層水體的光譜信息,所以由遙感影像反演只能得到表層水體的泥沙質(zhì)量濃度。而在實(shí)際應(yīng)用(例如海岸侵蝕沉積研究)中水底的泥沙質(zhì)量濃度往往是很有用的指標(biāo)。要想通過表層水體的泥沙質(zhì)量濃度來獲取深層水體的泥沙質(zhì)量濃度是一個需要進(jìn)行進(jìn)一步研究的課題。

[1] 劉志國, 周云軒, 蔣雪中, 等．近岸Ⅱ類水體表層懸浮泥沙濃度遙感模式研究進(jìn)展[J]．地球物理學(xué)進(jìn)展,2006, 21(1)： 321-326.

[2] 樊輝, 黃海軍, 唐軍武.黃河口水體光譜特性及懸沙濃度遙感估測[J]. 武漢大學(xué)學(xué)報(信息科學(xué)版), 2007,32(7)： 601-604.

[3] Forget P, Ouillon S, Lahet F,et al. Inversion of reflectance spectra of nonchlorophyllous turbid coastal waters[J]. Remote Sensing of Environment, 1999, 68(3)：264-272.

[4] Doxaran D, Froidefond J, Lavender S,et al. Spectral signature of highly turbid waters, application with SPOT data to quantify suspended particulate matter concentrations[J]. Remote Sensing of Environment,2002, 81(1)： 149-161.

[5] Dekker A G, Malthus T J, Wijnen M M,et al. The effect of spectral band width and positioning on the spectral signature analysis of inland waters[J]. Remote Sensing of Environment, 1992, 41(2/3)： 211-222.

[6] Doxaran D, Froidefond J, Castaing P. A reflectance band ratio used to estimate suspended matter concentrations in sediment-dominated coastal waters[J]. International Journal of Remote Sensing, 2002, 23(23)：5079-5085.

Based on remote sensing spectral reflectance inversion of suspended sediment concentration model of surfase water at the Yellow River Estuary

FAN Yan-guo1, MENG Zhi-he2, ZHANG Lei1, LIU Fu-sheng1, ZHU Hao1
(1. College of Geo-Resource and Information in China University of Petroleum, Qingdao 266555, China; 2. Geotechnical Engineering Office in East China Investigation and Design Institute, China National Petroleum Corporation, Qingdao 266071, China)

Feb., 3, 2010

the Yellow River Estuary; suspended sediment concentration (SSC); remote sensing; spectral reflectance

Based on the experimental data of spectral reflectance at the Yellow River Estuary, a regression model was established, combining spectral data of the water surface, which can be used to estimate suspended sediment concentrations and suspended sediment figures can be accessed. It is of great significance to analysis of changes in estuarine and coastal erosion and deposition process.

TP79

A

1000-3096(2010)09-0060-04

2010-02-03;

2010-06-27

山東省科技攻關(guān)項目(2008GG10009018); 山東省科技計劃項目(J08LD55)

樊彥國(1965-), 男, 河北望都人, 博士, 碩士導(dǎo)師, 從事3S技術(shù)在數(shù)字國土、城市及海岸帶方向的教學(xué)與研究工作, E-mail：ygfan@upc.edu.cn

(本文編輯： 劉珊珊)