淄河源区浅层地下水化学特征及主要离子来源研究

游 京，齐跃明，邵光宇，马 超，杨雅琪，裴毅峰

(1.中国矿业大学 资源与地球科学学院，江苏 徐州 221116；2.上海市地矿工程勘察院，上海 200072；3.山东省地矿工程勘察院，山东 济南 250014)

地下水的化学成分是地下水与外界环境长期相互作用的产物[1]，开展地下水化学特征及主要成分来源的研究有助于阐明地下水的来源与形成过程，揭示地下水环境状况，同时也为地下水开发与保护提供科学合理的依据。岩溶地下水是淄博市博山区石马镇主要供水水源，在保障供水、支持社会经济发展及维持生态平衡等方面具有重要地位。近年来全国关于地下水化学成分及水质方面的研究较多，刘伟江等[2]研究表明贵阳市地下水水化学演化主要受水-岩作用的控制，碳酸盐矿物以及第四系硅酸岩和蒸发岩的风化和溶滤是区域地下水化学组分的主要来源，为贵阳市地下水污染治理提供了依据。王珺瑜等[3]综合运用水化学和多元统计方法，分析了济南泉域地下水水化学特征，探讨了不同区域水质影响因素及影响强度，结果表明人类活动对当地地下水造成了一定程度的污染，对济南泉域岩溶水资源开发利用及生态环境保护具有十分重要的意义。董维红等[4]为改善鄂尔多斯盆地地下水水质恶化、水污染蔓延发展等问题，采用Schoeller图和Piper三线图两种方法进行地下水类型分类，得出鄂尔多斯白垩系地下水盆地地下水化学类型具有明显的水平分带规律，为深入研究区域地下水的循环特征和水动力场特征奠定基础，对区域地下水资源可持续开发利用和综合管理具有指导作用。李贵恒等[5]调查分析了泰莱盆地地下水水化学特征及水质现状，结果表明其水化学类型以 HCO3·SO4-Ca型水为主；pH值整体呈弱碱性，TDS从大到小顺序为孔隙水、岩溶水、裂隙水。袁建飞等[6]研究了西昌市南部岩溶地下水系统，得出了该区域岩溶水、基岩裂隙水和地表水水化学组分、氢氧同位素和氚同位素特征，并探讨了这3类水的成因和溶质演化过程。王广才，沈照理等[7]运用质量平衡法研究平顶山矿区岩溶水的组分，得到了岩溶水化学成分的特征及形成机理。齐跃明等[8]通过中美交流互访、水资源研究现状调查和文献统计分析，对比研究了中美水资源研究的特点、发展现状及展望，为未来我国水资源研究和利用决策提供了科学依据。Zhou等[9]运用K-XAJ耦合模型模拟漓江流域的降雨-径流过程，结果较新安江模型更为准确。Kalhor等[10]综述了岩溶地下水中水质评价方法和修复技术，探讨了地表水和地下水之间的交互作用，为地下水利用和管理提供了研究基础。Sappa等[11]基于环境同位素 (18O、2H)和水化学调查技术，建立了拉丁美洲南部的碳酸盐岩含水层的水文地质概念模型，确定了水流的补径排条件以及控制地下水演化的化学作用过程，为当地水资源管理和保护提供了大量水文地球化学信息，拓展了研究区所欠缺的水文地球研究。

总体来说，前人对各自研究区域的水文地球化学(如浅层地下水特征、运移规律、水-岩相互作用及地下水污染治理等方面)都做了较好的研究工作，为之后更加深入的研究或者相关工程实施奠定了基础。然而，对于淄河源区来说，基础性、系统性、较多组分的岩溶水水化学组分特征研究工作未见报道，而这个工作恰恰对于源区的水资源合理开发利用及当地水环境、生态环境变迁的研究具有较大的理论和实际意义。

鉴于此，本文在山东省淄博市博山区石马镇谢家店村拟建水源地前期水文地质勘查和水质分析的基础上，结合淄河源区当地的地质和水文地质条件，通过采集不同区域岩溶水样品，综合运用水文地球化学(Piper三线图、离子比例系数、相关分析等)分析方法，初步确定区域地下水化学特征及主要离子来源，为之后的岩溶水化学分区提供参考，也有助于确定水污染分析背景场，为未来规划应急地下水源地、防止水污染以及合理开发与保护岩溶地下水提供科学依据。

1 淄河源区概况

1.1 自然地理条件

淄河源区[12]位于山东省淄博市博山区南部，北抵莱芜青石关村—崮山—源泉镇一线，东南及南部至博山区与沂源区行政边界，即地表分水岭，东北部为博山区与淄川区行政边界，西到莱芜市变质岩山区地表分水岭(图1所示)。行政区划涉及莱芜市莱城区与淄博市博山区。地理坐标：北纬36°15′43″～36°26′58″；东经117°48′44″～118°12′44″，面积605 km2。源区水文地质调查和水源勘探井工程主要位于博山区博山镇的北博山村—谢家店—邀兔崖村一带，属温带半湿润大陆性季风气候，多年平均气温12 ℃，平均降水量720.6 mm。地貌类型为周边中低山，中部为淄河河谷地形。区内地表水主要有石马水库及季节性河流淄河。淄河西北支流发源于鲁山山脉，西南支流发源于禹王山麓，两支流在源泉镇泉河村汇流，并入淄河主干道，继续向北径流。

