环境科学  2024, Vol. 45 Issue (6): 3196-3204   PDF    
安徽涡河流域水化学与同位素特征及水体转化关系
许金昭, 刘桂建, 司雯, 刘锐佳, 王冠宇     
中国科学技术大学地球和空间科学学院, 合肥 230026
摘要: 以安徽涡河流域为研究区域, 通过采集涡河地表水、浅层地下水和中深层地下水样品, 运用Gibbs图、离子比例和MixSIAR模型等方法对各类水体的水化学参数和氢氧同位素进行分析, 揭示并量化了各类水体之间的转化关系. 结果表明, 研究区地表水和地下水主要呈现中性至弱碱性, 地表水的水化学类型以Cl·SO4·HCO3-Na和Cl·SO4-Na型为主, 浅层地下水以HCO3-Ca·Mg和HCO3-Mg·Na型为主, 中深层地下水主要为Cl·HCO3-Na型. 各类水体的水化学特征受到岩石风化、蒸发浓缩及正向阳离子交换等多重因素共同影响. 地表水和地下水的δ18O和δ2H分布特征表明大气降水是该区域水体的主要补给来源, 且地下水δ18O和δ2H值与K+、Na+、Cl-、SO42-和NO3-浓度存在显著相关性. 根据MixSIAR模型分析结果, 地表水接受大气降水和浅层地下水的贡献率分别为46.5%和53.5%;浅层地下水的补给来源分别为大气降水(57.4%)和地表水(42.6%);中深层地下水的补给主要来自上游地下水的侧向径流补给.
关键词: 涡河流域      氢氧稳定同位素      水化学特征      转化关系      MixSIAR模型     
Hydrochemical and Stable Isotopic Characteristics and Transformation Relationship of Water Bodies in the Guohe River Basin, Anhui Province
XU Jin-zhao , LIU Gui-jian , SI Wen , LIU Rui-jia , WANG Guan-yu     
School of Earth and Space Sciences, University of Science and Technology of China, Hefei 230026, China
Abstract: The Guohe River Basin in Anhui Province was selected as the research area for this study. By collecting surface water, shallow groundwater, and middle-deep groundwater samples, various hydrochemical parameters and stable isotopes of water in different water bodies were analyzed using methods such as the Gibbs diagram, ion ratios, and MixSIAR model to reveal and quantify the transformation relationships between these water bodies. The results indicated that both surface water and groundwater in the study area were predominantly neutral to weakly alkaline. The hydrochemical types of surface water were mainly characterized by Cl·SO4·HCO3-Na and Cl·SO4-Na types, whereas the shallow groundwater exhibited HCO3-Ca·Mg and HCO3-Mg·Na types, and the middle-deep groundwater was of the Cl·HCO3-Na type. The hydrochemical characteristics of various water bodies were influenced by multiple factors such as rock weathering, evaporation concentration, and positive cation exchange. The distribution characteristics of δ18O and δ2H values in surface water and groundwater indicated that atmospheric precipitation was the main water source. The δ18O and δ2H in groundwater were significantly correlated with K+, Na+, Cl-, SO42-, and NO3-. According to the analysis using the MixSIAR model, the contribution of atmospheric precipitation to surface water was 46.5%, whereas the contribution from shallow groundwater was 53.5%. The sources of shallow groundwater were identified as atmospheric precipitation (57.4%) and surface water (42.6%), and the main source of supply for middle-deep groundwater was lateral flow from upstream groundwater.
Key words: Guohe River Basin      hydrogen and oxygen stable isotopes      hydrochemical characteristics      transformation relationships      MixSIAR model     

水资源是地球生物生存的关键要素, 但由于经济的快速发展、生活方式的显著改变以及工业、农业和生活用水需求的持续增长, 全球范围内的水资源危机日益加剧[1, 2]. 在此背景下, 深入研究流域水体的转化过程和相关机制, 对于流域水资源的科学管理和有效保护具有重要意义[3~5].

涡河, 位于淮北平原, 是淮河流域的第二大支流. 地表水和地下水是安徽省涡河流域农业、工业和生活用水的主要水源, 在区域生态文明建设、自然资源利用和经济发展方面发挥了重要作用. 然而, 该流域各类水体的水循环关系由于气候变化和人类活动的增加而发生了显著变化. 流域内地表水和地下水之间的转化关系是水循环过程的核心环节. 尽管此关系至关重要, 但目前的研究还存在不足, 尤其是缺乏不同水体间转化关系的定量化研究. 现有研究主要侧重于涡河流域的水质评价[6, 7]和污染物成因分析[8], 对水化学来源分析和转化关系的深入探讨较少.

近年来, 氢氧同位素技术已经成为研究流域水循环过程的有效工具[9~13]. 结合水化学特征分析, 这一技术有助于精确反映区域水文的形成和演化机制[14~16]. 作为一种天然的示踪剂, 氢氧同位素被广泛应用于水循环过程研究, 可以揭示不同水体间的转化关系, 评估蒸发等因素的影响[17~20].

因此, 本研究采用水文地球化学方法和氢氧同位素技术, 对安徽涡河流域内不同水体的水化学和同位素特征进行分析, 探讨地表水和地下水的补给来源, 揭示不同水体间的转化关系, 以期为区域水资源的开发、利用和保护提供建议和科学数据.

1 材料与方法 1.1 研究区概况

