年 月 文章编号： 水 利 学 报 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.58m/s～0.74m/s，小潮期为0.20m/s～0.34m/s。 收稿日期： 基金项目：项目名称（项目号） ；项目名称（项目号） 作者简介：第一作者姓名（1973-），男，籍贯，职称，博士生，主要从事水文水资源工作。E_mail:1@2.net 

深圳湾具有重要的航运功能，湾内有蛇口港、东角头港等主要港口。在湾顶附近的福田国家红树 林自然保护区和米埔自然保护区分布着大面积的泥滩和红树林，是华南地区最重要的湿地生态系统。 该区域为许多珍稀和濒危生物物种提供了栖息地，尤其是以水鸟和候鸟为主的鸟类，其中包括一些在 全球范围内受到威胁的种类［8］。 计划在深圳湾北岸进行的大规模填海造地工 程，不可避免地会改变泥沙的运动规律，并对港口、 航道和自然保护区造成影响。为减少填海对港口、 航道的负面影响和避免生态环境恶化，有必要对深 圳湾的填海工程的设计方案进行全面规划和综合论 证。拟议中的深圳湾北岸填海工程岸线方案共有五 种(见图1)，方案一至方案五的填海面积分别为8.8、 13.2、10. 95、5.88、8.40km2。文献［1］利用水动力 学数学模型，预测了深圳湾填海后纳潮量、流场、 潮位等海湾水动力条件的变化，并讨论了流速、潮 位变化对防洪、港口航道和红树林自然保护区的影 响。本文将在水动力学数学模型研究的基础上，进 一步探讨填海前后泥沙运动特征的变化及其影响。 3 泥沙模型 图 1 深圳湾填海方案示意 3.1 基本方程 根据实测资料［7］分析，冬季深圳湾各测站实测平均含沙量在0.0073～0.205kg/m3之间， 而夏季在0.0137～0.2385kg/m3之间。由于深圳湾水体中含沙量不高，泥沙对水流的影响极小，泥沙数 学模型与水动力学模型可以采用非耦合的方式求解。深圳湾的床沙组成为砂层及亚黏土层，上覆海相 淤泥层，底质细、中值粒径为0.004～0.006mm［9,10］，泥沙运动以悬移质运动为主，推移质运动可以忽 略。由于缺乏足够的悬移质 泥沙级配资料，本研究将悬沙按均匀沙处理。悬移质非平衡输沙方程和河 床变形方程分别如下［11］ S  t   uS  x   vS  y   ( H S  )  x     SK  x  x     y     SK  y  y      ( t  m  S  ) (1) (2) 式中：S为垂线平均含沙量；u，v分别为沿x和y方向的垂线平均流速；α为恢复饱和系数；ω为悬沙沉 降速度;H为实际水深；φ为水流挟沙力；KxKy分别为沿x和y方向的泥沙扩散系数；η为河床变形厚度(η <0为冲刷)；γm为泥沙干容重。水流挟沙力公式采用 常见的半经验公式［12］，即  m 3  VA  Hg      (3) 式中：A和m为经验系数；V为流速大小；g为重力加速度。 深圳湾的地形如图1。泥沙模型的计算范围从深圳河口到深圳湾口赤湾烂角嘴一线，共有168×132 个网格，网格大小为100m×100m，模型中的流速和实际水深，直接引用水动力学模型计算结果。差分 格式采用ADI（Alternative Direction Implicit）法［13］构造，即计算时段分为两个半时间步长，在前半个 时间步长内，只对空间x方向作变量运算，而后半个步长内，只对空间y方向作变量运算。差分格式中 对流项采用迎风格式。时间步长经调试取100s。计算域开边界条件根据泥沙浓度的监测资料获得；对 2 

