研究区为淄河流域（图1），其东发源于鲁山北麓，西发源于禹王山东麓，自西南向东北流经博山、淄川、临淄三区，干流全长179.7 km，流域面积1397 km²，是淄博市主要的水系之一. 该区域属暖温带大陆型季风气候，多年平均气温12.2℃，区内多年平均降水量672.68 mm（1959—2018年），多年平均蒸发量2008 mm. 最大年降水量1118.7 mm（1964年），最小年降水量319.5 mm（2006年），降水量年内分布不均.

淄河流域地势南高北低，东西两侧分别为与弥河及孝妇河地表分水岭，南部为低山丘陵，向北逐渐降低为山前倾斜平原，地形低洼向北呈倾伏状. 中部分布南北向淄河断裂带^[11]，总体走向NE35°，主断裂倾向SE，倾角60°—70°. 自淄河上游太河水库建成蓄水以来，下游河道来水量锐减，绝大多数年份均断流干涸.

该流域水文地质分区属于淄河地堑及太河水库上游单斜构造水文地质小区，区内寒武—奥陶纪碳酸盐岩广泛分布，岩性以中厚层灰岩及白云岩为主，南部分水岭一带分布岩浆岩. 淄河河床由卵砾石、砂砾石等松散堆积物组成，渗漏能力极强^[12]，素有“淄河十八漏”之说. 沿淄河一带下伏奥陶系灰岩受断裂构造影响，裂隙岩溶发育，具有极强的渗透性，主要接受大气降水入渗和淄河河道渗漏补给，东、西、南三面地下水、地表水汇集后经淄河断裂带源源不断渗入地下，迅速补给岩溶地下水. 自淄河上游谢家店至下游大武一带分布多个具有强富水性的岩溶型地下水水源地，单井涌水量一般大于1000 m³·d⁻¹，局部如大武、刘征水源地单井涌水量均大于5000 m³·d⁻¹. 该流域地下水主要富水类型为碳酸盐岩类裂隙岩溶水，为本次取样研究主要目的含水层.

调查采样于2018年6月（枯水期）、9月（丰水期），共采集地下水样品70个（图1）. 采样现场用水位计测定水位，用便携式HACH水质分析仪测定地下水的pH、电导率（ EC） 、溶解氧（DO）. ${\mathrm{N}\mathrm{a}}^{+}$、${\mathrm{K}}^{+}$、${\mathrm{C}\mathrm{a}}^{2+}$、${\mathrm{M}\mathrm{g}}^{2+}$等阳离子测定通过电感耦合等离子体发射光谱仪（ICP-OES，OPTIMA 7000DV型），${\mathrm{C}\mathrm{l}}^{-}$、${{\mathrm{S}\mathrm{O}}_{4}}^{2-}$、${{\mathrm{N}\mathrm{O}}_{3}}^{-}$浓等阴离子通过CIC-D120离子色谱仪完成，${{\mathrm{C}\mathrm{O}}_{3}}^{2-}$、${{\mathrm{H}\mathrm{C}\mathrm{O}}_{3}}^{-}$、TDS采用EDTA滴定法测定，通过主离子的浓度总和减去1/2${{\mathrm{H}\mathrm{C}\mathrm{O}}_{3}}^{-}$计算得到溶解性总固体（ TDS）. 所有样品在山东省地矿工程勘察院实验室完成检测.

测试完毕后，对水化学数据进行分析整理. 利用分析统计软件SPSS25.0对地下水水化学离子进行特征统计及因子分析；利用AQQA软件绘制piper三线图；利用excel软件绘制Gibbs图及离子比例关系图；利用R统计软件PerformanceAnalytics包进行相关性可视化分析.

