２０１５年４月　水利水电快报　ＥＷＲＨＩ　第３６卷第４期　文章编号：１００６－００８１（２０１５）０４－００５４－０３　

贵州乌江梯级水库入库流量计算分析　及改进方法　

王俊莉　（贵州乌江水电开发有限责任公司，贵州贵阳　５５０００２）　

摘要：针对乌江梯级水库小时入库流量“锯齿”波动及负值、梯级区间流量负值等现象，对坝址入库和　区间流量计算的影响因素进行分析。以洪家渡、东风、乌江渡及构皮滩４座水库为例，提出了相应的改进计　算方法，并将七点滑动平均法应用于乌江梯级，明显提升了入库流量过程的平稳性，成效显著。从改进后的　实际应用情况看，只要使用方法得当，就可以提高入库流量计算精度，以便－为作好乌江梯级水库防洪决策　及调度服务。　关键词：梯级水库；入库流量；流量计算；乌江　中图法分类号：ＴＶ６９７．２１　文献标志码：Ａ　

水调自动化系统（以下简称“水调系统”）是一　个复杂且专业的信息自动化系统工程。自水调系统　建设运行以来，已成为贵州乌江水电开发有限责任　公司水电站远程集控中心（以下简称“贵州乌江公　司集控中心”）生产运行、调度管理的核心系统，为　乌江公司的相关业务开展提供了数据和技术支持。　随着该公司生产管理的不断提升，水调系统也越来　越受到重视。受发电、泄流、库水位测量和计算方法　等因素的影响，入库流量的计算结果偏差较大，过程　线也大多呈“锯齿”状剧烈波动，甚至出现负值，计　算结果失真，从而对区间实时流量过程计算造成一　定影响，无法确定汇流时间，对防洪决策、指挥调度　以及发电计划编制产生了较大影响…。因此，应对　入库流量加以分析。　

１　入库流量特性及作用　

１．１入库流量特性　入库流量是河川径流的一种特殊状态，其发生、　发展有特定规律，并表现出以下几点特性。　（１）时变性。受自然或人为因素的影响，人库流　量随时间而不断改变。　

（２）连续性。入库流量的变化连续不断，并遵循　质量守恒定律。　（３）非线性和非恒定性。与自然河道水流一样，　入库流量的变化不是均匀、线性的。　（４）不重复性。虽然河川径流有其本身的周期　性，但不会重复同样的变化轨迹。　（５）非模型化。不能通过建立计算模型来推求　入库流量。　（６）不可测量性。与河道断面流量的可测量性　不同，人库流量只能在发生后，利用公式进行反推，　才能得出“实际”值。　

１．２　入库流量在水库调度中的作用　入库流量是指导水库运行最直接有效的指标之　一。

水库优化调度、洪水调度以及洪水预报等都需要　

准确的入库流量作为支持，特别是在水库调度自动　化程度日益提高的今天，若入库流量计算精度偏低，　不但会限制水调自动化功能的实用性，也会影响到　电网端对电厂水调数据的使用。可以说，对入库流量　的准确推求，不仅能促进水库的安全经济运行，还可　促进电网的优化调度，有利于实现水资源的优化配　置，促进国民经济可持续发展。　

收稿日期：２０１５．０４—０９　作者简介：王俊莉，女，贵州乌江水电开发有限责任公司，高级工程师　

・５４・　王俊莉　贵州乌江梯级水库入库流量计算分析及改进方法　２　入库流量计算误差　２．１计算基本原理　目前，依据水量平衡原理和水库调度管理规范　来计算入库流量，水库的入库流量、出库流量及梯级　区间流量的通用计算方法如下。　出库流量＝各机组发电流量＋　各泄流设备弃水流量＋航运流量＋损失流量　（１）　入库流量＝出库流量＋　（时段末库容一时段初库容）／时段长　梯级区间流量＝下级电站入库流量一　

