雷达降水估测中变分方法的改进研究.docxVIP

雷达降水估测中变分方法的改进研究?慕熙昱 1 雷达降水估测中变分方法的改进研究? 慕熙昱 1 徐琪 2 刘国庆 3 1.江苏省气象科学研究所 210008 南京 2.民航华东空中交通管理局江苏分局 210000 南京 3.南京工业大学 210000 南京 摘要：雷达作为定量测量降水的工具，在探测降水分布和结构上比较准确，但是在精度上还是低于雨 量计。本文利用变分法，综合多普勒天气雷达数据及雨量计数据进行定量降水估测，其本质是通过调整模 式解以达到模式解与观测数据之间最小二乘条件。最开始在利用变分方法校正雷达降雨场的过程中，对变 分方程中的系数是经验选取的，但是这种经验关系有很大的局限性。在对变分方程求解过程中还有一个问 题就是边界的处理。在对变分方程求解的过程中，采用交叉验证法。交叉验证( cross－validation) 理论 是：对某一点来说,用其余若干个观察值对该点进行估计,每一点都重复这个过程, 用估计值与真实值进行 对比得误差参量,最后进行处理,使这些误差参量达到要求。该方法基于定理：椭圆型偏微分方程中，区域 内的解连续依赖于边界。首先利用变分方程可以求出雨量与边界值关系的表达式；然后根据已有的雨量计 数据，利用最优化方法可以求得边界值。将边界值代入变分方程，则方程可解，就求得每个离散点上的雨 量。这种方法的优点在于不需要对边界条件进行假定，改进了以往运用变分法进行雷达定量降水估测时存 在的问题。通过对交叉验证方法对变分方程求解，以及将求解后的雨量进行拟合，得到新的 Z ? R 关系。利 用实测雨量计数据对新的 Z ? R 关系进行检验。结果表明交叉验证方法能够在不假定边界条件的情况下， 对变分方程求解，解的状况比较符合实际情况。对层状云降水来说，经典的 Z ? R 关系能够很好的反映降 水的情况，拟合后的结果与直接应用经典 Z ? R 关系得到的雨量没有较大差异，两者与实际情况的符合程 度都很高；对于积层混合云来说，经过拟合后的降水结果明显优于直接用 Z ? R 关系计算的结果。 关键词：变分方法，降水估测，交叉验证 一、引言 雷达作为定量测量降水的工具，在探测降水分布和结构上比较准确，但是在精度上还是 低于雨量计。尽管双偏振雷达的出现使雷达估测降水精度大大增加，但是由于国内目前双偏 振雷达研究还处于刚起步阶段，并且目前国内布网雷达基本上不具有双偏振功能，因此，如 何有效地利用雨量计校正雷达估测降水来提高雷达测雨精度是非常必要的。为了用雨量计校 正雷达降雨，国内外科研工作者提出了很多种方法，如平均法（Wilson，1979），变分方法 （Ninomiya and Akiyama，1978；Sun Shouxiang et al.,1993;赵坤等,2001）、卡尔曼滤波方法 （Ahnert，1986；DakSyangLin and Krajewski，1989）。 Daley 曾明确指出伴随变分和 Kalman 滤波是大气数据同化的未来发展方向。如今他的 预见已被越来越多的工作所验证。变分方法通过非线性模式解与不同时次观测资料集的全局 调整达到同化的目的。 二、研究方法的提出 在雷达定量估测降水领域,近年来的主要发展在变分分析方案，其本质是通过调整模式 解以达到模式解与观测数据之间最小二乘条件。在对变分方程求解过程中的一个问题就是边 界的处理。已往的研究中对于边界的处理方法是假定边界。这种方法从理论上来说存在误差。 为解决这个问题，本文尝试用交叉验证法对变分方程求解，以改善雷达定量测量降水结果。 三、详细算法 首先考虑任一矩形区域 Ω ，有边界 ? ，在 Ω 区域内有 Ng 个雨量站，即为集合? 。在 Ω ? 资助项目：江苏省自然科学基金 BK2010599。本文部分内容投稿至其它学术期刊。 1 区域内建立如下变分：?????? ? ?22R ? ? ? ?R?2 ? ? ?R ? R dxdy?2? ? ? ? R ? RJ ? ? ??? ?Ω ????? ?x ? ? ?y ?rg????? ?????其中 Rr 表示由经典 Z－R 关系估测出来的降雨量， Rg 是地面面雨量站的实测降雨量。第一项约束表示 R 在 Ω 区 区域内建立如下变分： ?? ?? ?? ? ?2 2 R ? ? ? ? R ?2 ? ? ?R ? R dxdy ?2 ? ? ? ? R ? R J ? ? ? ?? ? Ω ?? ? ? ? ?x ? ? ?y ? r g ??? ? ? ? ?? ?? 其中 Rr 表示由经典 Z－R 关系估测出来的降雨量， Rg 是地面面雨量站的实测降雨量。 第一项约束表示 R 在 Ω 区域内的光滑性，第二项约束表示 R 在整个区域内与