CN108008392A - 一种基于船载高频地波雷达的海洋表面风场测量方法 - Google Patents
一种基于船载高频地波雷达的海洋表面风场测量方法 Download PDFInfo
- Publication number
- CN108008392A CN108008392A CN201711173938.3A CN201711173938A CN108008392A CN 108008392 A CN108008392 A CN 108008392A CN 201711173938 A CN201711173938 A CN 201711173938A CN 108008392 A CN108008392 A CN 108008392A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mfrac
- msup
- ocean
- 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.)
- Granted
Links
- 238000000691 measurement method Methods 0.000 title claims abstract description 13
- 238000001514 detection method Methods 0.000 claims abstract description 24
- 108010003272 Hyaluronate lyase Proteins 0.000 claims abstract description 14
- 238000001228 spectrum Methods 0.000 claims abstract description 14
- 238000000034 method Methods 0.000 claims description 23
- 230000008569 process Effects 0.000 claims description 14
- JXASPPWQHFOWPL-UHFFFAOYSA-N Tamarixin Natural products C1=C(O)C(OC)=CC=C1C1=C(OC2C(C(O)C(O)C(CO)O2)O)C(=O)C2=C(O)C=C(O)C=C2O1 JXASPPWQHFOWPL-UHFFFAOYSA-N 0.000 claims description 2
- 230000001133 acceleration Effects 0.000 claims description 2
- 230000005540 biological transmission Effects 0.000 claims description 2
- 239000004576 sand Substances 0.000 claims 1
- 230000003595 spectral effect Effects 0.000 claims 1
- 238000005259 measurement Methods 0.000 abstract description 14
- 238000000605 extraction Methods 0.000 abstract 1
- 238000012423 maintenance Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 7
- 230000009286 beneficial effect Effects 0.000 description 4
- 238000009826 distribution Methods 0.000 description 3
- 238000002592 echocardiography Methods 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000004083 survival effect Effects 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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
- G01S13/95—Radar or analogous systems specially adapted for specific applications for meteorological use
- G01S13/956—Radar or analogous systems specially adapted for specific applications for meteorological use mounted on ship or other platform
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Ocean & Marine Engineering (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
一种基于船载高频地波雷达的海洋表面风场测量方法，本发明涉及海洋表面风场测量方法。本发明解决了传统海洋表面风场测量工具部署维护困难、无法大面积测量和观测效率低的问题；以及岸基高频地波雷达探测海域与测量性能受限的问题。一：选取某一距离门对应的一阶海面回波多普勒展宽谱；二：确定被探测方位的海洋单元回波的入射方位角对应的多普勒频率；三：建立风向与扩展因子的对应关系；四：提取无模糊风向和对应的扩展因子；五：建立风速与扩展因子的对应关系；六：重复执行二到五，依次得到该距离门内不同探测方位海洋单元对应的风场；七：重复执行N遍一到六，得到雷达探测范围内整个海域的海面风场。本发明用于海洋表面风场测量领域。
Description
技术领域
本发明涉及海洋表面风场测量方法。
背景技术
在海洋表面动力学参数中，海面风是研究大气和海洋之间能量和物质交换的重要因素，是影响海浪和水流等要素的活跃因子，是海洋学研究的重点对象之一，对海洋工程、海上资源开发、气象预报以及海洋学的研究等方面具有重要影响。然而传统的海洋表面风场测量工具部署维护困难、无法大面积测量及观测效率低。
高频地波雷达的工作频率为3-30MHz，高频地波雷达利用垂直极化的电磁波沿海面传播绕射性好、衰减小的特性，可以突破地球曲率的限制，探测到视线以下的海域，已经成为海洋表面动力学参数遥感的重要工具。与传统的海洋表面动力学监测工具相比，高频地波雷达能够提供一个全天时、全天候、大面积以及超视距的监测环境。
按照雷达系统的部署方式，高频地波雷达可以分为岸基型和船载型。岸基高频地波雷达利用海面对电磁波的一阶散射机制，可以从海面回波提取风向、海流等信息。然而，岸基高频地波雷达需要大型的天线接受阵列，利用数字波束形成算法或超分辨算法进行波达方向估计，这使得岸基高频地波雷达不但造价高昂、系统复杂，而且覆盖海域固定，限制了其应用范围。船载高频地波雷达除了保留岸基高频地波雷达的优点外，还具有机动灵活性，不但扩展了雷达的探测范围，而且在复杂的海洋环境中增加了生存能力。由于船载平台的前向运动，高频地波雷达一阶多普勒谱被展宽，使得来自不同方位的海面回波对应不同的多普勒频率，因此船载高频地波雷达可以利用单接收阵元进行高方位分辨力的风场测量。
发明内容
本发明的目的为了解决传统海洋表面风场测量工具部署维护困难、无法大面积测量和观测效率低的问题；以及岸基高频地波雷达探测海域、测量性能受限的问题，而提出一种基于船载高频地波雷达的海洋表面风场测量方法。
一种基于船载高频地波雷达的海洋表面风场测量方法具体过程为：
步骤一：选取某一距离门对应的一阶海面回波多普勒展宽谱；
步骤二：根据雷达系统的距离分辨力等间距划分被探测海域，然后基于多普勒波束锐化技术对被探测海域进行方位划分，从而得到被探测海域的网格图，其中每个网格称作一个海洋单元；
基于步骤一，确定被探测方位的海洋单元回波的入射方位角对应的多普勒频率；
步骤三：基于步骤二，利用正、负Bragg(布拉格)峰的相对强度建立该海洋单元风向与扩展因子的对应关系；
步骤四：基于步骤三，结合相邻海洋单元风向与扩展因子的对应关系，提取该探测海洋单元的无模糊风向和对应的扩展因子；
步骤五：基于步骤四，建立风速与扩展因子的对应关系；
步骤六：对不同探测方位的海洋单元重复执行步骤二到步骤五，得到该距离门内不同探测方位海洋单元对应的风场；
风场包括无模糊风向和风速；
步骤七：重复执行N遍步骤一到步骤六，得到雷达探测范围内整个海域(所有距离门内不同方位)的海面风场；
N为距离门个数，取值为正整数。
本发明的有益效果为：
本发明解决了传统风场测量工具部署维护困难、无法大面积测量和观测效率低的问题，能够进行全天时、全天候和大面积海域的风场测量，观测效率高；
本发明解决了岸基高频地波雷达探测海域、测量性能受限的问题，船载高频地波雷达利用船载平台的机动性和灵活性不但扩展了雷达的探测范围，而且增大了在复杂海洋环境中的生存能力；
本发明利用单个接收阵元进行海面风场测量，无需利用整个接收阵列，克服了岸基高频地波雷达接收阵列规模较大的局限性，不但节约了雷达系统成本，而且更加适用于甲板狭小的船载平台，为中小船进行海面风场测量提供了可能性；
本发明利用展宽的一阶海面回波多普勒谱得到的扩展因子和风速的关系进行风速提取，克服了利用二阶海面回波多普勒谱进行风速测量的局限性，增加了风速的测量距离。
如图7a和图7b分别为船载高频雷达在探测海域内测量的风向和风速统计图。统计次数大于10的样本用来验证本发明的有益效果，有41.89％雷达测量的风向落在当地气象预报的范围内，有71.79％雷达测量的风速落在当地气象预报的范围内；雷达测量的平均风向和均方根误差分别为15.81°和15.96°，雷达测量的平均风速和均方根误差分别为12.33m/s和2.6m/s。因此利用船载高频地波雷达进行风场测量具有可行性。
附图内容
图1为一种基于单站船载高频地波雷达的海洋表面风场测量方法流程图；
图2a为某一距离单元上相邻海洋单元示意图；
图2b为船载高频地波雷达一阶多普勒展宽谱图；
图3为船载高频地波雷达风向分布图；
图4为同时获取无模糊风向和扩展因子结果图；
图5为确定风速结果图；
图6为船载高频雷达在探测海域内测量的风场分布图，Wind speed表示风速，m/s表示米/秒，km表示千米；
图7a为船载高频雷达在探测海域内测量的风向统计图；
图7b为船载高频雷达在探测海域内测量的风速统计图。
具体实施方式
具体实施方式一：结合图1说明本实施方式，本实施方式所述的一种基于船载高频地波雷达的海洋表面风场测量方法具体过程为：
步骤一：选取某一距离门对应的一阶海面回波多普勒展宽谱；
步骤二：根据雷达系统的距离分辨力等间距划分被探测海域，然后基于多普勒波束锐化技术对被探测海域进行方位划分，从而得到被探测海域的网格图，其中每个网格称作一个海洋单元；
基于步骤一，确定被探测方位的海洋单元回波的入射方位角对应的多普勒频率；
步骤三：基于步骤二，利用正、负Bragg(布拉格)峰的相对强度建立该海洋单元风向与扩展因子的对应关系；
步骤四：基于步骤三，结合相邻海洋单元风向与扩展因子的对应关系，提取该探测海洋单元的无模糊风向和对应的扩展因子；
步骤五：基于步骤四，建立风速与扩展因子的对应关系；
步骤六：对不同探测方位的海洋单元重复执行步骤二到步骤五，得到该距离门内不同探测方位海洋单元对应的风场；
风场包括无模糊风向和风速；
步骤七：重复执行N遍步骤一到步骤六，得到雷达探测范围内整个海域的海面风场(所有距离门内不同方位)；
N为距离门个数，取值为正整数。
具体实施方式二：本实施方式与具体实施方式一不同的是：所述步骤一中选取某一距离门对应的一阶海面回波多普勒展宽谱，具体过程为：
1)利用移动的船载平台及单个接收阵元获得海面回波距离-多普勒谱，船载平台的最大移动速度为其中g是重力加速度，λ是电磁波波长；
2)根据船载高频地波雷达的系统距离分辨力将被探测海域等间距划分为N个距离门，其中ρmax是船载高频地波雷达最大探测距离，ΔR是雷达的系统距离分辨力，是向下取整符号；
3)选择第i个距离门对应的一阶海面回波多普勒展宽谱σi(fd)，其中fd是多普勒频率，i＝1,2,...,N。
其它步骤及参数与具体实施方式一相同。
具体实施方式三：本实施方式与具体实施方式一或二不同的是：所述步骤二中确定被探测方位的海洋单元回波的入射方位角对应的多普勒频率，具体过程为：
设φA是探测海洋单元A处回波的入射方位角，为船载平台运动方向与海洋单元回波的夹角，其中φA∈[0,π]，通过关系式得到探测海洋单元A处回波的入射方位角φA对应的正、负多普勒频率，分别为和
其中是岸基高频地波雷达一阶Bragg频率，v是船载平台的速度。
其它步骤及参数与具体实施方式一或二相同。
具体实施方式四：本实施方式与具体实施方式一至三之一不同的是：所述步骤三中利用正、负Bragg峰的相对强度建立该海洋单元风向与扩展因子的对应关系，具体过程为：
1)定义正、负Bragg峰的相对强度为
其中，
和分别是多普勒频率和对应的Bragg峰强度，为有向波高谱方向因子，ξ为逆风回波与顺风回波的强度比，s为扩展因子，θA是海洋单元A处回波入射方向与船载平台运动方向法线方向的夹角，是海洋单元A处风向与船载平台运动方向法线方向的夹角；χ为或
因此，公式(1)可以写成
2)定义
y为中间变量；
将公式(3)代入公式(2)，获得海洋单元A处的风向与扩展因子的对应关系
其中，海洋单元A处风向与船载平台运动方向法线方向的夹角随着扩展因子s变化，±表示风向的模糊性。
其它步骤及参数与具体实施方式一至三之一相同。
具体实施方式五：本实施方式与具体实施方式一至四之一不同的是：所述逆风回波与顺风回波的强度比ξ为0.004。
其它步骤及参数与具体实施方式一至四之一相同。
具体实施方式六：本实施方式与具体实施方式一至五之一不同的是：所述步骤四中结合相邻海洋单元风向与扩展因子的对应关系，提取该探测海洋单元的无模糊风向和对应的扩展因子，具体过程为：
1)设B为海洋单元A的相邻海洋单元，类似于获得海洋单元A处风向的过程，可以得到海洋单元B处的风向与扩展因子的对应关系
θB是海洋单元B处回波入射方向与船载平台运动方向法线方向的夹角，为海洋单元B处风向与船载平台运动方向法线方向的夹角；
2)对于充分发展的海域，相邻海洋单元之间的风向认为是慢变或不变的，通过比较与扩展因子s相关的的曲线和与扩展因子s相关的的曲线，则两条曲线的交点对应海洋单元A处无模糊风向及唯一的扩展因子s。
其它步骤及参数与具体实施方式一至五之一相同。
具体实施方式七：本实施方式与具体实施方式一至六之一不同的是：所述步骤五中建立风速与扩展因子的对应关系，具体过程为：
1)扩展因子s与风速的关系通过能量传动因子μ得到
其中，μ＝(CD)1/2(4π/gλ)1/2U/κ，U是风速，κ为冯卡曼常数，κ＝0.4，CD＝(0.8+0.065U)×10-3是阻力系数；
公式(6)可以表示为
其中，U*为由雷达工作频率确定的值；
α1＝0.065×10-3；
β1＝0.8×10-3；
γ1＝0；
y1、p1、q1、β1、α1、γ1、δ1为中间变量；
2)观察公式(7)，当U＞U*时，风速U可以由扩展因子s唯一确定，即
其中，
α2＝0.065×10-3；
β2＝0.8×10-3；
γ2＝0；
y2、q2、p2、β2、α2、γ2、δ2为中间变量。
其它步骤及参数与具体实施方式一至六之一相同。
采用以下实施例验证本发明的有益效果：
实施例一：
实验中采用的数据来源于哈尔滨工业大学电子所：雷达发射频率为6.45MHz，带宽为50kHz，探测距离范围为6-120km，探测方位范围为相对于船载平台航行方向53.5°-151.1°，距离分辨率为3km，相干积累时间为129s，平均航速为4.67m/s，平均航向为北偏西168°。由当地的气象预报统计，在船载高频地波雷达从北到南的探测海域内，风向由东北偏北风27.5°到北风10.6°慢变，风速由13.8m/s到8m/s慢变。如图2a、图2b、图3、图4、图5；
图6为船载高频雷达在探测海域内测量的风场分布图。从图中可以看出大部分的风向是东北偏北风或北风；且在从北到南的船载高频雷达的探测海域内，风向由东北偏北风到北风慢变，风速逐渐减小。因此船载高频地波雷达的测量结果展示了与当地气象预报良好的一致性。
图7a和图7b分别为船载高频雷达在探测海域内测量的风向和风速统计图。统计次数大于10的样本用来验证本发明的有益效果，有41.89％雷达测量的风向落在当地气象预报的范围内，有71.79％雷达测量的风速落在当地气象预报的范围内；雷达测量的平均风向和均方根误差分别为15.81°和15.96°，雷达测量的平均风速和均方根误差分别为12.33m/s和2.6m/s。因此利用船载高频地波雷达进行风场测量具有可行性。
本发明还可有其它多种实施例，在不背离本发明精神及其实质的情况下，本领域技术人员当可根据本发明作出各种相应的改变和变形，但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。
Claims (7)
1.一种基于船载高频地波雷达的海洋表面风场测量方法，其特征在于，所述方法的具体过程为：
步骤一：选取某一距离门对应的一阶海面回波多普勒展宽谱；
步骤二：根据雷达系统的距离分辨力等间距划分被探测海域，然后基于多普勒波束锐化技术对被探测海域进行方位划分，从而得到被探测海域的网格图，其中每个网格称作一个海洋单元；
基于步骤一，确定被探测方位的海洋单元回波的入射方位角对应的多普勒频率；
步骤三：基于步骤二，利用正、负Bragg峰的相对强度建立该海洋单元风向与扩展因子的对应关系；
所述正、负Bragg峰为正、负布拉格峰；
步骤四：基于步骤三，结合相邻海洋单元风向与扩展因子的对应关系，提取该探测海洋单元的无模糊风向和对应的扩展因子；
步骤五：基于步骤四，建立风速与扩展因子的对应关系；
步骤六：对不同探测方位的海洋单元重复执行步骤二到步骤五，得到该距离门内不同探测方位海洋单元对应的风场；
风场包括无模糊风向和风速；
步骤七：重复执行N遍步骤一到步骤六，得到雷达探测范围内整个海域的海面风场；
N为距离门个数，取值为正整数。
2.根据权利要求1所述一种基于船载高频地波雷达的海洋表面风场测量方法，其特征在于：所述步骤一中选取某一距离门对应的一阶海面回波多普勒展宽谱，具体过程为：
1)利用移动的船载平台及单个接收阵元获得海面回波距离-多普勒谱，船载平台的最大移动速度为其中g是重力加速度，λ是电磁波波长；
2)根据船载高频地波雷达的系统距离分辨力将被探测海域等间距划分为N个距离门，其中ρmax是船载高频地波雷达最大探测距离，ΔR是雷达的系统距离分辨力，是向下取整符号；
3)选择第i个距离门对应的一阶海面回波多普勒展宽谱σi(fd)，其中fd是多普勒频率，i＝1,2,...,N。
3.根据权利要求2所述一种基于船载高频地波雷达的海洋表面风场测量方法，其特征在于：所述步骤二中确定被探测方位的海洋单元回波的入射方位角对应的多普勒频率，具体过程为：
设φA是探测海洋单元A处回波的入射方位角，为船载平台运动方向与海洋单元回波的夹角，其中φA∈[0,π]，通过关系式得到探测海洋单元A处回波的入射方位角φA对应的正、负多普勒频率，分别为和
其中是岸基高频地波雷达一阶Bragg频率，v是船载平台的速度。
4.根据权利要求3所述一种基于船载高频地波雷达的海洋表面风场测量方法，其特征在于：所述步骤三中利用正、负Bragg峰的相对强度建立该海洋单元风向与扩展因子的对应关系，具体过程为：
1)定义正、负Bragg峰的相对强度为
<mrow>
<msub>
<mi>R</mi>
<mi>A</mi>
</msub>
<mo>=</mo>
<mn>10</mn>
<msub>
<mi>log</mi>
<mn>10</mn>
</msub>
<mrow>
<mo>(</mo>
<mfrac>
<msubsup>
<mi>B</mi>
<mi>A</mi>
<mo>+</mo>
</msubsup>
<msubsup>
<mi>B</mi>
<mi>A</mi>
<mo>-</mo>
</msubsup>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
其中，
和分别是多普勒频率和对应的Bragg峰强度，为有向波高谱方向因子，ξ为逆风回波与顺风回波的强度比，s为扩展因子，θA是海洋单元A处回波入射方向与船载平台运动方向法线方向的夹角，是海洋单元A处风向与船载平台运动方向法线方向的夹角；χ为或
公式(1)可以写成
<mrow>
<msub>
<mi>R</mi>
<mi>A</mi>
</msub>
<mo>=</mo>
<mn>10</mn>
<msub>
<mi>log</mi>
<mn>10</mn>
</msub>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>&xi;</mi>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>&xi;</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>sin</mi>
<mrow>
<mn>2</mn>
<mi>s</mi>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>&theta;</mi>
<mi>A</mi>
</msub>
<mo>-</mo>
<msubsup>
<mi>&alpha;</mi>
<mi>A</mi>
<mo>*</mo>
</msubsup>
</mrow>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&xi;</mi>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>&xi;</mi>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>&theta;</mi>
<mi>A</mi>
</msub>
<mo>-</mo>
<msubsup>
<mi>&alpha;</mi>
<mi>A</mi>
<mo>*</mo>
</msubsup>
</mrow>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mi>s</mi>
</msup>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
2)定义
<mrow>
<mi>y</mi>
<mo>=</mo>
<msup>
<mi>sin</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>&theta;</mi>
<mi>A</mi>
</msub>
<mo>-</mo>
<msubsup>
<mi>&alpha;</mi>
<mi>A</mi>
<mo>*</mo>
</msubsup>
</mrow>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
y为中间变量；
将公式(3)代入公式(2)，获得海洋单元A处的风向与扩展因子的对应关系
<mrow>
<msubsup>
<mi>&alpha;</mi>
<mi>A</mi>
<mo>*</mo>
</msubsup>
<mo>=</mo>
<msub>
<mi>&theta;</mi>
<mi>A</mi>
</msub>
<mo>&PlusMinus;</mo>
<mn>2</mn>
<mi>a</mi>
<mi>r</mi>
<mi>c</mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msqrt>
<mi>y</mi>
</msqrt>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
其中，海洋单元A处风向与船载平台运动方向法线方向的夹角随着扩展因子s变化，±表示风向的模糊性。
5.根据权利要求4所述一种基于船载高频地波雷达的海洋表面风场测量方法，其特征在于：所述逆风回波与顺风回波的强度比ξ为0.004。
6.根据权利要求5所述一种基于船载高频地波雷达的海洋表面风场测量方法，其特征在于：所述步骤四中结合相邻海洋单元风向与扩展因子的对应关系，提取该探测海洋单元的无模糊风向和对应的扩展因子，具体过程为：
1)设B为海洋单元A的相邻海洋单元，得到海洋单元B处的风向与扩展因子的对应关系
<mrow>
<msubsup>
<mi>&alpha;</mi>
<mi>B</mi>
<mo>*</mo>
</msubsup>
<mo>=</mo>
<msub>
<mi>&theta;</mi>
<mi>B</mi>
</msub>
<mo>&PlusMinus;</mo>
<mn>2</mn>
<mi>a</mi>
<mi>r</mi>
<mi>c</mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msqrt>
<mi>y</mi>
</msqrt>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
θB是海洋单元B处回波入射方向与船载平台运动方向法线方向的夹角，为海洋单元B处风向与船载平台运动方向法线方向的夹角；
2)通过比较与扩展因子s相关的的曲线和与扩展因子s相关的的曲线，则两条曲线的交点对应海洋单元A处无模糊风向及唯一的扩展因子s。
7.根据权利要求6所述一种基于船载高频地波雷达的海洋表面风场测量方法，其特征在于：所述步骤五中建立风速与扩展因子的对应关系，具体过程为：
1)扩展因子s与风速的关系通过能量传动因子μ得到
<mrow>
<mi>s</mi>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>0.2</mn>
<msup>
<mrow>
<mo>(</mo>
<mi>&mu;</mi>
<mo>-</mo>
<mn>0.1</mn>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>&mu;</mi>
<mo>></mo>
<mn>0.1</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>2</mn>
</mtd>
<mtd>
<mrow>
<mi>&mu;</mi>
<mo>&le;</mo>
<mn>0.1</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
其中，μ＝(CD)1/2(4π/gλ)1/2U/κ，U是风速，κ为冯卡曼常数，κ＝0.4，CD＝(0.8+0.065U)×10-3是阻力系数；
公式(6)可以表示为
<mrow>
<mi>s</mi>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>0.2</mn>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mi>U</mi>
<mi>&kappa;</mi>
</mfrac>
<msqrt>
<mrow>
<mfrac>
<mrow>
<mn>4</mn>
<mi>&pi;</mi>
</mrow>
<mrow>
<mi>g</mi>
<mi>&lambda;</mi>
</mrow>
</mfrac>
<msub>
<mi>C</mi>
<mi>D</mi>
</msub>
</mrow>
</msqrt>
<mo>-</mo>
<mn>0.1</mn>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>U</mi>
<mo>></mo>
<msup>
<mi>U</mi>
<mo>*</mo>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>2</mn>
</mtd>
<mtd>
<mrow>
<mi>U</mi>
<mo>&le;</mo>
<msup>
<mi>U</mi>
<mo>*</mo>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
其中，U*为由雷达工作频率确定的值；
<mrow>
<msub>
<mi>y</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mroot>
<mrow>
<mo>-</mo>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msqrt>
<mrow>
<msubsup>
<mi>q</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>p</mi>
<mn>1</mn>
<mn>3</mn>
</msubsup>
</mrow>
</msqrt>
</mrow>
<mn>3</mn>
</mroot>
<mo>+</mo>
<mroot>
<mrow>
<mo>-</mo>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msqrt>
<mrow>
<msubsup>
<mi>q</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>p</mi>
<mn>1</mn>
<mn>3</mn>
</msubsup>
</mrow>
</msqrt>
</mrow>
<mn>3</mn>
</mroot>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>9</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>&beta;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&alpha;</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<mfrac>
<msub>
<mi>&gamma;</mi>
<mn>1</mn>
</msub>
<mrow>
<mn>3</mn>
<msub>
<mi>&alpha;</mi>
<mn>1</mn>
</msub>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>27</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>&beta;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&alpha;</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>&beta;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&gamma;</mi>
<mn>1</mn>
</msub>
</mrow>
<mrow>
<mn>6</mn>
<msubsup>
<mi>&alpha;</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<msub>
<mi>&delta;</mi>
<mn>1</mn>
</msub>
<mrow>
<mn>2</mn>
<msub>
<mi>&alpha;</mi>
<mn>1</mn>
</msub>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
α1＝0.065×10-3；
β1＝0.8×10-3；
γ1＝0；
<mrow>
<msub>
<mi>&delta;</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mn>0.01</mn>
<mfrac>
<mrow>
<msup>
<mi>&kappa;</mi>
<mn>2</mn>
</msup>
<mi>g</mi>
<mi>&lambda;</mi>
</mrow>
<mrow>
<mn>4</mn>
<mi>&pi;</mi>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
y1、p1、q1、β1、α1、γ1、δ1为中间变量；
2)当U＞U*时，风速U可以由扩展因子s唯一确定，即
<mrow>
<mi>U</mi>
<mo>=</mo>
<msub>
<mi>y</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<mfrac>
<msub>
<mi>&beta;</mi>
<mn>2</mn>
</msub>
<mrow>
<mn>3</mn>
<msub>
<mi>&alpha;</mi>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
其中，
<mrow>
<msub>
<mi>p</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>9</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>&beta;</mi>
<mn>2</mn>
</msub>
<msub>
<mi>&alpha;</mi>
<mn>2</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<mfrac>
<msub>
<mi>&gamma;</mi>
<mn>2</mn>
</msub>
<mrow>
<mn>3</mn>
<msub>
<mi>&alpha;</mi>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>q</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>27</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>&beta;</mi>
<mn>2</mn>
</msub>
<msub>
<mi>&alpha;</mi>
<mn>2</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>&beta;</mi>
<mn>2</mn>
</msub>
<msub>
<mi>&gamma;</mi>
<mn>2</mn>
</msub>
</mrow>
<mrow>
<mn>6</mn>
<msubsup>
<mi>&alpha;</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<msub>
<mi>&delta;</mi>
<mn>2</mn>
</msub>
<mrow>
<mn>2</mn>
<msub>
<mi>&alpha;</mi>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
α2＝0.065×10-3；
β2＝0.8×10-3；
γ2＝0；
<mrow>
<msub>
<mi>&delta;</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<msup>
<mi>&kappa;</mi>
<mn>2</mn>
</msup>
<mi>g</mi>
<mi>&lambda;</mi>
</mrow>
<mrow>
<mn>4</mn>
<mi>&pi;</mi>
</mrow>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>5</mn>
<mi>s</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mn>0.1</mn>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
y2、q2、p2、β2、α2、γ2、δ2为中间变量。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711173938.3A CN108008392B (zh) | 2017-11-22 | 2017-11-22 | 一种基于船载高频地波雷达的海洋表面风场测量方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711173938.3A CN108008392B (zh) | 2017-11-22 | 2017-11-22 | 一种基于船载高频地波雷达的海洋表面风场测量方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108008392A true CN108008392A (zh) | 2018-05-08 |
| CN108008392B CN108008392B (zh) | 2021-07-06 |
Family
ID=62053064
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711173938.3A Expired - Fee Related CN108008392B (zh) | 2017-11-22 | 2017-11-22 | 一种基于船载高频地波雷达的海洋表面风场测量方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108008392B (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109884337A (zh) * | 2019-03-05 | 2019-06-14 | 哈尔滨工业大学 | 一种利用高频地波雷达探测海面风向的方法 |
| CN110109079A (zh) * | 2019-05-15 | 2019-08-09 | 哈尔滨工业大学 | 海洋回波一阶Bragg峰位置检测方法 |
| CN112782691A (zh) * | 2020-12-24 | 2021-05-11 | 西安空间无线电技术研究所 | 一种基于环扫雷达的海表面风浪流联合探测方法 |
| CN115267717A (zh) * | 2022-07-29 | 2022-11-01 | 哈尔滨工业大学 | 风向夹角反演方法、解模糊方法及存储介质和设备 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4633255A (en) * | 1983-05-16 | 1986-12-30 | The United States Of America As Represented By The Secretary Of The Navy | Method for sea surface high frequency radar cross-section estimation using Dopler spectral properties |
| WO1992021041A1 (en) * | 1991-05-16 | 1992-11-26 | The Commonwealth Of Australia | Oceanographic and meteorological data |
| CN102681033A (zh) * | 2012-04-27 | 2012-09-19 | 哈尔滨工程大学 | 一种基于x波段航海雷达的海面风场测量方法 |
| CN102879766A (zh) * | 2011-07-11 | 2013-01-16 | 哈尔滨工业大学 | 一种高频地波雷达弱目标检测跟踪方法和装置 |
| CN103616690A (zh) * | 2013-12-11 | 2014-03-05 | 哈尔滨工业大学 | 基于船载高频地波超视距雷达的海面风向提取方法 |
| CN104749563A (zh) * | 2015-03-26 | 2015-07-01 | 武汉大学 | 从高频地波雷达海洋回波一阶Bragg谐振峰中提取浪高的方法 |
| US20160157073A1 (en) * | 2013-05-13 | 2016-06-02 | Mitsubishi Electric Corporation | Tsunami monitoring system including transmitting antenna for radiating transmission signal for detecting tsunami as radio wave toward sea |
-
2017
- 2017-11-22 CN CN201711173938.3A patent/CN108008392B/zh not_active Expired - Fee Related
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4633255A (en) * | 1983-05-16 | 1986-12-30 | The United States Of America As Represented By The Secretary Of The Navy | Method for sea surface high frequency radar cross-section estimation using Dopler spectral properties |
| WO1992021041A1 (en) * | 1991-05-16 | 1992-11-26 | The Commonwealth Of Australia | Oceanographic and meteorological data |
| CN102879766A (zh) * | 2011-07-11 | 2013-01-16 | 哈尔滨工业大学 | 一种高频地波雷达弱目标检测跟踪方法和装置 |
| CN102681033A (zh) * | 2012-04-27 | 2012-09-19 | 哈尔滨工程大学 | 一种基于x波段航海雷达的海面风场测量方法 |
| US20160157073A1 (en) * | 2013-05-13 | 2016-06-02 | Mitsubishi Electric Corporation | Tsunami monitoring system including transmitting antenna for radiating transmission signal for detecting tsunami as radio wave toward sea |
| CN103616690A (zh) * | 2013-12-11 | 2014-03-05 | 哈尔滨工业大学 | 基于船载高频地波超视距雷达的海面风向提取方法 |
| CN104749563A (zh) * | 2015-03-26 | 2015-07-01 | 武汉大学 | 从高频地波雷达海洋回波一阶Bragg谐振峰中提取浪高的方法 |
Non-Patent Citations (2)
| Title |
|---|
| G.L.TYLER等: ""Wave directional spectra from synthetic aperture observations of radio scatter"", 《DEEP SEA RESEARCH AND OCEANOGRAPHIC ABSTRACTS》 * |
| 孙明磊: ""船载高频地波雷达海面电磁散射及风向提取研究"", 《中国博士学位论文全文数据库信息科技辑》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109884337A (zh) * | 2019-03-05 | 2019-06-14 | 哈尔滨工业大学 | 一种利用高频地波雷达探测海面风向的方法 |
| CN110109079A (zh) * | 2019-05-15 | 2019-08-09 | 哈尔滨工业大学 | 海洋回波一阶Bragg峰位置检测方法 |
| CN112782691A (zh) * | 2020-12-24 | 2021-05-11 | 西安空间无线电技术研究所 | 一种基于环扫雷达的海表面风浪流联合探测方法 |
| CN115267717A (zh) * | 2022-07-29 | 2022-11-01 | 哈尔滨工业大学 | 风向夹角反演方法、解模糊方法及存储介质和设备 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108008392B (zh) | 2021-07-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108008392B (zh) | 一种基于船载高频地波雷达的海洋表面风场测量方法 | |
| Barrick | HF radio oceanography—A review | |
| Fujii et al. | An overview of developments and applications of oceanographic radar networks in Asia and Oceania countries | |
| Stewart et al. | Radio measurements of oceanic winds at long ranges: An evaluation | |
| Heron et al. | On the application of HF ocean radar to the observation of temporal and spatial changes in wind direction | |
| CN106886021B (zh) | 高分辨率星载sar成像质量提升方法 | |
| Shearman | Radio science and oceanography | |
| Hauser et al. | Overview of the CFOSAT mission | |
| Xie et al. | Measuring ocean surface wind field using shipborne high-frequency surface wave radar | |
| CN102736073B (zh) | 一种通用模式下星载sar距离向模糊度的计算方法 | |
| CN102565772B (zh) | 基于sar子孔径序列图像的海洋动态信息提取方法 | |
| Xie et al. | Ocean surface wind direction inversion using shipborne high-frequency surface wave radar | |
| Liu et al. | A new Doppler model incorporated with free and broken-short waves for coherent S-band wave radar at near-grazing angles | |
| Wan et al. | Research on the applicability of the E spectrum and PM spectrum as the first guess spectrum of SAR wave spectrum inversion | |
| CN104977583A (zh) | 一种基于经验正交分解的x波段雷达海浪反演方法 | |
| Chen et al. | Ocean wave directional spectrum measurement using microwave coherent radar with six antennas | |
| CN102176012A (zh) | 一种二维浅海水下地形合成孔径雷达影像仿真的方法 | |
| CN112162282A (zh) | 一种基于合成孔径雷达海表流速反演方法 | |
| Wyatt | Wave mapping with HF radar | |
| CN114325625A (zh) | 海洋-电离层混沌动力学特征提取方法、装置、介质、应用 | |
| Li et al. | Ocean surface wind and wave monitoring at Typhoon Fung-Wong by HFSWR OSMAR 071 | |
| Nghiem et al. | Wind fields over the Great Lakes measured by the SeaWinds scatterometer on the QuikSCAT satellite | |
| McGregor et al. | Ocean surface currents obtained from microwave sea‐echo Doppler spectra | |
| Shen et al. | Spatial distribution characteristics of surface tidal currents in the southwest of Taiwan Strait | |
| Chu et al. | Algorithm to eliminate the wind direction ambiguity from the monostatic high‐frequency radar backscatter spectra |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210706 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |