里山取水口水质达饮用水水质的方案探讨

里山取水口水质达饮用水水质的方案探讨

杭州的自来水水源从池塘中提取85%，并将于每年7月至11月受到洪水的影响。如果河水中含有25010-6，则不能用作供水水源。随着城市的发展,用水量日益增加。目前,作为主要水源的珊瑚沙水库及其取水口的含氯度指标只有80%的保证率,将不能满足日益增长的用水需求。为提高供水保证率,计划将取水口再移至上游12km处的里山(图1)。据调查,钱塘江两岸均设有排污口,在规划年(2010年)里山位置的水质能否达到饮用水水质的标准是公众关心的问题。本文通过收集水文水质资料,借助Mike21水质模型,预测规划期里山位置的水质情况,并对上游径流量、上游来水水质、下游潮汐条件、污染负荷的改变等情况作敏感性分析。1枯水期总磷指标水质监测资料包括为取水工程专门所测的资料和常规测验站的资料。从收集的资料来看,里山河段枯水期NH3-N为II-III类,劣于平水期(总体II类);枯水期总磷总体为II类,优于平水期(II-III类);CODMn和CODCr在2个水期总体基本相当,为I-II类;平水期总氮指标较高,为V-劣V类。分析原因可能是平水期下泄流量大,水环境稀释扩散作用优于枯水期,因此平水期NH3-N水质优于枯水期,而平水期总氮和总磷指标较高,可能是径流量增大,带入农业面源增加所造成。2最大连带以6102q3.2xq为中心的内部水网络拓扑结构采用丹麦水力学研究所研制的平面二维数值模型Mike21研究工程区域的潮流和污染物对流扩散运动。该模型计算软件具有算法可靠、计算稳定、界面友好、前后处理方便等突出优点,并在70多个国家得到应用。潮流模型包括1个连续性方程和2个动量方程。基本方程为:∂ζ∂t+∂p∂x+∂q∂y=QsAs;(1)∂ζ∂t+∂p∂x+∂q∂y=QsAs;(1)∂p∂t+∂∂x(p2h)+∂∂y(pqh)+gh∂ζ∂x-Ωq+gp√p2+q2C′2h2=E(∂2p∂x2+∂2p∂y2);(2)∂p∂t+∂∂x(p2h)+∂∂y(pqh)+gh∂ζ∂x−Ωq+gpp2+q2√C′2h2=E(∂2p∂x2+∂2p∂y2);(2)∂q∂t+∂∂x(pqh)+∂∂y(q2h)+gh∂ζ∂y+Ωp+gq√p2+q2C′2h2=E(∂2q∂x2+∂2q∂y2)。(3)∂q∂t+∂∂x(pqh)+∂∂y(q2h)+gh∂ζ∂y+Ωp+gqp2+q2√C′2h2=E(∂2q∂x2+∂2q∂y2)。(3)对流扩散方程:∂∂t(hC)+∂∂x(uhC)+∂∂y(vhC)=∂∂x(hDx∂C∂x)+∂∂y(hDy∂C∂y)+Qs⋅CsAs。(4)∂∂t(hC)+∂∂x(uhC)+∂∂y(vhC)=∂∂x(hDx∂C∂x)+∂∂y(hDy∂C∂y)+Qs⋅CsAs。(4)式中:ζ为潮位,即水面到某一基准面的距离;p,q分别为x,y方向上的垂线平均流量分量;u,v分别为x,y方向上的垂线平均流速分量;h为水深;g为重力加速度;Ω为柯氏力参数;C′为谢才系数;E为涡动粘性系数;t为时间;C为沿程平均污染物浓度;Dx,Dy分别为x,y方向上的扩散系数;Qs为污水排放流量;Cs为污水排放浓度;As为排放口位置的流体微团面积,在数值计算中取x,y向网格宽度的乘积。上述方程组的初始条件:ζ(x,y)|t=0=ζ0(x,y);p(x,y)|t=0=p0(x,y);q(x,y)|t=0=q0(x,y);C(x,y)|t=0=C0(x,y)。ζ(x,y)|t=0=ζ0(x,y);p(x,y)|t=0=p0(x,y);q(x,y)|t=0=q0(x,y);C(x,y)|t=0=C0(x,y)。边界条件:(1)水边界。ξ(x,y,t)=ζ*(x,y,t),“*”表示已知值;C(x,y,t)=C*(x,y,t),流入计算域;∂(hC)∂t+∂(uhC)∂x+∂(vhC)∂y=0∂(hC)∂t+∂(uhC)∂x+∂(vhC)∂y=0,流出计算域。(2)陆边界。Qn=0,法线方向流量为零;∂C∂n=0∂C∂n=0。根据上述定解条件,Mike21软件采用交替方向隐格式(ADI)求解二维浅水潮波方程,方程矩阵采用双向消除(DoubleSweep)算法求解。3水质验证结果二维上边界取在上游窄溪,下边界选在钱塘江闸口站。模型采用矩形网格,空间步长为50m×100m,时间步长5s,江道地形采用2004年实测地形。枯水期验证采用2004年1月的资料,平水期验证采用2004年9月的资料。现状条件下,模型范围内污水总量为95万t/d,CODCr排放量为95t/d,富阳渔山江段NH3-N排放量为1.73t/d。计算中,CODCr与CODMn的换算值取2.5。各监测断面的水质验证结果见表1。通过水流、水质的数值模拟和验证表明,本文所建立的数学模型基本能反映钱塘江河道的实际情况,验证结果合理可靠,率定的参数可用于设计潮型下的水质预测。4规划期间的水质和敏感分析4.1污染源污染物排放量按照国家规范,水环境评价必须要在上游流量85%～90%设计保证率的水文条件下,其水质仍能满足水质标准。统计上游富春江电站多年径流量,保证率在90%的径流流量为300m3/s,但近两年,水资源日趋紧张,小于200m3/s的流量频频出现,针对具体情况,本着慎重的态度,认为对于杭州市的重要水源保护区宜采用保证率95%的流量为妥,其值为200m3/s。小潮稀释水体较少,污染物排放口形成的混合区浓度高,影响较大。因此,取小潮汛作为下边界是最不利的工况条件。规划期,富阳至小沙江段南北岸的负荷之比约为8∶1。江北片污水总量为15.7万t/d,CODCr排放量为9.5t/d,较现状污染负荷增幅为25%;江南片污水总量为76万t/d,CODCr排放量为75.2t/d,较现状污染负荷增幅22%。江北片NH3-N排放量为2.35t/d,江南片NH3-N排放量为4t/d。对于污染物CODMn,模型范围内所有污染源都有其对应的规划期负荷,因此可以计算绝对值。对于污染物NH3-N,仅收集到富阳渔山江段的规划期负荷,缺乏其他众多污染源的NH3-N负荷,因此NH3-N的计算采用“增量法”,集中反映富阳城镇污水排放对于杭州市饮用水源保护区的“贡献值”。表2为上游枯水流量遭遇下游小潮时在里山取水口水源保护区的CODMn和NH3-N增量全潮平均值。可见,规划期保护区上游断面(取水口上游4km)的污染物CODMn值为4.00mg/L,保护区下游断面(取水口下游3km)的CODMn值为3.98mg/L,总体接近Ⅱ类水体的临界值;NH3-N全潮浓度增量为0.12mg/L。4.2染负荷控制情况里山取水口的水质与上游径流量、上游来水水质、下游潮汐条件、污染负荷控制情况有一定的关系,但具体的敏感性如何,需要方案计算进行确定。本文拟定了7种工况,以CODMn指标为例,分析不同条件下的水质敏感性,具体的工况设置和计算的成果见表3。4.2.1南片和北片codn指标比较通过方案1,2,3的比较,分析污染负荷削减对取水口水质的影响。方案1为规划期负荷的工况,和现状相比,江北片CODcr负荷增幅25%;江南片CODcr负荷增幅22%,在这种工况下,保护区CODMn指标已经接近Ⅱ类水体临界值。方案2将富阳南片的负荷维持在现状水平,北片为规划排污量,在这种工况下,保护区的CODMn平均值有所下降,基本保持在Ⅱ类水体的标准值内。如果进一步削减南片的污染物,即方案3削减江南片污染负荷为现状的75%,江北片为规划排污量,则保护区的CODMn平均值将对应下降,但幅度不是很大。4.2.2黄竹片与北片的污染负荷变化对取水口水质的影响通过方案3,4的比较,分析上游来水水质对保护区的影响。随着上游城镇和经济的发展,污染物排放量增加,使上游来水的水质变差,其浓度的改变将对取水口水质产生影响。在富阳南片污染负荷为现状的75%,北片为规划排污量的负荷下,方案4上游边界浓度在方案3的基础上,增加15%。可知,尽管富阳南片的污染负荷在现状基础上削减了25%,但只要其上游的污染物没有加以控制,保护区的CODMn平均值仍满足不了Ⅱ类水体的水质要求。因此,应严格控制上游的污染物排放,其总量应保持现状的水平。4.2.3现有两种流速对保护区codn的影响通过方案1,5,6的比较,分析径流量改变对保护区水质的影响。在方案1富阳城镇污染负荷规划工况的基础上,考虑将上游径流量进一步削减为150m3/s(方案5)和抬升为300m3/s(方案6)2种工况。通过比较发现,径流量的改变直接影响保护区的CODMn平均值,在径流量削减为150m3/s时,保护区CODMn平均值超过Ⅱ类水体的水质要求;而将上游径流量提高到300m3/s,则保护区CODMn平均值有较大幅度的下降,可满足Ⅱ类水体的水质要求。4.2.4边界条件的确定通过方案1,7的比较,分析下游潮汐条件改变对保护区水质的影响。一般情况下,小潮期污染物稀释扩散条件差,而大潮期流速大,污染物迁移、扩散作用强于小潮,污染物推移上溯的距离会长一些。因此,在规划工况的基础上,考虑大潮的边界条件(方案7)。通过比较发现,大潮工况对于保护区CODMn平均值的改变量是明显的,它并没有将下游排放的污染物推移至影响保护区,而是有效地稀释了保护区的污染物浓度。5富宁县范围内水质分析(1)为了确保杭州市取水口的安全,必须设置一级保护区来加以保护,并严格遵守国家“饮用水水源保护区污染防治管理规定”中一级保护区有关禁止事项,严格加以保护。根据具体情况在里山取水口的上游4km、下游3km设定为一级保护区是必