EN中文

地理学报 >

2020 , Vol. 75 >Issue 5: 1053 - 1064

DOI: https://doi.org/10.11821/dlxb202005012

气候变化与地表过程

南海珊瑚礁人工岛淡水透镜体形成过程及影响因素

  • 韩冬梅 1, 2 ,
  • 曹国亮 3 ,
  • 宋献方 1, 2
展开
  • 1.中国科学院地理科学与资源研究所 中国科学院陆地水循环及地表过程重点实验室,北京 100101
  • 2.中国科学院大学资源与环境学院,北京 100101
  • 3.中国水利水电科学研究院,北京 100083

韩冬梅(1978-), 女, 新疆焉耆人, 博士, 副研究员, 研究方向为滨海地下水水文过程和海水入侵。E-mail: handm@igsnrr.ac.cn

收稿日期: 2019-06-03

  要求修回日期: 2020-02-10

  网络出版日期: 2020-07-25

基金资助

中国科学院战略性先导科技专项(XDA13010303)

版权

版权所有，未经授权，不得转载、摘编本刊文章，不得使用本刊的版式设计。
收起

Formation processes and influencing factors of freshwater lens in artificial island of coral reef in South China Sea

  • HAN Dongmei 1, 2 ,
  • CAO Guoliang 3 ,
  • SONG Xianfang 1, 2
Expand
  • 1.Key Laboratory of Water Cycle & Related Land Surface Processes, Institute of Geographic Sciences and Natural Resources Research, CAS, Beijing 100101, China;
  • 2.College of Resources and Environment, University of Chinese Academy of Sciences, Beijing 100101, China
  • 3.China Institute of Water Resources and Hydropower Research, Beijing 100083, China

Received date: 2019-06-03

  Request revised date: 2020-02-10

  Online published: 2020-07-25

Supported by

Strategic Priority Research Program of the Chinese Academy of Sciences(XDA13010303)

Copyright

Copyright reserved © 2020
Fold

摘要

中国目前在珊瑚礁上建设人工岛的目标是建成生态宜居岛屿,人工岛地下海水在降雨入渗驱替作用下逐渐淡化,形成的淡水透镜体是支撑岛屿生态系统的重要水源。本文针对南海珊瑚礁普遍存在的上土下岩的二元地质结构,考虑不同测量方法得到的南海人工岛填筑所用珊瑚砂渗透系数存在较大差异,建立了珊瑚礁人工岛淡水驱替海水形成淡水透镜体的动态模拟模型,对南海某珊瑚礁人工岛开展了系列数值模拟研究。结果表明,人工岛礁二元地质结构不同的渗透性对人工岛淡水开始形成时间无明显影响,填岛完成2年左右的时间可观测到淡水,估计（或理论计算）在10~20年内可形成相对稳定的淡水透镜体;在人工岛面积达到一定规模的前提下,全新世珊瑚砂沉积物渗透性是控制人工岛礁淡水透镜体的主要因素,能够形成的淡水透镜体体积随渗透系数增大而呈对数减小,并且淡水透镜体体积随全新世沉积物渗透性增大而减小的速率较淡水透镜体厚度更加显著;淡水透镜体越过不整合面进入更新世礁灰岩后,由于含水层渗透系数突然增大,淡水透镜体向水平方向扩展,在礁灰岩顶界形成“髭”状淡水体。

关键词： 人工岛礁; 二元地质结构; 淡水透镜体; 珊瑚砂渗透性; 南海

本文引用格式

韩冬梅 , 曹国亮 , 宋献方 . 南海珊瑚礁人工岛淡水透镜体形成过程及影响因素[J]. 地理学报, 2020 , 75(5) : 1053 -1064 . DOI: 10.11821/dlxb202005012

Abstract

The goal of constructing artificial islands on coral reefs in China is to build ecological habitable islands. The seawater in the artificial islands is gradually replaced by rainfall infiltration, and the freshwater lens formed is an important water source supporting the island ecosystem. The double-layer geological structure is widespread in coral reefs in the South China Sea. Considering the fact that there exists a great large difference in the permeability coefficients obtained using different measurement methods, this study established numerical simulation models for a coral reef island in the South China Sea to simulate the seawater refreshing process. The results show that different permeabilities of the double-layer geological structure of the artificial islands have no obvious influence on the formation time of fresh water on the artificial islands. The permeability of the Holocene coral sand sediment is the main factor controlling the freshwater lens of the artificial reef. The volume of the forming freshwater lens decreases logarithmically with the increase of the permeability coefficient, and the volume of the freshwater lens decreases more significantly with the increase of the permeability of the Holocene sediment than the thickness of the freshwater lens. After the freshwater lens crossed the unconformity surface and entered the Pleistocene reef limestone, due to the sudden increase of the permeability coefficient of the aquifer, the freshwater lens expanded to the horizontal direction and formed the mustached-shape freshwater body in the reef limestone.

Key words： artificial island; double-layer geological structure; freshwater lens; coral sand permeability; South China Sea

1 引言

