污染物输运数值模拟方法.docxVIP

污染物输运数值模拟方法

适用范围

本附录规定了污染物输运数值模拟的原则、方法、内容及要求。

属于宽浅型水域且潮混合较强烈、各要素垂向分布较均匀的近岸海域或河口、海湾，可采用二维数值模型近似描述污染物的三维运动；潮混合不充分、各要素垂向分布不均匀的海域，海域水文条件较复杂或模拟分辨率要求较高等状况时，应采用三维污染物输运数值模拟方法。

模型计算域的确定及网格剖分

模型计算域确定

模型计算域的与潮流模拟区域相同。

网格剖分

网格剖分应符合：网格大小应不小于潮流网格，并能够反映污染物的输运规律和特征，满足海洋水质环境评价内容的预测需求。

污染物输运数值模拟

基本资料

用于污染物输运数值模拟的基本资料应按照本规范相应的规定获取，选取现状调查的资料作为背景值。

基本方程

二维污染物输运方程

二维污染物输运的基本方程为：

………(H.1)

式中：

——为水深(m)；

——为某种污染物浓度(mg/L)；

，——分别为向及向紊动扩散系数()；

——为污染物源(汇)项。

三维污染物输运方程

三维污染物输运的基本方程为：

………(H.2)

式中：

C——为某种污染物浓度（mg/L）；

，，——分别为向、向和向紊动扩散系数（m2/s）；

Q——污染物的源或汇。

计算模式

应根据计算域地形特征、项目布置方案等具体情况，可采现行计算模式进行计算。

初始条件和边界条件

初始条件

二维初始条件按下列公式确定：

……(H.3)

三维初始条件按下列公式确定：

………………(H.4)

边界条件

当计算域内存在大面积潮间浅滩时，宜采用动边界技术处理露滩问题。

固边界

固边界应按下列方法确定。

法向污染物通量为零：

…………………(H.5)

式中：

——为固边界法方向。

水边界

污染物按入流和出流情况分别控制。

二维水边界应按下列方法确定。

入流时：

……(H.6)

出流时：

…………(H.7)

式中：

——水边界；

——已知污染物浓度；

——法向流速（m/s）。

三维水边界的边界条件可按下列方法确定：

入流时：

…………………(H.8)

出流时：

………(H.9)

式中：

——水边界；

——已知污染物浓度；

——边界法向流速。

水边界

三维水面边界条件可按下列各式确定：

…………………(H.10)

基本参数

二维污染物输运紊动扩散系数和可取与相应的水流紊动粘性系数、相同数值。

三维污染物输运紊动扩散系数、宜采用试验或经验公式确定，可取50m3／s～500m3／s。宜采用试验或经验公式确定，也可取的值。

验证计算及精度控制

验证计算应通过参数的调整，满足模拟计算结果与实测结果基本相符的要求，同时应满足验证计算精度的要求。

计算成果

通过项目建设前后的模拟计算，应给出各方案的污染物输运的浓度分布特征等模拟结果，并附以相应的图表。

非保守污染物输运数值模拟

上述模型仅适用于保守性物质。对于非保守性物质，进入海洋后会发生一系列的生物、化学过程转化，从而使其浓度发生变化，因此，可根据现行研究成果，结合不同的污染物的特性适当增加方程的内容。