永定河（北京段）河道生态补水效益分析与方案评估

孙冉; 潘兴瑶; 王俊文; 任宇; 杜鹏; 马尧; 邢渊; Ran SUN; Xing-yao PAN; Jun-wen WANG; Yu REN; Peng DU; Yao MA; Yuan XING

中国农村水利水电 ›› 2021 ›› (6) : 19-24.

水环境与水生态

永定河（北京段）河道生态补水效益分析与方案评估

作者信息 +

An Analysis and Evaluation of Ecological Water Replenishment Benefit of Yongding River （Beijing Section）

Author information +

稿件信息 +

摘要

为了更加科学地制定北方断流河道的生态补水方案，选取常年断流河道永定河北京段为研究区，基于构建的河道水力学模型，模拟分析研究区2019年和2020年两次生态补水过程，量化了生态补水过程中蒸发和渗漏损失量约占20%。官厅水库补水流量越大，下游各断面峰值流量的削减率越小，为减小流量损失，提高补水流量利用率，在水源相对充足的条件下，建议补水流量介于40~80 m³/s之间，同时采取小流量和大流量交替补水的方法，将流域生态补水效益实现最优化。

Abstract

In order to formulate the ecological water replenishment plan of northern cut-off river more scientifically，this paper takes the Beijing section of the Yongding River as the research area. Based on the established channel hydraulics model， it simulates and analyzes the ecological water replenishment process of the study area in 2019 and 2020， and quantifies the amount of evaporation and leakage loss during the ecological replenishment process accounting for about 20%. The larger the replenishment flow of Guanting Reservoir， the lower the reduction rate of the peak flow of each downstream section. In order to reduce the flow loss and improve the utilization rate of the replenishment flow， it is suggested that the replenishment flow should be between 40~80 m³/s under the condition of relatively abundant water source， and the method of alternating small flow and large flow replenishment should be adopted at the same time. Under the condition of limited water source， the alternations of high-flow and low-flow water replenishment can be considered proper to realize the optimization of ecological water replenishment benefit in the basin.

关键词

河道断流 / 生态补水 / 水力学模型 / 效益评估 / 永定河

Key words

channel cut-off / ecological replenishment / hydraulics model / benefit assessment / Yongding River

基金

国家水专项课题项目(2018ZX07109-004)

国家自然科学基金项目(41730749)

引用本文

EndNote

Ris (Procite)

Bibtex

导出引用

孙冉 , 潘兴瑶 , 王俊文 , 任宇 , 杜鹏 , 马尧 , 邢渊. 永定河（北京段）河道生态补水效益分析与方案评估[J].中国农村水利水电, 2021(6): 19-24

Ran SUN , Xing-yao PAN , Jun-wen WANG , Yu REN , Peng DU , Yao MA , Yuan XING. An Analysis and Evaluation of Ecological Water Replenishment Benefit of Yongding River （Beijing Section）[J].China Rural Water and Hydropower, 2021(6): 19-24

0 引言

流域水循环过程改变一方面受气候变化等自然节律影响，同时还受水库、堤坝工程建设和水资源利用等人为因素影响，使得河道下游水量、水质、泥沙等过程发生改变^［1］。进入21世纪以来，我国北方很多流域气象、水文要素变化显著，降雨量明显减少、上游来水也持续减少，导致下游河道出现断流引发系列生态环境问题，如河道生态流量缺失、平原荒漠化、地表水与地下水补排关系失衡、甚至有些河道由于常年干枯断流，河床、河漫滩土壤干化、沙化，有些已经形成“沙龙”，生物环境遭到破坏，众多本底水生物种消失^［2，3］。

生态补水是解决断流河道水资源短缺的有效途径，科学的生态补水需要针对性的结合河道水文地理条件和水文地质条件，制定科学有效的补水方案^［4］。一般认为对断流河道进行生态补水需要科学确定断流河道所需的生态补水总量，通过补水过程分析，量化补水过程的主要损耗（包括蒸发量和渠道入渗量等）^［5，6］，以确保补水的生态效益最优。最终实现加强水资源节约、水环境保护和水生态修复的目标。

