基于第三代浅水波浪模型的海区泥沙骤淤数值预报及工程应用.pdfVIP

十、河流模拟技术 ’1491’ 基于第三代浅水波浪模型的海区 泥沙骤淤数值预报及工程应用 白玉川1 张胤祺1 纪 平2 吴朝安31 (1．天津大学建筑工程学院河流海岸工程泥沙研究室天津 300072； 北京 100044；3．香港大学机械工程系 2．中国水利水电科学研究院 香港) 摘要：简要介绍现有近岸波浪数值模型，对第三代浅水波浪数值预报模型进行研究和应用，通过数值 计算得到黄程山局部区域有效波高分布，与骤淤机理相结合，数值预报近岸区域骤淤量随时间的变化， 以及讨论了模型未来发展的趋势。 缓坡方程 第三代浅水波浪数值预报模型 骤淤 关键词：Boussinesq方程 近岸波浪对近岸海域工程的建设与规划、能源开发、海岸地形迁移变化有着重大的 影响，与人类实践活动息息相关，因此确定近岸区域波浪场是十分重要的。目前预报波 浪的近岸海浪平面二维数值模型主要分为经典的缓坡方程[1-3](椭圆型、抛物型、双曲型)、 J和能量平衡方程(SwAN)。 Boussinesq方程【4q 缓坡方程和Boussinesq方程由于模型控制方程和数值求解的不同，具有各自的特点 和局限性。缓坡方程是模拟近岸波浪传播变形的高效数学模型，可以考虑波浪由外海向 内海传播变形过程中所发生的波浪浅水变形、折射、绕射、反射等效应，同时也可以考 虑底摩擦和破碎引起的能量损耗。其局限性表现为：数值求解相对比较复杂，计算工作 量大，计算效率低；缓坡方程至少要求给定入射和侧边界条件，将导致计算结果不准确； 未能考虑自浪损耗、风能摄入等过程；不适用于地形变化急剧的复杂区域。Boussinesq 方程可以考虑波浪浅水变形、折射、绕射和反射等机理，其局限性主要表现为：Boussinesq 方程为二维非线性浅水波方程，对水深的适用范围有所限制[61；方程求解计算量大，只 适用于小范围区域；要求给定闭合边界条件，而实际工程中很难准确实现，极易导致计 算失真；在波能损耗和波与波之间相互作用的处理上比较不足。 SWAN【7q】波浪数值模型为适用于近岸、湖泊以及河口地区的第三代浅水波浪数值 预报模型，模型可较全面合理地考虑了波浪传播过程：多向不规则波的传播，波浪折射 和反射作用、浅水变形、障碍物的阻挡或部分传播作用。模型还可考虑波浪的成长和消 减过程：风成浪作用，白浪的削减作用，水深引起的破碎作用，海底摩擦作用，3～4个 波一波的非线性相互作用。边界条件给定也非常灵活，只要计算域的边界选择合适，即使 基金项目：国家自然科学基金资助项目50279030)。 作者简介：白玉川，1967年出生，男，山西神池人，博士，教授。 ．1492． 第六届全国泥沙基本理论研究学术讨论会论文集 只能确定迎浪边界条件也可获得较准确的计算结果。模型用不规定谱型的方向谱表示随 机波浪，使得海浪模拟更加真实准确。其局限性表现为：模型对波浪反射和绕射模拟不 佳，不适宜于模拟反射和绕射效果明显的近岸波浪。事实上，通过现场观测和试验证明， 当入射波方向谱不太窄时，不规则短峰波的绕射效应是相当小的。 1基本控制方程 在SWAN波浪数值模型中，以二维动谱密度N(o-，∞表示随机波浪。动谱密度为 能谱密度E(仃，0)与相对频率砣比，0为谱分布方向空间。因为考虑到水流作用，动 谱密度守恒，而能谱密度不守恒。 在Cartesian坐标系下，动谱平衡方程形式为 鲁Ⅳ+去GⅣ+杀。Ⅳ+杀GⅣ+嘉0CoN=孑S (·， 方程左边第一项表示动谱密度随时间的变化率；第二项和第三项表示动谱密度在地 理坐标空间上的传播；第四项表示由于水深和流场所导致的动谱密度随相对频率仃的变 化；第五项表示动谱密度随谱分布方向0的变化，亦表示因水深和流场产生的波浪折射。 方程右边的S为源汇项，表示因风能输入、底摩擦、白浪、波浪破碎等所产生的能量损 耗和波与波之间所产生的非线性相互作用影响。 2数值模型中物理过程的处理 2．1 风能输入 风能输入项根据Cavaleri和Malanott