楚曉亮張 杰王曙曜紀(jì)永剛王祎鳴
①(國(guó)家海洋局第一海洋研究所 青島 266061)
②(中國(guó)海洋大學(xué)信息科學(xué)與工程學(xué)院 青島 266100)
高頻地波雷達(dá)風(fēng)速直接反演的經(jīng)驗(yàn)?zāi)P?/p>
楚曉亮*①②張 杰①王曙曜②紀(jì)永剛①王祎鳴①
①(國(guó)家海洋局第一海洋研究所 青島 266061)
②(中國(guó)海洋大學(xué)信息科學(xué)與工程學(xué)院 青島 266100)
利用風(fēng)浪經(jīng)驗(yàn)?zāi)P蛷母哳l地波雷達(dá)(HFSWR)的回波譜數(shù)據(jù)中反演風(fēng)速需要有效波高等先驗(yàn)信息,因此風(fēng)速的反演受有效波高反演精度的影響。該文基于風(fēng)浪經(jīng)驗(yàn)?zāi)P?,利用風(fēng)速和高頻地波雷達(dá)海面回波二階譜與一階譜能量之比的關(guān)系,發(fā)展了無(wú)需波高信息的風(fēng)速直接反演的經(jīng)驗(yàn)?zāi)P?。將風(fēng)速反演經(jīng)驗(yàn)?zāi)P蛻?yīng)用到高頻地波雷達(dá)風(fēng)速的反演中,對(duì)兩部不同頻率的雷達(dá)在不同海域獲得的數(shù)據(jù)進(jìn)行了比較分析,結(jié)果表明,該文中采用的經(jīng)驗(yàn)?zāi)P湍軌蛴行У貙?duì)風(fēng)速進(jìn)行反演,其中三參數(shù)模型的結(jié)果略好于雙參數(shù)模型。
高頻地波雷達(dá);風(fēng)速反演;經(jīng)驗(yàn)?zāi)P?/p>
高頻地波雷達(dá)在海洋探測(cè)中具有大范圍、全天候及實(shí)時(shí)性的優(yōu)勢(shì),由于雷達(dá)的回波譜中含有豐富的海態(tài)信息,因此從回波譜中可以提取海面流場(chǎng)、風(fēng)場(chǎng)和海浪等信息。目前,利用高頻地波雷達(dá)進(jìn)行風(fēng)速反演時(shí),通常是先獲取有效波高后,再借助風(fēng)浪模型計(jì)算出來(lái)。有效波高的反演最初是基于文獻(xiàn)[1]以高頻電磁波海面雷達(dá)散射理論為基礎(chǔ)而建立的反演經(jīng)驗(yàn)?zāi)P?。利用此模型與實(shí)測(cè)數(shù)據(jù)擬合分析獲得相應(yīng)的擬合參數(shù),從而求得波高,此模型在工程上適用性較強(qiáng)。在此基礎(chǔ)上,文獻(xiàn)[2-4]對(duì)模型進(jìn)行了修正,文獻(xiàn)[5]采用這一經(jīng)驗(yàn)?zāi)P屠枚囝l雷達(dá)對(duì)有效波高進(jìn)行了反演。另外一些學(xué)者則提出了幾種通過(guò)反演出浪高譜而獲得有效波高的方法[6?12],通過(guò)這些方法獲得有效波高及相應(yīng)的峰波頻率后,利用風(fēng)浪充分發(fā)展條件下的風(fēng)浪模型即可求得風(fēng)速[13]。還有一些學(xué)者通過(guò)尋找有效波高與雷達(dá)回波多普勒譜特征值間的關(guān)系來(lái)獲得風(fēng)速,如文獻(xiàn)[14]在風(fēng)浪充分發(fā)展的條件下, 通過(guò)建立風(fēng)速與二階譜峰值頻率位置之間的關(guān)系來(lái)反演風(fēng)速,并通過(guò)實(shí)驗(yàn)數(shù)據(jù)驗(yàn)證了該方法的實(shí)用性,但這種方法對(duì)獲得二階譜峰值點(diǎn)位置的精確度要求較高。文獻(xiàn)[15]通過(guò)人工神經(jīng)網(wǎng)絡(luò)方法對(duì)風(fēng)速進(jìn)行了反演,取得了不錯(cuò)的結(jié)果。文獻(xiàn)[16]建立了海面風(fēng)速與有效波高的經(jīng)驗(yàn)?zāi)P停ㄟ^(guò)對(duì)實(shí)測(cè)數(shù)據(jù)的反演對(duì)比驗(yàn)證了該模型的穩(wěn)定性。在實(shí)際應(yīng)用中,經(jīng)驗(yàn)?zāi)P头椒ū容^簡(jiǎn)單實(shí)用,但在求解風(fēng)速時(shí)都增加了反演波高這一環(huán)節(jié),即在對(duì)高頻地波雷達(dá)反演風(fēng)速進(jìn)行定標(biāo)時(shí)需要先獲得有效波高的反演值或?qū)崪y(cè)值。本文根據(jù)風(fēng)浪經(jīng)驗(yàn)?zāi)P?,直接建立了二階譜與一階譜能量之比和風(fēng)速關(guān)系的雙參數(shù)經(jīng)驗(yàn)?zāi)P停瑴p少了反演有效波高的環(huán)節(jié)。在此基礎(chǔ)上,發(fā)展了風(fēng)速反演的三參數(shù)經(jīng)驗(yàn)?zāi)P?,通過(guò)高頻地波雷達(dá)反演風(fēng)速與風(fēng)速計(jì)數(shù)據(jù)的對(duì)比,表明雙參數(shù)模型和三參數(shù)經(jīng)驗(yàn)?zāi)P投寄芎芎玫剡M(jìn)行風(fēng)速的反演,并且三參數(shù)經(jīng)驗(yàn)?zāi)P头囱莸慕Y(jié)果略好于雙參數(shù)模型。
根據(jù)風(fēng)浪經(jīng)驗(yàn)?zāi)P?SMB)(SMB:是由Sverdrup, Munk and Bretschneider 3個(gè)人名的首字母組成),以風(fēng)浪為主導(dǎo)的海面風(fēng)速與波高滿足[13]
式中,g是重力加速度,Hs是有效波高,V10是海面上方10 m處風(fēng)速,fp是峰值頻率。式(1)不宜直接求解,只能采用迭代等方法。SMB關(guān)系式可以簡(jiǎn)化成的經(jīng)驗(yàn)公式為[17]
文獻(xiàn)[16]基于式(2)建立了海面風(fēng)速與有效波高的經(jīng)驗(yàn)?zāi)P?,風(fēng)速可表示為
式中,a和b是待定系數(shù)。通過(guò)經(jīng)驗(yàn)?zāi)P头囱莩鲇行Рǜ?,但是利用此模型進(jìn)行風(fēng)速反演在很大程度上依賴于海面有效波高的反演[16]。
有效波高的反演可根據(jù)文獻(xiàn)[1]的經(jīng)驗(yàn)?zāi)P颓蠼狻?/p>
式中,k0為雷達(dá)波數(shù),ω為多普勒頻率,σ(2)和σ(1)分別為二階和一階雷達(dá)散射截面,W(ω)表示與兩列海浪波矢之間耦合作用相關(guān)的權(quán)函數(shù)。式(4)可以改寫成
Rw為二階譜能量與一階譜能量的加權(quán)之比。由于加權(quán)系數(shù)只依賴于多普勒頻率,對(duì)于一定的區(qū)間可視為常數(shù),文獻(xiàn)[3]利用了二階譜能量與一階譜能量的無(wú)權(quán)重比值R,將式(5)簡(jiǎn)化為
文獻(xiàn)[2]則將式(5)推廣為
式中,a和b仍然為待定系數(shù)。
將式(6)代入式(3),可以得到高頻地波雷達(dá)海面回波二階譜能量與一階譜能量之比和風(fēng)速關(guān)系的雙參數(shù)模型,即
式(8)中參數(shù)a包含了雷達(dá)波數(shù)k0。由此可見(jiàn),對(duì)于在風(fēng)作用下充分成長(zhǎng)的海面來(lái)說(shuō),可以不用通過(guò)求解Hs,而直接建立比值R與風(fēng)速的關(guān)系,從而對(duì)風(fēng)速進(jìn)行反演求解。式(8)與有效波高計(jì)算式(7)形式上基本相同,在統(tǒng)計(jì)分析時(shí),需要對(duì)a, b參數(shù)進(jìn)行擬合。為使風(fēng)速V與R更加有效地?cái)M合,本文中對(duì)式(8)進(jìn)行改進(jìn),加入第3個(gè)參數(shù)c以加強(qiáng)擬合曲線的上、下偏移的調(diào)節(jié),形成三參數(shù)模型,即
本文采用中船重工鵬力(南京)大氣海洋系統(tǒng)有限公司研發(fā)的高頻地波雷達(dá)系統(tǒng)分別在福建(2013年10月份)和濱海(2013年2月份)獲得的實(shí)驗(yàn)數(shù)據(jù)對(duì)模型進(jìn)行驗(yàn)證。雷達(dá)系統(tǒng)的發(fā)射頻率分別為7.815 MHz和10.750 MHz,該雷達(dá)系統(tǒng)采用了收發(fā)分開共站體制,其工作帶寬為30 kHZ,發(fā)射天線為三元八木天線,接收天線為24元雙排陣型,反演風(fēng)浪時(shí)采用了數(shù)字波束合成技術(shù),具體的雷達(dá)系統(tǒng)參數(shù)見(jiàn)表1。
福建海域的高頻雷達(dá)數(shù)據(jù)是2013年10月份獲得的,風(fēng)速反演區(qū)域距離雷達(dá)站點(diǎn)75 km,與正北方向夾角為90o,頻率為7.815 MHz。將反演區(qū)域內(nèi)布放浮標(biāo)所提供的風(fēng)速數(shù)據(jù)作為比測(cè)數(shù)據(jù)。圖1給出了雷達(dá)反演風(fēng)速結(jié)果與浮標(biāo)數(shù)據(jù)風(fēng)速結(jié)果的散點(diǎn)圖,根據(jù)式(9)三參數(shù)模型,擬合得到a, b, c為46.67, 0.35和-15.29,相關(guān)系數(shù)為0.73,均方根誤差為1.73m/s。圖2給出了擬合后的雷達(dá)反演結(jié)果與浮標(biāo)數(shù)據(jù)的比較圖,從圖中看出二者結(jié)果總體符合較好。根據(jù)式(8)雙參數(shù)模型對(duì)雷達(dá)數(shù)據(jù)和浮標(biāo)數(shù)據(jù)擬合得到a和b分別為51.62和0.98,均方根誤差為1.82 m/s,相關(guān)系數(shù)為0.71,數(shù)據(jù)散點(diǎn)圖和比較結(jié)果如圖3和圖4所示,從擬合結(jié)果上看,式(9)模型反演結(jié)果要略好于式(8)模型。比較圖2和圖4,在浮標(biāo)風(fēng)速在18 m/s左右情況下,三參數(shù)模型要顯著好于雙參數(shù)模型。對(duì)于風(fēng)速較小(<6 m/s)的情況下,兩種模型得到的結(jié)果與實(shí)測(cè)結(jié)果差別都較大。由于這兩種模型都是風(fēng)浪經(jīng)驗(yàn)?zāi)P?,都需要滿足風(fēng)浪充分發(fā)展的海態(tài)條件。因此,風(fēng)速較小情況下,海浪無(wú)法得到足夠的能量,不能充分發(fā)展,因而得到的誤差較大。
表1 高頻地波雷達(dá)系統(tǒng)指標(biāo)
濱海海域采用的是從2013年02月01號(hào)8點(diǎn)到2013年02月18號(hào)23點(diǎn)的數(shù)據(jù),雷達(dá)頻率為10.75 MHz,選取參考點(diǎn)位置距離雷達(dá)站點(diǎn)大約30 km,與正北方向夾角為85o。風(fēng)速實(shí)測(cè)數(shù)據(jù)選取濱海岸基自動(dòng)氣象站風(fēng)速計(jì)同一時(shí)段的觀測(cè)數(shù)據(jù)。利用三參數(shù)模型得到的結(jié)果如圖5和圖6所示。圖5給出了風(fēng)速計(jì)與高頻地波雷達(dá)反演風(fēng)速的散點(diǎn)圖,擬合得到a=14.24, b=0.24, c=?4.09,相關(guān)系數(shù)為0.5,均方根誤差為1.99 m/s。相關(guān)系數(shù)較低,但均方根誤差可以接受。圖6給出了風(fēng)速計(jì)和高頻地波雷達(dá)風(fēng)速比較圖。從圖上看,雖然風(fēng)速計(jì)測(cè)量的風(fēng)速和高頻雷達(dá)反演的風(fēng)速符合不太理想,但還是能夠反映出風(fēng)速的變化趨勢(shì)。圖中部分對(duì)應(yīng)風(fēng)速計(jì)風(fēng)速較高的數(shù)據(jù),雷達(dá)反演結(jié)果誤差較大。這一方面由于風(fēng)速計(jì)距雷達(dá)反演區(qū)域約為30 km,會(huì)造成一定的誤差;另一方面,通過(guò)觀察雷達(dá)多普勒譜,發(fā)現(xiàn)相對(duì)應(yīng)的地波雷達(dá)二階譜較小,因此可以解釋為雖然風(fēng)速較大,但是由于海浪沒(méi)有充分發(fā)展,無(wú)法得到足夠的能量,從而導(dǎo)致地波雷達(dá)二階譜較小并且不穩(wěn)定,反演結(jié)果較差。另外,較大的誤差也有可能來(lái)自雷達(dá)系統(tǒng)或風(fēng)場(chǎng)的不均勻分布。
根據(jù)兩部雷達(dá)的數(shù)據(jù)得到的擬合結(jié)果如表2所示,綜合比較分析可以發(fā)現(xiàn):(1)雷達(dá)頻率較低(7.815 MHz),在風(fēng)浪充分成長(zhǎng)下,其反演高風(fēng)速較之低風(fēng)速要好;而對(duì)于雷達(dá)頻率較高時(shí)(10.750 MHz),反演低風(fēng)速較之低頻雷達(dá)要好一些,這符合高頻地波雷達(dá)的遙感機(jī)理;(2)綜合兩個(gè)頻率的雷達(dá)數(shù)據(jù)風(fēng)速反演結(jié)果,三參數(shù)模型反演的結(jié)果略好于雙參數(shù)模型的結(jié)果;(3)風(fēng)速經(jīng)驗(yàn)?zāi)P蛯?duì)于充分成長(zhǎng)的海浪,反演結(jié)果較好,而對(duì)于不充分成長(zhǎng)的海浪,對(duì)風(fēng)速反演結(jié)果誤差較大,有一定的局限性。
表2 不同模型的風(fēng)速擬合參數(shù)
圖1 福建雷達(dá)三參數(shù)模型反演風(fēng)速結(jié)果與浮標(biāo)數(shù)據(jù)散點(diǎn)圖
圖2 福建雷達(dá)三參數(shù)模型反演風(fēng)速結(jié)果與浮標(biāo)比較圖
圖3 福建雷達(dá)雙參數(shù)模型反演風(fēng)速結(jié)果與浮標(biāo)數(shù)據(jù)散點(diǎn)圖
圖4 福建雷達(dá)雙參數(shù)模型反演風(fēng)速結(jié)果與浮標(biāo)比較圖
圖5 濱海雷達(dá)三參數(shù)模型反演風(fēng)速結(jié)果與風(fēng)速計(jì)數(shù)據(jù)散點(diǎn)圖
圖6 濱海雷達(dá)三參數(shù)模型反演風(fēng)速結(jié)果與風(fēng)速計(jì)比較圖
本文在風(fēng)浪經(jīng)驗(yàn)?zāi)P偷幕A(chǔ)上,建立了風(fēng)速和高頻地波雷達(dá)海面回波二階譜能量與一階譜能量之比的關(guān)系,發(fā)展了無(wú)需波高信息的高頻地波雷達(dá)風(fēng)速直接反演的經(jīng)驗(yàn)?zāi)P?。反演風(fēng)速經(jīng)驗(yàn)?zāi)P瓦m合以風(fēng)浪為主導(dǎo),并充分發(fā)展的海面。文中采用兩部不同頻率的雷達(dá)在不同海域的數(shù)據(jù)進(jìn)行了比較分析,結(jié)果表明,文中采用的經(jīng)驗(yàn)?zāi)P湍軌蛴行У貙?duì)風(fēng)速進(jìn)行反演,并且三參數(shù)模型反演的結(jié)果略好雙參數(shù)模型。不過(guò)從反演結(jié)果上看,對(duì)于風(fēng)浪還沒(méi)有充分成長(zhǎng)的海面,該模型應(yīng)用受到限制。另外,本文共采用了48天的數(shù)據(jù),數(shù)據(jù)量相對(duì)較少,三參數(shù)模型相對(duì)于雙參數(shù)模型的優(yōu)勢(shì)有待于獲取更多數(shù)據(jù)進(jìn)行分析。
致謝:感謝中船重工鵬力(南京)大氣海洋信息系統(tǒng)有限公司周濤研究員提供的高頻地波雷達(dá)數(shù)據(jù)及比測(cè)數(shù)據(jù)。
[1] Barrick D and Weber B. On the nonlinear theory for gravity waves on the ocean's surface. Part II: interpretation and applications[J]. Journal of Physical Oceanography, 1977, 7(1): 11-21.
[2] Maresca J W Jr and Georges T. Measuring rms wave height and the scalar ocean wave spectrum with HF skywave radar[J]. Journal of Geophysical Research, 1980, 85(C5): 2759-2771.
[3] Heron S and Heron M. A comparison of algorithms for extracting significant wave height from HF radar ocean backscatter spectra[J]. Journal of Atmospheric and Oceanic Technology, 1998, 15(5): 1157-1163.
[4] 吳雄斌, 李倫, 李炎. 高頻地波雷達(dá)海面有效波高探測(cè)實(shí)驗(yàn)研究[J]. 海洋與湖沼, 2012, 43(2): 210-216. Wu Xiong-bin, Li Lun, Li Yan, et al.. Experimental research on significant wave height detecting with HFSWR OSMAR071[J]. Oceanologia et Limnologia Sinica, 2012, 43(2): 210-216.
[5] Zhao Chen, Chen Ze-zong, Jiang Yan-ni, et al.. Exploration and validation of wave-height measurement using multifrequen HF radar[J]. Journal of Atmospheric and Oceanic Technology, 2013, 30(9): 2189-2202.
[6] Lipa B and Barrick D. Extraction of sea state from HF radar sea echo: mathematical theory and modeling[J]. Radio Science, 1986, 21(1): 81-100.
[7] Wyatt L. The measurement of the ocean wave directional spectrum from HF radar Doppler spectra[J]. Radio Science, 1986, 21(3): 473-485.
[8] Howell R and Walsh J. Measurement of ocean wave spectra using narrow beam HF radar[J]. IEEE Journal of Oceanic Engineering, 1993, 18(3): 296-305.
[9] Hisaki Y. Nonlinear inversion of the integral equation to estimate ocean wave spectra from HF radar[J]. Radio Science, 1996, 31(1): 25-29.
[10] Green J and Wyatt L. Row-action inversion of the Barrick-Weber equations[J]. Journal of Atmospheric and Oceanic Technology, 2006, 23(3): 501-510.
[11] 李倫, 吳雄斌, 龍超, 等. 基于正則化方法的高頻地波雷達(dá)海浪方向譜反演[J]. 地球物理學(xué)報(bào), 2013, 56(1): 219-229.
Li Lun,Wu Xion-bin, Long Chao, et al.. Regularization inversion method for extracting ocean wave spectra from HFSWR sea echo[J]. Chinese Journal of Geophysics, 2013, 56(1): 219-229.
[12] 李倫, 吳雄斌, 徐興安. 基于優(yōu)化理論的高頻地波雷達(dá)海浪參數(shù)反演[J], 物理學(xué)報(bào), 2014, 63(3): 441-451.
Li Lun, Wu Xiong-bin, and Xu Xing-an. Optimization method for extracting ocean wave parameters from HFSWR[J]. Acta Physica Sinica, 2014, 63(3): 441-451.
[13] 文必洋, 黃為民, 王小華. OSMAR2000探測(cè)海面風(fēng)浪場(chǎng)原理與實(shí)現(xiàn)[J]. 武漢大學(xué)學(xué)報(bào)(理學(xué)版), 2001, 47(5): 642-644. Wen Bi-yang, Huang Wei-min, and Wang Xiao-hua. The principle and realization of sounding ocean's wave height and wind with OSMA2000[J]. Journal of Wuhan University (Nature Science Edition), 2001, 47(5): 642-644.
[14] Green D, Gill E, and Huang W. An inversion method for extraction of wind speed from high-frequency ground-wave radar oceanic backscatter[J]. IEEE Transactions on Geoscience and Remote Sensing, 2009, 47(10): 3338-3346.
[15] Shen W, Gurgel K, and Voulgaris G. Wind-speed inversion from HF radar first-order backscatter signal[J]. Ocean Dynamics, 2012, 62: 105-121.
[16] 李倫, 吳雄斌, 徐興安, 等. 高頻地波雷達(dá)風(fēng)速反演經(jīng)驗(yàn)?zāi)P蚚J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2012, 37(9): 1096-1098.
Li Lun, Wu Xiong-bin, Xu Xing-an, et al.. An empirical model for wind speed inversion by HFSWR[J]. Geomatics and Information Science of Wuhan University, 2012, 37(9): 1096-1098.
[17] 文圣常, 余宙文. 海浪理論與計(jì)算原理[M]. 北京: 科學(xué)出版社, 1984: 469-473.
Wen Sheng-chang and Yu Zhou-wen. Ocean Wave Theory and Calculation Principle[M]. Beijing: Science Press, 1984: 469-473.
楚曉亮: 男,1977年生,博士,講師,研究方向?yàn)閄波段導(dǎo)航雷達(dá)和高頻地波雷達(dá)海態(tài)反演.
張 杰: 男,1963年生,博士,研究員,研究方向?yàn)楹Q筮b感.
王曙曜: 男,1988年生,碩士生,研究方向?yàn)楦哳l地波雷達(dá)海洋參數(shù)反演.
紀(jì)永剛: 男,1977年生,博士,副研究員,研究方向?yàn)楦哳l地波雷達(dá)海態(tài)反演及目標(biāo)探測(cè).
王祎鳴: 男,1977年生,碩士,助理研究員,研究方向?yàn)楦哳l地波雷達(dá)信號(hào)處理及雜波抑制.
An Empirical Model for Wind Speed Inversion Directly from High Frequency Surface Wave Radar Sea Echo
Chu Xiao-liang①②Zhang Jie①Wang Shu-yao②Ji Yong-gang①Wang Yi-ming①
①(First Institute of Oceanography of National Bureau of Oceanography, Qingdao 266061, China)
②(College of Information Science and Engineering, Ocean University of China, Qingdao 266100, China)
The information of the significant wave height is needed in the wind speed inversion from the sea echo of High Frequency Surface Wave Radar (HFSWR) by using the empirical model of wind waves. Therefore, the accuracy of the significant wave height has an effect on the wind speed inversion. Based on the empirical model of wind waves, a wind speed inversion empirical model for wind retrieval without the information of wave height is developed, which uses the relationship between the wind speed and ratio of second-order to first-order spectrum energies. The inversion model is applied to the wind speed extraction. And the data obtained from two radars with different frequency in different detecting area are analyzed. The results show that the proposed model can be used to extract the wind speed from the HFSWR sea echo and the results of which three-parameter model are better than the two-parameter model.
High Frequency Surface Wave Radar (HFSWR); Wind speed invesion; Empirical model
TN958
: A
:1009-5896(2015)04-1013-04
10.11999/JEIT140850
2014-06-26收到,2015-01-06改回
國(guó)家自然科學(xué)基金重點(diǎn)項(xiàng)目(61032011)和中國(guó)博士后科學(xué)基金(2013M531559)資助課題
*通信作者:楚曉亮 xlchu@ouc.edu.cn