一种毫米波云雷达三维风场反演方法及系统的制作方法
【技术领域】
[0001] 本发明涉及雷达技术领域，尤其涉及一种毫米波云雷达三维风场反演方法及系 统。
【背景技术】
[0002] 风场信息是大气探测的主要气象因子之一。大型活动的气象保障、机场气象保障、 天气预报等领域都对风场的精细观测有着强烈需求。除了风速风向仪、风廓线雷达等仪器 直接测风外，利用空中漂浮的示踪颗粒(如云滴、雨滴、烟雾等）间接测量风场信息也是一类 重要的方法。正因如此，目前业务运行的天气雷达已将风场反演信息作为一个常规的输出 产品。
[0003] 然而天气雷达所反演三维风场存在许多不足之处，主要表现在:第一，天气雷达反 演风场所依赖的示踪物是雨滴及较强的云滴粒子，这些粒子在天空中出现的概率并不是十 分频繁，这就导致了在大部分时间里有效测风数据的空白；第二，由于观测的对象不同，天 气雷达通常的俯仰角设置较低，这使得圆周扫描时水平半径达到数百公里，很难满足线性 均匀的风场反演假设，造成较大的风场反演误差;第三，受失踪对象的限制，天气雷达所能 够反演风场的高度范围较小。
[0004] 相比天气雷达，近些年蓬勃发展的波长更短的毫米波云雷达能够探测到空中更弱 的粒子，因此在三维风场反演方面拥有更好的潜力。毫米波云雷达主要用于云、雾、弱降水 等目标的探测，利用小粒子对电磁波的散射作用，连续测量站点上空气象目标的回波信号， 获取高时空分辨率的气象目标信息，具有全方位、全天候的观测能力，广泛应用于大气科学 研究、人工影响天气、云自动化观测、机场气象保障、军事气象保障等方面和领域。
[0005] 由于自身体制(波束窄、相关性弱)及被测目标的特性(强度弱、对环境及杂波扰动 敏感），毫米波雷达的探测数据规模大、单个数据的随机性强、连续性及稳定性差。这些特点 也给三维风场反演的求解方法设置及反演准确性造成了一定的困难。
【发明内容】
[0006] 本发明所要解决的技术问题是针对现有技术的不足，提供一种毫米波云雷达三维 风场反演方法及系统，解决了无法有效且高精度地利用毫米波云雷达得到三维风场产品的 问题。本发明解决上述技术问题的技术方案如下：
[0007] -种毫米波云雷达三维风场反演方法，包括以下步骤：
[0008] 步骤1:控制毫米波云雷达按预设俯仰角的PPI扫描模式进行气象目标探测，并生 成目标探测数据；
[0009] 步骤2:对所述目标探测数据进行信号处理并输出原始径向速度数据；
[0010] 步骤3:采集所述信号处理模块输出的原始径向速度数据，并将其进行打包处理；
[0011] 步骤4:将打包处理之后的原始径向速度数据进行保存；
[0012] 步骤5:读取保存的原始径向速度数据并对其进行用以消除毫米波雷达回波信号 的不稳定性的预处理，得到可供三维风场反演的径向速度数据；
[0013] 步骤6:读取所述可供三维风场反演的径向速度数据，并在预设求解准则下构建反 演方程以对对应的所述目标探测数据进行三维风场反演，得到三维风场速度数据。
[0014] 本发明的有益效果是:利用毫米波云雷达进行了三维风场的反演，充分利用了毫 米波段雷达对弱目标粒子探测能力强的优点，增强了测风数据的可靠性和有效性，丰富了 毫米波测云雷达的数据产品类型。另外，从数据预处理和求解准则两个方面建立反演方程， 从而保证了利用原始径向数据反演三维风场的有效性与准确性。
[0015] 在上述技术方案的基础上，本发明还可以做如下改进。
[0016] 在上述技术方案的基础上，本发明还可以做如下改进。
[0017] 进一步，毫米波云雷达的水平扫描速度为1~2度/秒，预设俯仰角的最大值不超过 82度。
[0018] 采用上述进一步方案的有益效果是合理的扫描速度既保证了波束对同一目标有 足够的驻留时间，又不会使反演过程所需的时间过长。同时俯仰角的合理设置则既提高了 对较高空处风场反演的能力，又避免了因俯仰角过大而水平风场的测量误差较大。
[0019] 进一步，所述步骤5的预处理具体实现为：
[0020] 步骤5.1:对所述原始径向速度数据按方位角的大小进行升序排列；
[0021]步骤5.2:设置一个径向速度数据的预设阈值以去除所述原始径向速度数据中的 奇异点；
[0022]步骤5.3:使用5X1的滤波窗口按方位角排序并对去除了奇异点的原始径向速度 数据进行平滑滤波，得到相关参数数据。
[0023]采用上述进一步方案的有益效果是通过预处理将原始径向速度数据中的有用信 息提取出来，增强其可用性，保证了后续的求解方程可解及结果准确。
[0024]进一步，所述步骤5.3具体实现为：
[0025] 步骤5.3.1:设扫描控制模块进行气象目标探测得到的总观测点数为M，总迭代次 数为N;
[0026] 步骤5.3.2:从1开始，对于第i点，计算其周围5个观测点的径向速度数据的均值m;
[0027] 步骤5.3.3:若其中一个所述观测点与所述均值m的绝对值之差> =0.5dBZ，则去除 所述第i点，否则保留所述第i点；
[0028] 步骤5.3.4:计算所述观测点的下一个点，直至计算到Μ;
[0029] 步骤5.3.5:进行下一次迭代，直至迭代到次数Ν。
[0030] 采用上述进一步方案的有益效果是通过利用单个径向速度数据在其方位上相邻 数据的信息，对因谱识别误差、环境杂波扰动等原因导致的误差较大的数据进行去除，保证 了后续风场反演的精度。
[0031]进一步，所述步骤6具体实现为：
[0032]步骤6.1:假设扫描控制模块进行气象目标探测得到的不同观测点的方位角分别 为扮，&，......，&，Μ为观测点的总个数，建立Μ行的反演矩阵：
[0033]
[0034]其中，VR为降水粒子径向速度，VA为降水粒子径向速度的水平分量，VE为降水粒子 径向速度的垂直分量，水平分量进一步可分解为平行于正北方向和垂直于正北方向的两个 分量VAN和ν ΑΕ，β〇为水平风向与正北方向的夹角，α为天线仰角，根据关系式XV = F，即可求出 三维风场速度V。
[0035] 步骤6.2:建立反演方程：
[0036]
[0037]
[0038] 其中，VR(&)为方位角为氏的实测径向速度，VR(&)'为经最小二乘法拟合后得到的 在氏处的径向速度；
[0039] 求解方程(1)后，即可解得Van、Vae及VE此三个方向上的风速分量。
[0040]采用上述进一步方案的有益效果是充分挖掘了预处理后得到的所有径向速度中 蕴含的风场信息，且考虑所有单个径向速度数据在预期风场假设下的综合误差最小，使得 所反演得到的风场精度较高。
[0041 ]本发明解决上述技术问题的另一技术方案如下：
[0042] 一种毫米波云雷达三维风场反演系统，包括扫描控制模块、信号处理模块、信息采 集模块、主控模块和数据反演模块，其中，所述数据反演模块包括预处理模块和数据求解模 块；
[0043]所述扫描控制模块，用于控制毫米波云雷达按预设俯仰角的PPI扫描模式进行气 象目标探测，并将产生的目标探测数据传输给信号处理模块；
[0044]所述信号处理模块，用于对所述目标探测数据进行信号处理并输出原始径向速度 数据；
[0045] 所述信息采集模块，用于采集所述信号处理模块输出的原始径向速度数据，并将 其打包处理之后传输给主控模块；
[0046] 所述主控模块，用于将打包处理之后的原始径向速度数据进行保存；
[0047] 所述预处理模块，用于读取保存在所述主控模块中的原始径向速度数据并对其进 行用以消除毫米波雷达回波信号的不稳定性的预处理，得到可供三维风场反演的径向速度 数据；
[0048]所述数据求解模块，用于读取所述可供三维风场反演的径向速度数据，并在预设 求解准则下构建反演方程以对对应的所述目标探测数据进行三维风场反演，得到三维风场 速度数据。
[0049] 本发明的有益效果是:利用毫米波云雷达进行了三维风场的反演，充分利用了毫 米波段雷达对弱目标粒子探测能力强的优点，增强了测风数据的可靠性和有效性，丰富了 毫米波测云雷达的数据产品类型。另外，从数据预处理和求解准则两个方面建立反演方程， 从而保证了利用原始径向数据反演三维风场的有效性与准确性。
[0050] 进一步，所述预处理模块包括：
[0051] 排序单元，用于对所述原始径向速度数据按方位角的大小进行升序排列；
[0052] 删除单元，用于设置一个径向速度数据的预设阈值以去除所述原始径向速度数据 中的奇异点；
[0053] 滤波单元，用于使用5X1的滤波窗口按方位角排序并对去除了奇异点的原始径向 速度数据进行平滑滤波，得到相关参数数据。
[0054]采用上述进一步方案的有益效果是通过预处理将原始径向速度数据中的有用信 息提取出来，增强其可用性，保证了后续的