荒漠区平原水库多目标生态调度能力与潜力分析.pdf

D O I:1 0.1 9 3 2 2/j.c n k i.i s s n.1 0 0 6-4 7 1 0.2 0 2 3.0 1.0 1 0h t t p s:/x u e b a o.x a u t.e d u.c n.引文格式:牟琳辉,白涛,加沙尔蒙古拜,梁宇航.荒漠区平原水库多目标生态调度能力与潜力分析J.西安理工大学学报,2 0 2 3,3 9(1):9 6-1 0 4.M O UL i n h u i,B A IT a o,M E N G G U B A I J i a s h a e r,L I A N GY u h a n g.M u l t i-o b j e c t i v ee c o l o g i c a l o p e r a t i o nc a p a c i t ya n dp o t e n t i a la n a l y s i so f p l a i nr e s e r v o i r i nd e s e r t a r e aJ.J o u r n a l o fX ia nU n i v e r s i t yo fT e c h n o l o g y,2 0 2 3,3 9(1):9 6-1 0 4.收稿日期:2 0 2 2-0 4-1 4;网络首发日期:2 0 2 2-1 1-0 2网络首发地址:h t t p s:/k n s.c n k i.n e t/k c m s/d e t a i l/6 1.1 2 9 4.N.2 0 2 2 1 1 0 2.0 9 2 3.0 0 4.h t m l基金项目:国 家 自 然 科 学 基 金 资 助 项 目(5 2 1 7 9 0 2 5,5 1 8 7 9 2 1 3);陕 西 省 自 然 科 学 基 础 研 究 计 划 项 目(2 0 1 9 J LM-5 2,2 0 2 1 J LM-4 4);中国博士后科学基金资助项目(2 0 1 9 T 1 2 0 9 3 3,2 0 1 7 M 6 2 3 3 3 2 X B)第一作者:牟琳辉,男,硕士生,研究方向为水文学及水资源。E-m a i l:2 2 8 3 7 9 2 4 1 8q q.c o m通信作者:白涛,男,博士,副教授,硕导,研究方向为水资源系统工程。E-m a i l:b a i t a o x a u t.e d u.c n荒漠区平原水库多目标生态调度能力与潜力分析牟琳辉1,白 涛1,加沙尔蒙古拜2,梁宇航1(1.西安理工大学 西北旱区生态水利国家重点实验室,陕西 西安7 1 0 0 4 8;2.新疆维吾尔自治区伊犁哈萨克自治州伊犁水文勘测局,新疆 伊犁8 3 5 0 0 0)摘要:针对干旱荒漠区水资源匮乏、区域水资源分配不均等问题,提出一种优化调度方案。本文以塔里木河为例,建立水库多目标生态调度模型,采用融入人工经验的自迭代模拟优化算法求解。通过水库多目标生态调度能力与优化潜力的对比分析,最大程度地满足了枯水期河道内、河道外生态流量以及补湖水量,挖掘了荒漠区平原水库显著的生态调度潜力。结果表明:优化后大西海子水库多年下泄水量变化较实测值更加均匀;实测的大西海子下游3年断流,优化后仅有1年断流,枯水期断流天数大幅减少;在胡杨漂种的关键生态期,河道外可供生态水量较实测值提高了6 6%;补湖水量较实测值减少了2 8 9.5万m3,有效控制了台特玛湖补湖水量及其变化过程。研究结果优化了季节性河流的河道内、河道外以及入湖等多要素生态供水过程,对于干旱荒漠区内陆河水资源的集约利用、河岸植被的生态保护与修复具有重要的理论意义和应用价值。关键词:河道内、河道外生态流量;补湖水量;融入人工经验;生态调度潜力;耗散型河段中图分类号:T V 2 1 3.9 文献标志码:A 文章编号:1 0 0 6-4 7 1 0(2 0 2 3)0 1-0 0 9 6-0 9M u l t i-o b j e c t i v e e c o l o g i c a l o p e r a t i o nc a p a c i t ya n dp o t e n t i a l a n a l y s i so fp l a i nr e s e r v o i r i nd e s e r t a r e aMOUL i n h u i1,B A IT a o1,ME NG GU B A I J i a s h a e r2,L I ANGY u h a n g1(1.S t a t eK e yL a b o r a t o r yo fE c o-h y d r a u l i c s i nN o r t h w e s tA r i dR e g i o no fC h i n a,X ia nU n i v e r s i t yo fT e c h n o l o g y,X ia n7 1 0 0 4 8,C h i n a;2.K a z a kA u t o n o m o u sP r e f e c t u r eo fX i n j i a n gU y g u rA u t o n o m o u sI l iR e g i o nB u r e a uo fH y d r o l o g yI n v e s t i g a t i o n,I l i 8 3 5 0 0 0,C h i n a)A b s t r a c t:A i m i n ga t t h e s h o r t a g eo fw a t e r r e s o u r c e s a n d t h eu n e v e nd i s t r i b u t i o no f r e g i o n a lw a t e rr e s o u r c e s i na r i dd e s e r ta r e a s,a no p t i m a ls c h e d u l i n gs c h e m ei sp r o p o s e d.I nt h i sp a p e r,t a k i n gT a r i mR i v e ra sa ne x a m p l e,ar e s e r v o i rm u l t i-o b j e c t i v ee c o l o g i c a lo p e r a t i o nm o d e l i se s t a b l i s h e da n ds o l v e db yas e l f i t e r a t i v es i m u l a t i o no p t i m i z a t i o na l g o r i t h mi n c o r p o r a t i n gh u m a ne x p e r i e n c e.T h r o u g ht h ec o m p a r a t i v ea n a l y s i so fm u l t i-o b j e c t i v ee c o l o g i c a l r e g u l a t i o nc a p a c i t ya n do p t i m i z a-t i o np o t e n t i a lo f t h er e s e r v o i r,t h ee c o l o g i c a l f l o wi na n do u to f t h er i v e rc h a n n e l a n dt h ew a t e rv o l u m e f o r l a k er e p l e n i s h m e n t i nt h ed r ys e a s o na r em e t t ot h eg r e a t e s te x t e n t,w i t ht h es i g n i f i-c a n t e c o l o g i c a lr e g u l a t i o np o t e n t i a lo fp l a i nr e s e r v o i r si nt h ed e s e r ta r e ae x p l o i t e d.T h er e s u l t ss h o wt h a t:a f t e ro p t i m i z a t i o n,t h ev a r i a t i o no fd r a i n a g eq u a n t i t y i nD a x i h a i z iR e s e r v o i r i sm o r eu n i f o r mt h a nt h em e a s u r e dv a l u e;T h em e a s u r e df l o ww a sc u to f f f o r3y e a r s i nt h e l o w e r r e a-69 西安理工大学学报J o u r n a l o fX ia nU n i v e r s i t yo fT e c h n o l o g y(2 0 2 3)V o l.3 9N o.1 c h e so fD a x i h a i z i,b u t o n l y f o r 1y e a r a f t e r o p t i m i z a t i o n;I n t h e c r i t i c a l e c o l o g i c a l p e r i o do f p o p-u l u se u p h r a t i c ad r i f t i n gs p e c i e s,t h ee c o l o g i c a lw a t e ra v a i l a b l eo u t s i d et h ec h a n n e l i n c r e a s e db y6 6%c o m p a r e dw i t ht h em e a s u r e dv a l u e;C o m p a r e dw i t ht h em e a s u r e dv a l u e,t h ea m o u n to fl a k ew a t e rs u p p l e m e n t i n gd e c r e a s e db y2.8 91 06m3,w h i c he f f e c t i v e l yc o n t r o l l e dt h ea m o u n to fl a k ew a t e rs u p p l e m e n t i n ga n d i t sc h a n g i n gp r o c e s so fT e t m aL a k e.T h e r e s e a r c hr e s u l t so p t i m i z et h em u l t i-f a c t o r e c o l o g i c a lw a t e r s u p p l yp r o c e s s e so f s e a s o n a l r i v e r s,s u c ha s i na n do u t o f t h e r i v-e ra n d i n t o t h e l a k e,w h i c h i s o f t h e o r e t i c a l s i g n i f i c a n c e a n da p p l i c a t i o nv a l u e f o r t h e i n t e n s i v eu t i-l i z a t i o no fw a t e r r e s o u r c e s i n i n l a n dr i v e r s a n d t h e e c o l o g i c a l p r o t e c t i o na n dr e s t o r a t i o no f r i p a r i a nv e g e t a t i o n i na r i dd e s e r t a r e a s.K e yw o r d s:e c o l o g i c a l f l o wi na n do u to f r i v e rc o u r s e;l a k ew a t e rr e c h a r g e;i n c o r p o r a t ea r t i f i c i a le x p e r i e n c e;e c o l o g i c a l o p e r a t i o np o t e n t i a l;d i s s i p a t i v er e a c h 水资源是人类社会经济和生态向前发展不可或缺的基础资源1。随着如何实现生态系统的保护与修复成为热点问题,许多学者针对塔里木河流域(简称“塔河”)生态保护与修复开展了深入研究。邓铭江等2结合塔里木河下游生态输水的特点及方式、生态输水后的生态环境变化,对塔河下游生态调度策略与模式进行了初步探讨。王光焰等3指出塔河流域内用水存在大面积超载,导致中下游水量严重减少,引起生态问题,并提出了生态调度、水权配置和生态补偿、跨流域调水等未来生态修复方向。李丽琴等4针对内陆干旱区城市发展和生态环境保护的强烈互斥性,在整体识别内陆干旱区水循环与生态演变耦合作用机理的基础上,以塔河为例,构建了基于生态水文阈值调控的水资源多维均衡配置模型。李均力等5利用L a n d s a t8O L I、S e n t i n e l-2 A等遥感数据,分析了2 0 1 32 0 1 8年荒漠河岸带植被的时空变化特征,并结合地下水位数据分析了荒漠河岸带植被对生态输水的响应。Z h a n g等6利用生态系统服务价值(E S V)评价模型和灰色关联度,揭示了新疆生态系统服务价值的时空变化特征及其与人类活动的相互作用。白涛等7基于荒漠区耐旱性植被胡杨特殊的生物节律,提出了汊渗轮灌的灌溉理念,搭建了汊渗轮灌的基本框架。X i n等8计算了近3 0年来台特玛湖水域周边的植被覆盖率,分析了生态输水前后植被覆盖率与湖区的演变和驱动因素。综合以往塔河流域生态保护的研究成果,主要侧重于是否满足“大西海子水库年均下泄水量达到3.5亿m3、有水到达台特玛湖(简称台湖)”的目标,而对于如何高效利用生态水资源研究不足,缺乏在水量刚性约束下对下游生态水资源的集约利用,无法充分发挥大西海子水库对下游河道内生态流量、河道外生态补水、补湖水量等多项用水的联合调控作用,导致事先制定的水库调节体系难以应对水库下游生态环境新形势,生态系统的脆弱性无法得到根本性的改变。因此,本文以塔河流域为研究对象,结合干旱荒漠区自然属性,改变传统以单因子调度为主的思维定式,构建面向生态系统的多保护对象的水库中长期生态优化调度模型,采用融入人工经验的自迭代模拟优化算法求解,通过水库实测调度能力与优化潜力的对比分析,最大限度地满足枯水期河道内、河道外的生态流量以及补湖水量,创建一种新型集约调度模式,充分挖掘塔河下游显著的生态调度潜力。研究结果对干旱荒漠区耗散型河段的水资源集约利用、河岸植被的生态保护与修复具有一定的指导意义。1 研究区域概况西北干旱区具有降水稀少、蒸发强烈、有水就有绿洲、无水皆为荒漠的典型干旱特征,资源性缺水与生态环境的脆弱性使得流域生态环境保护面临严峻挑战9,1 0。塔河是我国最长的内陆河流,全长24 8 6k m,流域面积为2 5.8 6万k m21 0,1 1。自1 9 7 2年大西海子水库投运以来,下游3 5 0k m河道断流近3 0年,塔里木河也成为中国乃至世界生态最为脆弱的流域之一7,1 0,1 2。研究区位于东经3 9.4 5 4 0.5 5,北纬8 7.5 4 8 8.4 2 之间,如图1所示。为遏制塔河下游生态环境恶化,保障阻隔两大沙漠合拢的塔河下游生态防护体系,自2 0 0 02 0 2 1年,向大西海子水库下游河道实施了2 2次生态输水1 3,1 4。生态输水后,随着地下水位的抬升,天然植被得到拯救和恢复,植被的响应范围也逐渐扩大,但距实现保护塔河下游绿色走廊的目标尚有很大差距7。河道枯水期的季节性断流和粗犷的调度方式导致下游河道内生态流量的年际变化过大,生态系统碎片化严重;生态供水过程无法满足植被需水要求,导致在关键生态期无法为胡杨漂种、着床提供充足的生态水;缺乏水库、生态闸等控制性水利设施,导致大量生态水注入尾闾湖台湖,蒸发耗散严重,区域水资源分配严重不均7,1 5。塔河下游生态保护与修复任重道远。79 牟琳辉,等:荒漠区平原水库多目标生态调度能力与潜力分析 图1 塔里木河下游研究区示意图F i g.1 S c h e m a t i cd i a g r a mo f t h e l o w e rT a r i m R i v e r s t u d ya r e a 鉴于此,为优化区域水资源配置,实现旱区水资源集约利用,亟待创建一种新型水库优化调度模式。通过发挥大西海子水库对下游河道内生态流量、河道外生态补水、补湖水量等多项用水的联合调控作用,将水库生态调度与河道外保护性植被的生长节律相结合,构建水库新型集约调度模式,优化水库调度潜力,将生态水资源最大限度用于河道外植被生长、发育,减少补湖水量,提高河道内生态流量的保证程度,以最大效率合理利用塔河下游宝贵的生态水资源。2 生态调度模型的建立因塔河下游生态环境脆弱,多次放水实践证实,在没有水库调节的情况下,天然来水难以顺利流达台特玛湖,大西海子水库成为了一座专门用于生态供水的水库(无农业灌溉供水需求)。由于大西海子水库入库的资料较少,本文基于恰拉水文站不同保证率的年径流量成果,借用阿拉尔水文站到恰拉水文站区间单位长度河道衰减率,将恰拉水文站不同保证率的设计 年径流量 推算至 大 西 海 子 水 库 断面1 6。结果表明,受冰川融雪洪水影响,塔河下游的大西海子水库来水从每年4月开始,5月中下旬至6月下旬出现季节性断流,到7月下旬来水陡增,至9月中、下旬达到峰值,1 1月中下旬以后来水急剧减少,随后进入季节性断流时段。由此可见,大西海子水库来水主要集中在每年71 0月。同时,1 0月份为胡杨和柽柳的集中落种繁育期1 7。从生态保护与修复的角度看,生态调度输水时段应与天然植物落种时间相契合1 8,因此,本文以1 0月作为河道外胡杨林生态保护的关键期。根据塔河一期规划,大西海子水库年平均下泄水量不得低于3.5亿m3,但由于塔河干流缺水严重,在连续枯水年无法保证3.5亿m3的生态补水,需要到丰水年再对其补充,以保证年均下泄水量达3.5亿m3。本文将通过对长系列水文资料的计算,制定精细化的水库调度方案,最大限度保证每年的下泄水量更均衡,水资源分配更加合理。2.1 目标函数1)由于塔河为季节性河流,没有生态基流,本文将有水到达台特玛湖的下泄流量作为河道内生态流量,以提高河道内生态流量保证程度为目标,最大限度减少河道断流时间,避免周围生态环境遭到不可恢复的破坏。1N TNi=1Tj=1(q i,jt)-Ws(1)式中:N为调度期年数;T为一年内总月份;q i,j为调度期内第i年第j月河道内生态流量,m3/s;Ws为允许的河道内最小生态需水量,亿m3;t为月时段差;为计算精度,本文取0.0 1亿m3。1NNi=1PiP(2)式中:Pi为第i年河道内生态流量保证程度,%;P为允许的河道内最小生态流量保证程度,%。2)在合理的湖面面积条件下,保证合适的入湖水量,避免因湖面面积过大而无效蒸发。01N TNi=1Tj=1(q i,jt)W1(3)式中:q i,j为第i年第j月入湖流量,m3/s;W1为允许的最小多年平均入湖水量,亿m3。3)将节约的水资源最大限度地应用于河道外生态屏障的构筑,增加胡杨灌溉面积。m a xW=Ni=1Wi(1in)(4)式中:W为河道外可供生态水量,亿m3;Wi为第i年调度期内河道外可供生态水量,亿m3。89 西安理工大学学报(2 0 2 3)第3 9卷第1期 2.2 约束条件1)水库水量平衡约束Vt=Vt-1+It+It-12t-qt+qt-12t(5)式中:Vt为t时刻水库库容,亿m3;It为t时刻入库流量,m3/s;qt为t时刻下泄流量,m3/s;t为t时刻、t-1时刻的时段长。2)库容约束Vm i nVtVm a x(6)式中:Vm i n为死库容,亿m3;Vm a x为兴利库容,亿m3。3)水位约束Zm i nZtZm a x(7)式中:Zm i n为死水位,m;Zm a x为正常蓄水位,m。4)下泄流量约束Qm i ntQtQm a xt(8)式中:Qt为第t时刻流量;Qm i nt和Qm a xt为水库第t时刻的最小、最大限制流量,m3/s。5)下泄总水量约束1N TNi=1Tj=1(qi,jt)We(9)式中:qi,j为第i年第j月下泄生态流量,m3/s;We为允许的最小多年平均下泄生态水量,亿m3。6)变量非负约束等。3 模型求解3.1 模型处理一般来说,处理水库多目标生态调度主要有两种方法:一种是通过确定最重要的目标,将多个目标转化为一个主目标,然后将剩余的子目标作为约束条件1 9;另一种则是将各子目标合并为一个总体目标,一般通过子目标的权重系数2 0进行合并。求解子目标权重值的方法有很多,如层次分析法2 1-2 2、交互式多目标决策方法2 3-2 4等。本文选用第一种方法求解多目标模型。以河道外生态水量最大为主目标,将式(1)中的Ws转化为河道内生态水量约束,将式(2)中的P转化为生态流量保证程度约束,将式(3)中的Wl转换为补湖水量约束,从而将多目标问题转化为单目标问题,并在模型中融入人工经验,对关键问题进行求解。由于塔里木河存在季节性断流,水资源匮乏严重,本文将最大限度满足目标约束,以输出最优调度结果。3.2 融入人工经验的调度规则大西海子水库为年调节水库,本文以月为调度时段,年为调度周期,进行长系列调度模拟。塔河属于季节性河流,大西海子水库无法保证全年1 2个月放水,其调度规则如下:1)水库于当年4月起调,在蓄水过程中保证下泄水量不低于最小生态流量,若来水过小,无法满足下泄要求,则停止下泄;2)7月下旬起来水陡增,水库蓄水至正常蓄水位,若来水过小,无法保证水库于1 0月蓄满,则停止下泄;若水量较大,则保证下泄流量不低于河道内最小生态流量;3)生态关键期1 0月,水库以最大生态调度流量集中下泄,以实现1 0月河道外生态输水最大化以及输水入台湖的目标;4)若水库蓄水较多,生态关键期1 0月水库未放至死水位,则多余水量在1 1月至来年3月(枯水期)以最小生态流量均匀下泄至死水位,尽可能保证枯水期河道不断流,维系塔河下游生态系统的完整性。3.3 约束条件值的确定1)汛期为每年8月1日至9月3 0日,控制水位不超过8 4 9.2 5 m,水库在1 0月集中下泄至合理水位8 4 6.6 1 m(见表1)。表1 水库特征水位及库容T a b.1 C h a r a c t e r i s t i cw a t e r l e v e l a n ds t o r a g e c a p a c i t yo f r e s e r v o i r特征水位特征值/m特征库容特征值/(万m3)死水位8 4 6.6 1死库容6 6 0正常蓄水位8 4 9.2 5兴利库容80 3 1.1 2汛限水位8 4 8.6 0防洪库容58 7 0.3 0设计洪水位8 4 9.5 7拦洪库容96 9 6.7 6校核洪水位8 4 9.6 0调洪库容98 7 02)根据塔河一期综合治理规划,大西海子水库年均下泄水量不得小于3.5亿m3。3)根据下游河道实测断面,运用曼宁公式计算得到塔河下游最大生态调度流量为4 0m3/s,结合历年生态输水流量,保证塔河一期规划有水到达台特玛湖,据此确定塔河下游河道内最小生态流量为1 5m3/s2。4)年均河道内最小生态调度流量保证程度为8 5%。5)多目标生态需水量:河道内年最小生态需水量为2.0 41 08m32,两河区间年生 态需水量为1.3 91 08m32,维系台湖适宜湖面面积的生态需水量约为31 07m32。3.4 求解步骤本文采用自迭代模拟优化算法,在优化模型中融入人工经验进行求解,算法技术路线如图2,求解步骤如下:99 牟琳辉,等:荒漠区平原水库多目标生态调度能力与潜力分析 图2 模型求解步骤F i g.2 M o d e l s o l v i n gs t e p S t e p 1:输入基本资料,包括大西海子水库特征参数及来水过程等数据;S t e p 2:拟定初始下泄水量;S t e p 3:根据融入人工经验的调度规则,合理分配下泄水量;S t e p 4:根据长系列计算指标及目标函数,输出计算结果;S t e p 5:判断输出是否满足补湖水量约束,若不满足则反馈修正,重新拟定下泄水量;若满足,进入下一步;S t e p 6:判断输出是否满足年均河道内水量约束,若不满足则反馈修正,重新拟定下泄水量;若满足,进入下一步;S t e p 7:判断输出结果是否满足年均河道内生态流量保证程度约束,若不满足则反馈修正,重新拟定下泄水量;若满足,进入下一步;S t e p 8:输出最优调度结果,计算结束。4 实例应用4.1 调度范围及周期本次水库生态调度的范围为从大西海子水库至尾闾湖台湖。中长期调度以2 0 0 0年4月至2 0 2 0年3月共2 0年作为调度期,以月为计算时段。4.2 生态输水方式及水库调度网络节点1)生态输水方式。单通道输水,在来水量较少的情形下,利用其文阔尔河或老塔里木河其中的一条通道轮换输水;双通道输水,利用其文阔尔河和老塔里木河两条通道同时输水;汊河输水,在主河道输水过程中,利用两岸自然分布的汊河实施生态输水,扩大天然植被受水范围;面状输水,在河道上修建壅水坝或节制闸,营造河水漫溢条件,扩大河道沿线受水范围2 5。2)水库调度网络节点。大西海子水库至阿拉干段,主河道分为东西两支,西支为老塔里木河,东支为其文阔尔河,它们于阿拉干站汇合。地下水监测断面包括英苏、阿拉干、依干布及麻、台湖等。结合塔里木河下游耗散型河段水资源开发利用现状,绘制涵盖大西海子水库、水文站、3个重要控制断面、河道生态闸和湖泊等节点在内的水库调度网络节点图,如图3所示。通过优化调度,实现了控制补湖水量、保证河道内生态流量、提高两河区间生态供水的多生态目标。图3 塔里木河下游大西海子水库生态优化调度节点F i g.3 E c o l o g i c a l o p t i m i z a t i o ns c h e d u l i n gn o d eo fD a x i h a i z iR e s e r v o i r i nt h e l o w e r r e a c h e so fT a r i mR i v e r 001 西安理工大学学报(2 0 2 3)第3 9卷第1期 4.3 多目标结果分析为挖掘荒漠区平原水库的生态调度潜力,本文通过自主研发的漂流浮球勘测下游耗散型河段各断面的生态流量及其流达过程,以2 0年来大西海子水库生态输水后塔河下游各区域耗水过程作为实测能力,以转变新型调度模式后各区域水量变化过程作为优化潜力值,通过实测值、优化值及目标值三者的对比分析,揭示荒漠区平原水库调度能力和调度潜力之间的量化关系,主要包括水库调度过程分析、河道内外生态流量分析及补湖水量分析。4.3.1 水库调度过程分析在优化调度分析中,以大西海子水库长系列水位变化过程来验证模型及其算法的合理性。长系列水位变化过程见图4,下泄水量实测、优化调度过程见图5。1)经计算,2 0 0 0年4月至2 0 2 0年3月共2 4 1个月的中长期优化调度过程中,大西海子水库多年平均来水量5.2 5亿m3,多年出库下泄水量3.8 7亿m3,入库、出库水量差额与大西海子水库库容、水位变化一致,蓄放水过程满足水量平衡,验证了优化调度结果的合理性、可靠性。2)大西海子水库的水位在死水位8 4 6.6 1 m至正常蓄水位8 4 9.2 8 m之间波动,符合水位约束条件,多年平均水位为8 4 7.9 9 m;2 0年的中长期优化调度过程中,水库蓄满7次,库满率1 0.3%,放空3 3次,库空率4 6.6%,满足了水库年平均下泄水量3.5亿m3、有水到达台湖的基本目标,且大幅度提高了水库的运行效率。3)在水库下泄水量的刚性约束下,实测出库水量3年为零,出库水量变化呈现年内、年际分布极不均匀的特点,且出库水量主要集中在71 0月,导致枯水期河道断流,生态系统呈碎片化状态;从优化调度后出库水量的变化过程来看,水库4月起调,蓄水至1 0月集中下泄,枯水期均匀下泄,出库水量仅有1年为零,尽最大可能保障了河道内生态流量不为0,实现了年调节水库的“蓄丰补枯”作用。图4 长系列水位变化过程F i g.4 L o n gs e r i e so fw a t e r l e v e l c h a n g e s 图5 大西海子水库下泄水量优化与实测调度过程对比F i g.5 C o m p a r i s o nb e t w e e no p t i m i z a t i o no fd i s c h a r g ew a t e rv o l u m eo fD a x i h a i z iR e s e r v o i ra n da c t u a l o p e r a t i o np r o c e s s 101 牟琳辉,等:荒漠区平原水库多目标生态调度能力与潜力分析 4.3.2 河道内生态流量分析针对塔河下游河段耗散性严重、季节性断流且全年无稳定生态流量等实际问题,本文通过调节大西海子水库的调度模式,将塔河下游分为四个区段:大西海子英苏段、英苏阿拉干段、阿拉干依干布及麻段、依干布及麻台湖段。对比实测及优化后各区段年平均耗水量的差异性(见图6),分析塔河下游各区段耗水过程。图6 长系列实测、优化调度下各河段年平均耗水量对比F i g.6 C o m p a r i s o no f a n n u a l a v e r a g ew a t e r c o n s u m p t i o no f e a c hr i v e r s e c t i o nu n d e r l o n g-t e r mm e a s u r e da n do p t i m i z e do p e r a t i o n 由图6可知。1)实测调度满足生态流量的月份占比5 0.0%,生态流量不为0的年份占比8 5.0%,河道内生态流量及不稳定。优化后,满足生态流量的月份占比提高1 3.2%,年生态流量保证程度提高1 5.0%,减少断流天数超过15 0 0天,且均满足设计保证率要求,验证了本文优化调度模式的有效性,大幅提高了下游河道两岸生态系统的完整性,维护了河岸保护性植被的稳定性。2)相较实测中大量生态水资源在河道内蒸发耗散的现状,优化后河道内各区段的蒸发损耗水量占比减小,有效保障了河道生态流量不为0,为两河区间植被的生长发育奠定了坚实的水量基础。4.3.3 补湖水量分析研究表明,随着台湖面积的持续扩大,将严重危及区域生态环境健康7。因此,本节将通过优化调度,使台湖面积控制在合理区间。计算长系列优化调度前后的补湖水量过程,如图7所示。图7 长系列台湖补湖水量的实测、优化调度过程F i g.7 A c t u a lm e a s u r e m e n t a n do p t i m a l o p e r a t i o np r o c e s so f l a k er e c h a r g e w a t e rv o l u m eo f t h e l o n gs e r i e so fT a i h uL a k e 201 西安理工大学学报(2 0 2 3)第3 9卷第1期 1)实测补湖水量多年平均33 7 6万m3,优化后补湖水量多年平均30 8 6.4 5万m3,较实测值减少了2 8 9.5万m3,有效控制了台湖面积,减少了塔河下游生态水资源的无效蒸发。2)实测台湖补湖水量在1 4 8万m31 08 8 2万m3之间波动,且连续3年补湖水量为0,补湖水量变化极不均匀。优化后补湖水量趋于均匀化,补湖水量于2 0 0 22 0 0 3年出现最低值,为6 9 8万m3,2 0 1 02 0 1 1年出现最高值,为72 3 9万m3,仅一年入湖水量为零,在节省了大量生态水资源的同时,大大提升了补湖水量的稳定性,为稳定台湖周边生态环境提供了有力保证。4.3.4 河道外生态水量分析由4.3.3节可知,优化调度后补湖水量减少至一定范围有利于保障台湖周边生态稳定性。将可供水量最大限度应用于河道外生态修复,发挥塔河下游生态水资源的最大价值。综合考虑河岸保护性植被的生长节律,选定胡杨漂种的生态关键期1 0月令水库集中下泄,尽可能保障胡杨等保护性植被的飘种、着床、萌发、生长与发育。优化调度前后每年1 0月两河区间的供水过程如图8所示。图8 1 0月河道外可供水量对比图F i g.8 C o m p a r i s o nm a pf o ra v a i l a b l ew a t e rs u p p l yo u t s i d er i v e r i nO c t o b e r 1)在胡杨漂种的关键生态期,实测两河区间可供水量多年平均为92 7 7万m3,优化后两河区间可供水 量 平 均 增 加61 2 7万m3,较 实 测 值 提 高 了6 6%。大幅提高了河岸保护性植被的受水面积,抬高了两河区间地下水位,为建立“地下生态水银行”打下了坚实的水量基础。2)实测河道外供水量变化明显,且连续3年可供水量为0。优化后供水变化趋于稳定,在3年无入库水量的情况下,仅1年无法供水,为重建两河区间生态系统提供了重要保障。综上所述,通过采用新型调度模式,水库优化潜力较实测能力得到显著提高,最大限度满足了塔河下游多目标生态要求,有效缓解了塔河下游区域水资源供需矛盾,验证了大西海子水库显著的生态调度潜力。在水资源刚性约束下,有效提高了生态水资源的集约利用效率,提高了河道周边植被多样性和群落稳定性;减少补湖水量,重构区域水资源配置关系。5 结 论1)以河道外可供生态水量最大为目标函数,建立塔河下游多目标生态调度模型,采用自迭代模拟优化算法并融入人工经验进行求解,获得了水库中长期优化调度计算结果,验证了模型的合理性、可靠性。2)在满足塔河一期综合治理规划目标的前提下,结合耗散型河段的时空演变规律,综合考虑两河区间生态系统变化的响应关系,通过构建新型集约调度模式,优化了塔河下游水资源配置,保证在水资源刚性约束下实现塔河下游的生态效益最大化,验证了荒漠区平原水库显著的生态调度潜力。3)量化了水库调度运行中耦合径流过程与生态系统变化的响应关系,在采用新型调度模式后,优化调度结果满足了两河区间植被生态需水、补湖水量、河道内生态流量保证程度的要求,体现了干旱荒漠区水资源集约利用的治水方针。301 牟琳辉,等:荒漠区平原水库多目标生态调度能力与潜力分析 本文仅考虑了塔河下游大西海子水库调度下的生态水资源优化配置,下一步将考虑补湖水量与河道外生态水量的胁迫关系,以寻求适宜的台湖面积对应的入湖水量与河道外最大生态供水量的阈值关系,从而进一步提高水资源集约利用效率,助力塔河生态体系修复。参考文献:1刘哲杰,白涛,高凡,等.面向生态系统多对象保护与修复的水库优化调度J.干旱区研究,2 0 2 2,3 9(2):4 1 0-4 1 8.L I UZ h e j i e,B A IT a o,G AOF a n,e t a l.