钦州湾三墩作业区排污阈值的数值模拟研究*

裴木凤，高劲松

(1.国家海洋局北海海洋环境监测中心站，广西北海 536000；2.南宁师范大学北部湾环境演变与资源利用教育部重点实验室，广西南宁 530001)

0 引言

钦州湾位于广西壮族自治区南部钦州市以南，是广西北部湾重点开发的海湾之一，随着向海经济的发展，钦州湾大面积海域被围填成陆，如2005-2015年间，钦州湾围垦规模近3 000 hm2[1]，钦州湾岸线发生较大改变。与之对应的是近年来钦州湾近岸海域局部水质下降和海洋生态环境发生剧烈变化[2-4]，尤其是钦州市茅尾海水质长期达不到功能区目标。导致钦州湾近岸海域水质下降的原因，主要是入海河流携带大量污染物入海、非法及不合理设置入海排污口或不达标废水排放入海，使得污染物的入海量超出河口海湾的环境承载能力。为逐步改善钦州市近岸海域环境质量，钦州市政府部门提出钦州市入海排污口环境综合整治和污水外迁至三墩排放口的方案，并确立钦州湾海域水质“近期有所改善，远期全部达标”的目标。

排污容量是水污染防治和污染减排工作的重要依据，其研究对保护海洋生态环境有重要意义。早期排污量主要基于实测数据进行计算[5]，随着计算机技术的发展，使用数值模型计算区域环境容量和排污量取得了长足发展[6,7]。目前，关于钦州湾的数值模拟有较多的研究成果[6,8-11]，但关于钦州湾排污量的研究仍然少见于研究文献中，其中夏华永等[6]采用二维潮流数值模型和水质模型探讨了钦州湾污染物的扩散规律，利用1993年海域现场观测数据计算了钦州湾大规模填海之前的海水交换率、平均半更换期及环境容量；郭森等[8]采用多源线性迭加的方法计算了钦州湾金鼓江东、西两侧两个排污口的排污量。

在夏华永等[6]的研究中，排污总量是基于整个海湾整体纳污能力的计算，无法为实际排污口的污染减排提供依据，而郭森等[8]的研究则主要考虑金鼓江东、西两侧两个排污口的排污情况，对钦州湾的水污染防治参考有限，此外，两者的研究均没有考虑近年来钦州湾大规模填海造成的岸线及地形的变化。三墩作业区位于钦州湾外湾中部，已规划为广西钦州石化产业园和中国(广西)自由贸易试验区的重要组成部分，将与金鼓片区形成国家级石化产业基地，区域的环保设施服务水平亟需提升，而该区域的排污管理也是改善钦州湾水质的重要部分。本研究基于钦州湾的现状海岸线及地形情况(图1),建立钦州湾二维数值模型，结合实测资料，对三墩作业区的排污口进行选划及最大排污量计算，可为管理部门提供决策依据，保护钦州湾的海洋生态环境。

图1 钦州湾水深地形

1 材料与方法

1.1 EFDC数值模型

数值模型采用Environmental Fluid Dynamics Code (EFDC)模型。该模型是由美国弗吉尼亚海洋科学研究所开发，并由数家科研单位后期维护下发展起来的，可用于近岸海洋及陆地地表水水动力-水质模拟的水环境数学模型[12,13]，该模型也是美国最大日负荷总量(TMDL)等环境保护计划主要使用的模型。

1.1.1 水动力数学模块

水动力数学模型基于静水压假定和Boussinesq近似，对不可压缩的自由表面水体进行描述，其控制方程如下：

(1)

(2)

(3)

1.1.2 水质模块

水质变量的质量控制方程由物理输运、平流扩散和动力学过程组成，公式如下：

(4)

其中C为物质浓度，u、v为x、y方向的速度，Kx、Ky分别为x、y方向上的扩散系数，K为一阶降解速率，R为源汇项。

1.2 模型配置与验证

1.2.1 模型配置

模型计算区域北至茅尾海顶部21.95°N，南至钦州湾湾口21.37°N，西至108.389°E，东至108.856°E(图1)，水平网格为520行×411列，有效网格数共104 875个，网格依据地形及工程区域而大小不一，三墩作业区附近海域的网格分辨率约为60 m×100 m，动边界的处理采用干湿节点法。

海底地形数据由海图水深插值处理得到；潮汐开边界采用调和常数(M2、S2、N2、K2、K1、O1、P1、Q1、Sa)确定；模型的初始条件，包括流速初始场和水位场(开边界除外)均为0，预热运行30 d后达到稳定状态的结果作为正式计算的初始条件；时间步长为5 s。

1.2.2 模式验证

国家海洋局北海海洋环境监测中心站于2012年3月3日、4日、6日、7日于钦州湾布设了4个站位进行海流观测(观测仪器：亚力克AEM-USB)。本研究使用该实测数据对钦州湾的潮流进行验证，此外还使用国家海洋局北海海洋环境监测中心站在钦州站2012年3月的实测潮汐资料进行潮位验证。各站位分布见图1。

