唐清竹，徐宗学*，王京晶，赵 刚，李 鹏

(1. 北京师范大学 城市水循环与海绵城市技术北京市重点实验室，北京 100875; 2. 北京师范大学 水科学研究院，北京 100875; 3. 布里斯托大学 地理科学学院，布里斯托 bs8 1th)

随着全球气候的持续变暖和城市化进程的快速推进，城市洪涝问题频发。目前，城市洪涝灾害已经成为影响城市经济社会发展和人民生命财产安全的主要自然灾害之一，防灾减灾问题已经引起了中央和地方各级政府的广泛关注和高度重视。在中国目前大力建设海绵城市以缓解城市内涝的大背景下，城市雨洪过程模拟逐渐成为一种有效的防灾减灾措施与雨水管理非工程措施。swmm模型(storm water management model)是由美国环境保护署(epa)于1971年开发的一款暴雨洪水管理模型，随着几十年的不断完善，其功能愈发强大，可以进行动态的降雨-径流模拟计算，并因其开源免费、功能强大、简单易学等特点被广泛应用于城市防汛排洪及雨水污水管道系统等的规划、设计和评价。模型参数是swmm模型的重要组成部分，对参数敏感性进行准确的识别，是保证模型模拟精度和提高参数率定效率的重要前提。

参数敏感性分析是模型模拟中重要的环节，可以对不同情况下模型的参数率定与模型应用提供一定依据。敏感性分析可以识别不同参数对模拟结果的影响程度，从而筛选出对输出值高敏感和不敏感的参数。高敏感的参数在模型率定时需重点关注，不敏感的参数则可以直接选取经验值，这不仅可以提高模型参数率定效率，还能提高模拟结果准确性。目前，参数敏感性分析主要分为全局敏感性分析和局部敏感性分析。①全局敏感性分析是综合考虑参数整体作用对模型输出结果的总体影响，从而得到全局最优解。常用的全局敏感性分析方法主要有普适似然不确定性估计方法(glue)、sobol方差分析法、区域灵敏度分析法(hsy)以及多元逐步回归分析法等，其中glue方法可以更明确地认识降雨径流模型的局限性，通过提供模型参数的后验分布，定性分析出模型结果不确定性范围，能够很好地解释“异参同效”现象，避免单一最优参数组合造成的洪水预警预报风险。②局部敏感性分析也叫单参数分析，其原理是控制其他参数不变，依次对某一个参数在一定范围内按照固定步长进行人工扰动，分析单个参数对模拟结果的影响。目前国内应用最多的局部敏感性分析方法是修正morris筛选法，该方法较于原始morris筛选法来说，操作简便、容易计算且计算结果精度更高，可定量筛选出对模型结果影响显著的参数。全局敏感性分析和局部敏感性分析的原理差异较大，目前大多数研究主要着重关注某一种方法的参数敏感性分析结果，而对比研究两种方法的参数敏感性分析结果相对较少。

本文从全局和局部角度出发，分别采用glue方法和具有明显优势的修正morris筛选法对构建的swmm模型中11个参数进行敏感性分析，得出对模型模拟结果影响显著的参数并进行排序，对比两种方法的敏感性分析结果，分析总结这两种方法的优缺点及适用条件，为山东省济南市主城区及其他相似区域洪水模拟和流域水资源管理提供科学依据。

山东省济南市主城区面积约为322 km，属温带季风气候，四季分明、雨热同期，年平均气温约为13.8 ℃，多年平均降雨量约为647 mm。降水时空分布极不均匀，6月至9月为主汛期，降水量约占全年降水量的75%。市域水量丰沛，小清河在中北部横穿城市，流域内地形复杂。济南市属于典型的山前平原型城市，地势南高北低，南北高差百余米，其数字高程模型(dem)如图1所示。济南市主城区河流湖泊众多，主要包含黄河、小清河两大水系和大明湖、白云湖等湖泊。研究区主要属于小清河水系，是城区唯一的排洪干道，流域内共有5座雨量站，流域唯一出口断面设置有黄台桥水文站，用于监测降水径流等水文数据。降水数据用于驱动swmm模型，实测流量数据用于校验模型，流域水系如图1所示。

