1. 引言
水资源是人类赖以生存和发展所必需的资源,自然界中地表水与地下水之间的密切作用,直接导致两者的量与质的变化,并对人类生活及社会发展产生影响。因此,地表水与地下水之间的相互作用一直是地下水研究领域的一个重要课题 [1]。系统地研究地表水与地下水之间的转化关系,对水资源合理开发利用、水资源调控及水资源污染防治等都具有重要作用。
水化学是研究各水体成分的形成机理、发展趋势、演变规律及水文循环的基础工具和重要手段。自然状态下水体中的氢氧同位素具有稳定的化学性质,是研究天然状态下水运动过程的理想示踪剂 [2]。地表水与地下水相互作用影响着各自水化学成分的分布和演变规律以及氢氧稳定同位素数值的变化,因此国内外许多学者应用水化学方法和同位素技术研究流域内地表水与地下水之间的转化关系 [3] - [9],并利用同位素来定量计算转化程度 [10] [11]。
前人对挠力河流域内地表水与地下水的水化学及同位素的特征进行过分析 [12],但对于两者间的转化关系多利用水位资料和模拟模型来确定 [13] [14],利用水化学和同位素手段在不同时空条件下分析转化关系的研究较少。因此,本文在前人研究基础上,通过水化学和氢氧稳定同位素的手段,对挠力河干流在不同区段、不同时期下的河水与地下水的转化关系进行分析,并定量计算转化程度,为合理利用挠力河流域水资源等提供科学依据。
2. 研究区概况
挠力河发源于长白山脉北延支脉——完达山,汇入乌苏里江,全长596 km,为雨雪混合补给的河流,平均年径流量为19.8 × 108 m3。三江平原挠力河流域位于黑龙江省东北腹地,属中温带湿润大陆性季风气候区,年内降雨分配不均,降雨多集中于夏季,多年平均降雨量为534.5 mm,多年平均蒸发量约为732.2 mm。地形主要呈西南高,东北低,上游为山丘区,河道坡度较大,河道曲折;中游为平原区,河道坡度极缓,河道呈现极度弯曲状,中游区域河曲、牛轭湖、湖泊及湿地较多;下游为平洼区,主槽坡降在0.07‰~0.03‰。
挠力河流域,平原广袤农业发达,地下水开发利用程度较高,是农业用水的主要水源,地下水处于超量开采状态,部分地段地下水位埋深持续增加 [15],地表水与地下水的联合利用和科学管理问题亟待解决。
3. 数据与方法
3.1. 采样与分析
2018年4、7和10月采集不同河段河水、近岸地下水及远岸地下水 [16] 各14个水样,采样点如图1所示。水化学数据由吉林大学测试中心及中国科学院东北地理与农业生态研究所测定,同位素数据由美国Beta Analytic公司进行分析,δ18O和δD的精度达到0.3‰,测量结果以相对维也纳标准海水(VSMOW)的千分偏差表示。
![](//html.hanspub.org/file/5-2410877x10_hanspub.png)
Figure 1. Distribution map of Naoli River basin and water samples
图1. 挠力河流域及水样点分布图
3.2. 分析方法
整理水样点的水化学和同位素数据并绘制Piper图和δ18O-δD关系图,绘制沿挠力河流向的TDS(总溶解固体)、δ18O的变化曲线和各区段各监测期水样点的δ18O、Cl−分布图,考虑地形及水文地质条件等因素,分析不同区段、时期下挠力河干流河水与地下水之间的转化关系。根据质量守恒原理,利用同位素做示踪剂,定量估算河水与地下水之间的转化程度,公式如下:
(1)
(2)
(3)
在计算河水入渗补给地下水时:式中
为河水流量,
为远岸地下水流量,
为近岸地下水流量;
为河水中同位素浓度,
为远岸地下水中同位素浓度,
为近岸地下水中同位素浓度,
为河水入渗地下水的补给比例 [17]。
在计算地下水补给河水时:式中
为地下水流量,
为上游河水流量,
为下游河水流量;
为地下水中同位素浓度,
为上游河水中同位素浓度,
为下游河水中同位素浓度,
为地下水对河水的补给比例 [18]。
4. 结果
4.1. 水化学组成
挠力河流域内,各水体的水化学数据分布如表1所示,全部水样点的pH值均分布在弱酸、弱碱性范围,TDS均小于1 g/L。挠力河流域的浅层含水层主要由砂砾石和细砂组成,上覆约10 m厚的粉质黏土 [19],流域内地表植被丰富,腐殖层较厚,降水入渗过程中的水岩作用复杂,导致地下水TDS含量明显高于河水。
![](Images/Table_Tmp.jpg)
Table 1. Statistics of chemical and isotopic data of samples from Naoli River
表1. 挠力河水样点水化学和同位素数据统计
由Piper图(图2所示)看出,
为河水和地下水的优势阴离子,占总阴离子的40.23%~50%以上;阳离子以Ca2+为主。除
和Ca2+外,河水水样点
和Na++K+含量相对较高,地下水Cl−和Mg2+含量相对较高。
4.2. 同位素组成
挠力河流域全部水样点均分布在当地大气降水线附近,表明大气降水是各水体的主要来源之一,河水水样点大多数分布在降水线的右下侧,且气温较高、蒸发更为强烈的7月同位素更为富集,10月次之,4月最为贫化,反映出河水在受当地降水补给过程中,受到更为显著的蒸发作用影响。通过图3可以看出,三种水体水样点分布区域大致相同,且河水水样点分布最为紧密。因近岸地下水较远岸地下水与河水联系更加紧密,故近岸地下水同位素分布较为分散,远岸地下水则较为紧密。
![](//html.hanspub.org/file/5-2410877x33_hanspub.png)
Figure 2. Piper diagram of all water samples
图2. 全部水样点Piper图
![](//html.hanspub.org/file/5-2410877x34_hanspub.png)
Figure 3. δ18O-δD diagram of all water sample points
图3. 全部水样点δ18O-δD关系图
5. 讨论
5.1. 不同区段河水与地下水转化关系
挠力河流域河水与地下水的转化关系具有空间变异特性。根据水化学和同位素数据,在同一监测期,沿挠力河流向绘制TDS和δ18O的变化曲线,如图4所示。
4月,上游区段,起始断面处河水的δ18O数值高于地下水,沿河流流向,河水的δ18O下降明显,表明受到了同位素较为贫化的地下水的补给;近岸地下水受到离子浓度含量较低的远岸地下水的补给,TDS出现下降趋势。中游区段,同位素方面,近岸地下水呈现增长趋势,河水受到较贫化的冰雪融水和大气降水的补给,弱化蒸发作用的结果,使河水的δ18O值保持不变;水化学方面,因地形原因,近岸地下水受远端贫化地下水补给较少,TDS值出现较大幅度的增加,同时侧向补给河水使其呈现增长趋势,因此本区段为地下水补给河水。下游区段,地下水的δ18O和TDS数值迅速下降,河水变化不大,表明河水大量补给地下水。
7月,上游区段地下水TDS始终高于河水且呈增长趋势,受到蒸发作用和降水入渗补给,地下水同位素值出现升高趋势,而受到蒸发作用大于降水补给的河水,其同位素值未出现上升趋势,表明河水受到较贫化的地下水补给。中上游区段,河水的同位素富集作用强烈,远大于近岸地下水,近岸地下水受到远端贫化的地下水补给出现下降趋势;水化学方面,矿化度较高的地下水受矿化度较低的河水补给,出现减小并趋近河水变化的趋势。中下游区段,受到δ18O值较高的七星河补给,挠力河河水同位素出现上升趋势,两种水体水化学变化趋势不大,考虑夏季人均用水量增加及左岸广阔的农耕区,水田灌溉面积占总灌溉面积的98%,地下水开采量达54.56 × 108 m3 [20],并参考地形及水位曲线,判断为河水补给地下水。
10月水样点分布与7月相同,河水与地下水的转化关系也具有相似性。上游区段,地下水TDS受灌区农作物追肥等人类活动导致离子浓度增加,且高矿化度的河水受到相对贫化的地下水补给,出现下降趋势。中上游区段,两种水体的TDS变化趋势不大,同位素较富集的河水,其δ18O随地下水出现下降的趋势,说明地下水补给河水。中下游区段,考虑地形和用水情况,为河水补给地下水。
![](//html.hanspub.org/file/5-2410877x36_hanspub.png)
![](//html.hanspub.org/file/5-2410877x35_hanspub.png)
![](//html.hanspub.org/file/5-2410877x37_hanspub.png)
Figure 4. TDS andδ18O curves of river and inshore groundwater during each monitoring period
图4. 各监测期河水与近岸地下水TDS和δ18O沿程变化曲线
5.2. 不同时期河水与地下水转化关系
挠力河流域内河水与地下水的转化关系具有时间动态变化性。在确定不同区段的转化关系后,绘制各监测期各断面水样点δ18O和Cl−分布图,如图5所示。
上游区段,7月河水受蒸发作用远大于大气降水和地下水补给作用,δ18O变化不大,以富集为主。4、10月河水受到同位素较贫化的冰雪融水、大气降水和地下水的补给程度大于蒸发作用,同位素出现贫化的趋势;水化学方面,河水离子浓度与地下水相类似,出现贫化的现象,说明受到了地下水的补给。表明4、7和10月上游区段为地下水补给河水,且河水左右两岸地下水都向河水补给。
中上游区段,4月地下水受水岩作用及施肥等人类活动导致Cl−较为富集;同位素方面,受到融雪水和地下水补给,上游同位素较富集的河水出现贫化趋势。7月,δ18O从河水向两侧地下水逐渐减小,河水Cl−含量最高。挠力河左岸地下水同位素和水化学指标较高于右岸,且有靠近河水的趋势。考虑近岸地下水与远岸地下水受山地丘陵地势的影响,水力联系相对不畅,使得近岸地下水各项指标最低,推测转化关系为挠力河右岸补给河水,河水补给挠力河左岸。10月同位素及水化学分布与7月呈现一致性,δ18O指标上河水最高,远岸地下水最低;Cl−指标上河水最高,近岸地下水最低。但河水同7月相比,与地下水在数值上较为相近,不仅由蒸发作用减弱导致,还与受到地下水补给有关。
中下游区段,4月中上游断面处δ18O和Cl−数值较低的河水受到数值较高的地下水补给,在本区段河水的同位素和水化学都呈现增长趋势;7、10月中上游断面处δ18O和Cl−数值较低的地下水皆因受到河水的补给,各项指标趋近与河水水体。
下游区段,仅涉及4月水样点,河水的Cl−最高,δ18O最小,地下水与河水呈现相反的趋势。挠力河下游区域,因城市人口密集,农业水平发达,地下水用水量较大,考虑降落漏斗的出现,本区段为河水补给地下水。
![](//html.hanspub.org/file/5-2410877x38_hanspub.png)
Figure 5. δ18O and Cl- distribution diagram of water samples in each section and monitoring period
图5. 各区段各监测期水样点δ18O和Cl-分布图
5.3. 转化量计算
在确定转化关系的基础上,根据质量守恒原理,选择δ18O数据,利用公式(3)计算,可得各监测期各区段下河水与地下水的转化关系及补给比例,如表2所示。
4月上游、中游区段,以地下水补给河水为主,补给比例为77.27%和66.67%。因4月为非汛期,流域降水量甚少,河道基流主要由地下水补给,计算结果与正常认知相同。下游地下水开采量大,导致河水补给地下水,其补给比例为38.43%。7月为汛期,降水量较多,河流水位普遍高于地下水,因此中上、下游河水对地下水进行补给,补给比例为15.04%和11.82%。10月降水量略有减少,补给情况受地形条件的影响,故中上游地下水补给河水,补给比例为15.95%,中下游以河水补给地下水为主,补给比例为16.18%。
![](Images/Table_Tmp.jpg)
Table 2. Transformation relationship and recharge ratio of every section in every monitoring period of Naoli River
表2. 挠力河各监测期各区段转化关系及补给比例
6. 结论
挠力河干流河水与地下水水化学类型相近,优势阴、阳离子都为
、Ca2+;同位素数据均分布在当地大气降水线附近。流域内湖泊湿地众多,河道中游及中下游坡度极缓,呈现极度弯曲状,河水与地下水交互频繁。两者的转化关系主要受地形地势及水文地质条件的影响,同时也受到补给源和补给量的影响,在非汛期转化关系主要由地势和水文地质条件决定,在汛期河水受大气降水及冰雪融水的大量补给,水位上升,改变了原有的转化关系。
1) 河水与地下水的转化关系具有时间和空间特性:4月以地下水补给河水为主,7月以河水补给地下水为主,10月两者相近;上游区段为地下水补给河水,中游区段转化关系与季节相关,下游因地下水开采量大,为河水补给地下水。
2) 4月上游及中游区段地下水补给河水,地下水补给比例为77.27%和66.67%,下游河水补给地下水,其补给比例为38.43%;7月上游地下水补给河水,中上游及中下游河水补给地下水,河水补给比例为15.04%和11.82%,同时挠力河受到七星河补给,其补给比例为32.50%;10月上游及中上游地下水补给河水,地下水补给比例为15.95%,中下游河水补给地下水,河水补给比例为16.18%,七星河补给挠力河,补给比例为18.99%。
基于挠力河干流河水与地下水转化关系的确定,对当地水资源合理利用提出如下建议:河流上游河道受流水侵蚀作用较强,地下水补给河水,应注意植被保护,防止水土流失;中游转化关系受降雨等影响,可建立地下水–河水联合利用模式,并加强统一管理调度;下游因农业和生活用水,开采地下水导致降落漏斗,建议适度减少地下水开采、同时加强地表水灌溉规模,维持水资源与农业和生态环境的可持续发展。
基金项目
国家重点研发计划“大面积农业灌溉的地表水与地下水联合调控”(2017YFC0406002)。