永定河北京段生态补水流量及模式.pdf
《永定河北京段生态补水流量及模式.pdf》由会员分享,可在线阅读,更多相关《永定河北京段生态补水流量及模式.pdf(10页珍藏版)》请在咨信网上搜索。
1、永定河北京段生态补水流量及模式徐佳琪1,吴悦1,仇文顺2,李述2,王梦瑶2,时洋1,魏加华1(1.清华大学水圈科学与水利工程全国重点实验室,水利部水圈科学重点实验室,北京,100084;2.北京市水资源调度管理事务中心,北京100089)摘要:以永定河北京段为例,建立水文学和水力学相结合的水流演进模型,利用生态补水实测监测数据率定并验证模型。利用建立的模型研究永定河北京段的生态补水规律及模式,结果表明:在持续稳定流量补水条件下,固安(出境断面)通水后 11d 左右达到稳定状态,停止放水后可维持河道 1213d 有水;在持续稳定流量下,随着补水流量的增大,河道首次贯通的时间和损失水量占比均降低,
2、且降低速率逐渐减缓,官厅补水下泄流量以 3035m3/s为宜;在生态补水总水量受限条件下,先以大流量下泄再调整为小流量的补水方式,对快速实现全线通水、河道渗漏回补地下水等效果更佳;平原段同步加入再生水和南水北调水,将进一步缩短全线通水时间;在官厅水库和平原段再生水、南水北调水向河道补水 2.24 亿4.14 亿 m3条件下,全线通水时间约需要 15d,固安断面出境水量约1.05 亿2.22 亿 m3,蒸发和渗漏水量 1.19 亿1.93 亿 m3,且补水水量越大,下渗和蒸发水量的占比越低;卢沟桥-六环路河段的入渗能力最强。本研究对永定河北京段生态补水水量、下泄过程及补水方式的确定具有实践意义。
3、关键词:永定河;生态补水;马斯京根;MIKE11HD;水量平衡中图分类号:TV213文献标志码:ADOI:10.13476/ki.nsbdqk.2023.0066为解决永定河河道断流、生态系统退化、地下水位下降严重等问题,水利部、海河水利委员会、北京市和河北省大力推进永定河流域的生态修复与环境治理,2019 年以来开展了多次生态补水工作,尤其是华北地区陆续开展了生态补水的尝试,如永定河引黄生态补水、滹沱河生态补水等1-2。魏健等3通过环境需水量和生态需水量计算,提出永定河官厅山峡段不同阶段河流生态修复的生态需水量,未来 3 个典型枯水年(对应来水频率 75%、90%和 95%)最低生态补水量分
4、别为 4.88 亿、6.11 亿和6.37 亿m3;邵慧芳等4基于水量平衡模型对比分析永定河北京段不同可配置水源、不同入渗条件下的生态补水效果,提出首次通水由官厅水库 4070m3/s大流量集中补水,首次通水后官厅水库 56m3/s、平原南段分段补水,全年 270d 不断流;孙冉等5基于河道水力学模型量化生态补水中的损失水量,提出流量介于 4080m3/s、同时采取小流量和大流量交替补水的建议。现阶段,我国华北和西北地区进行生态补水,主要是统筹安排一定的生态水量,利用工程调节补充河道因缺水而受损的生态系统。目前已开展的生态补水多为应急补水,规模和时间均有限,生态补水理论和调控方法研究不多。尤其
5、是华北平原,南水北调东中线工程正式通水以来,为受水区河湖复苏、地表-地下水协同修复提供了条件,但在水资源不足、生态需水得不到充分满足的条件下,如何合理确定生态补水过程、流量及模式,具有重要的理论和应用价值。本文探索开展基于河道水流贯通、地下水涵养和水生态系统指示物种用水需求下的补水流量及模式研究,深化生态补水规律认识。2019 年之前永定河北京段常年断流,实测径流时间序列有限,为收稿日期:2023-02-23修回日期:2023-07-16网络出版时间:2023-07-19网络出版地址:https:/ 3110km2,干流河道长 172km。自官厅水库起流经官厅山峡,进入北京市,于三家店出山峡进
6、入平原段,后流入河北省廊坊市。20192022 年,永定河北京段开展了 4 次生态补水(表 1)。2019 年春季首次补水为小流量“以水探路”;2020 年两次较大流量的集中补水分别实现了“以水开路”和“以水验路”;2022 年继续小流量、长历时补水,完成永定河连山通海目标。目前永定河生态补水还处于探索阶段。表1永定河历次补水情况Tab.1PreviouswaterreplenishmentsituationoftheYongdingRiver补水情况2019年春季2020年春季2020年秋季2022年春季补水时长/d(起止日期)86(3月14日6月7日)25(4月20日5月14日)19(10
7、月14日11月10日)76(4月10日6月15日)官厅水库累计下泄水量/亿m32.231.640.672.3官厅水库下泄流量变化过程/(m3s-1)10-25-30-4040-10040-80-34-2615-60-51-40水头到达位置黄良铁路桥下游800m天津武清永定河新河大桥下游400m河北廊坊固安G106国道上游屈家店枢纽注:“10-25-30-40”代表水库下泄流量从10m3/s变化到25m3/s、30m3/s和40m3/s,表中“40-100”“40-80-34-26”“15-60-51-40”含义相同。研究利用河道地形、生态补水逐日下泄流量、断面流量-水位等 20192022 年
8、的实测数据,主要流量断面包括官厅水库、雁翅(山峡段测站)、三家店(山峡与平原分界断面)、六环路(平原段测站)、河北固安(出境断面),固安断面作为下游水位边界,见图 1。233234234235235236236237020 40 60 80 100 x/m343638404244460150 300 450 600 750 x/mD-D:六环路252729313335370100 300 500 700z/mz/mz/mz/mz/mx/mE-E:固安57585960616263640 50 100150200250 x/mC-C:卢沟桥98100102104106108110050 10015
9、0200250300 x/mA-A:雁翅B-B:三家店AABBCCDDEE固安出流官厅水库入流南水北调水入流再生水入流(b)河道断面(a)永定河流域水系概化图N图例高程最大值:2 245最小值:39永定河流域研究断面市界干流支流水库图1永定河流域河道研究断面概化图Fig.1SchematicdiagramofwatersystemandtheresearchsectioninYongdingRiverbasin第21卷第4期南水北调与水利科技(中英文)2023年 8 月6702研究方法与模型2.1研究方法2.1.1水文学方法马斯京根(Muskingum)法7根据上游水情推求下游的河道流量,可满
10、足河道补水中上游水库放水后的下游断面河道流量过程演算的需要。洪水过程一般水量大、历时短,地下水补给和排泄可以忽略不计8,而生态补水流量小且历时长,尤其是长期干涸的河道,沿线地下水埋深大,补水过程中河水下渗回补地下水的占比较大,河道渗漏不可忽略。采用包含渗漏项的水量平衡公式,方程为dWdt=IQq(1)q=f I(2)W=kQ=kxI+(1x)Q(3)x=12l2L0.5 x 0.5式中:I、Q 分别为河段上游入流量和下游出流量,m3/s;q 为渗漏率,m3/s;f 为渗漏系数;Q为示储流量,m3/s;W 为河段的槽蓄量,m3;k 为槽蓄系数,h;x 为流量比重因子。其中,k 值等于在蓄水量 W
11、 下稳定流状态下的河段传播时间,x 的大小反映河段的调蓄作用大小以及水面线的形状。有学者9得到,其中 l 和 L 分别为河段的特征河长和河段长度,据此推断 x 的取值范围为10。将水量平衡方程和槽蓄方程联立求解即可得到马斯京根流量演算的方程为Q2=C0I2+C1I1+C2Q1(4)式中:I1、I2分别为计算时段始、末的入流量,m3/s;Q1、Q2分别为计算时段始、末的出流量,m3/s;t 为计算时段长,s。其中:C0=kx+(1 f)t/2k(1x)+t/2C1=kx+(1 f)t/2k(1x)+t/2C2=k(1x)t/2k(1x)+t/2(5)通过遗传算法率定参数 k、x、f 的值11,就
12、能够确定 C0、C1、C2,再代入流量演算方程,即可由上游河道断面流量演算得到下游河道断面流量。2.1.2水力学方法圣维南提出非恒定流偏微分方程,描述一维非恒定流运动规律。圣维南方程组由连续方程以及动量守恒方程组成12-14:BSht+Qx=qQt+gA(hx+Q|Q|K2)+x(aQ2A)=0(6)K=CAR式中:Bs为河道的断面宽度,m;h 为河道断面水位,m;Q 为河道断面流量,m3/s;t 为时间坐标,s;x 为空间坐标,m;q 为河道旁侧入流量,m3/s;A 为河道过水断面面积,m2;g 为重力加速度,m/s2;K 为河道断面流量模数,C 为河道断面谢才系数;R为断面水力半径,m;为
13、动量校正系数。MIKE11HD 模型利用六点中心 Abbott-Ionescu隐式差分法离散圣维南方程,利用“追赶法”求解差分方程,从而求解一维非恒定流的河道15-16。计算时水位 h 和流量 Q 在同一时间步长之下交替计算,因此计算水位 h 和计算流量 Q 无法同时在某点得到。水位点 h 布设在断面位置处,相邻的两个计算水位点之间布设一个计算流量点 Q,距离不一。该方法为无条件稳定,可以节省计算时间17。2.1.3模型精度评价方法选用均方根误差(ERMS)和纳什效率系数(ENS)评估模型精度18。ERMS用于描述观测值和模拟值时序数据间的偏差19,值越小表示模型性能越好;ENS用于评估模型的
14、预测能力20,值越接近最大值 1 模拟精度越高。ERMS和 ENS计算公式为ERMS=vuuuuutNi=1(PiOi)2N(7)ENS=1Ni=1(PiOi)2Ni=1(OiO)2(8)PiOiO式中:、分别代表模拟值和实测值;为实测值的平均值;N 为样本总数。参考水文情报预报规范21洪水预报精度评定采用的洪峰流量、洪峰出现时间、洪量和洪水过程等指标,生态补水水流演进过程用峰值流量、峰值出现时间和径流量三者的相对误差 评价模拟效果,计算见式(9),小于 20%时认为是合格预报;用流量过程确定性系数 R2表示预报和实测过程之间的吻合程度见式(10)。=PiOiOi(9)徐佳琪,等永定河北京段生
15、态补水流量及模式671R2=1Ni=1(PiOi)2Ni=1(OiO)2(10)2.2河道一维水文-水力学模型构建水文学方法模型简单且计算量小,流量仅取决于时间,忽略系统内部的不均匀性和河道水流运动规律22,即使没有河段的断面资料和粗糙度、比降等数据,水文学方法仍然适用。水力学方法中流量被看作是时间和空间的函数,可同时演算断面水位与流量,且能反映断面形状及糙率等河道特征的影响,演算结果更精确。因此,水文学和水力学方法之间存在准确性和复杂性之间的权衡23。从应用角度看,在缺少断面资料或对演算结果要求不严格的情况下适宜选用水文学方法,在需要同时演算流量和水位数据或对演算精度要求较高的情况下适宜选用
16、水力学方法。永定河北京段以三家店为界。上游为山峡段,河段地形较复杂,缺乏实测的断面资料,且系统内部的水流运动过程对本研究的影响不大,采用水文学模型模拟三家店流量作为下游模型的输入;水文学模型的输入数据为 20192022 年永定河北京段补水的实测流量数据,采用遗传算法率定。下游为平原段,河道宽度及断面面积更大,河道特征会影响流量,采用水力学模型模拟沿程流量及水位关系。建立永定河平原段的 MIKE11HD 模型,输入数据为 20192022 年补水实测流量、水位数据,和水文学模型输出的三家店流量数据。根据河道实测地形资料确定三家店至固安的河网文件、断面文件、边界文件以及参数文件,率定其糙率系数及
17、渗漏系数。2.3模型精度评价2.3.1马斯京根模型精度考虑到地层岩性导致的渗漏系数差异,将永定河山峡段划分官厅水库-雁翅和雁翅-三家店两个河段建立马斯京根模型。将官厅水库和雁翅 20192020 年补水期(共 338d)流量数据用作模型率定,20212022 年补水期流量数据用作模型验证(共148d)。将雁翅和三家店 20192021 年 11 月中旬补水期(共 186d)流量数据用作模型率定,2021 年11 月中旬2022 年补水期流量数据用作模型验证(共 62d)。官厅-雁翅段率定和验证时段的 ENS分别为 0.98 和 0.92,ERMS分别为 3.18 和 5.66;雁翅-三家店段率
18、定和验证时段的 ENS分别为 0.90 和 0.89,ERMS分别为 7.90 和 6.61。这是因为雁翅-三家店段的实测数据相对较少,导致模拟结果和实测数值之间的偏差更大,精度低于官厅-雁翅段。山峡段模型实测和模拟流量过程见图 2(a)。参数组合取值见表 2。0204060801001201316191121151181255565流量/(m3s1)补水天数/d实测流量模拟流量01020304050605 10152025303540455055 2 7 121722273237流量/(m3s1)补水天数/d实测流量模拟流量率定期验证期率定期 验证期(a)山峡段(b)平原段图2模型流量模拟效
19、果Fig.2Theflowsimulationeffectofthemodel表2马斯京根模型参数取值Tab.2ParametervaluesofMuskingummodel参数xkfC0C1C2官厅-雁翅0.4990.3270.0010.6690.3400.010雁翅-三家店0.0180.3840.1860.4640.4490.122以峰值流量、峰值出现时间、径流量和流量过程评价生态补水模型的模拟效果(表 3)。在两段模型的交集模拟时间内,整个山峡段模型的峰值流量 为 0.25,峰值出现时间均在同一天内;山峡段模拟的总出流量为 4.18 亿m3,实际总出流量为 3.90 亿m3,模拟出流水量
20、略大于实测,径流量 为 0.07,R2为0.90。认为山峡段模型峰值流量的模拟能力相对不足,对水量和流量过程模拟更精确。表3永定河水流演进模型模拟精度Tab.3SimulationaccuracyoftheYongdingRiverwaterflowevolutionmodel评价指标山峡段平原段峰值流量相对误差0.250.04峰值出现时间误差/d00径流量相对误差0.070.02流量过程确定性系数R20.900.93第21卷第4期南水北调与水利科技(中英文)2023年 8 月6722.3.2MIKE11HD 模型精度永定河平原段河道的入流水源除三家店外,还有再生水、南水北调中线工程向永定河平
21、原段补水(以下简称“南水”)。将三家店流量、再生水及南水补水流量作为研究河段入流流量,固安水位作为出流断面水位用于模型率定和验证。将 2021 年秋季补水期(共 58d)用作模型率定,2022 年春季(共36d)用作模型验证。率定时段的 ENS为 0.91,ERMS为 1.51;验证时段的 ENS为 0.67,ERMS为 4.98。糙率三家店-六环路河段为 0.20,六环路-固安河段为0.18。平原段率定及验证精度不如山峡段精度高,可能原因是实测水位数据量的限制和三家店拦河闸及卢沟桥拦河闸的人工闸门调控影响,导致出流断面的实际流量过程波动更为复杂。实测和模拟流量过程见图 2(b)。参数率定和验
22、证时段内,平原段模型的峰值流量 为0.04,峰值出现时间在同一天内;平原段模拟的总出流量为1.17 亿m3,实际总出流量为1.14 亿m3,模拟出流水量略大于实测,径流量 为 0.02,R2为0.93。平原段模型的预报与实测过程之间吻合程度较高,见表 3。整体上,水文学和水力学模型均能较为准确地模拟流量变化趋势及数值大小,可以准确反映补水过程及水量均衡要素,表明所建立的模型及其参数能够满足对生态补水流量过程的模拟需求。3生态补水规律研究利用建立的水流演进模型,探索生态补水的流量-水头演进时间、水量损失规律,并结合研究区域实际,提出适宜的补水方式,确定生态补水适宜水量和流量过程,为今后制定永定河
23、河道生态补水方案提供科学参考。3.1生态补水适宜流量补水流量与通水时间、稳定流量、损失水量等之间存在复杂的相关关系。为明确不同补水流量下的水头演进速度、损失水量占比等,模拟官厅水库不同下泄流量下的流量过程及均衡要素,绘制拟合曲线,进而选择时间和损失水量均在理想范围内的流量区间,获得适宜补水流量并推算补水周期(即一次补水自开始至终止放水所经历的时间)。水头演进时间为自官厅放水至固安断面有水流通过所用的天数(d)。损失水量包括渗漏水量、蒸发损失和槽蓄变化量,以渗漏水量为主,其中河道渗漏水量主要回补地下水,对地下水恢复有积极意义。在持续稳定流量下,通水后约 11d 达到稳定状态。停止放水后河道将维持
24、有水约 1213d,之后固安断面流量将减小为 0。流量与水头演进时间、达到稳定后日均损失水量占比之间的关系见图 3。随着流量的增大,通水时间和损失水量占比均降低,且降低速率逐渐减缓。因此,为将北京段全线通水时间控制在较小的天数(25d 以内),同时实现有效的地下水回补(损失占比大于 0.45%),将官厅补水流量控制在 3035m3/s 为宜。010203040506020406080水头演进时间/d补水流量/(m3s1)0.30.40.50.60.70.8020406080损失水量占比/%补水流量/(m3s1)(a)水头演进时间随补水流量变化(b)损失水量占比随补水流量变化图3不同补水流量的补
25、水效果Fig.3Waterreplenishmenteffectofdifferentreplenishmentflow补水起止时间主要依据最适宜物候、生境的时间确定。近年来永定河通过万家寨引黄和官厅、南水以及再生水等补水,河道已经形成了局部全年有水的水生生物栖息地。河道生态补水时机应尽量满足生物繁殖对物理环境的需求。根据孙涛等的研究,淡水鱼繁殖期为 46 月,其鱼卵为半浮性卵24,沉入水底过久会死亡,生态补水脉冲过程使鱼卵更不容易沉底,能够提升繁殖率。另外,鸟类繁殖高峰期大约在 57 月,34 月植物生长需水量旺盛。综合考虑动植物生长需求、河湖复苏、结冰期、汛前防洪库容预留等因素,选定生态补
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 永定河 北京 生态 水流 模式
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。