人工岛是填海造陆的一种方式,是指人工建造而非自然形成的岛屿[1]。与大陆沿岸的填海造陆不同,南海人工岛主要是在珊瑚礁上建造的,珊瑚礁沉积为二元地质结构,表层为全新世无胶结的珊瑚碎屑沉积,下层为更新世已成岩的礁灰岩[2]。环礁是南海诸岛中最为常见的一种珊瑚礁地貌类型,其上建设人工岛采取“自然仿真”的技术路线,绞吸或挖取澙湖中的珊瑚砂,向礁坪上吹填形成潮上陆域即人工岛[1]。人工岛面积达到一定规模后,由于入渗雨水比重小,浮于岛内海水之上,形似透镜体,即“淡水透镜体”[3,4]。淡水透镜体的形成,对珊瑚礁人工岛土壤淡化—岛屿绿化过程起到关键作用,其形成时间和规模对仿真造岛进程中岛屿生态系统的形成有重要影响。
珊瑚岛淡水透镜体的形成和影响因素已有大量研究工作,淡水透镜体的体积规模与岛屿的地质结构、沉积物渗透性以及岛屿气候条件密切相关[5]。国外对淡水透镜体研究集中在环礁基础上自然形成的岛屿—灰沙岛[6],对淡水透镜体的形成特征、影响因素和保护管理等方面进行了详细研究[3]。国内对淡水透镜体的研究工作集中在西沙永兴岛[4],利用原位观测、物理模型、数值模型等对永兴岛淡水储量以及淡水透镜体形态的时空变化、降雨和开采对淡水透镜体的影响进行了研究[7,8,9,10,11,12,13]。永兴岛淡水透镜体的数值模拟研究中主要针对达到动态稳定条件下的淡水透镜体形态,即最大可能形成的淡水透镜体规模,以及降雨量季节变化对其影响。对于人工岛,淡水透镜体形成是一个从无到有的渐变过程,对其形成时间和规模的影响因素分析是深入认识和合理规划利用南海珊瑚岛礁淡水资源的重要基础工作。
珊瑚礁的水文地质条件以及含水层渗透性是影响淡水透镜体形成的重要因素[5]。全新世碎屑沉积物的渗透性通常比礁灰岩低1~2个数量级。上下两层地质体的渗透性差异一方面造成地下水流线在非整合面上发生折射,另一方面造成礁灰岩中的咸淡水混合,两种作用往往造成淡水透镜体在非整合面处发生尖灭[6]。西沙永兴岛的不整合面位于海平面以下20 m处[7, 14],文献显示永兴岛珊瑚砂淡水透镜体的数值模型历经几次修正,早期模型中全新世珊瑚砂渗透系数根据粒径由经验公式给出,取值为50 m/d [15]。后据室内实验和抽水实验将渗透系数调整为70 m/d,淡水透镜体最大厚度为16.9 m [7]。永兴岛珊瑚砂非完整井抽水实验得到渗透系数平均值为110 m/d [16],利用岛上观测井地下水位对模型参数进行识别后采用渗透系数为100 m/d,并且认为含水层是各向同性的（水平向和垂向渗透系数相同）,模拟得到淡水透镜体最大厚度为15.0 m [14]。可见,全新世珊瑚砂渗透性是影响淡水透镜体形成的重要因素。
中国在珊瑚礁人工岛建设过程中,针对珊瑚砂渗透性已开展大量实验室和原位测定工作[17,18,19,20,21,22,23,24,25],对珊瑚砂颗粒粒径、级配以及微生物固化等因素对珊瑚砂的渗透性进行了系统研究。前人对岛屿淡水透镜体的研究多针对天然岛屿,多采用稳态模型模拟淡水透镜体达到稳定状态后可能形成的最大淡水透镜体体积及其影响因素,就珊瑚砂渗透性对天然海岛淡水透镜体最大深度的影响已进行了比较系统的研究[26]。但对于人工岛,填岛完工后多长时间能够开始形成淡水透镜体?淡水透镜体的最大规模如何?这些仍是人工岛内渗地下水淡化研究需要回答的问题。为合理预测人工岛淡水透镜体规模和形成时间,本文在前人研究基础上,以南海永暑礁人工岛为研究对象,构建人工岛淡水透镜体形成演化模型,重点研究珊瑚砂和礁灰岩渗透性对淡水透镜体形成时间和规模的影响,为中国加速生态岛礁建设进程提供理论支撑。

2 南海珊瑚礁地貌地质特征

中国南海除少数基岩岛以外大多数为珊瑚岛,南沙群岛是由一群珊瑚礁及部分礁坪上岛屿所组成,珊瑚礁类型主要为环礁[27]。环礁是热带海洋中的一种巨大珊瑚礁,能够独自在海洋中形成巨大礁体,并且在礁体上形成众多低平岛屿——珊瑚岛[28]。环礁一般分布在珊瑚易于生长的太平洋和印度洋的热带和亚热带海洋中,礁体呈一圈环状,中间为封闭或半封闭的澙湖,礁外缘多成陡坡。环礁在海面附近形成广大礁盘,其上可形成沙洲和灰沙岛。环礁在南海分布广泛、类型众多,自20世纪80年代开始前人已对南海珊瑚礁地貌进行了考察和系统研究[29]。中国南海环礁从其地貌形态可分为典型环礁、残缺环礁、沉没环礁、开放环礁和环礁链等六大类[28]。
永暑礁是中国南沙群岛西北部的一座开放型环礁。永暑礁从水深2000 m的海底拔起,礁顶长轴约25 km,呈NEE-SWW向,短轴约6 km,呈NW-SE向,面积约110 km2。礁顶中部为一开放式澙湖,面积105 km2,水深20~30 m。该珊瑚礁礁坪非常平坦,外礁坪略高于内礁坪,礁体西南向海坡坡度几近90°,呈垂直墙状态[30],但礁体陡峭的向海坡不是从礁缘外开始的,而是在礁缘外存在水深6~8 m的第一级浅水平台和水深36~38 m的第二级浅水平台,越过第二级浅水平台外缘,向海坡坡度变陡,水深陡然增大（图1）。礁坪在高潮时淹没于水下,低潮时露出。
图1 永暑礁地形[31]

注：等值线为海水等深线(m)。

Fig. 1 Topography of Yongshu Reef with bathymetric contours

