水资源学报.doc

水利学报 年 月 SHUILI XUEBAO 第 卷 第 期 文章编号： 题目 姓名1,2，姓名1,2，姓名2 (1.单位，省 市 邮编； 2. 单位，省 市 邮编) 摘要：泥沙运动特征的变化及其对港口航道的影响是大型海岸工程实施时需要考虑的主要问题。由于海岸生态系统往往对泥沙冲淤的变化非常敏感，生态评估也逐渐成为工程规划及方 案可行性论证研究的重要内容。本文利用泥沙数学模型，以深圳湾为例研究了五种填海方案造成的海湾来沙量、含沙量和泥沙冲淤分布的变化，讨论了泥沙淤积速率的变化对港口航道和红树林生态系统的影响，并通过岸线方案的对比指出了在分别考虑各种泥沙评价指标时对 应的五种填海方案的优劣。定量分析表明，不同评价指标对应的方案优劣次序有所差异，协调多种指标要求的综合优化岸线可以通过比较各种方案后给出。 关键词：泥沙数学模型；填海工程；岸线；深圳湾 中图分类号：TV148 文献标识码：A 1 研究背景 填海对当地或周边地区潮波特性的影响，已经在许多填海工程的岸线确定中得到充分重视［1］。近年来，填海工程造成的区域泥沙运动特征的改变以及对港口、航道和海岸生态系统的影响也愈来愈受到广泛的关注。泥沙淤积带来的问题尤其复杂。例如，香港维多利亚港附近由于填海造地使得一些区域潮流作用减弱，形成了有利于泥沙沉降的环境，沉积速率明显增大［2］；厦门西海域的海岸工程和围垦造地使纳潮量减少了67%左右，潮汐对泥沙的输运能力因此大大减弱，使厦门港的淤积加重、通航条件受到影响，每年必须花800多万元费用清理航道；此外，沉积环境的改变还导致珍贵的文昌鱼几乎绝迹，海豚减少［3］。 过去关于填海后泥沙运动特征变化及其对港口航道影响的研究较多［4～6］，对海岸生态系统影响的研究却较少。本文以深圳湾填海工程为例，尝试利用泥沙数学模型对填海前后海湾的来沙量、冲淤分布规律、港口航道和生态保护区的影响进行综合考虑，通过对比分析典型填海岸线方案的优劣为填海工程的优化提供依据。 2 研究区域概况 深圳湾是深圳和香港的交界水域，为珠江口东岸的一个半封闭型浅水湾(如图1)，面积约为90km2，平均水深2.9m。注入深圳湾的河流中，深圳河是最大的河流。受南海不规则半日混合潮影响，海湾潮流方向为往复流，年均潮差1.37m。根据近来的实测资料［7］和水动力学模型计算结果①，涨潮时深圳湾内大潮期流速为0.39m/s～0.65m/s，小潮期为0.31m/s～0.40m/s；落潮时，深圳湾内大潮期流速为0 收稿日期： 基金项目：项目名称（项目号） ；项目名称（项目号） 作者简介：第一作者姓名（1973-），男，籍贯，职称，博士生，主要从事水文水资源工作。E_mail:1@2.net 深圳湾具有重要的航运功能，湾内有蛇口港、东角头港等主要港口。在湾顶附近的福田国家红树林自然保护区和米埔自然保护区分布着大面积的泥滩和红树林，是华南地区最重要的湿地生态系统。该区域为许多珍稀和濒危生物物种提供了栖息地，尤其是以水鸟和候鸟为主的鸟类，其中包括一些在全球范围内受到威胁的种类［8］。 图1深圳湾填海方案示意计划在深圳湾北岸进行的大规模填海造地工程，不可避免地会改变泥沙的运动规律，并对港口、航道和自然保护区造成影响。为减少填海对港口、航道的负面影响和避免生态环境恶化，有必要对深圳湾的填海工程的设计方案进行全面规划和综合论证。拟议中的深圳湾北岸填海工程岸线方案共有五种(见图1)，方案一至方案五的填海面积分别为8.8、13.2、10. 95、5.88、8.40km2。文献［ 图1深圳湾填海方案示意 3 泥沙模型 3.1 基本方程 根据实测资料［7］分析，冬季深圳湾各测站实测平均含沙量在0.0073～0.205kg/m3之间，而夏季在0.0137～0.2385kg/m3之间。由于深圳湾水体中含沙量不高，泥沙对水流的影响极小，泥沙数学模型与水动力学模型可以采用非耦合的方式求解。深圳湾的床沙组成为砂层及亚黏土层，上覆海相淤泥层，底质细、中值粒径为0.004～0.006mm［9,10］，泥沙运动以悬移质运动为主，推移质运动可以忽略。由于缺乏足够的悬移质 泥沙级配资料，本研究将悬沙按均匀沙处理。悬移质非平衡输沙方程和河床变形方程分别如下［11 (1) (2) 式中：S为垂线平均含沙量；u，v分别为沿x和y方向的垂线平均流速；α为恢复饱和系数；ω为悬沙沉降速度;H为实际水深；φ为水流挟沙力；KxKy分别为沿x和y方向的