安徽涡河流域地处淮北平原, 行政区划包括安徽省蒙城县、涡阳县和谯城区. 流域总长约173 km, 流域面积约4 039 km2. 该地区位于中国南北气候过渡地带, 属于暖温带半湿润大陆性季风气候区, 夏秋季节降水集中, 容易引发洪涝灾害, 冬春季节则干旱少雨. 多年平均降水量为600~900 mm, 但年内分布不均, 以6~9月为主要降水期. 降水量年际变幅较大, 最大年降水量约为最小年降水量的4倍. 多年平均蒸发量为1 668.8 mm. 年平均径流深度为60~150 mm, 北部地区较小, 南部地区较大. 涡河流域地势平坦, 海拔范围为21.8~88.0 m, 呈现西北高、东南低的趋势[21]. 该流域地表岩性主要由棕红色、灰黄色粉质黏土、黄色粉质黏土和粉土构成. 总体地貌为典型的黄淮堆积型地貌, 包括河谷平原、河间地块和残丘. 地区内无基岩出露, 地表被厚厚的晚新生代松散层所覆盖.

根据地下水赋存条件、含水层岩性和水力特征等条件, 研究区内地下水可分为3个含水层:浅层、中层和深层孔隙水含水岩组[22]. 浅层孔隙水含水岩组处于开放状态, 与大气降水联系密切, 地下水位在这个含水岩组中存在强烈的季节性波动, 地下水年龄多小于5 × 10³ a. 浅层孔隙水含水岩组对应的浅层地下水主要赋存于第四纪晚更新世含水层中, 含水层底板埋深小于50.0 m, 厚度在10.0~35.0 m之间, 富水性较好. 中层孔隙水含水岩组处于半封闭状态, 随着深度增加, 与大气降水的联系逐渐减弱, 地下水年龄主要在10 × 10³~20 × 10³ a之间. 中层孔隙水含水岩组对应的中层地下水主要赋存于第四纪中更新世和早更新世含水层中, 含水层不连续, 底板埋深为50.0~150.0 m, 厚度为17.5~34.0 m, 富水性较差. 深层孔隙水含水岩组大多处于封闭状态, 地下水年龄超过20 × 10³ a. 深层孔隙水含水岩组对应的深层地下水主要赋存于新近纪中新世含水层中, 含水层底板埋深为150.0~480.0 m, 厚度为42.4~118.5 m.

1.2 样品采集与分析方法

2022年8~11月期间对安徽涡河流域地表水和地下水进行系统采样, 共采集涡河流域地表水样品32组、地下水样品27组, 其中浅层孔隙地下水样品15组, 中深层孔隙地下水样品12组. 采样点空间分布如图 1所示. 平均每隔5 km采集一个涡河地表水样品, 浅层孔隙地下水样品全部来源于涡河流域周边农户自用水井, 井深为10.0~50.0 m;中深层孔隙地下水取自涡河旁的乡镇集中式水源地水井和地质环境监测水井, 井深为50.0~400.0 m.

图 1 涡河流域采样点空间分布与高程叠加 Fig. 1 Spatial distribution of sampling locations in the Guohe River Basin with elevations superimposed

利用便携式多参数水质分析仪(WTW multi 3400i, Germany)对水样中的总溶解固体(TDS)和pH等参数进行现场测量. 所有水样用聚乙烯瓶收集, 取样前用待测水(地表水或地下水)润洗瓶子3次以上, 放置于便携式冰箱中4℃下保存, 并在24 h内送往中国科学技术大学实验室进行检测. 所有水样分别过0.45 μm和0.22 μm微孔滤膜进行主要离子和氢氧同位素值测试. 主要阳离子(K+、Ca2+、Na+和Mg2+)采用电感耦合等离子体质谱法(NexION 1000, PerkinElmer, USA)分析, 主要阴离子(SO42-、NO3-和Cl-)采用离子色谱法(ICS-900, Thermo Fischer Scientific, USA)测定. 通过酸碱中和滴定法测定HCO3-含量. 所有水样主要离子的电荷平衡误差均低于± 5%, 保证了数据的准确性. 利用液态水同位素分析仪(LWIA-24d, USA)测定水样δ18O和δ2H值, δ18O和δ2H的测定精度分别为± 0.1‰和± 0.4‰.

1.3 数据处理

本文采用描述性统计方法对各类水体的水化学参数进行分析, 研究其化学组成成分以及氢氧同位素的空间分布特征. 绘制Piper图阐明水体的水化学类型, 并结合Gibbs图和主要离子比值图, 以确定水中化学成分来源. 氯碱指数(CAI)和矿物饱和指数(SI)的计算公式如下[23]

(1)
(2)
(3)

式中, 各离子浓度的单位为meq·L-1, IAP为离子活度积, K为溶解平衡常数. SI指数通过PHREEQC软件(版本:3.7.3-15968)计算得出.

SPSS 26.0用于氢氧同位素和水中主要离子的相关性分析, 贝尔斯混合模型(MixSIAR)是基于R构建的带有马尔可夫链蒙特卡洛(MCMC)的贝叶斯框架, 它包含与多种来源有关的判别因子, 能够辨别各类水体的转化关系, 并更加准确地计算水源贡献率[24]. 通过将水样δ18O和δ2H的平均值和标准差代入MixSIAR, 在R语言环境下运行模型, 最终提取和分析模型结果. 所有图例均使用Arcgis 10.8和Origin 2022绘制.

