新安江模型介绍
- 格式:docx
- 大小:752.45 KB
- 文档页数:5
第二章新安江流域水文模型60年代初,河海大学(原华东水利学院)水文系赵人授等开始研究蓄满产流模型,配合一定的汇流计算,将模型应用于水文预报和水文设计。
1973年,他们在对新安江水库做人库流量预报的工作中,把他们的经验归纳成一个完整的降雨径流流域模型——新安江模型。
模型可用于湿润地区和半湿润地区的湿润季节径流模拟和计算。
最初的新安江模型为两水源模型,只能模拟地表径流和地下径流。
80年代初期,模型研制者将萨克拉门托模型与水箱模型中,用线性水库函数划分水源的概念引入新安江模型,提出了三水源新安江模型,模型可以模拟地面径流、壤中流、地下径流。
1984至1986年,又提出了四水源新安江模型,可以模拟地面径流、壤中流、快速地下径流和慢速地下径流。
三水源新安江模型一般应用效果较好,但模拟地下水丰富地区的日径流过程精度不够理想。
在新安江三模型中增加慢速地下水结构就成为四水源新安江模型。
当流域面积较小时,新安江模型采用集总模型,当面积较大时,采用分块模型。
分块模型把流域分成许多块单元流域,对每个单元流域做产、汇计算,得到单元流域的出口流量过程。
再进行出口以下的河道洪水演算,求得流域出口的流量过程。
把每个单元流域的出流过程相加,就求得了流域出口的总出流过程。
划分单元流域的主要目的是处理降雨分布的不均匀性,因此单元流域应当大小适当,使得每块面积上的降雨分布比较均匀.并有一定数目的雨量站。
其次尽可能使单元流域与自然流域相一致,以便于分析与处理问题,并便于利用已有的小流域水文资料。
如果流域内有大中型水库,则水库以上的集水面积即应作为一个单元流域。
因为各单元流域的产汇、流计算方法基本相同,以下只讨论一个单元流域的情况。
2.1新安江两水源模型1.模型结构和参数新安江两水源模型的产流子模型采用蓄满产流模型,蒸发计算采用三层蒸发计算模型。
利用稳定下渗率FC将径流划分为地面径流和地下径流两种水源。
地面径流采用单位线汇流,地下径流采用一次线性水库汇流。
《流域水文模拟》结课报告新安江模型的原理、结构及应用、发展历程The principle, structure, application and development process of Xin anjiang Model作者姓名:孔旭学科、专业:水文学及水资源学号:21506149指导教师:王国利完成日期:2016年8月30日大连理工大学Dalian University of Technology摘要新安江模型是河海大学提出的一个概念性降雨径流模型,具有原创性,是我国为数不多的被国际上广泛认可的水文模型。
新安江水文模型在我国湿润与半湿润地区广为应用,取得了良好的效果。
经过近50年的发展,新安江模型已经从最初的专门从事水库入库洪水预报的单一功能模型发展为适合用于水文预报、水资源管理、水土资源评价、面源污染预测、气候变化和人类活动影响研究的多功能的水文模型;其部分参数已从靠经验率定发展为可以进行物理推求。
总之,新安江模型是一个不断发展的模型体系。
本文主要由三部分构成。
第一部分为新安江模型简介,回顾了新安江模型产生的历史背景和发展历程,介绍了新安江模型的基本原理和结构体系;第二部分讲述了新安江模型参数的物理意义及其率定;第三部分为新安江水文模型在英那河流域防洪规划编制当中的应用。
关键词:水文模型;新安江模型;洪水预报The principle, structure, application and development process of Xinanjiang ModelAbstractXin anjiang Model originally proposed by Hehai University is a conceptual rainfall runoff model and is also one of the few widely recognized international hydrological model in China. Xin anjiang hydrological model was widely used in our humid and semi-humid areas, and achieved good results.After nearly 50 years study, Xin anjiang model has been developed from the single-function of reservoir flood forecasting into multi-purpose model including hydrological forecasting, water resources management, water and soil resources evaluation, non-point source pollution prediction, climate change and human activities versatile hydrological model studies. And part of its parameters can be acquired through physical calculation instead of experience. In short, Xin anjiang model is an evolving model system.This paper consists of three parts. The first part is about the brief introduction of Xin anjiang model, which recalls the historical background and the development, as well as introduces the basic principles and architecture; the second part describes the physical meaning of Xin anjiang model parameters and calibration; the third part is about the application of Xin anjiang model in Ying Na River Basin flood control planning.Key Words: hydrological model; Xin anjiang model; Flood forecasting目 录摘 要 .............................................................. I Abstract .. (II)引 言 (1)1 新安江模型简介 (2)1.1 新安江模型起源 ............................... 错误!未定义书签。
《流域水文模拟》结课报告新安江模型的原理、结构及应用、发展历程The principle, structure, application and development process of Xin anjiang Model作者姓名:孔旭学科、专业:水文学及水资源学号:21506149指导教师:王国利完成日期:2016年8月30日大连理工大学Dalian University of Technology摘要新安江模型是河海大学提出的一个概念性降雨径流模型,具有原创性,是我国为数不多的被国际上广泛认可的水文模型。
新安江水文模型在我国湿润与半湿润地区广为应用,取得了良好的效果。
经过近50年的发展,新安江模型已经从最初的专门从事水库入库洪水预报的单一功能模型发展为适合用于水文预报、水资源管理、水土资源评价、面源污染预测、气候变化和人类活动影响研究的多功能的水文模型;其部分参数已从靠经验率定发展为可以进行物理推求。
总之,新安江模型是一个不断发展的模型体系。
本文主要由三部分构成。
第一部分为新安江模型简介,回顾了新安江模型产生的历史背景和发展历程,介绍了新安江模型的基本原理和结构体系;第二部分讲述了新安江模型参数的物理意义及其率定;第三部分为新安江水文模型在英那河流域防洪规划编制当中的应用。
关键词:水文模型;新安江模型;洪水预报The principle, structure, application and development process of Xinanjiang ModelAbstractXin anjiang Model originally proposed by Hehai University is a conceptual rainfall runoff model and is also one of the few widely recognized international hydrological model in China. Xin anjiang hydrological model was widely used in our humid and semi-humid areas, and achieved good results.After nearly 50 years study, Xin anjiang model has been developed from the single-function of reservoir flood forecasting into multi-purpose model including hydrological forecasting, water resources management, water and soil resources evaluation, non-point source pollution prediction, climate change and human activities versatile hydrological model studies. And part of its parameters can be acquired through physical calculation instead of experience. In short, Xin anjiang model is an evolving model system.This paper consists of three parts. The first part is about the brief introduction of Xin anjiang model, which recalls the historical background and the development, as well as introduces the basic principles and architecture; the second part describes the physical meaning of Xin anjiang model parameters and calibration; the third part is about the application of Xin anjiang model in Ying Na River Basin flood control planning.Key Words: hydrological model; Xin anjiang model; Flood forecasting目 录摘 要 .............................................................. I Abstract .. (II)引 言 (1)1 新安江模型简介 (2)1.1 新安江模型起源 ............................... 错误!未定义书签。
一、模型的结构与参数三水源新安江模型的流程图如图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所示的关系。
新安江模型介绍:
三水源新安江模型蒸散发计算采用三层模型;产流计算采用蓄满产流模型;用自由水蓄水库结构将总径流划分为地表径流、壤中流和地下径流三种;流域汇流计算采用线性水库。
模型结构:
模型计算:
在新安江模型中,流域蒸散发计算没有考虑流域内土壤含水量在面上分布的不均匀性,而是按土壤垂向分布的不均匀性将土层分为三层,用三层蒸散发模型计算蒸散发量。
参数有流域平均张力水容量WM(mm),上层张力水容量UM(mm),下层张力水容量LM(mm),深层张力水容量DM(mm),蒸散发折算系数KC和深层蒸散发扩散系数C。
具体计算为
若P+WU>=EP,则EU=EP,EL=0,ED=0;
若P+WU<EP,则EU=P+WU;
若WL>C*LM,则WL=(EP-EU)WL/LM,ED=0;
若WL<C*LM且WL>=C*(EP-EU),则EL=C*(EP-EU),ED=0;
若WL<C*LM且WL<C*(EP-EU),则EL=WL,ED=C*(EP-EU)-WL;
水源划分中,本小组采用的是三水源划分。
三水源用自由水蓄水库结构解决水源划分问题,自由水蓄水库结构考虑了包气带的垂向调蓄作用,按蓄满产流模型计算出总径流量R,先进入自由水蓄水库调蓄,再划分水源。
模型参数调整:
1蒸散发能力折算系数KC
KC是影响产流量计算最为重要和敏感的参数,产流计算中KC控制着水量平衡,因此,对水量计算是最重要的。
KC主要反映流域平均高程与蒸发站高程之间差别的影响和蒸发皿蒸散发于陆面蒸散发间差别的影响。
在实际模拟计算中KC值往往变化很大,最后需经模型调试并验证后确定。
2流域平均张力水容量WM
流域平均张力水容量WM表示流域干旱程度,分为UM,LM,DM。
根据经验,南方湿润
地区WM约为120~150mm,半湿润地区WM约为150~200mm。
3流域蓄水容量—面积分布曲线指数B
B值反映划分单元流域张力水蓄水分布的不均匀程度。
在一般情况下其取值与单元流域面积有关。
在山丘地区,若单元流域面积较小,B=0.1左右;若单元流域面积中等,有几百到一千平方千米,B=0.2~0.3;若单元流域面积有几千平方千米,B=0.4左右。
4不透水面积占全流域面积的比例IM
IM值可由大比例尺的地形图通过GIS现代技术测量出来。
在天然流域,IM=0.01~0.02。
5深层蒸散发扩散系数C
C值主要取决于流域内深根植物的覆盖面积,南方多林地区C=0.15~0.20;在北方半湿润地区C=0.09~0.15。
6自由水蓄水容量SM
SM反映表土蓄水能力,其值受降雨资料时段均化的影响明显。
当以日作为时段长时,在图层很薄的山区,其值为10mm或更小;在土深林茂透水性很强的流域,其值可取50mm 或更大;一般流域在10~20mm。
7自由水蓄水容量—面积分布曲线指数EX
EX反映流域自由水蓄水分布的不均匀程度,一般EX=1.0~1.5
8自由水蓄水库对地下水和壤中流的日出流系数KG+KI
KG的大小反映基岩和深层土壤的渗透性,KI的大小反映表层土的渗透性。
9地下水消退系数CG
CG可根据枯季地下径流退水规律来推求。
若以日作为计算时段长,CG=0.950~0.998
10壤中流消退系数CI
若无壤中流CI=0,若壤中流丰富,则CI=0.9
本小组经过程序的编写,用C++语言编写了新安江模型的计算界面,通过界面填写各种参数的数值调试计算结果,从而达到参数的简单率定。
部分程序:
//三层蒸发模型的时段蒸发计算
for(i=0;i<N;i++)
{
W[i]=WU[i]+WL[i]+WD[i];
if(WU[i]+P[i]>=EP[i])
{
EU[i]=EP[i];
EL[i]=0;
ED[i]=0;
}
else
{
EU[i]=WU[i]+P[i];
if(WL[i]>=C*WLM)
{
EL[i]=(EP[i]-EU[i])*WL[i]/WLM;
ED[i]=0;
}
else if(WL[i]>=C*(EP[i]-EU[i]))
{
EL[i]=C*(EP[i]-EU[i]);
ED[i]=0;
}
else
{
EL[i]=WL[i];
ED[i]=C*(EP[i]-WU[i])-EL[i];
}
}
//时段总蒸发量的计算
E[i]=EU[i]+EL[i]+ED[i];
//时段净雨量的计算
PE[i]=P[i]-E[i];
//时段初蓄水容量值的计算
a[i]=WMM*(1.0-pow((1.0-W[i]/WM),(1.0/(1.0 +B))));
//时段产流量的计算
if(PE[i]>0)
{if(a[i]+PE[i]<=WMM)
R[i]=PE[i]+W[i]-WM+WM*pow((1.0-(PE[i]+a[i] )/WMM),(B+1.0));
else
R[i]=PE[i]-(WM-W[i]);
}
//时段末土壤各层蓄水量的计算
if(i+1==N)
break;
else
{if(PE[i]<=0)
{if(WU[i]+P[i]>=EP[i])
{
WU[i+1]=WU[i]+P[i]-E[i];
WL[i+1]=WL[i];
WD[i+1]=WD[i];
}
else if(WL[i]-EL[i]>0)
{
WU[i+1]=0;
WL[i+1]=WL[i]-EL[i];
WD[i+1]=WD[i];
}
else
if(WD[i]-ED[i]>0)
{
WU[i+1]=0;
WL[i+1]=0; WD[i+1]=WD[i]-ED[i];
}
}
else
{
WU[i+1]=WU[i]+PE[i]-R[i];
WL[i+1]=WL[i];
WD[i+1]=WD[i];
if(WU[i+1]>WUM)
{
WL[i+1]=WL[i]+WU[i+1]-WUM;
WU[i+1]=WUM;
if(WL[i+1]>WLM)
{
WD[i+1]=WD[i]+WL[i+1]-WLM;
WL[i+1]=WLM;
WU[i+1]=WUM;
}
else
WD[i+1]=WD[i];
}
}
}
}
在参数的率定过程中,KC,SM,KG,KI,CI,CG都属于敏感参数,而UM,LM,C,WM,B,IM,EX都是不敏感参数
本小组采用的资料为
模型设计总结:
小组经过将大量的数据资料编辑整理到电脑后,开始进行程序代码的编写。
我们采用的是C++语言进行模型的构造,因为运用C++可以设计出相应的程序界面,方便进行各种参数的调试和计算结果的比较。
因为曾经进行过水文预报课程设计的程序编写,所以在程序设计方面刚开始进展较快,主要的难点在于三水源划分。
因为曾经进行水文预报课程设计的时候,我们程序中采用的是二水源划分,所以在这块内容上,我们需要进行重新摸索和设计,以满足这次任务的要求。
因为对相应的过程和参数的不了解,所以在三水源划分这块程序设计上进展较慢,在听取了学长对其模型的解析和讲解,然后经过于川对第一组的模型设计的讲解,让我们也有了新的认识,对于我们本组的程序设计工作也有很大的帮助。
经过一段时间的摸索,我们组也完成了程序设计部分的工作,然后通过设计的人机对话界面,进行各种参数的调整调试,从而完成了本次实验任务的相关内容。
水文第七组成员:刘俊、冯远、曹胜荣、杨春智
数据输入:刘俊、冯远、曹胜荣、杨春智
程序编写:杨春智
总结报告:杨春智、曹胜荣、冯远、刘俊
小组成员:杨春智、曹胜荣、冯远、刘俊
2010年11月08日。