黄茅海围垦工程潮流泥沙变化数值模拟.pdfVIP

泥沙研究 2008年4月 JoumalofSedimentResearch 第2期 黄茅海围垦工程潮流泥沙变化数值模拟 汤立群1，梁建林2，刘大滨1 (1．中国水利水电科学研究院。北京100044；2．黄河水利职业技术学院，河南开封475001) 摘要：简述了黄茅海海域平面二维潮流泥沙数学模型的基本方程、基本问题的处理、初边界条件和网格划分。 中水大潮情况下，着重对选定的治导线围垦工程方案ElW2进行了工程前后潮流泥沙动力条件变化模拟计算 分析，结果表明：围垦工程实施后，对海区尤其是口门内河道沿程不会产生防洪压力，对崖门水道和虎跳门水 道排涝的影响不大，围垦工程不会引起水流动力轴线的明显变化，对围堤和口门内河岸的稳定不会构成威胁。 关键词：围垦工程；潮流泥沙；数值模拟；方案计算 中图分类号：Tvl48 文献标识码：A 文章编号：0468．155x12008)02．0009—0r7 l。前言 黄茅海是珠江三角洲八大口门中最西面的崖门和 虎跳门水道人海处的河口湾，崖门上游为潭江，人海段 为银洲湖，银洲湖河段又汇合了由西江分流的江门河 和虎坑水道水量。位于东经113咖7一11305l7，北纬2lo 50’。22。13 7，北起湾顶崖门、虎跳门出口，宽1．8km，下 延约30km处，有大芒岛、三角山岛、南水尾的湾口，宽 约20km，浅海水域面积为403km2。再下延10km，外侧 有大襟岛、荷包岛和高栏岛，并将湾口分隔成三个口 门，即大襟岛以西至台山市陆域，为西口门；大襟岛与 荷包岛间峡口为中口门，称为西槽；荷包岛与高栏岛间 峡口为东口门，称为东槽。湾顶至湾口长约40km，中 图1 黄茅海河口湾示意图 湾宽约17km，三个湾口总宽15km，海湾面积537．7km2， 11le of sea 如图l。 Fig．1 estuarybayHuangmaoIlai 在黄茅海治理规划与围海造田工程中，根据滩槽 水下地形的格局，在东西滩地拟定了5条治导线，E。、E2、w。、W2和W3(见图1)。根据造价低、投入少、 效果佳、对社会生态环境影响小的原则，综合比较分析后推荐ElW2方案作为黄茅海海区治导线规划方 案。在E1W2推荐方案下，黄茅海海域的潮流泥沙动力条件如何变化，是否会对环境产生不利影响，特 别是潮位的变化对于崖门水道和虎跳门水道的防洪排涝产生怎样的影响，应该进行重点模拟与预测分 析。本文首先对平面二维潮流泥沙数学模型进行验证，其次，利用验证后的数学模型进行围垦工程方案 收稿日期：2007．09．“ 基金项目：中国水利水电科学研究院专项基金(泥集05KY01)；教育部回国人员科研启动基金 作者简介：汤立群(1962一)，男，教授级高工，博士后，主要从事河口海岸动力过程和泥沙输移、土壤侵蚀及水土保持研究。 E-mil：tmIglq@iwllr．咖 9 计算，重点对比分析方案前后潮位、纳潮量、流速流向、含沙量以及冲淤地形等的变化，最后，对推荐方案 作出评价。 2数学模型 2．1正交曲线坐标系下水流泥沙基本方程 平面二维正交曲线坐标下沿水深平均的零方程紊流模型控制方程为 连续方程 。r (1)”7 萼+去『掣I掣1：o a善 J一” at’CfC口L a叩 动量方程 ， 掣+矗【委c