李建秋 曾俏俏 苏跃涵 李恒翔 卢耀斌 骆海萍 刘广立

摘要：海水抽水蓄能電站上水库温度升高可能对下水库海域海水水温产生影响 ， 文章对珠海市大万山岛拟建的海水抽水蓄能电站温排水进行模拟分析 ， 结果表明 ， 海水抽水蓄能电站在设计水量12.3m3/s的条件下 ， 下水库排水8h后 ， 不会对所在海域海流流场产生显著影响;当排放温度为30.2℃时 ， 春季温排水影响的面积为1.76km2 ， 最大温升为1.35℃ ， 最大温升面积为0.006km2;夏

季影响面积约1.58 km2 ， 最大温升为0.7℃ ， 最大温升面积为0.013 km2;秋季影响面积约为1.66km2 ， 最大温升为0.8℃ ， 最大温升面积约为0.085km2;冬季影响面积约为1.95 km2 ， 最大温升为1.9℃ ， 最大温升面积约为0.030km2。周边海域最大温升随着排水温度增加逐渐增加 ，但对应

的最大温升面积变化较小。在排水温度为32.2℃时 ， 夏季最高温升约0.9℃ ， 冬季最高温升约2.3℃;排水温度为34.2℃ ， 夏季最高温升约1.15℃ ， 冬季最高温升约2.6℃。最大温升面积在0.019～0.027km2。排水流量从10.3m3/s逐渐增加到20.3m3/s， 周边海域的最大温升值增加 ， 但最大温升面积增幅较小 ， 当排水流量达到20.3m3/s 时 ， 最大温升值达到最大值2.25℃ ， 最大温升面积为0.007km2。文章的研究结果为未来大万山岛海水抽水蓄能电站的环境影响评价及周边海域生态管理提供了科学依据。

关键词：海水抽水蓄能电站;大万山岛;温排水;最大温升

中图分类号：P76;X55      文献标志码：A文章编号：1005-9857（2022）03-0009-07

Simulation and Ana;lysis on the Influence of Thermal Dischargeof Seawater Pumped Storage Power Station on SurroundingSea Area in Dawanshan Island of 'Zhuhai City

LI Jianqiu'，ZENG Qiaoqiao2，SU Yuehan3，LI Hengxiang'，LU Yaobin3，LUO Haiping3，LIU Guangli3

（1.CSG Power Generation Company，Guangzhou 510630.China;2.Guangdong Institute of water conservancy andPower Survey and Design Co.，Ltd，Guangzhou 510635，China;3.School of Environmental Science and Engineer-ing，Sun Yat-sen University，Guangzhou 510006，China;4. South China Sea Institute of Oceanology，Chinese A-cademy of Sciences ，Guangzhou 510301， China）

Abstract：The temperature rise of the upper reservoir of the seawater pumped storage power station mayaffect the seawater temperature of the lower reservoir.In this paper，the thermal discharge of the seawa-ter pumped storage power station proposed to be built on Dawanshan Island in Zhuhai City was simulatedand analyzed. The results showed that the seawater pumped storage power station would not have a sig-nificant effect on the local sea flow field under the condition of 12.3 m3/s within 8 h drainage from thelower reservoir. When the discharge temperature was 30.2 ℃，the affected area was 1.76 km ， and themaximum temperature rise was 1.35℃ with the maximum temperature rise area of 0.006 km2in spring.In summer，the affected area was 1.58 km2 ， and the maximum temperature rise was 0.7 ℃ with themaximum temperature rise area of 0.013 km" . In autumn，the affected area was 1.66 km ， and the max-imum temperature rise was 0.8 ℃ with the maximum temperature rise area of 0.085 km . In winter，theaffected area was 1.95 km2 ， and the maximum temperature rise was 1.9 ℃ with the maximum tempera-ture rise area of 0.030 km2 . The maximum temperature rise in the surrounding sea area increased gradur-ally with the increase of drainage temperature，but the corresponding maximum temperature rise areachanged little. When the drainage temperature was 32.2℃，the maximum temperature rise was 0.9℃and 2.3℃ in summer and winter，respectively.When the drainage temperature was 34.2℃，the maxi-mum temperature rise was 1.15℃ and 2.6℃ in summer and winter，respectively. The maximum temper-ature rise area was in a range of 0.019～0.027 km . With the drainage flow increasing from 10.3 to 20.3m3/s，the maximum temperature rise in the surrounding sea area increased but the maximum tempera-ture rise area was almost kept stable.The maximum temperature rise reached 2.25 ℃ with the maximumtemperature rise area of 0.007 km2 at the drainage flow of 20.3 m'/s. The results should provide a scien-tific basis for environmental impact assessment and ecological management of seawater pumped storagepower station in Dawanshan Island in the future.