永定河是北京市重要的水源地，同时也是全国四大重点防洪江河之一，永定河作为“京西绿色生态走廊与城市西南的生态屏障”，对北京的生态环境品质提升和区域可持续发展意义重大。当前永定河河道生态用水长期得不到保障，地表干枯、河床沙化、植被退化，河流水陆生态系统严重退化，已严重影响了永定河生态服务功能的发挥，通过生态补水进行生态修复刻不容缓。水利部分别于2019和2020年启动了永定河生态补水工作，有效缓解了永定河流域缺水的紧张局面。

1 数据与方法

1.1 研究区概况

永定河上游是阴山和太行山支脉所包围的高原，有桑干河、洋河两大支流^［7］，总流域面积4.7 万km²，永定河流域多年平均降水量在370~650 mm之间。其中，永定河北京段全长169.7 km（官厅水库坝下至梁各庄出境），高程介于22.3~425 m之间。永定河北京段通常分为三段，其中从官厅水库至三家店为官厅山峡段，长度92 km^［8］，高差340 m。三家店至南六环定水桥为平原城市河段，河道长37 km；南六环至梁各庄出境段为平原郊野河段，长度为41 km。沿程主要建有：珠窝水库、落坡岭水库、三家店拦河闸、卢沟桥拦河闸、滞洪水库等（如图1所示）。

图1 永定河流域研究区示意图

Fig.1 Study area of Yongding River Basin

Full size|PPT slide

1.2 研究数据

本研究通过构建永定河一维河道模型定量模拟河道生态补水过程及其影响。模型覆盖清水河、永定河干流，模型构建所需数据主要包括空间地理基础数据和补水过程监测数据，其中空间地理基础数据包括：永定河北京段流域边界、数字高程图（DEM，Digital Elevation Model）、土壤类型图、土地利用图、河网图等（如表1所示）。补水过程监测数据主要结合永定河生态补水，选取沿线5个关键断面进行过程监测：官厅水库（坝下）、雁翅站、陇驾庄站、三家店拦河闸和卢沟桥拦河闸，监测指标主要包括水位、流量过程数据，以进行生态补水过程与效益分析。

表1 数据来源

Tab.1 Data sources

参数	来源	属性（备注）
流域边界	由30 m精度DEM提取得到	1 km×1 km网格
地形	由30 m精度DEM提取得到	1 km×1 km网格
河网水系	由30 m精度DEM提取得到	1 km×1 km网格
土壤类型	HWSD土壤数据库	Shape文件
植被类型	基于Landsat TM的土地利用图	网格文件
雨量站分布	泰森多边形法	1 km×1 km网格
降水	气象观测	13站
蒸发	气象观测+彭曼公式	13站
流量	水文观测	5站

1.3 模型构建

本研究通过构建永定河一维河道水动力学模型来模拟河道水流过程，一维水动力学模型可以快速、准确地模拟复杂河网的水位、流量过程，模型在无支流交汇河段采用圣维南方程组的动力波演算方法^［9］，基于垂向积分的水量和动量守恒方程^{［10，11］}。

质量守恒方程（连续方程）：

\frac{\partial Q}{\partial x} + \frac{\partial A}{\partial t} = q

（1）

动量方程：

\frac{\partial Q}{\partial t} + \frac{\partial}{\partial x} (\frac{Q^{2}}{A}) + g A \frac{\partial h}{\partial x} + g \frac{Q}{C^{2}} \frac{|Q|}{A R} = q \frac{\partial Q}{\partial x} + \frac{\partial A}{\partial t} = q

（2）

式中：A为过水断面面积，m²；Q为流量，m³/s；x为距离坐标，m；t为时间坐标，s；h为水位，m；q为旁侧入流流量，m³/s，；c为谢才系数；n为河床糙率系数，R为水力半径，m；g为重力加速度，m/s²。

本研究以实测2019年和2020年两次生态补水流量过程为模型输入。具体生态补水流量过程线如图2所示。

图2 生态补水流量过程线（官厅水库断面）

Fig.2 Process line of ecological replenishment flow （Guanting Reservoir section）

