岩溶水是人类重要的供水水源, 在我国北方分布面积广且资源丰富[1, 2], 在保障城乡居民生活和支持经济社会可持续发展等方面发挥着十分重要的作用[3, 4]. 济南趵突泉泉域作为北方典型岩溶泉域, 分为间接补给区、直接补给区和汇集出露区这3个区域, 岩溶水是泉城济南经济社会发展不可缺少的重要水源[5, 6]. 但是随着社会经济的快速发展和人口数量的增加, 济南市岩溶地下水的开发利用与城市的发展产生了巨大的矛盾, 20世纪70年代以来济南岩溶地下水开采量猛增, 城区岩溶水位不断下降导致泉水断流和泉城失色[7 ~ 9], 同时岩溶水中Cl-、NO3-、SO42-和总硬度等指标持续上升, 水体污染问题日益突出[10, 11]. 探明济南岩溶地区地下水水化学特征及形成机制, 对岩溶水资源可持续利用及保泉工作具有重要的指导意义.
近年来, 已有诸多学者着重于岩溶发育规律、岩溶水运移规律以及岩溶水水质形成与演化的研究, 如孙斌等[12, 13]、邢立亭等[14]和李严等[15]采用泉流量相关分析、示踪试验、动态水文时间序列和水化学统计等方法对趵突泉泉域水文地质条件、地下水补给来源、水循环特征及水环境等进行大量研究；李常锁等[16]、李传生等[17]和Wu等[18]对济南岩溶地下水的岩溶发育条件、地下水水动力特征和水化学特征等进行了探讨；高宗军等[19]和殷秀兰等[20]对济南地区泉水、岩溶地下水中主要阴阳离子、微量元素及稳定环境同位素进行分析, 研究了济南地区地下水补给来源和比例等问题；孟祥玲等[21]和管清花等[22]对济南泉域岩溶水水质变化规律及其影响因素进行了详细分析. 前人关于趵突泉泉域地下水化学特征的研究多采用结合水文地质条件、水文地球化学理论及人类活动影响的方法, 主要集中在汇集出露区(排泄区)、径流区等局部地下水系统的水化学成因分析. 但是在趵突泉泉域尺度上, 对包含间接补给区、直接补给区和汇集出露区这3个功能区的岩溶水系统水化学特征空间差异性分析、控制因素及形成机制的研究较少. 此外, 单一的传统水文地球化学的研究方法如水化学图示、离子比值及数理统计分析等难以对复杂多变的水环境系统进行评价与分类研究[23], 为了研究结果的可靠性需要综合多种方法分析, 这也是当前水文地球化学研究的趋势[24].
因此, 本文在分析趵突泉泉域水文地质条件的基础上, 着眼于趵突泉泉域间接补给区、直接补给区和汇集出露区3个功能区且沿着补径排路径进行了样品的采集, 综合运用水化学参数统计、多元统计、自组织神经网络(SOM)分析、水化学图示和离子比例系数等方法, 系统研究了趵突泉泉域不同功能区地下水水化学特征空间差异性及其形成机制, 解释了不同区域水化学成分的控制因素, 以期为济南岩溶地下水资源可持续利用和水污染防治等提供理论依据.
1 研究区概况
济南趵突泉泉域位于山东省中部、泰山山脉的西北麓, 北临黄河与鲁西北平原相接, 整体地势南高北低. 研究区属暖温带大陆性半湿润季风气候区, 多年(1971~2015年)平均气温为12.8℃, 多年平均降雨量约为646.55 mm, 每年6~9月降水较多, 12月至次年5月降水量较少, 降水量在空间上自东南向西北递减.
研究区属华北地层大区、鲁西地层分区, 区内地层发育较为齐全, 分布有基底岩系和盖层, 属二元结构. 地层自南向北、从老至新依次出露太古代泰山群、古生代寒武系、奥陶系和新生代第四系(图 1).
图 1
Fig. 1
图 1 泉域地质及取样点分布示意
Fig. 1 Geological sketch map of spring area and distribution of sampling sites
结合研究区岩溶水的补给、径流和排泄特征, 可将泉域划分为间接补给区、直接补给区和汇集出露区这3个功能区[22, 25], 以区分泉水与其母体岩溶地下水形成过程中, 不同地段所发挥的不同作用与功能(图 2). 其中间接补给区地层岩性为寒武纪炒米店组灰岩以下至变质岩, 位于岩溶水子系统上游, 与市区泉水所属不同的含水岩组, 岩溶地下水径流通过断裂构造进入和补给泉水岩溶含水系统；直接补给区地层岩性主要为寒武系九龙群炒米店组上段、三山子组及奥陶系马家沟群灰岩、白云岩, 岩溶地下水能够直接接受大气降水及地表水渗漏补给, 且与市区泉水为同一含水岩组, 有直接水力联系和具有统一流场和水位, 对出露区地下水产生直接的水源补给；汇集出露区指整个岩溶泉域水系统下游地下水汇聚、储存和排泄的地区, 为岩溶水的承压排泄区.
图 2
Fig. 2
图 2 泉域水文地质剖面
Fig. 2 Hydrogeological profile of spring area
2 样品采集与分析
本次样品采集点结合研究区的地质、水文地质条件, 在地下水补给、径流和排泄区域进行了布设. 共采集地下水全分析样品81组, 取样点分布如图 1所示, 其中间接补给区取样点42个, 直接补给区取样点18个, 汇集出露区样品取样点21个, 采集时间为2022年6月初(枯水期). 各样品的采取、保存和送检均按《水质采样、样品的保存和管理技术规定》(HJ 493-2009)相关规定执行. 所有样品送至山东省地矿工程勘察院实验室进行测试, 水样在测试前都要通过0.45 μm的过滤器进行过滤. pH、电导率等易变参数采用Clean多参数测试计现场测试；主要阴离子(Cl-、SO42-、NO3-和SiO2)测试仪器为Thermo Fisher的ICS-600离子色谱仪；主要阳离子(Mg2+、Ca2+、Na+和Sr2+)测试仪器为PerkinElmer的电感耦合等离子体发射光谱仪(ICP-OES)；HCO3-采用盐酸滴定法测定, 相关检测标准依据HJ 776-2015水质32种元素的测定. 所有样品的阴阳离子平衡误差均在±5%以内.
3 结果与讨论
3.1 地下水化学特征空间差异
3.1.1 水化学参数统计分析
基于四分位法, 选取Na+、Ca2+、Mg2+、K+、Cl-、SO42-、HCO3-、NO3-、Sr2+、SiO2、TDS和pH共12个水化学参数绘制了箱线图(图 3). 从图 3中可以看出, 3个功能区地下水样品都表现出典型的碳酸盐矿物溶解呈碱性, pH值介于7.58~8.0之间；直接补给区水样具有相对较高的TDS、Ca2+、Na+、Mg2+、Sr2+、SO42-、HCO3-和NO3-, 可能与该区域地下水径流路径较长有关；间接补给区水样各项参数值均相对比较低, 表明该区域地下水径流时间较短, 水岩作用不够充分；汇集出露区各参数的取值范围相对较大, 可能与该区处于地下水径流的排泄区, 同时接受多水源的补给, 以及城区和郊区人为因素的影响存在差异有关. 此外, 3个功能区的ρ(NO3-)的均值均大于20 mg·L-1, 表明人为活动对趵突泉泉域地下水水化学成分的改变比较明显.
图 3
Fig. 3
图 3 泉域不同功能区地下水化学参数箱线图
Fig. 3 Box-line diagram of chemical parameters of groundwater in different functional areas of spring area
3.1.2 地下水化学类型
从趵突泉泉域不同功能区岩溶地下水水化学Piper图看出(图 4), 研究区几乎所有的地下水样点均落在1区, 表明碱土金属离子(Ca2+和Mg2+)大于碱金属离子(Na+和K+), 大部分样点落在5区, 表明碳酸盐毫克当量百分比大于50%, 即各功能区阴离子以HCO3-和SO42-为主, 阳离子以Ca2+为主, Na++K+和Cl-含量不足总离子含量的30%；地下水水化学类型以HCO3-Ca和HCO3·SO4-Ca型为主, 从间接补给区到汇集出露区部分地下水样点水化学类型由HCO3-Ca型向HCO3·SO4-Ca型转变. 直接补给区地下水中的SO42-浓度高于间接补给区. 汇集出露区水样点大体分布于直接补给区和间接补给区之间, 结合不同功能区的TDS可知, 汇集出露区岩溶地下水可能同时受到直接补给区和间接补给区的补给, 此外, 该区Cl-浓度相对较高, 如落在7区的水样点, 水化学类型为Cl·HCO3-Na型, 该点位于城区基岩裸露区, 受人类活动影响较大.
图 4
Fig. 4
图 4 不同功能区地下水化学Piper图
Fig. 4 Piper diagram of groundwater chemistry in different functional areas
3.2 地下水化学组分来源分析
3.2.1 相关分析
通过地下水化学各参数的相关分析可以判断离子之间的来源关系[26, 27], 由趵突泉泉域不同功能区地下水化学参数Pearson相关性分析发现(表 1), TDS与Ca2+、Mg2+、Cl-、SO42-、HCO3-和NO3-的相关系数均在0.6以上, 说明以上离子对TDS贡献较大, 其中和Ca2+的相关关系最为显著, 相关系数达到了0.972；Cl-与Na+、Ca2+、Mg2+和TDS存在明显的相关关系, 相关系数均在0.6以上, 表明Cl-与Na+、Ca2+和Mg2+可能有同一来源, 即来源于岩盐溶解或蒸发岩的风化；HCO3-与Ca2+和Mg2+的相关系数达到了0.6左右, 表明HCO3-、Ca2+和Mg2+来源于碳酸盐岩(方解石和白云石)的溶解；此外, SO42-与Ca2+和TDS的相关系数分别为0.73和0.741存在明显的相关关系, 表明含水层中可能发生了石膏的溶解.
表 1
(Table 1)
表 1 趵突泉泉域地下水化学参数Pearson相关性1)
Table 1 Pearson correlation of groundwater chemical parameters of the Baotu Spring area
Na+
Ca2+
Mg2+
K+
Cl-
SO42-
HCO3-
NO3-
Sr2+
SiO2
TDS
pH
Na+
1
Ca2+
0.444**
1
Mg2+
0.353**
0.707**
1
K+
0.121
-0.288**
-0.159
1
Cl-
0.721**
0.694**
0.596**
-0.0317
1
SO42-
0.569**
0.730**
0.453**
-0.147
0.490**
1
HCO3-
0.308**
0.696**
0.583**
-0.376**
0.342**
0.469**
1
NO3-
0.253*
0.755**
0.644**
-0.053
0.537**
0.360**
0.236*
1
Sr2+
0.023
0.083
0.023
-0.159
-0.054
0.176
0.094
0.008
1
SiO2
0.259*
0.550**
0.339**
-0.199
0.503**
0.293**
0.309**
0.424**
0.074
1
TDS
0.568**
0.972**
0.764**
-0.181
0.761**
0.741**
0.616**
0.817**
0.071
0.532**
1
pH
-0.371**
-0.594**
-0.370**
-0.002
-0.515**
-0.335**
-0.411**
-0.499**
-0.117
-0.582**
-0.602**
1
1)*表示在0.05水平上显著相关；**表示在0.01水平上显著相关
表 1 趵突泉泉域地下水化学参数Pearson相关性1)
Table 1 Pearson correlation of groundwater chemical parameters of the Baotu Spring area
3.2.2 自组织神经网络(SOM)分析
自组织神经网络(简称SOM)作为一类无监督学习的神经网络模型, 提供了一种将高维数据投射到低维空间(如二维空间)的可视化方法, 可以保留原始数据的拓扑结构和度量关系, 这使其成为非线性降维和无监督分类的理想选择, 适用于对复杂多变的水环境系统进行评价与分类研究[28, 29]. 本文采用MATLAB软件中的SOM工具箱建立算法程序, 将研究区81组地下水样本的9种水化学指标的标准化数据(均值为0, 方差为1)通过SOM训练后, 分别得到各指标对应的64个六边形(神经元)组成的平面映射图[图 5(a)], 每个小六边形表示1个神经元, 每个平面映射上的神经元的颜色表示其对应参数值的大小, 其中蓝色表示该属性值较小, 红色则相反, 通过颜色深浅的变化可以直观地反映输入的水化学指标之间的关系. 由图 5(a)可以看出, 研究区地下水中9种水化学指标分布差异明显. 其中HCO3-、Ca2+、Mg2+和SO42-浓度与TDS变化规律相似, 表现为从上到下逐渐增加, 表明以上指标受水岩相互作用(如矿物溶解)的影响；Na+和Cl-与NO3-浓度呈现对角线梯度, 其中Na+与Cl-从右上角向左下角递增, NO3-从左上角向右下角递增, 表明以上离子具有相似的来源, 如人为来源等；K+的浓度与其他离子的相似度较低.
图 5
Fig. 5
图 5 泉域地下水9个变量的SOM映射图与聚类结果空间分布
Fig. 5 SOM mapping of nine variables of spring area groundwater and the spatial distribution of clustering results
由于自组织神经网络不能自动确定最佳聚类数, 因此可将SOM的结果作为K-means聚类分析的输入值, 对映射对应的样本进行聚类, 聚类结果如图 5(b)所示, 可以看出, 研究区地下水样品可分为4类, 且各聚类结果具有较好的空间差异性[图 5(c)]. 第1类地下水样品主要分布在间接补给区及汇集出露区与直接补给区交界处, 该类地下水样品具有较高浓度的K+, 浓度平均值为2.15 mg·L-1；第2类地下水样品主要分布在间接补给区与汇集出露区的西北角, 该类地下水样品除HCO3-浓度较高之外, 其余离子浓度均较低；第3类地下水样品主要分布在直接补给区与汇集出露区的东北角, 具有较高的Na+浓度, 其余各项水化学指标为中等水平；第4类地下水样品在3个区域均有分布, 空间分布规律性不强, TDS值相对较高, 且具有较高浓度的Ca2+、Mg2+、Cl-、SO42-、HCO3-和NO3-, 其中ρ(NO3-)范围为24.3~112 mg·L-1, 平均值为41.1 mg·L-1, 表明该类地下水受人类活动影响较大[23].
3.2.3 因子分析
因子分析是基于降维的思想, 在不损失或者少损失原始数据信息的情况下, 将多变量聚合成少数几个独立的公共因子, 以提高数据分析的效率与准确性[30]. 因此, 为进一步分析地下水各指标的内在联系, 探讨地下水中离子组成及控制因素[31], 选取了10项水化学指标(K+、Na+、Ca2+、Mg2+、Cl-、SO42-、HCO3-、NO3-、可溶性SiO2和TDS)进行因子分析[32], 所有数据均通过了KMO和巴特利特检验, 显著性水平接近于0. 基于特征值大于1的筛选要求, 采用最大方差法对成分矩阵进行旋转, 分别识别出影响水化学指标的两个主成分因子(F1和F2), 累计贡献率68.181%, 各因子提取结果与旋转成分矩阵见表 2.
表 2
(Table 2)
表 2 趵突泉泉域地下水样品旋转成分矩阵1)
Table 2 Rotational composition matrix of water samples of the Baotu Spring area
参数
F1
F2
K+
0.026
-0.863
Na+
0.751
-0.323
Ca2+
0.873
0.420
Mg2+
0.733
0.315
Cl-
0.872
-0.084
SO42-
0.727
0.152
HCO3-
0.507
0.611
NO3-
0.724
0.156
SiO2
0.539
0.266
TDS
0.947
0.270
特征值
5.495
1.323
贡献率/%
54.947
13.285
累计贡献率/%
54.947
68.181
1)黑体字表示高荷载
表 2 趵突泉泉域地下水样品旋转成分矩阵1)
Table 2 Rotational composition matrix of water samples of the Baotu Spring area
第一主成分因子F1的贡献率为54.947%(表 2), 与Na+、Ca2+、Mg2+、Cl-、SO42-、NO3-和TDS显著正相关, 其中NO3-为人类活动的特征因子, 其他离子与盐岩、蒸发岩等矿物有关, 结合相关分析(表 1), 推测主成分因子F1受盐岩、蒸发岩等矿物的风化和人类活动影响；第二主成分因子F2的贡献率为13.285%(表 2), 与HCO3-正相关, 与K+、Na+和Cl-负相关, 由于HCO3-与碳酸盐岩矿物有关, 推测主成分因子F2可能受碳酸盐岩矿物风化的影响. 据主成分因子的荷载矩阵及对应变量的相关系数绘制趵突泉泉域地下水样品离子的主成分分析(图 6), 由图 6(a)可知, Ca2+、Mg2+、SO42-、HCO3-、NO3-和可溶性SiO2离子浓度主要分布在第Ⅰ象限, 说明以上离子浓度同时受主成分因子F1和F2的影响；K+、Na+和Cl-离子浓度分布在第Ⅳ象限, 说明以上离子浓度主要受主成分因子F1的影响. 由图 6(b)可知, 3个区域地下水样品的离子浓度分布规律性不强, 差异性不明显, 并且有少数地下水样品分布在第Ⅲ象限, 说明除主成分因子F1和F2之外, 有其他因素影响离子浓度.
图 6
Fig. 6
图 6 趵突泉泉域地下水样品离子的主成分分析
Fig. 6 Principal component analysis of ions of groundwater samples of the Baotu Spring area
3.2.4 水化学离子起源
对于水中离子起源的自然影响因素, Gibbs根据世界河流、湖泊及主要海洋水TDS与离子(Na+、Ca2+、Cl-和HCO3-)的关系确定了天然水化学成分的主要来源[33, 34], 因此Gibbs图可用于识别研究区地下水中各离子的化学组成及其主要来源, 从而判断地下水样品所经受的不同作用, 即大气降水补给、岩石矿物的风化水解和矿物蒸发-结晶等. 从趵突泉泉域不同功能区地下水化学Gibbs图可以发现(图 7), 3个功能区地下水取样点基本位于岩石风化作用区域内, 且大部分点落在图中部左侧, 表明各功能区含水层中溶解盐岩类矿物受到风化作用是地下水中离子成分组成的主导因素, 即地下水中离子成分来源于成岩矿物的风化水解, 地下水径流过程中水动力条件较好, 溶滤作用比较充分, 受降水作用和蒸发-结晶作用的影响不明显. 此外, 有少部分汇集出露区和间接补给区取样点位于Gibbs图右侧中部, ρ(TDS)中等(约300 mg·L-1), Na+/(Na++Ca2+)的值约为0.8, Cl-/(Cl-+HCO3-)的值接近于0.5, 表明可能存在阳离子交换作用或这些取样点周围的地下水可能存在外源输入的Na+和Cl-.
图 7
Fig. 7
图 7 不同功能区地下水化学Gibbs图
Fig. 7 Gibbs map of groundwater chemistry in different functional areas
3.2.5 水化学离子比
地下水化学组分、各离子的比值关系与地下水环境相关, 不同的水环境下发生不同的水岩相互作用, 因此可以通过地下水中阴、阳离子的种类和比例关系, 来揭示地下水中溶质的来源成因及其演化过程[35 ~ 37]. 为了识别不同成因或不同条件下形成的地下水, 分析不同功能区离子比值的差异, 选取Ca2+、Mg2+、HCO3-、SO42-、Na+和Cl-等离子之间的比例系数来探讨各区域地下水水化学组分的来源、地下水水化学形成特征及其差异性.
地下水中Ca2+、Mg2+、HCO3-和SO42-如果全部来源于含水层中矿物的溶解, 如方解石、白云石(碳酸盐矿物)和石膏(蒸发盐矿物)的溶解, 那么(Ca2++Mg2+)/(HCO3-+SO42-)的值应该为1[38]. 然而, 研究区水样的(Ca2++Mg2+)/(HCO3-+SO42-)值大部分接近于1, 直接补给区地下水样点最靠近1∶1线, 表现出明显的碳酸盐矿物主导的特点[图 8(a)]. 此外, 从图 8(b)中可以看出(Ca2++Mg2+)与(HCO3-)的比值整体在1∶1线上方偏移, 表明除碳酸盐矿物溶解之外, Ca2+与SO42-还存在别的来源. 由于研究区部分含水层夹有石膏层[9], 结合(Ca2++Mg2+)/(HCO3-)和(Ca2++Mg2+)/(HCO3-+SO42-)值的差异性, 推测可能存在少量的石膏溶解作用生成了Ca2+与SO42-. 此外, 汇集出露区和间接补给区地下水样点有不少样本点偏离1∶1线且大部分比值大于1, 即Ca2++Mg2+多于HCO3-+SO42-, 表明阴离子有HCO3-和SO42-之外的其他来源构成.
图 8
Fig. 8
图 8 (Ca2++Mg2+)-(HCO3-+SO42-)与(Ca2++Mg2+)⁃HCO3-散点图
Fig. 8 (Ca2++Mg2+)-(HCO3-+SO42-)and (Ca2++Mg2+)-(HCO3-)scatter plots
自然状态下, 如果地下水中的Na+和Cl-仅来源于岩盐的溶解, 那么样品中Na+和Cl-的毫克当量比值为一个固定值1[9]. 从图 9(a)可以看出, 研究区大部分水样品的Na+和Cl-的相关性相对较好, 两者的比值接近于1, 且两者的相关系数达到了0.721(表 1), 表明研究区地下水中存在盐岩的溶解[39]. 从图 9(b)中可以看出, 大部分地下水样的Na+/Cl-值都不同程度地小于1, Cl-相对于Na+来说存在明显的过量, 表明可能存在阳离子交换从而使得地下水中的Na+浓度降低, 或者Cl-还有其他的来源. 此外, 间接补给区中的Cl-浓度相对较低, 直接补给区与汇集出露区水样具有相对较高的Cl-浓度(图 3), 并且表现出随着Cl-浓度的增加, Na+浓度有较明显的增长趋势[图 9(a)], 可能是汇集出露区含水层位于地下水排泄区, 随着地下水径流, 从间接补给区到汇集出露区地下水中的Na+和Cl-浓度逐渐累积, 同时也可能存在人为活动等外源输入.
图 9
Fig. 9
图 9 Na+-Cl-相关关系及其比值
Fig. 9 Correlation and ratio of Na+ and Cl-
3.3 控制因素分析
3.3.1 矿物溶解平衡
矿物饱和指数是水文地球化学研究的一个重要指标, 可以用来确定地下水与矿物的存在状态, 其公式为：
(1)
式中, SI为矿物饱和指数, IAP为矿物溶解反应中的离子活度积, Ksp为溶解反应的平衡常数. 当SI = 0时, 水与矿物达到溶解平衡状态；当SI > 0时, 矿物处于过饱和状态, 矿物可能会沉淀析出；当SI < 0时, 矿物处于非饱和状态, 矿物可能会继续溶解. 由于研究区含水层岩性以碳酸盐岩为主, 局部地层中夹有石膏层, 因此本文采用Phreeqc软件计算了地下水样品中不同矿物的饱和指数(图 10). 从图 10(a)中可以看出, 除少量间接补给区和汇集出露区水样品的方解石和白云石饱和指数为负值外, 其余样品的饱和指数均为正值, 即研究区方解石和白云石基本处于平衡或饱和状态. 此外, 所有的地下水样品的盐岩和石膏的饱和指数均小于0, 说明研究区盐岩和石膏均处于不饱和状态, 仍有继续溶解的趋势, 尤其是盐岩远远未饱和, 可能与含水层中盐岩矿物较少有关.
图 10
Fig. 10
图 10 不同矿物的饱和指数
Fig. 10 Saturation indices of different minerals
3.3.2 阳离子交换
阳离子交换可以通过氯-碱指数即Scholler(舍勒)指数来分析[40, 41], Schoeller(舍勒)指数CAI-I和CAI-II的计算公式为(离子浓度单位为meq·L-1)：
(2)
(3)
如果Schoeller指数为负值, 说明含水层成岩矿物中的Na+和K+被水中的Ca2+和Mg2+置换出来, Ca2+和Mg2+从溶液中析出附着在矿物中；如果得到的是正值, 说明发生了相反的阳离子交换, 即含水层成岩矿物中的Ca2+和Mg2+被Na+和K+取代. 根据公式(2)和公式(3)计算的Schoeller指数如图 11(a)所示, 可以看出, 仅有少部分地下水取样点的Schoeller指数为负值, 说明泉域内少部分地区存在Na+和K+取代水溶液中Ca2+和Mg2+的阳离子交换作用；几乎全部汇集出露区地下水样品的Schoeller指数为正, 说明可能存在Ca2+和Mg2+取代水溶液中Na+和K+的阳离子交换作用, 使得汇集出露区内Na+浓度降低. 通过(Ca2++Mg2+-HCO3--SO42-)⁃(Na++K+-Cl-NO3-)的关系可知[图 11(b)], 3个区域水样点均分布于斜率为-1(d = -1∶1)的线附近, 即地下水中增加(或减少)的Ca2+和Mg2+与增加(或减少)的Na+和K+的毫克当量相等, 且大多位于左上第二象限, 也可同样说明趵突泉泉域各功能区存在阳离子交换为主导的水文地球化学作用[42], 据前文分析可知, 研究区水样中的Ca2+和Mg2+除了碳酸盐岩和石膏矿物溶解之外有额外的来源, 推测可能是含水层矿物中的Ca2+和Mg2+与地下水中的Na+和K+发生了交换. 因此, 大部分地下水样的Na+/Cl-值小于1是由于阳离子交换使得地下水中Na+浓度降低, Cl-浓度相对升高所致.
图 11
Fig. 11
图 11 Scholler指数及(Ca2++Mg2+-HCO3-SO42-)⁃(Na++ K+-Cl--NO3-)关系
Fig. 11 Scholler index value and (Ca2++Mg2+-HCO3--SO42-)-(Na++K+-Cl--NO3-)relations plots
3.3.3 人类活动影响分析
人类活动与地下水化学组分的演化息息相关, 尤其是NO3-、Cl-和SO42-与人类活动产生的污染物质相关的组分, 因此可利用(SO42-/ Ca2+)⁃(NO3-/ Ca2+)和(Cl-/ Na+)⁃(NO3-/ Na+)的关系研究人类活动对地下水化学组分的影响[43, 44]. 由(SO42-/ Ca2+)-(NO3-/ Ca2+)的关系可知[图 12(a)], 部分汇集出露区和间接补给区地下水样品具有较高的NO3-/ Ca2+值, 说明这两个区域局部地下水受农业活动影响相对较大；直接补给区NO3-/ Ca2+值较低, 受农业活动影响相对较小；大部分样品的SO42-/ Ca2+值均高于NO3-/ Ca2+值, 说明研究区受工矿业活动的影响大于农业活动[45]. 此外, 间接补给区的大部分样品的SO42-/ Ca2+值高于汇集出露区, NO3-/ Ca2+值高于直接补给区, 说明间接补给区局部地区受人类活动的综合影响最为严重.
图 12
Fig. 12
图 12 (SO42-/ Ca2+)-(NO3-/ Ca2+)与(Cl-/ Na+)-(NO3-/ Na+)散点图
Fig. 12 (SO42-/ Ca2+)-(NO3-/ Ca2+)and(Cl-/ Na+)-(NO3-/ Na+)scatter plots
从(Cl-/ Na+)⁃(NO3-/ Na+)的关系中可以看出[图 12(b)], 研究区地下水样点多位于生活污水、降水和农业活动端元内部, 同样也可说明3个区域地下水均受人类活动影响较大, 其中部分间接补给区和直接补给区的地下水受生活污水的影响最为强烈, 并且间接补给区地下水样点有向农业活动影响区域迁移的趋势.
4 结论
(1) 趵突泉泉域各功能区地下水均为碱性, 表现出典型的碳酸盐矿物溶解特征, SOM分析结果显示各类地下水中的水化学参数具有较好的空间差异性. 由于地下水径流路径较长, 直接补给区水样具有相对较高的TDS、Ca2+、Na+、Mg2+、Sr2+、SO42-、HCO3-和NO3-；间接补给区地下水径流时间较短, 水岩作用不够充分, 水样各项参数值均比较低；汇集出露区水化学各参数的取值范围相对较大, 与该区处于地下水径流的排泄区、同时接受多水源的补给有关.
(2) 趵突泉泉域地下水阴离子以HCO3-和SO42-为主, 阳离子以Ca2+为主；从间接补给区到汇集出露区水化学类型由HCO3-Ca型向HCO3·SO4-Ca型转变, 含水层内存在少量的石膏溶解作用生成了Ca2+与SO42-.
(3) 趵突泉泉域地下水中离子主要来源于含水层岩石矿物的溶滤作用, 同时受阳离子交换作用、矿物溶解-平衡及人类活动的综合影响. 从间接补给区到汇集出露区地下水中的Cl-浓度逐渐累积, 可能是由于反向阳离子交换所致, 同时也存在人为活动等外源输入. 间接补给区地下水受人类活动综合影响较大, 人类活动的输入改变了泉域地下水水化学的演化.
感谢中国科协和济南城乡水务局对本研究工作的大力支持.