生态补水对白洋淀流域地表水和地下水水化学特征的影响*

崔 旭，张 兵，何明霞，夏文雪，王义东，赵 勇

(1:天津师范大学天津市水资源与水环境重点实验室，天津 300387)(2:天津师范大学地理与环境科学学院，天津 300387)(3:中国水利水电科学研究院水资源研究所，北京 100038)

湿地被誉为“地球之肾”，它是由水陆生态系统相互作形成的生态系统，在维系地区生态平衡，推动经济社会发展中具有不可替代的地位[1-3]. 近年来，在人类活动与全球气候变化的综合影响下，湿地面临着生态环境恶化、面积持续萎缩以及水体污染严重等危机[4]. 充足的水源是保障湿地水环境及生态服务功能的基础. 跨流域调水可以改善水资源时空分布不均，缓解区域水资源供需矛盾[5]. 通过跨流域调水向湿地生态水网补水，可以修复、改善其受损的生态系统结构与功能，使湿地自我调节能力得到恢复，对维持湿地生态系统生态平衡起到至关重要的作用.

人类活动会对湿地水环境产生深远影响. 目前，人类工、农业活动对白洋淀水生态产生影响的研究已较为深入，主要包括：水质时空变化[6-7]、湿地景观格局变化[8]、水体富营养化监测[9-10]、污染物的特征与溯源[11-12]以及水生生物对人类活动的响应[13]等方面. 水化学是水体的基本要素之一，分析水体水化学特征可以表征其环境化学特点、水质状况、水体离子来源、分布及迁移转化规律. 生态补水在增加受水区地表水水量、涵养和修复地下水资源、提升地表水与地下水水位[14-15]、改善淀区水动力状况[16]和水质[17]的同时，也会使区域地表水与地下水的水化学类型发生改变[18-19]. 研究生态补水条件下水体水化学变化，可为受水区水环境及生态的修复或进行合理的农业灌溉提供理论支撑.

白洋淀是华北平原最大的淡水湿地生态系统. 自1960年以来，气候、水文等自然因素，以及工农业与生活用水量加大、污水排放、上游建设水利工程设施以及地下水的过度开采等人为因素，导致区域水环境逐渐恶化[6,20]. 白洋淀湿地景观破碎化、干化问题严重[21]，水体富营养化和沼泽化问题并存，污染物总量居高不下，水质长期处于Ⅳ～劣Ⅴ类之间[22-23]. 2017年4月1日，中共中央、国务院设立雄安新区，白洋淀在其建设过程中具有调节区域气候、优化地区生态环境的重要功能. “建设雄安新区，一定要把白洋淀修复好、保护好”. 水是白洋淀生态之基，向淀区补水以恢复白洋淀“华北之肾”功能的任务迫在眉睫. 截至2020年，共向白洋淀补水约197.6×108m3. 近年来国内外学者对白洋淀生态补水量[24-25]、配水水源路径[26]及补水对区域水动力的影响[16]等方面进行了大量研究，也对补水方案产生的生态效益进行了评价[27-28]. 作为生态补水的典型受水区域，有关多水源长效补水机制对白洋淀地表水与地下水水化学组成的影响方面的研究仍需深入.

2018年4月16日-6月22日，“南水北调”中线工程首次通过瀑河及北易水退水闸向白洋淀进行生态补水，补水量为1×108m3. 同年4月5日-6月5日，王快水库和西大洋水库向白洋淀补水2.8×107m3. 本文通过野外采样和室内测试，分析2018年白洋淀淀水、河水及地下水生态补水前、后水化学时空分布特征及离子来源，阐释生态补水对地表水和地下水水化学特征的影响机理，探讨生态补水对未来白洋淀水环境可能产生的影响，以期为白洋淀水体生态环境改善及雄安新区生态核心区的保护与治理提供科学依据.

1 研究区概况

