金宇豪

（1.中山大學(xué) 地理科學(xué)與規(guī)劃學(xué)院，廣東 廣州 510275）

3S 技術(shù)應(yīng)用

基于GOCI的輻射沙脊群海表流場遙感分析

金宇豪1

（1.中山大學(xué) 地理科學(xué)與規(guī)劃學(xué)院，廣東 廣州 510275）

針對江蘇沿海尤其是輻射沙脊群海域水沙動力環(huán)境復(fù)雜，難以實時更新實地海洋流場調(diào)查數(shù)據(jù)的問題，根據(jù)2014年1 月韓國高分遙感影像GOCI與現(xiàn)場同步潮汐數(shù)據(jù)資料，運用最大交叉相關(guān)系數(shù)法，遙感模擬分析了輻射沙脊群海表瞬時流場。結(jié)果表明，在輻射沙脊群海域瓊港-洋口港一帶，近岸海流場變化明顯，海流自沿岸灘涂帶流向外海；同時遠岸瞬時流場的渦旋及流場瞬時方向變化亦較清晰。通過遙感分析法觀測海表流場，在一定程度上解決了基礎(chǔ)資料更新難的問題，實現(xiàn)了實時檢測輻射沙脊群海表流場變化的目的。

海表流場；遙感分析；MCC法；輻射沙脊群

南黃海輻射沙脊群海域規(guī)模龐大﹑水沙動力環(huán)境復(fù)雜，在現(xiàn)有技術(shù)條件下，難以實現(xiàn)實時更新實地海洋流場調(diào)查數(shù)據(jù)?，F(xiàn)有的實測資料已無法滿足對輻射沙脊群進行科學(xué)研究﹑工程建設(shè)和環(huán)境保護的需求，嚴重制約了江蘇沿海開發(fā)的進度[1]。因此，利用遙感分析法，并結(jié)合數(shù)學(xué)建模方法獲取輻射沙脊群海表流場信息，對科學(xué)研究及沿海開發(fā)至關(guān)重要。關(guān)于南黃海海表流場的變化，曾有多種不同的觀點。Yanagi T[2-3]等根據(jù)其診斷模型的結(jié)果認為，冬季南黃海流場存在一個上下方向一致的順時針環(huán)流。Takahashi S[4]等采用一個風作用下的簡化地形模型進行研究，結(jié)果表明在冬季風作用下將在黃海西側(cè)產(chǎn)生逆時針環(huán)流，在東側(cè)產(chǎn)生順時針環(huán)流，兩個環(huán)流結(jié)構(gòu)都以陸架波的形式傳播，順時針環(huán)流將在南黃海占據(jù)主導(dǎo)地位。Lee S H[5]等對南黃海表層潮流場進行了模擬，得到輻射沙脊群表層潮流向北的結(jié)論。LIU G[6]等對南黃海M2分潮進行了研究，得出蘇北沿岸夏季潮余流為東北向的結(jié)論。

高分海洋遙感技術(shù)是近年來新興的探測技術(shù)。韓國的第一顆地球靜止軌道海洋水色衛(wèi)星GOCI具有高時間分辨率﹑全天候﹑大范圍的連續(xù)觀測能力，其在南黃海水動力環(huán)境﹑海表流場監(jiān)測中具有其他傳統(tǒng)傳感器無法比擬的優(yōu)勢。本文基于2014年1月GOCI衛(wèi)星遙感數(shù)據(jù)和同步海域潮汐資料，運用數(shù)學(xué)模型計算并分析了南黃海輻射沙脊群表層流場的特征。

1 研究區(qū)與數(shù)據(jù)來源

1.1 研究區(qū)概況

南黃海是東中國海的重要組成部分，南連東海寬闊海域，為太平洋海表潮波入口；北接北黃海和渤海，為口袋型海灣[7]。黃海北﹑南連接成“呂”字形海灣，平均水深為44 m，總面積為4.2×105km2。輻射沙脊群海域是南黃海的一部分，分布于蘇北海岸帶外側(cè)﹑黃海南部陸架海域，范圍從射陽港一帶以南到長江口一帶以北的蒿枝港：南北范圍介于32°00'N～33°48'N之間，長度為199.6 km；東西范圍介于120°40'E～122°10'E之間，寬度為140 km[8]。

1.2 遙感數(shù)據(jù)

輻射沙脊群海表流場分析的遙感數(shù)據(jù)選用2014年1月韓國GOCI衛(wèi)星遙感數(shù)據(jù)，以實時高時間分辨率遙感影像為主[9]。遙感數(shù)據(jù)質(zhì)量較好，影像色調(diào)充實，云量較少，空間信息非常清晰。GOCI影像是由COMS傳感器獲取，其軌道為地球靜止同步軌道，運行周期為23 h 56 min 4 s。GOCI具有高時間分辨率，時間分辨率為1 h，地面像元分辨率為500 m。

