王曉晨，萬劍華，張振華，孟俊敏，范陳清

(1. 中國石油大學(xué)(華東) 地球科學(xué)與技術(shù)學(xué)院，山東 青島 266580；2. 國家海洋局第一海洋研究所，山東 青島 266061；3. 北京遙測技術(shù)研究所，北京 100076)

?

小入射角下海浪后向散射系數(shù)不對(duì)稱性和各向異性分析驗(yàn)證

王曉晨1，2，萬劍華1，張振華3，孟俊敏2*，范陳清2

(1. 中國石油大學(xué)(華東) 地球科學(xué)與技術(shù)學(xué)院，山東 青島 266580；2. 國家海洋局第一海洋研究所，山東 青島 266061；3. 北京遙測技術(shù)研究所，北京 100076)

基于北京遙測技術(shù)研究所機(jī)載波譜儀實(shí)測數(shù)據(jù)以及ECMWF提供的風(fēng)場數(shù)據(jù)，分析了不同風(fēng)速條件下小入射角海浪平均后向散射系數(shù)不對(duì)稱性與各向異性隨入射角的變化特性，并對(duì)導(dǎo)致這一現(xiàn)象的原因進(jìn)行了初步分析。結(jié)果表明：小入射角下，海浪平均后向散射系數(shù)隨入射角增大呈現(xiàn)遞減趨勢并在18°時(shí)達(dá)到某一穩(wěn)定值，之后基本維持不變；海浪后向散射系數(shù)不對(duì)稱性與各向異性均隨入射角增大呈現(xiàn)遞增趨勢，同樣在18°增大至某一穩(wěn)定值，之后基本維持不變；低風(fēng)速條件下，海浪不對(duì)稱性和各向異性與風(fēng)速有關(guān)，且風(fēng)速越大不對(duì)稱性與各向異性越明顯。

機(jī)載波譜儀；海浪；小入射角；不對(duì)稱性；各向異性

1 引言

海浪作為一種復(fù)雜的風(fēng)生重力波，按照成因可以分為重力波、毛細(xì)重力波和毛細(xì)波，其中重力波波長從幾十厘米到幾百米，主要是重力作用主導(dǎo)的大尺度波；毛細(xì)重力波波長約為幾厘米，是由表面張力和重力混合作用，毛細(xì)波波長為毫米級(jí)，主要是表面張力作用，二者都為小尺度波。重力波描述的是宏觀大尺度范圍內(nèi)海浪的分布特性，而毛細(xì)重力波和毛細(xì)波描述的是微觀小尺度空間內(nèi)海浪微細(xì)結(jié)構(gòu)，研究海浪空間分布特性以及微細(xì)結(jié)構(gòu)，對(duì)于理解海浪風(fēng)成機(jī)制和內(nèi)部動(dòng)力過程，提高人類對(duì)于海洋的認(rèn)識(shí)具有重要意義。海浪生成機(jī)制主要基于共振和平行流不穩(wěn)定性理論，前者適用于海浪生成階段，后者主要適用于海浪發(fā)展階段[1—2]，由于波與波之間的相互作用，海浪中顯著波峰附近波面前側(cè)出現(xiàn)短重力波和毛細(xì)重力波，Gent與Talyor認(rèn)為這些短波的分布增大了海表面的粗糙度改變了海面能量分布[3]，海浪長波的存在又表現(xiàn)在對(duì)短波的調(diào)制，因而可以通過探測海面后向散射系數(shù)得到海面調(diào)制長波以及短波的分布。目前，通過微波手段探測海浪主要包括合成孔徑雷達(dá)(SAR)、高度計(jì)、散射計(jì)、波譜儀等。其中SAR采用大入射角斜視工作方式，散射計(jì)采用大入射角360°掃描工作方式，波譜儀采用小入射角360°掃描工作方式，近年來由于波譜儀具有波數(shù)截?cái)嘈?yīng)小、散射模型簡單便于計(jì)算、不需要額外的風(fēng)場數(shù)據(jù)輸入等優(yōu)勢[4]，逐漸成為海浪探測的重要手段。