较小时，坝前水面线基本呈渐变状态，变化比较稳　定；当出库流量变化较大时，近坝区的水面线则有明　显的升跌变化，这会直接导致坝前水位的大幅波动，　进而影响人库流量的计算精度。　（４）出库流量计算。出库流量主要包括发电流　量、各泄流设施泄流、船闸流量等，各部分计算均可　能出现误差。　（５）静库容曲线。水库建成后，受水流冲刷、泥　沙淤积等影响，其河道地形必然发生变化，这将导致　库容曲线与建库前不一致，且建库后不同运行时期　的库容曲线也可能不一致。　

（２）　３　入库流量计算改进方法及分析　

上级电站出库流量　（３）　其中，式（１）中的航运流量仅存在于有独立航道的　水库，难以量化损失流量，在大多数系统中直接置　零。　

２．２计算误差来源　入库流量计算误差主要来源于以下几个方面。　（１）库水位。目前很多水库仅在坝前才有库水　位观测站，且以坝前水位作为整个水库的代表水位，　这对水库总库容的计算影响较大。坝前的机组启闭、　闸门开关等操作均会严重影响到当时的瞬时水位。　由于操作时间短，水库库容并未产生对应变化，从而　造成库容变化的假象。因此，库水位的代表性不强是　导致入库流量计算不准确的最主要因素。　（２）动库容影响。水库库容曲线是假定水库水　面为平面的情况下得到的，这只有当水流速度很小　时才可成立。这种静库容曲线在一般的湖泊型水库　或较短的河道型水库计算中基本可满足精度要求，　但对于较长的河道型水库而言，其水库末端总会存　在较大回水，从而使水面形成向上翘的曲线，而非水　平线。若不考虑动库容影响，必将影响人库流量的计　算精度，水库动静库容如图１所示。　图１　水库动静库容示意　（３）坝前闸门控泄。当相邻时段出库流量变化　目前，河道型水库水务计算的小时入库流量普　遍存在“锯齿”波动、入库流量负值、梯级区间流量　负值等现象属专业难题，故在水调自动化系统中，不　能采用简易方法直接处理。针对这一现状，三峡集团　等流域梯级公司均开展了专题研究，包括动库容分　析、泄流曲线校正、机组ＮＨＱ关系曲线修正、梯级水　量平衡关系率定、卡尔曼滤波校正、七点滑动平均校　正等，效果均十分显著。　本文仅针对贵州乌江公司所辖的乌江梯级水电　站（洪家渡、东风、乌江渡、构皮滩、索风营、思林、沙　沱、大花水、格里桥９级），其中前４座电站水库库容　大、调节性能良好，其余５座则为日调节电站。针对　任一水库，假设某ｔ时段内库水位变幅为△ｚ，对应的　水库蓄水量变幅为△　，则ｔ时段流量变幅为ａＱ＝　△Ｖｚ／ａｔ。由此可见，水库库容越大，相同水位变幅△ｚ　产生的ａＱ流量变幅就越大。　因此，对于乌江梯级水库人库流量计算改进方　法和分析仅以洪家渡、东风、乌江渡及构皮滩４座水　库（２０１３年１月１—５日）为例加以说明，与之相关　的改进算法可推广至全梯级应用。研究表明，平滑后　入库流量过程平稳性显著提升。　

３．１坝址入库流量　坝址入库流量计算改进方式通常有分段库容叠　加法、库容差修正法、卡尔曼滤波法和七点滑动平均　法等。　（１）分段库容叠加法。对于河道地形资料齐全　且库区遥测水位站点多的水库，可对河道进行分段。　首先按大断面资料积分，计算各分段库容曲线，然后　以各分段遥测水位查分段库容曲线得各分段库容，　叠加后即为总库容，最后再根据水量平衡方程计算　

・５５・　２０１５年４月　水利水电快报ＥＷＲＨＩ　第３６卷第４期　入库流量。在该法中，充分考虑了库区回水末端及沿　程附加比降的影响，是精度最高的入库流量算法，目　前仅应用于三峡水库，其余多数水库不具备基础资　料条件，难以推广。　（２）库容差修正法。假定计算入库偏差与库容　差之间存在一定相关性，建立计算入库偏差与库容　差的相关关系，进而计算修正入库流量。运用该法需　要入库流量偏差与库容差具有较好相关性，经分析　计算，乌江梯级各水库按该法计算的小时入库流量　过程改进效果欠佳，仍存在较大的“锯齿”波动和负　值现象。　（３）卡尔曼滤波校正法。将入库流量的　锯齿”　跳变视为随机干扰，运用离散状态空间方程，以“预　测一实测一修正”的顺序递推，从而实现入库流量　滤波和预测。在该法中，需要实测流量作为滤波依　据。若乌江梯级各调节型水库附近有代表性较好的　流量站或水文站，且具有小时流量实测数据，可尝试　采用该法进行人库流量修正。　（４）七点滑动平均法。设某入库流量序列为Ｑ　，　Ｑ：，…，Ｑ　．．，Ｑ　，其中Ｔ＞７　０则其七点滑动平均流量　序列计算公式为：　