2 结果与分析 2.1 水化学特征

研究区地表和地下水的理化性质参数如表 1所示. 涡河地表水、浅层地下水和中深层地下水的pH值分别为8~8.4、7.4~8.3和7.1~8.1, 整体呈中性至弱碱性. 地表水、浅层地下水和中深层地下水的ρ(TDS)平均值分别为1 077.4、771.2和532.3 mg·L-1, 水质总体为低矿化度[ρ(TDS) < 1 000 mg·L-1]. 地表水中主要阳离子浓度大小顺序为:Na+ > Ca2+ > Mg2+ > K+, 浓度平均值分别为307.3、54.5、43.8和17.1 mg·L-1;主要阴离子浓度大小顺序为:Cl- > HCO3- > SO42- > NO3-, 浓度平均值分别为270.1、243.8、219.2和1.7 mg·L-1. 浅层地下水和中深层地下水中阴阳离子浓度大小顺序均一致, 分别为:Na+ > Ca2+ > Mg2+ > K+和HCO3- > Cl- > SO42- > NO3-, 对应离子在浅层地下水中的浓度平均值分别为135.0、55.9、51.7、1.2、390.6、118.0、103.4和4.7 mg·L-1, 在中深层地下水中的浓度平均值分别为271.2、43.6、30.8、1.8、391.8、186.1、183.9和0.5 mg·L-1. 地表水和地下水中主要阳离子浓度大小顺序一致, 而主要阴离子浓度顺序存在差异, 主要是由地表水中的Cl-浓度较高导致, 研究区周边主要为农业用地, 农药和化肥中含有大量Cl-, 农田排水进入河流可能引起地表水Cl-浓度增高.

表 1 涡河流域地表水及地下水理化性质参数1) Table 1 Physical and chemical properties of surface water and groundwater in the Guohe River Basin

Piper图是研究水化学类型的重要方法, 被广泛用于分析水体水化学组成特征[25, 26]. 如图 2所示, 研究区水体中的阳离子以Na+、Ca2+和Mg2+为主, 阴离子主要为HCO3-和Cl-. 在地表水中, Na+、Ca2+和Mg2+的毫克当量分别占阳离子总量的66.22%、13.51%和18.09%, HCO3-、Cl-和SO42-的毫克当量分别占阴离子总量的24.70%、47.04%和28.23%;在浅层和中深层地下水中, Na+、Ca2+和Mg2+的毫克当量分别占阳离子总量的44.78%、21.33%和32.88%以及71.09%、13.14%和15.49%, HCO3-、Cl-和SO42-的毫克当量分别占阴离子总量的53.88%、27.98%和18.12%以及41.45%、33.83%和24.72%. 地表水的水化学类型主要为Cl·SO4·HCO3-Na和Cl·SO4-Na型, 浅层地下水以HCO3-Ca·Mg和HCO3-Mg·Na型为主, 中深层地下水主要表现为Cl·HCO3-Na型.

图 2 涡河流域水体Piper三线图 Fig. 2 Piper plot of water bodies in the Guohe River Basin

2.2 氢氧同位素特征

氢氧稳定同位素能在区域尺度上为水资源的量化分析和蒸发规律研究等提供重要的信息支持[27, 28]. 从表 1可知, 涡河地表水的δ18O和δ2H值范围分别为-8.0‰~-6.5‰和-60.5‰~-53.5‰, 平均值分别为-7.3‰和-55.7‰. 浅层地下水的δ18O和δ2H值范围分别为-10.9‰~-7.3‰和-73.8‰~-50.3‰, 平均值分别是-8.5‰和-56.5‰. 中深层地下水的δ18O和δ2H值范围分别为-10.4‰~-7.6‰和-77.1‰~-61.6‰, 平均值分别为-9.8‰和-72.8‰. 研究区各类水体中δ18O和δ2H值的大小顺序均为:涡河地表水 > 浅层地下水 > 中深层地下水, 呈现出随着水深的增加, δ18O和δ2H值逐渐递减的趋势.

3 讨论 3.1 基于水化学特征的不同水体成因解析 3.1.1 Gibbs图分析

Gibbs图能够用于识别水体水化学特性的主要控制机制, 将控制水体水化学特性的因素主要分解为降水、蒸发浓缩和岩石风化这3个自然因素[29]. 研究区大部分水样点均分布在Gibbs图的可识别区间内, 只有部分阳离子例外(图 3). 涡河流域地表水样品的Na+ /(Na+ + Ca2+)和Cl- /(Cl- + HCO3-)值分别为0.80~0.87和0.42~0.63, 表明研究区地表水中阳离子主要受蒸发浓缩的影响, 阴离子受岩石风化和蒸发浓缩的共同作用. 浅层地下水样品的Na+ /(Na+ + Ca2+)和Cl- /(Cl- + HCO3-)值分别为0.29~0.97和0.02~0.64, 说明浅层地下水中阳离子除了受岩石风化和蒸发浓缩的影响, 可能还受人类活动或阳离子交换的影响, 而阴离子主要受岩石风化的控制. 中深层地下水样品的Na+ /(Na+ + Ca2+)和Cl- /(Cl- + HCO3-)值分别为0.30~0.98和0.07~0.56, 其水化学特征的影响因素和浅层地下水的一致. 从Gibbs图中可以看出研究区地下水中阳离子变化的幅度大于阴离子. 研究区地下含有丰富的硅酸盐矿物, 以上矿物在水岩相互作用过程中会释放大量Na+, Na+又会进一步和Ca2+、Mg2+进行阳离子交换, 导致地下水阳离子波动大[30].