Full size|PPT slide

1.4 模型率定与验证

为了提升模型模拟精度，利用2019年（3月13日至6月16日）生态补水监测数据进行参数率定，利用2020年（4月20日至5月20日）生态补水数据进行模型验证，主要率定的参数为：河道入渗系数、河道糙率系数。

山峡段河道渗漏主要集中在下落坡岭水库至三家店段^［12］，河道入渗试验结果表明落坡岭至三家店河道平均入渗强度约0.2 m/d。根据水量平衡原理，基于河道监测断面实测的流量数据，逐段计算河道入渗能力，设置河道渗漏参数。根据永定河防洪规划的糙率取值，主河槽糙率n＝0.035，有杂草的滩地n=0.045，有树林滩地n=0.06，结合实测水位流量数据，采用试错法设置和率定糙率值。

本研究选取4个判断指标来量化模拟精度：绝对误差（用于表示模拟值和实测值的差值，该值越小表示模拟精度越高），相对误差（绝对误差与实测值的比值），纳什效率系数（NSE，Nash Sutcliffe Efficiency，取值范围为负无穷至1，值接近于1表示模拟精度越高；值接近于0表示模拟结果接近观测值的平均值），均方根误差（RMSE， root-mean-square error，也称作标准误差，用来衡量观测值同真值之间的偏差，能够很好地反映出测量的精密度）。

（1）参数率定。参数率定过程中，三家店实测水位波动较大，这是由于人为调度闸门所致。三家店拦河闸通过拦蓄上游来水形成一定的河道调蓄容积（设计水位107.51 m、校核水位109.35 m），调蓄能力最大为200 万m³。三家店拦河闸位于永定河出山口附近，具有重要的拦蓄洪水的作用。拦河闸上游左岸的进水闸就是永定河引水渠的起点，三家店拦河闸蓄水后可通过永定河引水渠向城区河湖生态供水，2019年生态补水期选取时段为3月13日至6月20日（图3），历时100 d，期间向城区补水4 726 万m³，最大补水量为10 m³/s。综合受到2019年生态补水期间上游来水量小（官厅出库最大仅为40 m³/s），同时兼顾三家店水库要维持一定水面面积、确保下泄量和为城区河湖生态供水等目标，管理单位对三家店拦河闸进行频繁调度，水位波动受人为影响波动相对明显，介于106.66~108.19 m之间。三家店断面是山峡段的下游出口边界，且山峡段地形落差高达336.2 m，因此，三家店断面水位波动对山峡段整体计算的水位、流量结果影响都较小。经过参数率定，模拟水位和实测水位一致性较好，模型可以用于补水过程分析。

图3 2019年生态补水模拟水位定结果

Fig.3 Simulation results of water level for ecological replenishment in 2019

Full size|PPT slide

表2 2019年生态补水模型计算水位与实测对比

Tab.2 Comparison between calculated water level and measured water level of ecological replenishment model in 2019

断面

计算水位值

（均值）/m

实测水位值

（均值）/m

绝对误差/m

相对误差/%

NSE

RMSE

雁翅

233.892

233.894

-0.002

0.001

0.863

0.039

陇驾庄

134.464

134.484

-0.020

0.015

0.782

0.120

三家店拦河闸

107.570

107.910

-0.340

0.320

0.580

0.300

永定河（北京段）自1993年来基本没有洪水，河床沙粒径变化不大，河床质为砂卵石，跟永定河目前河床形态及植被状况，糙率介于0.03左右。构建的山峡段河道模型为一维河道模型，糙率系数在合理取值范围内主要通过影响河水流速进而影响流量，入渗系数主要影响河道渗漏补给量。参数率定结果表明，2个参数分别在合理取值范围时对模拟结果相对敏感，当取值超出合理范围后模拟结果变化不大。基于2019年生态补水过程数据，率定的河道综合糙率系数在合理取值范围内，介于0.015~0.05之间，并分段确定了河道的入渗系数（如表3所示）。

表3 河道分段糙率赋值表

Tab.3 Roughness assignment of channel segments

