河流最小生态基础流量分析计算

河流最小生态基本流量计算措施研究（张新华 李红霞 肖玉成 赵少华）摘要：为了可以计算几何断面复杂、泥沙淤积、河宽敞、水深浅的河道生态基本流量，本文以渭河关中地区的生态基本流量计算为例，通过综合水力学中的湿周法和R2CROSS法建立了一套简朴、适合管理的最小生态基本流量措施综合法。该措施计算出的5个站的生态基本流量占近年平均流量的比值都在10.8%14.9%之间，满足生态基流的合理规定。通过与Tennant措施比较，成果表白：综合法计算拟定的生态基本流量更为合理、是一种在管理上更值得推广应用的措施。核心词：生态需水量；湿周法；斜率法；曲率最大法；R2CROSS法中图分类号：X171文献标记码：A文章编号：1672-3031（）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 profiles，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 indicate 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：ecological water demand；wetted perimeter method；slope method；curvature method and R2CROSS1研究背景河流系统涉及河流、湖泊及其相邻洪泛区。河流系统功能一般涉及生态功能、环境功能和资源功能1。然而，在开发水资源的时候重要考虑到满足工业、农业和城乡生活用水等需求，忽视了河流系统功能得以正常发挥并维持其基本功能时其自身对水量的规定，从而导致河段内的天然流量逐渐减少、甚至浮现断流的状况。此外，经济社会的高速发展导致大量的污水排放，使得有些河流虽然水量上没问题，但水质污染十分严重导致水生生物种群及其数量急剧下降，水体严重富营养化，浮现藻类或水华爆发等现象。为了维持河流系统健康，保障其基本功能能正常发挥，国内外学者为此开展了大量的研究工作，重要集中在河流生态环境需水上并获得了一定的研究成果。河流生态环境需水量计算内容重要涉及：（1）生态基本流量；（2）输沙需水量；（3）自净需水量；（4）蒸发需水量；（5）防海水入侵的河道最小水量等2。其中的河流生态基本流量是维持河道功能不退化、特性水生生物物种能正常繁衍生息所必须维持的河道最小流量，本文重点分析讨论这一河道最小生态基本流量。然而，国内外针对这一流量的研究尚未形成统一结识，基本理论和措施都尚有待进一步完善。随着近几年的研究进一步，国内外学者提出了多种计算措施，概括起来重要涉及4类措施：（1）水文学法；（2）水力学法；（3）栖息地法；（4）地形构造法等。在这几种措施中，以水文学法与水力学法最为常用。因此，本研究措施的重点在分析讨论水力学法中的湿周法并结合R2CROSS法和Tennat法，提出一套简朴、实用、便于管理部门操作的措施综合法。2计算措施2.1湿周措施湿周法计算河道的最小生态需水量时要基于一种假设，即浅滩是临界的河流栖息地，而保护浅滩栖息地也能保护到其她的水生生物栖息地。通过河道的实测大断面和实测流量资料建立湿周流量关系曲线，在曲线中用斜率1法和曲率最大法确立变化点，该点相应的流量即为最小生态流量，这就是湿周法确立河道最小生态流量的基本原理3-5。根据水力学中计算明渠均匀流的谢齐公式：Q = ACRJ以及曼宁公式：C = 1/nR1/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、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）曲率法是选用曲线上曲率最大的一点相应的流量值作为河道内的最小生态流量。为求得曲率最大的一点，需对式（6）和（7）中的Q求一阶导数，一阶导数值等于0的一点即为曲率最大的一点。该点相应的流量值就作为河道内的最小生态流量。曲率最大法求得的最小生态流量公式为：对数函数：Q = a/2（10）幂函数：2.2 Tennant措施Tennant法是Tennant，D.L等人于1976年提出的河流生态需水量评估措施6。该措施是一种历史流量法，按照Tennant措施的原则，河道内的最小生态流量不能不不小于近年平均流量的10%。近年平均流量1 0是避免河道退化的临界条件；近年平均流量的20是保护水生栖息地的合适原则；在较小河流中，近年平均流量的30为最佳栖息地原则。2.3 R2CROSS措施R2CROSS措施最初是由美国森林委员会提出的，其目的是为了美国高原山区河流的水资源与水环保，维持生物多样性而拟定山区河流环境流量的措施。该措施一方面根据研究河段控制断面的河顶宽度，查表1得到环境流量所需的水力学参数：平均水深、湿周率和平均流速。然后再根据该断面建立的水深、湿周率和平均流速与流量的关系分别得到3个流量Q1、Q2和Q3，最后在3个流量中选出所需要的环境流量。措施是按季节规定：（1）如果是在夏季和秋季，那么平均水深、平均流速及湿周率必须所有满足，即生态流量为Q1、Q2和Q3中的最大值；（2）如果是在冬季和春季，3个水力参数满足两个即可，即为Q1、Q2和Q33个流量中的第二大值。表1 R2CROSS法拟定生态流量的原则3应用及存在的问题分析3.1渭河关中段河道概况近十几年来，由于渭河上游过渡开发，维持河流基本生态功能所需的最小流量即生态基流得不到保证，一系列河流生态问题已经日益突显，为此，本研究以渭河关中段为研究对象。渭河关中段全长504.2km，流域面积为6.17104km2。渭河关中段有5个重要的水文站，分别为林家村水文站、魏家堡水文站、咸阳水文站、临潼水文站以及华县水文站（见图1）。5个水文站的近年平均流量分别为84.8m3/s、117m3/s、156m3/s、236m3/s以及251m3/s。渭河关中段河道的水面宽普遍较大，正常状况下都能达到150300m之间，而水深比较小，因此渭河关中段河道属于宽浅型河道。图1研究区域与水文站3.2计算中的问题及分析讨论本文对建立湿周与流量关系中的问题以及R2CROSS法应用中浮现问题进行讨论。（1）与否需要建立控制站点完整的湿周流量关系。由渭河关中段河道特性可知，河道非常宽，沟槽、滩地较多，很不规整，如下图2所示。是选择整个大断面来建立完整的湿周流量关系，还是以枯水期的主河槽来建立湿周与流量关系。图2咸阳站的大断面根据澳大利亚学者Gippel的研究4，湿周流量关系的变化点（Breakpoint，该点相应的流量即为生态基本流量）常发生在河道断面有较大变化的地方。如图2所示咸阳站的大断面为典型的复式断面，若建立整个断面的湿周流量关系，则会浮现变化点相应的基本流量非常大（图2中河道几何形态变化最大发生在左侧滩地处，河道水位要很高时才会浮现），只有汛期漫滩时才干满足规定。因此，在建立计算生态基本流量的湿周流量关系时只能选择枯期有水的主河槽部分。（2）河道与否需要概化。在建立湿周与流量关系时，是直接根据实测大断面资料按均匀流计算湿周和流量的关系，还是需要将实测大断面绘出后将其概化成某一规则断面再计算湿周与断面的关系。为此，以林家村站为例（实测资料绘制的图形如图3所示，类似抛物线）进行了有关计算分析。图3 林家村水文站实测断面概化为抛物线。根据实测大断面资料（如图3所示），计算得到林家村水文站的抛物线方程为：y =1.8910-4x2（12）其中y为河底高程，x为起点距。湿周按如下公式计算其中a为抛物线系数，对于林家村水文站a=1.8910-4，c为水面宽。根据公式（13）可计算出湿周与流量关系（如图4所示）。从图4可知湿周与流量之间的关系为对数关系。直接根据大断面实测资料计算。然而，根据实测大断面资料直接计算出的实际湿周与流量的关系则为幂函数关系，如图5所示。图4林家村站概化后计算得到的P-Q关系图5根据林家村站实测资料计算出的P-Q关系从图4和图5对比可以发现：没有概化解决的湿周流量关系包涵的河道信息更完整、丰富，拟合出的关系更能代表实际状况，由此得到的生态基本流量应当更为合理。（3）R2CROSS法应用中浮现问题。从前面的简介可知R2CROSS措施只合用于河顶宽度为0.3m31m的山区性河流，不合用于大中型河流。而渭河的河顶宽度都在100m以上，因此直接运用R2CROSS措施有较大的难度。3.3措施改善因渭河关中段的河道宽阔，为典型的沟槽、滩多的复式河道。考虑到生态基本流量为河道健康和河道功能不退化的最小流量，因此，选择枯期有水的主河槽来建立湿周流量关系。对选择的主河槽，根据实测大断面资料直接计算不同水位相应的湿周和流量，然后建立湿周流量关系。根据拟合好的线型，结合湿周法拟定生态基本流量的措施计算生态基本流量。针对渭河河道特宽、多沙、缺水等特点，结合湿周法和R2CROSS法，考虑管理机构要根据不同季节河流健康规定、或不同河段中对水环境比较敏感的水生生物物种（涉及鱼类或其他批示物种）在不同生长阶段对水动力学条件的规定而提出的一套措施综合法。在该措施中将生态基本流量相应的水深、流速和湿周率分别定义为生态水深、生态流速和生态湿周率。该措施的具体分析环节如下。（1）根据研究河段控制点的大断面实测资料，分析计算出水深与流量，流速与流量，湿周与流量等有关关系；（2）拟定生态流速。可以根据河道特性从泥沙（河道泥沙问题突出的河流）或保护物种生长繁衍最合适的流速考虑。渭河关中段河道的泥沙浓度很高并且重要为很细的悬沙，可根据实验或经验公式拟定不冲不淤流速7：vm = h00.64（14）式中：h0为平均水深（m）；为与河道泥沙特性有关的系数，当沙比较粗时=0.600.70；当泥沙颗粒中档时=0.540.57；当沙比较细时=0.390.41。也可以根据研究河段保护的特性生物物种生长最合适的流速选定。（3）拟定生态水深。可以从控制站点枯水期近年平均水深的变化范畴来拟定。若有研究河段保护物种或批示物种生长繁殖对水深的规定，也可以根据其规定来拟定。（4）拟定生态湿周率。该指标可以根据生态水深结合研究河段的大断面综合拟定，一般选大断面的50%左右为宜。（5）拟定生态基本流量。根据拟定的生态流速、生态水深和生态湿周率，在相应的水深与流量，流速与流量和湿周与流量关系曲线上查出3个相应流量，然后综合拟定生态基本流量。4计算成果与分析4.1湿周法计算成果根据实测大断面资料，假定不同的水深，计算出各控制站点相应的水力学参数（如表2所示）。表2计算渭河林家村站的水力参数根据表2建立各控制站的湿周和流量关系（如图5所示）。按照斜率1法和曲率最大法得到各控制点（水文站）的最小生态基本流量。渭河关中段各水文站计算得到的最小生态流量成果见表3所示。表3渭河关中段各水文站最小生态需水量计算成果按照Tennant措施的原则，河道内的最小生态流量不能不不小于近年平均流量的10%，应在10%30%之间。渭河关中段的前4个水文站，即林家村水文站、魏家堡水文站、咸阳水文站和临潼水文站的最小生态流量占近年平均流量都在10%20%之间，满足河道不退化的条件。华县水文站的最小生态流量占近年平均流量的20%30%之间，在合适值与最佳原则值之间。由此可见上述湿周法计算渭河关中段最小生态流量成果都满足最小生态需水规定，并且最小生态流量呈现出从上游到下游增大的合理趋势。综合斜率1法和曲率最大法计算成果，再结合渭河关中地区水资源短缺的现实，湿周法拟定渭河关中段各水文站的最小生态基本流量及其相应水力参数如表4所示。表4渭河关中段最小生态流量成果（湿周法）从表4中可看出，渭河关中段各水文站最小生态流量占近年平均流量的10%20%时，平均水深在0.20.7m之间，断面流速在0.40.65m/s之间。4.2综合法计算成果分别对渭河关中段的平均生态流速、平均生态水深、平均生态湿周率进行计算，计算成果如下。（1）平均生态流速。根据湿周法计算的成果，5个站的生态基本流量相应的平均水深约为0.50m，由公式（14）计算得到的不淤最小流速约为0.25 m/s，又据水力学等8，河道中的不冲不淤的容许流速范畴为0.210.32m/s。因此，拟定渭河关中段的平均生态流速为0.3m/s。（2）平均生态水深。由前面的湿周法算出的各站的生态流量值，可知生态流量所相应的平均水深变化范畴较大，但最大值不会超过1m，因此，拟定平均水深的原则范畴为0.2 0.9m。（3）平均生态湿周率。由于渭河的断面较大，河顶宽度较长，湿周率定得太大后并没有足够的水量来满足规定，综合考虑取湿周率为研究主河道（河槽）湿周的50%。根据以上原则，由各站建立的水深流量关系，流速流量关系和湿周率流量的关系，就可以得到相应的3个流量，最后再结合渭河关中段的实际状况拟定各站的生态基本流量，最后成果如表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%14.9%之间。根据Tennant措施可知，综合法计算拟定的生态基本流量可以保证河道功能不退化、水生生物的生存不会受到威胁，并且在水资源短缺的关中地区更能得到保障。5结论针对湿周法和R2CROSS法在渭河关中段计算最小生态流量中遇到的问题提出了一套相对简朴、实用的综合法。该措施以平均生态流速-平均生态水深和平均生态湿周率为基本进行生态基本流量计算，与其她措施对比该措施的长处重要表目前：（1）可以充足考虑维持河道特性及其形态稳定需要（从不冲不淤生态流速上反映）；（2）可以最大限度的保障生物栖息地信息（根据实测断面资料建立多种关系）；（3）以便、灵活便于管理部分调控：可以根据研究河段需保护的水生生物物种的生活习性来拟定生态基流。但是，该措施在应用也也许受到研究河段的基本资料、水生态和水生生物多样性调查资料等限制。计算成果表白：该措施计算成果比湿周法的成果更合理，适合在水资源缺少地区推广应用。参照文献：1李丽娟，郑红星.海滦河流域河流系统生态环境需水量计算J.地理学报，55（4）：496-500 .2宋进喜，李怀恩，王伯铎.河流生态环境需水量研究综述J.水土保持学报，17（6）：95-97 .3吉利娜，刘苏峡，王新春.湿周法估算河道内最小生态需水量-以滦河为例J.地理科学进展，29（3）：287-291 .4Gippel C J，Stewardson M J . Use of wetted perimeter in defining minimum environmental flowsJ. Regul . River：Res . Mgmt .，1998，14：53-67 .5Shang S H，A multiple criteria decision-making approach to estimate minimum environmental flows based on wetted perimeterJ. River Research and ApplicationsJ. ，24：54-67 .6Tennant 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 .7Men 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，52（11）：3401-3405 .8吴持恭.水力学M.北京：高等教育出版社， .作者简介：张新华（1965-），男，副专家，重要从事水力学与河流动力学及水环境研究。