新安江模型报告
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
202蓄1/2水/11能力小于某个点蓄水能力的面积占12总面积之比
L新O安G江O模型
RunModle 运行新安江模型
将该曲线表示为蓄水量公式为
f 1(1 w'm )B
F
w'mm
(1)
（1）式将B赋值为1的话,正好反映了线性关系。 B为其他方次时，分别反映不同的概化曲线。
整个流域的平均蓄水容量应 对流域求和平均,即:
4.模型状态量
土壤湿度:W=Wu+Wl+Wd,程序中假定流域经过了长时间降雨 ，数值分别等于各自的蓄水容量。
流域面积：Area=537 km2 每一步长的小时数：DeltaT=2h U=Area/(3.6*DeltaT）一个将径流深转化为流量的转换系数
2021/2/11
5
L新O安G江O模型
Baidu Nhomakorabea
SetParameters设置模型参数
新安LO江G模O型
1
InitModel模型初始化
2 SetParameters设置模型参数
3 RunModle 运行新安江模型
4 SaveResults 结果输出
2021/2/11
1
L新O安G江O模型
概述
1
为对新安江水库 作入库流量预报 而提出的一种降 雨径流流量模型 。
2
主要特色在于认 为湿润地区为蓄 满产流。提出的 流域蓄水容量曲 线为模型核心。
2021/2/11
9
L新O安G江O模型
RunModle 运行新安江模型
▪ 蒸散发计算——算法的实现(PE<0时)
一、WU+PE>0(上层足够蒸发）
EU=K*E EL=ED=0 二、WU+PE<0(上层不足蒸发)
EU=WU程+序P(在降W雨L与=C上*W层L蓄M水处全还部有蒸发） 1.若WL>C一*个（分k*E界-E点U,没)（有下看层懂足!够蒸发） （解释）
(2) W'm
Wm
0
m(1F f )dW 'mW 1'm Bm
增加参数IMP后,形式变成:
2021/2/11
W'mm11IB MW Pm (3)
13
L新O安G江O模型
RunModle 运行新安江模型
▪ 产流计算——算法的实现
当天含水量折算成径流深 A=WMM*[1-(pow(1-W/WM,1/(B+1))]
2021/2/11
11
L新O安G江O模型
RunModle 运行新安江模型
▪ 产流计算——蓄水容量曲线
最大点蓄水容量
由于产流计算要根据蓄满产
流原理,故必须了解蓄水容
量。但流域内各点的蓄水容
量并不相等，故概化为如图
的曲线。
如果连最大点蓄水容量都已不 能继续存蓄,则全流域产流; 反之，则并非全流域产流，就 需要引入抛物线指数B来概化 考虑了。
EL=(K应*当E-有EU相) 应（的剩理余论的解全释由。下层）
ED=0 2.若WL<C*(K*E-U)（下层不够蒸发）
EL=WL （下层按最大能力蒸发） ED=C*(K*E-EU)-EL(正好是不足的部分）
2021/2/11
10
L新O安G江O模型
RunModle 运行新安江模型
▪ 产流计算——蓄满产流原理
3
为分散性模型， 按一定方法分块 。对每个单元流 域进行产汇流计 算，然后相加得 到总出流过程。
2021/2/11
2
L新O安G江O模型
概述
正演程序流程图
2021/2/11
3
L新O安G江O模型
InitModel模型初始化
1 步数:steps
原程序步长为1天,步数为731，即两年。 现依据《常用水文预报算法和计算程序》P21算例，将步长改为2h，时段 数为24。其余参数均按照该算例取值。
▪ 1.流域蒸散发能力与实测水面蒸发之比 K=0.65
▪ 2.流域不透水面积占全流域面积之比 IMP=0.
▪ 3.蓄水容量曲线的方次
B=0.3
▪ 4.上层蓄水容量
Wum=20mm
▪ 5.下层蓄水容量
Wlm=75mm
▪ 6.深层蓄水容量
Wdm=80mm
▪ 7.深层蒸散发系数
C=0.11
▪ 8.地下径流消退系数
2 驱动数据:
降雨量P 蒸散发量E
2021/2/11
4
数据自 《常用水文 预报算法 和计算程序》
算例
L新O安G江O模型
InitModel模型初始化
3.模型输出 蒸散发项:E=Eu+El+Ed 产流项： R=Rg+Rs 出流项： Q=Qrg+Qrs
该程序中产流项和出流项 均只分地表和地下 无壤中流部分
▪ 11.汇流计算参数
Kstor 汇流参数
▪ 12.稳定入渗率（mm/h） Fc=2.3mm/h
关于稳渗率的取值:因该程序采用二水源划分,需要此参数，算例在文 献中为三水源划分，不需要此参数。故从算例中得不到该参数，只能大 致得到一个值，以作参考。
2021/2/11
7
L新O安G江O模型
RunModle 运行新安江模型
▪ 蒸散发计算——三层蒸发模型
上层 下层 深层
按蒸发能力蒸发
上层不足蒸发 补足 按蒸发能力蒸发 下层不足蒸发
补足
2021/2/11
8
L新O安G江O模型
RunModle 运行新安江模型
▪ 蒸散发计算——“下层按蒸发能力”的解释
下层蒸发与剩余蒸散发能力（流域蒸散发量与上层 蒸发量之差）及下层含水量成正比,而与下层蓄水 容量成反比。所谓的下层按蒸发能力蒸发是指: 下层蒸发量与剩余蒸散发能力之比不小于深层蒸散 发系数C。若下层蓄水不足，则只好由深层补足。
KKG=0.99
注:该程序未划分壤中流,故不需要壤中流消退系数 KKSS; 该程序的水源划分为二水源，没用到两个出流系数KSS、KG
2021/2/11
6
L新O安G江O模型
SetParameters设置模型参数
▪ 9.流域平均蓄水容量 WM=WuM+WlM+WdM
▪ 10.流域内最大点蓄水容量 WMM=WM*(1+B)/(1-IMP)
W为流域的初始平均蓄水量（不是容量）,在上图中属于面积（点的累积）方面的量。 在下面的分段中，不等式右边为WMM（最大的点蓄水容量），在上图中
属于纵坐标方面(点）的量。故为了对应，需把WMM对应的纵坐标即A算出。
一、A+PE>=WMM（最大点蓄水容量） R=W+PE-WM
蓄满产流:土壤湿度达到田间持水量之前不产流,所有降雨被 土壤吸收成为张力水。达到田间持水量后，所有降雨都产流。
?为>:（全流域蓄满产流） R=土壤湿度+净降 雨量-土壤蓄水量
总共的水量_?_ 流域内最大点蓄水容量
?为<:（有部分未蓄满产流） R=土壤湿度+净降 雨量-土壤蓄水量+一个考 虑蓄水容量方次的部分
L新O安G江O模型
RunModle 运行新安江模型
将该曲线表示为蓄水量公式为
f 1(1 w'm )B
F
w'mm
(1)
（1）式将B赋值为1的话,正好反映了线性关系。 B为其他方次时，分别反映不同的概化曲线。
整个流域的平均蓄水容量应 对流域求和平均,即:
4.模型状态量
土壤湿度:W=Wu+Wl+Wd,程序中假定流域经过了长时间降雨 ，数值分别等于各自的蓄水容量。
流域面积：Area=537 km2 每一步长的小时数：DeltaT=2h U=Area/(3.6*DeltaT）一个将径流深转化为流量的转换系数
2021/2/11
5
L新O安G江O模型
Baidu Nhomakorabea
SetParameters设置模型参数
新安LO江G模O型
1
InitModel模型初始化
2 SetParameters设置模型参数
3 RunModle 运行新安江模型
4 SaveResults 结果输出
2021/2/11
1
L新O安G江O模型
概述
1
为对新安江水库 作入库流量预报 而提出的一种降 雨径流流量模型 。
2
主要特色在于认 为湿润地区为蓄 满产流。提出的 流域蓄水容量曲 线为模型核心。
2021/2/11
9
L新O安G江O模型
RunModle 运行新安江模型
▪ 蒸散发计算——算法的实现(PE<0时)
一、WU+PE>0(上层足够蒸发）
EU=K*E EL=ED=0 二、WU+PE<0(上层不足蒸发)
EU=WU程+序P(在降W雨L与=C上*W层L蓄M水处全还部有蒸发） 1.若WL>C一*个（分k*E界-E点U,没)（有下看层懂足!够蒸发） （解释）
(2) W'm
Wm
0
m(1F f )dW 'mW 1'm Bm
增加参数IMP后,形式变成:
2021/2/11
W'mm11IB MW Pm (3)
13
L新O安G江O模型
RunModle 运行新安江模型
▪ 产流计算——算法的实现
当天含水量折算成径流深 A=WMM*[1-(pow(1-W/WM,1/(B+1))]
2021/2/11
11
L新O安G江O模型
RunModle 运行新安江模型
▪ 产流计算——蓄水容量曲线
最大点蓄水容量
由于产流计算要根据蓄满产
流原理,故必须了解蓄水容
量。但流域内各点的蓄水容
量并不相等，故概化为如图
的曲线。
如果连最大点蓄水容量都已不 能继续存蓄,则全流域产流; 反之，则并非全流域产流，就 需要引入抛物线指数B来概化 考虑了。
EL=(K应*当E-有EU相) 应（的剩理余论的解全释由。下层）
ED=0 2.若WL<C*(K*E-U)（下层不够蒸发）
EL=WL （下层按最大能力蒸发） ED=C*(K*E-EU)-EL(正好是不足的部分）
2021/2/11
10
L新O安G江O模型
RunModle 运行新安江模型
▪ 产流计算——蓄满产流原理
3
为分散性模型， 按一定方法分块 。对每个单元流 域进行产汇流计 算，然后相加得 到总出流过程。
2021/2/11
2
L新O安G江O模型
概述
正演程序流程图
2021/2/11
3
L新O安G江O模型
InitModel模型初始化
1 步数:steps
原程序步长为1天,步数为731，即两年。 现依据《常用水文预报算法和计算程序》P21算例，将步长改为2h，时段 数为24。其余参数均按照该算例取值。
▪ 1.流域蒸散发能力与实测水面蒸发之比 K=0.65
▪ 2.流域不透水面积占全流域面积之比 IMP=0.
▪ 3.蓄水容量曲线的方次
B=0.3
▪ 4.上层蓄水容量
Wum=20mm
▪ 5.下层蓄水容量
Wlm=75mm
▪ 6.深层蓄水容量
Wdm=80mm
▪ 7.深层蒸散发系数
C=0.11
▪ 8.地下径流消退系数
2 驱动数据:
降雨量P 蒸散发量E
2021/2/11
4
数据自 《常用水文 预报算法 和计算程序》
算例
L新O安G江O模型
InitModel模型初始化
3.模型输出 蒸散发项:E=Eu+El+Ed 产流项： R=Rg+Rs 出流项： Q=Qrg+Qrs
该程序中产流项和出流项 均只分地表和地下 无壤中流部分
▪ 11.汇流计算参数
Kstor 汇流参数
▪ 12.稳定入渗率（mm/h） Fc=2.3mm/h
关于稳渗率的取值:因该程序采用二水源划分,需要此参数，算例在文 献中为三水源划分，不需要此参数。故从算例中得不到该参数，只能大 致得到一个值，以作参考。
2021/2/11
7
L新O安G江O模型
RunModle 运行新安江模型
▪ 蒸散发计算——三层蒸发模型
上层 下层 深层
按蒸发能力蒸发
上层不足蒸发 补足 按蒸发能力蒸发 下层不足蒸发
补足
2021/2/11
8
L新O安G江O模型
RunModle 运行新安江模型
▪ 蒸散发计算——“下层按蒸发能力”的解释
下层蒸发与剩余蒸散发能力（流域蒸散发量与上层 蒸发量之差）及下层含水量成正比,而与下层蓄水 容量成反比。所谓的下层按蒸发能力蒸发是指: 下层蒸发量与剩余蒸散发能力之比不小于深层蒸散 发系数C。若下层蓄水不足，则只好由深层补足。
KKG=0.99
注:该程序未划分壤中流,故不需要壤中流消退系数 KKSS; 该程序的水源划分为二水源，没用到两个出流系数KSS、KG
2021/2/11
6
L新O安G江O模型
SetParameters设置模型参数
▪ 9.流域平均蓄水容量 WM=WuM+WlM+WdM
▪ 10.流域内最大点蓄水容量 WMM=WM*(1+B)/(1-IMP)
W为流域的初始平均蓄水量（不是容量）,在上图中属于面积（点的累积）方面的量。 在下面的分段中，不等式右边为WMM（最大的点蓄水容量），在上图中
属于纵坐标方面(点）的量。故为了对应，需把WMM对应的纵坐标即A算出。
一、A+PE>=WMM（最大点蓄水容量） R=W+PE-WM
蓄满产流:土壤湿度达到田间持水量之前不产流,所有降雨被 土壤吸收成为张力水。达到田间持水量后，所有降雨都产流。
?为>:（全流域蓄满产流） R=土壤湿度+净降 雨量-土壤蓄水量
总共的水量_?_ 流域内最大点蓄水容量
?为<:（有部分未蓄满产流） R=土壤湿度+净降 雨量-土壤蓄水量+一个考 虑蓄水容量方次的部分