济南市洪涝灾害频发，暴雨高度集中、强度大、历时短、降水时空变化剧烈。老城区雨水管网系统建造时间较早，设计标准低，泄洪能力较差，同时由于近20年高度城市化导致下垫面发生剧烈改变，下渗能力低，汇水区自然调蓄能力弱。加之南高北低、落差大的地形特点，在遭遇短历时、高强度降水时，南部山区下泄洪水直接汇入主城区，导致雨洪猛涨缓落，极易在短时间内在低洼区域形成大面积严重积水，排泄不及，极易引发城市洪涝灾害。

图1 山东省济南市地形及水系概况fig.1 topography and water system maps of jinan city, shandong province

利用济南市主城区现有河道及排水管网数据，基于管网空间拓扑关系，利用arcgis软件对管道长度、流向、糙率和最大水深等数据进行处理和整理，将研究区内复杂的河道水系及管网资料进行合理概化。研究区排水管网经过概化共得到160个节点、1个排水口、156条管线。选用划分精度更高的泰森多边形法，结合人工修正方法划分各个子汇水区。首先，利用arcgis软件对整个研究区进行泰森多边形的自动划分，初步形成各个子汇水区域；然后，结合研究区流域水文特征与数字高程模型，依据研究区用地类型、街道建筑分布手动进行局部调整，从而进一步细划子汇水区。按照上述方法，共将济南市主城区流域排水系统划分为112个由河道和排水管网共同控制的子汇水区，流域出口位于黄台桥水文站。同时，提取swmm模型所需的参数数据(如子汇水区面积、子汇水区平均坡度和子汇水区特征宽度等)，将子汇水区、管道、节点文件以及各个属性数据通过inp. pins软件输入模型中，添加刘家庄雨量站、东红庙雨量站、兴隆雨量站、燕子山雨量站和黄台桥水文站，并添加处理好的降雨数据，输入各个子汇水区所对应的控制雨量站与排水节点，从而构建济南市主城区swmm模型，最终建立的swmm模型结构如图2所示。

图2 swmm模型概化fig.2 generalization of swmm model

swmm模型参数较多，根据获取方法可分为物理参数和过程参数两大类，其中物理参数主要依据实测资料或者通过arcgis软件进行计算与提取，而过程参数则需要通过城市雨洪模拟进行率定及验证得到。本研究中，参数敏感性分析对象为11个过程参数，通过参考相关文献结果及swmm模型用户手册确定其取值范围，相关过程参数取值范围见表1。

表1 swmm模型参数取值范围

glue方法是英国水文学家beven等于1992年提出的水文模型参数不确定性估计方法，其评估的是参数集，而不是参数的单一值，可通过设定敏感性阈值，筛选出多个高敏感性参数组合，使模型模拟结果达到最优。本研究中，全局敏感性分析步骤如下：①首先，采用拉丁超立方法在确定的11个参数取值范围内进行10 000次均匀随机抽样，将参数组代入swmm模型得出模拟值；②然后，选择合适的似然目标函数，本文以纳什效率系数(nash-sutcliffe efficiency coefficient,nse)为目标函数估算模拟值与实测值的拟合程度，得出似然度值；③接着，根据贝叶斯理论指定参数空间的先验分布为均一分布，以计算的参数似然度值为后验分布，每个参数集实现的后验似然度值将反映所有参数之间的相互作用；④最后，筛选出符合一定似然度(纳什效率系数大于0)的结果进行分析，将散点数进行统计，根据各个参数的取值范围合理划分单位区间，统计各单位区间中有效参数组合个数，计算各参数空间分布概率及累计概率。

morris筛选法是任意选取某个参数，在合理范围内进行随机扰动，得出相应输出结果()=(,,,…,)，用影响值判断输出值对各个参数的响应。其中，=1,2,…,。影响值计算公式为

=(-)

(1)

式中：为参数变化后的模型输出值；为参数变化前的模型输出值；为参数的变化幅度值。

本文采用修正morris筛选法对参数进行特定百分率的扰动，运行模型得到参数的敏感性系数序列，序列的平均值即为参数敏感性的判别因子。敏感性的判别因子()计算公式为

(2)

式中：为模型运行次数；为模型第次运行输出值；为初始参数运行得到的初始输出值；为第次运行参数相对于初始参数值的变化百分率。

参数的敏感性可根据的大小分为4类：Ⅰ类，||≥1，参数极敏感；Ⅱ类，0.2≤||