赫尔辛基城市暴雨洪水预报系统介绍

合集下载

科技成果——2D城市洪涝与流域汛情风险预警评估系统

科技成果——2D城市洪涝与流域汛情风险预警评估系统

科技成果——2D城市洪涝与流域汛情风险预警评估系统技术开发单位Wallingford Software(英国)主要应用领域洪水灾害风险预测、防汛指挥调度成果简介1、2D城市洪涝与流域汛情风险预警评估系统,可连接实时水文与气象数据,集成复杂的水文及水力模型,实时对流域汛情以及城市洪涝灾害进行预报预警,并结合地理信息,评估洪灾的风险、洪灾的影响范围及其影响程度,帮助管理者快速准确地制定流域水资源规划、城市防汛调度应急措施,为公共报警及事故应急措施提供支持。

2、系统实现了2D洪涝灾害的预测及评估,针对防洪区域的漫堤式淹没及溃堤式淹没洪水,系统集成了一维和二维地面洪水模拟引擎;可预测城市降雨(无论是来自于河道还是城市的地下排水系统)产生的洪涝灾害,结合地面数字高程模型如DEM,预测洪水的发生地点,淹没范围,流向及淹没水深;同时系统会综合考虑地理相关信息以及城市复杂的地形、建筑物、街道以及其他设施的影响,对防洪减灾、洪水风险分析和洪涝灾害的评估都具有重要的意义。

3、系统采用了先进的雷达测量降雨技术,将雷达采集到的实时降雨数据与水文水力学模型进行结合,提高了暴雨产流的预报精度,延长了预报周期,为河流及排水系统快速响应的洪水预警,提供一种非常有效的预报机制。

4、系统具有强大的网络发布与瘦客户远程驱动功能,远程的用户一旦具有相应的权限,通过网络浏览器(如IE)能够远程驱动服务器进行水情预测,大大提高了不受地域限制的远程调度决策的能力。

主要性能指标1、1D与2D洪水引擎的完美结合,可真实地反映地形地貌要素,精确而有效地模拟洪水:可根据土地利用情况设置不同的下垫面粗糙系数,根据城市地面建筑物定义其不同的透水程度;根据地形地貌自动调节2D格网的大小,大大提高了计算的速度;更准确、真实地确定洪水的淹没范围、水深及流速。

2、为洪水风险分析及灾害评估提供准确的评判依据,无缝连接GIS系统,使洪水风险评估与地理信息紧密结合;快速、强劲的模拟与数据分析功能;可方便而充分地利用数据接口将洪水风险预测结果提供给各类灾害评估系统。

基于深度学习的雷达降雨临近预报及洪水预报

基于深度学习的雷达降雨临近预报及洪水预报

