舟山近海域懸浮泥沙濃度遙感反演研究

錢揚風， 郭碧云， 黃 東， 羅鳳云， 韓錄維

(浙江海洋學院, 浙江舟山 31600)

?

舟山近海域懸浮泥沙濃度遙感反演研究

錢揚風， 郭碧云*， 黃 東， 羅鳳云， 韓錄維

(浙江海洋學院, 浙江舟山 31600)

摘要[目的]采用衛星遙感方法反演舟山海域懸浮泥沙的分布，并進一步探索適合的反演模型，為以后舟山海域懸浮泥沙遙感研究提供參考依據。[方法]利用實測數據、遙感數據和實驗室分析的方法，對舟山海域不同水體進行光譜分析，建立水體遙感反射率與懸浮泥沙濃度之間3種不同的反演模型。[結果]遙感參數反演模型的反演效果最好；對數反演模型次之，對數模型下TM4波段的相關系數最大，TM3波段的相關系數最??；線性反演模型最差，線性模型下TM2波段的相關系數最大，TM3波段的相關系數最小。[結論]在舟山海域采用遙感參數模型可以獲得較高的反演精度。

關鍵詞舟山海域；懸浮泥沙；遙感反射率；反演模型

Remote Sensing Retrieval of Suspended Sediment Concentration in Near Sea Area of Zhoushan

QIAN Yang-feng, GUO Bi-yun*, HUANG Dong et al(Zhejiang Ocean University, Zhoushan, Zhejiang 31600)

Abstract[Objective] The distribution of suspended sediment in Zhoushan sea area was retrieved by using satellite remote sensing method, and the suitable inversion model was further explored, which will provide reference for remote sensing study of suspended sediment in Zhoushan sea area.[Method] Using field data, remote sensing data and laboratory analysis method, spectral analysis of different water bodies in Zhoushan sea area was carried out, and three kinds of inversion models between water body remote sensing reflectance and suspended sediment concentration were established.[Result] The results showed that the inversion effect of remote sensing parameters model was the best.The second was Logarithmic inversion model.The correlation coefficient was the largest in TM4 band and the smallest in TM3.The linear inversion model was the worst.The correlation coefficient was the largest in TM2 band and the smallest in TM3.[Conclusion] The higher inversion accuracy can be obtained by using the remote sensing parameters model in Zhoushan sea area.

Key wordsZhoushan sea area; Suspended sediment; Remote sensing reflectance; Inversion model

水體中的懸浮泥沙含量是重要的水質參數之一。懸浮泥沙濃度對研究水流挾沙能力和河口水流泥沙運動規律等具有重要作用，對研究海岸帶的水質、生態、環境以及海岸工程、港口建設等具有重要意義[1]。傳統懸浮泥沙調查方法是利用船只在相關海域進行逐點采樣分析，速度慢，周期長，只能獲得在時間和空間分布上都很離散的少量數據，而且這種傳統的調查方法存在著很多實際操作上的困難。隨著衛星遙感技術的發展，衛星遙感監測逐漸從定性發展到定量，遙感精度逐漸提高。利用衛星遙感大范圍、周期短的優點，可以實現大面積水域懸浮泥沙動態、連續的觀測，為近海海域懸浮泥沙研究提供了新的有效手段。

關于我國海區懸浮泥沙遙感的研究較多，但比較集中在一些重點河口和港灣，如長江口、珠江口、杭州灣、渤海灣等，對舟山海域懸浮泥沙遙感研究較少，建立的反演模型不夠完善。舟山海域緊鄰長江入?？?，島嶼星羅棋布，近岸海流運動復雜，影響泥沙運動的水動力機制多樣，其特殊的地理位置、島嶼的緊密分布以及人類活動對舟山海域懸浮泥沙分布有較大影響。對舟山海域懸浮泥沙分布的研究可以為舟山未來的海洋生態建設、海洋工程建設與海洋資源的開發利用提供參考依據。筆者在前人研究的基礎上，探索適合舟山海域的懸浮泥沙反演模型。

1研究區概況

舟山群島位于浙江省東北部，長江口南側，杭州灣外緣的東海洋面上，地理坐標為121°31′~123°25′E、29°32′~31°04′N，總面積22 000 km2，其中海域面積20 800 km2，陸域面積1 440 km2。舟山群島由1 390個大小島嶼組成，是全國第一大群島(圖1)。舟山群島為海島丘陵地貌，群島呈西南—東北走向，南部大島較多，海拔較高，排列密集；北部以小島為主，地勢漸低，分布稀疏。舟山群島屬北亞熱帶南緣季風海洋型氣候，四季分明，冬暖夏涼，年溫適中，光熱充裕，蒸發量大，降雨偏少[2]。2011年6月30日，國務院正式批準設立浙江舟山群島新區，也是首個以海洋經濟為主題的國家級新區。研究舟山近海域懸沙分布對舟山建設國家級新區，發展海洋經濟具有重要意義。

