水文模型参数敏感性分析方法研究评述
关键 词 : 感性 分析 : 敏 局部 分析 法 ; 全局 分 析 法 ;H- L OAT分析 法
中图分 类号:3 49 P3. 2 文献 标识 码:A 文章 编 号:0 0 0 5 (0 00 — 0 9 0 1 0 — 8 22 1 )4 0 0 — 4
1 引 言
过 度参 数 化 是 水 文 模 型 中 的 难 点 问 题 【 特 别 是 ” . 对 于 分 布式 水 文 模 型 而 言 。 为 复 杂 。 在 实 际 应 用 尤 中 ,集合 了众 多参 数 的分 布式 水文 模 型还需 要 与模 型 的预报 作对 比 。 因此 , 数敏 感性 分 析 的方法 是 十分 必 参 要 的 , 须 能 适应 多 参数 , 必 同时考 虑 多 输 出变量 [ 3 1 。参 数 敏感 性分 析 方法 的 目的是 尽 量减 少参 数率 定 时需要
第 3 卷 第4 0 期
2 1年 8 0 0 月
水
文
V01 0 No4 . 3 .
Au . 01Biblioteka g,2 0J OURNAL O HI F C NA HYDROL Y OG
水 文模 型参数敏 感性分析方 法研 究评述
田 雨 , 雷 晓辉 , 云钟 , 一 , 蒋 李 薇
21 局 部 分 析 法 .
后 也成 为敏 感性 分 析 的方法 。
一
次 二 阶矩法 的步 骤是 先 选 择 理论 分 布族 g y = ()
g Y ) 近 M= ( t :… , ) (, 逼 f e , , e 的概 率分 布 , e 然后 用 泰勒 公 式 将 在 ( e, , ) 均值 处 展 开 , e,:… e 的 舍去 二 次 以 上 的 高 阶项 , 这样 近似 求 得 的 肘 的二 阶矩 。 而估 计 进 出参 数 。泰 勒公 式 的一 次项 需要 计算 某 一点 上 的模 型 输 出 . 确定 其导 数 : 并
(. 1天津 大学 建工 学 院 , 津 3 0 7 ; 天 0 0 2 2 中国水利 水 电科 学研 究 院水资 源研 究所 , 京 10 3 ; .水利 部水 文局 , . 北 0 0 83 北京 1 0 5 ) 0 0 3
摘
要: 参数 敏 感性 分析 一 直被认 为 是 用来确 定模 型 关键 参数 、 制模 型效 率 非 常有效 的过 滤 工具 。 同 控
式 中 :=e -, ) 与均 值 =F 一, 有关 的 随机 输 e{ - 为 e I ) 入 变量 ; 为 e在均 值 的导数 值 。 跏 处 输 出 函数 的均值 和标 准 差 o 可 以通 过 下式 计 r
S =
M (l…e, ,。 e , i… e)
垒
e
调 整 的参 数 个 数 1 识 别 哪 些 参 数 在 某 个 具 体 流域 对 4 ] .
模 拟精 度具 有 重要影 响 。
式 中 : 为模 型 输 出 代 表 不 同 的模 型参 数 ; 代 表 对 e参 数 的微 小扰 动 。这种 方 法要 求模 型运 行 的次数 并 不多 , 早期 研究 中很 流行 。 在
一
2 参 数 敏 感 性 分 析 方 法 评 述
参 数 的多少 是对 现 存参 数敏 感性 分 析方 法分 类 的
一
次 二 阶矩法 通过 计算 某 一点 的导 数估 计 出模 型
输 出的平 均 值 ( 次 二 阶原 点 矩 ) 标 准 差 ( 阶 中心 一 和 二
矩 ) 。该方 法 最 初是 为 解 决 不确 定 性 问题 而提 出的 , 同
个重 要依 据 ,一般 可 以分 为局 部分 析 法和 全局 分 析
法[ 5 1 部 分 析法 关 注 的是 局 部 参数 对 模 型 输 出 的 影 。局
响 , 方 法分 析 的是特 定 参 数 ( 值 、 省 值 或 者 最优 该 均 缺 值 ) 化 时 的影 响 , 之 相反 , 变 与 全局 分 析 法 是 同时 分析 所 有参 数对 模 型输 出 的影响 。
时, 参数 敏 感性 分 析 可 以帮助 理 解模 型 结 构 , 乃至 发现 模 型 结构 缺 陷 , 而改善 模 型 结构 。本 文 详 细归 从
纳 了典 型 的 参数 敏 感 性 分析 方 法及 其 求 解 过程 , 论 了各 个 方 法 的利 弊 , 对参 数 敏 感 性 分析 方 法在 讨 并
敏感 性通 常是 通 过 计 算 与参 数 相 关 的输 出 函 数 的 偏
导数 获得 。参 数敏 感 度 S的计 算 方 法见 公 式 ( )S可 1, 以通 过 在保 持 其 他参 数 不 变 的情 况 下 .改 变 △ e值 得
到嘲 :
丝(!: 垒 }:旦2丝(3: : 堡:2± 竺! : 二 璺:! : :曼 : : 2 :
局 部分 析法 是对 多维 参数 空 间 中的某 一点 进 行参 数 敏感性 分 析 , 点 通过 缺省 值或 人 工粗 略率 定赋 值 , 该
收 稿 日期 : 0 9 1 — 5 20 — 0 2 基金项 目: 国家 自然 科 学 基 金 创 新 研 究 群 体 科 学 基 金 项 目 (0 2 0 6 ; 十 一 五 ” 家 科 技 支 撑 计 划 (0 6 A o A 72o B 2 B 8 ; 5 7 10 ) “ 国 2 o B B 4 o ,0 8 AB 9 0 ) 国家 重 点 基 础 研 究 发 展 计 划 (7 9 3计 划 )20 C 4 30 (0 6 B 0 4 4) 作者简介 : 田雨 (9 5 )女 , 津 人 , 士研 究 生 , 18一 , 天 博 主要 研 究 方 向 为水 文 水 资 源 。
1 0
水 文
第 3卷 0
(: (+ e ∑ ( ) ) ee i m)
il U =
( 2 )
计算 所求 参数 的统计 特征 , 最后 给 出所 求解 的近 似值 , 解 的精确 度可用 估计 值 的标准误 差来 表示 。 ( ) 计分 析法 。Mo t C r 2统 ne al o法通 常采 用两 种方
基于SWAT模型的水文模拟参数区域差异研究
基于SWAT模型的水文模拟参数区域差异研究刘倩【摘要】基于SWAT模型对区域水文参数值及参数敏感区域差异进行研究,有利于减少水文过程模拟的不确定性,从而提高SWAT模型在不同区域水文过程参数区域差异模拟中的应用水平.基于此,文章主要针对SWAT模型,对该研究区流域进行合理划分,采用定性分析法对该研究区水文过程参数差异及空间分布敏感性和相关影响因素进行分析.研究结果表明,土壤可利用量及径流曲线系数和研究区地表水文过程,在各个不同区域都较为敏感.在此研究过程中发现,由于研究区植被及土壤和气候等条件存在差异,因此该流域水文过程模拟的参数存在一定的区别.故各地区在社会实践过程中,应充分考虑不同气象因素及地理因素,尽量减少不同气象因素对人类生产实践活动的影响,结合区域水文过程模拟差异,因地制宜展开业生产实践.%The research of SWAT model on regional hydrological parameters and the sensitivity of parameters based on regional differences helps to reduce the uncertainty of hydrological processes,so as to improve the level of application of SWAT model in different areas of regional differences in parameters in the simulation of hydrological process.Based on this,this paper focuses on the SWAT model,a reasonable division of the basin in the study area,the distribution of sensitivity and related differences in the study area hydrological process parameters and spatial factors were analyzed using the method of qualitative analysis.The results show that the soil available capacity and runoff curve coefficient and surface hydrological process in the study area are sensitive to different regions.In the course of the study,it is found that there are some differences between thevegetation and soil and climate,so there are some differences in the parameters of hydrological process simulation.Therefore,all regions in the process of social practice,should fully consider the different meteorological factors and geographic factors,minimize the impact of different meteorological factors on the production of human practice,combined with the regional hydrological process simulation to show differences,business practice.【期刊名称】《黑龙江水利科技》【年(卷),期】2016(044)012【总页数】4页(P8-11)【关键词】SWAT模型;水文模拟参数;区域差异;研究【作者】刘倩【作者单位】辽宁省锦州水文局,辽宁锦州 121000【正文语种】中文【中图分类】P333随着自然科学和地质科学不断发展进步,在对人类社会活动及气候变化、流域内产汇流、水资源时空演变等自然和社会规律进行研究时,通常采用分布式水文模型进行分析模拟。
SWMM敏感性参数的分析方法探讨
做 了很多研究 , 但仍 没有得出统一的结 论。为进一步探究 S WMM 的敏感参数 , 本文首先介绍 了国内采用 的敏 感性
分析方法 ; 再者 , 根据 S WMM 的特点、 结合其它领域应 用的理论 , 提 出 了两 种适用 于 S WMM 的分析 方法 ; 最后 , 对 当前研究 中需要注意 的问题和研究方 向给出了建议 。
S WMM 敏 感 性 参 数 的分 析 方 法 张 敏
( 昆明t l- r 大学建筑工程学院, 昆明 6 5 0 0 5 1 )
【 摘
要J准确分析 S WM M的敏感 参数 , 能极大 的提 高城市洪 涝灾 害的模拟 效率与 准确性 。国内学者 虽然
第3 9卷 第 1期 2 0 1 7年 1月
低 温 建 筑 技 术
L 0W TE MP ERATURE ARC Hr r E CT URE TE CHN0L OGY 1 23
建筑科学管理
D O I : 1 0 . 1 3 9 0 5 / j . c n k i . d w j z . 2 0 1 7 . 0 1 . 0 4 3
【 关键词 】S WM M; 敏感参数 ; 全局分析法 ; 局部分析法 【 中图分类号 】T U 8 2 3 . 3 【 文献标 识码 】A 【 文章编号 】1 0 0 1 — 6 8 6 4 ( 2 0 1 7 ) 0 1 — 0 1 2 3— 0 3
DI S CUS S I ON 0N THE M ETHoD oF ANALYZ I NG THE
S c i e n c e a n d T e c h n o l o g y , K u n m i n g 6 5 0 0 5 1 , C h i n a )
SWAT模型参数率定和验证
泥沙校核
? 调整 HRU/ 子流域负荷:
?调整USLE 方程中的作物管理因子USLE_P (.sub 文件 中)
?调整USLE 方程中的坡长因子SLSUBBSN (.sub 文件中 )
?调整HRU 的坡度(.sub 文件中) ?调整土地利用作物因子USLE_C ( crop.dat 文件中) ?调整耕作措施(*.mgt 文件中 )、作物残留系数
?调整P渗透系数PPERCO (.bsn) ?调整土壤 P比例系数 PHOSKD( .bsn)
?河道内可溶性P运移过程校核
?调整藻类生物量比率 AI2(.wwq)
水质校核
?有机P校核 ?子流域内有机P校核
?调整土壤中营养物的初始浓度 SOL_ORGP ( .sol)
?确定并调整土壤表层的施肥率 FRT_LY1 ( .mgt)
模型率定
?分布真实式世界水质模型误观测差-模源型
误差
模型率定
?参数自动率定结果的优劣指标
?Nash-Suttcliffe 效率系数 ?相关系数 ?相对误差 ?模拟和实测数据的平均误差和标准差 ?斜率、截距、回归系数等
?比较目标:
利用界面进行参数率定
水量校核
?检查水量平衡
?sub.dbf文件中的PET、ET、SNOMELT 、SW、 GWQ、SURQ、WLYD 、 rch.dbf 文件中的 FLOW_IN、FLOW_OUT
RSDCO 和作物混合有效系数BIOMIX (. bsn) ?调整N渗透系数NPERCO (.bsn)
?河道内矿物性N运移过程校核
?调整藻类生物量比率AI1(.wwq)
流体动力学模型的参数优化与敏感性分析方法
流体动力学模型的参数优化与敏感性分析方法引言流体动力学模型是研究和预测流体行为的有效工具。
然而,模型的准确性和可靠性取决于参数的选择。
为了优化模型并评估参数对模型输出的影响,研究者们开发了各种参数优化和敏感性分析方法。
本文将介绍流体动力学模型的常见参数优化和敏感性分析方法,并讨论它们的优缺点及应用领域。
参数优化方法1. 遗传算法(Genetic Algorithm)遗传算法是一种模拟生物进化过程的优化算法。
它通过选择、交叉和变异操作来生成新的参数组合,并根据适应度函数评估它们的优劣。
适应度函数通常是模型输出与实测数据之间的均方根误差或相关系数。
遗传算法可以帮助研究者在搜索参数空间时快速找到最优解。
2. 粒子群优化算法(Particle Swarm Optimization)粒子群优化算法是另一种基于群体智能的优化算法。
它模拟了鸟群或鱼群中个体之间的协作和信息传递。
算法中的每个粒子代表一个参数组合,并根据其自身的历史最优解和群体的最优解进行更新。
粒子群优化算法具有全局搜索能力和快速收敛性。
3. 贝叶斯优化算法(Bayesian Optimization)贝叶斯优化算法是一种基于贝叶斯统计的优化算法。
它利用模型输出与实测数据之间的先验知识来指导参数搜索过程。
算法通过更新参数的后验分布来逐步收敛到最优解。
贝叶斯优化算法适用于对高代价函数进行优化的问题,并能够在少量样本点的情况下取得较好的效果。
敏感性分析方法1. 全局敏感性分析(Global Sensitivity Analysis)全局敏感性分析通过评估模型输出与各参数之间的敏感性来确定参数对模型输出的影响程度。
常见的全局敏感性分析方法包括Sobol指数、Morris元分析方法和FAST(Fourier Amplitude Sensitivity Test)等。
全局敏感性分析可以帮助研究者识别对模型输出影响最大的参数,并优先进行优化。
2. 局部敏感性分析(Local Sensitivity Analysis)局部敏感性分析旨在评估模型在参数变化时的灵敏程度。
SWAT模型参数率定和验证
12
Sol_z
4 BIOMIX
10
SOL_ORG N
自动校核
?单击Tools菜单下的auto-calibration and uncertainty,弹出下图所示对话框
自动校核
?根据敏感性分析结果选择参数 ?选择子流域14进行自动校核
13 7
1
Out 4 3 6 1 5 7 11 4
9
Out 5 4 9 1 2 7 6 23
5
Out 6 4 5 2 3 6 10 8
15
Out 7 3 6 2 5 7 10 4
1
Out 8 4 5 1 6 7 10 3
2
Out 9 参数
5 ALPHA_BF
6
ESCO
1
CN2
3 SOL_AWC
7
surlag
?径流:地表径流和地下径流
?调整CN2(.mgt) ?若地表径流仍然不符合要求,则调整SOL_AWC
(.sol)或ESCO(.bsn or .hru) ?基流值太高:增加 GW_REVAP(.gw) ,减小
REVAPMN(.gw) ,增加 GWQMN(.gw),
?蒸发:土壤蒸发、植被蒸发
?EPCO、EPSO (.bsn or .hru文件中)
模型率定
?分布真实式世界水质模型误观测差-模源型
误差
模型率定
?参数自动率定结果的优劣指标
?Nash-Suttcliffe 效率系数 ?相关系数 ?相对误差 ?模拟和实测数据的平均误差和标准差 ?斜率、截距、回归系数等
?比较目标:
利用界面进行参数率定
水量校核
?检查水量平衡
?sub.dbf文件中的PET、ET、SNOMELT 、SW、 GWQ、SURQ、WLYD 、 rch.dbf 文件中的 FLOW_IN、FLOW_OUT
基于USLE和GIS的水土流失敏感性空间分析以河北太行山区为例
基于USLE和GIS的水土流失敏感性空间分析以河北太行山区为例一、本文概述本文旨在利用通用土壤流失方程(USLE)和地理信息系统(GIS)工具,对河北太行山区的水土流失敏感性进行空间分析。
通过结合这两种强大的工具,我们希望能够更准确地评估该地区的水土流失风险,为土地管理、环境保护和可持续发展提供决策支持。
河北太行山区是我国重要的生态屏障,但由于地形复杂、气候多变、人类活动频繁,水土流失问题日益严重。
因此,开展水土流失敏感性空间分析对于该地区的生态保护和可持续发展具有重要意义。
本文首先介绍了USLE的基本原理和计算过程,包括降雨侵蚀力因子、土壤可侵蚀性因子、坡长坡度因子和植被覆盖与管理因子的确定方法。
然后,我们详细阐述了如何利用GIS技术进行空间数据处理和分析,包括地形数据的提取、土壤数据的分类、降雨数据的空间插值等。
在方法部分,我们将详细介绍如何结合USLE和GIS进行水土流失敏感性空间分析的具体步骤。
包括建立水土流失敏感性评价模型、确定评价因子权重、计算敏感性指数等。
我们将以河北太行山区为例,展示水土流失敏感性空间分析的结果,并分析其空间分布特征和影响因素。
通过本文的研究,我们期望能够为该地区的水土流失防治工作提供科学依据和技术支持。
二、研究方法和数据来源本研究采用修正的通用土壤流失方程(USLE)与地理信息系统(GIS)相结合的方法,对河北太行山区的水土流失敏感性进行空间分析。
USLE作为一种广泛使用的土壤侵蚀预测模型,其通过综合考虑降雨侵蚀力、土壤可蚀性、坡长坡度、植被覆盖及管理措施等因素,为水土流失的定量评估提供了有效工具。
结合GIS的空间分析功能,我们能够更加精准地识别和评估水土流失敏感区域,为区域水土保持和生态环境建设提供科学依据。
在具体操作中,我们首先利用GIS平台,对河北太行山区的地形地貌、土壤类型、植被分布等基础数据进行整合与预处理。
然后,基于USLE模型,构建水土流失敏感性评估模型,将各项影响因子进行量化并整合到GIS平台中。
基于改进蝴蝶优化算法的新安江模型参数率定研究
基于改进蝴蝶优化算法的新安江模型参数率定研究目录一、内容概览...............................................21.1 研究背景及意义.........................................2 1.2 国内外研究现状分析.....................................3 1.3 新安江模型介绍.........................................41.4 蝴蝶优化算法概述.......................................5二、理论基础与预备知识.....................................62.1 新安江模型的数学描述...................................6 2.2 蝴蝶优化算法原理解析...................................82.3 参数率定理论框架.......................................9三、改进蝴蝶优化算法设计..................................103.1 IBFOA算法流程设计.....................................11 3.2 IBFOA算法参数设置.....................................123.3 IBFOA算法性能评估指标.................................13四、新安江模型参数率定实验设计与实施......................144.1 实验环境搭建..........................................15 4.2 实验数据准备..........................................16 4.3 参数率定实验步骤......................................16五、新安江模型参数率定结果分析............................185.1 实验结果展示..........................................195.2 参数敏感性分析........................................205.3 参数最优解确定........................................21六、改进策略与方法........................................226.1 对比分析IBFOA与传统BFOA...............................236.2 提出新的参数调整策略..................................246.3 算法改进效果验证......................................26七、结论与展望............................................277.1 研究结论总结..........................................277.2 研究局限性与不足......................................287.3 未来研究方向建议......................................29一、内容概览1.研究背景与意义新安江模型作为水文预报和水资源管理中的关键工具,其参数率定的准确性直接影响到预报结果的可靠性。
swat模型
• 大型流域模拟过程中,一般采用流域-子流域-水文响应单元的空间离散 方法:
以分水岭上的分水线为界进行分割,可以把一个流域分成若干个子流域,每个子流域内 部有相应的河道,这些河道又将一个个的子流域连成一个整体。 当流域面积较大时,在划分出来的子流域内部,依然分布着多种土壤类型和多种土地利用 方式。为了反映子流域内部不同的土地利用和不同的土壤类型引起的蒸散发、表面径流、 入渗水、农业管理措施等其他水文条件和人类经济活动引起的差异性,可以在子流域内部 进一步划分水文响应单元。
• 优点:确保所有参数在其取值范围内均被采样,并且明确地确定哪一个 参数改变了模型的输出,减少了需要调整的参数数目,提高了计算效率 。
• 原理: LH-OAT先执行LH采样,然后执行OAT采样,见下图。
模型参数敏感性分析
• 首先,每个参数被划分为N个区间,在每个区间内取一个采样点(LH采 样)。然后,一次改变一个采样点(OAT)。
• 水文循环是最主要的流域过程,也是流域内部过程的主要驱动力, 包括:降雨、植被截留、蒸散发、地表积水、入渗、坡面流及河道 水流、土壤水侧向运动(壤中流)、深层渗水等。
SWAT模型 的流域水文
过程
水循环的 陆面部分 控制每个子流域内主河道的水、沙、
(产流、坡面汇流) 营养物质和化学物质等的输入量
水循环的 水面部分
(河道汇流)
决定水、沙等物质从河网向流域 出口的输移运动
SWAT模型的流域地理过程
•陆相水文循环
模拟过程:气候模拟(降水、融雪、土壤温度等)、水文模拟(植物截 留、入渗、壤中流、蒸散发、表面径流、地下径流等)、作物/植 被生长模拟、土壤侵蚀与产沙模拟、养分(氮、磷)输移模拟、杀虫 剂模拟、农业管理操作模拟。对应于SWAT模型的几个功能模块。
水文模型心得
水文模拟的相关认识水文模拟是通过水文模型对自然系统的动态描述,进而揭示复杂的水文现象实质。
水文模型可以分为随机模型、确定模型、随机—耦合模型;也可以分为集总式模型、分布式模型;还可以分为物理模型、数学模型等。
现今,数学模型是目前最简单、应用最为广泛、发展速度最快的模型。
水文模拟可追溯到1858年Mulvany建立的推理公式方法、1932年Sherman的单位线、1933年Horton的入渗方程、1948年Penman 的蒸发公式等,这些标志着水循环各种过程模拟的进展。
随后,水文模拟经历了三次“革命”。
第一次革命发生在19世纪60年代,当时的计算机数字革命给予了水文模拟更加方便使用。
第一类是研制的用于洪水预报的概念性模型,如SAC模型、日本的水箱模型、中国的新安江模型、瑞典的包夫顿模型。
第二类是基于物理基础的分布式水文模型。
Freeze和Harlan提出了分布式水文模型的构造框架。
第二次革命是建立在卫星和雷达的是大量信息的提供,水文模拟所需的数据精度大大提高,并且逐渐迈向实用领域。
21世纪,第三次革命的到来。
随着人们需求的不断增长,水文模拟不单单只是研究产汇流机制,也注重水文循环和大气圈、植被、土地利用等之间的相互作用。
此时,GIS和数据库技术向水文模拟注入了新的血液。
水文模型是水循环研究、水文机理研究、水资源管理、水环境保护的基本工具。
但是在我们建立水文模型和水文模拟过程中存在很多问题:1.模型问题。
如今我们大多数情况下使用的SWAT模型、MIKE模型、AGNPS模型、Stanford模型等,都是从国外引进的。
尽管在我国某些区域,相应的模型与当地实际情况拟合比较好,但是由于这些模型内部结构是以国外的环境情况为基础做的。
所以,我们需要从模型内部结构、其方案、算法上出发建立出适合我国国情的模型2.数据瓶颈。
无论是哪种模型,其对数据要求精度及数量特别高,一般涉及水文数据、土地利用、土壤类型、地形、植被等。
常见水文模型参数率定
3/8/2011
三、常用洪水预报模型
降雨径流相关图法
P—Pa—R关系曲线
Pa:
前期雨量计算,也称前期影
响雨量,是反映土壤湿度的参
数。
通常从主汛期开始计算。
P:
把时段雨量序列变成累积雨
量序列。
R:
由累积净雨转化成时段净雨
3/8/2011
二、模型参数的性质与约值
三水源蓄满产流模型参数
K: 蒸发皿实测的蒸发量乘上K就是流域蒸散发能力。 K的取值控制总的水量平衡,是非常重要的参数。 K的参数率定通常将目标函数固定为多年水量平衡再计算而得。
3/8/2011
二、模型参数的性质与约值
三水源蓄满产流模型参数
B:
7 //P~R曲线条数N//
0 10 20 30 40 50 60 //每条P—R曲线对应的Pa值//
10 //每条P—R曲线的节点数M//
1
0
0
0
0
0
0
0
0
2
20
3.7
5
6.3
7.6
9
12
15
。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。
3/8/2011
10
300 205 215 225 235 245 255 265
常用洪水预报模型
三水源蓄满产流模型 三水源滞后演算模型 马斯京根河道分段连 X 续演算模型
MP
3/8/2011
一、常用洪水模型概述
三水源蓄满产流模型参数
WM(WUM、WLM、WDM):流域平均蓄水容量,反映干旱程度的指标。