第34卷第5期2023年9月㊀㊀水科学进展ADVANCES IN WATER SCIENCEVol.34,No.5Sep.2023DOI:10.14042/ki.32.1309.2023.05.003基于深度学习的雷达降雨临近预报及洪水预报李建柱,李磊菁,冯㊀平,唐若宜(天津大学水利工程仿真与安全国家重点实验室,天津㊀300350)摘要:为探究深度学习的雷达降雨临近预报在流域洪水预报中的适用性,采用U-Net㊁嵌入注意力门的Attention-Unet 和添加转换器的多级注意力TransAtt-Unet 开展雷达降雨临近预报,将预报降雨作为HEC-HMS 水文模型的输入,对柳林实验流域进行洪水预报㊂结果表明:1h 预见期时,Attention-Unet 对短时强降雨预报结果较好,TransAtt-Unet 预报降雨模拟的洪峰流量和径流量相对误差小于20%,各深度学习模型对量级较大的降雨和洪水预报精度较高;2h 预见期的预报降雨强度㊁降雨总量㊁洪峰流量和径流量存在显著低估,U-Net 能取得相对较好的降雨预报结果㊂基于深度学习的1h 预见期雷达降雨临近预报及洪水预报可为流域防洪减灾提供科学依据㊂关键词:雷达降雨临近预报;降雨定量估计;深度学习;洪水预报;柳林实验流域中图分类号:P333㊀㊀㊀文献标志码:A㊀㊀㊀文章编号:1001-6791(2023)05-0673-12收稿日期:2023-05-19;网络出版日期:2023-09-19网络出版地址:https :ʊ /urlid /32.1309.P.20230918.1718.002基金项目:国家自然科学基金资助项目(52279022)作者简介:李建柱(1981 ),男,河北沧州人,教授,博士,主要从事水文水资源方面研究㊂E-mail:lijianzhu@近年来,极端降雨事件导致洪涝灾害频发,对洪水预报精度和时效性的要求越来越高[1]㊂准确的降雨预报是洪水预报的关键,可为防洪减灾工作提供重要的科学依据[2]㊂传统的洪水预报主要以地面雨量站实测降雨输入水文模型,导致洪水预报的有效预见期短,且雨量站实测降雨无法反映其空间分布特性[3]㊂天气雷达具有高时空分辨率㊁高精度和可靠性等特点[4]㊂将雷达降雨临近预报结果作为水文模型的输入进行洪水预报,能在一定程度上延长洪水预见期[5],是水文预报领域的主要发展趋势,但降雨临近预报是尚未解决的重要科学难题[6]㊂以光流法㊁质心跟踪法㊁交叉相关法为主的传统雷达回波外推方法[7],能在较短预见期内对缓慢变化的回波过程取得较好的外推效果,但无法准确描述迅速变化的回波过程[8]㊂近年来大量研究将深度学习的方法引入到水文气象领域,并取得显著成效[9-10]㊂Zhang 等[6]提出一种具有预测误差优化的神经网络模型NowcastNet,基于中国和美国的雷达资料开展降雨预报,显著提高了极端强降雨的预报精度㊂Ritvanen 等[11]提出一种基于卷积神经网络U-Net 的拉格朗日模型,改善了强降雨的临近预报效果㊂作为一种应用广泛的基础模型,U-Net 在临近预报领域得到推广[12]㊂Han 等[13]将U-Net 模型用于雷达回波外推,并与循环神经网络的TrajGRU 和交叉相关法的预报结果进行对比,表明U-Net 模型在时空序列预测问题的适用性㊂预报降雨的重要用途之一是作为水文模型的输入进行洪水预报[14-15]㊂Heuvelink 等[16]使用确定性和概率性方法进行降雨预报,并将结果作为WALRUS 水文模型的输入进行流量预报,在小流域取得较好的预报效果;Nguyen 等[17]将雷达回波外推和数值模式天气预报结合,提高了降雨预报精度,并指出降雨预报与分布式水文模型结合使用的优势;包红军等[18]构建了临近降雨集合预报的中小河流洪水预报模型,延长了洪水预报的预见期㊂但目前使用深度学习的方式开展雷达降雨临近预报,从而进行洪水预报的精度还有待提高,尤其是在半干旱半湿润地区的小流域㊂U-Net 结构简单,可根据目标灵活调整和添加模块,为实现更加精确的雷达降雨临近预报和洪水预报提供了更多可能㊂因此,本研究采用深度学习的U-Net,并尝试使用嵌入注意力门的Attention-Unet(简写为Att-674㊀水科学进展第34卷㊀Unet)和添加转换器的多级注意力TransAtt-Unet进行柳林实验流域典型降雨过程的1h和2h预见期雷达回波外推,利用动态雷达反射率因子和降雨强度关系计算外推回波对应的逐小时降雨,作为半分布式水文模型HEC-HMS的输入进行洪水预报,对比分析预报降雨与雨量站实测降雨模拟的洪水精度差异,探讨基于深度学习的雷达降雨临近预报在小流域洪水预报中的适用性㊂1㊀研究区域与数据1.1㊀流域概况选择河北省邢台市内丘县柳林实验流域为研究区㊂如图1所示,流域出口位于114ʎ21ᶄE㊁37ʎ17ᶄN,布设有柳林水文站进行水位和流量观测㊂流域面积为57.4km2,主河道长13.2km,流域坡度为30.9ɢ㊂地处半干旱半湿润气候区,季节变化分明,径流的年内和年际分配不均,降雨多发生在6 9月,多年平均降水量为594.5mm,多年平均径流深为80.7mm㊂雷达资料来源于中国新一代多普勒天气雷达监测网河北省石家庄市的Z9311站㊂雷达位于114ʎ42ᶄ50ᵡE㊁38ʎ21ᶄ00ᵡN,采用VP21体扫模式,扫描半径为230km,时间分辨率为6min,能够完成9个不同仰角的扫描㊂柳林实验流域在距离雷达120km范围内,可保证雷达回波数据的质量㊂图1㊀柳林实验流域位置和雷达监测范围Fig.1Location of Liulin experimental watershed and radar monitoring range1.2㊀雷达数据预处理将2018 2020年降雨时段的雷达基数据,经过编码转换㊁杂波抑制㊁衰减订正㊁地物遮挡订正㊁坐标转换后,形成反射率混合扫描数据图,用于深度学习模型的训练㊂由于柳林实验流域面积较小,将回波图裁剪至流域周围128行128列的范围(113ʎ36ᶄ36ᵡE㊁37ʎ30ᶄ36ᵡN到114ʎ53ᶄ24ᵡE㊁36ʎ47ᶄ24ᵡN)㊂以训练集中20180521场次降雨为例,图2为经过预处理与未处理回波对比图,经过质量控制后的回波剔除了杂波干扰并显著减少波束遮挡,具有较高的可靠性㊂训练集和验证集共包含20000帧回波图,按照8ʒ2的比例划分进行模型训练和验证㊂测试集选取2012年㊁2016年和2021年的4场典型降雨过程对应的雷达数据㊂1.3㊀典型降雨洪水过程降雨洪水资料来源于河北省邢台水文勘测研究中心㊂洪水资料为柳林水文站汛期实测逐小时流量数据,㊀第5期李建柱,等:基于深度学习的雷达降雨临近预报及洪水预报675㊀图2㊀20180521场次降雨未处理与预处理后回波对比Fig.2Comparison of unprocessed and preprocessed echoes of20180521case暴雨资料为流域内菩萨岭㊁神头㊁任庄㊁安上和柳林5个雨量站汛期逐小时雨量,筛选出与雷达回波时段对应的典型降雨洪水过程对降雨和洪水预报精度进行评价㊂降雨洪水信息如表1所示㊂由于20211006场次降雨过程在汛期后发生,柳林水文站未对其洪水过程进行观测㊂表1㊀典型降雨洪水过程信息Table1Information of typical rainfall and flood processes降水场次降水量/mm最大小时雨量/mm降雨时长/h洪水径流深/mm洪峰流量/(m3㊃s-1) 2012072675.958.45 5.147.020160719251.942.718129.2368.020210721149.426.92929.993.32021100625.58.25//2㊀研究方法2.1㊀雷达回波外推的深度学习模型2.1.1㊀U-NetU-Net网络由4层编码器-解码器组成,是一种卷积神经网络㊂如图3(a)所示,网络的左边是编码器,应用最大池化和双重卷积来减小图像大小和加倍特征映射的数量㊂编码器之后的右侧为解码器,通过双线性插值进行上采样操作,使特征图大小增加1倍㊂每层编码器和解码器之间通过1个跳跃连接保存来自较浅层的细尺度信息㊂完成上述采样操作之后,模型通过一个1ˑ1的卷积,输出代表网络预测值的单个特征图㊂2.1.2㊀Att-UnetAtt-Unet(图3(a))在U-Net的解码器前添加注意力门,以此过滤跳跃连接传播的特征,再将编码器的特征与解码器中相应的特征进行拼接,有效抑制无关区域的激活,减少编码器中无关信息的跳跃连接,达到改善预测效果的目的[19]㊂本研究将Att-Unet模型调整为时间序列预测模型进行回波外推㊂676㊀水科学进展第34卷㊀图3㊀U-Net㊁Att-Unet㊁TransAtt-Unet和TSA㊁GSA模块结构Fig.3Structure of U-Net,Attention-Unet,TransAtt-Unet,TSA and GSA modules2.1.3㊀TransAtt-UnetTransAtt-Unet将多层次引导注意和多尺度跳跃连接联合嵌入U-Net,如图3(b)㊂将变换器自注意力(TSA)和全局空间注意力(GSA)嵌入到网络中,同时在解码器中使用多尺度跳跃连接来聚合不同语义尺度的特征,从而有效减少卷积层叠加和连续采样操作造成的细节损失㊂如图3(c)和图3(d)所示,TSA将特征嵌入到Q㊁K㊁V3个矩阵中,在Q和K的转置之间采用Softmax函数进行归一化运算,形成注意力图,再与V 矩阵相乘得到注意力权重㊂GSA对特征进行卷积转置映射为W㊁M㊁N,对M㊁N采用Softmax函数进行归一化运算得到位置注意力信息,再与W相乘得到位置特征㊂在解码器部分采用残差多尺度跳跃连接的方式[20],输入特征图通过双线性插值向上采样到输出的分辨率,然后与输出特征图进行级联,作为后续块的输入㊂2.2㊀深度学习模型训练1h外推数据集取前1h间隔6min共10帧回波图(反射率因子),预测后1h共10帧回波图(2h外推为前20帧预测后20帧)㊂深度学习模型基于Pytorch环境,初始学习率设置为0.001,批处理大小设置为8,损失函数采用均方根误差(E MS)[13],在NVIDIA Geforce RTX3050上采用Adam优化器训练200个轮次㊂当损失函数在4个周期内没有增加时,学习率调节器将学习率自动减小10%㊂采用二元评价指标命中率(Proba-bility of detection,D PO)㊁虚警率(False alarm ratio,R FA)㊁临界成功指数(Critical success index,I CS)和准确率(Accuracy,A)进行雷达回波外推精度的评价[21]㊂2.3㊀雷达降雨临近预报多普勒天气雷达采用Z R关系描述雷达反射率因子(Z)和降雨强度(R)的幂指数关系[22]㊂中国的多普勒雷达普遍采用Z=aR b(a=300,b=1.4)进行降雨定量估计,但仅适用于平均情况㊂基于实测资料动态调整的Z R关系,可以实现更加精确的降雨估计[23]㊂殷志远等[24]采用4种不同的Z R关系开展雷达降雨定量估计,并将结果用于水文模拟,表明动态Z R关系的降雨定量估计精度最高,洪水模拟效果最好㊂动态Z R关系建立在逐小时快速更新资料的基础上,通过动态调节参数a和b,使逐小时雷达估测降雨与对应㊀第5期李建柱,等:基于深度学习的雷达降雨临近预报及洪水预报677㊀的雨量站观测降雨的最优判别函数δ达到最小[23],从而确定适用于逐小时雷达定量降雨估计的多组Z R关系参数㊂为保证参数a和b的取值合理,限定a和b数值调节范围分别为[150.00,400.00]㊁[0.80, 2.40],调整间隔分别为10和0.05㊂经过上述步骤最终确定出每场降雨过程的动态Z R关系参数如图4所示㊂采用相关系数(Correlation coefficient,C C)㊁平均偏差(Mean bias,B M)和平均绝对误差(Mean absolute error,E MA)进行降雨预报精度评价[14]㊂图4㊀动态Z R关系参数Fig.4Parameters of dynamic Z R relationship2.4㊀水文模型李建柱等[25]研究了地形数据源和分辨率对柳林实验流域洪水模拟精度的影响,结果表明,基于无人机三维倾斜摄影构建的1m分辨率DEM能反映流域真实地形的变化,在此基础上构建的HEC-HMS模型能较好地模拟流域洪水过程㊂本研究采用作者基于无人机三维倾斜摄影构建的1m分辨率HEC-HMS模型,将雷达降雨临近预报结果作为HEC-HMS水文模型的输入,进行柳林实验流域洪水预报㊂采用洪峰流量相对误差(E RP)㊁径流量相对误差(E RV)㊁峰现时差(ΔT)和纳什效率系数(E NS)进行洪水预报精度评价[25]㊂3㊀结果及分析3.1㊀回波外推结果分析以反射率20dBZ和30dBZ为阈值计算二元评价指标㊂1h回波外推评价结果见表2,Att-Unet对20120726场次降雨的回波外推效果相对较好,注意力门加强了Att-Unet对强回波的识别和外推效果,同时抑制弱回波或杂波产生的干扰,但对弱回波或中等回波的外推效果较差;TransAtt-Unet采用的多尺度跳跃连接和注意力机制使模型能综合不同尺度的图像特征,提高模型精度和稳定性,因此,该模型对于持续时间较长㊁降雨过程变化丰富的20160719场次回波过程取得了较好的外推效果;20210721和20211006场次降雨的过程回波总体偏弱,各模型的1h预见期回波外推精度差异并不显著㊂2h回波外推评价结果见表3,U-Net 模型对30dBZ阈值的回波外推效果优于其他模型,其原因是雷达回波外推需要预测每个像素的精确值;Att-Unet和TransAtt-Unet所采用的注意力门或者多尺度跳跃连接结构仅增强局部特征的学习,而忽略随时间动态变化的信息,因此导致预测时效性的不足㊂U-Net模型尽管结构简单,但以往研究表明其在时间序列预测中具有一定适用性[26],对不同等级回波信息具有一定的泛化能力[27],因此,尽管U-Net的1h预见期回波外推效果略差于添加注意力机制的模型,但能在更长预见期的回波外推中保持相对较好的效果㊂总体来看, 3种模型对中等强度回波外推效果均好于强回波,1h预见期回波外推效果好于2h预见期㊂国内外研究主要依靠天气雷达外推实现1h预见期降雨预报[5]㊂与传统的雷达外推方法相比,深度学习678㊀水科学进展第34卷㊀对回波和降水的演变趋势具有更好的预报效果,更适用于剧烈变化的降雨过程[28]㊂曹伟华等[29]使用基于U-Net网络搭建的RainNet模型开展雷达降雨临近预报,并与交叉相关的外推结果进行对比,指出了深度学习模型对降雨消亡过程的时空演变趋势和强度变化范围具有更好的预报效果,而交叉相关法更适合于稳定降雨的预报㊂本研究预报的4场典型降雨过程,除20211006场次持续时间短㊁降雨强度较小外,其余场次降雨过程变化较为剧烈,回波过程变化较为迅速,因此,采用深度学习的方法进行雷达降雨临近预报更为合适㊂表2㊀1h回波外推结果评价指标值Table2Evaluation index value of1h echo extrapolation results降雨场次模型D PO R FA I CS A20dBZ30dBZ20dBZ30dBZ20dBZ30dBZ20dBZ30dBZ20120726U-Net0.830.370.070.190.780.330.820.56 Att-Unet0.820.520.120.170.760.440.800.65 TransAtt-Unet0.750.350.060.170.720.340.770.5820160719U-Net0.710.370.140.440.610.210.610.31 Att-Unet0.790.460.150.450.680.240.690.33 TransAtt-Unet0.870.470.140.490.750.240.750.3120210721U-Net0.740.420.140.210.670.360.810.71 Att-Unet0.700.490.260.440.540.280.640.51 TransAtt-Unet0.720.520.260.460.550.290.650.5320211006U-Net0.450.240.560.760.190.090.740.64 Att-Unet0.440.210.550.790.240.080.750.61 TransAtt-Unet0.430.280.560.720.210.090.690.63表3㊀2h回波外推结果评价指标值Table3Evaluation index value of2h echo extrapolation results降雨场次模型D PO R FA I CS A20dBZ30dBZ20dBZ30dBZ20dBZ30dBZ20dBZ30dBZ20120726U-Net0.660.370.140.160.600.350.640.63 Att-Unet0.620.180.140.170.590.180.650.48 TransAtt-Unet0.660.200.110.130.590.200.630.5120160719U-Net0.590.200.420.630.380.120.430.47 Att-Unet0.480.120.440.640.310.070.370.47 TransAtt-Unet0.570.150.360.600.410.100.420.4820210721U-Net0.690.480.140.330.620.330.770.62 Att-Unet0.600.360.250.390.460.220.590.51 TransAtt-Unet0.650.380.290.460.480.220.570.4820211006U-Net0.500.170.480.770.320.080.670.58 Att-Unet0.420.180.670.810.210.020.600.54 TransAtt-Unet0.410.170.610.800.200.040.640.52㊀第5期李建柱,等:基于深度学习的雷达降雨临近预报及洪水预报679㊀3.2㊀降雨预报精度分析利用动态Z R关系计算各降雨场次真实回波和外推回波的逐小时降雨,精度和相关性评价指标见表4㊂20160719㊁20210721和20211006场次真实回波的反演降雨与实际降雨有较高的相关性,但20120726场次真实回波反演降雨与实际降雨的相关性较弱,但该场次回波外推取得较高的技能评分,其原因可能是2012年Z9311雷达为单偏振雷达,所采集的原始回波数据存在一定的数值和发生时间的系统误差,且该场降雨过程变化迅速,导致预报降雨与实际产生误差㊂TransAtt-Unet对20160719和20210721场次降雨在1h预见期内具有较小的误差,Att-Unet在20160719场次降雨的预报中具有最高的相关性㊂3种模型对20211006场次降雨预报效果差别并不显著,其原因是该场降雨强度较小,深度学习模型对此类降雨预报性能较接近㊂在2h 预见期降雨预报中,U-Net模型的降雨预报效果优于其他模型,这与雷达回波外推结果相对应,表明了U-Net在较长预见期降雨预报中的适用性㊂表4㊀降雨相关性评价指标值Table4Rainfall correlation evaluation index value降雨场次模型E MA/mm C C B M/mm1h2h1h2h1h2h20120726U-Net22.5722.49-0.28-0.41-1.96-2.32 Att-Unet23.9921.28-0.22-0.32-0.66-3.08 TransAtt-Unet24.2022.66-0.28-0.32-1.76-2.66真实回波14.780.31 2.9320160719U-Net8.868.840.480.420.22-6.21 Att-Unet8.7010.710.590.170.25-7.05 TransAtt-Unet7.8810.520.580.13-0.74-6.77真实回波 6.020.77-0.0620210721U-Net 4.57 4.780.520.380.31-1.90 Att-Unet 4.85 5.580.320.19-0.93-1.31 TransAtt-Unet 4.40 5.650.470.21-0.22-1.29真实回波 3.090.67-0.4120211006U-Net 2.81 2.650.840.61-2.81-2.39 Att-Unet 2.61 2.790.790.63-2.61-2.79 TransAtt-Unet 2.48 2.620.820.65-2.48-2.62真实回波 1.670.88-1.67㊀㊀预报降雨过程如图5所示㊂3种模型在1h预见期时,预报的20120726场次降雨峰值与实际较为一致,但出现1h时差;对20160719和20210721场次降雨过程预报结果出现部分异常值,这与回波外推过程较大的虚警率有关,但总体上能反映降雨过程变化和雨强峰值;20211006场次降雨则存在少量低估,但能预报出该场降雨峰值出现的时间,这与动态Z R算法对较弱降雨的系统性低估有关㊂2h预见期降雨可以一定程度预报降雨过程的变化,但对各场降雨的峰值存在显著低估㊂目前,小流域降雨临近预报效果普遍较差㊂Heuvelink等[16]采用拉格朗日持续性方法在一个40km2的680㊀水科学进展第34卷㊀流域上对一场强降雨进行预报,产生了50%相对误差,发现面积越小的流域对过程变化迅速的降雨越容易产生误报;石毅[30]采用Farneback光流法和ConvLSTM在柳林实验流域进行降雨预报,结果表明光流法对回波演变的敏感性相对较低,ConvLSTM对强回波存在显著的均化趋势,导致1h预见期强降雨存在严重的低估㊂本研究采用的深度学习方法,在1h预见期内对不同类型的降雨均取得相对较好的回波外推效果,且能较准确的预报出强降雨峰值和变化过程,尽管2h预见期的预报精度相对较差,但能预报出降雨变化过程㊂预报结果存在的误差与定量降雨估计方法的系统误差和小流域上有限的雷达回波信息相关㊂动态Z R关系在定量降雨估计中具有相对较高的精度,但在降雨预报的业务化应用中仍然具有优化的空间,如Mihulet等[31]使用机器学习的方法改善了定量降雨估计的效果㊂另外,由于流域面积较小,随着预见期的延长降雨发生的实际位置也许出现在雷达图之外,使得深度学习的方法对剧烈变化的降雨产生较大的误差,Heuvelink等[16]也指出面积越小的流域对降雨发生的位置敏感性越高㊂图5㊀预报逐小时降雨过程Fig.5Forecasted hourly rainfall process3.3㊀洪水预报精度分析表5为洪水预报精度评价结果㊂3场实测降雨模拟的洪峰流量均小于实测洪峰流量,但E NS均达到了0.7以上㊂预报降雨模拟的20120726和20210721场次峰量较小的洪水,E NS均小于0.3,但20120726场次洪水1h预见期径流量相对误差小于20%,20210721场次洪水的洪峰流量预报效果也好于实测降雨模拟结果㊂对于20160719场次峰量较大的洪水,1h预见期预报的洪水E NS均能达到0.7以上且预报洪峰流量相对误差均小于20%,满足预报的精度要求㊂2h预见期洪水预报效果显著变差,洪峰流量和径流量的预报也存在较大误差㊂1h预见期时,TransAtt-Unet对于3场洪水的预报洪峰流量和径流量相对误差均小于20%,且20160719场次洪水E NS达0.78;Att-Unet则较准确地预报出20210721场次洪水的洪峰流量,相对误差仅为-0.9%㊂由于2h预见期预报的降水量存在显著低估,使得预报洪峰流量显著小于实测值,但U-Net模型对3场洪水预报的E NS为3个深度学习模型的最优值,且预报的20160719场次洪水E NS达0.52㊂㊀第5期李建柱,等:基于深度学习的雷达降雨临近预报及洪水预报681㊀表5㊀洪水预报精度评价结果Table5Accuracy evaluation of forecasted floods洪水场次模型预见期Q S/(m3㊃s-1)E RP/%V S/mm E RV/%ΔT/h E NS20120726U-Net1h31.3-33.4 4.4-13.610.23 2h19.2-59.2 4.7-8.010.16Att-Unet1h40.3-24.5 5.68.810.06 2h17.8-62.1 3.9-43.630.08TransAtt-Unet1h38.0-19.2 5.3 2.210.23 2h19.0-59.6 4.3-17.110.15实测降雨模拟36.7-21.97.0537.200.7220160719U-Net1h372.4 1.2162.525.810.76 2h198.4-46.189.3-30.930.52Att-Unet1h407.310.7165.928.410.77 2h185.3-49.779.3-38.640.34TransAtt-Unet1h354.9-3.6154.819.910.78 2h191.3-48.183.9-35.030.37实测降雨模拟343.3-6.7167.729.810.8320210721U-Net1h123.232.033.010.400.25 2h32.3-65.47.7-74.300.14Att-Unet1h92.5-0.917.9-40.110.04 2h41.6-55.411.6-61.31-0.09TransAtt-Unet1h101.08.325.7-14.000.19 2h57.0-38.912.8-57.210.02实测降雨模拟76.2-18.327.7-7.300.79㊀㊀预报洪水过程线如图6所示㊂3种模型1h预见期洪水变化过程与实际较为一致,预报的20120726场次洪水峰现时间和洪水涨落时间较实际滞后1h;20160719场次预报洪水与实测降雨模拟的峰现时间均较实际滞后1h,预报洪峰与实际较为接近,但径流量存在一定的高估;对于20210721场次洪水预报,U-Net预报的洪峰流量较实际偏大,Att-Unet预报的峰现时间较实际滞后1h,TransAtt-Unet的预报结果与实际更为接近,3种模型均能预报出该场洪水的涨落过程㊂深度学习的方法对于剧烈变化的降雨引发的洪水,1h预见期的预报洪水E NS较低,但能较好地预报出洪水的变化过程和洪峰流量,对量级较大的洪水能取得较高的E NS,且能在准确预报洪水变化过程的基础上,较准确地预报出洪峰流量㊁径流量和峰现时间㊂2h预见期降雨虽然可以预报出洪峰形成过程,但对洪峰流量和径流量存在显著低估㊂糜佳伟等[32]在梅溪流域(面积约956km2)进行降雨预报和洪水预报,指出1h预见期降雨预报结果能满足中小流域洪水预报需求㊂本研究在降雨径流响应时间更快的柳林实验流域进行洪水预报,尽管预报洪水E NS较小,但能在1h预见期对不同类型降雨引发的洪水取得较为准确的洪峰流量和径流量预报效果,预报的20160719场次大洪水的洪峰流量和径流量相对误差小于实测降雨模拟洪水结果,因此,1h预见期洪水预报效果具有一定的准确性,为流域的防洪减灾工作争取了更长的时间㊂未来可在更多流域开展雷达降雨临近预报和洪水预报研究,以验证本文采用的深度学习方法在其他流域的适用性㊂682㊀水科学进展第34卷㊀图6㊀模拟和预报洪水过程线Fig.6Simulated and forecasted flood hydrographs4㊀结㊀㊀论采用深度学习的U-Net㊁Att-Unet和TransAtt-Unet进行雷达回波外推,通过动态雷达反射率因子和降雨强度关系实现雷达降雨临近预报,将降雨预报的结果输入HEC-HMS水文模型对柳林实验流域典型洪水过程进行预报,得到以下主要结论:(1)1h预见期时Att-Unet对强回波过程外推效果较好,TransAtt-Unet对变化更丰富的回波过程外推效果较好;2h预见期时U-Net外推效果更稳定㊂(2)深度学习模型在1h预见期对短时强降雨存在时间上的误差,对持续时间较长的降雨存在少量预报异常值,但均能较准确地预报降雨强度和过程;2h预见期降雨存在显著低估和较大误差㊂(3)3种模型的1h预见期预报的洪水能反映实际变化过程,TransAtt-Unet预报的洪峰流量和径流量误差更小,Att-Unet能对部分场次洪水取得较准确的洪峰预报效果㊂U-Net在2h预见期洪水预报效果精度最高㊂参考文献:[1]雍斌,张建云,王国庆.黄河源区水文预报的关键科学问题[J].水科学进展,2023,34(2):159-171.(YONG B, ZHANG J Y,WANG G Q.Key scientific issues of hydrological forecast in the headwater area of Yellow River[J].Advances in Water Science,2023,34(2):159-171.(in Chinese))[2]金君良,舒章康,陈敏,等.基于数值天气预报产品的气象水文耦合径流预报[J].水科学进展,2019,30(3):316-325.(JIN J L,SHU Z K,CHEN M,et al.Meteo-hydrological coupled runoff forecasting based on numerical weather prediction products[J].Advances in Water Science,2019,30(3):316-325.(in Chinese))[3]IMHOFF R O,BRAUER C C,van HEERINGEN K J,et rge-sample evaluation of radar rainfall nowcasting for flood early warning[J].Water Resources Research,2022,58(3):e2021WR031591.[4]IMHOFF R O,de CRUZ L,DEWETTINCK W,et al.Scale-dependent blending of ensemble rainfall nowcasts and numerical weather prediction in the open-source pysteps library[J].Quarterly Journal of the Royal Meteorological Society,2023,149(753):㊀第5期李建柱,等:基于深度学习的雷达降雨临近预报及洪水预报683㊀1335-1364.[5]刘佳,邱庆泰,李传哲,等.降水临近预报及其在水文预报中的应用研究进展[J].水科学进展,2020,31(1):129-142.(LIU J,QIU Q T,LI C Z,et al.Advances of precipitation nowcasting and its application in hydrological forecasting[J]. Advances in Water Science,2020,31(1):129-142.(in Chinese))[6]ZHANG Y C,LONG M S,CHEN K Y,et al.Skilful nowcasting of extreme precipitation with NowcastNet[J].Nature,2023, 619(7970):526-532.[7]SOKOL Z,SZTURC J,ORELLANA-ALVEAR J,et al.The role of weather radar in rainfall estimation and its application in mete-orological and hydrological modelling:a review[J].Remote Sensing,2021,13(3):351.[8]HAN L,ZHAO Y Y,CHEN H N,et al.Advancing radar nowcasting through deep transfer learning[J].IEEE Transactions on Geoscience and Remote Sensing,2022,60:1-9.[9]BI K F,XIE L X,ZHANG H H,et al.Accurate medium-range global weather forecasting with3D neural networks[J].Nature, 2023,619(7970):533-538.[10]李步,田富强,李钰坤,等.融合气象要素时空特征的深度学习水文模型[J].水科学进展,2022,33(6):904-913.(LI B,TIAN F Q,LI Y K,et al.Development of a spatiotemporal deep-learning-based hydrological model[J].Advances in Water Science,2022,33(6):904-913.(in Chinese))[11]RITVANEN J,HARNIST B,ALDANA M,et al.Advection-free convolutional neural network for convective rainfall nowcasting[J].IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing,2023,16:1654-1667. [12]AYZEL G,SCHEFFER T,HEISTERMANN M.RainNet v1.0:a convolutional neural network for radar-based precipitation now-casting[J].Geoscientific Model Development,2020,13(6):2631-2644.[13]HAN L,LIANG H,CHEN H N,et al.Convective precipitation nowcasting using U-net model[J].IEEE Transactions on Geo-science and Remote Sensing,2022,60:1-8.[14]熊立华,刘成凯,陈石磊,等.遥感降水资料后处理研究综述[J].水科学进展,2021,32(4):627-637.(XIONG LH,LIU C K,CHEN S L,et al.Review of post-processing research for remote-sensing precipitation products[J].Advances in Water Science,2021,32(4):627-637.(in Chinese))[15]刘家宏,梅超,刘宏伟,等.特大城市外洪内涝灾害链联防联控关键科学技术问题[J].水科学进展,2023,34(2):172-181.(LIU J H,MEI C,LIU H W,et al.Key scientific and technological issues of joint prevention and control of river flood and urban waterlogging disaster chain in megacities[J].Advances in Water Science,2023,34(2):172-181.(in Chi-nese))[16]HEUVELINK D,BERENGUER M,BRAUER C C,et al.Hydrological application of radar rainfall nowcasting in the Netherlands[J].Environment International,2020,136:105431.[17]NGUYEN H M,BAE D H.An approach for improving the capability of a coupled meteorological and hydrological model for rain-fall and flood forecasts[J].Journal of Hydrology,2019,577:124014.[18]包红军,曹勇,曹爽,等.基于短时临近降水集合预报的中小河流洪水预报研究[J].河海大学学报(自然科学版),2021,49(3):197-203.(BAO H J,CAO Y,CAO S,et al.Flood forecasting of small and medium-sized rivers based on short-term nowcasting and ensemble precipitation forecasts[J].Journal of Hohai University(Natural Sciences),2021,49(3): 197-203.(in Chinese))[19]GAO Y B,GUAN J P,ZHANG F H,et al.Attention-unet-based near-real-time precipitation estimation from Fengyun-4A satel-lite imageries[J].Remote Sensing,2022,14(12):2925.[20]FANG J,YANG C,SHI Y T,et al.External attention based TransUNet and label expansion strategy for crack detection[J].IEEE Transactions on Intelligent Transportation Systems,2022,23(10):19054-19063.[21]TREBING K,STA CZYK T,MEHRKANOON S.SmaAt-UNet:precipitation nowcasting using a small attention-UNet architec-ture[J].Pattern Recognition Letters,2021,145:178-186.[22]ZOU H B,WU S S,TIAN M X.Radar quantitative precipitation estimation based on the gated recurrent unit neural network andecho-top data[J].Advances in Atmospheric Sciences,2023,40(6):1043-1057.[23]GOU Y B,CHEN H N,CHANDRASEKAR V.A dynamic approach to quantitative precipitation estimation using multiradar mul-tigauge network[J].IEEE Transactions on Geoscience and Remote Sensing,2020,58(9):6376-6390.。