白洋淀(38°44′～38°59′N，115°45′～116°07′E)地处河北省雄安新区，总面积为366 km2，由143个淀泊与沟壕组成，是华北平原最大的淡水湿地[20](图1a). 流域地处暖温带半湿润大陆季风气候区，夏季炎热多雨，冬季寒冷干燥，多年年平均气温为12.8℃，多年年平均降水量为524.8 mm，多年年平均水面蒸发量达1369 mm[21]. 气温、降水量与蒸发量年内变化较大，夏季(6-9月)月均温度均大于20℃，其中7月平均温度高达26.9℃，约67%的降雨集中于6-9月. 流域地势呈现西北高、东南低的特征. 地处海河流域大清河水系九河下梢的白洋淀，承接南支潴龙河、孝义河、唐河，西侧府河、漕河、瀑河、萍河以及北侧的白沟引河水，并经由枣林庄枢纽排入赵王新河，在当地具有蓄水行洪的重要作用(图1a).

图1 白洋淀流域水系(a)及淀区补水线路、采样点分布(b)Fig.1 River system of Lake Baiyangdian(a), hydrating line and distribution of sampling sites(b)

以水文地质要素为依据，将该地区4个含水层组划分为浅层地下水系统和深层地下水系统(图2). 浅层地下水主要存在于第Ⅰ和第Ⅱ含水层组内，为全新统(Qh)和上更新统(Qp3)含水层. 第Ⅰ、Ⅱ含水层组属潜水和浅层承压水，水位埋深为5～20 m，含水层厚度小于40 m，底界深度为150～180 m，以中、细、粉砂为主. 浅层地下水受大气降水、农业灌溉、地表水入渗及侧向径流补给. 深层地下水主要赋存于第Ⅲ和第Ⅳ含水层组内，为中更新统(Qp2)和下更新统(Qp1)含水层. 第Ⅲ和第Ⅳ含水层组主要为承压水，水位埋深为25～35 m，含水层厚度为100～180 m，低界深度为350～600 m，岩性以粗、中、细砂为主，主要受侧向径流补给[29].

图2 研究区水文地质剖面[29]及地下水水化学组成Fig.2 Hydrogeology profile in the study area[29] and hydrochemistry of groundwater

1960年以来，受气候变化与人类活动的影响，入淀水量逐年减少，除白沟引河、府河及孝义河外的大部分入淀河流基本断流，白洋淀由天然过水型湖泊逐渐变为人工调蓄型湖泊[30]. 与此同时，水位下降、水体污染、湿地萎缩、淀区沼泽化以及生物多样性退化等一系列问题亟待解决. 2017年雄安新区建立，作为区域生态屏障，白洋淀在雄安新区建设中有“以淀兴城、城淀共荣”的作用. 为保障区域生态水源，1981-2018年政府相继采取了上游水库补水(王快水库、西大洋水库、安格庄水库、岳城水库)、“引黄济淀”、“引岳济淀”、“南水北调”中线工程补水等一系列措施(图3). 上游水库补水及跨流域调水的多水源补水措施，使白洋淀淀区水质逐步优化，湿地面积增加，区域水环境得到进一步改善，遏制了白洋淀生态环境的进一步退化.

图3 白洋淀生态补水信息Fig.3 Information of ecological water replenishment to the Lake Baiyangdian

补水线路如图1b所示. 王快水库水流经沙河总干渠、月明河、孝义河向白洋淀补水；西大洋水库水通过唐河总干渠流至白洋淀；安格庄水库由中易水、南拒马河、白沟引河向白洋淀补水. 2004年3月1日至6月29日“引岳济淀”工程以岳城水库为起点，利用滏阳河、老漳河、滏东排河、小白河等现有河渠向白洋淀补水1.59×108m3[31]. 2006-2011年自“引黄济淀”工程向白洋淀补水5次，从黄河位山闸引水，经位山三干渠、穿卫工程，沿清临干渠、清凉江、滏东排河、小白河、人文干渠进入淀区，补水量为4.52×108m3[27]. 2012年5月30日，“两库(王快水库-西大洋水库)连通”工程竣工，王快水库水跨沙河流域至唐河流域，经三会河流入西大洋水库，再经保定市区“大水系输水”工程，由白洋淀应急输水工程及府河向白洋淀补水. 2017年11月16日“引黄入冀”补淀工程自卫河、老漳河、滏东排河、滹沱河、古洋河、小白河和任文干渠输水至白洋淀. 2018年4月16日“南水北调”中线工程首次通过瀑河退水闸、北易水退水闸、郑家佐分水口，向白洋淀上游河道及淀区补水，入淀水量达1×108m3.

