1、 1 流域地表径流系数的计算方法研究 摘要 ：径流系数是描述降雨和径流关系的重要参数 ,在雨洪控制利用系统的理论研究、 规划、 设计计算中应用广泛 ,在流域或区域的雨水径流总量、 径流峰流量、 流量过程线以及非点源污染物总量、 各设施规模的计算中也起着极其重要的作用。由于 径流系数有 着 不同的含义 ,其相应 的统计计算方法、适用条件、应用目的和取值 不尽相同。而且要获得流域的径流系数通常是比较困难的，在一些特殊流域基本上很难获得能满足要求的径流实测资料，尤其在多年平均径流量的计算中实测数据资料往往相当缺乏，在这样的情况下有必 要利用一些特殊的方法去满足工程建设对水文数据的需求。本文综合了大量
2、的数据以及列举了多个例子，详细地介绍了不同情况下径流系数的推求方法，并在此基础上研究总结提出了过程中发现的一些问题和心得。 关键词 ：流域 径流量 降雨量 径流系数 一 引言 流域径流系数是指 同一流域面积、同一时段内径流量与降水量的比值，以小数或百分数表示。 计算式为： R/P，式中 为径流系数， R 为径流深度， P 为降水深度。 值变化于 0 1 之间，湿润地区 值大，干旱地区 值小。我国台湾地区河流年平均径流系数 0.7，表明径流十分丰富 ；径流贫乏的海滦河平原，年平均径流系数仅有 0.1。 根据计算时段的不同，可分为 瞬时雨量径流系数、雨量径流系数、 年径流系数、多年平均径流系数 等
3、。 径流系数综合反映流域内自然地理要素对降水 径流关系的影响。 瞬时雨量径流系数是指某一特定的流域或汇水面上 ,降雨期间随时间变化的径流厚度和降雨厚度之间的瞬时变化关系 ,是一个动态的变量 ,这个意义上的径流系数就是瞬时雨量径流系数。雨量径流系数是指降雨时 ,在某一汇水面上产生的径流量 (厚度 )和降雨量 (厚度 )的比值 ,一般用于估计一场降雨在某一汇水区域内单位面积产 生的平均径流厚度。年径流系数和多年平均径流系数反映了流域降雨厚度和径流厚度长时间的关系 ,是一个累积结果。在各种径流系数中应用较为广泛的是年径流系数和多年平均径流系数。径流系数的计算主要是要计算流域相应时间段内径流量与降雨量
4、。 二 径流量的计算 （一） 年径流量的计算 流域年降雨次数为 n 次，且每次降雨所产生的径流量均有实测数据资料，则流域的年径流量可按下式计算。 Q=n1i Qi（ 1） 式中 Q 流域年径流总量（ mm）； Qi 第 i 次降雨产生的径流量（ mm）。 （二） 多年平均径流量的计算 1.有长期实测资料的多年平均径流量的计算 所谓的有长期实测资料，是指实际观测的年数 n 在 20 年以上。它包括有丰、平、枯水年的观测资料，由它计算的径流量多年平均值基本上是稳定的。在这种情况下，可以由下式（ 2）计算径流量的多年平均值，以此值代表多年平均径流量，即： 2 Q 平均 =n1 n1i Qi（ 2）
5、式中 ： Q 平均 流域多年平均径流量（ mm）； Qi 序列号为 i的年份流域径流径流量（ mm）。 2有短期实测资料的多年平均径流量的计算 若实测系列长度小于 20a，用上述方法计算误差会太大，为了提高计算的精度，保证计算结果的可靠性，就必须对实测径流系列进行相关展延，然后用上述方法计算。目前常用的方法有两种。 （ 1） 利用径流资料插补延展系列。 该方法的关键是找到和流域自然地理条件相似 的邻近的参照流域，且参照流域具有长序列实测年径流序列。利用相关分析验证研究区流域和参照流域是否具有相关关系，如果两流域相关关系密切，则可以用两者的相关趋势函数插补延展研究区的短期年径流量实测资料序列。当
6、年径流量资料序列很短，不足以建立年相关时，也可以先建立月相关序列，插补延展月径流量，然后计算年径流量，但用月径流量相关来插补延展年径流量时，会使误差累积，精度较低。 （ 2） 利用降雨量资料插补延展序列。 如果研究区的上下游或邻近地区找不到具有长序列的径流量资料的参照流域，而在研究区流域或邻近参照流域有较长序列 的年降水量资料，则可以选择降水量作为参照变量与研究区短期实测径流量资料进行相关分析，如果两个序列之间相关关系密切，则可以利用较长实测年降水量资料进行插补延展研究区的径流量序列。 如下表列举了贵州省大煤矿流域 1967 年到 2004 年的降雨量实测系列，而该流域只有1967年到 198
7、4年及 1991年到 2004年的实测径流序列，因此需要用降雨量序列插补延展 1984年到 1991 年径流量序列序列。插补结果如下表。 表 1 贵州省大煤矿流域多年雨量径流量数据 年份 雨量 /mm 径流深 /mm 年份 雨量 /mm 径流深 /mm 19671968 1035.7 846.8 19681969 1006.0 470.3 19691970 974.3 457.9 19701971 1133.4 532.7 19711972 850.5 399.7 19721973 1006.0 470.3 19731974 992.9 466.7 19741975 987.6 464.2 1
8、9751976 974.6 458.1 19761977 1013.9 476.5 19771978 973.0 457.3 19781979 958.6 450.5 19791980 1026.5 482.5 19801981 1006.5 473.1 19811982 864.8 406.5 19821983 1082.7 508.9 19831984 1025.3 481.9 19841985 880.2 413.5 19851986 973.0 457.0 19861987 893.9 419.9 19871988 895.1 420.5 19881989 860.2 404.1 19
9、891990 869.3 408.4 19901991 766.1 360.0 19911992 1056.7 496.6 19921993 827.9 389.1 19931994 901.9 423.9 19941995 966.8 454.4 19951996 1009.5 474.5 19961997 1048.8 492.9 19971998 824.0 387.3 19981999 946.7 444.9 20002001 1068.3 502.1 20012002 1038.3 488.0 20022003 843.1 396.3 20032004 783.1 368.1 3 y
10、 = 0.4691x + 0.5978R2= 0.9993350370390410430450470490510530550800 900 1000 1100 1200降雨量( m m )径流量(mm)图 1 径流量与降雨量相关曲线图 3无实测资料时河流径流量推求方法 （ 1）等值线图法 把相同数值的点连接起来的线叫等值线。某一流域的水文特征值的等值线图即可反映出该流域水文特征值的地理分布规律。闭合流域多年径流量的主要影响因素是气候因素，而气候因素有地区性，即降雨量与蒸发量具有地理分布规律，同理，受降雨量和蒸发量影响的多年平均年径流量也具有地理分布规律。因此可利用这一特点绘制多年平均年径流量的
11、等值线图，并可以用它来推算无实测资料流域的多年平均年径流量。 应用等值线图推求多年平均年径流深时，先在图上勾绘出 研究流域的分水线，再找出流域的形心，而后根据等值线内插读出形心处的多年平均年径流深值。如果流域面积较大或地形复杂，等值线分布不均匀，也可用加权平均法推算，即： Y =（ y1f1 +y2f2 + +ynfn） /F （ 3） 式中 ： y1 相邻两径流深等值线的平均值； f1 相邻两等值线间面积； F 流域总面积。 （ 2） 水位比拟法 如前所述，水文现象具有地区性，如果某几个流域处在相似的自然地理条件下，则其水文现象具有相似 的发生、发展、变化规律和相似变化特点。与研究流域有相似
12、自然地理特征的流域称为相似流域（即参证流域）。水文比拟法就是以流域间的相似性为基础，将相似流域的水文资料移用至研究流域的一种简便方法。其中移用相似流域研究资料的方法较多，如选择相似流域的径流模数、径流深度、径流量、径流系数以及降水径流相关图等。但是，地球上不可能有两个流域完全一致，或多或少都存在一些差异，倘若相似流域与研究流域之间仅在个别因素上有些差异时，可以考虑不同的修正系数加以修正。 若研究流域与相似流域的气象条件和下垫面因素基本相似，仅流域面积有所 不同，这时只考虑面积的影响，则研究流域的正常年径流量有如下关系式： Q 研 /F 研 =Q 相 /F 相 。如果使用径流深或径流模数，则不需
13、要修正即可使用。 若两流域的年降水量有不同时，则 Q 研 /P 研 =Q 相 /P 相 。 4 式中 ： P 研 和 P 相 分别为研究流域、相似流域的年降水量。 （ 3） 径流系数法 当小流域内（或附近）有年降水量资料，且降水量与径流关系密切时，可利用多年平均降雨量与径流量间的定量关系计算年径流量，即利用年降雨量的多年平均值乘以径流系数推求多年平均径流量间的流量，可由下式计算之： W=1000 C P F （ 4） 式中： W 多年平均径流总量， m； C 该地区年径流系数，与研究区植被、地形，地质、主河道长度等因素有关，可通过调值并参考省、地水文手册确定； P 研究地区多年平均降雨量（ m
14、m），可从省、地区的水文手册查出，或向附近水文站、雨量站查询； F研究流域的集水面积（ km ）。 （ 4） 水文查勘法 对于完全没有资料，也找不到相似流域的小河或间歇性河流，此时可进行水文查勘，收集水文资料， 进行正常年径流量的估算。这项任务一般是通过野外实地查勘访问，了解多年期间典型水位过程线，河道特性，建立水位流量关系曲线，从而推算出近似的流量过程线，并估算其正常年径流量。水文查勘工作，不仅对完全无资料的小河有必要，就是对有资料的水流域也是不可缺少的。 （ 5） 经验公式法 经验公式都是根据各地实测资料分析得出的，这些径流公式一般可以在当地的水文手册中查得。如： Q=KFn （ 5） 式
15、中： Q 多年平均径流量（ m3）； F 流域面积（ km3）； K、 n 分别为地区性参数，取值可查水文手册。 三 降水量的计算 （一） 一次降水过程中降水量的计算方法 1 算术平均法 算术平均法是将该流域各站测得的同期雨量相加后 ,除以总站数 ,即为流域面雨量。其数学表达式为 : P 平均 = n1 n1i Pi（ 6） 式中： n 总站数； Pi 为各站同期雨量（ m）。 2 等雨量线法 根据流域内各测站实测的雨量资料绘出等雨量线 ,然后用求积仪或其它方法求各相邻两等雨量线间的面积 ,再分别乘以各相邻两等雨量线雨深的平均值 ,即得该面积上的降水总量。将各面积上的降水总量相加 ,除以全流域
16、的总面积 ,即得流域的面雨量。用数学关系式表示为 : P =A1 n1i 2PP i1-i ai （ 7） 式中： Pi-1， Pi 第 i-1 和第 i 条等雨量线代表的降雨量， i=1,2,3， n； ai 第 i-1 和第 i 条等雨量线之 间的面积， i=1,2,3， n； A 总面积； 5 n 等雨量线的条数。 表 2 邵阳地区 2000 年各站年降水量 站名 东经 北纬 降雨量 / ( ) / ( ) /mm 站名 东经 北纬 降雨量 / ( ) / ( ) /mm 汪家田 110. 58 26. 52 1 744. 3 武 冈 110. 60 26. 72 1 391. 2 鸦槎
17、铺 110. 80 26. 75 1 313. 8 武 阳 110. 32 26. 73 1 514. 2 瓦屋塘 110. 33 26. 93 1 582. 7 红 岩 110. 50 26. 90 1 548. 3 高 沙 110. 68 26. 97 1 326. 9 黄 桥 110. 85 27. 03 1 369. 2 洗 马 110. 47 27. 32 1 768. 0 栗山界 110. 45 27. 15 1 693. 0 长 塘 110. 52 27. 13 1 563. 4 金屋塘 110. 35 27. 07 1 982. 3 金盆形 110. 45 27. 03 1 6
18、58. 7 洞 口 110. 55 27. 07 1 541. 5 山 门 110. 68 27. 23 1 495. 8 苏家洞 110. 80 27. 32 1 396. 1 匡家铺 110. 93 27. 53 1 423. 4 隆 回 110. 98 27. 13 1 329. 9 枯 桑 111. 03 26. 93 1 276. 6 麻 林 110. 63 26. 47 1 623. 6 新 宁 110. 87 26. 45 1 599. 4 回龙寺 111. 10 26. 75 1 264. 2 三门江 111. 23 26. 93 1 425. 8 罗家庙 111. 25 27
19、. 02 1 333. 7 檀木塘 111. 85 27. 12 1 223. 6 邵 东 111. 73 27. 23 1 301. 1 官桥铺 111. 62 27. 18 1 292. 2 五丰铺 111. 47 26. 88 1 429. 9 诸家亭 111. 58 27. 03 1 284. 9 短陂桥 111. 57 27. 27 1 330. 5 茅 坪 111. 47 27. 23 1 390. 6 刘付冲 111. 47 27. 43 1 946. 2 高 坪 111. 12 27. 48 1 316. 7 岩口铺 111. 23 27. 22 1 360. 9 大河滩 11
20、1. 38 27. 50 1 243. 6 双 林 110. 93 27. 73 1 762. 7 白茅坪 110. 37 26. 27 1 402. 9 界 溪 110. 13 26. 38 1 880. 3 黄 桑 110. 05 26. 43 1 691. 0 黄土头 109. 98 26. 50 1 867. 3 朝 仪 109. 97 26. 58 1 863. 9 丰 家 110. 07 26. 57 1 659. 4 党 坪 110. 08 26. 60 1 570. 3 小沙江 110. 75 27. 52 1 551. 8 图 2 用距离倒数加权法绘等值线 3 数值法 6 通
21、常数值法是指森 泰多边行法或三角形法 ,它实际上是以权重系数决定的计算方法。三角形法比森泰多边行法计算方便。三角形法的具体做法是先在流域内自先设定均匀分布的网格点 ,再将各网格点置于相应的三角形的重心上 ,以此求出代表此三角形面积上的平均雨深。然后各三角形区域的平均雨深乘上三角形面积 ,得出三角形区域的降水总量 ,再将各三角形区域降水总量相加 ,除以流域总面积 ,即得流域面雨量 ,用数学关系式表示为 : P =n1i/AaP ii （ 8） 式中 : Pi 任意三角形区域平均雨深 ; ai 任意三角形面积 ; A 流域总面积。 443322 11aaaaPPPP图 3 泰森多边形 (二 ) 年
22、降雨量的计算方法 如果流域一年中所有场次的降雨量均已测出，则流域的年降雨总量可按下式计算： P=n1i iP（ 9） 式中： Pi 第 i 次降雨的雨量； n 一年中的降雨次数。 （三） 多年平均降雨量的计算 若流域有长期年降雨量资料，则流域的多年平均降雨量为： P =n1 n1i iP，一般要求实际观测年数 n 在 20a 以上。 若流域只有短期年降雨量实测资料，则可以用插补法延长实测序列，通常所用的插补法有： 1 对选用的 雨量站 , 若本站与参证站某时段 (如 : 年、月 ) 降水量数值较为接近 , 则说明成雨条件大体一致 , 降水量在面上分布比较均匀 ,这时 , 各参证站降水量算术平均
23、值作为插补站相应时段降水量。 2 非汛期降水量较少 , 各年变化不大时 , 采用本站同月降水量的历年平均值插补缺测月份。 3 年降水量缺测时 , 若参证站与插补站有较长的平行观测资料 , 且降水成因一致 , 根据相临站的长系列资料用相关分析法插补或延长 , 也可用降水量等值线图内插。 四 总结和建议 7 1径流系数是受降雨过程、 土壤性质、 流 域坡度、 土地利用类型以及先前湿润状况等因素综合影响的代表径流和降雨之间关系的系数 ,以不同的含义及大小具体体现在降雨过程、 产流过程和汇流过程中。 2径流系数的计算关键是对径流量与降雨量的计算，而径流量与降雨量的计算关键是要获得实测数据资料，因此实测
24、径流数据和降雨数据是非常重要的。 3场 (次 )雨量径流系数可用于计算单场降雨的径流总量。在应用中 ,一方面可对开发前、 后设计重现期的降雨事件的雨水径流总量进行计算分析和制定控制对策 ;也可以粗略计算单场降雨的雨水资源量 ;峰流量径流系数可用于计算设计重现期 降雨事件的径流峰流量 ,除了用于传统的管渠设计计算 ,也用于对开发前后径流峰流量的控制。 4年均径流系数一方面可用于粗略计算年均径流总量或雨水资源总量 ,另一方面可间接估算非点源污染物年均总量。 5为准确反映降雨和径流之间的关系 ,首先需要明确不同径流系数的概念 ,对汇水面的特征参数、 气候因素、 不同条件下的产流和汇流进行监测和综合分
25、析 ,以区别和修正用于不同条件的不同径流系数。 6传统的城市雨水排放理念和方法主要针对一定设计标准下的峰流量排放 ,因此主要用峰流量径流系数。但现代城市新型雨洪控制利用 理念和方法则必须针对径流污染控制、 雨水径流总量和峰值的控制 (限制排放 )、 滞蓄利用等 ,因而需要根据目标和具体要求 ,采用不同的径流系数。对不同的雨水径流系数概念及其主要影响因素进行梳理和总结 ,主要为城市雨洪控制利用系统的理论分析、 规划、 设计计算中径流系数的合理确定提供参考。但是降雨和径流之间的关系非常复杂、 多变 ,在应用中要根据实际情况分辨主次因素 ,选择不 同的径流系数并确定其大小 ,其中一些问题 (如参数选
26、取、 确定等 )还需进一步深入研究。 参考文献 ： 1 徐恒力 .水资源开发与保 护 .地质出版社 .北京 .2001 年 8 月 .124-125. 2 芮孝芳 .水文学原理 .中国水利水电出版社 .北京 .2004 年 8 月 .52-54. 3 张红梅 .水文统计讲义 .水资源与环境学院 .北京 .2010 年 3 月 . 20-25. 4 何玉敬 丁国梁 . 无实测资料时河流径流量推求方法的适用性分析 . 新疆 昌吉 . 昌吉水文水资源勘测局 . 2009 年 9 月（上旬）第 08 卷第 25期 第 186 期 . 5 李学美 刘新 华任伯帜 . 无资料非闭合小流域矿区径流来水量计算
27、方法 . 湖南科技大学 土木工程学院 . 湖南 湘潭 . 2009 年 9 月第 24 卷第 3 期 . 6 于丽玲 卢昌 卢 昊 . 山区小流域径流系数计算方法探讨 . 伊春水文局 .黑龙江 . 2005 年 第 6 期 (第 33 卷 ). 7 雷祖培 彭国才 . 无资料地区径流计算中儿个技术问题的探讨 . 韶关市水利水电勘测设计院 . 广东 . 2000 年第二期 . 8 郭天亮 . 沟尾溪涝区整治工程地表径流系数计算 . 潮州市水利水电勘测设计院 . 2009 年 8 月 . 总第 219 期 . 9 唐宁远 车伍 潘国庆 . 城市雨洪控制利用的雨水径流系数分析 . 北京建筑工程学院 城市雨水与水环境研究中心 . 北京 . 2009 年 11 月第 25 卷第 22 期 . 10 范银贵 . 空间插值方法在绘制降水量等值线中的应用 . 湖南省邵阳水文水资源勘测局 . 湖南 邵阳 . 2002 年 6 月第 22 卷第 3期 .
