河流水质模型综合衰减系数确定的探讨.doc
- 文档编号:29189315
- 上传时间:2023-07-21
- 格式:DOC
- 页数:7
- 大小:137KB
河流水质模型综合衰减系数确定的探讨.doc
《河流水质模型综合衰减系数确定的探讨.doc》由会员分享,可在线阅读,更多相关《河流水质模型综合衰减系数确定的探讨.doc(7页珍藏版)》请在冰豆网上搜索。
河流水质模型综合衰减系数确定的探讨
李慧珑第一作者:
李慧珑,女,1933年生,教授,主要从事水文水资源专业的教学与科研工作。
诸晓华1包宝华2韩伯成2
(1.扬州大学水利科学与工程学院,江苏扬州225009;
2.江苏省水文水资源勘测局南京分局,江苏南京210008)
摘要河流水质模型综合衰减系数对计算水体的纳污能力、预测未来水质、制订污染物控制方案有着重要的影响。
采用实测水文数据、污染物初始浓度及设定的综合衰减系数代入水质模型,率定COD和TN的综合衰减系数,并进行合理性分析,作为同类河道应用时的参考。
关键词综合衰减系数纳污能力水质模型CODTN
Discussiononthedeterminationofcomprehensiveattenuationcoefficientinriverwaterqualitymodel
LiHuilong1,ZhuXiaohua1,BaoBaohua2,HanBocheng2.(1.HydraulicScienceandCivilEngineeringCollege,YangzhouUniversity,YangzhouJiangsu225009;2.NanjingBranchofJiangsuHydrologyandWaterResourceSurveysBureau,NanjingJiangsu210008)
Abstract:
Thecomprehensiveattenuationcoefficientinriverwaterqualitymodelhasanimportantimpactoncalculatingthepermissiblepollutionbearingcapacityofwaterbodies,predictingthefuturewaterquality,establishingthepollutantcontrolplan.Theauthorsusedthemeasuredhydrologicaldataandtheinitialconcentrationofpollutants,putthesupposedcomprehensiveattenuationcoefficientintowaterqualitymodel,educedagroupoftheCODandTNofcomprehensiveattenuationcoefficient,andhadareasonableanalysis.Itcanbeusedasareferenceforsimilarriver.
Keywords:
comprehensiveattenuationcoefficient;permissiblepollutionbearingcapacity;waterqualitymodel;COD;TN
有机污染物进入河道在随水流输移过程中,由于物理、化学与生物的作用污染物发生降解,其降解速度随河流的水文条件,如水深、流量、流速、水温、泥沙和污染物含量等因素而异,也与河道的形态有关。
在各类水质模型中,常将上述影响因素概化为综合衰减系数k,作为综合反映其导致的污染物降解强弱程度。
在实际计算水体纳污能力、预测污染物浓度、制定污染物排放控制方案工作中,综合衰减系数的合理性与准确性有着重要的影响。
然而,在不同河流不同时期影响因素千差万别、复杂多变,又使得k难以确定。
《水域纳污能力计算规程》(SL348—2006)提出了确定k的3种方法:
分析借用法、实测法和经验公式法[1]。
在实用中,分析借用法和经验公式法有时限于不同河段影响因素的差别,精度难以保证,而实测法要求的实验条件严格,较难操作。
笔者介绍了在进行南京市六合区水环境规划中,根据河流的水文情势选择水质模型,从实用性和可操作性出发,将与实测水位资料相应的河段流量、流速或蓄水量、水质污染物初始浓度及设定的k等代入水质模型,逐日模拟污染物浓度变化并与实测浓度值比较。
经过多次调整对k的取值,视模拟过程与实测值的拟合程度,选取拟合较好时的系数,再经过合理性分析后确定采用值。
1河流纳污能力计算模型的选择
(1)对于污染物在河段横断面上基本上能均匀混合的中、小型河段,采用一维恒定流水质模型以模拟污染物沿河段纵向的迁移,流经河段长L后的污染物浓度计算式[1,2]如下:
(1)
式中:
为河段L末端的污染物质量浓度,mg/L;为河段初始断面的污染物质量浓度,mg/L;为污染物综合衰减系数,1/s或1/d;u为过水断面平均流速,m/s。
若将沿程入河的点源与非点源污染物概化为自一个排污口入河,入河口位于河段L中间,则式
(1)成为:
(2)
式中:
为污染物入河速率,g/s或t/d;为河段流量,m3/s。
(2)对于水流基本上处于停滞状态的河段和中、小型水库、湖泊,可作为一个均匀混合的水体,其水质与污染物的入河位置无关,采用零维水质模型计算如下:
(3)
式中:
、为起始时刻与t时刻水体污染物质量浓度,mg/L;V为水体容积,m3。
2实例
滁河源于安徽,入江苏境曲折东流经南京市的浦口、六合区进入长江。
滁河干流六合段全长73.4km,其中的主要水功能区有3个,即孔湾至沿河段、沿河至铁路桥段和铁路桥至红山窑闸段,长度分别为16.7、3.5、12.5km。
前者为保留区,后两者分别为工业用水区和农业用水区。
3个功能区的河段上、下游有控制枢纽工程。
为保证区内灌溉、航运的需要,全年大部分时间关闸蓄水,仅在洪水期按工程控制调度预案有计划启闸泄洪。
关闸期间河水基本上不流动,水质较差,只有在短时间内抽、引长江水入滁河,或由滁河输水入灌区和支流径流汇入等才引起河道水位变化。
水功能区的水质目标已确定为Ⅳ类,但时有Ⅴ类及劣Ⅴ类的情况发生。
2.1污染物入河量年内分配
以滁河沿河至铁路桥段水功能区为例说明计算过程。
该区承受河段上游及区间支流的来水,随区间径流入河的非点源污染物主要来自水土流失、农田施肥、畜禽养殖、居民生活污水及城乡地面废弃物等,点源污染物则主要来自工业废水的排放。
经对2005年区间面积内各项污染物的COD及TN分别估算出年入河量后,再将其分配为月、日的入河量见表1、表2。
表12005年沿河至铁路桥水功能区COD入河量计算
月份
1月
2月
3月
4月
5月
6月
非点源COD量小计/t
140.89
140.89
140.89
245.59
245.59
245.59
点源COD入河量/t
252.36
252.36
252.36
252.36
252.36
252.36
COD合计入河量/t
393.25
393.25
393.25
497.95
497.95
497.95
COD日入河量/(t·d-1)
12.69
14.04
12.69
16.6
16.06
16.6
月份
7月
8月
9月
10月
11月
12月
非点源COD量小计/t
371.23
371.23
371.23
140.89
140.89
140.89
点源COD入河量/t
252.36
252.36
252.36
252.36
252.36
252.36
COD合计入河量/t
623.59
623.59
623.59
393.25
393.25
393.25
COD日入河量/(t·d-1)
20.12
20.12
20.79
12.69
13.11
12.69
表22005年沿河至铁路桥段COD计算表
日期
水位
/m
(m)
蓄水量
/万m3
(104m3)
COD
/(mg·L-1)
COD/t
备注
1月1日
6.86
432.46
20
86.49
水位为六合站实测值。
C0=20mg/L根据实际监测数据确定。
关闸期设k=0.07。
1月2日
6.86
432.46
20.76
89.76
1月3日
6.87
434.1
21.36
92.72
1月4日
6.87
434.1
21.98
95.4
…
…
…
…
…
8月6日
7.28
504.24
29.55
149
开闸期C0来自上一水功能区。
Q=276m3/s,u=0.267m/s。
设k=0.25。
8月7日
7.71
583.19
27.34
159.44
8月8日
7.96
631.64
25.37
160.25
8月9日
7.37
520.3
23.64
123
8月10日
7.32
511.35
22.11
113.06
关闸期设k=0.07。
8月11日
7.51
545.78
22.08
120.53
8月12日
7.45
534.79
23.8
127.29
…
…
…
…
…
9月2日
7.4
525.71
26.97
141.78
Q=209m3/s,u=0.194m/s。
设k=0.25。
9月3日
7.49
542.1
24.37
132.11
9月4日
8.77
801.44
22.18
177.76
9月5日
8.63
770.68
20.33
156.68
…
…
…
…
…
9月22日
6.95
447.39
10.86
48.59
9月23日
6.85
430.81
10.77
46.4
…
…
…
…
…
关闸期设k=0.07。
12月29日
6.56
384.41
31.51
121.12
12月30日
6.56
384.41
31.5
121.1
12月31日
6.58
387.53
31.24
121.08
2.2河段内污染物浓度变化模拟
根据《水域纳污能力计算规程》,滁河六合段属15m3/s<Q<150m3/s的中型河段,Q为多年平均流量。
(1)关闸期。
河段内水流基本上处于停滞状态,采用零维水质模型式(3)计算。
(2)开闸期。
采用一维恒定流式
(2)计算,其中河段流量Q以开闸期泄洪流量计。
按实际工况应用不同水质模型逐日计算沿河至铁路桥段COD见表2。
2005年8月6日至9日及9月2日至23日滁河段开闸泄洪,其余时间关闸。
将表3模拟的逐日COD与实测值绘于同一图中,视模拟与实测两者的拟合程度调整k,直至拟合较好时的k即为所取,见图1,k取0.07(关闸期),k取0.25(开闸期)。
同理通过实测TN与模拟的逐日过程拟合,见图2,得k为0.05与0.25。
采用同样方法可以得到孔湾至沿河及铁路桥至红山窑闸两河段的拟合过程,如图3至图6。
归纳三河段在开闸与关闸期的综合衰减系数见表3。
关闸期k=0.07,关闸期k=0.25
COD/(mg·L-1)
图1沿河-铁路桥段COD模拟与实测对照
关闸期k=0.05,关闸期k=0.25
TN/(mg·L-1)
图2沿河—铁路桥段TN模拟与实测对照
关闸期k=0.04,关闸期k=0.25
COD/(mg·L-1)
图3孔湾-沿河段COD模拟与实测对照
关闸期k=0.04,关闸期k=0.25
TN/(mg·L-1)
图4孔湾-沿河段TN模拟与实测对照
关闸期k=0.04,关闸期k=0.25
COD/(mg·L-1)
图5铁路桥-红山窑闸段COD模拟与实测对照
关闸期k=0.04,关闸期k=0.25
TN/(mg·L-1)
图6铁路桥-红山窑闸段TN模拟与实测对照
表3合区滁河干流COD、TN综合衰减系数k1/d
水功能区
COD
TN
关闸期
开闸期
关闸期
开闸期
孔湾至沿河段
0.04
0.25
0.04
0.25
沿河至六合铁路桥段
0.07
0.25
0.05
0.25
六合铁路桥至红山窑闸段
0.04
0.25
0.04
0.25
平均值(采用)
0.05
0.25
0.04
0.25
3合理性分析
图1~图6表明,模拟与实测值的拟合情况并不十分理想,尤其是TN的偏离更大。
这可能一方面是实际监测的数据较少,另一方面是污染物入河量的估算及年内分配、污染物入河在时空上的不确定性、k本身受各方面因素制约而不稳定以及水质监测误差等诸多原因所致。
为此,还需要对以上成果作合理性分析才能确定采用值。
(1)从表3可知,COD与TN在开闸期的k比关闸期大得多,这是由于关闸期间水体基本上不流动,污染物的降解主要依赖水体自身的物理、生化作用,缺少水流的迁移稀释;开闸后,河水具有一定的流速,与大气接触面大,曝气充分,有利于污染物的氧化和降解,其自净能力也比静水强得多。
(2)影响k的因素主要有:
水体的地形、微生物种类与数量、复氧能力、水温、流量、流速和污染物浓度等,这些因素在3个相邻河段中应大致相近,但沿河至六合铁路桥段流经城区,生活污水及工业废水的排放量都较大,该河段长度又较其余两段短得多,因而水体受污染的程度较其他两河段严重,从实际的监测数据也可得到证明。
据河南省水利厅郭金巨等在其省内13条河流所作的,有关COD、NH3-N与k关系的大规模实验表明:
当水体中的污染物较少时,水体的自净能力往往不能充分发挥,使k偏小,随着污染物浓度增大,在一定范围内k也随之增大至最大值,之后由于好氧反应减弱,k缓慢递减并趋于稳定。
由此认为沿河至铁路桥段的COD与TN的k略大于其他两段可能也是以上原因。
(3)中国环境规划院在《全国地表水水环境容量核定技术复核要点》(2004年)提出了水质降解系数参考值见表4。
表4水质降解系数参考值
水质及水生态环境状况
水质降解系数/d-1
COD
NH3-N
一般河道
湖泊水库
一般河道
湖泊水库
优(相应水质Ⅱ~Ⅲ类)
0.18~0.25
0.06~0.10
0.15~0.20
0.06~0.10
中(相应水质Ⅲ~Ⅳ类)
0.10~0.18
0.03~0.06
0.10~0.15
0.03~0.06
劣(相应水质Ⅴ类或劣Ⅴ类)
0.05~0.10
0.01~0.03
0.05~0.10
0.01~0.03
本次分析的河段水体水质状况属于Ⅳ~Ⅴ类,在关闸时的水文情势近于湖泊水库,所得综合衰减系数k在上表推荐的参考值范围内,说明成果具有一定的合理性。
4结语
本文以滁河实测水质数据率定模型参数,将2005年实测滁河六合站水位资料、水质污染物初始浓度及设定的综合衰减系数代入水质模型,逐日模拟污染物浓度变化并与实测浓度值比较。
经过多次对综合衰减系数的不同取值,取拟合程度较好的综合衰减系数为采用值。
所得综合衰减系数k在相关技术文献推荐的参考值范围内,说明成果具有一定的合理性,可以在滁河六合段及其支流采用,也可作为同类河道应用时的参考。
本法对监测断面水质资料依赖性较强,当实测资料有偏差时,会导致拟合结果的偏差而失真,因此应用本法进行河流综合衰减系数估值时,应充分重视水质断面的监测和资料的收集,要增加监测频率,提高监测精度。
参考文献
[1]SL348-2006,水域纳污能力计算规程[S].
[2]史晓新,朱党生,张建永.现代水资源保护规划[M].北京:
化学工业出版社,2005.
责任编辑:
贺锋萍(收到修改稿日期:
2008-04-10)
©版权所有《环境污染与防治》杂志社
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 河流 水质 模型 综合 衰减系数 确定 探讨