泰莱盆地地下水化学特征分类及成因分析

田明刚，徐 建，赵 耘，郑梦琪，徐秋林

（1.中国冶金地质总局山东正元地质勘查院，济南250000；2.山东正元建设工程有限责任公司，济南250000）

地下水化学特征空间分布规律蕴含了地下水在循环过程中的赋存环境、径流途径，能量质量交换等重要信息。地下水在演变过程中受到蒸发、降雨、水岩相互作用、人为因素等影响，其组分在不断的发生变化。分析水化学特征，揭示控制因素，可为地下水资源保护和可持续开发提供科学依据[1-4]。

泰莱盆地位于鲁中地区，区域内工农业发达，矿产资源丰富，但水资源相对短缺，成为制约当地经济发展的重要因素。随着环保力度不断加大，合理利用水资源对地区可持续发展具有重要意义。盆地四周多火成岩山地，地下水以基岩裂隙水为主，盆地中部深层岩溶发育，地下水以岩溶水和裂隙岩溶水为主，近地表为第四系孔隙水[5]。相比地表水，盆地内地下水可利用性更高。2017年，李贵恒等利用丰水期水化学数据，结合数理统计、piper 三线图等方法，对泰莱盆地地下水进行了特征分类和水质评价，指出该地区地下水具有很高的利用价值[6]。冯建国等针对该地区孔隙水进行了细化研究，分析得出水化学组分主要受水-岩作用控制，以硅酸盐岩溶解为主，并伴随碳酸盐溶解。在南、北部基岩区还受裂隙水和岩溶水补给影响[7]。该地区农业发达，地下水除受岩性控制以外还受种植施肥、喷洒农药等人为因素的控制[8]，为减少控制因素，此研究在枯水期取样测试，研究盆地各类地下水水化学特征及其控制因素。

1 研究区概况

1.1 自然概况

研究区位于山东省中南部，包括泰安市、莱芜市、新泰市及下属区县，面积8 173 km2，属温带季风大陆性气候，平均气温在11～13 ℃之间。多年平均降水量为765.3 mm，多年平均蒸发量为1 664.2～1 927.0 mm。研究区北依泰山，南至蒙山，主要地貌为山地、丘陵、盆地、平原。地势四周高，中间低，成盆地状，盆地整体向东南方倾斜。区域内多水系，主要为汇河、康王河、大汶河及附属支流等，最后由东南注入东平湖[9]，如图1所示。

图1 交通位置及取样分布图Fig.1 Traffic location and sampling distribution

1.2 水文地质条件

研究区内地下水类型主要为第四系孔隙水、碳酸岩类岩溶裂隙水、碳酸岩类岩溶水、岩浆岩与变质岩类裂隙水。研究区四周被高山丘陵环绕，尤其在北、东、南部地区，有大范围泰山群火成岩和变质岩出露，风化程度弱，裂隙发育一般，富水性差，地下水类型为岩浆岩与变质岩类裂隙水。地下水来源主要是大气降雨，排泄方式主要是蒸发和向岩溶水、孔隙水径流。受地势影响，地下水由高向低径流，但岩层渗透性差，运移能力弱。在研究区的西部和中间深部为寒武—奥陶系地层，岩性以灰岩为主，其中掺杂白云岩、页岩、砂岩等。该区域岩溶发育，富水性强，地下水类型为碳酸盐岩类岩溶水和岩溶裂隙水。地下水补给主要来自降雨入渗、河流渗漏、孔隙水入渗，排泄方式主要是蒸发、向西南径流、越流补给孔隙水和人工开采。该区域发育的岩溶和断层对局部地下水流场有一定影响，但整仍受地形控制，由高往低径流。研究区中部浅层是坡积、残坡积、洪坡积形成的第四系沉积物，其成分主要是片麻岩、角闪岩、灰岩碎屑与泥沙形成的混合物，地下水类型是第四系孔隙水。在河谷两侧，含水层厚度50 m左右，富水性强，水位埋深在1～7 m，年变幅为1～3 m。随着距河流越远，水位埋深也逐渐增大。孔隙水主要受大气降雨、河流渗漏、岩溶越流补给。地下水流运移主要受地形控制，整体呈现由地势高向地势低洼的河谷地区径流，最后排泄到大汶河和其支流的下游[5]。

2 样品采集与测试

为揭示泰莱盆地地下水类型及成因，避开了受降雨影响较大的雨季，选取《山东省地下水水资源地调查与评价——泰莱盆地》项目中枯水期有效数据。如图1所示，取样点均匀分布于泰莱盆地，其中地下水点95 个（G01-G95），河水点5个（H1-H5）。井水在采样前，需先抽水数分钟，得到含水层中最新的地下水。在取样过程中，样品瓶先用去离子水清洗三次，再用样品原液润洗3次，保证所取水样能代表含水层中水体的真实水化学特征。对本次所取水样进行测试项目包括地下水常规离子（Na+、K+、Ca2+、Mg2+、HCO3-、CO32-、Cl-、SO42-）和pH 值、温度、电导率，其中pH 值、温度、电导率、HCO3-和CO32-为野外现场测试；Na+、K+、Ca2+、Mg2+测试方法为电感耦合等离子体原子发射光谱法(ICP AES-6000)；Cl-、SO42-测试方法为离子色谱法(ICS-900)。