1.2 构造

自燕山运动以来，由于多次构造作用的结果，区内的构造行迹主要表现为：断裂构造较为发育，褶皱比较少见，特别是张性断裂发育广泛，规模亦较大。区内较大断裂主要有青龙山断裂、龙头山断裂、石马断裂、盆泉—北博山断裂、甘泉断裂等。

1.3 水文地质条件

研究区地下水的补给几乎全部来自于大气降水，地下水流场形态由地表地势和断裂构造控制，由于北面石马断裂为一阻水断层，地下水径流方向是由研究区的西南部、南部向淄河河谷方向径流，淄河断裂带(盆泉—北博山断裂)为该区地下水的强径流带，地下水向淄河河谷方向汇流之后继续沿淄河地堑向北东方向径流。地下水的排泄方式主要有：(1)人工开采；(2)河道排泄；(3)自流井、泉及地下水溢出排泄等。

2 浅层地下水化学特征

2.1 水样采集与测定

2.2 水化学特征分析

表1 淄河源区岩溶水水化学分析结果统计

本文采用舒卡列夫分类法[13]得出HCO3·SO4-Ca·Mg型水为研究区主要水化学类型，占比为总水样数的71.9%，水化学类型演化总体表现为由HCO3-Ca型水转化为HCO3·SO4-Ca型，再转变为HCO3·SO4-Ca·Mg型的过程。第四系松散岩类孔隙水含水岩组地下水的水化学类型以HCO3·SO4-Ca·Mg型为主，块状岩类风化裂隙水以HCO3·SO4-Ca型为主，碳酸盐岩岩溶裂隙水以HCO3·SO4-Ca·Mg型为主。本文将不同类型水质成分分别投影到Piper三线图(图2)中，可以发现3种地下水类型基本相似，主体区分为HCO3·SO4-Ca·Mg型和HCO3·SO4-Ca型水，结合研究区大型断裂构造较多的情况，本文推测该地区各个含水层之间存在着较好的水力联系。

图2 研究区岩溶地下水水化学Piper图

3 主要离子来源分析

3.1 离子浓度比值分析

图3 研究区地下水离子浓度比值

碳酸盐岩与硫酸盐岩溶解化学反应式如式(1)、式(2)所示：

(1)

(2)

3.2 相关性分析

地下水化学成分由于形成来源、流动途径、沿途所发生的各种水文地球化学作用等的差异，而有较大差异。而同一来源、流动途径相似、水化学作用相似的地下水则往往具有相似的水化学成分特征。根据这一特点，我们可通过地下水化学成分之间的相关性分析推断地下水在补给来源、流动途径或水文地球化学作用等方面的差异性。本文运用SPSS软件计算了研究区地下水11项水化学指标间的Pearson相关系数[18]，并在0.05的显著水平上对相关系数进行显著性双尾检验，得到相关矩阵如表2所示。

表2 岩溶水各水化学组分间的相关系数

相关性检验时，通常认为：相关系数r=0时，表示无相关性；0 < |r| < 0.4，表示低度线性相关；0.4 < |r| < 0.7时，表示显著线性相关；0.7 < |r| < 1，表示高度线性相关；|r|=1，表示完全相关[19]。

4 结论

3)地下水水质类型和低TDS值表明，研究区地下水溶滤作用强烈，浅层地下水水质主要与地下水流动条件及碳酸盐岩裸露或浅埋藏有关。区域地下水水质较好，源区生态环境保护较好，满足建设应急水源地的基本水质条件。

5 讨论

1)由于现实原因，第四系松散岩类孔隙含水层和块状岩类风化裂隙含水层所取得的水样偏少，使得两含水层在水化学特征分析结果中代表性可能不足。并且由于打一口基岩井成本高，采样点选取需依赖当地勘探井和民用井来进行取水工作，难以在空间上实现理论上的均匀布置，研究区东北、东南、南及西南端的水样与岩样数据缺乏，在条件允许的基础上，今后将加大东北及东南区域的水质水样的搜集工作。另一方面，还可在区域同位素分析、特征微量元素分析上展开进一步的研究工作，以从多角度认识区域地下水化学特征，全面展现研究区地下水环境的真实面貌。

2)大采样时间还可考虑增加不同时间(如枯水期与丰水期，夏季与冬季)的水样采取工作，以便进行对比分析，分析地下水化学演化的控制因素。另一方面，在地下水合理利用方面应该考虑该区域雨季和旱季分明、淄河源区岩溶地下水库容较大的特点，利用补偿疏干法等来评价当地的水资源量，更好地改善当地的水资源开发利用问题。