2數據來源與處理

2.1實測數據 該研究采用的測量數據是2015年7月2日搭乘浙江海洋學院浙?？?號科考船出海測量所得，觀測點共6個(圖2)。采用合肥儀思特光電技術有限公司生產的ISI921VF系列野外地物光譜輻射計對不同海域水體的光譜值進行測量，獲得該海域水體的各種光譜數據，同時利用采水器獲得不同觀測點位置的水樣。

圖2　測量點位置示意Fig.2　Position of observation point

2.2數據處理

2.2.1等效遙感反射率計算。建立懸浮泥沙實際濃度與遙感反射率之間的關系，需將各種實地測量的光譜值轉化為遙感反射率。由遙感反射率Rs的定義可知：

Rs=Lw/Es

(1)

式中，Lw為離水輻亮度；Es為水面入射輻照度。這2個值都不是直接的測量值,需對測量數據處理分析才能得到。海水的離水輻亮度Lw與海水輻射率Lseawater、天空漫散射光信號Lsky的關系為：

Lseawater= Lw+ rLsky+ Lg

(2)

式中，海水輻射率Lseawater、天空漫散射光信號Lsky是直接的測量值，r是水-氣界面反射率, 也稱為菲涅耳反射系數。Lg為太陽的鏡面輻射率，通常情況下這一項很難測量，但如果采用適當的測量方位角和儀器觀測角度，此項的作用可以避免。在避開太陽直接反射、忽略或避開白帽的情況下，可以忽略太陽鏡面反射率Lg。式(2)化簡可得到海水離水輻射率：

Lw=Lseawater- rLsky

(3)

水面入射輻照度Es與標準板的反射信號Lp及標準板的反射率ρP有關,兩者均是直接的測量值，計算公式為：

Es= π Lp/ρP

(4)

推算遙感反射率的目的是建立其與懸浮泥沙濃度的相關關系，但是反演泥沙分布使用的是TM衛星影像，所以在進行分析時,要實現光譜儀測得的光譜值與TM數據之間的等效遙感反射率的計算[3]。TM衛星影像指美國陸地衛星Landsat4-5號專題制圖儀(thematicmapper)所獲取的多波段掃描影像，分為7個波段：TM1(450~520nm)藍光波段、TM2(520~600nm)綠光波段、TM3(630~0.69nm)紅光波段、TM4(760~900nm)近紅外波段、TM5(1 550~1 750nm)中紅外波段、TM6(1 040~1 250nm)熱紅外波段、TM7(2 080~2 350nm)中遠紅外波段。該研究選用前4個波段，將遙感反射率等效計算為TM影像上1~4波段的反射率，結果見表1。

表1　等效計算各TM波段下對應的遙感反射率

2.2.2懸浮泥沙濃度測量。在實驗室采用過濾稱重法對在不同觀測點采集的水樣泥沙濃度進行測量。先定量量取200 mL搖勻水樣并置于漏斗中，用抽濾機將漏斗中的海水快速抽離，得到覆蓋懸浮泥沙的濾紙片。然后將覆蓋懸浮泥沙的濾紙片中的水分烘干，用電子天平稱量得到烘干后濾紙與懸浮泥沙總質量，減去提前測出的濾紙片質量即可得到懸浮泥沙質量，最后除以水樣體積即為懸浮泥沙含量。懸浮泥沙含量計算公式為：

懸沙含量=懸沙質量/水樣體積

(5)

2015年7月2日各取樣點表層含沙量見表2。

3結果與分析

3.1懸浮泥沙水體的光譜特征分析運用MATLAB繪制出遙感反射率光譜曲線，結果見圖3。不同波段電磁波對海水的穿透能力存在差異，不同含沙量水體對太陽輻射光譜的反射率明顯不同。從圖3可以看出，隨著水體含沙量的增加，各波段的反射率都普遍增大，但其增幅不同。懸浮泥沙光譜反射率普遍具有雙峰特征，第一反射峰位置在500～600 nm波段，第二反射峰位置在750～850 nm波段。當水體中懸浮泥沙含量增加時，波譜上的反射峰由短波向長波方向移動，即所謂的“紅移現象”。隨著水體中懸浮泥沙濃度改變，各波段上反射率的改變存在差異。有的波段改變大，說明該波段對懸浮泥沙變化敏感，懸浮泥沙信息比較豐富；而有的波段基本不改變，即該波段基本不含有懸浮泥沙信息，因此在反演過程中要選擇對懸浮泥沙濃度敏感的波段[4-5]。對此，考慮建立不同的反演模型來探求舟山海域懸浮泥沙的適合反演模型。

