吉林市暴雨强度公式推求和设计雨型
- 格式:pdf
- 大小:1.27 MB
- 文档页数:8
第39卷第1期2021年2月干旱气彖 Vol. 39 No. 1Journal of Arid Meteorology February ,2021邱译萱,壬晓峰,杨雪艳,等•吉林市暴雨强度公式推求和设计雨型[J].干旱气象,2201,29(1):151 - 158, [ QIU Yixuan, REN Xiaofena , YANG Xueyaa , et al. Rainstorm Intensity Formula and Design Rainstorm Pattern for Jilia City[ J]. Journal of And Meteoroloyy , 2221 , 39(1):151 - 158], DOI:12.31755/j. issn. 1006 -7639(2201) -21 -2151
吉林市暴雨强度公式推求和设计雨型
邱译萱1,任晓峰2,杨雪艳1,姜忠宝1,李宇凡1(.吉林省气候中心,长白山气象与气候变化吉林省重点实验室,吉林 长春 10062;2.吉林省吉林市气象局,吉林吉林132213)摘 要:利用吉林市1961—2217年分钟降雨量数据,分别采用年最大值法和年多个样法选样,P -皿 型、指数型和Gumbei分布拟合频率分布曲线,根据误差最小原则选择最佳取样及拟合方法,采用最小 二乘法求解暴雨强度公式参数,并将新旧暴雨强度公式进行比较,再利用推求出的公式参数模拟吉林 市短历时芝加哥雨型。结果表明:通过年最大值法选样和p- m型频率分布曲线拟合得到的暴雨 强度公式精度最高。5 ~22 ma历时新编暴雨强度公式计算的暴雨强度值与旧编公式相当或偏 小,22 - 182 ma历时新编暴雨强度公式计算的暴雨强度值偏大。在相同重现期下,随着降雨历 时的延长,雨强变化率逐渐加大。吉林市短历时雨型的综合雨峰位置系数为2.389。32 -182 mm历 时雨型形态均为单峰型,各历时瞬时雨强峰值接近,雨峰位置位于偏整场降雨过程的1/2处之前。累 计雨量的变化特征与设计暴雨雨型形态一致。关键词:暴雨;暴雨强度公式;雨型设计;吉林市文章编号:1026 -7639(2221)21 -2151 -08 D0I:12.11755/j. issn. 1026 -7639(2221) -21 -2151 中图分类号:P468.3 + 22 文献标志码:A引言目前全球变暖趋势仍在进一步持续,在气候变 暖背景下极端降水事件更加严重和频繁[1]0自22 世纪66年代以来,中国极端降水频次呈增多趋 势[2],至21世纪末中国极端降水事件强度仍将增 强、频率仍将增多,暴雨洪涝风险逐渐加大J]o短 历时强降雨是诱发城市洪涝的直接因素,造成交通 中断、电力故障、城市污染等问题,严重威胁国民经 济和人民群众生命财产安全⑷。吉林市是吉林省 第二大城市,地形属于低山丘陵,因松花江干流流经 吉林市区,水资源丰富,对防汛要求较高。吉林市近 12 n来发生暴雨16次,其中2215年和2217年的暴 雨降水量大、时段集中,造成严重洪涝灾害,带来了 巨大损失J ] o加快提高城市防洪排涝能力,保证城 市排水管网安全,已成为当前城市水利规划工作亟 待解决的问题。暴雨强度公式及暴雨雨型是城市排水标准计 算、城市规划设计等工作的基础[],其可靠性关乎 城市的安全。我国雨水流量设计计算普遍采用的是 1697年修订的《室外排水设计规范》J ] (GBJ14 - 87)规定的推理公式,该版公式较为陈旧,已无法满 足当前现实需要,从科学角度考量,应根据最新的统 计资料进行完善修编。近年来,国内外对城市暴雨 强度公式和设计雨型的相关研究主要集中在暴雨样 本筛选、频率拟合方法及雨型模型的选取。筛选暴 雨样本常用的方法主要有年最大值法和年多个样 法[8"12 ],基于这两种方法的重现期对应转换关系已 有相关研究[I5-15 ] o随着国内城市降雨资料观测年 限的增加,年最大值法选取的样本独立性增强,且在 数据收集和统计上相对容易,应用最为广泛。目前 较为流行的频率拟合方法为理论频率曲线适线法, 即对实测暴雨样本进行适线调整,得出重现期(P)、 降雨强度()和降雨历时()三者的经验数据表,以 此为基础计算暴雨强度公式。一般常用的频率分布收稿日期:2222 -29 -06;改回日期:222。-2 -31基金项目:国家级项目气象风险预警基本业务支撑能力建设(1410090012016028)和吉林省气象局科研课题“吉林省暴雨强度公式优化 和雨型设计研究”(2。1605)共同资助作者简介:邱译萱(1999—),女,硕士,工程师,主要从事气候变化及气候资源研究.E-maii: *****************通信作者:任晓峰(1977—),男,本科,高级工程师,主要从事天气预报及人工影响天气研究• E - maii: rrtte@naa. com
。152干 旱 气 象39卷曲线包括P -皿型、Gumbel型、指数型和对数正态 分布曲线[16_15]0 ALHASSOUN[15]在估算利雅得地 区的降雨强度公式时,对P -皿型、Gumbel型和对 数正态分布曲线3种频率拟合方法进行了对比,认 为理论频率曲线对暴雨样本拟合效果最好的,计算 得到的暴雨强度公式不一定精度最高。因此在计算 暴雨强度公式时,需要综合考虑选样方法、当地实际 暴雨分布规律、频率曲线误差和公式误差来确定最 适宜的暴雨强度频率曲线。暴雨雨型反映了暴雨过 程在时间尺度上的分布变化,常用的雨型有芝加哥 雨型、Huff雨型、Pilgrim & Cordery雨型、三角形雨 型、均匀雨型、同频率分析方法雨型等[20'23]o其中 芝加哥雨型的雨强过程线的计算方法相对简单,同 时雨峰位置系数可以根据暴雨强度公式计算得出, 因此在相关研究中应用较多。因极端降水事件特征不断变化,本研究利用 1991-2017年的逐分钟降雨量资料重新计算吉林 市暴雨强度公式,并首次进行短历时雨型设计,对吉 林市当前暴雨特性进行研究,以期为城市排水系统 设计及海绵城市规划提供参考。1资料与方法1 -1资料来源及合理性分析以吉林市城郊气象站1791—2217年逐分钟降 雨量为基础数据,其中1791—2204年分钟降雨数据 为降水自记纸记录经数字化处理及人工审核或修正 后得到的资料;2242—2417年分钟降雨数据从自动 站资料中获取。资料均来自吉林省气象档案馆。吉林市城郊气象站为吉林市区仅有的】个国家 级气象观测站,具有长期降水观测° 196—247 年,吉林市城郊气象站经历两次迁站,观测场海拔由 2043年第一次迁站前的784.】m变为798. 2叫 2417年第二次迁站海拔变为198. 8 m,与市区平均 海拔接近,具有区域代表性。对逐年降雨量资料采 用标准正态方法(SNHT)进行迁站前后时段的显著 性检验(显著性水平0.05),检验结果均不显著,说 明迁站前后的降水资料具有良好的完整性和均 一性01-0方法7 2.1暴雨样本序列建立及选择采用逐分钟滑动统计法,按照“不漏场次、不漏 大值”的原则,挑取 5、10、15、20、30、45、60、90、122、 172,178 min共7个历时的所有场次降雨数据。年最大值法:选取逐年各历时最大雨量,各历时 样本数为57个,作为建立短历时暴雨强度公式的统 计样本0年多个样法:将逐年各历时前8场最大雨量 (各历时基础样本数为57 X 8 =456个)从大到小 进行排序,并按从大到小选择资料年份的4倍 (各历时最终样本数为57 X4 =228个)数据,作 为统计样本0HO短历时强降水变化极值及均值( 7)短历时强降水的极值将各历时以年最大值法选取的统计样本从大到 小进行排序,选取最大的 70个降水极值,将资料年 份均分为前后两段,通过统计前后两个时段的降水 极值样本数,分析短历时强降水的极值分布特征0(2)短历时强降水的均值统计各历时每年10场最大雨量,对比不同年代 的多年暴雨量平均值,分析短历时强降水的均值分 布特征07.2-2暴雨强度公式(7样本频率和重现期计算及理论频率分布曲 线拟合将统计样本按照降序排列,计算样本的经验频 率,重现期P与经验频率P互为倒数。P计算公式 如下:
式中:m为样本排序数;为样本总数。基于选取的统计样本,采用理论频率曲线对样 本的经验频率数据序列进行趋势性拟合,根据暴雨 强度公式的重现期重点关注区间(2 ~22 a)和范围 要求(2-in a[25],对频率分布曲线进行适线调整 和外延。从P-皿型、Gumbel型和指数型3种分布曲线 中,选取拟合效果较好的理论频率曲线进行暴雨强 度公式拟合0P-皿型分布即伽玛分布,其概率密度函数如 下:/()= -久广甘心"0) (2)『(a)式中:a、和S为P -皿型分布包含的参数;r(a) 为a的伽马函数。3个原始参数经过换算,可用3个统计参数力、 C、和 CB 表示:a = 4/Cs2, = 2/xCvCs,a0 =斤(- 2Q/C") ° Cv为离差系数,C为偏差系数,"为均 值。这3个统计参数可通过矩法进行确定。Gumber型分布又称为极值I分布,其概率密度 函数如下
:第/期邱译萱等1吉林市暴雨强度公式推求和设计雨型159/() = aexp(-厂。0) (3)式中-y = a(x -b) ,a为分布的尺度参数,为分布 的位置参数。Gumbei型分布是P -皿分布参数固定的一 种特殊情况,Cs固定为1. 140,所以只有 2个参数。指数型分布概率密度函数如下:/()) = (6)式中:a为离散程度(a >0) ;2为分布的下限(对 城市短历时暴雨,2>0)o 2个原始参数经换算也 可用)、C、Cs 表示:a = 1/斤C,2 =)(1 -CJ,Cs = 2o(2)暴雨强度公式拟合依据《室外排水设计规范)[22](简称《规范》), 暴雨强度公式如下:q167 缶(1 + C(+b)"⑸式中:[L - (s - hm2)-1]为暴雨强度;P(a)为重现期;(min)为降雨历时;缶为雨力参数;;为雨力变 动参数;为降雨历时校正参数;为暴雨衰减指数, 与重现期有关由3种理论频率分布曲线确定P-i-t三联 表,采用最小二乘法计算参数,得到暴雨强度公 式。(3)精度检验精度检验重点为重现期2-22 a[25],在雨强一般的地区,暴雨强度平均绝对均方根误差bR均不 宜超过0. 35 mm - mid "1 ,在雨强较大的地区,平均相对均方根误差bq均不宜超过5% o 平均绝对均方根误差:(6)平均相对均方根误差:b =槡:打xl()2% ⑺ 式中:R:、R»(mm - mid"1)分别为暴雨强度公式计算 的理论分钟降雨量及频率拟合分钟降雨量(即三联 表对应的值- (s - hm2)-]分别为暴雨 强度公式计算的理论暴雨强度值及频率曲线拟合的 暴雨强度值;为样本数。1.2.4新、旧版暴雨强度公式对比(1)时段暴雨量极值与均值对比在两版公式对应资料年限内,统计各历时暴雨 量极大值,分析极值差异;统计各历时每年1场最 大雨量,求取多年平均值,分析均值差异。(2)暴雨强度值变化率分析对两版公式计算的暴雨强度值变化率进行分 析,其计算公式如下:△ q = qa _ 弘 x 102% (8)q。式中:Aq为变化率;q.和q。分别为新、旧公式计算 的雨强。1.2.2芝加哥法设计雨型(1)建立样本序列。挑取30、62、90、120、150、152 min共6个历时的所有场次降雨数据,记录每个 过程的起止时间及过程内逐分钟雨量值、累计雨量 值,将各历时的过程累计雨量值从大到小进行排序, 选择大约的资料年份数(各历时样本数约为57个) 的数据,作为统计样本。将过程累计雨量最小值,作 为该历时的暴雨雨量阈值。(2) 雨峰位置系数将确定的各历时统计样本以 5 min间隔为步长,滑动计算5 min降水量,挑出各 样本雨量最大的时刻,记为雨峰位置,雨峰位置出现 时间除以该样本历时,得到雨峰位置系数;计算各样 本雨峰位置系数的平均值,得到各历时的雨峰位置 系数。以历时为权重,对各历时的雨峰位置系数求 平均,得到芝加哥雨型的综合雨峰位置系数。(3) 芝加哥降雨过程线模型确定结合1.2.3中计算得到的新公式参数,计算雨峰前、后瞬时降雨强度。计算公式如下: 雨峰发生前(上升段):A [ (1 - n) itr + b] (b/r +b)"+ 雨峰发生后(下降段):A [1 _ n)//r1 —) + b] (//(1 -) +b)"+(9)