雷达预估信息(0~2 h)、数值天气预报产品(0~72 h)及基于高空和地面大气探测资料的综合预报信息(0~24 h)在不同预见期的降雨预报各有差异,预报效果上各有优劣,而针对山洪的预报存在准确率低、预警有效时间短等问题,研究面向山洪灾害防治区的数值预报模式的选取、模式最优物理参数化组合方案、同化方案、降水偏差订正技术以及基于大数据的自分型雷达估测降水最优化算法。以湖南临湘市为试验区进行验证,结果表明:WRF中尺度数值模式适用于山洪灾害降雨预报,其WSM6云微物理过程、Grell-Devenyi ensemble对流参数化方案和YSU边界层参数方案对致洪山洪暴雨过程模拟较好,基于频率(或面积)匹配的降水偏差订正方法能显著改善模式降水预报中雨量和雨区范围的系统性偏差,大数据分析方法应用于雷达估测降雨能显著提高准确度。在此研究基础上,提出了基于高空和地面大气探测、数值预报模式的0~72 h短期定量降雨预报和基于雷达的0~2 h临近定量降雨预报的多元信息耦合预报方法,将对致灾山洪的预报时间由2 h延长至72 h,并可显著提高了山洪预报的精度。 The precipitation forecasts from radar and satellite cloud pictures (with a lead time of 0 - 2 h), from upper air and ground surface atmospheric sounding (with a lead time of 0 - 24 h) and from numerical forecasting mode (with a lead time of 0 - 72 h) have different forecast effect in different forecast pe-riods. Considering the short lead time and low accuracy forecast for early warning of mountain tor-rents, study on the numerical forecasting model, the model’s optimized combination of parameters, data assimilation, correcting the errors on forecasted rainfall in flash-flood affected region and an op-timized calculating method based on big data for radar monitoring precipitation. Linxiang City of Hu-nan Province was selected as experimental area and the results show that the WRF numerical forecast is more suitable for the forecast of disasters including mountain torrents. Among them, WSM6 cloud microphysical process, Grell-Devenyi ensemble convective parameterization scheme and YSU boun-dary layer parameter configuration scheme can greatly simulate the process of mountain torrents and storm causing a flood, precipitation deviation correction method based on frequency (or area) match can considerably improve the systematic deviation in the rainfall and scope of rain area in the preci-pitation forecast, the big data analysis method used in radar rainfall forecast can remarkably improve the accuracy. A forecasting disastrous flash floods method has been proposed by coupling upper air and ground surface atmospheric sounding, the 0 - 72 h short-term quantitative precipitation forecasts (QPF) from numerical forecasting mode and the 0 - 2 h nowcast QPF based on radar, which helps to increase the lead time of flash flooding forecasts from 2 h to 72 h and enhance significantly the fore-casting accuracy.
熊明1,杨文发1,李俊2,周北平3,訾丽1
1长江水利委员会水文局,湖北 武汉
2中国气象局武汉暴雨研究所,湖北 武汉
3南京信息工程大学,江苏 南京
收稿日期:2017年3月20日;录用日期:2017年4月11日;发布日期:2017年4月14日
雷达预估信息(0~2 h)、数值天气预报产品(0~72 h)及基于高空和地面大气探测资料的综合预报信息(0~24 h)在不同预见期的降雨预报各有差异,预报效果上各有优劣,而针对山洪的预报存在准确率低、预警有效时间短等问题,研究面向山洪灾害防治区的数值预报模式的选取、模式最优物理参数化组合方案、同化方案、降水偏差订正技术以及基于大数据的自分型雷达估测降水最优化算法。以湖南临湘市为试验区进行验证,结果表明:WRF中尺度数值模式适用于山洪灾害降雨预报,其WSM6云微物理过程、Grell-Devenyi ensemble对流参数化方案和YSU边界层参数方案对致洪山洪暴雨过程模拟较好,基于频率(或面积)匹配的降水偏差订正方法能显著改善模式降水预报中雨量和雨区范围的系统性偏差,大数据分析方法应用于雷达估测降雨能显著提高准确度。在此研究基础上,提出了基于高空和地面大气探测、数值预报模式的0~72 h短期定量降雨预报和基于雷达的0~2 h临近定量降雨预报的多元信息耦合预报方法,将对致灾山洪的预报时间由2 h延长至72 h,并可显著提高了山洪预报的精度。
关键词 :致灾山洪,雷达,数值天气预报,耦合
Copyright © 2017 by authors and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
山洪灾害诱发原因较多,但暴雨量绝对是最主要的诱因。目前,在山洪灾害预报方面均依据地面实测雨量或雷达估测降水量信息 [
20世纪70年代以来,雷达估测降水方法发展迅速,但大多考虑降水类型和地理位置等影响 [
随着数值天气预报技术的发展,降水预报的精度不断提高,为增长山洪灾害预报的精度和预见期提供了可能 [
将雷达估测降水和数值天气模式结合,已有较多的研究成果 [
卫星遥感、天气雷达、高空和地面大气探测、地面自动雨量站等信息均能用于降水预报,且基于各类信息的降水预报模型也较多。目前,使用较广泛的两类模型,主要还是采用天气雷达测量信息的雷达估测降水模型以及采用高空和地面大气探测信息的数值天气预报模型。
自20世纪70年代以来,国内外对雷达估测降水方法开展了大量的研究,并在许多地区建设了天气雷达,能实现未来0~2 h降雨预估,但0~2 h的预警则显得时间太短。目前,数值模式已能实现天气预报的定时、定点和定量,不同模式、不同尺度的数值预报系统预见期可达72 h以上。如能将数值预报技术推广至山洪易发地区进行降水预报,则能使山洪灾害的有效预警时效大大延长。
研究表明 [
如果将雷达和卫星云图(0~2 h)、高空和地面大气探测(0~24 h)、数值预报产品(0~72 h)等信息,应用于山洪易发地区的山洪灾害预报,将能解决山洪预报时效性短的问题。为此,本文提出了基于雷达、高空、地面大气探测以及模式产品等多源信息融合的无缝滚动降雨预报方法,将预警时效由不足2 h提升到72 h。多源信息融合的无缝滚动降雨预报原理见图1。
数值模式和雷达外推方法的预报能力,随着预报时效的增加存在一个交叉点,交叉点以外,数值模式的预报能力要强于雷达外推方法。如何推求交叉点,许多科学家作了大量的对比分析,结论大都位于1~3 h之间。针对山洪预警需求,为简化计算工作量,本文交叉点统一取为2 h。同时为与高空、地面探测信息相配合,每日08 h、20 h开展数值模式预报,当各种时段的预报雨量值接近发生山洪临界雨量值时,开展并加密雷达外推预报频次。
图1. 多源信息融合的无缝滚动降雨预报原理图
山洪灾害预报的具体流程以及加密雷达外推预报频次的触发条件为:
1) 每日08 h、20 h开展数值天气预报,预测较大区域的72 h甚至更长预见期降雨过程;
2) 每日08 h、20 h利用最近高空、地面探测信息经过数据同化后,开展0~48 h较大区域逐小时精细化预报;
3) 将0~48 h预报成果划分为若干区域,与已有的各区域发生山洪灾害的时段(如1 h、6 h、12 h、24 h、48 h)临界雨量进行对比,若预报接近或超过临界雨量时,作为触发加密雷达外推预报频次的条件;
4) 针对可能发生山洪灾害的区域,利用雷达信息开展0~2 h临近降雨预估,特别是反演较为准确的未来30 min降雨,为山洪灾害预警提供决策依据。
数值天气预报模式在面积较小的山区应用,需要解决受局部地形影响的降水不均匀性,以及数值模式初值和模式自身存在的误差问题。数值模式初值和模式自身存在的误差主要通过选择合适的模式和预报区域实时的资料同化解决,受局部地形影响的降水不均匀性则需要利用地面雨量观测信息进行拟合率定。
由于造成山洪灾害的暴雨过程发生的时空尺度小,故采用数值天气预报模式也只能是中小尺度的。目前,可用于山洪灾害预警预报的模式主要有:中国气象局武汉暴雨研究所与中国科学院大气物理研究所的AREM模式系统、中国气象科学研究院的GRAPES-MESO模式系统以及美国的WRF模式系统。三个主要模式系统的性能指标如表1。
从模式系统的动力框架、可支持的最高分辨率、物理过程、同化系统等指标的对比可见:因静力平衡假设在中小时空尺度下不适用,需要非静力平衡的动力框架,同时还需要更完善的模式物理过程,以及更强的资料同化能力;此外,中小尺度模式需要高分辨率数值模拟和并行计算效率。
通过对以上三种模式2012年5月1日至2012年11月25日华中区域(24.5~36.53˚N, 104~117.23˚E)降水(小雨、中雨、大雨、暴雨、大暴雨)、2 m温度均方根误差、500 hPa高度场均方根误差、850 hPa风场均方根误差、850 hPa温度场均方根误差等9个要素对比评估,认为WRF模式的实际预报能力最优,因此,开展山洪灾害的数值预报试验选用WRF模式。
中小尺度数值模式的初值/边界条件一般为前6 h或者12 h全球模式的预报场,预报场的背景场往往与预报区域的观测值有一定的差异,两者之间资料同化能提高数值预报模式的精度。从资料同化的角度来说,同化可信的资料越多越能改善背景场的精度,有利于提高预报水平;而同化次数增多又会导致模式起报时间提前,模式积分时间增多,预报水平降低。如何找到这二者的合理配置,需要进行大量的试验。背景场资料一般与探空资料和地面观测资料进行同化,同化方案则考察非循环同化、不同循环同化等多种方案。表2为华中某区域同化试验设计方案。
试验Gfs0820DA0-Gfs0820DA3为考察模式不同的循环同化方案对降水的影响。Gfs0908DA0和Gfs0908DA1考察非循环同化的影响。结果表明,不管循环同化与否,同化观测资料能够改善24 h累积降水量的模拟,而循环同化中,以12 h的循环同化间隔模拟效果最好。资料同化后,即可开展累积雨量和逐时雨量分析计算。
由于数值模式初值和模式自身存在的误差,导致数值预报产品存在着一定的误差,因此,需采用适当的方法对数值模式直接输出的产品进行处理,从而改进预报效果。一般而言,降水偏差会随模式预报时效、降水不
模式 | 动力框架 | 分辨率 | 物理过程 | 同化系统 | 嵌套 | 成熟度 |
---|---|---|---|---|---|---|
AREM | 静力框架 | 公里级 | 较完善 | 客观分析/LAPS | 支持 | 发展中 |
GRAPES | 非静力框架 | 公里级 | 较完善 | 3Dvar | 支持 | 发展中 |
WRF | 非静力框架 | 米级 | 完善 | 3/4Dvar | 支持 | 较成熟 |
表1. 三个主要模式系统的性能技术指标
试验名称 | 起报时间 | GFS资料 | 同化资料 | 同化时刻 | 降水时段 | 模拟时段 |
---|---|---|---|---|---|---|
Gfs0820DA0 | 2011060820 | 00~36 h | / | / | 0908~1008 | 12~36 h |
Gfs0820DA1 | 2011060820 | 00~36 h | 探空+地面 | 0820 | 0908~1008 | 12~36 h |
Gfs0820DA2 | 2011060820 | 00~36 h | 探空+地面 | 0902,0908 | 0908~1008 | 12~36 h |
Gfs0820DA2p | 2011060820 | 00~36 h | 探空+地面 | 0820,0908 | 0908~1008 | 12~36 h |
Gfs0820DA3 | 2011060820 | 00~36 h | 探空+地面 | 0820, 0902,0908 | 0908~1008 | 12~36 h |
Gfs0908DA0 | 2011060908 | 12~36 h | / | / | 0908~1008 | 00~24 h |
Gfs0908DA1 | 2011060908 | 12~36 h | 探空+地面 | 0908 | 0908~1008 | 00~24 h |
表2. 同化试验设计方案
同临界值(阈值)和具体天气过程的不同而不同,而降水预报的准确性直接影响山洪灾害预报的准确性。
降水偏差订正可通过统计数值模式预报值和预报区域地面观测值在不同阈值条件下降水出现频率的两组数据(如图2),采用多项式插值的方法对它们进行曲线拟合、获得不同阈值模式降水预报的转换系数,即将模式降水预报的频率分布曲线调整到与观测降水频率分布曲线一致,这种保持出现频率一致的方法称为“频率匹配法”,如从空间上的分布来理解(图2(b),图2(c)),频率的大小实际上就是空间范围的大小,因为此法考虑了预报和实况雨区面积的一致性,所以“频率匹配法”也可称为“面积匹配法”,这样利用所获得的转换系数从而达到订正预报降水偏差的目的。
20世纪70年代以来,雷达估测降水技术得到长足的发展,雷达估测降水主要考虑降水类型和地理位置等的影响。由于山体对雷达电磁波的阻挡,雷达估测降水很少考虑海拔高度因素,山区随着海拔高度的变化,降水会呈现不同的变化趋势,而山区正是山洪灾害发生的高频区。为实现山洪易发区雷达估测降水,则须建立海拔高度和降水量的模型,并采用地面降水资料对雷达估测降水订正。
复杂地形雷达定量测量降水的误差来源有:1) 固定静止或近似静止的地物杂波具有较高的反射率值,导致降水量的高估,需剔除;2) 雷达不能完全将降水回波和非降水回波区分,如何评价数据质量的好坏仍然是个难题;3) 距离雷达较远的降水系统,当雷达波束从降水系统顶上扫过时,无法完整的探测到降水粒子,造成降水的过低估测;4) 不同降水类型、地区和时间的Z-I关系都存在着很大的差异,使用不合适的Z-I关系时导致降水估测误差。
因此,本文提出以Z-I最优化算法为基础,利用大数据提取数据属性间隐含信息,减少雷达估测中的主观因素,通过相异度函数检测海量雷达数据,杂波过滤清洗伪数据,DBSCAN算法聚类分析,偏差值动态修订雷达估测降雨等方法,构成基于大数据的自分型雷达估测降水最优化算法。
采用相异度函数对海量雷达数据进行清洗,并将相异度 > 10作为杂波过滤。
图2. 降水的频率和面积匹配方法原理图
式中: 为各簇内部的格点降水值; 为各簇的特征值。
采用最小化函数求聚类分析各簇的偏差值:
式中:
对雷达定量降水估计订正的基本思路是认为质量控制后的自动雨量站降水量是客观的,能代表一个区域的平均降水,而雷达定量降水估计能够反映降水场的结构,但存在量值上的系统偏差,仅仅使用未经订正的雷达回波资料会带来较大的误差。
雷达预估雨量分析流程图,如图3所示。
最优法估测降水。所有自动雨量站观测的1小时降雨量记录总数为N,观测值用Gn表示,
对降水强度进行时间积分就可以得到雷达测量的1小时降水,记为Rn,
式中:Wnm是资料所代表的时间权重系数,根据雷达观测时间确定。
图3. 雷达预估雨量分析流程图
利用自动雨量站资料,最优化估测降水方法中采用最佳判别函数CTF [
式中:Ii为雷达估计的各个样本的雨强值,Ri为自动雨量站测量的雨强值。不断的调整Z-I关系中的参数A和b的值,当CTF值为最小时,会使得雷达估测值最逼近实测值。
中国湖南省临湘市,位于北纬29˚10'~29˚52',东经113˚15'~113˚45'之间,地处东亚亚热带季风湿润气候区,属中亚热带向北亚热带过渡的边缘,年均降水量1469.1 mm,4~8月为雨季,降水量占全年的70%以上。临湘市地势南高北低,滑坡泥石流发育,由暴雨引发的山洪危害严重,引起人员伤亡,冲毁厂矿、房屋、道路、水库、良田,致使通讯、电力中断。1983年至2011年4~9月期间,直接经济损失1000万元以上且受灾人口超过10万人的暴雨山洪灾害有20次左右。
1) 数值预报模式选择
临湘市各类模式晴雨预报效率评分见图4,其中24 h AREM晴雨预报效率评分比GRAPES模式低,WRF/3d、WRF/L比GRAPES模式高。48 h AREM、WRF/3d晴雨预报效率评分比GRAPES模式低,WRF/L比GRAPES模式高。因此,该地区山洪灾害的数值预报试验选用WRF模式。
2) 资料同化及最优化参数组合方案
在暴雨的数值预报中,模式微物理过程、积云参数化方案对模拟结果的影响最为显著,为了获得山洪数值预报中最佳的物理参数组合,选取了在东亚地区应用效果较好的两组微物理过程,两组积云参数化方案,并通过4组模拟试验(其它物理参数设置一致),得出面向山洪灾害防治区的数值预报最优化参数组合方案。
采用WRFV3.4 (The Weather Research and Forecasting)模式三层双向嵌套方案(见表3),模式中心为研究区域,地图投影采用兰勃特投影,三层模式水平分辨率分别为27 km、9 km和3 km,水平格点数分别为:322*202、322*253和472*322,垂直层数为51层。使用的地形数据分别是2 m、30 s。
3) 降水偏差订正
图5为2014年5月21日08时至22日08时24小时累积观测降水和累积预报降水。从图中可以看出方案对整个雨带及其暴雨中心的落区都模拟的较好,分别将研究区域南部的暴雨中心模拟出来,特别是对临湘附近的降雨强度和落区模拟较好,但是预报强度较实况观测偏大。
从订正前后24小时、48小时预报累积降水量看,模式降水量一般存在系统性的偏差,订正后24小时预报的偏差波动范围在正负2毫米左右,48小时预报的偏差波动范围要稍大一些,降水预报的系统误差基本得到矫正。至于降水落区位置预报的准确性,通过TS和ETS评分看,0~24 h和24~48 h订正后小雨和大暴雨量级的评分都有改善,其中小雨量级的预报改善较明显。改进主要表现在小雨量段和大雨量段的原因是两端降水的偏差都较大,其中小雨量降水段有大片虚报的面积,而大雨量降水段虽然面积本身范围也许并不大,但雨量偏差较大,所以对两端订正的效果要明显一些。
图4. GRAPES、AREM15、WRF/3d与WRF/L模式24 h及48 h晴雨效率评分
微物理参数 | 对流参数 | 边界层参数 | 辐射方案 | 路面模式 |
---|---|---|---|---|
WSM 6-class graupel | Grell-Devenyi ensemble | YSU | RRTM长波、Goddard短波,15分钟计算一次 | Noah LSM |
三层嵌套;内层向外层网格反馈方案;三维变分12 h循环同化间隔模拟;基于频率(或面积)匹配的降水偏差订正。 |
表3. 山洪灾害防治区降雨数值模拟最优化参数组合方案
经过降水偏差订正后,即可判别是否达到该地区发生山洪的临界雨量,决定是否起用雷达预估短期降水。
临湘市属于岳阳市,岳阳市天气雷达台站位置见图6,以雷达站作为中心,半径150 km范围内的区域作为研究区域。
1) 雷达回波处理
采用S波段的多普勒天气雷达最大探测距离是230 km,降水产品分辨率是1 km,通过PPI扫描,可以得到不同高度上的降水强度。因为雷达探测到的原始数据存在着很多噪声,雷达回波很容易受到山脉的干扰,形成地物杂波,严重影响雷达的实测值。在复杂地形的雨量估测之前,应对雷达资料进行数据处理,主要包括孤立
图5. 2014年5月21日08时至22日08时24小时累积观测和预报降水量对比图(单位:mm)
图6. 岳阳雷达周边地形分布图
点剔除、畸异回波检测、阻挡订正、构成复合平面,其方法是通过平均值替换和分箱去噪技术进行数据清洗与数据聚类分析。
2) 雷达外推的定量降水预报
读取的雷达数据,首先经过数据清洗与数据聚类分析后,即可按照最优法估测累计降水量。自动雨量站数据的处理采用平均校准法,在整个测量区域内用一个标准订正因子去乘各点雨量的测量值,得到订正后的降水分布。读取估测偏差库的相应值作为偏差阈值,当偏差值超出偏差阈值时,通过调整Z-I反演关系直至偏差值小于偏差阈值。调整3次后若偏差值仍大于偏差阈值,则做为可疑值抛弃。
图7为雷达2014年5月21日10点45分TREC矢量图,外推30分钟后(11点15分)的雷达反射率因子图与相应时刻的回波实况图。对比分析图可知,30分钟外推回波的大体位置与实况较为吻合,雷达站西侧和北侧的强回波中心位置和范围也与实况对应较好,说明TREC能够较为真实的反映回波的移动情况,对于强对流天气提前预警具有重大的指示作用。
通过对2014年3~9月的中大型降水过程进行外推估测,得到30 min外推后的1 h累积降水量产品,并将30分钟外推的1 h累积降水量产品与12个具有代表性的雨量站降雨实况进行对比,见表4所示。
图7. 2014年5月21日10点45分TREC矢量外推的30分钟后与相应时刻的回波实况图
站号 | 平均相对误差 | 站号 | 平均相对误差 |
---|---|---|---|
P3698 | 29.8% | P2680 | 27.9% |
P3697 | 33.8% | P3666 | 32.0% |
P3696 | 40.4% | P3654 | 41.3% |
P3692 | 35.6% | P3647 | 32.0% |
P3690 | 31.6% | P3627 | 29.8% |
P3689 | 35.7% | P3583 | 34.1% |
表4. 12个雨量站降水与对应时次的外推估测产品检验
结果与实况比较表明:外推预测结果与实况相比较,在强度中心方面与实况吻合程度高,能很好的体现出降雨强度。在回波的范围、位置方面有较好的吻合,对于面积范围较大且移动稳定和变化不明显的回波在范围与位置方面,外推预测效果好于范围较小和变化较明显的回波,对流成熟阶段的外推预测正确率高于发展和消亡阶段。
1) 针对依据实测降雨开展山洪预报预警时效性短的问题,提出了基于雷达、高空、地面大气探测以及模式产品等多源信息融合的无缝滚动降雨预报方法,充分地利用雷达和卫星云图(0~2 h)、高空和地面大气探测(0~24 h)、数值预报产品(0~72 h)实现逐小时不间断无缝滚动预报,预警时效由2 h提升到72 h。
2) 研究认为WRF数值预报模式较为适用山洪灾害预报。其中,WSM6云微物理过程、Grell-Devenyi ensemble对流参数化方案和YSU边界层参数配置方案,对致洪山洪暴雨过程模拟较好。
3) 基于雷达数据速率高、周期短、数据量非常大等大数据特性,提出的通过相异度函数检测海量雷达数据,杂波过滤清洗伪数据,DBSCAN算法聚类分析,偏差值动态修订雷达估测降雨等大数据分析方法,能显著提高雷达短时预测降水的准确性。
4) 数值降水预报和雷达预测降水均需偏差修正。其中,基于频率(或面积)匹配的降水偏差订正方法,能显著改善模式降水预报中雨量和雨区范围的系统性偏差,订正后降水预报的范围和平均强度与实况更加接近,且偏差越大订正效果越好。
国家重点研发计划(批准号:2016YFC0402708)、水利行业专项经费项目(项目编号:201201063)。
熊 明,杨文发,李 俊,周北平,訾 丽. 多元信息耦合的致灾山洪降雨预报方法 Forecasting of Flash-Flood-Producing Precipitation by Coupling Multiple Sources of Information[J]. 水资源研究, 2017, 06(02): 91-102. http://dx.doi.org/10.12677/JWRR.2017.62013