1. 引言
全球可供人类利用的淡水资源已严重不足,且水质不断恶化,水污染造成的水体功能丧失进一步加重水资源危机的形势,全世界每年排入河流和湖泊的废水使全球水资源总量14%以上的水源受到不同程度的污染。严重的水质污染问题,不仅使水资源无法使用,也使农产品受到污染,影响了人民健康和农民的收入。因此,我国水污染问题己经成为不亚于洪灾、旱灾甚至更为严重的灾害。合理地利用环境资源,控制和减少环境污染已成为各国经济及其发展应当考虑的主要问题之一。迅速有效地控制水环境污染,对保护水资源质量,以水资源的可持续开发利用支持社会经济的可持续发展具有极其重要的意义。2009年我国提出以水资源配置、节约和保护为重点,实行更为严格的水资源管理制度,其中“明确水功能区限制纳污红线,严格控制入河排污总量”是三条“红线”之一。随着人口的增加和经济的发展,人类生存的环境日益受到破坏,尤其是水环境破坏倍受世人关注,纳污能力理论研究成为当前的研究热点之一。
渭河流域是陕西省重要的工农业科研和生产基地,人口多,然而渭河干、支流水质日益加剧恶化,严重影响着城市居民生产、生活用水,水污染治理是当务之急。为保证进入黄河的水质,从环保目标和管理需求出发,分析预测渭河流域各河流的水域纳污能力和进入河流污染物控制研究,对渭河水污染控制、水环境管理与水资源保护规划具有重要的意义 [1]。因此,对渭河流域水环境纳污能力计算研究和入河污染物控制研究能够为渭河流域综合治理提供依据,为关中地区提供有限的水资源,缓解水资源紧缺的现状,促进经济发展。
水域纳污量的概念多是以水质管理和水环境承载力的方式提出。美国环境保护局于1972年最先提出最大日负荷总量(total maximum daily loads, TMDL)概念及最大年负荷总量(total maximum yearly loads, TMYL)概念。污染负荷量可以表示为单位时间的质量、毒性和其他适合测定的指标。TMDL计划总的目标就是用来识别具体污染区和土地利用状况,并且考虑对这些具体区域点源和非点源污染物浓度和数量提出控制措施,从而引导整个流域执行最好的流域管理计划 [2]。TMDL计划在美国广泛实施,在点源和非点源综合控制方面成效显著 [3]。1983年12月8日,美国正式立法,实施以水质限制为基点的排污控制路线。
水域纳污量或水环境容量的计算方法是水域纳污量研究的一个重要组成部分,其适当与否直接影响计算结果的准确性。但由于我国地域广大,水体特性分异明显,加深了水域纳污量计算方法研究的难度。水域纳污量计算的研究方法主要有解析法 [4]、模型试错法 [5]、系统分析法 [6] 和概率稀释模型法 [7] [8] [9] [10]。制定一个能够使用于不同水体水环境容量计算的理论体系,从而推动水资源保护工作的深化,是一个重要课题。
因此,本文主要研究以渭河干流陕西河段,收集各种资料,开展入河排污口调查和水质监测,根据渭河水体水功能区划和水质保护目标,研究渭河开发利用的各个二级水功能区的纳污能力,由于水功能区水体纳污能力受多种因素的影响,尤其是水文设计条件对水功能区的纳污能力的影响至关重要,计算不同保证率下的允许纳污量,为水功能区限制纳污红线管理及水质监督管理提供技术支持。
2. 研究区域
渭河流域关中地区(图1)地处陕西中部,西安、宝鸡、咸阳、渭南等大中城市均坐落于内,是陕西省政治、经济、文化、金融及信息中心。但是由于人口增加、经济发展、用水剧增和排污量不断增大,生态环境恶化、污染加重的双重压力,渭河水环境严重超载。
根据黄河流域有关水利部《全国水功能区划技术大纲》中提供的方法和环保部门先后制定的有关黄河流域的水环境功能区划,区划结果按5类水环境功能区执行《地表水环境质量标准》(GB3838-2002)水质标准,可得到渭河干流陕西段的水功能区划,见表1。
渭河干流重点水功能区(陕西段)概化图,如图1。
![](Images/Table_Tmp.jpg)
Table 1. Water function zoning of the Shaanxi section of the main stream of the Weihe River
表1. 渭河干流陕西段水功能区划
![](//html.hanspub.org/file/15-1150783x9_hanspub.png?20210928113224262)
Figure 1. Generalized map of key water function areas (Shaanxi section) of the main stream of the Weihe River
图1. 渭河干流重点水功能区(陕西段)概化图
3. 纳污能力集散模型及参数的确定方法
3.1. 污染源概化
在河流纳污能力计算中,污染源的概化方式对河段纳污能力有较大影响,应尽量使污染源的概化符合实际情况,以准确反映功能区河段在现有污染源分布方式下的纳污能力。实际操作时,需对排入功能区河段的污染源进行必要的调查统计,对排污量大或集中排放、排入位置固定的点污染源按其实际情况进行概化,对排污量较小且分布较散的点污染源和面污染源采用均匀分布的概化方式,以此计算的河段纳污能力能较准确地反映实际情况。
3.2. 纳污能力计算模型
适用于污染物在横断面上均匀混合的中小型河段(Q ≤ 150 m3/s),其计算模型如下:
1) 河段的污染物浓度按式(1)计算:
(1)
式中:
Cx:流经x距离后的污染物浓度,单位为毫克每升(mg/L);
X:沿河段的纵向距离,单位为米(m);
U:设计流量下河道断面的平均流速,单位为米每秒(m/s);
K:污染物综合衰减系数,单位为负一次方秒(1/s);
其余符号意义同前。
2) 相应的水域纳污能力按式(2)计算:
(2)
式中符号意义同前。
3.3. 参数确定
3.3.1. 流速确定
本课题利用各水文站实测流量、流速资料,建立流量–流速关系曲线,分析计算渭河干流各水文条件下90%、75%、50%的保证率设计流量对应的流速。
根据已有的资料,拟合后的流量(Q)流速(u)关系为:
1) 林家村
(0 < Q < 50) (3)
2) 咸阳
(0 < Q < 2000) (4)
3) 华县
(0 < Q < 2000) (5)
可以通过这样的相关关系由已知的设计流量去推求流速。对于其他断面可根据附近区域的流量流速关系估计设计流速,以此,可计算出各断面不同设计水文条件下设计流量相对应的流速见下表所示表2。
![](Images/Table_Tmp.jpg)
Table 2. Design flow velocity of each section in the driest month (unit: m/s)
表2. 最枯月各断面的设计流速(单位:m/s)
3.3.2. 初始断面污染物浓度值Co的确定
目前渭河干流污染严重,除上游部分河段外基本上不存在河流水体的天然状态,而且由于河流水体的连续传递性,自下河沿断面以下河段,上游河段的污染均在一定程度上影响到下游水体的水质,若Co值仍沿用通常取值法显然是不尽合理的,因此在进行水域纳污能力审定时,为体现“上游河段污染不影响下游河段”,初始断面污染物浓度Co选用上游计算河段的水功能区水质目标,这既能够强化上游水质较好河段的水资源保护,又对下游河段是公平的,容易被各计算河段或水功能区所接受 [11]。
水质目标值是根据黄河流域总体规划目标要求,参考《西安市城市水资源保护规划报告》、《陕西省渭河流域综合治理五年规划(2008年~2012年)》,依据《全国水资源保护规划技术大纲》和黄河流域水污染特点,采用CODcr和氨氮作为污染物控制参数。CODcr和氨氮标准采用《地表水环境质量标准》(GB3838-2002),标准值见表3。
![](Images/Table_Tmp.jpg)
Table 3. Limits of surface water environmental quality standards (unit: mg/L)
表3. 地表水环境质量标准限值(单位:mg/L)
3.3.3. 降解系数K的确定
污染物降解系数反映了污染物在水中的降解速率,是计算纳污能力的关键参数。参数能否准确确定,直接影响到纳污能力的计算结果。综合上述方法,渭河干流COD和NH3-N的K值分别为0.43/d和0.30/d [12]。
4. 纳污能力计算结果
4.1. 纳污能力计算步骤
根据渭河干流陕西段地表水各功能区的划分与水质目标值,并结合2008陕西省水文水资源勘测局提供的渭河流域入河排污口排污资料,以及1951年到2010年的水文资料,进行本论文针对不同保证率下对于COD和NH3-N的纳污能力计算。
本课题利用国家标准模型(一维模型)进行计算,具体计算步骤如下:
步骤1:按式(1)计算第一个水功能区的污染物浓度Cx;
步骤2:按式(2)计算第一个水功能区的纳污能力M1;
步骤3:重复步骤1,2,算出所有断面的纳污能力并求和,即为渭河干流陕西段的纳污量。
4.2. 纳污能力计算结果
根据前两章计算的最枯月设计流量和设计流速,以及确定的其他参数,利用上面的步骤以此计算最枯月的COD和NH3-N河流的纳污量。结果见表4,在计算中出现了负值的现象,是与Cs的选择有关,由于本次毕设选择本功能区的水质类别下限值作为目标浓度Cs,下一断面的目标浓度值Cs小于上一断面的目标浓度值Cs,利用国标模型计算公式计算才会出现这种情况,这不是计算错误导致的,它不影响本文的结论,本次毕设是由计算的结果来反映总体的变化的。从结果中,也可以看出对于COD的纳污量会高于NH3-N的纳污量。
![](Images/Table_Tmp.jpg)
Table 4. Calculation of pollutant holding capacity in each functional area based on the driest month (unit: t/a)
表4. 基于最枯月各功能区的纳污量计算(单位:t/a)
5. 结论
本文在不同保证率下对于COD和NH3-N的纳污能力进行了计算,在年最枯月的90%条件下计算的纳污能力,在年内的纳污能力是一个定值,计算结果比较严格,应用于管理也比较安全。本论文根据我国2009年提出的水资源管理制度中三条红线之一——“明确水功能区限制纳污红线,严格控制入河排污总量”,结合各水功能区的水质管理目标,计算了年最枯月的90%条件下计算的纳污能力,从结果可以看出,对渭河干流陕西段污染严重的功能区必须采取一定的措施来有效治理,防止水质的进一步恶化。以便缓解水资源短缺的状况。
水域的纳污能力是排污总量控制的基础,其定量评价对于有效地保护水资源具有重要的现实意义。根据纳污能力制定科学的改善策略,具有一定的前瞻性和可操作性,不仅可以改善水环境,还可以节约人力、物力和财力,实现人与自然和谐相处和社会经济的可持续发展。
基金项目
长安大学省重点实验室开放基金(2019-JC08);陕西省土地工程建设集团内部项目(DJNY2021-25)。