图 3 涡河流域水体Gibbs图 Fig. 3 Gibbs plot of water bodies in the Guohe River Basin

3.1.2 离子比例分析

毫克当量浓度比例系数能进一步确定水体水化学特征的影响因素[31]. Cl-/Na+比例可以指示水体中Na+的主要来源:如果水样点位接近等值线, 说明水体中Na+主要来源于溶解的蒸发盐岩;位于等值线的下方, 表明水体中的Na+可能来自于阳离子交换作用和硅酸盐矿物的风化[32].(SO42- + HCO3-)/(Ca2+ + Mg2+)比例可以揭示水体水化学成分是否受硅酸盐、碳酸盐和硫酸盐矿物溶解的影响:水样点位位于(SO42- + HCO3-)/(Ca2+ + Mg2+)比例等值线上, 表明碳酸盐和硫酸盐矿物溶解起主导作用;位于等值线下方或附近, 说明水体水化学成分还会受到硅酸盐矿物溶解影响[33]. Ca2+/HCO3-比例位于1∶1和1∶2之间, 则表明水体中Ca2+和HCO3-主要来自于碳酸盐风化[34].(Na+ + K+ - Cl-)/(SO42- + HCO3- - Ca2+ - Mg2+)毫克当量浓度比例可以用来确定水体是否存在阳离子交换:水样点位接近等值线, 说明水化学成分受到阳离子交换作用的影响[35]. 此外, 氯碱指数(CAI)可以进一步判断水体中是否存在正向或反向的阳离子交换过程:CAI < 0表明水体存在正向的阳离子交换过程, CAI > 0表明水体存在反向的阳离子交换[36]. 同时, 水体矿物饱和指数(SI)能够指示水体矿物饱和度:SI > 0和SI < 0分别表示矿物过饱和和不饱和状态. 从图 4(a)可以看出, 涡河流域地表水、浅层地下水和中深层地下水大部分点位的Cl-/Na+比例均位于等值线以下, 说明涡河流域各类水体都受到硅酸盐矿物风化和阳离子交换作用的影响, 且地下水点位相比于地表水更加远离等值线, 证明地下水岩相互作用会进一步导致地下水中Na+浓度增加. 研究区地表水、浅层地下水和中深层地下水大多数点位的(SO42- + HCO3-)/(Ca2+ + Mg2+)比例均位于等值线上方和附近, 进一步说明地表和地下水中硅酸盐矿物的风化溶解是主要影响因素[图 4(b)]. 如图 4(c)所示, 地表水点位均位于y = xy = 2x之间, 表明Ca2+和HCO3-主要来源于碳酸盐矿物风化, 而浅层和中深层地下水点位只有部分位于y = xy = 2x区间内, 表明在地下水中Ca2+除了来源于碳酸盐风化, 还可能源自于阳离子交换作用或含钙硅酸盐矿物的溶解. 在(Na+ + K+ - Cl-)/(HCO3- + SO42- - Ca2+ - Mg2+)比例二元图中, 涡河流域地表水和地下水样品大多数分布在y = x附近, 说明涡河流域不同类型水体水化学受到了阳离子交换作用的影响[图 4(d)]. 图 4(e)显示涡河流域地表水和中深层地下水的全部水样点位和浅层地下水的大部分点位均落在CAI < 0的区域, 说明正向阳离子交换作用的影响占主导地位. 此外, 研究区所有水样中石膏和盐岩的SI值都小于0[图 4(f)], 表明地表水和地下水中这些矿物质均处于不饱和状态. 其中, 方解石和白云石SI值在地表水和中深层地下水中都为负值, 证明方解石和白云石在这两类水体中处于持续溶解的状态;而在少部分浅层地下水点位中方解石和白云石SI值大于0, 表明部分区域浅层地下水中方解石和白云石处于过饱和状态.

图 4 涡河流域各类水体主要离子关系 Fig. 4 Major ionic relationships of various water bodies in the Guohe River Basin

3.2 地下水氢氧同位素分布和相关性研究

通过Arcgis 10.8的反距离权重法绘制了地下水氢氧同位素值的空间分布[图(5)]. 结果显示, 浅层和中深层地下水δ18O和δ2H值的空间分布规律高度一致, 证明了地下水δ18O和δ2H值之间存在很强的相关性. 研究区地下水δ18O和δ2H值的空间分布上均呈现从西北向东南方向逐渐增大的趋势, 这种趋势与涡河流域的地理形态特征和地下水流动方向有显著的一致性. 地下水在其流动过程中会持续与含水层的岩石产生相互作用, 且蒸发过程中重同位素会出现富集现象, 这两种因素共同推动了上述δ18O和δ2H分布趋势的形成. 在涡河流域中间区域, 水体δ18O和δ2H值出现高值, 主要是因为涡阳枢纽作为一个大型的水闸位于该区域. 在通常情况下, 该水闸处于关闭状态, 使得上游地表水在此处滞留, 进而导致地下水中δ18O和δ2H逐渐富集.

图 5 涡河流域地下水氢氧同位素空间分布 Fig. 5 Spatial distribution of hydrogen and oxygen isotopes in groundwater in the Guohe River Basin

