四川地区精细化降水预报融合订正试验及检验.pdf

暴雨灾害TORRENTIAL RAIN AND DISASTERSVol.43 No.2Apr.2024第43卷 第2期2024年4月Experiment and verification of fine gridded precipitationforecast fusion correction in SichuanZHANG Wulong1,2,CHEN Chaoping1,2,YANG Kangquan1,2(1.Sichuan Provincial Meteorology Observatory,Chengdu 610072;2.Heavy Rain and Drought-Flood Disasterin Plateau and Basin Key Laboratory of Sichuan Province,Chengdu 610072)Abstract:Fine-scale quantitative precipitation forecast is a key issue and challenge in weather forecasting services.In this study,based onhourly precipitation from the 1 km 1 km resolution Southwest China WRF-based Intelligent Numeric Grid forecast System(SWC-WINGS),a fusion-corrected experiment was conducted using time lag and probability matching methods.The fusion-corrected forecast of hourly precipitation was then verified utilizing the CMA Multi-source Precipitation Analysis System(CMAPS)three-source merged precipitation observation grid data from 1 July to 31 August 2022 in Sichuan.Finally,the fusion-corrected method was applied to a short-term heavy precipitation process over the western Sichuan Basin.The results show that:(1)Compared with the model precipitation forecasts,the time-lagged ensemble forecast was over-optimistic for small-scale precipitation and over-conservative for large-scale precipitation.(2)However,the fusion-corrected method by time lag and probability matching methods overcame the above difficulties and showed significant improvement inthe TS score,particularly in 12 h nowcast time.The TS score for hourly precipitation exceeding 0.1 mm,5 mm,10 mm,and 20 mm were increased on average by 7.2%,17.2%,28.3%,and 36.3%,respectively.(3)A case studies also showed that the fusion-corrected method hadgood improvement and correction capabilities on the hourly precipitation forecast,especially for large-scale precipitation forecasts.Key words:SWC-WINGS model;probability matching;time lag;fusion correction张武龙,陈朝平,杨康权.2024.四川地区精细化降水预报融合订正试验及检验J.暴雨灾害,43(2):176-184.ZHANG Wulong,CHEN Cha-oping,YANG Kangquan.2024.Experiment and verification of fine gridded precipitation forecast fusion correction in Sichuan J.TorrentialRain and Disasters,43(2):176-184(in Chinese).doi:10.12406/byzh.2023-077四川地区精细化降水预报融合订正试验及检验张武龙1,2，陈朝平1,2，杨康权1,2(1.四川省气象台，成都 610072；2.高原与盆地暴雨旱涝灾害四川省重点实验室，成都 610072)摘要：精细化定量降水网格预报是天气预报业务的重点和难点，基于西南区域智能数值网格模式预报系统(SouthwestChina WRF-based Intelligent Numeric Grid forecast System，SWC-WINGS)1 km 1 km分辨率的小时降水预报，利用时间滞后和概率匹配方法开展融合订正试验，再利用中国气象局多源融合降水系统(CMA Multi-source Precipitation AnalysisSystem，CMPAS)三源融合降水实况格点数据，对2022年78月四川地区的小时降水预报融合订正结果进行检验，并在四川盆地西部一次短时强降水天气过程中进行应用，结果表明：(1)时间滞后集合降水预报相较于模式降水预报，存在小量级预报过度，大量级预报过于保守的问题；(2)时间滞后结合概率匹配的降水预报融合订正方法有效提升了各量级降水预报的TS评分，尤其12 h预报时效提升显著，小时雨量超过0.1 mm、5 mm、10 mm和20 mm量级的TS评分平均提升率分别为7.2%、17.2%、28.3%和36.3%；(3)一次短时强降水天气过程的应用结果表明，时间滞后结合概率匹配的融合订正方法对模式小时降水预报有较好的改进效果，尤其对大量级降水预报有较强的订正能力。关键词：SWC-WINGS模式；概率匹配；时间滞后；融合订正中图法分类号：P456文献标志码：ADOI：10.12406/byzh.2023-077收稿日期：2023-04-25；定稿日期：2023-10-30资助项目：国家重点研发计划项目(2021YFC3000900)；四川智能网格预报创新团队项目(SCQXCXTD-202201)；中国气象局气象能力提升联合研究项目重点专项（22NLTSZ006）；四川省重点实验室项目(SCQXKJZD202101,SCQXKJYJXMS202112)；四川省重点研发项目(2022YFS0542,2022YFS0540)；西南区域创新团队项目(XNQYCXTD-202202)；中国气象局创新发展专项(CXFZ2023J016,CXFZ2024J013)第一作者：张武龙，主要从事灾害性天气和智能网格预报研究。E-mail:zwl_通信作者：陈朝平，主要从事智能网格预报研究。E-maill: Editorial Office of Torrential Rain and Disasters.OA under CC BY-NC-ND 4.0第2期张武龙，等：四川地区精细化降水预报融合订正试验及检验引 言随着社会经济高速发展，气象理论和科技不断进步，天气业务进入精细化发展阶段(金荣花等，2019)。降水是多尺度天气系统相互作用的产物，易受地形影响，且不同区域降水机理复杂，预报难度非常大(孙继松等，2012；马玉坤等，2022；韦惠红等，2022)。因此，高时空分辨率的定量降水预报是无缝隙精细化网格天气预报业务发展的重点，也是难点(代刊等，2016；金荣花等，2019；曹勇等，2021)。随着物理过程、数据同化、动力框架等方面完善和提高，数值预报技术发展迅速，特别是高时空分辨率区域模式预报性能的提升使得快速滚动更新的精细化定量降水预报成为可能(衡志炜和程晓龙，2018；黄丽萍等，2022)。然而，模式自身的系统性误差，初值场的不确定性，使得模式降水预报不可避免存在偏差。不同预报时效产生的高频次预报结果也影响着模式降水产品的实际应用效果。因此，需要一些模式后处理技术，以减小模式预报偏差，提升模式应用水平(代刊等，2018)。不同于气温、气压等气象要素，降水为非连续变量，呈非正态分布，采用的模式后处理技术也更为复杂，包括时间滞后(Hoffman and Kalnay，1983)、频率匹配(周迪等，2015)、最优百分位(张华龙等，2022)、概率匹配(Ebert，2001)、邻域法(Theis et al.，2005)等偏差订正集成技术，其中时间滞后和概率匹配是两种较为成熟的偏差订正技术。傅娜等(2013)基于华东区域模式开展时间滞后集合试验，发现小到大雨量级集合平均优于集合成员，而暴雨量级反之。唐文苑和郑永光(2019)基于中国气象局全球/区域同化与预报系统的中尺度快速分析数值预报子系统(Global/Regional Assimilation and Prediction System-Rapid Analysis andForecast System，GRAPES-RAFS)，利用时间滞后集合技术，得到订正后的晴雨TS评分显著高于模式最新时次预报的结果。张武龙等(2021)利用西南区域中尺度数 值 预 报 系 统(South West Center-WRF ADAS Real-time Modeling System，SWCWARMS)不同时效的降水预报，通过时间滞后构建集合成员，结果表明集合平均对于晴雨TS评分提升明显，而对于大量级降水订正能力有限。马申佳等(2018)在开展对流尺度集合预报试验时发现，传统集合平均对小量级降水预报效果较好，概率匹配对极端降水预报有明显优势，可改进集合平均对大量级降水预报能力偏弱的问题。庞玥等(2022)利用欧洲中期天气预报中心(European Centrefor Medium-Range Weather Forecasts，ECWMF)集合预报产品对重庆暴雨预报进行检验分析，结果表明概率匹配对暴雨强度和落区预报均有较好的订正效果。刘扬等(2022)针对“217”河南暴雨开展灾害风险评估，研究发现概率匹配可有效综合各集合成员的预估结果，优于集合平均的简单集成。不难发现，前人的工作多是基于时间滞后或概率匹配的方法得到模式降水订正结果，而对于结合二者开展融合订正试验方面的研究尚不多见。西南区域智能数值网格模式预报系统(SouthwestChina WRF-based Intelligent Numeric Grid forecast System，SWC-WINGS)是中国气象局成都高原气象研究所最新研发的31 km分辨率快速更新循环同化预报模式，本文基于SWC-WINGS模式降水产品开展融合订正试验，利用时间滞后构建集合成员，并得到集合平均预报结果；再利用概率匹配对集合平均场进行优化，提升大量级降水预报订正能力，得到逐小时滚动更新的1 km网格降水预报融合订正结果，并通过检验分析，考察时间滞后结合概率匹配融合订正方法对模式预报的改进效果，以期提升快速更新模式降水预报能力。1 SWC-WINGS模式介绍SWC-WINGS使用美国国家环境预报中心全球预报 系 统(National Centers for Environmental PredictionGlobal Forecast System，NCEP GFS)全球分析预报资料作为模式背景场和边界条件，同化地面站、探空、飞机报、船舶、浮标、雷达反射率和卫星资料等。1 km模式垂直方向为 51 层(模式层顶 10 hPa)，格点数 11801030，每日逐小时滚动预报，预报时效30 h，可提供覆盖四川地区(97.25108.65E，2634.45N)的高时空分辨率的气象要素产品。实际业务中模式产品可在滞后4 h以内提供至服务器，供业务使用。2 资料与方法2.1 资料来源本文所用资料包括模式预报资料和实况资料，资料时段为 2022 年 7-8 月，资料范围为四川地区(97.25108.65E，2634.45N)，同模式预报范围一致。模式预报资料为SWC-WINGS模式每日逐小时滚动预报的小时降水产品，预报时效30 h，空间分辨率1 km 1 km，时间分辨率1 h。实况资料选取基于中国气象局多源融合降水系统(CMA Multi-source Precipitation Analysis System，CMPAS)的三源融合降水实况格点数据(师春香等，2019)，空间分辨率1 km，时间分辨率1 h，与模式资料一致。CMPAS实况降水融合产品的适用性评估工作(许冠宇等，2020；龙柯吉等，2020)较多，结果表明1 km产品优于5 km产品，三源产品优177第43卷暴雨灾害于二源产品，CMPAS的1 km三源融合小时降水实况格点数据与站点数据较为一致，具有很高的可信度(龙柯吉等，2020)。2.2 研究方法2.2.1 融合订正方法(1)时间滞后集合方法。时间滞后集合方法是利用不同起报时刻对某一时刻预报结果进行集合的方法(Hoffman and Kalnay，1983)。图 1 给出本文时间滞后集合预报生成示意图。对于H时刻的预报而言，假设当前能获得的模式最新起报时刻为h，则从h时刻开始搜索H时刻预报的集合成员，由于SWC-WINGS模式预报时效为30 h，则理论上最多可搜索到h-30h，那么最大集合成员数为30。集合成员数N的具体计算方式如下N=30-(H-h)+1(1)显然，当H-h=1时，N取最大值为30。然而，实际业务中考虑资料的传输运算时间所导致的滞后效应，最大集合成员数必然小于30，随着预报时刻H增加，集合成员数N逐渐减小。再将H时刻所有集合成员进行集合平均，则可得到H时刻的时间滞后集合预报。计算公式(2)中Tle为时间滞后集合预报(Time-Lagged Ensemble，TLE)，Ri为第i个集合成员的预报。Tle=i=1NRiN（2）hh-1h-2h-3h-30130 h时效H时刻时间滞后集合预报H图1 时间滞后集合示意图Fig.1 Time-lagged ensemble configuration(2)概率匹配订正法。概率匹配通常用于混合两种具有不同时空分布的数据类型，其中一种类型具有较好的空间表示，而另一种类型则可提供更高的精度信息(Ebert，2001)。研究表明，集合平均可使降水落区变得更加准确，但也会使集合成员降水预报强度被低估，损失掉量级预报精度，降低强降水的预报能力(马申佳等，2018；代刊等，2016)。因此，集合平均可体现出更好的空间分布，而集合成员则具有更高的量级精度信息。概率匹配方法原理演示如图2所示，假设通过时间滞后构建了三个集合成员Member1、Member2和Member3，各成员包含有33的网格降水预报。首先，对三个成员进行集合平均，并按照平均值大小进行排序，网格中最大值为51，标记为“Rank 1”，依次类推。其次，对三个成员的所有网格预报值从大到小进行排序，每三个一组，“Rank 1”、“Rank 2”、“Rank 3”为第1组(“Team 1”)，依次类推。最后，取每组的中位数，赋给集合平均中相同序号的网格点，如“Team 1”的中位数45，赋给集合平均中标记为“Rank 1”的网格点，依次类推。最终，通过概率匹配对时间滞后集合平均场的优化，得到时间滞后结合概率匹配的融合订正结果(Fusion Correction Ensemble，FCE)。图2 概率匹配示意图Fig.2 Probability matching configuration2.2.2 检验方法检验评估方法包括计算准确率PC(记为Pc)、空报率FAR(记为Rfa)、漏报率PO(记为Po)、TS评分(记为Ts)四种方法，具体计算公式如下Pc=Na+NdNa+Nb+Nc+NdRfa=NbNa+NbPo=NcNa+NcTs=NaNa+Nb+Nc(3)式(3)中Na为某一量级降水预报与实况均发生的格点数，Nb为某一量级降水预报发生而实况未发生的格点178第2期张武龙，等：四川地区精细化降水预报融合订正试验及检验数，Nc为某一量级降水预报未发生而实况发生的格点数，Nd为某一量级降水预报与实况均未发生的格点数。3 试验方案本文按照实际业务需求开展试验，需在每个整点前0.5 h形成新的124 h时效的1 km网格降水预报，运算时间15 m左右，试验方案如图3所示。首先，基于SWC-WINGS模式0.01 0.01分辨率的小时降水预报，利用时间滞后构建集合成员，再通过集合平均得到时间滞后集合预报(TLE)。然后，利用概率匹配对时间滞后集合平均降雨率进行重新分配，得到时间滞后结合概率匹配的融合订正结果(FCE)。以08 00(北京时，下同)起报为例，07 15前能获取的SWC-WINGS模式最新起报时次为 03 00，那么 08 00 起报的 1 h(2 h，3 h，.，24 h)时效对应的模式最新时次预报为模式03 00起报6 h(7 h，8 h，.，29 h)时效，集合成员数N为25(24，23，.，2)。将08 00起报的1 h时效对应的25个集合成员进行集合平均，则可得到该时效的时间滞后集合预报。然后，利用概率匹配的方法，整合该时效的25个集合成员网格降水预报信息并进行排序，再与时间滞后集合平均网格排序结果一一匹配，最终得到该时效的时间滞后结合概率匹配的融合订正结果，其他时效依此类推。研究表明，集合成员数的多少直接影响集合订正效果，集合成员数越多，预报准确率越高(唐文苑和郑永光，2019)。此外，考虑每隔1 h滚动更新，得到新的124 h集合订正场，更关注临近时效的预报结果。因此，本文主要针对2022年78月逐小时起报的112 h时效的TLE、FCE和相应的SWC-WINGS模式小时降水预报，利用逐小时1 km实况格点数据，采用“点对点”的方式进行晴雨准确率、空报率、漏报率、TS评分检验对比分析，在集合成员充足的情况下考察融合订正的提升效果。4 结果分析与检验研究表明，针对降水的模式后处理技术，可能对晴雨预报质量提升较大，但对大量级降水订正能力有限(马申佳等，2018；张武龙等，2021)。因此，本文分为晴雨检验和分级检验两个部分，对时间滞后结合概率匹配的融合订正结果进行分析。4.1 晴雨检验为了考察时间滞后和概率匹配对模式晴雨预报的改进效果，对小时降水超过0.1 mm量级进行准确率PC、空报率FAR、漏报率PO和TS评分检验。图4给出了2022年7月1日8月31日逐小时起报112 h时效的TLE、FCE和相应的SWC-WINGS模式小时降水预报的晴雨(小时降水超过0.1 mm)准确率PC，从图4中可见，模式预报在各时效晴雨准确率较为一致，112 h时效的平均预报准确率为86.6%。TLE的晴雨准确率在各时效明显低于模式预报，平均预报准确率为78.4%。进一步使用了概率匹配优化之后，各时效晴雨准确率较TLE有明显提升，且分布较为一致。FCE晴雨准确率在各时效均高于模式预报，平均预报准确率为87.2%。预报时效/h0.900.860.820.780.740.70准确率123456789101112模式TLEFCE图4 2022年7月1日8月31日逐小时起报112 h时效的TLE、FCE和相应的SWC-WINGS模式最新时次小时降水预报的晴雨准确率Fig.4 Sunny or rainy accuracy for hourly precipitation forecasts derived fromthe 1-h runs,112 h period validity TLE,FCE and the up-to-dateSWC-WINGS model from 1 July to 31 August in 2022图5给出了2022年7月1日8月31日逐小时起报112 h时效的TLE、FCE和相应的SWC-WINGS模式小时降水预报的晴雨空报率FAR、漏报率PO和TS评分以及FCE相较于模式TS评分的提高率。从图5a可以看到，模式降水预报的空报率随着预报时效的延长逐渐增加，平均空报率为0.669。FCE的空报率随时效的变化趋势与模式一致，但其空报率明显低于模式，平均空报率为0.646。TLE的空报率则明显高于模式，随着预报时效变化趋势不明显，平均空报率为时间滞后与概率匹配融合(FCE)检验分析做检验概率匹配时间滞后集合(TLE)集合平均Member1Member2Member3Member4Member构建集合成员SWC-WINGS实况格点图3 时间滞后结合概率匹配的融合订正试验方案Fig.3 Experiment scheme of the fusion correction by time lag andprobability matching methods179第43卷暴雨灾害0.732。从三者晴雨漏报率的对比来看(图5b)，模式预报的漏报率随着预报时效的延长逐渐增加，平均漏报率为0.624。FCE的漏报率随时效的变化趋势与模式一致，但其漏报率低于模式，平均漏报率为0.608。TLE 的漏报率也明显低于模式，平均漏报率为0.286。由此可见，TLE相较于模式预报，虽然减小了漏报率，但增加了空报率，从而使其晴雨准确率反而不如模式预报。FCE相较于模式预报，减小了漏报率的同时，也减小了空报率，从而使其晴雨准确率高于模式预报。0.750.720.690.660.630.66空报率模式TLEFCE123456789101112预报时效/h(a)0.70.60.50.40.30.2123456789101112预报时效/h模式TLEFCE(b)漏报率0.260.240.220.200.180.16TS评分预报时效/h123456789101112模式TLEFCETS提高率(c)图5 2022年7月1日8月31日逐小时起报112 h时效的TLE、FCE和相应SWC-WINGS模式最新时次小时降水预报的晴雨空报率(a)、漏报率(b)和TS评分以及FCE相较于模式TS评分的提高率(c)Fig.5 Sunny/rainy(a)false alarm rate,(b)missing alarm rate,(c)TS score and theincreased percentage of TS score between FCE and model for hourly precipitationforecasts derived from the 1-h runs,112 h period validity TLE,FCE and theup-to-date SWC-WINGS model from 1 July to 31 August in 2022从图5c中可见，模式预报随着预报时效的延长，TS评分逐渐降低，平均TS评分为0.213。三者对比来看，TLE的TS评分最高，且TS评分随着预报时效的延长变化不明显。根据时间滞后集合的计算原理，当某一个集合成员预报有雨的情况下，集合平均结果则必然预报有雨，从而使得空报率增加，漏报率减小，导致过度预报的问题。因此，TLE的TS评分较模式预报提高明显，是在增加空报率的基础上实现的，它的综合预报效果是不好的。FCE的TS评分随时效的变化趋势与模式一致，并且其TS评分在各时效均高于模式，平均TS评分为0.228。FCE在同时降低空报率和漏报率的情况下，有效的提升了模式晴雨TS评分，12 h时效的平均提高率为7.2%，112 h时效的平均提高率为7.0%。4.2 分级检验为了考察时间滞后和概率匹配对模式各量级小时降水预报的改进效果，分别对小时降水超过5 mm、10 mm、20 mm量级进行空报率FAR、漏报率PO和TS评分检验。图6给出2022年7月1日8月31日逐小预报时效/h123456789101112TS提高率/%0.070.060.050.040.030.020.010TS评分201816141210模式TLEFCETS提高率(c)0.960.950.940.930.920.910.90空报率123456789101112预报时效/h模式TLEFCE(a)0.920.900.880.860.840.820.80漏报率123456789101112预报时效/h模式TLEFCE(b)图6 同图5，但为小时降水超过5 mmFig.6 Same as Fig.5,but for hourly precipitation exceeding 5 mm180第2期张武龙，等：四川地区精细化降水预报融合订正试验及检验时 起 报 112 h 时 效 的 TLE、FCE 和 相 应 的SWC-WINGS模式预报的5 mm以上量级小时降水空报率、漏报率、TS评分以及FCE相较于模式TS评分的提高率。图 6a 表明，各预报时效 TLE 空报率小于FCE，也小于模式预报，且三者均表现为随预报时效增加而增加的趋势。模式预报、TLE和FCE的112 h时效平均空报率分别为0.947、0.919和0.937。从三者漏报率对比可以看到(图6b)，FCE漏报率小于模式预报，小于TLE，整体上也表现为随预报时效增加而增加的趋势。模式预报、TLE和FCE的112 h时效平均漏报率依次为0.850、0.878和0.837。时间滞后集合方法可以降低5 mm以上量级降水的空报率，但却增加了其漏报率，这是由于时间滞后集合对于大量级降水的平滑作用，导致其对于大量级降水预报过于保守。进一步使用概率匹配优化之后，可以消除时间滞后集合方法造成的平滑作用，突出单个集合成员的大量级降水预报，使得漏报率明显降低(图6b)。分析图6c可知，三者TS评分均表现为随预报时效增加而减小的趋势。TLE的112 h时效平均TS评分为0.051，较模式预报提升明显，但是是在增加漏报率的基础上实现的，存在预报过于保守的问题。FCE的112 h时效平均TS评分为0.048，虽略低于时间滞后集合，但高于模式预报，并且是在同时降低空报率和漏报率的情况下实现的，预报效果更好。FCE在12 h时效的5 mm以上量级降水的TS评分平均提高率为17.2%，112 h时效平均提高率为16.1%。2022年7月1日8月31日逐小时起报112 h时效TLE、FCE和相应的SWC-WINGS模式预报的小时降水超过10 mm量级的空报率、漏报率和TS评分对比与5 mm量级类似(图略)。TLE在10 mm以上量级降水的空报率小于模式起报，漏报率大于模式起报，依然存在预报过于保守的问题。经过概率匹配优化之后，空报率和漏报率同时减小，TS评分也有较明显提升。FCE的112 h时效10 mm以上量级降水平均TS评分为 0.025，12 h 时效的 TS 评分平均提高率为28.3%，112 h时效平均提高率为17.4%。三者的20 mm以上量级降水空报率、漏报率对比检验结果与5 mm和10 mm量级类似(图7a，b)，不再赘述。然而，三者20 mm以上量级降水的TS评分结果却与5 mm和10 mm量级有明显差异。从图7c中可见，TLE 的 TS 评分要明显低于模式预报，更低于FCE。在5 mm和10 mm以上量级降水的TS评分检验中，TLE在增加漏报率的情况下，TS评分是高于其余二者的。这说明相较于 5 mm 和 10 mm 以上量级降水，TLE对于20 mm以上量级降水的预报更加保守。FCE的TS评分在大多数时效内优于模式预报，特别是在临近时效有很好的改进效果。FCE的112 h时效的20 mm以上量级降水的平均TS评分为0.01，在12 h时效平均提高率为36.3%，112 h时效平均提高率为12.2%。图7 同图5，但为小时降水超过20 mmFig.7 Same as Fig.5,but for hourly precipitation exceeding 20 mm1.000.990.980.970.960.950.940.930.920.910.90空报率123456789101112预报时效/h1.000.990.980.970.960.950.940.930.920.910.90漏报率123456789101112预报时效/h模式TLEFCE模式TLEFCE(a)(b)0.0140.0120.0100.0080.0060.0040.0020TS评分6040200-20-40-60TS提高率/%123456789101112预报时效/h模式TLEFCETS提高率(c)由此可见，时间滞后集合预报的晴雨(小时降水超过0.1 mm)漏报率较SWC-WINGS模式预报有所减小，但空报率有所增加，存在预报过度的问题。对于5 mm、10 mm和20 mm以上量级降水空报率较模式有所减小，但漏报率有所增加，存在预报过于保守的问题。经概率匹配优化之后，在同时降低空报率和漏报率的情况下，时间滞后结合概率匹配的融合订正预报有效的提升了各量级的TS评分，尤其在12h时效提升显著，小时降水超过0.1 mm、5 mm、10 mm和20 mm量级的TS评分平均提升率分别为7.2%、17.2%、28.3%和36.3%。181第43卷暴雨灾害99102105108E34N3230282699102105108E99102105108E99102105108E34N3230282650201052.50.1mm(a)(b)(c)(e)(f)(h)(d)(g)5 降水预报个例订正检验通过以上分析得知，经时间滞后结合概率匹配融合订正后，预报质量较模式预报有明显提升。但是上述结果是基于 2022 年 78 月长时间的批量检验评分，对降水预报的订正效果还需要一些典型的强降水个例进一步检验。本文选取四川盆地西部一次强降水个例，采用融合订正技术改进模式降水预报并检验分析。2022年7月12日20 0013日08 00，四川盆地西部出现了一次暴雨到大暴雨过程，并伴有明显的短时强降水。雷达组合反射率拼图显示(图略)，7月12日23 0013日00 00，盆地西部有对流初生，到13日01 0002 00，逐渐发展为线状对流，选取这2个代表性时段的实况和预报降水进行评估分析。根据业务试验流程，7月13日00 00的降水实况对应的是7月12日23：00起报TLE和FCE的第1 h时效预报，模式最新时次预报则是18 00起报模式的第6 h时效预报。7月13日00 00的降水实况显示，在盆地西部沿山出现了南北两个强降水中心(图8a)；模式预报在整个盆地西部沿山均出现了分散的20 mm以上强降水(图8b)；TLE预报能够将南北两个降水中心突显出来，但量级偏小，且川西高原弱降水有明显空报(图8c)；经图8 2022年7月13日00 00(a)和02 00(e)1 h降水实况、SWC-WINGS模式7月12日18 00(b)和7月12日20 00(f)起报第6 h时效雨量预报、TLE 7月12日23 00(c)和7月13日01 00(g)起报第1 h时效雨量预报和FCE 7月12日23 00(d)和7月13日01 00(h)起报第1 h时效雨量预报Fig.8 Precipitation observation at(a)00 00 BT and(e)02 00 BT on 13 July,the 6-h forecast from SWC-WINGS model initialized at(b)18 00 BT and(f)20 00 BT on 12 July,the 1-h forecast from TLE initialized at(d)23 00 BT on 12 July and(h)01 00 BT on 13 July,the 1-h forecast from FCE initialized at(d)23 00 BT on 12 July and(h)01 00 BT on 13 July,2022图9 TLE、FCE 2022年7月12日20 0013日08 00逐小时起报第1 h时效和SWC-WINGS模式最新时次预报的空报率(a)、漏报率(b)及TS评分和FCE相较于模式的提高率(c)Fig.9 The(a)false alarm rate,(b)missing alarm rate,(c)TS score and theincreased percentage of TS score between FCE and model for the 1-hforecasts from TLE,FCE and the up-to-date SWC-WINGS model from20 00 BT on 12 July to 08 00 BT on 13 July,20221.00.80.60.40.20空报率0.151020降水量级/mm1.00.80.60.40.20漏报率0.151020降水量级/mm模式TLEFCE模式TLEFCE(b)(a)TS提高率/%0.50.40.30.20.10TS评分0.151020降水量级/mm6050403020100模式TLEFCETS提高率(c)182第2期张武龙，等：四川地区精细化降水预报融合订正试验及检验概率匹配优化之后，南北两个降水中心的预报雨量明显增大，与实况较为一致(图8d)。至7月13日02 00，南北两个实况强降水中心逐渐连接形成线状强降水对流雨带(图8e)；模式预报在盆地西部沿山的强降水落区与实况接近，但空间分布较为离散，且对流雨带右侧有空报的分散强降水(图8f)；TLE的强对流雨带量级偏弱，并存在大片弱降水的空报(图8g)；经概率匹配优化之后，强对流雨带更加清楚，降水量级增加，并且对于弱降水进行了明显的消空处理(图8h)。该过程中TLE、FCE的第1 h时效与模式预报的空报率、漏报率和TS评分如图9所示。可见在0.1 mm以上量级降水检验中，TLE 在增加空报率的情况下，TS 评分优于FCE，更优于模式预报。在较大量级检验中，FCE在同时降低空报率和漏报率的情况下，TS 评分均高于TLE，更高于模式预报，5 mm、10 mm和20 mm以上量级降水的 TS 评分提高率分别可达 45.7%、43.7%和34.0%。由此可见，个例订正试验与批量订正检验结果基本一致。无论从空间对比，还是TS评分检验，时间滞后结合概率匹配的融合订正方法对模式预报均有明显改善。6 结论与讨论本文基于SWC-WINGS模式1 km 1 km分辨率的小时降水产品，采用时间滞后和概率匹配方法在四川地区开展融合订正试验，得到逐小时滚动更新的1 km网格降水预报，再利用CMPAS三源融合降水实况格点数据，对2022年78月四川地区的小时降水融合订正试验结果进行检验，并应用于四川盆地西部一次短时强降水过程的模式预报订正，得到以下主要结论：(1)时间滞后集合预报的晴雨(小时降水超过0.1 mm)漏报率较SWC-WINGS模式预报有所减小，但空报率有所增加，存在预报过度的问题。对于5 mm、10 mm和20 mm以上量级降水空报率较模式有所减小，但漏报率有所增加，存在预报过于保守的问题。