适应水文非一致条件的水库防洪优化调度

　　摘要：为完善变化环境下水库防洪调度方案，建立了集水文非一致性诊断、水文模型还现、设计洪水计算与水库调度规则优化于一体的水库适应性调控方法体系，并应用于沂沭泗水系仕阳水库.结果表明，仕阳水库入库洪水变异点发生在1975年;利用新安江模型对水文序列还现处理，得到满足一致性的设计洪水，其洪峰、最大24和最大72洪量均值，较原设计分别增大了19.7%、14.5%和10.0%;建立了以防洪效益最大为目标函数的水库优化调度模型，采用动态规划进行求解.提出的适应性调度方案在20年、100年、000年和5000年一遇设计洪水条件下，较水库现行设计方案最高洪水位分别降低了0.32、0.11、0.19和0.11m.

　　关键词：防洪减灾;水库优化调度;非一致性;动态规划;仕阳水库

　　水库是在降雨时空分布不均匀地区，缓解防洪与兴利矛盾，保障防洪、供水和生态安全的重要工程措施[12].截至201年底，我国已建成9882座水库，总库容707亿，仅大中型水库近690座，总库容超过243亿[3].可以说，在我国已几乎不存在完全意义上的天然流域了.近年来，气候变化和人类活动加剧，流域水文情势发生了较大变化，很多水库(湖泊)其入库洪水序列已经不再满足水文一致性假设，设计洪水、水库设计及运行调度的原有计算方法面临严峻挑战，据此给出的水库调度方案，或因高估设计洪水、导致水库功能发挥不充分，或因低估设计洪水、增大水库安全风险.因此，在诊断水文非一致特性的基础上，对水文序列进行一致性处理，并建立适应非一致性条件、充分发挥水库设计效益的调度方法，是当前亟待解决的重要研究课题.自2008年Milly等在Science上发表论文“Climatechange.Stationarityisdead:whitherwatermanagement?”以来，有关水文非一致性诊断方面的研究引起了国内外学者的广泛关注.先后提出了多种检验方法，主要有MannKendall、Spearman秩次相关检验法和贝叶斯方法等10.