211219467_基于机器学习快速预报模型的城市洪涝预报预警系统研究及应用

211219467_基于机器学习快速预报模型的城市洪涝预报预警系统研究及应用

In the process of terrain generation,there are artificial corrections,and urban construction has a certain impact on local terrain.There⁃fore,it is necessary to survey and correct local terrain abnormal areas,and optimize the model based on actual waterlogging data.The ur⁃ban flood forecasting and decision-making system can achieve rapid and accurate urban flood forecasting,providing support for urban flood early warning and emergency management.Key words:urban flooding;GAST hydrodynamic numerical model;machine learning algorithms;forecast and early warning0 前 言气候变化结合城市 热岛效应”使得城市大气层结构不稳,对流性极端降水频发,近年来极端暴雨导致城市洪涝灾害频发[1]㊂如北京2012年 7㊃21”特大暴雨[2]㊁武汉2016年 7㊃6”大暴雨[3]㊁广州2020年 5㊃22”特大暴雨[4]以及河南郑州2021年特大暴雨灾害[5],尤其是2021年宁波 烟花”台风造成的城市大面积洪涝灾害[6]㊂随着智慧城市㊁孪生城市概念的提出与发展,全国各大城市相继开展了智慧城市建设,为城市管理构建一个智慧大脑,帮助城市更加健康㊁智能㊁高效的发展㊂其中,智慧城市综合管控平台建设核心之一是智慧防汛管理[7]㊂城市智慧防汛应急管理通常以城市管理传感器构建的物联网IOT㊁高速宽带有线无线互联网(5G)㊁数字孪生㊁大数据㊁云平台㊁地理信息系统GIS㊁预测预报㊁指挥调度等最新技术和排水防涝需求为支撑,形成 城市洪涝数字孪生体”,实现排水设施智慧管控及运维㊁洪涝灾害预报预警和防汛排涝智能调度,其中城市洪涝预报系统中的洪涝过程精确模拟计算过程是重点,模型的合理构建是关键㊂近年来,许多研究者针对城市洪涝预报开展研究[8-10],林珲等[11]指出城市洪涝预报系统的关键是具有能对整个暴雨洪涝过程进行精细化模拟的城市洪涝预报模型㊂孙永尚等[12]基于水位传感器等物联网技术构建北京市内涝积水监测预警系统;李海宏等[13]基于暴雨内涝灾情数据构建暴雨内涝灾情预测BP模型,对内涝灾情影响等级进行阈值划分,为暴雨内涝影响预报与风险预警业务㊁服务及灾害管理提供技术支撑㊂但上述模型均为经验模型,即利用实测监测和历史积涝数据,通过机器学习算法实现对未来降雨内涝情况的预测㊂显然,经验模型缺乏物理机制的表征,不能模拟城市暴雨㊁洪水导致的地表㊁河道及管网的全过程流动情况,更无法准确预测出内涝位置和淹没程度,尤其是针对历史未曾发生的极端暴雨事件,缺乏一个系统的可视化的城市洪涝快速预报预警系统,为城市雨洪管理者提供决策支撑㊂为实现对暴雨导致的城区洪涝全过程的精细预测,本文城市水循环规律深入研究的基础上,基于自主研发的城市洪涝水文水动力模型模拟暴雨-洪水-治涝过程,弥补历史降雨积涝数据的不足,通过设计各种降雨情景获取全面的降雨积涝数据,采用机器学习算法建立降雨-积涝关系[16],构建城市洪涝预报决策系统,实现城市洪涝快速地准确预报,为城市洪涝预警和应急管理提供支撑㊂1 系统构建城市洪涝预报预警平台以自主研发的城市洪涝水文水动力模型为核心,以精细化气象降雨数据驱动数学模型,结合人工智能快速预报算法搭建特征库,实现洪涝预报㊁洪涝预演㊁洪涝评估为核心功能为城市雨洪管理提供科学管理和辅助决策支持,形成洪涝管理信息智能综合管理㊂1.1 预报模型城市洪涝的快速预报模型主要依托于两个工具,一个是物理雨洪模型,另一个是机器学习算法㊂物理雨洪模型为机器学习提供可靠的训练以及测试数据,机器学习算法为快速预报模型提供基础支撑,使其在预报过程中不再反复迭代求解复杂方程组,而是通过算法所建立起来的联系直接得出结果㊂1.1.1 平原城区洪涝数值模拟方法城市内涝过程数值模型采用西安理工大学开发的GAST水动力数值模型[14-15,17]㊂模型控制方程为二维浅水方程(SWEs),采用Godunov格式求解,并使用高性能GPU并行技术进行加速㊂其二维非线性浅水方程守恒格式的矢量形式如下:əqət+əFəx+əGəy=S(1) qhq xqéëêêêêùûúúúúyF= uh uq x+gh2/2 uq yéëêêêêùûúúúúG= vh vq x vq x+gh2/éëêêêêùûúúúú231西北水电㊃2023年㊃第2期===============================================S = i -gh əz b /əx -gh əz b /əéëêêêêùûúúúúy + 0 -C f u u 2+v 2-C f v u 2+v éëêêêêùûúúúú2(2)式(1)~(2)中:q 为变量矢量,包括水深h ,两个方向的单宽流量q x 和q y ;g 为重力加速度,u 和v 分别为x ㊁y 方向的流速,m /s;F 和G 分别为x ㊁y 方向的通量矢量;S 为源项矢量;z b 为河床底面高程,m;C f 为谢才系数,n 为曼宁系数㊂1.1.2 机器学习快速预报模型预报模型通过采用机器学习算法(KNN 算法等)根据降雨的特征参数,建立预报降雨与情景库降雨的相关关系,根据降雨与内涝结果的对应关系实现城市洪涝快速预报,具体实现流程如图1所示㊂本研究采用的KNN 算法是一种常用的监督学习方法,根据待预测样本与训练集中样本数据的欧式距离获取其最靠近的k 个训练样本,然后通过这k个 邻居”的信息对待预测样本信息进行计算,具体原理如图2所示㊂图1 机器学习快速预报模型流程1.2 系统功能城市洪涝预报以气象预报信息㊁水文及下垫面资料为输入条件,运用数值模型定量模拟预测未来一段时间内的降水产汇流过程㊂城市洪涝预报可在灾害发生前预报出内涝积水点分布,水深变化㊁积水开始及结束时间等重要信息,可为决策人员赢得充足的思考决策时间㊂系统自动预报流程如图3所示㊂图2 K 值临近算法原理图3 城区洪涝自动预报(1)城区预报点位分析基于城区数字高程模型(DEM)数据,由洪涝模型计算推演出区域积水状况,并针对重要预报点位结果进行输出㊂系统在计算结果数据的支撑下,在平台进行可视化展示,提供城区预报点位的详细位置㊁具体坐标㊁积水程度㊁风险等级㊁积水时长,为城区防洪排涝提供支持㊂其中,预报点位风险等级的划分以积水深度为依据㊂(2)城区预报范围分析将城区预报范围内的重点关注对象予以管理,包括但不限于地下空间㊁医院㊁学校㊁加油站㊁高温熔融企业㊁化工医药企业㊁经营带存储企业㊁非煤矿山㊂以洪涝模型计算结果及其输出为数据支撑,前端借助GIS 平台进行预报范围自定义圈画,系统对圈画范围内的详细积水预报情况进行分析和计算,从而输出淹没区域面㊁积水点数量㊁影响房屋户数㊁影响企业数量㊁影响人员数量㊂41史超,王兴桦,吴新垒,等.基于机器学习快速预报模型的城市洪涝预报预警系统研究及应用===============================================(3)积水区域可视化展示积水区域可视化将宁波市海曙区城区范围的地理坐标信息进行网格划分,覆盖全域,在水动力模型积水模拟数据输出的基础上,通过实时渲染和更新,模拟区域内不同降雨程度下的积水退水动态过程㊂在低数据量情况下和大数据量情况下分别以不同的方案进行实现,达到最优展示的目的㊂在万级数据量下,可视化渲染流程如图4所示㊂图4 低数据量积水区域可视化展示方案其中,前端渲染基于WebGL(Web Graphics Li⁃brary)技术,可以将JavaScript 和OpenGL ES 2.0结合在一起,增加OpenGL ES 2.0的一个JavaScript 绑定,实现WebGL 为HTML5Canvas 提供硬件3D 加速渲染,借助系统显卡实现在浏览器中对积涝模拟场景的流畅展示,创建复杂的数据视觉化㊂(4)历史记录管理城区预报输入数据和预报输出数据做持久化处理,借助主流数据库平台MySql 进行大量记录数据的管理,具有记录浏览㊁查询㊁添加㊁删除等功能㊂记录浏览包含所有历史降雨数据的记录,可以查看显示历史预报的预报结果,展示历次积水点位;记录查询可根据记录日期㊁记录预报输入数据㊁记录起止时间范围进行快速定位㊂保存数据同时,也为后期预报精度和预报效率的提升提供数据支持㊂(5)汛情分析简报上传根据城区预报模型结果生成分析报告,支持上传编辑,直观反映海曙区汛情㊂系统界面详见图5㊂图5 城市洪涝快速预报预警平台界面2 应用实证应用区域为宁波海曙城区(见图6),海曙区由于自然地理㊁水文气象㊁河流水系㊁城市建设等原因,易受多种自然灾害威胁㊂特别是 烟花”台风期间,海曙区遭遇 风㊁暴㊁潮㊁洪”四碰头,多处干流㊁平原水位创历史新高,社会经济损失惨重㊂为确保新时代下防汛防台工作顺利开展,亟需建设海曙区 防台减灾智治”应用场景,为防汛人员在应急指挥时进行科学的决策,提供信息支撑,减少灾害损失㊂2.1 预报模型构建模型构建技术路线如图7所示,具体步骤如下:(1)获取研究区域资料,并根据研究区域的DEM 地形数据㊁土地利用情况㊁下垫面下渗条件等数据,构建相应输入文件驱动水动力模型,进而获取足够不同类型降雨条件下内涝灾害模拟数据;(2)进行降雨数据与模拟结果的相关性分析,确定与内涝灾害关系密切的特征参数,并基于py⁃51西北水电㊃2023年㊃第2期===============================================thon 语言,编程构建基于多机器学习算法的城市内涝快速模拟模型;图6 研究区域影像图7 模型构建技术路线(3)运用水动力模型模拟结果进行快速模拟模型训练,考虑运用网格搜索的方法构建参数优化程序,优化模型内部参数选取结构,提升模型模拟性能;(4)研究考虑水动力模型与机器学习造成的误差累积,构建误差校正程序,保障模型的整体模拟稳定性;(5)选择目标城市区域,建立洪城市内涝快速模拟模型,并运用实测及设计降雨资料,进行模型预报性能分析,评估所建快速模拟模型的模拟精度与模拟效率㊂2.1.1 物理模型构建及验证以海曙城区的5m 网格分辨率数字高程模型㊁土地利用数据等为输入数据构建城区洪涝数值模型,网格数量共1892865个㊂选取宁波市2021年台风 烟花”期间连续12h 降雨-内涝事件对物理模型进行验证,12h 累积降雨量为60mm,最大小时降雨量为16.5mm㊂模拟结果发现最大积涝时刻为降雨开始后10h(图),并与实测数据对比,根据6个易涝点位计算得出模型的相对误差为1.8%,故物理模型的精度较高㊂图8 城区淹没峰值水深61史超,王兴桦,吴新垒,等.基于机器学习快速预报模型的城市洪涝预报预警系统研究及应用===============================================表1 易涝点实测与模拟数据监测点名称实测最大水深/m模拟最大水深/m绝对误差/mA西湾路0.610.28-0.333B布政巷0.250.290.040C白云街0.330.400.066D苗圃路0.370.32-0.055E文昌街76号0.320.440.118F秀水社区0.640.660.019 2.2.2 机器学习模型构建及验证根据当地历史降雨资料和暴雨公式构建100种典型降雨模式作为输入,导入物理模型生成对应的内涝情况作为训练集㊂采用KNN算法训练生成城区快速内涝预报模型,并与多场实际降雨内涝结果和物理模型降雨内涝结果进行对比分析,其中 烟花”降雨内涝事件预报结果对比如图9所示,结果表明构建的机器学习模型预报精度较高,能满足本区域的城市内涝快速预报㊂图9 机器学习模型模拟结果与积涝实测对比表2 模型模拟准确度分析位置预报21:00水深/m复盘21:00水深/m实测水深/m预报相对误差复盘相对误差屠园小区0.110.250.25-56%0%广仁街46号0.370.500.323%67%西北街50号0.200.280.25-20%12%文昌街76号0.230.340.3-23%13%尚书街90号0.220.300.3-27%0%布政巷0.160.220.32-50%-31%新芝小区95弄13号0.360.420.6-28%-16%西湾路0.250.390.5-50%-22%咏归路0.190.21---三板桥-菱池街0.200.300.5-60%-40%大沙泥街0.170.260.26-35%0%白云街0.150.190.2-25%-5%前莫家巷0.220.22---新高路4弄0.140.220.2-30%10%体育场路30弄12号0.220.300.3-27%0%东恩中学涵洞0.380.50--56%0%大桥街与苍水街口0.320.370.5-36%-25%苗圃路0.180.210.2-9%3%云石小区0.150.160.2-23%-19%通达路0.180.190.2-8%-5%2.3 平台应用2022年 梅花”台风(见图10)于9月14日20时30分前后以强台风级在浙江舟山普陀沿海登陆,登陆时中心附近最大风力达14级(42m/s),中心最低气压为960hPa,是今年以来登陆我国华东的第一个台风,也是第一个正面袭击宁波市的台风㊂模型综合考虑强降雨㊁上游来洪和潮水的影响,快速预测未来3h和前期累计降雨导致的积水情况(见图11),做到降雨积水信息间隔3h滚动预报,为应急管理和市民出行提供有效信息㊂为综合评价预报模型的准确度,分别将实时预报结果㊁复盘结果与实测水深进行对比分析(见表1),提取翰香小学雨量监测站09:00 21:00(见图12),12h累计降雨量为209mm,最大1h(18:00)降雨量为30mm,最大3h(17:00 19:00)降雨量为78mm,实测18:00 21:00逐小时降雨量为25㊁21㊁11mm㊂根据9月14日21:00监测点数据分别对18:00进行的未来3h预报结果和09:00 21:00复71西北水电㊃2023年㊃第2期===============================================盘结果对模型准确性进行评价㊂预报结果准确率为57%,复盘结果准确率为70%㊂模型在降雨预报准确即复盘情景下,积水深度预测准确率能达到70%,能够为城市洪涝应急管理提供技术支撑㊂下一步将考虑降雨不确定性对城市洪涝预报的影响,提高实时预报的准确率㊂图10 台风 梅花”路径图11 洪涝平台预报预警结果81史超,王兴桦,吴新垒,等.基于机器学习快速预报模型的城市洪涝预报预警系统研究及应用===============================================图12 9:00 21:00逐小时实测降雨量3 结 论采用机器学习方法对大量由物理模型产生的降雨内涝数据进行驱动,构建了城市洪涝快速预报预警平台,并在宁波市海曙城区进行应用,经过对 梅花”台风的暴雨治涝事件的实时预报和复盘模拟,并与实际内涝情况进行对比分析,形成主要结论如下: (1)GAST水动力数值模型可准确模拟城市洪涝过程,为机器学习方法提供大量㊁可靠的训练数据㊂(2)采用的KNN机器学习算法可有效驱动降雨-内涝数据库,准确表征降雨径流过程,对城市内涝情况进行快速准确预报㊂(3)平台以预报模型为核心,采用气象降雨数据作为驱动进行内涝快速预报预警,可为交通出行和应急管理决策提供直观的可视化展示和数据支撑㊂(4)降雨预报的不确定性是导致预报结果偏差的一个重要因素,模型训练库需要持续更新和扩充,提高预报模型精度㊂参考文献:[1] 程晓陶,刘昌军,李昌志,等.变化环境下洪涝风险演变特征与城市韧性提升策略[J].水利学报,2022,53(07):757-768, 778.[2] 陈吉虎.北京市 7㊃21”特大暴雨对城市水土保持和监测工作的启示[J].水土保持通报,2013,33(06):313-316. [3] 陈昆仑,许红梅,李志刚,等.快速城市化下的武汉城市暴雨渍水空间特征及发生机制[J].经济地理,2020,40(05):129-135,19.[4] 靳奎峰,巢惟忐,陈虹杏,等.双偏振天气雷达在灾害性降水估测中的应用 以广州"5.22"特大暴雨灾害为例[J].灾害学,2021,36(03):140-146.[5] 刘家宏,裴羽佳,梅超,刘昌军.郑州 7㊃20”特大暴雨内涝成因及灾害防控[J].郑州大学学报(工学版),2023,44(02):38 -45.[6] 姚冰,余方顺,董其裕.宁波海曙区 烟花”台风复演及对策措施研究[J].水利技术监督,2022(06):66-69,74. [7] 周聂,侯精明,苏锋,等.基于陆气耦合的城市内涝高分辨率模拟预报方法[J].中国给水排水,2021,37(21):131-138. [8] 韩浩,侯精明,金钊.新型流域雨洪过程模拟方法研究[J].西北水电,2022(05):41-46.[9] 马广超.英国排水理念对中国海绵城市建设的启示[J].西北水电,2021(06):34-37.[10] 周思敏,侯精明,高徐军,等.老城区海绵改造降雨致涝及污染物输移过程数值模拟 以西安市小寨老城区为例[J].西北水电,2021(03):11-17.[11] 林珲,吴贤宇,潘家祎,邹海波.中国城市洪涝实时预报研究:现状与挑战[J].测绘学报,2022,51(07):1306-1316. [12] 孙永尚,刘扬,赵军合,王帅帅.北京市内涝积水监测预警系统设计与实现[J].地理空间信息,2015,13(04):117-119. [13] 李海宏,吴吉东,王强,等.基于机器学习方法的上海市暴雨内涝灾情预测模型研究[J].自然灾害学报,2021,30(01):191-200.[14] 侯精明,王润,李国栋,李桂伊.基于动力波法的高效高分辨率城市雨洪过程数值模型[J].水力发电学报,2018,37(03):40-49. [15] 侯精明,李桂伊,李国栋,等.高效高精度水动力模型在洪水演进中的应用研究[J].水力发电学报,2018,37(02):96-107. [16] 陈光照,侯精明,同玉,等.基于大气-水动力集成模型的城市内涝模拟预报研究[J].西安理工大学学报,2020,36(03):330-336.[17] 龚佳辉,侯精明,薛阳,等.城市雨洪过程模拟GPU加速计算效率研究[J].环境工程,2020,38(04):164-169,175. [18] 潘鑫鑫,侯精明,陈光照,等.基于K近邻和水动力模型的城市内涝快速预报[J/OL].水资源保护:1-17[2023-03-20]./kcms/detail/.20220825.1710.002.html.[19] 赵璐,袁立明,郝琨.多示例学习算法综述[J].计算机科学,2022,49(S1):93-99.91西北水电㊃2023年㊃第2期===============================================文章编号:1006 2610(2023)02 0020 04多布水电站鱼道工程设计符合性检验及改进措施研究高 繁1,3,崔 磊2,寇晓梅1,3,牛 乐1(1.中国电建集团西北勘测设计研究院有限公司,西安 710065;2.水电水利规划设计总院,北京 100120;3.陕西省四主体一联合河湖生态系统保护与修复校企联合研究中心,西安 710065)摘 要:鱼道作为主要过鱼设施类型被应用于缓解大坝对鱼类洄游的阻隔影响,近年来国内建成的鱼道工程大多采用竖缝式鱼道㊂由于施工工艺和方法等因素,竖缝式鱼道施工常出现设计符合性偏差问题,影响鱼道工程预期过鱼效果,进行设计符合性检验,分析工程施工与设计不相符导致的功能损失㊂通过以多布水电站鱼道工程为例,开展鱼道工程设计符合性检验,针对鱼道工程存在问题提出改进措施并分析效果,提升了多布水电站鱼道运行效果和管理能力,可为国内同类鱼道建设提供参考㊂关键词:鱼道;设计符合性;改进措施中图分类号:TV135 文献标志码:A DOI :10.3969/j.issn.1006-2610.2023.02.004Study on Design Compliance Inspection and Improvement Measures of Fishway at Duobu Hydropower StationGAO Fan 1,3,CUI Lei 2,KOU Xiaomei 1,3,NIU Le 1(1.Powerchina Northwest Engineering Corporation Limited ,Xian 710065,China ;2.China Renewable Energy Engineering Institute ,Beijing 100120,China ;3.Shaanxi Union Research Center of University and Enterprise for River and Lake Ecosystems Protection and Restoration ,Xian 710065,China )Abstract :As the main type of fish passage facilities ,fishways are used to alleviate the barrier effect of dams on fish migration.Due to factors such as construction technology and methods ,there is deviation in the design compliance of the vertical slot fishway ,which affects the expected passage effect of fishway.However ,carrying out the design compliance inspection can verify the functional loss caused by nonconformity between project construction and design scheme.This study takes the fishway at Duobu Hydropower Station as an example to carry out design compliance inspection.Aiming at the existing problems of the project ,the improvement measures are proposed and the effects after improvement are analyzed ,which improves the operation effect and management ability of the fishway at Duobu Hydropower Station.The study can provide reference for the construction of similar fishways in China.Key words :fishway ;design compliance ;improvement measures 收稿日期:2022-11-28 作者简介:高繁(1993-),男,陕西省西安市人,工程师,主要从事生态环境保护工程设计及科研工作.0 前 言水电工程的建设不仅会阻隔洄游鱼类的上溯通道,而且对半洄游性鱼类和非洄游性鱼类也有很强的阻隔效应,鱼类生境的片段化和破碎化将导致种群遗传多样性降低㊂因此,根据水电工程建设环境保护要求,减缓大坝对鱼类群体产生的阻隔影响,水电站应修建过鱼设施来连通河道生境,为洄游性鱼类留存洄游通道㊂鱼道作为目前中低坝过鱼设施的主要型式被广泛应用㊂鱼道根据型式的不同主要分为技术型鱼道和仿自然通道,其中技术型鱼道由于孔形不同又可分为水池式㊁竖缝式㊁丹尼尔式㊁溢流堰式等㊂因为竖缝式鱼道具备主流明确㊁流场结构02高繁,崔磊,寇晓梅,牛乐.多布水电站鱼道工程设计符合性检验及改进措施研究===============================================。