珊瑚礁的二元地质结构在世界范围内具有普遍性,其形成与第四纪以来的气候变化和海平面升降有关。末次冰期时海平面下降,珊瑚礁暴露在海平面以上接受降雨淋滤和风浪侵蚀作用,全新世以来,海平面上升,珊瑚礁继续增长,沉积了10~20 m厚的松散未胶结的碎屑沉积层,侵蚀面以下为更新世及以前经历成岩作用胶结较好的礁灰岩,形成上土下岩的二元地质结构[32]。礁灰岩岩体结构保留了原生生物骨架中的孔隙,并且在礁体暴露时期经受喀斯特化作用,溶蚀孔洞发育,具有孔隙率高,渗透性大的特点[33]。这种珊瑚礁的二元地质结构造成全新世松散沉积物和礁灰岩存在显著的渗透性差异,对于人工岛淡水透镜体的形成过程有着重要影响。
根据1990年中国科学院南沙综合科学考察队在永暑礁西南礁坪实施的“南永1井”钻孔,珊瑚礁全新世碎屑松散沉积物厚度17.3 m,为未胶结的珊瑚砾块和生物砾砂屑,以下为更新世礁灰岩[34]（图2）。全新世地层也存在分层现象,造成含水层渗透性存在明显的各向异性。人工岛建于西南礁坪,用珊瑚砂填筑沙岛,填筑高度约为3 m。人工岛礁吹填过程中,颗粒自重对钙质珊瑚砂不同粒径颗粒起到分选作用,每次吹填过程形成表层一定厚度的低渗透性的粉土层,多次间断进行的吹填作业便造成了多层粉土层呈间断分层[35]。人工岛浅层钙质珊瑚砂的垂向分层结构和细颗粒粉土层对淡水透镜体形成具有重要影响,人工岛特殊的土层结构是人工岛礁淡水透镜体形成机制的重要研究方面。
图2 永暑礁西南礁坪浅层结构示意图[1, 36]

Fig. 2 Schematic diagram of geological structure of the shallow layers in southwest reef flat of the Yongshu Reef

3 珊瑚礁淡水透镜体形成数值模拟

3.1 变密度水流数学模型

珊瑚岛淡水透镜体的形成实质上是降雨入渗后,由于岛屿和海平面间的地势差,入渗淡水在重力作用下驱替珊瑚礁中原有内渗海水的过程,淡水和海水之间的边界由于扩散作用是一个有一定宽度的咸淡水过渡带,过渡带内地下水盐分浓度随淡水向海水过渡深度增加而增大。因此,淡水透镜体形成过程的模拟需要考虑礁内水体密度变化。
采用SEAWAT软件建立三维变密度流模型,SEAWAT耦合了MODFLOW和MT3DMS软件,分别用于变密度流运动和溶质运移微分方程的求解。人工岛淡水透镜体形成模拟采用三维变密度流模拟,当渗透系数张量主方向与坐标轴方向一致时,变密度流运动微分方程为：
∂ ∂ x i ρ μ 0 μ K i ∂ h ∂ x i + ρ - ρ 0 ρ 0 ∂ z ∂ x i = ρ S s ∂ h ∂ t + θ ∂ ρ ∂ C ∂ C ∂ t - ρ s q s
式中：h为淡水相应的地下水头（L）;ρ为地下水密度（M/L3）;μ为地下水粘度（M/(L·T)）;μ0为淡水粘度（M/(L·T)）ρ0为淡水密度（M/L3）;Ki为渗透系数张量在x, y, z方向的分量（Kx, Ky, Kz; L/T）;Ss为贮水率（1/L）;θ为有效孔隙度;C为地下水盐分浓度（M/L3）;t为时间（T）;qs为源汇项（1/T）;ρs为源汇项流体密度（M/L3）;i = 1, 2, 3,对应于x, y, z坐标。
溶质运移微分方程为：
1 + ρ b K d θ ∂ ( θC ) ∂ t = ∂ x i θ D ii ∂ C ∂ x i - ∂ x i qC - q s C s
式中：ρb为含水层介质干密度（M/L3）;Kd为吸附常数（L3/M）;Dii为弥散系数张量分量;q为达西流速（L/T）,Cs为源汇项浓度（M/L3）;其他符号意义同上式。

3.2 模型范围及时空离散

模型模拟范围在礁缘向海外延500~1000 m,水平方向将第二级浅水平台包含在内。模型空间坐标系x轴与礁顶长轴NEE-SWW方向一致,西南近礁缘边界至东北澙湖边界间模型长4.2 km,y轴与礁顶短轴NW-SE方向一致,模型宽1.6 km,模型平面面积为6.72 km2,平面上等间距剖分为32行,84列,网格大小为50 m×50 m（图3）。垂向上从人工岛地平面向下延伸至海平面以下50 m,垂向上剖分为51列,海平面以上人工岛部分剖分为1层,以下剖分为50层。模型模拟50年内淡水透镜体的形成过程。
图3 人工岛模型网格与剖面线位置

Fig. 3 Model grid and location of the cross section

3.3 边界条件及模型参数

模型中顶部边界人工岛范围内接受降雨入渗补给,根据永暑礁海洋观测站降雨观测资料,1989—2000年多年平均降雨量1940 mm。太平洋珊瑚岛地下水补给量由于植被覆盖度不同,一般为降雨量的30%~60% [37],考虑到永暑礁人工岛建岛初期植被覆盖低,补给量按照降雨量的50%计算,即0.97 m/a。顶部边界人工岛以外部分以及模型周边侧向边界均处理成定水头（0 m,按海水计）和定浓度边界（TDS浓度35 g/L）。模型底部边界处理成隔水边界。
永暑礁珊瑚砂孔隙体积与固体颗粒体积之比,即孔隙比范围在0.75~1.26之间[38],但珊瑚砂的孔隙由颗粒间的外孔隙和颗粒本身的内孔隙两部分组成。本文中取珊瑚砂孔隙度为0.44。根据太平洋中珊瑚岛含水层贮水参数和弥散参数研究[6, 37],含水层贮水率取0.0001（1/m）,纵向弥散度取1 m,横向和垂向弥散度均取0.01 m。珊瑚岛含水层弥散参数通常采用数值模型识别方法得到,但永暑礁现有监测数据尚不足以用于模型参数识别,本研究重点分析珊瑚礁二元地质结构的渗透性对淡水透镜体的影响。

