北仑港三期工程废水排放对比仑海域水质影响的数值分析.pdfVIP

维普资讯 第 ll卷 第 3期 海 洋 通 报 Vo1．14．No．3 1995年 6月 M ARINE SCIENCE BULLETIN Jua．，19g5 一 ， ／ 北仑港三期工程废水排放对北仑 海域水质影响的数值分析 芏 姚炎明 耿兆铨 … 哺 斯 … z， 4／ 一。尹 麈术排敌对仑海譬 醣和藕 f ／． 美…|调皇黥苎 登登 查查匿髟髟响 锻缀 曰 f 沿海工程废水排放引起的海域水环境问题 已逐渐引起 们的重视 。废水中的污染物在海 水中的 自净过程十分复杂 ，主要通过污染物的扩散、迁移、降解和衰减等进行。而污染物的扩 散、迁移过程与海域的流场变化、潮位变化和海水本身的水质等因素有关 。因此，采用垂向平均 的二维数学模型来计算和预测污染物排海后的扩散、迁移规律是一项有益的工作，可为海洋环 境保护和污染治理提供依据 。 本文选用平面二维非恒定流的显式迎流有限元方法模拟计算了北仑海域的潮流，并在验 证的基础上 ，分别计算了北仑三期工程生活污水达标排放 、未经处理排放和船舶含油废水未经 处理排放时北仑海域的水质变化，达到预测北仑三期工程对海域水质影响程度和范围的目的。 1 模型建立 1．1 计算域的概化 北仑港是我国最大的深水 良港之一，港区海域在盒塘水道 。受峡道两侧地形的制约 ，水面 宽窄变化较大。如港区西面的穿山港南北水面宽约 8．0kni，而北仑 山与鳄鱼角两码头之间水 面又骤然缩窄到 3．Ikm，至水域东 口放宽到 4．5kni。水域东西两侧有大榭、黄蟒、中柱等岛屿 盘宣 (图 1)。水域面积约 80kni。，潮型属不正规半 日潮，潮流呈往复流特性。 水域 内潮强流急，温、盐分布较均匀，水体垂 向混合 良好 ，故宜选用平面二维的有限元方法 和三角形计算网格 ，既可节省机时和空闻，又可克服矩型网格对复杂边界的模拟困难 。 流场计算中计有 310个三角单元和 187个节点，最小空间步长约为 220m(图2)。计算浓 收稿 日期t1995—03—1O 维普资讯 海 洋 通 报 l4卷 度场时，在排放 El附近进一步加密 。流场计算 时间步长为 30s，浓度场计算时时间步长为 90s。 1．2 基本方程 数学模型控制方程采用经垂线积分含二 阶涡动粘滞项的浅水潮波方程和相应的二维 对流扩散方程 其 中潮流计算不考虑风应力 ， 浓度场计算不考虑污染物的蒸发、光合作用 和生物降解，以及悬浮或沉降等因素。即： 害一_缸 +。丛譬av：0 …(1) 围1计算域概况 搴+“磊9u,9u+塞+鐾一如 一 c塞，+c考 (2) 警+“+”舅+