CN107064929B - 一种利用s波段多普勒雷达探测海面浪高的方法 - Google Patents
一种利用s波段多普勒雷达探测海面浪高的方法 Download PDFInfo
- Publication number
- CN107064929B CN107064929B CN201710408429.8A CN201710408429A CN107064929B CN 107064929 B CN107064929 B CN 107064929B CN 201710408429 A CN201710408429 A CN 201710408429A CN 107064929 B CN107064929 B CN 107064929B
- Authority
- CN
- China
- Prior art keywords
- spectrum
- wave
- radar
- formula
- velocity
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种利用S波段多普勒雷达探测海面浪高的方法,包括步骤1,利用S波段多普勒雷达获取海面回波多普勒谱,对得到的N帧多普勒谱提取中心频率,获得中心频率序列,根据多普勒效应,得到雷达测量的径向速度序列;步骤2,对雷达测得的径向速度序列去除平均值,接着作离散傅里叶变换和取模平方得到速度谱;步骤3,根据线性波动理论,得到速度谱与海浪谱的转换关系,获得海浪谱;步骤4,根据有效波高与海浪谱之间的数学关系,利用谱矩法获得有效波高;步骤5,用零阶矩去奇异值法,使异常的有效波高恢复正常。本发明为利用微波多普勒雷达获取海面浪高提供了一种新的探测手段,具有较高的探测精度和浪场空间分辨率,结果简单清晰。
Description
技术领域
本发明属于雷达探测领域,尤其涉及一种利用S波段多普勒雷达探测海面浪高的方法。
背景技术
近岸处海浪对于海运和港口安全有着非常重要的作用,传统的测量手段如浮标、潜标、海洋调查船等只能提供局部有效的海洋数据,难以满足实际需求。随着计算机等新技术的快速发展,无线电磁波逐步应用于海洋监测领域。S波段多普勒雷达能够探测近岸处海面浪高信息,且由于其体积小、受环境干扰少、测量精度高等特点,目前已逐步发展成为一种用于监测海洋环境的新手段。对于S波段多普勒雷达,当入射角较大时,海面与电磁波信号发生Bragg散射,谐振的海浪波长与电磁波长处于同一量级,被称为Bragg毛细波。根据复合表面散射理论,毛细波跨骑在重力波上,并受到重力波的调制作用,这些调制包括轨道调制、倾斜调制、阴影调制和流体动力学调制等,因此Bragg毛细波的变化中携带了重力波的特征信息。从雷达回波中提取出重力波的速度信息,利用线性波动理论,速度与海浪谱存在着转换关系,进而获取海浪谱。因此,S波段多普勒雷达发射的电磁波与海面发生相互作用,其幅度和相位受到海浪、海流等参数调制,对调制后的回波信号经过一系列复杂的信号处理和数据分析,可从中提取出海面浪高信息。
发明内容
本发明针对上述问题,提出了一种利用S波段多普勒雷达探测海面浪高的方法。该方法的核心在于通过雷达回波多普勒快速、高精度地实时获取海面浪高信息。该方法不仅适用于S波段多普勒雷达,还可应用于其他波段的微波多普勒雷达。
本发明的技术方案是一种利用S波段多普勒雷达探测海面浪高的方法,包括以下步骤:
步骤1,利用S波段多普勒雷达获取海面回波多普勒谱,对得到的N帧多普勒谱提取中心频率,获得中心频率序列,根据多普勒效应公式,得到雷达测量的径向速度序列v(t)。
在该步骤中,中心频率的获取采用基于矩运算的能量均衡法:
式中:P(f)为多普勒功率谱密度函数。
步骤2,对雷达测得的径向速度序列v(t)去除平均值,接着作离散傅里叶变换和取模平方得到速度谱v2(f)。
步骤3,根据线性波动理论,得到速度谱与海浪谱的转换关系,获得海浪谱S(f);
在该步骤中,速度谱与海浪方向谱的转换关系如下:
式中:v2(f)表示速度谱,k和ω分别是水质点毛细波对应的波数和频率,Δf表示频率间隔,Δθ表示角度间隔,d是水深,α为雷达掠射角;
海浪谱通过下式获得:
式中:i表示第i根天线的海浪方向谱,雷达共有6根天线;
步骤4,根据有效波高与海浪谱之间的数学关系,利用谱矩法获得有效波高Hs。
在该步骤中,有效波高的计算公式如下:
式中,Hs为海浪的有效波高,m0为海浪谱的零阶矩,可通过下式得到:
m0=∫S(f)df (4)
步骤5,有效波高受地形、破碎波以及山体遮挡等因素影响,偶尔出现异常。在海浪谱零阶矩序列中应用零阶矩去奇异值法,使得异常的有效波高恢复正常。
在所述步骤2和步骤3中,由径向速度序列推导海浪谱的过程如下:
雷达测得的径向速度为:
式中:v(t)为雷达测得的径向速度,α为雷达掠射角,θ为雷达天线指向;对径向速度序列做平均,去除散射体相速度和海流径向流速;
接着,对去除平均值后的径向速度序列进行离散傅里叶变换和取模平方后得到速度谱:
式中:v2(f)表示速度谱;
在式(7)中,毛细波幅度A(f,θ)与海浪方向谱S(f,θ)存在以下的关系式:
式中:Δf表示频率间隔,Δθ表示角度间隔;
将式(8)代入式(7)就得到海浪方向谱与速度谱之间的转换关系,两者的转换公式如下所示:
海浪谱通过下式获得:
式中:i表示第i根天线的海浪方向谱,雷达共有6根天线。
在所述步骤5中,零阶矩去奇异值法的计算步骤如下:
a.利用转换公式得到海浪谱对应的零阶矩序列;
b.计算该序列平均值,记为A,设置k A为阈值,其中k为常数,取4~6;
c.检测高于阈值的位置并剔除,得到新的零阶矩序列;
d.计算新序列平均值,记为A1,设置k1A1为新阈值,其中k1为常数,取8~10;
e.检测原始零阶矩序列中高于新阈值的部分,标记此部分为异常位置;
f.重置异常位置的零阶矩为A1,其他位置不变,即得到剔除异常值后的零阶矩序列。
在所述步骤1和步骤3中,N的取值为2n,n是不小于7的正整数。
与现有技术相比,本发明具有如下优点:
1、本发明利用S波段多普勒雷达,为海面有效波高的获取提供了一种新的方法。
2、本发明利用S波段多普勒雷达作为探测工具,该款雷达体积小,受外部环境干扰少,且可全天候工作,不受雨雾影响。
3、本发明计算方法快速简单,满足雷达实时工作要求。
4、本发明能够连续、实时地观测有效波高,每间隔数分钟就能获取海洋表面有效波高分布,高效快速。
5、本发明针对异常的有效波高提出了一种简单快速的处理方法,使得有效波高恢复正常。
本发明不仅可适用于S波段多普勒雷达,还可应用于其他波段的微波多普勒雷达系统中。
附图说明
图1,本发明的算法流程图;
图2,提取的有效波高随距离元变化的示意图;其中,图2(a)为原始的有效波高结果;图2(b)为应用零阶矩去奇异值法后的有效波高结果;
图3,应用本发明方法后提取出的有效波高与浮标对比结果。
具体实施方式
下面以具体实例结合附图,对本发明作进一步说明。
海浪是海面起伏形状的传播,是水质点离开平衡位置,作周期性振动,并向一定方向传播而形成的一种波动。海面上水质点的速度可通过下式计算得到:
u(t)=U0cos(ωt)+vc+vs (5)
式中:u(t)和w(t)分别为水质点速度的水平分量和垂直分量,vc是海流的径向流速,vs是有效散射体的相速度,k和ω分别是水质点毛细波对应的波数和频率,d是水深,C是长重力波的相速度,一般远大于vs和vc。U0是轨道速度,由下式获得:
U0=ωA(f,θ)coth(kd) (7)
其中,A(f,θ)是毛细波对应的幅度。
S波段多普勒雷达直接对来自海面的回波信号进行两次傅里叶变换,得到回波多普勒谱,由于多普勒效应,海面水质点速度将会对雷达回波产生作用。因此雷达测得的径向速度为:
式中:v(t)为雷达测得的径向速度,α为雷达掠射角。上式第一项是均值为零的正弦信号,周期一般在十秒左右,因此对径向速度序列做几分钟的平均,可去除散射体相速度和海流径向流速。
具体实施步骤参见图1:
步骤1,利用S波段多普勒雷达获取海面回波多普勒谱,对得到的N帧多普勒谱提取中心频率,获得中心频率序列,根据多普勒效应公式,得到雷达测量的径向速度序列v(t)。
在该步骤中,中心频率的获取采用基于矩运算的能量均衡法:
式中:P(f)为雷达回波多普勒的功率谱密度函数。
步骤2,对雷达测得的径向速度序列v(t)去除平均值,接着作离散傅里叶变换和取模平方得到速度谱v2(f)。
雷达测量的径向速度序列v(t)如式(8)所示,对其去除平均值后,第二项和第三项的散射体相速度和海流径向流速被去除。经过离散傅里叶变换和取模平方后:
式中:v2(f)表示速度谱。
步骤3,根据线性波动理论,得到速度谱与海浪谱的转换关系,获得海浪谱S(f);
在该步骤中,速度谱与海浪方向谱的转换关系如下:
式中:v2(f)表示速度谱,k和ω分别是水质点毛细波对应的波数和频率,Δf表示频率间隔,Δθ表示角度间隔,d是水深,α为雷达掠射角;
海浪谱可通过下式获得:
式中:i表示第i根天线的海浪方向谱,雷达共有6根天线。
步骤4,根据有效波高与海浪谱之间的数学关系,利用谱矩法获得有效波高Hs。
在该步骤中,有效波高的计算公式如下:
式中,Hs为海浪的有效波高,m0为海浪谱的零阶矩,可通过下式得到:
m0=∫S(f)df (14)
步骤5,有效波高受地形、破碎波以及山体遮挡等因素影响,偶尔出现异常。在海浪谱零阶矩序列中应用零阶矩去奇异值法,使得异常的有效波高恢复正常。
由公示(13)可知,有效波高由海浪谱的零阶矩决定,零阶矩的异常导致了有效波高出现异常值,因此在海浪谱中对异常零阶矩进行检测并处理。基于海浪的短时平稳特征,提出了零阶矩去奇异值法。它是一种实际的工程处理方法,具体计算步骤如下:
a.利用转换公式得到海浪谱对应的零阶矩序列;
b.计算该序列平均值,记为A,设置k A为阈值,其中k为常数,由于正常序列波动很小,因此k取4~6就能检测出异常值;
c.检测高于阈值的位置并剔除,得到新的零阶矩序列;
d.计算新序列平均值,记为A1,设置k1A1为新阈值,其中k1为常数,考虑到实际海况中有效波高通常低于3倍平均值,k1一般取8~10;
e.检测原始零阶矩序列中高于新阈值的部分,标记此部分为异常位置;
f.重置异常位置的零阶矩为A1,其他位置不变,这样就得到了剔除异常值后的零阶矩序列。
图2为经过步骤5处理前后的有效波高结果,经过零阶矩去奇异值法处理后,异常的有效波高恢复正常。
图3为利用本发明方法提取的有效波高和浮标结果对比,两者之间均方根误差为0.13m,相关系数为0.92,在观测期间,海面波高变化幅度较大,雷达测量的结果始终与浮标高度一致,说明本发明方法获得的有效波高具有很高的准确性。
Claims (4)
1.一种利用S波段多普勒雷达探测海面浪高的方法,其特征在于:包括以下步骤:
步骤1,利用S波段多普勒雷达获取海面回波多普勒谱,对得到的N帧多普勒谱提取中心频率,获得中心频率序列,根据多普勒效应公式,得到雷达测量的径向速度序列v(t);
在该步骤中,中心频率f的获取采用基于矩运算的能量均衡法:
式中:P(f)为多普勒功率谱密度函数;
步骤2,对雷达测得的径向速度序列v(t)去除平均值,接着作离散傅里叶变换和取模平方得到速度谱v2(f);
步骤3,根据线性波动理论,得到速度谱与海浪谱的转换关系,获得海浪谱S(f);
在该步骤中,速度谱与海浪方向谱的转换关系如下:
式中:v2(f)表示速度谱,k和ω分别是水质点毛细波对应的波数和频率,Δf表示频率间隔,Δθ表示角度间隔,d是水深,α为雷达掠射角;
海浪谱通过下式获得:
式中:i表示第i根天线的海浪方向谱,雷达共有6根天线;
步骤4,根据有效波高与海浪谱之间的数学关系,利用谱矩法获得有效波高Hs;
在该步骤中,有效波高的计算公式如下:
式中,Hs为海浪的有效波高,m0为海浪谱的零阶矩,通过下式得到:
m0=∫S(f)df (5)
步骤5,在海浪谱零阶矩序列中应用零阶矩去奇异值法,使异常的有效波高恢复正常。
2.如权利要求1所述的一种利用S波段多普勒雷达探测海面浪高的方法,其特征在于:在所述步骤2和步骤3中,由径向速度序列推导海浪谱的过程如下:
雷达测得的径向速度为:
式中:v(t)为雷达测得的径向速度,α为雷达掠射角,θ为雷达天线指向;对径向速度序列做平均,去除散射体相速度和海流径向流速;
接着,对去除平均值后的径向速度序列进行离散傅里叶变换和取模平方后得到速度谱:
式中:v2(f)表示速度谱;
在式(7)中,毛细波幅度A(f,θ)与海浪方向谱S(f,θ)存在以下的关系式:
式中:Δf表示频率间隔,Δθ表示角度间隔;
将式(8)代入式(7)就得到海浪方向谱与速度谱之间的转换关系,两者的转换公式如下所示:
海浪谱通过下式获得:
式中:i表示第i根天线的海浪方向谱,雷达共有6根天线。
3.如权利要求1所述的一种利用S波段多普勒雷达探测海面浪高的方法,其特征在于:在所述步骤5中,零阶矩去奇异值法的计算步骤如下:
a.利用转换公式得到海浪谱对应的零阶矩序列;
b.计算该序列平均值,记为A,设置k A为阈值,其中k为常数,取4~6;
c.检测高于阈值的位置并剔除,得到新的零阶矩序列;
d.计算新序列平均值,记为A1,设置k1A1为新阈值,其中k1为常数,取8~10;
e.检测原始零阶矩序列中高于新阈值的部分,标记此部分为异常位置;
f.重置异常位置的零阶矩为A1,其他位置不变,即得到剔除异常值后的零阶矩序列。
4.如权利要求1或2所述的一种利用S波段多普勒雷达探测海面浪高的方法,其特征在于:在所述步骤1和步骤3中,N的取值为2n,n是不小于7的正整数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710408429.8A CN107064929B (zh) | 2017-06-02 | 2017-06-02 | 一种利用s波段多普勒雷达探测海面浪高的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710408429.8A CN107064929B (zh) | 2017-06-02 | 2017-06-02 | 一种利用s波段多普勒雷达探测海面浪高的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107064929A CN107064929A (zh) | 2017-08-18 |
CN107064929B true CN107064929B (zh) | 2020-02-14 |
Family
ID=59616130
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710408429.8A Active CN107064929B (zh) | 2017-06-02 | 2017-06-02 | 一种利用s波段多普勒雷达探测海面浪高的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107064929B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108802724A (zh) * | 2018-06-28 | 2018-11-13 | 北京豪恩科技有限责任公司 | 一种测试海浪高度的方法 |
CN112130122B (zh) * | 2020-09-01 | 2022-09-13 | 武汉大学 | 一种空基高频雷达海面散射系数的估计方法 |
CN113030894B (zh) * | 2021-03-02 | 2022-06-28 | 南京信息工程大学 | 一种利用快速扫描的相参雷达图像提取海浪参数的方法 |
CN113156393B (zh) * | 2021-03-29 | 2022-05-27 | 山东科技大学 | 一种机载激光测深破碎风浪海面模型构建方法 |
CN113900096B (zh) * | 2021-09-09 | 2024-04-26 | 武汉大学 | 一种双基地高频地波雷达海洋回波反演无向海浪谱的方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102879766A (zh) * | 2011-07-11 | 2013-01-16 | 哈尔滨工业大学 | 一种高频地波雷达弱目标检测跟踪方法和装置 |
-
2017
- 2017-06-02 CN CN201710408429.8A patent/CN107064929B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102879766A (zh) * | 2011-07-11 | 2013-01-16 | 哈尔滨工业大学 | 一种高频地波雷达弱目标检测跟踪方法和装置 |
Non-Patent Citations (3)
Title |
---|
"Doppler Spectra of Microwave Scattering Fields From Nonlinear Oceanic Surface at Moderate and Low-Grazing Angles";Yunhua Wang et al.;《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》;20120430;第50卷(第4期);第1102-1116页 * |
"基于非线性海面微波段多普勒回波的建模算法";陈泽宗 等;《电子技术应用》;20160531(第5期);第97-100页 * |
"微波多普勒雷达海浪参数提取算法";范林刚 等;《华中科技大学学报(自然科学版)》;20121031;第40卷(第10期);第21-24页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107064929A (zh) | 2017-08-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107064929B (zh) | 一种利用s波段多普勒雷达探测海面浪高的方法 | |
CN107678003B (zh) | 一种地波雷达海杂波背景下的目标检测方法 | |
Menelle et al. | Full digital high frequency surface wave radar: French trials in the Biscay bay | |
Lund et al. | Wind retrieval from shipborne nautical X-band radar data | |
CN108226933B (zh) | 一种基于条纹干涉结构的深海宽带目标深度估计方法 | |
CN111006743B (zh) | 一种基于平面雷达水位计的水位测量排除干扰的方法 | |
CN111665489B (zh) | 一种基于目标特性的线谱提取方法 | |
Lyzenga et al. | Real-time estimation of ocean wave fields from marine radar data | |
Thornhill et al. | Real time local sea state measurement using wave radar and ship motions | |
Xie et al. | Measuring ocean surface wind field using shipborne high-frequency surface wave radar | |
Braun et al. | Sea-surface current features observed by Doppler radar | |
Bouferrouk et al. | Field measurements of surface waves using a 5-beam ADCP | |
US20120078516A1 (en) | Water quality measuring method, water quality measuring device, and water quality measuring program for marine surface layers | |
Liu et al. | A new Doppler model incorporated with free and broken-short waves for coherent S-band wave radar at near-grazing angles | |
Zhao et al. | Inversion of wave parameters from time-Doppler spectrum using shore-based coherent S-band radar | |
Green et al. | An inversion method for extraction of wind speed from high-frequency ground-wave radar oceanic backscatter | |
Park et al. | Estimation of significant wave heights from x-band radar using artificial neural network | |
Sabinin et al. | Results of using acoustic doppler current profilers for studying the spatial structure of the marine environment | |
Deng et al. | Wave height and wave period measurements using small-aperture HF radar | |
CN114296046B (zh) | 基于人工神经网络的hfswr多海况有效波高提取方法与装置 | |
Liu et al. | Acoustic spectrum and signature analysis on underwater radiated noise of a passenger ship target based on the measured data | |
CN111948607B (zh) | 一种深海目标被动定位及粒子滤波跟踪方法 | |
Guo et al. | Modeling and simulation of water-surface vibration due to acoustic signals for detection with terahertz radar | |
Govea et al. | Motion compensation of opean-ocean acoustic Doppler velocimeter measurements for offshore turbulence assessment | |
Tian et al. | Wave height field measurement using a compact dual-frequency HF radar |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |