利用动态ZI关系雷达回波定量估测降水.doc

钤伟妙 联系电话基于动态Z-I关系雷达回波定量估测降水方法研究* 陈静1 钤伟妙1 韩军彩1 连志鸾2 1石家庄市气象局 石家庄 050081 2 河北省气象局 石家庄 050021 摘要：选取石家庄市新一代天气雷达2009—2012年29个降水过程雷达二次产品数据（低层5个仰角及组合反射率因子）和加密自动站逐分钟降水资料，依据与雷达回波强度最相近及与自动站位置距离最近的原则，建立随时空变化的动态Z-I关系，来定量估测无自动站观测格点的降水量，从而实时获取高分辨率雷达估测降水资料；对雷达定量估测降水与自动站观测降水的空间分布进行了对比，利用交叉检验方法分析了雷达定量估测降水的误差，结果表明：该方法反演降水的平均误差较小，为-0.6 mm·h-1，能较好的表现出降水的时空分布特征，显示出不同等级降水的面积和演变，具有反演40 mm·h-1以上的强降水中心的能力，表明该方法具有业务应用价值，对强降雨灾害的评估具有指导作用。 关键词：雷达回波，二次产品，加密自动站，逐分钟降水，动态Z—I关系，误差分析 引言 近年来，利用新一代多普勒天气雷达反演降水技术已成为雷达资料应用研究的重点，也是发展雷达技术的主要研究方向之一（李柏等，2013）。通过雷达反演降水可及时获取大面积高分辨率降水信息，弥补降水观测站点稀疏的不足。雷达反演降水是精细化气象服务的基础资料，对天气气候的监测和分析、强降雨灾害的评估、数值预报模式的检验和校准以及水文分析等具有很高的应用价值（姚燕飞等，2007；白先达等，2011；万玉发等，2013）。 雷达定量估测降水是利用雷达回波资料建立雷达回波与降水的Z—I关系，间接获得降水资料。雷达回波是对空中雨滴大小及密度的一种电磁学间接测量，由于雨滴粒径分布及雨滴下落末速度的变化性很大，不同性质的降水适用不同的Z—I关系，即使相同的回波强度所对应的降水强度也可能不同（仲凌志等，2007）。雨滴谱的自然可变性是雷达测量降水不确定性的一个重要原因（万玉发等，2013）。2006年全国自动站资料开始实时上传（任芝花等，2010），降水观测的时间和空间密度加大，加密降水资料与雷达回波时空分布有很好的对应关系（刘雨佳等，2014；Luo et al.，2013），使雷达估测降水的精度不断提高并稳定，尤其有利于提高积云降水过程降水反演的精度（东高红等，2012b）。目前建立Z—I关系的方法主要分为两种，一种是统计学方法，如采用最优化法（袁晓清等，2010；东高红等，2012a）、概率配对法（郑媛媛等，2004；王丽荣等，2008）等统计多个降水过程资料建立经验Z—I关系。此外，一些研究还根据降水天气形势、不同季节、雷达探测距离和不同降水类型（仲凌志等，2007；陈小华等，2013）等对雷达回波进行分类，建立本地化的Z—I关系。为进一步提高估测降水精度，研究者还利用平均校准法（张亚萍等，2013）、变分同化法（马慧等，2008）及卡尔曼滤波校准法（邵月红等，2009；王春青等，2010）等进行订正，这些均可在一定程度上改善估测降水的精度，但校正方法一般计算复杂，消耗时间长。另一种方法是动态分级Z—I关系法（汪瑛等，2011），该方法建立在快速更新资料的基础上，实时频繁变动Z—I关系中的A、b值，不依赖于长期历史资料，只考虑很短时间(6分钟)内降水强度和降水性质的延续性。动态分级Z—I关系方法将上一时次雷达回波反射率因子按强度分级，根据最优化方法选取最优的一组A、b值，将该最优分级Z—I关系用于下一时刻估测降水中，该方法能改进对短时强降水的估测，尤其是强降水低估严重的问题（胡胜等，2012）。 本文在以前研究成果的基础上，提出了利用日常业务应用的雷达二次产品回波反射率因子和加密自动站降水资料，确定随时间和空间变化的动态Z—I关系，通过自动站分钟降水数据和雷达回波数据建立动态Z—I关系，来定量估测无自动站观测格点的降水量。与汪瑛等（2011）不同，本文不仅考虑了Z—I关系时间上变化，还考虑了其空间上的动态变化，利用最优化方法得到某时次自动站位置上的A、b值，同时依照与雷达回波强度最接近及与自动站距离最接近的原则，将得到的A、b值延伸应用到该时次无自动站观测的雷达回波格点上，从而实时得到具有业务应用价值的高分辨率定量估测降水资料，为精细化预报服务及强降雨积涝风险评估提供支撑。 1 资料和方法 1.1 资料 本文以石家庄市为研究区域，包括石家庄市区和周边17个县（市），石家庄位于太行山东麓，地形比较复杂，西部为太行山区，东部为平原。降水资料为石家庄市191个自动站观测的逐分钟雨量数据，雷达数据为石家庄新一代天气雷达CINRAD—SA逐6分钟分辨率的雷达二次产品数据，包括组合反射率因子和低层5个仰角（0.5