新安江模型研究(完整资料).doc
- 格式:doc
- 大小:217.19 KB
- 文档页数:5
东南大学交通学院桥涵水文资料整理指导老师:许崇法姓名:郭赵元学号:21710131目录第一章新安江模型 (3)1.1 新安江模型简介 (3)1.2 新安江模型的基本原理 (3)1.3 新安江模型结构 (4)第二章陕北模型 (6)2.1陕北模型简介 (6)2.2 陕北模型结构 ............... .. (7)2.3 模型评述 (8)第一章新安江模型1.1新安江模型简介新安江模型始建于 1973 年,采用蓄满产流的概念,以土壤含水量达到田间持水量后才产流,是个分布式的概念性模型,30 多年来在我国湿润与半湿润地区有广泛应用,并发展改进为三水源的以及其他多水源的模型。
原华东水利学院的赵人俊教授于1963年初次提出湿润地区以蓄满产流为主的观点,主要根据是次洪的降雨径流关系与雨强无关,而只有用蓄满产流概念才能解释这一现象。
上个世纪70年代国外对产流问题展开了理论研究,最有代表性的著作是1978年出版的《山坡水文学》,它的结论与赵人俊先生的观点基本一致:传统的超渗产流概念只适用于干旱地区,而在湿润地区,地面径流的机制是饱和坡面流,壤中流的作用很明显。
20世纪70年代初建立的新安江模型采用蓄满概念是正确的。
但对于湿润地区,由于没有划出壤中流,导致汇流的非线性程度偏高,效果不好。
80年代初引进吸收了山坡水文学的概念,提出三水源的新安江模型。
1.2新安江模型的基本原理新安江模型是分散性模型,可用于湿润地区与半湿润地区的湿润季节。
当流域面积较小时,新安江模型采用集总模型,当面积较大时,采用分块模型。
它把全流域分为许多块单元流域,对每个单元流域作产汇流计算,得出单元流域的出口流量过程。
再进行出口以下的河道洪水演算,求得流域出口的流量过程。
把每个单元流域的出流过程相加,就求得了流域的总出流过程。
该模型按照三层蒸散发模式计算流域蒸散发,按蓄满产流概念计算降雨产生的总径流量,采用流域蓄水曲线考虑下垫面不均匀对产流面积变化的影响。
第二章新安江流域水文模型60年代初,河海大学(原华东水利学院)水文系赵人授等开始研究蓄满产流模型,配合一定的汇流计算,将模型应用于水文预报和水文设计。
1973年,他们在对新安江水库做人库流量预报的工作中,把他们的经验归纳成一个完整的降雨径流流域模型——新安江模型。
模型可用于湿润地区和半湿润地区的湿润季节径流模拟和计算。
最初的新安江模型为两水源模型,只能模拟地表径流和地下径流。
80年代初期,模型研制者将萨克拉门托模型与水箱模型中,用线性水库函数划分水源的概念引入新安江模型,提出了三水源新安江模型,模型可以模拟地面径流、壤中流、地下径流。
1984至1986年,又提出了四水源新安江模型,可以模拟地面径流、壤中流、快速地下径流和慢速地下径流。
三水源新安江模型一般应用效果较好,但模拟地下水丰富地区的日径流过程精度不够理想。
在新安江三模型中增加慢速地下水结构就成为四水源新安江模型。
当流域面积较小时,新安江模型采用集总模型,当面积较大时,采用分块模型。
分块模型把流域分成许多块单元流域,对每个单元流域做产、汇计算,得到单元流域的出口流量过程。
再进行出口以下的河道洪水演算,求得流域出口的流量过程。
把每个单元流域的出流过程相加,就求得了流域出口的总出流过程。
划分单元流域的主要目的是处理降雨分布的不均匀性,因此单元流域应当大小适当,使得每块面积上的降雨分布比较均匀.并有一定数目的雨量站。
其次尽可能使单元流域与自然流域相一致,以便于分析与处理问题,并便于利用已有的小流域水文资料。
如果流域内有大中型水库,则水库以上的集水面积即应作为一个单元流域。
因为各单元流域的产汇、流计算方法基本相同,以下只讨论一个单元流域的情况。
2.1新安江两水源模型1.模型结构和参数新安江两水源模型的产流子模型采用蓄满产流模型,蒸发计算采用三层蒸发计算模型。
利用稳定下渗率FC将径流划分为地面径流和地下径流两种水源。
地面径流采用单位线汇流,地下径流采用一次线性水库汇流。
一、模型的结构与参数三水源新安江模型的流程图如图1所示。
图1 三水源新安江模型流程图图1 中输入为实测雨量P ,实测水面蒸发EM ;输出为流域出口流量Q ,流域蒸散发E 。
方框内是状态变量,方框外是参数变量。
模型结构及计算方法可分为以下四大部分。
1. 蒸散发计算用三个土层的模型,其参数为上层张力水容量UM ,下层张力水容量LM ,深层蒸散发系数C ,蒸散发折算系数K ,所用公式如下:当上层张力水蓄量足够时,上层蒸散发EU 为EM E EU ⨯=当上层已干,而下层蓄量足够时,下层蒸散发EL 为LM WL EM K EL /⨯⨯=当下层蓄量亦不足,要触及深层时,蒸散发ED 为EM K C ED ⨯⨯=2. 产流量计算据蓄满产流概念,参数为包气带张力水容量WM ,张力水蓄水容量曲线的方次B ,不透水面积的比值IM ,所用公式为)1/()1(IM B WM WM -+⨯=))/1(1()1/(1B WM W MM A +--=当0≤⨯-EM K P ,则R=0不然,则当MM A EM K P <+⨯-,B MM A EM K P WM W WM EM K P R ++⨯--⨯++-⨯-=1)/)(1(不然,则W WM EM K P R +-⨯-=式中 R ——产流量;MM ——流域最大点蓄水容量。
3. 分水源计算分三种水源,即地面径流RS 、地下径流RG 和壤中流RI 。
参数为表层土自由水蓄水容量SM ,表层自由水蓄水容量曲线的方次EX ,表层自由水蓄量对地下水的出流系数KG 及对壤中流的出流系数KI ,所用公式为SM EX MS ⨯+=)1())/1(1()1/(1EX SM S MS AU +--⨯=)/())((EM K P EM K P IM R FR ⨯-⨯-⨯-=FR KG S RG ⨯⨯=FR KI S RI ⨯⨯=当 0,0=≤⨯-RS EM K P不然,当MS AU EM K P <+⨯-,则FR MS AU EM K P SM S SM EM K P RS EX ⨯+⨯--⨯++-⨯-=+))/)(1((1 当MS AU EM K P ≥+⨯-,则FR SM S EM K P RS ⨯-+⨯-=)(4. 汇流计算地下径流用线性水库模拟,其消退系数为CG ,出流进入河网。
8新安江模型第八章新安江模型8.1 概述新安江模型是由原华东水利学院(现为河海大学)赵人俊教授等(赵人俊,1984)提出来的。
从降雨径流经验相关图研究开始(华东水利学院水文系,1962),投入了水文预报教研室的十余位教师、研究生和上百的本科生前后经历了约20年才形成了蓄满产流概念、理论及其二水源新安江模型。
之后提出三水源新安江模型(赵人俊,1984),并开始在水情预报和遥测自动化的实时洪水预报系统中开始大量应用,通过对模型的结构、考虑的因素不断改进和完善,发展至今已形成了理论上具有一定系统性、结构较为完善、应用效果较好的流域水文模型,并被联合国教科文组织列为国际推广模型而广为国内外水文学家所了解和应用。
新安江模型研究概括起来可以分为二水源新安江模型、三水源新安江模型和新安江模型改进研究三个阶段。
8.2 二水源新安江模型二水源新安江模型包括直接径流和地下径流,产流计算用蓄满产流方法,流域蒸发采用二层或三层蒸发,水源划分用的是稳定下渗法,直接径流坡面汇流用单位线法,地下径流坡面汇流用线性水库,河道汇流采用马斯京根分河段演算法。
8.2.1 前期研究降雨径流相关图是径流估计最早使用的方法之一。
考虑前期气候指数的降雨径流相关图是蓄满产流概念形成的基础,见图8-1。
图中为降雨量,为径流深, P为前期气候指PRa,0数。
在实际应用中,要计算一次降雨所产生的洪水径流总量,为配合汇流计算,还需求出逐时段的净雨量。
利用上述相关图推求时段净雨量的具体步骤如下。
P(1)求本次降雨开始时的; a,0(2)按逐时段累积降雨量在关系图上查得累积径流量;图8-1 时段净雨量推求(3)由相邻时段的累积径流量之差得时段净雨量。
在这相关图应用过程中发现两个问题,一是前期气候指数不是一个物理量,二是关系不满足水量平衡方程。
为此,提出由土壤含水量来反应前期气候的W 湿润情况,点关系图,经大量的实践发RfPW,(,)现,在湿润地区曲线簇的上段均接近45?直线,若W点绘成与关系(是扣除雨期蒸发后的PEW,RPE净雨量),则呈现如图8-2所示的关系。
新安江模型
新安江模型的流程图
图中输入为实测降雨P和实测蒸散发能力EM,输出为流域出口断面流量Q和流域蒸散发量E 。
方框内是状态变量,方框外是常数常量。
模型主要由四部分组成,即蒸散发计算、产流计算、水源划分和汇流计算。
变量符号说明:
P,实测降雨(已知)
E,流域蒸散发量(第一步计算已知)
PE,降雨减去蒸发量(产流计算的输入)
R,流域产流量(产流计算的输出)
W,流域时段末土壤平均蓄水量,产流计算的输出
KC,蒸发折算系数
C,深层蒸散发折算系数
WM,流域总蓄水容量
B,抛物线指数
1蒸发量计算
各层蒸散发的计算原则是,上层按蒸散发能力蒸发,上层含水量蒸发量不够蒸发时,剩余蒸散发能力从下层蒸发,下层蒸发与蒸散发能力及下层含水量成正比,与下层蓄水容量成反比。
要求计算的下层蒸发量与剩余蒸散发能力之比不小于深层蒸散发系数C。
否则,不足部分由下层含水量补给,当下层水量不够补给时,用深层含水量补。
输入:
(1)KC ,蒸发折算系数
(2)上、下、深三层的蓄水容量WUM 、WLM 和WDM
(WM=WUM+WLM+WDM )
(3)C ,深层蒸散发折算系数
输出:
(1)上、下、深各层的流域蒸散发量EU 、EL 和ED (E=EU+EL+ED ) -- 时段值
(2)上、下、深三层的含水量WU 、WL 和WD (W=WU+WL+WD ) -- 时段值
注:以上的WM 、E 、W 分别表示总的流域蓄水容量、蒸散发量、土壤含水量。
计算公式:
EP=K×E 0(E 0,实测水面蒸发量,划分为时段值),WU 0日模拟合而来。
一阶段:当 P+WU ≧EP 时,EU=EP 、EL=0、ED=0,此时WU(t+1)=P(t)+WU(t)-EU(t) 二阶段:当 P+WU<EP 时, EU=P+WU
若WL>C×WLM 则 EL=(EP -EU)×WL/WLM 、ED=0
若WL<C×WLM 且 WL ≧C×(EP -EU) 则
EL=C×(EP -EU),ED=0
三阶段:若WL<C×WLM 且 WL<C×(EP -EU) 则
EL=WL ,ED=C×(EP -EU)-WL
E = EU + EL + ED
PE = P - E
注:一阶段上层够蒸;二阶段上层不够下层够蒸;三阶段上下层都不够蒸。
2产流计算
土壤蓄满之前不产流,当土壤湿度饱和时,所有降雨(减去同期蒸发量)都产流。
输入:
PE ,去除蒸发量后的降雨;
f/F ,f 为流域内包气带蓄水容量小于或等于WM’的面积,F 为全流域面积;
W’mm ,流域最大的点蓄水容量;
W’M ,流域内某一点的蓄水容量;
B ,反映流域包气带蓄水容量分布的不均匀性,常数,越小越均匀;B=0是均匀不变 W 0,流域初始土壤含水量
输出:
R ,降雨产生的径流
计算公式:
一般说来,流域内各点的蓄水容量并不相同,新安江三水源模型把流域内各点的蓄水容量概化成一条抛物线,即
B mm
m W W F f )1(1''--= W’mm 即流域最大的点蓄水容量;W’M 为流域内某一点的蓄水容量;f 为蓄水容量≤W’M 值时的流域面积;F 为流域面积;B 为抛物线指数。
根据上面公式可求得流域平均蓄水容量为 ⎰'+'='-=mm
W mm m B W W F f WM 01
d )1( 与流域初始平均蓄水量0W 相应的纵坐标(A )为])1(1[110+--'=B mm WM
W W A 当0>-=E P PE 时,则产流;否则不产流。
产流时,
当mm
W A PE '<+,则B mm
W A PE WM W WM PE R +'+-++-=10)1( 当mm
W A PE '≥+,则)(0W WM PE R --= 3水源划分
产流量R 进入自由水水库内,分为地面径流RS ,壤中流RSS 和地下径流RG 。
输入:
R ,时段产流量;(产流计算而来)
PE ,去除蒸发量后的降雨;(蒸发量计算而来)
FR0,上一时段的产流面积;
S0,产流面积上的平均自由水深;
SM ,流域的平均自由水容量;
EX ,流域自由水蓄水容量曲线的指数;
KSS ,自由水蓄水容量对壤中流的出流系统;
KG ,自由水蓄水容量对地下水的出流系统。
中间过程变量:
SMMF ,产流面积FR 上最大一点的自由水蓄水容量;
输出:
RS ,时段地面径流;
RSS ,时段壤中流;
RG ,时段地下径流;
S ,产流面积上的平均自由水深。
计算公式:
PE R FR /=
FR FR S t S /00)(⋅=
)1(EX SM SMM +⋅=
])1(1[/1EX FR SMM SMMF --⋅=
)1/(EX SMMF SMF +=
])1(1[11
EX SMF
S SMMF AU +--=
当SMMF AU PE ≥+时,则
))1((SMF t S PE FR RS --+=;FR KSS SMF RSS ⋅⋅=
FR KG SMF RG ⋅⋅=; FR RG RSS SMF t S /)()(+-=
当SSMF AU PE <+<0时,则 ])1()1([1++-
+-+-⋅=EX SMMF AU PE SMF t S SMF PE FR RS )/)1((FR RS t S PE FR KSS RSS --+⋅=
)/)1((FR RS t S PE FR KG RG --+⋅=;
FR RG RSS RS PE t S t S /)()1()(++-+-=
4汇流计算
汇流分为三个阶段进行:坡地汇流阶段、河网汇流阶段和河道汇流。
经过水源划分得到地面径流直接进入河网,成为地面径流对河网的总入流(QS ),壤中流和地下径流分别有消退系数CI 、CG ,分别形成壤中流、地下水对河网的总入流QI 、QG 。
时段河网总入流就是三者时段值之和。
河网汇流计算采用滞后演算法。
输入:
RS ,时段地面径流;
RSS ,时段壤中流(RI )
U ,单位转换系数,可将径流深转换成流量(F 为流域面积)t F U ∆=
6.3
CI ,深层水消退系数;
CG ,地下径流消退系数;
CS ,河网水流消退系数;
L ,滞后时间。
输出:
Q ,河网汇流结果,流量过程。
计算公式:
(1)地面径流的坡地汇流时间不计,直接进入河网: U t RS t QS ⋅=)()(
(2)表层自由水以KI 侧向出流后成为壤中流,进入河网,但如果土层较厚,表层自由水尚可渗入深层土,经过深层土的调蓄作用,才进入河网。
深层自由水用线性水库模拟,其消退系数为CI ,计算公式为
U t RSS CI t QI CI t QI ⨯⨯-+-⨯=)()1()1()(
(3)地下径流的坡地汇流用线性水库模拟,其消退系数为CG ,出流进入河网。
表层自由水以KG 向下出流后,再向地下水库汇流的时间不另计,包括在CG 之内,计算公式为
U t RG CG t QG CG t QG ⨯⨯-+-⨯=)()1()1()(
(4)单元面积的河网汇流是指水流有坡面进入河槽后,继续沿河网的汇集过程。
在河网汇流阶段,汇流特性受制于河槽水力学条件,各种水源是一致的。
河网汇流采用滞后演算法,其计算公式为
)()()()(t QG t QI t QS t QT ++=
)()1()1()(L t QT CS t Q CS t Q -⨯-+-⨯=。