对枯水期和丰水期地下水水化学成分及浓度进行特征统计，结果如表1所示. 总体而言，淄河流域岩溶地下水主要离子在枯水期、丰水期呈现一定的差异性. 地下水pH值处于7.1—8.3之间，呈弱碱性，变异系数较小，均值相差不大，季节差异不显著. 枯水期和丰水期TDS平均值分别为498.54 mg·L⁻¹ 、484.09 mg·L⁻¹，丰水期TDS值略小于枯水期. 地下水中${\mathrm{C}\mathrm{l}}^{-}$、${\mathrm{M}\mathrm{g}}^{2+}$、${{\mathrm{H}\mathrm{C}\mathrm{O}}_{3}}^{-}$、${{\mathrm{N}\mathrm{O}}_{3}}^{-}$平均质量浓度均表现为丰水期>枯水期，一般来说，降水进入含水层后会对地下水有稀释作用，而淄河流域相反，丰水期离子含量略高于枯水期，表明大气降水可能加快了岩石的溶滤作用，例如：碳酸盐岩的溶解会使得丰水期${{\mathrm{H}\mathrm{C}\mathrm{O}}_{3}}^{-}$离子含量偏高；地表硝酸盐污染物随雨水入渗地下水中，可能导致${{\mathrm{N}\mathrm{O}}_{3}}^{-}$含量升高. 且由于区域水文地质条件差异性使得不同采样点水质具有较大的变化范围.

Piper三线图可以用来分析岩溶地下水主要离子组成的总体特征及演化规律，该方法的优点是不受人为因素的影响，可以看出各离子的相对含量^[13-14]. 由图2可知，该流域岩溶地下水以${\mathrm{H}\mathrm{C}\mathrm{O}}_{3}-\mathrm{C}\mathrm{a}$型及${\mathrm{H}\mathrm{C}\mathrm{O}}_{3}·{\mathrm{S}\mathrm{O}}_{4}-\mathrm{C}\mathrm{a}$型为主，不论枯、丰水期优势阴离子均为${{\mathrm{H}\mathrm{C}\mathrm{O}}_{3}}^{-}$，分别占全部阴离子总量的45.65%—75.28%，50.25%—76.37%，其次为${{\mathrm{S}\mathrm{O}}_{4}}^{2-}$. 就阳离子而言，枯水期和丰水期的优势阳离子均为${\mathrm{C}\mathrm{a}}^{2+}$离子，分别占全部阳离子总量的55.60%—81.71%，48.25%—81.40%，说明碳酸盐岩风化溶解是水化学组分的主要来源^[15]. 究其原因是由于研究区岩溶发育较强，在灰岩的风化和溶滤作用下，致使该区地下水中${{\mathrm{H}\mathrm{C}\mathrm{O}}_{3}}^{-}$、${\mathrm{C}\mathrm{a}}^{2+}$普遍偏高. 根据前人研究，${{\mathrm{S}\mathrm{O}}_{4}}^{2-}$主要来源于淄河上游寒武系灰岩石膏层、芒硝等的溶解及石马镇一带局部煤系地层黄铁矿（${\mathrm{F}\mathrm{e}\mathrm{S}}_{2}$）的氧化^[16]，近年来，${{\mathrm{S}\mathrm{O}}_{4}}^{2-}$在地下水中占比逐渐增高^[17]，推测原因为人类开采活动加速了水岩作用过程.

从研究区TDS等值线分布图（图3）可知，淄河流域岩溶地下水TDS自上游至下游呈增高的趋势，在博山区北博山镇、源泉镇、崮山镇一带及下游的临淄区南王镇、金岭回族镇等一带，居民及工农业密集，且分布多个地下水水源地，地下水TDS明显增高，主要影响离子为总硬度、${\mathrm{C}\mathrm{l}}^{-}$、${{\mathrm{S}\mathrm{O}}_{4}}^{2-}$、${{\mathrm{N}\mathrm{O}}_{3}}^{-}$. 根据大武水源地供水水文地质勘探报告，淄河下游大武水源地一带在开采前的1977年TDS约200—360 mg·L⁻¹，${{\mathrm{N}\mathrm{O}}_{3}}^{-}$浓度仅为4—6.8 mg·L⁻¹；而至2018年枯、丰水期，该区域岩溶地下水中TDS含量已增加到510—650 mg·L⁻¹，${{\mathrm{N}\mathrm{O}}_{3}}^{-}$含量已达38—55 mg·L⁻¹，临淄区堠皋一带水质较差，TDS超过1000 mg·L⁻¹、总硬度、${\mathrm{C}\mathrm{l}}^{-}$、${{\mathrm{S}\mathrm{O}}_{4}}^{2-}$含量均明显高于流域其它地区，这固有地下水从补给区到排泄区TDS自然增高的因素，但更主要是由于该区域石油化工企业密集，工业活动较频繁^[18]，人类活动已经对地下水水质产生显著影响.

