高頻地波雷達風速直接反演的經(jīng)驗?zāi)Ｐ?/h1>
2015-07-12 13:59:59楚曉亮張杰王曙曜紀永剛王祎鳴
電子與信息學(xué)報 2015年4期
楚曉亮張 杰王曙曜紀永剛王祎鳴
①(國家海洋局第一海洋研究所 青島 266061)
②(中國海洋大學(xué)信息科學(xué)與工程學(xué)院 青島 266100)
高頻地波雷達風速直接反演的經(jīng)驗?zāi)Ｐ?/p>
楚曉亮*①②張 杰①王曙曜②紀永剛①王祎鳴①
①(國家海洋局第一海洋研究所 青島 266061)
②(中國海洋大學(xué)信息科學(xué)與工程學(xué)院 青島 266100)
利用風浪經(jīng)驗?zāi)Ｐ蛷母哳l地波雷達(HFSWR)的回波譜數(shù)據(jù)中反演風速需要有效波高等先驗信息，因此風速的反演受有效波高反演精度的影響。該文基于風浪經(jīng)驗?zāi)Ｐ?，利用風速和高頻地波雷達海面回波二階譜與一階譜能量之比的關(guān)系，發(fā)展了無需波高信息的風速直接反演的經(jīng)驗?zāi)Ｐ汀L速反演經(jīng)驗?zāi)Ｐ蛻?yīng)用到高頻地波雷達風速的反演中，對兩部不同頻率的雷達在不同海域獲得的數(shù)據(jù)進行了比較分析，結(jié)果表明，該文中采用的經(jīng)驗?zāi)Ｐ湍軌蛴行У貙︼L速進行反演，其中三參數(shù)模型的結(jié)果略好于雙參數(shù)模型。
高頻地波雷達；風速反演；經(jīng)驗?zāi)Ｐ?/p>
1 引言
高頻地波雷達在海洋探測中具有大范圍、全天候及實時性的優(yōu)勢，由于雷達的回波譜中含有豐富的海態(tài)信息，因此從回波譜中可以提取海面流場、風場和海浪等信息。目前，利用高頻地波雷達進行風速反演時，通常是先獲取有效波高后，再借助風浪模型計算出來。有效波高的反演最初是基于文獻[1]以高頻電磁波海面雷達散射理論為基礎(chǔ)而建立的反演經(jīng)驗?zāi)Ｐ汀＠么四Ｐ团c實測數(shù)據(jù)擬合分析獲得相應(yīng)的擬合參數(shù)，從而求得波高，此模型在工程上適用性較強。在此基礎(chǔ)上，文獻[2-4]對模型進行了修正，文獻[5]采用這一經(jīng)驗?zāi)Ｐ屠枚囝l雷達對有效波高進行了反演。另外一些學(xué)者則提出了幾種通過反演出浪高譜而獲得有效波高的方法[6?12]，通過這些方法獲得有效波高及相應(yīng)的峰波頻率后，利用風浪充分發(fā)展條件下的風浪模型即可求得風速[13]。還有一些學(xué)者通過尋找有效波高與雷達回波多普勒譜特征值間的關(guān)系來獲得風速，如文獻[14]在風浪充分發(fā)展的條件下， 通過建立風速與二階譜峰值頻率位置之間的關(guān)系來反演風速，并通過實驗數(shù)據(jù)驗證了該方法的實用性，但這種方法對獲得二階譜峰值點位置的精確度要求較高。文獻[15]通過人工神經(jīng)網(wǎng)絡(luò)方法對風速進行了反演，取得了不錯的結(jié)果。文獻[16]建立了海面風速與有效波高的經(jīng)驗?zāi)Ｐ?，通過對實測數(shù)據(jù)的反演對比驗證了該模型的穩(wěn)定性。在實際應(yīng)用中，經(jīng)驗?zāi)Ｐ头椒ū容^簡單實用，但在求解風速時都增加了反演波高這一環(huán)節(jié)，即在對高頻地波雷達反演風速進行定標時需要先獲得有效波高的反演值或?qū)崪y值。本文根據(jù)風浪經(jīng)驗?zāi)Ｐ?，直接建立了二階譜與一階譜能量之比和風速關(guān)系的雙參數(shù)經(jīng)驗?zāi)Ｐ?，減少了反演有效波高的環(huán)節(jié)。在此基礎(chǔ)上，發(fā)展了風速反演的三參數(shù)經(jīng)驗?zāi)Ｐ?，通過高頻地波雷達反演風速與風速計數(shù)據(jù)的對比，表明雙參數(shù)模型和三參數(shù)經(jīng)驗?zāi)Ｐ投寄芎芎玫剡M行風速的反演，并且三參數(shù)經(jīng)驗?zāi)Ｐ头囱莸慕Y(jié)果略好于雙參數(shù)模型。
2 基本模型
根據(jù)風浪經(jīng)驗?zāi)Ｐ?SMB)(SMB：是由Sverdrup, Munk and Bretschneider 3個人名的首字母組成)，以風浪為主導(dǎo)的海面風速與波高滿足[13]
式中，g是重力加速度，Hs是有效波高，V10是海面上方10 m處風速，fp是峰值頻率。式(1)不宜直接求解，只能采用迭代等方法。SMB關(guān)系式可以簡化成的經(jīng)驗公式為[17]
文獻[16]基于式(2)建立了海面風速與有效波高的經(jīng)驗?zāi)Ｐ?，風速可表示為
式中，a和b是待定系數(shù)。通過經(jīng)驗?zāi)Ｐ头囱莩鲇行Рǜ撸抢么四Ｐ瓦M行風速反演在很大程度上依賴于海面有效波高的反演[16]。
有效波高的反演可根據(jù)文獻[1]的經(jīng)驗?zāi)Ｐ颓蠼狻?/p>
式中，k0為雷達波數(shù)，ω為多普勒頻率，σ(2)和σ(1)分別為二階和一階雷達散射截面，W(ω)表示與兩列海浪波矢之間耦合作用相關(guān)的權(quán)函數(shù)。式(4)可以改寫成
Rw為二階譜能量與一階譜能量的加權(quán)之比。由于加權(quán)系數(shù)只依賴于多普勒頻率，對于一定的區(qū)間可視為常數(shù)，文獻[3]利用了二階譜能量與一階譜能量的無權(quán)重比值R，將式(5)簡化為
文獻[2]則將式(5)推廣為
式中，a和b仍然為待定系數(shù)。
將式(6)代入式(3)，可以得到高頻地波雷達海面回波二階譜能量與一階譜能量之比和風速關(guān)系的雙參數(shù)模型，即
式(8)中參數(shù)a包含了雷達波數(shù)k0。由此可見，對于在風作用下充分成長的海面來說，可以不用通過求解Hs，而直接建立比值R與風速的關(guān)系，從而對風速進行反演求解。式(8)與有效波高計算式(7)形式上基本相同，在統(tǒng)計分析時，需要對a, b參數(shù)進行擬合。為使風速V與R更加有效地擬合，本文中對式(8)進行改進，加入第3個參數(shù)c以加強擬合曲線的上、下偏移的調(diào)節(jié)，形成三參數(shù)模型，即
3 風浪經(jīng)驗?zāi)Ｐ偷膽?yīng)用
本文采用中船重工鵬力(南京)大氣海洋系統(tǒng)有限公司研發(fā)的高頻地波雷達系統(tǒng)分別在福建(2013年10月份)和濱海(2013年2月份)獲得的實驗數(shù)據(jù)對模型進行驗證。雷達系統(tǒng)的發(fā)射頻率分別為7.815 MHz和10.750 MHz，該雷達系統(tǒng)采用了收發(fā)分開共站體制，其工作帶寬為30 kHZ，發(fā)射天線為三元八木天線，接收天線為24元雙排陣型，反演風浪時采用了數(shù)字波束合成技術(shù)，具體的雷達系統(tǒng)參數(shù)見表1。
福建海域的高頻雷達數(shù)據(jù)是2013年10月份獲得的，風速反演區(qū)域距離雷達站點75 km，與正北方向夾角為90o，頻率為7.815 MHz。將反演區(qū)域內(nèi)布放浮標所提供的風速數(shù)據(jù)作為比測數(shù)據(jù)。圖1給出了雷達反演風速結(jié)果與浮標數(shù)據(jù)風速結(jié)果的散點圖，根據(jù)式(9)三參數(shù)模型，擬合得到a, b, c為46.67, 0.35和-15.29，相關(guān)系數(shù)為0.73，均方根誤差為1.73m/s。圖2給出了擬合后的雷達反演結(jié)果與浮標數(shù)據(jù)的比較圖，從圖中看出二者結(jié)果總體符合較好。根據(jù)式(8)雙參數(shù)模型對雷達數(shù)據(jù)和浮標數(shù)據(jù)擬合得到a和b分別為51.62和0.98，均方根誤差為1.82 m/s，相關(guān)系數(shù)為0.71，數(shù)據(jù)散點圖和比較結(jié)果如圖3和圖4所示，從擬合結(jié)果上看，式(9)模型反演結(jié)果要略好于式(8)模型。比較圖2和圖4，在浮標風速在18 m/s左右情況下，三參數(shù)模型要顯著好于雙參數(shù)模型。對于風速較小(＜6 m/s)的情況下，兩種模型得到的結(jié)果與實測結(jié)果差別都較大。由于這兩種模型都是風浪經(jīng)驗?zāi)Ｐ?，都需要滿足風浪充分發(fā)展的海態(tài)條件。因此，風速較小情況下，海浪無法得到足夠的能量，不能充分發(fā)展，因而得到的誤差較大。
表1 高頻地波雷達系統(tǒng)指標
濱海海域采用的是從2013年02月01號8點到2013年02月18號23點的數(shù)據(jù)，雷達頻率為10.75 MHz，選取參考點位置距離雷達站點大約30 km，與正北方向夾角為85o。風速實測數(shù)據(jù)選取濱海岸基自動氣象站風速計同一時段的觀測數(shù)據(jù)。利用三參數(shù)模型得到的結(jié)果如圖5和圖6所示。圖5給出了風速計與高頻地波雷達反演風速的散點圖，擬合得到a=14.24, b=0.24, c=?4.09，相關(guān)系數(shù)為0.5，均方根誤差為1.99 m/s。相關(guān)系數(shù)較低，但均方根誤差可以接受。圖6給出了風速計和高頻地波雷達風速比較圖。從圖上看，雖然風速計測量的風速和高頻雷達反演的風速符合不太理想，但還是能夠反映出風速的變化趨勢。圖中部分對應(yīng)風速計風速較高的數(shù)據(jù)，雷達反演結(jié)果誤差較大。這一方面由于風速計距雷達反演區(qū)域約為30 km，會造成一定的誤差；另一方面，通過觀察雷達多普勒譜，發(fā)現(xiàn)相對應(yīng)的地波雷達二階譜較小，因此可以解釋為雖然風速較大，但是由于海浪沒有充分發(fā)展，無法得到足夠的能量，從而導(dǎo)致地波雷達二階譜較小并且不穩(wěn)定，反演結(jié)果較差。另外，較大的誤差也有可能來自雷達系統(tǒng)或風場的不均勻分布。
根據(jù)兩部雷達的數(shù)據(jù)得到的擬合結(jié)果如表2所示，綜合比較分析可以發(fā)現(xiàn)：(1)雷達頻率較低(7.815 MHz)，在風浪充分成長下，其反演高風速較之低風速要好；而對于雷達頻率較高時(10.750 MHz)，反演低風速較之低頻雷達要好一些，這符合高頻地波雷達的遙感機理；(2)綜合兩個頻率的雷達數(shù)據(jù)風速反演結(jié)果，三參數(shù)模型反演的結(jié)果略好于雙參數(shù)模型的結(jié)果；(3)風速經(jīng)驗?zāi)Ｐ蛯τ诔浞殖砷L的海浪，反演結(jié)果較好，而對于不充分成長的海浪，對風速反演結(jié)果誤差較大，有一定的局限性。
表2 不同模型的風速擬合參數(shù)
圖1 福建雷達三參數(shù)模型反演風速結(jié)果與浮標數(shù)據(jù)散點圖
圖2 福建雷達三參數(shù)模型反演風速結(jié)果與浮標比較圖
圖3 福建雷達雙參數(shù)模型反演風速結(jié)果與浮標數(shù)據(jù)散點圖
圖4 福建雷達雙參數(shù)模型反演風速結(jié)果與浮標比較圖
圖5 濱海雷達三參數(shù)模型反演風速結(jié)果與風速計數(shù)據(jù)散點圖
圖6 濱海雷達三參數(shù)模型反演風速結(jié)果與風速計比較圖
4 結(jié)論
本文在風浪經(jīng)驗?zāi)Ｐ偷幕A(chǔ)上，建立了風速和高頻地波雷達海面回波二階譜能量與一階譜能量之比的關(guān)系，發(fā)展了無需波高信息的高頻地波雷達風速直接反演的經(jīng)驗?zāi)Ｐ?。反演風速經(jīng)驗?zāi)Ｐ瓦m合以風浪為主導(dǎo)，并充分發(fā)展的海面。文中采用兩部不同頻率的雷達在不同海域的數(shù)據(jù)進行了比較分析，結(jié)果表明，文中采用的經(jīng)驗?zāi)Ｐ湍軌蛴行У貙︼L速進行反演，并且三參數(shù)模型反演的結(jié)果略好雙參數(shù)模型。不過從反演結(jié)果上看，對于風浪還沒有充分成長的海面，該模型應(yīng)用受到限制。另外，本文共采用了48天的數(shù)據(jù)，數(shù)據(jù)量相對較少，三參數(shù)模型相對于雙參數(shù)模型的優(yōu)勢有待于獲取更多數(shù)據(jù)進行分析。
致謝：感謝中船重工鵬力(南京)大氣海洋信息系統(tǒng)有限公司周濤研究員提供的高頻地波雷達數(shù)據(jù)及比測數(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] 吳雄斌, 李倫, 李炎. 高頻地波雷達海面有效波高探測實驗研究[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] 李倫, 吳雄斌, 龍超, 等. 基于正則化方法的高頻地波雷達海浪方向譜反演[J]. 地球物理學(xué)報, 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)化理論的高頻地波雷達海浪參數(shù)反演[J], 物理學(xué)報, 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探測海面風浪場原理與實現(xiàn)[J]. 武漢大學(xué)學(xué)報(理學(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] 李倫, 吳雄斌, 徐興安, 等. 高頻地波雷達風速反演經(jīng)驗?zāi)Ｐ蚚J]. 武漢大學(xué)學(xué)報(信息科學(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] 文圣常, 余宙文. 海浪理論與計算原理[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年生，博士，講師，研究方向為X波段導(dǎo)航雷達和高頻地波雷達海態(tài)反演.
張 杰： 男，1963年生，博士，研究員，研究方向為海洋遙感.
王曙曜： 男，1988年生，碩士生，研究方向為高頻地波雷達海洋參數(shù)反演.
紀永剛： 男，1977年生，博士，副研究員，研究方向為高頻地波雷達海態(tài)反演及目標探測.
王祎鳴： 男，1977年生，碩士，助理研究員，研究方向為高頻地波雷達信號處理及雜波抑制.
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改回
國家自然科學(xué)基金重點項目(61032011)和中國博士后科學(xué)基金(2013M531559)資助課題
*通信作者：楚曉亮 xlchu@ouc.edu.cn
猜你喜歡
坐在同一條船上青年文學(xué)家(2023年14期)2023-06-20 12:08:46 反演對稱變換在解決平面幾何問題中的應(yīng)用中等數(shù)學(xué)(2022年5期)2022-08-29 06:07:38 風浪干擾條件下艦船航向保持非線性控制系統(tǒng)艦船科學(xué)技術(shù)(2022年10期)2022-06-17 06:27:54 基于Kmeans-VMD-LSTM的短期風速預(yù)測電機與控制應(yīng)用(2021年12期)2021-02-28 07:55:52 基于最優(yōu)TS評分和頻率匹配的江蘇近海風速訂正海洋通報(2020年5期)2021-01-14 09:26:54 不怕風浪的小船小雪花·小學(xué)生快樂作文(2020年3期)2020-10-13 09:48:18 East Side Story漢語世界(The World of Chinese)(2019年3期)2019-07-01 02:37:48 基于低頻軟約束的疊前AVA稀疏層反演石油地球物理勘探(2017年4期)2017-12-18 07:14:55 基于自適應(yīng)遺傳算法的CSAMT一維反演石油地球物理勘探(2017年2期)2017-11-23 06:02:04 基于GARCH的短時風速預(yù)測方法西南交通大學(xué)學(xué)報(2016年4期)2016-06-15 20:29:37
楚曉亮張 杰王曙曜紀永剛王祎鳴
①(國家海洋局第一海洋研究所 青島 266061)
②(中國海洋大學(xué)信息科學(xué)與工程學(xué)院 青島 266100)
高頻地波雷達風速直接反演的經(jīng)驗?zāi)Ｐ?/p>
楚曉亮*①②張 杰①王曙曜②紀永剛①王祎鳴①
①(國家海洋局第一海洋研究所 青島 266061)
②(中國海洋大學(xué)信息科學(xué)與工程學(xué)院 青島 266100)
利用風浪經(jīng)驗?zāi)Ｐ蛷母哳l地波雷達(HFSWR)的回波譜數(shù)據(jù)中反演風速需要有效波高等先驗信息，因此風速的反演受有效波高反演精度的影響。該文基于風浪經(jīng)驗?zāi)Ｐ?，利用風速和高頻地波雷達海面回波二階譜與一階譜能量之比的關(guān)系，發(fā)展了無需波高信息的風速直接反演的經(jīng)驗?zāi)Ｐ汀L速反演經(jīng)驗?zāi)Ｐ蛻?yīng)用到高頻地波雷達風速的反演中，對兩部不同頻率的雷達在不同海域獲得的數(shù)據(jù)進行了比較分析，結(jié)果表明，該文中采用的經(jīng)驗?zāi)Ｐ湍軌蛴行У貙︼L速進行反演，其中三參數(shù)模型的結(jié)果略好于雙參數(shù)模型。
高頻地波雷達；風速反演；經(jīng)驗?zāi)Ｐ?/p>
1 引言
高頻地波雷達在海洋探測中具有大范圍、全天候及實時性的優(yōu)勢，由于雷達的回波譜中含有豐富的海態(tài)信息，因此從回波譜中可以提取海面流場、風場和海浪等信息。目前，利用高頻地波雷達進行風速反演時，通常是先獲取有效波高后，再借助風浪模型計算出來。有效波高的反演最初是基于文獻[1]以高頻電磁波海面雷達散射理論為基礎(chǔ)而建立的反演經(jīng)驗?zāi)Ｐ汀＠么四Ｐ团c實測數(shù)據(jù)擬合分析獲得相應(yīng)的擬合參數(shù)，從而求得波高，此模型在工程上適用性較強。在此基礎(chǔ)上，文獻[2-4]對模型進行了修正，文獻[5]采用這一經(jīng)驗?zāi)Ｐ屠枚囝l雷達對有效波高進行了反演。另外一些學(xué)者則提出了幾種通過反演出浪高譜而獲得有效波高的方法[6?12]，通過這些方法獲得有效波高及相應(yīng)的峰波頻率后，利用風浪充分發(fā)展條件下的風浪模型即可求得風速[13]。還有一些學(xué)者通過尋找有效波高與雷達回波多普勒譜特征值間的關(guān)系來獲得風速，如文獻[14]在風浪充分發(fā)展的條件下， 通過建立風速與二階譜峰值頻率位置之間的關(guān)系來反演風速，并通過實驗數(shù)據(jù)驗證了該方法的實用性，但這種方法對獲得二階譜峰值點位置的精確度要求較高。文獻[15]通過人工神經(jīng)網(wǎng)絡(luò)方法對風速進行了反演，取得了不錯的結(jié)果。文獻[16]建立了海面風速與有效波高的經(jīng)驗?zāi)Ｐ?，通過對實測數(shù)據(jù)的反演對比驗證了該模型的穩(wěn)定性。在實際應(yīng)用中，經(jīng)驗?zāi)Ｐ头椒ū容^簡單實用，但在求解風速時都增加了反演波高這一環(huán)節(jié)，即在對高頻地波雷達反演風速進行定標時需要先獲得有效波高的反演值或?qū)崪y值。本文根據(jù)風浪經(jīng)驗?zāi)Ｐ?，直接建立了二階譜與一階譜能量之比和風速關(guān)系的雙參數(shù)經(jīng)驗?zāi)Ｐ?，減少了反演有效波高的環(huán)節(jié)。在此基礎(chǔ)上，發(fā)展了風速反演的三參數(shù)經(jīng)驗?zāi)Ｐ?，通過高頻地波雷達反演風速與風速計數(shù)據(jù)的對比，表明雙參數(shù)模型和三參數(shù)經(jīng)驗?zāi)Ｐ投寄芎芎玫剡M行風速的反演，并且三參數(shù)經(jīng)驗?zāi)Ｐ头囱莸慕Y(jié)果略好于雙參數(shù)模型。
2 基本模型
根據(jù)風浪經(jīng)驗?zāi)Ｐ?SMB)(SMB：是由Sverdrup, Munk and Bretschneider 3個人名的首字母組成)，以風浪為主導(dǎo)的海面風速與波高滿足[13]
式中，g是重力加速度，Hs是有效波高，V10是海面上方10 m處風速，fp是峰值頻率。式(1)不宜直接求解，只能采用迭代等方法。SMB關(guān)系式可以簡化成的經(jīng)驗公式為[17]
文獻[16]基于式(2)建立了海面風速與有效波高的經(jīng)驗?zāi)Ｐ?，風速可表示為
式中，a和b是待定系數(shù)。通過經(jīng)驗?zāi)Ｐ头囱莩鲇行Рǜ撸抢么四Ｐ瓦M行風速反演在很大程度上依賴于海面有效波高的反演[16]。
有效波高的反演可根據(jù)文獻[1]的經(jīng)驗?zāi)Ｐ颓蠼狻?/p>
式中，k0為雷達波數(shù)，ω為多普勒頻率，σ(2)和σ(1)分別為二階和一階雷達散射截面，W(ω)表示與兩列海浪波矢之間耦合作用相關(guān)的權(quán)函數(shù)。式(4)可以改寫成
Rw為二階譜能量與一階譜能量的加權(quán)之比。由于加權(quán)系數(shù)只依賴于多普勒頻率，對于一定的區(qū)間可視為常數(shù)，文獻[3]利用了二階譜能量與一階譜能量的無權(quán)重比值R，將式(5)簡化為
文獻[2]則將式(5)推廣為
式中，a和b仍然為待定系數(shù)。
將式(6)代入式(3)，可以得到高頻地波雷達海面回波二階譜能量與一階譜能量之比和風速關(guān)系的雙參數(shù)模型，即
式(8)中參數(shù)a包含了雷達波數(shù)k0。由此可見，對于在風作用下充分成長的海面來說，可以不用通過求解Hs，而直接建立比值R與風速的關(guān)系，從而對風速進行反演求解。式(8)與有效波高計算式(7)形式上基本相同，在統(tǒng)計分析時，需要對a, b參數(shù)進行擬合。為使風速V與R更加有效地擬合，本文中對式(8)進行改進，加入第3個參數(shù)c以加強擬合曲線的上、下偏移的調(diào)節(jié)，形成三參數(shù)模型，即
3 風浪經(jīng)驗?zāi)Ｐ偷膽?yīng)用
本文采用中船重工鵬力(南京)大氣海洋系統(tǒng)有限公司研發(fā)的高頻地波雷達系統(tǒng)分別在福建(2013年10月份)和濱海(2013年2月份)獲得的實驗數(shù)據(jù)對模型進行驗證。雷達系統(tǒng)的發(fā)射頻率分別為7.815 MHz和10.750 MHz，該雷達系統(tǒng)采用了收發(fā)分開共站體制，其工作帶寬為30 kHZ，發(fā)射天線為三元八木天線，接收天線為24元雙排陣型，反演風浪時采用了數(shù)字波束合成技術(shù)，具體的雷達系統(tǒng)參數(shù)見表1。
福建海域的高頻雷達數(shù)據(jù)是2013年10月份獲得的，風速反演區(qū)域距離雷達站點75 km，與正北方向夾角為90o，頻率為7.815 MHz。將反演區(qū)域內(nèi)布放浮標所提供的風速數(shù)據(jù)作為比測數(shù)據(jù)。圖1給出了雷達反演風速結(jié)果與浮標數(shù)據(jù)風速結(jié)果的散點圖，根據(jù)式(9)三參數(shù)模型，擬合得到a, b, c為46.67, 0.35和-15.29，相關(guān)系數(shù)為0.73，均方根誤差為1.73m/s。圖2給出了擬合后的雷達反演結(jié)果與浮標數(shù)據(jù)的比較圖，從圖中看出二者結(jié)果總體符合較好。根據(jù)式(8)雙參數(shù)模型對雷達數(shù)據(jù)和浮標數(shù)據(jù)擬合得到a和b分別為51.62和0.98，均方根誤差為1.82 m/s，相關(guān)系數(shù)為0.71，數(shù)據(jù)散點圖和比較結(jié)果如圖3和圖4所示，從擬合結(jié)果上看，式(9)模型反演結(jié)果要略好于式(8)模型。比較圖2和圖4，在浮標風速在18 m/s左右情況下，三參數(shù)模型要顯著好于雙參數(shù)模型。對于風速較小(＜6 m/s)的情況下，兩種模型得到的結(jié)果與實測結(jié)果差別都較大。由于這兩種模型都是風浪經(jīng)驗?zāi)Ｐ?，都需要滿足風浪充分發(fā)展的海態(tài)條件。因此，風速較小情況下，海浪無法得到足夠的能量，不能充分發(fā)展，因而得到的誤差較大。
表1 高頻地波雷達系統(tǒng)指標
濱海海域采用的是從2013年02月01號8點到2013年02月18號23點的數(shù)據(jù)，雷達頻率為10.75 MHz，選取參考點位置距離雷達站點大約30 km，與正北方向夾角為85o。風速實測數(shù)據(jù)選取濱海岸基自動氣象站風速計同一時段的觀測數(shù)據(jù)。利用三參數(shù)模型得到的結(jié)果如圖5和圖6所示。圖5給出了風速計與高頻地波雷達反演風速的散點圖，擬合得到a=14.24, b=0.24, c=?4.09，相關(guān)系數(shù)為0.5，均方根誤差為1.99 m/s。相關(guān)系數(shù)較低，但均方根誤差可以接受。圖6給出了風速計和高頻地波雷達風速比較圖。從圖上看，雖然風速計測量的風速和高頻雷達反演的風速符合不太理想，但還是能夠反映出風速的變化趨勢。圖中部分對應(yīng)風速計風速較高的數(shù)據(jù)，雷達反演結(jié)果誤差較大。這一方面由于風速計距雷達反演區(qū)域約為30 km，會造成一定的誤差；另一方面，通過觀察雷達多普勒譜，發(fā)現(xiàn)相對應(yīng)的地波雷達二階譜較小，因此可以解釋為雖然風速較大，但是由于海浪沒有充分發(fā)展，無法得到足夠的能量，從而導(dǎo)致地波雷達二階譜較小并且不穩(wěn)定，反演結(jié)果較差。另外，較大的誤差也有可能來自雷達系統(tǒng)或風場的不均勻分布。
根據(jù)兩部雷達的數(shù)據(jù)得到的擬合結(jié)果如表2所示，綜合比較分析可以發(fā)現(xiàn)：(1)雷達頻率較低(7.815 MHz)，在風浪充分成長下，其反演高風速較之低風速要好；而對于雷達頻率較高時(10.750 MHz)，反演低風速較之低頻雷達要好一些，這符合高頻地波雷達的遙感機理；(2)綜合兩個頻率的雷達數(shù)據(jù)風速反演結(jié)果，三參數(shù)模型反演的結(jié)果略好于雙參數(shù)模型的結(jié)果；(3)風速經(jīng)驗?zāi)Ｐ蛯τ诔浞殖砷L的海浪，反演結(jié)果較好，而對于不充分成長的海浪，對風速反演結(jié)果誤差較大，有一定的局限性。
表2 不同模型的風速擬合參數(shù)
圖1 福建雷達三參數(shù)模型反演風速結(jié)果與浮標數(shù)據(jù)散點圖
圖2 福建雷達三參數(shù)模型反演風速結(jié)果與浮標比較圖
圖3 福建雷達雙參數(shù)模型反演風速結(jié)果與浮標數(shù)據(jù)散點圖
圖4 福建雷達雙參數(shù)模型反演風速結(jié)果與浮標比較圖
圖5 濱海雷達三參數(shù)模型反演風速結(jié)果與風速計數(shù)據(jù)散點圖
圖6 濱海雷達三參數(shù)模型反演風速結(jié)果與風速計比較圖
4 結(jié)論
本文在風浪經(jīng)驗?zāi)Ｐ偷幕A(chǔ)上，建立了風速和高頻地波雷達海面回波二階譜能量與一階譜能量之比的關(guān)系，發(fā)展了無需波高信息的高頻地波雷達風速直接反演的經(jīng)驗?zāi)Ｐ?。反演風速經(jīng)驗?zāi)Ｐ瓦m合以風浪為主導(dǎo)，并充分發(fā)展的海面。文中采用兩部不同頻率的雷達在不同海域的數(shù)據(jù)進行了比較分析，結(jié)果表明，文中采用的經(jīng)驗?zāi)Ｐ湍軌蛴行У貙︼L速進行反演，并且三參數(shù)模型反演的結(jié)果略好雙參數(shù)模型。不過從反演結(jié)果上看，對于風浪還沒有充分成長的海面，該模型應(yīng)用受到限制。另外，本文共采用了48天的數(shù)據(jù)，數(shù)據(jù)量相對較少，三參數(shù)模型相對于雙參數(shù)模型的優(yōu)勢有待于獲取更多數(shù)據(jù)進行分析。
致謝：感謝中船重工鵬力(南京)大氣海洋信息系統(tǒng)有限公司周濤研究員提供的高頻地波雷達數(shù)據(jù)及比測數(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] 吳雄斌, 李倫, 李炎. 高頻地波雷達海面有效波高探測實驗研究[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] 李倫, 吳雄斌, 龍超, 等. 基于正則化方法的高頻地波雷達海浪方向譜反演[J]. 地球物理學(xué)報, 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)化理論的高頻地波雷達海浪參數(shù)反演[J], 物理學(xué)報, 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探測海面風浪場原理與實現(xiàn)[J]. 武漢大學(xué)學(xué)報(理學(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] 李倫, 吳雄斌, 徐興安, 等. 高頻地波雷達風速反演經(jīng)驗?zāi)Ｐ蚚J]. 武漢大學(xué)學(xué)報(信息科學(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] 文圣常, 余宙文. 海浪理論與計算原理[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年生，博士，講師，研究方向為X波段導(dǎo)航雷達和高頻地波雷達海態(tài)反演.
張 杰： 男，1963年生，博士，研究員，研究方向為海洋遙感.
王曙曜： 男，1988年生，碩士生，研究方向為高頻地波雷達海洋參數(shù)反演.
紀永剛： 男，1977年生，博士，副研究員，研究方向為高頻地波雷達海態(tài)反演及目標探測.
王祎鳴： 男，1977年生，碩士，助理研究員，研究方向為高頻地波雷達信號處理及雜波抑制.
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改回
國家自然科學(xué)基金重點項目(61032011)和中國博士后科學(xué)基金(2013M531559)資助課題
*通信作者：楚曉亮 xlchu@ouc.edu.cn
猜你喜歡
青年文學(xué)家(2023年14期)2023-06-20 12:08:46
中等數(shù)學(xué)(2022年5期)2022-08-29 06:07:38
艦船科學(xué)技術(shù)(2022年10期)2022-06-17 06:27:54
電機與控制應(yīng)用(2021年12期)2021-02-28 07:55:52
海洋通報(2020年5期)2021-01-14 09:26:54
小雪花·小學(xué)生快樂作文(2020年3期)2020-10-13 09:48:18
漢語世界(The World of Chinese)(2019年3期)2019-07-01 02:37:48
石油地球物理勘探(2017年4期)2017-12-18 07:14:55
石油地球物理勘探(2017年2期)2017-11-23 06:02:04
西南交通大學(xué)學(xué)報(2016年4期)2016-06-15 20:29:37