Keywords;Seawater pumped storage power station，Dawanshan Island，Thermal discharge，Max-imum temperature rise

0 引言

21世纪是海洋的世纪 ，海洋战略在世界上所有沿海国家中都是重中之重[1-4]。我国海洋经济发展迅猛 ，随着海岛开发步伐的加快 ，海岛用电装机容量不足的问题日益凸显 ，成为制约海洋开发的重要“瓶颈”[5-6]。建设海水抽水蓄能电站 ，有利于增强海岛地区的电力保障 ，推动能源-海水淡化联产 ，为维护海岛安全和经济发展提供有力保障[7-11]。海水抽水蓄能电站是利用电网中负荷低谷时的剩余电能从海洋（下水库）向上水库抽海水 ，将电能转化为势能存储起来;当电力系统需要时 ，从上水库向海洋泄放海水发电 ，再将势能转化为电能的一种电站。我国海水抽水蓄能电站资源站点丰富[12-13] ，其中广东省珠海地区岛屿资源总量达到310MW[14-16]。

大万山岛[17-18]位于珠海市正南39 km 处 ， 面积约8.1km2 ， 经勘察海水抽水蓄能资源站点为大万山岛东北部的推船湾 ，上水库到下水库（海洋）的水平距离约560 m ，高差可达170 m ，距高比为3.3， 装机容量估计20MW 左右 ，有望成为较为理想的电站建设选址。大万山岛属亚热带海洋性气候 ， 光照辐射强 ，海水提升至上水库后 ， 上水库海水最高水位约110 m ，其接受太阳辐射能量后 ，水温升高的程度要高于下水库的海洋温度。研究表明 ， 太阳辐射是海洋热量的重要来源 ， 太阳辐射中的红外部分可被海水表层迅速吸收转化为热量 ， 而可见光部分能够穿透一定深度的海水 ， 如果海水较为清澈 ， 透明度高 ，其穿透深度可达60 80 m ， 由此引起的海水混合层平均水温升高可达1.5℃ 2.0℃ ， 而表层海水的水温升高可达到6℃ 7℃[19-20]。海洋中的鱼类对温度的变化极为敏感 ， 当海水水温变化超过0.1℃ ，都会引起鱼类行动的变化[21]。因此 ， 在电站的正常运行中 ， 升高的上水库排水形成的温排水 ， 进入下水库周边海域 ， 有可能对周边海域造成不良的影响 ， 形成热污染。《广东省海洋功能区划（2011—2020年）》将万山群岛及周边海域的功能区划为旅游娱乐、港口航运、农渔业 ， 执行二类海水水质标准 ，人为造成的海水温升夏季不超过当时当地1℃ ，其他季节不超过2℃。

对海岸发电厂的研究表明 ， 其利用海水作为冷却水运行中 ， 向海域排放的温排水 ， 可导致局部海域升温 ，对海洋物种生存造成恶劣影响。如温排水会干扰鱼类的产卵、孵化和生长 ， 影响海洋养殖业的正常运行[22-24];温排水也可对海洋底栖生物或者藻类种群产生不利影响。大亚湾核电站的温排水 ，其过高的温度可抑制浮游植物繁殖 ， 甚至致其死亡[25]。美国比斯开湾在温排水影响升温3℃ 4℃后 ，生物数量急剧减少。海洋温度上高 ， 不仅使溶解氧含量降低 ， 影响生物生存繁衍 ， 还会造成富营养化、水质恶化等问题[26-28]。丹麦水利研究所研发的 MIKE系列水动力软件 ， 广泛应用于海洋、河流、海湾等水域的模拟预测分析 ， 其数值模拟结果的可靠性已被世界公认。广西来宾电厂[29]、马来西亚沙巴电厂[30]、浙江苍南电厂[31]等温排水对附近海域的影响均采用 MIKE模型进行模拟 ， 分析潮型、流量及季节下的最大温升包络面积和取水口的温升变化 ， 为电厂的环境影响评价及管理提供依据[32-34]。