洪水预警系统迅速响应洪水威胁

洪水预警系统迅速响应洪水威胁

洪水预警系统迅速响应洪水威胁洪水是一种极具破坏性的自然灾害,经常给人们的生命财产安全带来严重威胁。

为了及时预警和减少洪水带来的损失,洪水预警系统应运而生。

本文将介绍洪水预警系统的重要性以及如何迅速响应洪水威胁。

1. 洪水预警系统的重要性洪水预警系统的建立与实施对于保障人们的生命财产安全至关重要。

首先,洪水预警系统能够提前警示人们,使他们有充足的时间采取应对措施。

其次,系统能够提供准确可靠的洪水信息,帮助相关部门制定科学合理的应急计划。

不仅如此,洪水预警系统也对于减少洪水带来的经济损失和人员伤亡起到至关重要的作用。

2. 洪水预警系统的三个主要组成部分洪水预警系统由三个主要组成部分构成,分别是感知系统、传输系统和接收系统。

- 感知系统负责监测洪水可能发生的地区,并收集相关数据。

例如,感知系统可以包括水位传感器和雨量计,用于监测水位和降雨情况。

这些感知设备能够实时采集数据,并将其传输到传输系统。

- 传输系统负责将感知系统收集到的数据传输给接收系统。

传输系统通常采用无线通信技术,能够快速可靠地将数据传输到接收系统。