3 聚类分析

研究区地质条件复杂，地下水类型行丰富，传统的分组方法会受主观干扰，此次研究利用Q 型聚类分析法研究泰莱盆地地区地下水水化学参数[10-13]。首先，在聚类分析之前需要对样本数据（Na+、K+、Ca2+、Mg2+、HCO3

-、CO32-、Cl-、SO4

2-、pH 值）进行“对称性”检验，经过检验发现样本除pH值以外均有较大偏度。为满足数据接近对称性分布，利用零均值法对样本数据进行了数学变换[4,15]。最后对100个样本的9个数据参数进行聚类分析，得到分析结果如图2所示。为了直观地对样品进行分组，方便不同组之间的分析，此次研究将现象线设定为9和21，当现象线为9时，取样点被分为ABC三大组，当现象线为21 时，取样点被分为10 小组（A1、AR、A3、B1、B2、B3、B4、C1、C2、C3）。ABC 三组的位置分布和盆地岩性分布如图3所示，A 组除河水AR 外，其余分布在泰山群火成岩地区，该地区地下水以基岩裂隙水为主；B组主要分布在碳酸盐出露及下伏地区，地下水类型以碳酸盐岩溶水为主；C组分布于第四系较厚的河谷地区，地下水类型以第四系孔隙水为主。由此可知聚类具有一定的分组效果，可用于进一步分析研究。

图2 聚类结果树状图及化学成分Stiff图Fig.2 Clustering result dendrogram and chemical composition Stiff diagram

图3 地质岩性及分组位置图Fig.3 Geological lithology and component distribution map of grouping location

4 结果与讨论

4.1 水化学类型

Stiff图是一种可以将地下水中主要离子组分的含量进行直观表示的水化学成分图示法[16]。对每组离子平均值绘制stiff图，用于区分不同组分之间的水化学类型，如图2(b)所示。在A 组中A1 和A2 组分相似，相同的是阳离子运移Ca2+为主要组分，不同的是A1 阴离子以HCO3-为主，A2 以SO42-为主。与之相比AR 各离子含量明显高于A1、A2，阳离子以Na++K+为主，阴离子以SO42-和Cl-为主。地下水演化过程中，地下水由中主要阴离子HCO3-转变为SO42-和Cl-[17]，表明SO4—Ca 类型的A2的演化过程晚于HCO3—Ca 类型的A1。SO4·Cl—Na 类型AR 为河水，蒸发作用和溶滤作用远远大于A1 和A2。结合图3，A1、A2 全部位于裂隙水分布区域，溶滤作用和蒸发作用较弱。B 组中B1、B3、B4 组分十分相似，阳离子和阴离子均以为主，为HCO3—Ca 类型水，只是矿化度略有不同。与之相比B2 略有差异，阴离子以和为主，为HCO3·SO4—Ca 类型水。B 组内部差异不明显。结合图3，B 组通常位于岩溶水或下伏岩溶水分布区，水化学类型主要受溶滤作用控制，受蒸发影响较弱。C 组中C1、C2、C3 间组分差异大，其主要体现在阴离子含量上，C1、C2、C3 的阴离子主要成分分别为：Cl-和阳离子主要成分均为Ca2+，其水化学类型分别为Cl·SO4—Ca，HCO3·SO4—Ca，HCO3—Ca。由图3 可知C 组主要分布在第四系较厚地区，地下水来源复杂，受蒸发作用和溶滤作用强烈，导致了复杂的水化学类型。

4.2 控制因素分析

研究区内的地下水主要来自于大气降水，但在长期演化过程中受多重水化学作用的影响。Gibbs 图可以有效反映出蒸发控制、岩石控制、降雨控制的作用程度[18]。如图4所示，A1、A2 位于岩石控制区间内，水化学成分主要受溶滤作用控制，AR 介于岩石控制和蒸发控制之间，这主要是由于河水蒸发更加强烈所致；如图5所示，B 组均位于岩石控制区域，B1～B4之间区分并明显，这由于B组主要以岩溶水和裂隙岩溶水为主，埋深大，受蒸发和降雨作用弱；相比B 组，C 组具有偏向蒸发控制的趋势（图6），这是由于C 组以第四系潜水为主，相比深层岩溶水，受蒸发作用影响更为强烈。通过C2 位置可以得出，C2 具有远高于C1 和C3 的矿化度，但并不是由蒸发造成的，而是由蒸发或溶滤沉积岩共同导致。本次研究未了避免降雨干扰，选取了研究区枯水期水样，所以三组样品受降雨控制均较弱。

