应用MIKE21之BW模型进行日照港波浪数值模拟研究.docVIP

应用MIKE21之BW模型进行日照港波浪数值模拟研究 应用ＭＩＫＥ２１之ＢＷ模型进行日照港波浪数值模拟研究　 张娜 郭科　 （天津港湾工程研究所，天津，　３００２２２）　 摘要：ＭＩＫＥ２１是由丹麦ＤＨＩ公司在２０多年来大量工程试验基础上开发的平面二维河口与海岸模拟软件，其中ＢＷ模型是用于波浪数值模拟的一个重要模块，也是该软件发展最为成熟和完善的一个模块，在世界各国波浪数学模型试验中已有十分广泛的应用。本文结合日照港工作船防浪措施工程波浪数值模拟研究，系统阐述了ＭＩＫＥ２１之ＢＷ模型在实际工程试验中的适用性、具体方法以及使用过程中的一些关键性问题。尤其针对复杂地形及建筑物环境下的边界条件处理及计算参数的选择进行了深入的探讨。目前试验成果已为设计采用，受到业主及设计单位好评。　 关键词：ＭＩＫＥ２１；ＢＷ；波浪；数值模拟；日照港　　 　 １．ＭＩＫＥ２１之ＢＷ（Ｂｏｕｓｓｉｎｅｓｑ　Ｗａｖｅ）模块简介：　 该模型通过求解沿垂向积分的Ｂｏｕｓｓｉｎｅｓｑ方程获得沿水深平均的流速、水位变化以及波高等物理量。ＭＩＫＥ２１－ＢＷ模型可用来模拟任意地形或平面布置条件下的波浪运动过程，包括波浪的绕射、折射、反射、浅水变形以及它们之间的相互作用。在模拟反射时，既可以模拟完全反射、也可以模拟部分反射情况。它具有在指定区域吸收（即消除）波浪的功能，从而可以消除不希望的二次反射的影响。此外，它还具有模拟内部生波的功能。　 ＭＩＫＥ２１－ＢＷ不同于一般的波浪数学模型，它是在指定时段和指定区域内模拟波浪运动的物理过程，然后根据记录下来的水面波动过程，统计波高等特征值。从这个意义上讲，用ＭＩＫＥ２１－ＢＷ进行数学模型试验与通常的物理模型试验方法及步骤完全相同，所以它对现场波况的描述更加真实可靠。　 与此同时，ＭＩＫＥ２１之ＢＷ模块同样存在着一些不足：　 （１）关于波浪破碎一项，２００５版以前没有加入该参数，进行近岸浅水区波浪数值推算或者模拟直立式建筑物附近波浪反射现象较明显水域时，由于无法模拟波浪破碎，常常会使结果失真（大多数情况会导致结果偏大），或者致使模型溢出，无法模拟完全。　 （２）２００５升级版虽然加入了波浪破碎参数项，但使用起来比较繁琐，模拟破碎区域的网格需要划分得非常细，应为原模型网格大小的五分之一，故模拟时间也会相应延长，大大影响了计算效率。　 （３）航道的边坡开挖功能不健全，需要人工计算和手工输入，这样浪费了大量的人力 和时间，建议开发一下此项功能，给定航道边坡比就能沿航道自动进行边坡插值。　 ＭＩＫＥ２１之ＢＷ模型建立主要包括三个主要部分：（１）基本参数：需要设置地形数据、方程类型以及模拟时间等信息。（２）校准参数：涉及到诸如水位之类的边界条件、以及类似底摩阻之类的物理参数的选取等问题。（３）结果输出：可预置输出结果类型，以及数据采集时间和步长等信息。　 ２．基本方程　连续方程：?ξ?p?q++=0　?t?x?y p2q2pg+??p2???p.q???2p?2p??p??ξhhｘ方向动量方程：+?　?+?E??+gh?+?y???x+?y???t?x?hh?xc?????? =1??3p?3q?Dh?+??x2?t?x?y?t3?? ｙ方向动量方程：p2q2qg2+2???2q?2q??q??q2???p.q??ξhhh+?E?+?x?h?+gh?y+2??x2+?y2??　?t?yhc?????? 1??3q?3p?=Dh?+3??y?t?x?y?t?? 式中：　 ｘ，ｙ　－－－　　水平座标（ｍ）　 ｔ 　－－－　　时间（ｓ）　 ξ －－－　　高出平均水位的水面高度（ｍ）　 p?q　－－－　　ｘ，ｙ方向流量密度　 ｈ 　－－－　　水深（ｍ）　 Ｄ －－－　　平均水深（ｍ）　 ｃ 　－－－　　谢才阻力系数，　ｃ＝Ｍ?ｈ１／６　 Ｍ －－－　　曼宁系数，缺省值为　３２　（ｍ１／３／ｓ）　 Ｅ 　－－－　　紊动“涡粘”系数　 Ｇ 　－－－　　重力加速度　 　 ３．计算实例　 ３．１　工程概况　 为了给日照港的工作船停靠泊位提供良好的 防浪条件，在港内考虑了三个可供工作船停靠的 位置，它们分别是在东港区泊位处（ＡＢＣＤ所在点 位处），称为方案一；二期集装箱码头南端（ＥＦＧ图１　　日照港工作船防浪措施工程总平面及测点布置所在点位处），称为方案二；西港区南护岸的南端（ＨＪＫ所