基于人工甜味剂时序特征的白洋淀对地下水补给量研究

马 荣， 张 冰， 周晓妮

(1.南京水利科学研究院 水文水资源与水利工程科学国家重点实验室， 江苏 南京 210098；2.中国地质科学院 第四纪年代学与水文环境演变重点实验室， 河北 石家庄 050061)

1 研究背景

设立雄安新区，是重大的历史性战略选择，是千年大计、国家大事[1]。白洋淀是雄安新区内重要的生态要素之一，与大气圈、生物圈和岩石圈之间进行着密切的物质循环和能量转递，是各圈层相互作用的重要连接点，对雄安新区生态水城的建设具有重要意义。白洋淀内纵横交错的3 700余条沟壕将全淀分割成143个大小不等又相互连通的淀泊，均与地下水存在着密切的水力联系。

目前已对白洋淀流域水生态与水环境等开展了大量研究，如杨薇等[2]对2000-2018年关于白洋淀生态需水量的研究成果进行了归纳总结；王欢欢等[3]对1988-2016年白洋淀水质时空演化特征进行了研究；白洁等[4]对白洋淀流域水环境承载力开展了定量评价，但对白洋淀与地下水的相互作用模式，尤其是其对区域地下水补给量等，尚未见到广泛报道。国内外对湖水与地下水转化关系也开展了大量研究[5-7]，研究人员通过物理、化学、数值模拟等多种方法对二者间的补给模式和补给量进行刻画[8]。如Wilson等[9]以氡同位素为示踪剂，结合遥感影像数据探究了爱尔兰劳夫湖与地下水的交互补给作用模式，研究结果表明地下水主要通过北侧和东侧的石灰岩补给湖泊，在西部的砂岩和混合花岗岩中并未发现补给通道；Yang等[10]根据热力学模型精确计算了地下水对云南阳宗海的补给量；Su等[11]根据同位素和热力学模型对内蒙古达克湖与地下水作用模式进行了研究，计算结果表明地下水对湖水的补给速率在春、夏、秋3季为2×10-6～3×10-6m/s，在冬季由于湖水-地下水潜流带结冰，阻隔了地下水的补给；Liao等[12]利用氢、氧稳定同位素和镭同位素识别出地下水对鄱阳湖的补给区；尹立河等[13]以湖水与地下水位监测数据为基础，结合水化学数据，对鄂尔多斯高原乞盖淖与地下水的动态补排关系进行了研究。上述研究均取得良好成果，但同时也面临以下问题：当镭、氡等同位素在地下水和地表水间未有显著差异时，对二者间的交互作用模式便难以精确刻画[9]；利用水位或热力学模型进行研究时，需投入较大工作量完成对水位或温度的长期监测，且研究成果多以定性为主[14]。基于此，结合雄安新区实际的水化学特征和土地利用类型，本研究以人工甜味剂(artificial sweetners，AS)为示踪剂，通过构建其在白洋淀中的时序特征，来探究典型淀区对地下水的补给量。

AS是一类广泛应用于食品、饮料、药物和个人护理品的人工合成或半合成的代替蔗糖的有机化合物[15]，由于大部分AS几乎不被人体转化，因此被称为无热量的糖。自20世纪末通过美国食品和药物管理局许可以来，全球每年的AS消耗量为15.9×104t，我国是AS第一大消费与制造国，约占全球消耗量的32%[16]。部分AS，如安赛蜜(acesulfame，ACE)、三氯蔗糖(sucralose，SUC)等具有较高的稳定性与持久性[17]，研究人员在大气粉尘[18]、地表水[19-21]、地下水[22]、土壤[23]、地表植被[24]等不同环境要素中均发现其广泛存在。我国对AS的研究尚处于起步阶段，目前以厘清其在各生态要素单元内的空间分布特征为主，如马晓雁等[17]对浙江省某市水处理及供水系统内ACE、糖精(saccharin，SAC)、阿斯巴甜(aspartame，ASP)、纽甜(neotame，NEO)的分布特征进行了研究；Gan等[24]对天津市大气粉尘、土壤、地表水中AS的空间分布进行了调查，对比干、湿季节下大气粉尘中AS含量的变化，其浓度范围介于0.02～1 940 pg/m3；Yang等[25]对珠江支流东江流域内地表水和地下水中AS含量的时空分布进行了研究，查明其主要来源于东江沿岸生活污水。但充分利用AS较强水溶性、较弱吸附性、单一来源等特性[18]，将其作为示踪剂探究人类活动的影响范围与强度、精确刻画地表水与地下水交互补给作用模式等方面的研究较少[26]。

