城市面源污染特性及污染负荷预测模型的研究.docVIP

城市面源污染特性及污染负荷预测模型研究1 叶 闽 杨国胜 肖 彩 长江水资源保护科学研究所 武汉 430051 摘 要 本文通过分析武汉市汉阳地区城市面源污染的来源和污染特征，并针对城市土地利用和下垫面类型特点，建立城市面源污染预测模型，研究探讨城市降雨径流及污染负荷预测模拟和迁移规律。所建模型已成功地应用于武汉市汉阳地区城市面源污染控制技术研究，获得了满意的结果。该模型可为城市水环境规划、面源控制和管理研究提供科学手段。 关键词 城市面源 污染特性 预测模型 1 前 言 随着我国城市化水平的发展进程加快，可渗透地面的面积比例越来越小，由暴雨径流产生的突发性的、冲击性强的城市面源污染已成为受纳水体水环境污染的主要来源之一。城市暴雨径流中包含有大量的污染物，使得其初期径流产生的污染负荷远高于城市生活污水。因此，由面源引起的水环境污染已成为当前城市水环境综合治理中亟待解决的主要问题之一。 近十年来，武汉汉阳地区的城市面源污染加剧，大量污染物排入水体，水质急剧下降，由城市面源污染和其他污染引起的水环境问题已经严重制约了城市的可持续发展。本文通过分析武汉市汉阳地区城市面源污染的种类、来源和污染物特征；并针对城市土地利用和下垫面类型特点，建立城市面源污染预测模型，研究探讨城市降雨径流及污染负荷预测模拟和迁移规律。为有效控制城市面源污染提供强有力的科学手段。 2 汉阳地区城市面源污染来源及特性 武汉市汉阳地区位于武汉市西南部，长江、汉江交汇处的三角地带。为武汉三镇一方。以汉阳区汉阳大道和鹦鹉大道为轴区，两面临江（东北临汉江，东南临长江），硬质下垫面硬质下垫面质下垫面SWMM）和城市暴雨径流模型（STORM），化学物质径流负荷与流失模型（CREAMS）。另外，在此之后开发的MIKE－SWMM模型，将一维非恒定流模型与暴雨径流的分布式模型结合，模拟暴雨水和废水排放系统中管道有压流、排水渠和地表漫流等。目前，模型的发展重点是将地理信息系统GIS与原有的面源污染模型相结合，评价面源污染的发生区域，污染物的运移变化、控制措施的影响以及受纳水体的响应等。如土壤水体评价模型（SWAT，Soil and Water Assessment Tool），已经在美国18条主要河流系统得到了广泛应用。 本文结合汉阳地区城市面源污染的特性，将城市地表径流产、汇流联系在一起，基于单元网格产流产污、汇流模型，以及排水管网的水动力学模型和污染迁移模型，建立了分布式城市面源污染模型，模拟暴雨径流过程中城市面源污染特征及其变化规律。 3.2 建立城市面源径流模型 根据城市区域特点，其降雨径流产生呈三种形式，地表径流、屋面径流和排水管网径流，排水管网径流将滞后前两种形式，而这三种模拟方法也是不同的。 （1）地表径流模型 根据研究区域特点，将不同的地表类型归纳为三种类型的单元，即不透水有滞蓄A1、不透水无滞蓄A3和透水有滞蓄单元A2。并根据地表子单元的特征为具有相同的坡度，其水文参数如曼宁值，洼地蓄水，下渗为常量。通过联合曼宁方程的连续性方程建立非线性滞蓄方程。 对一个子单元，连续性方程为： （1） 其中，V=A·d—单元中的水量，m3；d—水深，m；t—时间，s；A—单元表面积，m2；i*—净雨=雨强–蒸发–下渗，m/s；Q—出流，m3/s。 出流用曼宁公式表示： （2） 其中，W—子集水区宽度，m；n—糙率；dp—滞蓄水深，m；S—子集水区坡度。 （3） 其中， (4) 方程（3）的求解是在时间步长内通过简单的有限区间的方法。将方程（3）可近似表示为： （5） 其中，△t=时间步长，s. 然后方程（5）的求解可用Newton - Raphson 迭代方法，得出h2。然后则可计算时间t+△t 的瞬时出流流量Q。Q也是排水系统的入口流量。 （2）屋顶产流模型 屋顶的特征为不可透水，屋顶净雨量为雨强I、屋顶面积S以及时间步长△t三者的积。当下泄水量小于排水管最大设计能力时，产流量为净雨量除以时间步长，即: (6) (7) 其中，Q0=净雨量，m3；Q=屋顶产流量，m3/s；I=雨强，m/s；S=屋顶面积，m2；△t=时间步长，s。 假定屋顶排水管的最大设计流量为Qs，则有如下关系： ， （8） ，