　　由于单一检验方法存在不稳定性，2010年谢平等[1提出了融合多种方法的水文变异诊断系统，应用结果表明与单一诊断方法相比，综合诊断系统结果更加可靠，现已成为趋势性和跳跃性分析的常用方法.对已发生变异的水文序列进行还原或还现，再对其进行频率分析是非一致分布水文序列频率分析的重要内容[1.目前，还原或还现的方法主要包括：变异点前后序列与某一参数的关系分析法、时间序列的分解与合成方法和水文模型方法.

　　其中，水文模型法利用模型参数的变化反映流域下垫面变化带来的影响，通过不同时期的降雨资料与某一时期(变异点前或后)的参数结合，实现水文序列的还原与还现,1.由于该方法具有一定的物理意义和较强的解释能力，是目前较为常用的方法.王国庆等[1以黄河中游三川河流域为例，采用流域水文模型模拟了1970后的天然径流过程，定量评估了气候变化人类活动对流域径流的影响.王忠静等[1针对变化环境条件下水资源评价中的“还原失真”和“还原失效”问题，基于分布式水文模型提出了变化环境下水资源评价新方法.韩瑞光等[1应用河北雨洪模型对阜平站1980前后的洪水序列进行了研究，把洪水序列统一修正到80后下垫面条件下的洪水序列.

　　环境变化下如何更好地对水库进行适应性调度受到国内外学者的广泛关注[1.张玮等[1针对气候变化对水库调度带来的不利影响，提出一种基于DS理论的水库适应性调度规则.吴书悦等20评估了气候变化对径流和新安江水库调度的影响，并结合气候模型与新安江模型绘制了水库适应性调度图.为在一定约束条件下，实现水库防洪与兴利效益最大化，学者们提出了多个数学规划和启发式算法21.数学规划方法一般有严密的数学理论支撑，具有收敛速度快、约束条件明确的优点，虽在计算效率和灵活性上有时不如启发式算法，但在目标函数和水库数量较少的情况下优势明显.

　　其中，动态规划(ynamicprograming，DP)在利用递归公式逐时段寻优的同时，保证了全局最优性，是水库优化调度的常用方法，应用广泛3].综上所述，国内外学者分别在水文非一致性诊断、一致性处理与水库适应性调度等方面开展了系列研究，取得了丰富的研究成果，但鲜见适应非一致性条件的水库优化调度方法体系.为此，本文建立集水文非一致性诊断、水文模型还现、设计洪水计算与面向防洪安全的水库调度规则优化于一体的水库适应性调控方法体系，并在沭河流域仕阳水库开展应用研究，旨在定量解析人类活动对流域洪水影响程度的基础上，通过优化流域内现有水利工程的调度方式，降低环境变化对入库洪水改变带来的不利影响，提出现状条件下最大化防洪效益的水库调度方案，保障水库预期效益的发挥.

　　1适应水文非一致条件的水库优化调度方法

　　适应水文非一致条件的水库优化调度方法体系主要包括：(1)水文非一致性诊断.目的是判断水文序列是否存在非一致性变化，特别是变异特性，并确定变异点;(2)水文一致性处理.对于确定存在变异特性的水文序列，利用水文模型进行还现，得到满足一致性的水文序列;然后进行频率分析，计算满足一致性的水库设计洪水;(3)水库调度方案优化.建立最大化防洪效益的水库优化调度模型，通过模型求解得到变化条件下水库适应性调度方案.

　　(1)水文序列变异诊断：水文序列一致性审查是工程设计的前提[.首先采用线性滑动平均法、Hurst系数法及累积模比系数法对水文序列进行初步诊断，判断其序列是否可能存在变异;如发生变异，就要进一步诊断和分析突变性和趋势性.借鉴综合诊断法[1，本研究中采用的趋势性诊断方法包括非参数Speraman秩次相关检验法和线性趋势相关系数;变异诊断方法包括采用annendall(MK)检验法、有序聚类检验法、LeeHeghinian检验法、滑动检验法、滑动秩和检验法及MWP检验法.通过上述非参数统计方法获得数值分析结果后，根据需要再结合物理成因，确定最可能的变异点.(2)水文序列一致性处理及设计洪水推求：对于诊断出具有非一致性的水文序列，首先利用水文模型方法进行还现，获得满足一致性的水文序列，然后进行频率分析.主要步骤是：首先建立符合现状下垫面条件(变异点后)的水文模型;然后将不同时期降雨、蒸发等资料，输入具有良好精度的水文模型，计算得到现条件下满足一致性要求的径流过程，并根据需要遴选水文特征序列，例如洪峰、时段洪量;最后采用III概率分布对还现后的特征序列进行频率分析，进而基于同频率放大法对水库入库典型洪水过程进行放大，得到不同频率的设计洪水过程线.

　　(3)最大化防洪效益的水库优化调度模型：水文序列一致性条件发生改变，直接影响设计洪水成果的大小，从而影响水库预期功能的发挥.若变异后设计洪水变小，即相当于水库工程原设计时高估了设计洪水，则会导致水库功能发挥不充分;反之，若变异后设计洪水变大，导致水库规模一定条件下，水库运行安全风险加大.但无论如何变化，对于建成后的水库而言，最有效的办法是完善水库调度方案，实现水库预期功能最大化.

　　对于前一种情况，通过抬升退水期的运行水位，适量增加水库蓄量，一方面能减轻下游河道或者水库的防洪压力;另一方面由于水库防洪库容偏大，也能保障水库自身安全.对于第二种情况，则需要在洪水起涨前，适当降低水库的运行水位，为超过原设计的洪水腾空更多的防洪库容.为此，本研究构建了以防洪效益最大化(水库最高水位最小)为目标的水库优化调度模型.

　　2应用实例

　　2.1研究区概况与数据

　　仕阳水库位于沂沭泗流域中的沭河流域，控制流域面积281km，总库容1.2亿³，兴利库容0.6亿³，是一座具有多年调节能力的大型水利枢纽，以防洪、灌溉为主，兼顾发电与养鱼.水库2010年经除险加固后，满足00年一遇洪水设计、000年一遇洪水校核要求.仕阳水库保护下游1.63万hm耕地和30万人口，还有206国道、335省道、日东高速公路、胶新铁路等重要交通设施.下游小仕阳村处河道安全下泄流量为400m/s，招贤镇断面河道安全行洪流量910m/s.

　　仕阳水库汛限水位为153.00，设计洪水位为155.55m，校核洪水位为158.11m.水库现状防洪调度规则为：遭遇年一遇设计洪水，153.00起调，控泄400m/s，保护小仕阳村，计算最高雍水位155.00;00年一遇设计洪水，153.0起调，水位上涨到155.0前控泄400m/s，155.0后敞泄，计算最高雍水位155.33.此外，本研究还收集了1961—2013年间水库入库洪水和库区逐小时降雨等资料.并对入库洪水年最大洪峰流量序列、年最大24h洪量序列与年最大72h洪量序列进行初步分析，个序列存在增大趋势，尤其是年最大洪峰增加趋势较为明显.

　　2.2计算结果分析

　　2.2.1变异诊断

　　经等方法诊断仕阳水库1961—2013年实测入库洪水年最大洪峰流量序列、年最大24h洪量序列与年最大72h洪量序列.如表所示，各序列变异点大致有965、975、1976、1998、1999、2000和001年个变异点，其中1976年是1975年的延迟变异点，999—2001年是1998年的延迟变异点，故上述序列存在965、975和1998年个可能变异点.为了从上述个变异点中识别最可能变异点，从气候变化和人类活动角度开展了如下分析：

　　(1)仕阳水库实测洪水序列年限为1961—2013年，故从统计角度而言，靠近序列首端的1965年作为变异点不太可靠.(2)借助降雨量分析，由于遭遇“1997•”暴雨，故诊断出的998年延迟变异点不能认为是由于流域下垫面变化导致的变异点，而归结于气候所导致的变异点又过于片面，还需针对长系列水文气象资料进一步分析确认，多种诊断方法中998年出现的次数很少，暂不考虑这个变异点.

　　年最大洪峰序列分析中，取置信度水平为.05时，有序聚类、LeeHeghinian、滑动检验的值都在975年达到了最大，为.68;年最大4洪量检验中，有序聚类、LeeHeghinian、滑动检验的值都在1975年达到了最大，为3.45，滑动秩和检验和MWP分析检验值同样最大，分别为.8和.04，均为显著跳跃;年最大72h洪量检验结果与最大72h洪量检验结果类似.综合各诊断方法的诊断结论，仕阳水库年最大洪峰序列、年最大24h洪量序列、年最大72h洪量序列的变异点确定为1975年.

　　2.3讨论

　　为分析水文非一致性对水库防洪安全的影响，基于仕阳水库现行调度方式，分别对上述不考虑变异和考虑变异两种入库设计洪水进行了调度.由图可以看出，如不改变水库调度方式，水文变异导致的设计洪水增大，显著抬升水库最高壅水位，特别是对防洪标准内的洪水.其中年和00年一遇设计洪水下，水库最高水位分别增加了和.31m.此外，虽然水库最大下泄流量没有增加，但泄洪时长均显著增加.事实上，这些出现概率相对较大的洪水，如本例中的年一遇和00年一遇洪水更容易造成洪水灾害，是流域防御的重点对象[2.

　　因此，针对水文变异导致的防洪风险增加，水库有必要进行对应的改变.本文提出适应性调度方法的基本思想是，变异后的设计洪水经过优化调度后水库的最高洪水位与最大下泄流量均不超过现状调度下的特征值.然而，对于大多数水库，特别是单库，现行的削峰调度规则对于防洪来说已然是非常优秀的调度规则，其优化空间很小[3031.

　　为实现洪水适应性调度的目的，优化模型实质上是通过在洪峰到来前进行预泄，腾空了部分库容.可以预见的是，若不进行预泄，现状条件下由于洪水量级变化某些水库可能会出现调洪最高水位高于常规调度的最高水位，这更突显了进行适应性调度研究的重要性.此外，通过优化调度得到的解最优解是水库实际调度能够达到的上限，预报调度只有在预见期足够长、预报精度非常高的情况下才能接近该效果.

　　对于设计标准外的稀遇洪水，即使水文变异增大了设计洪水的洪峰与洪量，但由于原设计洪水已经很大，即使调度规则不改变，水库最高水位变化不明显，本例中000年一遇和000年一遇分别增大了.06和0.22m.000年一遇洪水下最高洪水位虽然变化不大，但洪水下泄提前，持续时间增大;而000年一遇设计洪水下，洪水下泄过程几乎相同，说明常规调度规则已没有进一步调节洪水的能力.本文提出适应性调度方法的结果表明，在超标准洪水下，仕阳水库仍能以控泄形式对洪水进行调节，实现保护水库自身和下游保护目的. 从仕阳水库的运用结果来看，本文提出的非一致性条件下水库适应性调度方法相较于水库现行调度技术能够更加充分地发挥水库预期功能.

　　在实际运用中，对于入库洪水序列不满足一致性假设的水库，根据本文提出的方法计算得到现状条件下设计洪水，然后与原设计洪水比较，当设计洪水确有改变时应重新制定水库调度规则.当设计洪水增大时，可在现状调度规则的基础上，降低水库的起调水位，来确保水库安全;当设计洪水变小时，可在现行调度规则的基础上，适当减小最大下泄流流量以减小下游防洪压力，或者增加水库蓄量，更好地发挥水库兴利效益.由于不同水库的工程调控能力和现行调度规则不同，在运用本文所提出的方法时，适应性调度规则的提出需要决策者通过优化或者多方案比较来最终确定.相关研究27指出未来气候变化可能导致淮河流域中小洪水(20年一遇)的洪量设计值较历史时期增大.

　　因此，结合气候模式与本文提出的适应性调控技术框架，预测未来降雨变化，研究未来情景下的水库防洪调度方案是保障仕阳及类似水库防洪安全的可行措施之一.结论受气候变化和人类活动的强烈影响，绝大部分水库入库洪水序列不再满足一致性假设，继续采用水库设计时形成的调度规则，往往难以充分发挥水库设计功能与效益，甚至可能带来防洪或供水风险.为此，本文提出了水文非一致性诊断→水文模型还现→设计洪水计算→面向防洪安全的水库优化调度于一体的水库适应性调控方法体系，并将其运用于仕阳水库.取得的主要结论有：

　　1)仕阳水库年最大洪峰序列、年最大24h洪量序列、年最大72h洪量序列的变异点均为1975年.考虑一致性影响后，仕阳水库设计洪水明显增大，特别是洪峰，5000年一遇、1000年一遇、100年一遇和20年一遇设计洪水洪峰分别增加了2065、1595、942和512m/s.2)如继续采用水库原调度方式，水文变异导致的设计洪水增大，将显著抬升水库最高壅水位，特别是对防洪标准内的洪水.其中，20年和100年一遇设计洪水下，仕阳水库最高水位分别增加了0.2和0.31m.3)针对满足一致性要求的设计洪水，通过优化调度，水库最高水位下降量在0.110.32m之间，20年一遇设计洪水情景水位下降最为明显;水库最大下泄流量均未超过原设计值，能满足保护下游小仕阳村的要求，且在超标准洪水情况下，全程未出现敞泄情况，较现状调度规则，更好地实现了对下游的防洪.

　　水文论文范例：考虑水文非一致性变异的河道基本生态流量

　　本研究中仕阳水库属于水文变异导致设计洪水增大、防洪风险增大的案例，适应性调度的主要目的是降低由水文变异增加的防洪风险，这对于防洪为主的流域具有很好的借鉴意义.可以预见的是，对于水文序列具有明显下降趋势的流域，如黄河流域32与海河流域33等北方缺水地区，适应性调度的主要目标应是在水文变异导致水资源不足的情况下，如何通过适应性调度增加洪水资源的利用量，有待进一步研究.

　　参考文献

　　[1]WangZZ,ZhangLL,ChengLetal.OptimizingoperatingrulesforareservoirsysteminNorthernChinaconsideringecologicalflowrequirementsandwaterusepriorities.JournalofWaterResourcesPlanningandManagement,2020,146(7):04020051.DOI:10.1061/(asce)wr.19435452.0001236.

　　[2]ang,ZhouCeds.ResearchondynamiccontrolofreservoirNormalelevationinfloodseasonandapplications.Beijing:NationalHydraulicandHydroelectricPublicationHouse,2006.王本德周惠成.水库汛限水位动态控制理论与方法及其运用.北京:中国水利水电出版社,2006.]

　　[3]MinistryofWaterResourcesofthePeople'sRepublicofChina.Regulationforcalculatingdesignfloodofwaterresourcesandhydropowerprojects:SL44—2006.Beijing:ChinaWater&PowerPress,2006.[中华人民共和国水利部.中国水利统计年鉴2018.北京:中国水利水电出版社,2018.]

　　[4]ZhangJY,WangGQ,JinJLetal.EvolutionandvariationcharacteristicsoftherecordedrunoffforthemajorriversinChinaduring1956—2018.AdvancesinWaterScience,2020,31(2):153161.张建云,王国庆,金君良等.1956—2018年中国江河径流演变及其变化特征.水科学进展,2020,31(2):153161.]

　　[5]WangJ,XuYP,WangYFetal.Impactsofanthropogenicactivityontheresponseofwaterleveltorainfallintheurbanizedplainrivernetwork:AcasestudyintheHangzhouJiaxingHuzhouregionofTaihuBasin.JLakeSci,2019,31(3):779787.DOI:10.18307/2019.0316.王杰,许有鹏,王跃峰等.平原河网地区人类活动对降雨水位关系的影响——以太湖流域杭嘉湖地区为例.湖泊科学,2019,31(3):779787.]

　　作者：王宗志1**，贺雨晖1,2，王坤1,3，王卫光2