于闭边界，假设  ，0  n n  S  / ，为闭边界外法线方向的单位矢量；泥沙方程的初始值由实测含沙量 给定。 3.2 模 型 验证 模型参数的率定采用1998年1月的实测水文资料［7］进行。在深圳湾的水文流监测点位 中(如图1)，VS3和VS4接近湾口开边界，利用这两点的泥沙浓度时间序列，通过线性插值得到湾口开边 界。 （a）VS1 （b）VS2 图 2 泥沙含量计算值与测量值比较 假设枯水期(或丰水期)每个月的来沙量都等于模型预测的枯水期(或丰水期)月来沙量，以枯水期5 个月、丰水期7个月将12个月的来沙量累加，可以近似得到深圳湾的年来沙量。如表2所示，除方案二 对应的年来沙量减少外，其它方案的年来沙量都略有增加，其中方案四增加得最多，而方案三的年来 沙量变化最小。因此，单就进出深圳湾泥沙变化最小的要求来看，各方案的优劣顺序为：方案三、一、 五、二、四。 表2 填海前后深圳湾月来沙量和年来沙量变化(单位：103t) 方案 现 状 方案一 方案二 方案三 枯水期(30d) 丰水期(30d) 年(365d) 填海后年内沙量的增-减比例(%) 14.5 18.5 19.5 19.0 58.4 56.4 53.8 55.5 481.3 487.3 474.1 483.5 1.25 -1.50 0.46 对各种填海岸线方案的影响进行的分析表明：不同的考虑对象要求的泥沙评价指标及大小不同， 因而对应岸线方案的优劣次序排列也有所差异(如表4)，所以不能凭单个指标来确定填海方案。综合考 虑各种泥沙评价指标，可以认为方案三、四较优。 方案 方案一 方案二 方案三 年来 沙量 2 4 1 泥沙冲 淤分布 5 2 4 表4 方案优劣比较 港口航道、自然保护区冲淤变化 蛇口港 东角头港 内湾航道 外湾航道 福田 米埔 3 1 2 3 5 4 4 1 2 3 5 4 4 3 2 2 5 3 水动力学模型和泥沙数学模型为评估填海前后水动力学条件和泥沙运动特征变化提供了基础。填 海工程对海湾水质、环境容量和污染物排放方式方面的影响需要利用其它模型。 3 

4 结论 参 考 文 献： [ 1 ] 林秉南，赵雪华，施麟宝. 河口建坝对毗邻海湾潮波影响的计算(二维特征线理论法)[J]. 水利学报，1980， (3):16-25. [ 2 ] 潘少明，施晓冬，王建业，Peter Tanner，Lai Shing Leong. 围海造地工程对香港维多利亚现代沉积作用 的影响[J]. 沉积学报，2000，18，(1):22-28. [ 3 ] 洪华生，陈宗团. 海岸带综合管理中面临的科学问题. 海洋管理，1998， 1：28-31. [ 4 ] 松林，丁平兴. 湛江湾沿岸工程冲淤影响的预测分析，Ⅱ. 冲淤的数值计算[J]. 海洋学报，1997，19， (1):64-72. [ 5 ] Han, Z. C., and Cheng H. P.. Two-dimensional sediment mathematical model of Hangzhou Bay. Proc. 3rd Intern. Symp. on River Sedimentation[J]. The University of Mississippi, USA, 1986, 463-471. [ 6 ] Chien, N, and Wan, Z.H. Mechanics of Sediment Transport[D]. American Society of Civil Engineers, Virginia: [ 7 ] ASCE Press, 1999. Leendertse, J. J.. Aspects of SYMSYS2D: a system for two Corporation, dimensional flow computation. The Rand [ 8 ] 隋淑珍，张乔民. 华南沿海红树林海岸沉积物特征分析[J]. 热带海洋，1999，18(4):18-23. [ 9 ] 倪晋仁，杨小毛，王光谦. 深港交界带经济开发过程中泥沙对生态环境的影响[J]. 地理学报，1998， 53(4):350-355. [ 10 ] 张乔民，于红兵，陈欣树，等. 红树林生长带与潮汐水位关系的研究[J]. 生态学报，1997, 17(3)：258-265. Determination of Post-reclamation Coastline Based on Sediment Transport Model 姓名1,2，姓名1,2，姓名2 （注：全英文） their impacts on (1.单位，市 邮编，China； 2. 单位，市 邮编，China)（注：全英文） Abstract: Attention should be paid on variations of sediment transport characteristics and ports and navigation channels before reclaiming a bay. In general, a bay system would be sensitive to the change of sediment erosion or deposition, and thus the ecological assessment in terms of sediment models could be important for the feasibility study of the proposed options. By taking the reclamation of Shenzhen Bay as an example, this paper used a 2-D sediment transport model to predict the variations of net sediment inflow, sediment concentration, sediment erosion or deposition for five designated reclamation options. Studies were made for potential impacts on ports, navigation channels and ecosystems of mangrove according to modeling results. Then, the priority order of the five options for all the concerned issues was discussed. It was concluded that although the priority order for each of the issues is different, a harmonized optimal post reclamation coastline could be obtained from viewpoint of integrated impact of the reclamation options to all the concerned issues. Keywords: sediment transport model; reclamation; coastline 4