U-Net模型在京津冀临近降水预报中的应用和检验评估_徐成鹏.pdf
《U-Net模型在京津冀临近降水预报中的应用和检验评估_徐成鹏.pdf》由会员分享,可在线阅读,更多相关《U-Net模型在京津冀临近降水预报中的应用和检验评估_徐成鹏.pdf(12页珍藏版)》请在咨信网上搜索。
1、第 42 卷第 6 期气象科学Vol 42,No62022 年 12 月Journal of the Meteorological SciencesDec,2022徐成鹏,曹勇,张恒德,等U-Net 模型在京津冀临近降水预报中的应用和检验评估气象科学,2022,42(6):781-792XU Chengpeng,CAO Yong,ZHANG Hengde,et al Application and test evaluation of U-Net model in Beijing-Tianjin-Hebeiprecipitation nowcasting Journal of the Mete
2、orological Sciences,2022,42(6):781-792UNet 模型在京津冀临近降水预报中的应用和检验评估徐成鹏曹勇张恒德刘海知梅双丽(国家气象中心,北京 100081)摘要利用 20182019 年期间 10 min 定量降水估计(Quantitative Precipitation Estimation,QPE)实况观测,构建基于 U-Net 的分钟级临近降水预报模型,实现了京津冀地区未来 02 h 逐 10min 降水量滚动预报。以 TS、BIAS、POD、S、FA 作为评价指标,通过检验 2020 和 2021 年 69月长序列以及分析 2020 年 8 月 12
3、 日和 2021 年 7 月 1 日两次强降水个例,表明 U-Net 模型预报接近实况,局部伴随着一定程度的空报,相较光流法、持续性预报及 CMA-MESO 模式预报效果有明显提升。具体表现为:当分钟级降水预报不超过10 mm/(10 min)时,U-Net 模型明显优于光流法和持续性预报;当小时预报不超过 25 mmh1,U-Net 模型优于 CMA-MESO 模式和光流法。然而,当降水强度超过 10 mm/(10 min)或 25 mmh1时,U-Net 模型存在预报偏弱的情况,可能与强降水样本较少有关。关键词深度学习;定量降水;临近预报分类号:P456doi:1012306/2022jm
4、s0078文献标识码:A收稿日期(eceived):2022-06-02;修改稿日期(evised):2022-06-27基金项目:国家重点研发计划项目(2021YFC3000903);风云卫星应用先行计划(2022)FY-APP-20220113通信作者(Corresponding author):曹勇(CAO Yong)caoyong cmagovcnApplication and test evaluation of UNet model inBeijingTianjinHebei precipitation nowcastingXU ChengpengCAO YongZHANG Hen
5、gdeLIU HaizhiMEI Shuangli(National Meteorological Center,Beijing 100081,China)AbstractIn order to strengthen the application of deep learning in precipitation nowcasting over theBeijing-Tianjin-Hebei region,the 10-minute QPE observation during 2018-2019 was used to construct aminute-level precipitat
6、ion nowcasting model based on U-Net,which realizes 10-minute rolling precipitationforecast in the future 0-2 hours By verifying the long series from June to September in 2020 and 2021 andanalyzing two cases of heavy precipitation on both August 12,2020 and July 1,2021 based on evaluationindicators s
7、uch as TS,BIAS,POD,S and FA,results show that the prediction of U-Net model isclose to observation accompanied by false alarms to some extent and its forecasting effect is significantlyimproved compared with the optical flow method,persistent forecast and CMA-MESO model Specificperformance as follow
8、s:when the minute-level precipitation forecast does not exceed 10 mm/(10 min),the U-Net model is significantly better than the optical flow method and persistence forecast;when thehourly forecast does not exceed 25 mmh1,the U-Net model is significantly better than the CMA-MESOmodel and the optical f
9、low method However,when the precipitation intensity exceeds 10 mm/(10 min)or 25 mmh1,U-Net has a weak forecast,which may be related to fewer samples of heavy precipitationKey wordsdeep learning;quantitative precipitation;nowcasting引言京津冀地区地处山地平原过渡地带,西毗太行山脉,东临渤海,受复杂下垫面、海陆热力差异及中低纬系统相互作用等多因素共同影响,天气系统复杂多
10、变,造成该地区降水时空分布不均,常出现短时强降水1-4,易造成重大人员伤亡和社会经济损失。京津冀强降水具有明显局地性、强突发性、短预见期的特点,数值预报模式在短中期时效内难以准确预报,亟待加强。传统短临降水预报技术是在雷达回波外推基础上发展起来的,应用最广泛的方法是单体质心法5-6 和交叉相关法7-8。单体质心法通过三维雷暴跟踪技术对雷暴单体识别和运动路径跟踪,比较适合强对流风暴单体的跟踪和临近预报,在一般对流性天气或层状云降水天气的临近预报方面应用较差9-10。交叉相关法通过计算相邻时刻雷达回波位置间的最大相关系数来跟踪和外推回波,既可以跟踪对流降水系统,也可以跟踪层状云降水,但在强对流系统
11、快速发展时风场外推质量迅速下降预报偏差变大9-11。为改进雷达回波外推效果,短临预报引入了计算机视觉领域应用较广的光流法,通过计算雷达回波的光流场代替雷达回波的运动矢量场,并基于该运动矢量场对雷达回波进行外推12-16。在预报对流性降水和稳定性降水方面,光流法都明显优于交叉相关法和质心法17-18,具有较好的业务应用价值。尽管如此,短临预报在应用中仍存在许多需提升的地方,譬如对流系统的生消和发展如何考虑,预报时效如何延长,极端降水如何反映等19-20。近年,人工智能技术迎来第三次发展浪潮,以深度学习为代表的人工智能技术在图像处理、语音识别和决策分析等领域取得了巨大成功。气象数据体量巨大,本身具
12、有天然的数据优势,基于深度学习对海量数据进行挖掘,能得到样本更深层次的内在物理特征和内在规律。当前针对特定应用场景学术界已经发展出许多成熟神经网络模型,成功应用于天气预报领域21,同时为当前迫切提升短临降水预报精度提供了新思路和新技术。短临预报的深度学习应用大致可分为两类22:(1)围绕 NN 网络开发,它显式地模拟时间变化;(2)使用 CNN 网络,将输入图像转化成输出图像。SHI,et al23 最早将短临预报归类为图像的时空序列问题,构建卷积长短记忆网络(ConvLSTM)外推02 h 雷达回波,模型预测精度较传统光流算法具有明显优势。之后,SHI,et al24 引入 Trajecto
13、ryGU 弥补记忆网络无法学习空间变化关系的缺陷,又进一步提升模型预测性能。郭瀚阳等25 将卷积GU 模型成功应用于京津冀地区雷达回波预测,结果表明模型预测未来 30 min 和 60 min 雷达回波的准确率相比传统临近预报提升明显。而另一方面,有关基于 CNN 的短临预报模型也有不少的研究进展。Ayzel,et al26 仅仅用全卷积神经网络来预测雷达短时降水,取得与光流法相当的效果。Lebedev,et al27 利用 U-Net 对卫星图像进行反演降水,继而用光流法预测未来降水量。而 Kevin,et al28 更进一步提出基于 U-Net 预测未来 30 min 的雷达回波,通过引入
14、注意力机制和深度可分离卷积,既保持较高的预测准确率又极大减少了计算参数,提高运算效率。基于 U-Net 在图形学的分割处理能力,本文尝试开展 U-Net 在短临预报方向上应用,提升短临降水预报准确性。本文设计的 U-Net 深度学习模型直接将观测的降水量作为输入量,预测京津冀地区 02 h 逐 10 min 降水量,并与光流法、持续性预报以及CMA-MESO 模式结果进行对比,给出长时间序列的检验和个例检验效果评估。1资料和方法11资料本研究区域覆盖京津冀地区、山西省大部、山东省和河南省等部分地区(3543N,111119E)。以 1 km 为单位,将研究区域大致划分成800800 个网格。基
15、于国家气象信息中心提供的 20182020年 QPE 实况资料,水平分辨率为 1 km,其中 20182019 年数据用于训练深度学习模型参数,2020 年数据用于评估模型性能。12方法深度学习通常以图像方式处理数据,对于气象格点数据,可将格点与图像中的要素点一一对应起来。深度学习分割网络 U-Net 是一种端对端的全卷积神经网络,在生物医学图像分割领域具有很好的应用效果29。如图 1 所示,U-Net 由编码阶段和解码阶段两部分组成,整体呈 U 字型对称结构:左侧为要素特征提取的编码结构,右侧为要素特征还原的解码结构。编码阶段包括卷积层和池化层,实现特征图提取和图像降维,表达气象数据的高级特
16、征。卷积层采用二维卷积提取特征,若输入尺寸为287气象科学42 卷6800800 的图片,经过 16 个 33 卷积核作用,图片尺寸变为 16800800。池化层负责降低特征图像维度,若输入特征图为 16800800,经池化后变为 16400400。解码阶段包括上采样和通道复制拼接操作,起到图像特征复原的作用,既能保留气象数据高级特征又可重建位置信息。上采样层采用二维转置卷积提高特征图分辨率,若特征图大小为 5122525,经上采样后变为 2565050。复制拼接可将编码器与解码器的的特征图做通道拼接,如编码器 2565050 的特征图与解码器 2565050 拼接后可获得的大小为 51250
17、50 的新特征图。经过重复数次上采样和复制拼接操作,最底层特征图(5122525)由低分辨率逐步还原高分辨率特征图(16800800),最后经过网络尾部卷积层输出最终预报结果(12800800)30。图 1U-Net 网络结构Fig1The network structure of U-Net相比于一般 U-Net 结构,本文构建的 U-Net 落到模型具体实现上有一些不同30:(1)在选择模型深度上,使用了五次下采样和通道拼接方式,一定程度上提升模型复杂度;(2)特征图包含的特征高级程度正比于特征图的通道数,根据网络深度特征图通道数分别设置为 16、32、64、128、256 和 512;(
18、3)卷积层后接激活函数采用 LeakyeLU 函数,解决了 eLU 的神经元死亡问题,即当输入接近于零或负值时,神经网络无法继续学习;(4)模型输入输出均是具有特定长度的时间序列,实现了基于卷积层的多对多预报模式。考虑不同量级降水样本分布极为不均衡,训练时损失函数会偏向样本多的小量级降水,造成样本稀少且较大量级降水预报效果较差。因此,在训练模型过程中使用选取加权平均 L1 函数作为损失函数,即:L1(y,y)=mi=0iy(i)y(i),(1)其中:y(i)为实际观测值;y(i)为预测值;m 为样本个数;i为样本 i 的权重系数。样本权重参考 SHI,etal28 的赋值:当降水强度小于 0.
19、5 mm/(10 min)时,权重系数设为 1.0;当降水强度位于 0.5 1 mm/(10 min)时,权重系数设为 2.0;当降水强度位于12 mm/(10 min)时,权重系数设为 5;当降水强度位于 25 mm/(10 min)时,权重系数设为 10;当降水强度大于 5 mm/10 min,权重系数设为 30。通过赋予大量级降水更高的权重,模型能更多地关注大量级降水,提高对强降水预报能力。另外,降水格点数据进入模型前经过最大最小归一化处理。训练模型最大次数设定不超过 200次。但当损失函数值在最近 10 次训练周期都没有增长,则提前停止训练。初始学习率设定为 0.001,梯度下降采用
20、Adam 优化算法。本文使用的检验评分主要包括 TS,BIAS,命中率 POD,成功率 S,空报率 FA:TS=hit/(hit+miss+false),(2)BIAS=(hit+false)/(hit+miss),(3)POD=hits/(hit+miss),(4)S=hit/(hit+false),(5)3876 期徐成鹏,等:U-Net 模型在京津冀临近降水预报中的应用和检验评估FA=false/(hit+false)。(6)其中:hit 是命中样本数;false 是空报样本数;miss 是漏报样本数。2预报效果评估2.1评分指标统计检验本文使用 2020 年和 2021 年 69 月期
21、间 10min 降水观测数据对 U-Net 模型、LK 光流法、持续性预报的三家预报结果进行对比分析(图 2、3)。U-Net模型总体表现要优于光流法和持续性预报。其中,对于 2 mm/(10 min)以上降水量,U-Net 模型的 TS 在0120 min 各时效均显著高于后两种方法;对于 5 mm/(10 min)以上降水量,U-Net 模型在 80min 左右时效依然具有一定优势,但一旦超过这个时效,则不具有相对优势,而此时三种方法的 TS 已经非常小,比较意义不大;对于降水量超过 10 mm/(10 min)的强降水,三种方法都表现较差,三种的方法 TS 在第 30 min 以后几乎接
22、近于 0,其中 U-Net 模型预报 BIAS 明显偏小,表明预报强度和范围较其他方法显著偏弱偏小。通常,华北地区降水强度超过 10 mm/10 min 属于较为极端事件,在收集样本中占比极小,因此 U-Net 模型很难学习到。22典型过程预报分析评估为进一步考察 U-Net 模型在不同天气背景下的适用性,选取 2020 年 8 月 12 日和 2021 年 7 月 1 日的华北地区强降水个例进行分析。两次降水过程分别受到不同天气系统影响,前者对应暖区内对流性强降水,后者对应冷锋面上稳定性强降水。图 22020 年 69 月逐 10 min 降水量预报 TS(ad)和 BIAS(eh)评分Fi
23、g2(a-d)TS and(e-h)BIAS scores of 10-minute precipitation forecast from June to September in 2020图 32021 年 69 月逐 10 min 降水量预报 TS(ad)和 BIAS 评分(eh)Fig3(a-d)TS and(e-h)BIAS scores of 10-minute precipitation forecast from June to September in 2021487气象科学42 卷22.12020 年 8 月 12 日实况预报对比分析2020 年 8 月 12 日,受位于华
24、北南侧暖式切变线和副高外围气流的共同影响,配合低层水汽输送,河北大部、北京、天津、山东西北部等地出现大到暴雨、局地大暴雨,形成了一条西南东北向的雨带(图 4)。12 日 03 时(北京时,下同),位于河南北部和山东西部的对流单体随着中低层水汽汇集和对流不稳定能量增加,有组织地发展为成片对流云团。08 时,该地区出现较强的降水中心,10min 累积降水量普遍在 13 mm/(10 min),部分地区可达 5 10 mm/(10 min),甚至 10 20 mm/(10 min),小时雨量在 10 40 mmh1、局地最大可达 68 mm h1。受西南风气流影响,强雨带之后逐渐北抬,先后影响河北、
25、北京、天津等大部分区域。以河北任县站为例,12 日 08、09 时的小时雨量分别为 35.7 mm 和 16.1 mm,持续降水长达 4 h(图 5)。表 1 为光流法、U-Net 模型和 CMA-MESO模式小时降水预报结果。可以看出,U-Net 模型预报 08 时小时降水较实况偏弱,而预报 09 时小时降水与实况接近。图 42020 年 8 月 12 日 08 时 850 hPa 风场和 500 hPa 位势高度(a)和 12 日 08 时13 日 08 时的 24 h 降水实况(b)Fig4(a)The 850 hPa wind field and 500 hPa geopotentia
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- Net 模型 京津冀 临近 降水 预报 中的 应用 检验 评估 徐成鹏
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。