4 全新世珊瑚碎屑沉积物及更新世礁灰岩渗透性

4.1 全新世珊瑚碎屑沉积物渗透性

珊瑚礁全新世松散沉积物的渗透性通常与环礁所在地的盛行风向有关,环礁迎风向通常沉积较粗的沉积物,渗透性较大,反之在背风向易沉积较细的沉积物,渗透性较差。并且在岛内部存在渗透性差异,澙湖一侧沉积物较海洋一侧往往颗粒较细,渗透性较低。这种全新世沉积物在岛内部的渗透性差异甚至可以达到两个数量级,造成淡水透镜体形态上明显不对称[6]。国外文献中全新世珊瑚的水平向渗透系数一般为10~100 m/d量级[6]。中国以往对珊瑚砂的研究集中在其工程性质上,近年来开始对其渗透性质进行研究,主要采用实验室测定和原位抽水实验方法获得渗透系数,尚未有岛内部珊瑚砂渗透性明显变化的报道（表1）。
表1 珊瑚砂渗透系数实验成果表

Tab. 1 Hydraulic conductivity of coral sands

测定方法 试样个数/实验组数 渗透系数范围(m/d) 渗透系数平均值(m/d) 数据来源
室内实验 常水头渗透实验 7 535.3~1334.1 706.6 [39]
常水头渗透实验 5 0.3~1.5 0.9 [22]
常水头渗透实验 1 1023.8 [19]
常水头渗透实验 1 58.6 [17]
原位实验 双环渗透实验 12 1.0~52.0 [21]
多孔稳定流抽水实验 3 0.7~1.6 1.1 [20]
多孔稳定流抽水实验 4 103.1~114.4 110.0 [16]
珊瑚砂的渗透性与颗粒粗细、形状和分选性有关。实验室测定过程中主要分析珊瑚砂颗粒级配与渗透性的关系,一般首先对珊瑚砂进行筛分后,绘制颗粒级配曲线,从级配曲线读取珊瑚砂的有效粒径,并计算级配曲线的不均匀系数和曲率系数。珊瑚砂渗透系数一般与孔隙比、曲率系数、不均匀系数和有效粒径的二次方有很好的相关性,前人据此已建立多个经验公式,使用不均匀系数、曲率系数和有效粒径计算渗透系数[21,22]。前人利用室内渗透实验获取的南沙人工岛礁填筑用珊瑚砂的渗透系数差异较大,这可能与珊瑚砂中细颗粒粒径和含量有关。珊瑚砂中细颗粒对渗透性影响实验表明,当最小粒径≤ 0.075 mm时,渗透系数量级为10-2 cm/s（10 m/d量级）;当最小粒径介于0.075~0.5 mm时,渗透系数量级为10-1 cm/s（100 m/d量级）[17]。当细颗粒含量较小时,其在水力渗透作用下挣脱束缚而发生运移、丢失,使土体结构重组,渗透系数一般呈缓慢增大后趋于稳定的变化特征[17]。
珊瑚礁原位抽水实验工作较少,新建人工岛浅层抽水实验测得的渗透系数在1 m/d量级,而天然的成岛历史久远的永兴岛浅层抽水实验测定渗透系数在100 m/d量级。室内实验和原位实验测定的珊瑚砂渗透系数差异较大,说明珊瑚砂渗透性与当地地层结构特征有关,人工岛珊瑚砂渗透性可能随雨水入渗淡化过程发生变化。因此,预测某具体人工岛能够形成的淡水透镜体规模时,有必要对珊瑚砂渗透性的影响进行定量分析。

4.2 更新世礁灰岩渗透性

前人对珊瑚岛淡水透镜的模拟研究中,礁灰岩的渗透系数多采用500~1000 m/d [6],并且一般采用经验值直接用于模型或者通过模型校正确定礁灰岩渗透系数取值。但对于新建人工岛,已有监测数据尚不足以用于模型的校正和参数识别。因此,需要有礁灰岩渗透性的实测成果用于人工岛淡水透镜体的模拟预测。中国目前在海底油气藏研究中,对礁灰岩渗透性进行了一定研究,确定礁灰岩储层渗透性的方法主要包括室内岩芯测定法、测井法和地震反演法,这3种方法测定的渗透率反映的孔隙介质的流体渗流性质,或者为孔隙渗透率[40]。对于孔隙、裂隙和溶洞可能并存的礁灰岩,上述3种方法测定的渗透率都不能反映礁灰岩的真实渗透率。中国南海流花11-1油田礁灰岩储层渗透性研究表明,礁灰岩渗透率在1~10 μm2范围[40]。
本文建立32个系列模型,分别针对更新世礁灰岩渗透系数（KHP）在500 m/d和1000 m/d尺度下,全新世沉积物渗透系数（KHH）在5~500 m/d范围内的量级变化对淡水透镜体（地下水总溶解固体浓度,TDS ≤ 1 g/L）的影响。垂向渗透系数按照水平渗透系数的1/5取值。

5 人工岛礁淡水透镜体形成控制因素

5.1 含水层渗透性对淡水透镜体形成时间影响

全新世珊瑚砂和更新世礁灰岩不同渗透系数组合情况下的淡水透镜体开始形成时间、最大厚度和体积如表2所示。人工岛淡水透镜体的形成是一个随时间逐渐增大直至稳定,淡水体与周围海水逐渐平衡的过程（图4）。全新世珊瑚碎屑沉积物渗透系数对开始出现淡水体的时间影响不大,填岛完成2年左右即可观测到淡水,形成微小的淡水透镜体（图4）。人工岛淡水体形成初期,淡水体的形成主要与地下水补给这一外部驱动力有关。太平洋中珊瑚岛遭受风暴潮淹没后的淡水透镜体的恢复过程的监测和模拟研究发现,海水淹没对珊瑚礁淡水透镜体的破坏是暂时的,淡水透镜体能够在6~17个月时间恢复[41,42],说明珊瑚岛淡水透镜体的形成时间很快。永暑礁人工岛吹填工程2015年基本完毕,经过2年左右降雨淋滤,2017年地下水氯离子浓度监测显示在珊瑚礁体中心浅层位置已形成圆台型淡水体[20]。模型模拟的淡水透镜体开始形成时间与南海人工岛礁地下水盐分变化监测结果基本一致。人工岛淡水透镜体形成时间较大陆含水层浅层潜水补给更新时间尺度非常短暂。这说明人工岛淡水透镜体在合理开发利用下,对保证人工岛生活和生态用水具有重要意义的。
表2 全新世珊瑚砂和更新世礁灰岩不同渗透系数组合下淡水透镜体形成时间和体积

