安阳河动态纳污能力分析

来源 :人民黄河 | 被引量 : 0次 | 上传用户:estone
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
  摘 要:在不同保证率和不同时间尺度的设计水文条件下,动态计算出安阳河市区段2016年水体纳污能力和年内分配,并对计算结果进行了对比分析,结果表明:安阳河市区段整体处于污染状态,3个水功能区中仅安阳河安阳市排污控制区污染物排放不超标,其余水功能区COD和NH3-N排放量均超过纳污能力上限,且NH3-N超标情况更为严重;安阳河市区段汛期月均纳污能力与非汛期月均纳污能力相差较小;安阳河市区段纳污能力1—5月整体呈下降趋势,5—7月呈逐渐上升趋势,至8月纳污能力达到最大,随后又呈现波动下降趋势,5月为全年最小;纳污能力的大小与时间尺度成负相关关系,不同时间尺度下确定的年纳污能力月尺度>水期尺度>年尺度。
  关键词:纳污能力;保证率;时间尺度;安阳河
  中图分类号:X522;TV212.5 文献标志码:A
  doi:10.3969/j.issn.1000-1379.2021.07.021
  引用格式:耿文斌,赵姗,侯俊山,等.安阳河动态纳污能力分析[J].人民黄河,2021,43(7):110-113.
  Abstract:Under the designed hydrological conditions of different guarantee rates and different time scales, the paper dynamically calculated the sewage capacity and annual distribution of the urban section of Anyang River in 2016. The calculation results were comparatively analyzed. The results show that the pollutant discharge in Anyang sewage control area is not exceeding the standard, the COD and NH3-N discharge of other water function areas exceed the upper limit of the capacity and the NH3-N is more serious. The monthly average water environmental capacity of the river in flood season is smaller than that in non-flood season, the water environmental capacity of the river in July to October is larger, the water environmental capacity of August is the maximum in the whole year and the water environmental capacity of November to May is gradually reduced, in May, the water environmental capacity is the lowest in the whole year. There is a negative correlation between the size of the capacity and the time scale and the annual capacity is the monthly scale > the water scale > the annual scale.
  Key words: water environmental capacity; guarantee rate; time scale; Anyang River
  納污能力是指一定条件下水体的最大纳污量[1-2]。纳污能力是制定污染物限排量的重要依据,能直接反映水体的水质状况[3-4]。目前我国纳污能力大多按照《水域纳污能力计算规程》(GB/T 25173—2010)计算,其中设计水文条件选取90%保证率最枯月平均流量或近10 a最枯月平均流量,这种条件下得到的纳污能力为一固定值,且设计水文条件较为严格,以此作为标准应用于水资源管理和保护必然安全[5]。但实际纳污能力与水文过程关系密切,水文过程在自然条件下是动态变化的,纳污能力也是一个动态变化的量[6-7]。将一个确定的值作为控制指标难以揭示水体的实际污染负荷,不利于保护和改善水质,也不利于水资源的合理开发利用,因此开展水体动态纳污能力评价很有必要[8]。
  我国对动态纳污能力的研究众多,但多以不同保证率、不同水文条件来计算水体动态纳污能力。李乔臻[9]以95%、90%、75%、50%、20%保证率下的流量为设计流量,选取不同来水水质分析了琶江河段3个水功能区的纳污能力;周洋等[10]运用一维稳态水质模型和水环境容量模型,根据水质、水文监测资料和沿岸排污情况,计算了渭河陕西段丰水年、平水年、枯水年COD的水环境容量。这些研究都只分析了不同保证率下年纳污能力的动态变化特征,并未结合年内水文特性进行分析,只能确定水功能区全年的污染物限排量,年内的污染物限排量分配仍不清晰,以此为水资源管理的依据会出现全年污染物排放量不超标,但某些月份因污染物排放量过多而造成水功能区水质变差的现象。笔者基于安阳水文站48 a(1969—2016年)长系列资料,在不同保证率条件下,根据不同时间尺度确定设计水文条件,以COD和NH3-N为控制因子,建立河流一维水质模型[11-13],计算2016年安阳河市区段年、水期和月纳污能力,以期为安阳河确定污染物限排总量、年内污染物限排量分配提供数据支撑,为安阳河管理与保护提供科学依据。   1 研究区概况
  安阳河是海河流域漳卫南运河水系的第二大支流,发源于林州市清泉寺,自西向东流经林州市、安阳县、安阳市区、内黄县(见图1),在内黄县李大晁入卫河,全长约164 km,流域面积1 920 km2。上游穿过林州盆地后,在横水镇附近进入山区,于彰武水库以下进入平原。干流上建有彰武水库和小南海水库。安阳河的主要补给水源包括岩溶泉水(泉主要有小南海泉和珍珠泉,泉域面积935 km2)和大气降水,据安阳水文站多年观测资料统计,最大流量可达2 060 m3/s(1982年8月2日),最小流量为0.29 m3/s(1982年7月8日),多年平均流量为10.73 m3/s[14]。
  安阳河市区段共包括3个二级水功能区,分别为安阳河彰武水库饮用水源区、安阳河安阳县蒋村用水区和安阳河安阳市排污控制区,总长约49.1 km,水质目标分别为Ⅱ类、Ⅲ类和Ⅳ类。
  2 纳污能力计算模型
  2.1 计算模型
  本次计算单元为彰武水库大坝至入卫河口,共包含3个水功能区,计算单元示意见图2。
  安阳河总体上具有流速小、河道宽、水浅等特点,主要河道宽深比小,入河污染物在较短河段内混合均匀,污染物横向变化不大。因此,可采用概化一维水质模型进行纳污能力计算,计算公式如下:
  M=31.536(QP+Q)(CseKx86.4u-C0e-KL86.4u)(1)
  式中:M为水功能区纳污能力,t/a;Q为水功能区设计流量,m3/s;QP为废污水排放量,m3/s;Cs为水功能区水质控制目标质量浓度,mg/L;C0为水功能区初始断面质量浓度,mg/L;K为污染物综合衰减系数,1/d;x为概化点距控制断面的距离,m;L为水功能区长度,km;u为设计流量下断面平均流速,m/s。
  根据研究区排污口分布以及资料掌握情况[15],将各个水功能区内多个排污口概化为一个集中排污口,且集中排污口位于河段中点处,因此纳污能力计算公式改为
  M=31.536(QP+Q)(CseKL2×86.4u-C0e-KL86.4u)(2)
  2.2 模型参数
  2.2.1 设计水文条件的确定
  水体纳污能力动态变化受多因素共同影响,而设计水文条件是纳污能力计算的重要前提,对计算结果影响较大。基于安阳水文站多年监测数据,分别确定不同时间尺度和不同保证率下的设计水文条件,并对纳污能力的变化特征进行分析。
  (1)基于年尺度的设计水文条件的确定。选取系列资料中逐年最枯月平均流量,采用P-Ⅲ型曲线参数估计法绘制最佳经验频率曲线,确定50%、75%、90%保证率下的设计流量分别为2.94、1.76、1.03 m3/s,根据河道大断面状况分析得到对应的设计流速。
  (2)基于水期尺度的设计水文条件的确定。将年内月份分为汛期(6—10月)和非汛期(11月至翌年5月),选取逐年汛期和非汛期内的最枯月平均流量,采用P-Ⅲ型曲线参数估计法,分别绘制汛期和非汛期的最佳经验频率曲线,确定75%保证率下的设计流量分别为2.06、1.88 m3/s,90%保证率下的分别为1.16、1.07 m3/s,并根据河道大断面状况分析得到对应的设计流速。
  (3)基于月尺度的设计水文条件的确定。选取逐年各月平均流量,采用P-Ⅲ型曲线参数估计法分别绘制1—12月的最佳经验频率曲线,确定90%保证率下1—12月的设计流量分别为2.13、2.08、1.88、1.52、1.37、1.48、1.98、2.46、2.24、2.33、2.16、2.12 m3/s,再根据河道大断面状况分析得到对应的设计流速。
  根据经验频率曲线分析结果可知,相同保证率下,时间尺度越小,对应的设计流量越大。主要原因是设计水文条件的确定方法不同,年尺度下选择逐年最枯月平均流量进行计算,水期尺度则分别选择逐年汛期最枯月平均流量和非汛期最枯月平均流量进行计算,而安阳河属于常年河,枯水年和平水年没有明显的汛期特征,甚至个别年份汛期内月平均流量小于非汛期内月平均流量,因此汛期和非汛期确定的设计流量更大。月尺度分别选择对应的月平均流量作为系列值进行分析,确定的设计流量更大。
  2.2.2 水质浓度及综合衰减系数的确定
  水功能区初始断面质量浓度C0采用实测数据;水功能区水质目标质量浓度Cs以水质类别为控制指标。根据《全国地表水环境容量核定技术复核要点》中提供的河道水质综合衰减系数参考值,结合各水功能区优劣状况和水文特征选取综合衰减系数。由经验法,COD的综合衰减系数计算公式为
  K=0.05+0.68u(3)
  式中:u为设计流速。
  NH3-N的综合衰减系数计算公式为
  K=0.061+0.551u(4)
  2.3 数据来源
  安阳河水文数据来源于安阳水文站长系列监测数据;污染物排放数据来源于安阳市水利局和安阳水文水资源勘测局共同出版的《安阳市入河排污口监测报告》,水功能区污染物入河量统计见表1。
  3 纳污能力计算分析
  3.1 年尺度纳污能力计算分析
  安阳河市区段年尺度納污能力计算结果见表2,由表2可以看出,安阳河市区段3个水功能区中安阳河安阳市排污控制区纳污能力最大,安阳河彰武水库饮用水源区纳污能力最小;当上游来水水质优于水功能区水质目标,其他因素一定时,纳污能力与保证率的大小成负相关关系,即纳污能力随保证率的增大而变小。
  3.2 水期尺度纳污能力计算分析
  安阳河市区段水期尺度纳污能力计算结果见表3,由表3可以看出,安阳河市区段汛期月均纳污能力与非汛期月均纳污能力相差较小。主要原因是安阳河为常年河,在丰水年,径流量较为集中,主要集中在6—10月,占年径流量的60%~80%,平水年、枯水年年内径流分配较为均匀,汛期特征不明显,特枯水年基本无明显的汛期特征,且汛期安阳河上游进行灌溉、水库蓄水等,导致汛期月均流量与非汛期月均流量相差较小。   3.3 月尺度纳污能力计算分析
  安阳河市区段月尺度纳污能力计算结果见表4,安阳河市区段COD和NH3-N年纳污能力分别为1 393.4 t和91.5 t。由表4可知,安阳河市区段逐月纳污能力呈波动变化趋势,并未与流量季节变化保持完全一致,说明安阳河纳污能力不仅受自然水文条件的影响,还受到人类活动的影响,且人类活动对纳污能力的影响较大,如6月为汛期,但受人类活动影响(上游灌溉、水库蓄水等),其纳污能力较5月增大不明显。安阳河市区段纳污能力1—5月整体呈下降趋势,5—7月呈逐渐上升趋势,至8月纳污能力达到最大,随后又呈现波动下降趋势。
  3.4 综合对比分析
  由不同尺度(年、水期、月)纳污能力的对比分析可知,纳污能力与时间尺度成负相关关系,不同时间尺度下确定的年纳污能力月尺度>水期尺度>年尺度。
  主要原因是,年尺度纳污能力计算分析时考虑多年的不利水文条件,以此为全年的水文条件进行纳污能力计算,这种计算方法虽然对水环境的改善有促进效果,但计算尺度较大,确定的设计水文条件不具备代表性,水体水质状况大多数时段优于设计水文条件,确定的年纳污能力较小。随着计算尺度的减小,设计水文条件以水期、月为尺度确定,更接近水体实际情况,计算得到的纳污能力更接近水体污染承载力。
  由确定的纳污能力与年污染物入河量进行对比分析,安阳河市区段整体处于污染状态,3个水功能区仅安阳河安阳市排污控制区污染物排放不超标,其余水功能区COD和NH3-N排放量均超过纳污能力上限,且NH3-N超标更严重。主要原因是,水域周围农业生产及养殖活动较多,废水未经处理直接排入河道,因此研究区在控制污染物排放时,对NH3-N排放的限制应更严格。
  4 结 语
  纳污能力与水体水文过程关系密切,是一个动态变化的量,笔者在不同保证率和不同时间尺度下对安阳河市区段动态纳污能力进行了计算,分析了不同保证率和不同时间尺度对纳污能力的影响。考虑保证率和时间尺度的变化,可以使纳污能力的计算结果与水体实际状况更相符。目前水功能区存在监测时间间隔过长的问题,且以特定的监测点水质分析结果来判断整个水功能区的水质状况,存在极大的偶然性,因此建议以污染物排放量与纳污能力比值判断水功能区是否达标,并分别从月、水期和年尺度进行纳污能力考核,以避免传统考核模式存在的偶然性问题。本文仅进行了年、水期和月尺度的纳污能力动态分析,今后可考虑更小时间尺度(日、时)的计算,以了解水体的实际状况,及时对水资源管控做出调整,提高水功能区纳污能力。
  参考文献:
  [1] 李如忠,汪家权,王超,等.不确定性信息下的河流纳污能力计算初探[J].水科学进展,2003,14(4):359-363.
  [2] 赵鑫,黄茁,李青云.我国现行水域纳污能力计算方法的思考[J].中国水利,2012(1):29-32.
  [3] 王林,宓辰羲,宓永宁.辽阳市水功能区纳污能力计算研究[J].人民黄河,2017,39(12):80-84.
  [4] 张秀菊,杨凯,蔡爱芳,等.不确定性参数对水体纳污能力的影响分析[J].中国农村水利水电,2012(1):13-17.
  [5] 中华人民共和国质量监督检验检疫总局,中国国家标准化委员会.水域纳污能力计算规程:GB/T 25173—2010[S].北京:中国水利水电出版社,2011:1-4.
  [6] 张海欧,闵涛,罗军刚,等.渭河干流纳污能力动态计算模型系统开发[J].人民黄河,2013,35(6):70-72.
  [7] 张晓,罗军刚,杜金龙,等.水功能区纳污能力动态计算主题模式研究及应用[J].水利信息化,2016(2):24-28.
  [8] 史晓新,禹雪中,马巍.湖泊纳污能力动态特征分析及计算[J].中国水利水电科学研究院学报,2008,6(2):105-110.
  [9] 李乔臻.琶江河段动态纳污能力分析[J].水文,2009,29(增刊1):131-133.
  [10] 周洋,周孝德,冯民权.渭河陕西段水环境容量研究[J].西安理工大学学报,2011,27(1):7-11.
  [11] 张川,吕谋.河道纳污模型研究与应用[J].青岛理工大学学报,2015,36(1):74-79,113.
  [12] 韓龙喜,朱党生,蒋莉华.中小型河道纳污能力计算方法研究[J].河海大学学报(自然科学版),2002,30(1):35-38.
  [13] 路雨,苏保林.河流纳污能力计算方法比较[J].水资源保护,2011,27(4):5-9,47.
  [14] 范兰池,石政华,苗利芳.安阳河流域水文特性分析[J].海河水利,2009(6):23-24.
  [15] 赵嵩林,杨扬,王国涛,等.安阳市入河排污口监测报告(2016年)[R].安阳市水利局,2017:5-8.
  【责任编辑 吕艳梅】