關(guān)于海浪后向散射系數(shù)特性的研究，Simon在2000年利用雙極化Ku波段機(jī)載散射計(jì)數(shù)據(jù)分析了HH、VV極化下和隨入射角變化特性，并認(rèn)為在中等以上入射角，Bragg散射會(huì)導(dǎo)致正的不對(duì)稱性，而小入射角下的幾何光學(xué)散射會(huì)導(dǎo)致負(fù)的不對(duì)稱性[5]。Mouche等在2006年利用STORM(Système de Télédétection pour l’ Observation Radar de la Mer)數(shù)據(jù)，通過引入波與波之間非線性交互作用，分析了大入射角下短重力波和毛細(xì)重力波對(duì)雷達(dá)信號(hào)各向異性的影響[6]。儲(chǔ)小青在2011年利用Ku波段TRMM(Tropical Rainfall Mapping Mission)降雨雷達(dá)數(shù)據(jù)分析了小入射角下和的分布特性[2]。然而由于現(xiàn)有遙感手段工作模式限制，目前國內(nèi)外學(xué)者對(duì)于海浪后向散射系數(shù)的研究多為大入射角下情況下，波譜儀作為唯一一種在小入射角下進(jìn)行方位向360°掃描海浪的設(shè)備，為分析小入射角下海浪后向散射系數(shù)提供了數(shù)據(jù)支持。

本文利用機(jī)載波譜儀實(shí)測數(shù)據(jù)分析了小入射角下海浪后向散射系數(shù)不對(duì)稱性和各向異性。

2 實(shí)驗(yàn)數(shù)據(jù)

本文所選擇的數(shù)據(jù)為北京遙測技術(shù)研究所機(jī)載波譜儀掛飛實(shí)驗(yàn)數(shù)據(jù)和ECMWF提供的風(fēng)場數(shù)據(jù)。

本次實(shí)驗(yàn)機(jī)載波譜儀是由北京遙測技術(shù)研究所研制的一種脈沖雷達(dá)，其工作原理與STORM基本相同，均是在小角度下通過旋轉(zhuǎn)天線進(jìn)行360°掃描。不同的是，北京遙測研究所機(jī)載波譜儀采用脈沖技術(shù)，調(diào)頻連續(xù)波技術(shù)需要兩根天線用于發(fā)射接收信號(hào)，故系統(tǒng)設(shè)計(jì)較復(fù)雜，體積大，而脈沖技術(shù)可以克服以上缺點(diǎn)。另外，二者具體的雷達(dá)參數(shù)也有差異，北京遙測研究所機(jī)載波譜儀雷達(dá)參數(shù)如下表1所示。

表1 北京遙測研究所機(jī)載波譜儀雷達(dá)參數(shù)

北京遙測技術(shù)研究所波譜儀掛飛實(shí)驗(yàn)于2014年6月在黃海海域進(jìn)行，期間飛行了多個(gè)架次，分別獲取了海面回波功率、航向、天線方位角、飛機(jī)姿態(tài)以及GPS時(shí)間等數(shù)據(jù)，現(xiàn)場同步布放了波浪浮標(biāo)進(jìn)行驗(yàn)證，浮標(biāo)位置為36°37′18″N，122°43′40″E，本文選擇兩個(gè)架次實(shí)驗(yàn)數(shù)據(jù)，對(duì)應(yīng)實(shí)驗(yàn)區(qū)域與浮標(biāo)具體位置如圖1。

圖1 航飛試驗(yàn)區(qū)域以及浮標(biāo)位置Fig.1 Experimental flight and buoy location

歐洲中期天氣預(yù)報(bào)中心(European Centre for Medium-Range Weather Forecasts，ECMWF)為全球用戶提供的各類再分析數(shù)據(jù)，均可免費(fèi)下載。其中本文用到的風(fēng)場數(shù)據(jù)其空間分辨率為0.125°×0.125°，時(shí)間分辨率為6 h。

3 數(shù)據(jù)處理與結(jié)果分析