河流分段	糙率系数	入渗系数/（m $∙$ d^-1）
官厅水库-珠窝水库	0.035	0.017
珠窝水库-雁翅站	0.029	0.302
雁翅站-落坡岭水库	0.020	0.029
落坡岭水库-下苇甸	0.015	0.051
下苇甸-三家店拦河闸	0.050	0.051
三家店拦河闸-卢沟桥拦河闸	0.045	0.010
卢沟桥拦河闸-黄良铁路	0.025	0.017

（2）模型验证。基于2020年生态补水实测数据（图2）对模型进行验证，如图4和表4。可见雁翅站、陇驾庄站、三家店模型参数质量较好，模型可信度较高。2020年度生态补水期间，为了实现在下游开展脉冲实验蓄水的目标，要求三家店水库采用“平进平出”的调度模式。因此，补水期间三家店水库水位相对保持稳定，此次生态补水期历时较短，人为调度影响较弱，模型模拟结果与实测结果的一致性更好。

图4 2020年生态补水模型水位验证

Fig.4 Water level verification of ecological replenishment model in 2020

Full size|PPT slide

表4 2020年生态补水模型计算水位与实测对比

Tab.4 Comparison between calculated water level and measuredwater level of ecological replenishment model in 2020

断面

计算水位值（均值）/m

实测水位值（均值）/m

绝对误差/

相对误差/

NSE

RMSE

雁翅

234.080

234.140

-0.06

0.026

0.622

0.21

陇驾庄

134.653

134.797

-0.144

0.110

0.763

0.22

三家店拦河闸

108.050

108.106

-0.056

0.052

0.778

0.12

2 结果分析

2.1 2019年生态补水过程及其水量平衡分析

模拟结果表明2019年3月14日至6月6日永定河山峡区间春季补水总量为2.2 亿m³。官厅水库至雁翅段长度为59 km，损失水量是1 313.21 万m³，雁翅至陇驾庄段长度为30 km，损失水量737.05 万m³，陇驾庄至三家店拦河闸段长度为11.5 km，损失水量2 025.46 万m³，合计渗漏损失水量为4 075.72 万m³，河段综合损失率为18.5%。与实测渗漏损失水量4 399.7 万m³相差323.98万m³，则山峡段补水损失相对误差为7.4%。

模型计算各关键断面的流量过程如图5所示，官厅、雁翅、陇驾庄和三家店峰值流量分别为40、37.13、34.44和35.80 m³/s，官厅至雁翅、雁翅至陇驾庄、陇驾庄至三家店峰值削减率分别为7.2%、7.3%、9.2%，峰现时间雁翅、陇驾庄、三家店分别比官厅水库延迟36、42、54 h。

图5 2019年各监测断面流量过程

Fig.5 Flow process of each monitored section in 2019

Full size|PPT slide

2.2 2020年生态补水过程及其水量平衡分析

模拟结果表明2020年4月20日至5月15日永定河山峡区间春季补水总量1 6264.28 万m³。官厅水库至雁翅、雁翅至陇驾庄、陇驾庄至三家店拦河闸的损失水量分别为349.35万、630.38万和2 176.9 万m³。模型计算得出官厅水库至三家店拦河闸渗漏水量为3 156.63 万m³，河段综合损失率为19.4%。与实测渗漏水量3 630.01 万m³相差473.4 万m³，则山峡段补水损失相对误差为13%。

模型计算各关键断面的流量过程如图6所示，官厅、雁翅、陇驾庄和三家店峰值流量分别为100、96.26、94.23和79.89 m³/s，官厅至雁翅、雁翅至陇驾庄、陇驾庄至三家店峰值削减率分别为3.74%，2.11%，15.2%，峰现时间雁翅、陇驾庄、三家店分别比官厅水库延迟28、40、45 h。

图6 2020年各监测断面流量过程

Fig.6 Flow process of each monitored section in 2020

Full size|PPT slide

2.3 不同生态补水效益对比分析

基于模型分别模拟官厅水库分别以5，10，20，40，80，100，120和150 m³/s等8种出库情景，对比分析下游各断面峰值流量和和达到峰值所需时间（图7，图8）。

