CN107831481A - 动态地球坐标系下三部多普勒雷达三维风场反演方法 - Google Patents
动态地球坐标系下三部多普勒雷达三维风场反演方法 Download PDFInfo
- Publication number
- CN107831481A CN107831481A CN201711041700.5A CN201711041700A CN107831481A CN 107831481 A CN107831481 A CN 107831481A CN 201711041700 A CN201711041700 A CN 201711041700A CN 107831481 A CN107831481 A CN 107831481A
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- msup
- mfrac
- coordinate system
- 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.)
- Pending
Links
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
- 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
-
- 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
-
- 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
本发明涉及一种动态地球坐标系下三部多普勒雷达三维风场反演方法,其包括以下步骤:确定动态地球坐标系;采用与动态地球坐标系相应的等经纬度、等海拔高度网格作为反演网格;读取三部多普勒雷达的基数据,将三部多普勒雷达的径向速度和回波强度分别插值到等经纬度、等海拔高度的网格反演网格点上;利用插值后的三部多普勒雷达探测到的目标点的三个方向的径向速度和回波强度计算得到目标点处的三维风速分量,获得降水回波区域的三维风场。本发明采用三部多普勒雷达在地球坐标系下进行三维风场的反演,能够扩大反演区域,提高反演精度。
Description
技术领域
本发明属于大气遥感监测分析技术领域,具体涉及一种动态地球坐标系下三部多普勒雷达三维风场反演方法。
背景技术
多普勒雷达应用多普勒频移原理,通过测量降水粒子反射回来的电磁波与发射电磁波间频率的变化,来测定降水粒子相对于雷达移近或移远的速度。这种速度称为径向速度。由于无法测量降水粒子的真实矢量速度,因此单雷达在观测风场方面受到很大的限制。为了突破单雷达的限制,国内外科研人员发展出双多普勒雷达风场反演技术,即将处于一定距离的两部多普勒雷达所观测的径向速度进行合成,并借助相关假定和方程反演出降水系统内部的三维风场。
目前,国内外普遍采用的是笛卡尔坐标系下的双多普勒雷达风场反演方法。该方法采用笛卡尔坐标系,坐标原点O设置在双雷达中的某一雷达处,或者两部雷达连线的中点上,X轴指向正东方,Y轴指向正北方,Z轴指向天顶,反演网格采用的是等X、Y和Z坐标点。这种方法主要存在以下缺点:1)由于地球是近似球体,地面近似为球面,因此在现有反演技术采用的笛卡尔坐标系下, X-O-Y平面与地球表面不平行,目标点距离原点O越远,X-O-Y平面距离地面就越高。采用该方法反演出的网格点与气象业务科研常用的等经纬度、等海拔高度的网格点不一致,虽然这对小范围天气系统三维风场结构分析影响不大,但风场反演产品难以与数值预报等其他气象信息进行综合分析,直接影响业务应用效果。2)随着我国多普勒天气雷达网的进一步完善,雷达风场拼图分析显得越来越迫切。对于不同的双雷达组合,利用笛卡尔坐标系下的反演技术就需要设置不同的坐标系,导致不同双雷达反演的三维风场产品无法组合到统一坐标系里进行组网分析应用,这在一定程度上影响了雷达观测网建设效益的发挥。
为克服笛卡尔坐标系下双多普勒雷达风场反演方法存在的缺点,有人提出了动态地球坐标系下双多普勒雷达风场反演方法,该方法采用动态地球坐标系,以及等经纬度和等海拔高度的网格。这种方法既便于将反演风场与其他产品进行综合分析,又便于将不同双雷达反演的风场进行组网分析,从而能够有效地克服笛卡尔坐标系下反演方法存在的缺点。然而,这种方法仍存在以下缺点:1) 双多普勒雷达反演的区域较小,且由于双多普勒雷达连线上的区域误差很大,因此无法对双多普勒雷达连线上的区域进行反演。2)双多普勒雷达反演时采用的连续方程近似,垂直速度等存在误差。另外,目前笛卡尔坐标系下的双多普勒雷达和三多普勒雷达风场反演方法由于反演的水平面和地面是不平行的,因此反演出的结果影响业务应用效果。
发明内容
为了解决现有技术存在的上述问题,本发明提供了一种动态地球坐标系下三部多普勒雷达三维风场反演方法。
本发明所采用的技术方案为:动态地球坐标系下三部多普勒雷达三维风场反演方法包括以下步骤:
确定动态地球坐标系;
采用与动态地球坐标系相应的等经纬度、等海拔高度网格作为反演网格;
读取三部多普勒雷达的基数据,将三部多普勒雷达的径向速度和回波强度分别插值到等经纬度、等海拔高度的网格反演网格点上;
利用插值后的三部多普勒雷达探测到的目标点的三个方向的径向速度和回波强度计算得到目标点处的三维风速分量,获得降水回波区域的三维风场。
进一步地,所述动态地球坐标系的坐标原点O设置在地球球心,Z轴随反演目标点的变化而变化。
更进一步地,所述动态地球坐标系中,Z轴由原点O通过目标点P指向雷达站天顶方向,X轴指向目标点P的正东方向,Y轴指向目标点P的正北方向。
进一步地,所述反演网格的经度采用0.01°,纬度采用0.01°,海拔高度采用500m。
进一步地,所述将三部多普勒雷达的径向速度和回波强度分别插值到等经纬度、等海拔高度的网格反演网格点上的步骤中,水平方向采用Cressman距离权重插值方法进行插值,垂直方向采用线性插值方法进行插值。
进一步地,所述步骤获得降水回波区域的三维风场的具体过程为:
在动态地球坐标系下,建立目标点处三维风场的方程组;
利用插值后的回波强度计算得到降水粒子在静止大气中的下落末速度;
利用三部多普勒雷达探测到的目标点的三个方向的径向速度、三部多普勒雷达到目标点的各自距离以及降水粒子在静止大气中的下落末速度,求解目标点处三维风场的方程组,得到目标点处的三维风速分量,获得降水回波区域的三维风场。
进一步地,所述动态地球坐标系下目标点处三维风场的方程组为:
设目标点P(0,0,z)的三维风速分量分别为u、v和w,则三维风场的方程组为:
其中,vr1、vr2和vr3分别表示三部多普勒雷达探测到的目标点P的三个方向的径向速度;(x01,y01,z01)、(x02,y02,z02)和(x03,y03,z03)分别表示三部多普勒雷达的坐标位置;R1、R2和R3分别表示三部多普勒雷达到目标点P的距离,
更进一步地,所述降水粒子在静止大气中的下落末速度wt为:wt=3.8I0.072,其中,I表示插值后得到的回波强度。
由于采用以上技术方案,本发明的有益效果为:本发明动态地球坐标系下多普勒雷达三维风场反演方法采用三部多普勒雷达在地球坐标系下进行三维风场的联合反演,将常用的笛卡尔坐标系变换为动态地球坐标系,将风场反演网格点由等XYZ网格改进为等经纬度和等海拔高度网格。三部多普勒雷达反演能够解决双多普勒雷达连线区域观测自由度不足的问题,可以使得双多普勒雷达连线上的区域得以反演,也可以使得其他部分区域尝试进行三部多普勒雷达反演。三部多普勒雷达反演时不需要连续方程近似,所以反演的垂直速度会比双多普勒雷达更为可靠。本发明能够扩大反演区域,提高反演精度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明一实施例提供的一种动态地球坐标系下三部多普勒雷达三维风场反演方法的流程图;
图2是本发明一实施例提供的动态地球坐标系下三部多普勒雷达三维风场反演示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将对本发明的技术方案进行详细的描述。显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所得到的所有其它实施方式,都属于本发明所保护的范围。
本发明动态地球坐标系下三部多普勒雷达三维风场反演方法采用三部多普勒雷达进行联合反演,将常用的笛卡尔坐标系变换为动态地球坐标系,将风场反演网格点由等XYZ网格改进为等经纬度和等海拔高度网格。三部多普勒雷达反演能够解决双多普勒雷达连线区域观测自由度不足的问题,可以使得双多普勒雷达连线上的区域得以反演,也可以使得其他部分区域尝试进行三部多普勒雷达反演。三部多普勒雷达反演时不需要连续方程近似,所以反演的垂直速度会比双多普勒雷达更为可靠。
如图1所示,本发明提供了一种动态地球坐标系下三部多普勒雷达三维风场反演方法,其包括以下步骤:
S1、确定动态地球坐标系。
动态地球坐标系的坐标原点O设置在地球球心,Z轴随反演目标点的变化而变化。具体地,Z轴由原点O通过目标点P指向雷达站天顶方向,X轴指向目标点P的正东方向,Y轴指向目标点P的正北方向。
S2、采用与动态地球坐标系相应的等经纬度、等海拔高度网格作为反演网格。具体地,反演网格的经度采用0.01°,纬度采用0.01°,海拔高度采用500m。
S3、读取三部多普勒雷达的基数据,将三部多普勒雷达的径向速度和回波强度分别插值到等经纬度、等海拔高度的网格反演网格点上。其中,水平方向采用Cressman距离权重插值方法进行插值,垂直方向采用线性插值方法进行插值。
S4、利用插值后的三部多普勒雷达探测到的目标点的三个方向的径向速度和回波强度计算得到目标点处的三维风速分量,获得降水回波区域的三维风场,其具体过程为:
S41、在动态地球坐标系下,建立目标点处三维风场的方程组;
如图2所示,在步骤S1确定的动态地球坐标系中,设目标点P(0,0,z) 的三维风速分量分别为u、v和w,则三维风场的方程组为:
其中,vr1、vr2和vr3分别表示三部多普勒雷达探测到的目标点P的三个方向的径向速度;(x01,y01,z01)、(x02,y02,z02)和(x03,y03,z03)分别表示三部多普勒雷达的坐标位置;R1、R2和R3分别表示三部多普勒雷达到目标点P的距离,其分别为:
S42、利用插值后的回波强度计算得到降水粒子在静止大气中的下落末速度;
降水粒子在静止大气中的下落末速度wt为:wt=3.8I0.072。
式中,I表示插值后得到的回波强度。
S43、利用三部多普勒雷达探测到的目标点的三个方向的径向速度、三部多普勒雷达到目标点的各自距离以及降水粒子在静止大气中的下落末速度,求解目标点处三维风场的方程组,得到目标点处的三维风速分量,获得降水回波区域的三维风场。
将三部多普勒雷达探测到的目标点P的三个方向的径向速度vr1、vr2和vr3,三部多普勒雷达到目标点P的距离R1、R2和R3,以及降水粒子在静止大气中的下落末速度wt带入三维风场的方程组中,求解得到目标点P(0,0,z)的三维风速分量u、v和w,获得降水回波区域的三维风场。
本发明动态地球坐标系下三部多普勒雷达三维风场反演方法采用三部多普勒雷达在地球坐标系下进行三维风场的反演,能够扩大反演区域,提高反演精度,特别是增加水平散度的精确性。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。
Claims (10)
1.动态地球坐标系下三部多普勒雷达三维风场反演方法,其特征在于,它包括以下步骤:
确定动态地球坐标系;
采用与动态地球坐标系相应的等经纬度、等海拔高度网格作为反演网格;
读取三部多普勒雷达的基数据,将三部多普勒雷达的径向速度和回波强度分别插值到等经纬度、等海拔高度的网格反演网格点上;
利用插值后的三部多普勒雷达探测到的目标点的三个方向的径向速度和回波强度计算得到目标点处的三维风速分量,获得降水回波区域的三维风场。
2.如权利要求1所述的动态地球坐标系下三部多普勒雷达三维风场反演方法,其特征在于,所述动态地球坐标系的坐标原点O设置在地球球心,Z轴随反演目标点的变化而变化。
3.如权利要求2所述的动态地球坐标系下三部多普勒雷达三维风场反演方法,其特征在于,所述动态地球坐标系中,Z轴由原点O通过目标点P指向雷达站天顶方向,X轴指向目标点P的正东方向,Y轴指向目标点P的正北方向。
4.如权利要求1所述的动态地球坐标系下三部多普勒雷达三维风场反演方法,其特征在于,所述反演网格的经度采用0.01°,纬度采用0.01°,海拔高度采用500m。
5.如权利要求1或2或3或4所述的动态地球坐标系下三部多普勒雷达三维风场反演方法,其特征在于,所述将三部多普勒雷达的径向速度和回波强度分别插值到等经纬度、等海拔高度的网格反演网格点上的步骤中,水平方向采用Cressman距离权重插值方法进行插值,垂直方向采用线性插值方法进行插值。
6.如权利要求5所述的动态地球坐标系下三部多普勒雷达三维风场反演方法,其特征在于,所述步骤获得降水回波区域的三维风场的具体过程为:
在动态地球坐标系下,建立目标点处三维风场的方程组;
利用插值后的回波强度计算得到降水粒子在静止大气中的下落末速度;
利用三部多普勒雷达探测到的目标点的三个方向的径向速度、三部多普勒雷达到目标点的各自距离以及降水粒子在静止大气中的下落末速度,求解目标点处三维风场的方程组,得到目标点处的三维风速分量,获得降水回波区域的三维风场。
7.如权利要求1或2或3或4所述的动态地球坐标系下三部多普勒雷达三维风场反演方法,其特征在于,所述步骤获得降水回波区域的三维风场的具体过程为:
在动态地球坐标系下,建立目标点处三维风场的方程组;
利用插值后的回波强度计算得到降水粒子在静止大气中的下落末速度;
利用三部多普勒雷达探测到的目标点的三个方向的径向速度、三部多普勒雷达到目标点的各自距离以及降水粒子在静止大气中的下落末速度,求解目标点处三维风场的方程组,得到目标点处的三维风速分量,获得降水回波区域的三维风场。
8.如权利要求6所述的动态地球坐标系下三部多普勒雷达三维风场反演方法,其特征在于,所述动态地球坐标系下目标点处三维风场的方程组为:
设目标点P(0,0,z)的三维风速分量分别为u、v和w,则三维风场的方程组为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>r</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>01</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
</mfrac>
<mi>u</mi>
<mo>+</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>01</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
</mfrac>
<mi>v</mi>
<mo>+</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>01</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
</mfrac>
<mrow>
<mo>(</mo>
<mi>w</mi>
<mo>-</mo>
<msub>
<mi>w</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>r</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>02</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
</mfrac>
<mi>u</mi>
<mo>+</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>02</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
</mfrac>
<mi>v</mi>
<mo>+</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>02</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
</mfrac>
<mrow>
<mo>(</mo>
<mi>w</mi>
<mo>-</mo>
<msub>
<mi>w</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>r</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>03</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>R</mi>
<mn>3</mn>
</msub>
</mfrac>
<mi>u</mi>
<mo>+</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>03</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>R</mi>
<mn>3</mn>
</msub>
</mfrac>
<mi>v</mi>
<mo>+</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>03</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>R</mi>
<mn>3</mn>
</msub>
</mfrac>
<mrow>
<mo>(</mo>
<mi>w</mi>
<mo>-</mo>
<msub>
<mi>w</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
其中,vr1、vr2和vr3分别表示三部多普勒雷达探测到的目标点P的三个方向的径向速度;(x01,y01,z01)、(x02,y02,z02)和(x03,y03,z03)分别表示三部多普勒雷达的坐标位置;R1、R2和R3分别表示三部多普勒雷达到目标点P的距离,
<mrow>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>01</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>01</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>01</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
</mrow>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
</msup>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>02</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>02</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>02</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
</mrow>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
</msup>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>R</mi>
<mn>3</mn>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>03</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>03</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>03</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
</mrow>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
</msup>
<mo>.</mo>
</mrow>
9.如权利要求8所述的动态地球坐标系下三部多普勒雷达三维风场反演方法,其特征在于,所述降水粒子在静止大气中的下落末速度wt为:wt=3.8I0.072,其中,I表示插值后得到的回波强度。
10.如权利要求7所述的动态地球坐标系下三部多普勒雷达三维风场反演方法,其特征在于,所述动态地球坐标系下目标点处三维风场的方程组为:
设目标点P(0,0,z)的三维风速分量分别为u、v和w,则三维风场的方程组为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>r</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>01</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
</mfrac>
<mi>u</mi>
<mo>+</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>01</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
</mfrac>
<mi>v</mi>
<mo>+</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>01</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
</mfrac>
<mrow>
<mo>(</mo>
<mi>w</mi>
<mo>-</mo>
<msub>
<mi>w</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>r</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>02</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
</mfrac>
<mi>u</mi>
<mo>+</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>02</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
</mfrac>
<mi>v</mi>
<mo>+</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>02</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
</mfrac>
<mrow>
<mo>(</mo>
<mi>w</mi>
<mo>-</mo>
<msub>
<mi>w</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>v</mi>
<mrow>
<mi>r</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>03</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>R</mi>
<mn>3</mn>
</msub>
</mfrac>
<mi>u</mi>
<mo>+</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>03</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>R</mi>
<mn>3</mn>
</msub>
</mfrac>
<mi>v</mi>
<mo>+</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>03</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>R</mi>
<mn>3</mn>
</msub>
</mfrac>
<mrow>
<mo>(</mo>
<mi>w</mi>
<mo>-</mo>
<msub>
<mi>w</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
其中,vr1、vr2和vr3分别表示三部多普勒雷达探测到的目标点P的三个方向的径向速度;(x01,y01,z01)、(x02,y02,z02)和(x03,y03,z03)分别表示三部多普勒雷达的坐标位置;R1、R2和R3分别表示三部多普勒雷达到目标点P的距离,
<mrow>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>01</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>01</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>01</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
</mrow>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
</msup>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>02</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>02</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>02</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
</mrow>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
</msup>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>R</mi>
<mn>3</mn>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>03</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>03</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>03</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
</mrow>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
</msup>
<mo>.</mo>
</mrow>
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711041700.5A CN107831481A (zh) | 2017-10-30 | 2017-10-30 | 动态地球坐标系下三部多普勒雷达三维风场反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711041700.5A CN107831481A (zh) | 2017-10-30 | 2017-10-30 | 动态地球坐标系下三部多普勒雷达三维风场反演方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107831481A true CN107831481A (zh) | 2018-03-23 |
Family
ID=61650158
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711041700.5A Pending CN107831481A (zh) | 2017-10-30 | 2017-10-30 | 动态地球坐标系下三部多普勒雷达三维风场反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107831481A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108693534A (zh) * | 2018-03-27 | 2018-10-23 | 南京恩瑞特实业有限公司 | Nriet x 波段雷达协同组网分析方法 |
CN110288856A (zh) * | 2019-06-21 | 2019-09-27 | 中国民用航空总局第二研究所 | 基于风的精细化预报的航班动态监控系统及方法 |
CN110361722A (zh) * | 2019-07-17 | 2019-10-22 | 四川川大智胜软件股份有限公司 | 一种多雷达信号融合方法及系统 |
CN113009490A (zh) * | 2021-02-20 | 2021-06-22 | 江苏省气象台 | 一种基于高分辨率模式动力约束的雷达三维风场反演方法 |
CN114660591A (zh) * | 2022-05-25 | 2022-06-24 | 成都远望科技有限责任公司 | 一种基于多部天气雷达三维组网生成方法 |
CN114996937A (zh) * | 2022-06-01 | 2022-09-02 | 浙江省气象台 | 一种基于动力统计的台风精细化风场反演方法及其应用 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102508222A (zh) * | 2011-11-18 | 2012-06-20 | 中国科学技术大学 | 一种中高层大气风场反演方法 |
KR101354522B1 (ko) * | 2012-08-23 | 2014-01-22 | 부경대학교 산학협력단 | 변분법과 초기추정법을 이용한 고분해능 레이더 바람장 산출방법 |
CN105467393A (zh) * | 2016-01-06 | 2016-04-06 | 北京无线电测量研究所 | 一种毫米波云雷达三维风场反演方法及系统 |
CN106324601A (zh) * | 2016-08-08 | 2017-01-11 | 北京无线电测量研究所 | 一种基于风廓线雷达组网的三维风场反演方法 |
-
2017
- 2017-10-30 CN CN201711041700.5A patent/CN107831481A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102508222A (zh) * | 2011-11-18 | 2012-06-20 | 中国科学技术大学 | 一种中高层大气风场反演方法 |
KR101354522B1 (ko) * | 2012-08-23 | 2014-01-22 | 부경대학교 산학협력단 | 변분법과 초기추정법을 이용한 고분해능 레이더 바람장 산출방법 |
CN105467393A (zh) * | 2016-01-06 | 2016-04-06 | 北京无线电测量研究所 | 一种毫米波云雷达三维风场反演方法及系统 |
CN106324601A (zh) * | 2016-08-08 | 2017-01-11 | 北京无线电测量研究所 | 一种基于风廓线雷达组网的三维风场反演方法 |
Non-Patent Citations (2)
Title |
---|
罗昌荣等: ""三雷达、双雷达反演降雹超级单体风暴三维风场结构特征研究"", 《气象学报》 * |
罗昌荣等: ""双雷达反演台风外围强带状回波风场结构特征研究"", 《大气科学》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108693534A (zh) * | 2018-03-27 | 2018-10-23 | 南京恩瑞特实业有限公司 | Nriet x 波段雷达协同组网分析方法 |
CN110288856A (zh) * | 2019-06-21 | 2019-09-27 | 中国民用航空总局第二研究所 | 基于风的精细化预报的航班动态监控系统及方法 |
CN110361722A (zh) * | 2019-07-17 | 2019-10-22 | 四川川大智胜软件股份有限公司 | 一种多雷达信号融合方法及系统 |
CN113009490A (zh) * | 2021-02-20 | 2021-06-22 | 江苏省气象台 | 一种基于高分辨率模式动力约束的雷达三维风场反演方法 |
CN114660591A (zh) * | 2022-05-25 | 2022-06-24 | 成都远望科技有限责任公司 | 一种基于多部天气雷达三维组网生成方法 |
CN114996937A (zh) * | 2022-06-01 | 2022-09-02 | 浙江省气象台 | 一种基于动力统计的台风精细化风场反演方法及其应用 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107831481A (zh) | 动态地球坐标系下三部多普勒雷达三维风场反演方法 | |
CN107843895B (zh) | 一种双多普勒雷达三维风场反演方法 | |
CN109298388B (zh) | 基于方位信息的超视距目标地理坐标直接估计方法 | |
CN104765032B (zh) | 岸基阵列gnss反射信号潮汐与多波浪参数综合探测系统 | |
CN109100723B (zh) | 基于多普勒天气雷达数据的高空风反演方法 | |
CN109001846B (zh) | 一种复杂地形条件下s波段与x波段雷达组网测雨方法 | |
CN102928860B (zh) | 基于局部定位信息提高gps定位精度的方法 | |
CN106019217A (zh) | 基于aoa的二维无线传感器网络半定规划定位方法 | |
CN103841641B (zh) | 一种基于到达角度和Gossip算法的无线传感器网络分布式协作定位方法 | |
CN109884730B (zh) | 一种基于深海水下滑翔机平台的虚拟阵被动探测方法 | |
CN204556822U (zh) | 岸基阵列gnss反射信号潮汐与多波浪参数综合探测系统 | |
Tiryakioğlu | Geodetic aspects of the 19 May 2011 Simav earthquake in Turkey | |
CN109612441A (zh) | 一种基于最优分类和矢量平均的波浪方向浮标主波向计算方法 | |
CN105758405A (zh) | 一种适用于高海况的uuv跟踪母船航迹规划方法 | |
Itay et al. | Lagrangian kinematic criterion for the breaking of shoaling waves | |
CN110991087A (zh) | 一种基于多入射角组网sar卫星数据的风场反演方法及系统 | |
CN112782705B (zh) | 一种超声声层析流场成像方法及其装置 | |
CN105388470B (zh) | 一种估计目标运动参数的方法 | |
CN107871323A (zh) | 一种基于全固态雷达的海冰运动信息获取方法 | |
CN206459611U (zh) | 一种基于北斗卫星的超高层建筑高精度变形监测系统 | |
CN108061666A (zh) | 一种输电塔损伤识别方法 | |
CN107942316A (zh) | 一种基于多波束声呐波束输出信号的水中浓悬浮体运动速度估计方法 | |
CN104297753A (zh) | 一种基于自适应缩减算子的导航雷达图像反演海面风向方法 | |
CN109683201A (zh) | 一种基于轴频电场的船舶跟踪定位方法 | |
CN207215165U (zh) | 一种船舶数据共享系统 |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180323 |