R在水文建模中的应用初探课件

,Clique para editar o estilo do ttulo mestre,Clique para editar os estilos do texto mestre,Segundo nvel,Terceiro nvel,Quarto nvel,Quinto nvel,R,在水文建模中的应用初探,东北林业大学 生态学,06,级硕士 王化儒,2008.11.14,汇报的主要内容,水文过程与模拟模型,R,在水文模拟中的应用,实例演示,水文过程与模拟模型,水文过程,坡面水文过程,水文模型,y = f (,| x),y,径流,x,降水,参数集,水文模型,f,（）,降水,x,径流,y,Flow,Time,时间,水文模型的种类,集总式水文模型,：不包含任何空间信息，参数的物理意义不明显，容易估计,半分布式水文模型,：通过水文响应单元隐式包含空间信息，参数具有明显的物理意义，容易获取,分布式水文模型,：显式地包含空间信息，参数具有明显的物理意义，难以获取,个人观点：,当前的分布式水文模型严格地说都是半分布式水文模型，真正的分布式水文模型是不存在的,TOPMODEL,模型,半分布式水文模型代表,TOPMODEL,模型是英国,Lancaster,大学,Beven,和,Kirkb,于,1979,年,提出的一个基于地形的半分布式水文模型,该模型的理论基础是变动源面积，核心是地形指数。经历了三十年的发展，,TOPMODEL,已然成为了半分布式水文模型的代表，在国内外许多地区得到了很好的应用，受到了广泛的赞誉。,http:/,www.es.lancs.ac.uk/hfdg/research/hfdg_res_cat.htm,Prof. Keith,Beven,径流形成机制,P,整个流域内超渗产流,部分区域超渗产流,蓄满产流,P,P,q,r,q,s,q,o,P,P,P,q,o,f,P,P,P,q,o,f,f,TOPMODEL,中采用的是产流机制是蓄满产流，适用于湿润和半湿润地区，不适用于干旱地区,变动源面积,-,TOPMODEL,中的产流机制,变动源面积的含义是坡面流只在整个流域的在降雨事件中地下水位上升至地表的饱和区域产生，此时降落到这个饱和区域上的降水将形成直接径流。在整个降水过程中，源面积是不断变化的，流域源面积的位置受流域地形和土壤水力特性的影响。饱和区域土壤水力传导性差、坡面平缓，一般位于河道附近，随着下渗的持续，饱和区域向河道两边的坡面延伸，在一定意义上，变动源面积可看作是河道系统的延伸。在变动源面积的概念中，饱和区域上形成的径流有两种方式：饱和坡面流和壤中流。,P,P,q,r,q,s,q,o,流域水文表面分析,首先要识别那些周围高、中间低的凹陷点，因为凹陷点的水流方向不能用代表流向的特征码来表示。,而后进行流向分析，有单流向算法和多流向算法，前者包括,D8,算法、,Rho8,算法等。这里采用,D,算法。,凹陷点,D,算法,80,74,63,69,67,56,60,52,48,80,74,63,69,67,56,60,52,48,Tarboton, D. G., (1997), A New Method for the Determination of Flow Directions and Contributing Areas in Grid Digital Elevation Models, Water Resources Research, 33(2): 309-319.) (,http:/,www.engineering.usu.edu/cee/faculty/dtarb/dinf.pdf,),32,16,8,64,4,128,1,2,地形指数,-,TOPMODEL,中的水文响应单元,1,式中：,A:,特定上坡面积；,：坡度角。,地形是影响径流响应的重要因子，它决定了重力作用下流域中水分运动的趋势。,TOPMODEL,提出了地形指数作为水文响应单元来反映径流运动的分布规律，不论空间位置的水文单元，只要有着相同地形指数，其水文学特性就相似。地形指数体现了土壤湿度、地表饱和度的空间分布和径流生成的过程，它是,TOPMODEL,的核心，其计算公式如下。,地形数据易获取,地形图,免费的,DEM,：如,SRTM,（,90m,），,GTOPO30,（,1 km,）,http:/,srtm.usgs.gov/srtmimagegallery/index.html,http:/,srtm.csi.cgiar.org/SELECTION/inputCoord.asp,Santa Barbara, California,SRTM,数据,http:/www1.gsi.go.jp/geowww/globalmap-gsi/gtopo30/gtopo30.html,泾河流域,DEM,http:/, package in R,help(package,=,RHydro,),topmodel,package in R,help(package,=,topmodel,),Topographical analysis using,RHydro,dem,=,read.table(“dem.txt”,skip,=6),#,读入,DEM,数据,dem,=,as.matrix(dem,),demdem,=-9999=NA,#DEM,数据预处理,ARCGIS,导出的,GRID,格式的,DEM,dem.filled,-,sinkfind(dem,cellsize,=25, degree=0.1),#,凹陷点填充,topidx,-,atb(dem.filled,cellsize,=25),#,计算地形指数,topidx$atb,和上源面积,topidx$area,outlet(topidx$area,c(28,8),2),#,定义出水口,catchment -,subcatch(dem.filled, c(29,8),#,进行流域分割,flowlength,-,flowlength(dem.filled, c(29,8),#,计算河网长度,rivers -,river(dem.filled,topidx$atb,topidx$area,cellsize,=25,thatb,=12.35,tharea,=10000),#,划分河流网络,catchmentcatchment,= 0 - NA,rivers-rivers*catchment,#,提取流域内的河网,DEM,地形指数,上源面积,分割的流域,流域的河网,流域的地形指数,data(Huagrahuma,),#,没有观测流量数据情况下的模拟,Qsim,- topmodel(parameters,topidx,delay,rain,ET0),#,这个命令会返回,Qobs,qo,qs,和,S:,topmodel(parameters,topidx,delay,rain,ET0,verbose = TRUE),plot(Qobs,);,points(Qsim,col,=red, type=l),# Monte,carlo,模拟生成一个随机的样本参数集,:,runs-10;,qs0 -,runif(runs,)*4e-5;,lnTe,-,runif(runs,)*3-2;,m-,runif(runs,)*0.2,Sr0 -,runif(runs,)*0.02;,Srmax,-,runif(runs,)*2;,td -,runif(runs,)*3-3,vch,- 1000;vr - 100+runif(runs)*2400;,k0 -,runif(runs,)*0.01,CD -,runif(runs,)*5;dt - 0.25,parameters-cbind(qs0,lnTe,m,Sr0,Srmax,td,vch,vr,k0,CD,dt),result-topmodel(parameters,topidx,delay,rain,ET0,Qobs =,Qobs,),Hydrological simulation using,topmodel,模型模拟及效果评价,模拟效果采用,Nash-Sutcliffe,效率进行评价，它是评价水文,模型的常用指标。,E,表示,Nash-Sutcliffe,效率，,Qo,是,t,时刻观测的径流，,Qm,是,t,时刻模拟的径流。,E,的范围在,-,到,1,之间；当,E,为,1,时，模拟结果与实测数据能够完美地匹配；当,E,越接近于,1,，模拟效果越好。,Nash-Sutcliffe,效率为,0.828,讨论：,模型模拟的不确定性,数据与模型的融合，数据同化（,data assimilation,）,?,R,为多学科交叉提供一种理念和实现途径,不确定性分析,Generalised,Sensitivity Analysis (GSA),HSY Generalised Sensitivity Analysis,Generalised,Likelihood Uncertainty Estimation (GLUE),Monte Carlo Markov Chain methods,http:/www.uncertain-future.org.uk/RSoftware.htm,气候数据：,ACACIA Regional Climate Data Access System,Climate Research Unit (CRU) Global Climatology Data Set (10 minute Resolution),National Center for Atmospheric Research (NCAR),IPCC Data Distribution Centre Climate Model Data Portal,NASA Global Change Master Data Set Directory,National Climactic Data Center (NCDC) Archives,NCAR Climate and Global Dynamics Division,NOAA Earth System Research Laboratory /Physical Sciences Division (ESRL/PSD),Union of Concerned Scientists: Global Environment IPCC Page,World Climate Research Program/ Coupled Model,Intercomparison,Project (WCRP / CMIP3),R,中水文气候数据处理方法,http:/cran.r-project.org/web/views/Spatial.html,Acquire and process GRIB format GCM climate data files,Locate, read and analyze NetCDF format climate data files,实例演示,利用,R,中的,TOPMODEL,包和,RHydro,包在流域水文表面分析的基础上流域水文过程进行模拟，设置是否进行填洼处理来考察其对地形指数等水文响应单元的影响，并展示水分地形分析和水文模拟的结果。,谢谢！,