遙感影像數(shù)據(jù)預(yù)處理，主要包括幾何配準﹑圖像增強等，用以提高遙感影像質(zhì)量，減少誤差因素對數(shù)據(jù)的影響。

1.3 實測潮汐數(shù)據(jù)

輻射沙脊群潮汐數(shù)據(jù)選用兩組江蘇沿海環(huán)境監(jiān)測站洋口港實測潮汐數(shù)據(jù)：一組為2014年1月1日潮高值，另一組為2014年1月2日潮高值，數(shù)據(jù)采樣間隔為1 h，與遙感影像時間分辨率相吻合，數(shù)據(jù)質(zhì)量已得到驗證（見表1）。

表1 輻射沙脊群洋口港實測潮汐數(shù)據(jù)

2 研究方法與結(jié)果

2.1 流場示蹤水團追蹤方法

本文通過最大交叉相關(guān)系數(shù)（MCC）法建立模型，在高時間分辨率遙感影像中追蹤輻射沙脊群海表流場水團位移的情況[10]。MCC法[11]，即選取10 像素×10 像素的示蹤水團，在40像素×40像素的目標影像模塊中，進行逐一搜索匹配，并計算MCC值。設(shè)MCC系數(shù)R(m,n)為：

式中，f (i, j)為示蹤水團模塊的DN值；(i, j)為示蹤水團模塊的行列號；(m,n)為行列變量位移變化。

本文對示蹤水團的追蹤，主要是對輻射沙脊群海域懸浮泥沙水團移動變化情況的跟蹤，即按照模板窗口大小將大面積示蹤水團逐一細化，并在搜索窗口中對照追蹤定位的過程。本文所選用的template-window模塊的像素大小為10×10，search-window模塊的像素大小為40×40，即在search-window中通過MCC系數(shù)尋找到相關(guān)性最高且具有相同模板尺寸的templatewindow，如圖1所示?；谏鲜隼碚?，對研究區(qū)海洋表面流場水團進行逐一追蹤，以獲得研究區(qū)域完整的海洋表面流場瞬時流向情況。

圖1 MCC法定位海流位移示例（9:00～11:00）

通過MCC法遙感追蹤示蹤水團的優(yōu)勢在于，可精確地從微觀上定位水團下一時刻的具體流向情況；在計算海流場流速時，因MCC法能更好地考慮微觀水團的位移變化，消除誤差，所以其對移動水團的計算誤差相對較小。

2.2 海表流場分析結(jié)果

輻射沙脊群海表流場反演分析分別選取2014年1月1日上午及下午兩個研究時間段，處理范圍為121°E～122°E，32°N～33°N，以該時刻示蹤水團在下一相鄰序列遙感中得到的最大相關(guān)示蹤水團的追蹤方向為流場流向。圖2﹑3 分別為未作任何條件約束情況下上午時和下午時MCC相關(guān)系數(shù)流場反演結(jié)果，用以反映冬季表層流場情況。由圖2可知，輻射沙脊群瓊港-洋口港一帶，近岸上午時流場變化明顯，呈顯著落潮趨勢，海流自沿岸灘涂帶流向外海，且遠岸上午時海表流場的渦旋趨勢較清晰。由圖3可知，洋口港一帶，近岸下午時流場變化呈漲潮趨勢，海流自外海流向沿岸輻射沙脊群灘涂帶；遠岸下午時海表流場的渦旋變化呈順時針趨勢，比上午時流場流向更加明顯。由冬季南黃海輻射沙脊群海表流場可以看出，近岸潮汐上午時和下午時分別成落潮和漲潮趨勢，反演結(jié)果較清晰，遠岸海流呈渦旋趨勢，與實測潮汐規(guī)律基本吻合。

圖2 輻射沙脊群上午時海表流場圖

圖3 輻射沙脊群下午時海表流場圖

3 研究結(jié)果分析

3.1 對比分析表層流場

本文選取輻射沙脊群海域上午時和下午時海表流場反演結(jié)果，對比分析當MCC系數(shù)臨界值為0.6和0.8時，表層流場的海流流向情況。

從圖4可以看出，對于上午時流場，當Rmcc>0.6時（圖 4a），近岸海表流場中不合理的瞬時海流流向基本被消除，而遠岸海表流場海流整合篩選情況較少，部分不合理海流分布并未消除；當Rmcc>0.8時（圖4b），遠岸海表流場中不合理的瞬時海流流向亦被完全篩去，同時獲得了近海和遠岸完整合理的海表流場。對于下午時流場，當Rmcc>0.6時（圖4c），近岸海表流場呈清晰落潮趨勢，不合理流場影響較小，而遠岸海表流場海流呈順時針趨勢，不合理海流分布未得到較好的整合篩選；當Rmcc>0.8時（圖4d），遠岸海表流場海流流向清晰，基本已去除不合理流場干擾，獲得較為準確的海表流場分析結(jié)果。

