栅格型新安江模型的参数估计及应用.
栅格型新安江模型的参数估计及应用
栅格型新安江模型的参数估计及应用姚成;纪益秋;李致家;刘开磊【摘要】Based on the theory of the Xin'anjiang model and one-dimensional diffusion wave model, the grid-based Xin anjiang model (the Grid-Xin' anjiang model) was developed. This model regarded the grid DEM as the computational element, and the generated runoff of each element was divided into three components, the surface runoff, interflow, and groundwater runoff, which were routed using the cell-to-cell diffusion wave flow routing method according to the computational order of each grid cell in the DEM. During the calculation of the runoff generation and runoff concentration, the water exchange between grids and the influence of the river drainage network were taken into consideration. Based on the underlying surface conditions, including the topography, and the types of soil and vegetation in the study basin, a method for parameter estimation of the model was proposed and verified. The model was applied to the Tunxi Basin in Anhui Province for flood simulation and produced favorable results.%根据新安江模型与一维扩散波模型理论,构建了基于栅格的新安江( Grid-Xin' anjiang)模型.该模型以DEM栅格为计算单元,将每个单元的产流量划分为地表径流、壤中流以及地下径流3种水源,最后再根据栅格间的汇流演算次序,利用扩散波汇流方法依次将各种水源演算至流域出口.模型在进行产汇流计算时,考虑了栅格之间的水量交换以及河道排水网络的影响.以流域地貌特征、土壤及植被类型等下垫面条件为基础,研究了模型参数的估计方法,并对其进行了验证.将模型用于安徽省屯溪流域的洪水模拟,取得了良好的应用效果.【期刊名称】《河海大学学报（自然科学版）》【年(卷),期】2012(040)001【总页数】6页(P42-47)【关键词】栅格型新安江模型;新安江模型;扩散波模型;参数估计;洪水模拟;屯溪流域【作者】姚成;纪益秋;李致家;刘开磊【作者单位】河海大学水文水资源学院,江苏南京210098;南通市港闸区幸福街道办事处,江苏南通226012;河海大学水文水资源学院,江苏南京210098;河海大学水文水资源学院,江苏南京210098【正文语种】中文【中图分类】P333.1基于栅格的分布式水文模型是目前水文研究的热点之一.1969年,Freeze等[1]发表的“一个具有物理基础数值模拟的水文响应模型的蓝图(FH69蓝图)”被认为是分布式水文模型研究的开始.过去受制于技术的限制,许多流域信息、特征不能及时准确地获得,地理信息无法及时准确地获得和使用,使得分布式水文模型发展缓慢.而数字高程模型(DEM)的出现为分布式水文模型的构建创造了条件,分布式水文模型取得了飞速发展.比较具有代表性的分布式和半分布式流域水文模型有SHE[2],TOPKAPI[3]和TOPMODEL[4]模型等.随着信息、遥感、计算机等技术的发展,基于DEM的分布式水文模型以其充分考虑水文要素和各种参数空间变化的特点正成为水文模型的发展趋势,而新安江模型[5]作为一个半分布式的概念性模型,自研制以来,就得到了广大学者的认同,在中国的洪水预报中也得到了广泛成功的应用,因此研究栅格型新安江模型具有一定的现实意义[6-7].本研究利用DEM,提取了水文模拟与参数估计所需要的流域地貌特征,并根据新安江模型与扩散波模型理论,构建了栅格型新安江(Grid-Xin'anjiang)模型.同时,本文结合土壤类型、植被覆盖及利用等信息,对空间变化模型参数的估计方法进行了一定的研究与验证,并将Grid-Xin'anjiang模型应用于安徽省屯溪流域洪水模拟计算.1 Grid-Xin'anjiang模型原理Grid-Xin'anjiang模型以流域内每个DEM栅格作为计算单元,并假设在栅格单元内的降雨和地貌特征、土壤类型以及植被覆盖等下垫面条件空间分布均匀,模型只考虑各个要素在不同栅格之间的变异性.Grid-Xin'anjiang模型先计算出每个栅格单元的植被冠层截留量、河道降水量和蒸散发量,然后再计算出栅格单元的产流量,并采用自由水蓄水库结构对其进行水源划分,即划分为地表径流、壤中流以及地下径流3种水源,最后再根据栅格间的汇流演算次序,依次将各种水源演算至流域出口.在进行次洪模型计算时,采用一维扩散波模型,且假设在原来的坡地栅格内也存在一个“虚拟河道”,每个栅格单元的壤中流和地下径流直接流入河道,采用河道汇流的计算方法,则栅格间的扩散波汇流演算只包括坡面汇流和河道汇流2种方式.考虑到日径流模拟时对汇流演算精度要求不高且时间步长较大,为了增加模型的运行效率,Grid-Xin'anjiang模型采用简便的线性水库法和滞后演算法进行栅格间汇流计算.1.1 植被冠层截留植被冠层截留是指降雨在植被冠层表面的吸着力、承托力及水分重力、表面张力等作用下储存于其表面的现象[8].在一次降雨过程中,植被冠层对降雨的累积截留量[9]可表示为式中:I cum——植被冠层的累积截留量;f lc——植被覆盖率;P cum——累积降雨量;C vd——植被密度的校正因子;S cmax——植被冠层的截留能力,即植被冠层的最大截留量.1.2 河道降水假设在每一个河道栅格内,降雨分布均匀,河道形状不发生改变,则可以根据栅格单元内河道部分所占的面积比例进行河道降水的计算.河道降水I ch可表示为式中:L ch——河道长度;W ch——河道断面最大过水面积所对应的水面宽;A gr——栅格单元的面积;P——时段降雨量.其中,L ch W ch表示河道所占面积.对于流域内的坡地栅格,即使有“虚拟河道”存在,也不考虑河道降水的影响,即在坡地栅格上I ch=0.1.3 蒸散发将每个栅格单元内的土壤分为上层、下层和深层3层,每一层对应的张力水蓄水容量分别为 W UM,W LM和W DM.在栅格单元实际蒸散发计算时,冠层截留量按蒸散发能力蒸发,当截留水量小于蒸散发能力时,则采用3层蒸散发模型[10].1.4 单元产流及分水源土壤蓄满表示的是土壤含水量达到田间持水量,而不是饱和含水量.新安江模型引进张力水蓄水容量分布曲线来考虑土壤含水量面上分布不均的问题,而在Grid-Xin'anjiang模型中,由流域地貌特征以及土壤、植被等下垫面条件来确定任意一个栅格单元的张力水蓄水容量W M,而在栅格单元内暂不考虑张力水含水量分布不均的问题.将计算时段内栅格单元的实测降雨先扣除相应时段的蒸散发、植被冠层截留、河道降水后,再考虑上游入流是否补足当前单元的土壤含水量,即可得到实际用于产流计算的时段雨量P e.当P e≤0或P e+W0≤W M时,R=0;当 P e+W0＞W M时,R=P e+W0-W M.其中:R为时段产流量;W0为栅格单元实际的张力水含量.在Grid-Xin'anjiang模型中,任意栅格单元内的产流量R均被划分为3种水源:地面径流R s、壤中流R i和地下径流R g.所需参数包括表层土自由水容量S M、表层自由水含量对壤中流的出流系数K i、表层自由水含量对地下水的出流系数K g.1.5 汇流演算Grid-Xin'anjiang模型在进行短时段栅格间扩散波汇流演算时,假设任意栅格单元都由坡地和河道组成,即原来的坡地栅格上也存在一个“虚拟河道”,地下径流与壤中流都直接汇入河道或“虚拟河道”中,因此栅格间的汇流就由坡面汇流及河道汇流组成,均采用扩散波模型.在模型进行产汇流计算时,考虑了栅格间的水量交换问题以及河道排水网络的影响.即:如果当前计算单元的土壤含水量处于未蓄满状态,则上游栅格的入流量首先补充当前栅格的土壤含水量,直至其蓄满为止;如果当前栅格有河道存在,则地表径流的出流量将先按一定的比例汇入河道,然后再汇至下游栅格. Grid-Xin'anjiang模型的坡面水流以及河道水流运动均利用一维扩散波方程组来描述,在进行栅格间汇流演算时,坡面水流运动方程需在每个栅格单元上进行离散,其中的连续性方程为式中:h s——坡面水流的水深;t——时间;A gc——栅格单元的面积;Q s——栅格单元的地表径流流量;Q sout——栅格单元的地表径流出流量;Q sup——上游栅格入流量.本文采用基于两步MacCormack算法[11-12]的二阶显式有限差分格式进行坡面与河道水流扩散波方程组的求解.2 研究流域概况选用湿润地区的屯溪流域作为研究流域.该流域位于安徽省境内皖南山区,流域面积为2692.7 km2.屯溪流域地势西高东低,最大、最小以及平均海拔高程分别为1398m,116m,380m,相对高差较大.该流域雨量充沛,多年平均降雨量约为1800mm,降水在年内年际分配极不均匀,汛期内的降雨量一般占年总雨量的60%以上.屯溪流域内植被良好,主要包括常绿针叶林、落叶阔叶林、混合林、森林地、林地草原、牧草地与作物地,土壤类型主要为黏壤土.本研究采用USGS提供的精度为1km(30″)的DEM高程数据提取了屯溪流域的地形地貌特征.屯溪流域水系及站点分布见图1.图1 屯溪流域水系及站点分布Fig.1 River network and distribution of gauging stations in Tunxi Basin3 Grid-Xin'anjiang模型参数估计3.1 蒸散发参数在湿润流域,平均的上、下层张力水容量一般可以取20mm与60mm,而整个包气带张力水容量常以120mm作为估值[5,10].以此为基础,对与屯溪流域W UM与W LM的估计分别取W UM≈0.167W M;W LM≈0.5W M.参数C与栅格单元的植被覆盖率有关,在植被密集地区可取0.18,因此可令C=0.18f lc,f lc可通过相关公式求得.本研究暂未考虑蒸散发折算系数K在不同栅格单元的高程修正问题,而是认为它在流域内空间分布均匀,因此该参数主要与测量水面蒸发所用的蒸发器有关.对于国内普遍采用的E-601蒸发器而言,可取K≈1.3.2 产流及分水源参数单元产流及分水源计算中W M与S M[13]可表示为式中:W M——栅格单元张力水蓄水容量;S M——栅格单元自由水蓄水容量;θs——饱和含水量;θfc——田间持水量;θwp——凋萎含水量;L a——包气带厚度;L h——腐殖质土层厚度.式(4)与式(5)中的θs,θfc,θwp均可以根据栅格单元的土壤类型通过查土壤参数统计表获取,因此只要知道每个栅格单元的L a与L h即可获得W M与S M在流域的空间分布.在自然界中,影响包气带厚度的因素较多,很难进行直接推求.L a与 L h 可通过与地形指数及土壤类型对应的土壤水分常数进行估算,可假定地形指数大的地方包气带较薄而地形指数小的地方包气带较厚.这与实际情况基本相符.一般而言,在湿润地区地形指数大的地方大多位于河道附近,而这些区域的地下水埋深较浅,包气带相对较薄;相反,地形指数小的地方基本位于流域的上游山坡,远离河道,包气带相对较厚.因此,可以假设流域上地形指数最大的栅格单元对应的张力水蓄水容量最小,而地形指数最小的栅格单元对应的张力水蓄水容量最大.图2为估计的屯溪流域包气带厚度以及腐殖质土厚度的空间分布.图2 屯溪流域土壤厚度分布Fig.2 Distribution of thickness of soil in Tunxi BasinK i与K g这2个参数属于并联参数,其和K i+K g代表的是自由水出流的快慢,应与单元的土壤类型有关,而自由水指的是饱和含水量与田间持水量之间那部分可以在重力作用下自由流动的水,因此可以将θs与θfc作为衡量自由水出流快慢的指标.K i/K g表示的是壤中流与地下径流的比,此比值可以通过θwp来反映[14].具体计算公式为式中:m oc——自由水出流综合影响因子;m r——自由水出流校正系数,根据已有研究成果,在估计模型参数时,可以取m r=1.m oc可直接通过新安江模型分水源中的结构性约束(K i+K g=0.7)[5]选取.具体步骤是:先给m oc赋予一初值,然后根据每个栅格单元的土壤类型由土壤参数统计表和式(6)确定栅格单元对应的K i+K g,再统计出流域内所有栅格K i+K g的均值并与0.7进行比较,由此对m oc进行调整,使其尽量满足该结构性约束.3.3 汇流参数当Grid-Xin'anjiang模型进行扩散波汇流演算时,模型参数的估计主要是基于流域地貌特征以及河道断面信息.其中,坡面汇流的曼宁糙率系数n h根据栅格单元的植被类型由植被参数统计表确定.而每个栅格单元河道(包括虚拟河道)汇流的曼宁糙率系数n c[15]由下式进行计算:式中:A d——栅格单元的上游汇水面积,km2;S oc——河道坡度,可由DEM提取得到;n0——系数,n0可以由流域出口点的糙率及坡度代入式(7)反算出.地表径流汇入河道的比例也是在生成的研究流域水系基础上,采用面积比例法进行计算[7].在判断河道形状时,主要是根据流域内实测站点的断面数据,反演出每一栅格单元的河道形状.以屯溪流域为例,定义某一断面宽度指数α,考虑到随着上游汇水面积的增加,越到流域下游其河道过水断面面积应该越大,且变化比较明显,因此可认为α在流域内应当是变化的,即河道断面的尺寸是空间变化的.根据流域内已有的实测断面资料分析,即可获得α的空间分布.4 应用结果分析屯溪流域具有1982—2001年以及2006年的日降雨资料,有40场洪水的时段降雨资料,屯溪出口站有对应的流量资料.本文将1982—1995年的日资料和对应的次洪资料用于模型参数的率定,剩余资料用于模型参数的检验,对于模型参数估计方法的验证采用率定期资料进行,对于屯溪流域降雨资料的空间插值采用距离平方反比法.采用率定期内的资料对模型估计参数进行验证.其中,日模型模拟的径流深以及洪峰相对误差的合格率分别为100%和79%,平均相对误差水平分别为3.7%和12.6%,确定性系数的均值为0.91.次洪模型模拟的径流深、洪峰相对误差以及峰现时差的合格率分别为91%,48%和83%,平均相对误差水平分别为9.1%,19.9%和2.2h,确定性系数的均值为0.87.图3为应用估计参数模拟的2次洪水过程线比较.图3 Grid-Xin'anjiang次洪模型估计参数模拟结果比较Fig.3 Comparison ofestimated parameters of Grid-Xin'an jiang model for floods由估计参数的模拟结果可以看出,根据土壤、植被及地貌等下垫面条件估计出的Grid-Xin'anjiang日模型参数值较为合理,对于屯溪流域率定期内的日径流模拟精度较高,因此直接将日模型的参数估值用于该流域检验期内的日径流模拟,检验期内模拟的径流深与洪峰相对误差合格率均为91%,平均相对误差水平分别为10.1%和7.0%,确定性系数的均值为0.90.对于Grid-Xin'anjiang次洪模型而言,需要率定检验的主要是扩散波汇流方法中的坡面糙率系数n h,率定方法选择人工试错法.对于次洪模型中其他参数则不再进行率定.通过率定扩散波汇流方法的糙率系数,对于率定期内洪水而言,Grid-Xin'anjiang次洪模型模拟的径流深、洪峰和峰现时间的平均相对误差水平分别为7.7%,7.9%和1.8 h;合格率分别为96%,96%和91%,确定性系数的均值为0.91.对于检验期内洪水而言,Grid-Xin'anjiang次洪模型模拟的径流深、洪峰和峰现时间的平均相对误差水平分别为8.0%,9.7%和2.6h;合格率均为94%,确定性系数的均值为0.90,模拟结果较好.图4为应用率定参数模拟的2次洪水过程线比较.图4 Grid-Xin'anjiang次洪模型率定参数模拟结果比较Fig.4 Com parison of calibrated parameters of Grid-Xin'an jiang model for floods5 结语以新安江蓄满产流模型与一维扩散波汇流模型为基础,构建了适用于湿润流域的分布式Grid-Xin'anjiang模型,并在模型进行产汇流计算时,考虑了栅格间的水量交换问题以及河道排水网络的影响.根据流域的地貌特征、土壤类型及植被覆盖等下垫面特性,对基于物理机制的模型参数估计方法进行了一定的研究与验证.从Grid-Xin'anjiang模型在屯溪流域的应用情况可以看出,模型的应用效果较好,模拟精度较高.为了使Grid-Xin'anjiang模型可以更好地应用于不同类型的流域,还应当加入超渗产流、融雪模型等其他水文过程的计算方法,对于栅格间壤中流与地下径流的汇流演算方法也有待进一步改进.为了使模型更好地用于缺乏资料地区的洪水模拟及预报,还需要进一步加强对模型的应用研究,开展模型在其他流域以及与其他模型应用的对比分析,同时也需要进一步完善模型参数的估计方法.参考文献:【相关文献】[1]FREEZE R,HARLAN R.Blueprint for a physically-based,digitally-simulated hydrologic response model[J].Journal of Hydrology,1969,9:237-258.[2]ABBOTT M,BATHURST J,CUNGE J,et al.An introduction to the European Hydrological System-Système Hydrologique Européen,“SHE”,1:history and philosophy of a physically-based,distributed modelling system[J].Journal of Hydrology,1986,87:45-59. [3]TODINI E,CIARAPICA L.The TOPKAPI model[M]//SINGH V.Mathematical models of large watershed hydrology.Colorado:Water Resources Publications,2001:471-504.[4]BEVEN K,KIR KBY M.A physically based variable contributing area model of basin hydrology[J].Hydrology Science Bulletin,1979,24(1):43-69.[5]赵人俊.流域水文模拟[M].北京:水利电力出版社,1984.[6]李致家,姚成,汪中华.基于栅格的新安江模型的构建和应用[J].河海大学学报:自然科学版,2007,35(2):131-134.(LI Zhijia,YAO Cheng,WANG Zhong-hua.Development and application of grid-based Xin'anjiang model[J].Journal of Hohai University:Nature Sciences,2007,35(2):131-134.(in Chinese))[7]YAO Cheng,LI Zhi-jia,BAO Hong-jun,et al.Application of a developed Grid-Xin'anjiang model to Chinese watersheds for flood forecasting purpose[J].Journal of Hydrologic Engineering,2009,14(9):923-934.[8]芮孝芳.水文学原理[M].北京:中国水利水电出版社,2004.[9]ASTON A.Rainfall intercep tion by eight small trees[J].Journal ofHydrology,1979,42:383-396.[10]赵人俊,王佩兰.新安江模型参数的分析[J].水文,1988(6):2-9.(ZHAO Ren-jun,WANG Pei-lan.Analysis of the parameters of the Xin'anjiang model[J].Journal of China Hydrology,1988(6):2-9.(in Chinese))[11]MACCORMACK R.Numerical solution of interaction of shock wave with laminar boundary layer[J].Lecture Notes in Physics,1971,8:151-163.[12]WANG M,HJELMFELT A.DEM based overland flow routing model[J].Journal of Hydrologic Engineering,1998,3(1):1-8.[13]李致家.水文模型的应用与研究[M].南京:河海大学出版社,2008.[14]KOREN V,SMITH M,WANG D,et e of soil property data in the derivation of conceptual rainfall-runoff model parameters[C]//Proceedings of the 15th Conference on Hydrology.Long Beach:AMS,2000:103-106.[15]KOREN V,REED S,SMITH M,et al.Hydrology laboratory research modeling system(HL-RMS)of the national weather service[J].Journal of Hydrology,2004,291:297-318.。
栅格型新安江模型的参数估计及应用
关键 词 ： 格 型新 安 江模 型 ； 安 江模 型 ； 散 波模 型 ； 栅 新 扩 参数 估计 ； 洪水模 拟 ； 屯溪流域
中图分 类 号 ：３ ３ １ Ｐ３． 文 献标 志码 ： Ａ 文章 编号 ：００ １８ （０ ２ ０ — ０ ２ ０ １０ —９０ ２ １ ） １０ ４ — ６
Ｄ ：０ ３ ７ ／ ． ｓ ．０ ０１８ ．０ ２ ０ ．０ ＯＩ １ ．８ ６ ｊｉ ｎ １０ －９ ０ ２ １ ． １０ ９ ｓ
栅 格 型 新 安 江 模 型 的参 数 估 计 及 应 用
姚 成 纪 益秋２李 致 家 刘 开 磊 ， ， ，
（． １河海大学水文水 资源学院 ， 苏 南京 江 ２０９ ； ． １０８ ２ 南通市港 闸区幸福街 道办事处 ， 江苏 南通 ２６ １ ） ２０２
ｃ ｍｐ ｔｔ ｎ ｌ ｒ ｅ ｏ ｅ ｃ ｇ ｉ ｃ ｌ ｉ ｔ ｅ ｏ ｕ ａｉ ａ ｏ ｒ ｆ ａ ｈ ｒｄ ｅｌ ｎ ｈ ＤＥＭ ． Ｄｕ ｎ ｔ ｅ ａｃ ｌｔｏ ｏ ｔｅ ｕ ｎ ｆ ｇ ｎ ｒｔｏ ａ ｄ ｕ ｏ ｏ ｄ ｉ ｒ ｇ ｈ ｃ ｕ ａｉｎ ｆ ｈ ｒ ｏ ｌ ｅ ｅ ａｉｎ ｎ ｒ ｎ ｆ ｃ ｎ ｅ ｔａｉｎ ｔ ｅ ｏ ｃ ｎｒ ｔ ， ｈ ｗａｅ ｅ ｃ ａ ｇ ｂ ｔ ｅ ｇ ｄ ａ ｄ ｈ ｉｆｕ ｎ ｅ ｆ ｔ ｅ ｉｖ ｒ ｒ ｉａ ｅ ｅ ｏ ｋ ｏ ｔｒ ｘ ｈ ｅ ｅｗｅ ｎ ｒ ｓ ｎ ｔｅ ｎ ｅ ｃ ｏ ｒ ｅ ｄ ａｎ ｇ ｎ ｔ ｒ ｗｅ ｔｋｅ ｉ ｔ ｎ ｉ ｌ ｈ ｗ ｅ ｒ ａ ｎ ｎｏ ｃ ｎ ｉ ｅａｉｎ．Ｂ ｓ ｄ ｏ ｅ ｕ ｄ ｒｙｎ ｕ ｆｃ ｏ ｄｔｏ ｓ ｎ ｌ ｄ ｎ ｅ ｔｐ ｇａ ｈ ｎ ｄ ｔ ｅｔｐ ｓｏ ｏｌａ ｄ ｖ ｇ ｔｔ ｎ ｏ ｓｄ ｒ ｔ ｏ ａ ｅ ｎ ｔ ｎ ｅ ｌｉ ｇｓ ｒａ ｅｃ ｎ ｉ ｎ ，ｉｃｕ ｉ ｇｔ ｏ ｏ ｒ ｐ ｙ，ａ ｙ ｅ ｆｓ ｉ ｎ ｅ ｅａｉ ｈ ｉ ｈ ｈ ｏ ｉ ｈ ｔ ｄ ａ ｉ ｎ ｔ ｅ ｓｕ ｙ ｂ ｎ，ａｍｅ ｏ ｏ ａａ ｔｒｅ ｔ ｔ ｎ ｏ ｈ ｄ ｅ ｓ ｒ ｐ ｓ ｄ ａ ｄ ｖ ｒｆ ｄ．Ｔｈ ｄ ｅ ｓ ｐ ｌ ｄ ｔ ｓ ｔｄ ｆｒｐ ｒｍｅｅ ｓｉ ｉ ｆｔｅ ｍｏ ｌＷａ ｐ ｏ ｅ ｎ ｅ ｉ ｈ ｍａ ｏ ｏ ｉ ｅ ｅ ｍｏ ｌＷａ ａ ｐ ｉ ｏ ｅ ｔｅ Ｔｕ ｘ ａ ｉ ｎ Ａｎ ｉＰｒｖｎ ｅ ｆ ｒｆｏ ｄ ｓｍｕａｉｎ ａ ｄ ｐ ｏ ｕ ｅ ａ ｏａ ｅ ｒｓ ｌ ｈ ｎ ｉＢ ｓｎ ｉ ｈｕ ｏ ｉｃ ｏ ｌｏ ｉ ｌｔ ｒｄ ｃ ｄ ｆｖ ｒｂｌ ｕｔ ｏ ｎ ｅ ｓ．
清华大学《高等水文学》L06-新安江模型_95380393
R = (P − E) − (W2 − W1)
(3)
大量资料表明，WWM～f/F有如下关系：
1 − f = (1 − WWM )B
F
WWMM
或 f =1 − (1 − WWM )B
F
WWMM
∫ = WM
1
= WWMd(f / F)
WM = WUM + WLM + WDM
W = WU + WL + WD
(12)
WUM,WLM,WDM----上层,下层,深层土壤蓄水容量
WU, WL, WD----上层,下层深, 层土壤蓄水量
当则WU > EM,
EU = EM
EL = 0
(13)
当则WU < EM,
EU = WU
EL = (EM-EU) WL WLM
二层模型：WU>EM EU=EM WU<=EM EU=WU EL=(EM-EU)WL/ WLM
二层模型在久旱之后，WL已很小，如仍不下雨，计算出 的蒸发EL会很小，这与实际不符。因为由于植物根系的 作用，深层的水分会通过植物腾发到大气。故宜采用三 层模型。
16
三层模型：
模型中流域蓄水容量和WM流域蓄水量都是W上层下层和,深层之和即 , :
n+1 j +1
=
C1Q
n j
+
C2Q
n+1 j
+
C3Q
n j
+1
C1
=
kx + 0.5∆t k(1− x) + 0.5∆t
新安江模型参数的分析
一、模型的结构与参数三水源新安江模型的流程图如图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 W M W MM A +--=当0≤⨯-EM K P ，则R=0不然，则当MM A EM K P <+⨯-,B MM A EM K P W M W W M 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 ，出流进入河网。
新安江模型在资料匮乏的长江中下游山区中小流域洪水预报应用
ke SC)湖泊科学)，2021,33(2)：581-594DOI10.18307/2021.0223©2021by Jofnal of Lake Sciences新安江模型在资料匮乏的长江中下游山区中小流域洪水预报应用**龚f夫1**，陈红兵2,朱芳2,李永凯2,崔明2,李致家1(1:河海大学水文水资源学院#南京210098)(2：宜昌市水文水资源勘测局，宜昌443000)摘要：水文资料匮乏流域的洪水预报(PUBs)是水文科学与工程中一个尚未解决的重大挑战.中国湿润山区中小流域大多是水文资料匮乏的流域，在此地区进行洪水预报的重要手段之一就是水文模型参数的估计•对基于参数物理意义的估算方法(以下简称物理估算法)及两种区域化方法进行了研究，将其用于新安江模型参数的估算及移植•皖南山区的29个中小流域被选作水文资料丰富的测量流域，鄂西山区的3个中小流域被视为水文资料匮乏的目标流域#目的是研究目标流域与测量流域空间位置较远但物理条件相似时，区域化等方法是否可以有效估计模型参数.结果表明，即使目标流域与测量流域空间距离较远，区域化及物理估算法也能一定程度上减少参数估计导致的模型效率损失，且在研究区的最优参数估计方案为单流域物理相似法结合回归法及物理估算法.为长江中下游资料匮乏的山区中小流域提出了可行的新安江模型参数估计方案，为该地区的洪水预报提供指导•关键词：长江中下游;山区中小流域&无资料;新安江模型;参数估计;区域化Applicahon of Xio'njiong Modei io the flow pradiction of ungauged smal i-and mediumsized catchments io the middle and lower raaches of the Yangtze River Basio*Gong Junfu1**,Chen Hongbing2,Zhu Fang2,Li Yongkai2,Cui Ming2&Li Zhijia1(1:College of Hydrologd and Water Resources#Hohai University#Nanjing210098#P.R.China)(2：Hydrologd and Water Resources Surved Bureru of YicCang#YicCang443000#P.R.Rhina)Abstrach:FU w prediction in ungauged basins(PUB)J as impoUant as it J ddicult.PUB J a major unsolved chW U nge in scientific and engineeJng hydrology.So far,empiwcal and stochastic methods have mainly been used for PJ pujose.Many smal l catchments in humid mountainous areas of Chine are with little os no hydrological data.The principal method of Jood forecasting in this area is the parameter estimation of the hydrological model.Xin'anjiang Model was used to predict Jow in this study,and Xin'anjiang Model parameters were estimated based on tUe physical meaning of Pem and regionabzationi Regression and physical similaJty approach were included in regionalization metUods in this study.The29small basins in mountainous areas in southern Anhui Province were selected for gauged catchments,and three small basins in mountainous regions in westem Hubei Province were treated as ungauged catchments.The main goal of this approach was to research whePes regionalization and oPes methods couad e o e ctoeeayestomatemodeapaameteswhen thespatoaaaocatoon oothegauged catchmentsand theungauged catchmentsweCe far away,but the physical conditions were similes.It is shown that:(1)Xin'anjiang Model has high accuracy in Jow prediction of study areas•(2)Even J the space distance of gauged catchments and ungauged catchments J far,regionalization and physicalbased estimation method can also reduce the e/iJency loss caused by parameter es/mabon to some extent.(3)WiP the increase in the number of donor catchments,the etiJency of the physical similaJty approach J reduced.(4)The optimal parameter estimation scheme in the study area is Pe ingenious combination of physical similaJty approach(only one donor catchment),regression approach and physical-based estimation method.A feasible parameter estimation scheme of the Xin'anjiang Model J proposed in this study,which provides guidance for Jow forecasting in smal l ungauged catchments in humid mounWinous areas among Pe middle and Uwes reaches of the Yangtze RJeSi*2019-10-10收稿&2020-06-23收修改稿•国家自然科学基金项目(51679061)和国家重点研发计划项目(2018YFC1508103)联合资助.**通信作者;E-mail:gon/junfu321@.ke Sd(湖泊科学)，2021,33(2)Keywords：The middle and Uwes reaches of the Yangtze Rives；smal l and medium catchmenb in the mountainous area;ungauged；Xin'anjiang Model¶meter estimation®ionalization可靠的洪水预报是防洪决策中的重要依据［1］，是减少甚至避免洪水破坏的第一道防线.在水文学研究及实际工程应用中,相较于水文资料充足的地区，更常见且更有挑战性的是对资料匮乏流域进行洪水预报.国际水文学会（IAHS）曾在20世纪初提出了水文资料匮乏流域的洪水预报（PUBs）10年计划［2］,期间涌现了大量研究成果+ 3A,-资料匮乏流域的洪水预报研究一般被概括为以下两部分［6］:（1）直接针对径流或径流指标进行区域化；！2）针对水文模型参数进行区域化.两者都是基于回归方法或目标流域与测量流域之间的某种距离度量.首先，简单的回归模型显示出较高的应用价值，如Maz/weei等在非洲52个流域建立了径流指标与流域属性的多元线性回归模型，以期通过年平均降雨、流域平均坡度、河网密度和土地覆盖等流域属性估计无资料流域的年平均径流量、基流指数、历时曲线和月平均径流量等径流指标.Marsh—等⑷基于英国土壤类型水文计划（HOST）的土壤数据，使用了一个简单的回归方法较精确地预测了未测量流域的基流指数等径流指标.也有部分研究基于更复杂的回归模型，如Laahe等［10-11］揭示了根据区域特点定制回归模型而不是使用全局模型的好处，其研究表明基于按季节性分组的流域回归模型提供了鲁棒性更强的区域化结果.同样地,地质统计学方法也被用于估计资料匮乏流域的水文变量，比如Sk~ien等阳提出了拓扑克里金法，可以用于解释水动力及地貌特征的空间离散性，该方法避免了输入数据的误差和参数可识别性问题.部分研究［13］认为该方法提供了比回归模型更稳健的估计.另一方面，研究者们越来越认识到，流域空间上的接近不一定导致流域产流模式的相似［14］，当使用更具有水文学意义上的距离量度时，基于距离的方法的估计效果会显著提高［15］-如Mew等［16］将拓扑克里金法与流域特征相结合，提高了该方法的预测性能. Domes等®认为只要保证流域之间的水文相似性，水文模型的参数甚至可以传递数千千米.时至今日，资料匮乏地区的水文预报获得了长足进步，但仍有大量难题无法解决［18］-我国有防洪任务的重点中小河流（汇流面积小于3000km2）有5186条问.近年来，受气候变化影响，由局地强降水造成的中小河流突发性洪水频繁发生，已成为造成人员伤亡的主要灾种.据统计,我国中小流域涝灾害和山洪地质灾害损失约占全国洪涝灾害经济损失的70%~80%，对人民群众的生命财产安全构成了严重的威胁［20］-由于大部分中小河流源短流急，洪水具有历时短、上涨快的特点，加之水文资料短缺，传统的预报方案通常难以对这些地区的洪水做出有效的预报［21］-区域化方法是解决此问题的思路之一，目前针对无资料中小流域的区域化研究较少,姚成等［22］使用区域化方法（空间临近法和回归法）将API模型和新安江模型应用于大别山区及皖南山区的29个中小流域，得到了较高的预报精度.但是在现有的区域化研究中，测量流域与目标流域大多空间分布均匀且空间距离较近.在这些地区,基于空间距离的区域化方法（以下称为空间临近法）由于操作简单且所需资料少而多被使用.同时，基于水文属性距离的区域化方法（以下称为物理相似法）和回归方法也被用来提高预测性能.在自然条件与历史因素的限制下［23］，设立在我国湿润山区中小流域的水文测站分布稀疏,可能导致某一具体区域内的测量流域数量较少.是否可以通过其他资料充足但相隔较远地区的测量流域作为参证流域进行参数估计？为了回答此问题,本研究关注于目标流域与测量流域空间位置较远时,区域化方法是否可以有效估计参数.此时，测量流域与目标流域之间显著的空间差异导致空间临近法无法使用，物理相似法及回归方法是否适用，需要接下来进一步的研究.本研究将皖南山区与鄂西山区分别作为已观测区及未观测区，选用新安江模型进行洪水模拟，并验证区域化方法（物理相似法及回归法）在研究区的适用性-正如PUetc等［24］所言，单独使用参数区域化有诸多弊端，本研究基于新安江模型参数的物理意义先估算出部分参数，剩余参数使用区域化方法（物理相似法及回归法）估计,并比较不同参数估计方案的效果，提出在研究区内最为适用的参数估计方案，为长江中下游资料匮乏的山区中小流域水预供参.1研究区概况及资料来源研究区域为长江中下游地区的32个山区中小流域,其中位于皖南山区的29个流域为有观测资料的测量流域，鄂西山区的3个流域被视为资料匮乏的目标流域，研究区域见图1.各流域面积均不超过1000km2,龚f夫等：新安江模型在资料匮乏的长江中下游山区中小流域洪水预报应用583均为自然闭合流域，受水工建筑物的影响较小，流域平均海拔多在400-1000m,其中雾渡河流域平均海拔最高，为1109mi各流域年平均降水量均在1200-1500m叫年平均蒸发量均在600-800m叫年平均气温14-1601研究流域大多地势，密集，平均坡度大在10-30°,密度大多在00~10km/km2.研究区内型和黏土为主，饱和渗透系数大多在3~12mm/h.土地利用（LUCC）以林地为主，大流域林地面积占比70%以上，甚至达到100%i研究区内植被类型为主，大为常绿针叶林、常阔叶林及灌木i图1研究区地理位置FigO Locations of the research catchments降雨、径流及蒸发资料来自当地水文局，时间步长均为1hi29个测量流域均具有至少30年的观测资料.3个目标流域中，有20082017间共11场观水，有20162018年间共5场观水，下仅有2016共2场观测洪水.DEM来自中国科学院地理空间数，分辨率为30m，用取流域地形、地貌属性;土地利用数据来自中国科学院遥感的2010年LUCC图，分辨率为1km，用取流域土地利用属性；数自于HWSD（世界数）的型图，分辨率为1km，用取流域土属.2研究的研究思路大致如下：在量流域进行新安型的校准，得到量流域的新安型校参数集.然流域物理属性的筛选，筛的流域物理属，运用区域化（回归法、物理相似法），将量流域的校准参数转移至目标流域.与此同时，使用参数物理意义的（以下简称物理），参数•最后，将理与区域化结合，针目标流域的参数案，并不同参数案的效果，得到在研究区内最为适用的参数案•本研究的线图图％22新安江模型及物理估算法2.1.1新安江模型新安江模型是赵人俊+25］于1984年前后提岀的概念性模型,在我国湿润及半湿润地区具有较高的预报精度，已到广泛应用+ 26气巴该模型共有16参数,各参数物理意［29］，其中不敏感ke SO.(湖泊科学)，2021,33(2)参数10个(WM、WUM、WLM、B、C、CP、EX、CC、KE、XE$，敏感参数7个(KC、SM、KI(KG)、CG、CS、L)-根据赵人俊等[3°]提岀的客观优选法，不敏感参数不参与率定，直接赋予固定值•敏感参数的率定使用人工优选结合SCE-NA算法[31],这样既可以快速稳定的优选参数，又能有效避免自动优选参数可能不合理的问题•选用洪峰相对误差、峰现时间误差,纳什效率系数(NSE)对参数率定结果进行评价.依据《水文情报预报规范》两，洪峰相对误差在20%以内，峰现时间误差在3h以内为合格-在目标流域，不敏感参数同样直接赋予固定值，敏感参数需要进行估计.在敏感参数中K与KG满足线性约束，只需估计K即可,因此本研究在目标流域需要估计的参数为6个，见表1-在新安江模型的待估计参数中，壤中流岀流系数K和河网水消退系数CS可以基于参数的物理意义进行估算(物理估算法)-表1目标流域待估计参数Tab.1Estimated parameters of ungauged basins参数参数意义参数取值范围KC蒸散发折算系数一般小于1SM自由水容量一般为10〜50mm,在土层较厚、透水性强的流域可取50mm或更大K/壤中流出流系数0〜06之间CG地下水消退系数0〜1之间CS水流系数0〜1之间L河网汇流滞时大于等于0,在山区中小流域一般小于10h276壤中流与地下水出流系数KC与KG土壤质地可以用来推断土壤的性质，进而估算自由蓄水量对壤中流及地下水的岀流系数(新安江模型中的K与KG).K和KG之和决定了自由水的排水速率，赵人俊等认为KI+KG的值可固定为0.5&KI/KG决定了壤中流与地下水岀流的比例，姚成等⑶〕在GXM模型中提岀了针对栅格的K/KG估计公式，将其修改成适用于新安江模型的形式：KC KG=0"K/=1+2(1_九)(1)KG表2不同质地土壤凋萎含水量Tab.2Wilting moisture content of various soil textures 土壤质地/p土壤质地&呷砂土0.04砂质黏壤土0.16壤砂土0.05粉质黏壤土0.24砂0.09黏0.21粉壤土0.15砂质黏0.21粉土0.11粉黏0.28壤土0.14黏0.28式中，m c为自由水岀流校正系数，根据前人的研究成果a35〕，在湿润流域m c可取为1;/为土壤凋萎含水量，以体积含水量表示，可根据土壤质地取值#见表2[36]-2.1.3河网水消退系数CS李致家等⑶〕针对二叉树结构河网，对河链进行编号，链接河源的河链编号记为厂=1，下游河链的编号等于该条河链的上游总链数+1•并根据矩形河道线性蓄泄关系(式！2))#引入一个无量纲时间i+：vQ(ftt=/(ft(2)式中，Q(ft为河链厂的岀流量,”3/s；W(ft)为河链厂的蓄量,m3&v为平均河道流速，”/s；L为平均河链长，”i在资料匮乏地区一般没有实测流速资料，可以使用经验公式估算.RodOguez-tturbe等+38〕基于运动波理论提岀如下流速公式：v=0.665-弓(p/A)(4)龚f夫等：新安江模型在资料匮乏的长江中下游山区中小流域洪水预报应用585tU—2\J丿o.D式中，P为次洪的平均雨强,c m/h；2为流域面积,km2;D为平均河宽,m；S。
基于栅格的分布式新安江模型构建与分析(硕士论文)
关键词：数字高程模型、新安江模型、Grid-Xinanjiang 模型、TOPMODEL 模型、GTOPMODEL 模型、密赛流域、洪水模拟
Abstract
无 河海大学
论文中文题名
基于栅格的分布式新安江模型构建与分析
论文中文副题名
无
论文英文题名
Development and Application of Grid-Based Distributed Xinanjiang Model
论文英文副题名
None
论文语种 汉语 论文摘要语种 汉、英 论文页数 74论文字数 3.1 (万)
Grid-Xinanjiang model takes every DEM unit grid as a sub-basin, using the initial Xinanjiang model to calculate runoff-yield of the current grid, then, according to the calculation order among the DEM grids, the flow of the outlet of each grid to the outlet of the whole basin is routed by the grid-by-grid method of Muskingum. When calculating the runoff-yield and runoff-concentration in every unit gird, the water balance and influence of the river drainage network are taken into consideration. In other words, in the calculation of the runoff-concentration, first, considering whether the current grid is saturated. If it is not, the soil water content of this grid will increases due to the outflow of the upstream grids before the soil moisture of it reaches water storage capacity. Moreover, if there is a channel in unit grids and because of the effect of the drainage network, not all the overland inflow and interflow flow into the downstream grid, part of them will flow into the channel.
新安江水文模型简介
《流域水文模拟》结课报告新安江模型的原理、结构及应用、发展历程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 新安江模型起源 ............................... 错误!未定义书签。
第三章新安江模型
（2）三水源
由于饱和坡面流的产流面积是不断变化的，所以在
产流面积上自由水蓄水容量分布是不均匀的。三水源
水源划分结构是采用类似于流域蓄水容量面积分配曲
线的流域自由水蓄水容量面积分配曲线来考虑流域内
自由水蓄水容量分布不均匀的问题。所谓流域自由水
蓄水容量面积分配曲线是指：部分产流面积随自由水
蓄水容量而变化的累计频率曲线。流域自由水蓄水容
马斯京根 或
滞后演算 法
SM、EX、KG 、KI
UH或 CS、CI、CG
KE、XE或L
3.2 模型计算 1 蒸散发计算
具体计算为蒸散发计算采用三层模型，其参数有上 层张力水蓄水容量UM，下层张力水蓄水容量LM，深层张 力水蓄水容量DM，流域平均张力水蓄水容量WM，蒸散发 折算系数KC，深层蒸散发扩散系数C，计算公式为：
新安江模型各层次结构功能、计算方法和相应参数
层次 第一层次 第二层次
第三层次
第四层次
功能
蒸散发 计算
产流 计算
水源 划分
二水源 三水源
汇流 计算
坡面汇流 河道汇流
计算 方法
三层 模型
蓄满 产流
参数
KC、UM、 LM、DM、
C
WM、B、IM
稳定 下渗
fC
自由水 蓄水库
单位线或 线性水库或 滞后演算法
存在的主要问题：
③对许多流域资料的分析表明，即使是同一流域，
各次洪水所分析出的 fc也不尽相同，而且有的时候
变化还很大，很难进行地区综合和在时空上外延， 应用时其任意性大，常造成较大误差。
（2）三水源
三水源的水源划分结构应用了 山坡水文学的概念，去掉了 fc， 用自由水蓄水库结构解决水源 划分问题。
扩散波与马斯京根法在栅格汇流演算中的应用比较
扩散波与马斯京根法在栅格汇流演算中的应用比较姚成;章玉霞;李致家【摘要】以数字高程模型(DEM)为基础,对分布式水文模型中的汇流演算方法进行研究.以皖南山区呈村流域为例,分别采用扩散波模型与马斯京根法进行栅格单元间的汇流过程模拟,对2种汇流演算方法进行应用比较.结果表明,采用2种方法模拟的洪水过程精度基本相当,均能取得较好的应用效果.通过进一步对比分析后发现,两者在描述参数与变量的空间分布、运算效率等方面存在一定的差异.%The flow routing method in the distributed hydrological model was investigated based on the digital elevation model ( DEM). The cell-to-cell diffusion wave and Muskingum routing methods were applied to the flood simulation in the Chengcun watershed, located in the southern mountainous area of Anhui Province. The results show that the performances of the two methods are comparable and can produce desirable simulation results. Further comparison and analysis show that the two methods have differences in representation of the spatial distribution of parameters and variables, as well as their computational efficiency.【期刊名称】《河海大学学报（自然科学版）》【年(卷),期】2013(041)001【总页数】5页(P6-10)【关键词】分布式水文模型;汇流演算;扩散波;马斯京根法;洪水模拟;呈村流域【作者】姚成;章玉霞;李致家【作者单位】河海大学水文水资源学院,江苏南京210098;河海大学水文水资源学院,江苏南京210098;河海大学水文水资源学院,江苏南京210098【正文语种】中文【中图分类】P338水文模型作为一种研究流域内复杂水文现象的有效工具，在水文学研究领域中一直发挥着重要作用。
新安江模型的应用
新安江模型的应用张利茹河海大学水文水资源学院，南京（210098）摘要：新安江降雨径流模型应用在梁辉水库上，采用2002年至2006年五年的降雨和蒸发资料对该流域进行日模和次模的模拟，得出的结果还比较满意。
为了找出新安江模型的敏感性参数，本文在其他研究人员的基础上，选出公认的比较敏感的参数，把它们的值分别变成初始值的80%、90%和110%（CG除外）后进行模拟计算，得出的结果证实了学者们的说法。
关键词：新安江模型，梁辉水库，敏感性分析1. 新安江模型简介新安江模型始建于1973年，采用蓄满产流的概念，以土壤含水量达到田间持水量后才产流，是个分布式的概念性模型，30多年来在我国湿润与半湿润地区有广泛应用，并发展改进为三水源的以及其他多水源的模型【1】。
几十年来，很多专家和学者都致力于新安江模型的应用和发展上，发表了数以百计篇文章(像赵仁俊，1992；程等人，2002)，但很少有用一个实际例子来研究新安江模型参数的敏感性问题的，实际上，新安江模型参数的命感性分析会有助于该模型的更广泛的应用，例如，对于无资料的地区或是资料不全的地区，参数的敏感性分析将显得更加有用。
2. 新安江模型结构新安江模型是分散性的模型，常按泰森多边形法把全流域分成许多单元流域，产流部分采用蓄满产流模型，另增加了流域不透水面积占全流域面积之比的参数IMP。
蒸发部分采用三水源蒸散发模式。
河道洪水演算采用马斯京根法。
地面径流的汇流采用经验单位线，并假定每个单元流域上的无因次单位线相同，简化结构。
地下径流的汇流采用线性水库。
对每一个单元流域作汇流计算，求得单元流域出口流量过程。
再进行出口以下的河道洪水演算，得出流域出口的流量过程。
把每个单元流域的出流过程相加，就求得了流域出口的总出流过程[2]。
新安江模型流程图如图1。
基于概念型降雨径流蓄满产流的新安江模型，其参数可大致划分为四种类型，如下述：（1）蒸散发。
此部分的参数包括K、C、WUM、WLM。
河海大学教师获得2020年度中国航海学会科学技术进步特等奖
52第1期姚 成,等 基于栅格型新安江模型的中小河流精细化洪水预报[25]李致家,黄鹏年,张建中,等.新安江海河模型的构建与应用[J].河海大学学报(自然科学版),2013,41(3):5⁃11.(LI Zhijia,HUANG Pengnian,ZHANG Jianzhong,et al.Construction and application of Xin’anjiang⁃Haihe model[J].Journal of Hohai University(Natural Sciences),2013,41(3):5⁃11.(in Chinese))[26]姚成,纪益秋,李致家,等.栅格型新安江模型的参数估计及应用[J].河海大学学报(自然科学版),2012,40(1):42⁃47.(YAO Cheng,JI Yiqiu,LI Zhijia,et al.Parameter estimation and application of grid⁃based Xin’anjiang model[J].Journal of Hohai University(Natural Sciences),2012,40(1):42⁃47.(in Chinese))[27]姚成,李致家,刘晋,等.无资料地区洪水预报模型研究与应用[M].南京:河海大学出版社,2018.[28]YAO Cheng,YE Jinyin,HE Zhixin,et al.Evaluation of flood prediction capability of the distributed Grid⁃Xin’anjiang modeldriven by weather research and forecasting precipitation[J].Journal of Flood Risk Management,2019,12:e12544. [29]姚成,李致家,章玉霞.DEM分辨率对分布式水文模拟的影响[J].水利水电科技进展,2013,33(5):11⁃14.(YAOCheng,LI Zhijia,ZHANG Yuxia.Effects of DEM resolution on distributed hydrologic simulation[J].Advances in Science and Technology of Water Resources,2013,33(5):11⁃14.(in Chinese)))(收稿日期:20201120 编辑:王 芳 ㊃简讯㊃河海大学教师获得2020年度中国航海学会科学技术进步特等奖2020年12月25日,2020年度中国航海学会科学技术进步奖揭晓,河海大学作为主要完成单位参与的‘基于波流边界层特性的海湾水沙运动模拟关键技术研究与应用“获评特等奖㊂该成果由水利部交通运输部国家能源局南京水利科学研究院㊁河海大学㊁中国科学院力学研究所㊁天津大学㊁交通运输部天津水运工程科学研究所等单位共同研究完成㊂河海大学龚政教授㊁姚鹏副教授㊁苏敏副研究员㊁张茜博士后作为技术骨干,共同参与科研攻关㊂该成果从波流边界层角度出发,阐明和揭示水沙动力过程的内在物理机制,进而推广应用到宏观动力地貌模拟技术,形成了基于波流边界层特性的海湾水沙运动模拟关键技术研究与应用的创新成果㊂该成果已推广应用到渤海湾曹妃甸㊁珠江口㊁江苏沿海㊁泉州港㊁瓯江口等港口航道㊁桥隧基槽㊁滩涂保护㊁河口治理等多个方面,对于促进水沙科学与港口航道工程学等发展具有重要的理论意义和实践价值,产生了显著的社会经济效益㊂中国航海学会科学技术进步奖设立于2002年,由中国航海学会主办,是中国水运行业领域高级别的科技奖项,对促进我国水运科技发展具有重要意义㊂(本刊编辑部供稿)。
新安江模型参数的局部灵敏度分析
新安江模型参数的局部灵敏度分析
薄会娟;董晓华;邓霞
【期刊名称】《人民长江》
【年(卷),期】2010(041)001
【摘要】新安江模型是我国应用广泛的一种降雨径流流域模型,适用于湿润和半湿润地区.新安江模型参数较多,人工率定时需要首先率定对模拟结果影响最为敏感的参数,因此对参数进行灵敏度分析是必要的.采用水文模拟法和摩尔斯分类筛选法,对清江流域1989年9月30日至10月5日次洪水进行径流模拟,并让参数变化一定的步长,用Nash系数对模拟结果的精度进行评定,从而反映出参数的灵敏度.结果表明,最灵敏的参数是FC和Cg;Kc和WM是灵敏参数;b对峰值流量是高灵敏参数,而对径流总量来说是灵敏参数;WUM,WLM,WDM,C和IMP是不灵敏参数.
【总页数】4页(P25-28)
【作者】薄会娟;董晓华;邓霞
【作者单位】三峡大学,土木水电学院,湖北,宜昌,443002;三峡大学,土木水电学院,湖北,宜昌,443002;三峡大学,土木水电学院,湖北,宜昌,443002
【正文语种】中文
【中图分类】TV131.4;P338.9
【相关文献】
1.SWMM模型中参数率定及局部灵敏度分析 [J], 官奕宏;吕谋;王烨;杨婷婷
2.高承压含水层水文地质参数局部和全局灵敏度分析 [J], 赵春虎
3.城市降雨径流模型的参数局部灵敏度分析 [J], 黄金良;杜鹏飞;何万谦;欧志丹;王浩昌;王志石
4.新安江模型系统响应参数优化方法在邵武流域的应用 [J], 赵丽平;邢西刚;吴楠;刘云;周艳先
5.新安江模型参数线性化率定研究 [J], 王红艳;王新龙;周倩;杨宇
因版权原因，仅展示原文概要，查看原文内容请购买。
新安江流域水文模型
世界上第一个流域水文模型－Stanford模型出现在20世纪60年代，目前全世界已提出数 以百计的流域水文模型。主要包括由美国天气局V. T. Sitten提出的API模型、N. H. Crawford和R. K. Linsley提出的斯坦福模型以及R. J. C. Bernash等提出的萨克拉门托模 型，日本国立防灾科学研究中心菅原正已教授提出的水箱模型，丹麦技术大学提出的 NAM模型，以及原华东水利学院赵人俊教授提出的新安江模型。这些概念性水文模型 对流域的降雨径流过程进行了较为细致的模拟。由于这些模型具有较好的结构形式和 良好的模拟预报精度，因此在洪水实时预报中得到广泛地应用。本文主要介绍国内应 用最为广泛的新安江三水源模型。
－2002.12－
新安江模型的结构
蒸散发计算原理
各层蒸散发的计算原则是，上层按蒸散发能力蒸发，上层含水量蒸发量不够 蒸发时，剩余蒸散发能力从下层蒸发，下层蒸发与蒸散发能力及下层含水量 成正比，与下层蓄水容量成反比。要求计算的下层蒸发量与剩余蒸散发能力 之比不小于深层蒸散发系数。否则，不足部分由下层含水量补给,当下层水量 不够补给时，用深层含水量补充。(上层以蒸散发能力蒸发,直到上层水分耗尽, 才蒸发下层;下层土壤蒸散发量与剩余蒸散发能力(流域蒸散发能力与上层蒸散 发量之差)及下层土壤实际含水量成正比。)
－2002.12－
流域蓄水容量曲线
新安江模型的结构－产流量计算
－2002.12－
新安江模型的结构－产流量计算
据此可求得流域平均蓄水容量为
WM
Wmm (1
0
f F
)dWm
Wmm B 1
与流域初始平均蓄水量W0 相应的纵坐标A(
)为 A
基于SoilGrids的栅格新安江模型参数空间分布估算
基于SoilGrids的栅格新安江模型参数空间分布估算
童冰星;李致家;姚成
【期刊名称】《水科学进展》
【年(卷),期】2022(33)2
【摘要】为实现中小流域降雨径流过程精细化模拟,合理估算水文模型参数的空间分布具有重要意义。
基于新版全球数字土壤制图系统(SoilGrids)构建栅格新安江模型(GXM)参数化方案,对陕西省陈河流域2003—2012年16场洪水进行模拟,与新安江模型计算结果进行对比,开展基于洪水过程划分的自由水蓄水容量敏感性及空间分布特征量化分析。
结果表明:GXM模拟的峰现时间误差水平降低约0.31 h,洪峰和洪量模拟精度较高,模型能够对土壤水饱和度等水文要素的动态空间分布进行较合理的模拟;自由水蓄水容量参数对洪峰和涨洪过程的确定性系数以及涨洪段的洪量相对误差影响较大,对退水过程影响小;自由水蓄水容量在陈河流域河谷和山脊附近较大,坡段中部较小。
【总页数】8页(P219-226)
【作者】童冰星;李致家;姚成
【作者单位】河海大学水文水资源学院
【正文语种】中文
【中图分类】TV122
【相关文献】
1.基于栅格的新安江模型与GTOPMODEL模型
2.基于分形理论的新安江模型参数空间尺度效应分析
3.基于分形理论的新安江模型参数空间尺度效应分析
4.栅格新安江-地表地下双人工调蓄分布式水文模型
5.基于栅格的分布式新安江模型在三岔河流域的洪水预报研究
因版权原因，仅展示原文概要，查看原文内容请购买。
《河海大学学报(自然科学版)》被Scopus数据库收录
（ Ｄ Ａ Ｉ Ｊ ｉ ａ ｎ ｎ ａ ｎ ， Ｌ Ｉ Ｚ ｈ ｉ ｊ ｉ ａ ，ＨＵ Ａ Ｎ Ｇ Ｐ ｅ ｎ ｇ ｎ ｉ ａ ｎ ，ｅ ｔ １．Ｕ ａ ｎ ｃ ｅ ｒ ｔ ａ ｉ ｎ ｔ ｙ ａ ｎ ｌ ａ ｙ ｓ ｉ ｓ ｏ ｆ Ｘ ｉ ｎ’ ｎｊ ａ ｉ ａ ｎ ｇ ｍ ｏ ｄ ｅ ｌ ｐ ａ ｒ ａ ｍ ｅ ｔ ｅ ｒ ｓ［ Ｊ ］ ．Ｊ ｏ ｕ ｎａ ｒ ｌ ｏ ｆ
ｓ ｉ ｍ ｐ ｌ ｅ ｘ ｍ ｅ ｔ ｈ ｏ ｄ［ Ｊ ］ ．Ｅ ｎ ｇ ｉ ｎ ｅ ｅ ｉ ｒ ｎ ｇ Ｊ ｏ ｕ ｒ ｎ ｌ ａ ｏ ｆ Ｗｕ ｈ ａ ｎ Ｕ ｎ ｉ ｖ ｅ ｒ ｓ ｉ ｔ ｙ ， ２ ０ ０ ９ ， ４ ２ （ １ ） ： ６ － ９ ． （ ｉ ｎ Ｃ ｈ ｉ ｎ ｅ ｓ ｅ ） ）
（ Ｙ Ａ Ｏ Ｃ ｈ ｅ ｎ ｇ ， Ｊ ｌ Ｙ ｉ ｑ ｉ ｕ ， Ｌ Ｉ Ｚ ｈ ｉ ｊ ｉ ａ ， ｅ ｔ １．Ｐ ａ ａ ｒ ａ ｍ ｅ ｔ ｅ ｒ ｅ ｓ ｔ ｉ ｍａ ｔ ｉ ｏ ｎ ａ ｎ ｄ印ｐ ｌ ｉ ｃ ａ ｔ ｉ ｏ ｎ ｏ ｆ ｄ — ｂ ａ ｓ ｅ ｄ Ｘ ｉ ｎ ’ ａ ｎ ｊ ｉ ａ ｎ ｇ ｍ ｏ ｄ ｅ ｌ［ Ｊ ］ ．Ｊ ｏ ｕ ｒ ｎ ａ ｌ ｏ ｆ
Ｈ ｏ ｈ ａ ｉ Ｕ ｎ ｉ ｖ ｅ ｒ ｓ ｉ ｔ ｙ ： Ｎ ａ ｔ ｕ ｒ ｌ ａ Ｓ ｃ ｉ ｅ ｎ ｃ ｅ ｓ ， ２ ０ １ ２ ， ４ ０ （ １ ） ： ４ ２ － ４ ７ ． （ ｉ ｎ Ｃ ｈ ｉ ｎ ｅ ｓ ｅ ） ）
基于栅格的新安江模型的构建和应用
基于栅格的新安江模型的构建和应用
李致家;姚成;汪中华
【期刊名称】《河海大学学报（自然科学版）》
【年(卷),期】2007(035)002
【摘要】根据分布式模型和新安江模型的最新理论,建立了基于栅格的分布式新安江水文模型,该模型能够考虑网格之间的水量交换.将建立的模型用于钱塘江支流密赛流域的洪水模拟计算,并与新安江模型的计算结果进行比较,可以看出,基于栅格的分布式新安江水文模型可以很好地应用于该流域进行洪水模拟,其模拟精度于原新安江模型的精度相当.
【总页数】4页(P131-134)
【作者】李致家;姚成;汪中华
【作者单位】河海大学水文水资源学院,江苏,南京,210098;河海大学水文水资源学院,江苏,南京,210098;济宁市水文水资源勘测局,山东,济宁,272045
【正文语种】中文
【中图分类】P333
【相关文献】
1.栅格新安江模型在陕西省大河坝流域中的构建与应用 [J], 贺成民;顾钊
2.基于栅格型新安江模型的中小河流精细化洪水预报 [J], 姚成;李致家;张珂;朱跃龙;刘志雨;黄迎春;龚珺夫;张锦堂;童冰星
3.基于栅格的新安江模型在滦河流域的构建与应用 [J], 厉治平;范辉;徐嘉
4.基于栅格新安江模型的秦淮河流域洪水模拟及实时校正研究 [J], 孙文宇;王容;姚
成;何健;李致家;龚珺夫;黄鹏年
5.基于栅格的分布式新安江模型在三岔河流域的洪水预报研究 [J], 范火生
因版权原因，仅展示原文概要，查看原文内容请购买。
网格新安江模型在于桥水库洪水预报中的应用
网格新安江模型在于桥水库洪水预报中的应用赵考生【摘要】根据于桥水库流域的地理概况和水文气象特性,利用网格新安江模型GRID-XAJ对流域历史资料进行降雨径流模拟及参数率定,研制洪水预报方案.从模拟结果上看,网格新安江模型能有效反映流域下垫面变化及降水空间分布的离散对降水径流的影响.【期刊名称】《水科学与工程技术》【年(卷),期】2017(000)005【总页数】4页(P30-33)【关键词】于桥水库;洪水预报;网格新安江模型;参数率定【作者】赵考生【作者单位】天津市水文水资源勘测管理中心,天津300061【正文语种】中文【中图分类】TV122根据于桥水库流域的地理概况和水文气象特性，考虑人类活动对下垫面的影响，应用ArcView对于桥流域的数字高程模型DEM进行识别和相应的处理，选择网格新安江模型GRID-XAJ对流域历史资料进行降雨径流模拟及参数率定，模拟各子流域降雨、洪水过程，预报全流域暴雨之后的入库水量、水位、流量及峰现时间等，建立洪水预报方案，提高洪水预报精度。
水文资料主要包括降水量、径流量、水库水位及流域地下水位等，资料年限为1960～2012年。
本研究采用于桥水库的库面蒸发资料作为于桥水库流域的蒸发值处理，以流域的DEM数据为基础，根据生成的集水面积和水系流向，客观充分利用地形资料，并结合人工经验、流域内雨量站、水文站分布情况等划分自然子流域。
于桥水库流域的控制面积2060km2，分为前毛庄、水平口、淋河桥、水库区间4个一级子流域，其中每个一级子流域再用自然子流域方法来划分二级子流域，划分原则是尽可能利用水文站作为子流域的控制断面，根据流域内雨量站和水文站的布设情况及自然流域的边界，用Arcgis软件把流域（包括于桥水库库区水域）划分为17个二级子流域。
采用1960～2012年各子流域内平均降水量、实测径流量（或反推入库流量）、天然径流量和用水量、流域面积等数据系列，做降雨径流相关图。
基于栅格的新安江模型的构建和应用
基金项目： 国家自 然科摹基金资助项目 ５４ ０７． （ ７１ ０９ ）
作者简介 ： 李致家（９ ２ ）男 ， １ 一 ， 山西运城 人 ， ６ 教授 ， 博士生导师 ， 主要从事水文模 型与水 文预报研 究
维普资讯
型能够考虑 网格之间的水量交换 ． 将建立的模型用于钱塘江支流密进 行 比较 ， 以看 出， 于栅格 的分布 式新 安 江水 文模 型 可以很好 地 应 用 于该 可 基
流域进 行 洪水模 拟 ， 其模 拟精 度 于 原新安 江模 型 的精 度相 当 ． 关键 词 ： 分布 式 水文模 型 ； 新安 江模 型 ； ｒ 一 ｉ ｊ Ｉ 模 型 ； Ｇ ｉ ）嗍 ｉ ｌ ｄ【 ａｇ 洪水模拟
１２ ３
河 海 大 学 学 报 （自 然 科 学 版 ）
１ 基 于 网格 的 新 安 江 模 型 构 建
１ １ 新安 江模 型 ．
新安江模型是一个分散参数的概念性模型． 对于较大流域 ， 根据流域下垫面的水文 、 地理情况将其分为
若 干 个单元 面 积 ， 每 个单元 面 积预 报 的流量 过程 演算 到 流域 出 口， 后叠 加起来 即为整 个流 域 的预报 流量 将 然 过程 ． 三水源新安江模型由 ４ 个模块组成 ， 分别为蒸散发计算 、 产流量计算 、 分水 源计算 、 汇流计算Ｅ９． ８３ － １２ ． 】瑚 ｉ 模 型 ！ 【 ｉ ａ ‘
匀以及水流 的横向运动 ， 因此单元流域的产流机制 属于蓄满产 流 有学者采用新安江 、Ｏ Ｍ Ｄ Ｌ Ｈ Ｃ ６． ＴＰ Ｏ Ｅ 和 Ｅ 这 ３个适 用 于湿 润地 区 的模型 对半 干旱 的海 河 流域 进 行模 拟 ， 果差 别 不 大 ， 明 在 目前 资 料 的 情 况下 ， 结 说 蓄 满产流一类 的模型可以扩展到半干旱地区【 新安江模型是一个半分布式的概念性模型 ， 引． 在中国的洪水预报
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
再考虑上游入流是否补足当前单元的土壤含水量，即可得到实际用于产流计算的时段雨量 (; ! 当 (;!? 或 (; B *?! *& 时，, > ?；当 (; B *? C *& 时，, > (; B *? D *& ! 其中：, 为时段产流量；*? 为栅格单元实际
的张力水含量 !
在 ’()*+,)-’.-/).-0 模型中，任意栅格单元内的产流量 , 均被划分为 1 种水源：地面径流 ,E、壤中流 ,) 和地下径流 ,0 ! 所需参数包括表层土自由水容量 %&、表层自由水含量对壤中流的出流系数 -)、表层自由水 含量对地下水的出流系数 -0 !
收稿日期：."!!!"!!!. 基金项目：中央高校基本科研业务费专项（."!"/""!!0）；国家博士点基金（.""-""-0!!"""1） 作者简介：姚成（!-%. —），男，江苏扬州人，讲师，博士，主要从事水文模型与水文预报研究 # 2345*6：7589:;,<= ::># ;?># 9,
第:期
将每个栅格单元内的土壤分为上层、下层和深层 1 层，每一层对应的张力水蓄水容量分别为 *@&，*A& 和 *$& ! 在栅格单元实际蒸散发计算时，冠层截留量按蒸散发能力蒸发，当截留水量小于蒸散发能力时，则 采用 1 层蒸散发模型［:?］!
!!. 单元产流及分水源
土壤蓄满表示的是土壤含水量达到田间持水量，而不是饱和含水量 ! 新安江模型引进张力水蓄水容量分
该模型以 C2D 栅格为计算单元，将每个单元的产流量划分为地表径流、壤中流以及地下径流 $ 种
水源，最后再根据栅格间的汇流演算次序，利用扩散波汇流方法依次将各种水源演算至流域出口 #
模型在进行产汇流计算时，考虑了栅格之间的水量交换以及河道排水网络的影响 # 以流域地貌特
征、土壤及植被类型等下垫面条件为基础，研究了模型参数的估计方法，并对其进行了验证 # 将模型
栅格型新安江模型的参数估计及应用
姚 成!，纪益秋.，李致家!，刘开磊!
（! # 河海大学水文水资源学院，江苏 南京 .!""-%；. # 南通市港闸区幸福街道办事处，江苏 南通 ..’"!.）
摘要：根据新安江模型与一维扩散波模型理论，构建了基于栅格的新安江（@A*?3B*,’5,)*5,<）模型 #
89" :;(-3!，<# 8+4=+>.，?# @;+47+%!，?#A B%+40(+! （! ! "#$$%&% #’ ()*+#$#&) ,-* .,/%+ 0%1#2+3%1，(#4,5 6-57%+15/)，8,-95-& .!""-%，"45-,；
. ! :5-&’2 ;/+%%/ <&%-3) 5- =,-&>4, ?51/+53/ #’ 8,-/#-& "5/)，8,-/#-& ..’"!.，"45-,）
假设在每一个河道栅格内，降雨分布均匀，河道形状不发生改变，则可以根据栅格单元内河道部分所占
的面积比例进行河道降水的计算 ! 河道降水 "5=可表示为
"5= #
)5= *5= ( + 0(
（2）
式中：) 5=———河道长度；* 5=———河道断面最大过水面积所对应的水面宽；+ 0(———栅格 单 元 的 面 积；( ———时 段降雨量 ! 其中，)5= *5=表示河道所占面积 ! 对于流域内的坡地栅格，即使有“虚拟河道”存在，也不考虑河道 降水的影响，即在坡地栅格上 "5= > ? ! !!- 蒸散发
! "#$%&’$(’)(*$)(+ 模型原理
’()*+,)-’.-/).-0 模型以流域内每个 $%& 栅格作为计算单元，并假设在栅格单元内的降雨和地貌特征、
土壤类型以及植被覆盖等下垫面条件空间分布均匀，模型只考虑各个要素在不同栅格之间的变异性 # ’()*+
,)-’.-/).-0 模型先计算出每个栅格单元的植被冠层截留量、河道降水量和蒸散发量，然后再计算出栅格单元
姚 成，等 栅格型新安江模型的参数估计及应用
F1
现实意义［!!"］# 本研究利用 $%&，提取了水文模拟与参数估计所需要的流域地貌特征，并根据新安江模型与 扩散波模型理论，构建了栅格型新安江（’()*+,)-’.-/).-0）模型 # 同时，本文结合土壤类型、植被覆盖及利用等 信息，对空间变化模型参数的估计方法进行了一定的研究与验证，并将 ’()*+,)-’.-/).-0 模型应用于安徽省屯 溪流域洪水模拟计算 #
用于安徽省屯溪流域的洪水模拟，取得了良好的应用效果 #
关键词：栅格型新安江模型；新安江模型；扩散波模型；参数估计；洪水模拟；屯溪流域
中图分类号：E$$$ F !
文献标志码：G
文章编号：!"""!!-%（" ."!.）"!!""0.!"’
$%&%’()(& (*)+’%)+,- %-. %//0+1%)+,- ,2 3&+.45%*(. 6+-’%-7+%-3 ’,.(0
的求解 !
日径流模拟时对汇流演算精度要求不高且时间步长较大，为了增加模型的运行效率，’()*+,)-’.-/).-0 模型采
用简便的线性水库法和滞后演算法进行栅格间汇流计算 #
!#! 植被冠层截留
植被冠层截留是指降雨在植被冠层表面的吸着力、承托力及水分重力、表面张力等作用下储存于其表面
的现象［3］! 在一次降雨过程中，植被冠层对降雨的累积截留量［4］可表示为
的产流量，并采用自由水蓄水库结构对其进行水源划分，即划分为地表径流、壤中流以及地下径流 1 种水源，
最后再根据栅格间的汇流演算次序，依次将各种水源演算至流域出口 # 在进行次洪模型计算时，采用一维扩
散波模型，且假设在原来的坡地栅格内也存在一个“虚拟河道”，每个栅格单元的壤中流和地下径流直接流入
河道，采用河道汇流的计算方法，则栅格间的扩散波汇流演算只包括坡面汇流和河道汇流 2 种方式 # 考虑到
B(C D,&.*：<A*?3L5+;? B*,’5,)*5,< 48?;6；B*,’5,)*5,< 48?;6；?*II>+*8, J5K; 48?;6；M5A54;H;A ;+H*45H*8,；I688? +*4>65H*8,；N>,O* /5+*,
基于栅格的分布式水文模型是目前水文研究的热点之一 # !-’- 年，QA;;R; 等［!］发表的“一个具有物理基 础数值模拟的水文响应模型的蓝图（QS’- 蓝图）”被认为是分布式水文模型研究的开始 # 过去受制于技术的 限制，许多流域信息、特征不能及时准确地获得，地理信息无法及时准确地获得和使用，使得分布式水文模型 发展缓慢 # 而数字高程模型（C2D）的出现为分布式水文模型的构建创造了条件，分布式水文模型取得了飞速 发展 # 比较具有代表性的分布式和半分布式流域水文模型有 TS2［.］，NUEVGE［W $］和 NUEDUC2X［0］模型等 # 随着 信息、遥感、计算机等技术的发展，基于 C2D 的分布式水文模型以其充分考虑水文要素和各种参数空间变化 的特点正成为水文模型的发展趋势，而新安江模型［1］作为一个半分布式的概念性模型，自研制以来，就得到 了广大学者的认同，在中国的洪水预报中也得到了广泛成功的应用，因此研究栅格型新安江模型具有一定的
!" + !#
$
,［ % *-
&
+./
’
&+
(
& +0.1］
（2）
式中："+———坡 面 水 流 的 水 深；#———时 间；%*-———栅 格 单 元 的 面 积；&+———栅 格 单 元 的 地 表 径 流 流 量； &+0.1———栅格单元的地表径流出流量；&+ቤተ መጻሕፍቲ ባይዱ/———上游栅格入流量 !
本文采用基于两步 3(-40"5(-6 算法［,,!,7］的二阶显式有限差分格式进行坡面与河道水流扩散波方程组
" 567
#
( ) $ 85 %57.9
: & ;&
’ <*
( 567 % 57.9
（:）
式中：" 567———植被冠层的累积截留量；$ 85———植被覆盖率；( 567———累积降雨量；’ <*———植被密度的校正因 子；%57.9———植被冠层的截留能力，即植被冠层的最大截留量 ! !!, 河道降水
量交换问题以及河道排水网络的影响 ! 即：如果当前计算单元的土壤含水量处于未蓄满状态，则上游栅格的 入流量首先补充当前栅格的土壤含水量，直至其蓄满为止；如果当前栅格有河道存在，则地表径流的出流量
将先按一定的比例汇入河道，然后再汇至下游栅格 !
!"#$%&#’’(’)#(’* 模型的坡面水流以及河道水流运动均利用一维扩散波方程组来描述，在进行栅格间汇 流演算时，坡面水流运动方程需在每个栅格单元上进行离散，其中的连续性方程为
95*)&%1)：/5+;? 8, H:; H:;8A7 8I H:; B*,’5,)*5,< 48?;6 5,? 8,;3?*4;,+*8,56 ?*II>+*8, J5K; 48?;6，H:; <A*?3L5+;? B*,’ 5,)*5,< 48?;6（H:; @A*?3B*,’5,)*5,< 48?;6）J5+ ?;K;68M;? # N:*+ 48?;6 A;<5A?;? H:; <A*? C2D 5+ H:; 984M>H5H*8,56 ;6;4;,H，5,? H:; <;,;A5H;? A>,8II 8I ;59: ;6;4;,H J5+ ?*K*?;? *,H8 H:A;; 984M8,;,H+，H:; +>AI59; A>,8II，*,H;AI68J，5,? <A8>,?J5H;A A>,8II，J:*9: J;A; A8>H;? >+*,< H:; 9;663H839;66 ?*II>+*8, J5K; I68J A8>H*,< 4;H:8? 5998A?*,< H8 H:; 984M>H5H*8,56 8A?;A 8I ;59: <A*? 9;66 *, H:; C2D# C>A*,< H:; 9569>65H*8, 8I H:; A>,8II <;,;A5H*8, 5,? A>,8II 98,9;,HA5H*8,，H:; J5H;A ;O9:5,<; L;HJ;;, <A*?+ 5,? H:; *,I6>;,9; 8I H:; A*K;A ?A5*,5<; ,;HJ8AP J;A; H5P;, *,H8 98,+*?;A5H*8, # /5+;? 8, H:; >,?;A67*,< +>AI59; 98,?*H*8,+，*,96>?*,< H:; H8M8<A5M:7，5,? H:; H7M;+ 8I +8*6 5,? K;<;H5H*8, *, H:; +H>?7 L5+*,，5 4;H:8? I8A M5A54;H;A ;+H*45H*8, 8I H:; 48?;6 J5+ MA8M8+;? 5,? K;A*I*;? # N:; 48?;6 J5+ 5MM6*;? H8 H:; N>,O* /5+*, *, G,:>* EA8K*,9; I8A I688? +*4>65H*8, 5,? MA8?>9;? I5K8A5L6; A;+>6H+ #