1. 引言
红海湾是位于我国南海北部的一个半封闭的海湾,地处汕尾市东部,东临碣石湾,西依大亚湾,三面环海。随着近年红海湾内众多工程的建设,有必要摸清湾内的动力特征。
南海的环流复杂[1],风场与沿岸地形所形成的强风应力旋度可能是其主要的驱动机制[2]。粤东沿岸流的流向随季风的转换而改变,在冬季,粤东沿岸流流向SW,夏季,粤东沿岸流流向NE[3]。
针对以上特点,于2015年在红海湾海域进行了冬、夏季全潮水文测验[4],并采用数学模型进行了模拟分析。
2. 红海湾水文观测及分析
红海湾海域于2015年夏、冬季分别进行过一次全潮测验,观测站位相同,图1~4分别为夏季大、小潮和冬季大、小潮的潮流矢量图。观测结果表明:
1) 红海湾的潮波并不是单纯的驻波或者前进波,而是介于驻波与前进波之间,兼有驻波与前进波的特征,为不规则半日潮流主导的海区。
2) 湾内潮差小,潮流弱,流态复杂。
3) 夏季中潮和小潮期间红海湾西岸的海流由东北向流向湾内,东岸的海流东南向流向湾外,中部的海流以ENE向为主,外侧的海流沿等深线向东,整个海湾呈现顺时针旋转特征。大潮期间整个海湾西岸海流由西南向流出湾外,其他方位的海流比较混杂,不是很一致,总体来看,整个观测区域呈现逆时针旋转特征。
4) 冬季海湾外部的海流都是沿着等深线从东向西流。中潮和小潮期间,海湾内的海流为西岸东北向流入湾内,东岸东南向流向湾外,海湾内的海流呈现顺时针旋转特征。大潮期间,海湾内的海流为西岸西南向流向湾外,东岸西北向流向湾内,海湾内的海流呈现逆时针旋转特征。
3. 数值模拟方法
海流计算采用Mike系列软件中的二维水动力模块(MIKE21 FM模块) [5]。该软件由丹麦水工所开发,可以应用于海洋、海岸、河口区域的二、三维水动力计算。FM模块采用三角形网格,在处理潮流动边界、负责工程建筑物边界方面具有强大的功能,在国内外许多工程项目研究中得到了广泛应用。
![](//html.hanspub.org/file/23990x7_hanspub.png)
Figure 1. perpendicular average sea current vector graph during spring tide in summer
图1. 夏季大潮期垂线平均海流平面分布矢量图
![](//html.hanspub.org/file/23990x8_hanspub.png)
Figure 2. Perpendicular average sea current vector graph during neap tide in summer
图2. 夏季小潮期垂线平均海流平面分布矢量图
![](//html.hanspub.org/file/23990x9_hanspub.png)
Figure 3. perpendicular average sea current vector graph during spring tide in winter
图3. 冬季大潮期垂线平均海流平面分布矢量图
![](//html.hanspub.org/file/23990x10_hanspub.png)
Figure 4. perpendicular average sea current vector graph during neap tide in winter
图4. 冬季小潮期垂线平均海流平面分布矢量图
3.1. 控制方程
连续方程:
(3.1)
运动方程:
(3.2)
(3.3)
温度输运方程:
(3.4)
式中:
t为时间,x, y为笛卡尔坐标;
为水位,d为静水深,
为总水深;u和v为x和y方向的水深平均速度;f为柯氏力参数,
为旋转角速度,
为地理纬度;g为重力加速度;
为水的密度,
为底部应力,
为水面风应力,
为参考水密度; 为大气压力,
水面上空10 m处风速;
为源汇项。
3.2. 模型设置及参数的选取
1) 模型设置
模型计算范围约为81 km*66 km,模型计算范围及水深情况见图5,计算网格见图6,区域内有51978个结点,这些结点将计算区域分成26653个互不重合的三角形单元。空间步长最小为20 m,最大空间步
![](//html.hanspub.org/file/23990x25_hanspub.png)
Figure 5. Calculation range and verifying point position
图5. 数模计算范围及验证点位置图
长约2 km,计算时间步长从0.01~20 s自动调节。计算范围内的水深及岸线采用中国航海图书出版社出版发行的海图“田尾角至大昆山”。
2) 开边界设置
本研究开边界由南海区域环流预测结果提取得到[6]。因红海湾内流态复杂,除潮差影响因素外,还受到环流、局部地形引起的沿岸流等影响,因此只考虑潮流动力的模型不符合本海域的动力特征。而南海区域环流预测模型充分考虑了季风、温度、盐度等综合因素的作用,以此作为开边界可以较准确的模拟红海湾海域内的海流。
3) 温度模型参数设置
本次模型研究依据厂址周边的汕尾气象站夏季7月和冬季12月的实际测量气温、湿度、风速、风向、降雨、蒸发等资料进行计算。
4. 水动力模型验证计算及结果分析
4.1. 模型验证
为了验证海流模型的合理性,采用2015年7月与2015年12月水文观测资料对模型潮位、流速、流向进行了验证。共有13个潮流站(图5中V1~V13)和4个潮位站(图5中小漠、汕尾-2、汕尾、马官)。限于篇幅,在红海湾内选取一个代表点,图7、图8为夏季和冬季实测潮位与计算结果的比较,图9、图10分别为夏季大、小潮和冬季大、小潮代表站位的流速、流向实测与计算结果的比较。可见,模型计算潮位与实测潮位拟合较好,流速、流向在连续的变化过程中与实测值比较接近,流态合理,模型计算结果与实际测量的变化趋势基本一致,本研究所建立的海流模型较全面的反映了红海湾海域的流动情况。
4.2. 模拟结果与分析
图11~14为夏季大、小潮期间涨急与落急流场。夏季大潮涨急时刻,涨潮流自红海湾东侧向西进入湾内,顺岸线逐渐向西,最后由西侧向西流出红海湾,计算的流场与实测的流向基本一致。夏季大潮落急时刻,落潮流则由红海湾西南侧进入湾内,在湾东部沿东岸线方向流出红海湾,湾东部的流速大于西部,计算的流场与实测主流方向基本相符。从流速值分布看,夏季大、小潮期间红海湾流速分布相同:湾外的流速均大于湾内,湾内东部流速略大于西部,计算的结果与实测值分布趋势也是一致的;夏季小潮期,涨急时刻流场表现为涨潮流从红海湾西侧进入,在湾内顺岸线逐渐向东,最后从东侧向东流出红海湾;落急时流向与涨急时相同。流速分布亦相同:湾外流速大于湾内流速,湾内西部流动弱于东部。计算的流场和实测流速分布表现为较一致的趋势。
图15~18为冬季大、小潮期间涨急与落急流场。大潮期间,海湾内的海流由东岸西北向流向湾内,西岸西南向流向湾外,海湾内的海流呈现逆时针旋转特征。小潮期间,海湾内的海流由西岸东北向流入湾内,东岸东南向流向湾外,海湾内的海流呈现顺时针旋转特征。
从以上夏季和冬季的计算结果来看,模拟计算所给出的红海湾流场与实测海流分布规律一致,本模型的计算方法是合理可行的。
5. 温排水数值模拟及其对环境的影响分析
5.1. 预测方案介绍
红海湾内广东太平岭核电项目拟建设2台百万千瓦级核电机组,循环水系统采用直流供水系统,以厂址周围海水作为冷却水源,循环冷却水量为130 m3/s,取排水温升7.83℃。工程取水工程及排放口位置见图19。
![](//html.hanspub.org/file/23990x27_hanspub.png)
![](//html.hanspub.org/file/23990x28_hanspub.png)
Figure 7. Verification curve of tidal level in summer
图7. 夏季潮位验证曲线
![](//html.hanspub.org/file/23990x29_hanspub.png)
![](//html.hanspub.org/file/23990x30_hanspub.png)
Figure 8. Verification curve of tidal level in winter
图8. 冬季潮位验证曲线
![](//html.hanspub.org/file/23990x34_hanspub.png)
Figure 9. Verification curve of flow velocity and flow direction in summer
图9. 夏季流速、流向验证曲线
![](//html.hanspub.org/file/23990x38_hanspub.png)
Figure 10. Verification curve of flow velocity and flow direction in winter
图10. 冬季流速、流向验证曲线
![](//html.hanspub.org/file/23990x39_hanspub.png)
Figure 11. Velocity vector of flood tide at spring tide in summer
图11. 夏季大潮涨潮流场图
![](//html.hanspub.org/file/23990x40_hanspub.png)
Figure 12. Velocity vector of ebb tide at spring tide in summer
图12. 夏季大潮落潮流场图
![](//html.hanspub.org/file/23990x41_hanspub.png)
Figure 13. Velocity vector of flood tide at neap tide in summer
图13. 夏季小潮涨潮流场图
![](//html.hanspub.org/file/23990x42_hanspub.png)
Figure 14. Velocity vector of ebb tide at neap tide in summer
图14. 夏季小潮涨潮流场图
![](//html.hanspub.org/file/23990x43_hanspub.png)
Figure 15. Velocity vector of flood tide at spring tide in winter
图15. 冬季大潮涨潮流场图
![](//html.hanspub.org/file/23990x44_hanspub.png)
Figure 16. Velocity vector of ebb tide at spring tide in winter
图16. 冬季大潮落潮流场图
![](//html.hanspub.org/file/23990x45_hanspub.png)
Figure 17. Velocity vector of flood tide at neap tide in winter
图17. 冬季小潮涨潮流场图
![](//html.hanspub.org/file/23990x46_hanspub.png)
Figure 18. Velocity vector of ebb tide at neap tide in winter
图18. 冬季小潮落潮流场图
![](//html.hanspub.org/file/23990x47_hanspub.png)
Figure 19. Geographical location of warm discharge and sensitive target map
图19. 温排口地理位置和敏感目标图
太平岭核电厂址东北方向约7 km为华润海丰电厂厂址,近期建设2台百万千瓦机组,冷却水量为72m3/s,取排水温升7.5℃。
5.2. 夏季温升计算结果与影响分析
图20为夏季全潮最大温升包络线图,表1为夏季温升最大影响分布。由图表可见,广东太平岭电厂温排水主要沿东北-西南方向扩散,温升范围主要围绕排水口附近的水域分布,夏季温排水造成的温升超过
1 ℃
的面积最大约为77.6 ha,对生态红线区域内造成的最大温升为
0.91 ℃
,夏季
1 ℃
最大温升包络位于近岸海域三类环境功能区内,未进入考洲洋和生态红线区。
海丰火电厂温排水主要沿岸边呈南-北向扩散,影响范围主要位于排放口附近的水域,夏季温排水造成的温升超过
1 ℃
的面积最大约为335.2 ha,温升超过
1 ℃
的范围对生态红线区域内造成的最大温升为0.76℃,夏季
1 ℃
最大温升包络位于近岸海域三类环境功能区内,未进入考洲洋和生态红线区。
5.3. 冬季温升计算结果与影响分析
图21为冬季温排水全潮最大温升范围包络线图,表2为冬季温升最大影响分布。由图表可见,广东太平岭电厂温排水扩散规律与夏季类似,主要沿东北-西南方向扩散,冬季温排水造成的温升超过
2 ℃
的面积最大约为27.0ha,对生态红线区域内造成的最大温升为
1.61 ℃
。冬季温升超过
3 ℃
的面积最大约为1.4ha,超过
4 ℃
的面积最大约为0.4 ha。冬季
2 ℃
最大温升包络位于近岸海域三类环境功能区内,未进入考洲洋和生态红线区。
海丰火电厂温排水主要沿岸边呈南-北向扩散,冬季温排水造成的温升超过2℃的面积最大约为455.9 ha,
![](Images/Table_Tmp.jpg)
Table 1. Thermal drainage discharge environmental impact analysis on water environment in summer
表1. 夏季电厂温排水最大温升对水环境的影响
![](//html.hanspub.org/file/23990x48_hanspub.png)
Figure 20. The influence range of the maximum temperature rise of the whole tide in summer
图20. 夏季全潮最大温升包络影响范围
![](Images/Table_Tmp.jpg)
Table 2. Thermal drainage discharge environmental impact analysis on water environment in winter
表2. 冬季电厂温排水最大温升对水环境的影响
对生态红线区域内造成的最大温升为1.95℃,冬季温升超过
3 ℃
的面积最大约为173.1 ha,超过
4 ℃
的面积最大约为128.6 ha。冬季2℃最大温升包络位于近岸海域三类环境功能区内,未进入考洲洋和生态红线区。
5.4. 温排水对海洋环境的影响
综上所述,红海湾内电厂温排水主要受风、海流等动力作用,在排水口附近海域输移扩散,形成一定区域的温升区域。温排水扩散区域及形态基本上与海流特征一致,太平岭电厂温升包络主要呈东北和西南向扩散特征,海丰电厂温升包络主要呈南北向顺岸扩散,向外海的扩散相对较弱。夏季
1 ℃
和冬季
2 ℃
最大温升包络位于近岸海域三类环境功能区内,未进入考洲洋和生态红线区。
![](//html.hanspub.org/file/23990x49_hanspub.png)
Figure 21. The influence range of the maximum temperature rise of the whole tide in winter
图21. 冬季全潮最大温升包络影响范围
6. 结论
通过对红海湾的海流观测资料进行分析,并进行了红海湾二维海流与温排水数值模拟,主要得到以下结论:
1) 红海湾海流符合南海海域的总体特征:潮差小,潮流弱,流态复杂,季节差异性大。红海湾海流除潮差影响因素外,还受到环流、局部地形引起的沿岸流等影响。
2) 红海湾二维海流数值模拟开边界由南海区域环流预测结果提供是合理的。南海区域环流预测考虑了季风、温度、盐度等综合因素的作用,计算结果与实际测量结果一致。
3) 红海湾在夏季和冬季以及不同潮型情况下的流态差异较大,测量及计算应分别考虑不同季节(至少夏季和冬季)、不同潮型的情况。
4) 红海湾内电厂温排水主要受风、海流等动力作用,温排水扩散区域及形态基本上与海流特征一致。太平岭电厂温升包络主要呈东北和西南向扩散特征,海丰电厂温升包络主要呈南北向顺岸扩散,向外海的扩散相对较弱。夏季
1 ℃
和冬季
2 ℃
最大温升包络位于近岸海域三类环境功能区内,未进入考洲洋和生态红线区。
资金项目
中央级公益性科研院所基本科研业务费专项资金项目(TKS160227, TKS160209)。