涡河流域地下水氢氧同位素值与TDS和主要阴阳离子浓度间的相关性分析能够进一步揭示地下水中氢氧同位素的分布成因. 地下水氢氧同位素值与TDS和离子浓度均呈显著相关, 则说明地下水氢氧同位素值主要受到蒸发效应的影响, 蒸发分馏会导致氢氧同位素数值增加, TDS和离子浓度也会因为随着水体蒸发浓缩作用而提高;反之, 说明地下水受到补给来源的影响大于蒸发效应[37]. 根据表 2所示的Pearson相关性分析结果, 涡河流域地下水TDS与δ18O和δ2H之间没有显著相关性, 表明地下水较蒸发效应而言主要受到补给来源的影响. 同时, 研究区地下水δ18O、δ2H值与K+、Na+、Cl-和SO42-浓度呈现显著负相关, 这可能是由于地下水在流动过程中会与土壤和岩石发生相互作用, 这些离子在地表水补给地下水的过程中会被吸附于土壤和岩石中或与之发生化学反应, 导致其在地下水中的浓度相应降低. 该现象反映了地下水中离子浓度会在地层过滤过程中降低. 此外, 地下水中δ18O、δ2H值与NO3-浓度呈显著正相关. NO3-主要来源于农业施肥和畜禽养殖等人类活动, 在地下水δ18O和δ2H值增加的同时, 受到了较多人类活动废水的补给, 这些水源可能含有较高浓度的NO3-, 使得地下水中NO3-浓度同步增加, 反映出了人为活动对地下水化学组成的影响. 由此可以推断, 安徽涡河流域地下水除了受到大气降水的补给外, 地表水以及人类活动产生的废水也对地下水产生一定补给作用. 此外, 该区域地下水δ18O与δ2H值之间呈极显著正相关, 相关性系数为0.90. 因此, 地下水资源的保护和科学管理意义重大, 需要考虑到多种补给源的影响, 并加强对区域人类活动污染的监管, 以确保地下水资源的可持续利用.

表 2 涡河流域地下水氢氧同位素皮尔逊相关性分析1) Table 2 Pearson correlation analysis of hydrogen and oxygen isotopes in groundwater of the Guohe River Basin

3.3 基于氢氧同位素的各类水体补给来源和转化关系分析

通过氢氧同位素技术可以确定不同水体的补给来源. 南京和郑州与亳州在地理位置、年平均降水量和年平均蒸发量等方面具有相似性[38], 因此在缺乏亳州当地的同位素观测数据的情况下, 本研究采用南京和郑州2个国际原子能机构观测站1985~1992年的112个降水氢氧同位素数据来计算亳州的当地大气降水线(LMWL). 计算得出的亳州LMWL为δ2H = 7.34 δ18O + 5.59(n = 112, R2 = 0.89, 图 6). 亳州处于暖温带半湿润季风气候, 年蒸发量常大于降水量, 因此该区域降水氢氧同位素存在一定分馏, 使得亳州LMWL的斜率低于全球大气降水线[39](GMWL, 方程为δ2H = 8 δ18O +10). 对涡河流域地表水和地下水的δ18O与δ2H关系进行分析(图 6), 发现涡河地表水的所有水样点位均位于LMWL右下方, 且分布较为集中, 表明研究区地表水的主要补给来源之一为大气降水[40]. 浅层和中深层地下水样品点位大多位于LMWL附近, 少数点位处于LMWL上部, 反映了大气降水也是研究区地下水的主要补给来源. 相较于地表水, 涡河流域地下水氢氧同位素更加贫化并且分布更散乱, 说明地下水在经历长期的水岩作用和转化过程中发生了氢氧同位素分馏. 地表水的δ18O与δ2H线性关系分别为δ2H = 1.45 δ18O - 45.18(n = 32, R2 = 0.19). 浅层和中深层地下水的δ18O与δ2H线性关系分别为δ2H = 6.24 δ18O - 3.72(n = 15, R2 = 0.80)和δ2H = 6.82 δ18O - 6.16(n = 12, R2 = 0.93), 两者斜率均低于LMWL, 且该斜率在浅层地下水中小于中深层地下水, 进一步验证了地下水存在蒸发分馏效应. 具体而言, 浅层地下水距离地表近, 蒸发作用更为明显, 而中深层地下水埋深较大, 封闭性高, δ18O和δ2H值的影响因素少, 因此蒸发作用较不明显, 斜率也更接近亳州LMWL.

图 6 涡河流域各类水体δ18O与δ2H关系 Fig. 6 Relationships between δ18O and δ2H of various water bodies in the Guohe River Basin

进一步地, 利用MixSIAR模型分析涡河流域各类水体的δ18O和δ2H数据, 以揭示并量化各水体间的转化关系. 研究区中深层地下水与浅层地下水和地表水的水力联系小, 将涡河流域(河南段)的中深层地下水氢氧同位素值作为研究区中深层地下水的侧向补给来源数据[41]. 模型结果显示, 在安徽涡河流域, 大气降水和浅层地下水对地表水的补给贡献率分别为46.5%和53.5%, 大气降水和地表水对浅层地下水的补给贡献率分别为57.4%和42.6%, 浅层地下水和上游地下水侧向补给对中深层地下水的补给贡献率分别为9.0%和77.0%. 由此表明, 涡河流域地表水和浅层地下水主要由大气降水补给, 中深层地下水主要依赖上游地下水的侧向径流补给.