因此 ，本研究基于 MIKE模型 ，对珠海市大万山岛拟建海水抽水蓄能电站的温排水 ，开展其对周边海域水温影响的研究 ，为未来大万山岛海水抽水蓄能电站的环境影响评价及环境生态管理提供科学依据。

1 大万山岛拟建海水抽水蓄能电站位置及地质概况

大万山岛距珠海市直线距离约48km ， 西北距澳门35km ，东北距香港59km ， 岛上年平均气温为22.6℃ ， 气温变幅不大 ， 年平均温差为4.3℃。大万山岛所在地属于亚热带季风气候 ， 多年平均最大风速为17.3 m/s， 冬季多东北向风 ， 夏季多偏西南向风。大万山岛在地质上属北东向的万山隆起带（滨海大断裂） ， 为构造基本稳定区。拟建海水抽水蓄能电站位于山坳处 ， 库区所处位置三面环山 ， 为天然谷地 ，呈“V”型 ，远离岛上居民区。库区地貌属海岛地貌 ，库区地形起伏较大 ， 以构造侵蚀地貌为主;未发现较大滑坡、潜在不稳定地质体、泥石流及可溶岩地区喀斯特地貌等的分布。

2 模型建立

2.1 模型的边界条件及网格划分

如图1所示 ，研究模拟区域为大万山岛周围海域 ，平均水深10 28 m ，该模型基于非结构网格建立 ， 网格生成最大单元面积为20万 m2 ，最小允许角度为31°，最大节点数为50000， 然后对网格取迭代次数为10进行平滑处理 ，插值方法为自然邻点插值法。模型由11514个网格节点、20667个三角形网格单元构建而成。设置干水深、淹没水深和湿水深分别为0.005 m、0.05 m 和0.1 m ， 涡黏系数取为0.28，底床摩擦力取曼宁系数70 m1/3/s[35]。

海水抽水蓄能电站的运行参考淡水抽水蓄能电站的运行方式 ，按8h排水计算 ，取排水流量设置多个梯度（10.3 m3/s、12.3 m3/s、14.3 m3/s、16.3 m3/s、18.3 m3/s、20.3 m3/s）;背景溫度参考南海海面温度[36]及大万山岛当地温度值 ，取大万山岛海域春季水温为28.4℃ ， 夏季水温为29.3℃ ， 秋季水温为28.9℃ ，冬季水温为26.7℃。由于上下水库高差达170 m ，下水库的水被抽到上水库后 ，接收光辐射增强导致最终排水温度升高 ，但目前尚未有明确可参考的温升范围[29]。参考太阳辐射引起的海水水温升高的情况 ， 本研究在初步选取电站温排水温度30.2℃进行模拟计算后 ， 对32.2℃～40.2℃的排水温度下 ，最大温升值及最大温升面积进行了对比分析。设定的排水温度为40.2℃时 ， 大万山岛海域的水温变化范围为0.9℃～3.5℃。

2.2 模型验证

潮位验证数据实测值采用国家海洋信息中心东澳岛监测站（113°41'33E，22°0'39N）和大万山监测站（113°42'40E，21°55'46N）公开发布的潮汐表数据 ，分别选取春季、夏季、秋季、冬季不同时间段24h内的潮位与模型模拟值进行比较。验证结果如图2和图3所示 ，可以看出模拟结果与东澳岛、大万山岛监测站的实测数据吻合较好 ，平均相对误差低于10% ， 由此认为所建立的模型是合理的 ，可以满足预测要求。

3 结果与讨论

3.1 电站排水对周边海域潮流场的影响

按电站排水水量12.3 m3/s，排水与周边海域无温差计算 ，对拟建电站工程实施前后海域潮流场进行模拟分析 ，结果如图4所示。由图可知 ， 电站排水前 ，下水库附近的潮流自西北沿大万山岛顺岸流动。当电站开始排水后 ， 排水在入海口处流速较快 ，汇入外海时形成小型旋涡 ， 之后逐渐向外海扩散 ， 由于电站排水流量较小 ， 流速迅速降低 ， 对潮流方向无明显影响。因此抽水蓄能电站排水不会对所在海域海流的流场产生显著影响。

3.2 不同季节电站排水对周边海域温度的影响分析