这样,接收系统能够及时获得关键信息,以便进行分析和响应。

- 接收系统是整个洪水预警系统的核心部分,负责接收和处理传输系统传输过来的数据。

接收系统可以将数据分析处理,并生成预警信息。

该系统可以由专门的人员操作,也可以自动化运行。

3. 洪水预警系统的迅速响应洪水预警系统在接收到数据后,需要能够快速响应洪水威胁。

具体而言,它应该具备以下几个特点:- 迅速性:系统需要能够及时处理数据,并在最短的时间内生成预警信息。

这样,相关部门和个人就能够及时采取应对措施,以减少洪水给人们的生命财产安全带来的威胁。

- 可靠性:洪水预警系统应该提供准确可靠的预警信息,以便人们能够进行有效的应对。

系统需要基于科学和专业的数据分析,以确保预警信息的准确性。

- 实时性:接收系统应该能够实时监测洪水情况,并及时更新预警信息。

只有在获取到及时的数据后,人们才能够做出明智的决策。

基于 ARIMAX 的开封“7·20”特大暴雨城市内涝预报研究

基于 ARIMAX 的开封“7·20”特大暴雨城市内涝预报研究

基于 ARIMAX 的开封“7·20”特大暴雨城市内涝预报研究作者:喻谦花冯峰罗福生霍继超冯跃华张鹏飞来源:《人民黄河》2022年第10期摘要:针对暴雨引发的城市内涝问题,将多元时间序列ARIMAX 模型应用于内涝预报,考虑到积水变化过程是一个时间序列,其主要受到降水过程的影响,用该模型可更好地描述降水量与积水深度之间的关系,更精准地预测城市内涝积水深度及變化过程。

