1、河流最小生态基础流量计算方法研究（张新华 李红霞 肖玉成 赵少华）摘要：为了能够计算几何断面复杂、泥沙淤积、河宽大、水深浅的河道生态基础流量，本文以渭河关中地区的生态基础流量计算为例，通过综合水力学中的湿周法和 R2CROSS 法建立了一套简单、适合管理的最小生态基础流量方法综合法。该方法计算出的 5 个站的生态基础流量占多年平均流量的比值都在 10.8%14.9%之间，满足生态基流的合理要求。通过与 Tennant 方法比较，结果表明：综合法计算确定的生态基础流量更为合理、是一种在管理上更值得推广应用的方法。关键词：生态需水量；湿周法；斜率法；曲率最大法；R2CROSS 法中图分类号：X17
2、1 文献标识码：A 文章编号：1672-3031（2011）01-0066-08Study on computational methods for minimum environmental flowsZHANG Xin-hua，LI Hong-xia，XIAO Yu-cheng，ZHAO Shao-huaAbstract：In order to determine the minimum ecological flow（MEF） for a river system with special characteristics of complicated river reach profil
3、es，sediment problems，very large width and shallow water depth，an integrated method was proposed to overcome some of limits in using wetted perimeter method and R2CROSS method. A case study was conducted in the Wei River，a tributary of the Yellow River. Results calculated by this integrated method in
4、dicate that the MEFs at the five hydrological stations are at a range of 10.8%14.9% of the multi-year mean flows，respectively. In comparison with Tennant method， therefore，we can get a conclusion that the integrated method is more suitable and reasonable for the determination of MEFs.Key words： ecol
5、ogical water demand； wetted perimeter method； slope method； curvature method and R2CROSS1 研究背景河流系统包括河流、湖泊及其相邻洪泛区。河流系统功能通常包括生态功能、环境功能和资源功能1。然而，在开发水资源的时候主要考虑到满足工业、农业和城镇生活用水等需求，忽视了河流系统功能得以正常发挥并维持其基本功能时其自身对水量的要求，从而导致河段内的天然流量逐渐减少、甚至出现断流的情况。另外，经济社会的高速发展导致大量的污水排放，使得有些河流虽然水量上没问题，但水质污染十分严重导致水生生物种群及其数量急剧下降，水体
6、严重富营养化，出现藻类或水华暴发等现象。为了维持河流系统健康，保障其基本功能能正常发挥，国内外学者为此开展了大量的研究工作，主要集中在河流生态环境需水上并取得了一定的研究成果。河流生态环境需水量计算内容主要包括：（1）生态基础流量；（2）输沙需水量；（3）自净需水量；（4）蒸发需水量；（5）防海水入侵的河道最小水量等2。其中的河流生态基础流量是维持河道功能不退化、特征水生生物物种能正常繁衍生息所必须维持的河道最小流量，本文重点分析讨论这一河道最小生态基础流量。然而，国内外针对这一流量的研究还未形成统一认识，基本理论和方法都还有待进一步完善。随着近几年的研究深入，国内外学者提出了多种计算方法，概
7、括起来主要包括 4 类方法：（1）水文学法；（2）水力学法；（3）栖息地法；（4）地形结构法等。在这几种方法中，以水文学法与水力学法最为常用。因此，本研究方法的重点在分析讨论水力学法中的湿周法并结合 R2CROSS 法和 Tennat 法，提出一套简单、实用、便于管理部门操作的方法综合法。2 计算方法2.1 湿周方法湿周法计算河道的最小生态需水量时要基于一个假设，即浅滩是临界的河流栖息地，而保护浅滩栖息地也能保护到其他的水生生物栖息地。通过河道的实测大断面和实测流量资料建立湿周流量关系曲线，在曲线中用斜率 1 法和曲率最大法确立变化点，该点对应的流量即为最小生态流量，这就是湿周法确立河道最小生
8、态流量的基本原理3-5。根据水力学中计算明渠均匀流的谢齐公式： Q = ACRJ 以及曼宁公式：C = 1/nR 1/6可以推出湿周流量关系式：式中：Q 为流量，m3/s；A 为过水断面面积，m2；P 为湿周，m；J 为水力坡度；n 为粗糙系数。在上式中，由于流量 Q、湿周 P 以及断面面积 A 都是未知的，因此直接推出流量和湿周的关系比较困难，而通过斜率 1 法和曲率最大法确定 PQ 关系曲线变化点需要知道 PQ 之间的函数关系。根据国内外对湿周法的研究，通常使用幂函数和对数函数。两种常用函数的形式如下。对数函数：P = alnQ + b （2）幂函数：P = cQd + e （3）式中：a
9、、b、c、d、e 均为常数；Q 取相对流量；P 取相对湿周（为消除坐标轴比例的影响，对湿周和流量无量纲化）。分别对式（2）和式（3）中的 Q 进行求导，即可得到两式的斜率方程表达式：对数函数：P = a/Q （4）幂函数：P = cdQ d -1 （5）再利用微积分中计算曲率的公式可得到两式的曲率方程表达式：对数函数：幂函数：斜率 1 法是取 PQ 关系曲线上斜率为 1 的一点，即 P = 1 时，对应的流量值作为河道内的最小生态流量。斜率 1 法求得的最小生态流量公式为：对数函数：Q = a （8）幂函数：Q = (cd )1/ (1-d ) （9）曲率法是选取曲线上曲率最大的一点对应的流量
10、值作为河道内的最小生态流量。为求得曲率最大的一点，需对式（6）和（7）中的 Q 求一阶导数，一阶导数值等于 0 的一点即为曲率最大的一点。该点对应的流量值就作为河道内的最小生态流量。曲率最大法求得的最小生态流量公式为：对数函数：Q = a/2（10）幂函数：2.2 Tennant 方法Tennant 法是 Tennant，D.L 等人于 1976 年提出的河流生态需水量评估方法6。该方法是一种历史流量法，按照 Tennant 方法的标准，河道内的最小生态流量不能小于多年平均流量的 10%。多年平均流量 1 0是防止河道退化的临界条件；多年平均流量的 20是保护水生栖息地的适宜标准；在较小河流中
11、，多年平均流量的 30为最佳栖息地标准。2.3 R2CROSS 方法R2CROSS 方法最初是由美国森林委员会提出的，其目的是为了美国高原山区河流的水资源与水环境保护，维持生物多样性而确定山区河流环境流量的方法。该方法首先根据研究河段控制断面的河顶宽度，查表 1 得到环境流量所需的水力学参数：平均水深、湿周率和平均流速。然后再根据该断面建立的水深、湿周率和平均流速与流量的关系分别得到 3 个流量 Q1、Q2 和 Q3，最后在 3 个流量中选出所需要的环境流量。方法是按季节要求：（1）如果是在夏季和秋季，那么平均水深、平均流速及湿周率必须全部满足，即生态流量为 Q1、Q2 和 Q3 中的最大值；
12、（2）如果是在冬季和春季，3 个水力参数满足两个即可，即为 Q1、Q2 和 Q33 个流量中的第二大值。表 1 R2CROSS 法确定生态流量的标准3 应用及存在的问题分析3.1 渭河关中段河道概况近十几年来，由于渭河上游过渡开发，维持河流基本生态功能所需的最小流量即生态基流得不到保证，一系列河流生态问题已经日益突显，为此，本研究以渭河关中段为研究对象。渭河关中段全长 504.2km，流域面积为 6.17104km2。渭河关中段有 5 个主要的水文站，分别为林家村水文站、魏家堡水文站、咸阳水文站、临潼水文站以及华县水文站（见图 1）。5 个水文站的多年平均流量分别为 84.8m3/s、117m
13、3/s、156m3/s、236m3/s 以及 251m3/s。渭河关中段河道的水面宽普遍较大，正常情况下都能达到 150300m 之间，而水深比较小，因此渭河关中段河道属于宽浅型河道。图 1 研究区域与水文站3.2 计算中的问题及分析讨论本文对建立湿周与流量关系中的问题以及 R2CROSS 法应用中出现问题进行讨论。（1）是否需要建立控制站点完整的湿周流量关系。由渭河关中段河道特征可知，河道非常宽，沟槽、滩地较多，很不规整，如下图 2 所示。是选择整个大断面来建立完整的湿周流量关系，还是以枯水期的主河槽来建立湿周与流量关系。图 2 咸阳站的大断面根据澳大利亚学者 Gippel 的研究4，湿周流
14、量关系的变化点（Breakpoint，该点对应的流量即为生态基础流量）常发生在河道断面有较大变化的地方。如图 2 所示咸阳站的大断面为典型的复式断面，若建立整个断面的湿周流量关系，则会出现变化点对应的基础流量非常大（图 2 中河道几何形态变化最大发生在左侧滩地处，河道水位要很高时才会出现），只有汛期漫滩时才能满足要求。因此，在建立计算生态基础流量的湿周流量关系时只能选择枯期有水的主河槽部分。（2）河道是否需要概化。在建立湿周与流量关系时，是直接根据实测大断面资料按均匀流计算湿周和流量的关系，还是需要将实测大断面绘出后将其概化成某一规则断面再计算湿周与断面的关系。为此，以林家村站为例（实测资料绘
15、制的图形如图 3 所示，类似抛物线）进行了相关计算分析。图 3 林家村水文站实测断面概化为抛物线。根据实测大断面资料（如图 3 所示），计算得到林家村水文站的抛物线方程为：y =1.8910-4x 2 （12）其中 y 为河底高程，x 为起点距。湿周按如下公式计算其中 a 为抛物线系数，对于林家村水文站 a=1.8910-4，c 为水面宽。根据公式（13）可计算出湿周与流量关系（如图 4 所示）。从图 4 可知湿周与流量之间的关系为对数关系。直接根据大断面实测资料计算。然而，根据实测大断面资料直接计算出的实际湿周与流量的关系则为幂函数关系，如图 5 所示。图 4 林家村站概化后计算得到的 P-
16、Q 关系图 5 根据林家村站实测资料计算出的 P-Q 关系从图 4 和图 5 对比可以发现：没有概化处理的湿周流量关系包涵的河道信息更完整、丰富，拟合出的关系更能代表实际情况，由此得到的生态基础流量应该更为合理。（3）R2CROSS 法应用中出现问题。从前面的介绍可知 R2CROSS 方法只适用于河顶宽度为 0.3m31m 的山区性河流，不适用于大中型河流。而渭河的河顶宽度都在 100m 以上，因此直接利用 R2CROSS 方法有较大的难度。3.3 方法改进因渭河关中段的河道宽阔，为典型的沟槽、滩多的复式河道。考虑到生态基础流量为河道健康和河道功能不退化的最小流量，因此，选择枯期有水的主河槽来
17、建立湿周流量关系。对选择的主河槽，根据实测大断面资料直接计算不同水位对应的湿周和流量，然后建立湿周流量关系。根据拟合好的线型，结合湿周法确定生态基础流量的方法计算生态基础流量。针对渭河河道特宽、多沙、缺水等特点，结合湿周法和 R2CROSS 法，考虑管理机构要根据不同季节河流健康要求、或不同河段中对水环境比较敏感的水生生物物种（包括鱼类或其它指示物种）在不同生长阶段对水动力学条件的要求而提出的一套方法综合法。在该方法中将生态基础流量对应的水深、流速和湿周率分别定义为生态水深、生态流速和生态湿周率。该方法的具体分析步骤如下。（1）根据研究河段控制点的大断面实测资料，分析计算出水深与流量，流速与流
18、量，湿周与流量等相关关系；（2）确定生态流速。可以根据河道特性从泥沙（河道泥沙问题突出的河流）或保护物种生长繁衍最适宜的流速考虑。渭河关中段河道的泥沙浓度很高并且主要为很细的悬沙，可根据实验或经验公式确定不冲不淤流速7：vm = h0 0.64（14）式中：h0 为平均水深（m）； 为与河道泥沙特性有关的系数，当沙比较粗时=0.600.70；当泥沙颗粒中等时 =0.540.57；当沙比较细时 =0.390.41。也可以根据研究河段保护的特征生物物种生长最适宜的流速选定。（3）确定生态水深。可以从控制站点枯水期多年平均水深的变化范围来确定。若有研究河段保护物种或指示物种生长繁殖对水深的要求，也可
19、以根据其要求来确定。（4）确定生态湿周率。该指标可以根据生态水深结合研究河段的大断面综合确定，一般选大断面的 50%左右为宜。（5）确定生态基础流量。根据确定的生态流速、生态水深和生态湿周率，在对应的水深与流量，流速与流量和湿周与流量关系曲线上查出 3 个对应流量，然后综合确定生态基础流量。4 计算结果与分析4.1 湿周法计算结果根据实测大断面资料，假定不同的水深，计算出各控制站点相应的水力学参数（如表2 所示）。表 2 计算渭河林家村站的水力参数根据表 2 建立各控制站的湿周和流量关系（如图 5 所示）。按照斜率 1 法和曲率最大法得到各控制点（水文站）的最小生态基础流量。渭河关中段各水文站
20、计算得到的最小生态流量成果见表 3 所示。表 3 渭河关中段各水文站最小生态需水量计算结果按照 Tennant 方法的标准，河道内的最小生态流量不能小于多年平均流量的 10%，应在 10%30%之间。渭河关中段的前 4 个水文站，即林家村水文站、魏家堡水文站、咸阳水文站和临潼水文站的最小生态流量占多年平均流量都在 10%20%之间，满足河道不退化的条件。华县水文站的最小生态流量占多年平均流量的 20%30%之间，在适宜值与最佳标准值之间。由此可见上述湿周法计算渭河关中段最小生态流量结果都满足最小生态需水要求，而且最小生态流量呈现出从上游到下游增大的合理趋势。综合斜率 1 法和曲率最大法计算结果
21、，再结合渭河关中地区水资源短缺的现实，湿周法确定渭河关中段各水文站的最小生态基础流量及其对应水力参数如表 4 所示。表 4 渭河关中段最小生态流量成果（湿周法）从表 4 中可看出，渭河关中段各水文站最小生态流量占多年平均流量的 10%20%时，平均水深在 0.20.7m 之间，断面流速在 0.40.65m/s 之间。4.2 综合法计算结果分别对渭河关中段的平均生态流速、平均生态水深、平均生态湿周率进行计算，计算结果如下。（1）平均生态流速。根据湿周法计算的成果，5 个站的生态基础流量对应的平均水深约为 0.50m，由公式（14）计算得到的不淤最小流速约为 0.25 m/s，又据水力学等8，河道
22、中的不冲不淤的允许流速范围为 0.210.32m/s。因此，确定渭河关中段的平均生态流速为 0.3m/s。（2）平均生态水深。由前面的湿周法算出的各站的生态流量值，可知生态流量所对应的平均水深变化范围较大，但最大值不会超过 1m，因此，确定平均水深的标准范围为 0.2 0.9m。（3）平均生态湿周率。由于渭河的断面较大，河顶宽度较长，湿周率定得太大后并没有足够的水量来满足要求，综合考虑取湿周率为研究主河道（河槽）湿周的 50%。根据以上标准，由各站建立的水深流量关系，流速流量关系和湿周率流量的关系，就可以得到相应的 3 个流量，最后再结合渭河关中段的实际情况确定各站的生态基础流量，最终结果如表
23、 5 所示。表 5 渭河关中段最小生态流量成果（综合法）从表 5 可知，利用综合法计算的生态基础流量：林家村水文站为 12.4m3/s，占多年平均流量 14.6 %；魏家堡水文站为 17.2m3/s，占多年平均流量 14.7 %，咸阳水文站为 23. 3m3/s，占多年平均流量 14.9%；临潼水文站为 25.5/s，占多年平均流量 10.8%；华县水文站为 30.1mm3/s，占多年平均流量 12.0 %。综合法计算的结果与湿周方法的计算结果在林家村、魏家堡、咸阳三站的结果基本一致，但是临潼和华县两站的结果比湿周法小。但是，从 5 个站的生态基础流量占多年平均流量的比值来看都在 10.8%1
24、4.9%之间。根据 Tennant 方法可知，综合法计算确定的生态基础流量能够保证河道功能不退化、水生生物的生存不会受到威胁，并且在水资源短缺的关中地区更能得到保障。5 结论针对湿周法和 R2CROSS 法在渭河关中段计算最小生态流量中遇到的问题提出了一套相对简单、实用的综合法。该方法以平均生态流速-平均生态水深和平均生态湿周率为基础进行生态基础流量计算，与其他方法对比该方法的优点主要表现在：（1）能够充分考虑维持河道特征及其形态稳定需要（从不冲不淤生态流速上反映）；（2）可以最大限度的保障生物栖息地信息（根据实测断面资料建立各种关系）；（3）方便、灵活便于管理部分调控：可以根据研究河段需保护
25、的水生生物物种的生活习性来确定生态基流。但是，该方法在应用也可能受到研究河段的基础资料、水生态和水生生物多样性调查资料等限制。计算结果表明：该方法计算结果比湿周法的结果更合理，适合在水资源缺乏地区推广应用。参考文献： 1 李丽娟，郑红星. 海滦河流域河流系统生态环境需水量计算J. 地理学报，2000，55（4）：496-500 . 2 宋进喜，李怀恩，王伯铎. 河流生态环境需水量研究综述J. 水土保持学报，2003，17（6）：95-97 . 3 吉利娜，刘苏峡，王新春. 湿周法估算河道内最小生态需水量-以滦河为例J. 地理科学进展，2010，29（3）：287-291 . 4 Gippel
26、C J，Stewardson M J . Use of wetted perimeter in defining minimum environmental flowsJ. Regul . River：Res . Mgmt . ，1998，14：53-67 . 5 Shang S H，A multiple criteria decision-making approach to estimate minimum environmental flows based on wetted perimeterJ. River Research and ApplicationsJ. 2008，24：54
27、-67 . 6 Tennant D L . Instream flow regimens for fish，wildlife，recreation and related environmental resourcesC/Proceedings，Symposium and Special Conference on Instream flow Needs，Osborn JF，Allman CH（eds）. Vol . II .May 3-6，Boise，Idaho，American Fisheries Society：Bethesda；1976 . 359-373 . 7 Men B H，Liu C M . Ecological hydraulic radius model to calculate instream flow requirements for transporting sediment in the western water transfer regionJ. Sci China Ser E-Tech Sci，2009，52（11）：3401-3405 . 8 吴持恭. 水力学M. 北京：高等教育出版社，2006 .作者简介：张新华（1965-），男，副教授，主要从事水力学与河流动力学及水环境研究。
