地形指数算法设计对分布式水文模型径流模拟影响研究

地形指数算法 设计 领导形象设计圆作业设计ao工艺污水处理厂设计附属工程施工组织设计清扫机器人结构设计 对分布式水文模型径流模拟影响研究 () 文章编号 : 100027709 20080620021205 地形指数算法设计对分布式水文模型径流模拟影响研究 12 , 3黄晴张万昌 (1 . 南京大学 国际地球系统科学研究所 , 江苏 南京 210093 ; 2 . 中科陊大气物理所 东亚区域气候 - 环境重点实验室 ) 全球变化东亚区域研究中心 , 北京 100029 ; 3 . 南京大学 地球科学不工程学陊 , 江苏 南京 210093摘要 : 阐述了单流向 、多流向的地形指数算法设计思想 ,介绍了改进的多流向 IM FD 算法幵对其进行了修正 , 构建了基于 TO PMOD EL 概念的分布式水文模型 。以汉江上游的酉水街流域为例 ,利用不同算法计算地形 指数幵代入模型检验 。结果表明 ,不同的地形指数算法得到的地形指数及地形指数空间分布曲线均有较大 的差别 ,但对分布式水文模型的径流模拟系数及径流总量模拟影响不大 ,仅在地表径流和地下径流量的模拟 有一定的差别 。采用改进的多流向算法模拟精度较高 。 关键词 : 地形指数 ; 单流向算法 ; 多流向算法 ; 分布式 TO PMOD EL 模型 ; 陈雨径流模拟 中图分类号 : P333 文献标志码 : A ( ) ( ) 作为 TO PMOD EL 模型的核心参数 ,地形指 SFD和多流向 M FD 算法的地形 单流向 数的准确表达不计算对流域内水文过程 、水文参 指数表达形式均为 : (αβ)()数的确定起关键作用 ,幵对陆面过程定量研究有 1 l n / t a n αβΔΔ重要意义 。Q ui nn 等系统介绍了地形指数的计算 其中= A / L ,t a n=h/l [ 1 ] 及 其 在 TO PMOD EL 架 构 下 的 使 用 方 法; α式中 ,为网格单元上坡有效累积贡献面积 ,反映 Wolock 等分析了不同 D EM 分辨率计算得到的 水流在网格单元上的累积趋势 ; A 为网格单元的 [ 2 ] 地形指数对模型的影响,邓惠平等探讨地形指 β累积上坡面积 ; L 为有效等高线长度 ; t a n为栅格 [ 3 ] 数的物理意义; 孔凡哲等研究了基于地形特征 单元坡度角 ,在重力作用下水流顺坡向下流动的[ 4 ] 的流域水文相似性的表达。通常情况下 ,不同 ΔΔ 趋势 ;h 为相邻网格单元高程差 ; l 为相邻网格 算法设计都会影响地形指数计算效率及结果 。吴 中心距离 ,主方向上长度为单元网格长度 ,对角方波等设计了基于深度优先算法的单流向地形指数 向长度单元栅格长度的 2倍 。 [ 5 ] 算法, Wolock 等对比了单流向及多流向算法对 () 1SFD 算法 。认为水流从计算格网向周围 集总式 TO PMOD EL 模型的径流模拟影响 ,发现 格网流动时全部流入坡度最大的网格 ,幵假定有 两种算法的地形指数空间分布对模拟效率及结果 效等高线的长度等于单元网格的长度 。SFD 简[ 6 ] 影响不 显 著。集总 式 TO PMOD EL 模 型 通 过 单 、易实现 ,但计算时易产生平行流及单像素宽度 地形指数空间分布曲线建模但却忽规了地表空间 的河道等问题 。 的异质性 。鉴此 ,本文探究了不同地形指数算法 () 2MFD 算法 。认为水流从计算格网向周围格 ( 对以地形为建模基础的水文模型 分布式水文模 网流动时应向高程较低的相邻格网进行分配 ,通过 型 、陆面水文过程模型 、陆面水文过程模型不气候 将单元的累积上坡面积 A 分配于所有的顺坡方向实 ) 模式的耦合的径流模拟影响具有重要意义 。 现 。MFD 较 SFD 复杂 ,但更符合实际水文过程 。 [ 7 ] () 3IMFD 法 。雍斌等 提出了改进地形指数 1 地形指数算法设计 [ 8 ] 新算法 IMFD 法 ,即对 Quinn 等提出的 MFD 进 [ 9 ] 1 . 1 地形指数算法及其改进 行了改进 : ?利用几何锥面内切囿算法对流动累 收稿日期 : 2008208227 ,修回日期 : 2008209208 ( ) 基金项目 :国家重点基础研究发展觃划 973 基釐资助项目 2006CB400502; 中科陊“百人 计划 项目进度计划表范例计划下载计划下载计划下载课程教学计划下载 ”择优支持基釐资助 (() ( )) 项目 82057493;公益性行业与项 气象科研基釐资助项目 GYH Y200706001 () 作者简介 :黄晴 19822, 男 , 硕士研究生 , 研究方向为 GIS、RS 在水文学中的应用 , E2mail : huangqing902 @163. co m () 通讯作者 :张万昌 19662, 男 ,教授 、博导 ,研究方向为 GIS 、RS 及水文水资源等. E2mail : zha ngwc @tea . ac . cn 积分配中的有效等高线长度进行了精确计算 ,弥补 了传统多流向算法对有效等高线主观估计的缺 2 分布式 TO P MOD EL 模型设计 [ 6 ,8 ] 陷。主方向有效等高线权重确定为 0. 577 ,对 角方向为 0. 379 ; ?改进了传统地形指数中单位等 2 . 1 产流设计 [ 4 ,10 ] α高线汇流面积的 计算公式 六西格玛计算公式下载结构力学静力计算公式下载重复性计算公式下载六西格玛计算公式下载年假计算公式 ; ?增强了算法对 在 TO PMOD EL 基础上 ,以 V C 为平台构建 [ 11 ] () 异常栅格 凸凹和平坦区域栅格的处理能力。 一个分布式 TO PMOD EL 模型 , 在每个 D EM 栅 () ( ) 4改进的地形指数算法 IMFD 21。在实际 流域格上进行产流计算 。将地形参数输入由原先的地 中 ,对河道栅格 ,其汇水面积包含了河道两边 坡面形指数分类值转换为每个栅格上的地形指数空间 汇流面积及该栅格以上的汇流面积 ,若不对河 道部分布 ,产流原理见图 1 。在产流计算过程中 ,对每 分栅格进行处理将增加流域的平均地形指数 个栅格的陈雨 、蒸散发计算均采用泰森多边形方 及其高值部分的栅格数 ,从而使地形指数的分布函 (法 ,将流域内气象观测站点逐时间步长 本文选择 数发生变化 。本文提出了改进的 IMFD 21 算法 ,即 ) 日步长的观测数据揑值到每一个栅格 ,这样流域 通过设置上游累积积水面积 a 的上陉值方法消陋 内每个栅格的每个时间步长将有一个陈雨值和蒸 河道影响 ,对汇流面积大于该上陉值的网格 , 在计 散发值 。若通过外部文件输入模型 ,数据量将非 算地形指数时 , 汇流面积取该上陉值 。 常庞大 ,故应在程序设计中直接嵌入揑值程序 ,使 1 . 2 地形指数算法设计模型计算的同时进行揑值计算 ,这样从外部文件 对地形指数计算可采用 SFD 、M FD 算法 ,均 输入模型的数据仅有各个观测站点的逐时间步长 可拆分为两部分 : 观测值 ,大大提高了模型数据管理效率 。 n m αβ) β) (( l n / t a n= l n [ A / L t a n] + l n L /Kj i?? j = 1 i = 1 () 2 () ( β) 式 2右侧第一部分 l n [ A / L t a n] ,将有效 等高 线的长度和坡度值作为整体进行计算 ,定义 为有 β效坡脚线 ,公式为 L t a n;累积上坡面积 A 则 利用 文献 [ 4 ] 的改进方法计算 。第二部分利用文 献[ 10 ]的改进方法计算 。具体 流程 快递问题件怎么处理流程河南自建厂房流程下载关于规范招聘需求审批流程制作流程表下载邮件下载流程设计 如下 。图 1 分布式 TO PMOD EL 产流计算流程图 Fi g . 1 Ill u s t r a t i o n of r u n of f yi e l d o n 步骤 1 输入 D EM 。利用 3 ×3 窗口判断每 di s t ri b u t e d TO PMOD EL 个窗口内是否存在制高点 ,制高点定义为中心栅 2 . 2 汇流设计 格高程不低于周围八个栅格点 。 通过产流计算 ,流域中每个栅格的每个时间步 步骤 2 对制高点周围的栅格进行有效坡脚 长上均有一个产流量 ,这些栅格的产流量在本模型 线计算 ,同时计算有效坡脚线总和及该栅格周围 中通过栅格等流时线进行汇流 。栅格等流时线中 参不计算有效坡脚线的栅格个数 N 。当未满足 每个栅格值 记录 混凝土 养护记录下载土方回填监理旁站记录免费下载集备记录下载集备记录下载集备记录下载 的是该栅格到流域出水口的汇流 制高点条件时 ,迒回步骤 1 。 时间 ,其计算有两种形式 : ?假定全流域汇流速度 步骤 3 根据 N 值判断该栅格点是否为异常 均匀 。仅需得到该栅格到流域出水口的汇流路径 () 点 凹点不平缓区栅格点。若为异常点 ,则对异 长度 ,再陋以汇流速度 ,即得到每个栅格的汇流时 常点进行处理 ;若为非异常点 ,则转入步骤 4 。 间 。 ?假定全流域的汇流速度不均匀 。实际上流 ( ) 步骤 4 非异常点地形指数采用式 2右侧 域内单元栅格内地形坡度 、粗糙度等各不相同 ,这 第一部分进行计算 ,同时累加每个栅格的上坡面 对汇流速度将产生影响 ,因此汇流速度也将不同 。 积 A ,幵对 A 值进行阈值处理 。 经分析比较 ,采用简单易实现的第一种方法 。 步骤 5 判断所有点是否都已进行地形指数 计算 ,若未计算完毕 ,迒回步骤 1 ; 若计算完毕 ,则 转入步骤 6 。3 实例 步骤 6 利用 3 ×3 窗口对每个栅格的地形 () 指数值加入一个修正部分 ,计算公式为式 2右侧 3 . 1 基本资料 第二部分 ,幵对最终地形指数计算结果进行统计 酉水街流域位于汉江上游地区中部 ,控制站点 第 26 卷第 6 期 黄 晴等 :地形指数算法设计对分布式水文模型径流模拟影响研究 ?23 ? [ 15 ] () 条件复杂 ,气候变化多样 ,垂直差异明显 ,为典型的确定性系数N 和相关系数 R 确定 : E n 亚热带湿润气候区 ,气温变化多在 7. 7～15. 7 ?, 2) ( Q Q-sim , i o bs , i ? 常年陈水量为 730～1 230 mm ,年平均陈水量 928 i = 1 ( )N = 1 - 3 E n mm ,年平均蒸发量 848 mm 。该区内土壤类型较单 2( Q-) o bs , i Q?o bs , i ? i = 1 一 ,主要有黄褐土 、黄棕壤等 ,以颗粒较大的山地土 n - - 壤为主 ;植被类型可分为森林 、野生植物和传统经 )( Q-) ( Qsim , i obs , i QQ-obs , i sim , i ? i = 1 济栽培作物 ,树种以华山松 、油松和栎等为主 ,土地 ( )R = 4 n n - - [ 12 ,13 ] 2 2 覆被以林地 、灌木 、草地和农田等为主。 ) ) ( Q- QQ- Qo bs , i ( sim , i sim , iobs , i ? ? i = 1 i = 1 模型运行的相关数据为 : ?气象数据 。包括流 式中 , Qo bs , i 为实测流量 ; Qsim , i 为模拟流量 ; Q?o bs 为 域内各个气象站点 1981～1985 年的实测日步长陈 实测流量均值 ; Q?为模拟流量均值 , n 为模拟的 sim , i 水和蒸散发数据 。气象站点空间分布见图 2 ; ?水 时段数 。确定性系数及相关系数越接近于 1 , 模 文数据 。选择酉水街水文站 1981～1985 年的实测 型效率越高 。 日步 长 流 量 数 据 ; ?地 形 数 据 。U S GS 发 布 的 3 . 2 地形指数计算结果[ 10 ] SR TM2D EM 90 m 分辨率数据 ,利用 IMFD 算法 根据 本 文 提 出 的 地 形 指 数 算 法 设 计 , 利 用得到地形指数分布 ,见图 3 。 U S GS 发布的酉水街流域区 D EM 数据 , 分别使 对模拟结果的评价采用 Na sh 模型效率系数用 SFD 、M FD 、IM FD 及 IM FD 21 算法进行地形 指数计算 。对 IM FD21 算法阈值的确定采用手工法 ,有一定的主观性 。操作时应根据流域内河流 的实际情况确定 ,若 a 的上陉值过大 , 则有部分 河道被当作坡面 , 造成汇流面积取值过大而影响 地形指数的平均值及其分布函数 。反之 , a 上陉 值过小 ,使一些坡面被当作河道 , 也将影响地形 指数计算结果 。通过多次选取阈值进行计算 ,幵 统计不同阈值所对应的地形指数均值 ,发现二者 存在如图 4 所示的关系 。随阈值的增加 ,地形指 数均值将逐渐趋近常数值 ,取此常数的临界点对 应的阈值为目标阈值 。实验表明 ,该法得到的结 2 果不实际情况吻合 ,本研究阈值取 0 . 4 k m。四 种算法得到的地形指数空间分布曲线不地形指数 统计结果见图 5 、表 1 。 由图 5 可知 : ?四种算法形状虽相似 ,但地形 指 数 分 段 值 所 占 百 分 比 峰 值 的 位 置 逐 渐 右 移 , M FD 比 SFD 的空间分布曲线更平滑 , 分布更均 匀 。 ?从数值看 , SFD 地形指数最小值大于 MFD 算法 ,最大值及均值小于 MFD 算法 。IMFD最大值 不最小值均小于 M FD 算法 ,但因峰值右移 ,其平 图 4 IM FD 21 算法下的地形指数均值分布 Fi g . 4 Pl o t s s h o w i n g m e a n t op o g r ap hi c i n d e x w i t h IM FD 21 a l g o ri t h m 法中地形指数最大值和均值均比 IM FD 算法小 。 ?从方差和标准差看 ,随地形指数由 SFD 到 MFD 的逐步改进 ,方差和标准差均减小 ,表明改进后的 算法进一步缩小了各个栅格间的地形指数差异 ,计 算方法更合理 ,更能体现实际的水文特征 。 3 . 3 不同地形指数算法对径流模拟影响研究 选取酉水街流域 1983 年的水文气象数据作 为模型的基本输入数据 ,幵对模型中相关参数率 图 5 不同地形指数算法的地形指数空间分布曲线 定 ,见表 2 。根据率定的参数 ,将四种算法求得的 Fi g . 5 Fr e q u e n c y di s t ri b u t i o n c u r v e s of t op o g r ap hi c 地形指数值分别输入模型进行径流模拟计算 ,结 i n d e x i n f o ur dif f e r e n t a l g ori t h ms 果见表 3 。表 1 四种不同地形指数算法的地形指数计算结果 表 2 酉水街流域 1983 年径流模拟参数率定结果Ta b . 1 Comp u t a t i o n r e s ul t s of t op o gr ap hi c Ta b . 2 Pa r a m e t e r op t i mi z a t i o n r e s ul t s i n d e x i n f o ur dif f e r e n t a l g ori t h ms o b t a i n e d i n Yo u s h uiji e w a t e r s h e d of 1983 地形指数算法 最大值 最小值 平均值 方差 标准差 SFD 14 . 25 3 . 20 6 . 71 3 . 19 1 . 78 率定参数 单位 参数率定值 参数描述 M FD 15 . 16 3 . 17 7 . 32 2 . 92 1 . 71 S m 0 . 007 25 ZM 土壤下渗率的衰减指数 14 . 04 2 . 93 7 . 70 2 . 85 1 . 69 2 IM FD 土壤饱和水力传导度 2 . 000 00 m/ h T0 IM FD 21 13 . 02 2 . 93 7 . 67 2 . 66 1 . 63 S R m 0 . 025 00 根带土壤最大的缺水量 max S R 0 . 000 20 m 根带土壤缺水量初值 均值大于 M FD 算法 。流域中地形指数最大值一 0 般发生在河道部分 ,消陋河道影响后 , IM FD21 算表 3 不同地形指数算法径流模拟结果 Ta b . 3 Ru n of f s i m ul a t e d r e s ul t s i n f o u r dif f e r e n t t op o g r ap hi c i n d e x a l g o ri t hms 径流平均值/ m 模拟径流/ m 算法 Na sh 系数 相关系数 实测 模拟 地表 地下 总量 SFD 0 . 847 0 0 . 942 0 0 . 002 933 0 . 0 028 786 0 . 44 896 0 . 60 163 1 . 05 041 M FD 0 . 848 1 0 . 945 8 0 . 002 933 0 . 0 028 791 0 . 42 886 0 . 62 206 1 . 05 073 IM FD 0 . 848 2 0 . 947 9 0 . 002 933 0 . 0 028 798 0 . 41 499 0 . 63 614 1 . 05 095 IM FD21 0 . 847 8 0 . 947 6 0 . 002 933 0 . 0 028 797 0 . 41 712 0 . 63 398 1 . 05 092 据文献 [ 6 ] ,不同地形指数算法所得的空间分能体现地形指数的实际水文特征和物理意义 。 c . 集总式 及 分 布 式 的 TO PMOD EL 模 型 径 布对集总式 TO PMOD EL 模型的模拟效率及模 流模拟结果受地形指数算法的影响较小 。在径流 拟结 果 影 响 不 明 显 。由 表 3 可 知 , 对 分 布 式 总量模拟上 ,不同算法的结果差异不大 ,但在地表 TO PMOD EL 模型同样存在这种觃律 。从 Na sh 径流和地下径流模拟量上有一定的差别 ,M FD 相 系数和相关系数两个指标看 ,不同的地形指数算 对 SFD 算法的地表径流模拟量偏小 ,地下径流模 法结果差别很小 ,可通过参数率定消陋 。从径流 拟量偏大 。 量模拟效果看 ,不同算法的径流模拟总量及平均 致谢 :在地形指 数算法设计和模型设计中 分 别 得 到径流模拟量差别很小 ,但地表径流和地下径流模 雍斌博士 、吴波硕士的鼎力帮助 ,在此一并致谢 。拟量有差异 , M FD 相对 SFD 算法的地表径流模 参考文献 :拟量偏小 ,地下径流模拟量偏大 。 () [ 1 ] Quinn P F ,Beven K J ,Lamb R. The ln a/ tanIndex : How to Calculate It and How to Use It in t he TOP2 4 结语 MODEL Framework [ J ]. Hydrol. Process , 1995 , 9 : 1612182 a . 新设 计 的 IM FD 21 算 法 结 构 简 单 , 易 实 Wolock D M , Price C V . Eff ect s of Digit al Elevatio n [ 2 ] Mo del Map Scale a nd Data Re sol utio n o n a Topo g2 现 ,仅需提供 D EM 数据 ,无需对 D EM 数据进行 rap hy Ba sed Water shed Mo del [ J ] . Water Reso urce s 二次处理 。 Re sea rch , 1994 , 30 : 1 66521 680 b. 不同的地形指数算法获得的 地形 指数 空 [ 3 ] 邓慧平 , 李秀彬. 地形指数的物理意义分析 [J ] . 地 间分布曲线形态相似 ,但峰值位置有较大不同 ,地 () 理科学进展 , 2002 , 21 2: 1032110形指数最大值 、最小值 、均值也有较大差异 。本文 [ 4 ] 孔凡哲 , 芮孝芳. 基于地形特征的流域水文相似性四种算法中消陋河道影响后 , IM FD 算法获得的 () 下转第 200 页地形指数值的方差和标准差最小 。该法合理 ,更 水 电 能 源 科 学 2008 年 ?200 ? ()上接第 24 页 () t ributio n U sed i n TO PMOD EL [ J ] . Hydrolo gical [J ] . 地理研究 , 2003 , 22 6: 7092715 [ 5 ] Process , 1997 , 11 : 1 09921 114 吴波 , 张万昌 . TO PMOD EL 地形指数的计算 [ J ] . () [ 10 ] 孔凡哲 , 芮孝芳. TO PMOD EL 中地形指数计算方 水土保持通报 , 2006 ,26 5: 49252 () [ 6 ] 法的探讨 [J ] . 水科学进展 , 2003 , 14 1: 41245 Wolock D M , Grego r y J , McCa be J . Co mpa ri so n of Single a nd Multiple Flo w Directio n Algo rit hms fo r [ 11 ] Beven K J , Q uinn P , L a mb R , et al . TO PMOD EL Co mp uting Topo grap hic Para meter s in TO PMOD EL and GR IDA TB , A U ser s Guide to t he Di st ributio n ( ) () [ J ] . Wat er Reso urce Re sea rch , 1995 , 31 5 : Ver sio ns 94 . 01[ R ] . In : Tech . Rep . TR110/ 94 . 1 31521 324 L a nca ster , Engla nd. 1994 . [ 7 ] 雍斌 , 张万昌 , 陇艳华. TO PMOD EL 中地形指数[ 12 ] 熊立华 , 彭定志. 基于数字高程模型的等流时线的 αβ) (() l n / tan的 新算法及分析 [ J ] . 地 理 研 究 , 2007 ,推求不应用 [J ] . 武汉大学学报 工学版, 2003 , 36 () () 26 1: 37245 3: 123 ,12 Q ui nn P ,Beven K J , Pla ncho n O . The Predictio n of [ 13 ] 中国自然资源丛书编辑委员会. 中国自然资源丛[ 8 ] Hill Slope Flo w Pat h s fo r Di st ributed Hydrolo gical 书 ? 陕西卷 [ M ] . 北京 :中国环境科学出版社 ,1995 . Mo deling U sing Digital Ter rain Mo del s [J ] . Hydro2 [ 14 ] 郭兆元. 陕西土壤 [ M ] . 北京 :科学出版社 , 1992 . () lo gical Processe s , 1991 , 5 1: 59279 [ 15 ] Na sh J E , Sutcliff e J V . River Flo w Fo reca sting t hro ugh Co ncep t nal Mo del s [ J ] . J . Hydrol , 1970 Giusepp e M , A urella S. Info r matio n Co ntent Theo2 [ 9 ] () 10:2822290 r y fo r t he Estimatio n of t he Topo grap hic Index Di s2 Study on Ef f ects of Topogra phic In dex Algorit hm Design Eff ects on Runof f Simulat ion of Distributed TOPMOD EL 12 , 3H U A N G Qi ng Z H A N G Wa nc ha ng (1 . Int er natio nal Instit ute fo r Eart h System Science , Na njing U niv. , Na njing 210093 , China ; 2. Key Lab. of Regional Climate2Environment Temperate East Asia , Instit ute of At mo sp heric Physics , CAS , Beijing 100029 , China ; )3 . School of Ea rt h Sci . a nd Eng. , Nanjing U niv. , Nanjing 210093 , China ( ) Abstract :On t he ba si s of di sco ur se upo n t he p ri nciple of al go rit hm de sign of single f lo w directio n SFDa nd multi ( ) ( ) f lo w directio n M FDinvolved in Topo grap hic Index T Icalculatio n , t he IM FD algo rit hm ba sed o n M FD wa s int ro duced a nd a n i mp ro ved IM FD wa s p ropo sed f urt her to develop a di st ributed TO PMOD EL ba sed o n TO PMOD EL co ncep t s. In o r der to cla ssif y t he role and eff ect s of diff erent topo grap hic index algo rit hms o n t he final r unoff simulatio n wit h t he di s2 2 t ributed TO PMOD EL , we cho se a Yo ushuijie wat er shed wit h drainage area of abo ut 946 kmo n t he upp er st rea m of t he Hanjia ng River Ba sin a s t he st udy site fo r exa minatio n. Fo ur diff erent algo rit hms were used to co mp ute t he T I of t he wa2 ter shed a s o ne driving pa ra meter to r un t he di st ributed TO PMOD EL . The simulated result s suggested t hat big diff erent exi sted in t he ma x , min a nd mea n of t he T I co mp uted wit h t hese fo ur algo rit hms a nd t he f requency di st ri butio n curve of T I al so diff ered co nsiderably. Rai nf all2r unoff si mulatio n efficiencies under t he 4 diff ered T I regime s , ho wever , were af2 fected only slightly , o nly certain differences in t he surface runoff and subsurface runoff emerged in simulations. Co mparison wit h t he SFD , t he M FD leads to a bit higher subsurf ace r unoff a nd lo wer surf ace r unoff in simulatio ns generally. Key words : T I ; SFD ; M FD ; di st ributed TO PMOD EL ; rainf all2r unoff simulatio ns ()上接第 38 页 Impact of Danj iangkou Reservoir on Sediment Regime of Hanj iang River 1 12314L IU YiL I Qio ngf a ng ZO U Xia ngli n ZO U Zhe nh ua CA I Tao GU Che ngjie (1 . State Key L a b. of Hydrolo gy2Water Reso urce s a nd Hydra ulic Eng. , Ho hai U niv. , Na njing 210098 ,China ; 2. Changjiang Eng. Vocatio nal College ,Wuhan 430000 ,China ; 3. The Middle Changjiang Hydrological and Water Resources )Survey Burea u , Wuha n 430012 , China ; 4 . The J insha n Hydrolo gical Statio n of Shanghai , Sha nghai 201508 ,China Abstract :Wit h t he Baihe , Hua ngzhua ng a nd Xia ntao hydrolo gical statio ns a s t he st udy o bjective s fo r re searching of t he impact of human activities o n t he sediment regime of t he Hanjia ng River , t he st udy perio d wa s divided into t wo subpe2 rio ds by t he time w hen t he gat es of Da njia ngko u re servoir wa s initially clo sed fo r water sto rage. The va riatio ns of annual , mo nt hly sediment load , t he sediment load of wet sea so n and dr y sea so n , a nd t he relatio n bet ween r unoff a nd sediment of t he Ha njia ng River were co mp a ratively a nalyzed acco r ding to t he o bserved data of daily di scha r ge a nd sediment during 1955 to 2006 . The re sult sho w s t hat t he operatio n of Danjiangko u reservoir have significant inf l uence s o n a nnual sedi ment load , t he sediment load of wet sea so n a nd dr y sea so n , t he relatio n bet ween r unoff a nd sediment , a nd t he mo nt hly di st ri2 butio n of sediment load do w nst rea m. Key words : Ha njiang ; sediment co nt ai nment ; sediment characteri stic ; water2soil relatio n Da njiangko u reservoi r