《水利学报》论文样板doc.doc

水利学报年 月 SHUILI XUEBAO 第 卷 第 期文章编号：文章编号： 题目姓名1,2，姓名1,2，姓名2(1.单位，省 市 邮编； 2. 单位，省 市 邮编)摘要：摘要：泥沙运动特征的变化及其对港口航道的影响是大型海岸工程实施时需要考虑的主要问题。由于海岸生态系统往往对泥沙冲淤的变化非常敏感，生态评估也逐渐成为工程规划及方 案可行性论证研究的重要内容。本文利用泥沙数学模型，以深圳湾为例研究了五种填海方案造成的海湾来沙量、含沙量和泥沙冲淤分布的变化，讨论了泥沙淤积速率的变化对港口航道和红树林生态系统的影响，并通过岸线方案的对比指出了在分别考虑各种泥沙评价指标时对 应的五种填海方案的优劣。定量分析表明，不同评价指标对应的方案优劣次序有所差异，协调多种指标要求的综合优化岸线可以通过比较各种方案后给出。 关键词：关键词：泥沙数学模型；填海工程；岸线；深圳湾中图分类号：中图分类号： 文献标识码：文献标识码：A 1 研究背景填海对当地或周边地区潮波特性的影响，已经在许多填海工程的岸线确定中得到充分重视1。 近年来，填海工程造成的区域泥沙运动特征的改变以及对港口、航道和海岸生态系统的影响也愈来愈 受到广泛的关注。泥沙淤积带来的问题尤其复杂。例如，香港维多利亚港附近由于填海造地使得一些 区域潮流作用减弱，形成了有利于泥沙沉降的环境，沉积速率明显增大2；厦门西海域的海岸工程 和围垦造地使纳潮量减少了67%左右，潮汐对泥沙的输运能力因此大大减弱，使厦门港的淤积加重、 通航条件受到影响，每年必须花800多万元费用清理航道；此外，沉积环境的改变还导致珍贵的文昌 鱼几乎绝迹，海豚减少3。 过去关于填海后泥沙运动特征变化及其对港口航道影响的研究较多4-6，对海岸生态系统影响 的研究却较少。本文以深圳湾填海工程为例，尝试利用泥沙数学模型对填海前后海湾的来沙量、冲淤 分布规律、港口航道和生态保护区的影响进行综合考虑，通过对比分析典型填海岸线方案的优劣为填 海工程的优化提供依据。2 研究区域概况深圳湾是深圳和香港的交界水域，为珠江口东岸的一个半封闭型浅水湾(如图1)，面积约为 90km2，平均水深2.9m。注入深圳湾的河流中，深圳河是最大的河流。受南海不规则半日混合潮影响， 海湾潮流方向为往复流，年均潮差1.37m。根据近来的实测资料7和水动力学模型计算结果，涨潮 时深圳湾内大潮期流速为0.390.65m/s，小潮期为0.310.40m/s；落潮时，深圳湾内大潮期流速为 0.580.74m/s，小潮期为0.200.34m/s。收稿日期： 基金项目：项目名称（项目号） ；项目名称（项目号）作者简介：第一作者姓名（1973-），男，籍贯，职称，博士生，主要从事水文水资源工作。E_mail:12.net2深圳湾具有重要的航运功能，湾内有蛇口港、东角头港等主要港口。在湾顶附近的福田国家红树 林自然保护区和米埔自然保护区分布着大面积的泥滩和红树林，是华南地区最重要的湿地生态系统。 该区域为许多珍稀和濒危生物物种提供了栖息地，尤其是以水鸟和候鸟为主的鸟类，其中包括一些在 全球范围内受到威胁的种类8。 计划在深圳湾北岸进行的大规模填海造地工程， 不可避免地会改变泥沙的运动规律，并对港口、航 道和自然保护区造成影响。为减少填海对港口、航 道的负面影响和避免生态环境恶化，有必要对深圳 湾的填海工程的设计方案进行全面规划和综合论证。 拟议中的深圳湾北岸填海工程岸线方案共有五种 (见图1)，方案一至方案五的填海面积分别为 8.8、13.2、10. 95、5.88、8.40km2。文献1 利用水动力学数学模型，预测了深圳湾填海后纳潮 量、流场、潮位等海湾水动力条件的变化，并讨论 了流速、潮位变化对防洪、港口航道和红树林自然 保护区的影响。本文将在水动力学数学模型研究的 基础上，进一步探讨填海前后泥沙运动特征的变化 及其影响。3 泥沙模型 3.1 基本方程 根据实测资料7分析，冬季深圳湾各测站实测平均含沙量在0.00730.205kg/m3之 间，而夏季在0.01370.2385kg/m3之间。由于深圳湾水体中含沙量不高，泥沙对水流的影响极小， 泥沙数学模型与水动力学模型可以采用非耦合的方式求解。深圳湾的床沙组成为砂层及亚黏土层，上 覆海相淤泥层，底质细、中值粒径为0.0040.006mm9-10，泥沙运动以悬移质运动为主，推移质运 动可以忽略。由于缺乏足够的悬移质 泥沙级配资料，本研究将悬沙按均匀沙处理。悬移质非平衡输 沙方程和河床变形方程分别如下11(1) ySKyxSKxSHyvS xuS tSyx)(2)(Stm式中：S为垂线平均含沙量；u、v分别为沿x和y方向的垂线平均流速；为恢复饱和系数；为悬沙 沉降速度;H为实际水深；为水流挟沙力；KxKy分别为沿x和y方向的泥沙扩散系数；为河床变形厚 度(<0为冲刷)；m为泥沙干容重。水流挟沙力公式采用常见的半经验公式12，即(3)mHgVA 3式中：A和m为经验系数；V为流速大小；g为重力加速度。 深圳湾的地形如图1。泥沙模型的计算范围从深圳河口到深圳湾口赤湾烂角嘴一线，共有 168×132个网格，网格大小为100m×100m，模型中的流速和实际水深，直接引用水动力学模型计算结 果。差分格式采用ADI（Alternative Direction Implicit）法13构造，即计算时段分为两个半时间步长， 在前半个时间步长内，只对空间x方向作变量运算，而后半个步长内，只对空间y方向作变量运算。差图 1 深圳湾填海方案示意3分格式中对流项采用迎风格式。时间步长经调试取100s。计算域开边界条件根据泥沙浓度的监测资料获得；对于闭边界，假设，为闭边界外法线方向的单位矢量；泥沙方程的初始值由实nnS，0/测含沙量给定。 3.2 模型验证 模型参数的率定采用1998年1月的实测水文资料7进行。在深圳湾的水文流监测点 位中(如图1)，VS3和VS4接近湾口开边界，利用这两点的泥沙浓度时间序列，通过线性插值得到湾口 开边界。假设枯水期(或丰水期)每个月的来沙量都等于模型预测的枯水期(或丰水期)月来沙量，以枯水期 5个月、丰水期7个月将12个月的来沙量累加，可以近似得到深圳湾的年来沙量。如表2所示，除方案 二对应的年来沙量减少外，其它方案的年来沙量都略有增加，其中方案四增加得最多，而方案三的年 来沙量变化最小。因此，单就进出深圳湾泥沙变化最小的要求来看，各方案的优劣顺序为：方案三、 一、五、二、四。表2 填海前后深圳湾月来沙量和年来沙量变化 (单位：103t)方 案枯水期(30d)丰水期(30d)年(365d)填海后年内沙量的增-减比例/%现 状方案一方案二方案三14.518.519.519.058.456.453.855.5481.3487.3474.1483.51.25-1.500.46对各种填海岸线方案的影响进行的分析表明：不同的考虑对象要求的泥沙评价指标及大小不同， 因而对应岸线方案的优劣次序排列也有所差异(如表4)，所以不能凭单个指标来确定填海方案。综合 考虑各种泥沙评价指标，可以认为方案三、四较优。表4 方案优劣比较港口航道、自然保护区冲淤变化 方 案年来沙量泥沙冲淤分布蛇口港东角头港内湾航道外湾航道福田米埔方案一方案二方案三241524312354412354432253水动力学模型和泥沙数学模型为评估填海前后水动力学条件和泥沙运动特征变化提供了基础。填 海工程对海湾水质、环境容量和污染物排放方式方面的影响需要利用其它模型。（a）VS1 （b）VS2图 2 泥沙含量计算值与测量值比较44 结论参参 考考 文文献：献： 1 林秉南，赵雪华，施麟宝. 河口建坝对毗邻海湾潮波影响的计算(二维特征线理论法)J. 水利学报，1980(3):16-25. 2 潘少明，施晓冬，王建业，等. 围海造地工程对香港维多利亚现代沉积作用的影响J. 沉积学报，2000，18(1):22-28. 3 洪华生，陈宗团. 海岸带综合管理中面临的科学问题. 海洋管理，1998（1）：28-31. 4 松林，丁平兴. 湛江湾沿岸工程冲淤影响的预测分析. 冲淤的数值计算J. 海洋学报，1997，19(1):64-72. 5 Han Z C, Cheng H P. Two-dimensional sediment mathematical model of Hangzhou Bay. Proc. 3rd Intern. Symp. on River SedimentationJ. The University of Mississippi, USA, 1986, 463-471. 6 Chien N, Wan Z H. Mechanics of Sediment TransportD. American Society of Civil Engineers, Virginia: ASCE Press, 1999. 7 Leendertse J J. Aspects of SYMSYS2D: a system for two dimensional flow computationM. The Rand Corporation, 1990. 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（注：全英文）(1.单位，市 邮编，China； 2. 单位，市 邮编，China)（注：全英文）Abstract: Attention should be paid on variations of sediment transport characteristics and their impacts on 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.4Key words: sediment transport model; reclamation; coastline