热带海洋学报 ›› 2016, Vol. 35 ›› Issue (2): 10-18.doi: 10.11978/2015034CSTR: 32234.14.2015034

• 海洋水文学 • 上一篇    下一篇

基于FVCOM的虎门水道盐水入侵三维数值模拟

刘祖发1, 2, 关帅1, 2, 张淦濠1, 2, 丁波1, 2, 林颖妍1, 2, 查悉妮1, 2   

  1. 1. 中山大学水资源与环境研究中心, 广东 广州 510275; 2. 华南地区水循环与水安全广东省普通高校重点实验室, 广东 广州 510275
  • 收稿日期:2015-03-05 修回日期:2015-08-11 出版日期:2016-02-29 发布日期:2016-02-29
  • 作者简介:刘祖发(1961—), 男, 河南省潢川市人, 副教授, 博士, 主要从事水文与水资源研究。E-mail: eeslzf@mail.sysu.edu.cn
  • 基金资助:
    国家自然科学基金项目(41301627)

Three-dimensional numerical simulation of saltwater intrusion into the Humen estuary based on FVCOM

LIU Zufa1, 2, GUAN Shuai1, 2, ZHANG Ganhao1, 2, DING Bo1, 2, LIN Yingyan1, 2, ZHA Xini1, 2   

  1. 1. Center for Water Resources and Environment, Sun Yat-sen University, Guangzhou 510275, China; 2. Key Laboratory of Water Cycle and Water Security in Southern China of Guangdong Higher Education Institutes, Guangzhou 510275, China
  • Received:2015-03-05 Revised:2015-08-11 Online:2016-02-29 Published:2016-02-29
  • Supported by:
    The National Natural Science Foundation of China (41301627)

摘要: 采用不规则网格有限体积近岸海洋三维模型FVCOM (finite volume coastal ocean model), 对2001年虎门水道至伶仃洋的盐水入侵现象进行数值模拟, 并采用多个站的实测潮位、盐度对模型进行了率定和验证, 然后对其水动力特性、盐度的时空变化规律及其垂向分层现象进行分析。结果表明, 水位的模型预测值与实测值之间的相关系数达到0.93以上, 绝对平均误差小于0.1m; 盐度模拟的相对误差能控制在10%以内; 大虎站表底层的流向较为一致, 盐度表现出明显的分层现象; 而漳澎站多次出现表底层流向相反的情况, 盐度的垂向混合较充分; 文章构建的模型能够较好地模拟盐水入侵的最大距离, 虎门水道大潮时入侵距离和强度均大于小潮时, 盐水楔活动受潮流强弱变化影响显著。

关键词: 盐水入侵, 数值模拟, 有限体积近岸海洋三维模型(FVCOM), 虎门

Abstract: An unstructured grid finite-volume coastal ocean model (FVCOM) was applied to simulate the saltwater intrusion from Lingdingyang into the Humen estuary. Measurements of water level and salinity at multiple stations obtained in February 2001 were used for model calibration and verification, aiming at the hydrodynamic characteristics and spatiotemporal variation of salinity. The results showed that 1) the correlation coefficients between the predicted and measured water level exceeded 0.93 with a mean absolute bias less than 0.1 m, while the relative error of the simulated salinity was within 10%; 2) at Station Dahu, the flow direction in the surface layer was consistent with that in the bottom layer most of the time, but there was apparent salinity gradient in the vertical. However, for Station Zhangpeng, opposite flow directions between the surface and bottom layers occurred repeatedly while the salinity was more uniform vertically; and 3) the maximum distance of saltwater intrusion was well simulated. The intrusion distance as well as the strength at spring tide was greater than that at neap tide in the Humen estuary, and saltwater wedge movement was significantly affected by the tidal strength.

Key words: saltwater intrusion, numerical simulation, finite-volume coastal ocean model (FVCOM), Humen