Gibbs（1970）模型能够清楚地表示出控制地下水化学成分的自然来源以及水化学形成机制，其中形成机制分为降水作用机制、岩石风化机制以及蒸发浓缩机制^[19-20].

从研究区岩溶地下水的Gibbs图4可知，枯、丰水期地下水样品TDS介于335—1495 mg·L⁻¹、294—1370 mg·L⁻¹，${\mathrm{N}\mathrm{a}}^{+}/({\mathrm{N}\mathrm{a}}^{+}+{\mathrm{C}\mathrm{a}}^{2+})$比值分别在0.01—0.74、0.05—0.88范围内，${\mathrm{C}\mathrm{l}}^{-}/({\mathrm{C}\mathrm{l}}^{-}+{\mathrm{H}\mathrm{C}\mathrm{O}}_{3}^{-})$比值分别在0.03—0.51、0.02—0.38范围内. 枯水期、丰水期采样点在均分布在岩石风化控制区附近，并远离降水控制型区域. 这表明该流域岩溶地下水离子组分主要受岩石风化作用控制. 该区分布多个地下水水源地，开采量较大，地下水水位埋深较深，受大气降雨和蒸发的作用较小，但受构造条件的影响，裂隙岩溶发育，为岩石的风化提供了有利条件. 由于地下水离子来源的广泛性和复杂性，需要通过其他方法进一步深入探讨.

相关分析是研究不同变量相关性的数理统计方法，将相关性可视化可清楚地看出变量之间的相关性强弱，从而有效的揭示地下水主要组分来源的一致性和差异性^[21-22]. 相关性分析可视化结果如图5所示，图中对角线是变量自身的分布；对角线右上方表示两个属性的相关性值，星号表示显著程度，*在0.05 级别（双侧），相关性显著，**表示在 0.01 级别（双侧），相关性显著；对角线的左下方是两个属性的散点图.

TDS与各离子相关关系能较好反映地下水的成因. 由图5可知，除${\mathrm{K}}^{+}$外，TDS与各离子之间均呈显著相关. 在枯水期，${\mathrm{C}\mathrm{a}}^{2+}$与TDS相关性最强（r=0.90），其次为${{\mathrm{S}\mathrm{O}}_{4}}^{2-}$、${\mathrm{C}\mathrm{l}}^{-}$、${\mathrm{M}\mathrm{g}}^{2+}$（r＞0.7），${\mathrm{N}\mathrm{O}}_{3}^{-}$、${\mathrm{H}\mathrm{C}\mathrm{O}}_{3}^{-}$相关系数分别为0.48、0.39；在丰水期，${\mathrm{C}\mathrm{l}}^{-}$与TDS相关性最强（r=0.90），其次为${\mathrm{C}\mathrm{a}}^{2+}$、${{\mathrm{S}\mathrm{O}}_{4}}^{2-}$，表明研究区地下水主要离子的来源可能与广泛分布的灰岩、白云岩等可溶性碳酸盐岩的溶解有关. 在枯水期主要以碳酸盐岩溶解为主，到丰水期地下水位抬升，增加了岩盐、硫酸盐等矿物的溶解.

