余真真，闫 莉，王瑞玲，马秀梅，张世坤

(黄河水资源保护科学研究院，河南 郑州 450004)

以往关于自净需水的研究[1]仅从量的方面考虑，一般采用河流多年平均天然径流量的10%及近10年最枯月平均流量或90%保证率最枯月平均流量，或以历史水文数据为基础的7Q10法、Tennant法[2]等计算得到，大都缺乏对研究河段水体功能和生态保护目标的识别。在黄河流域相关研究中，郝伏勤等[3]研究了黄河干流2003年现状纳污、污染可控及目标控制3种纳污模式下的宁蒙、龙三、小花河段的自净水量；黄锦辉等[4]在河流水质及污染源调查的基础上，分析河流水质和污染源分布及排放规律之间的关系，建立了湟水、渭河、沁河、伊洛河自净需水模型，提出保障水功能区水质目标所需的水量；张建军等[5]建立了黄河干流自净需水量模型，分别对黄河花园口以上河段以及利津断面的自净需水量进行了计算。以上研究偏重以水功能区水质目标为约束条件，较少考虑生态系统功能发挥对水质的需求，以及生态系统变化和流量、水质的响应关系。本文建立基于河道边界条件水动力学与水质相耦合模型，结合伊洛河生态系统特点及功能性需水组成，以河流水功能区水质目标为协控因子，进一步考虑伊洛河河口黄河鲤繁殖期对水质的要求进行自净需水流量计算，同时考虑不同河流纳污水平，将入河排污口作为分散点源，实现水质、水量及污染源同步输入情景模拟。

1 伊洛河生态环境状况

1.1 生态环境特征

受水沙条件、水体物理化学性质及流域气候、地理条件等因素影响，黄河水生生物种类和数量相对贫乏，生物量较低，鱼类种类相对较少，但许多特有土著鱼类具有重要保护价值，是国家水生生物保护和鱼类物种资源保护的重要组成部分。伊洛河是小浪底水利枢纽以下黄河一级支流，流域面积 1.89万km2，流经陕西、河南两省。随着小浪底水库“调水调沙”运行，黄河下游频繁出现“流鱼”现象，同时水库低温水下泄也对鱼类资源破坏极大。黄河下游支流水环境相对稳定，生物行为学特性研究表明，鱼类在不适宜生境条件下，可以依靠其侧线管等感应系统的判断，转移到环境相对适宜的水域。目前黄河鱼类有自黄河下游干流迁徙到连通支流的趋势，因此对作为鱼类主要生境的伊洛河的保护和修复显得尤为重要。

伊洛河流域位于我国暖温带和北亚热带的分界区，体现出南北过渡、东西交汇的特征，流域上游是自然植被主要分布区，中下游地区人类活动和干扰比较严重，以农田植被为主。上游水质优良，为Ⅱ类，河流经过城镇后，水质变差，为Ⅲ类或Ⅳ类。干流灵口以上河段水环境相对良好，其他河段均存在水电无序开发现象，河流纵向连通性遭到严重破坏，枯水期及平水期多个河段存在断流、脱流现象，水生生物合适的生境数量锐减、面积萎缩且片段化。

1.2 水生生物及生态习性

黄河水生生态系统中，鱼类处于水生生物群落食物链的顶层，对其他类群的存在和丰度有着重要作用。伊洛河流域曾经是我国著名经济鱼类“洛鲤伊鲂”的出产地，是多种地方土著鱼类的重要分布区，共有鱼类4目8科36种，其中鲤科鱼类25种，占总数的69%；鳅科鱼类4种，占总数的13%；鮠科、鲿科、鲶科、鰕虎鱼科、合鳃鱼科、攀鲈科、鳢科鱼类各1种，各占3%[4]。鲤科鱼类是伊洛河标志性鱼种群，具有重要遗传与生态保护价值。

伊洛河入黄口河段由于水量充足、营养丰富，一直是流域的鱼类传统产卵场，具有较佳的栖息生境，目前被划为黄河鲤国家级水产种质资源保护区，其中伊洛河核心区位于巩义市境内，保证黄河鲤繁殖生长期的水量是该河段水生态保护的主要目标。黄河鲤属鲤形目鲤科鲤亚科，是黄河流域著名经济鱼种之一，以洛河和黄河中下游最为盛产，洛河内鲤鱼又称“洛鲤”，伊洛河口成为维持黄河鲤渔业资源的主要集中产卵场，也是黄河中下游其他鱼类繁殖、育幼阶段最为重要的栖息场所之一。鲤鱼属分批产卵鱼类，繁殖期延续时间较长，黄河鲤产卵时间为4—6月，要求有较为适宜的产卵条件(依据繁殖期黄河鲤栖息地适宜度曲线，黄河鲤产卵期适宜流速为0.15～0.4 m/s，适宜水深为0.8～1.5 m，并对水质有较高要求)。

2 标志性鱼种群繁育对水质的要求