图4 A区水化学Gibbs图Fig.4 Gibbs diagram of water chemistry in area A

图5 B区水化学Gibbs图Fig.5 Gibbs diagram of water chemistry in area B

图6 C区水化学Gibbs图Fig.6 Gibbs diagram of water chemistry in area C

通过Gibbs 图可知，枯水期研究区地下水主要受岩石控制，进一步分析控制岩石的类型可利用Ca2+/Na+、Mg2+/Na+与离子比值端元图[19]。由图7(a)、图7(b)可知，A1A2主要受硅酸盐控制，AR 受蒸发岩控制，A1A2 主要位于泰山群火成岩地区，硅酸盐成分占主导，AR 为河水，蒸发作用强烈，受蒸发岩控制；由图7(c)、图7(d)可知，相比A 组，B 组地下水受碳酸盐控制更多，B1～B3 受碳酸盐岩控制更加明显，这是由于B组地下水位于碳酸岩岩溶发育区；以第四系潜水为主的C组，并未表现出强烈的蒸发岩控制，而是更倾向于硅酸盐控制。

图7 离子比值端元图Fig.7 Endmember diagram of ion ratio

4.3 基于质量平衡模拟研究水岩相互作用

研究区内地下水来源是大气降水，其水化学类型主要受岩石控制，该地区主要发生溶解作用的矿物有方解石、白云石、芒硝、岩盐和石膏。为进一步确定矿物成分对地下水成因的影响，此次研究选取研究区雨水作为端源，在经过矿物溶解作用后形成现有地下水化学组分，在此基础上建立质量平衡模型[20,21]，定量分析水化学成因，如表1所示。陈晨等在研究中发现研究区农业发达，化肥和农药用量较高，城镇化建设，水泥路铺设都会导致地下水中Cl-含量明显提高[6]。结合此次研究，尤其是在B、C组分布区域，增设Cl-其他来源。同时自然条件下，石膏和芒硝具有共生和相同形成条件，我们默认A、B、C 组地区石膏和芒硝的比值相同，比值由A 组计算结果近似得出。

表1 A、B、C三组质量平衡计算 mmol/LTab.1 Mass balance calculation of A，B，C three groups

经计算A、B、C 组的平衡质量模型，得到由雨水演变成现有地下水的过程中，各物质的溶解量见表2。根据计算结果，A 组矿物溶解矿物较B、C 组少，A1 和A2 相比，在形成过程中A1 溶解了更多的方解石，吸收了更多的CO2，但是A2却溶解了更多的岩盐、石膏和芒硝，其溶解量将近A1 的一倍。由于A 组所处的位置均为泰山群组火成岩地区，不应存在自生沉积岩，固推测岩盐、石膏、芒硝等矿物来自所在区域的风积物，而A2溶解的岩盐、石膏、芒硝更多，是因为A2所在的区域海拔较低，靠近河谷，风积物较多。B 组相比A 组白云石、方解石的溶解量和CO2的吸收量明显增高，与离子端元分析结果一直，其主要原因是B组位于碳酸盐岩溶区，碳酸岩溶解量较大，白云石最大达到1.42 mmol/L。B1～B4 之间各部分溶解的量各有不同，其中B2石膏、芒硝和其他来源Cl-溶解量偏高。C组受蒸发作用影响大，但其水化学成分仍来自于矿物溶解。相比B 组，C 组岩盐、石膏、芒硝以及其他来源Cl-溶解量明显增高，C2 的白云石、岩盐、石膏、芒硝、其他来源Cl-的溶解量均为最大值，但CO2的吸收量是零，这表第四系潜水对方解石的溶解减弱，对易溶盐溶解量较高。C组水具有老水的特征，在枯水期第四系潜水更多来自下伏含水层的越流补给。

表2 矿物质溶解计算结果 mmol/LTab.2 Calculation results of mineral dissolution

5 结论

（1）利用聚类分析对泰莱盆地枯水期地下水化学成分进行分类，可排除主观干扰得到A、B、C 三组，其分别对应泰山群火成岩裂隙水、碳酸盐岩溶水、第四系潜水。在此基础可细分成10 小组，便于样品数量多、地下水类型复杂地区水化学特征分类及成因分析研究。

（2）通过Gibbs图分析枯水期ABC 三组地下水主要受岩性控制，其中C 组和地表水AR 也受蒸发作用影响。通过离子比值端元图分析得出A 组主要受硅酸盐控制，其中地表水受蒸发岩影响；B 组受硅酸盐和碳酸盐控制；C 组主要受硅酸盐控制。

（3）通过质量平衡模拟计算出10小组的不同矿物溶解量。该地区地下水溶滤过程中，除方解石、白云石矿物外，蒸发岩类溶解量占有一定比例。离子比值端元图中未显示出蒸发岩具有明显的控制作用，是因为蒸发岩溶解过程中石膏和芒硝占据相当大的比例，岩盐并不占据对优势。除此之外，研究区内人类活动导致地下水Cl-增加明显。