N波传播与爬高的数值模拟及其特征研究.pdfVIP

N 波传播与爬高的数值模拟及其特征研究1) * * *,2) 赵 曦，王本龙 ，刘 桦 * （上海交通大学工程力学系，上海 200240 ） 摘要：大多数以断层错动为特征的海底地震激发的海啸波为N 波。N 波是指由一个波峰和一个波谷组成的形状如 “N ”的海洋长波，其在海洋中传播以及近岸爬高与以往作为海啸波模型的孤立波有所不同。本文基于完全非线 性高阶色散性的Boussinesq 方程，模拟了N 波在平底上长距离传播及在缓坡上浅化，并对N 波在斜坡上爬高速度 和能量转化问题进行了细致的研究。 关键词：海啸；N 波；爬高；Boussinesq 方程 引 言 最初的海啸研究认为海啸波的形状为孤立波，这可以解释实际观测到的海啸波经过跨洋传播后能量衰 减很小和波形基本不变的现象。但是从更多的地震海啸(如1993年9月1 日的尼加拉瓜海啸，1994年爪哇 岛海啸等)实际观测中发现多数海啸波是由一个大波峰和一个大波谷组成，形状像 “N”，Tadepalli Synolakis将这类波称为N-wave(N波)，并最先给出了N波的一般表达式[1]、[2] 。从许多关于海啸的报道中 可以发现，在海啸到来之前多数目击者称曾看到海岸线后退的现象，但经典的孤立波理论无法解释这一现 [3] 象。对于深海地震来说，特别是沿着断层面上下错动为特征的 “倾滑型”地震，形成的海面波动为N波 ， 因此在多数情况下，N波比孤立波更能很好地模拟海啸波的形态与特征。 本文基于完全非线性高阶色散性的Boussinesq 水波模型，模拟N 波在平底地形上长距离传播，研究 了频率色散性的累计效果，并得到海啸波跨洋传播以及在缓坡上浅化时的波高变化规律。研究了N 波在到 达近岸时，岸线移动速度的变化问题。通过Boussinesq 模型重构水波速度场，计算N 波特征水层面上的 速度，由水面变形和流场内的速度计算波动势能和动能分布，研究爬高时的能量转化问题。 1 控制方程 考虑不可压流体的无旋流动，自由表面运动学边界条件和动力学边界条件分别为：    w  u 0 (1) t x  1 2 2 2 U g   [U w (1 )] 0 (2) t x x x 2 ~ 其中，(x,t ) 为水面波动，t 为时间，g 为重力加速度，U 为辅助速度变量，由下式表示： 1) 国家自然科学基金 和博士点基金 （20060248046 ）资助项目 2) E-mail ：hliu@sjtu.edu.cn    U u wx (3) 底面边