${\mathrm{N}\mathrm{a}}^{+}$与${\mathrm{C}\mathrm{l}}^{-}$相关性最好（r≈0.71），其次为${{\mathrm{S}\mathrm{O}}_{4}}^{2-}$，表明${\mathrm{N}\mathrm{a}}^{+}$的来源与${\mathrm{C}\mathrm{l}}^{-}$具有相似性，可能与岩盐矿物的溶解有关；${\mathrm{C}\mathrm{a}}^{2+}$与${{\mathrm{S}\mathrm{O}}_{4}}^{2-}$相关性最好（r≈0.76）、其次为${\mathrm{C}\mathrm{l}}^{-}$、${\mathrm{M}\mathrm{g}}^{2+}$、${\mathrm{N}\mathrm{O}}_{3}^{-}$，说明除白云岩水岩作用、石膏矿物溶解以外，部分${\mathrm{C}\mathrm{a}}^{2+}$的来源与含${\mathrm{C}\mathrm{l}}^{-}$、${\mathrm{N}\mathrm{O}}_{3}^{-}$的污染物入渗有关.

${\mathrm{M}\mathrm{g}}^{2+}$和${{\mathrm{S}\mathrm{O}}_{4}}^{2-}$相关性（r=0.63、0.40）表明存在脱白云石化作用. 由于石膏溶解未达到饱和，水中${\mathrm{C}\mathrm{a}}^{2+}$浓度增加，并产生含${{\mathrm{S}\mathrm{O}}_{4}}^{2-}$溶液，白云石中${\mathrm{M}\mathrm{g}}^{2+}$被${\mathrm{C}\mathrm{a}}^{2+}$置换产生方解石沉淀，发生脱白云石化.