4 结论

(1)安徽涡河流域地表水及地下水整体呈现中性至弱碱性. 涡河地表水水化学类型以Cl·SO4·HCO3-Na和Cl·SO4-Na型为主, 而浅层地下水为HCO3-Ca·Mg和HCO3-Mg·Na型, 中深层地下水以Cl·HCO3-Na型为主.

(2)岩石风化(包括盐岩、硅酸盐、碳酸盐和硫酸盐矿物)、蒸发浓缩及正向阳离子交换对研究区内地表水和地下水的主要阴阳离子分布影响显著. 大部分水体中的石膏、盐岩、方解石和白云石均呈现不饱和状态.

(3)各类水体中δ18O和δ2H平均值的大小顺序为:涡河地表水 > 浅层地下水 > 中深层地下水. δ18O和δ2H值在水平空间分布上呈现西北低、东南高的趋势. 此外, δ18O和δ2H的线性关系表明大气降水是涡河流域水体的主要补给来源. 地下水δ18O、δ2H值与K+、Na+、Cl-和SO42-浓度呈显著负相关, 与NO3-呈显著正相关, 阐明了研究区地层过滤过程中地下水中离子浓度的降低以及人为活动对地下水化学组成的影响.

(4)MixSIAR模型揭示了涡河流域地表水的主要来源为大气降水(贡献率为46.5%)和浅层地下水(贡献率为53.5%);浅层地下水的主要补给来自大气降水(贡献率为57.4%)和涡河地表水(贡献率为42.6%);中深层地下水主要依赖上游地下水的侧向径流补给.