2 材料与方法

2.1 样品采集与测试

2.2 数据分析

对各离子浓度进行描述性统计分析，运用ArcGIS软件绘制白洋淀流域水系、淀区补水线路、采样点分布及水体EC值图；用 Microsoft Excel软件绘制Stiff图，水文地质剖面图采用CorelDRAW软件绘制. 利用AquaChem软件制作Piper图以分析水体水化学类型. 为明晰水体离子来源，应用Origin软件绘制Gibbs模型图.

3 结果与分析

3.1 水化学特征

表1 2018年各水体补水前后水化学指标Tab.1 Physicochemical parameters of water samples before and after ecological water replenishment in 2018

3.2 水化学类型

地表水水化学类型时空变异性较大. 补水前淀水水化学类型主要为Na-Cl·SO4、Na-Cl·HCO3，补水后淀水水化学类型种类繁多，主要为Na·Mg-HCO3·Cl、Na·Ca-HCO3·Cl、Ca·Na-HCO3·Cl、Ca·Mg-HCO3·SO4. 补水前后河水水化学类型差异较大. 补水前府河沿程流向水化学类型为：Na-Cl·SO4、Na-Cl·HCO3、Na-Cl，瀑河水化学类型为Na-HCO3·Cl，白沟引河河水化学类型为Na-Cl·SO4，赵王新河水化学类型为Na-SO4·Cl；补水后府河、瀑河、白沟引河以及南河干渠的水化学类型分别为Na·Ca-Cl·HCO3、Mg·Ca-HCO3、Ca·Mg-HCO3·SO4和Na·Ca-HCO3.

补水前后深层地下水水化学类型较稳定，浅层地下水水化学类型变化较复杂. 淀区及周边深层地下水水化学类型以Na-HCO3为主，补水后小部分深层地下水水化学类型发生改变. 位于瀑河周边的深层地下水G14w水化学类型为Na·Mg-HCO3，G3w水化学类型为Na·Mg-HCO3·Cl，G11w水化学类型为Na-HCO3·Cl. 浅层地下水水化学类型变化显著. 补水前位于唐河附近的浅层地下水水化学类型为Na-SO4，补水后浅层地下水水化学类型为Na·Ca-HCO3·Cl、Na·Mg-HCO3、Mg·Na-Cl·HCO3.

位于瀑河入淀口的深层地下水G14w(400 m)与浅层地下水G15w(100 m)的水化学类型在瀑河河水采样点R4w与R6d之间，且浅层地下水G15w(100 m)水化学类型更接近于补水后瀑河水水化学类型R4w；位于白沟引河入淀口的地下水G1w(100 m)水化学类型在R1w与R7d之间且与R1w水化学类型基本一致；府河入淀口的地下水G16w(50 m)水化学类型为Na·Ca-HCO3·Cl与府河河水R6w水化学类型Na·Ca-Cl·HCO3较为接近；地表水(R1w、R4w、R6w)与其附近的浅层地下水(G1w、G15w、G16w)水化学类型类似. 这些现象都表明白洋淀地表水与浅层地下水存在密切的水力联系，而深层地下水与浅层地下水水力联系较弱. 补水后地表水与地下水水化学类型位于补水前地表水、太行山区地下水及水库水水化学类型之间(图4)，表明补水改变了白洋淀水体水化学类型，且地下水受地表水、太行山区侧向径流的共同补给.

图4 白洋淀不同时期各水体Piper图Fig.4 Piper plots of different waters of Lake Baiyangdian in different periods

3.3 水化学组分影响因素

图5 白洋淀湿地不同时期各水体Gibbs图Fig.5 Gibbs diagram of different waters of Lake Baiyangdian in different periods

