Delft3D对污染物的扩散输运计算.pdfVIP

Delft3D 对污染物的扩散输运计算 朱宇新 (河海大学海洋学院，江苏 南京 210098) 摘要：应用水动力数学模型对江苏省如东县小洋口附近海域水质因子 COD 进行扩散场的模拟。模型采用 正交曲线网格对整个计算区域进行剖分，对大、小潮情况下的流速、潮位及扩散场进行了研究，并利用排 放口附近的现场实测水文资料进行验证。结果表明：水动力模拟结果与实测资料吻合较好，在大、小潮时 各个阶段中，流速与流向对 COD 的迁移扩散起着决定作用，且小潮落憩时刻污水的影响范围最大。 关键词：正交曲线网格；COD ；扩散输运；深度平均 在我国海洋水质监测中，通常以 COD 代表海水中有机物的耗氧量，本文选用此因子进 行模型的构建，用潮流场作用下污染物的对流-扩散方程进行污染物输运、扩散过程的研究。 研究区域地处江苏省中部，东经 120°56，北纬 32°36附近，东临黄海，南近长江口。 排污口设在小洋口附近，排海污水日排放量约为 2×104m3[1] 。由于排污口所在海域形成的混 合区，不能影响邻近功能区水质，排放口海域其水质需维持《海水水质标准》（GB3097-1997 ） 中第Ⅲ类标准之内，所以需要预报污染物排放的分布情况并分析污水排放对纳污水体的影 响。 1 数值模拟方法 水动力模型采用荷兰Delft水力研究院的Delft3D模型。该模型可以模拟在海岸、河流、河 口地区的浪、潮、流、水质、生态以及泥沙输运过程，并包含有前处理模块（网格生成 [2] （Post-Processing ，简称GPP ）。 RGFGRID 、数据插值QUICKIN ）和数据的后处理模块 1.1 平面二维水流、水质模型 为了研究污水排放中污染物的扩散及COD 的分布情况，使用Delft3D模型的二维水动力 模块对南黄海的流速和潮位进行模拟。采用东中国海大范围数学模型导出开边界处的潮位， 通过插值给出本文网格开边界的边界条件。水动力模块建立在Navier-Stokes方程的基础上， [2] 应用了浅水简化，采用交替方向法（ADI ）对该坐标系下的控制方程进行离散求解 。在正 交曲线坐标系ξ −η下，二维潮流运动连续性方程、动量方程和质量守恒方程以及对流-扩散 方程可分别表示如下： 沿水深积分的连续方程 [( ) ∂ +] ξ [( ) ∂ +] ξ d U G d V G ∂ς 1 ηη 1 ξξ + + 0 (1) ∂t ∂ξ ∂η G G G G ξξ ηη ξξ ηη 1 η ξ和 方向的动量方程 ∂U UU ∂V ∂U U ω ∂UV ∂ Gξξ + + + + ∂t ∂ ∂ d + ∂ ∂ G ξ G η ξ σ G G η ξξ ηη