Tab. 2 Formation time and maximum volume of fresh groundwater affected by different hydraulic conductivity of Holocene coral sand and Pleistocene limestone

全新世珊瑚砂渗透系数(KHH)(m/d) 更新世礁灰岩渗透系数(KHP)为500 m/d 更新世礁灰岩渗透系数(KHP)为1000 m/d
淡水透镜体开始形成时间(a) 淡水透镜体最大厚度(m) 淡水透镜体体积(万m3) 淡水透镜体开始形成时间(a) 淡水透镜体最大厚度(m) 淡水透镜体体积(万m3)
5 1.64 22.5 4774.75 1.72 19.5 4463.25
10 1.47 22.5 3969.25 1.84 20.5 3103.5
20 1.7 20.5 2850.25 1.75 19.5 2670.75
30 1.54 21.5 2492.25 1.64 19.5 2378.75
40 1.52 20.5 2257.25 1.53 19.5 2134
50 1.51 20.5 2081.5 1.53 18.5 1964.5
60 1.51 20.5 1963 1.52 18.5 1829.75
70 1.47 19.5 1835.5 1.48 18.5 1769.75
80 1.49 19.5 1715.25 1.49 19.5 1699.75
90 1.46 18.5 1599 1.46 18.5 1587.75
100 1.43 17.5 1502.25 1.43 17.5 1487.25
150 1.44 14.5 1148.75 1.43 14.5 1137.75
200 1.42 12.5 958 1.45 12.5 945.75
300 1.43 8.5 633.5 1.4 7.5 554.75
400 1.43 6.5 387.75 1.43 6.5 389.25
500 1.45 4.5 194.5 1.42 5.5 179.5

注：a淡水入渗2年后,开始形成淡水透镜体;b淡水入渗8年,淡水透镜体的咸淡水界面进入礁灰岩;c淡水入渗15年,淡水透镜体达到最大体积;剖面位置见图3。

图4 淡水透镜体体积随时间变化规律

Fig. 4 Variation of fresh groundwater volume with time

5.2 含水层渗透性对淡水透镜体形成体积影响

淡水透镜体最大厚度随全新世珊瑚碎屑沉积物渗透系数增大呈指数衰减（图5）,与前人研究发现规律一致[16]。永暑礁全新世碎屑松散沉积物与更新世礁灰岩不整合面深度为17.3 m,当全新世珊瑚砂渗透系数小于100 m/d时,淡水透镜体最大深度越过不整合面进入礁灰岩。当淡水透镜体进入礁灰岩后,礁灰岩渗透系数增大同样造成淡水透镜体体积一定程度的减小（表2）。当全新世珊瑚砂渗透系数大于100 m/d时,礁灰岩渗透性对淡水透镜体的最大厚度和体积无显著影响。
图5 淡水透镜体最大厚度随全新世珊瑚砂渗透系数变化

Fig. 5 Variation of calculated maximum thickness of fresh groundwater lens with hydraulic conductivities of Holocene sediments

模拟各情景中,淡水透镜体进入礁灰岩深度不超过2 m（表2）,礁灰岩渗透性对淡水透镜体体积影响作用有限,说明南海珊瑚礁二元地质结构中,对淡水透镜体起控制作用的是全新世珊瑚砂的渗透性。
相比淡水透镜体最大厚度随渗透系数衰减,淡水透镜体形成体积随渗透系数的衰减更加显著。淡水透镜体形成体积随全新世珊瑚碎屑沉积物渗透系数增大呈对数衰减（即渗透系数取对数后,呈线性衰减趋势）（图6）。若以永兴岛浅层抽水实验测定渗透系数（100 m/d量级）为基准,能够形成的淡水透镜体体积在1500万m3左右,这对于人工岛生活用水保障和生态岛屿建设有重要意义。
图6 淡水透镜体体积随全新世珊瑚砂渗透系数变化

Fig. 6 Variation of calculated volume of fresh groundwater with hydraulic conductivities of Holocene sediments

5.3 珊瑚礁二元地质结构对淡水透镜体形态影响

虽然岛礁的二元地质结构中全新世珊瑚砂对淡水透镜体最终能够形成的体积起控制作用,但当全新世珊瑚砂渗透系数在10 m/d量级以下时,淡水透镜体的最大深度浓度均超过了全新世珊瑚砂和更新世礁灰岩的不整合面,更新世礁灰岩对淡水透镜体的形成过程会造成影响。当淡水透镜体的咸淡水界面越过不整合面进入更新世礁灰岩后,由于含水层渗透系数突然增大,淡水透镜体的咸淡水界面向水平向扩展,形如“髭”状[43]（图7）。珊瑚礁二元地质结构不仅造成最终形成的淡水透镜体在不整合面发生尖灭,对于淡水透镜体的动态演化过程也有重要影响。
图7 沿A-A'剖面不同时间淡水透镜体分布特征

Fig. 7 Spatial distribution of salinity concentration along cross section A-A' at selected times during the formation process of the fresh groundwater lens.

