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

基于多变量统计分析和水化学特征的海水入侵特征研究以青岛市崂山区为例刚什婷1,2，吕明荟1,2，卢茜茜1,2，高帅1,2，赵志强1,2，陈奂良1,2，彭同强1,2，王玺1,2，邢立亭3，李莉霞1,2（1.山东省地质矿产勘查开发局八一水文地质工程地质大队(山东省地矿工程勘察院),山东济南250014；2.山东省地下水环境保护与修复工程技术研究中心,山东济南250014；3.济南大学水利与环境学院,山东济南250022）SO24摘要：沿海地区地下水环境问题日益突出，进行地下水水化学特征及演化规律的研究，能够更有效地开展地下水环境的监测和保护。以青岛市崂山区地下水为研究对象，综合运用统计分析、主成分分析、Piper 图解法、HFE-D 图解法、Chadhas 矩形图法等方法，对研究区海水入侵特征与地下水化学特征演化进行分析，探究地下水水化学特征及演化规律，并进一步评价了海水入侵现状。结果表明，研究区地下水以 Na+、Ca2+、Cl、为主要优势离子，地下水化学类型多为 ClSO4-Na 型和 SO4Cl-CaMg 型。地下水中 Cl浓度变化幅度较大，且其均值超出了有无海水入侵的分界值（250mgL1），地下水可能发生一定程度的海水入侵；青岛市崂山区地下水呈中性至弱碱性（pH 均值=7.08.0），是沿海地区长期的水文地球化学过程的影响；地下水化学变化主要受自然因素（岩石与水的相互作用）或人为因素（农业和家庭活动）的控制；采用反距离加权（IDW）方法，结合地理信息系统（GIS），进行海水入侵位置的空间映射，研究结果表明崂山区海水入侵主要分布于江家土寨东浦里社区北入侵段，王哥庄港西港东入侵段、仰口湾入侵段、登瀛村栲栳岛入侵段。关键词：地下水；沿海含水层；水化学特征；地下水化学演化；海水入侵；青岛市崂山区中图分类号：P641.3文献标识码：A文章编号:10014810（2023）05098213开放科学（资源服务）标识码（OSID）：0引 言海水入侵是一种在海平面上升等自然因素下或者大量开采内陆地下水等人类活动影响下，地下淡水与海水之间的平衡状态遭到破坏，引起海水和地下淡水交界面向内陆含水层推进的现象13。海水入侵是全世界沿海含水层面临的主要环境问题45。海水入侵会破坏淡水资源、改变当地生态地质环境、腐蚀地下设施、造成土壤盐渍化等问题。为防止海水入侵，国内外学者在海水入侵机理、影响因素、预测分析和防治措施等领域开展了大量的研究68。沿海含水层的咸化可能是由于简单的、直接的基金项目：山东省地矿局八一水文地质工程地质大队（暨山东省地下水环境保护与修复工程技术研究中心）基金项目（801KY2021-4）第一作者简介：刚什婷（1990），女，硕士研究生，主要从事地下水污染控制与数值模拟研究。E-mail：。通信作者：吕明荟（1982），女，工程师，主要从事水文地质研究工作。E-mail：。收稿日期：20230420第42卷第5期中国岩溶Vol.42No.52023年10月CARSOLOGICASINICAOct.2023刚什婷，吕明荟，卢茜茜，等.基于多变量统计分析和水化学特征的海水入侵特征研究以青岛市崂山区为例J.中国岩溶，2023，42（5）：982-994.DOI：10.11932/karst20230509海水入侵造成的，也可能涉及一系列复杂的地球化学过程，例如，岩石风化、蒸发浓缩、补径排、离子交换、氧化还原作用和人类活动等，这些过程以不同的方式影响水质变化9。因此利用水化学指标来分析水化学特征，探究水化学演化规律及影响因素是研究沿海含水层咸化的化学过程以及盐度的来源最常见有效的方法10。目前，对于水化学过程研究，最广泛使用的水文地球化学技术是主要、次要和微量元素的离子比率、离子比值法、混合计算和地球化学模型1116。青岛市的海(咸)水入侵起源于 20 世纪 70 年代后期，到 20 世纪 80 年代中期，青岛地区发生连续多年干旱，地下水严重超采，海水入侵形成迅速蔓延趋势，导致区域水资源贫乏，尤其基岩地区贫水最甚，海水入侵严重威胁着青岛市崂山区弥足珍贵的地下水资源1725。本研究以青岛市崂山区地下水为主要研究对象，研究：研究区水文地球化学特征；取样井的主要离子化学特性及影响因素；地下水的水文地球化学演变过程；对海水入侵现状进行评价。旨在为研究区及类似沿海地区地下水水化学特征研究以及生态环境保护提供科学依据。1研究区概况研究区位于青岛市崂山区，北起土寨河北岸，与即墨市海岸接壤；南至麦岛，与市南区海岸相接；东至海岸线；自海岸线向内陆延伸 35km 为西部边界，研究区位置见图 1。研究区地貌主要为游冲积平原及山前冲洪积平原，部分为滨海养殖滩涂地带，面积约 245km2。属温带大陆季风气候区，四季分明，雨量较为丰富，多年平均降水量 828.3mm，多年平均蒸发量 1461.1mm，多年平均气温 12.1。境内河流水系主要为沿海近缘水系，共有 23 条大小河流，呈放射状展布。河流特点是河长短、流速快、多直流入海、属于季节性河流。研究区出露的地层较简单，从老至新依次为：元古界、中生界白垩系和新生界第四系。根据研究区含水介质类型和埋藏条件，将地下水划分为三大类：松散岩类孔隙水、碎屑岩类孔隙裂隙水和基岩裂隙水（图 2，图 3）。第四系孔隙含水岩组为主要含水层，主要分布在研究区河流中下游的河谷平原，该区域渗透能力较好，含水量随含水层厚度的增加而增大，其单井出水量在 100500m3d1之间。碎屑岩类孔隙裂隙水仅赋存于研究区东南的黄山青山村一带，含水层透水性、富水性均较弱。2研究方法2.1样品采集与测试+4HCO3SO24NO3本研究采用的数据均来自青岛市崂山区水文地质调查数据（20162017 年），所有水样检测分析均在山东省地矿工程勘察院实验室完成，测试内容主要包括 K+、Na+、Ca2+、Mg2+、NH、Cl、可溶性总固体（TDS）、总碱度（TA）、总硬度（TH）、pH 等 指 标。采 样 点 位 置 如 图 4。采 用Origin2017 软件对水化学参数进行描述性统计分析，绘制了 Piper 三线图、Chadhas 矩形图、Chadhas 矩形图，并运用主成分分析法讨论研究区地下水主要影响因素。借助 ACRGIS 软件空间分析功能对研究区浅层地下水海水入侵现状的空间分布规律进行分析，结合自然条件及人类活动背景资料，揭示研究区地下水化学特征及成因。2.2统计和水化学特征分析主成分分析法经常被用来确定地下水质量的主要影响因子2627。这种方法在许多研究中被证实是非常有效的2832。利用水化学分析方法研究地下水特征和演化有助于识别地下水的水文地球化学过程。Piper三线图是运用最广泛的水化学分析方法之一，但该方法只能表示地下水的化学组成，不能表示地下水化学演化。GimnezForcada 提出了水化学相演化图HFED，根据海水阳离子/阴离子的组成，可以有效地检查海水入侵的位置和淡水/海水的混合百分比3334。研究区范围0510 km图1研究区位置图Fig.1Locationofthestudyarea第42卷第5期刚什婷等：基于多变量统计分析和水化学特征的海水入侵特征研究以青岛市崂山区为例983HCO3SO24HCO3HFED 的计算过程涉及的阳离子有 Na+、K+、Ca2+、Mg2+，阴离子有、CO32、Cl、。海水入侵的水文地球化学过程复杂，HFED 在简化复杂过程时忽略了一些次要过程，因此 HFED 的绘制过程只考虑占比最大的 Na+、K+、Ca2+、Cl、SO24HCO3SO24HCO3。其中 K+、Mg2+不在图中展示，仅在计算中使用。另外和都可以代表淡水，在图中一个指标表示即可。本文选取作为淡水中的特征阴离子。HFED 中横纵坐标均用各阴阳离子在总阴、总阳离子的浓度百分比表示。HFED 把含水层发生海0124 km剖面线河流研究区边界地下水类型喷出岩类孔洞裂隙水块状岩类裂隙水层状岩类裂隙水松散岩类孔隙水碎屑岩类孔洞裂隙水图2青岛崂山区水文地质略图Fig.2HydrogeologicalsketchofLaoshandistrict,Qingdao984中国岩溶2023年水入侵时的水化学演化分为侵入期和恢复期，分别表示咸化过程和淡化过程。HCO3HCO3HCO3SO24一般来说，Na+和 Cl是海水/盐水的主要离子、而 Ca2+和通常是淡水的主要离子35。因此，沿海地下水中高浓度的 Na+和 Cl离子可能表明海水混合的显著影响和盐水的出现，而大量的 Ca2+和主要反映了水岩石相互作用的贡献36。通过绘制Chadhas 矩形图可以确定地下水盐化的来源。Chadhas 矩形图中阴阳离子的值都为相对质量浓度，阴阳离子的和分别为 100%。大矩形可以分为四个大区，分别代表补给水；反向离子交换水；海水入侵；离子交换水。每个大区分为 4 个次级的矩形，总计16 个次一级的矩形，每一个矩形都以c(Ca2+Mg2+)c(Na+)为 X 轴，c()c(+Cl)为 Y 轴。HCO3SO24海水入侵地下水质量指数37GQISWI由 Na+、K+、Ca2+、Mg2+，、Cl、七种离子浓度计算得到，该指数实现了把海水入侵过程中一系列复杂的水文化学过程量化。在 GQISWI图中可以把 Piper 三线图划分成 I 到 VI 区（图 5），并给各分区分别赋值，实现把 Piper 三线图结果量化。3结果与讨论3.1水化学指标统计结果分析对地下水样品检测结果的主要指标进行数理统计分析。表 1 中给出各水化学指标的最小值、最大值、平均值和标准偏差。表中各数据，尤其是 TDS、SO24HCO3TH、Na+、Cl等水质指标，变化较大，表明研究区内地下水化学组分有很大差异。地下水中阳离子以Na+为主，浓度均值关系为 Na+Ca2+Mg2+K+；阴离子浓度 Cl。pH 是一个可以反映水文地球化学平衡信息的重要因子，结果表明，研究区地下水 pH 均值处于 7.08.0 之间，变异系数为 0.04，说明该地区的地下水呈中性至弱碱性，在区域变化较小。表中各数据，尤其是 TDS、TH、Na+、Cl等水质指标，变化较大，表明研究区内地下水化学组分有很大差异。地下水中 Cl浓度变化幅度较大，且其均值超出了地下水质量标准（GB/T148482017）38的限值（250mgL1）和海水入侵的分界值（250mgL1），说明地下水可能发生一定程度的海水入侵，且基岩裂隙水中 Cl浓度远低于第四系孔隙水浓度，说明降雨使地下水 Cl浓度降低。TDS 质量浓度变化比较大，最大值为 17138.1mgL1，根据 TDS 含量分类39，地下水存在淡水、微咸水、咸水，部分地区地下水为盐水。SO24在进行主成分分析时，有 9 个变量被考虑在内：Na、Mg、Ca、K、Cl、SO4、HCO3、NO3、pH。由TotalVarianceExplained（主成分特征根和贡献率）可知，特征根 1=4.76，特征根 2=1.725，前两个主成分的累计方差贡献率达 72.06%，即涵盖了大部分信息。表明前两个主成分能够代表最初的 9 个指标来分析地下水水化学情况，故提取前两个主成分即可。主成分分别记作 PC1、PC2。正如预期的 Na+、K+、Ca2+、Mg2+、Cl、在 PC1 有较高的荷载（表 2，图 6），相75502502575海拔高程海拔高程mm75土寨河LS012QSQXKQLZK01养虾池入侵锋线黄海5025000.51.0 km25粉质黏土含砾混砾砂淤泥中粗粒碱长花岗岩黏土质砂细砂、中细砂1Tk35图3土寨河流域水文地质剖面图 A-AFig.3Hydrogeologicalprofile(A-A)oftheTuzhairiverbasin第42卷第5期刚什婷等：基于多变量统计分析和水化学特征的海水入侵特征研究以青岛市崂山区为例985HCO3NO3NO3NO3NO3关性强，反映的是水岩作用对水化学的影响，如蒸发、岩石与水的相互作用，以及可能与海水入侵有关。为证实这一假设，研究 Piper 三线图、HFED 图、Chadhas 矩形图等其他相关指标，pH、在第二主成分上有较高载荷，相关性强，说明发生了碳酸盐矿物的溶解；为负值（代表了农业活动和碱性矿物肥料的主要特征），反映了农业活动对地下水的影响，也可能与降水直接补给有关。由于农业和人类活动对青岛市地下水质量的影响前人已有大量研究40,25，本文只关注海水入侵对地下水的影响，因此对不做过多研究。3.2水化学特征分析3.2.1Piper 三线图解析水质演化Piper 三线图能反映地下水主要离子组成变化和HCO3HCO3SO24水化学特征，可以用来划分水化学类型。对于阳离子（图 7 左下三角形），大部分数据中主导阳离子是Ca2+，平均毫克当量百分比为 75%，碱土金属元素（Ca2+Mg2+）浓度大于碱金属元素（K+Na+），仅在少数水样中碱金属元素（K+Na+）大于碱土金属元素（Ca2+Mg2+）。对于阴离子（图 7 右下三角形），大部分数据主导阴离子是，平均毫克当量百分比为75%，弱酸根离子（）浓度大于强酸根离子（Cl+）。研究区地层中常见的矿物有方解石（CaCO3）、白云石（CaMg(CO3)2）等碳酸盐矿物，石膏（CaSO42H2O）等硫酸盐矿物，由此初步判断地下水化学组分受碳酸盐矿物和硫酸盐矿物溶滤作用较强。按照 GQISWI中对 Piper 图中水化学类型的划分，崂山区地下水多为 ClSO4-Na 型和 SO4Cl-CaMg 型地下水，部分地区为 HCO3-Ca 型，还有少数监测点地下水为 Cl-Na 型水组成。HCO3-Ca 型是淡水的主要泉水水库河流024 km 研究区范围机民井取样点图4取样监测点位置图Fig.4Locationsofsamplingandmonitoringsites986中国岩溶2023年特征，Cl-Na 型是海水的主要特征，Cl-CaMg 和HCO3-Na型是咸淡水之间过渡的水化学类型。图 7 表明取样点位置只有少部分地区地下水是 HCO3-Ca 型淡水，多数地区地下水中 Cl含量较大，并且有部分区域的地下水化学类型为 Cl-Na 型，上述研究结果表明，部分地区已经发生了严重的海水入侵现象，且多数地区地下水化学类型处于咸淡水之间的过渡阶段。HCO3地下水化学类型分布较一致，河流中上游阳离子主要以 Ca2+和 Mg2+为主，阴离子主要以为主，由陆地向海洋，主要阳离子逐渐向 Na+变化，阴离子向 Cl变化。图中取样点水化学类型由地势高的地方向海洋趋近，水化学类型由 HCO3-Ca 型逐渐向 Cl-Na 型演化，可能沿径流路径发生了阳离子交换，以及海水不同程度的混合作用。如在石人河，LS015、LS011、LS009、LS010 从上游到海岸边依次排列，水化学类型分别由 HCO3Cl-CaMgSO4Cl-CaMgCl-CaMgClSO4-Na 变化；在土寨河，LS012、LS005、LS001、LS003、LS010 从上游到海岸边依次排列，水化学类型分别由 SO4-CaMgSO4Cl-NaCaMgSO4Cl-CaMgSO4Cl-NaMgSO4Cl-Na 变 化；在王哥庄河 LS019、LS020、LS016、LS033、LS035 从上游到海岸水化学类型由 HCO3Cl-CaMgHCO3Cl-NaCaMgHCO3Cl-NaCaMgCl-CaMgCl-NaMg 变 化；在 南 部 的 大 河 东小 河 东 LS085、LS081、LS082、LS092、LS089、LS090 从上游到海岸边依次排列，水化学类型分别由 HCO3Cl-NaCaMgHCO3-NaCaCl-NaCaSO4Cl-NaCaMgSO4Cl-NaCaMgCl-NaCaMg 变化。根据各取样点离子浓度关系，绘制了地下水七大离子浓度占比图（图 8）。由图 8 可知，地下水中 Cl含量占主导地位，阳离子以 Ca2+、Na+为主。3.2.2HFED 图解译水文地球化学过程崂山区地下水 HFED 如图 9 所示，显示了研表 1地下水水化学参数统计特征值（单位：mgL1，pH 除外）Table1Statisticsofhydrochemicalparametersofgroundwater(unit:mgL1,exceptforpH)分区项目pHTDSTHCa2+Mg2+K+Na+HCO3SO24ClNO3基岩裂隙水Mean7.34998.81380.6094.2334.711.78199.04104.33135.62385.9435.94SD0.272032.16670.59200.6846.612.49512.0182.42218.811103.4532.57Cv0.042.031.762.131.341.402.570.791.612.860.91Min6.9067.2146.1212.933.360.132.177.653.0521.530Max8.208877.643064.85910.97191.908.822125.00306.08778.464812.71101.35第四系孔隙水Mean7.201435.83554.2688.3581.0412.92277.95108.77235.23560.8369.69SD0.383211.69802.7263.12163.2736.36927.5675.61403.371722.6367.98Cv0.052.241.450.712.012.813.340.701.713.070.98Min6.5074.8339.208.314.480.176.6715.3010.6813.130.36Max8.5017138.103597.12263.18728.01200.005000.00369.841617.989160.70244.75地表水Mean7.30301.99138.7140.888.904.4330.6985.0756.1140.9123.25SD0.32302.26127.0635.809.185.7541.8290.2662.8153.3833.27Cv0.041.000.920.881.031.301.361.061.121.301.43Min6.8091.6948.5415.432.310.164.4327.9819.217.185.52Max7.60814.84351.78100.8024.3112.59103.60243.46165.70132.7882.60注：Min为最小值，Max为最大值，Mean为平均值，SD为标准差，Cv为变异系数。Note:Minrepresentsminimumvalue;Maxrepresentsmaximumvalue;Meanrepresentsaveragevalue,SDrepresentsstandarddeviation;Cvrepresentsvariationcoefficient.1007550VCa-ClCa2+Mg2+Ca-Mg-ClSO42+ClCa-HCO3Na+K+Ca-Na-HCO3HCO3Na-HCO3Na-ClGQlpiper/mixGQlpiper/dom2525507510000图5GQISWI对 Piper 图的分区Fig.5DomainsofGQISWinpiperdiagram第42卷第5期刚什婷等：基于多变量统计分析和水化学特征的海水入侵特征研究以青岛市崂山区为例987究区地下水部分采样点为 Cl-Na 型海水，表明河流中下游沿海地区地下水发生了海水入侵。部分采样点位于混合线以下，表明该区域地下水正处于海水侵入阶段，并且多数地区地下水已经完成由HCO3-Ca 型向 Cl-Ca 型的演化。HFED 图隐藏了坐标轴 033.3%段各离子浓度相对均衡的区域，用“Mix”表示图中主导离子浓度占比小于 50%的水化学类型33。3.2.3Chadhas 矩形图解析水质演化根据 Chadhas 图（图 10），受海水入侵影响的钻孔（LS010、LS029、LS082 号样品）位于海岸线上，分别对应江家土寨东浦里社区北入侵段、王山口南表 2主成分分析法组成矩阵Table2MatrixformedbyprincipalcomponentanalysisPC1PC2pH0.007700.64180Ca2+0.270270.10640Mg2+0.431210.05697K+0.396810.00419Na+0.444050.00605HCO30.099170.48226SO240.414960.07280Cl0.447430.01636NO30.063990.579080.80.60.40.2700.4310.3970.4440.4150.4820.4470.5790.642PC1PC20.200.20.40.6CaMgKNa SO4HCO3ClNO3pH图6主成分分析结果图：PC1、PC2Fig.6Resultsofprincipalcomponentanalysis:PC1andPC21008060Ca+Mg40200100806040SO4+Cl200100806040Mg20010080604020002040Na+K60SO4801000204060HCO3801001008060CaCl40200020406080100孔隙水裂隙水地表水.Ca-HCO3.Na-Cl.Ca-Na-HCO3.Ca-Mg-Cl.Ca-Cl.Na-HCO3图7地下水水化学 Piper 图Fig.7Piperdiagramofgroundwaterhydrochemsitry988中国岩溶2023年SO24桑园社区北王哥庄港东入侵段、登瀛村南栲栳岛入侵段。区域、表示 c(Ca2+Mg2+)c(Na+)，其中c(Na+)c(Na+)和 c()c(+Cl)，并且 c(Na+)和 c()相对质量浓度均小于 25%，这种水有永久硬度，用于灌溉时，不会形成 NaCO3盐沉积，地下水化学类型为 SO4Cl-CaMg 型水（采样点 LS001、LS089、LS111、LS033）。3.3海水入侵现状评价SO24HCO3采用基于熵权的集对分析法和 GIS软件结合，选取了 Cl、矿化度、Cl/、SAR 五项化学特征指标作为评价因子，评价因子等级划分指标见表 3，对崂山区 53 个监测井点海水入侵程度分别进行了评价，并绘制了海水入侵现状图（图 11）。由图 11 可以看出，青岛市崂山区地下水综合质量较好，受海水入侵影响较少，海水入侵主要分布于沿海冲洪积平原、河流下游及河道两侧区域，尤其是河流入海口。自沿海地区向内陆地区海水入侵现状逐渐减缓，海水入侵主要分布于江家土寨东浦里社区北入侵段，王哥庄港西港东入侵段、仰口湾入侵段、登瀛村栲栳岛入侵段。4结 论本文以 20162017 年地下水水质调查数据为背景，以地下水中主要离子浓度为水化学指标，对青岛市崂山区地下水化学特征及海水入侵状况进行分析，得出以下结论：（1）水化学指标统计分析结果显示 TDS、TH、Na+、Cl等水质指标变化较大，表明研究区内地下水化学组分有很大差异。地下水中 Cl浓度变化幅度较大，且其均值超出了地下水质量标准（GB/T148482017）的 限 值（250mgL1）和 海 水 入 侵 分 界 值020406080100采样点编号LS019LSl10LSl11LSl04LSl02LS097LS096LS092LS090LS089LS085LS084LS083LS082LS081LS077LS076LS075LS074LS073LS071LS063LS060LS051LS048LS046LS043LS041LS040LS035LS033LS030LS029LS027LS026LS024LS023LS022LS020LS016LS015LS012LS011LS010LS009LS005LS003LS001ClSO42HCO3Na+K+Mg2+Ca2+图8采样点地下水七大离子浓度占比图Fig.8Percentagesofconcentrationsofthesevenionsingroundwaterofthesamplingsites第42卷第5期刚什婷等：基于多变量统计分析和水化学特征的海水入侵特征研究以青岛市崂山区为例989表 3海水入侵指标的等级划分（单位/mgL1）Table3Classificationofindexesofseawaterintrusion(unit/mgL1)特征因子IIIIIIIV无入侵轻度入侵中度入侵严重入侵Cl25060015001500SO2420045012001200M1000200030003000SAR23.551010Cl/HCO30.51.06.66.61005050Cl%33.3%1005033.3%50100505033.3%100HCO3NaMix HCO3NaHCO3MixNaMixHCO3MixNaHCO3MixCaMixClMixNaMixHCO3MixCaMixHCO3CaMixClCaClCaHCO3CaMixClMixCaClMixCaClMixNaMixClNaClNa5050Na%+K%Ca%33.3%1123456789101112131415162678123459131014151611淡化过程咸化过程HCO3%+CO3%图9崂山区地下水 HFEDFig.9EvolutionofhydrochemicalfaciesofgroundwaterinLaoshandistrict7550251001075505075100252525507510000981576541431613211121c(Ca2+Mg2+)c(Na+)c(HCO3)c(SO42+Cl)图10Chadha 矩形水化学类型图Fig.10Chadharectanglehydrochemicaldiagram990中国岩溶2023年（250mgL1），说明地下水可能发生一定程度的海水入侵。根据主成分分析可知控制研究区地下水水质的因素主要有大气降水和人为活动。（2）Piper 三线图、HFED、Chadhas 矩形图结果显示，取样点位置只有少部分地区地下水是 HCO3-Ca 型淡水，多数地区地下水中 Cl含量较大，并且有部分区域的地下水化学特征是 Cl-Na 型，表明部分地区已经发生了严重的海水入侵现象，多数地区地下水化学类型处于咸淡水之间的过渡阶段。研究区地下水海水入侵有向内陆扩展的趋势。（3）根据海水入侵评价结果显示，海水入侵水主要分布于江家土寨东浦里社区北入侵段，王哥庄港西港东入侵段、仰口湾入侵段、登瀛村栲栳岛入侵段。参考文献AndersenMS,NyvangV,JakobsenR,PostmaD.Geochemicalprocessesandsolutetransportattheseawater/freshwaterinter-1faceofasandyaquiferJ.GeochimicaetCosmochimicaActa,2005,69:39793994.WernerAD,SimmonsCT.Impactofsea-levelriseonseawaterintrusionincoastalaquiferJ.Groundwater,2009,47：197-204.2LiuY,JiaoJJ,LiangW,KuangX.Hydrogeochemicalcharac-teristicsincoastalgroundwatermixingzoneJ.AppliedGeo-chemistry,2017,85：49-60.3KazakisN,PavlouA,VargemezisG,VoudourisKS,SouliosG,PliakasF,TsokasG.Seawaterintrusionmappingusingelectri-calresistivitytomographyandhydrochemicaldata.Anapplica-tioninthecoastalareaofeasternThermaikosGulf,GreeceJ.ScienceoftheTotalEnvironment,2016,543：373-387.4Alfarrah N,Walraevens K.Groundwater overexploitation andseawater intrusion in coastal areas of arid and semi-aridregionsJ.Water,2018,10（2）：143.5BovoloCI,ParkinG,SophocleousM.Groundwaterresources,climate and vulnerabilityJ.Environmental Research Letters,2009,4（3）：35001.6郭明华,王敬,戴长国,蒋雷,张浩清,柳超.文登区浅层地下水化学演化与海水入侵研究J.海洋科学,2021,45（7）：57-65.GUO Minghua,WANG Jing,DAI Changguo,JIANG Lei,7024 km河流水库研究区范围海水入侵程度值严重入侵无入侵图11青岛市崂山区海水入侵空间分布Fig.11SpatialdistributionofseawaterintrusioninLaoshandistrict,QingdaoCity第42卷第5期刚什婷等：基于多变量统计分析和水化学特征的海水入侵特征研究以青岛市崂山区为例991ZHANGHaoqing,LIUChao.Studyonthehydrochemicalevo-lutionofgroundwaterandseawaterintrusionintheshallowlayerofWendengdistrictJ.MarineScience,2021,45（7）：57-65.赵倩.锦州小凌河扇地海水入侵过程中水文地球化学作用及影响因素研究D.廊坊：防灾科技学院,2021.ZHAOQian.Hydrogeochemistryandinfluencingfactorsintheprocess of seawater intrusion in Xiaoling river fan,JinzhouCityD.Langfang:InstituteofDisasterPrevention,2021.8MondalNC,SinghVP,SinghVS,SaxenaVK.Determiningtheinteraction between groundwater and saline water throughgroundwater major ions chemistryJ.Journal of Hydrology,2010,388,100-111.9王昊,王中良.天津北大港水库水质咸化的地球化学机理分析J.天津师范大学学报(自然科学版),2016,36（6）：29-34.WANG Hao,WANG Zhongliang.Geochemical analysis onwater salinization mechanism of the Tianjin Beidagang reser-voirJ.JournalofTianjinNormalUniversity(NaturalScienceEdition),2016,36（6）：29-34.10辛祥.大沽河下游地区地下水氯离子时空变化特征与调控的数值模拟D.青岛：青岛大学,2019.XINXiang.Numericalsimulationofspatialandtemporalchar-acteristics and control of chloride ion in groundwater in thedownstreamareaofDaguriverD.Qingdao:QingdaoUniver-sity,2019.11林曼曼.青岛近海海域灾害地质特征研究D.石家庄：石家庄经济学院,2014.LINManman.Thestudyonthecharacteristicsofgeo-hazardfac-torsinQingdaooffshoreD.Shijiazhuang:UniversityofEco-nomics,2014.12陆求裕.长乐滨海地区海水入侵现状及成因分析J.地质灾害与环境保护,2020,31（4）：101-105.LUQiuyu.ThepresentsituationoftheseawaterintrusionalongthecoastalareaofChangleJ.JournalofGeologicalHazardsandEnvironmentPreservation,2020,31（4）：101-105.13梁彬,关碧珠,崔光中,揭江,唐祚旺.广东湛江开发巨型盆地防治海水入侵研究J.中国岩溶,2001,20（1）：77.14秦凯凯.管状海水入侵形成机理及影响因素研究D.大连：大连海事大学,2017.QINKaikai.Studyonformationmechanismandinfluencingfac-torsoftubularseawaterintrusionD.Dalian:DalianMaritimeUniversity,2017.15崔震,陈广泉,徐兴永,王东亮,喻龙.北长山岛海水入侵成因机理及现状评价J.海洋环境科学,2015,34（6）：930-936.CUI Zhen,CHEN Guangquan,XU Xingyong,WANGDongliang,YULong.MechanismandassessmentofseawaterintrusioninthenorthernChangshanislandJ.MarineEnviron-mentalScience,2015,34（6）：930-936.16徐铭霜.某沿海地区海水入侵动态变化规律及防治措施研究D.青岛：青岛理工大学,2022.XUMingshuang.Researchonthedynamicvariationruleofsea-waterintrusionandcontrollingmeasuresinacoastalareaD.17Qingdao:QingdaoUniversityofTechnology,2022.王蜜蕾,窦衍光,邹亮,薛碧颖,胡睿,岳保静,徐刚,林曦,李林森.青岛崂山周边地下水化学特征与矿泉水成因分析J.海洋地质前沿,2021,37（9）：17-24.WANGMilei,DOUYanguang,ZOULiang,XUEBiying,HURui,YUEBaojing,XUGang,LINXi,LILinsen.Hydrochemi-calcharacteristicsofgroundwaterandgenesisofmineralwateratLaoshanmountainandsurroundingareas,QingdaoJ.MarineGeologyFrontiers,2021,37（9）：17-24.18胡慧君.青岛市大沽河典型区非饱和饱和带地下水耦合模型D.南京：南京大学,2019.HU Huijun.An integrated modeling of variably saturatedgroundwater flow system in a typical area of Dagu riverbasinD.Nanjing:NanjingUniversity,2019.19刘久潭,周丹,高宗军,王敏,马媛媛,张洪英,时孟杰,董杰.青岛西海岸新区地下水水化学特征及水质评价J.山东科技大学学报(自然科学版),2019,38（2）：14-24,43.LIU Jiutan,ZHOU Dan,GAO Zongjun,WANG Min,MAYuanyuan,ZHANGHongying,SHIMengjie,DO