图7 不同官厅水库补水情景下游各断面峰值流量

Fig.7 Peak flow of each section downstream of different Guanting reservoirs under different water replenishment scenarios

Full size|PPT slide

图8 不同情景下游各断面达到峰值流量所需时间

Fig.8 The time required for each section downstream of different scenarios to reach the peak flow

Full size|PPT slide

官厅水库补水流量越大，下游各断面峰值流量的削减率越小。当补水量为5 m³/s时，三家店断面峰值流量为2.18 m³/s，流量削减率达到56%。当补水量达到150 m³/s时，三家店断面峰值流量为130.1 m³/s，流量削减率达到13%。在情景5当官厅水库补水流量为80 m³/s时，下游断面流量削减率也维持在较低水平达到15%，因此，在水量水源相对充足的条件下，为了减少沿途流量损失补水流量适宜维持在80 m³/s。

官厅水库下游各断面流量达到峰值流量所需时间随着补水流量的增大而缩短，不同官厅水库生态补水情景下三家店达到峰值流量所需时间介于33~261 h之间。水量传输主要受传输距离和高程落差影响，雁翅、陇驾庄、三家店与官厅水库之间的距离分别为63、97.4和108.3 km，高差分别为205.6、305.2和336.2 m。在情景1补水量为5 m³/s时，官厅水库到三家店的水量传输时间主要消耗在官厅水库到雁翅之间约占60%。当补水流量大于10 m³/s时，传输时间主要消耗在陇驾庄至三家店之间。当补水流量介于5~40 m³/s之间，陇驾庄至三家店之间随着流量增加水量传输时间显著减少。当流量大于80 m³/s时，3处断面达到峰值流量所需时间相近。因此，从水量传输时间效率分析，在水源相对充足的条件下，补水流量介于40~80 m³/s之间效率最高。

2.4 补水方案优化与建议

根据情景1模拟结果，当官厅水库补水量为5 m³/s时，三家店断面峰值流量为2.18 m³/s，则官厅水库至三家店河道入渗能力约为3 m³/s，进一步模拟结果表明当官厅水库以3 m³/s下泄进行为期20 d的生态补水时，恰好无水量到达三家店，生态补水全部用于官厅水库至三家店全段渠道渗漏损失。表明在补水期官厅水库下泄流量保持不小于10 m³/s、非补水期官厅水库下泄流量保持不小于3 m³/s进行精准补水，这是恢复流域水生态的最低补水阈值。

生态补水流量越小，流量损失率越大，且随着补水流量的减小，流量损失率的增幅越大。因此，生态补水时，为减小流量损失，提高补水流量利用率，在水源相对充足的条件下，建议补水流量介于40~80 m³/s之间。在水源有限条件下，可适当考虑大流量补水和小流量补水交替进行，可实现减少过程损失，同时还能实现全线贯通的目的。

此外，干旱断流地区蒸发量较大，可以考虑中水回用、人工补水等非常规水源。通过景观营造强化河流的节水和保水措施，降低补水用水过程中的消耗，也是维系河流水文平衡的重要方法，如永定河的“五湖一线”工程^［15］等，如对河床采用防渗膜或防渗墙等做防渗处理，尽量减少渗透量，达到蓄存河道水量的效果。

3 结语

北方地区河道断流现象普遍，带来一系列生态环境问题，实施生态补水是恢复河道基本生态功能的有效途径。本研究永定河流域北京段为研究区，通过构建河道一维水动力学模型，基于实测2019年和2020年生态补水数据资料，对模型进行验证，并基于模型分析了永定河两次生态补水的过程和实际效益。

结果表明生态补水是解决断流河道水资源短缺的有效途径，补水过程中蒸发和渗漏损失量约占20%。官厅水库补水流量越大，下游各断面峰值流量的削减率越小，当官厅水库补水量为5 m³/s时，三家店峰值流量削减率为56%，如果补水量小于3 m³/s时，则三家店流量为0。生态补水时，若水源相对充足，建议补水流量介于40~80 m³/s之间；若水源有限时，可适当考虑大流量补水和小流量补水交替进行。