永暑礁人工岛氯离子浓度监测显示当入渗淡水进入渗透性较大地层后向水平向扩展,形成“帽沿状”淡水体[20]。人工岛氯离子监测虽然目前仍限于全新世松散沉积物,但渗透性的分层差异对淡水透镜体形成过程的影响观测和模拟结果是一致的。淡水体在人工岛下部主要沿垂向向下扩展,进入礁灰岩顶界后沿水平向有所扩大,但这种水平向的扩展距离应该是有限的。
模拟结果和盐分浓度监测说明不仅珊瑚礁上土下岩的地质结构分层对淡水透镜体的形成过程和形态有重要影响,全新世沉积物和人工岛吹填土层的分层同样对淡水透镜体有重要影响。当考虑含水层存在内部分层造成的渗透性各向异性时,淡水透镜体的形成过程会更加复杂。

6 讨论及结论

对南海珊瑚礁人工岛淡水透镜体的模拟预测是在目前认识和资料条件下的、针对珊瑚礁普遍存在的二元地质结构,建立了珊瑚礁人工岛淡水驱替海水形成淡水透镜体的动态模拟模型,系统分析了地层渗透性对淡水透镜体体积的影响,给出了人工岛礁淡水形成时间和淡水透镜体可能形成的体积规模,根据珊瑚礁二元地质结构不同渗透系数组合下淡水驱替海水的非稳态模拟结果,得出如下结论：
（1）人工岛礁二元地质结构、上层全新世沉积物和下层礁灰岩渗透性对淡水开始形成时间无明显影响,人工岛填岛完成2年左右可观测到淡水,在10~20年内可形成相对稳定的淡水透镜体。由于南海珊瑚礁热带气候下的强降雨量,人工岛淡水透镜体的形成时间较大陆含水层地下水的补给更新时间尺度是非常短暂的,其对人工岛生活和生态用水是有重要保证意义的。
（2）全新世沉积物渗透性是控制人工岛礁淡水透镜体体积的主要因素,能够形成的淡水透镜体最大体积随渗透系数增大而呈对数减小,并且淡水透镜体体积随全新世沉积物渗透性增大而减小的速率较淡水透镜体厚度更加显著。渗透系数小于100 m/d时,淡水透镜体最大深度将越过不整合面,但进入礁灰岩深度有限（< 2 m）,淡水透镜体的体积进而受礁灰岩渗透性影响,随礁灰岩渗透系数增大而减小。基于数值模型对人工岛淡水透镜体形成体积的可靠预测,一方面需要对含水层渗透性空间分布进行深入调查,另一方面也需要对地下水盐分浓度进行详细观测,能够利用其对模型参数进行针对性校正。
（3）淡水透镜体越过不整合面进入更新世礁灰岩顶界后,由于含水层渗透系数突然增大,淡水透镜体咸淡水界面向水平向扩展,形成“髭”状淡水体。人工岛地层的渗透性分层差异对淡水透镜体形成过程的影响值得深入研究,关键在于对珊瑚砂含水层垂向渗透系数沿深度变化的详细准确刻画。

参考文献

原文顺序 | 文献年度倒序 | 文中引用次数倒序
[1]
Zhao Huanting, Wang Lirong . Construction of artificial islands on coral reef in the South China Sea Islands. Tropical Geography, 2017,37(5):681-693.

DOI

[ 赵焕庭, 王丽荣 . 南海诸岛珊瑚礁人工岛建造研究. 热带地理, 2017,37(5):681-693.]

[2]
Zhao Huanting, Wang Lirong, Song Chaojing . Geomorphological model of coral reefs in the South China Sea. Acta Oceanologica Sinica, 2014,36(9):112-120.

DOI

[ 赵焕庭, 王丽荣, 宋朝景 . 南海珊瑚礁地貌模型研究. 海洋学报, 2014,36(9):112-120.]

[3]
Zhao Huanting, Wang Lirong . Review on the study of freshwater lens in the coral reef island. Tropical Geography, 2015,35(1):120-129.

[ 赵焕庭, 王丽荣 . 珊瑚礁岛屿淡水透镜体研究综述. 热带地理, 2015,35(1):120-129.]

[4]
Zhao Huanting, Wang Lirong, Song Chaojing . Review on freshwater lens of lime-sand island in Nanhai Zhudao. Marine Science Bulletin, 2014,33(6):601-610.

[ 赵焕庭, 王丽荣, 宋朝景 . 南海诸岛灰沙岛淡水透镜体研究述评. 海洋通报, 2014,33(6):601-610.]

[5]
Cao Yingjie . Numerical simulation of the formation and evolution of freshwater lens in islands[D]. Nanjing: Hohai University, 2008.

[ 曹英杰 . 海岛淡水透镜体形成和演化规律的数值模拟研究[D]. 南京: 河海大学, 2008.]

[6]
Werner A D, Sharp H K, Galvis S C , et al. Hydrogeology and management of freshwater lenses on atoll islands: Review of current knowledge and research needs. Journal of Hydrology, 2017,551:819-844.

DOI

[7]
Zhou Congzhi, He Li, Yang Qin , et al. Three-dimensional numerical simulation of freshwater len in coral islands. Journal of Hydraulic Engineering, 2010,41(5):560-566.

[ 周从直, 何丽, 杨琴 , 等. 珊瑚岛礁淡水透镜体三维数值模拟研究. 水利学报, 2010,41(5):560-566.]

[8]
Fang Zhendong, Zhou Congzhi, Liang Hengguo , et al. Research on influential factors upon the up-cone of the freshwater lens in a coral island. Journal of Logistical Engineering University, 2012,28(4):57-66.

[ 方振东, 周从直, 梁恒国 , 等. 珊瑚岛淡水透镜体抽水倒锥影响因素研究. 后勤工程学院学报, 2012,28(4):57-66.]

[9]
Zhou Congzhi, Fang Zhendong, Liang Henguo , et al. Influence of rainfall on freshwater lens in a coral island. China Water & Wastewater, 2006,22(1):53-57.