相比于传统的水化学类型分析方法，分析水化学组分的比例关系特征更能有效的判断地下水的成因、水化学成分的来源及演化过程^[23]. 假设岩溶地下水中${\mathrm{C}\mathrm{a}}^{2+}$、${\mathrm{M}\mathrm{g}}^{2+}$主要来自方解石、白云石、石膏的溶解及蒸发岩的溶解，所以可用${(\mathrm{C}\mathrm{a}}^{2+}+{\mathrm{M}\mathrm{g}}^{2+})$与$({\mathrm{H}\mathrm{C}\mathrm{O}}_{3}^{-}$+${{\mathrm{S}\mathrm{O}}_{4}}^{2-}$）摩尔浓度关系来判断${\mathrm{C}\mathrm{a}}^{2+}$、${\mathrm{M}\mathrm{g}}^{2+}$的主要来源. 方解石、白云石和石膏的溶解方程如下：

若地下水中相关离子主要来源于方解石、白云石和石膏的溶解，则水中${(\mathrm{C}\mathrm{a}}^{2+}+{\mathrm{M}\mathrm{g}}^{2+})$与$({\mathrm{H}\mathrm{C}\mathrm{O}}_{3}^{-}$+${{\mathrm{S}\mathrm{O}}_{4}}^{2-}$）的毫克当量比应该接近1:1^[24]，由图6（a）可以看出，研究区岩溶地下水在枯、丰水期绝大部分采样点位于1:1直线上侧，说明地下水中溶解组分季节差异较小，主要与地质背景条件有关，地下水中${\mathrm{C}\mathrm{a}}^{2+}$和${\mathrm{M}\mathrm{g}}^{2+}$组分主要来源于碳酸盐岩的溶解作用. 位于1:1直线下侧的点位主要分布在太河乡、太河水库下游附近，说明可能伴随着硅酸盐岩和蒸发岩的溶解作用^[25].

${\mathrm{K}}^{+}$、${\mathrm{N}\mathrm{a}}^{+}$和${\mathrm{C}\mathrm{l}}^{-}$主要来源于大气降水、生活废水、农业施肥和盐岩的溶解^[26]. 当盐岩溶解起主导作用时，${(\mathrm{K}}^{+}+{\mathrm{N}\mathrm{a}}^{+})$与${\mathrm{C}\mathrm{l}}^{-}$的比值应为1:1^[27]. ${(\mathrm{K}}^{+}+{\mathrm{N}\mathrm{a}}^{+})$与${\mathrm{C}\mathrm{l}}^{-}$的关系如图6（b）所示，可见${\mathrm{K}}^{+}$、${\mathrm{N}\mathrm{a}}^{+}$和${\mathrm{C}\mathrm{l}}^{-}$受季节影响较大. 枯水期采样点大部分位于1:1直线上方，说明枯水期时该区主要存在盐岩溶解作用和铝硅酸盐矿物的溶滤. 丰水期采样点大部分位于1:1直线下方，说明丰水期时地下水中${\mathrm{C}\mathrm{l}}^{-}$可能有其他来源，如生活污水和农业面源污染等.

通常用$\left[\right({\mathrm{C}\mathrm{a}}^{2+}+{\mathrm{M}\mathrm{g}}^{2+})-({\mathrm{H}\mathrm{C}\mathrm{O}}_{3}^{-}+{\mathrm{S}\mathrm{O}}_{4}^{2-}\left)\right]$与$\left[\right({\mathrm{N}\mathrm{a}}^{+}+{\mathrm{K}}^{+})-{\mathrm{C}\mathrm{l}}^{-}$]的关系来研究地下水中阳离子交换作用，$\left[\right({\mathrm{C}\mathrm{a}}^{2+}+{\mathrm{M}\mathrm{g}}^{2+})-({\mathrm{H}\mathrm{C}\mathrm{O}}_{3}^{-}+{\mathrm{S}\mathrm{O}}_{4}^{2-}\left)\right]$表示除石膏、方解石和白云石溶解外所带来的${\mathrm{C}\mathrm{a}}^{2+}$和${\mathrm{M}\mathrm{g}}^{2+}$的增加和减少量，$\left[\right({\mathrm{N}\mathrm{a}}^{+}+{\mathrm{K}}^{+})-{\mathrm{C}\mathrm{l}}^{-}$]表示除岩盐溶解外所带来的${\mathrm{N}\mathrm{a}}^{+}$的增加或减少量^[28]. 如果$\left[\right({\mathrm{C}\mathrm{a}}^{2+}+{\mathrm{M}\mathrm{g}}^{2+})-({\mathrm{H}\mathrm{C}\mathrm{O}}_{3}^{-}+{\mathrm{S}\mathrm{O}}_{4}^{2-}\left)\right]$与$\left[\right({\mathrm{N}\mathrm{a}}^{+}+{\mathrm{K}}^{+})-{\mathrm{C}\mathrm{l}}^{-}$]呈显著负相关关系，那么各采样点均落在斜线-1的直线附近，表明地下水在径流过程中发生了阳离子交换作用.

由图6（c）可见，大部分水样点落在斜线-1的直线附近及上方，这说明研究区地下水中存在阳离子交换作用，水中${\mathrm{C}\mathrm{a}}^{2+}$和${\mathrm{M}\mathrm{g}}^{2+}$会把岩土中吸附的${\mathrm{N}\mathrm{a}}^{+}$和${\mathrm{K}}^{+}$置换出来，再结合由相关性的可视化图3中可知，${\mathrm{C}\mathrm{a}}^{2+}$与${\mathrm{N}\mathrm{a}}^{+}$、${\mathrm{M}\mathrm{g}}^{2+}$与${\mathrm{N}\mathrm{a}}^{+}$的相关系数分别为0.26和 0.18，相关性不显著，说明阳离子交换作用在地下水化学组分构成中起到不可忽视的作用.

有研究表明${\mathrm{N}\mathrm{O}}_{3}^{-}$、${\mathrm{C}\mathrm{l}}^{-}$来源于含水层介质的溶解较少，参与大规模的离子交换反应也不多，而主要来自人类活动产生的生活污水、农业无机化肥的使用及人畜粪便^[29]. 由图6（d）以及相关性可视化图可知，${\mathrm{N}\mathrm{O}}_{3}^{-}$与${\mathrm{C}\mathrm{l}}^{-}$相关性不显著，说明它们有各自不同来源，可能来自生活污水和农业活动等不同的污染物.

影响地下水水化学组分的变量很多，变量研究已经广泛用于地下水水化学组成及影响因素，而因子分析能用几个变量反映多个变量的内在的联系和本质特征^[30-31]. 通过对研究区岩溶地下水中主要阴阳离子（${\mathrm{N}\mathrm{a}}^{+}$、${\mathrm{K}}^{+}$、${\mathrm{M}\mathrm{g}}^{2+}$、${\mathrm{C}\mathrm{a}}^{2+}$、${\mathrm{C}\mathrm{l}}^{-}$、${\mathrm{S}\mathrm{O}}_{4}^{2-}$、${\mathrm{H}\mathrm{C}\mathrm{O}}_{3}^{-}$、${\mathrm{N}\mathrm{O}}_{3}^{-}$）进行KMO 检验和Bartlett 球形检验统计，检验结果显示，枯、丰水期KMO 检验取值分别为0.516和0.521，Bartlett 球形检验统计量的 Sig值均为零（Sig=0＜0.01），有研究表明 KMO统计量在 0.5 以上，Bartlett 球形检验统计量的Sig=0（＜0.01） 时，因子分析法可用于分析初始变量内在联系. 因此选取研究区地下水中主要的阴阳离子（${\mathrm{N}\mathrm{a}}^{+}$、${\mathrm{K}}^{+}$、${\mathrm{M}\mathrm{g}}^{2+}$、${\mathrm{C}\mathrm{a}}^{2+}$、${\mathrm{C}\mathrm{l}}^{-}$、${\mathrm{S}\mathrm{O}}_{4}^{2-}$、${\mathrm{H}\mathrm{C}\mathrm{O}}_{3}^{-}$、${\mathrm{N}\mathrm{O}}_{3}^{-}$）进行主成分分析和因子分析，通过分析结果进一步把握研究区自然条件和人为活动等因素对地下水水化学组成的影响^[32]. 将主成分分析和因子分析进行整理，其分析结果如表2、3所示.

以特征值大于1为提取原则，符合要求的共三个主成分，提取分析结果共同解释了原数据74.567%的方差. 丰水期时，第1主成分解释了37.21%的方差，该成分的主要贡献者是${\mathrm{C}\mathrm{a}}^{2+}$、${\mathrm{S}\mathrm{O}}_{4}^{2-}$、${\mathrm{H}\mathrm{C}\mathrm{O}}_{3}^{-}$，反映了以方解石和石膏等碳酸盐溶解对岩溶地下水组分的影响. 第2主成分的解释了21.38%的方差，该成分的主要贡献者是${\mathrm{N}\mathrm{a}}^{+}$、${\mathrm{K}}^{+}$、${\mathrm{C}\mathrm{l}}^{-}$、${\mathrm{N}\mathrm{O}}_{3}^{-}$，反映了人为活动和大气输出的影响. 第3主成分解释了15.98%的方差，该成分的主要贡献者是${\mathrm{M}\mathrm{g}}^{2+}$，反映了白云岩溶解以及脱白云石化等水岩作用. 故可认为研究区地下水水化学组分和演化进程主要受到这三种因素的影响.

结合因子的得分图7和采样点的分布情况可以直观解释第1因子、第2因子与原变量的关系，较好的反映采样点离子的来源与影响因素. 总体而言，枯水期、丰水期的地下水在图中分布均匀，但分布在左下角较多，枯、丰水期均占到34.21%. 左下角地下水离子的因子的第1、2因子得分为负值，说明第1、2因子影响较小，结合采样图可知采样位于淄河支流与淄河交汇的上游、太河水库及仁河水库下游，在交汇的河流段长度不长，水岩作用时间短，加之水库周围保护措施较好，人为活动影响较小；左上角第1因子为负值，第2因子为正值，说明地下水化学组分是由第2因子代表的人为活动和大气输出的影响所控制的，采样点分别位于陡沟村、拔子村以及淄河下游的南仇村至前李官庄村河段. 枯、丰水期分布于右下角水样较少，只分别占到18.42%、13.16%，右下角第1因子为正，第2因子为负，则主要由第1因子代表的碳酸盐溶解的影响所控制，采样点分布也较散，无规律可循. 右上角第1因子为正，第2因子为正，则由碳酸盐溶解、人为活动和大气输出影响共同控制，这些采样点位于淄河上游，上游有谢家店水源地、源泉水源地等.

（1）研究区岩溶地下水枯水期和丰水期的TDS平均值分别为498.54 mg·L⁻¹ 、484.09 mg·L⁻¹，丰水期TDS值略小于枯水期. 在地下水中${\mathrm{C}\mathrm{l}}^{-}$、${\mathrm{M}\mathrm{g}}^{2+}$、${\mathrm{H}\mathrm{C}\mathrm{O}}_{3}^{-}$、${\mathrm{N}\mathrm{O}}_{3}^{-}$平均质量浓度均表现为丰水期>枯水期. 一般来说，降水进入含水层后会对地下水有稀释作用，但淄河流域相反，丰水期离子含量略高于枯水期，表明降水作用加快了岩石的溶滤作用.

（2）分析piper三线图可知，淄河流域岩溶地下水水化学类型以${\mathrm{H}\mathrm{C}\mathrm{O}}_{3}-\mathrm{C}\mathrm{a}$型及${\mathrm{H}\mathrm{C}\mathrm{O}}_{3}·{\mathrm{S}\mathrm{O}}_{4}-\mathrm{C}\mathrm{a}$型为主. TDS自上游至下游呈增高的趋势，符合地下水演化自然规律. 但在乡镇一带居民及工农业聚集区，TDS明显增高，主要影响离子为总硬度、${\mathrm{C}\mathrm{l}}^{-}$、${\mathrm{S}\mathrm{O}}_{4}^{2-}$、${\mathrm{N}\mathrm{O}}_{3}^{-}$. 淄河下游临淄区堠皋一带，水质较差，反映化工企业对该区域地下水水质造成影响.

（3）运用Gibbs模型及相关性分析等方法得出，淄河流域岩溶地下水离子组分主要受岩石风化作用控制，与广泛分布的灰岩、白云岩等可溶性碳酸盐岩的溶解有关. 在枯水期主要以碳酸盐岩溶解为主，到丰水期地下水位抬升，增加了岩盐、硫酸盐等矿物的溶解. 地下水中存在脱白云石化作用.

（4）运用离子比例分析可知，岩溶地下水中${\mathrm{C}\mathrm{a}}^{2+}$和${\mathrm{M}\mathrm{g}}^{2+}$组分主要来源于碳酸盐岩的溶解作用；${\mathrm{C}\mathrm{l}}^{-}$、${\mathrm{N}\mathrm{O}}_{3}^{-}$受人类活动影响较大，可能来自生活污水和农业活动等不同的污染物，并且丰水期加速了污染物质的溶解；阳离子交换作用在地下水化学组分构成中起到不可忽视的作用.

（5）因子分析法提取了三个主因子，能够较好地解释淄河流域岩溶水的主要影响因素. 第1主成份主要贡献者是${\mathrm{C}\mathrm{a}}^{2+}$、${\mathrm{S}\mathrm{O}}_{4}^{2-}$、${\mathrm{H}\mathrm{C}\mathrm{O}}_{3}^{-}$，反映了方解石和石膏等碳酸盐溶解对岩溶地下水组分的影响，贡献率37.21%；第2主成分主要贡献者是${\mathrm{N}\mathrm{a}}^{+}$、${\mathrm{K}}^{+}$、${\mathrm{C}\mathrm{l}}^{-}$、${\mathrm{N}\mathrm{O}}_{3}^{-}$，反映了人为活动和大气输出的影响，贡献率21.38%；第3主成分主要贡献者是${\mathrm{M}\mathrm{g}}^{2+}$，反映了白云岩溶解以及脱白云石化等水岩作用对岩溶地下水组分的影响，贡献率15.98%.