Hydrogeochemical characteristics and evolution mechanism of karst groundwater in the catchment area of the Sangu Spring
-
摘要: 三姑泉域岩溶地下水对晋东南地区居民生活及煤炭基地建设起着重要支撑作用。受煤炭资源大规模开采及人类活动加剧的影响,岩溶地下水水化学环境演变剧烈。通过分析三姑泉域的125 个岩溶地下水、14 个地表水及14个雨水样品,综合利用统计分析、变异系数分析、氢氧稳定同位素、Gibbs模型、离子相关关系、矿物饱和指数及因子分析,确定岩溶地下水的补给来源,分析三个子系统间岩溶水水化学特征、差异及形成演化机制。结果表明:岩溶水δD值为−77‰~−42‰,δ18O值为−10.6‰~−4.5‰,且沿大气降雨线分布,显示岩溶地下水以大气降水入渗补给为主。子系统Ⅰ与子系统Ⅱ和子系统Ⅲ具有明显不同的水化学特征和水文地球化学演化进程。子系统Ⅰ岩溶水为低TDS软水,73%为HCO3—Ca(Mg)型水;子系统Ⅱ、子系统Ⅲ为低TDS软水—高TDS硬水,水化学类型复杂,36%~40%为HCO3·SO4—Ca·Mg型水,24%~45%为SO4·HCO3—Ca型水。因子分析表明,区域岩溶地下水化学形成演化主要受控于水-岩作用、人类活动、地表水和局部裂隙水渗漏的影响。水-岩作用使得岩溶水化学组分以
HCO−3 、SO2−4 、Ca2+、Mg2+为主,采矿活动导致SO2−4 含量增加,农业活动导致NO−3 含量增加,污水排放和裂隙水渗漏导致Cl−和Na+含量增加。研究结果可为岩溶水水资源开发和管理提供科学依据。Abstract: The karst groundwater in the catchment area of the Sangu Spring plays an important supporting role in the residents' living and coal base construction in southeastern Shanxi Province.With the large-scale coal mining and the intensification of human activities, the hydrochemistry of the regional karst groundwater has changed drastically. In this study, 125 karst groundwater, 14 surface water and 14 rain water samples were collected and tested. The descriptive statistical analysis, coefficient of variation analysis, hydrogen and oxygen stable isotopes, Gibbs model, ion correlation, mineral saturation index and factor analysis are comprehensively used to determine the replenishment sources of the karst groundwater, the differences of hydrochemical characteristics and evolution mechanism among the subsystems. The results show that the δD and δ18O values of the karst groundwater range from −77‰ to −42‰ and from −10.6‰ to −4.5‰, respectively, and the data points fall near the local meteoric water line (LMWL), indicating that the karst groundwater are mainly recharged by infiltration from precipitation. The hydrochemical characteristics and hydrogeochemical evolution processes between subsystems Ⅰ and Ⅱ-Ⅲwere obviously different. The karst groundwater of subsystem Ⅰhas low salinity and is soft water, and 73% of the groundwater is of HCO3—Ca(Mg) type. However, the karst groundwater in subsystems Ⅱ and Ⅲ changes from low salinity and soft water to high salinity and hard water, and 36% to 40% of the water is of HCO3·SO4—Ca·Mg type and 24% to 45%, of SO4·HCO3—Ca type. Factor analyses show that the hydrochemical evolution processes of the regional karst groundwater are mainly controlled by water-rock interactions, human activities, leakage of surface water and local fissure water. Hydrochemical compositions ofHCO−3 ,SO2−4 , Ca2+ and Mg2+are mainly controlled by water-rock interactions. A higher concentration ofSO2−4 is controlled by mining activities,NO−3 is controlled by agriculture activities and Cl−and Na+are controlled by sewage discharge and leakage of fissure water. The results will be instructive for improving the management and utilization of groundwater resources for the local government. -
世界上约有1/4人口在利用岩溶地下水资源,其对人类生活和社会发展的重要性不言而喻[1-2]。由于岩溶地区具有特殊的地表、地下二元结构体系,地表水和污染物极易通过裂隙发育段、岩溶天窗、落水洞等裂隙、管道直接进入岩溶地下水中,地下水水质迅速恶化[3-5]。水质型缺水已是制约岩溶区经济发展、保障饮用水安全的重大问题之一[6]。
区域地下水水化学特征反映了地下水的演化特征、地质及水文地质背景和影响因素,是水文地质及水文地球化学研究的主要内容[7-10]。然而,地下水水化学组成的形成演化除受地质特性及水文地质条件等自然因素控制外,对人类活动的响应也是快速且敏感的[11-13]。如何识别影响和控制地下水水化学组分的主要因素、区分自然因素和人为因素的影响,是水文地球化学领域研究的热点[14-15]。
三姑泉域位于山西省东南部,是晋东煤炭基地晋城矿区分布区,岩溶水是当地居民生活和煤炭工业发展的重要水源。受采煤活动、城镇生活污水及化工企业污水排放的影响,区域岩溶地下水环境已发生不同程度的恶化。前人研究主要集中在三姑泉域水环境质量评价、地下水资源评价、煤炭开采对地下水环境的影响[16-20],缺乏区域上水化学特征及形成演化机制的系统研究。本文以此为背景,选取三姑泉域岩溶地下水为研究对象,采用统计分析、稳定同位素、Piper三线图、Gibbs模型、离子相关关系和因子分析等方法,对区内岩溶地下水补给来源、水化学特征及形成演化机制进行了系统研究。成果有助于深化理解北方岩溶水水化学特征及水文地球化学形成演化机制,预测污染物来源,为该区地下水污染防治及岩溶水资源的合理开发利用和管理提供科学依据。
1. 研究区概况
研究区位于山西省东南部边缘(图1),太行山脉的南端,总体地势是四面高山环绕,中部地势低且呈波状起伏,形成宽阔的山间盆地和丘陵。三姑泉出露于山西省晋城市泽州县河西乡孔庄村东北5 km处的丹河河谷西岸,现已被青天河水库淹没,泉域面积3214 km2。区内出露地层有寒武系、奥陶系、石炭系、二叠系、第四系。
研究区属于暖温带半湿润大陆性季风气候。年平均气温最低值出现在1956年,为9.0 ℃;最高值出现在1989年,为12.8 ℃。晋城近50年气温呈上升趋势,1986年以后,气温增温明显。近50年来,晋城市降水量总体呈下降趋势,多年平均降雨量为593.66 mm,多年平均蒸发量为1633.95 mm[21]。
区内水系分属于黄河流域和海河流域,以丹株岭、金泉山为界,北部为海河流域,南部为黄河流域,以丹河水系为主,海河流域的漳河水系仅分布在工作区的北部边缘地带。区内最大的河流为丹河,是沁河的一级支流。丹河主要支流有东仓河、许河、东大河、巴公河、东丹河和白水河。
碳酸盐岩类岩溶水是研究区最主要的地下水资源,奥陶系中统马家沟组和寒武系中统张夏组为主要含水层,奥陶系下统和寒武系上统为相对隔水层(图2)。岩溶水在泉域的北部、东部灰岩裸露、半裸露区直接接受大气降水入渗补给,在灰岩分布区的河谷和水库区接受地表水体的渗漏补给,在覆盖和埋藏区接受上覆含水层的越流补给;地下水自北东向南西运动;最后以泉、人工开采和潜流的方式排泄。根据区域构造和水文地质条件,三姑泉域可以进一步划分为三个子系统:东丹河子系统(Ⅰ)、任庄子系统(Ⅱ)、高平—晋城子系统(Ⅲ)[21]。
图 2 研究区水文地质剖面(A-A’剖面)Figure 2. Hydrogeological profile (along line A-A’ in Fig.1) of the study area2. 样品测试与研究方法
2015年5月—2016年11月,在研究区采取岩溶地下水样共125 个。其中高平—晋城子系统地下水样69 个样品,主要分布在泉域西部高平市—巴公镇—晋城市一带,以机民井点为主,采样深度为195~620 m,由北向南采样深度逐渐减小,取水层位主要为奥陶系马家沟组;任庄子系统地下水样41 个样品,主要分布于泉域中部陵川县—任庄水库—郭壁泉一带,以泉点和机民井点为主,井点采样深度为34~500 m,东北部泉域边缘的地下水位埋深相对较大,取水层位主要为奥陶系马家沟组,排泄带部分为寒武系张夏组地层;东丹河子系统15 个样品,分布于泉域东部路城镇—台北村—三姑泉一带,以泉点和机民井点为主,井点采样深度为80~295 m,北部补给区含水层主要为奥陶系马家沟组,南部排泄带含水层主要为寒武系张夏组。样品用于常量组分和微量组分测试,采样点分布图见图1。岩溶地下水样品中,选取105组进行氢氧同位素(δD、δ18O)测试,地表水14组,大气降雨14组用于当地大气降雨线分析。采样瓶选用500 mL聚乙烯棕色瓶,采样前分别用蒸馏水和水样润洗3 次,水样采用0.45 μm的滤膜过滤,一份加入纯HNO3至pH值小于2,用于阳离子分析;另外两份不做任何处理,用于氢氧稳定同位素和阴离子测试。
水样的温度、pH值、电导率(Ec)、溶解性总固体(TDS)等易变性指标采用便携式HANNA多参数水质测试仪,碱度在24 h内滴定。其他离子指标和氢氧同位素送至中国地质科学院水文地质环境地质研究所地下水与矿泉水监测中心进行测试。测试水样阴阳离子电荷平衡误差均在5%以内,氢同位素和氧同位素测定仪器的误差分别为±2‰和±0.1‰。
采用SPSS 21统计分析软件,对离子组分的极大值、极小值、平均值和变异系数等参数进行统计分析;选取岩溶水主要水化学参数进行相关性分析和主成分分析;借助PHREEQC软件对岩溶水中方解石、白云石、石膏、硬石膏等矿物饱和指数进行分析。
3. 结果分析
采集的岩溶地下水和地表水样品水化学统计参数见表1。
表 1 岩溶水和地表水水化学参数Table 1. Statistical parameters of the dissolved chemical components of karst groundwater and surface water项目 岩溶地下水 东丹河子系统(Ⅰ) 任庄子系统(Ⅱ) 晋城-高平子系统(Ⅲ) 地表水 最小值 最大值 均值 Cv 最小值 最大值 均值 Cv 最小值 最大值 均值 Cv 最小值 最大值 均值 Cv 最小值 最大值 均值 Cv 总硬度 175.10 1324.00 474.34 0.47 175.10 431.90 296.48 0.21 264.70 1315.00 414.60 0.47 246.00 1324.00 548.51 0.41 159.10 983.80 444.61 0.54 TDS 218.70 2781.00 751.21 0.60 218.70 510.20 338.81 0.23 301.70 2781.00 731.46 0.66 371.40 2762.00 852.59 0.51 218.90 1540.00 759.06 0.54 pH 7.08 8.57 7.74 0.03 7.43 7.89 7.62 0.02 7.08 8.57 7.82 0.04 7.20 8.53 7.72 0.03 7.42 8.61 7.96 0.04 Ca2+ 49.37 455.00 132.88 0.56 49.37 148.60 88.27 0.26 74.87 394.00 118.08 0.56 53.60 455.00 151.38 0.53 38.30 250.40 121.19 0.53 Mg2+ 12.59 80.10 34.49 0.40 12.59 24.98 18.45 0.19 16.90 80.10 29.05 0.35 18.02 77.90 41.21 0.30 9.19 86.95 34.49 0.59 K+ 0.15 16.56 1.33 1.34 0.40 2.74 1.12 0.51 0.35 7.42 1.09 1.05 0.15 16.56 1.51 1.45 1.55 33.54 9.78 0.87 Na+ 3.32 304.00 51.72 1.35 3.32 25.27 7.28 0.74 4.42 304.00 67.77 1.20 5.04 284.00 51.84 1.28 4.30 165.60 76.84 0.62 Cl− 0.36 224.00 25.40 1.16 2.50 27.07 9.91 0.67 0.36 118.70 22.26 1.10 5.36 224.00 30.64 1.10 6.07 239.30 74.27 0.82 SO2−4 21.86 1731.00 269.23 0.98 21.86 107.50 50.10 0.45 26.28 1731.00 253.36 1.18 53.92 1322.00 326.30 0.75 66.08 850.50 281.55 0.73 HCO−3 146.00 441.10 299.80 0.15 159.00 321.00 270.22 0.18 146.00 428.00 294.30 0.13 208.00 441.10 309.49 0.15 134.50 422.00 239.28 0.35 NO−3 0.20 209.50 19.41 1.17 3.92 85.60 16.31 1.23 3.76 209.50 21.40 1.46 0.20 98.30 18.90 0.87 1.88 124.40 27.14 1.21 F− 0.10 2.10 0.54 0.68 0.10 0.52 0.23 0.54 0.18 1.53 0.51 0.71 0.26 2.10 0.63 0.59 0.50 1.70 0.85 0.41 耗氧量 0.20 10.40 1.49 1.09 0.37 6.40 1.35 1.36 0.29 4.80 1.27 0.81 0.20 10.40 1.65 1.13 1.08 24.00 5.18 1.15 注:pH、变异系数(Cv)为无量纲;其余单位为mg/L;总硬度以CaCO3计。 3.1 主要离子组分特征
岩溶地下水pH值为7.08~8.57,整体上属于中—弱碱性水。TDS及总硬度变化范围较大,分别为218.7~2781 mg/L和175.1~1324 mg/L,均值分别为751.21 mg/L和474.34 mg/L属低TDS软水—高TDS硬水。TDS均值显著高于世界平均值100 mg/L,显示了三姑泉域岩溶地下水剧烈的水-岩作用。对于三个岩溶水子系统而言,东丹河子系统(Ⅰ)与任庄子系统(Ⅱ)和晋城—高平子系统(Ⅲ)相比,存在明显差异。子系统Ⅰ的TDS、总硬度和其他离子组分与子系统Ⅱ和Ⅲ相比明显较小,TDS为218.7~510.2 mg/L,总硬度为175.1~431.9 mg/L,为低TDS软水。子系统Ⅱ和Ⅲ的水化学统计参数分布范围及均值较为相近。
岩溶地下水特征阳离子主要为Ca2+,含量为49.37~455 mg/L,均值为132.88 mg/L;其次为Na+,含量为3.32~304 mg/L,均值为51.72 mg/L。特征阴离子主要为
HCO−3 ,含量为146~441.1 mg/L,均值为299.8 mg/L;其次为SO2−4 ,含量为21.86~1731 mg/L,均值为269.23 mg/L。子系统Ⅰ的Na+、Cl−、SO2−4 组分均值与子系统Ⅱ、Ⅲ存在明显差异,子系统Ⅱ、Ⅲ的Na+是子系统Ⅰ的9.3倍和7.1倍;Cl−分别为2.2倍和3.1倍;SO2−4 分别为5.1倍和6.5倍。其他离子组分无明显差异。变异系数(Cv)可说明地下水化学组分形成及演化的影响因素复杂程度,通常Cv≥1,为强变异性;0.1<Cv<1,为中等变异性;Cv≤0.1,为弱变异性。岩溶地下水各统计参数中K+、Na+、Cl−、
NO−3 和耗氧量为强变异性,pH值为弱变异性,其他离子组分属中等变异性。3.2 水化学类型
将岩溶地下水和地表水投影至Piper三线图[22-23],其中岩溶地下水按照子系统分别进行投影,结果见图3。
子系统Ⅰ岩溶水与Ⅱ、Ⅲ在Piper图中具有明显不同的分布格局。子系统Ⅰ岩溶水点集中分布在菱形左部,子系统Ⅱ、Ⅲ岩溶水点则相对比较分散,分布于菱形上半部。地表水与子系统Ⅱ、Ⅲ的岩溶水有相似的分布特征。
子系统Ⅰ岩溶水的水化学类型主要为HCO3—Ca(Mg)型,占比73%;其次为HCO3·SO4—Ca·Mg型,占比为27%。子系统Ⅱ岩溶水的水化学类型主要为HCO3·SO4—Ca·Mg型(36%),其次为SO4·HCO3—Na·Ca型(24%)和HCO3—Ca·Mg型(24%),部分岩溶水为SO4—Ca·Na型和SO4·HCO3—Ca·Na型。子系统Ⅲ岩溶水的水化学类型主要为SO4·HCO3—Ca型(45%),其次为HCO3·SO4—Ca·Mg型(40%),部分岩溶水为HCO3—Ca·Mg型、SO4—Ca·Na型和SO4·HCO3—Na·Ca型。
地表水的水化学类型相对比较复杂,HCO3·SO4型水占比为29%;SO4型水占比达71%,主要为SO4·HCO3—Ca·Na型、SO4·HCO3—Ca·Mg型,部分为SO4·HCO3·Cl—Ca·Na型、SO4—Ca·Mg型和SO4·Cl·HCO3—Ca·Na型,说明地表水已受到不同程度的污染。
4. 讨论
4.1 岩溶水补给来源
氢氧同位素广泛应用于地下水补给来源示踪、不同水体混合比例及地下水循环特征等研究中[22]。不同水体的氢氧同位素值分布见图4。
全球大气降雨线为y=8x+10,当地雨水线为y=8.42x+13.35。岩溶水趋势线为y=6.12x−11.59,其斜率和截距均小于全球和当地降雨线。这是由于降雨过程中的蒸发作用产生氢氧同位素的不平衡分馏,进而造成氢氧同位素的富集。研究区空气干燥,年平均蒸发量为降雨量的3倍左右,蒸发作用较为强烈,导致其斜率和截距相对较低。
区内岩溶水δD值为−77‰~−42‰,δ18O值为−10.6‰~−4.5‰;地表水的δD值为−70‰~−42‰,δ18O值为−9.1‰~−3.9‰。岩溶水点分布与大气降雨线趋势一致,说明大气降水是岩溶水的主要补给来源。岩溶水氢氧同位素特征明显,Ⅰ组岩溶水主要为浅部和表层岩溶水,直接接受大气降雨入渗和地表水渗漏补给,循环路径短,水位埋藏浅,蒸发作用较强,出现明显的氧漂移现象。Ⅱ组岩溶水主要分布在径流区,氢氧同位素较为贫乏,循环缓慢且路径较长;排泄区泉点也分布在该区,这是因为排泄区为大面积碳酸盐岩裸露区,在泉点排泄周边接受大气降水入渗的直接补给;该区分布有地表水点,与Ⅰ组不同,代表了河流下游岩溶水向地表水排泄的过程,使得地表水表现出与该区岩溶水同位素一致的特征。Ⅲ组岩溶水主要分布在深埋滞缓区,位于泉域北部及西北部,上覆厚层石炭、二叠系地层,地下水位埋藏深,径流滞缓,氢氧同位素最为贫乏。
4.2 Gibbs模型
Gibbs模型能够清楚地表示地下水中各种离子的大气输入、水-岩作用及蒸发浓缩效应等起源机制和变化的趋势过程[23-24]。该模型将天然水化学组分的控制因素分为三类:降雨控制型、水-岩作用型和蒸发浓缩作用型。
根据研究区内岩溶地下水的Gibbs图(图5),子系统Ⅰ的全部水点和子系统Ⅱ、Ⅲ的大部分岩溶水点落在Gibbs模型内部及水-岩作用控制区,说明区内岩溶地下水水化学组分主要受控于水-岩相互作用;子系统Ⅱ、Ⅲ的小部分水样点落在蒸发作用控制区,说明蒸发浓缩作用对地下水化学组分有一定的控制作用,这部分水点主要为浅层及表层岩溶水。大部分水样点的c(Cl−)/c(Cl−+
HCO−3 )和c(Na+)/c(Na++Ca2+)值小于0.5,说明总体上岩溶水的c(HCO−3) 相对于c(Cl−),c(Ca2+)相对于c(Na+)占有优势。随着c(Na+)/c(Na++Ca2+)值的增加,TDS含量无明显差异,这可能是因为存在阳离子交换作用,因为1 mmol/L的Ca2+(40 mg/L)与2 mmol/L的Na+(46 mg/L)交换,其质量浓度不会发生明显变化。4.3 阳离子交换作用
众多学者采用γ[(Ca2++Mg2+)−(
HCO−3 +SO2−4 )]和γ(Na+−Cl−)的关系分析地下水中的阳离子交换作用[25-27]。该关系图反映了方解石、白云石和石膏的溶解或沉淀产生的剩余的Ca2+和Mg2+,与岩盐的溶解或沉淀产生的剩余的Na+之间的相关关系。图6所示,岩溶地下水点分布于斜率大致为−1的直线两侧,表明了剩余Na+与Ca2+和Mg2+之间呈负相关,即γ[(Ca2++Mg2+)−(
HCO−3 +SO2−4 )]随着γ(Na+−Cl−)的增加而减小,指示了在地下水流过程中,Ca2+和Mg2+与吸附于岩石表面的Na+发生交换作用,Ca2+和Mg2+发生沉淀导致剩余量减小。4.4 矿物溶解作用
岩溶含水层中方解石、白云石、石膏的溶解方程为:
CaCO3+CO2+H2O=Ca2++2HCO−3 (1) CaMg(CO3)2+2CO2+2H2O=Ca2++Mg2++4HCO−3 (2) CaSO4=Ca2++SO2−4 (3) 方解石和白云石共同溶解方程为:
CaCO3+CaMg(CO3)2+3CO2+3H2O=2Ca2++Mg2++6HCO−3 (4) 式(1)、式(2)、式(4)生产的γ(Ca2+)/
γ( {\rm{HCO}}_3^- ) 摩尔比值为1∶2、1∶4和1∶3,可用其判别主要矿物溶解类别;γ(Ca2+−{\rm{SO}}_4^{2-} )和γ( {\rm{HCO}}_3^- ) 的关系可用来表示非石膏来源的Ca2+;γ(Ca2+−0.33{\rm{HCO}}_3^- )和γ( {\rm{SO}}_4^{2-} ) 的关系可用来识别石膏溶解所贡献的Ca2+含量(图7)。当岩溶地下水水化学组分受控于碳酸盐岩和石膏共同溶解时,地下水点的γ(
0.5{\rm{HCO}}_3^- +{\rm{SO}}_4^{2-} )和γ(Ca2+ + Mg2+)摩尔浓度比应接近1∶1关系线。由图7(a)可以看出,岩溶水点沿1∶1关系线分布,说明岩溶水中的Ca2+、Mg2+、{\rm{HCO}}_3^- 和{\rm{SO}}_4^{2-} 组分的来源主要为碳酸盐岩和石膏矿物的溶解。但大部分岩溶水点位于1∶1线的右下方,说明尚有其他矿物的溶解导致γ(0.5{\rm{HCO}}_3^- +{\rm{SO}}_4^{2-} )相对于γ(Ca2+ + Mg2+)过剩。由图7(b)可以看出,子系统Ⅰ岩溶水点大致沿1∶2关系线分布,说明子系统Ⅰ岩溶水方解石为主要溶解矿物。子系统Ⅱ和Ⅲ的岩溶水点一部分分布在1∶2和1∶4关系线之间,指示了方解石和白云石矿物的共同溶解;一部分分布于1∶4关系线的下部,其形成过程可能存在两个方面的原因,一是导致Ca2+浓度减小的过程,如碳酸盐岩矿物的沉淀或阳离子交换作用,使得Ca2+−{\rm{SO}}_4^{2-} 偏小;二是除石膏等硫酸盐矿物溶解外,还存在有其他{\rm{SO}}_4^{2-} 的来源,使得Ca2+−{\rm{SO}}_4^{2-} 偏小。图7(c)表明非碳酸盐岩来源的Ca2+与{\rm{SO}}_4^{2-} 的关系,可以看出,子系统Ⅰ基本沿1∶1关系线分布,表明硫酸盐矿物的溶解作用。子系统Ⅱ和Ⅲ的大部分岩溶水点分布于1∶1关系线的右下方,说明{\rm{SO}}_4^{2-} 有其他来源供给。图7(d)中,子系统Ⅰ岩溶水和子系统Ⅱ、Ⅲ的部分岩溶水,其Cl−与Na+大致沿1∶1关系线分布,指示了岩盐矿物的溶解作用;子系统Ⅱ、Ⅲ的其他岩溶水点多分布于1∶1关系线的上方,Na+摩尔浓度显著高于Cl−,说明对于子系统Ⅱ、Ⅲ而言,除岩盐矿物溶解外,还存在有其他来源供给了相当一部分量的Na+。4.5 矿物饱和指数
矿物饱和指数反映了矿物在地下水中的溶解/沉淀趋势。方解石、白云石、石膏、硬石膏、岩盐和萤石的饱和指数见图8。各子系统岩溶水中方解石和白云石矿物接近饱和状态,仅子系统Ⅰ中部分岩溶水白云石未达到饱和状态。其余矿物如石膏、硬石膏、岩盐和萤石均未达到饱和状态。区内大气降水pH值为5.65~6.4,TDS为20.35~46.4 mg/L[21],低TDS的弱酸性雨水在入渗过程中,与碳酸盐岩迅速反应,使得岩溶水中方解石和白云石矿物迅速达到饱和。
矿物饱和指数与TDS关系表明,区内岩溶水中的方解石、白云石的饱和指数与TDS无明显相关性;石膏、硬石膏和岩盐矿物的饱和指数与TDS具有较好的相关性,呈指数正相关;萤石与TDS具有一定的正相关关系。说明了岩溶地下水盐度的升高主要受控于石膏、硬石膏和岩盐矿物的溶解作用。
4.6 因子分析
为了解区域上岩溶地下水化学形成演化机制,对研究区内的125个样品,选取总硬度、TDS、pH、Ca2+、Mg2+、K+、Na+、Cl−、
{\rm{SO}}_4^{2-} 、{\rm{HCO}}_3^- 、{\rm{NO}}_3^- 、F−、SI-CaCO3、SI-CaMg(CO3)2、SI-CaSO4·2H2O、SI-CaSO4、SI-NaCl、SI-CaF2指标进行主成分分析,经KMO度量及Bartlett球型度检验,满足因子分析的要求,并选取特征值大于1的为主因子。本次共提取出了4个主因子,其累积贡献率达82.729%(表2),即4个主因子可对82.729%的样点数据进行解释,能够反映样本数据的基本信息。对每个主因子叙述如下:表 2 岩溶水各参数因子载荷矩阵Table 2. Rotated component matrix and extraction sums of squared loadings项目 F1 F2 F3 F4 总硬度 0.863* −0.363 −0.046 0.240 TDS 0.943* −0.052 0.103 −0.167 pH 0.197 0.867* 0.396 −0.142 Ca2+ 0.820* −0.425 0.086 0.214 Mg2+ 0.504 0.635* −0.469 0.239 K+ 0.218 −0.240 −0.212 −0.378 Na+ 0.500 0.226 0.227 0.699* Cl− 0.512 −0.484 0.468 0.642* {\rm{SO}}_4^{2-} 0.926* −0.031 −0.025 −0.205 {\rm{HCO}}_3^- 0.295 0.810* 0.008 0.311 {\rm{NO}}_3^- 0.268 −0.529 0.526* 0.054 F− 0.409 0.394 0.649* −0.083 SI−CaCO3 0.455 0.609* 0.415 0.313 SI−CaMg(CO3)2 0.386 0.738* 0.283 0.296 SI−CaSO4·2H2O 0.949* −0.044 −0.106 0.013 SI−CaSO4 0.949* −0.044 −0.106 0.013 SI−NaCl 0.512 −0.050 0.243 0.336 SI−CaF2 0.456 0.263 0.594* −0.009 特征值 8.410 2.845 2.054 1.583 贡献率/% 46.721 15.808 11.409 8.792 累积贡献率/% 46.721 62.529 73.938 82.729 注:*为主因子中具有较高载荷。 第一主因子F1的贡献率为46.721%,其中总硬度、TDS、Ca2+、
{\rm{SO}}_4^{2-} 、SI-CaSO4·2H2O、SI-CaSO4表现出较高的载荷,与F1呈正相关。F1得分高值点主要分布在泉域西北部、晋城市北义城镇以北,为岩溶水径流区及深埋滞缓区。F1表征了岩溶水较为充分的水-岩作用,主要是石膏及硬石膏矿物的溶解。此外,该区为采煤活动较为强烈的地区,除硫酸盐矿物溶解外,{\rm{SO}}_4^{2-} 的其他来源主要有两方面:一是采煤活动导致煤系地层中的含硫矿物氧化溶解而产生{\rm{SO}}_4^{2-} ;二是采煤过程中矿坑水疏排至地表水中,并渗漏段补给至岩溶水中。地表水中{\rm{SO}}_4^{2-} 含量均值为281.55 mg/L,最大为850.50 mg/L,因此地表水对{\rm{SO}}_4^{2-} 的输入贡献不容忽视。故F1因子代表了硫酸盐矿物的溶解以及采煤活动对地下水化学组分的影响。第二主因子F2的贡献率为15.808%,其中pH、Mg2+、
{\rm{HCO}}_3^- 、SI-CaCO3和SI-CaMg(CO3)2表现出较高的载荷。F2得分高值区主要分布在泉域中部,为地下水径流区及主径流带分布区。F2表征了方解石和白云石的溶解作用。故F2因子代表了碳酸盐矿物的溶解对地下水化学组分的影响。第三主因子F3的贡献率为11.409%,其中
{\rm{NO}}_3^- 、F−和SI-CaF2表现出较高的载荷。F3得分高值点主要分布在任庄水库以南丹河及白水河周边,泉域东北部及东部碳酸盐岩半裸露区,上覆第四系黄土,为地表水影响区和降雨入渗直接补给区,同时也是主要的农业活动区,使得地下水中的{\rm{NO}}_3^- 浓度较高,最高达125.3 mg/L。{\rm{NO}}_3^- 与农业活动过程中农药和化肥的施用密切相关,也与人类及动物粪便有关。区内有14%的岩溶水属高氟地下水(F−>1.0 mg/L),这主要与萤石矿物的溶解有关;地表水中F-含量介于0.5~1.7 mg/L之间,均值为0.85 mg/L,因此地表水的输入也不容忽视。故F3因子代表了人类活动、农业活动及萤石矿物的溶解对地下水化学组分{\rm{NO}}_3^- 、F−的影响。第四主因子F4的贡献率为8.792%,其中Na+和Cl−表现出较高的载荷。F4得分高值点主要分布在高平以南的丹河、白水河附近,呈点状分布,为地表水强烈影响区;泉域北部岩溶水深埋区,为岩溶水与裂隙水沟通区。在地表水影响区,城镇及化工企业排污是岩溶地下水Cl-和Na+整体偏高的主要原因;在岩溶水深埋区,裂隙水的渗漏补给是Na+组分偏高的主要原因,裂隙水中含钠矿物的溶解,使其Na+含量最高达407.4 mg/L[21],经由构造或钻孔与岩溶水发生联系。因此,区内岩溶水Na+和Cl−组分主要受控于污水排放和裂隙水渗漏补给的影响。
5. 结论
(1)区内岩溶地下水的主要补给来源为大气降水入渗补给,在地表水渗漏段岩溶水接受河流渗漏的补给。浅层岩溶水和表层岩溶水、径流区岩溶水及深部滞缓区岩溶水其氢氧同位素表现出明显不同的特征。
(2)岩溶水子系统Ⅰ和子系统Ⅱ、Ⅲ相比具有明显不同的水化学特征。子系统Ⅰ为低TDS软水,水化学类型主要为HCO3—Ca(Mg)型;而子系统Ⅱ、Ⅲ岩溶水总硬度和TDS变化范围大,为低TDS软水—高TDS硬水,水化学类型复杂,以HCO3·SO4—Ca·Mg型和SO4·HCO3—Ca型水为主,反应了子系统Ⅱ、Ⅲ具有强烈的水文地球化学过程。
(3)岩溶水子系统Ⅰ水化学组分特征主要受控于碳酸盐岩、硫酸盐矿物的溶解作用及阳离子交换作用的影响。而子系统Ⅱ、Ⅲ水化学形成机制复杂,除受到水-岩作用和阳离子交换作用的影响,还有其他水化学过程供给相对量的
{\rm{SO}}_4^{2-} 和Na+,使岩溶水中{\rm{SO}}_4^{2-} 和Na+明显偏高。岩溶地下水盐度的升高主要受控于石膏、硬石膏和岩盐矿物的溶解作用。(4)因子分析表明,区内岩溶地下水化学组成及演化主要受水-岩作用(碳酸盐岩、硫酸盐岩、萤石的溶解及含硫矿物氧化、阳离子交换作用)、人类活动(农业活动、人畜粪便、城镇和企业污水排放及采煤活动)、地表水和裂隙水渗漏等控制因素的影响。
-
图 2 研究区水文地质剖面(A-A’剖面)
Figure 2. Hydrogeological profile (along line A-A’ in Fig.1) of the study area
表 1 岩溶水和地表水水化学参数
Table 1 Statistical parameters of the dissolved chemical components of karst groundwater and surface water
项目 岩溶地下水 东丹河子系统(Ⅰ) 任庄子系统(Ⅱ) 晋城-高平子系统(Ⅲ) 地表水 最小值 最大值 均值 Cv 最小值 最大值 均值 Cv 最小值 最大值 均值 Cv 最小值 最大值 均值 Cv 最小值 最大值 均值 Cv 总硬度 175.10 1324.00 474.34 0.47 175.10 431.90 296.48 0.21 264.70 1315.00 414.60 0.47 246.00 1324.00 548.51 0.41 159.10 983.80 444.61 0.54 TDS 218.70 2781.00 751.21 0.60 218.70 510.20 338.81 0.23 301.70 2781.00 731.46 0.66 371.40 2762.00 852.59 0.51 218.90 1540.00 759.06 0.54 pH 7.08 8.57 7.74 0.03 7.43 7.89 7.62 0.02 7.08 8.57 7.82 0.04 7.20 8.53 7.72 0.03 7.42 8.61 7.96 0.04 Ca2+ 49.37 455.00 132.88 0.56 49.37 148.60 88.27 0.26 74.87 394.00 118.08 0.56 53.60 455.00 151.38 0.53 38.30 250.40 121.19 0.53 Mg2+ 12.59 80.10 34.49 0.40 12.59 24.98 18.45 0.19 16.90 80.10 29.05 0.35 18.02 77.90 41.21 0.30 9.19 86.95 34.49 0.59 K+ 0.15 16.56 1.33 1.34 0.40 2.74 1.12 0.51 0.35 7.42 1.09 1.05 0.15 16.56 1.51 1.45 1.55 33.54 9.78 0.87 Na+ 3.32 304.00 51.72 1.35 3.32 25.27 7.28 0.74 4.42 304.00 67.77 1.20 5.04 284.00 51.84 1.28 4.30 165.60 76.84 0.62 Cl− 0.36 224.00 25.40 1.16 2.50 27.07 9.91 0.67 0.36 118.70 22.26 1.10 5.36 224.00 30.64 1.10 6.07 239.30 74.27 0.82 {\rm{SO}}_4^{2-} 21.86 1731.00 269.23 0.98 21.86 107.50 50.10 0.45 26.28 1731.00 253.36 1.18 53.92 1322.00 326.30 0.75 66.08 850.50 281.55 0.73 {\rm{HCO}}_3^- 146.00 441.10 299.80 0.15 159.00 321.00 270.22 0.18 146.00 428.00 294.30 0.13 208.00 441.10 309.49 0.15 134.50 422.00 239.28 0.35 {\rm{NO}}_3^- 0.20 209.50 19.41 1.17 3.92 85.60 16.31 1.23 3.76 209.50 21.40 1.46 0.20 98.30 18.90 0.87 1.88 124.40 27.14 1.21 F− 0.10 2.10 0.54 0.68 0.10 0.52 0.23 0.54 0.18 1.53 0.51 0.71 0.26 2.10 0.63 0.59 0.50 1.70 0.85 0.41 耗氧量 0.20 10.40 1.49 1.09 0.37 6.40 1.35 1.36 0.29 4.80 1.27 0.81 0.20 10.40 1.65 1.13 1.08 24.00 5.18 1.15 注:pH、变异系数(Cv)为无量纲;其余单位为mg/L;总硬度以CaCO3计。 表 2 岩溶水各参数因子载荷矩阵
Table 2 Rotated component matrix and extraction sums of squared loadings
项目 F1 F2 F3 F4 总硬度 0.863* −0.363 −0.046 0.240 TDS 0.943* −0.052 0.103 −0.167 pH 0.197 0.867* 0.396 −0.142 Ca2+ 0.820* −0.425 0.086 0.214 Mg2+ 0.504 0.635* −0.469 0.239 K+ 0.218 −0.240 −0.212 −0.378 Na+ 0.500 0.226 0.227 0.699* Cl− 0.512 −0.484 0.468 0.642* {\rm{SO}}_4^{2-} 0.926* −0.031 −0.025 −0.205 {\rm{HCO}}_3^- 0.295 0.810* 0.008 0.311 {\rm{NO}}_3^- 0.268 −0.529 0.526* 0.054 F− 0.409 0.394 0.649* −0.083 SI−CaCO3 0.455 0.609* 0.415 0.313 SI−CaMg(CO3)2 0.386 0.738* 0.283 0.296 SI−CaSO4·2H2O 0.949* −0.044 −0.106 0.013 SI−CaSO4 0.949* −0.044 −0.106 0.013 SI−NaCl 0.512 −0.050 0.243 0.336 SI−CaF2 0.456 0.263 0.594* −0.009 特征值 8.410 2.845 2.054 1.583 贡献率/% 46.721 15.808 11.409 8.792 累积贡献率/% 46.721 62.529 73.938 82.729 注:*为主因子中具有较高载荷。 -
[1] KOLLARITS S, KUSCHNIG G, VESELIC M, et al. Decision-support systems for groundwater protection: innovative tools for resource management[J]. Environmental Geology,2006,49(6):840 − 848. DOI: 10.1007/s00254-006-0179-3
[2] GHASEMIZADEH R, HELLWEGER F, BUTSCHER C, et al. Review: Groundwater flow and transport modeling of karst aquifers, with particular reference to the North Coast Limestone aquifer system of Puerto Rico[J]. Hydrogeology Journal,2012,20(8):1441 − 1461. DOI: 10.1007/s10040-012-0897-4
[3] 何守阳, 朱立军, 董志芬, 等. 典型岩溶地下水系统地球化学敏感性研究[J]. 环境科学,2010,31(5):1176 − 1182. [HE Shouyang, ZHU Lijun, DONG Zhifang, et al. Study on geochemical susceptivity of groundwater system in representative karstic regions[J]. Environmental Science,2010,31(5):1176 − 1182. (in Chinese with English abstract) [4] 刘再华, 袁道先. 中国典型表层岩溶系统的地球化学动态特征及其环境意义[J]. 地质论评,2000,46(3):324 − 327. [LIU Zaihua, YUAN Daoxian. Features of geochemical variations in typical epikarst systems of China and their environmental significance[J]. Geological Review,2000,46(3):324 − 327. (in Chinese with English abstract) DOI: 10.3321/j.issn:0371-5736.2000.03.014 [5] 徐尚全, 杨平恒, 殷建军, 等. 重庆雪玉洞岩溶地下河地球化学敏感性研究[J]. 环境科学,2013,34(1):77 − 83. [XU Shangquan, YANG Pingheng, YIN Jianjun, et al. Research on the sensitivity of geochemical of underground river in Chongqing Xueyu cave[J]. Environmental Science,2013,34(1):77 − 83. (in Chinese with English abstract) [6] 张超, 张保祥, 张吉圣, 等. 肥城市岩溶水水化学特征及形成机制[J]. 中国岩溶,2018,37(5):698 − 707. [ZHANG Chao, ZHANG Baoxiang, ZHANG Jisheng, et al. Analysis of hydrochemical characteristics and formation mechanism of karst water in Feicheng City[J]. Carsologica Sinica,2018,37(5):698 − 707. (in Chinese with English abstract) [7] 霍俊伊, 于奭, 张清华, 等. 湘西峒河流域水化学特征及无机碳通量计算[J]. 水文地质工程地质,2019,46(4):64 − 72. [HUO Junyi, YU Shi, ZHANG Qinghua, et al. Hydrochemical characteristics and estimation of the Dissolved inorganic carbon flux in the Donghe River Basin of western Hunan[J]. Hydrogeology & Engineering Geology,2019,46(4):64 − 72. (in Chinese with English abstract) [8] 殷晓曦, 陈陆望, 谢文苹, 等. 采动影响下矿区地下水主要水-岩作用与水化学演化规律[J]. 水文地质工程地质,2017,44(5):33 − 39. [YIN Xiaoxi, CHEN Luwang, XIE Wenping, et al. Main water-rock interactions and hydrochemical evolution in the aquifers under the mining-induced disturbance in a mining district[J]. Hydrogeology & Engineering Geology,2017,44(5):33 − 39. (in Chinese with English abstract) [9] 黄奇波, 覃小群, 程瑞瑞, 等. 左江中游岩溶峰林区河流交互带水化学特征与控制因素[J]. 水文地质工程地质,2019,46(5):1 − 8. [HUANG Qibo, QIN Xiaoqun, CHENG Ruirui, et al. Hydrochenmical characteristics and control factors of karst hyporheic zones in the karst peak forest region of the middle reaches of the Zuo River[J]. Hydrogeology & Engineering Geology,2019,46(5):1 − 8. (in Chinese with English abstract) [10] 刘英学, 梁韵, 乔梁. 邢台百泉岩溶地下水系统水-岩反应特征分析[J]. 南水北调与水利科技,2009,7(4):63 − 66. [LIU Yingxue, LIANG Yun, QIAO Liang. Analysis on water-rock reaction characteristics of groundwater system in Xingtai-Baiquan springs[J]. South-to-North Water Transfers and Water Science & Technology,2009,7(4):63 − 66. (in Chinese with English abstract) DOI: 10.3969/j.issn.1672-1683.2009.04.018 [11] 邓启军, 李方红, 李伟, 等. 蒲阳河流域地下水水化学及同位素特征[J]. 水文地质工程地质,2017,44(2):8 − 14. [DENG Qijun, LI Fanghong, LI Wei, et al. The chemical and isotopes characteristic of groundwater in Puyang river basin[J]. Hydrogeology & Engineering Geology,2017,44(2):8 − 14. (in Chinese with English abstract) [12] 张军以, 王腊春, 苏维词, 等. 岩溶地区人类活动的水文效应研究现状及展望[J]. 地理科学进展,2014,33(8):1125 − 1135. [ZHANG Junyi, WANG Lachun, SU Weici, et al. Status and prospect of the hydrological effects of human activities in the Karst area[J]. Progress in Geography,2014,33(8):1125 − 1135. (in Chinese with English abstract) DOI: 10.11820/dlkxjz.2014.08.013 [13] 刘绍华, 郭芳, 姜光辉, 等. 桂林市峰林平原区岩溶水文地球化学特征[J]. 地球与环境,2015,43(1):55 − 65. [LIU Shaohua, GUO Fang, JIANG Guanghui, et al. Hydrogeochemical characteristics of peak forest plain in Guilin City, China[J]. Earth and Environment,2015,43(1):55 − 65. (in Chinese with English abstract) [14] 黄科云, 刘德深, 马祖陆, 等. 云南鹤庆西山岩溶地下水主要离子雨季和旱季对比及来源分析[J]. 地球与环境,2015,43(2):183 − 189. [HUANG Keyun, LIU Deshen, MA Zulu, et al. Major ion chemistry and their sources of karstic ground water from the Heqing west mountain, China during flood and dry seasons[J]. Earth and Environment,2015,43(2):183 − 189. (in Chinese with English abstract) [15] 蔡月梅, 蔡五田, 刘金巍, 等. 河南省某大型水源地岩溶水水化学及同位素特征[J]. 水文地质工程地质,2018,45(5):41 − 47. [CAI Yuemei, CAI Wutian, LIU Jinwei, et al. Chemical and isotopic characteristics of the karst groundwater in the wellfield in Henan Province[J]. Hydrogeology & Engineering Geology,2018,45(5):41 − 47. (in Chinese with English abstract) [16] 李慧清. 晋城市三姑泉域地下水生态现状分析及修复对策[J]. 山西水利科技,2018(2):61 − 64. [LI Huiqing. Analysis on the groundwater ecological status in the Sanguquan spring-feeding area of Jincheng City and restoration countermeasures[J]. Shanxi Hydrotechnics,2018(2):61 − 64. (in Chinese with English abstract) DOI: 10.3969/j.issn.1006-8139.2018.02.020 [17] 王振兴, 李向全, 侯新伟, 等. 煤炭开采条件下三姑泉域岩溶含水层保护评价[J]. 中国岩溶,2019,38(1):28 − 39. [WANG Zhenxing, LI Xiangquan, HOU Xinwei, et al. Karst aquifer protection evaluation to the Sangu spring basin under the condition of coal mining[J]. Carsologica Sinica,2019,38(1):28 − 39. (in Chinese with English abstract) [18] 刘晓红, 张永波. 山西晋城三姑泉域地下水资源评价[J]. 太原科技大学学报,2009,30(3):261 − 263. [LIU Xiaohong, ZHANG Yongbo. Groundwater resources evaluation on Sangu spring region of Jincheng of Shanxi Province[J]. Journal of Taiyuan University of Science and Technology,2009,30(3):261 − 263. (in Chinese with English abstract) DOI: 10.3969/j.issn.1673-2057.2009.03.019 [19] 刘晓红, 张永波. 山西晋城三姑泉域岩溶地下水水质评价[J]. 山西建筑,2009,35(5):173 − 174. [LIU Xiaohong, ZHANG Yongbo. On the karst underwater water quality assessment of Sanguquan in Jincheng of Shanxi[J]. Shanxi Architecture,2009,35(5):173 − 174. (in Chinese with English abstract) DOI: 10.3969/j.issn.1009-6825.2009.05.108 [20] 许志峰, 张志祥, 刘晓霞. 三姑泉域水环境质量评价及水污染控制对策[J]. 地下水,2012,34(4):87 − 90. [XU Zhifeng, ZHANG Zhixiang, LIU Xiaoxia. Water environmental quality assessment and strategy for water pollution control in Sangu spring basin[J]. Ground Water,2012,34(4):87 − 90. (in Chinese with English abstract) [21] 侯新伟, 李向全, 王振兴, 等. 晋东能源基地水文地质环境地质调查报告[R]. 石家庄: 中国地质科学院水文地质环境地质研究所, 2016: 45-189. HOU Xinwei, LI Xiangquan, WANG Zhenxing, et al. Hydrogeological and environment geological survey report of Jindong energy base[R]. Shijiazhuang: Institute of Hydrogeology and Environmental Geology, 2016: 45-189. (in Chinese)
[22] 黄平华, 陈建生, 宁超, 等. 焦作矿区地下水水化学特征及其地球化学模拟[J]. 现代地质,2010,24(2):369 − 376. [HUANG Pinghua, CHEN Jiansheng, NING Chao, et al. Hydrochemical characteristics and hydrogeochemical modeling of groundwater in the Jiaozuo mining district[J]. Geoscience,2010,24(2):369 − 376. (in Chinese with English abstract) DOI: 10.3969/j.issn.1000-8527.2010.02.023 [23] PIPER A M. A graphic procedure in the geochemical interpretation of water-analyses[J]. Eos, Transactions American Geophysical Union,1944,25(6):914 − 928. DOI: 10.1029/TR025i006p00914
[24] 韦虹, 吴锦奎, 沈永平, 等. 额尔齐斯河源区融雪期积雪与河流的水化学特征[J]. 环境科学,2016,37(4):1345 − 1352. [WEI Hong, WU Jinkui, SHEN Yongping, et al. Hydrochemical characteristics of snow meltwater and river water during snow-melting period in the headwaters of the Ertis River, Xinjiang[J]. Environmental Science,2016,37(4):1345 − 1352. (in Chinese with English abstract) [25] CAROL E, MAS-PLA J, KRUSE E. Interaction between continental and estuarine waters in the wetlands of the northern coastal plain of Samborombón Bay, Argentina[J]. Applied Geochemistry,2013,34:152 − 163. DOI: 10.1016/j.apgeochem.2013.03.006
[26] HUANG G X, SUN J C, ZHANG Y, et al. Impact of anthropogenic and natural processes on the evolution of groundwater chemistry in a rapidly urbanized coastal area, South China[J]. Science of the Total Environment,2013,463/464:209 − 221. DOI: 10.1016/j.scitotenv.2013.05.078
[27] LIU F, SONG X F, YANG L H, et al. The role of anthropogenic and natural factors in shaping the geochemical evolution of groundwater in the Subei Lake basin, Ordos energy base, Northwestern China[J]. Science of the Total Environment,2015,538:327 − 340. DOI: 10.1016/j.scitotenv.2015.08.057
-
期刊类型引用(38)
1. 刘海,黄健敏,宋阳,赵国红,王旭东,魏伟. 太湖县花凉亭水库流域地下水水化学特征及水质评价. 环境科学. 2025(02): 809-820 . 百度学术
2. 杨永强,李帆,杜倩倩,韩志强,韩自强. 山西朔州源子河流域地下水水化学特征及其影响因素. 地球与环境. 2025(01): 99-111 . 百度学术
3. 金爱芳,殷秀兰,李长青,李文娟,庞菊梅,金晓媚. 张家口地区枯水期地下水水化学特征及其成因机制分析. 环境科学. 2024(02): 826-836 . 百度学术
4. 郑慧铭,张丰,史启朋,仝路,孟甲. 煤炭资源型城市浅层地下水化学特征及成因分析. 人民长江. 2024(02): 51-56 . 百度学术
5. 刘海,宋阳,李迎春,魏伟,赵国红,王旭东,黄健敏. 长江流域安庆段浅层地下水水化学特征及控制因素. 环境科学. 2024(03): 1525-1538 . 百度学术
6. 吕晓立,郑跃军,刘可,李春燕,赵伟,韩占涛. 兰州不同城镇功能区地下水氟赋存特征及影响因素. 水文地质工程地质. 2024(02): 215-226 . 本站查看
7. 陶兰初,邹祖建,涂春霖,和成忠,马一奇,杨坤,刘振南,张七道. 滇东富源地区岩溶地下水水化学特征及控制因素. 地球与环境. 2024(02): 253-265 . 百度学术
8. 黄勇,胡稳,朱世超,程冠全,梁菊. 高原盆地区域地下水无机物污染、源解析及健康风险评价. 环境科学研究. 2024(06): 1389-1400 . 百度学术
9. 潘登,王帅,郭红,王琳,李屹田. 黄河下游(河南段)潜水水文地球化学特征及补给来源识别. 安全与环境工程. 2024(04): 170-180 . 百度学术
10. 周智强,黄奇波,汪玉松,罗飞,梁建宏,熊江俣. 典型岩溶矿区地表水和地下水补给来源及水化学演化机制. 环境科学. 2024(09): 5264-5276 . 百度学术
11. 张春来,杨慧,黄芬,邱成,朱同彬. 亚热带季风区边缘坡立谷河流水化学特征及外源水的增汇效应:以广西马山县清波河流域为例. 地学前缘. 2024(05): 377-386 . 百度学术
12. 曾蛟,李海波,徐华源,褚宏宪,赵正鹏,庞国涛,王丽鑫,孔令号. 山东莱州-平度一带水化学特征成因分析及水质评价. 海洋地质前沿. 2024(11): 15-24 . 百度学术
13. 陈吉吉,陶蕾,刘保献,杨庆,席玥,徐蘇士,陈圆圆,荆红卫. 北京市平原区地下水铁锰分布特征及成因分析. 水文地质工程地质. 2024(06): 198-207 . 本站查看
14. 纪明宇,贾德彬,郝玉胜,刘婷,杨丽娜,李晓艳,吕晨歌,尚紫琴. 闪电河流域正蓝旗段水化学和同位素特征及水体转换关系. 应用生态学报. 2024(11): 3149-3156 . 百度学术
15. 王继华,于福荣,田鹏州,潘登,李涛,孙淼. 南太行山前冲积平原潜水水化学特征及控制因素分析. 干旱区资源与环境. 2023(05): 130-138 . 百度学术
16. 贾洋,崔素丽. Cu~(2+)对重塑黄土饱和渗透系数的影响研究. 水文地质工程地质. 2023(03): 93-103 . 本站查看
17. 王震威,周施阳,李振,杨帆,万军. 天台盆地浅层地下水水化学特征及成因分析. 资源环境与工程. 2023(04): 438-445 . 百度学术
18. 魏世博,韩金辉,毕俊杰,张福存,徐蓉桢. 保定市一亩泉漏斗区地下水水化学特征演变及驱动机制解析. 地质论评. 2023(05): 1845-1858 . 百度学术
19. 陶小虎,王赫生,龚建师,许乃政,檀梦皎,胡小雨,戴飞虎,孟显超,朱春芳,周锴锷,李亮. 苏北阜宁县地下水水化学特征及成因分析. 华东地质. 2023(03): 357-366 . 百度学术
20. 柳浩然,张文强,刘文,马雪莹,关琴,张海林. 基于长期监测的济南趵突泉水文地球化学演化规律研究. 中国岩溶. 2023(05): 1061-1073 . 百度学术
21. 张文强,滕跃,唐飞,王金晓,许庆宇,张海林. 山东省肥城断块岩溶水系统地下水水化学特征及演化分析. 中国岩溶. 2023(05): 1047-1060+1084 . 百度学术
22. 黄静,肖先煊,许模,李兆峰. 车马碧引水隧洞近场区水文地球化学特征及成因浅析. 四川地质学报. 2023(04): 638-647 . 百度学术
23. 赵宝峰,吕玉广. 新上海一号井田水文地球化学特征及控制因素. 干旱区资源与环境. 2022(03): 92-98 . 百度学术
24. 岑雷,杨苏诗,曹飞龙,沈龙梅,廖柱坤. 粤港澳大湾区滨海区域地下水质量特征及成因浅析——以江门市新会区为例. 地下水. 2022(01): 64-68 . 百度学术
25. 李琪,赵小二,武周虎,武桂芝,张成. 溶潭容积对岩溶管道穿透曲线的影响实验. 吉林大学学报(地球科学版). 2022(03): 979-991 . 百度学术
26. 向超,王运生,杜春阳,赵方彬,苏毅,潘元贵. 川藏铁路察雅段低温温泉钙华成因及其工程意义. 吉林大学学报(地球科学版). 2022(04): 1244-1254 . 百度学术
27. 倪泽华. 新会区地下水水化学特征及成因. 陕西水利. 2022(07): 124-128 . 百度学术
28. 吴涛. 煤层深埋矿区地下水化学组分含量随含水层埋深变化特征与指示意义——以榆横矿区大海则井田为例. 勘察科学技术. 2022(04): 33-38 . 百度学术
29. 李玉山,李惠平,王虎,周彦兆,何泽林. 河西堡化工园区地下水化学特征与高氟水成因机制. 干旱区资源与环境. 2022(12): 119-126 . 百度学术
30. 杨海娇,魏加华,任倩慧. 柴达木盆地典型流域地表水-地下水转化关系及水化学特征. 干旱区研究. 2022(05): 1543-1554 . 百度学术
31. 查学芳,吴攀,李学先,陈世万,黄家琰,李清光,陈思睿. 基于水化学与硫同位素的卡林型金矿区岩溶水文地球化学特征及控制因素. 环境科学. 2022(11): 5084-5095 . 百度学术
32. 卢丽,陈余道,邹胜章,樊连杰,林永生,王喆. 岩溶区典型工业型城市地下水水化学特征及成因机制. 中国岩溶. 2022(04): 588-598 . 百度学术
33. 李晓波,李杭,杨宝萍,赵新村. 泰安市旧县水源地水化学特征及成因分析. 环境工程技术学报. 2022(06): 2002-2010 . 百度学术
34. 龚亚兵,龚绪龙,许书刚,唐鑫,苏东,吴夏懿. 苏南地区地下水化学特征及演化分析. 地质论评. 2022(06): 2207-2218 . 百度学术
35. 朱刚. 贵州省德江县高位隐蔽性地质灾害形成条件分析. 工程技术研究. 2022(19): 190-192 . 百度学术
36. 韩庆定,罗锡宜. 广东佛山市高明区李家村岩溶塌陷群成因机理分析. 中国地质灾害与防治学报. 2021(04): 56-64 . 百度学术
37. 李向全,马剑飞,张春潮,王振兴,付昌昌,白占学. 川藏铁路格聂山和察雅段构造岩溶发育规律及岩溶地下水循环模式研究. 水文地质工程地质. 2021(05): 34-45 . 本站查看
38. 马剑飞,李向全,张春潮,付昌昌,白占学,王振兴. 西藏波密冰川覆盖区大型河流与断裂带地下水转化关系. 水文地质工程地质. 2021(05): 23-33 . 本站查看
其他类型引用(7)