参考文献
[1] Kadam A, Wagh V, Jacobs J, et al. Integrated approach for the evaluation of groundwater quality through hydro geochemistry and human health risk from Shivganga river basin, Pune, Maharashtra, India[J]. Environmental Science and Pollution Research, 2022, 29(3): 4311-4333. DOI:10.1007/s11356-021-15554-2
[2] Zhou Y H, Li P Y, Chen M J, et al. Groundwater quality for potable and irrigation uses and associated health risk in southern part of Gu'an County, North China Plain[J]. Environmental Geochemistry and Health, 2021, 43(2): 813-835. DOI:10.1007/s10653-020-00553-y
[3] 刘鑫, 向伟, 司炳成. 汾河流域浅层地下水水化学和氢氧稳定同位素特征及其指示意义[J]. 环境科学, 2021, 42(4): 1739-1749.
Liu X, Xiang W, Si B C. Hydrochemical and isotopic characteristics in the shallow groundwater of the Fenhe River Basin and indicative significance[J]. Environmental Science, 2021, 42(4): 1739-1749.
[4] 朱金峰, 刘悦忆, 章树安, 等. 地表水与地下水相互作用研究进展[J]. 中国环境科学, 2017, 37(8): 3002-3010.
Zhu J F, Liu Y Y, Zhang S A, et al. Review on the research of surface water and groundwater interactions[J]. China Environmental Science, 2017, 37(8): 3002-3010. DOI:10.3969/j.issn.1000-6923.2017.08.024
[5] 李双慧, 黄选明, 朱宏军, 等. 基于水化学和环境同位素的准格尔煤田地下水循环特征[J]. 南水北调与水利科技(中英文), 2021, 19(3): 561-571.
Li S H, Huang X M, Zhu H J, et al. Characterization of groundwater circulation in the Jungar coalfield based on water chemistry and environmental isotopes[J]. South-to-North Water Transfers and Water Science & Technology, 2021, 19(3): 561-571, 589.
[6] 于品正. 基于水污染指数法的涡河(蒙城段)水质评价分析[J]. 水利科技与经济, 2022, 28(8): 67-70.
Yu P Z. Evaluation and analysis of water quality of Guohe (Mengcheng section) based on water pollution index method[J]. Water Conservancy Science and Technology and Economy, 2022, 28(8): 67-70. DOI:10.3969/j.issn.1006-7175.2022.08.015
[7] 胡铮琪. 近5年皖北地表水水质变化趋势及治理措施—以涡河亳州境内水质变化为例[J]. 绿色科技, 2021, 23(2): 96-99.
Hu Z Q. Changing trend and control measures of surface water quality in northern Anhui province in the past five years-taking the change of Water Quality in Guohe River, Bozhou city as an example[J]. Journal of Green Science and Technology, 2021, 23(2): 96-99. DOI:10.3969/j.issn.1674-9944.2021.02.034
[8] 叶永红, 龚建师, 许乃政, 等. 涡河流域河南段浅层高碘地下水分布及成因[J]. 华东地质, 2023, 43(3): 292-299.
Ye Y H, Gong J S, Xu N Z, et al. Distribution and genesis of high iodine shallow groundwater in Henan section of Guohe Basin[J]. East China Geology, 2023, 43(3): 292-299.
[9] 黄静, 周忠发, 丁圣君, 等. 双河洞流域氢氧同位素特征及其指示意义[J]. 地球与环境, 2022, 50(4): 516-525.
Huang J, Zhou Z F, Ding S J, et al. Hydrogen and oxygen isotope characteristics and their indicative significance in Shuanghe Cave Basin[J]. Earth and Environment, 2022, 50(4): 516-525.
[10] 成玉婷, 李鹏, 徐国策, 等. 丹江流域氢氧同位素变化特征[J]. 水土保持学报, 2014, 28(5): 129-133.
Cheng Y T, Li P, Xu G C, et al. Characteristics of hydrogen and oxygen isotopes in Danjiang watershed[J]. Journal of Soil and Water Conservation, 2014, 28(5): 129-133.
[11] Kong X L, Wang S Q, Liu B X, et al. Impact of water transfer on interaction between surface water and groundwater in the lowland area of North China Plain[J]. Hydrological Processes, 2018, 32(13): 2044-2057. DOI:10.1002/hyp.13136
[12] 雷米, 周金龙, 张杰, 等. 新疆博尔塔拉河流域平原区地表水与地下水水化学特征及转化关系[J]. 环境科学, 2022, 43(4): 1873-1884.
Lei M, Zhou J L, Zhang J, et al. Hydrochemical characteristics and transformation relationship of surface water and groundwater in the plain area of Bortala River Basin, Xinjiang[J]. Environmental Science, 2022, 43(4): 1873-1884.
[13] 苏小四, 林学钰. 银川平原地下水循环及其可更新能力评价的同位素证据[J]. 资源科学, 2004, 26(2): 29-35.
Su X S, Lin X Y. Cycle pattern and renewablity evaluation of groundwater in Yinchuan Basin: isotopic evidences[J]. Resources Science, 2004, 26(2): 29-35. DOI:10.3321/j.issn:1007-7588.2004.02.005
[14] 王忠亮, 郭春艳, 张彦鹏. 涞源北盆地地下水氢氧同位素特征及北海泉形成模式[J]. 水文地质工程地质, 2021, 48(1): 27-35.
Wang Z L, Guo C Y, Zhang Y P. Characteristics of hydrogen and oxygen isotopes in the groundwater and formation mode of the Beihai springs in the northern Laiyuan Basin[J]. Hydrogeology & Engineering Geology, 2021, 48(1): 27-35.
[15] 翟婧雅, 金彦香, 金鑫. 巴音河流域水化学与氢氧同位素特征研究[J]. 灌溉排水学报, 2022, 41(11): 101-106.
Zhai J Y, Jin Y X, Jin X. Hydrochemical property and hydrogen and oxygen isotopes in the Bayin River Basin[J]. Journal of Irrigation and Drainage, 2022, 41(11): 101-106.
[16] 杨楠, 苏春利, 曾邯斌, 等. 基于水化学和氢氧同位素的兴隆县地下水演化过程研究[J]. 水文地质工程地质, 2020, 47(6): 154-162.
Yang N, Su C L, Zeng H B, et al. Evolutional processes of groundwater in Xinglong County based on hydrochemistry and hydrogen and oxygen isotopes[J]. Hydrogeology & Engineering Geology, 2020, 47(6): 154-162.
[17] 蒲俊兵. 重庆岩溶地下水氢氧稳定同位素地球化学特征[J]. 地球学报, 2013, 34(6): 713-722.
Pu J B. Hydrogen and oxygen isotope geochemistry of karst groundwater in Chongqing[J]. Acta Geoscientica Sinica, 2013, 34(6): 713-722.
[18] Evaristo J, Jasechko S, Mcdonnell J J. Global separation of plant transpiration from groundwater and streamflow[J]. Nature, 2015, 525(7567): 91-94. DOI:10.1038/nature14983
[19] Bowen G J, Good S P. Incorporating water isoscapes in hydrological and water resource investigations[J]. WIREs Water, 2015, 2(2): 107-119. DOI:10.1002/wat2.1069
[20] Jiang W J, Wang G C, Sheng Y Z, et al. Isotopes in groundwater (2H, 18O, 14C) revealed the climate and groundwater recharge in the Northern China[J]. Science of the Total Environment, 2019, 666: 298-307. DOI:10.1016/j.scitotenv.2019.02.245
[21] 马涛, 张朝君. 安徽涡河地区浅层地下水化学特征分析[J]. 地下水, 2021, 43(5): 22-24, 31.
Ma T, Zhang C J. Chemical characteristics of shallow groundwater in Guo River area, Anhui province[J]. Ground Water, 2021, 43(5): 22-24, 31.
[22] 郑涛, 焦团理, 胡波, 等. 涡河流域中部地区地下水化学特征及其成因分析[J]. 环境科学, 2021, 42(2): 766-775.
Zheng T, Jiao T L, Hu B, et al. Hydrochemical characteristics and origin of groundwater in the central Guohe River Basin[J]. Environmental Science, 2021, 42(2): 766-775.
[23] 张志雄, 王仕琴, 张依章, 等. 雄安新区唐河污水库残留污染物对地下水水化学动态的作用机制[J]. 环境科学, 2021, 42(11): 5312-5321.
Zhang Z X, Wang S Q, Zhang Y Z, et al. Dynamic mechanisms of groundwater quality by residual contaminants of the tanghe wastewater reservoir in Xiong'an New Area[J]. Environmental Science, 2021, 42(11): 5312-5321.
[24] Wang J, Lu N, Fu B J. Inter-comparison of stable isotope mixing models for determining plant water source partitioning[J]. Science of the Total Environment, 2019, 666: 685-693. DOI:10.1016/j.scitotenv.2019.02.262
[25] 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
[26] 肖春艳, 刘怡凡, 赵同谦, 等. 小浪底水库水沙调控期滨河湿地地下水与河水转化关系[J]. 环境科学, 2022, 43(12): 5499-5508.
Xiao C Y, Liu Y F, Zhao T Q, et al. Transformation relationship of groundwater and river water in riparian wetland during water and sediment regulation of xiaolangdi reservoir in yellow river[J]. Environmental Science, 2022, 43(12): 5499-5508.
[27] 金可, 张乾柱, 卢阳, 等. 巴丹吉林沙漠湖泊群水体氢氧同位素和水化学特征[J]. 人民长江, 2022, 53(4): 65-72.
Jin K, Zhang Q Z, Lu Y, et al. Research on stable isotopes and hydrochemical features of lakes water in Badain Jaran Desert[J]. Yangtze River, 2022, 53(4): 65-72.
[28] 谢姝, 孟玉川, 刘国东, 等. 大渡河丰水期水体氢氧同位素特征及其影响因素[J]. 人民长江, 2022, 53(3): 55-60.
Xie S, Meng Y C, Liu G D, et al. Variations of hydrogen and oxygen isotopes in waterbody of Dadu River during flood period and their influencing factors[J]. Yangtze River, 2022, 53(3): 55-60.
[29] Gibbs R J. Mechanisms controlling world water chemistry[J]. Science, 1970, 170(3962): 1088-1090. DOI:10.1126/science.170.3962.1088
[30] 孔晓乐, 杨永辉, 曹博, 等. 永定河上游地表水-地下水水化学特征及其成因分析[J]. 环境科学, 2021, 42(9): 4202-4210.
Kong X L, Yang Y H, Cao B, et al. Hydrochemical characteristics and factors of surface water and groundwater in the upper Yongding River Basin[J]. Environmental Science, 2021, 42(9): 4202-4210.
[31] 李政, 范成博, 尚琳群, 等. 石家庄平原区浅层地下水化学特征及成因分析[J]. 河北地质大学学报, 2023, 46(3): 71-79.
Li Z, Fan C B, Shang L Q, et al. Chemical characteristics and genesis analysis of shallow groundwater in the plain area of Shijiazhuang[J]. Journal of Hebei GEO University, 2023, 46(3): 71-79.
[32] 艾力哈木·艾克拉木, 周金龙, 张杰, 等. 伊犁河谷西北部地下水化学特征及成因分析[J]. 干旱区研究, 2021, 38(2): 504-512.
Ailihamu A K L M, Zhou J L, Zhang J, et al. Chemical characteristics and genesis analysis of groundwater in northwest Yili River Valley[J]. Arid Zone Research, 2021, 38(2): 504-512.
[33] 冯建国, 鲁统民, 高宗军, 等. 新泰市地下水水化学特征及成因探讨[J]. 山东科技大学学报(自然科学版), 2020, 39(1): 11-20.
Feng J G, Lu T M, Gao Z J, et al. Hydrochemical characteristics and causes of groundwater in Xintai city[J]. Journal of Shandong University of Science and Technology (Natural Science), 2020, 39(1): 11-20.
[34] Li P Y, Wu J H, Qian H. Hydrochemical appraisal of groundwater quality for drinking and irrigation purposes and the major influencing factors: a case study in and around Hua County, China[J]. Arabian Journal of Geosciences, 2016, 9(1). DOI:10.1007/s12517-015-2059-1
[35] Xiao J, Jin Z D, Wang J, et al. Hydrochemical characteristics, controlling factors and solute sources of groundwater within the Tarim River Basin in the extreme arid region, NW Tibetan Plateau[J]. Quaternary International, 2015, 380-381: 237-246. DOI:10.1016/j.quaint.2015.01.021
[36] Liu J T, Peng Y M, Li C S, et al. Characterization of the hydrochemistry of water resources of the Weibei Plain, Northern China, as well as an assessment of the risk of high groundwater nitrate levels to human health[J]. Environmental Pollution, 2021, 268. DOI:10.1016/j.envpol.2020.115947
[37] 汪敬忠, 吴敬禄, 曾海鳌, 等. 内蒙古河套平原水体同位素及水化学特征[J]. 地球科学与环境学报, 2013, 35(4): 104-112.
Wang J Z, Wu J L, Zeng H A, et al. Characteristics of water isotope and hydrochemistry in Hetao plain of inner Mongolia[J]. Journal of Earth Sciences and Environment, 2013, 35(4): 104-112.
[38] 郑涛, 侯香梦, 管后春, 等. 氢氧同位素在亳州东南部地区水循环研究中的应用[J]. 安徽地质, 2021, 31(3): 280-285.
Zheng T, Hou X M, Guan H C, et al. Application of hydrogen and oxygen isotopes in the study of water cycle in southeastern Bozhou[J]. Geology of Anhui, 2021, 31(3): 280-285.
[39] Craig H. Isotopic variations in meteoric waters[J]. Science, 1961, 133(3465): 1702-1703.
[40] Guo Q L, Yang Y S, Han Y Y, et al. Assessment of surface–groundwater interactions using hydrochemical and isotopic techniques in a coalmine watershed, NW China[J]. Environmental Earth Sciences, 2019, 78(3). DOI:10.1007/s12665-019-8053-2
[41] 李潇. 涡河流域(河南段)地下水循环特征研究[D]. 郑州: 华北水利水电大学, 2020.
Li X. Study on groundwater circulation characteristics in Guohe River Basin (Henan section)[D]. Zhengzhou: North China University of Water Resources and Electric Power, 2020.