2003年江淮梅雨期一次特大暴雨中尺度三维动力结构的双多普勒雷达反演探究.pdfVIP

七、新技术应用 ·399· 2003年江淮梅雨期一次特大暴雨中尺度三维 动力结构的双多普勒雷达反演研究’ 周海光1王玉彬2 (1中国气象科学研究脘，北京1000812北京市气象局，北京，100089) 摘要：2003年6月下旬至7月上旬，淮河流域出现了持续强暴雨。本文使用地基多部多普勒雷达综合和连续调 究表明，中低层辐合线和中尺度气旋是造成此次暴雨的重要原因。 关■词：特大暴雨辐合线双多普勒天气雷达 三维风场风场反演 强度之大为历史上少见。研究其形成机理，具有重要意义。 国家重大基础研究“973中国暴雨”项目2003年在安徽省组织实施了双多普勒雷达同步观测，本 文使用7月4～5日合肥和马鞍山雷达探测到的特大暴雨的双多普勒雷达资料，用MUSCAT反演技 术对其三维动力结构进行了反演研究。 双多普勒雷达联合探测可以明显地提高三维风场的反演精度，从20世纪60年代末国外就开始 了这方面的研究工作，先后发展了一系列双雷达风场反演算法，随着变分技术的引入，明显地提高了 反演的可靠性和精度，从而加深了对暴雨、飑线等灾害性天气的认识。 最近，我们也在双、三多普勒雷达三维风场反演技术，使用双、三多普勒雷达研究暴雨中尺度三维 动力结构等方面做了一些研究。 1 MUSCAT三维风场反演技术 地基MUSCAT双、三多普勒雷达三维风场反演技术采用变分法，定义泛函： F—l (1) EA(“，，叫)4-B(“，口，)4-C(“，口，训)]dxdy Js 令一阶导数为零，便可求得三维速度场，即： 堡一0，。aF一0，婴一0 (2)“’ a—ui一’瓦一0’丽一0 方程(1)中A为数据调整项，是多普勒雷达探浏到的径向速度与笛卡儿坐标系中三维速度的最小 二乘法表示，常规的反演技术是将雷达探测点的径向速度插值到反演网格点的径向速度，做为单独的 一个步骤来完成，也即先通过插值算子得到反演网格点上的径向速度，然后再通过风场反演方程组得 到网格点上的三维速度场。但由于径向速度是矢量，探测点径向速度与网格点的径向速度的方向不 一致，这样插值得到的径向速度会存在一些误差，从而影响了反演精度。函数A是通过最小二乘法拟 合调整技术在泛函中实现插值，避免了插值问题带来的误差，这种在目标函数中直接完成插值计算的 技术，在基于变分法的单、双多普勒天气雷达风场反演技术中使用得比较广泛，这种方案可以提高反 演精度。 B项为质量连续方程的最小二乘法表达式，连续方程的平方项作约束条件，使积分连续方程所造 成的误差累积有了改善，而这正是反演技术中最难解决的。 C项为二次微分约束，目标函数C的作用不仅仅是一个各向同性低通滤波器，滤去风场的不合理 的小尺度扰动，使反演得到的三维风场更加合理；正如Sasaki所论述的，低通滤波器对改善数据稀疏 *本文耨判科技都2002年鏖科研院所社会公釜研究专项蠢全项日“突发性强灾害天气璜警系统的研制”的赘助， W够／ ’400‘新世纪气象科技创新与大气科学发展——03．7淮河大水的水文气象学问题 或者数据噪音较大区域的分析结果很有益处。 数值表明·MUSCAT算法具有较高的反演精度，而且其最佳反演区比常规迭代算法的反演区要 大得多。 2双多普勒雷达同步观测 2003年6～7月，“973中国暴雨”项目外场试验在安徽省设置了一个双多普勒雷达同步探测区， 需要指出，试验中所指的同步观测并不是严格意义上的时间同步观测，参与同步观测的两部雷达以北 京时间为参考时间，在整点时刻开始各自的体积扫描进行数据采集，由于两部雷达完成一个体扫所需 的时间并不完全相同，因此有部分体扫资料采集时所处的时间段并不一致，我们选用时间段比较一致 的双多普勒雷达体扫数据用于三维风场反