将开封市12个内涝监测站数据进行筛选和分析,对数据进行预处理、时间序列平稳性和白噪声检验,通过绘制降水量与积水深度的自相关系数和偏自相关系数图,判断其是否拖尾,同时对时间序列进行差分操作,选用相应的 ARIMAX 模型。

利用贝叶斯信息准则( BIC)确定模型阶数,确定了12个内涝点对应的 ARIMAX 模型的最终形式。

为了检验模型的精度,将基于12个内涝点建立的ARIMAX 模型进行内涝预测,并计算预测值与实际内涝积水过程的 MAE(平均绝对误差)、MSE(均方误差)、RSME(机器学习误差)。

将开封市“7·20”特大暴雨的降水过程作为输入时间序列,进行内涝积水预测,并与实测值进行对比,认为2 h 以内、50 cm 以下积水深度的预报效果较好。

通过“7·20”特大暴雨的实例验证认为构建的城市内涝 ARIMAX 预报模型适用性较强,进行更多的样本学习训练会取得更精准的预报结果。

关键词:城市内涝;ARIMAX 模型;积水预报;开封市中图分类号:TV122+.1 文献标志码:A doi:10.3969/j.issn.1000-1379.2022.10.008引用格式:喻谦花,冯峰,罗福生,等.基于 ARIMAX 的开封“7·20”特大暴雨城市内涝预报研究[J].人民黄河,2022,44(10):40-45,50.Research on Urban Waterlogging Forecast of“July 20”Rainstorm in Kaifeng City Based on ARIMAXYU Qianhua1,FENG Feng2,3,LUO Fusheng4,HUO Jichao1,FENG Yuehua4,ZHANG Pengfei2,3(1.Kaifeng Key Laboratory of Meteorological Disaster Prevention and Mitigation,Kaifeng 475004,China;2.Yellow River Conservancy Technical Institute ,Kaifeng 475004,China;3.Henan Engineering Technology Center forWater Resources Conservation and Utilization in the Middle and Lower Reaches of the Yellow River,Kaifeng 475004,China;4.Eastern Henan Water Conservancy Administration Bureau ,Kaifeng 475004,China)Abstract:For urban waterlogging caused by heavy rain forecast issue ,multivariate time series ARIMAX model is applied to waterlogging forecast,considering the water change process is time series ,it is mainly influenced by precipitation process ,use the model can better de⁃ scribe the relationship between the precipitation and accumulated water,can be more effective ,more accurate to predict urban waterlogging water depth and the change ing the monitoring data of 12 stations in Kaifeng to screen and analyze ,pre⁃process the data,test the stationarity of time series and white noise ,through mapping the autocorrelation coefficient of precipitation and water depth and partial autocor⁃relation coefficient diagram to determine its tail and truncation the difference for the time series of operation to choose the correspondingARIMAXmodel.Bayesian information criterion ( BIC ) was used to determine the model order,and the final form of ARIMAX model corre⁃sponding to 12 waterlogging stations was determined.In order to test the accuracy of the model the ARIMAX model established by 12 water⁃logging stations was used to predict waterlogging,and MAE ( mean absolute error),MSE ( mean square error) and RSME ( machine learning error) between the predicted value and the actual waterlogging process were calculated.Taking the precipitation process of the“ July 20”heavy rain in Kaifeng as the input time series the prediction of water logging was carried out and compared with the measured value.Itwasbelieved that the prediction effect of water depth within 3 hours and below 50 cm was better.Through the example of the“ July 20”rainstormit is verified that the constructed ARIMAX forecast model of urban waterlogging has strong applicability.More sample learning and trainingwill obtain more accurate forecast results providing data basis and scientific basis for the development of waterlogging warning plan andemer⁃gency and disaster reduction plan for the city in response to the rainstorm.Key words:urban waterlogging;ARIMAX model;waterlogging forecast;Kaifeng City2021年7月20日,河南省遭遇極端强降雨。

