Please wait a minute...
浙江大学学报(工学版)
浙江大学学报(工学版)  2019, Vol. 53 Issue (2): 347-354    DOI: 10.3785/j.issn.1008-973X.2019.02.018
水资源与海洋工程     
采用Sobol方法的暴雨径流管理模型参数灵敏度分析
周云峰1(),周永潮1,郑春华2,翁献明2,马晓宇2,柳景青1,*()
1. 浙江大学 建筑工程学院，浙江 杭州 310058
2. 温州市排水有限公司，浙江 温州 325000
Sensitivity analysis of parameters of storm water management model with Sobol method
Yun-feng ZHOU1(),Yong-chao ZHOU1,Chun-hua ZHENG2,Xian-ming WENG2,Xiao-yu MA2,Jing-qing LIU1,*()
1. College of Civil Engineering and Architecture, Zhejiang University, Hangzhou 310058, China
2. Wenzhou Drainage Co. Ltd, Wenzhou 325000, China
 全文: PDF(881 KB)   HTML
摘要：

为了定量探究典型城市暴雨径流管理模型（SWMM）中水文水力参数的全局灵敏度，利用基于方差分解的Sobol方法，计算不同雨强下SWMM中各参数对径流量、峰值流量以及峰现时间的灵敏度，简要分析结果的收敛性. 研究结果表明，在不同雨强下，SWMM中水文水力模块的灵敏参数及相关排序有所差异. 在大雨强条件下，不渗透地区表面的曼宁糙率对3类输出变量都是最灵敏参数，总阶灵敏度均大于0.6；混凝土管段的曼宁糙率系数对峰现时间的灵敏度远大于对其他2类输出变量的灵敏度. 在小雨强条件下，无洼地蓄水不渗透性百分比为灵敏参数，不渗透性地面洼地储蓄为最灵敏参数，第1阶灵敏度指数为0.396. 在模型前63 000次运行中，计算得到的Sobol指数上下波动，在运行84 000次之后，Sobol指数逐渐趋于收敛.
关键词： 灵敏度;  Sobol方法;  暴雨径流管理模型(SWMM);  模型参数;  定量分析    
Abstract:

The Sobol method based on variance decomposition was used to calculate the sensitivity of parameters in the storm water management model (SWMM) to runoff, peak flow, and peak time under two different rainfall conditions, and the convergence of the results was briefly analyzed, in order to quantitatively explore the global sensitivity of hydro-hydraulic parameters in the typical urban SWMM. Results showed that the sensitive parameters and relative ranking of hydro-hydraulic modules in SWMM were different under different rainfall intensities. The impervious areas’ Manning roughness was the most sensitive parameter for all three kinds of output variables under the heavy rainfall intensity condition, and the total sensitivity values were all greater than 0.6. The Manning roughness coefficient of the concrete pipe section was much more sensitive to the peak time than to the other two types of output variables. Percent of impervious area with no depression storage was sensitive parameter, and the impervious ground surface store was the most sensitive parameter with the first-order sensitivity index of 0.396 in the small rainfall intensity condition. The calculated Sobol indexes fluctuated in the first 63 000 runs of the model and tended to converge after 84 000 runs.
Key words: sensitivity    Sobol method    storm water management model (SWMM)    model parameter    quantitative analysis
收稿日期: 2018-01-08 出版日期: 2019-02-21
CLC:  X 143  
通讯作者: 柳景青     E-mail: zhouyunfeng1993@zju.edu.cn;liujingqing@zju.edu.cn
作者简介: 周云峰(1993—)，男，硕士生，从事排水管网模型优化校核研究. orcid.org/0000-0001-8552-3211. E-mail: zhouyunfeng1993@zju.edu.cn
服务　　
把本文推荐给朋友
加入引用管理器
E-mail Alert
作者相关文章　　
周云峰
周永潮
郑春华
翁献明
马晓宇
柳景青

引用本文:

周云峰,周永潮,郑春华,翁献明,马晓宇,柳景青. 采用Sobol方法的暴雨径流管理模型参数灵敏度分析[J]. 浙江大学学报(工学版), 2019, 53(2): 347-354.