ｇ￡＝　：１）　一　、’　一，　旦孚　２）　Ｑ１＋Ｑ２＋Ｑ３＋Ｑ４×２＋Ｑ５＋Ｑ６　７　，　（ｔ＝３）　

Ｑ　一３＋Ｑ　一２＋Ｑ　一１＋Ｑ　＋Ｑ　１＋Ｑ　＋２＋Ｑ　＋３　７　（３＜ｔ＜Ｔ一２）　Ｑ　一５＋Ｑ　一４＋ＱＴ一３×２＋ＱＴ一２＋Ｑｆ一１十Ｑｒ　７　（ｔ＝Ｔ一２）　Ｑ　一４＋Ｑ　一３×２＋Ｑｒ一２×２＋Ｑｒ一１＋Ｑｒ　７　（ｔ＝Ｔ一１）　

（　＝　）　（４）　式中，　为当前时段，ｈ；　为计算时段，ｈ；ｑ　为计算时　段入库流量，ｍ。／ｓ；Ｑ　为第一个时段入库流量，　ｍ　／ｓ；Ｑ　为第二个时段入库流量，ｍ。／ｓ；下同。Ｑ¨为　计算时段ｔ向前推一个时段对应的入库流量，ｍ　／ｓ；　Ｑ　为计算时段ｔ向前推二个时段对应的入库流量，　

・５６・　

ｍ’／ｓ；Ｑ川为当前时段　向前１　ｈ对应的入库流量，　ｍ。／ｓ；Ｑ　为当前时段　向前推２　ｈ对应的入库流　量，ｍ　／ｓ，下同。　对于ｔ＝Ｔ一２、Ｔ一１、Ｔ这３个时段，当前计算的　是暂时值。在未来时段滚动过程中，这３个时段的数　值会根据未来的变化趋势滚动修正。例如当时段ｔ　＝Ｔ＋１时，需重算ｔ＝Ｔ一２、Ｔ一１、Ｔ这３个时段，　其中Ｔ一２时段转为稳定值，以后无需继续修正，但　

一１、　时段还需滚动修正，后续时段以此类推。平　滑后若存在负值，则置为零。该法的物理意义明确，　具有良好的适应性，能较好地模拟小时入库流量过　程变化趋势，且无需任何参照依据和辅助资料，原理　简单、实现便捷，但在该法中，每次只能确定３　ｈ前　的入库流量过程，３　ｈ前至当前的３个入库流量值会　被滚动修正，实时性较差。　

３．２　区间入库流量　由于出库流量计算相对稳定，因此，区间入库流　量计算精度主要依赖于坝址入库流量和区间传播时　间。在坝址入库流量方面，可采用平滑入库流量作为　计算依据，确定区间传播时间则是计算区间入库流　量的难点。从直观上看，上级电站出库流量越大，区　间传播越快，反之则传播越慢。这一定性认识存在一　定局限性，且无法量化。事实上，根据水力学理论，河　道水流速度的计算公式应表达如下。　流速＝（水力半径的２／３次方）×　（坡度的平方根）／河床糙率　（５）　水力半径＝断面面积／除水面以外的断面周长　（６）　式（５）、（６）中，坡度为比值，水力半径单位为ｍ，流　速单位为ｍ／ｓ。　河床糙率与河床粗糙程度有关，对计算结果影　响较大，可参考表１。　表１　河床糙率表　

由此可见，影响区间传播的本质因素是水力半　径、坡度和河床糙率，而非流量大小。流量对区间传　（下转第７４页）