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.　Key　words:　sediment　transport　model;　reclamation;　coastline(英文关键词:5号,Times　New　Roman)　1　研究背景(一级标题:小四号宋体,占3行,左齐)　(正文:5号宋体,通栏,数字和单位不空,排方程断行时+-+÷在上一行行尾)填海对当地或周边地区潮波特性的影响,已经在许多填海工程的岸线确定中得到充分重视[1]。近年来,填海工程造成的区域泥沙运动特征的改变以及对港口、航道和海岸生态系统的影响也愈来愈受到广泛的关注。泥沙淤积带来的问题尤其复杂。例如,香港维多利亚港附近由于填海造地使得一些区域潮流作用减弱,形成了有利于泥沙沉降的环境,沉积速率明显增大[2];厦门西海域的海岸工程和围垦造地使纳潮量减少了67%左右,潮汐对泥沙的输运能力因此大大减弱,使厦门港的淤积加重、通航条件受到影响,每年必须花800多万元费用清理航道;此外,沉积环境的改变还导致珍贵的文昌鱼几乎绝迹,海豚减少[3]。　过去关于填海后泥沙运动特征变化及其对港口航道影响的研究较多[4-6],对海岸生态系统影响的研究却较少。本文以深圳湾填海工程为例,尝试利用泥沙数学模型对填海前后海湾的来沙量、冲淤分布规律、港口航道和生态保护区的影响进行综合考虑,通过对比分析典型填海岸线方案的优劣为填海工程的优化提供依据。　2　研究区域概况　深圳湾是深圳和香港的交界水域,为珠江口东岸的一个半封闭型浅水湾(如图1),面积约为90km2,平均水深2.9m。注入深圳湾的河流中,深圳河是最大的河流。受南海不规则半日混合潮影响,海湾潮流方向为往复流,年均潮差1.37m。根据近来的实测资料[7]和水动力学模型计算结果,涨潮时深圳湾内大潮期流速为0.39~0.65m/s,小潮期为0.31~0.40m/s;落潮时,深圳湾内大潮期流速为0.58~0.74m/s,小潮期为0.20~0.34m/s。　深圳湾具有重要的航运功能,湾内有蛇口港、东角头港等主要港口。在湾顶附近的福田国家红树林自然保护区和米埔自然保护区分布着大面积的泥滩和红树林,是华南地区最重要的湿地生态系统。该区域为许多珍稀和濒危生物物种提供了栖息地,尤其是以水鸟和候鸟为主的鸟类,其中包括一些在全球范围内受到威胁的种类[8]。　计划在深圳湾北岸进行的大规模填海造地工程,不可避免地会改变泥沙的运动规律,并对港口、航道和自然保护区造成影响。为减少填海对港口、航道的负面影响和避免生态环境恶化,有必要对深圳湾的填海工程的设计方案进行全面规划和综合论证。拟议中的深圳湾北岸填海工程岸线方案共有五种(见图1),方案一至方案五的填海面积分别为8.8、13.2、10.　95、5.88、8.40km2。文献[1]利用水动力学数学模型,预测了深圳湾填海后纳潮量、流场、潮位等海湾水动力条件的变化,并讨论了流速、潮位变化对防洪、港口航道和红树林自然保护区的影响。本文将在水动力学数学模型研究的基础上,进一步探讨填海前后泥沙运动特征的变化及其影响。　3　泥沙模型　3.1　基本方程(二级标题:5号黑体,占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方向的泥沙扩散系数;η为河床变形厚度(η<0为冲刷);γm为泥沙干容重。水流挟沙力公式采用常见的半经验公式[12],即(3)　式中:A和m为经验系数;V为流速大小;g为重力加速度。　深圳湾的地形如图1。泥沙模型的计算范围从深圳河口到深圳湾口赤湾烂角嘴一线,共有168×132个网格,网格大小为100m×100m,模型中的流速和实际水深,直接引用水动力学模型计算结果。差分格式采用ADI(Alternative　Direction　Implicit)法[13]构造,即计算时段分为两个半时间步长,在前半个时间步长内,只对空间x方向作变量运算,而后半个步长内,只对空间y方向作变量运算。差分格式中对流项采用迎风格式。时间步长经调试取100s。计算域开边界条件根据泥沙浓度的监测资料获得;对于闭边界,假设,为闭边界外法线方向的单位矢量;泥沙方程的初始值由实测含沙量给定。　3.1.1　三级标题　(5号楷体,占1行,左齐,与正文接排(若下面接标题则另起一行))　模型验证　模型参数的率定采用1998年1月的实测水文资料[7]进行。在深圳湾的水文流监测点位中(如图1),VS3和VS4接近湾口开边界,利用这两点的泥沙浓度时间序列,通过线性插值得到湾口开边界。　假设枯水期(或丰水期)每个月的来沙量都等于模型预测的枯水期(或丰水期)月来沙量,以枯水期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　(表格:6号宋体,居中排,小数点对位)　对各种填海岸线方案的影响进行的分析表明:不同的考虑对象要求的泥沙评价指标及大小不同,因而对应岸线方案的优劣次序排列也有所差异(如表4),所以不能凭单个指标来确定填海方案。综合考虑各种泥沙评价指标,可以认为方案三、四较优。　表4　方案优劣比较　方案年来　沙量　泥沙冲　淤分布　港口航道、自然保护区冲淤变化　蛇口港　东角头港　内湾航道　外湾航道　福田　米埔　方案一　方案二　方案三　2　4　1　5　2　4　3　1　2　3　5　4　4　1　2　3　5　4　4　3　2　2　5　3　水动力学模型和泥沙数学模型为评估填海前后水动力学条件和泥沙运动特征变化提供了基础。填海工程对海湾水质、环境容量和污染物排放方式方面的影响需要利用其它模型。　4　结论　参考文献:　林秉南,赵雪华,施麟宝.　河口建坝对毗邻海湾潮波影响的计算(二维特征线理论法)[J].　水利学报,1980(3):16-25.　潘少明,施晓冬,王建业,等.　围海造地工程对香港维多利亚现代沉积作用的影响[J].　沉积学报,2000,18(1):22-28.　　洪华生,陈宗团.　海岸带综合管理中面临的科学问题.　海洋管理,1998(1):28-31.　　松林,丁平兴.　湛江湾沿岸工程冲淤影响的预测分析Ⅱ.　冲淤的数值计算[J].　海洋学报,1997,19(1):64-72.　　Han　Z　C,　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.　　Chien　N,　Wan　Z　H.　Mechanics　of　Sediment　Transport[D].　American　Society　of　Civil　Engineers,　Virginia:　ASCE　Press,　1999.　　Leendertse　J　J.　Aspects　of　SYMSYS2D:　a　system　for　two　dimensional　flow　computation[M].　The　Rand　Corporation,　1990.　隋淑珍,张乔民.　华南沿海红树林海岸沉积物特征分析[J].　热带海洋,1999,18(4):18-23.　倪晋仁,杨小毛,王光谦.　深港交界带经济开发过程中泥沙对生态环境的影响[J].　地理学报,1998,53(4):350-355.　张乔民,于红兵,陈欣树,等.　红树林生长带与潮汐水位关系的研究[J].　生态学报,1997,　17(3):258-265.　(参考文献:小五号,宋体)