某洪水气象预报系统设计及实现研究

某洪水气象预报系统设计及实现研究

某洪水气象预报系统设计及实现研究洪水是常见的自然灾害之一,针对其发生,气象预报可起到极大的作用。

本文将介绍某洪水气象预报系统的设计和实现研究。

一、需求分析为有效应对洪水灾害,需要建立一个洪水气象预报系统,能够对未来洪水灾害发生的情况进行预测和预警。

通过预测和预警,我们可以及时采取措施,防止洪水灾害的发生和蔓延,最大程度地保护人民的生命、财产和生态环境。

二、系统设计1. 数据源洪水的发生和蔓延涉及多种因素,如降雨量、水位、风速、气压等。

因此,需要通过气象观测站、卫星遥感等多种数据源获得相关数据,以确保预测的准确性。

2. 数据处理为了更好地获得连续准确的预测结果,需要将多个数据源进行融合处理。

在这个过程中,需要考虑多种算法,并进行有效的评估和比较。

3. 模型选择根据前期的数据处理结果进行模型选择,以达到最好的预测结果。

通常,常用的模型包括FloodR、HBV、HEC-HMS、SWMM 等。

4. 预测和预警通过选定的模型,对全国各地的洪水灾害情况进行预测和预警。

预测结果将会包括降雨量、水位、洪水范围和灾害程度等信息。

5. 可视化与应用在预测和预警完成后,需要将结果进行可视化,并将结果提供给相关的单位和个人使用。

同时,在开发洪水气象预测系统时,需要考虑手机APP、微信公众号等工具,以便用户能够充分了解相关预测信息。

三、系统实现1. 基础技术系统实现需要使用相关的技术,包括Python、Matlab、Java、R 等。

为加快系统的开发和降低开发难度,可以使用基于云计算的开发平台,如Google Cloud Platform、Azure等。

2. 算法开发与优化选定算法后,需要进行相关的算法开发与优化。

在开发过程中,需要考虑到算法的稳定性和可扩展性,以达到更好的效果。

3. 模型建立根据前期的数据处理结果及算法开发与优化结果,进行模型的建立。

在建模过程中,需要不断地调整参数和算法,以获得更好的预测结果。

4. 评估和比较建立完毕的模型需要进行评估和比较。

科技成果——基于网格化的大数据洪水预报系统

科技成果——基于网格化的大数据洪水预报系统

科技成果——基于网格化的大数据洪水预报系统
技术开发单位
北京艾力泰尔信息技术股份有限公司
成果简介
基于网格化的大数据洪水预报系统采用适合区域水文气象特征的分布式水文模型,实现流域断面水位、流量预报。

并发展了大数据驱动的产汇流模拟理论,创新了产汇流参数、地理规律分析方法及创新了大数据支撑下新安江模型系列,增强了其理论基础,建立了大数据驱动的智能型实时预报水文模型系列、创新了大数据驱动的洪水预报和实时校正方法,开发出云计算洪水预报平台。

主要性能指标
(1)10秒内实现地面站数据的处理,得到小流域的相关数据参数。

(2)实现以10分钟(30分钟/60分钟)为降雨径流过程模拟时间步长的实时连续计算。

(3)平移和缩放速度小于0.3秒,多图形包含、重叠、相交计算小于3秒,数据名称查询、定位查询速度小于1秒,客户端响应时间小于3秒。

适用范围
该系统应用在对大江大河、小流域等洪水预报中,在水利行业、农业方面广泛应用。

暴雨洪水管理模型-概述说明以及解释

暴雨洪水管理模型-概述说明以及解释

暴雨洪水管理模型-概述说明以及解释1.引言1.1 概述概述暴雨洪水是一种常见的自然灾害,其突发性和破坏力对人类的生命和财产安全构成了严重威胁。

随着气候变化和城市化进程的加速,暴雨洪水事件的频率和强度也呈现出逐渐增加的趋势。

因此,如何有效地管理和应对暴雨洪水成为了亟待解决的问题。

本文将探讨暴雨洪水管理模型,旨在通过科学合理的方法和手段来提高对暴雨洪水的应对能力和管理水平。

该模型将涵盖暴雨洪水的危害、现有的管理方法以及需要改进的地方,以期为相关领域的研究人员和政策制定者提供参考和借鉴。

暴雨洪水的危害主要表现在以下几个方面:一是对人类生命安全的威胁,包括溺水、水质污染和疾病传播等;二是对城市基础设施和建筑物的破坏,如道路、桥梁、房屋等;三是对农田和农作物的影响,导致农业生产的损失;四是对环境的破坏,包括土壤侵蚀、水资源浪费等。

针对这些危害,目前存在一些管理方法,如建设防洪堤坝、整治河道、提高排水能力等。

然而,面对不断变化的暴雨洪水形势,现有方法仍然存在一些问题和不足之处,需要进一步改进和完善。

因此,本文将重点关注暴雨洪水管理模型的构建和应用。

通过分析暴雨洪水的形成机制和发展规律,结合相关的地理、气象和水文数据,建立起一套科学、精确的模型,以指导暴雨洪水的应对和管理工作。

该模型将基于先进的技术手段,如遥感、地理信息系统和数值模拟等,为预测洪水的发生和演变提供准确的数据支持。

同时,该模型还将考虑不同地区和城市的特点,制定相应的管理策略和措施,以实现对暴雨洪水的有效管理和防控。

综上所述,暴雨洪水管理模型的建立具有重要的理论和实践意义。

通过深入研究暴雨洪水的危害、现有的管理方法以及需要改进的地方,可以为相关领域的决策者提供科学依据,推动暴雨洪水管理工作的持续发展。

与此同时,需要进一步加强科学研究,不断改进和完善暴雨洪水管理模型,以应对不断变化的自然环境和社会需求,为构建安全、可持续的社会发展环境作出贡献。

