离岸人工岛渗流和防渗数值模拟研究.docVIP

离岸人工岛渗流和防渗数值模拟研究 　　摘 要：采用FEFLOW模型，对海南海花岛内海堤坝及海花岛围堤无防渗和有防渗措施下的渗流量、渗透范围及对内海水位变化的影响进行了计算和分析。结果显示：（1）实施防渗措施后，在正常潮位情况下内海水位波动减小0.03m～0.13m；在极端情况下（50年极高水位、50年极低水位）分别使内海潮位波动减小0.29m～0.39m；其渗透坡降均小于允许渗透坡降值。（2）内海堤坝及其两侧为渗流发生的主要区域，建议对渗流重点区域进行适当范围内进行防渗处理；根据计算结果，给出了防渗的大致范围。（3）在渗流影响下，人工岛成陆区，浸润线（饱和度1）的位置，基本在高程1.1m～1.5m（国家85基面）之间。为了人工岛成陆区土质和植被的生长环境，建议对成陆区适当深度进行反渗处理。 　　关键词：海花岛；人工岛；渗流；防渗措施；数值模拟 　　1 概述 　　近年来，我国围海造地的建设项目迅速发展。在围填海整体布局上，多采取人工岛式、连片分离式等平面形态，形成一定面积的内河（内海）水域。在离岸人工岛建设工程的关键技术问题上，前人将主要研究方向集中在平面形态、护岸结构、地基处理和生态环保等问题上，而对于渗流问题、渗流引起的内河（内海）水位变化的研究却极少[1][2]。在现有的文献中，戚蓝等[3]采用Ansys有限元数值模拟方法，对福建某人工岛渗流与外海、内河水位变动之间的相互影响进行了计算与分析，并提出了防渗措施和建议；庄志农[2]针对目前围海造陆工程所形成的内河（水域）的防渗要求，以福建某典型工程为例，采用Ansys有限元数值模拟方法对无防渗措施情况下的渗透坡降、流速、流量进行了计算与分析，表明各闸门两侧为渗流发生的主要区域，建议对渗流重点区域，采用混凝土垂直防渗墙的防渗措施，并给出了最佳方案。 　　文章采用基于有限元方法的FEFLOW水流和物质迁移模型，以海花岛工程为例，对其内海堤坝和海花岛围堤无防渗、有防渗措施下的渗流量、渗透范围及对内海水位变化的影响进行计算和分析，对工程设计提供一定参考。 　　海南海花岛项目位于儋州市排浦港与洋浦湾之间的海湾区域，南起排浦镇，北至马井镇，总跨度6.8km，拟在工程区域新建护岸，吹填造陆，共建3个人工岛[4]，人工岛之间以联系桥连接，总圈围形成陆地面积7.53km2，形成海岸线39.44km，见图1。三个人工岛与内陆之间建设四条堤坝形式的桥梁，与外海隔断（利用箱涵开关在涨落潮期间进行内外水体自然交换），使3个岛屿与东南侧内陆之间形成内海，面积约3.7km2，回填陆地标高为+7.0m～+7.5m（85国家高程，下同），内海水域设计底标高为-3.0m，水面常水位高程为+1.5m。 　　2 工程区地质特征与渗流风险 　　整个人工岛调查区域钻探揭露的地表下50m范围揭露的地层综合划分为6个岩性单元层，自上而下各土层为：①1含细粒土粉砂（Q4m）、①2珊瑚碎屑（Q4m）、②低液限粘土（Q4m）、③海滩岩（Q1m）、④低液限粘土（Q1m）、⑤低液限粘土（Q1m）、⑤1含细粒土粗砂（Q1m）和⑥低液限粘土（N2m）。 　　根据地质勘察结果显示，海花岛工程周围均存在珊瑚碎屑强透水层；而且对已形成陆域的1#岛的勘察，拟建场地浅层地下水含水层系赋存于①流泥、①-1中砂混淤泥、②珊瑚碎屑中的孔隙型潜水及各粘土层中的少量孔隙型潜水；以海水补给为主，水量丰富，场地地下水与海水水力作用联系较紧密。而且内海与外海之间存在一定的水位差，当多种条件汇聚在一起，满足形成渗流的水力条件，进而产生一定量的渗流量和渗流流速，对内海水位波动产生一定影响。 　　3 渗流模型与计算 　　3.1 模型介绍 　　FEFLOW（Finite Element Subsurface Flow System）是德国 WASY水资源规划和系统研究所于20世纪70年代末开发的数值模拟软件，是迄今为止功能最为齐全的地下水流和物质迁移过程的模拟软件，已在地下水污染防控和修复、地表水/地下水耦合模拟、大坝渗流计算、下渗计算、环境影响评价、地下水监测规划等多个领域得到广泛应用[5-8]。在人工岛渗流方面的研究文章尚属首例。模型原理见参考文献[9][10]，在此不再赘述。 　　3.2 模型的建立与处理 　　整体三维计算模型按照海花岛实际尺寸和周边实测地形建立，整个模型范围，南北8.9km，东西4.5km，内海底高程-3.0m。根据实际地质勘察情况对模型中地质进行了分层，共分了6层（流泥、含细粒土粉砂、珊瑚碎屑、粘土粉砂、粗砂和粘土粉砂），如图2。 　　具体堤坝三维模型，按照内海围堤4条堤坝实际尺寸建立三维模型，其中1#堤坝三维结构图如图3。 　　3.2.1 参数设置 　　各土层材料参数主要为渗透系数K，护岸抛填石料K=1×1