• 全国中文核心期刊
  • 中国科技核心期刊
  • 中国科学引文数据库收录期刊
  • 世界期刊影响力指数(WJCI)报告来源期刊
  • Scopus, CA, DOAJ, EBSCO, JST等数据库收录期刊

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于多变量统计分析和水化学特征的海水入侵特征研究

刚什婷 吕明荟 卢茜茜 高帅 赵志强 陈奂良 彭同强 王玺 邢立亭 李莉霞

刚什婷,吕明荟,卢茜茜,等. 基于多变量统计分析和水化学特征的海水入侵特征研究−以青岛市崂山区为例[J]. 中国岩溶,2023,42(5):982-994 doi: 10.11932/karst20230509
引用本文: 刚什婷,吕明荟,卢茜茜,等. 基于多变量统计分析和水化学特征的海水入侵特征研究−以青岛市崂山区为例[J]. 中国岩溶,2023,42(5):982-994 doi: 10.11932/karst20230509
GANG Shenting, LYU Minghui, LU Qianqian, GAO Shuai, ZHAO Zhiqiang, CHEN Huanliang, PENG Tongqiang, WANG Xi, XING Liting, LI Lixia. Characterization of seawater intrusion based on multivariate statistical analysis and water chemistry characteristics: A case study of Laoshan district, Qingdao City[J]. CARSOLOGICA SINICA, 2023, 42(5): 982-994. doi: 10.11932/karst20230509
Citation: GANG Shenting, LYU Minghui, LU Qianqian, GAO Shuai, ZHAO Zhiqiang, CHEN Huanliang, PENG Tongqiang, WANG Xi, XING Liting, LI Lixia. Characterization of seawater intrusion based on multivariate statistical analysis and water chemistry characteristics: A case study of Laoshan district, Qingdao City[J]. CARSOLOGICA SINICA, 2023, 42(5): 982-994. doi: 10.11932/karst20230509

基于多变量统计分析和水化学特征的海水入侵特征研究——以青岛市崂山区为例

doi: 10.11932/karst20230509
基金项目: 山东省地矿局八〇一水文地质工程地质大队(暨山东省地下水环境保护与修复工程技术研究中心)基金项目(801KY2021-4)
详细信息
    作者简介:

    刚什婷(1990-),女,硕士研究生,主要从事地下水污染控制与数值模拟研究。E-mail:gangshenting@163.com

    通讯作者:

    吕明荟(1982-),女,工程师,主要从事水文地质研究工作。E-mail:lmhyoxiang@163.com

  • 中图分类号: P641.3