Yun-feng ZHOU,Yong-chao ZHOU,Chun-hua ZHENG,Xian-ming WENG,Xiao-yu MA,Jing-qing LIU. Sensitivity analysis of parameters of storm water management model with Sobol method. Journal of ZheJiang University (Engineering Science), 2019, 53(2): 347-354.

链接本文:

http://www.zjujournals.com/eng/CN/10.3785/j.issn.1008-973X.2019.02.018    或    http://www.zjujournals.com/eng/CN/Y2019/V53/I2/347

图 1  研究区SWMM排水管网模型
降雨
日期 		降雨历
时/min 		降雨
编号 		降雨记录
间隔/min 		降雨峰值强度/
（mm/5 min） 		降雨总
量/mm
2017-07-08 135 Rain1 5 7.40 43.30
2017-06-28 20 Rain2 1 1.32 3.26
表 2  暴雨径流管理模型灵敏度分析所选降雨场次
分类类别 标签 地表属性 IL/% AL/hm2 PL/%
屋顶 A 混凝土，瓷砖 95 35.112 6 62.12
公用街道 B 沥青，水泥 95 1.845 3 3.27
绿化人行道 C 绿色植被，沥青 70 2.714 7 4.80
绿地 D 绿色植被 10 3.540 4 6.26
道路及旁
绿化带 		E 绿色植被，沥青，
水泥 		85 6.991 4 12.37
绿色屋顶 F 混凝土，瓷砖，
微小植被 		80 6.320 7 11.18
总计 ─ ─ 85.6 56.525 100.00
表 1  研究区土地利用情况
参数名称 物理意义 单位 所属模块 最小值 最大值 初始估计值
si 不渗透性地面洼地储蓄 mm 水文径流模块 0.18 2.54 1.27
sp 渗透性地面洼地储蓄 mm 水文径流模块 2.54 5.08 3.78
Ni 不渗透地区表面曼宁粗糙系数 ? 水文径流模块 0.010 0.050 0.025
Np 渗透性地区表面曼宁粗糙系数 ? 水文径流模块 0.020 0.400 0.131
C1 混凝土管段曼宁糙率系数 ? 动态演算模块 0.011 0.017 0.015
C2 塑性管段曼宁糙率系数 ? 动态演算模块 0.011 0.015 0.013
C3 PVC管段曼宁糙率系数 ? 动态演算模块 0.008 0.012 0.010
Rmax 霍顿模型最大下渗系数 mm/h 水文下渗模块 50.8 101.6 76.2
Rmin 霍顿模型最小下渗系数 mm/h 水文下渗模块 7.6 15.2 12.7
Dc 霍顿模型衰减系数 h?1 水文下渗模块 2 7 4
Ks 子汇水区坡度系数 ? 水文径流模块 0.67 1.67 1.00
Kw 子汇水区特征宽度系数 ? 水文径流模块 0.6 1.2 0.8
PA 屋顶不渗透百分比 % 水文径流模块 90 100 95
PB 公用街道不渗透百分比 % 水文径流模块 90 100 95
PC 绿化人行道不渗透百分比 % 水文径流模块 65 85 70
PD 绿地不渗透百分 % 水文径流模块 0 15 10
PE 道路及旁绿化带不渗透百分比 % 水文径流模块 75 95 85
PF 绿色屋顶不渗透百分比 % 水文径流模块 75 90 80
PO 无洼地蓄水不渗透性百分比 % 水文径流模块 10 35 25
表 3  暴雨径流管理模型水文水力模块中主要参数及取值范围、初始估计值
图 2  不同降雨强度下参数灵敏度指数比较
径流量灵敏度
参数 Si 参数 $\scriptstyle{S_i^{\text{T}}}$
si 0.396 si 0.504
Ni 0.139 Ni 0.178
PO 0.067 PO 0.103
Kw 0.023 Kw 0.089
Ks 0.017 Ks 0.082
PA ?0.002 C1 0.079
C1 ?0.002 PA 0.073
PD ?0.003 PD 0.069
Rmax ?0.003 PE 0.061
$ \cdots $ $ \cdots $ $ \cdots $ $ \cdots $
PB ?0.005 Np 0.051
$ \cdots $ $ \cdots $ $ \cdots $ $ \cdots $
Dc ?0.012 PF 0.045
总计 0.560 总计 1.733
表 5  Rain2下暴雨径流管理模型参数灵敏度分析
径流量灵敏度 峰值流量灵敏度 峰现时间灵敏度
参数 Si $\scriptstyle{S_i^{\text{T}}}$ 参数 Si $\scriptstyle{S_i^{\text{T}}}$ 参数 Si $\scriptstyle{S_i^{\text{T}}}$
Ni 0.502 0.682 Ni 0.414 0.609 Ni 0.482 0.635
Kw 0.138 0.259 Kw 0.144 0.263 Kw 0.111 0.223
Ks 0.051 0.142 Ks 0.058 0.169 C1 0.102 0.233
C1 0.011 0.092 C1 0.037 0.140 Ks 0.047 0.132
si 0.010 0.081 Dc 0.012 0.069 si 0.011 0.077
PA ?0.006 0.056 si 0.010 0.062 Dc 0.007 0.081
PF ?0.007 0.056 PA ?0.001 0.072 PA 0.003 0.075
PD ?0.008 0.056 PD ?0.001 0.051 PO 0.002 0.073
Rmax ?0.008 0.056 PB ?0.002 0.052 PD ?0.002 0.072
$ \cdots $ $ \cdots $ $ \cdots $ $ \cdots $ $ \cdots $ $ \cdots $ $ \cdots $ $ \cdots $ $ \cdots $
PB ?0.009 0.056 sp ?0.006 0.055 PB ?0.006 0.065
$ \cdots $ $ \cdots $ $ \cdots $ $ \cdots $ $ \cdots $ $ \cdots $ $ \cdots $ $ \cdots $ $ \cdots $
Dc ?0.009 0.058 PF ?0.010 0.053 PE ?0.011 0.062
总计 0.594 2.051 总计 0.597 2.029 总计 0.686 2.266
表 4  Rain1下暴雨径流管理模型参数灵敏度分析
图 3  径流量计算所得参数灵敏度指数收敛性
1 SCHUELER T R The importance of imperviousness[J]. Watershed Protection Techniques, 1994, 1 (3): 100- 111
2 ARNOLD C L, GIBBONS J Impervious surface coverage: the emergence of a key environmental indicator[J]. Journal of the American Planning Association, 1996, 62 (2): 243- 258
doi: 10.1080/01944369608975688
3 JAKEMAN A J, LETCHER R A, NORTON J P Ten iterative steps in development and evaluation of environmental models[J]. Environmental Modelling and Software, 2006, 21 (5): 602- 614
doi: 10.1016/j.envsoft.2006.01.004
4 ROSSMAN L A. Storm water management model user's manual [EB/OL]. [2018-01-08] . https://csdms.colorado.edu/w/images/Epaswmm5_user_manual.pdf.
5 JIANG L, CHEN Y, WANG H. Urban flood simulation based on the SWMM model [C]// 3rd Remote Sensing and Hydrology Symposium and the 3rd International Conference of GIS/RS in Hydrology, Water Resources and Environment. Guangzhou: PIAHS, 2014: 186–191.
6 WU X, WANG Z, GUO S, et al Scenario-based projections of future urban inundation within a coupled hydrodynamic model framework: a case study in Dongguan City, China[J]. Journal of Hydrology, 2017, 547: 428- 442
doi: 10.1016/j.jhydrol.2017.02.020
7 KNIGHTON J O, WALTER M T Critical rainfall statistics for predicting watershed flood responses: rethinking the design storm concept[J]. Hydrological Processes, 2016, 30 (21): 3788- 3803
doi: 10.1002/hyp.v30.21
8 SONG X, ZHANG J, ZHAN C, et al Global sensitivity analysis in hydrological modeling: review of concepts, methods, theoretical framework, and applications[J]. Journal of Hydrology, 2015, 523: 739- 57
doi: 10.1016/j.jhydrol.2015.02.013
9 黄金良, 杜鹏飞, 何万谦, 等 城市降雨径流模型的参数局部灵敏度分析[J]. 中国环境科学, 2007, 27 (4): 549- 553
HUANG Jin-liang, DU Peng-fei, HE Wan-qian, et al Local sensitivity analysis for urban rainfall runoff modelling[J]. China Environmental Science, 2007, 27 (4): 549- 553
doi: 10.3321/j.issn:1000-6923.2007.04.024
10 TSAI L Y, CHEN C F, FAN C H, et al Using the HSPF and SWMM models in a high pervious watershed and estimating their parameter sensitivity[J]. Water, 2017, 9 (10): 780
doi: 10.3390/w9100780
11 ZEYNEP A, BASAK G Assessing the sensitivity of SWMM to variations in hydrological and hydraulic parameters: a case study for the city of Istanbul[J]. Global Nest Journal, 2016, 18 (4): 831- 841
doi: 10.30955/gnj.001717
12 史蓉, 庞博, 赵刚, 等 SWMM模型在城市暴雨洪水模拟中的参数敏感性分析[J]. 北京师范大学学报: 自然科学版, 2014, (5): 456- 460
SHI Rong, PANG Bo, ZHAO Gang, et al Parameter sensitivity analysis of SWMM model in urban storm flooding simulation[J]. Journal of Beijing Normal University: Natural Science, 2014, (5): 456- 460
13 SUN N, HONG B, HALL M Assessment of the SWMM model uncertainties within the generalized likelihood uncertainty estimation (GLUE) framework for a high-resolution urban sewershed[J]. Hydrological Processes, 2014, 28 (6): 3018- 3034
14 王荫茵, 陶涛, 吕永鹏 城市雨水管网模型的水文参数灵敏度分析[J]. 城市道桥与防洪, 2016, (2): 159- 162
WANG Yin-yin, TAO Tao, LV Yong-peng, et al Analysis of hydrological parameters sensitivity in urban rainwater pipe network model[J]. Urban Roads Bridges and Flood Control, 2016, (2): 159- 162
15 高颖会, 沙晓军, 徐向阳, 等 基于Morris的SWMM模型参数敏感性分析[J]. 水资源与水工程学报, 2016, 27 (3): 87- 90
GAO Ying-hui, SHA Xiao-jun, XU Xiang-yang, et al Sensitivity analysis of SWMM model parameters based on Morris method[J]. Journal of Water Resources and Water Engineering, 2016, 27 (3): 87- 90
16 SARRAZIN F, PIANOSI F, WAGENER T Global sensitivity analysis of environmental models: convergence and validation[J]. Environmental Modelling and Software, 2016, 79: 135- 152
17 YU P S, YANG T C, CHEN S J Comparison of uncertainty analysis methods for a distributed rainfall-runoff model[J]. Journal of Hydrology, 2001, 244 (1/2): 43- 59
doi: 10.1016/S0022-1694(01)00328-6
18 KREBS G, KOKKONEN T, VALTANEN M, et al A high resolution application of a stormwater management model (SWMM) using genetic parameter optimization[J]. Urban Water Journal, 2013, 10 (6): 394- 410
doi: 10.1080/1573062X.2012.739631
19 ZHANG C, CHU J, FU G Sobol′s sensitivity analysis for a distributed hydrological model of Yichun River Basin, China[J]. Journal of Hydrology, 2013, 480: 58- 68
20 齐伟, 张弛, 初京刚, 等 Sobol’方法分析TOPMODEL水文模型参数敏感性[J]. 水文, 2014, 34 (2): 49- 54
QI Wei, ZHANG Chi, CHU Jing-gang, et al Sensitivity analysis of TOPMODEL hydrological model parameters Based on Sobol’ method[J]. Journal of China Hydrology, 2014, 34 (2): 49- 54
doi: 10.3969/j.issn.1000-0852.2014.02.010
21 张质明, 王晓燕, 李明涛 基于全局敏感性分析方法的WASP模型不确定性分析[J]. 中国环境科学, 2014, 34 (5): 1336- 1346
ZHANG Zhi-ming, WANG Xiao-yan, LI Ming-tao Uncertainty analysis of WASP based on global sensitivity analysis method[J]. China Environmental Science, 2014, 34 (5): 1336- 1346
22 GHASEMIZADE M, BARONI G, ABBASPOUR K, et al Combined analysis of time-varying sensitivity and identifiability indices to diagnose the response of a complex environmental model[J]. Environmental Modelling and Software, 2017, 88: 22- 34
23 SOBOL I M Global sensitivity indices for nonlinear mathematical models and their Monte Carlo estimates[J]. Mathematics and Computers in Simulation, 2001, 55 (1?3): 271- 280
doi: 10.1016/S0378-4754(00)00270-6
24 SALTELLI A Sensitivity analysis for importance assessment[J]. Risk Analysis, 2002, 22 (3): 579- 590
doi: 10.1111/risk.2002.22.issue-3
25 NOSSENT J, ELSEN P, BAUWENS W Sobol’ sensitivity analysis of a complex environmental model[J]. Environmental Modelling and Software, 2011, 26 (12): 1515- 1525
doi: 10.1016/j.envsoft.2011.08.010
26 BARCO J, WONG K M, STENSTROM M K Automatic calibration of the U. S. EPA SWMM model for a large urban catchment[J]. Journal of Hydraulic Engineering, 2008, 134 (4): 466- 474
doi: 10.1061/(ASCE)0733-9429(2008)134:4(466)
27 李春林, 胡远满, 刘淼, 等 SWMM模型参数局部灵敏度分析[J]. 生态学杂志, 2014, 33 (4): 1076- 1081
LI Chun-ling, HU Yuan-man, LIU Miao, et al Local sensitivity analysis of parameters in storm water management model[J]. Chinese Journal of Ecology, 2014, 33 (4): 1076- 1081
[1] 黄腾逸,周瑾,徐岩,孟凡许. 基于多场耦合分析的磁流变阻尼器建模与结构参数影响[J]. 浙江大学学报(工学版), 2020, 54(10): 2001-2008.
[2] 程锋,张栋,唐硕. 高超声速运载器宽速域气动/推进耦合建模与分析[J]. 浙江大学学报(工学版), 2019, 53(5): 1006-1018.
[3] 周瑾,高天宇,陈怡,席鸿皓. 采用动力测试的双转子卧螺离心机模型修正[J]. 浙江大学学报(工学版), 2019, 53(2): 241-249.
[4] 潘秋萍, 汪震, 甘德强, 谢欢, 李尚远. 基于数值微分的阻尼灵敏度方法比较[J]. 浙江大学学报(工学版), 2018, 52(1): 174-183.
[5] 刘凯, 曹毅, 周睿, 葛姝翌, 丁锐. 一移两转平板折展柔性铰链的建模及优化[J]. 浙江大学学报(工学版), 2017, 51(12): 2399-2407.
[6] 臧勇, 赵泽波, 秦勤, 于洋. 奥氏体晶粒尺寸对工艺参数的灵敏度方程及应用[J]. 浙江大学学报(工学版), 2015, 49(4): 662-669.
[7] 吕云腾, 祝长生. 基于温漂补偿的高温电涡流位移传感器[J]. 浙江大学学报(工学版), 2015, 49(4): 749-753.
[8] 齐虎,李云贵,吕西林. 混凝土弹塑性损伤本构模型参数及其工程应用[J]. 浙江大学学报(工学版), 2015, 49(3): 547-554.
[9] 刘杰, 王宇平. 一种基于单纯形法的改进中心引力优化算法[J]. 浙江大学学报(工学版), 2014, 48(7): 2-.
[10] 张进锋, 杨强, 颜文俊. 基于模型参数聚类的分解炉温度建模[J]. 浙江大学学报(工学版), 2014, 48(12): 2139-2143.
[11] 刘杰, 王宇平. 一种基于单纯形法的改进中心引力优化算法[J]. 浙江大学学报(工学版), 2014, 48(12): 2115-2122.
[12] 尹喜珍, 马成炎, 叶甜春, 肖时茂, 于云丰. 高灵敏度GNSS接收机频率合成器设计[J]. J4, 2013, 47(1): 70-76.
[13] 张杨 周晓军, 杨辰龙, 李雄兵, 郑慧峰. 变厚度曲面构件超声检测灵敏度补偿[J]. J4, 2013, 47(1): 116-121.
[14] 乔龙学, 王伟刚, 何姗, 鲁阳. 双悬臂梁直径测量传感器性能分析[J]. J4, 2012, 46(6): 1060-1066.
[15] 杨剑,王智化,张彦威,陈云,周俊虎,岑可法. 热化学硫碘开路循环制氢系统的设计与模拟[J]. J4, 2011, 45(5): 869-877.