4 讨论

4.1 生态补水对地表水环境的影响

生态补水增加了入淀水量，使淀区水位上升，水面积扩大，水质得到一定改善. 1965年前，白洋淀主要由降水补给，水量丰富，1957年水位高达9.5 m；1965年后，上游水库的修建及区域用水量的增加，导致入淀水量减少，1983-1988年出现连续5年干淀的现象[20]，1990-2000年，白洋淀水质呈下降趋势[23]. 为化解干淀危机，改善淀区水质，遏制白洋淀生态退化，政府多次采取生态补水及淀区污染治理措施. 1988-1996年，白洋淀水位大幅度上升，2007年后水质状况明显改善(图6). 但在补水过程中仍旧存在诸多问题，如：河道渗漏蒸发，沿岸居民及工厂非法取水、排放污染物，引水工程存在安全隐患且缺乏较为完善的长效补水机制[32]. 1997-2002年白洋淀水位持续下降，2000-2008年再次出现干淀现象，白洋淀“水危机”并没有从根本上得到解决.

图6 1950-2020年白洋淀水位、水面积与水质等级[21,23] Fig.6 Water level, water area and water quality grade from 1950 to 2020[21,23] in Lake Baiyangdian

2018年4月发布的《河北雄安新区规划纲要》明确指出，“建立多水源补水机制，坚持控源-截污-治河系统治理”. 通过改进补水方案、整治排污企业、优化淀中村废水处理体系以及规范水产养殖，淀区水环境好转，水质逐步优化[33]. 截至2019年底，淀区水位始终保持在6.5 m以上，水面积增长至约250 km2[34](图6). 2019年颁布的《白洋淀生态环境治理和保护规划(2018-2035年)》提出：到2022年，淀区生态环境质量将初步恢复，2035年综合治理全面完成，淀区水位将维持在6.5～7 m，面积将扩大至360 km2，水质将达到国家地表水环境质量Ⅲ～Ⅳ类标准[34]. 白洋淀将逐步恢复其“华北之肾”的功能，维护雄安新区乃至华北平原的生态平衡.

图7 补水前、后地表水的水化学特征及EC值Fig.7 The hydrochemical characteristics and EC value of surface water before and after ecological water replenishment

4.2 生态补水对地下水环境的影响

白洋淀地区地下水与地表水存在着密切的水力联系[38]，地表水渗漏补给地下水贡献率高达52%[39]，周边地区浅层地下水平均更新速率为4.4%/a[40]. 近几十年来，随着经济社会的发展，区域用水量不断增加，白洋淀地下水超采现象严重，地下水位大幅下降. 1974-2007年, 白洋淀流域地下水埋深以平均每年0.51 m的速率下降, 在淀区附近形成地下水降落漏斗群[41]. 在2000年后，随着补水量不断加大，淀区周边地下水资源涵养修复效果显著，地下水降落漏斗群面积减小，水位有所上升[14,42]. 2015-2018年间，白洋淀流域大部分地区地下水储量呈上升趋势[43]. 预计到2030年，白洋淀地区浅层地下水位将较2016年上升20～30 m，淀区南部或将上升35 m[44].

表2 浅层地下水与深层地下水不同时期水化学指标Tab.2 Physicochemical parameters of water samples in the surface and deep groundwater

5 结论

采集2018年白洋淀枯、丰水期(补水前、补水后)地表水与地下水水样，分析水体水化学组成，主要结论为：

2) 深层地下水与浅层地下水存在不透水层，因此受生态补水影响较小，水化学类型以Na-HCO3为主. 枯水期(补水前)水体离子来源受蒸发结晶作用影响较大，丰水期(补水后)岩石风化作用影响更为强烈.

水化学是反映湿地水环境的重要计量指标，本研究主要依据2018年生态补水前后及2010年数据进行解析，未来将融合多学科、新方法，增加监测时间、监测点及生态监测数据，更全面地分析生态补水后由地下水位上升、水体在含水层中发生阳离子交换和水岩相互作用可能引发的一系列生态环境问题.