[ 周从直, 方振东, 梁恒国 , 等. 雨量变化对珊瑚岛礁淡水透镜体的影响. 中国给水排水, 2006,22(1):53-57.]

[10]
Sheng Chong, Xu Hehua, Zhang Wentao . Numerical simulation of the effect of geomorphologic changes on freshwater lens in the Yongxing Island. Hydrogeology and Engineering Geology, 2018,45(6):7-14.

[ 盛冲, 许鹤华, 张文涛 . 地貌变化对永兴岛淡水透镜体影响的数值模拟. 水文地质工程地质, 2018,45(6):7-14.]

[11]
Zhen Li . Physical simulation of formation and evolution of freshwater lens in island[D]. Nanjing: Hohai University, 2008.

[ 甄黎 . 海岛淡水透镜体形成与演化规律的物理模拟[D]. 南京: 河海大学, 2008.]

[12]
Li Daning, Xu Jiexin, Zhao Hui , et al. Temporal variations and influencing factors of freshwater lens in Yongxing Island. Journal of Guangdong Ocean University, 2019,39(2):67-74.

[ 黎大宁, 许洁馨, 赵辉 , 等. 永兴岛淡水透镜体的时间变化规律及其影响因素. 广东海洋大学学报, 2019,39(2):67-74.]

[13]
Zhao Jun . Analysis on the evolution law of freshwater lens in island under mining condition[D]. Nanjing: Hohai University, 2007.

[ 赵军 . 开采条件下海岛淡水透镜体演变规律分析[D]. 南京: 河海大学, 2007.]

[14]
Zhou Congzhi, Qiao Hua, Du Rong . Simulation and exploitation of the freshwater lens in coral island. Journal of Logistical Engineering University, 2016,32(3):1-10.

[ 周从直, 谯华, 杜蓉 . 珊瑚岛淡水透镜体的模拟与开发. 后勤工程学院学报, 2016,32(3):1-10.]

[15]
Zhou Congzhi, Fang Zhendong, Liang Hengguo , et al. Numerical simulation of freshwater lens on coral island. Marine Sciences, 2004,28(11):77-80.

[ 周从直, 方振东, 梁恒国 , 等. 珊瑚岛礁淡水透镜体的数值模拟. 海洋科学, 2004,28(11):77-80.]

[16]
Zhou Congzhi, Fang Zhendong, Wei Ying , et al. Development and Utilization of Freshwater Lenses on Coral Reefs. Chongqing: Chongqing University Press, 2017.

[ 周从直, 方振东, 魏营 , 等. 珊瑚岛礁淡水透镜体的开发利用. 重庆: 重庆大学出版社, 2017.]

[17]
Hu Mingjian, Cui Xiang, Wang Xinzhi , et al. Experimental study on the effect of fine particles on permeability of the calcareous sand. Rock and Soil Mechanics, 2019,40(8):1-6.

[ 胡明鉴, 崔翔, 王新志 , 等. 细颗粒对钙质砂渗透性的影响试验研究. 岩土力学, 2019,40(8):1-6.]

[18]
Zhao Xiaowei . Experimental study on the effect of freshwater seawater mixing on the distribution and permeability of coral sand [D]. Beijing: China University of Geosciences (Beijing), 2018.

[ 赵孝伟 . 咸淡水混合对珊瑚砂溶蚀及渗透特性影响试验研究[D]. 北京: 中国地质大学(北京), 2018.]

[19]
Ren Yubin, Wang Yin, Yang Qing . Effects of particle size distribution and shape on permeability of calcareous sand. Rock and Soil Mechanics, 2018,39(2):491-497.

[ 任玉宾, 王胤, 杨庆 . 颗粒级配与形状对钙质砂渗透性的影响. 岩土力学, 2018,39(2):491-497.]

[20]
Jiang Hanghai . Study on the permeability and characteristics of the groundwater quality of corel reef[D]. Nanning: Guangxi Univeristy, 2018.

[ 蒋航海 . 珊瑚礁砂渗透性与地下水质变化特征研究[D]. 南宁: 广西大学, 2018.]

[21]
Hu Mingjian, Jiang Hanghai, Zhu Zhangqi , et al. Discussion on permeability of calcareous sand and its influencing factors. Rock and Soil Mechanics, 2017,38(10):2895-2900.

[ 胡明鉴, 蒋航海, 朱长歧 , 等. 钙质砂的渗透特性及其影响因素探讨. 岩土力学, 2017,38(10):2895-2900.]

[22]
Qian Kun, Wang Xinzhi, Chen Jianwen , et al. Experimental study on permeability of calcareous sand for islands in the South China Sea. Rock and Soil Mechanics, 2017,38(6):1557-1564, 1572.

[ 钱琨, 王新志, 陈剑文 , 等. 南海岛礁吹填钙质砂渗透特性试验研究. 岩土力学, 2017,38(6):1557-1564, 1572.]

[23]
Ren Yubin . Experimental study on the permeability characteristics of calcareous sand in South China Sea[D]. Dalian: Dalian University of Technology, 2016.

[ 任玉宾 . 南海钙质砂渗透特性试验研究[D]. 大连: 大连理工大学, 2016.]

[24]
Ou Yixi, Fang Xiangwei, Zhang Nan , et al. Influence of solution salinity on microbial biocementation of coral sand. Journal of Logistical Engineering University, 2016,32(1):78-82.

[ 欧益希, 方祥位, 张楠 , 等. 溶液盐度对微生物固化珊瑚砂的影响. 后勤工程学院学报, 2016,32(1):78-82.]

[25]
Ou Yixi, Fang Xiangwei, Shen Chunni , et al. Influence of particle sizes of coral sand on biocementation. Journal of Water Resources and Architectural Engineering, 2016,14(2):35-39.

[ 欧益希, 方祥位, 申春妮 , 等. 颗粒粒径对微生物固化珊瑚砂的影响. 水利与建筑工程学报, 2016,14(2):35-39.]