基于此，本研究首先对雄安新区水环境中AS空间分布特征进行研究，查明人类活动对雄安新区水环境的影响；在此基础上通过构建AS在典型淀区内的时序特征，结合气象资料与补排特征，精确计算白洋淀内的金龙淀对地下水的补给量，为雄安新区水资源科学调控提供支撑。

2 研究区概况

雄安新区位于华北腹地，主要包括安新、雄县、容城三县及周边区域，面积约1 760 km2[1]，区域内地势平坦，新区土地利用分为农业种植区、生态林地、交通运输地、建成区、水体5种类型，其中农田面积最大，2019年遥感影像资料显示农业种植区面积占比达70.4%[27]，主要分布于新区西部和南部，作物以玉米和小麦为主。研究区属暖温带大陆性季风气候，具有夏季高温多雨、冬季干旱寒冷的特点，1965-2019年平均气温为7.3～12.7 ℃，年平均降水量为553.1 mm，70%集中于6-8 月份，降水年际变化大。区域内有华北地区最大的湖泊湿地白洋淀，具有完备的沼泽和水域生态系统，在区域气候、地下水补给、洪水调蓄以及生态系统稳定性等方面发挥着重要作用。白洋淀淀区总面积为366.4 km2，多年平均明水面积约占41.0%[28]，从北、南、西3个方向接纳9 条入淀河流。

研究区地势由西北向东南逐渐降低，地面海拔高程多在15～50 m，坡降小于2‰，第四纪沉积以冲洪积、冲湖积、冲积为主，地下水以孔隙水为主。根据华北平原含水层组划分原则，可分为浅层与深层地下水两大类，其中大部分地区浅层地下水水位埋深在5～15 m，单井涌水量在1 000～3 000 m3/d[1]，大气降雨和河湖入渗是其主要补给来源，地面蒸发、农业灌溉、生活用水是其主要排泄方式。

3 样品采集与测试

2019年7月在雄安新区采集42组浅层地下水和6组白洋淀地表水样品，测试ACE、SUC、SAC、甜蜜素(cyclamate，CYC)4种在我国最为常见的人工甜味剂[29]，所有样品均由自然资源部地下水科学与工程重点实验室测试完成，通过测试得出的雄安新区地下水位及人工甜味剂空间分布特征见图1。在此基础上，选择白洋淀南部的金龙淀作为试点，金龙淀紧临由中国地质科学院水文地质环境地质研究所建立的雄安新区地球多圈层交互带监测站[30]。为恢复白洋淀流域生态功能，河北省利用南水北调与引黄工程在每年3-7月对白洋淀实施生态补水，故在8-9月内，金龙淀补给来源仅为大气降雨，排泄渠道为蒸发与地下水渗漏，无外界水源对淀区内的AS含量造成干扰。本次研究从2019年8月1日起，每间隔7 d在金龙淀水域中心取地表水样品，截止至2019年9月26号，累积获取8组样品，结合不同AS半衰期，选择半衰期最长的ACE作为示踪剂[21]，构建其在金龙淀内的时序特征，计算该淀对地下水的补给量。

4 结果与分析

相较于其他天然示踪剂，AS最大的特点在于其具有较强极性、不易被土壤吸附、极易溶于水，且在自然界中无背景值[31]。目前AS被作为一种理想人工示踪剂，广泛用于定量刻画人类活动对生态系统的影响[32-34]。图1表明在雄安新区内，地下水中ACE含量为0.000 5～1.34 μg/L，平均值为0.11 μg/L，在白洋淀中其浓度为0.95～1.40 μg/L，平均值为1.16 μg/L；地下水中SUC浓度为0.003～0.93 μg/L，平均值为0.14 μg/L，在白洋淀中其浓度为1.14～3.16 μg/L，平均值为2.17 μg/L。SAC和CYC在地下水中浓度较低，SAC在地下水中浓度最大值为0.32 μg/L，平均值为0.04 μg/L,在白洋淀中其浓度最大值为0.27 μg/L，平均值为0.21 μg/L;CYC在地下水中浓度最大值为0.30 μg/L，平均值为0.05 μg/L,在白洋淀中其浓度最大值为0.50 μg/L,平均值为0.38 μg/L。因SAC和CYC在水中半衰期较短，且在自然条件下极易发生生物降解，故其含量远低于ACE和SUC[20]。

雄安新区水环境中的AS含量远高于德国[35]、瑞士[22]、加拿大[34]、新加坡[36]以及中国天津等地[37]，但与珠江流域东江河内的AS浓度相似，这充分表明人类活动对雄安新区与东江流域的影响强度较为相似，且远大于其他区域。从总体空间分布特征来看，上述4种AS在雄安新区西南部含量显著高于其他区域，这表明西南部区域受人类活动的影响更为强烈，结合新区目前的产业结构，人类活动主要为生活污水排放和农业灌溉，这些活动导致大量的AS进入白洋淀内，在白洋淀对地下水补给作用下，AS又进一步进入含水层；在研究区北部，地下水以山前侧向补给为主，地下水受白洋淀影响较小，故其AS浓度较低。

通过位于金龙淀南部的地球多圈层交互带监测站，采集到2019年8-9月金龙淀的日降雨和蒸发量及数据ACE时序分布特征，见图2。由图2可见，2019年8-9月金龙淀内降雨量较小，日降雨量为0～19.46 mm，累计降雨量为62.75 mm，蒸发量较大，日蒸发量为0.13～8.33 mm，累计蒸发量为235.8 mm，累计蒸发量远大于累计降雨量，在强烈蒸发和渗漏作用下，金龙淀内AS浓度持续升高。

图2 2019年8-9月金龙淀内降雨量、蒸发量及ACE浓度时序分布特征

由于ACE在水中的半衰期长达15 a[21]，其在短期内衰减量基本可忽略不计，ACE在淀内时空演化主要与大气降雨、蒸发、地下水渗漏有关, 金龙淀与浅层地下水作用模式示意图见图3。

图3 金龙淀与浅层地下水作用模式示意图

根据公式(1)～(5)即可反演出湖水对地下水的渗漏量。

(1)

(2)

(3)

(4)

(5)

式中：wi、wi-1分别为第i期、第i-1期金龙淀中ACE的浓度(i=1，2，…，8)，μg/L；v为湖水体积，m3；s为第i-1到i期间湖水净蒸发量，即蒸发量与降雨量之差，m3；z为湖水向地下水渗漏量，m3。

根据实地调查，金龙淀的水域面积为8.947 km2，结合净蒸发量数据和淀内ACE浓度的时序数据(图2)，根据公式(1)～(5)，确定金龙淀每7 d内平均向地下水渗漏量为5.25×104m3，则年渗漏量为273.75×104m3。

5 结 论

(1)研究结果表明，人类活动对雄安新区水环境的影响大于天津等地，尤其是在研究区南部和白洋淀内受人类活动影响最为严重，在后期的水生态环境保护和治理中需对该区域进行重点研究。ACE和SUC在地下水中的浓度平均值分别为0.11和0.14 μg/L，在白洋淀中的平均浓度分别为1.16和2.17 μg/L。在农业灌溉与生活污水的影响下，AS首先在地表水中富集，在其从湖水向地下水渗透过程中，在包气带内发生吸附与衰减，导致其在白洋淀内浓度远高于地下水。

(2)本研究以AS中ACE为示踪剂，构建其在地表水中的时序特征，在整个监测时间段内ACE浓度介于0.73～1.31 μg/L，在强烈蒸发和渗漏作用下，金龙淀内ACE浓度持续升高。结合地球多圈层交互带监测站降雨与蒸发数据，通过迭代计算，查明金龙淀向浅层地下水的渗漏量为0.75×104m3/d，则年渗漏总量为273.75×104m3。白洋淀是雄安新区浅层地下水系统主要补给来源，其对区域内地下水量和水质均有重要影响，精确刻画白洋淀与地下水的交互补给作用，对雄安新区水资源科学调控具有重要意义。