CN104035096B - 一种基于多普勒天气雷达的垂直风廓线非线性反演方法 - Google Patents
一种基于多普勒天气雷达的垂直风廓线非线性反演方法 Download PDFInfo
- Publication number
- CN104035096B CN104035096B CN201410254985.0A CN201410254985A CN104035096B CN 104035096 B CN104035096 B CN 104035096B CN 201410254985 A CN201410254985 A CN 201410254985A CN 104035096 B CN104035096 B CN 104035096B
- Authority
- CN
- China
- Prior art keywords
- radar
- vertical wind
- doppler
- rvd
- data
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 54
- 230000008569 process Effects 0.000 claims abstract description 7
- 230000009021 linear effect Effects 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 abstract description 7
- 238000004422 calculation algorithm Methods 0.000 description 10
- 238000004458 analytical method Methods 0.000 description 8
- 230000008859 change Effects 0.000 description 6
- 238000013178 mathematical model Methods 0.000 description 6
- 239000011159 matrix material Substances 0.000 description 4
- 230000006870 function Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 230000009022 nonlinear effect Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 241000700199 Cavia porcellus Species 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000002459 sustained effect Effects 0.000 description 1
- 238000003786 synthesis reaction Methods 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
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
-
- 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)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及一种基于多普勒天气雷达的垂直风廓线非线性反演方法,步骤如下:任意选取某一高度的实际水平风场,并用泰勒多项式展开的形式来表示;将雷达的速度观测用rVd的形式来表示;将步骤1中多项式代入步骤2中的rVd的表达式中,得到非线性风场下场rVd的理论表达式;计算雷达每个观测点的多普勒数据Vd和距离雷达径向距离r的乘积,得到rVd的实际物理量;将步骤4中观测得到的实际物理量rVd代入步骤3的rVd的理论表达式中,对这一数学表达式进行多项式拟合从而计算出该高度上的水平风速;选取不同高度的数据,重复步骤1至步骤5的过程,即可得到垂直风廓线。该方法能够更好地进行水平风场以及垂直风廓线的反演,改变了以往方法由于资料分布不均匀而造成地计算不稳定的问题。
Description
技术领域
本发明涉及一种基于多普勒天气雷达的垂直风廓线反演方法,尤其是一种针对实际风场包含大量缺测以及强非线性成分时的垂直风廓线反演方法,属于大气科学(遥感资料分析)研究领域,用于多普勒天气雷达的垂直风廓线反演及临近预报。
背景技术
垂直风廓线(水平风速随高度的变化)的反演是多普勒天气雷达最重要的业务应用之一。以往的反演方法一般包括以下两步:
①针对实际风场的空间分布特征选择合理的数学模型对实际风场进行建模;
②使用多普勒天气雷达观测到的径向速度来拟合相应数学模型,从而得到需要的水平风信息。使用相同的方法对不同高度的数据进行计算,即可得到垂直风廓线。
最早基于多普勒天气雷达的垂直风廓线反演算法(VAD,距离方位显示)由(Browning and Wexler1968)在1968年提出,其方法的核心是假定实际风场在同一高度平面上呈线性变化。在这样的线性假设下,雷达观测到的多普勒速度Vd在一个距离圈上是方位角α的三角函数:Vd=usinαcosβ+vcosαcosβ+wsinβ (1)
通过选取某一距离圈(同一仰角和径向距离,但方位角不同的数据)上的数据进行计算,即可得到水平风u和v,再通过计算不同高度上的距离圈的数据,即可得到垂直风廓线。(Waldteufel and Cobin1979)对以上的方法进行了改进,不再使用单距离圈的资料进行拟合,而是使用一个雷达体积扫描的所有数据来进行相应的计算。(Tabary et al.2001)进一步使用使用径向速度在方位方向的梯度来进行拟合,从而克服了多普勒速度折叠的问题。
虽然VAD方法从提出以来得到了以上的改进,但其方法的基本假设,即实际风场在空间呈线性变化,并没有改变。这一假设虽然在一般的情况下可以近似成立,但是实际天气过程的风场可能包含很强的非线性性。因此基于以上线性假设而计算得到的垂直风廓线可能存在较大的误差。(Caya and Zawadzki1992)在VAD方法的基础之上,提出了一个可以分析非线性风场的垂直风廓线反演方法(NVAD)。这个方法的实现分成两个步骤,第一步的分析和VAD的方法相似,即同样使用三角函数来对一个距离圈上的多普勒速度数据进行建模。唯一的区别是和实际风场对应的数学模型了包含了非线性项的贡献。在包含了非线性成分的贡献后(这里以3阶的非线性风场为例),多普勒速度在一个距离圈上的表达式为:
在式(2)里,三角函数的系数不再与雷达所在位置的水平风的直接对应,而同时依赖与水平距离r。因此,NVAD第二步分析以第一步计算的结果为基础,对不同距离圈的cosα和sinα对应的系数进行一个多项式的拟合来得到雷达上空的水平风场变化。(Caya andZawadzki1992)的方法虽然从理论上解决了雷达风廓线反演中的非线性污染问题。但在实际应用当中还有很大的局限性。
(Caya and Zawadzki1992)的方法第一步使用了和VAD计算相同的算法,当计算的距离圈上出现大范围资料缺测时,由于“过拟合”问题该方法容易出现严重的计算偏差。同时,由于该方法计算的第二个步骤以第一个步骤的结果为基础,第一个计算步骤中出现的错误结果会进一步污染第二步的计算结果。另外,由于NVAD的计算需要对雷达观测到的每一个距离圈独立的进行,因而难以通过对中间计算结果进行分析的方法来去除出现严重计算偏差的结果。最后,从(2)可以直观地看出,NVAD方法所包含的数学模型风场复杂,在实时的雷达分析应用中并不方便。
发明内容
本发明针对现有基于天气雷达反演垂直风廓线技术的不足,提出一种新的垂直风廓线的反演方法。其首要技术目的是改变现有技术在包含缺测和噪音的观测中不能有效处理非线性风场的问题。
本发明的方法不以雷达直接观测到的多普勒速度Vd来进行计算,而是使用多普勒速度Vd和对应数据点到雷达径向距离r的乘积rVd为基本的分析量来进行建模和反演。通过使用这一新的合成观测量,实际的风场可以简单地表示成一个二元的高次多项式。雷达所在位置的水平风速u和v可以通过对以上的多项式分析而得到。同时,由于新的分析方法只有一个拟合的过程,因而在拟合当中的数据量变得很大,这样可以显著提高拟合的稳定性。当实际的观测包含较大的缺测时,仍然可以提供准确的计算结果。
为了实现以上的技术目的,本发明将采取以下的技术方案:将基于多普勒天气雷达的垂直风廓线反演分成“非线性风场建模”和“反演”两个部分。
本发明的基于天气雷达的垂直风廓线反演方法,具体步骤如下:
1.任意选取某一高度的实际水平风场,并用泰勒多项式展开的形式来表示;
2.将雷达的速度观测用rVd的形式来表示;
3.将步骤1中多项式代入步骤2中的rVd的表达式中,得到非线性风场下场rVd的理论表达式;
4.计算雷达每个观测点的多普勒数据Vd和距离雷达径向距离r的乘积,得到rVd的实际物理量;
5.将步骤4中观测得到的实际物理量rVd代入步骤3的rVd的理论表达式中,对这一数学表达式进行多项式拟合从而计算出该高度上的水平风速;
6.选取不同高度的数据,重复步骤1至步骤5的过程,即可得到垂直风廓线。
上述步骤1至步骤3为“非线性风场建模”步骤,步骤4至步骤5为“反演”步骤。
本发明的基于天气雷达的垂直风廓线反演方法较现有的反演方法,其优点主要有如下三个方面:
1、本发明采用多普勒速度和径向距离的乘积来分析雷达的速度观测。使用这一分析量使得实际的非线性风场可以用简洁的数学模型来表达。与现有技术相比,这种新的数学模型能够更好地进行水平风场以及垂直风廓线的反演。
2、本发明所述的风场非线性估计方法,可以后验的对实际风场中包含的非线性性进行估计,从而选择和实际风场最为吻合的数学模型。
3、本发明采用的风廓线数学模型可以使用一定垂直高度内的雷达资料进行数学分析,因而从根本是改变了以往方法由于资料分布不均匀而造成地计算不稳定的问题。
附图说明
图1为用于检验本方法所使用的在不同高度的模拟风场图;
图2为理想情况下,虚拟雷达在图1所示的风场里观测到的多普勒速度图;
图3为理想情况下,现有反演算法与本发明算法的效果对比图;
图4为在图2理想多普勒数据之上,增加了噪音和缺测数据后的多普勒速度图;
图5为模拟实际情况下(数据包含缺测和噪音),现有反演算法与本发明算法的对比图;
图6为一次实际天气过程的现有反演算法与本发明算法的对比图。
具体实施方式
以下将结合附图详细地说明本发明的技术方案。
一、建模
实际风场的水平分量在雷达中心处的泰勒展开可以表示为:
多普勒径向速度Vd和径向距离r的乘积rVd表示为:
rVd=ux+vy+wz (4)
将(3)带入(4)式,rVd可以和实际风场的对应关系为:
上式可以简化为:
对比(6)式中各项系数cij与(5)式中的实际风场的参数可以得到:
c10=u0,c01=v0 (7)
通过求解式(6)中的系数,即可得到不同高度的水平风从而得到垂直风廓线。
二、反演
为了方便说明,式(6)可以写成如下的矩阵的形式:
Ac+∈=b (8)
其中矩阵A为系数矩阵,其数值通过数据点所在位置的水平坐标x,y和泰勒展开的阶数可以直接确定。向量b为所有数据点多普勒径向速度Vd和径向距离r的乘积rVd。由于雷达观测每给数据点的水平坐标x,y径向距离r和多普勒径向速度Vd都为已知。因此通过这些观测量计算出A,b,再使用标准的矩阵求解算法即可得到c。由式(7)可以看出,这样就求得了水平风u,v。再使用不同高度上的数据进行计算,即可得到垂直风廓线。
三、本发明反演方法的有效性论证
为了验证本发明内容的有效性和可行性,这里使用模拟和实际数据并以附图的形式来进行说明。由于实际天气观测时并不知道观测变量的真实值,因此使用了模拟的观测资料来检验新方法的正确性。同时,我们使用了业务布网的S波段多普勒雷达的实测资料来对比新旧方法在实际应用中的优劣。
图1显示了用于检验本方法所使用的在不同高度的模拟风场。其中十字标记表示模拟雷达的位置。图3中带圆点的线条(实线为东西风,虚线为南北风)表示了雷达所在位置不同高度水平风场的真实值。图2显示了理想情况下,虚拟雷达在图1所示的风场里观测到的多普勒速度。分别使用前面所示的NVAD方法和本专利提出的新方法计算得到的垂直风廓线如图3所示。其中左为NVAD的结果,右侧为本专利方法得到的结果。通过对比可以看出,在理想情况下,这两种方法是一致的,这说明了本专利方法在理论上的正确性。
实际雷达观测到的多普勒速度由于各种原因会包含各种缺测和噪音,以上例子里的理想数据在实际观测中是不可能存在的。为了模拟这些在实际观测中的问题,我们在图2中的理想多普勒数据之上,人为地增加了噪音和缺测数据。在增加了这些数据问题后的多普勒速度如图4所示。以图4所示的多普勒数据作为输入,两个不同方法计算得到的的结果如图5所示。从图中可以看出,NVAD方法计算出现了严重的偏差,而本发明提出的方法几乎没有受到噪音和缺测资料的影响。
图6显示了一个实际天气过程中两个方法计算结果的对比。其中上图为NVAD方法,下图为本专利提出的方法。图中重叠及缺少部分深色的风标为独立的探空观测作为验证。从图中可以看出,NVAD方法出现了明显的风场计算误差,而本专利的方法计算出的垂直风廓线显示出了很好的一致性,同时和独立的探空资料也非常的接近。
Claims (2)
1.一种基于多普勒天气雷达的垂直风廓线非线性反演方法,其特征在于包含如下步骤:
①任意选取某一高度的实际水平风场,并用泰勒多项式展开的形式来表示;
②将雷达的速度观测用rVd的形式来表示;
③将步骤1中多项式代入步骤2中的rVd的表达式中,得到非线性风场下场rVd的理论表达式;
④计算雷达每个观测点的多普勒数据Vd和距离雷达径向距离r的乘积,得到rVd的实际物理量;
⑤将步骤4中观测得到的实际物理量rVd代入步骤3的rVd的理论表达式中,对这一数学表达式进行多项式拟合从而计算出该高度上的水平风速;
⑥选取不同高度的数据,重复步骤1至步骤5的过程,即可得到垂直风廓线。
2.如权利要求1所述的一种基于多普勒天气雷达的垂直风廓线非线性反演方法,其特征在于所述步骤1中的泰勒多项式展开形式为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410254985.0A CN104035096B (zh) | 2014-06-06 | 2014-06-06 | 一种基于多普勒天气雷达的垂直风廓线非线性反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410254985.0A CN104035096B (zh) | 2014-06-06 | 2014-06-06 | 一种基于多普勒天气雷达的垂直风廓线非线性反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104035096A CN104035096A (zh) | 2014-09-10 |
CN104035096B true CN104035096B (zh) | 2017-05-03 |
Family
ID=51465933
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410254985.0A Expired - Fee Related CN104035096B (zh) | 2014-06-06 | 2014-06-06 | 一种基于多普勒天气雷达的垂直风廓线非线性反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104035096B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104569981B (zh) * | 2015-01-28 | 2017-02-22 | 中国科学院大气物理研究所 | 协同自适应观测方法 |
CN105022036B (zh) * | 2015-08-26 | 2018-02-23 | 成都信息工程大学 | 风廓线雷达风速测定方法 |
CN109100723B (zh) * | 2018-07-25 | 2022-05-27 | 南京信息工程大学 | 基于多普勒天气雷达数据的高空风反演方法 |
CN109283535B (zh) * | 2018-11-16 | 2022-12-20 | 北京无线电测量研究所 | 非均匀风场的探测方法 |
CN110780299A (zh) * | 2019-10-24 | 2020-02-11 | 上海眼控科技股份有限公司 | 散度场的获取方法、装置、计算机设备和存储介质 |
CN112965084B (zh) * | 2021-01-28 | 2021-10-22 | 中国人民解放军国防科技大学 | 基于激光雷达的机场风场特征探测方法、装置和设备 |
CN112748434B (zh) * | 2021-02-07 | 2023-04-11 | 江西师范大学 | 一种基于多项式拟合的多普勒天气雷达速度退模糊方法 |
CN118091666B (zh) * | 2024-04-29 | 2024-06-28 | 江西省气象科学研究所 | 深度学习风廓线雷达全波束的三维风场反演方法及系统 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101988963B (zh) * | 2010-04-19 | 2012-09-26 | 南京恩瑞特实业有限公司 | 利用风廓线雷达获取三维风场的方法 |
RU2477490C1 (ru) * | 2011-07-06 | 2013-03-10 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Московский государственный университет приборостроения и информатики" (МГУПИ) | Способ определения вертикального профиля ветра в атмосфере |
CN102426352B (zh) * | 2011-10-17 | 2013-11-20 | 中国人民解放军理工大学气象学院 | 风廓线雷达风计算方法 |
-
2014
- 2014-06-06 CN CN201410254985.0A patent/CN104035096B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN104035096A (zh) | 2014-09-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104035096B (zh) | 一种基于多普勒天气雷达的垂直风廓线非线性反演方法 | |
Bechmann et al. | The Bolund experiment, part II: blind comparison of microscale flow models | |
CN108763825B (zh) | 一种模拟复杂地形的风场的数值模拟方法 | |
CN103353923B (zh) | 基于空间特征分析的自适应空间插值方法及其系统 | |
An et al. | Sensitivity of inflow boundary conditions on downstream wind and turbulence profiles through building obstacles using a CFD approach | |
CN105224715A (zh) | 一种山区地貌下强风三维脉动风场综合模拟方法 | |
CN104597430B (zh) | 机载气象雷达低空风切变三维回波仿真方法 | |
CN103472502B (zh) | 动态展示区域空气质量和气象场的方法 | |
CN108959742A (zh) | 大跨越输电塔线体系气动弹性模型设计方法 | |
CN105842676A (zh) | 一种雷达盲区分析方法及装置 | |
CN104535981A (zh) | 海杂波Pareto幅度分布参数的双分位点估计方法 | |
CN104573249A (zh) | 基于时变arma模型的非平稳风速模拟方法 | |
CN107885964A (zh) | 一种顾及复杂地形的风能cfd模拟方法 | |
CN103700091A (zh) | 基于多尺度低秩分解且结构信息敏感的图像显著性物体检测方法 | |
CN111859748A (zh) | 一种基于垂向混合坐标的海洋内波模拟方法 | |
CN103366045A (zh) | 基于格子Boltzmann的流体可视化仿真方法 | |
CN115329691A (zh) | 一种基于cfd与gis的超大城市风环境模拟方法 | |
CN101750616B (zh) | 植被风阻的测量方法及系统 | |
Burlando et al. | A simple and efficient procedure for the numerical simulation of wind fields in complex terrain | |
CN104951798A (zh) | 基于emd的lssvm非平稳脉动风速预测方法 | |
CN107085646A (zh) | 一种输电线路下地面场强测量方法 | |
CN103870699B (zh) | 基于双层异步迭代策略的水动力学洪水演进模拟方法 | |
Sikirić et al. | Coupling of the Regional Ocean Modeling System (ROMS) and wind wave model | |
Yudin | Comparison of FDM and FEM models for a 2D gravity current in the atmosphere over a valley | |
CN103678788A (zh) | 基于曲面论的空间数据插值与曲面拟合方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into 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: 20170503 |
|
CF01 | Termination of patent right due to non-payment of annual fee |