圖2 回波功率隨入射角變化Fig.2 Echo power as functions of incidence angles

本文對(duì)機(jī)載波譜儀回波信號(hào)進(jìn)行預(yù)處理得到海浪后向散射系數(shù)，入射角選取4°～18°，得到某方位向與入射角關(guān)系如圖2所示。5°～12°海浪后向散射系數(shù)隨著入射角增大呈線性減小，15°以后基本穩(wěn)定在某一定值。這是由于天線采用小角度入射，電磁波與海面作用可視為準(zhǔn)鏡面散射，海面回波較強(qiáng)，然而隨著入射角的增大，布拉格散射的影響越來越顯著，電磁波與海面的作用逐漸被布拉格散射主導(dǎo)。

由于波譜儀采用小入射角下旋轉(zhuǎn)掃描的工作方式，因而可以觀測到0°～360°整個(gè)方位向的分布，故可以用來分析海浪后向散射系數(shù)方位向不對(duì)稱性和各向異性，本文統(tǒng)計(jì)了其中一圈數(shù)據(jù)4°～18°不同入射角下0°～360°方位向分布，如圖3所示。

圖3 不同入射角下方位向分布，a-h分別代表入射角4°、6°、8°、10°、12°、14°、16°、18°Fig.3 As functions of azimuth in different incidence angles，a-h incidence angles equal 4°、6°、8°、10°、12°、14°、16°、18°

由圖3a所示，入射角為4°時(shí)，并沒有出現(xiàn)明顯的雙峰現(xiàn)象，說明此時(shí)σ0在方位向的分布并沒有明顯的風(fēng)向信息，這與儲(chǔ)小青等的結(jié)論基本一致[7]，與Tran和Chapron的結(jié)論不一致[8]，Tran和Chapron通過分析高度計(jì)Jason-1數(shù)據(jù)認(rèn)為在天底點(diǎn)存在較明顯的方向信號(hào)。當(dāng)入射角超過4°后，σ0方位向分布均出現(xiàn)明顯雙峰趨勢，其中出現(xiàn)雙峰位置方位向?yàn)轫橈L(fēng)、逆風(fēng)向，出現(xiàn)波谷位置方位向?yàn)閭?cè)風(fēng)向，說明海浪后向散射系數(shù)σ0存在明顯的各向異性。由ECMWF提供的風(fēng)向約為135°，實(shí)際風(fēng)向約為110°，可見順逆風(fēng)方向基本都是在波峰位置，這與散射計(jì)觀測結(jié)果以及儲(chǔ)小青等的結(jié)論一致，與Hesany等較不一致[9]，Hesany認(rèn)為回波信號(hào)峰值并不總出現(xiàn)在順逆風(fēng)風(fēng)向，二者會(huì)有一定偏差，當(dāng)入射角位于6°～12°之間時(shí)，如圖3b～d，此時(shí)順風(fēng)向σ0大于逆風(fēng)向σ0，當(dāng)入射角為12°時(shí)，如圖3e，此時(shí)順風(fēng)向和逆風(fēng)向σ0基本相近，當(dāng)入射角超過12°后，如圖3f～h，此時(shí)逆風(fēng)向σ0大于順風(fēng)向，這與散射計(jì)觀測結(jié)果基本一致[5]，兩個(gè)波峰峰值存在一定偏差，說明σ0存在明顯的海浪后向散射系數(shù)不對(duì)稱性，而且當(dāng)入射角小于12°時(shí)，σ0呈現(xiàn)負(fù)的不對(duì)稱性，大于12°時(shí)呈現(xiàn)正的不對(duì)稱性，這與Simon結(jié)論基本一致[5]。

為了更好地分析海浪后向散射系數(shù)特性，對(duì)σ0使用余弦函數(shù)展開：

(1)

3.1 平均后向散射系數(shù)

平均后向散射系數(shù)指海浪后向散射系數(shù)去除方位向特性后的結(jié)果，其表示海浪后向散射系數(shù)在整個(gè)方位向的平均值。依照式(1)得到不同風(fēng)速下平均后向散射系數(shù)A0隨入射角變化如圖4所示。

圖4 不同風(fēng)速下后向散射系數(shù)A0隨入射角變化Fig.4 A0 as functions of incidence angles in different wind speed

由圖可見，海浪平均后向散射系數(shù)值A(chǔ)0受入射角影響較明顯，A0隨入射角增大逐漸減小且減小趨勢越來越不明顯，當(dāng)入射角接近18°時(shí)，A0基本穩(wěn)定至某一定值，這主要是因?yàn)樵诮斓c(diǎn)下，準(zhǔn)鏡面反射占主導(dǎo)地位，因而回波值最高，隨著入射角的增大，布拉格散射的作用越來越明顯，當(dāng)入射角大于18°時(shí)，布拉格散射開始占主導(dǎo)地位，此時(shí)回波值基本為布拉格散射。而在3～5 m/s低風(fēng)速條件下，A0基本不隨風(fēng)速發(fā)生較大變化，這是因?yàn)轱L(fēng)對(duì)海面的影響主要是改變其粗糙度，即改變海面小尺度毛細(xì)重力波和短重力波的分布，總而言之，風(fēng)的作用僅僅是對(duì)海浪后向散射系數(shù)方位向特性的影響，因而去除了方位向特性的平均后向散射系數(shù)基本不隨風(fēng)速發(fā)生變化，風(fēng)對(duì)海面更多的是對(duì)海浪后向散射系數(shù)方位向特性(包括不對(duì)稱性和各向異性)的影響。

3.2 不對(duì)稱性

不對(duì)稱性指逆風(fēng)方向與順風(fēng)方向海浪后向散射系數(shù)值不同，二者不成對(duì)稱分布的現(xiàn)象，是各向異性的一種特殊現(xiàn)象。其產(chǎn)生的原因主要是由于風(fēng)吹拂海面時(shí)順風(fēng)向海面比較光滑，而逆風(fēng)向海面比較粗糙，即逆風(fēng)向海面小尺度毛細(xì)重力波和短重力波分布更多，在海浪后向散射系數(shù)數(shù)值上表示為順風(fēng)逆風(fēng)出現(xiàn)差異，即不對(duì)稱性。依照式(1)得到不同風(fēng)速下不對(duì)稱性A1隨入射角變化如圖5所示。

圖5 不同風(fēng)速下不對(duì)稱性A1隨入射角變化Fig.5 A1 as functions of incidence angles in different wind speed

由圖5可見，海浪不對(duì)稱性A1受風(fēng)速和入射角影響，入射角在5°～12°之間時(shí)，海浪不對(duì)稱性A1基本不隨風(fēng)速變化，當(dāng)入射角超過12°之后，風(fēng)速對(duì)A1的作用越來越明顯，其中風(fēng)速越大，平均后向散射系數(shù)A1越大。這與儲(chǔ)小青的結(jié)論并不一致[2]，儲(chǔ)小青認(rèn)為在低速風(fēng)條件下，不對(duì)稱性A1隨入射角呈遞減趨勢且全為負(fù)值，這可能與降雨雷達(dá)數(shù)據(jù)以及海況有關(guān)，降雨雷達(dá)數(shù)據(jù)僅僅覆蓋了海浪后向散射系數(shù)負(fù)不對(duì)稱性部分，沒有對(duì)正不對(duì)稱性部分進(jìn)行比較。對(duì)于某一風(fēng)速，A1隨著入射角增大呈現(xiàn)緩慢增長的趨勢，且增長趨勢越來越明顯，不對(duì)稱性越明顯，這說明增大入射角可以更好地觀測小尺度海浪騎行毛細(xì)重力波和短重力波在海表面的分布。

3.3 各向異性