Characterization of seawater intrusion based on multivariate statistical analysis and water chemistry characteristics: A case study of Laoshan district, Qingdao City

  • 摘要: 沿海地区地下水环境问题日益突出,进行地下水水化学特征及演化规律的研究,能够更有效地开展地下水环境的监测和保护。以青岛市崂山区地下水为研究对象,综合运用统计分析、主成分分析、Piper图解法、HFE-D图解法、Chadha’s矩形图法等方法,对研究区海水入侵特征与地下水化学特征演化进行分析,探究地下水水化学特征及演化规律,并进一步评价了海水入侵现状。结果表明,研究区地下水以Na+、Ca2+、ClSO24为主要优势离子,地下水化学类型多为Cl·SO4-Na型和SO4·Cl-Ca·Mg型。地下水中Cl浓度变化幅度较大,且其均值超出了有无海水入侵的分界值(250 mg·L−1),地下水可能发生一定程度的海水入侵;青岛市崂山区地下水呈中性至弱碱性(pH均值=7.0~8.0),是沿海地区长期的水文地球化学过程的影响;地下水化学变化主要受自然因素(岩石与水的相互作用)或人为因素(农业和家庭活动)的控制;采用反距离加权(IDW)方法,结合地理信息系统(GIS),进行海水入侵位置的空间映射,研究结果表明崂山区海水入侵主要分布于江家土寨东−浦里社区北入侵段,王哥庄−港西−港东入侵段、仰口湾入侵段、登瀛村−栲栳岛入侵段。

     

  • 海水入侵是一种在海平面上升等自然因素下或者大量开采内陆地下水等人类活动影响下,地下淡水与海水之间的平衡状态遭到破坏,引起海水和地下淡水交界面向内陆含水层推进的现象[1-3]。海水入侵是全世界沿海含水层面临的主要环境问题[4-5]。海水入侵会破坏淡水资源、改变当地生态地质环境、腐蚀地下设施、造成土壤盐渍化等问题。为防止海水入侵,国内外学者在海水入侵机理、影响因素、预测分析和防治措施等领域开展了大量的研究[6-8]

    沿海含水层的咸化可能是由于简单的、直接的海水入侵造成的,也可能涉及一系列复杂的地球化学过程,例如,岩石风化、蒸发浓缩、补径排、离子交换、氧化还原作用和人类活动等,这些过程以不同的方式影响水质变化[9]。因此利用水化学指标来分析水化学特征,探究水化学演化规律及影响因素是研究沿海含水层咸化的化学过程以及盐度的来源最常见有效的方法[10]。目前,对于水化学过程研究,最广泛使用的水文地球化学技术是主要、次要和微量元素的离子比率、离子比值法、混合计算和地球化学模型[11-16]

    青岛市的海(咸)水入侵起源于20世纪70年代后期,到20世纪80年代中期,青岛地区发生连续多年干旱,地下水严重超采,海水入侵形成迅速蔓延趋势,导致区域水资源贫乏,尤其基岩地区贫水最甚,海水入侵严重威胁着青岛市崂山区弥足珍贵的地下水资源 [17-25]。本研究以青岛市崂山区地下水为主要研究对象,研究:①研究区水文地球化学特征;②取样井的主要离子化学特性及影响因素;③地下水的水文地球化学演变过程;④对海水入侵现状进行评价。旨在为研究区及类似沿海地区地下水水化学特征研究以及生态环境保护提供科学依据。

    研究区位于青岛市崂山区,北起土寨河北岸,与即墨市海岸接壤;南至麦岛,与市南区海岸相接;东至海岸线;自海岸线向内陆延伸3~5 km为西部边界,研究区位置见图1。研究区地貌主要为游冲积平原及山前冲洪积平原,部分为滨海养殖滩涂地带,面积约245 km2。属温带大陆季风气候区,四季分明,雨量较为丰富,多年平均降水量828.3 mm,多年平均蒸发量1 461.1 mm,多年平均气温12.1 ℃。境内河流水系主要为沿海近缘水系,共有23条大小河流,呈放射状展布。河流特点是河长短、流速快、多直流入海、属于季节性河流。

    图  1  研究区位置图
    Figure  1.  Location of the study area

    研究区出露的地层较简单,从老至新依次为:元古界、中生界白垩系和新生界第四系。根据研究区含水介质类型和埋藏条件,将地下水划分为三大类:松散岩类孔隙水、碎屑岩类孔隙−裂隙水和基岩裂隙水(图2图3)。第四系孔隙含水岩组为主要含水层,主要分布在研究区河流中下游的河谷平原,该区域渗透能力较好,含水量随含水层厚度的增加而增大,其单井出水量在100~500 m3·d−1之间。碎屑岩类孔隙−裂隙水仅赋存于研究区东南的黄山−青山村一带,含水层透水性、富水性均较弱。

    图  2  青岛崂山区水文地质略图
    Figure  2.  Hydrogeological sketch of Laoshan district, Qingdao
    图  3  土寨河流域水文地质剖面图A-A′
    Figure  3.  Hydrogeological profile(A-A′) of the Tuzhai river basin

    本研究采用的数据均来自青岛市崂山区水文地质调查数据(2016—2017年),所有水样检测分析均在山东省地矿工程勘察院实验室完成,测试内容主要包括K+、Na+、Ca2+、Mg2+、NH+4HCO3、ClSO24NO3、可溶性总固体(TDS)、总碱度(TA)、总硬度(TH)、pH等指标。采样点位置如图4。采用Origin2017软件对水化学参数进行描述性统计分析,绘制了Piper三线图、Chadha’s矩形图、Chadha’s矩形图,并运用主成分分析法讨论研究区地下水主要影响因素。借助ACRGIS软件空间分析功能对研究区浅层地下水海水入侵现状的空间分布规律进行分析,结合自然条件及人类活动背景资料,揭示研究区地下水化学特征及成因。

    图  4  取样监测点位置图
    Figure  4.  Locations of sampling and monitoring sites

    主成分分析法经常被用来确定地下水质量的主要影响因子[26-27]。这种方法在许多研究中被证实是非常有效的[28-32]。利用水化学分析方法研究地下水特征和演化有助于识别地下水的水文地球化学过程。Piper三线图是运用最广泛的水化学分析方法之一,但该方法只能表示地下水的化学组成,不能表示地下水化学演化。Giménez−Forcada提出了水化学相演化图HFE−D,根据海水阳离子/阴离子的组成,可以有效地检查海水入侵的位置和淡水/海水的混合百分比[33-34]

    HFE−D的计算过程涉及的阳离子有Na+、K+、Ca2+、Mg2+,阴离子有HCO3、CO32–、ClSO24。海水入侵的水文地球化学过程复杂,HFE−D在简化复杂过程时忽略了一些次要过程,因此HFE−D的绘制过程只考虑占比最大的Na+、K+、Ca2+HCO3、ClSO24。其中K+、Mg2+不在图中展示,仅在计算中使用。另外HCO3SO24都可以代表淡水,在图中一个指标表示即可。本文选取HCO3作为淡水中的特征阴离子。HFE−D中横纵坐标均用各阴阳离子在总阴、总阳离子的浓度百分比表示。HFE−D把含水层发生海水入侵时的水化学演化分为侵入期和恢复期,分别表示咸化过程和淡化过程。

    一般来说,Na+和Cl是海水/盐水的主要离子、而Ca2+HCO3通常是淡水的主要离子[35]。因此,沿海地下水中高浓度的Na+和Cl离子可能表明海水混合的显著影响和盐水的出现,而大量的Ca2+HCO3主要反映了水−岩石相互作用的贡献[36]。通过绘制Chadha’s矩形图可以确定地下水盐化的来源。Chadha’s矩形图中阴阳离子的值都为相对质量浓度,阴阳离子的和分别为100%。大矩形可以分为四个大区,分别代表①补给水;②反向离子交换水;③海水入侵;④离子交换水。每个大区分为4个次级的矩形,总计16个次一级的矩形,每一个矩形都以c(Ca2++Mg2+)−c(Na+)为X轴,c(HCO3)−c(SO24+Cl)为Y轴。

    海水入侵地下水质量指数[37]GQISWI由Na+、K+、Ca2+、Mg2+HCO3、ClSO24七种离子浓度计算得到,该指数实现了把海水入侵过程中一系列复杂的水文化学过程量化。在GQISWI图中可以把Piper三线图划分成I到VI区(图5),并给各分区分别赋值,实现把Piper三线图结果量化。

    图  5  GQISWI对Piper图的分区
    Figure  5.  Domains of GQISW in piper diagram

    对地下水样品检测结果的主要指标进行数理统计分析。表1中给出各水化学指标的最小值、最大值、平均值和标准偏差。表中各数据,尤其是TDS、TH、Na+、Cl等水质指标,变化较大,表明研究区内地下水化学组分有很大差异。地下水中阳离子以Na+为主,浓度均值关系为Na+> Ca2+> Mg2+> K+;阴离子浓度Cl> SO24>HCO3。pH是一个可以反映水文地球化学平衡信息的重要因子,结果表明,研究区地下水pH均值处于7.0~8.0之间,变异系数为0.04,说明该地区的地下水呈中性至弱碱性,在区域变化较小。表中各数据,尤其是TDS、TH、Na+、Cl等水质指标,变化较大,表明研究区内地下水化学组分有很大差异。地下水中Cl浓度变化幅度较大,且其均值超出了《地下水质量标准》(GB/T 14848−2017)[38]的限值(250 mg·L−1)和海水入侵的分界值(250 mg·L−1),说明地下水可能发生一定程度的海水入侵,且基岩裂隙水中Cl浓度远低于第四系孔隙水浓度,说明降雨使地下水Cl浓度降低。TDS质量浓度变化比较大,最大值为17 138.1 mg·L−1,根据TDS含量分类[39],地下水存在淡水、微咸水、咸水,部分地区地下水为盐水。

    表  1  地下水水化学参数统计特征值(单位:mg·L−1,pH除外)
    Table  1.  Statistics of hydrochemical parameters of groundwater (unit: mg·L−1, except for pH)
    分区项目pHTDSTHCa2+Mg2+K+Na+HCO3SO24ClNO3
    基岩裂
    隙水
    Mean7.34998.81380.6094.2334.711.78199.04104.33135.62385.9435.94
    SD0.272032.16670.59200.6846.612.49512.0182.42218.811103.4532.57
    Cv0.042.031.762.131.341.402.570.791.612.860.91
    Min6.9067.2146.1212.933.360.132.177.653.0521.530
    Max8.208 877.643 064.85910.97191.908.822 125.00306.08778.464 812.71101.35
    第四系
    孔隙水
    Mean7.201 435.83554.2688.3581.0412.92277.95108.77235.23560.8369.69
    SD0.383 211.69802.7263.12163.2736.36927.5675.61403.371 722.6367.98
    Cv0.052.241.450.712.012.813.340.701.713.070.98
    Min6.5074.8339.208.314.480.176.6715.3010.6813.130.36
    Max8.5017 138.103 597.12263.18728.01200.005 000.00369.841 617.989 160.70244.75
    地表水Mean7.30301.99138.7140.888.904.4330.6985.0756.1140.9123.25
    SD0.32302.26127.0635.809.185.7541.8290.2662.8153.3833.27
    Cv0.041.000.920.881.031.301.361.061.121.301.43
    Min6.8091.6948.5415.432.310.164.4327.9819.217.185.52
    Max7.60814.84351.78100.8024.3112.59103.60243.46165.70132.7882.60
    注:Min为最小值,Max为最大值,Mean为平均值,SD为标准差,Cv为变异系数。
    Note: Min represents minimum value; Max represents maximum value; Mean represents average value, SD represents standard deviation; Cv represents variation coefficient.
    下载: 导出CSV 
    | 显示表格

    在进行主成分分析时,有9个变量被考虑在内:Na、Mg、Ca、K、Cl、SO4、HCO3、NO3、pH。由 Total Variance Explained(主成分特征根和贡献率)可知,特征根λ1=4.76,特征根λ2=1.725,前两个主成分的累计方差贡献率达72.06%,即涵盖了大部分信息。表明前两个主成分能够代表最初的9个指标来分析地下水水化学情况,故提取前两个主成分即可。主成分分别记作PC1、PC2。正如预期的Na+、K+、Ca2+、Mg2+、ClSO24在PC1有较高的荷载(表2图6),相关性强,反映的是水−岩作用对水化学的影响,如蒸发、岩石与水的相互作用,以及可能与海水入侵有关。为证实这一假设,研究Piper三线图、HFE−D图、Chadha’s矩形图等其他相关指标,pH、HCO3NO3在第二主成分上有较高载荷,相关性强,说明发生了碳酸盐矿物的溶解;NO3为负值(NO3代表了农业活动和碱性矿物肥料的主要特征),反映了农业活动对地下水的影响,也可能与降水直接补给有关。由于农业和人类活动对青岛市地下水质量的影响前人已有大量研究[40,25],本文只关注海水入侵对地下水的影响,因此对NO3不做过多研究。

    表  2  主成分分析法组成矩阵
    Table  2.  Matrix formed by principal component analysis
    PC1PC2
    pH0.007 700.641 80
    Ca2+0.270 27−0.106 40
    Mg2+0.431 21−0.056 97
    K+0.396 81−0.004 19
    Na+0.444 050.006 05
    HCO30.099 170.482 26
    SO240.414 96−0.072 80
    Cl0.447 43−0.016 36
    NO3−0.063 99−0.579 08
    下载: 导出CSV 
    | 显示表格
    图  6  主成分分析结果图:PC1、PC2
    Figure  6.  Results of principal component analysis: PC1 and PC2
    3.2.1   Piper三线图解析水质演化

    Piper三线图能反映地下水主要离子组成变化和水化学特征,可以用来划分水化学类型。对于阳离子(图7左下三角形),大部分数据中主导阳离子是Ca2+,平均毫克当量百分比为75%,碱土金属元素(Ca2++Mg2+)浓度大于碱金属元素(K++Na+),仅在少数水样中碱金属元素(K++Na+)大于碱土金属元素(Ca2++Mg2+)。对于阴离子(图7右下三角形),大部分数据主导阴离子是HCO3,平均毫克当量百分比为75%,弱酸根离子(HCO3)浓度大于强酸根离子(Cl+SO24)。研究区地层中常见的矿物有方解石(CaCO3)、白云石(CaMg(CO3)2)等碳酸盐矿物,石膏(CaSO4·2H2O)等硫酸盐矿物,由此初步判断地下水化学组分受碳酸盐矿物和硫酸盐矿物溶滤作用较强。

    图  7  地下水水化学Piper图
    Figure  7.  Piper diagram of groundwater hydrochemsitry

    按照GQISWI中对Piper图中水化学类型的划分,崂山区地下水多为Cl·SO4-Na型和SO4·Cl-Ca·Mg型地下水,部分地区为HCO3-Ca型,还有少数监测点地下水为Cl-Na型水组成。HCO3-Ca型是淡水的主要特征,Cl-Na型是海水的主要特征,Cl-Ca·Mg和HCO3-Na型是咸淡水之间过渡的水化学类型。图7表明取样点位置只有少部分地区地下水是HCO3-Ca型淡水,多数地区地下水中Cl含量较大,并且有部分区域的地下水化学类型为Cl-Na型,上述研究结果表明,部分地区已经发生了严重的海水入侵现象,且多数地区地下水化学类型处于咸淡水之间的过渡阶段。

    地下水化学类型分布较一致,河流中上游阳离子主要以Ca2+和Mg2+为主,阴离子主要以HCO3为主,由陆地向海洋,主要阳离子逐渐向Na+变化,阴离子向Cl变化。图中取样点水化学类型由地势高的地方向海洋趋近,水化学类型由HCO3-Ca型逐渐向Cl-Na型演化,可能沿径流路径发生了阳离子交换,以及海水不同程度的混合作用。如在石人河,LS015、LS011、LS009、LS010从上游到海岸边依次排列,水化学类型分别由HCO3·Cl-Ca·Mg→SO4·Cl-Ca·Mg→Cl-Ca·Mg→Cl·SO4-Na变化;在土寨河,LS012、LS005、LS001、LS003、LS010从上游到海岸边依次排列,水化学类型分别由SO4-Ca·Mg→SO4·Cl-Na·Ca·Mg→SO4·Cl-Ca·Mg→SO4Cl-Na·Mg→SO4·Cl-Na变化;在王哥庄河LS019、LS020、LS016、LS033、LS035从上游到海岸水化学类型由HCO3·Cl-Ca·Mg→HCO3·Cl-Na·Ca·Mg→HCO3·Cl-Na·Ca·Mg→Cl-Ca·Mg→Cl-Na·Mg变化;在南部的大河东—小河东LS085、LS081、LS082、LS092、LS089、LS090从上游到海岸边依次排列,水化学类型分别由HCO3·Cl-Na·Ca·Mg→HCO3-Na·Ca→Cl-Na·Ca→SO4·Cl-Na·Ca·Mg→SO4·Cl-Na·Ca·Mg→Cl-Na·Ca·Mg变化。根据各取样点离子浓度关系,绘制了地下水七大离子浓度占比图(图8)。由图8可知,地下水中Cl含量占主导地位,阳离子以Ca2+、Na+为主。

    图  8  采样点地下水七大离子浓度占比图
    Figure  8.  Percentages of concentrations of the seven ions in groundwater of the sampling sites
    3.2.2   HFE−D图解译水文地球化学过程

    崂山区地下水HFE−D如图9所示,显示了研究区地下水部分采样点为Cl-Na型海水,表明河流中下游沿海地区地下水发生了海水入侵。部分采样点位于混合线以下,表明该区域地下水正处于海水侵入阶段,并且多数地区地下水已经完成由HCO3-Ca型向Cl-Ca型的演化。HFE−D图隐藏了坐标轴0~33.3%段各离子浓度相对均衡的区域,用“Mix”表示图中主导离子浓度占比小于50%的水化学类型[33]

    图  9  崂山区地下水HFE−D
    Figure  9.  Evolution of hydrochemical facies of groundwater in Laoshan district
    3.2.3   Chadha’s矩形图解析水质演化

    根据Chadha’s图(图10),受海水入侵影响的钻孔(LS010、LS029、LS082号样品)位于海岸线上,分别对应江家土寨东−浦里社区北入侵段、王山口南−桑园社区北−王哥庄−港东入侵段、登瀛村南−栲栳岛入侵段。区域②、③表示c(Ca2++Mg2+)>c(Na+),其中c(Na+)<25%,25%≤c(SO24+ Cl)≤75%,此区内地下水化学类型为HCO3·SO4·Cl-Ca·Mg型或SO4·HCO3-Ca·Mg型水(采样点LS009、LS019、LS030、LS040、LS023、LS024、LS102、LS110、LS081、LS022、LS046、LS060、LS074、LS015)。

    图  10  Chadha矩形水化学类型图
    Figure  10.  Chadha rectangle hydrochemical diagram

    区域④表示c(Ca2++ Mg2+)>c(Na+)和c(HCO3)<c(SO24+ Cl),并且c(Na+)和c(HCO3)相对质量浓度均小于25%,这种水有永久硬度,用于灌溉时,不会形成NaCO3盐沉积,地下水化学类型为SO4·Cl-Ca·Mg型水(采样点LS001、LS089、LS111、LS033)。

    采用基于熵权的集对分析法和GIS 软件结合,选取了Cl、矿化度、SO24、γClHCO3、SAR五项化学特征指标作为评价因子,评价因子等级划分指标见表3,对崂山区53个监测井点海水入侵程度分别进行了评价,并绘制了海水入侵现状图(图11)。由图11可以看出,青岛市崂山区地下水综合质量较好,受海水入侵影响较少,海水入侵主要分布于沿海冲洪积平原、河流下游及河道两侧区域,尤其是河流入海口。自沿海地区向内陆地区海水入侵现状逐渐减缓,海水入侵主要分布于江家土寨东−浦里社区北入侵段,王哥庄−港西−港东入侵段、仰口湾入侵段、登瀛村−栲栳岛入侵段。

    表  3  海水入侵指标的等级划分 (单位/mg·L−1
    Table  3.  Classification of indexes of seawater intrusion (unit/mg·L−1)
    特征因子IIIIIIIV
    无入侵轻度入侵中度入侵严重入侵
    Cl ≤250 ≤600 ≤1 500 >1 500
    SO24 ≤200 ≤450 ≤1 200 >1 200
    M ≤1 000 ≤2 000 ≤3 000 >3 000
    SAR ≤2 ≤3.55 ≤10 >10
    γCl/γHCO3 ≤0.5 ≤1.0 ≤6.6 >6.6
    下载: 导出CSV 
    | 显示表格
    图  11  青岛市崂山区海水入侵空间分布
    Figure  11.  Spatial distribution of seawater intrusion in Laoshan district, Qingdao City

    本文以2016—2017年地下水水质调查数据为背景,以地下水中主要离子浓度为水化学指标,对青岛市崂山区地下水化学特征及海水入侵状况进行分析,得出以下结论:

    (1)水化学指标统计分析结果显示TDS、TH、Na+、Cl等水质指标变化较大,表明研究区内地下水化学组分有很大差异。地下水中Cl浓度变化幅度较大,且其均值超出了《地下水质量标准》(GB/T 14848−2017)的限值(250 mg·L−1)和海水入侵分界值(250 mg·L−1),说明地下水可能发生一定程度的海水入侵。根据主成分分析可知控制研究区地下水水质的因素主要有大气降水和人为活动。

    (2)Piper三线图、HFE−D、Chadha’s矩形图结果显示,取样点位置只有少部分地区地下水是HCO3-Ca型淡水,多数地区地下水中Cl含量较大,并且有部分区域的地下水化学特征是Cl-Na型,表明部分地区已经发生了严重的海水入侵现象,多数地区地下水化学类型处于咸淡水之间的过渡阶段。研究区地下水海水入侵有向内陆扩展的趋势。

    (3)根据海水入侵评价结果显示,海水入侵水主要分布于江家土寨东−浦里社区北入侵段,王哥庄−港西−港东入侵段、仰口湾入侵段、登瀛村−栲栳岛入侵段。

  • 图  1  研究区位置图

    Figure  1.  Location of the study area

    图  2  青岛崂山区水文地质略图

    Figure  2.  Hydrogeological sketch of Laoshan district, Qingdao

    图  3  土寨河流域水文地质剖面图A-A′

    Figure  3.  Hydrogeological profile(A-A′) of the Tuzhai river basin

    图  4  取样监测点位置图

    Figure  4.  Locations of sampling and monitoring sites

    图  5  GQISWI对Piper图的分区

    Figure  5.  Domains of GQISW in piper diagram

    图  6  主成分分析结果图:PC1、PC2

    Figure  6.  Results of principal component analysis: PC1 and PC2

    图  7  地下水水化学Piper图

    Figure  7.  Piper diagram of groundwater hydrochemsitry

    图  8  采样点地下水七大离子浓度占比图

    Figure  8.  Percentages of concentrations of the seven ions in groundwater of the sampling sites

    图  9  崂山区地下水HFE−D

    Figure  9.  Evolution of hydrochemical facies of groundwater in Laoshan district

    图  10  Chadha矩形水化学类型图

    Figure  10.  Chadha rectangle hydrochemical diagram

    图  11  青岛市崂山区海水入侵空间分布

    Figure  11.  Spatial distribution of seawater intrusion in Laoshan district, Qingdao City

    表  1  地下水水化学参数统计特征值(单位:mg·L−1,pH除外)

    Table  1.   Statistics of hydrochemical parameters of groundwater (unit: mg·L−1, except for pH)

    分区项目pHTDSTHCa2+Mg2+K+Na+HCO3SO24ClNO3
    基岩裂
    隙水
    Mean7.34998.81380.6094.2334.711.78199.04104.33135.62385.9435.94
    SD0.272032.16670.59200.6846.612.49512.0182.42218.811103.4532.57
    Cv0.042.031.762.131.341.402.570.791.612.860.91
    Min6.9067.2146.1212.933.360.132.177.653.0521.530
    Max8.208 877.643 064.85910.97191.908.822 125.00306.08778.464 812.71101.35
    第四系
    孔隙水
    Mean7.201 435.83554.2688.3581.0412.92277.95108.77235.23560.8369.69
    SD0.383 211.69802.7263.12163.2736.36927.5675.61403.371 722.6367.98
    Cv0.052.241.450.712.012.813.340.701.713.070.98
    Min6.5074.8339.208.314.480.176.6715.3010.6813.130.36
    Max8.5017 138.103 597.12263.18728.01200.005 000.00369.841 617.989 160.70244.75
    地表水Mean7.30301.99138.7140.888.904.4330.6985.0756.1140.9123.25
    SD0.32302.26127.0635.809.185.7541.8290.2662.8153.3833.27
    Cv0.041.000.920.881.031.301.361.061.121.301.43
    Min6.8091.6948.5415.432.310.164.4327.9819.217.185.52
    Max7.60814.84351.78100.8024.3112.59103.60243.46165.70132.7882.60
    注:Min为最小值,Max为最大值,Mean为平均值,SD为标准差,Cv为变异系数。
    Note: Min represents minimum value; Max represents maximum value; Mean represents average value, SD represents standard deviation; Cv represents variation coefficient.
    下载: 导出CSV

    表  2  主成分分析法组成矩阵

    Table  2.   Matrix formed by principal component analysis

    PC1PC2
    pH0.007 700.641 80
    Ca2+0.270 27−0.106 40
    Mg2+0.431 21−0.056 97
    K+0.396 81−0.004 19
    Na+0.444 050.006 05
    HCO30.099 170.482 26
    SO240.414 96−0.072 80
    Cl0.447 43−0.016 36
    NO3−0.063 99−0.579 08
    下载: 导出CSV

    表  3  海水入侵指标的等级划分 (单位/mg·L−1

    Table  3.   Classification of indexes of seawater intrusion (unit/mg·L−1)

    特征因子IIIIIIIV
    无入侵轻度入侵中度入侵严重入侵
    Cl ≤250 ≤600 ≤1 500 >1 500
    SO24 ≤200 ≤450 ≤1 200 >1 200
    M ≤1 000 ≤2 000 ≤3 000 >3 000
    SAR ≤2 ≤3.55 ≤10 >10
    γCl/γHCO3 ≤0.5 ≤1.0 ≤6.6 >6.6
    下载: 导出CSV
  • [1] Andersen M S, Nyvang V, Jakobsen R, Postma D. Geochemical processes and solute transport at the seawater/freshwater interface of a sandy aquifer[J]. Geochimica et Cosmochimica Acta, 2005, 69: 3979–3994.
    [2] Werner A D, Simmons C T. Impact of sea-level rise on seawater intrusion in coastal aquifer[J]. Groundwater, 2009, 47:197-204. doi: 10.1111/j.1745-6584.2008.00535.x
    [3] Liu Y, Jiao J J, Liang W, Kuang X. Hydrogeochemical characteristics in coastal groundwater mixing zone[J]. Applied Geochemistry, 2017, 85:49-60. doi: 10.1016/j.apgeochem.2017.09.002
    [4] Kazakis N, Pavlou A, Vargemezis G, Voudouris K S, Soulios G, Pliakas F, Tsokas G. Seawater intrusion mapping using electrical resistivity tomography and hydrochemical data. An application in the coastal area of eastern Thermaikos Gulf, Greece[J]. Science of the Total Environment, 2016, 543:373-387. doi: 10.1016/j.scitotenv.2015.11.041
    [5] Alfarrah N, Walraevens K. Groundwater overexploitation and seawater intrusion in coastal areas of arid and semi-arid regions[J]. Water, 2018, 10(2):143.
    [6] Bovolo C I, Parkin G, Sophocleous M. Groundwater resources, climate and vulnerability[J]. Environmental Research Letters, 2009, 4(3):35001. doi: 10.1088/1748-9326/4/3/035001
    [7] 郭明华, 王敬, 戴长国, 蒋雷, 张浩清, 柳超. 文登区浅层地下水化学演化与海水入侵研究[J]. 海洋科学, 2021, 45(7):57-65. doi: 10.11759/hykx20200708002

    GUO Minghua, WANG Jing, DAI Changguo, JIANG Lei, ZHANG Haoqing, LIU Chao. Study on the hydrochemical evolution of groundwater and seawater intrusion in the shallow layer of Wendeng district[J]. Marine Science, 2021, 45(7):57-65. doi: 10.11759/hykx20200708002
    [8] 赵倩. 锦州小凌河扇地海水入侵过程中水文地球化学作用及影响因素研究[D]. 廊坊: 防灾科技学院, 2021.

    ZHAO Qian. Hydrogeochemistry and influencing factors in the process of seawater intrusion in Xiaoling river fan, Jinzhou City[D]. Langfang: Institute of Disaster Prevention, 2021.
    [9] Mondal NC, Singh VP, Singh VS, Saxena VK. Determining the interaction between groundwater and saline water through groundwater major ions chemistry[J]. Journal of Hydrology, 2010, 388, 100-111.
    [10] 王昊, 王中良. 天津北大港水库水质咸化的地球化学机理分析[J]. 天津师范大学学报(自然科学版), 2016, 36(6):29-34.

    WANG Hao, WANG Zhongliang. Geochemical analysis on water salinization mechanism of the Tianjin Beidagang reservoir[J]. Journal of Tianjin Normal University (Natural Science Edition), 2016, 36(6):29-34.
    [11] 辛祥. 大沽河下游地区地下水氯离子时空变化特征与调控的数值模拟[D]. 青岛: 青岛大学, 2019.

    XIN Xiang. Numerical simulation of spatial and temporal characteristics and control of chloride ion in groundwater in the downstream area of Dagu river[D]. Qingdao: Qingdao University, 2019.
    [12] 林曼曼. 青岛近海海域灾害地质特征研究[D]. 石家庄: 石家庄经济学院, 2014.

    LIN Manman. The study on the characteristics of geo-hazard factors in Qingdao offshore[D]. Shijiazhuang: University of Economics, 2014.
    [13] 陆求裕. 长乐滨海地区海水入侵现状及成因分析[J]. 地质灾害与环境保护, 2020, 31(4):101-105. doi: 10.3969/j.issn.1006-4362.2020.04.019

    LU Qiuyu. The present situation of the seawater intrusion along the coastal area of Changle[J]. Journal of Geological Hazards and Environment Preservation, 2020, 31(4):101-105. doi: 10.3969/j.issn.1006-4362.2020.04.019
    [14] 梁彬, 关碧珠, 崔光中, 揭江, 唐祚旺. 广东湛江开发巨型盆地防治海水入侵研究[J]. 中国岩溶, 2001, 20(1):77.
    [15] 秦凯凯. 管状海水入侵形成机理及影响因素研究[D]. 大连: 大连海事大学, 2017.

    QIN Kaikai. Study on formation mechanism and influencing factors of tubular seawater intrusion[D]. Dalian: Dalian Maritime University, 2017.
    [16] 崔震, 陈广泉, 徐兴永, 王东亮, 喻龙. 北长山岛海水入侵成因机理及现状评价[J]. 海洋环境科学, 2015, 34(6):930-936. doi: 10.13634/j.cnki.mes.2015.06.021

    CUI Zhen, CHEN Guangquan, XU Xingyong, WANG Dongliang, YU Long. Mechanism and assessment of seawater intrusion in the northern Changshan island[J]. Marine Environmental Science, 2015, 34(6):930-936. doi: 10.13634/j.cnki.mes.2015.06.021
    [17] 徐铭霜. 某沿海地区海水入侵动态变化规律及防治措施研究[D]. 青岛: 青岛理工大学, 2022.

    XU Mingshuang. Research on the dynamic variation rule of seawater intrusion and controlling measures in a coastal area[D]. Qingdao: Qingdao University of Technology, 2022.
    [18] 王蜜蕾, 窦衍光, 邹亮, 薛碧颖, 胡睿, 岳保静, 徐刚, 林曦, 李林森. 青岛崂山周边地下水化学特征与矿泉水成因分析[J]. 海洋地质前沿, 2021, 37(9):17-24. doi: 10.16028/j.1009-2722.2021.126

    WANG Milei, DOU Yanguang, ZOU Liang, XUE Biying, HU Rui, YUE Baojing, XU Gang, LIN Xi, LI Linsen. Hydrochemical characteristics of groundwater and genesis of mineral water at Laoshan mountain and surrounding areas, Qingdao[J]. Marine Geology Frontiers, 2021, 37(9):17-24. doi: 10.16028/j.1009-2722.2021.126
    [19] 胡慧君. 青岛市大沽河典型区非饱和−饱和带地下水耦合模型[D]. 南京: 南京大学, 2019.

    HU Huijun. An integrated modeling of variably saturated groundwater flow system in a typical area of Dagu river basin[D]. Nanjing: Nanjing University, 2019.
    [20] 刘久潭, 周丹, 高宗军, 王敏, 马媛媛, 张洪英, 时孟杰, 董杰. 青岛西海岸新区地下水水化学特征及水质评价[J]. 山东科技大学学报(自然科学版), 2019, 38(2):14-24, 43. doi: 10.16452/j.cnki.sdkjzk.2019.02.002

    LIU Jiutan, ZHOU Dan, GAO Zongjun, WANG Min, MA Yuanyuan, ZHANG Hongying, SHI Mengjie, DONG Jie. Hydrochemical characteristics and water quality assessment of groundwater in Qingdao West Coast New District[J]. Journal of Shandong University of Science and Technology (Natural Science), 2019, 38(2):14-24, 43. doi: 10.16452/j.cnki.sdkjzk.2019.02.002
    [21] 尹子悦, 林青, 徐绍辉. 青岛市大沽河流域地下水水化学时空演化及影响因素分析[J]. 地质论评, 2018, 64(4):1030-1043. doi: 10.16509/j.georeview.2018.04.017

    YIN Ziyue, LIN Qing, XU Shaohui. Spatial-temporal variations and controlling factors of groundwater hydrochemical characteristics in the Dagu river basin[J]. Geological Review, 2018, 64(4):1030-1043. doi: 10.16509/j.georeview.2018.04.017
    [22] 康晓雨. 青岛市典型富水区地下水环境质量评价及演化规律[D]. 石家庄: 河北地质大学, 2017.

    KANG Xiaoyu. Groundwater quality evaluation and evolution law of typical rich water area in Qingdao Hebei[D]. Shijiazhuang: Hebei GEO University, 2017.
    [23] 肖菲, 周勇华, 陈小英. 青岛沿岸海(咸)水入侵灾害的成因与现状[J]. 海洋地质前沿, 2012, 28(5):54-58. doi: 10.16028/j.1009-2722.2012.05.011

    XIAO Fei, ZHOU Yonghua, CHEN Xiaoying. Research on the origin and current status of sea (salt) water intrusion in coastal areas of Qingdao[J]. Marine Geology Frontiers, 2012, 28(5):54-58. doi: 10.16028/j.1009-2722.2012.05.011
    [24] 姜兴钰. 青岛新河地区海水入侵综合研究[D]. 青岛: 国家海洋局第一海洋研究所, 2010.
    [25] 董少杰, 孟春霞, 王成见. 青岛市地下水源地水质评价及污染原因分析[J]. 水资源与水工程学报, 2006(6):54-57. doi: 10.3969/j.issn.1672-643X.2006.06.013

    DONG Shaojie, MENG Chunxia, WANG Chengjian. Water quality evaluation and reasons analysis of pollution for groundwater sources in Qingdao[J]. Journal of Water Resources and Water Engineering, 2006(6):54-57. doi: 10.3969/j.issn.1672-643X.2006.06.013
    [26] Shaheen A, Iqbal J, Hussain S. Adaptive geospatial modeling of soil contamination by selected heavy metals in the industrial area of Sheikhupura, Pakistan[J]. International Journal of Environmental Science and Technology, 2019, 16(8): 4447-4464.
    [27] Khan M K, Ayoub W, Saied S, Hussain M M, Masood S S, Siddique A, Khwaja H A. Statistical and geospatial assessment of groundwater quality in the megacity of karachi[J]. Journal of Water Resource and Protection, 2019, 11(3): 311-332.
    [28] Fatma B B, Emna B, Jalila M, Salem B. Evaluation of groundwater hydrogeochemical characteristics and delineation of geothermal potentialities using multi-criteria decision analysis: Case of Tozeur region, Tunisia[J]. Applied Geochemistry, 2020, 113: 104504.
    [29] 周爱红, 牛鹏飞, 袁颖, 黄虎城. 基于PCA-PSO-SVM的凡口铅锌矿地区岩溶地表塌陷危险性预测[J]. 中国岩溶, 2020, 39(4):622-628.

    ZHOU Aihong, NIU Pengfei, YUAN Ying, HUANG Hucheng. Prediction of karst surface subsidence risk in the Fankou lead-zinc mine area based on PCA-PSO-SVM[J]. Carsologica Sinica, 2020, 39(4):622-628.
    [30] 叶慧君, 张瑞雪, 吴攀, 李学先, 覃应机, 査学芳, 韩志伟. 基于主成分分析的岩溶水水化学组成及影响因素研究: 以贵州水城盆地为例[J]. 中国岩溶, 2017, 36(2): 215-225.

    YE Huijun, ZHANG Ruixue, WU Pan, LI Xuexian, QIN Yingji, ZHA Xuefang, HAN Zhiwei. Hydrochemical characterization of groundwater and surface water and their influencing factors based on principal component analysis: An example in the Shuicheng basin of Guizhou[J]. Carsologica Sinica, 2017, 36(2): 215-225.
    [31] 赵艳丽, 王敏, 李常锁. 山东省淄河流域地下水化学污染源解析[J]. 人民长江, 2022, 53(12):37-43, 49. doi: 10.16232/j.cnki.1001-4179.2022.12.006

    ZHAO Yanli, WANG Min, LI Changsuo. Analysis of chemical pollution source of groundwater in Zihe river basin of Shandong Province[J]. Yangtze River, 2022, 53(12):37-43, 49. doi: 10.16232/j.cnki.1001-4179.2022.12.006
    [32] 张洁, 梁杏, 刘延锋, 张鑫, 孙立群, 赵枫, 付鹏宇. 基于主成分的协克里金法对地下水砷空间分布预测[J]. 地球科学, 2023, 48(10): 3820-3831.

    ZHANG Jie, LIANG Xing, LIU Yanfeng, ZHANG Xin, SUN Liqun, ZHAO Feng, FU Pengyu. The cokriging method is based on principal components to predict the spatial distribution of arsenic in groundwater[J]. Earth Science, 2023, 48(10): 3820-3831.
    [33] Giménez Forcada E. Dynamic of sea water interface using hydrochemical facies evolution diagram[J]. Groundwater, 2010, 48:212-216. doi: 10.1111/j.1745-6584.2009.00649.x
    [34] Giménez Forcada E. Space/time development of seawater intrusion: A study case in Vinaroz coastal plain (Eastern Spain) using HFE-Diagram, and spatial distribution of hydrochemical facies[J]. Journal of Hydrology, 2014, 517:617-627. doi: 10.1016/j.jhydrol.2014.05.056
    [35] Chadha D K. A proposed new diagram for geochemical classification of natural waters and interpretation of chemical data[J]. Hydrogeology Journal, 1999(7):431-439. doi: 10.1007/s100400050216
    [36] Mondal N C, Singh V P, Singh V S, Saxena V K. Determining the interaction between groundwater and saline water through groundwater major ions chemistry[J]. Journal of Hydrology, 2010, 388:100-111. doi: 10.1016/j.jhydrol.2010.04.032
    [37] Nosair A M, Shams M Y, AbouElmagd L M, Hassanein A E, Fryar A E, Abu Salem H S. A predictive model for progressive salinization in a coastal aquifer using artificial intelligence and hydrogeochemical techniques: A case study of the Nile Delta aquifer, Egypt[J]. Environmental Science and Pollution Research, 2022, 29(6):9318-9340. doi: 10.1007/s11356-021-16289-w
    [38] GB/T 14848-2017. 地下水质量标准[S].

    GB/T 14848-2017. Standard for Groundwater Quality[S].
    [39] Wang Min, Yin Yueping, Wen Dongguang. Handbook of hydrogeology[M]. Beijing: Geology Press, 2012.
    [40] 寇馨月, 丁军军, 李玉中, 毛丽丽, 李巧珍, 徐春英, 郑欠, 庄姗. 青岛市农区地下水硝态氮污染来源解析[J]. 环境科学, 2021, 42(7): 3232−3241.

    KOU Xinyue, DING Junjun, LI Yuzhong, MAO Lili, LI Qiaozhen, XU Chunying, ZHENG Qian, ZHUANG Shan. Identifying the sources of groundwater NO3 -N in agricultural region of Qingdao[J]. Environmental Science, 2021, 42(7): 3232−3241.
  • 期刊类型引用(1)

    1. 邸俊强,孙大明,原晖,王主华,吕宝阔. 基于GQI_(SWI)模型的大连市地下水海水入侵水质评价. 水电能源科学. 2025(02): 63-66 . 百度学术

    其他类型引用(1)

  • 加载中
图(11) / 表(3)
计量
  • 文章访问数:  237
  • HTML浏览量:  20
  • PDF下载量:  47
  • 被引次数: 2
出版历程
  • 收稿日期:  2023-04-20
  • 录用日期:  2023-08-11
  • 修回日期:  2023-08-07
  • 刊出日期:  2023-10-01

目录

/

返回文章
返回