版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
6新安江模型2023/6/42目录概述二水源新安江模型三水源新安江模型新安江模型的改进新安江模型的应用2023/6/431、概述新安江模型简介一、新安江流域水文模型系列新安江模型是华东水利学院(河海大学)水文系1973年对新安江水库作入库流量预报时提出来的,是一个分块式的概念性流域降雨径流模型。可以用于湿润地区和半湿润地区的湿润季节。最初的新安江模型:二水源模型——地表径流、地下径流;编制新安江入库洪水预报方案80年代初:三水源模型——地面径流、壤中流、地下径流(引入了萨克拉门托模型与水箱模型中的用线性水库函数划分水源的概念);1984~1986年:提出四水源模型——地面径流、壤中流、快速地下径流、慢速地下径流。之后,其它改进。2023/6/44
二、模型的总结构小流域——集总模型大面积流域——分块模型分块模型把流域分成许多块单元流域,对每个单元流域做产汇流计算,可以得到单元流域的出口流量过程。再进行出口一下的河道洪水演算,并得流域的流量过程。把每个单元流域的出流过程相加,就求得了流域出口的总出流过程。划分单元的目的:处理降雨分布的不均匀性(用面雨量带来面积均化);其次是下垫面条件变化(水库等)。因此:单元流域面积要适中,使得在每块面积上降雨比较均匀,并有一定数目的雨量站;(泰森多边形)其次,尽可能是单元流域与自然流域相一致;若流域中有大中型水库,则水库以上的集水面积即可作为一个单元流域。2023/6/45ThiessenPolygonsDACBEA1A2A3A4A5测站ai(%)A24B21C37D8E10合计100RainfallAveragingMethods2023/6/462023/6/47模型结构为了考虑降水和流域下垫面分布不均匀的影响,新安江模型的结构设计为分散性的,分为:蒸散发计算,产流计算,分水源计算和汇流计算四个层次结构。
层次1层2层3层4层功能蒸发计算产流计算水源划分汇流计算二水源三水源坡面汇流河道汇流方法三层模型蓄满产流稳定下渗率自由水水库单位线或线性水库或滞后演算法马斯京根法或滞后演算法参数KCUMLMCWMBIMFCSMEXKGKIUH或CSCICGKEXE或L2023/6/482、二水源新安江模型一、二水源新安江模型基本结构(状态变量(模型参数))降雨P蒸发皿蒸发EI透水面积土壤湿度W上层WU下层WL深层WD径流REUELED蒸散发EWUMWLMC不透水面积IMPWMB地面径流RS地下径流RGFC地面径流过程地下径流过程单元流域出流过程UHKKGKEXE径流R2023/6/49二、二水源新安江模型的微结构(一)用超蓄产流(即“蓄满产流”)模型计算总径流R、地表径流RS及地下径流RG(1)超蓄产流模型概念超蓄产流模型是目前我国湿润地区的主要产流模型。“蓄满”,指含气层的土壤含水量达到田间持水量,而非土壤完全饱和;“超蓄产流”指土壤达到田间持水量以前不产流,所有降雨都被土壤吸收,成为薄膜水和张力水;而在土壤达到田间持水量以后,所有降雨(除去同期蒸发)都产流。这时土壤的下渗能力为稳定下渗率,稳定下渗量FC补充地下水,形成地下径流,而超渗的部分则形成地表径流。与“超渗产流”模型的区别:“超蓄产流”模型先计算R,在分成RS、RG;“超渗产流”模型先计算RS、RG,再合成R。2023/6/410(2)超蓄产流模型的结构
a)点模型以含气层缺水量为控制条件,就流域中某点而言:2023/6/411
b)流域蓄水容量曲线(超蓄产流模型的核心)WWM:流域蓄水容量WWMM:流域最大蓄水容量WM:流域平均蓄水容量2023/6/412利用流域蓄水容量曲线计算产流量(右图):
W:流域原有蓄水量,相应纵标A
W分布:(f/F)A左边蓄满,右边未蓄满,假定按水平分布。
以此时段为基础:降雨P,蒸散发E,径流量R,损失量L满足如下水量平衡关系(超蓄产流方程):大量资料表明,WWM~f/F有如下关系:2023/6/413则:c)流域产流计算P-E>0时,产流,否则不产流,产流时:产流计算特点:雨强对产量无影响,产流量取决于P-E与W。2023/6/414模型参数:WM与BWM:流域干燥时的缺水量,代表流域干旱情况,气候因素;
B:蓄水容量在流域上的分布不均匀性,B=0时分布均匀,愈大愈不均匀,决定于地形、地质条件。
d)地面、地下径流的划分(分水源)产流面积变化,则:
2023/6/415例6-1:超蓄产流模型产流量计算示例WM=120mm,B=0.3,FC=18mm/d年.月.日P-EAf/FRRGRSW1978.7.1716.61
184.2216.890.03380.160.16020.67
199.4021.110.04270.500.50029.57
20-5.9830.510.063200023.682160.3524.270.04957.462.225.2376.572254.2484.620.209117.615.8411.76113.202320.27138.850.484413.4711.961.51120.0024-2.79156.001.0000000117.212023/6/416IMP:不透水面积参数(新安江模型新增参数),流域不透水面积占总面积的百分比,增加后,需修改(6-5),(6-8)式,其它都不变。
尤其半湿润地区需要考虑2023/6/417(二)稳定下渗率fc的推求1、求一场洪水的RS、R、RG(1)据上图求RS(2)根据图求R(3)求RG=R-RS
(4)fc=RG/TT为净雨时间2023/6/418AEGBCHIDt(h)Q(m3/s)F本次降雨形成的径流过程C’D’B’直接径流地下径流N2023/6/419
2、用试算法求fc(三)、不透水面积上的直接径流2023/6/420
(四)、透水面积上的蒸散发模型因为不透水面积上没有蒸散发,因此,计算出来的蒸散发量要乘以透水面积所占比重,才是流域上的蒸散发量。1、蒸散发模型原理蒸散发能力(EM,mm/d)
新安江模型中,认为流域土壤含水量达到最大时,实际蒸散发量E=EM;当土湿很小时,蒸散发量几乎维持为一常数。2、模型结构该模型不考虑蒸散发在面上分布的不均匀性,以模拟土湿纵向分布。可以把土壤分成一、二、三层,现主要采用三层模型。2023/6/421三层模型,其参数有上层张力水蓄水容量
UM,下层张力水蓄水容量LM,深层张力水蓄水容量DM,流域平均张力水蓄水容量WM,蒸散发折算系数KC,深层蒸散发系数C,计算公式为:
WM=UM+LM+DMW=WU+WL+WDE=EU+EL+ED上层(Upperlayer)下层(Lowerlayer)深层(Deeplayer)2023/6/422
2023/6/423
2023/6/4242023/6/425Ifw(1)+p(i)>ep(i)Thene(1)=ep(i)e(2)=0e(3)=0Elsee(1)=w(1)+p(i)e(2)=(ep(i)-e(1))*w(2)/wm(2)Ifw(2)<=c*wm(2)Thene(2)=c*(ep(i)-e(1))e(3)=0Ifw(2)>=c*(ep(i)-e(1))Thene(2)=c*(ep(i)-e(1))e(3)=0Elsee(2)=w(2)e(3)=c*(ep(i)-e(1)-e(2))EndIfEndIfEndIf
w(1)=w(1)+p(i)-r-e(1)w(2)=w(2)-e(2)w(3)=w(3)-e(3)Ifw(1)>wm(1)Thenw(2)=w(1)-wm(1)+w(2)w(1)=wm(1)Ifw(2)>wm(2)Thenw(3)=w(3)+w(2)-wm(2)w(2)=wm(2)EndIfEndIf2023/6/426
C值取决于深根植物面积占流域面积的比重,同时也与(WUM+WLM)值有关,此值越大,C值越小。一般经验,江南湿润地区为0.15~0.20华北半湿润地区0.09~0.12
※蒸散发能力的推求1、多年平均值2、水面蒸发实测3、气象因素推算
2023/6/4273、模型的计算(1)计算蒸散发能力EM=K×EI(2)计算PEPE=P-E(3)利用PE按超蓄产流计算R,PE<0时不产流(4)计算WUi+1=WUi+P(5)计算EU,EL(6)计算ED(7)计算E=EL+EU+ED(8)计算WU(9)计算WL(10)计算WD4、确定K值:蒸发皿系数2023/6/428
5、单元流域汇流计算(1)地面径流采用经验单位线法无因次单位线相同
(2)地下径流汇流计算
QRG1,QRG2-时段初末的地下径流量;KKG-地下径流日退水系数-时段长度;F-单元流域面积;D-一日内时段数;RG-时段内地下径流产流量
2023/6/4296、河槽汇流计算特征河长法
马斯京根法(扩散波解的差分求解)滞后演算法线性扩散模拟法2023/6/4302023/6/431二水源新安江模型参数确定新二模型参数:K—流域蒸散发能力于蒸发皿蒸发量之比
C—深层蒸散发系数
IMP—不透水面积比重
WM—流域平均蓄水容量(指张力水)
WUM—流域平均上层蓄水容量
WLM—流域平均下层蓄水容量
WDM—流域平均深层蓄水容量
B—蓄水容量曲线指数
FC—稳定下渗量
KKC—地下水日退水系数
UH—无因次的地表径流单位线纵表
KE—单元河段的马斯京根K值
XE—单元河段的马斯京根X值2023/6/432新二模型参数初值确定:B—蓄水容量曲线指数,反映流域的不均匀性,流域越大,则B越大小300平方公里大0.10.2—0.30.3—0.4FC—稳定下渗量,各场雨不同,相差很大,需注意。
KKG—地下水日退水系数
UH—无因次的地表径流单位线纵表,单元流域的地面径流的单位线,无因次,可找相邻流域值作为初值
KE—单元河段的马斯京根K值
XE—单元河段的马斯京根X值
KE、XE可以根据河段特性,用水力学方法求出2023/6/433新二模型参数初值确定:K—流域蒸散发能力于蒸发皿蒸发量之比,E601可以作为初值,但要高程修正
C—深层蒸散发系数,决定于深根植物占流域面积的比数,同时和WUM+WLM有关,南方0.15—0.2,北方0.09—0.12.
IMP—不透水面积比重,干旱降小雨,有一个小洪水,此时径流系数就是IMP,也可以在地图上量出
WM—流域平均蓄水容量(指张力水),反映流域干旱程度,久旱下大雨的资料可以分析,雨前为0,雨后为WM
WUM—流域平均上层蓄水容量,20mm,差5-10mm
WLM—流域平均下层蓄水容量,60-90mm
WDM—流域平均深层蓄水容量,
2023/6/434新二模型参数确定的步骤:1、定初始值:取5年资料,以天为时段,进行计算2、比较多年径流:最基本的水量平衡,先改K,冬夏不同3、比较每年径流:干旱和湿润年份有无系统误差调WUM和WLM,以及C,WUM变小,雨季蒸发小,旱季影响不大,C加大,干旱季节蒸发加大4、年内干湿差比较:如洪水计算值偏大,调WUM,WLM和C,如W在久旱后出现负值,加大WM不改WUM和WLM5、比较枯季地下径流:如有系统偏差,调FC,快慢调KKG6、比较小洪水:可以调IMP和B,湿润区不敏感7、比较地面径流过程:降雨中心误差造成汇流偏早偏迟,调KE,仍有误差,调UH和XE2023/6/435存在的主要问题:①用FC划分水源是建立在包气带岩土结构为水平方向空间分布均匀的基础上,这假定往往与实际情况不符。②用FC划分水源没有考虑包气带的调蓄作用,在某些流域实际计算结果表明,壤中流的坡面调蓄作用有时比地面径流大得多;直接进入地下水库没有考虑坡面垂向调节作用,即包气带的调蓄作用;由于地表径流和壤中流的汇流规律和汇流速度不同,两者合在一起采用同一种方法进行计算,常会引起汇流的非线性变化。③对许多流域资料的分析表明,即使是同一流域,各次洪水所分析出的也不相同,而且有的时候变化很大,很难进行地区综合和在时空上外延,应用时任意性大,常造成较大误差。
Hortonoverlandflowdominateshydrograph;contributionsfromsubsurfacestormflowarelessimportantDirectprecipitationandreturnflowdominatehydrograph;subsurfacestormflowlessimportantSubsurfacestormflowdominateshydrographvolumetrically;peaksproducedbyreturnflowanddirectprecipitationAridtosub-humidclimate;thinvegetationordisturbedbyhumansHumidclimate;densevegetationSteep,straighthillslopes;deep,verypermeablesoils;narrowvalleybottomsThinsoils;gentleconcavefootslopes;widevalleybottoms;soilsofhightolowpermeabilityClimate,vegetationandlanduseTopographyandsoilsVariablesourceconceptHowdorunoffprocessesdependonClimate,Vegetation,Landuse,TopographyandSoils?(FromDunneandLeopold,1978)2023/6/437二水源依据霍尔顿概念.没有考虑壤中流造成的没有包气带的调蓄作用2023/6/438超蓄产流自由水蓄水库2023/6/439有积水的面积形成了界面上的产流面积。在此面积上,下渗的部分是地下径流,超渗的部分是直接径流。总产流量R增大,积水面积也增大。但这个过程与供水强度有关,因为积水要按FCB下渗。强度小,历时长,下渗就多,积水面积增加得慢。强度大,历时短,下渗量少,积水面积增加得快(但FCB比FCA相对很小,所以这个现象不容易显示出来)。当积水面积增加后,地下径流与直接径流都应增加,如图所示。但当R大到一定程度后,全流域面积都产流,则地下径流就不能再增加了,只能等于全流城面积的FCB,接近一常数。影响产流的不均匀因素:1)B层土壤稳定入渗FCB分布不均。2)在A、B界面上不均匀分布的洼蓄量。3)A层土壤供水强度的不均匀分布。2023/6/4403、新安江模型的改进一、三水源新安江模型1、概述
二水源模型由于没有考虑壤中流作用,故在壤中流丰富流域常常得不到好的模拟结果。三水源模型认为,土壤中水有张力水(田间持水量以下的水)和自由水(田间持水量以上的水)之分。二水源模型只考虑了张力水的调蓄作用,没有考虑自由水的调蓄作用。因此,新三模型中增加了一个自由水蓄水库,把总径流划分成三种水源:地面径流、壤中流、地下径流,代替新二模型中用FC划分水源的办法。2023/6/441自由水蓄水库①②③小单元2023/6/442
2、新三模型的基本结构
与新二模型基本相同,不产流时是一致的;产流后,在产流面积上R先进入自由蓄水库,形成RG、RSS,且只出现在产流面积上,去掉FC,增加参数:KG:地下水从自由蓄水库中的出流系数
KSS:壤中流出流系数
SM:自由蓄水库容量(即最大蓄量,mm)则划分水源的计算为:
底宽地下径流壤中流2023/6/443输入降雨P,蒸发皿蒸发量EI透水面积不产流面积1-FR产流面积FR产流量R张力水W上层WU下层WL深层WD自由水S不透水面积产流量RIMPWMBIMP输出蒸散发EEUEUEUWUMWLMC地面径流RS壤中流RSS地下径流RGSMEXKSSKG地下径流流域出口流量KKG地面径流及壤中流总入流地面径流及壤中流总入流UH输出总流量Q三水源新安江模型流程图2023/6/4443、自由水蓄水容量曲线F-全流域面积;FS-自由水蓄水量小于等于SS的面积;EX-抛物线指数;SSM-流域上自由水蓄水量最大的某点的蓄量值;S-流域平均自由水蓄水量;SM-流域平均自由水蓄水容量;AU-与自由水蓄水量S对应的蓄水容量曲线的纵标值RGRI2023/6/4454、利用自由水蓄水容量曲线计算产流量
FR=R/PE注:以上计算认为SS-FS/F曲线是按全流域(FR=1)计算产流,再缩至FR。RGRI2023/6/4465、汇流计算(1)二水源汇流计算①地面径流汇流地面径流汇流采用单位线法,计算公式为:
②地下径流汇流地下径流汇流可采用线性水库或滞后演算法模拟。当采用滞后演算法时,计算公式为:
CG为消退系数;U——单位转换系数,U=流域面积F(km2)/(3.6△t)2023/6/447
③单元面积河网总入流单元面积河网总入流为地面径流与地下径流出流之和,计算公式为:
④单元面积河网汇流单元面积河网汇流可采用线性水库或滞后演算法模拟。当采用滞后演算法时,计算公式为:
单元面积河网汇流计算在很多情况以简化。这是由于单元流域的面积一般不大而且其河道较短,对水流运动的调蓄作用通常较小。在单元流域面积较大或流域坡面汇流极其复杂的情况下,才考虑单元面积内的河网汇流。
河网消退系数2023/6/448
从单元面积以下到流域出口是河道汇流阶段。河道汇流计算采用马斯京根分段连续演算法。参数有槽蓄系数KE和流量比重因素XE,各单元河段的参数取相同值。为了保证马斯京根法的两个线性条件,每个单元河段取KE≈Δt,已知KE、XE和Δt,可求出C0、C1和C2,即可用下式进行河道演算。
(2)三水源汇流计算①地表径流汇流地表径流的坡地汇流可以采用单位线,也可以采用线性水库,采用单位线的计算公式与二水源相同,采用线性水库的计算公式为:
2023/6/449
②壤中流汇流壤中流汇流可采用线性水库或滞后演算法模拟。当采用线性水库时,计算公式为:
③地下径流汇流采用线性水库时,与二水源相同相同。④单元面积河网总入流
2023/6/450参数的物理意义
(参数可根据其物理意义,分为4类)
(1)蒸散发参数:K、WUM、WLM、CK为蒸散发能力折算系数,是指流域蒸散发能力与实测水面蒸发值之比。此参数控制着总水量平衡,因此,对水量计算是重要的。
WUM为上层蓄水容量,它包括植物截留量。在植被与土壤很好的流域,约为20mm;在植被与土壤颇差的流域,约为5~6mm。WLM为下层蓄水容量。可取60~90mm。C为深层蒸散发系数。它决定于深根植物占流域面积的比数,同时也与WUM+WLM值有关,此值越大,深层蒸散发越困难。一般经验,在江南湿润地区值C约为0.15~0.20左右,而在华北半湿润地区则在0.09~0.12左右。新安江模型的参数2023/6/451(2)产流量参数:WM、B、IMPWM为流域蓄水容量,是流域干旱程度的指标。找久旱以后下大雨的资料,如雨前可认为蓄水量为0,雨后可认为已蓄满,则此次洪水的总损失量就是WM。一般分为上层WUM、下层WLM和深层WDM。在南方约为120mm,北方半湿润地区约为180mm。
B为蓄水容量曲线的方次。它反映流域上蓄水容量分布的不均匀性。一般经验,流域越大,值也越大。在山丘区,很小面积(几平方公里)的为0.1左右,中等面积(300平方公里以内)的为0.2~0.3左右,较大面积(数千平方公里)的值为0.3~0.4左右。IMP为不透水面积占全流域面积之比。可找干旱期降小雨的资料来分析,这时有一很小洪水,完全是不透水面积上产生的。求出此洪水的径流系数,就是IMP。新安江模型的参数2023/6/452
-2002.12-
(3)水源划分参数:SM、EX、KSS、KG
SM为流域平均自由水蓄水容量,本参数受降雨资料时段均化的影响,当用日为时段长时,一般流域的SM值约为10~50mm。当所取时段长较少时,SM要加大,这个参数对地面径流的多少起着决定性作用,因此很重要。EX为自由水蓄水容量曲线指数,它表示自由水容量分布不均匀性。通常EX取值在1~1.5之间。KSS为自由水蓄水库对壤中流的出流系数,KG为自由水蓄水库对地下径流出流系数,这两个出流系数是并联的,其和代表着自由水出流的快慢。一般来说,KSS+KG=0.7,相当于从雨止到壤中流止的时间为3天。
新安江模型的参数2023/6/453
-2002.12-
(4)汇流参数:KKSS、KKG、CS、LKKSS为壤中流水库的消退系数。如无深层壤中流时,KKSS趋于零。当深层壤中流很丰富时,KKSS趋于0.9。相当于汇流时间为10天。KKG为地下水库的消退系数。如以日为时段长,此值一般为0.98~0.998,相当于汇流时间为50~500日。CS为河网蓄水消退系数,L为滞时,它们决定于河网地貌。新安江模型的参数2023/6/454
-2002.12-模型参数率定
新安江模型的参数按照物理意义分为4层,上面已作了介绍。参数的率定可以按照蒸散发~产流~分水源~汇流的次序进行,各类参数基本上是相互独立的。按照以下次序率定参数。新安江模型的参数2023/6/455日模型(1)定出各参数的初始值。(2)比较多年总径流。这是最基本的水量平衡校核。如有误差,要首先修改K值,K是影响蒸发计算最大的参数,对于某些北方河流,夏季植物茂盛,而冬季则有封冻。冬季蒸发不可能用E601观测,则应考虑把分为冬、夏各月定为不同的数值。(3)多年总水量基本平衡后,再比较每年的径流,看很干旱的年与湿润年份有无系统误差。如有应调整WUM、WLM和C。WUM减小将使少雨季节的蒸发减少,而对于很干旱的季节则无影响。WLM的作用与此相仿。加大C值将使很干旱的季节的蒸散发增大,而对于有雨季节则无此影响。在北方半湿润地区可以找到干旱年份与湿润年份之间的系统误差,而在南方湿润地区则不易找到。新安江模型的参数参数的率定2023/6/456(4)如上述差别并不明显,则应比较年内干旱季与湿润季之间的差别。在南方,主要是伏旱季的蒸散发计算是否正确的问题。如伏旱以后的初次洪水具有系统误差,例如,各年中这种洪水的计算值都偏大,则应调整WUM、WLM和值C,使基本符合。如果在计算中发现值在久旱后W出现负值,则应加大WM,不改变WUM和WLM。在计算中当为负值时以零处理是不对的,它破坏了产流量计算的前提。新安江模型是蓄满型,只要蒸散发计算基本正
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2024年度绿化工程承包合同
- 大班种子课件教学课件
- 2024山西劳动合同范本
- 2024年度J企业卫星通信技术服务合同
- 2024年店面续租协议:市中心
- 2024互联网销售涂料产品独家代理合同
- 2024年工程进度与安全合同
- 2024年建筑修正协议
- 2024年家用电器维修服务合同
- 2024双方关于影视制作与发行委托合同
- 业主业主委员会通用课件
- 了解金融市场和金融产品
- 南京理工大学2015年613物理化学(含答案)考研真题
- 初中数学应用题解题思路分享
- 安全生产科技创新与应用
- 人工智能在文化传承与遗产保护中的价值实现
- 2024年汽修厂开业计划书
- ISTA标准-2A、2B、2C系列解读(图文)
- 日间手术应急预案方案
- 退费账户确认书
- 幼儿园小班《汽车滴滴响》
评论
0/150
提交评论