[发明专利]动态地球坐标系下三部多普勒雷达三维风场反演方法

说明书

技术领域

本发明属于大气遥感监测分析技术领域，具体涉及一种动态地球坐标系下三部多普勒雷达三维风场反演方法。

背景技术

多普勒雷达应用多普勒频移原理，通过测量降水粒子反射回来的电磁波与发射电磁波间频率的变化，来测定降水粒子相对于雷达移近或移远的速度。这种速度称为径向速度。由于无法测量降水粒子的真实矢量速度，因此单雷达在观测风场方面受到很大的限制。为了突破单雷达的限制，国内外科研人员发展出双多普勒雷达风场反演技术，即将处于一定距离的两部多普勒雷达所观测的径向速度进行合成，并借助相关假定和方程反演出降水系统内部的三维风场。

目前，国内外普遍采用的是笛卡尔坐标系下的双多普勒雷达风场反演方法。该方法采用笛卡尔坐标系，坐标原点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距离权重插值方法进行插值，垂直方向采用线性插值方法进行插值。

进一步地，所述步骤获得降水回波区域的三维风场的具体过程为：

在动态地球坐标系下，建立目标点处三维风场的方程组；

利用插值后的回波强度计算得到降水粒子在静止大气中的下落末速度；

利用三部多普勒雷达探测到的目标点的三个方向的径向速度、三部多普勒雷达到目标点的各自距离以及降水粒子在静止大气中的下落末速度，求解目标点处三维风场的方程组，得到目标点处的三维风速分量，获得降水回波区域的三维风场。

进一步地，所述动态地球坐标系下目标点处三维风场的方程组为：