[26]
Underwood M R, Peterson F L, Voss C I . Groundwater lens dynamics of atoll islands. Water Resources Research, 1992,28(11):2889-2902.

DOI

[27]
Zhao Huanting, Yuan Jiayi . The landscape and regional palaeogeography of Nansha Islands. Acta Geographica Sinica, 1995,50(2):107-117.

[ 赵焕庭, 袁家义 . 南沙群岛景观及区域古地理. 地理学报, 1995,50(2):107-117.]

[28]
Ceng Zhaoxuan . Types of atolls in China. Marine Science Bulletin, 1982,4(4):43-50.

[ 曾昭璇 . 中国环礁的类型划分. 海洋通报, 1982,4(4):43-50.]

[29]
Zhao Liang, Zhang Zhengsheng, Nan Wenlong . Zeng Zhaoxuan's contributions to the geographical research of the South China Sea. Acta Geographica Sinica, 2016,71(3):515-523.

DOI

[ 赵亮, 张争胜, 南文龙 . 曾昭璇先生对南海地理研究的重要贡献. 地理学报, 2016,71(3):515-523.]

[30]
Yu Hongbing, Sun Zongxun, Mao Qingwen . Geomorphic features of upper seaward slope of Yongshu Reef, Nansha Islands. Marine Science Bulletin, 1999,18(3):49-54.

[ 于红兵, 孙宗勋, 毛庆文 . 南沙群岛永暑礁上部向海坡地貌特征研究. 海洋通报, 1999,18(3):49-54.]

[31]
Sun Zongxun, Zhan Wenhuan, Zhu Junjiang . Rockmass stability of Yongshu Reef, Nansha Islands. Marine Geology & Quaternary Geology, 2003,23(4):9-14.

[ 孙宗勋, 詹文欢, 朱俊江 . 南沙群岛永暑礁岩体稳定性评价. 海洋地质与第四纪地质, 2003,23(4):9-14.]

[32]
Liu Jingru, Wang Jianping, Guo Dong , et al. Three-dimensional morphology and stratigraphic structure of coral reef in Xisha and Nasha Islands. Advances in Geosciences, 2019,9(2):61-68.

DOI

[ 刘景儒, 王建平, 郭东 , 等. 西南沙珊瑚礁三维形态和地层结构特征研究. 地球科学前沿, 2019,9(2):61-68.]

[33]
Tang Guoyi, Zheng Jianguo . Engineering properties of reef calcareous rock in Southeast Asia. Geotechnical Investigation & Surveying, 2015,43(6):6-10.

[ 唐国艺, 郑建国 . 东南亚礁灰岩的工程特性. 工程勘察, 2015,43(6):6-10.]

[34]
Zhao Huanting, Zhu Yuanzhi, Nie Baofu , et al. Preliminary division of geological age and Quaternary strata of Yongshu reef. Chinese Science Bulletin, 1992,37(23):2165-2168.

[ 赵焕庭, 朱袁智, 聂宝符 , 等. 永暑礁地质年代和第四纪地层初步划分. 科学通报, 1992,37(23):2165-2168.]

[35]
Wang Xinzhi, Wang Xing, Hu Mingjian , et al. Study of permeability of calcareous silty layer of foundation at an artificial reclamation island. Rock and Soil Mechanics, 2017,38(11):3127-3135.

[ 王新志, 王星, 胡明鉴 , 等. 吹填人工岛地基钙质粉土夹层的渗透特性研究. 岩土力学, 2017,38(11):3127-3135.]

[36]
Sun Zongxun, Zhan Wenhuan, Zhu Junjiang . Rockmass structure and engineering geological zones of coral reefs, Nansha Islands of South China Sea. Journal of Tropical Oceanography, 2004,23(3):11-20.

[ 孙宗勋, 詹文欢, 朱俊江 . 南沙群岛珊瑚礁岩体结构特征及工程地质分带. 热带海洋学报, 2004,23(3):11-20.]

[37]
Post V E A, Bosserelle A L, Galvis S C , et al. On the resilience of small-island freshwater lenses: Evidence of the long-term impacts of groundwater abstraction on Bonriki Island, Kiribati. Journal of Hydrology, 2018,564:133-148.

DOI

[38]
Sun Zongxun . Engineering properties of coral sands in Nansha Islands. Tropic Oceanology, 2000,19(2):1-8.

[ 孙宗勋 . 南沙群岛珊瑚砂工程性质研究. 热带海洋, 2000,19(2):1-8.]

[39]
Shu Longcang, Zhou Congzhi, Zhen Li , et al. Measurement of the hydrological properties of coral sand in a laboratory. Journal of Hohai University (Natural Sciences), 2008,36(3):330-332.

[ 束龙仓, 周从直, 甄黎 , 等. 珊瑚砂含水介质水理性质的实验室测定. 河海大学学报(自然科学版), 2008,36(3):330-332.]

[40]
Hou Lianhua, Wu Xiling, Lin Chengyan , et al. Determination of permeability in ecologic reef reservoir. Acta Petrolei Sincia, 2003,24(5):67-70, 73.

[ 侯连华, 吴锡令, 林承焰 , 等. 礁灰岩储层渗透率确定方法. 石油学报, 2003,24(5):67-70, 73.]

[41]
Terry J P, Falkland A C . Responses of atoll freshwater lenses to storm-surge overwash in the Northern Cook Islands. Hydrogeology Journal, 2010,18(3):749-759.

DOI

[42]
Bailey R T . Quantifying transient post-overwash aquifer recovery for atoll islands in the western Pacific. Hydrological Processes, 2015,29(20):4470-4482.

DOI

[43]
Stoeckl L, Houben G J, Dose E J . Experiments and modeling of flow processes in freshwater lenses in layered island aquifers: Analysis of age stratification, travel times and interface propagation. Journal of Hydrology, 2015,529(6851):159-168.

DOI

Options
文章导航

/

〈 〉