北支咸潮倒灌对崇明岛东风西沙水库取水口水质的影响.docxVIP

北支咸潮倒灌对崇明岛东风西沙水库取水口水质的影响 崇明岛位于长江三角洲的交汇处，是世界上最大的沙岛。目前,崇明岛没有集中式供水水源地,以32座水厂供水为主, 这些供水厂分散于内河取水,规模偏小,且多数水厂建于20多a前,设备简陋,建筑物陈旧,不符合现代化水厂设计标准,因此供水能力不足的问题日益突出。 随着近年来崇明岛农业基础设施建设力度的加大,化肥和农药年使用量越来越大,内河水质受到严重污染。此外,随着长江口大通以下沿江引排水量的增大,长江口咸潮入侵呈加重态势,对崇明岛供水安全构成严重威胁。因此,崇明岛极有必要利用东风西沙夹泓兴建长江边滩水库以保障全岛供水安全。 鉴于水库取水水源的重要性和咸潮入侵的危害性,近年来不少学者对长江口咸水入侵开展了大量研究,对长江口咸潮入侵的途径、规律有了基本的认识。在此基础上,茅志昌等、王国峰等、肖成猷等就北支咸水倒灌南支对青草沙水库、陈行水库等水源地的影响进行了深入研究,顾莉等、刘树峰等也对近海水域的水环境质量进行了现状评价和模拟。由于水源地取水口所在河段属感潮河段,受径潮流影响的水动力条件复杂,因此研究的重点主要集中在咸潮入侵期间水库最长连续不可取水天数、可取淡水时间和盐度过程线等方面。由于针对北支水体倒灌南支影响取水口所在河段水质(除氯化物指标外)的研究较少,本研究以崇明岛东风西沙水库工程为例,在收集取水口、徐六泾断面同步水质监测资料和北支2010年洪枯季水质监测资料的基础上,采用delft-3D水质数学模型分析预测北支咸潮倒灌对取水口水质的影响。 1 模型和数值的选择 1.1 模型网格、网格划分及精度分析 模型研究范围包括大通以下长江下游河段、长江口、杭州湾和邻近海区,东边界到124.5°E附近,北边界到33°N附近,南边界到28°N附近,上游一直延伸至大通。采用Delft 3D软件自带的Rgfgrid工具生成曲线正交网格,网格生成过程中兼顾网格的正交性、拟合岸线,以及与河道走向尽可能一致。模型网格的分辨率较高,长江口内的分辨率为100～500 m,南北支分汊口区域的网格加密,在工程水域、北支进口及北支上段的分辨率为100 m左右,口外较疏,最大为10 km左右。本研究所用模型经过多次的水动力、氯(盐)度及水质要素的率定、验证,具有较高的精度,在长江口的多项研究工作中得到成功应用。限于篇幅,本文仅给出1999年1—4月崇头及陈行水库氯度观测站的验证效果,见图1。 1.2 扩散系数的计算 采用对流扩散方程研究北支咸潮倒灌对取水口水质的影响,二维对流扩散模型的基本方程为 ??t(hc)+??x(uhc)+??x(vhc)=??x(hDx?c?x)+??y(hDy?c?y)+fr(c,t)+S(1)??t(hc)+??x(uhc)+??x(vhc)=??x(hDx?c?x)+??y(hDy?c?y)+fr(c,t)+S(1) 式中:c为水质浓度;h为水深;t为时间;u、v分别为ξ方向、η方向的流速;Dx、Dy分别为ξ方向、η方向的扩散系数;S为源项;fr为生物学、细菌学、生态学、化学反应项。 1.3 计算边界条件 1.3.1 文化条件设计 以长江河口上游大通水文站P=97%保证率流量过程为基础,同时考虑三峡蓄水、南水北调、区间沿江引水影响,以减去三峡蓄水、南水北调、区间沿江引水影响后的流量过程为设计水文条件,以此水文条件为上游边界,下游外海边界选用对应时段的外海天文潮。上游大通流量边界条件见图2。 1.3.2 计算水质的限制 1.3.2. 同潮同潮水水质调查 a. 拟建取水口水质现状分析(枯水期)。于2008年12月15日(大潮期)、12月23日(小潮期)对取水口水域水质进行监测,结果见表1。大小潮期间,拟建取水口水域水质绝大多数指标达到Ⅱ类标准,主要超标项目为TP,其Ⅱ类达标率为45.8%;DO的Ⅱ类达标率为62.5%;Pb的Ⅱ类达标率为87.5%;石油类的Ⅱ类达标率为87.5%;其余项目的Ⅱ类达标率均为100%。 除DO、TP在大、小潮期的质量浓度相差一个水质类别外,其他水质指标在大、小潮期基本持平。 b. 北支口门、徐六泾断面、取水口断面的水质监测对比(平水期)。2009年4月27日(大潮)和5月5日(小潮)分别对徐六泾断面、北支口门断面、取水口断面进行了同步水质监测。监测结果显示,徐六泾、东风西沙水域总体水质较好,除TP、TN、Fe指标外,其他指标均能达到Ⅰ～Ⅱ类,见图3。TP为Ⅲ类标准,TN为Ⅴ类标准,Fe、Mn在个别潮期超过集中式生活饮用水地表水源地补充特定项目标准限值。 从空间分布情况来看,北支口门断面NH3-N、TP、TN、Fe、Mn、Zn、Cu等水质指标均高于徐六泾断面和东风西沙取水口断面。 c. 北支水质现状分析(枯水期和丰水期)。自20世纪60年代以来,长江径流在北支的分