表2不同觀測點表層水體含沙量

Table 2The sediment concentration of surface water in various observation sites

樣點標號Locationlabel濃度Concentrationmg/L經緯度Latitudeandlongitude1112.75122°20.463'E、29°53.908'N2108.25122°11.680'E、29°55.099'N3108.00122°15.191'E、29°51.302'N495.75122°18.132'E、29°54.248'N578.50122°17.880'E、29°50.638'N678.25122°13.547'E、29°52.513'N

圖3　不同濃度懸浮泥沙水體反射波譜曲線Fig.3　Reflection spectrum curve of suspended sediment in different concentration

3.2懸浮泥沙反演模型的建立

3.2.1線性模型。不同TM波段對懸浮泥沙的敏感程度不同，得到的反演模型精度也不同，據此建立不同TM單波段與遙感反射率的線性相關模型。首先在整個波段，如圖3所示，每個固定波長下存在6個觀測點遙感反射率，將這6個觀測點的遙感反射率與其對應的6個觀測點的懸浮泥沙濃度比較，對遙感反射率與懸浮泥沙濃度兩組數作線性回歸分析得出該波長下的線性相關系數，進而推廣至整個波段，計算出每個波長下采樣點遙感反射率與懸浮泥沙濃度的線性相關系數，結果見圖4。

圖4　波長與線性相關系數對應曲線Fig.4　The corresponding curves of wavelength and linear correlation coefficient

從圖4可以看出，在400~900 nm波段下相關系數均非常大，平均值在0.9以上。在900 nm之后的波段出現不規則大幅波動，平均相關系數不大，故不考慮900 nm之后的TM5、TM6、TM7波段。分別在TM1、TM2、TM3、TM4波段上，對不同TM波段等效遙感反射率與泥沙濃度做線性相關分析，并得到線性回歸方程與相關系數(表3)。從表3可以看出，TM2波段下的相關系數最大，TM1與TM4波段次之，TM3波段下的相關系數最小。

表3　線性回歸方程與相關系數

注：參數X代表遙感反射率；Y代表懸浮泥沙濃度。

Note: ParameterXstands for remote sensing reflectance;Ydenotes suspended sediment concentration.

3.2.2對數模型。2000年朱立俊等[6]根據不同海區的陸地衛星資料及地面同步實測數據，建立了水體懸浮泥沙對數反演模型，其表達式為：

X= - 5 714 + 0.595 2 log10S

(6)

式中，X為陸地衛星TM3波段與TM2波段(或MSS5波段與MSS4波段)反射率比值；S為水體懸浮泥沙濃度[6-8]。參考建立此對數模型的方法，建立適合舟山海域的對數模型。采用“3.2.1”同樣的數據分析方法，得出整個波長下遙感反射率與懸浮泥沙濃度的對數相關系數，結果見圖5。從圖5中可以得出與表3線性模型相似的結論，在400~920 nm波段下相關系數均非常高，平均相關系數在0.9以上。在920 nm之后的波段出現不規則大幅波動，不適宜取TM5、TM6、TM7波段，故略去。據此在TM1、TM2、TM3、TM4波段建立回歸方程，得到對應的相關系數(表4)。從表4可以看出，總體上4個波段相關系數都比較大，不同波段之間差別不明顯。其中TM4波段下的相關系數最大，TM1、TM2波段次之，TM3波段下的相關系數最小。

圖5　波長與對數相關系數對應曲線Fig.5　The corresponding curves of wavelength and logarithm correlation coefficient

Table 4 Logarithmic regression equation and correlation coefficient

波段Band∥nm回歸方程Regressionequation相關系數CorrelationcoefficientTM1450~520Y=22.435lnX+57.160.9814TM2520~600Y=25.548lnX+44.110.9832TM3620~690Y=23.091lnX+53.860.9814TM4760~900Y=12.648lnX+87.530.9867

注：參數X代表遙感反射率；Y代表懸浮泥沙濃度。

Note: Parameter X stands for remote sensing reflectance; Y denotes suspended sediment concentration.

