新安江模型参数局部灵敏度分析.docVIP

新安江模型参数局部灵敏度分析 第41卷第1期 2 0 1 0年1月 人 民 长 江Yangtze River Vo.l 41,No. 1Jan., 2010 收稿日期:2009-08-24 基金项目:国家自然科学基金 作者简介:薄会娟,女,硕士,主要从事水文水资源研究。 文章编号: 1001-4179(2010)01-0025-04 新安江模型参数的局部灵敏度分析 薄会娟,董晓华,邓 霞 (三峡大学土木水电学院,湖北宜昌443002) 摘要:新安江模型是我国应用广泛的一种降雨径流流域模型,适用于湿润与半湿润地区。新安江模型参数较 多,人工率定时需要首先率定对模拟结果影响最为敏感的参数,因此对参数进行灵敏度分析是必要的。采用 水文模拟法与摩尔斯分类筛选法,对清江流域1989年9月30日至10月5日次洪水进行径流模拟,并让参数 变化一定的步长,用Nash系数对模拟结果的精度进行评定,从而反映出参数的灵敏度。结果表明,最灵敏的 参数是FC与Cg;Kc与WM是灵敏参数; b对峰值流量是高灵敏参数,而对径流总量来说是灵敏参数;WUM, WLM,WDM,C与IMP是不灵敏参数。 关 键 词:参数;局部灵敏度分析;新安江模型;清江流域 中图法分类号:TV131. 4; P338. 9 文献标志码:A 新安江模型是河海大学赵人俊等于1973年在作 新安江水库入库径流预报工作时提出来的降雨径流流 域模型。最初的新安江模型为两水源模型, 20世纪80 年代初期,提出了三水源新安江模型,中期提出了四水 源新安江模型。 新安江模型把全流域分成若干单元流域,对每个 单元流域做产汇计算,得出单元流域的出口流量过程。 再进行出口以下的河道洪水演算,求得流域出口的流 量过程。把每个单元流域的出流过程相加,就求得了 流域出口的总出流过程[1]。本文采用两水源3层蒸发 模型,对次洪水进行模拟。 由于新安江模型的参数比较多,而且参数的灵敏 度不同,所以对参数的灵敏度进行分析有助于更好地 利用新安江模型。参数灵敏度分析的作用是定性或定 量地评价模型输入误差对模型输出结果的影响。灵敏 度分析主要分为局部灵敏度分析与全局灵敏度分 析[2, 3]。局部灵敏度分析检验单个参数的变化对模型 结果的影响程度;全局灵敏度分析考虑了参数在整个 取值空间的影响及参数之间的共同作用。 本文采用的是局部灵敏度分析的方法,分析单个 参数的灵敏度。局部灵敏度分析的优点在于其可操作 性。具有以下意义:1选出对模型输出结果影响大的 参数,对于这些参数,在模型分析过程中需要集中精力 尽可能地提高参数的准确度;而对于那些对模型结果 影响不大的不灵敏参数,只需选取其经验值,这会在很 大程度上减少模型参数率定与验证的工作量[4]。o 加深对模型的理解。不同参数的变化对模型的影响程 度与方式都不同,全面掌握参数对模型的影响程度与 方式,有助于在不同的模型使用条件下选择相应的敏 感参数以进行重点识别。?加深对所模拟系统行为 的理解。通过灵敏度分析可以发现所模拟系统对哪些 参数变化敏感,从而可以确定各模拟因子对所模拟系 统的影响程度。 1 新安江模型 新安江两水源模型的蒸发计算采用3层蒸发计算 模型。3层蒸发模型就是把整个流域划分为上层、下 层与深层3层,只有当上层蒸发完之后下层才开始蒸 发,下层蒸发完之后深层开始蒸发。 产流子模型采用蓄满产流模型,蓄满产流是指包 人 民 长 江2010年 气带土壤含水量达到田间持水量之前不产流,之前的 降雨全部补充包气带缺水量,蓄满后开始产流,下渗的 雨量形成地下径流,超渗的雨量形成地面径流。 利用稳定下渗率FC将径流划分为地面径流与地 下径流两种水源。地面径流采用单位线汇流,地下径 流采用一次线性水库汇流。模型把流域面积划分为透 水面积与不透水面积两部分,不透水面积上的降水在 满足蒸发后将直接转化为地面径流。透水面积上将发 生下渗,下渗的水量一部分存储于土壤层,后期耗于蒸 发;满足了流域土壤蓄水容量后下渗水量才能转化为 径流。不透水面积用参数IMP表示,它是用流域内不 透水面积占全流域面积的百分比表示的。 地面径流汇流所用的时段单位线的推求采用最小 二乘法,即:由降雨资料推求地面净雨过程,使得地面 净雨量和实测地面径流量相等。 新安江模型可分为4个层次:第一层是蒸散发,参 数有Kc,WUM,WLM,WDM,C;第二层是产流,参数包 括WM,b;第三层是水源划分,参数是FC;第四层是汇 流,参数是Cg。表1是新安江两水源模型主要参数的 物理意义。 表1 新安江模型主要参数及取值范围 参数物理意义取值范围 Kc流域蒸发折算系数 0~1 C深层蒸散发系数 0.10~0.20 WUM流域平均上层蓄水容量5~20