水污染是黄河鱼类保护面临的主要威胁因子。自1980年起，由于经济的快速发展，各种工业废水及生活污水大量排放，黄河流域水质总体呈现恶化的趋势[6]。2000年前后是全流域水质最差的时段，满足Ⅲ类水的河长比例仅占38.7%[6]。黄河鱼类资源自20世纪80年代以来大幅度减少，鱼类种质资源数量减少约42%，鱼类资源及种类减少除受不合理捕捞、大坝阻隔、气候变化等因素影响外，河川径流条件的变化和水质恶化是导致黄河鱼类资源减少的重要原因。

水质直接影响着黄河水生生物的繁殖及栖息，水体污染造成鱼类生长发育滞缓、生殖能力减弱；同时污染造成水体富营养化，引起水中藻类多样性减少，优势度提高，鱼类的可利用食料减少，影响到鱼类的多样性和数量。根据野外监测结果和实践经验，地表水Ⅲ类的标准可以满足和保证鱼类繁育，黄河水域的Ⅲ类水质对于鱼类胚胎发育的影响几乎接近阈值，水质恶化将影响鱼类产卵和孵化。GB 11607—1989《渔业水质标准》对鱼类水质要求是：淡水pH值为6.5～8.5，溶解氧质量浓度连续 24 h 中16 h以上必须大于5 mg/L(对应Ⅲ类标准)。GB 3838—2002《地表水环境质量标准》依据地表水水域环境功能和保护目标，规定Ⅱ类水适用于珍稀水生生物栖息地、鱼虾类产卵场、仔稚幼鱼的索饵场等，Ⅲ类水适用于鱼虾类越冬场、洄游通道、水产养殖区等渔业水域及游泳区。

伊洛河入黄口河段分布有国家级水产种质资源保护区，对水质和水量过程的要求较高。笔者认为维持伊洛河口水产种质资源保护区核心区的黄河鲤正常繁育所需的水质为Ⅲ类，维持鱼类良好产卵生存状态下所需的水质为Ⅱ类。

3 自净需水流量核算模型构建

3.1 控制方程及求解方法

建立机理性水质模型，根据质量守恒定律、能量守恒定律与动量守恒定律，应用模型方法建立数学方程，通过实测资料确定相关参数，对水环境中污染物在时间和空间上的分布做定量的分析和预测。

采用数值模拟方法研究伊洛河流域自净需水量。圣维南方程作为求解控制方程，由以下连续方程和运动方程联立组成：

(1)

式中：x、t分别为计算点空间坐标和时间坐标；A为过水断面面积；Q为过流流量；h为水位；q为旁侧入流流量或排污口流量；c为谢才系数；R为水力半径；α为动量校正系数；g为重力加速度。

此外，水质模型的控制方程为

(2)

式中：C为水质浓度；ux为纵向流速；Ex为扩散系数。

方程组利用Abbott-Ionescu六点隐式有限差分格式求解。该格式在每一个网格点按顺序交替计算水位和流量，分别称为h点和Q点。离散后的线性方程组用追赶法求解[7-8]。

3.2 计算河段选取

研究范围为伊洛河流域及小浪底以下黄河干流河段。小浪底以下黄河干流河段自小浪底至入海河口，全长895.7 km。伊洛河为伊河陆浑水库坝址断面至洛河交汇处，洛河为故县水库坝址断面至入黄口。伊洛河流域自净需水流量研究范围与河网构建详见图1和图2。

3.3 入河排污口调查与设置

伊洛河流域污染源达标排放后黄河不得不接纳的污染物量及其时空分布的量化标准，是计算伊洛河流域自净需水流量的前提[9]。为准确了解伊洛河纳污现状，调查流域入河排污口数据，发现伊洛河共接纳废污水2.8亿m3/a，COD 1.8万t/a，氨氮0.1万t/a。其中洛河纳污量占3/4，伊河占1/4。伊洛河的纳污特点是：纳污量80%以上集中在洛阳、偃师及巩义河段，超标因子为COD、BOD5和氨氮等，主要来源于生活排污。因此，伊洛河水质沿程变化趋势是河流经过城镇后，受人为排污影响水质变差。将调查到的37个入河排污口以点源形式设置于模型中，给定排污流量及浓度值。

图2 伊洛河流域自净需水流量模型河网构建Fig.2 River network building of self-purificationwater demand flow model in the Yiluo River

4 结果与分析

4.1 水功能区环境目标条件下的自净需水流量

分别以现状排污与控制排污两种情景模式，推算主要控制断面满足水质目标所需的自净需水流量，见表1。现状排污条件下，伊洛河入黄断面黑石关要确保满足汛期11.50 m3/s、非汛期 14.00 m3/s 的自净需水流量，才能保证入黄水质满足Ⅳ类水质目标要求，即ρ(COD)≤30 mg/L，ρ(NH3-N)≤1.5 mg/L。控制排污条件下(入河情景模拟所有排污口满足一级A达标排放标准，该条件下，经计算COD、氨氮入河控制量分别小于1.87万t、1 055 t，满足限制排污总量控制要求[10])，伊洛河入黄断面黑石关要确保满足汛期6.6 m3/s、非汛期8.0 m3/s的自净需水流量，才能保证入黄水质满足水质目标要求。

