崖门水道悬沙输运及底床冲淤特性研究 Study on Suspended Sediment Transport and Bed Changes at Tidal Channel of Yamen.pdfVIP

第18卷4期 应用基础与工程科学学报 V01．18．No．4 2010 2010年8月 AND August JOURNALOFBASICSCIENCEENGINEERING 文章编号：1005-0930 2010 04-0637-08中图分类号：TVl48文献标识码：A l doi：10．3969／j．issn．1005-0930．2010．04．01 崖门水道悬沙输运及底床冲淤特性研究 刘金贵， 李瑞杰， 张义丰， 董德信 1．国家海洋环境预报中心，北京100081；2．河海大学海岸灾害及防护教育部重点实验室，江苏南京210098 摘要：针对感潮水道独特的水动力特征和泥沙输移情况，同时考虑径流、潮流等 多种动力因素共同作用，建立了三维悬沙输运及床面冲淤数学模型．水平方向采 用非结构三角形网格，垂向采用盯坐标，利用有限体积法离散积分形式三维浅 水守恒方程，对每个控制体单元分别进行水量和动量平衡计算．利用模型对崖门 水道的悬沙分布和输运情况进行了数值模拟，对泥沙浓度场的分布特征进行了 分析，并对底床冲淤变化进行了定量计算． 关键词：感潮水道；有限体积法；悬沙输运；底床冲淤 崖门位于新会县崖南镇，是银洲湖的人海口门，它与虎跳门均位于黄茅海湾的头部， 是珠江八大口门中最西边的一个I I门，潭江水域的径流主要通过银洲湖从崖门出海．已有 学者对珠江口所在海域的潮流运动及泥沙迁移进行了数值模拟¨引，但未能针对崖门水 道水域作系统的分析．李瑞杰旧1等结合崖门出海航道工程对回淤情况进行了分析，对不 同工程方案下的回淤强度进行了估算．崖门水道为重要的通航水道，同时由于崖门外两侧 围垦，海域水下地形的滩槽处于不断演变中，近年来深泓局部地段淤积加重，已显著碍航， 因此对该水域的悬沙输运以及河床冲淤进行分析研究具有重要的意义，以期解决影响崖 门水道进一步开发的泥沙问题，同时用于河口、海岸地区污染物扩散、盐水入侵及泥沙输 运等问题的研究． 基于有限体积法的泥沙数学模型在近十年才逐步发展起来．赵棣华等M1基于有限体 积法，采用黎曼近似解通量差分裂 FDS 格式计算各单元边的含沙量法向数值通量，应用 相关的悬移质、推移质河床变形计算公式计算了感潮河段泥沙冲淤变化；余明辉等∞1基 于三角形网格划分，连续方程采用有限体积法进行计算，保证了水量的守恒，建立了平面 二维非均匀沙数值模型；张杰等¨1基于有限体积法，采用SIMPLE系列算法求解控制方 程，建立了二维水沙和温度及浓度场数学模型．而三维有限体积悬沙输运及底床冲淤模型 的研究仍处于探索阶段． 本文利用海洋模型FVCOM[7o对崖门水道三维水动力场数值模拟，模型采用非结构三 角形网格和有限体积法．非结构网格较早应用于有限单元法，能与边界及水下地形拟合较 好，并利于边界条件的处理，便于控制网格密度和调整其适应性．而有限体积法在非结构 收稿日期：2008—10437；修订日期：2009-06-09 基金项目：交通部港口舟道泥沙工程重点试验室开放基金；河海大学自然科学基金科学研究重点项目 座机电话号码1 作者简介：刘金贵 1981～ ，男，博士．E．口aail：Liujg@nmefc．gov．cn 万方数据 638 应用基础与工程科学学报 三角形单元内离散积分形式的控制方程，对每个控制体单元分别进行水量和动量平衡计 算，因变量的积分守恒对任意一组控制体积得到满足，对整个计算区域，自然也保证了守 恒律．在此基础上，考虑悬沙输运过程及底床冲淤变化，对崖门水道悬沙分布及输运情况 作了进一步的研究，并对底床侵蚀和淤积作了计算分析． 1数学模型 1．1基本方程 矿坐标系下三维非恒定浅水方程为 篮+些+丝+aw：0 1 df dx do