河道的生态环境恢复难度大影响因素复杂，需要合理采取科学的补水措施进行生态修复，逐步恢复河流自然的水文过程和自然功能，形成自然的河漫滩和河岸系统，改善断流河道地区常年缺水、河道萎缩的现状。

参考文献

原文顺序 | 文献年度倒序 | 文中引用次数倒序

1	崔建军，郑振华，张韬.河北省太行山区水库下游河道生态恶化特征及成因分析：以滹沱河下游河道为例［J］.中国集体经济，2015，21（3）：59-60. 本文引用 [1]

2	纪玉琨，武春侠，朱毕生.新时期永定河流域生态修复思考［J］.水利发展研究，2018（11）：31-34. 本文引用 [1]

3	李昂，刘定湘.河湖长制立法完善思考［J］.水利发展研究，2020（6）：9-11，31. 本文引用 [1]

4	徐震，赵进勇，李庆国，等.玉符河干流生态适宜流量与生态补水效应分析［J］.中国农村水利水电，2018（6）：79-83. 本文引用 [1]

5	杨卫，许明祥，李瑞清，等.面向生态环境的河湖连通引水调控方案研究［J］.武汉大学学报（工学版），2020，53（10）：861-870. 本文引用 [1]

6	李计生，于海超，张扬，等.基于Tennant改进法的讨赖河生态流量分析与评价［J］.节水灌溉，2020（8）：72-75. 本文引用 [1]

7	杨柠.永定河引黄生态补水长效机制初步探索［J］.水利发展研究，2020（2）：13-16. 本文引用 [1]

8	张琳.民国以来永定河流域北京段生态环境变迁研究［D］.北京：北京林业大学，2015. 本文引用 [1]

9	刘家宏，梅超，向晨瑶，等.城市水文模型原理［J］.水利水电技术，2017（5）：1-5，13. 本文引用 [1]

10	卢慧，董红霞，轩晓博.基于减水河MIKE11建库后对水环境的影响［J］.水科学与工程技术，2018，18（6）：20-22. 本文引用 [1]

11	潘兴瑶，李其军，陈建刚，等.城市地区流域洪水过程模拟：以清河为例［J］.水力发电学报，2015，34（6）：71-80. 本文引用 [1]

12	杨勇，谷健芬，李元春.永定河长时放水条件下的入渗能力及影响范围分析［J］.工程勘察，2019，37（6）：29-35. 本文引用 [1]

13	刘易庄，杨树清，蒋昌波.洞庭湖萎缩对湖内洪水影响［J］.水科学进展，2019，30（4）：485-495.

14	史源，章少辉，白美健，等.欧拉法融合拉格朗日法高效模拟灌溉二维地表水运动规律［J］.农业工程学报，2019，35（12）：150-156.

15	张敏秋.北京市永定河“五湖一线”综合治理思路与经验启示［J］.水利规划与设计，2016，17（2）：7-8，20. 本文引用 [1]

PDF(2573 KB)

1033

访问

引用

详细情况

段落导航

摘要
Abstract
关键词
Key words
基金
引用本文
0 引言
1 数据与方法
1.1 研究区概况
图1 永定河流域研究区示意图
1.2 研究数据
表1 数据来源
1.3 模型构建
图2 生态补水流量过程线（官厅水库断面）
1.4 模型率定与验证
图3 2019年生态补水模拟水位定结果
表2 2019年生态补水模型计算水位与实测对比
表3 河道分段糙率赋值表
图4 2020年生态补水模型水位验证
表4 2020年生态补水模型计算水位与实测对比
2 结果分析
2.1 2019年生态补水过程及其水量平衡分析
图5 2019年各监测断面流量过程
2.2 2020年生态补水过程及其水量平衡分析
图6 2020年各监测断面流量过程
2.3 不同生态补水效益对比分析
图7 不同官厅水库补水情景下游各断面峰值流量
图8 不同情景下游各断面达到峰值流量所需时间
2.4 补水方案优化与建议
3 结语
参考文献

收稿日期	出版日期
2020-08-14	2021-06-15
发布日期
2021-07-02

选择文件类型/文献管理软件名称

选择包含的内容