不同季节下 ， 电站排放流量为12.3m3/s，排放温度30.2℃ ， 连续排水8h后 ， 周边海域的温升结果模拟计算如图5所示。与周边海域海水温差超过0.1℃的区域 ， 即认为受到温排水影响的区域 ，计算温排水的影响面积。由图5可知 ， 春季温排水受潮流作用 ， 温升场逐渐向东南侧移动 ， 最远延伸至2.09km ， 影响面积约1.76 km2;最大温升为1.35℃ ，最大温升面积约为0.006 km2 ， 最大温升出现在距离排放口约0.5km 的海域。温排水主要呈南北方向移動 ， 东西方向移动距离较小。夏季温升场最远延伸至1.86 km ， 影响面积约1.58km2;最大温升为0.7℃ ， 最大温升面积约0.013km2 ，最大温升出现在距离排放口约0.6km 的海域。秋季温升场主要向排水口南侧移动 ， 最远延伸2 km ， 影响面积约1.66km2;最大温升为0.8℃ ，最大温升面积约为0.085km2。冬季温升场主要向西北方向移动 ，最远延伸至2.9km 的海域 ， 影响面积约1.95 km2;最大温升为1.9℃ ， 最大温升面积约为0.030km2。

研究表明 ，火力发电厂温排水对周围水域造成的温升可达6.11℃ ，这是由于火力发电站的装机容量大 ，且尾水温度较高 ， 温差可达到13.3℃ ， 因此对周围水域造成的温升大[29]。而海水抽水蓄能电站的排水温度按30.2℃计算时 ， 与周边海域的温差较低 ，夏季仅相差0.9℃ ， 因此对周围海域造成的温升较小。

3.3 电站温排水温度变化和流量对周边海域最大温升的影响分析

为进一步考察电站温排水对周边海域的影响 ， 在电站排水量为12.3 m3/s， 排水温度分别达到32.2℃、34.2℃的条件下 ，模拟计算夏季和冬季时周边海域温升的结果 ， 如图6和图7所示。在排水温度为32.2℃时 ，夏季最高温升约0.9℃ ， 出现在距离排放口约0.75 km的海域 ， 最大影响面积约为0.027km2;温排水总影响面积约为1.6km2。冬季最高温升约2.3℃ ， 最大温升面积约0.019 km2 ， 出现在距离排放口约0.37km 的海域 ， 温排水总影响面积约为1.81km2。当排水温度为34.2℃ ，夏季最高温升约1.15℃ ，最大温升面积约0.012km2 ， 出现在距离排放口约0.75km 的海域 ， 温排水总影响面积约为1.69 km2。冬季最高温升约2.6℃ ， 最大温升面积约0.02km2 ， 出现在距离排放口约0.37km 的海域 ，温排水总影响面积约为1.83km2。

选取电站温排水温度30.2℃ ， 周边海域海水水温为28.4℃的条件下 ，研究不同电站排水流量对海域温升的影响 ，结果如图8所示。随着排水流量从10.3 m3/s逐渐增加到20.3 m3/s， 周边海域的最大温升值增加 ， 但最大温升面积增幅较小 ， 当排水流量达到20.3 m3/s 时 ， 最大温升值达到最大值2.25℃ ，相应的温升面积为0.007 km2。这表明 ， 即使电站排水流量达到20.3 m3/s，也难以对周边海域潮流的运动产生显著影响 ， 导致最大温升值虽然增加 ，但最大温升面积并没有显著增大。可能是由于周边海域潮流影响起主要作用 ， 电站排水流量相对较小 ， 即使排水的温度与海水的温度温差较大 ， 进入海域中也可被很快稀释降温。

对大万山岛周边海域的海洋生物初步调查发现 ，鱼类主要有绯鲤、鲻科、多鳞和舌鳎等。绯鲤一般栖息于水深20～110 m ， 生殖期为2—4月（春季） ，在广东省上川岛和下川岛附近一带水深8～35 m处有密集的生殖群 ， 该区域距离大万山岛约140km 。鲻科鱼类以底藻和碎屑为食 ， 耐寒冷 ， 是亚热带和热带地区的沿岸经济性鱼类。多鳞为近海小型鱼类 ， 喜栖息在沙质海底中下层 ， 7—10月（夏秋季）为其活跃觅食期。舌鳎栖息于浅海底层泥沙底质海区 ，繁殖期为3—10月（春夏秋季） ，适宜温度为22℃～23℃。因此 ，海水抽水蓄能电站在运行中 ，上水库排出的温排水对周边海域鱼类造成的影响可能各不相同 ，绯鲤的影响较小 ， 而鲻科、多鳞和舌鳎类鱼类因在近岸海底觅食繁殖 ， 上层海水温度升高可能对此类鱼类的生长有一定影响;若在周边海域温升达到2.25℃的条件下 ， 则需要密切关注鲻科、多鳞和舌鳎类鱼类的生长繁殖 ，对海水抽水蓄能电站下水库排水的方式进行改进 ， 强化其扩散过程 ， 降低温排水对周边海域的温升效应。