1.2 文章结构文章结构部分的内容可以依照以下方式编写:文章结构:本文主要包括引言、正文和结论三个部分。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

赫尔辛基城市暴雨洪水预报系统介绍 摘要:城市暴雨洪水预警预报系统以暴雨洪水管理模型为基础,结合雷达观测技术以及水文测报技术,对赫尔辛基城市暴雨洪水进行分析预测,对预测得到的数据进行处理、展示,当某区域积水水位值超过预先设定的最高水位值时,触发报警系统,及时提醒相关部门采取应对措施,从而减小城市暴雨洪水对人民生命财产造成的损害。

关键词:暴雨洪水;预警预报;管理模型;水文测报 芬兰位于欧洲北部,被誉为“千湖之国”,水域占区域面积百分比如图1所示,该地区平均降雨量为500-800mm(其中30%-50%为降雪)。近年来,社会高速发展,带来社会经济效益增长的同时也对居住环境造成巨大压力,温室效应、热岛效应等引发的全球气候变暖,造成了诸如冰川融化、强降雨等一系列环境问题。受全球气候变化的影响,赫尔辛基城市降雨主要集中在夏季,随着全球气候变暖,夏季持续高温带来的强降水,引起城市径流量瞬间增大;同时由于城市化进程引起下垫面环境的变化,当产生的径流流量超过城市管网的排泄能力时,雨水在城市低洼处堆积而产生的城市内涝,严重影响城市居民的正常生活,并且暴雨洪水产生的雨洪污水汇入城市排水管网,最终排入城市河道,将改变水体环境,危及水体中生活的动植物。因此,采用先进的雨洪管理理念,结合现代水文自动测报技术以及雷达观测技术,采用暴雨洪水管理模型,开展洪水预警预报技术研究,对于实时监测赫尔辛基城市水文信息,把握城市暴雨洪水发展动向,及时采取应对措施减少生命财产损失具有重要意义。

一、城市暴雨洪水预警预报系统介绍 城市暴雨洪水预警预报系统(系统原理图如图2所示),以暴雨洪水管理模型为基础,结合水文自动测报技术、卫星雷达观测技术、通信技术、数据库技术、计算机技术等,对雨量数据进行采集、存储、传输、处理、展示以及预警预报。系统采用雨量站及雷达两种方式观测比对获得雨量数据,将该数据输入到洪水预报管理模型中,经处理后输出雨量数据预测值,同时计算出城市径流情况,用于手机端及电脑端洪水预警预报信息展示,各子模块介绍如下。

图1 水域占区域面积百分比 图2 城市暴雨预警预报系统工作原理图

1.系统网络架构总体部署 系统选用B/S(Browser/Server)架构,通过一定的权限控制实现多客户访问,交互性强,版本更新方便。客户端主要负责水文数据查询展示以及预警预报信息展示,服务器端主要负责数据处理、与数据库以及其他子系统交互。联网状态下,客户端和服务器端进行通信,实现数据处理、存储、展示与推送等过程。

2.雨量观测 暴雨监测是暴雨洪水预报的基础,过去采用遥测雨量站对暴雨数据进行采集,该方法存在站点布设密度低,站点分布不均匀,数据缺乏时空连续性且精度不高等局限性,特别是持续强降水情形下,数据错报、漏报等都会给当地环境部、农林部以及城市居民使用、查看过程中造成极大不便。雷达测报技术,对区域现有雨量站测得的数据进行校准比较,弥补水文测报系统的短板,从而提高数据的质量。该技术通过探测雨滴形状推定雨量大小,实现对地区强降雨信号的捕捉、监控、预报、传输,提高数据时空连续性以及提供实时的更高精度的雨量数据。数据采集间隔可设为5min、10min、1h等,获取得到雨量数据后,作为输入值输入到暴雨洪水管理模型,从而得出该地区雨量数据某时段内的预测值。

3.暴雨洪水管理模型 暴雨洪水管理模型(Storm Water Management Model,SWMM),1969至1971年间由美国环境保护署(Unite States Environment Protection Agency,EPA)提出,旨在解决日益严峻的城市暴雨洪水问题。该模型通过建立动态的降雨-径流模型,对城市的排水管网进行长期或者短期的跟踪、模拟,从而计算分析出任意时间点城市管网中任意径流的流量及水质情况,用于指导优化现有城市排水管网。当暴雨洪水来临时,对暴雨洪水进行优化调度,降低城市洪水发生的概率。选取赫尔辛基最近五年海量的雨量、径流数据,采用数据挖掘技术,对该地区数据进行选取、转换、清洗等操作,从而挖掘出有效的水文数据。同时选取M5模型树算法,输入挖掘得到的水文数据,不断对模型进行训练学习,得出适用于本系统的模型。该模型树将输入属性空间划分成互相平行的矩形区域,对每一个分区域,得出对应的回归模型,实现雨量的模拟及预测,M5模型树算法流程图如图3所示。通过M5模型树算法建立的洪水预报模型,操作方便、维护简单,并且能够提供高精度预测数据。

图3 M5模型树算法流程图 4.数据库设计 根据系统功能的需求,本系统的设计既支持测站数据、雨量数据、预警预报数据、用户数据等历史数据的查询,也支持实时数据的查询、存储与处理,历史数据通过查阅当地环境部以及农林部历史资料可得,实时数据可通过实时读取雨量站数据接口得到。如1.2节中介绍,通过暴雨洪水管理模型,生成城市暴雨洪水预测值,用以进行数据展示以及预警预报。本系统采用主流数据库SQL Server。系统部分E-R设计图如图4。其中,StationData对应测站数据表,包含站点精度、站点ID号等信息;RainFall对应降雨信息表,包含时间、测站降雨量等信息;User对应用户数据表,包含用户ID号、用户类型ID号等信息; WarningInformation对应预警信息表,包含预警信号ID、预警解除时间等信息。 图4 系统E-R图 对数据库表进行设计,表设计分别如下:(1)测站数据表。此表用来存放观测站的基本信息,如站点ID号、站点编号、经度、纬度、高程等(表1)。

表1 StationData表结构

字段名称 字段说明 字段类型 空/非空 是否为主键

StationID 站点ID号 Int NOTNULL 是 StationNumber 站点编号 Int NULL 否

StationName 站点名称 Nvarchar(50) NULL 否

Longtitude 经度 Float NULL 否 Latitude 纬度 Float NULL 否 Altitude 高程 Float NULL 否 … (2)降雨信息表。此表用来记录各个站点的逐日降雨量,表中包括时间、各个站点的雨量值等信息(表2)。

表2 RainFall表结构

字段名称 字段说明 字段类型 空/非空 是否为主键

Time 时间 Datatime NOTNULL 是 R1 站点1雨量 Float NOTNULL 否 R2 站点2雨量 Float NOTNULL 否 … … … … … (3)预警信息表。此表用于存放预警的信息,如预警信号ID、预警信号名称、预警信号类型、预警级别、预警方式、是否应用、预警word文档地址、是否解除预警及解除时间等(表3)。

表3 WarningInformation表结构

字段名称 字段说明 字段类型 空/非空 是否为主键

WarningInformationID 预警信号ID Int NOTNULL 是

WarningInformationTitle 预警信号名称 Nvarchar(200) NOTNULL 否

WarningInformationType 预警信号类型 Int NOTNULL 否

WarningEventLevel 预Int NOT否 警级别 NULL WarningMethod 预警方式 Int NOTNULL 否 isValid 是否应用 Int NOTNULL 否

WordURL 预警word文档地址 Nvarchar(60) NOTNULL 否 IsRelieved 是否解除 Int NULL 否 RelieveTime 解除时间 Datatime NULL 否 (4)用户信息表。此表用来存放用户的具体信息,如用户ID、登录账号、登录密码、账号类别、用户类型ID等(表4)。

表4 User表结构

字段名字段说字段类型 空/非是否为称 明 空 主键 UserID 用户ID Int NOTNULL 是 LoginName 登录账号 Nvarchar(30) NOTNULL 否

Password 登录密码 Nvarchar(30) NOTNULL 否

Type 账号类别 Nvarchar(30) NOTNULL 否 UserTypeID 用户类型ID Int NOTNULL 否

5.信息查询及预警模块介绍 信息检索查询模块用于收集水情信息,当用户需要查询雨量信息时,在数据查询区选择需要查询的测站编号以及查询时间间隔,系统即将实时查询到的结果显示在降水量预报记录框中,即可查看雨量信息。预警模块用于发布和搜索预报信息,按照设定的扫描时间间隔,对雨量站雨量数据进行扫描,将扫描得到的数据与提前录入模块中的雨量阈值进行比对,超过阈值时,触发预报语音、手机短信等推送服务,将预警预报信息推送给相应个人以及部门,同时记录下扫描结果、预警预报记录信息生成.txt文件,用于后期查阅。 二、赫尔辛基暴雨洪水模拟展示及应用 系统选定赫尔辛基地区作为模拟分析对象,结合ArcGIS技术,调用Google地图数据,确定赫尔辛基地区物理结构网络图,将该地区雨量计站点信息、雨量数据、子流域、节点、连接点、时间步长等信息输入到系统中(考虑到赫尔辛基地区冬季下雪,还需要录入该地区的气象数据,模拟融雪过程)。当暴雨洪水来临时,暴雨洪水管理模型对分析出的降雨量进行计算,结合该区域结构信息、水文数据、气象数据等进行径流模拟,分析出该地区未来特定时段的降雨-径流关系,得出积水情况,进而实现洪水模拟、预警预报。赫尔辛基城市暴雨洪水展示图如图4所示(该图片来源于芬兰VTT公司)。

图5 赫尔辛基城市暴雨洪水展示图 以赫尔辛基2020年雨量信息为例,该区域雨量-时间关系图如下(数据来源于芬兰环境部)。由图5可以直观看出,赫尔辛基夏季雨水充沛,城市雨量主要集中在6-8月,日均降雨量为15mm/d左右。其中,实线为该地区雨量数据实测结果,虚线为该地区预测结果,采用相对误差公式:

经过分析计算,得到每个月对应的相对误差,求平均值可得该地区2020年度洪水预报相对误差为18.02%,在洪水预报误差范围内,基本符合要求。

相关文档
最新文档