其他文献
摘 要:浅水方程由纳维-斯托克斯方程推导而来,是一种描述具有自由表面的浅水体在重力作用下流动的数学模型。实际浅水流动不可能是完全理想的情况,必须考虑底坡和摩阻等源项的存在,它们影响计算的稳定性和精确性,如果处理不当还会造成格式的不和谐,因此源项处理是求解浅水流动方程的关键,也是国内外学者研究的一个重要方向。介绍目前求解浅水方程存在的4个主要困难,介绍浅水方程的各种离散方法以及优缺点,阐述源项处理的
摘 要:采用黄河流域九省(区)2004—2018年统计资料,从水资源、社会、经济、生态环境4个方面构建水资源承载力评价指标体系,运用熵权法和层次分析法获得各指标权重,利用综合评价模型评价黄河流域水资源承载力,并利用障碍度函数定量诊断影响黄河流域水资源承载力提升的主要障碍因素。研究结果表明:2004—2018年,黄河流域水资源承载力呈现缓慢上升趋势,整体向良好方向发展,空间上呈现明显的波动特征且差异
摘 要:为揭示卫宁灌区地下水埋深对农业节水灌溉的响应规律,根据灌区2002—2017年地下水埋深、灌溉水量和种植结构资料,采用线性倾向估计法、累计距平法和灰色关联度法进行研究。结果表明:近年来受农业节水的影响,卫宁灌区地下水埋深呈增大趋势,2012年开始地下水埋深增大速率变大;受季节和农业灌溉的影响,地下水埋深年内呈周期性变化,分布呈W形,年内5—10月地下水埋深增幅最大,原因是该时期灌溉水量减少
摘 要:为分析气候变化背景下黄河流域蒸发量变化特征及其影响因素,采用线性回归、Mann-Kendall突变检验分析和经验正交分解等方法对1961—2019年黄河流域95个气象观测站蒸发皿蒸发量变化进行分析。结果表明:近59 a黄河流域蒸发量呈显著下降趋势,全流域平均气候倾向率为-60.4 mm/10 a;2004年以后,蒸发量下降趋势加剧(气候倾向率为-219.4 mm/10 a);夏季下降趋势相
摘 要:气候变化和人类活动加剧使得水文序列的一致性遭到破坏,在变化环境下传统频率分析方法不再适用。选择沁河下游武陟站1951—2016年共66 a实测年径流量数据为研究对象,使用Hurst系数法对年径流量序列进行初步诊断,再运用线性回归、Spearman秩次相关检验和Kendall秩次相关检验3种方法对序列进行趋势诊断,采用Bayesian变点分析、Pettitt检验等9种方法进行跳跃诊断,基于水
摘 要:黄河下游的洪水灾害历来是中华民族的心腹之患,黄河决溢、河道迁移等事件对下游人民的生命财产安全构成了极大的威胁。根据历史文献记载及对前人研究成果的梳理,探讨了明清时期河工建设与黄河洪水事件之间的关系。研究表明:明清时期河工建设明显响应于黄河洪水事件,河工与水患呈现正相关性,明代后期黄河水患的发生频率高于前期,河工建设也呈现相应趋势;明代前中期,河南地区水患频繁,河工建设也较多;明代后期和清代
快速准确地获取水体信息对于水资源管理利用以及灾害防治具有重要意义。利用不同分类方法探讨多种数据源在不同天气场景下水体提取的最优技术,结合Sentinel和Landsat系列数据,对红碱淖1973—2018年湖泊面积变化进行分析。结果表明:在无云的情况下,使用Sentinel-2结合最大似然分类方法提取效果最好,精度为99.30%;Sentinel-1利用面向对象分类精度最高,总体精度为95.70%
摘要:长江马鞍山河段主流摆动频繁,滩槽演变规律复杂,汊道分流比调整幅度较大,局部崩岸频发,影响防洪安全、河势稳定、航道尺度提升以及沿江涉水工程的正常运行,对长江经济带发展和黄金水道建设造成影响,不利于长江大保护工作开展,有必要对该河段演变规律开展进一步研究。基于前人研究成果,采用实测资料归纳总结了该河段河道演变规律,分析了河道整治的重点问题。研究结果表明:马鞍山河段历史上经历了洲滩合并靠岸、汊道减
摘 要:为探究黄河流域生态保护和高质量发展的研究热点及前沿趋势,以1992—2021年CNKI数据库为数据来源,采用Citespace文献计量软件,结合样本文献时间、期刊和学科分布统计图,进行作者及機构合作网络分析、关键词分析、关键词LLR聚类和Burst分析,并得出相应的知识图谱。结果显示:黄河流域生态保护和高质量发展的研究与政策颁布密切相关;研究热点及前沿趋势可划分为黄河流域生态环境、水资源污