圖4 海表流場系數(shù)控制結(jié)果圖

3.2 漂流的分布

輻射沙脊群長時間連續(xù)海表流場情況可通過在該海域放置漂流浮標進行探究，圖5為2013年12月14 日～2014年1月10日的輻射沙脊群漂流實驗，記錄了漂流浮標經(jīng)過竹根沙﹑苦水洋﹑黃沙洋﹑蔣家沙﹑太陽沙﹑爛沙洋﹑草米樹洋的連續(xù)漂流軌跡，并結(jié)合期間GOCI影像，對比分析了實驗區(qū)模擬情況。漂流軌跡表明，輻射沙脊群海表流場漂流特征為自北向南全局性漂移，與輻射沙脊群受蘇北沿岸海流影響的預(yù)知條件吻合。其中，海表流場南部爛沙洋為典型旋轉(zhuǎn)流，流速最強；北部草米樹洋為往復(fù)流，流速較小。此外，輻射沙脊群海域海表流場遠岸爛沙洋一帶為典型旋轉(zhuǎn)流，亦與流場遙感分析結(jié)果吻合。因此，南黃海輻射沙脊群漂流浮標實測數(shù)據(jù)清楚地反映了該時段流場變化的特點，可從宏觀上驗證輻射沙脊群海表流場的遙感分析結(jié)果。

圖5 研究區(qū)漂流軌跡

4 結(jié) 語

通過研究，本文得到以下結(jié)論：①將遙感技術(shù)與MCC法相結(jié)合對海表流場進行模擬與分析是切實可行的。本文綜合運用遙感技術(shù)﹑MCC法，得出南黃海輻射沙脊群海表流場的瞬時運動趨勢；并通過海流場漂流浮標實測數(shù)據(jù)﹑潮汐實測數(shù)據(jù)驗證了輻射沙脊群海流場瞬時運動情況的準確性。②通過遙感分析可知，冬季輻射沙脊群瓊港-洋口港一帶，近岸海流場變化明顯，上午時呈顯著落潮趨勢，海流自沿岸灘涂帶流向外海，下午時呈漲潮趨勢，海流沿灘涂帶流出；與此同時，遠岸瞬時流場的渦旋及流場瞬時方向變化亦較清晰，呈順時針趨勢。

[1] 江蘇省908專項辦公室.江蘇金海海洋綜合調(diào)查與評價總報告[R].北京:科學(xué)出版社,2012:337

[2] Yanagi T. Fundamental Study on the Tidal Residual Circulation, II[J]. Journal of Oceanography,1976,32(5):199-208

[3] Yanagi T, Takahashi S. Seasonal Variation of Circulations in the East China Sea and the Yellow Sea[J]. Journal of Oceanography,1993,49(5):503-520

[4] Takahashi S, Isoda Y, Yanagi T. A Numerical Study on the Formation and Variation of a Clockwisecirculation During Winter in the Yellow Sea[J]. Journal of Oceanography,1995,51(1):83-98

[5] Lee S H, Beardsley R C. Influence of Stratification on Residual Tidal Currents in the Yellow Sea[J]. Journal of Geophysical Research Atmospheres,1999,1 041(C7):15 679-15 702

[6] LIU G, WANG H, SUN S, et al. Numerical Study on the Velocity Structure Around Tidal Fronts in the Yellow Sea[J]. Advances in Atmospheric Sciences,2003,20(3):453-460

[7] 張長寬.江蘇省近海海洋環(huán)境資源基本現(xiàn)狀[M].北京:海洋出版社,2013:333-340

[8] 王穎.黃海陸架輻射沙脊群[M].北京:中國環(huán)境科學(xué)出版社,2002:119-228

[9] WANG M, Ahn J H, JIANG L, et al. Ocean Color Products from the Korean Geostationary Ocean Color Imager (GOCI)[J]. Optics Express,2013,21(3):3 835

[10] 項二玲.關(guān)于水環(huán)境監(jiān)測實驗質(zhì)量分析[J].環(huán)境與生活, 2014(14):70

[11] 楊耀中,馮衛(wèi)兵.大范圍波浪傳播數(shù)學(xué)模型[J].水道港口, 2010,31(2):83-89

P237

B

1672-4623（2017）09-0037-04

10.3969/j.issn.1672-4623.2017.09.013

2016-12-16。

項目來源：國家海洋公益性行業(yè)科研專項資助項目（201005006-3）；江蘇省基礎(chǔ)研究計劃（自然科學(xué)基金）資助項目（BK2012414）。

金宇豪，博士研究生，主要研究方向為地圖學(xué)與地理信息系統(tǒng)。