ELCIRC模型在番禺区台风暴潮数值模拟中的应用.pdfVIP

第8期 广东水利水电 No．8 2012年8月 GUANGDONG WATER RESOURCESAND HYDROPOWER Aug．2012 ELCIRC模型在番禺区台风暴潮数值模拟中的应用 黄武平 ，徐林春 ，张庭荣 (广东省水利水 电科学研究院，广东省水动力学应用研究重点实验室，广东广州 510635) 摘 要：基于水动力模型ELCIRC建立了珠江12'天文潮与风暴潮耦合模型，风场计算采用QSCAT／NCEP混合风场和台风 经验模型风场相加的方法，南中国海大范围模型提供外海边界，在对模型进行验证后 ，分别选取 了“0313号”台风 “杜 鹃”，“0601号”台风 “珍珠”和 “0814号”台风 “黑格比”，代表登陆型、北进型和西进型路径的台风进行番禺区台风暴潮 的数值模拟计算，计算结果表明能较准确地反映番禺区台风暴潮的增水值及出现时间。 关键词 ：ELCIRC；风暴潮 ；数值模拟 中图分类号：P731 文献标识码 ：B 文章编号：1008一O112(2012)08—00o4一o4 ELCIRC模型由美国俄勒冈州 OregonHealth&Sci— 区域剖分为2 ×2的网格 ，网格数为600×460。模型 enceUniversity的 OGISchoolofScience& Engineering 采用二维球坐标系下的潮波运动方程 ，开边界 由62个 开发，最初被用来模拟美国哥伦 比亚河 口海岸复杂动 分潮驱动。南中国海模型在数值计算 中采用球面坐标 ， 力环境，可用于计算流场、温度及盐度等。ELCIRC 以水位及流速分量作为方程的3个变量，计算采用交 模型采用了UnTrim模型的部分思想，它使用无结构网 替方向隐式格式 (ADI格式)。除以自然(陆域或岛屿) 格，在垂向上使用 z坐标，采用半隐格式的欧拉 一拉 岸线为陆边界外，开边界分别为：台湾海峡 由福建平 格朗 日有限体 有限差分法求解浅水方程，但 EL— 潭至台湾岛西岸 ，东部以吕宋海峡为控制边界 ，南部 CIRC和 UnTrim在处理切 向速度和温盐输运方面采用 以北纬 10。线为控制边界。 了不同的算法。此外，在该模型中还包含了多种紊流 本文珠江 口风暴潮模型基于 ELCIRC建立 ， 闭合格式以及对水气交换的详细描述。本文基于EL· 模型采用非结构混合网格，河道内采用正交曲线网格， CIRC模型建立的天文潮与风暴潮耦合模型，选取不同 在河道交汇处及 口门外复杂区域采用三角形网格。上 路径台风对番禺区的风暴潮进行数值模拟，研究各种 游边界分别取在潭江上游石咀、西江上游高要、北江 路径带来的风暴潮增水值及 出现时间，并与实测值进 上游石角、东江上游博罗、增江上游麒麟咀和广州片 行 比较。 区的老鸦 岗；模型外海开边界 由南中国海模型提供， 1 数学模型 取在珠江 口门外南海 30m等深线附近，包括整个伶仃 1．1 模型建立 洋和黄茅海。模型河道水深资料大部分为 1999～2005 QSCAT／NCEP混合风场是对高分辨率 的 Quik— 年测图，部分小河道为2009年的测图。模型采用非结 SCAT卫星散射计观测数据 (QSCAT风场数据 )和 构网格，河道 内采用贴体四边形网格，在河道交汇处 NCEP再分析数据在时间和空间上混合得到的覆盖全 及 口门外复杂区域采用三角形网格，网格的设计充分 球(88。S～88。N)10m高度处的风场…。QSCAT／NCEP 考虑到了贴合边界，网格数为99595个，网格节点为 混合风场资料始于 1999年7月，其时间分辨率为6h， 97583个，网格最小步长为20m左右。模型计算网格 空间分辨率为0．5。×0．5。，数据大概每 6个月更新一