3.2.3遙感參數模型。550 nm波段對葉綠素濃度較為敏感，而490 nm波段對葉綠素濃度也較為敏感。通過兩者的相關性分析，發現490 nm處的遙感反射率和550 nm處的遙感反射率高度相關，兩者相除，可以去除葉綠素對懸浮泥沙遙感信息的干擾，于是韓震等[9]提出采用遙感參數XS來建立與泥沙濃度的反演模型，遙感參數XS即：

XS=[R(λ5)+R(λ6) /R(λ4)]/R(λ5)

(7)

式中，R為各個波段下的遙感反射率，λ4、λ5、λ6分別對應490、550、670 nm波段。R(λ5)+R(λ6)可以較好地反映出懸浮泥沙的特征，R(λ4)/R(λ5)可以去除葉綠素對懸浮泥沙遙感信息的干擾[9-10]。利用公式(7)計算得到6個觀測點的遙感參數XS(表5)，與對應的懸浮泥沙濃度做相關分析，結果發現當遙感參數XS與懸浮泥沙濃度呈對數關系時相關系數最大。圖6為遙感參數XS與懸浮泥沙濃度的擬合曲線。圖6中反演模型的相關系數為0.993 4，是3個模型中最高的，可見遙感參數模型對懸浮泥沙濃度遙感反演結果最好。

表5　不同采樣點的遙感參數XS

4結論

懸浮泥沙分布研究對舟山海域的生態建設與經濟建設具有重要意義。該研究基于遙感技術在該海域不同海區進行現場測量，建立遙感反射率與懸浮泥沙濃度之間3種不同的反演模型，結果表明：

(1)線性模型是3個模型中相關系數最小的，在TM1~TM4波段中相關系數最大的在TM2波段，最小在TM3波段。

圖6　泥沙遙感參數與實際含沙量的相關性分析曲線Fig.5　 Correlation analysis curve between remote sensing parameters and actual sediment

從整個波段上來看，平均相關系數最好的波段是420~600 nm和730~830 nm,兩者的平均相關系數均在0.94以上。

(2)對數模型的相關性在3者中處于中間位置，相關系數最大的是在TM4波段，最小的是在TM3波段。

(3)遙感參數模型是3個模型中相關系數最大的，相關系數高達0.993 4。在以后對舟山海域懸浮泥沙遙感反演研究中可以首先考慮遙感參數模型，以取得較高的精度。

參考文獻

[1] 汪小欽，陳崇成． 遙感在近岸海洋環境監測中的應用[J].海洋環境科學，2000，19(4):72-76．

[2] 宋亞民.舟山群島水文特性[J]. 水文，2001，21(6):59-62.

[3] 鐘凱文，劉旭，攏解靚. 基于遙感方法反演珠江三角洲西江干流懸浮泥沙分布研究[J] 遙感應用，2009,12(3):49-59.

[4] 李四海,惲才興.河口表層懸浮泥沙氣象衛星遙感定量模式研究[J].遙感學報,2001,5(2):154-161.

[5] 李炎,李京.基于海面遙感器光譜反射率料率傳遞現象的懸浮泥沙遙感算法[J]. 科學通報,1999,44(17):1892-1897.

[6] 朱立俊，尤玉明.遼東灣綏中海岸演變及懸沙分布特征的遙感分析[J].海洋工程，2000，18(1)：66-69.

[7] 李洪靈,張鷹,姜杰. 基于遙感方法反演懸浮泥沙分布[J ].水科學進展,2006 ,17(2):242-2451

[8] 鄧孺孺,何執兼,陳曉翔. 珠江口水域水污染遙感定量分析[J ].中山大學學報,2002 ,41(3):99-103

[9] 韓震,惲才興,蔣雪中.懸浮泥沙反射光譜特性實驗研究[J].水利學報,2003(12):118-122.

[10] 黎夏 .遙感定量的統一模式及其在珠江口中的應用[J ].環境遙感,1992,7(2): 106-113.

收稿日期2015-12-16

作者簡介錢揚風(1994- )，男，安徽舒城人，本科生，專業：海洋科學。*通訊作者，副教授，博士，從事海洋遙感和地理信息科學研究。

基金項目浙江省大學生科技創新活動計劃暨新苗人才計劃項目(2014R4 11036)；國家自然科學基金項目(51479179)；浙江海洋學院大學生科技創新項目；浙江海洋學院科研啟動項目。

中圖分類號S 951.4

文獻標識碼A

文章編號0517-6611(2016)03-294-04