各向異性指海浪后向散射系數(shù)在不同方位向下其值不同的現(xiàn)象，其產(chǎn)生的原因主要是由于風(fēng)吹拂海面，使得不同方位向下海面粗糙度出現(xiàn)變化，從而影響小尺度毛細(xì)重力波和短重力波概率密度分布函數(shù)，在海浪后向散射系數(shù)數(shù)值上表示為不同方位向其值出現(xiàn)差異，即各向異性。依照式(1)得到不同風(fēng)速下各向異性A2隨入射角變化如圖6所示。

圖6 不同風(fēng)速下各向異性A2隨入射角變化Fig.6 A2 as functions of incidence angles in different wind speed

由圖6可見，海浪各向異性A2同樣受風(fēng)速和入射角影響，其中風(fēng)速越大，后向散射系數(shù)A2越大。對(duì)于某一風(fēng)速，當(dāng)入射角在5°～12°之間時(shí)，A2隨著入射角增大呈現(xiàn)緩慢增長的趨勢，當(dāng)入射角大于12°時(shí)，增長趨勢越來越不明顯，入射角增大至18°時(shí)，A2逐漸趨近一定值，此時(shí)各向異性最明顯。這說明增大入射角可以更好地觀測小尺度海浪騎行毛細(xì)重力波和短重力波在海表面的分布。

4 結(jié)論與總結(jié)

本文利用機(jī)載波譜儀掛飛試驗(yàn)數(shù)據(jù)，通過分析小入射角下海浪后向散射系數(shù)特性，得到以下結(jié)論：

(1)在低風(fēng)速條件下，海浪平均后向散射系數(shù)值基本不隨風(fēng)速變化，對(duì)于某一風(fēng)速條件下，散射值隨入射角增大逐漸減小直至某一定值，此后不再隨入射角變化。

(2)在低風(fēng)速條件下，風(fēng)速越大，海浪后向散射系數(shù)不對(duì)稱性越大，并且對(duì)于某一風(fēng)速條件下，不對(duì)稱性隨著入射角增大呈現(xiàn)緩慢增長的趨勢，隨著入射角變大增長趨勢越來越不明顯，最后穩(wěn)定在一定值。

(3)在低風(fēng)速條件下，風(fēng)速越大，海浪后向散射系數(shù)各向異性越大，并且對(duì)于某一風(fēng)速條件下，各向異性隨著入射角增大呈現(xiàn)緩慢增長的趨勢，且增長趨勢越來越明顯。

本文利用北京遙測技術(shù)研究所機(jī)載波譜儀掛飛實(shí)驗(yàn)數(shù)據(jù)，分析了小入射角下海浪后向散射系數(shù)特性，比較了不同風(fēng)速下海浪平均后向散射系數(shù)、不對(duì)稱性和各向異性，由于實(shí)驗(yàn)數(shù)據(jù)所限，本文并沒有針對(duì)中、高風(fēng)速條件下海浪散射系數(shù)特性分析，并且由于本次實(shí)驗(yàn)波譜儀沒有進(jìn)行定標(biāo)，因而只分析了海浪后向散射系數(shù)特性的相對(duì)變化，下一步有望針對(duì)中、高風(fēng)速條件下，分析海浪不對(duì)稱性以及各向異性，并分析海浪破碎對(duì)二者的影響。

致謝：感謝北京遙測技術(shù)研究所提供機(jī)載波譜儀數(shù)據(jù)集以及ECMWF提供的風(fēng)場數(shù)據(jù)。

[1]文圣常. 海浪理論與計(jì)算原理[M]. 北京：科學(xué)出版社， 1984．

Wen Shengchang. The Wave Theory and Calculation Principle[M]. Beijing：Science Press， 1984．

[2]儲(chǔ)小青. 海浪波譜儀海浪遙感方法及應(yīng)用基礎(chǔ)研究[D]. 青島：中國科學(xué)院海洋研究所， 2011．

Chu Xiaoqing. The algorithm and applied basic study on remote sensing of ocean wave spectrum for ocean wave spectrometer[D]. Qingdao：Institute of Oceanology，Chinese Academy of Science， 2011．

[3]Gent P， Taylor P A. A numerical model of the air flow above water waves[J]. Journal of Fluid Mechanics， 1976， 77(77):105-128.