表1 现状排污与控制排污条件下自净需水流量Table 1 Self-purification water demand flow undercurrent and control of sewage conditions m3/s

4.2 生态功能性需水条件下的自净需水流量

4.2.1基本维持鱼类正常繁育的自净需水流量

a. 模拟条件设置。现状洛河故县、长水和伊河陆浑水库断面水质满足Ⅱ类水质标准，背景浓度选取伊河、洛河Ⅱ类水质目标上限浓度，推算入河排污口达标排放浓度下，黑石关能够基本维持黄河鲤4—6月正常繁育期间，水质目标达到Ⅲ类(ρ(COD)= 20 mg/L、ρ(NH3-N)=1 mg/L)所需的自净需水流量。

b. 预测结果。经模拟预测，在达标纳污水平条件下，洛河故县水库坝址下泄断面COD、氨氮自净需水流量分别为4.62 m3/s、2.26 m3/s；伊河陆浑水库坝址下泄断面COD、氨氮自净需水流量分别为 3.08 m3/s、1.52 m3/s；入黄口黑石关断面COD、氨氮自净需水流量为13.62 m3/s、9.60 m3/s。按照偏严格确定，在控制排污条件下，伊洛河入黄断面黑石关4—6月应满足13.62 m3/s的自净需水流量，才能基本维持鱼类正常繁育需求所需的水质要求。

4.2.2维持鱼类良好生存状态下的自净需水流量

a. 模拟条件设置。背景浓度选取伊河、洛河Ⅱ类水质目标上限浓度，推算入河排污口达标排放浓度下黑石关断面维持黄河鲤等代表鱼类良好生存状态，4—6月正常繁育期间，水质目标达到Ⅱ类(ρ(COD)=15 mg/L、ρ(NH3-N)=0.5 mg/L)所需的自净需水流量。

b. 预测结果。经模拟预测，在达标纳污水平条件下，洛河故县水库坝址下泄断面COD、氨氮自净需水流量为10.20 m3/s、18.00 m3/s；伊河陆浑水库坝址下泄断面COD、氨氮自净需水流量为6.80 m3/s、12.00 m3/s；入黄口黑石关断面COD、氨氮自净需水流量为22.70 m3/s、35.80 m3/s。按照偏严格确定，在控制排污条件下，伊洛河入黄断面黑石关4—6月的自净需水流量应满足35.80 m3/s，这是能够维持鱼类繁殖期良好生存状态下的自净需水流量。

4.3 自净需水流量满足程度与可实现性分析

伊洛河现状入河排污口排污、满足限排总量控制要求下，黑石关断面达到水功能区水质目标Ⅳ类所需的自净需水流量分别为14.0 m3/s、8.0 m3/s，占黑石关多年平均流量的10.70%、6.18%。从考虑伊洛河水生态功能角度出发，维持伊洛河口黄河鲤等土著代表鱼类栖息生境对优良水质的需求，满足4—6月基本维持鱼类正常繁育水质要求的自净需水流量为13.62 m3/s，维持鱼类良好繁育状态下水质要求的自净需水流量为35.80 m3/s，分别占黑石关多年平均流量的10.79%、28.37%。

根据黑石关2009—2018年长系列实测流量数据，开展伊洛河自净需水流量满足程度与可实现性分析，结果发现，与多年实测流量相比，伊洛河现状入河排污状态下满足水功能目标所需自净需水流量满足程度为88.89%，控制排污状态下自净需水流量满足程度为96.80%。2009—2018年每年1月1日满足水功能水质目标条件下自净需水流量与实测流量对比见图3。

图3 满足水功能水质目标条件下自净需水流量与实测流量对比Fig.3 The comparison between self-purification andmeasured water flow under the condition ofmeeting water quality objectives

提取2009—2018年4—6月黑石关断面实测流量数据进行对比分析，发现基本维持鱼类正常繁育水质要求的自净需水流量满足程度为94.11%，维持鱼类良好繁育状态下水质要求的自净需水流量满足程度为61.28%(图4)。因此，2009—2018年的黑石关断面实际流量满足鱼类良好繁育状态所需的满足程度不高(特别是4月上旬至5月下旬春耕用水挤占了部分鱼类产卵期满足良好水质要求的水量)，凸显了人类用水与生态竞争性用水之间的矛盾。

图4 满足鱼类繁育水质要求条件下4—6月自净需水流量与实测流量对比Fig.4 The comparison between self-purification andmeasured water flow under the condition of meetingthe requirements of fish breeding water qualityfrom April to June

5 结 语

今后应加强重要生态单元的水质监测，提高水质生物监测能力，逐步改善生态保护重点河段水环境质量；实施流域多目标生态调度，在确保防汛安全的前提下，将伊洛河河道内满足鱼类生存的自净水量需求纳入黄河水资源统一调度体系，尽可能提高关键期黄河土著代表鱼类繁育期对水质的需求；把水资源作为最大的刚性约束，协调人类用水与生态竞争性用水关系，完善伊洛河流域用水总量控制指标，加大农业节水力度，推动用水方式从粗放向节约集约转变。