城市暴雨强度计算方法及系统、设备、存储介质转让专利
申请号 : CN202111381892.0
文献号 : CN113821939B
文献日 : 2022-03-15
发明人 : 包正铎 , 贺卫宁 , 吴羽璇 , 陈蕃 , 邓征宇
申请人 : 中机国际工程设计研究院有限责任公司
摘要 :
权利要求 :
1.一种基于时空分布特征的城市暴雨强度计算方法,其特征在于,包括以下内容:构建区域气象模型并设置目标城市的模拟区域和模拟网格单元大小,并获得适用于目标城市降雨事件模拟的最优模型配置方案;
利用获得的最优模型配置方案配置区域气象模型,采用配置好的区域气象模型执行历史降雨模拟和/或未来降雨模拟,分别获得各个模拟网格单元的历史降雨时序数据和/或未来降雨时序数据;
基于各个模拟网格单元的历史降雨时序数据计算目标城市网格化的历史暴雨强度,和/或,基于各个模拟网格单元的未来降雨时序数据计算目标城市网格化的未来暴雨强度;
所述构建区域气象模型并设置目标城市的模拟区域和模拟网格单元大小,并获得适用于目标城市降雨事件模拟的最优模型配置方案的过程具体包括以下内容:采用区域气象模型WRF模拟城市降雨动态演进和网格化降雨时空分布,并根据实际需要设置模拟区域和模拟网格单元大小;
获取目标城市的历史暴雨记录,筛选出多场暴雨及以上降雨等级的强降雨事件;
设计不同模型配置方案逐一模拟各场强降雨事件演化过程,获得逐场强降雨事件的模拟降雨数据;
将各场强降雨事件的历史降雨数据与对应的模拟降雨数据进行比较,统计分析两者之间的相关性和模拟偏差,并根据分析结果综合评价不同模型配置方案的模拟效果,筛选出模拟效果最优的模型配置方案。
2.如权利要求1所述的基于时空分布特征的城市暴雨强度计算方法,其特征在于,所述采用配置好的区域气象模型执行历史降雨模拟,获得各个模拟网格单元的历史降雨时序数据的过程具体包括以下内容:
收集目标城市连续多年的年最大日降雨量的发生日期,并按照预设顺序进行排序;
采用配置好的区域气象模型按照顺序逐一进行年最大日降雨量事件的网格化模拟,模拟结果记录时间间隔为1小时,提取出每个模拟网格单元的历史模拟降雨时序数据,重复上述内容,完成所有年最大日降雨量的网格化模拟;
针对每个年最大日降雨量事件,基于每个模拟网格单元的历史模拟降雨时序数据统计不同降雨历时的累计降雨量,并计算得到各个模拟网格单元对应不同降雨历时的历史最大降雨量,以构建各降雨历时的历史年最大降雨量空间矩阵。
3.如权利要求1所述的基于时空分布特征的城市暴雨强度计算方法,其特征在于,所述采用配置好的区域气象模型执行未来降雨模拟,获得各个模拟网格单元的未来降雨时序数据的过程具体包括以下内容:
获取全球气候模型的气候预测结果数据;
将全球气候模型的气候预测结果数据作为区域气象模型的边界条件和初始条件输入,进行连续的未来降雨模拟;
提取目标时段内各个模拟网格单元的未来模拟降雨时序数据,统计不同降雨历时的累计降雨量,并计算得到各个模拟网格单元对应不同降雨历时的未来最大降雨量,以构建各降雨历时的未来年最大降雨量空间矩阵。
4.如权利要求2所述的基于时空分布特征的城市暴雨强度计算方法,其特征在于,所述基于各个模拟网格单元的历史降雨时序数据计算目标城市网格化的历史暴雨强度具体包括以下内容:
基于各降雨历时的历史年最大降雨量空间矩阵,采用广义极值分布函数对模拟网格单元的各降雨历时对应降雨量极值进行分布拟合;
拟合得到广义极值分布函数对应各降雨历时的位置参数空间矩阵、尺度参数空间矩阵和形状参数空间矩阵;
采用基于广义极值分布函数推求的暴雨强度公式计算各个模拟网格单元的各重现期、各降雨历时对应的设计暴雨强度,以构建各重现期、各降雨历时对应的设计暴雨强度空间矩阵;
将各重现期、各降雨历时对应的设计暴雨强度空间矩阵导入ArcGIS系统,根据目标城市市域shp文件,提取目标城市市域的暴雨强度空间数据。
5.如权利要求1所述的基于时空分布特征的城市暴雨强度计算方法,其特征在于,模拟区域采用单层模拟网格单元或者多层嵌套模拟网格单元设置,目标层模拟网格单元的水平分辨率设置为1km 4km。
~
6.如权利要求1所述的基于时空分布特征的城市暴雨强度计算方法,其特征在于,还包括以下内容:
对历史暴雨强度和/或未来暴雨强度的计算结果进行可视化展示。
7.一种基于时空分布特征的城市暴雨强度计算系统,采用如权利要求1 6任一项所述~
的方法,其特征在于,包括:
模型构建模块,用于构建区域气象模型并设置目标城市的模拟区域和模拟网格单元大小,并获得适用于目标城市降雨事件模拟的最优模型配置方案;
网格化降雨模拟模块,用于利用获得的最优模型配置方案配置区域气象模型,采用配置好的区域气象模型执行历史降雨模拟和/或未来降雨模拟,分别获得各个模拟网格单元的历史降雨时序数据和/或未来降雨时序数据;
暴雨强度计算模块,用于基于各个模拟网格单元的历史降雨时序数据计算目标城市网格化的历史暴雨强度,和/或,基于各个模拟网格单元的未来降雨时序数据计算目标城市网格化的未来暴雨强度。
8.一种设备,其特征在于,包括处理器和存储器,所述存储器中存储有计算机程序,所述处理器通过调用所述存储器中存储的所述计算机程序,用于执行如权利要求1 6任一项~
所述的方法的步骤。
9.一种计算机可读取的存储介质,用于存储基于时空分布特征计算城市暴雨强度的计算机程序,其特征在于,所述计算机程序在计算机上运行时执行如权利要求1 6任一项所述~
的方法的步骤。
说明书 :
城市暴雨强度计算方法及系统、设备、存储介质
技术领域
背景技术
到城市排水防涝基础设施建设和海绵城市建设的规划和工程项目的设计建设,进而影响到
整个排水工程建设投资规模,以及系统的效能发挥和维护成本。我国暴雨强度公式编制经
过数十年的实践发展,暴雨强度公式的精度已经有较大提高,随着新理论、新技术和新方法
的提出和应用,暴雨强度公式的编制体系仍在完善。
大气环流、海陆因素以及气候变化等大尺度的环境影响,还受到地形地貌、城市化造成的
“热岛”效应、“雨岛”效应等多种区域因素的共同作用,而且在城市尺度上,降雨量的时空分
布也往往存在显著差异。然而,我国大部分城市的气象监测站点数量较少、建成时间短,降
雨数据记录存在时间间隔长、质量参差不齐以及数据缺失等问题,极大影响了传统暴雨强
度公式编制方法的适用性。目前,绝大部分城市仅编制有基于单个降雨监测站点的城市暴
雨强度公式,尚未在城市暴雨强度公式编制工作中考虑城市尺度上的降雨空间分布差异
性,仅北京、厦门等地根据城市暴雨特性进行了暴雨分区划分,但其仍以城市内的多个气象
站点降雨记录为基础,降雨监测站点数量和分布,降雨数据历史记录完整性、合理性和代表
性等因素对结果的影响较大;同时,在暴雨分区划分中,未考虑地形、下垫面等因素的影响,
极易造成在局部地区失真。此外,以暴雨强度公式为依据设计的城市排水防涝系统、海绵城
市系统在建成后通常需要长期服务,因此,在城市排水防涝系统、海绵城市系统规划设计中
应当包含未来的城市降雨特征,然而,传统暴雨强度公式编制均以历史降雨记录为基础开
展,仅能体现当前城市降雨特征,而缺乏对气候变化背景下城市未来降雨特征变化的预测。
发明内容
或未来降雨时序数据;
度。
选出模拟效果最优的模型配置方案。
重复上述内容,完成所有年最大日降雨量的网格化模拟;
最大降雨量,以构建各降雨历时的历史年最大降雨量空间矩阵。
建各降雨历时的未来年最大降雨量空间矩阵。
空间矩阵;
~
单元的历史降雨时序数据和/或未来降雨时序数据;
市网格化的未来暴雨强度。
法的步骤。
的步骤。
实际需求设计目标城市的模拟区域和模拟网格单元大小,同时,区域气象模型耦合了城市
地形、下垫面、植被等多种环境因素的共同作用,能够很好地模拟城市的历史降雨过程,特
别适用于丘陵、高原等地形复杂地区城市的暴雨强度计算,而传统采用多雨量计进行城市
尺度空间插值,无法考虑地形、下垫面分布、城市热岛效应等多种因素的影响。并且,还获得
了适用于目标城市降雨事件模拟的最优模型配置方案,保证了降雨模拟结果的准确度。然
后,利用配置好的区域气象模型执行历史降雨模拟和/或未来降雨模拟,以分别获得各个模
拟网格单元的历史降雨时序数据和/或未来降雨时序数据,最后基于模拟得到的历史降雨
时序数据和/或未来降雨时序数据分别计算历史暴雨强度和/或未来暴雨强度,基于区域气
象模型的模拟数据来进行暴雨强度计算,不仅可以很好地适用于城市气象监测站点数据
少、建站时间短、降雨记录中断/缺失、降雨记录尚不足以支撑暴雨强度分析的地区,而且能
够应用于预测气候变化背景下未来不同时期城市设计暴雨强度的变化过程,为城市涉水系
统的设计、城市防洪排涝、海绵城市建设等提供了精细化设计、评估的工具。
附图说明
示意图。
图。
意图。
具体实施方式
序数据和/或未来降雨时序数据;
暴雨强度。
而且可以根据实际需求设计目标城市的模拟区域和模拟网格单元大小,同时,区域气象模
型耦合了城市地形、下垫面、植被等多种环境因素的共同作用,能够很好地模拟城市的历史
降雨过程,特别适用于丘陵、高原等地形复杂地区城市的暴雨强度计算,而传统采用多雨量
计进行城市尺度空间插值,无法考虑地形、下垫面分布、城市热岛效应等多种因素的影响。
并且,还获得了适用于目标城市降雨事件模拟的最优模型配置方案,保证了降雨模拟结果
的准确度。然后,利用配置好的区域气象模型执行历史降雨模拟和/或未来降雨模拟,以分
别获得各个模拟网格单元的历史降雨时序数据和/或未来降雨时序数据,最后基于模拟得
到的历史降雨时序数据和/或未来降雨时序数据分别计算历史暴雨强度和/或未来暴雨强
度,基于区域气象模型的模拟数据来进行暴雨强度计算,不仅可以很好地适用于城市气象
监测站点数据少、建站时间短、降雨记录中断/缺失、降雨记录尚不足以支撑暴雨强度分析
的地区,而且能够应用于预测气候变化背景下未来不同时期城市设计暴雨强度的变化过
程,为城市涉水系统的设计、城市防洪排涝、海绵城市建设等提供了精细化设计、评估的工
具。
效果,筛选出模拟效果最优的模型配置方案。
~
平分辨率的大气科学研究和数值预报,已广泛地应用于云尺度到天气尺度等不同尺度的重
要天气特征模拟和预报,本发明采用WRF模型来模拟城市降雨动态演进和网格化降雨时空
分布,可以根据实际需求设计目标城市的模拟区域和网格单元大小。可选地,模拟区域采用
单层模拟网格单元或者多层嵌套模拟网格单元设置,目标层模拟网格单元的水平分辨率设
置为1km 4km,其中,目标层模拟网格单元指的是单层模拟网格单元,或者多层嵌套中最底
~
层的模拟网格单元,后续说明中的模拟网格单元均指的是目标层模拟网格单元。然后,根据
目标城市的历史暴雨记录,筛选出多场暴雨及以上降雨等级的强降雨事件,并获取其对应
的历史降雨记录。再参考WRF模型用户手册及文献报道的区域暴雨事件模拟的模型物理机
理参数化方案的推荐组合方案,建立一系列可用于区域暴雨事件模拟的模型物理机理参数
化方案,包括微气象方案、积云对流方案、长波辐射方案、短波辐射方案、边界层方案和陆面
过程方案,构建方案库,选择不同的组合配置方案,对筛选出的多场强降雨事件进行逐一模
拟,获得每场强降雨事件对应的模拟降雨数据。最后,将各场强降雨事件的历史降雨数据与
对应的模拟降雨数据进行比较分析,例如分析历史降雨数据与模拟降雨数据之间的相关性
2
(R)和均方根误差(RMSE),当相关性最高、均方根误差最小时,意味着模拟效果最佳,其中,
2
相关性(R)和均方根误差(RMSE)的计算公式均为公知技术,故在此不再赘述。基于历史降
雨数据与模拟降雨数据的比较分析结果综合评价区域气象模型在各个组合配置方案下的
模拟效果,筛选出模拟效果最佳的组合配置方案,将其作为最佳模型配置方案。
序数据,重复上述内容,完成所有年最大日降雨量的网格化模拟;
历时的历史最大降雨量,以构建各降雨历时的历史年最大降雨量空间矩阵。
理,核验无误后,筛选出逐年年最大日降雨量的发生日期,并按照时间顺序进行排序,并编
号。然后,选择编号1的年最大日降雨量事件,利用配置好的区域气象模型进行网格化模拟,
模拟结果记录时间间隔为1小时,然后提取出每个模拟网格单元的历史模拟降雨时序数据;
再选择编号2的年最大日降雨量事件,重复上述内容,不断重复,直至完成所有年最大日降
雨量的网格化模拟,并分别提取出对应的每个模拟网格单元的历史模拟降雨时序数据。然
后,针对每个年最大日降雨量事件,基于每个模拟网格单元的历史模拟降雨时序数据可以
统计得到不同降雨历时的累计降雨量,并计算得到各个模拟网格单元对应不同降雨历时的
历史年最大降雨量,例如2019年的年最大降雨量事件发生日期为5月13日,降雨历时为3小
时所对应的最大降雨量为该日00时 24时逐3小时降雨量累加,筛选累加值最大的值得到。
~
从而,构建出各降雨历时的历史年最大降雨量空间矩阵 ,其中,
表示网格坐标(m,n)在降雨历时为t时的历史最大降雨量,m和n为网格坐标的最大值,t表示
降雨历时。
对应的历史最大降雨量进行对比分析,分别计算每场年最大日降雨量事件对应的模拟降雨
校正系数,具体可以采用线性拟合的方式得到校正系数,或者直接采用步骤S1中计算得到
的相关性系数作为校正系数。再采用每场年最大日降雨量事件对应的模拟降雨校正系数对
步骤S23a中得到的各降雨历时的历史年最大降雨量空间矩阵中最大降雨量进行修正,进一
步确保了模型模拟结果的准确度。
雨量,以构建各降雨历时的未来年最大降雨量空间矩阵。
模型WRF的边界条件和初始条件输入,进行连续的未来降雨模拟。然后,根据实际需要提取
出目标时段内各个模拟网格单元的未来模拟降雨时序数据,与历史模拟降雨时序数据相同
的,统计不同降雨历时的累计降雨量,并计算得到各个模拟网格单元对应不同降雨历时的
未来最大降雨量,同样构建各降雨历时的 未来年最大降雨量空间矩阵
,其中, 表示网格(m,n)在降雨历时为t时的未来最大降
雨量,m和n为网格坐标的最大值,t表示降雨历时。
暴雨强度空间矩阵;
历时,可以获取N个年最大日降雨量事件对应的历史年最大降雨量空间矩阵,共N个矩阵,采
用GEV函数对各个模拟网格单元的降雨量极值进行分布拟合,具体的概率密度分布函数表
示为:
对应的GEV函数的位置参数、尺度参数和形状参数。
矩阵,具体为: 。
据。
各降雨历时的历史年最大降雨量空间矩阵替换为各降雨历时的未来年最大降雨量空间矩
阵即可,故具体过程在此不再赘述,参考上述内容即可。
分布。以长沙市城区作为模型模拟中心,在WRF模型中配置模型框架,设置模型模式为
“ARW”,嵌套层数为3层,分别为D01、D02和D03,平面网格分辨率比例为1:3:3,其中D01、D02
和D03的网格数分别为6708、4752和3240,地理数据分辨率分别设置为gtopo_5m+modis_
30s_lakes+5m+default、gtopo_2m+modis_30s_lakes+2m+default和gmted2010_30s+
modis_30s_lakes+30s+default,模型地图投影方式为lambert,中心经纬度为28.2和
113.0,truelat1=30,turelat2=60,stand_lon=113。
括微气象方案、积云对流方案、长波辐射方案、短波辐射方案、边界层方案和陆面过程方案,
构建方案库,具体如表2所示:
~
暴雨事件模拟期(约7 10天)及暴雨事件后模拟期(约1 3天),共计约13 20天,其中的年最
~ ~ ~
大降雨日对应为暴雨事件模拟期最后一天,如遇到暴雨事件时长超过预设时长时,需要根
据实际情况调整模拟期。
拟降雨量和实测记录降雨量的相关性(R)和均方根误差(RMSE),将其作为指标评估不同组
合方案对目标城市不同暴雨事件降雨量的模拟能力,筛选出最适用于长沙暴雨事件模拟的
方案组合,具体为如表4所示:
出各年度的年最大降雨日,根据日降雨量从大到小对降雨事件进行排序并记录对应的日降
雨量,具体如表5所示。
雨日事件数据集,根据降雨事件编号依次进行模拟,提取出D03网格年最大降雨日对应的降
雨值。根据长沙站记录的40场次最大日降雨事件的降雨数据和对应网格单元的模拟降雨
量,计算逐场降雨事件的模拟降雨校正系数。通过D03网格中各个模拟网格单元的模拟降雨
量与校正系数乘积得到校正后的模拟网格单元降雨量。
函数推求得到的设计暴雨强度公式计算D03网格24小时降雨历时对应的不同重现期暴雨强
度。在ArcGIS系统中,导入长沙市域shp文件,使用extract by mask命令,提取长沙市市域
内各个模拟网格单元的暴雨强度分布。计算结果如图7和图8所示,图7为2019年长沙市市域
24小时降雨历时对应的50年一遇的分钟暴雨强度分布图,图8为2019年长沙市市域24小时
降雨历时对应的50年一遇暴雨强度的空间分布差异图。
结果数据的2040年、2060年长沙市市域网格化各重现期24小时降雨历时的设计暴雨强度空
间分布图,具体结果如图9和图10所示。
单元的历史降雨时序数据和/或未来降雨时序数据;
市网格化的未来暴雨强度。
而且可以根据实际需求设计目标城市的模拟区域和模拟网格单元大小,同时,区域气象模
型耦合了城市地形、下垫面、植被等多种环境因素的共同作用,能够很好地模拟城市的历史
降雨过程,特别适用于丘陵、高原等地形复杂地区城市的暴雨强度计算,而传统采用多雨量
计进行城市尺度空间插值,无法考虑地形、下垫面分布、城市热岛效应等多种因素的影响。
并且,还获得了适用于目标城市降雨事件模拟的最优模型配置方案,保证了降雨模拟结果
的准确度。然后,利用配置好的区域气象模型执行历史降雨模拟和/或未来降雨模拟,以分
别获得各个模拟网格单元的历史降雨时序数据和/或未来降雨时序数据,最后基于模拟得
到的历史降雨时序数据和/或未来降雨时序数据分别计算历史暴雨强度和/或未来暴雨强
度,基于区域气象模型的模拟数据来进行暴雨强度计算,不仅可以很好地适用于城市气象
监测站点数据少、建站时间短、降雨记录中断/缺失、降雨记录尚不足以支撑暴雨强度分析
的地区,而且能够应用于预测气候变化背景下未来不同时期城市设计暴雨强度的变化过
程,为城市涉水系统的设计、城市防洪排涝、海绵城市建设等提供了精细化设计、评估的工
具。
如上所述的方法的步骤。
上所述的方法的步骤。
片(punch cards)、纸带(paper tape)、任何其余的带有洞的图案的物理介质、随机存取存
储器(RAM)、可编程只读存储器(PROM)、可抹除可编程只读存储器(EPROM)、快闪可抹除可编
程只读存储器(FLASH‑EPROM)、其余任何存储器芯片或卡匣、或任何其余可让计算机读取的
介质。指令可进一步被一传输介质所传送或接收。传输介质这一术语可包含任何有形或无
形的介质,其可用来存储、编码或承载用来给机器执行的指令,并且包含数字或模拟通信信
号或其与促进上述指令的通信的无形介质。传输介质包含同轴电缆、铜线以及光纤,其包含
了用来传输一计算机数据信号的总线的导线。
改、等同替换、改进等,均应包含在本发明的保护范围之内。