[4]Hauser D， Caudal G， Rijckenberg G J， et al. RESSAC: a new airborne FM/CW radar ocean wave spectrometer[J]. IEEE Transactions on Geoscience & Remote Sensing， 1992， 30(5):981-995．

[5]Yueh S H， West R， Li F K， et al. Dual-polarized Ku-band backscatter signatures of hurricane ocean winds[J]. IEEE Transactions on Geoscience & Remote Sensing， 2000， 38(1):73-88．

[6]Mouche A， Hauser D， Caudal G， et al. Use of Dual Polarization Radar Measurements to Understand the Azimuth Behavior of the Sea Surface Backscattered Signal[C]// Geoscience and Remote Sensing Symposium. Denver, Colorado, 2006：1867-1870.

[7]Chu Xiaoqing, He Yijun, Chen Gengxin. Asymmetry and anisotropy of microwave backscatter at low incidence angles[J]. IEEE Transactions on Geoscience & Remote Sensing， 2012， 50(10):4014-4024．

[8]Tran N， Chapron B. Combined wind vector and sea state impact on ocean nadir-viewing Ku-and C-band radar cross-sections[J]. Sensors， 2006， 6(3):193-207．

[9]Hesany V， Plant W J， Keller W C. The normalized radar cross section of the sea at 10° incidence[J]. Geoscience & Remote Sensing IEEE Transactions on， 2000， 38(1):64-72.

[10]Caudal G. A physical model for the narrowing of the directional sea wave spectra in the short gravity to gravity-capillary range[J]. Journal of Geophysical Research Atmospheres， 2002， 107(C10):12-1-12-9．

Analysis and validation of wave backscatter coefficients asymmetry and anisotropy at low incidence angle

Wang Xiaochen1，2，Wan Jianhua1，Zhang Zhenhua3， Meng Junmin2，F(xiàn)an Chenqing2

(1.SchoolofGeosciences，ChinaUniversityofPetroleum，Qingdao266580，China; 2.TheFirstInstituteofOceanography，StateOceanicAdministration，Qingdao266061，China; 3.ChinaAcademyofSpaceTechnology，Beijing100076，China)

Based on the China Academy of Space Technology airborne spectrometer test data and the ECMWF wind data， the average backscattering cross sections， asymmetry and anisotropy changed with the incidence angle under different wind speed conditions， and the reason of wave backscatter asymmetry and anisotropy were preliminary analyzed. The results showed that at low incidence angle， the average backscattering cross section decreased as the incidence angle increased and reached a stable value at 18°，after essentially unchanged；wave backscatter asymmetry and anisotropy increased as incidence angle increased， reached a stable value at 18° ， after essentially unchanged; under small winds speed conditions， waves asymmetry and anisotropy were related with wind speed， and the anisotropy and asymmetry were more obvious as wind speed increased．

airborne spectrometer; wave; low incidence angle; asymmetry; anisotropy

2015-11-06；

2016-03-27。

國家863計(jì)劃(SS2013AA091206)；國家自然科學(xué)基金青年基金(2013AA09A505，41406197)。

王曉晨(1990—)，男，山東省濰坊市人，主要從事波譜儀海浪反演研究。E-mail:15863014591@163.com

孟俊敏，男，研究員，主要從事海洋微波動(dòng)力過程研究。E-mail:mengjm@fio.org.cn

P731.22

A

0253-4193(2016)11-0043-06

王曉晨，萬劍華，張振華，等. 小入射角下海浪后向散射系數(shù)不對(duì)稱性和各向異性分析驗(yàn)證[J]. 海洋學(xué)報(bào)， 2016， 38(11): 43-48， doi:10.3969/j.issn.0253-4193.2016.11.004

Wang Xiaochen， Wan Jianhua， Zhang Zhenhua， et al. Analysis and validation of wave backscatter coefficients asymmetry and anisotropy at low incidence angle[J]. Haiyang Xuebao， 2016， 38(11): 43-48， doi：10.3969/j.issn.0253-4193.2016.11.004