1. 引言
在我国海洋水质监测中,通常以CODcr代表海水中有机物的耗氧量,本文选用此因子进行模型的构建,用潮流场作用下污染物的对流扩散方程进行污染物输运、扩散过程的研究。研究区域斯道拉恩索(广西)林浆纸厂排污口B3点,位于广西铁山港湾内[1] (见图1)。排污口设在铁山港航道的西北侧。
由于排污口所在海域形成的混合区,不能影响邻近功能区水质,排放口海域其水质需维持《海水水质标准》(GB3097-97)中第三类标准之内,所以需要预报污染物排放的分布情况并分析污水排放对纳污水体的影响。
![](//html.hanspub.org/file/23988x7_hanspub.png)
Figure 1. the comparison diagram of location, tide, tidal current verification point and the current velocity in research area
图1. 研究区地理位置、潮位、潮流验证点及流速流向对比点图
2. 数值模拟方法
2.1. 三
维水动力、水质模型
本模型采用三维河口海岸海洋模式(ECOMSED),该模型可以模拟在海岸、河流、河口地区的浪、潮、流、水质、生态以及泥沙输运过程。水动力模块建立在N~S方程的基础上,采用有限差分(半隐式)数值解法对σ坐标系下的控制方程进行离散求解。三维潮流运动连续性方程、动量方程和质量守恒方程以及对流扩散方程可见文献[2]。Ecom水质模型[3]是在水动力模拟的基础上,借助于质量守恒原理,考虑物质由于对流、紊动扩散[4]等因素计算污染物在水体中的输运、扩散以求得污染物浓度在水域分布的状况。
为提高潮流场计算的准确性,在计算域内涨潮淹没、退潮露滩的潮间带,利用干湿判别法的动边界处理技术,使计算区域内的潮间浅滩动态变化,真实地反映出潮涨、潮落潮过程中,潮间带水陆边界的移动过程。干湿网格的判别标准为水深 ≤ 0.05 m。在对流扩散中,规定流速为零的“干点”没有扩散;由于水动力的计算区域比较大,排污口距开边界较远,为了在计算过程中开边界不会对计算区域内的扩散有影响,所以规定开边界处没有扩散。
2.2. 模型及参数的选取
本文计算域设置、水深和岸界条件、模型水边界输入条件与文献[2]相同,其他模型参数设置如下:
1) 风场参数选取
废水扩散的季节性差别体现在季风风向的差别,根据北海地区冬季海面风的多年平均统计资料,冬季季风主要考虑为N向风的影响,夏季季风为S向风的影响,选取铁山港附近海域6级风速10.8 m/s条件下对污水排海后水质影响进行模拟。
2) 扩散系数
根据紊动扩散理论及对实测资料的修正,选取横向与纵向扩散系数[5];降解系数K与水流特性、海洋生物的种类和数量、海面气体交换速度、污染物浓度及水温等因素有关,参考国内外有关手册[5],模型实际采用的K一般取为常数[5],且水温范围在10℃~28℃时,降解系数K值介于0.02~0.07 d−1之间。根据北部湾沿海的具体情况,选择水温
21 ℃
时研究区域内的降解系数K为0.03 d−1。
3. 验证计算及结果分析
3.1. 潮位与流速验证
水文资料选择2007年7.2~2日大潮4个潮流站( #1~#4)及2个潮位站(A、B)实测的潮流资料。
根据上述资料和条件进行计算,对各测点流速和流向的实测结果和计算结果进行对比,结果列于图2~5。由图可见,三维模型的水位验证良好,三维模拟结果的流速呈现表、中层大,底层小的特点,而且流态也较为合理,基本上能反映出工程附近水域的水流状况,可以作为进一步分析计算的基础资料。
3.2. 风吹流场模拟结果分析
在工程区外侧区域选取代表点进行流速、流向对比说明,如图1所示。
3.2.1. 风吹流对流速的影响
图6~9是本海域风吹流的三维水动力数值模拟结果,选取了涨急和落急时刻受风应力影响大的表层的流场矢量图,中层和底层受风应力影响较小,由于本文篇幅有限,不给出矢量图。
1) 表层流速
N向风:在强风作用下,水体表面受风应力的作用,形成风吹流。风吹流与潮流相互作用,根据风
![](//html.hanspub.org/file/23988x9_hanspub.png)
Figure 3. surface flow velocities and direction validation curve
图3. 表层流速、流向验证曲线
向、流向的不同,流场显示不同的变化特性。计算海区涨潮时遇N向风,流向与风向相反,迎风逆流使得风吹流后的流速与无风方案相比涨潮流速减小,其中1#断面涨潮最大流速减小8%~43%,2#断面减小13%~30%,由河口向外海流速减小趋势加大。落潮时风向与流向相同,顺风顺流故流速增加5%~13%。
S向风:S向风对计算海区的表层流场影响与N向风相反,涨潮流速增加落潮流速减小。涨潮1#断面流速增加5%~42%,2#断面增加5%~29%,落潮1#断面减小25%~49%,2#断面减小35%~57%。
2) 中层流速
根据海水具有连续性和不可压缩性的特点,中层水体以对应表层流速变化相反的方式,形成一个闭合的循环。
N向风:1#断面涨潮流速增加5%~31%,2#断面增加8%~23%。落潮1#断面减小2%~10%,2#断面减小4%~7%。
![](//html.hanspub.org/file/23988x10_hanspub.png)
Figure 4. middle velocity and direction validation curve
图4. 中层流速、流向验证曲线
![](//html.hanspub.org/file/23988x11_hanspub.png)
Figure 5. bottom velocity and direction validation curve
图5. 底层流速、流向验证曲线
![](//html.hanspub.org/file/23988x12_hanspub.png)
Figure 6. velocity vector of surface-layer of N wind at the flood tide
图6. 表层N风作用下涨急时流场图
![](//html.hanspub.org/file/23988x13_hanspub.png)
Figure 7. velocity vector of surface-layer of N wind at the ebb tide
图7. 表层N风作用下落急时流场图
![](//html.hanspub.org/file/23988x14_hanspub.png)
Figure 8. velocity vector of surface-layer of S wind at the flood tide
图8. 规划后表层S风作用下涨急时流场图
![](//html.hanspub.org/file/23988x15_hanspub.png)
Figure 9. velocity vector of surface-layer of S wind at the ebb tide
图9. 规划后表层S风作用下落急时流场图
S向风:涨潮流速1#断面减小3%~27%,2#断面减小2%~14%。落潮1#断面、2#断面的增加值为9%~35%、26%~37%。
3) 底层流速
两个风向的风吹流计算后,底层流速不受影响,无变化。
3.2.2. 风吹流对流向的影响
1) 表层流向
N向风:计算海区北风为离岸风,涨、落潮流向受N向风顶托、压迫向东偏,1#断面偏转16˚~46˚,流向变化最大处在1-5#点附近,2#断面偏转19˚~25˚。落潮流向1#断面流向变化14°~39°,最大变化在1~7#点附近,2#断面流向偏东17˚~20˚。
S向风:S向风计算后流向的偏转方向为,涨、落潮流向西偏。1#、2#断面涨潮流向变化4˚~11˚、5˚~11˚;落潮流向变化4˚~16˚、8˚~12˚。
2) 中层流向
6级强风对计算海区中层水体的流速方向影响不大,N向风计算后,涨、落潮流向变化1˚~2˚,S向0˚~1˚。
3) 底层流向
N、S风向的风吹流计算后,底层流向不变。
综合上述,风吹流计算后的表层流速变化为:N向风(离岸风)涨潮流速减小,落潮流速增加;S向风(向岸风)涨潮流速增加,落潮流速减小。
4. 风吹流作用下污染物CODcr的扩散情况
4.1. N风作用下污染物扩散数值模拟
本次预测计算中CODcr与CODMn的换算关系考虑为2.5,该海域COD本底值为1.5 mg/L,排放量为109,158 m3/d,排放浓度为100 mg/L,污水中CODcr排放量为9.82 t/d。预测计算时分别以此作为计算源强。对COD扩散连续进行15d的计算,并将计算结果叠加COD本底值。
各层预测计算结果见图10及表1,其中各层的包络图为污染物扩散连续进行15 d的计算结果统计。结果显示污水从扩散器进入海水后,100 mg/L排放浓度下中层最大,表层最小。
4.2. S风作用下污染物扩散数值模拟
根据规划,本次预测按1种排放浓度进行预测计算( 计算中CODcr与CODMn的换算关系考虑为2.5,根据资料该海域COD本底值为1.5 mg/L),排放量为109158 m3/d,排放浓度为100 mg/L,污水中CODcr排放量为9.82 t/d。预测计算时分别以此作为计算源强。对COD扩散连续进行15 d的计算,并将计算结果叠加COD本底值。
各层预测计算结果见图11及表2,其中各层的包络图为污染物扩散连续进行15 d的计算结果统计。结果显示污水从扩散器进入海水后,100 mg/L排放浓度下中层最大,表层最小。
4.3. 污水中COD对海洋环境的影响
根据上述研究结果可知:污水排放后表层、中层和底层的污水总体上随涨、落潮沿岸线往复运动,垂直岸线方向的运动幅度很小,这有利于减少污染物对岸滩和敏感区的影响。从垂向看,排放污水在上升过程中浓度不断降低,在垂直方向上呈现排污口附近的中层浓度最高、表层浓度相对最低的特点。
![](//html.hanspub.org/file/23988x16_hanspub.png)
Figure 10. CODcr each layer’s influence envelope of 100 mg/L
图10. 规划后COD排放浓度为100 mg/L时各层的影响包络
![](Images/Table_Tmp.jpg)
Table 1. sewage disposal CODcr environmental impact analysis on water environment
表1. 污水排放COD对水环境影响分析
![](Images/Table_Tmp.jpg)
Table 2. sewage disposal CODcr environmental impact analysis on water environment
表2. 污水排放COD对水环境影响分析
![](//html.hanspub.org/file/23988x17_hanspub.png)
Figure 11. CODcr each layer’s influence envelope of 100 mg/L
图11. 规划后COD排放浓度为100 mg/L时各层的影响包络
总体来说,从垂向看排放口附近各特征水层的浓度存在一定区别,中层喷口附近的水域存在一个高浓度区,随着污染羽流上升,不断与周围海水掺混,垂向浓度梯度较大(浓度迅速降低)。
CODMn各层的超海水水质标准的影响范围都很小,不会对营盘海水养殖区、儒艮保护区及红树林保护区产生直接危害。
5. 小结
研究结果表明,研究区海域水动力条件和污染物扩散规律具备以下特点:
1) 风吹流计算后,垂向流速的变化特点为:强风作用于水面,风向与潮流向综合作用造成表层流速增、减,中层水体以对应表层流速变化相反的方式形成闭合的循环,底层流速几乎无变化。
2) 风吹流与潮流综合作用后的表层流场与无风潮流场相比,流速方向受风向影响有所变化,但流速方向仍由潮流向控制。中层流向变化很小,底层流向无变化。
3) 在N、S向风作用下使得表面水流加强,结果更有利于各层与层之间掺混,由预测结果可知对表层、中层污染物扩散产生一定影响。
4) 在垂直方向污染物规律与无风条件相似,随着不断与周围海水掺混,排放污水在上升过程中浓度不断降低,在垂直方向上呈现排污口附近的中层浓度高、表层浓度最低的特点,在有风条件下各浓度等值线较无风条件下向岸边靠近。
5) 在季风影响下,污染物扩散范围较无风条件下略大,但是差别不明显,从各层超海水水质标准的最大影响范围来看没有对敏感目标产生影响。
因此,该水质模型模拟结果可以更准确的确定环境保护目标是否受到海洋工程排污的影响,研究方法可以作为海域环境污染研究的可靠方法,研究成果有助于海洋环境保护的实施,可以作为环境保护部门的参考依据。
资金项目
中央级公益性科研院所基本科研业务费专项资金项目(TKS160227, TKS160209).