GISECOlecture5ext02.pptVIP

地理信息系统在资源环境研究中的应用 宋现锋 中国科学院研究生院资源与环境学院 Phone: 010-8825-6462 Email: xfsong@ VII 河流主干道作用过程（Main Channel Processes） 流域内的水流动类型分为：地表漫流（overland）与沟渠疏导（channelized）。二者的主要区别：沟渠疏导考虑了水道蓄水（water storage）及其对流量（flow rate）的影响。 SWAT建模的河流主干道作用过程包括：水的移动，沉积物的沉积过程，河网内部的氮循环，以及杀虫剂等的相变。其他过程还包括：河床下切侵蚀/侧切侵蚀（downcutting/widening）引起的河道空间尺寸随时间的动态变化。 1 Water Routing 明渠水流（ open channel flow）是指具有自由水面的水道，比如：河流的水流、部分充水的管道水流。SWAT使用曼宁方程定义流量与流速，使用可变槽蓄法或者马斯京根法 （Muskingum）模拟河网的水路演算（Water Routing）。二者均是运动波模式 （Kinematic Wave Model，Chow et al. 1988）的变种。 Open channel flow is defined as channel flow with a free surface, such as flow in a river or partially full pipe. SWAT uses Manning’s equation to define the rate and velocity of flow. Water is routed through the channel network using the variable storage routing method or the Muskingum river routing method. Both the variable storage and Muskingum routing methods are variations of the kinematic wave model. A detailed discussion of the kinematic wave flood routing model can be found in Chow et al. (1988). 1.1 河道的特征参数 SWAT中，假设主干河道、河段为梯形横断面形状。 用户输入参数：河道的宽度/深度（蓄水到达河堤顶部时），河道长度，河道坡度，河道曼宁系数 - n值。 ‘ Wbtm = Wbnkfull - 2·zch·depthbnkfull Wbtm 河道的底宽（米） Wbnkfull 河道蓄满水后的顶宽（米） zch 河道边坡坡度的倒数 depthbnkfull 河道蓄满水后的水深（米） SWAT中，假设zch=2，上式的河道底宽Wbtm可能会小于或者等于零！这种情况下，SWAT调整上式： Wbtm = 0.5·Wbnkfull zch = (Wbnkfull-Wbtm)/(2·depthbnkfull) 如果河道的水深depth，则河道的水平面宽度为： W = Wbtm + 2·zch·depth Wbtm 河道的底宽（米） W 河道蓄水的水平面宽度（米） zch 河道边坡坡度的倒数 depth 河道蓄水的水深（米） 河道水流的横断面积： Ach = (Wbtm + zch·depth ) ·depth Ach横断面积（米2） Wbtm 河道底宽（米） zch 河道边坡坡度的倒数 depth 河道水深（米） 河道的湿润周界： Pch = Wbtm + 2·depth·(1+zch2) 1/2 河道内的水流，在横断面上形成的湿润外缘的总长度 水力半径： Rch = Ach / Pch 河道水流的横断面面积与湿周的比值 河道的水量： Vch = 1000·Lch·Ach Vch 河道内的水量（m3） Lch 河道的长度（Km） Ach 河道水流的断面积（m2） 河道水量超过河道容量时，过剩的水将在漫滩（阶地）漫延。 SWAT模型的假设条件： 漫滩的底宽为河道顶宽的5倍。 漫滩的边坡坡度为0.25。 漫滩水流的深度、横断面积、湿周长度，计算如下： 1.2 流量与流速 均匀流体的曼宁方程，按照时间步长计算河段内的流量与流速： qch 河道水流的流量（m3/s） vch 河道水流的流速（m/s） Ach 河道水流的断面积（m2） Rch 河道的水力半径系数（m） slpch 河道的边坡坡度（m/m） n 河道的曼宁系数 SWAT模型的水流演算