对比模型计算结果与实测数据，模拟潮位与实测潮位基本吻合(图2)，小高潮时的潮位误差略大，验证时段潮位的平均误差约0.2 m，相关系数为0.97；流速、流向验证见图3，模拟的流速、流向与实测值的变化趋势大体一致，基本反映了区域涨落潮的往复流特性。总体上，模拟结果可代表模拟区域的流场状况，可以为区域污染物的扩散模拟提供实际的潮流场。

图2 钦州站潮位验证结果(2012年3月1日0：00-31日23：00)

图3 站位H1-H4流速、流向验证结果(2012年3月)

2 结果与分析

2.1 钦州湾流场分布

从图4可以发现钦州湾落潮流大于涨潮流，航道和深槽处流速最大，流向与航道、深槽走向基本一致，浅滩和岸边流速相对较小，流向多变，最大流速位于龙门湾颈处，涨潮流最大约1.8 m/s，落潮流最大约2.0 m/s。外海的西南向流在三墩作业区前沿转为北偏西进入钦州湾，通过3条主要深槽在青菜头汇合后以西北向直达龙门，在龙门口以西北向茅尾海发散，落潮流流向则基本与之相反，钦州湾的潮流呈往复流性质。涨、落潮时，保税港区、三墩作业区和三墩公路所围区域为相对弱流区，三墩公路东侧进、出鹿耳环江的潮流基本与公路平行，三墩作业区西北侧涨潮流为东北向、落潮流为西南向。本研究的钦州湾潮流场特征与前人的模拟结果相似[6,10,11]。

(a) 涨急流场，(b) 落急流场

2.2 排污口的设置

根据《广西壮族自治区近岸海域环境功能区划调整方案》(桂政办发〔2011〕74号)[17]，排污口拟选范围为三墩污水深海排放区(GX071DⅣ)，具体以108.679 167°E、21.603 056°N为中心，半径1 km的海域，海域主导功能为港口、工业、生活排污用海，水质目标为《海水水质标准》[18]第四类标准，此区域也是《钦州市海洋功能区划(2008-2020)》[19]划定的三墩污水深海排放区。

根据规划排污区、三墩规划填海区等情况，初步选择4个排污口位置，分别为规划排污区西端(A0)、中部北端(B0)、东端(C0)和中部南端(D0)，此后根据模拟结果再进行局部加密(A1、B1和C1)，规划排污区和排污口位置设置情况见表1和图5。

图5 三墩区拟选排污口位置示意

表1 排污口位置

在拟选排污口设置相同的污水排放量及出水标准(污水排放量2 800 m3/h，无机氮浓度40 mg/L)，模拟计算一个大小潮周期(15 d)，输出每个小时的计算结果进行统计，结果见表2。根据模拟结果，不同排污口的物质扩散方向基本为SW-NE向，其中西端即A0处排污的最大浓度值和物质扩散范围最小，东端C0处排污的最大浓度值和扩散范围最大，其余位置的相应值在两者之间，总体呈现西南区排污口的污染物扩散范围小于东北区排污口的污染物扩散范围，而不同排污口物质浓度增量较大(大于0.5 mg/L)的范围均近乎圆形。因此三墩排污口的最优位置可设置在A0处。

表2 不同排污口的无机氮扩散情况统计

2.3 三墩作业区的排污阈值

三墩作业区拟纳入广西钦州石化产业园，该产业园重点布置石化产业及石化物流仓储，因此污水排放标准以《污水综合排放标准》(GB 8978-1996)[20]一级和《石油化学工业污染物排放标准》(GB 31571-2015)[21]中的表2“水污染物特别排放限值”的低值为执行标准，即CODCr为60 mg/L，折算成CODMn为24 mg/L，总氮为40 mg/L，磷酸盐为0.5 mg/L。在拟定的排污点A0周围半径约1.0 km处设置水质控制点，水质控制点处的模拟计算值叠加现状值基本满足四类海水水质标准时，即为相应物质的最大允许排污量。

表3 三墩排污口的排污阈值

3 结论

本研究构建了钦州湾水动力模型，模拟值与实测值吻合较好。同时构建一阶降解系数的物质对流-扩散模型，研究三墩作业区不同排污口的污染物扩散情况，选取钦州湾三墩作业区的最优排放口，并依据环境水质标准计算该海区化学需氧量(CODCr)、总氮(TN)和总磷(TP)的排污阈值。研究结论如下：

钦州湾落潮流大于涨潮流，涨潮流最大约1.8 m/s，落潮流最大约2 m/s，外海的西南向潮流在三墩作业区前沿转为北偏西入湾，之后在青菜头以西北向通过龙门向茅尾海发散，落潮时流向则相反，三墩作业区西北侧涨潮流为东北向、落潮流为西南向。钦州湾流场呈往复流特征。

三墩排污区内不同位置排放的物质扩散方向为SW-NE向，总体呈现西南端排放的污染物扩散范围小于东北端排放的污染物扩散范围，最优排污口可设置在排污区的西南端(21.601 101°N，108.670 364°E)。

在CODCr、TN和TP的出水标准分别执行≤60 mg/L、≤40 mg/L和≤0.5 mg/L情况下，满足四类海水水质标准时拟选三墩深海排放口的CODCr、TN和TP的排污阈值依次为76 845.25，1 581.7和150.2 t/a。