4 结论

（1）海水抽水蓄能电站在设计水量为12.3 m3/s 的条件下 ，下水库排水8h后 ，不会对所在海域海流流场产生显著影响。

（2）电站排放流量为12.3 m3/s， 排放温度为30.2℃ ，连续排水8h后 ， 春季温排水影响的面积为1.76 km2 ， 最大温升为1.35℃ ， 最大温升面积为0.006km2;夏季影响面积约1.58 km2 ， 最大温升为0.7℃ ，最大温升面积为0.013km2;秋季影响面积约1.66 km2 ， 最大温升为0.8℃ ， 最大温升面积约为0.085km2;冬季影响面积约1.95 km2 ， 最大温升为1.9℃ ，最大温升面积约为0.030km2。

（3）电站排放流量为12.3 m3/s时 ，周边海域最大温升随着排水温度增加而逐渐增加 ，但对应的最大温升面积变化较小。在排水温度为32.2℃时 ， 夏季最高温升约0.9℃ ， 冬季最高温升约2.3℃;排水温度为34.2℃ ， 夏季最高温升约1.15℃ ， 冬季最高温升约2.6℃。最大温升面积在0.019～0.027 km2之间。排水流量从10.3m3/s 逐渐增加到20.3 m3/s时 ，周边海域的最大温升值增加 ， 但最大温升面积增幅较小 ， 当排水流量达到20.3 m3/s时 ， 最大温升值达到最大值2.25℃ ， 最大温升面积为0.007km2。

参考文献

[1]   李琳梅.发展海水淡化事业促进海洋强国建设[J].海洋开发与管理 ， 2012， 29（12）：59-61.

[2]   毕建国 ， 段志霞.我国海洋渔业生态环境污染及治理对策[J].中国渔业经济 ， 2008（2）：16-21.

[3]   CAVALLO M ， BORJA ?， ELLIOTT M ， etal. Impedimentsto achieving integrated marine management across borders： ThecaseoftheEUMarineStrategyFrameworkDirective[J].MarinePolicy， 2019， 103：68-73.

[4] 吴凯 ， 卢布.中国海洋产业结构的系统分析与海洋渔业的可持续发展[J].中国农学通报 ， 2007（1）：367-370.

[5] 刘哲.海水抽水蓄能技术在海岛供电中的应用探讨[J].电力与能源 ， 2017（5）：593-597.

[6]WANG Z，LIN X，TONG N，et al. Optimal planning of a l00%renewable energy island supply system based on the integration of aconcentrating solar power plant and desalination units[J]. Interna-tional Journal of Electrical Power 8. Energy Systems，2020 ，，l17;105707.

[7]钱钢粮.我国海水抽水蓄能电站站点资源综述[J].水电与抽水蓄能 ， 2017， 3（5）：1-6.

[8]石文辉 ， 查浩 ， 罗魁 ， 等.我国海水抽蓄电站发展初探[J].中国能源 ， 2015（12）：36-40.

[9]谭雅倩 ， 周学志 ， 徐玉杰 ， 等.海水抽水蓄能技术发展现状及应用前景[J].储能科学与技术 ， 2017， 6（1）：35-42.

[10]张皓天 ， 吴世东 ， 芮德繁 ， 等.我国海水抽水蓄能电站示范项目选择及开发研究[J].水电与抽水蓄能 ， 2017， 3（5）：7-10.

[11]侯平利 ， 柴建峰 ， 陈雄 ， 等.单库海水抽水蓄能电站技术研究[J].水电与抽水蓄能 ， 2017， 3（5）：19-21.

[12]張东.我国南部沿海地区海水抽水蓄能资源初步评价[J].低碳世界 ， 2015（1）：69-70.

[13]刘万久 ， 张皓 ， 黎俞琳 ， 等.浅谈海水抽水蓄能电站在中国的发展前景[J].四川水泥 ， 2018（3）：295-296.

[14]DUARTECM ， LOSADA IJ， HENDRIKSIE， etal. Theroleofcoastalplantcommunitiesforclimatechangemitigationandad- aptation[J]. NatureClimateChange， 2013， 3（11）：961-968.

