刘林平，刘维杰, 2，孙志林

(1. 浙江大学 海洋学院，浙江 舟山 316021；2. 南京水利科学研究院 水文水资源与水利工程科学国家重点实验室，江苏 南京 210029)

海啸波通常由海底地震、火山爆发或者山体滑坡等引起，并通过沿岸爬坡过程侵袭陆地，造成极大的破坏，因此研究海啸波在近岸地区的传播以及爬坡过程具有重要的科学和工程意义。尽管实际海啸波通常以瞬态非周期性的形式出现，但孤立波被认为可模拟海啸首波的一些重要特征，因此学者通常采用孤立波进行海啸波的物模试验以及数值模拟研究[1-5]。自印度洋海啸之后，珊瑚礁在抵御海啸袭击中的作用引起了人们的重视，但众多灾后调查发现珊瑚礁对于海啸灾害的防御程度不一[6]。因此，为了更好地了解珊瑚礁在抵御海啸侵袭过程中的作用，学者们进行了相应的物理模型试验和数值模拟来探究珊瑚礁不同地形和水动力参数对海啸灾害的影响。

在物理模型试验方面，Quiroga等[7]通过大比尺水槽试验研究了礁床粗糙度对孤立波在礁坪上传播的影响；Yao等[8]在水槽试验中采用圆柱体阵列近似模拟了珊瑚礁大糙率表面，并提出了预测珊瑚礁海岸孤立波爬高的经验公式。相比物理模型试验，数值模拟能更便捷地调整珊瑚礁地形参数(如礁坪宽度、礁前斜坡、礁后斜坡等)和水动力参数(如礁坪水深、入射波高、珊瑚礁糙率等)。过去的几十年中，Boussinesq方程因其计算精度和效率较高的优势，在珊瑚礁波浪水动力数值模拟研究中应用最为广泛[9-14]。……