深圳河口防洪能力数值研究.docx

１０００－７７０９２０１１１０－００２２－０４深圳河口防洪能力数值研究胡明军，毛献忠（清华大学深圳研究生院环境工程与管理研究中心，广东深圳５１８０５５）摘要：以深圳河为例，针对下游泥沙淤积严重、暴雨形成洪峰易导致洪灾问题，应用Ｄｅｌｆｔ３Ｄ模型区域分解技术构建深圳河湾二维水动力耦合模型，并基于实测洪水过程验证模型，模拟计算了防洪能力与清淤效果。关键词：深圳河口；洪水；防洪能力；清淤；数值模拟中图分类号：ＴＶ１４８文献标志码：Ｂ深圳河为深圳与香港的界河，由于受区域气候的影响多次发生大暴雨，并汇流形成洪峰，导致沿河低地洪灾［１，２］。深圳河泥沙淤积严重［３］，河床抬高及河口围垦［４］、填海工程［５］均给防洪带来隐患。深圳河作为城市感潮河流，与长江三角洲、珠江三角洲等河流具有相似的特点，研究深圳河的防洪可为类似情况提供借鉴。数学模型是研究洪水过程常用的方法［６，７］，Ｄｅｌｆｔ３Ｄ是一套功能强大的软件包，主要应用于地表水环境模拟，框架灵活，可模拟二维、三维的水流、波浪、水质、生态、泥沙运输等及各过程之间的相互作用，是最先进的水动力—水质模型之一［８］。鉴此，本文采用Ｄｅｌｆｔ３Ｄ模型模拟计算现有河道防洪能力，并进一步讨论了清淤效果。图１深圳河湾水系示意图Ｆｉｇ．１ＳｃｈｅｍａｔｉｃｄｉａｇｒａｍｏｆＳｈｅｎｚｈｅｎＲｉｖｅｒａｎｄＢａｙ个太阴日有两次涨潮和两次落潮，历时约２５ｈ，深圳湾口的烂角咀平均高潮位０．８６ｍ，罗湖水文站平均高潮位１．３３ｍ。深圳河受潮汐影响泥沙颗粒一般较小，床沙质泥沙中值粒径约为０．００４ｍｍ，悬移质泥沙中值粒径一般在０．００３～０．２００ｍｍ之间，属黏粒和粉砂。由于泥沙粒径较小，受咸水的影响易发生絮凝，从而极易沉降，造成泥沙淤积。深圳河干流含沙量一般在０．０２～０．３０ｋｇ／ｍ３之间，随潮汐周期性变化。深圳河口一般大潮涨潮输沙率大于落潮，导致泥沙不易往外扩散，也易造成淤积。区域概况1深圳河发源于梧桐山牛尾岭，全长３３．１ｋｍ，平均坡降１．４‰，流域面积３１２．５ｋｍ２，其中深圳境内１８７．５ｋｍ２，香港境内１２５．０ｋｍ２，图１为深圳河湾水系示意图。深圳河流域内气候湿润多雨，４～９月汛期降水量占全年８５％～９０％，降水强度大，暴雨多。近１０年来平均降水量２１８１．６ｍｍ，年最大降水量２９９５．６ｍｍ。流域内河流为雨源型河流，径流量主要受降水量控制，丰水年平均流量约１．６０×１０８ｍ３。深圳河干流全长１３．１６ｋｍ，受潮汐影响，为感潮河流。深圳湾潮汐为不规则半日混合潮，每深圳河淤积现状2收稿日期：２０１１－０４－１０，修回日期：２０１１－０５－２４基金项目：国家水体污染控制与治理科技重大专项基金资助项目（２００８ＺＸ０７４２３－００１）作者简介：胡明军（１９８６－），男，硕士研究生，研究方向为水环境数值模拟，Ｅ－ｍａｉｌ：ｈｍｊ０８＠ｍａｉｌ．ｔｓｉｎｇｈｕａ．ｅｄｕ．ｃｎ通讯作者：毛献忠（１９６８－），男，副研究员，研究方向为计算流体力学，Ｅ－ｍａｉｌ：ｍａｏｘｚ＠ｓｚ．ｔｓｉｎｇｈｕａ．ｅｄｕ．ｃｎ第２９卷第１０期胡明军等：深圳河口防洪能力数值研究·２３·为提高深圳河防洪能力，对深圳河进行了治理。工程按５０年一遇防洪标准，从深圳河口至平原河口约１３．５０ｋｍ河段进行裁弯取直、扩宽加深，三期工程相继完成［９］。前两期工程完成后，河床出现了严重泥沙淤积现象，深圳河的防洪能力大幅下降。基于深圳河历年汛前、汛后河床地形实测资料分析，２００５～２００６年淤泥总量基本维持于３００×１０４ｍ３，２００６～２００７年对河道进行清淤，约厚１．０ｍ。２００８年后河道回淤，至２００９年河道又基本达到动态平衡，图２为深圳河历年淤泥总量的变化。图３为２００９年汛后深圳河罗湖至深圳河口段平均高程沿程分布，对比设计高程，淤泥厚度基本在２．３～３．２ｍ之间，淤积严重。根据２００９年汛后实测地形资料计算，罗湖至深圳河口主河道上约淤积３００×１０４ｍ３泥沙。梧桐河口和布吉河口附近的河底平均高程分别达－１．３、－１．０ｍ，比邻近断面河底平均高程约高０．４ｍ，对该河段防洪造成了很大压力。平均１１０．２ｍ×１３９．４ｍ。③采用区域分解方法将深圳河与深圳湾计算域耦合计算。④整个区域面积８８．４４８ｋｍ２，共１０２６３个计算网格，见图４。图４计算网格及边界、验证点Ｆｉｇ．４Ｃｏｍｐｕｔａｔｉｏｎａｌｇｒｉｄ，ｂｏｕｎｄａｒｙａｎｄｖｅｒｉｆｉｃａｔｉｏｎｐｏｉｎｔｓ采用２００８年６月１１～１６日零点实测洪水过程所建的模型以验证罗湖及深圳河口两站实测数据，验证指标为潮位与流量；模拟区域地形数据采用２００８年汛前深圳河及２００７年汛后深圳湾实测地形数据。各支流采用实测流量数据计算其对干流的贡献。图５为罗湖与深圳河口潮位与流量验证结果。计算结果表明：①模型