浅层地下水安全开采动态控制水位确定.pdf
第 3 2卷第1 0期 2 0 1 4年1 0月 水 电 能 源 科 学 W a t e r R e s o u r c e s a n d P o w e r V o l . 3 2N o . 1 0 O c t . 2 0 1 4 文章编号1 0 0 0-7 7 0 9( 2 0 1 4)1 0-0 0 6 3-0 4 浅层地下水安全开采动态控制水位确定 刘 波,刘 革,束龙仓 ( 河海大学 水文水资源学院,江苏 南京2 1 0 0 9 8) 摘要为深入研究地下水安全开采问题, 提出了地下水安全开采动态控制水位概念, 具体包括控制水位上限、 控制水位下限、 上限预警水位和下限预警水位四个指标, 各指标值随时间变化。详细介绍了水均衡法确定动 态控制水位的基本原理和方法步骤, 并以阜阳市为例, 选取代表性观测站( 周棚站) 进行计算, 得到研究区年内 浅层地下水安全开采动态控制水位过程。结果表明, 年内控制水位上限、 控制水位下限、 上限预警水位和下限 预警水位均随时间动态变化, 为地下水实时动态监督管理提供了可靠依据。 关键词浅层地下水;水均衡法;动态控制水位;预警水位 中图分类号TV 2 1 1. 1+2 文献标志码A 收稿日期2 0 1 4-0 6-2 3, 修回日期2 0 1 4 -0 8-1 0 基金项目国家自然科学基金青年基金项目(4 1 3 0 1 0 1 7) ; 国家自然科学基金项目( 4 1 1 7 2 2 0 3,4 1 2 0 1 0 2 9) 作者简介刘波(1 9 8 0-) , 女, 讲师, 研究方向为地下水资源评价及管理, E-m a i lB o b o l 1 3 7 0 5@1 6 3. c o m 1 引言 地下水是重要的供水水源, 随着社会经济的 快速发展, 很多地区地下水开采量不断增长, 开采 深度和范围不断扩大, 导致地下水位下降, 产生降 落漏斗、 机井掉泵和地面沉降等危害生产、 生活的 环境地质问题。地下水位动态是判断地下水是否 出现问题的重要指标, 因此有关地下水控制水位 的研究非常必要。目前, 通过分析地下水位埋深 与作物生长、 土壤含盐量、 蒸发量等关系, 已提出 了一系列地下水位控制指标, 如临界深度[ 1,2]、 合 理地下水生态水位[ 3]、 地下水控制性关键水位[4]、 适宜水位和危急水位[ 5]等概念, 并分别采用水均 衡法、 数理统计法、 回归分析法、 人工神经网络法 和数值模拟法等确定了地下水位控制指标值。但 已有研究中各地下水位控制指标值均是一个固定 值, 而地下水位受补给、 迳流、 排泄条件的影响, 是 随时间动态变化的。因此, 本文提出了地下水安 全开采动态控制水位概念, 探讨了水均衡法确定 浅层地下水安全开采动态控制水位的方法, 并以 阜阳市区为例, 确定了阜阳市区浅层地下水安全 开采动态控制水位, 对地下水实时动态监督管理 具有指导作用。 2 地下水安全开采动态控制水位 地下水安全开采控制水位是指使地下水系统 保持健康状态, 水量可持续更新, 且不产生不良后 果时相应的地下水位。地下水位的变动直接或间 接影响地下水系统的稳定性, 由于水文气象和地 下水开采的时空分布不均, 地下水安全开采控制 水位不是一个固定值, 而是随年内、 年际地下水补 给量、 排泄量、 储存量等要素变化而动态变化, 因 此也称地下水安全开采动态控制水位, 简称动态 控制水位。 根据浅层地下水利用特点, 浅层地下水安全 开采动态控制水位包括控制水位上限、 控制水位 下限、 上限预警水位和下限预警水位四个指标, 各 指标均随时间变化, 各动态控制水位示意图见图 1。控制水位上限和控制水位下限是指在设计天 然补给条件下, 地下水在相应各时段的最高、 最低 控制水位。上限预警水位和下限预警水位是在改 变降水保证率后, 地下水位在相应时段的最高、 最 低控制水位。在实际应用过程中, 当地下水位高 于控制水位上限时, 需要加强地下水位监测, 适当 图1 动态控制水位示意图 F i g . 1 S c h e m a t i c d i a g r a m o f d y n a m i c c o n t r o l w a t e r l e v e l 减少地表水利用量, 增大地下水开采, 必要时需排 水除涝; 当地下水位在控制水位上、 下限之间时, 该地区地下水可以按照核定的安全开采量正常开 采; 当地下水位低于控制水位下限时, 则无法满足 地下水系统安全, 必须减小开采量, 重点监测, 必 要时禁采, 防止超采产生不良的环境地质问题。 3 动态控制水位确定方法 3. 1 基本原理 水均衡法也称水量平衡法或水量均衡法, 是 全面研究某一地区( 均衡区) 在一定时段内( 均衡 期) 地下水补给量、 排泄量和储存量之间数量转化 关系的平衡计算。根据质量守恒定律, 在任何均 衡区, 任一均衡期内, 含水层系统的补给量与排泄 量之差恒等于该系统中地下水储存量的变化量。 任一均衡期内水均衡表达式为 It-Qt-Q安=ΔWt ( 1) 式中, It为均衡期内均衡区的地下水平均总补给 量, 1 0 8 m 3; Ot为均衡期内均衡区的地下水平均 天然总排泄量, 1 0 8 m 3; Q安为均衡期内均衡区的 安全开采量, 1 0 8 m 3; ΔWt为均衡期内均衡区的地 下水储存量变化量, 1 0 8 m 3。 含水层作为天然储水系统, 其作用如同地下 蓄水工程, 也被广泛称为地下水库。地下水库能 将丰水期的补给水量储存起来, 通过地下水取水 工程的有计划开采, 控制和调节地下水库的排泄 量。当补给量大于排泄量时, 地下水库蓄水量增 加, 水位上升, 反之当补给量小于排泄量时, 地下 水库蓄水量减小, 地下水位下降。地下水库蓄水 量变化过程的计算类似于地表水库的径流调节计 算, 地下水库内蓄水量的变化将引起其地下水位 相应的改变, 二者之间的关系为 ΔWt=0. 0 1(Ht-Ht- 1) μA ( 2) 式中,Ht-1、 Ht分别为均衡期初、 末地下水位,m; μ 为均衡区含水层的给水度; A为均衡区面积,k m 2。 3. 2 计算步骤 地下水的历史开采过程和实际需水过程是确 定地下水开采量年内分配的依据。研究地下水库 的年内调节过程, 采用水均衡法进行年调节计算, 确定浅层地下水安全开采动态控制水位, 具体计 算步骤如下。 步骤1 选择水文气象资料、 水位监测资料 充足, 具有连续长系列年份资料的监测站点作为 代表性观测站, 进行水均衡分析。 步骤2 调查研究不同部门地下水的需水特 点, 确定年内典型地下水开采过程, 核定地下水安 全开采过程, 为地下水安全开采动态控制水位的 确定提供基础。 步骤3 根据式( 1) 、 (2) , 以供水期末T末 ( 供 水期末根据地下水位动态确定) 的水位H末为推 算起点, 由地下水月水均衡方程推算第T 末-1 个 时段的末水位。 步骤4 采用水均衡法, 重复步骤3逆时序 逐月推算水位至供水期初相应的地下水位, 逐年 逆时序推算年内地下水位过程, 得到各典型年全 年按照安全开采量开采时的逐月地下水位过程。 步骤5 将上述各典型年的逐月地下水位过 程绘制于同一张图上, 纵坐标为地下水位, 横坐标 为时间( 月份) , 去除特丰和特枯年份, 取其上包线 即得到地下水安全开采控制水位上限; 取下包线, 并经适当修正, 则得到地下水安全开采控制水位 下限。改变降水保证率, 重复上述计算, 取其上包 线和下包线分别为上限预警水位和下限预警水位。 4 实例 4. 1 研究区概况 阜阳市位于安徽省西北部, 淮河北岸, 属于暖 温带半湿润气候区, 有明显的过渡性特征, 年平均 气温1 4. 6 ℃, 多年平均(1 9 5 6~2 0 1 0年) 降水量 8 9 1. 3 mm, 降水年内分配极不均匀, 主要集中 在6~8月。地势平坦, 整个阜阳市地势总体自西 北向东南倾斜, 平均地面坡降为1∶8 0 0 0。区内 广泛沉积了巨厚的具多层结构的河湖相沉积物。 其中第四纪地层发育齐全, 分布全区, 地下水类型 为单一的松散岩类孔隙水, 浅层地下水赋存于埋 深小于5 0m 的全新世、 晚更新世地层中, 与大气 降水、 地表水关系密切, 浅层含水层组广泛分布于 全区, 岩性以灰黄、 棕黄色粉砂、 细砂和中细砂为 主, 部分地区存在亚砂土和亚粘土。含水层总体 渗透性和富水性好, 渗透系数在0. 5~1 5. 0m/d 之间, 大部分地区单井涌水量大于5 0 0m 3/ d; 地 下水位普遍较浅, 埋深1. 2 4~5. 8 4m, 埋深变化 与地形总趋势基本一致, 自西北至东南埋深逐渐变 小, 见图2。 4. 2 阜阳市区动态控制水位 阜阳市下辖颍州、 颍泉、 颍东三区和临泉、 太 和、 阜南、 颍上四县及省直辖县级市界首市, 阜阳 市区包含颍州、 颍泉、 颍东三区, 选取阜阳市区代 表性地下水位观测井周棚站, 进行地下水均衡计 算, 计 算 时 间 为 1 9 6 0~2 0 1 0 年。阜 阳 市 区 面 46 水 电 能 源 科 学 2 0 1 4年 第3 2卷第1 0期刘 波等 浅层地下水安全开采动态控制水位确定 图2 阜阳市2 0 1 0年平均浅层地下水位埋深分布 F i g . 2 T h e a v e r a g e s h a l l o w g r o u n d w a t e r d e p t h d i s t r i b u t i o n o f F u y a n g C i t y i n 2 0 1 0 积1 6 8 4k m 2, 扣 除不透水面积, 计 算面积 为 1 5 2 1k m 2, 经计算, 阜阳市区多年平均地下水总补 给量为3. 2 7 71 0 8 m 3。参考已有的安全开采量 研 究[ 6,7],综 合 确 定 阜 阳 市 区 安 全 开 采 系 数 为0. 4 2 9, 采用安全开采系数法确定地下水安全 开采量, 由此计算安全开采量为1. 4 0 61 0 8 m 3。 根据研究区不同部门地下水的需水特点, 确定地 下水开采过程, 采用等比例法拟合得到地下水安 全开采过程, 使得年内地下水开发利用实现最优 化, 阜阳市区年内地下水安全开采量过程见表1, 阜阳市区安全开采量为1. 4 0 61 0 8 m 3。浅层地 下水开采主要用于农业灌溉, 其次为生活用水和 少量业用水。开采主要集中在3~5月, 因此年最 低水位一般出现在4~5月,6~9月为雨季, 开采 量较少, 1 2月至次年2月农业灌溉基本停止, 少 量开采浅层地下水用于居民生活和工业生产。 表1 2 0 1 0年阜阳市区( 周棚站) 逐月水均衡过程 及水位推算表 T a b . 1 T h e m o n t h l y w a t e r b a l a n c e a n d w a t e r l e v e l p r e d i c t i o n t a b l e o f Z h o u p e n g S t a t i o n a t F u y a n g C i t y i n 2 0 1 0 月 份 面积 /k m2 总补 给量 天然 排泄量 安全 开采量 储存量变化 增减 水位变幅/m 升降 水位 /m 5 1 6 8 4 0. 3 1 5 0. 1 8 0 0. 2 0 2 0. 0 6 7 0. 1 1 2 4. 8 5 6 1 6 8 4 0. 8 8 4 0. 2 0 6 0. 0 9 0 0. 5 8 8 0. 9 7 2 5. 8 2 7 1 6 8 4 0. 3 8 1 0. 1 6 0 0. 1 1 8 0. 1 0 3 0. 1 7 2 5. 9 9 8 1 6 8 4 0. 2 4 9 0. 1 8 6 0. 1 1 8 0. 0 5 5 0. 0 9 2 5. 9 0 9 1 6 8 4 0. 3 6 0 0. 1 2 7 0. 1 1 8 0. 1 1 5 0. 1 9 2 6. 0 9 1 0 1 6 8 4 0. 1 5 0 0. 1 2 1 0. 0 9 0 0. 0 6 1 0. 1 0 2 5. 9 9 1 1 1 6 8 4 0. 2 5 6 0. 0 7 7 0. 0 8 8 0. 0 9 1 0. 1 5 2 6. 1 4 1 2 1 6 8 4 0. 0 2 5 0. 1 7 7 0. 0 6 0 0. 2 1 2 0. 3 5 2 5. 7 9 1 1 6 8 4 0. 0 7 6 0. 1 5 5 0. 0 6 0 0. 1 3 9 0. 2 3 2 5. 5 6 2 1 6 8 4 0. 0 8 4 0. 1 6 9 0. 0 6 0 0. 1 4 5 0. 2 4 2 5. 3 2 3 1 6 8 4 0. 3 0 1 0. 2 8 9 0. 2 0 0 0. 1 8 8 0. 3 1 2 5. 0 1 4 1 6 8 4 0. 0 8 5 0. 2 2 2 0. 2 0 2 0. 3 3 9 0. 5 6 2 4. 4 5 总计1 6 8 4 3. 1 6 6 2. 0 7 1 1. 4 0 6 0. 3 1 1 注 总补给量、 天然排泄量、 安全开采量、 储存量变化单位均为1 0 8 m3。 自供水期末开始, 按照步骤3逆时序推算不 同水平年的逐月地下水位过程, 以2 0 1 0年阜阳市 区( 周棚站) 的水均衡过程为例, 计算时段为月, 根 据研究区特点, 给水度( μ) 取经验值0. 0 3 6, 起点 水位取年内最低月份4月水位的平均值2 4. 4 5 m, 具体计算过程见表1, 得到年内地下水位过程。 对不同水平年代表性观测站, 按照步骤5统 计安全开采过程的逐月地下水位指示线; 去除特 丰和特枯年份, 取其上包线即得到地下水安全开 采控制水位上限; 取下包线, 并经适当修正, 则得 到地下水安全开采控制水位下限, 见图3。由图3 可知, 地下水安全开采动态控制水位随时间变化, 9月地下水控制水位上限和下限均达到最高点, 水位分别为2 6. 8 7、2 5. 4 3m。 图3 阜阳市区( 周棚站) 地下水安全开采动态控制水位过程 F i g . 3 T h e d y n a m i c g r o u n d w a t e r s a f e t y m i n i n g c o n t r o l l e v e l p r o c e s s o f Z h o u p e n g S t a t i o n i n F u y a n g C i t y 改变降水保证率, 取降水量值处于丰水年 ( 2 5%) 值和枯水年(7 5%) 值之间的水平年, 重复 上述水均衡计算过程, 得到所选水平年的年内地 下水位过程, 取各地下水位过程线的上包线, 即为 地下水安全开采控制水位的上限预警线; 各地下 水位过程线的下包线, 即为地下水安全开采控制 水位的下限预警线, 见图4。由图4可知, 根据确 定的动态控制水位, 制定不同的管理方案。当代 表性观测井的地下水位处于上限预警水位和下限 预警水位之间时, 则认为当前的开采方案合理, 此 时地下水开采以高效节水为主; 当地下水位处于 上限预警水位和控制水位之间, 或处于下限预警 水位和控制水位下限之间时, 则认为地下水开采 方案相对合理, 需要加强监测预警, 采取适当限采 或排水措施, 使得地下水位能够恢复至预警水位 区间内; 当地下水位高于控制水位上限, 或低于控 制水位下限时, 必须采取排水或禁采措施, 使地下 水位得以恢复。 图4 阜阳市区( 周棚站) 地下水安全开采动态 控制水位和预警水位 F i g . 4 T h e d y n a m i c g r o u n d w a t e r s a f e t y m i n i n g c o n t r o l l e v e l a n d w a r n i n g l e v e l o f Z h o u p e n g S t a t i o n i n F u y a n g C i t y 56 5 结论 a . 提出了地下水安全开采动态控制水位概 念, 并结合浅层地下水特点, 分析了水均衡法确定 动态控制水位的原理和步骤。 b . 采用水均衡法对阜阳市区( 周鹏站) 进行了 实例分析, 结果表明年内浅层地下水安全开采动 态控制水位和预警水位是随时间变化的, 为实现 地下水安全开采及动态监测提供了可靠依据。 参考文献 [1] 袁长极. 地下水临界深度的确定[J]. 水利学报,1 9 6 4 (3) 5 0-5 3. [2] 殷丹, 许春东, 束龙仓, 等. 淮北市岩溶地下水可持 续开采量及临界水位的确定[J].水电能源科 学, 2 0 1 2,3 0(7) 2 5-2 8. [3] 荣丽杉, 束龙仓, 王茂枚, 等. 合理地下水生态水位 的估算方法研究 以塔里木河下游为例[J].地下 水,2 0 0 9, 3 0(1) 1 2-1 5,4 3. [4] 谢新民, 柴福鑫, 颜勇, 等.地下水控制性关键水位 研究初探[J]. 地下水, 2 0 0 7,2 9(6) 4 7-5 0,6 4. [5] 陈文芳. 中国典型地区地下水位控制管理研究[D]. 北京 中国地质大学,2 0 1 0. [6] J I Y,S HU L,WANG Z,e t a l . N e w M e t h o d f o r E- v a l u a t i n g S a f e Y i e l d a n d E x p l o i t a t i o n P o t e n t i a l o f S h a l l o w G r o u n d w a t e r i n H u a i b e i P l a i n,A n h u i P r o v- i n c e[J]. I AH S-A I S H P u b l i c a t i o n,2 0 1 1,3 5 02 0 5-2 1 1 . [7] 王振龙, 鲁程鹏, 刘猛.地下水安全开采量的概念与 评价方法研究[J] . 水文, 2 0 1 0,3 0(2) 1 4-1 9,1 3. D e t e r m i n a t i o n o f S a f e t y M i n i n g D y n a m i c C o n t r o l L e v e l o f S h a l l o w G r o u n d w a t e r L I U B o,L I U G e,S HU L o n g c a n g (C o l l o g e o f H y d r o l o g y a n d W a t e r R e s o u r c e s,H o h a i U n i v e r s i t y,N a n j i n g 2 1 0 0 9 8,C h i n a) A b s t r a c tI n o r d e r t o s t u d y t h e s a f e t y m i n i n g o f g r o u n d w a t e r d e e p l y,t h i s p a p e r p r o p o s e s t h e c o n c e p t o f s a f e t y m i n i n g d y n a m i c c o n t r o l l e v e l o f g r o u n d w a t e r,i n c l u d i n g f o u r i n d e x e sc o n t r o l w a t e r l e v e l u p p e r l i m i t,c o n t r o l w a t e r l e v e l l o w e r l i m i t,u p p e r l i m i t w a r n i n g w a t e r l e v e l a n d l o w e r l i m i t w a r n i n g w a t e r l e v e l . E a c h i n d e x c h a n g e s w i t h t i m e . A n d t h e n i t i n- t r o d u c e s t h e b a s i c t h e o r y a n d p r o c e d u r e o f w a t e r b a l a n c e m e t h o d t o d e t e r m i n e d y n a m i c c o n t r o l w a t e r l e v e l i n d e t a i l . T a k - i n g F u y a n g C i t y a s a n e x a m p l e,t h e r e p r e s e n t a t i v e s t a t i o n s(Z h o u p e n g s t a t i o n)i s s e l e c t e d t o t a k e w a t e r b a l a n c e c a l c u l a t i o n a n d i t g e t s t h e s a f e t y m i n i n g d y n a m i c c o n t r o l l e v e l p r o c e s s o f s h a l l o w g r o u n d w a t e r i n s t u d y a r e a .T h e r e s u l t s s h o w t h a t t h e i n n e r c o n t r o l w a t e r l e v e l u p p e r l i m i t,c o n t r o l w a t e r l e v e l l o w e r l i m i t,u p p e r l i m i t w a r n i n g w a t e r l e v e l a n d l o w e r l i m i t w a r n i n g w a t e r l e v e l a l l c h a n g e s w i t h t i m e ,w h i c h c a n p r o v i d e a r e l i a b l e b a s i s f o r g r o u n d w a t e r d y n a m i c s u p e r v i s i o n a n d m a n a g e m e n t . K e y w o r d ss h a l l o w g r o u n d w a t e r;w a t e r b a l a n c e m e t h o d;d y n a m i c c o n t r o l w a t e r l e v e l;w a r n i n g w a t e r l e v e l 櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀櫀 ( 上接第5 9 页) [7] 黄金良, 杜鹏飞, 何万谦, 等.城市降雨径流模型的参 数局部灵敏度分析[J].中国环境科学, 2 0 0 7,2 7(4) 5 4 9-5 5 3. [8] L e w i s A R o s s m a n . S t o r m W a t e r M a n a g e m e n t M o d e l U s e r’s M a n u a l(V e r s i o n 5. 0) [M]. W a s h i n g t o n D C, U SE P A,2 0 0 9. [9] 赵冬泉, 王浩正, 陈吉宁, 等.城市暴雨径流模拟的参数 不确定性研究[ J].水科学进展,2 0 0 9,2 0(1) 4 5-5 1 . [1 0]上海市政工程设计研究总院, 北京市市政工程设计 研究总院, 中国市政工程东北设计研究院, 等.室外 排水设计规范(G B 5 0 0 1 4-2 0 0 6) [ S].北京 中国计划 出版社,2 0 0 6. [1 1]张磊.平原感潮河网区城市防洪规划中的水文计算 方法研究[D].南京 河海大学, 2 0 0 5. A p p l i c a t i o n o f S WMM M o d e l i n U r b a n D r a i n a g e a n d F l o o d C o n t r o l C a l c u l a t i o n o f T i d a l R i v e r N e t w o r k Z HANG B i n 1, D A I X i a n b o 2, XU X i a n g y a n g 1, Z HAO H u i 1, CHE N Q i n g s h a 1 (1. C o l l e g e o f H y d r o l o g y a n d W a t e r R e s o u r c e s,H o h a i U n i v e r s i t y,N a n j i n g 2 1 0 0 9 8,C h i n a; 2. Z h a n g j i a g a n g C i t y W a t e r C o n s e r v a n c y B u r e a u,Z h a n g j i a g a n g 2 1 5 6 0 0,C h i n a) A b s t r a c tF o r t h e r u n o f f c h a r a c t e r i s t i c s o f u r b a n t i d a l r i v e r n e t w o r k a r e a a n d c o n s i d e r i n g t h e i n f l u e n c e o f r i v e r d r a i n- a g e a n d d r a i n a g e p i p e,SWMM m o d e l i s u s e d t o g e n e r a l i z e t h e p i p e d r a i n a g e s y s t e m o f t h e s t u d y a r e a .U n i t h y d r o g r a p h m e t h o d o f d i m e n s i o n l e s s t i d e l e v e l i s u s e d t o d e t e r m i n e t h e t i d e l e v e l p r o c e s s o f t h e a r e a a l o n g t h e Y a n g t z e R i v e r . T a k i n g c a l c u l a t i o n r e s u l t o f r e g i o n a l f l o o d c o n t r o l a s c o m p u t i n g b o u n d a r y c o n d i t i o n s o f t h e o u t e r r i v e r a n d i n p u t t i n g s h o r t d u r a t i o n r a i n s t o r m p r o c e s s t h r o u g h t h e d e s i g n o f s a m e f r e q u e n c y s c a l i n g,t h i s p a p e r s i m u l a t e d a n d c a l c u l a t e d t h e u r b a n d r a i n a g e a n d w a t e r l o g c o n t r o l p r o c e s s . T h e r e s u l t s s h o w t h a t t h e s l u i c e e n g i n e e r i n g c a n c o n t r o l t h e f l o o d c o n t r o l w a t e r l e v e l o f t h e p a r t i t i o n e f f e c t i v e l y;f l o o d c o n t r o l w a t e r l e v e l m a y b e a d j u s t e d b y t h e r a t e o f r i v e r w a t e r i n s o m e c e r t a i n a n d t h e h i g h i n i- t i a l w a t e r l e v e l o f t h e r i v e r h a s i m p a c t o n t h e p i p e d r a i n a g e a n d i m p o u n d e d s u r f a c e w a t e r . K e y w o r d sSWMM m o d e l;t i d a l r i v e r n e t w o r k;s t o r m r u n o f f;f l o o d c o n t r o l a n d d r a i n a g e;p i p e d r a i n a g e 66 水 电 能 源 科 学 2 0 1 4年