[15]KATSAPRAKAKIS D A， CHRISTAKIS D G， STEFANAKIS I， etal. Technicaldetailsregarding the design， the construction and the operation ofseawater pumped storage systems[J].Energy， 2013， 55（1）：619-630.

[16]MORGAN LK， WERNER A D. Seawaterintrusionvulnera- bilityindicators for freshwater lenses in strip islands[J].JournalofHydrology， 2014， 508：322-327.

[17]廖静 ， 方琼玟.珠海万山岛：修复生态探索海岛休闲游[J].海洋与渔业 ， 2017（8）：34-35.

[18]余志 ， 蒋念东 ， 游亚戈.大万山岸式振荡水柱波力电站的输出功率[J].海洋工程 ， 1996（2）：78-83.

[19]刘志亮.穿透性太阳辐射和盐度作用对混合层深度的影响[D].青岛：中国海洋大学 ， 2004.

[20]刘童童.短波辐射传输模式对上层海温及混合的影响研究[D].厦门：厦门大学 ， 2017.

[21]李雪渡.海水温度与渔场之间的关系[J].海洋学报（中文版） ， 1982（1）：103-113.

[22]於凡 ， 张永兴 ， 曹颖.全球变暖背景下核电站温排水对海洋生态系统的影响[J].海洋湖沼通报 ， 2010（2）：155-161.

[23]  JAVED M ， AHMAD I， USMANIN， et al. Bioacumulation，oxidativestresandgenotoxicityin fish （Channa punctatus） ex- posedto a thermalpower plant efluent[J]. Ecotoxicology andEnvironmentalSafety， 2016， 127：163-169.

[24]  ANSAHYB， FRIMPONGEA， AMISAHS. BiologicalAses-mentofAquacultureEfects on Efluent-Receiving Streams in GhanaUsingStructuralandFunctionalCompositionofFish and MacroinvertebrateAsemblages[J]. EnvironmentalManagement，2012， 50（1）：166-180.

[25] 谢福武 ， 刘华雪 ， 黄洪辉 ， 等.大亚湾浮游植物粒级结构对温排水和营养盐输入的响应[J].热带海洋学报 ， 2018， 37（3）：55-64.

[26] 常小军.胶州湾温排水数值模拟[D].青岛：中国海洋大学 ， 2011.

[27]  PENK M R， WILLIAMSM A. Thermalefluentsfrom powerplantsboostperformanceoftheinvasiveclamCorbiculaflumineainIreland'slargestriver[J].ScienceofTheTotalEnvironment， 2019， 693（133546）.

[28]  JAVED M ， USMANIN. Stresresponseofbiomolecules（carbo-hydrate， protein andlipidprofiles） in fishChannapunctatusin- habiting riverpolutedbyThermalPowerPlantefluent[J]. SaudiJournalofBiologicalSciences， 2015， 22（2）：237-242.

[29] 王慶改 ， 戴文楠 ， 赵晓宏 ， 等.基于 Mike21FM 的来宾电厂扩建工程温排水数值模拟研究[J].环境科学研究 ， 2009， 22（3）：332-336.

[30] 严冰 ， 张娜 ， 赵洪波.沿岸往复流海岸电厂平面二维温排水数值模拟研究[J].水道港口 ， 2011， 32（4）：291-296.

[31] 张舒羽 ， 黄世昌 ， 韩海骞.浙江苍南电厂冷却水温排放的数值模拟[J].海洋学研究 ， 2009， 27（3）：61-66.

[32] 谢晓倩. MIKE21 FM 在洪水风险图制作中的应用研究[D].大连：大连理工大学 ， 2016.

[33] 袁雄燕 ， 徐德龙.丹麦 MIKE21模型在桥渡壅水计算中的应用研究[J].人民长江 ， 2006（4）：31-52.

[34] 张志林 ， 贾艾晨.基于 MIKE21 FM 模型的河道流场图绘制[J].东北水利水电 ， 2016， 34（10）：35-70.

[35] 张大茹.基于 Mike21FM 的山区小流域涉水工程防洪影响研究[D].北京：中国水利水电科学研究院 ， 2015.

[36] 孙伟富 ， 张杰 ， 孟俊敏 ，等.中国南海及邻近海域 SST时空分布和变化特征分析[J].海洋科学进展 ， 2018， 36（3）：402-411.