实时洪水预报（共8篇）
实时洪水预报 篇1
0 引言
20世纪80年代前,我国水文学家一直致力于水文物理规律的研究,研制了适合我国自然地理和水文特性的一系列水文模型,其中比较典型的是河海大学赵人俊教授研发的新安江流域模型[1],在中国得到了广泛的应用,是联合国教科文组织推荐使用的十大水文模型之一。但是,中国幅员辽阔,自然地理和水文特性差异巨大,用单一模型很难解决不同地区的水文问题。在此期间,国内也有一些学者开始探讨水文过程的实时校正问题。进入20世纪90年代后,实时校正技术的研究得到了广泛的关注。目前,水文水资源预测预报实际上就是水文数学模型加上实时校正。
在实际洪水预报作业中,实时校正方法主要有以下3种:1)滤波方法[2,3],包括卡尔曼、自适应滤波等;2)随机过程方法,如时间序列模型;3)人工智能方法[4],包括人工神经网络、进化算法等。实时校正的变量可以是预报模型的参数,也可以是预报误差或状态变量。
汤河水库洪水预报对汤河水库及其下游的防洪安全极其重要,为此要对汤河水库实时洪水预报和校正模型进行研究,以提高洪水预报精度,确保汤河水库流域人民生命财产的安全。
1 汤河水库概况
汤河水库坝址位于辽宁省辽阳市弓长岭区境内的太子河一大支流汤河干流上,是一座以防洪、工业和生活用水为主,兼顾灌溉、发电、养鱼等综合利用的国家大Ⅱ型水利枢纽工程。水库为多年调节,按百年一遇洪水设计,可能最大洪水校核。坝址以上的控制流域面积为1 228 km2,流域多年平均降雨量771.2 mm,多年平均径流量2.89亿m3。水库总库容7.07亿m3,调洪库容3.68亿m3,兴利库容3.59亿m3,多年调节水量2.15亿m3,是鞍山、辽阳2座城市生活和工业用水的主要来源。水库枢纽工程由土坝、溢洪道、泄洪洞、水电站、引水建筑物等部分组成。
汤河水库所属流域属温带半湿润和半干旱大陆性季风气候区,四季分明,气候干燥,日照时间长,蒸发量大,土壤地表干燥,植被较差,包气带较厚。降雨量年际不均匀,年最大降水量为1964年的1 155.7 mm,年最小降水量为1959年的501.0 mm;降水量年内分布极不均匀,7,8月份降水量占全年降水量的49.4%。降水历时短,降水量空间分布不均匀。
目前,汤河水库以上流域设有8个遥测雨量站和2个遥测水位站。
2 汤河水库预报系统结构
根据汤河水库坝址以上水系情况,将大坝以上流域概化为若干子系统,概化图如图1所示,图中q为单元面积出流;Qi为支流出流;Q为流域总出流。
汤河水库以上流域位于辽宁省中部,降水量中等偏旱,土层较薄,所以汤河水库模型采用双层蒸发结构。汤河水库流域水文模型结构实际上是产流采用两水源新安江模型,汇流采用经验单位线。该模型的特点是,在产流计算中将相关经验与数学模型结合,即在考虑前期影响雨量时,采用前期影响雨量指标Pa,在蒸发计算中引进两层蒸散发模型,所以这一模型充分考虑了汤河水库以上流域自身的自然地理和水文气象特性。考虑到降雨的空间分布不均匀,采用分散参数性的洪水预报模型,即分单元流域计算产汇流。以雨量站划分单元面积,即每个雨量站为1个产汇流单元,然后,对每个单元流域的出流进行河道汇流计算,最后在出口处进行叠加得到水库的总入库流量过程,即形成单元面积——河网系统。每个单元流域的洪水预报可采用不同的降雨径流模型,包括数学和经验模型。在汤河水库以上流域采用汤河水库流域水文模型。
3 汤河水库洪水预报实时校正模型
汤河水库降雨主要集中于汛期,流域包气带较薄,产流集中在地面径流和壤中流部分,洪水过程表现为有明显的涨落过程,所以通常的降雨径流模型不可能精确预报所有的洪水过程,同时考虑到误差之间存在一定的相依性,故采用人工神经网络方法构建实时校正系统。
人工神经网络是一种由大量简单的人工神经元广泛连接而成的,用以模仿人脑神经网络的复杂网络系统。人工神经网络具有高维性、并行分布处理性、自适应、自组织、自学习的优良特性[5]。这种网络依靠系统的复杂程度,通过调整内部大量节点之间相互连接的关系,达到信息处理的目的。人工神经网络具有自学习和自适应的能力,可以通过预先提供的一批相互对应的输入和输出数据,分析掌握两者之间潜在的规律,最终根据这些规律,用新的输入数据推算输出结果,这种学习分析的过程被称为“训练”。目前,人工神经网络模型很多,考虑到生产实际中既要满足预报精度,又要考虑资料的可获得性和计算的简捷方便性,为此采用前馈型网络误差反向传播(BP)算法。
在实时校正模型中,一般采用多输入单输出模型,输出为预报误差,输入为与预报误差有关的参数,通常取前若干时段的预报误差。利用历史资料对网络进行训练,求得各个连接的权重值即可建立人工神经网络模型。
在汤河水库预报校正模型中,输入采用前3个时段入库流量预报值与实测值的误差,输出为误差校正值。实时校正模型如图2所示,图中ωij为上一单元与下一单元间的权重系数,也即实时校正模型需要通过实测资料进行参数的校准。
4 汤河水库洪水预报实时校正模型参数校准
4.1 降雨径流模型参数校准
选用汤河水库2000—2010年总共11 a的资料对降雨径流模型的参数进行校准,并在这些资料中选取17次中大小洪水进行模拟计算。模型校准的准则是采用实测与模拟总径流差值最小、确定性系数最大,实测与模拟洪峰流量差值最小,并考虑洪峰滞时,最后进行综合考虑与平衡。
通过对实测资料模拟运行,实测与模拟年径流平均相对误差为2.30%,最大相对误差为15.71%;平均确定性系数为0.92;实测与模拟洪峰流量平均相对误差为5.40%,合格率为91.20%;洪水滞时合格率为92.30%。
4.2 人工神经网络模型的训练与学习
人工神经网络模型的训练是通过观测资料求各个网络连接的权重值。人工神经网络模型的学习是利用已经训练得出的各个网络连接的权重,输入实测数据计算出输出值,考察人工神经网络模型的模拟精度。
在校准时,首先用降雨径流模型对17场洪水进行模拟计算,然后对模拟与实测的结果进行比较,计算出模拟误差系列。将误差系列分为2个样本,其中3/4的样本数据用于人工神经网络模型的训练,1/4的样本数据用于人工神经网络模型的学习。
5 实例分析
以汤河水库2010年8月5日和2011年8月8日的次洪为例,对模型进行检验。2010年8月5日洪水是由8月5日2∶00—17∶00和8月8日5∶00—9日11∶00的2场暴雨形成的,8月5日流域平均降雨量为50.8 mm,点最大降雨量为麻屯81.5 mm,最小降雨量为吉洞22.5 mm;8月8日流域平均降雨量为83.8 mm,点最大降雨量为麻屯106.5 mm,最小降雨量为汤河73.9 mm。2011年8月8日2∶00—9日17∶00流域平均降雨量为95.7 mm,点最大降雨量为吉洞135 mm,最小降雨量为麻屯76.0 mm。
通过对2场洪水的模拟计算,并用校正模型进行实时校正,得到最终的预报结果,计算成果分别如表1、图3和4所示,图3和4中横坐标的标值点时刻均为0时。
由表1可见,第1场洪水经实时校正后洪峰流量相对误差由7.80%降到1.60%,确定性系数则由0.85提高到0.95;第2场洪水经实时校正后洪峰流量相对误差由21.60%降到4.70%,确定性系数则由0.81提高到0.91。由此,汤河水库水文模型加实时校正模型后,可有效地提高实时洪水预报的精度。
6 结语
汤河水库自身和下游的防洪对洪水预报精度要求很高,为了进一步提高预报精度,在汤河水库洪水预报模型中引进了基于人工神经网络的实时校正模型。人工神经网络模型具有较好的适应性和灵活性,各神经元的连接权重可以在计算中实时校准。用实测洪水资料对汤河水库预报模型进行校准与检验,从计算结果可见,校正后的预报结果令人满意。目前生产上所用的实时校正模型很多,人工神经网络模型在实际预报中的稳定性等还有待于进一步的研究。
参考文献
[1]赵人俊.流域水文模型——新安江模型与陕北模型[M].北京:水利电力出版社,1984:32.
[2]宋文饶.卡尔曼滤波[M].北京:科学出版社,1991:123.
[3]丛日凡,姜志群,张淑芬,等.水库洪水预报实时校正算法[J].水土保持应用技术,2010(4):31.
[4]谢新民,杨小柳.半干旱半湿润地区枯水水资源实时预测理论与实践[M].北京:水利电力出版社,1999:47.
实时洪水预报 篇2
BP网络优化结合峰值识别理论进行洪水预报
随着国民经济的.不断发展,各经济部门对水文预报的要求越来越高.及时准确地预报洪水显得非常重要.进一步提高水文预报精度一直是一个努力的方向,本着这个原则,比较详细地介绍了将自适应遗产算法对BP网络的初始权重的优化,并结合峰值识别理论在水文预报中的应用,将遗传算法和人工神经网络的优点提取出来,取长补短,取得了较好的预报效果.
作 者：陈南祥 呼唤 CHEN Nan-xiang HU Huan 作者单位：华北水利水电学院,资源环境学院,河南,郑州,450008 刊 名：人民长江 PKU英文刊名：YANGTZE RIVER 年，卷(期)：2008 39(2) 分类号：P338 关键词：洪水预报 BP网络 遗传算法 自适应实时洪水预报 篇3
水文系统是一个受多种因素影响的复杂系统, 这些影响因素大多具有不确定性的时变特征。尤其是人类活动对流域下垫面的影响, 包括流域内大量小水电的开发及引调水工程的建设, 使得洪水预报精度难以令人满意。
自20世纪70年代以来, 随着计算机的普及应用, 水文模型得以广泛应用, 实时预报越来越受到关注。但由于人类活动对流域下垫面的影响, 实时预报精度很难达到生产要求。故在实际洪水预报作业中, 实时校正是一个不可或缺的部分。实时校正从方法论上可分为3类:1) 滤波方法, 包括卡尔曼滤波[1]、自适应滤波等;2) 时间序列方法[2];3) 人工智能方法[3], 包括人工神经网络、进化算法等。3种方法各有优劣, 具体如下:卡尔曼滤波适用于任何线性随机系统, 如果预报模型选用合适, 预见期和预报精度都不会损失, 但该法使用条件多, 且有大量的矩阵计算, 许多流域难以满足, 实际应用时需要灵活处理;时间序列法以误差序列的自相关作为校正方程, 方法简单, 计算量小, 应用较多的是AR模型;人工智能方法目前实际应用较少, 尚处于试验阶段, 不能确定校正的实际效果。由于时间序列模型结构简单, 要求资料相对较少, 在生产上得到广泛的应用[4,5], 四川宝珠寺水库就采用时间序列AR模型进行实时洪水预报校正。
1 实时预报模型
宝珠寺水库以上流域多年平均降水量为820 mm, 多年平均入库径流量为106亿m3, 属于半湿润地区, 从水文气象特征及入库流量过程性态分析, 该地区预报模型可选用新安江模型。
新安江模型为分散参数的概念性模型, 结构简单, 参数少且具有物理意义, 模型的参数必须根据历史资料进行率定。新安江模型根据流域雨量站网的分布, 用泰森多边形法将全流域分成为许多单元子流域, 对每个单元子流域作产汇流计算, 得出单元子流域的出口流量过程, 再进行出口以下的河道洪水演算, 求得流域出口的流量过程。把每个单元子流域的出口流量过程相加, 即为流域出口的总出流过程。分散性主要考虑影响产流诸因素在流域面上分布的不均匀性, 诸如降雨分布及下垫面不均匀性的影响, 防止均化对预报结果的影响。在每个单元流域上, 应用流域蓄水容量曲线考虑土壤含水量面上分布的不均匀性对产流的影响。新安江模型总径流分为地面径流RS、地下径流Rg、壤中流RI等3部分。模型的蒸散发部分采用3层蒸发模型, 即上层蒸发Eu, 下层蒸发EL和深层蒸发ED。模型的汇流计算统一采用汇流曲线法。
新安江模型已广泛用于中国许多地区的洪水预报作业中, 应用表明, 新安江模型适用于中国湿润半湿润流域, 具有良好的精度。
2 实时校正模型
宝珠寺水库上游建有碧口水库, 几年前在碧口水库下游又建了河床式电厂, 对宝珠寺水库入库流量产生了人为的调节作用, 干扰了宝珠寺水库正常的预报过程, 使宝珠寺的预报与实际过程产生较大的偏差, 需要对预报过程进行实时校正。由于宝珠寺水库入库缺乏水文站控制, 无法通过入库站采用其它方法进行预报校正, 同时受到资料条件的限制, 所以本研究采用时间序列方法进行实时校正。
实时校正的变量可以是预报模型的参数, 也可以是预报误差或状态变量, 但大多数是对预报误差系列进行校正。对预报参数进行校正存在较大的风险, 因为预报模型参数是通过长系列历史洪水资料进行率定得到的, 基于历史的平均情况, 同时水文模型参数在某种意义上具有一定的物理意义, 对模型参数进行校正可能会因某一场特殊分布的洪水得出不符合参数物理意义的参数值。同时, 由于洪水预报误差序列具有拖尾的特性, 且为白噪声, 故洪水预报误差序列为平稳时间序列, 可采用时间序列AR模型进行实时校正。
AR模型为预报误差自相关模型, 反映了预报误差的前后相依特性, 误差序列的拖尾性反映了这种相依性。作业预报中, 根据前几个时段的实测流量和预报流量的误差, 由校正模型 (AR模型) 计算出预报时刻的误差量, 将其加到新安江模型的预报值上, 即为该时刻的预报流量。
AR模型的阶数由AIC准则确定, 模型的系数由最小二乘法确定。
AR (p) 模型误差计算公式为xt=Φtxt-1+Φ2xt-2+…+Φpxt-p+εt
式中:p为模型阶数;xt为误差项;{εt}为白噪声序列。
式中:N为资料样本总数;为资料样本标准差。
使得AIC信息量取值最小的p值, 即是模型理想的阶。AIC信息量由2部分构成:前一部分体现模型的拟合好坏, 后一部分表明模型参数的多少。显然希望模型拟合得越精确越好, 但过高的精度要求又会导致参数的增多及模型的复杂, 反而会影响模型的拟合效果, 因此, 实质上, AIC信息量就是对拟合精度和参数个数两者加以适当权重, AIC的最小值处对应着最佳模型的阶数。
3 参数校准
3.1 降雨径流模型参数校准
新安江模型参数包括蒸散发折减系数K;张力水容量WM, 分为上层WUM、下层WLM和深层WDM;张力水蓄水容量曲线的指数B;不透水面积的比例IMP;深层蒸散发系数C;表层土自由水容量SM;表层自由水蓄水容量曲线指数EX;表层自由水蓄水库对地下水和壤中流的出流系数KI和Kg;地下水库和深层壤中流的消退系数Cg和GI;滞后演算法中的滞后时间与河网蓄水消退系数L和CS。
新安江模型参数的率定是分层次进行的。首先用日降雨、日流量、日蒸发资料对模型水量平衡参数进行校准, 日模型校准的准则采用实测年径流与模拟年径流差最小, 确定性系数DC最大。
利用洪水过程资料对模型影响洪水过程的参数进行校准, 次洪模型校准的准则采用实测洪峰流量与模拟洪峰流量差最小, 且确定性系数最大, 并综合考虑洪峰滞时。
最后, 综合日模型和次模型率定成果, 统筹考虑确定新安江模型使用的参数。
3.2 AR (p) 模型的定阶和参数的确定
用新安江模型计算出模拟流量过程后, 与实测 (反推) 流量过程相减, 得到预报误差过程系列。然后对预报误差系列进行白化处理 (零均值化) , 即误差系列减去均值, 得到零均值化后的残差系列。计算后的残差系列均值为0, 方差为有限值。取不同的p值, 计算残差系列的AIC值, 取AIC最小值所对应的p, 即为该模型的阶数p。根据最小二乘法, 计算AR (p) 模型的系数。
4 实例研究
本实例选用宝珠寺水库实时洪水预报模型利用AR模型进行校正。宝珠寺水电站位于四川省广元市境内白龙江下游宝轮镇附近, 以发电为主, 兼有灌溉、防洪等综合效益。总装机容量70万k W, 设计年发电量23.5亿k W时, 大坝坝顶全长524.5 m, 最大坝高132.0 m, 水库总库容25.5亿m3, 调节库容13.4亿m3。电站具有不完全年调节能力, 承担着四川省电网调峰、调频和事故备用重任。宝珠寺水库坝址以上流域面积28 428 km2, 坝址多年来平均流量336 m3/s, 多年平均年径流量106亿m3。
由于宝珠寺水库所处海拔高程较大, 年蒸发量较大, 需对模型结构进行适当的修改, 主要对流域3层蒸发结构进行调整, 考虑分汛期和非汛期设置深层蒸发系数C。经率定, 汛期深层蒸发系数为0.2, 非汛期为0.1。
选用宝珠寺水库2001—2010年10年资料, 包括各雨量站逐日降雨量、逐日入库流量、逐日蒸发量、次洪入库流量过程及其对应的各雨量站降雨过程, 对新安江和AR模型参数进行率定。根据模拟与实测流量误差系列, 确定AR模型的阶数为2。用AR模型校正后, 校正前后的流量误差系列如图1所示。实时校正前后实测与模拟年径流误差统计如表1所示。由表1可见, 经校正后, 精度有了较大的提高, 平均绝对误差由1.10降低为-0.41, 平均相对误差由-1.33%降低为-0.64%, 最大相对误差由14.05%降到-8.80%, 平均确定性系数由校正前的0.84提高到校正后的0.94。
用2011年8月26日和2012年7月3日次洪对实时预报和校正模型进行检验, 结果如图2和3所示。
由图2和3可见, 实时预报模型加了实时校正模型后, 对预报洪水过程精度有明显的改进, 尤其是洪水过程的拟合更令人满意。2011年8月26日次洪实测洪峰流量为1 331和1 610 m3/s, 预报为1 197和1 579 m3/s, 校正后为1 443和1 760 m3/s;校正前预报相对误差为-10.07%和1.93%, 校正后预报相对误差为8.39%和-9.31%;校正前确定性系数为0.86, 校正后为0.97。2012年7月3日次洪实测洪峰流量为1 847 m3/s, 预报为1 877 m3/s, 校正后为1 865 m3/s;校正前预报相对误差为-1.64%, 校正后为-0.97%;校正前确定性系数为0.87, 校正后为0.99。
5 结语
人类活动极大地改变了流域下垫面特性, 使得单一洪水预报模型预报精度很难达到生产要求, 必须对洪水预报进行实时校正。由于洪水预报误差具拖尾特性和接近于白噪声, 属于平稳时间序列, 可用时间序列模型进行实时校正, 时间序列模型作为实时校正模型在实际生产中得到广泛应用。宝珠寺水库受人类活动干扰很大, 常规预报模型精度不足, 加入AR模型进行校正后, 精度明显提高, 特别是对洪水过程的拟合更为理想。由于AR实时校正模型本身是依赖于预报误差相依特性的, 其对洪峰附近 (曲线拐点处) 的流量和滞时可能校正不够, 影响预报精度。与人工智能等实时校正模型相比, 时间序列实时校正模型概念更加清楚, 在资料充分的情况下模型制作更加简单, 且使用更加方便。
参考文献
[1]宋文饶.卡尔曼滤波[M].北京:科学出版社, 1991:123.
[2]王必义.应用AR (P) 模型预测水文系列[J].河南理工大学学报:自然科学版, 1993 (3) :58-66.
[3]丛日凡, 姜志群, 张淑芬, 等.水库洪水预报实时校正算法[J].水土保持应用技术, 2010 (4) :31-33.
[4]赵超, 洪华生, 张格平.实时校正模型的抗差递推算法[J].中国科学院研究生院学报, 2008 (5) :666-670.
清河水库洪水自动预报系统研究 篇4
洪水自动预报[1,2,3,4], 就是根据定时遥测的水文数据, 预报出未来一定时期内入库洪水总量、洪峰、峰现时间、入库洪水过程等。所谓定时, 就是取固定的时间间隔。一般系统取1 h间隔, 即系统每过1 h, 到每个准点时刻就会自动作出洪水预报。系统实现此种功能, 不需要任何的人为操作, 是计算机全自动进行模拟预报。
1 清河水库洪水自动预报系统主要研究内容
清河水库洪水自动预报主要研究内容包括:水文资料系列摘录与分析整理;模型结构选择;模型参数率定;洪水预报实时修正;模型参数修正;模型中间变量初值估计。
洪水自动预报结果自动存入数据库, 不作任何人工修改, 其预见期范围内的结果作为预报方案优劣的考核依据。但由于水文预报的精度依赖于实时遥测的降雨、水位、蒸发等的精度, 如果这些资料不精确或缺测, 必然会导致定时洪水预报结果的误差, 在此种情况下做出的洪水预报不作为考核范围。同时还可能由于水库库容曲线、下泄流量关系线等误差, 造成反推的实测流量误差大时, 也不作为考核范围。考虑到实时修正的重要性, 系统特设置了人性化的实时修正栏, 可人工在自动预报软件界面进行设置。
2 自动预报系统主要功能
(1) 报警流量设置。在自动预报软件里面, 特别添加了自动报警功能。报警以流量为判别条件, 不同的用户, 报警流量可以不同, 可以自行设置。自动报警功能可作为洪水期值班人员辅助工具, 以进行实时提醒。
(2) 实时显示预报过程曲线。自动预报界面默认情况为“实时显示清河入库预报流量过程”。主要实时显示最近1 h的入库流量预报值。其显示的预报值可为操作人员提供更加简洁明了的信息。该流量过程预报的图形显示可以通过“软件设置”菜单底下进行设置是否需要显示。
(3) 预报延时设置。自动洪水预报, 需要预报时段的降雨、水位、蒸发等资料, 而这些资料是由另外一台遥测工作站接收后送入服务器数据库的。自动洪水预报必须等这些资料入库后, 才能使用这些资料。遥测工作站准点开始把所有的水文数据接收到写入数据库需要时间, 因此要求自动预报有等待时间, 特别添加了预报延时设置。通过预报延时设置大大降低了实时信息滞后所带来的影响。但同时也需注意:自动预报延时, 既要考虑遥测接收和写库时间, 还要考虑遥测工作站与预报软件所在计算机时间的误差。假如2台计算机时间误差大, 还要校正计算机时间。
(4) 预报时间设置。因为自动预报每做1次, 都会自动记录。原来未做的, 当自动预报软件执行时会自动补做。但如果上一预报时间与当前时间差超过3个月, 自动预报就不能补做, 导致软件不能运行。此时需要人工设置“上一预报时间”为自动预报补做限定一个时间。
(5) 重复时段数设置。自动预报软件, 只采用当前时段的雨量、水位等资料信息, 实时资料信息, 一经使用, 就被保存到历史信息库中。假如使用时, 资料信息有误, 有误的信息也被保存到历史信息库中, 且会不断地影响着以后的预报结果。考虑到遥测系统, 雨量、水位资料等问题不易实时发现, 有时要过几小时之后被发现, 为免除这些资料误差对未来预报的影响, 软件设置了重复时段数参数。如重复时段数设置为5, 在5 h内改正的资料, 计算机会自动改正5 h内的历史信息库。这样, 5 h内改正了误差, 就不会影响以后的预报结果。
(6) 加密预报设置。考虑到有些山区性小流域或短预见期的洪水, 1 h预报1次还嫌不够, 可通过加密预报设置, 设定15 min的降雨量阈值, 则系统每过15 min, 如雨量超过设定的阈值, 就会自动加密预报。自动加密预报的设置为历时短、强度大的暴雨提供了更加精确的预报。
(7) 预报时段数设置。预报时段数与流域面积有关。流域面积越大, 预报时段数也要设置的越大。对于一般的水库流域, 预报时段数设置为24或36就可。
3 自动预报可能出现的问题及解决办法
洪水自动预报是整个洪水预报调度系统软件中技术性最强、难度最大、对环境要求最高的技术支持软件。其运行过程中, 可能发生的问题有水文资料问题、模型中间变量问题等。现介绍几个自动预报疑难问题分析与解决方法。
(1) 水文资料问题分析与解决方法。洪水自动预报需要有实时的雨量、水位、蒸发等资料作为模型输入, 这些资料的问题会引起水文预报的误差或甚至使模型无法预报。雨量资料值的偏大、偏小或缺测, 会使预报的流量偏大、偏小。全流域的雨量站全无雨量观测时或雨量冒大数至不合理范围时, 会使模型无法预报而缺报。水位资料值的偏大或偏小误差, 主要影响由水位反推的入库流量或河道流量, 从而影响水文预报的实时修正。蒸发资料误差主要影响流域土壤蓄水量和雨间蒸发等, 从而也影响水文预报的结果。
水文资料问题可以通过实时信息查询, 了解其具体问题。但洪水自动预报结果不能恢复。如果有些洪水还想用模型计算水文数据恢复后的结果, 需要通过人工干预洪水预报, 来修正洪水自动预报结果。
(2) 模型中间变量问题分析与解决方法。模型中间变量, 就是指模型中随时间变化的状态变量初始值。水文模型常有的中间变量为流域土壤蓄水量、坡面退水流量、河网退水流量等。由于模型中间变量是时空变化的, 需要由水文资料作为模型输入, 时间上连续的递推估计。如果水文资料有误差, 或模型计算不连续, 都会引起模型中间变量的误差, 从而影响洪水预报的精度。例如, 某个时段的雨量观测值偏大或偏小, 就会使模型的中间变量偏大或偏小, 其结果保存下来会影响下个时段的洪水预报偏大或偏小。当由于某些原因, 使预报间断一个时段或数个时段或更长时间未作预报, 这时模型只能调出停报前的中间变量, 显然是存在很大误差的。每个模型都有中间变量, 在系统软件启动时, 其中间变量的初值需要人工估计。在系统软件运行过程中, 正常情况下随时间变化的中间变量模型软件会自动逐时段递推估计出并保存到计算机中。若遇非正常情况, 如中心站计算机设备故障、遥测设备故障等原因使模型中间变量遭受破坏后, 系统软件重新启动时, 需重新估计模型中间变量初值。
参考文献
[1]李海明, 赵文忠.山口水库水情自动预报与洪水调度自动化研究[J].黑龙江水利科技, 2005, 33 (2) :37-38.
[2]赵咸榕, 陈宜安.碛口水库水情自动测报系统若干问题探讨[J].西北水电, 1997 (1) :12-16.
[3]肖毅.浅谈深圳水库水情自动测报系统软件设计特点[J].广东水电科技, 1994 (4) :45-47.
周家水库洪水预报方法研究 篇5
1 洪水预报方法
水库短期洪水预报, 是根据已经发生的雨情、水情以及水库的降雨径流与洪水规律, 预报水库即将出现的洪水情况, 为水库的防汛和调度运用提供科学依据[1,2,3]。
1.1 研究资料
周家水库的降雨、水位和出流等项观测资料比较短缺, 同时精度较差, 不能较准确地计算入库流量过程和分析该水库的降雨径流和洪水规律。在这种情况下, 借用流域降雨洪水特性相似的三道岭水库所在区域的降雨径流关系和洪水规律分析资料, 来预报该水库的洪水过程。
周家水库洪水预报方案:产流方案采用降雨径流相关法;洪峰流量推求采用净雨深与洪峰流量相关即峰量关系相关图。采用资料:周家水库位于大清河上游, 属于无实测水文资料地区, 由于地理和自然条件同三道岭水库相近, 故借用三道岭水库预报方案相关资料及成果, 进行该水库的洪水预报方案的编制。
1.2 产流计算
1.2.1 流域平均降雨量计算。
预报时流域平均降雨量采用周家水库站降雨量。由于周家水库没有降水资料记载, 预报分析借用三道岭水库预报方案相关资料。
1.2.2 产流参数的确定。
该水库无蒸发资料, 最大初损值 (Im) 及折算系数均采用大清河三道岭站, Im=130 mm, 其中折算系数作了必要的修正 (表1) 。
1.2.3 前期影响雨量Pa的计算。计算公式:
式 (1) 中:Pa, t、Pa, t+1分别为t、t+1日前期影响雨量 (mm) ;Pt为t日降雨量 (mm) ;Rt为Pt所产生的径流量 (mm) ;K为日折减系数。
前期影响雨量Pa自5月1日开始计算, 初始值为Pa=0 mm。Pa计算以Im为控制, 最大初损Im采用大清河资料为130 mm, 折算系数K见表1。晴天Pt<5 mm;阴天5 mm≤Pt<15 mm;雨天Pt≥15 mm。
1.2.4 径流深的计算。
根据周家水库的降雨量, 借用三道岭水库预报的K和Im值计算Pa, 用 (P+Pa) 直接从降雨径流相关图中查R值, 即为周家水库预报的洪水径流深。对1974—2005年间的23次降雨径流关系点据和15次为2006—2009年报汛资料计算降雨径流点据进行分析, 根据平均降雨量P、Pa、径流深R, 建立P+Pa与R的相关关系, 点在P+Pa与R相关图上。以P+Pa为纵坐标, R为横坐标, 建立P+Pa与R (降雨径流) 的相关图。利用所求的P+Pa之和, 在P+Pa与R相关图中查得R值。
1.3 汇流计算
1.3.1 峰量相关图。
对水库资料计算分析, 共收集到22次峰量相关数据, 以Qm为横坐标, R为纵坐标, 点绘Qm与R相关图。对短缺实测水文资料的中、小型水库, 同样可借用流域特性等相似的水文站 (或水库) 的峰量关系, 但必须根据集水面积, 应用经验公式将借用水文站 (库) 的洪峰流量 (Q借m) 转化为本水库的入库洪峰流量 (Q入m) 。其转化公式为:
式 (2) 中:Q入m为本水库的入库洪峰流量 (m3/s) ;Q借m为借用水文站 (库) 的洪峰流量 (m3/s) ;F为本水库的集水面积 (km2) ;F借为借用水文站 (库) 的集水面积 (km2) ;2/3为经验指数。
利用三道岭水库站峰量关系移植到周家水库折算系数K= (F周/F三) 2/3=0.464, 洪峰流量Q周=0.464*Q三, R不变, 点绘出R与Qm相关图。
1.3.2 推求洪水历时。
(1) 用三角形法推求洪水历时。用求三角形底边的方法推求洪水历时, 已知R (或W) 和Q入m, 可根据三角形的面积值即为洪水总量W, 则:
式 (3) 中:T为洪水历时 (h) ;W为洪水总量 (万m3) ;Q入m为入库洪峰流量 (m3/s) ;5.56为单位换算系数。
曾建立R与T关系, 因为资料不完整, 相关性不明显, 主要是观测时段太长。
(2) 用经验公式法推求洪水历时。用经验公式法推求洪水历时, 根据有效降雨历时和流域集流时间, 推求洪水历时的经验公式为:
式 (4) 中:T为洪水历时 (h) ;t有效为有效降雨历时 (h) , 即水库水位起涨到降雨停止的历时;τ为流域集水时间 (h) , 可用经验公式求得:
式 (5) 中, L为流域河道总长 (km) , V为平均流速 (m/s) 。一般中、小型水库地处丘陵地带, 粗略计算时可取V=1.4 m/s, 则τ=0.2 L。如有河道比降资料也可用满宁公式计算V值。
应用时, 已知有效降雨历时 (t有效) , 加上流域集流时间τ, 即可推求出本次洪水历时 (T) 。
1.4 方案精度评定
借用三道岭水库分析成果, 用三道岭水库典型洪水资料12年20场洪水和周家水库2年2场洪水分析, 确定了P+Pa与R降雨径流相关图, 径流深方案合格率为68.18%。净峰流量预报采用三道岭水库12年20次典型洪水和周家水库2年2场洪水分析, 点绘相关图, 预报合格率为54.55%。该方案分析成果精度较差, 使用价值较低。
2 水库抗洪能力分析
水库抗洪能力, 就是水库在某一库水位情况下, 水库能安全抗御多大暴雨洪水的实际能力。当水库水位很低、空库容很大时, 水库抗洪能力就大, 能拦蓄较多洪水, 即使产生较大的降雨洪水, 也可以安全渡讯;当水库水位已经很高;空库容已经不多, 这时水库的抗洪能力就小, 即使出现一场不大的降雨, 水库也可能抗御不了, 因而水库的危险性较大, 必须及早采取防汛措施[4,5]。
编制水库抗洪能力图表的基本原理是根据水库水量平衡方程式和降雨径流关系进行相关计算。在有了水库集水面积、水位—库容关系曲线、降雨径流相关图, 并确定了水库所允许的最高洪水位之后, 可以计算水库各级水位与不同前期影响雨量情况下的允许最大降雨量, 即水库的防洪能力。
2.1 水位不超过溢洪道底的防洪能力计算
计算公式:
式 (6) 中:P1为现时水位到溢洪道底高程的允许最大降雨量 (mm) ;W1为现时水位以上到溢洪道堰底尚可蓄水量 (万m3) ;F为水库集雨面积 (km2) ;α为径流系数。
2.2 水位超过溢洪道底的防洪能力计算
计算公式:
式 (7) 中:W防为现时水位到允许最高水位之间还可以防御的水量;W1为现时水位到溢洪道底之间还能拦蓄的水量 (万m3) ;W2为溢洪道底以上各级水位到允许最高洪水位之间调洪库容 (万m3) ;△V调为溢洪道底以上各级水位至允许最高水位库容 (万m3) ;W泄为在某一段时间内溢洪道泄出的水量 (万m3) , 即W泄=1/2Q泄m×0.36△t=0.18 Q泄m△t;△t为泄流历时 (h) 。分别对周家水库不泄洪和泄洪的情况下进行抗洪能力计算 (表2、表3) 。
3 结语
水库洪水预报调度的目的, 就是要求事先预报出在一次暴雨洪水过程中, 水库可能发生的最高水位、最大泄洪流量、开始溢洪时间和最高水位及最大泄洪流量的出现时间。以便进行科学的调度运用, 正确处理好蓄洪与泄洪的关系, 以确保水库工程和下游河道的安全, 并合理拦蓄洪水, 充分发挥工程效益。周家水库缺少长期的水文观测资料积累, 这给后续水库运行管理、水利水文分析计算以及洪水预报调度方案的编制带来极大的困难。
注:F=42 km2;Im=130 mm;泄流历时ΔT=24 h;Pa=0 mm。
参考文献
[1]张改红, 周惠成.利用防洪预报调度调整汛限水位的风险分析[J].水力发电, 2010, 36 (1) :26-28, 82.
[2]彭涛, 张利平, 沈铁元, 等.现代水文气象预报技术研究进展[J].人民黄河, 2010, 32 (9) :29-32.
[3]白继平, 丁书红, 许东红, 等.洪水预报技术研究进展与展望[J].河南水利与南水北调, 2010 (8) :14-15, 18.
[4]丛日凡, 姜志群, 张淑芬, 等.水库洪水预报实时校正算法[J].水土保持应用技术, 2010 (4) :31-32.
南四湖洪水预报调度模型研究 篇6
南四湖为南阳、独山、昭阳、微山4湖的总称。南四湖入湖主要河流共53条,其中29条入上级湖,24条入下级湖,出湖河流共4条。入上级湖的河流中有控制站的主要有:湖东的白马河、四河、光府河、城河、北沙河;湖西的洙赵新河、梁济运河、万福河、复新河、东鱼河。入下级湖的河流中有控制站的主要有:十字河、不牢河、新薛河、郑集河、沿河。出湖河流为老运河、韩庄运河、伊家河、不牢河,自东向西分布于湖区南部。南四湖流域多年平均降水量695.2 mm,多年平均径流量29.6亿m3。
南四湖流域水系复杂,湖东为山洪河道,源短流急;湖西为平原坡水河道,集流入湖缓慢。在南四湖洪水预报时,没有比较实用的预报调度系统。在实时性和准确性方面,已经不能满足防汛调度对洪水预报信息的需求,因此,开展南四湖洪水预报调度模型研究具有重要的理论意义和实际应用价值。为了解决南四湖洪水预报调度中存在的问题,构建了新安江模型,率定有实测资料的11个子流域,对于主要控制站以下的无实测资料区间,考虑移用下垫面情况相似的相邻流域参数,使预报结果达到一定精度[1]。
1 新安江模型应用
1.1 模型概述
新安江模型是一个概念性流域水文模型[2]。根据降水和下垫面的水文、地理情况将流域划分为若干个单元,对每个单元流域作产汇流计算,再将每个单元面积预报的流量过程演算到流域出口后叠加起来即为整个流域的预报流量过程。河道洪水演算采用马斯京根法。单元流域的新安江模型由4个层次组成:三层蒸散发计算、蓄满产流计算、三水源划分,以及坡地、河网汇流计算。
1.2 模型验证与应用分析
根据流域水文站网的报汛情况和流域洪水的汇流时间,南四湖水库洪水预报采用6 h时段长度。各子流域均采用新安江模型进行洪水模拟,参数律定。限于篇幅,仅列举鱼城、孙庄及腾县3个子流域结果。20世纪70年代该流域是天然流域,利用历史资料进行验证,模型效果较好,90年代后,由于水利工程的建设,破坏了流域的天然特性,预报效果较差,为此系统采用20世纪90年代及以前的资料进行验证。
1.2.1 模型验证
鱼城流域内有5个雨量站,选取1971—1995年的27场洪水资料;孙庄流域内有3个雨量站,选取1971—1995年的29场洪水资料;腾县流域内有3个雨量站,选取1962—1995年的32场洪水资料。从洪量的模拟精度来看,鱼城、孙庄和腾县3个流域合格率分别为88.9%,75.9%和75.0%,其中鱼城精度达到甲级标准,孙庄和腾县达到乙级标准。图1为鱼城流域1973-07-14 T 8:00—08-05 T 14:00实测与模拟洪水过程线,图2为孙庄流域1973-07-14 T 20:00—07-19 T 8:00实测与模拟洪水过程线(图1、2中各纵向虚线均代表所标日期的零时),表1为各流域新安江模型参数。
通过对比鱼城、孙庄及腾县3个流域的参数值可以发现:1)由于KG+KI=1,EX=0,此时的新安江模型实际是二水源新安江模型[3];2)孙庄流域的上层张力水容量WUM较其他2个流域要小,表明此流域的植被和土壤发育状况不及另外2个流域;3)鱼城流域的河网水流消退系数CS比较大,参数CS决定于河网的地貌,CS值愈大,流域调节能力愈大[4]。腾县流域调蓄能力明显不及鱼城和孙庄流域,这一点由自由水容量参数SM也可看出,究其原因,腾县流域丘陵地区占70%,属山区地貌,而鱼城和孙庄流域属湖西平原地貌。此外,鱼城流域面积远大于另外2个流域,流域的调蓄能力自然也较另外2个流域要大。
1.2.2 实际应用
在20世纪90年代后期,南四湖流域部分河道内进行了水利工程建设,修筑了拦河闸坝拦蓄上游来水,工程的建设改变了天然河道特性。通过实际应用来看,鱼城站实际应用效果较好,孙庄和藤县2站受工程影响较为明显。详细情况如图3～5所示(此3图中各纵向虚线均代表所标日期的零时)。
1.3 无资料区间处理
上级湖流域控制流域面积为27 263 km2,湖区面积为602 km2,水文站以上控制面积为18 564 km2,这样湖东无资料地区786 km2,湖西7 004 km2,将湖东和湖西无资料区间流域各分为4块,湖东每块为196.5 km2,湖西每块为1 751 km2,湖区单独作为1块。下级湖流域控制流域面积为4 250 km2,湖区面积为666 km2。将下级湖区间分为4块,湖东湖西各2块,每块820 km2,湖区作为1块。
这样上下级湖无资料流域和湖面分为14块,由于没有直接的流量资料,这些子流域水文模型的参数无法直接率定,对于产流和分水源参数,在考虑东西2岸地理差异的基础上,从邻近相似流域直接移用,对于汇流参数,根据地形和汇流河段的长度,估算得出。
2 洪水调度模型
南四湖是1个浅水型河流堰塞湖,湖内芦苇、水草生长茂盛,加之群众在湖内围湖生产养鱼等的影响,导致湖内洪水下泄缓慢。2级坝段闸前和闸后芦苇和其他人工障碍的阻水十分严重。历史上南四湖地区就是洪水灾害的多发区,本次仅考虑南四湖防洪调度。
2.1 基本原理
水库调洪计算的直接目的,在于求出水库逐时段的蓄水、泄水变化过程,从而获得调节该次洪水后的水库最高洪水位和最大下泄流量,以供进一步防洪计算分析之用。
水库调洪演算有2种情况,一种为下泄流量受控制的调洪演算,其控制调节方式由水库防洪运行规则决定。这只能在有闸门控制的泄洪设备条件下才存在。
另一种情况为自由泄流条件下的调洪演算,无闸门溢洪道泄流、或设闸门的泄洪设备闸门开启程度一定的条件下泄流均属此种情况。对此种情况,若按静库容条件考虑,需联解水库水量平衡方程和相应水库蓄泄方程才能实现逐时段的调洪演算。对调洪中过程中任一∆t(∆t=t2-t1)时段,计算式可表示如下:
式中:Q1,Q2分别为时段初和末入库流量;q1i,q2i分别为时段初和末第i号闸门或泄流设施的下泄流量;q1j,q2j分别为时段初和末第j号机组的下泄流量;V1,V2分别为时段初和末水库蓄水量。
蓄泄方程qi=ƒ(Z)表示闸门开度不变的条件下水库水位与泄流量之间的关系[5]。
对于控制泄流的计算可采用列表法,方法非常简单,这里不加以叙述。对于自由泄流的情况,需要联解式(1)和(2)及水库库容曲线(Z～V)。这里解析解是不可能的,需要使用试算求解。对于任意的∆t时段,在南四湖流域不需考虑机组下泄流量(因为南四湖没有安装发电机组),其Q1,Q2,q1i,V1为已知,通过试算法可求解q2i和V2。
2.2 约束条件
在进行模型计算时,主要考虑以下4个约束条件:
1)水库最高水位约束
式中:Zt为t时刻水库水位;Zmax(t)为t时刻容许最高水位,通常采用防洪高水位。
2)水库泄流能力约束
式中:qt为t时刻的下泄量;q(Zt)为t时刻对应于水位Zt的下泄能力。
3)出库变幅约束
式中:|qt–qt-1|为相邻时段出库流量变幅;∆qm为相邻时段出库流量变幅的容许值。
4)下游防洪能力约束
式中:qt为t时刻的下泄量;qmax为下游防洪安全泄量。
2.3 调度原则
根据《沂沭泗洪水调度方案》,涉及到南四湖的洪水调度方案主要概括为以下几点:
1)当上级湖南阳站水位达到34.2 m并继续上涨时,2级坝枢纽开闸泄洪,视水情上级湖洪水尽量下泄。
预报南阳站水位超过36.5 m,2级坝枢纽敞泄。当南阳站水位超过36.5 m时,滨湖洼地滞洪。
2)当下级湖微山站水位达到32.5 m并继续上涨时,韩庄枢纽开闸,视南四湖、中运河、骆马湖水情,下级湖洪水尽量下泄。当预报微山站水位不超过36.0 m,如中运河运河站水位达到26.50 m或骆马湖水位达到25.0 m时,韩庄枢纽控制下泄。
预报微山站水位超过36.0 m,韩庄枢纽尽量泄洪,尽可能控制中运河运河站流量不超过5 500 m3/s。
当微山站水位超过36.0 m时,韩庄枢纽敞泄,滨湖洼地滞洪;在不影响徐州城市、工矿安全的前提下,蔺家坝闸参加泄洪。
本节中的水位与流量均为目前沂沭泗防洪调度方案中的控制水位与流量。
2.4 调度模式
根据调度原则,确定2种防洪调度模式,一种是按照调度规则放水的情况下推求上级湖和下级湖的蓄洪和水位变化过程;另一种是各个闸都按照指定的放水过程推求上级湖和下级湖的蓄洪和水位变化过程,整个放水过程由人工控制,充分考虑了调度员的经验,使得洪水调度更具灵活性。从调度结果可以发现,经过调度之后,上级湖与下级湖的水位均可控制到汛限水位,保证了实际洪水调度的灵活性。
3 结语
南四湖流域水系复杂,本次构建新安江模型进行洪水预报,并利用洪水调度模型分别对上、下级湖入湖洪水进行调度,基本满足南四湖防汛调度对洪水预报信息的需求。通过初步的研究和探索,对南四湖流域的水文特性有了一定的认识,调度结果为南四湖防灾减灾提供依据,人工干预的调度方式使调度员可以根据实际情况调整下泄流量,保证洪水调度的灵活性。但由于资料缺乏,仍有不少问题留待解决:
1)随着流域内水利工程的大规模建设,流域的下垫面条件发生较大变化,本次仍采用1套参数,可以考虑将已有洪水资料分为80年以前和80年以后分别率定以适应流域变化,提高模拟精度;
2)调洪演算的条件由于资料的原因不能完全满足,中运河运河镇或者骆马湖的水位不能及时取得,并且防洪高水位在调度方案中没有给出。因此,对于该模型的实际运行情况,目前只能按照普通的调洪演算进行。
参考文献
[1]李致家,孔凡哲,王栋,等.现代水文模拟与预报技术[M].南京:河海大学出版社,2010:247-248.
[2]赵人俊.流域水文模拟—新安江模型和陕北模型[M].北京:水利电力出社,1984:11-17.
[3]李致家.水文模型的应用与研究[M].南京:河海大学出版社,2008:34-37.
[4]王义民,黄强,畅建霞,等.基于面向对象技术的水库洪水调度仿真模型研究[J].2003(4):8-14.
降雨径流型洪水预报方法实例分析 篇7
关键词：中小流域,下垫面,产汇流影响,相关图,预报模型
1降雨径流预报介绍
针对某一流域, 降雨须经过蒸发、下渗、坡地产流和河槽汇流几个阶段才能形成一次洪水过程, 见图1。
一次降雨形成径流可分为地面径流 (图2中1、2部分) 和地下径流 (图2中3、4部分) , 地下径流分浅层和深层地下径流, 深层地下径流由前期降雨形成。降雨径流预报法就是根据已知降雨情况来推求其产生的洪水过程。
2降雨径流预报影响因素
自然地理、气象和土壤湿润程度是影响降雨在地面形成径流的主要因素。气象因素包括降雨量级、强度、降雨在时间和空间上的分布等;流域自然地理特性包括流域面积、植被覆盖率、土壤结构等;流域降雨前期土壤湿润情况是影响降雨产流的重要因素, 因为其直接决定着下渗损失雨量的大小。上述影响因素成因复杂, 均难以定量掌握。
3示例流域介绍
永翠河流域位于伊春市带岭区, 流域面积约677km2, 河长67km。流域地貌属小兴安岭山地, 植被良好, 地下水丰富。流域为黑龙江省暴雨中心之一。流域位于黑龙江省多年冻土分区南部边缘, 年内土壤11月结冻, 第二年4月解冻。一些山地阴坡或塔头地带, 年内始终存在小兴安岭岛状冻土, 最大冻土深约2.5m, 冻土总体分布自上游至下游面积、深度逐渐减小。流域冻土导致降雨下渗异常, 透水较差, 冻土融冻时吸收热量, 抑制蒸发, 实际蒸发小于正常蒸发, 土壤含水量大于无冻土区。该区域因冻土存在使降雨下渗情况更加区域复杂化, 难以用单一方法来预报其降雨产流过程。
4预报方法介绍
流域下垫面影响因素复杂, 其产汇流规律复杂多变, 根据流域特性, 采用降雨~洪峰流量相关图和拟合预报方程法来预报流域洪峰流量。
4.1相关图绘制
预报影响因素之间的经验相关图间接反应了因素对洪水的影响情况, 在明确主要影响预报要素因子时该方法是较实用的预报方法。经分析, 选择场次降雨量和流域前期土壤含水量与洪峰流量建立相关关系 (见图3) , 相关点据48个。图3点据分四个带状, 共定出4条相关线。1、2号线间部分点据降雨强度均小于0.5mm/h, 降雨中心多数在流域上游区域, 洪峰流量偏小。2、3号线间降雨强度在0.5-1.0mm/h之间, 降雨中心在流域中、下游。3、4号线间降雨强度在1.0-8.0mm/h之间, 洪峰流量与同级降雨相比偏大。
4.2回归方程模型法预报
选择上述实测点据, 利用降雨强度 (pi) 、流域冻土深度 (hi) 、平均雨量、前期土壤含水量 (Pa) 与洪峰流量 (qm) 建立相关拟合分析方程。预报方程为:
式中:X为系数, C为常数, 利用实测点据率定出方程系数后预报方程为:。洪水拟合见表1, 合格率为62%, 达到丙级预报标准。
5结束语
(1) 受冻土影响流域的下垫面环境复杂多变, 在降雨蓄渗、产流方面直接影响着流域洪水规律特征。由于冻土分布面积和埋藏深度的不确定性, 使得传统单一理论的洪水预报方法难以在该流域得以应用, 因此, 必须考虑采用多种方法结合使用来进行冻土影响流域的洪水预报。
(2) 文章所述洪水预报方法在实际工作中得到较好的应用效果。实际使用时首先根据流域的值、降雨强度和降雨中心位置等参数, 在相关图上选择适宜的线型初步确定预报洪峰流量, 然后再选择预报方程法进行洪峰流量预报计算。方程计算值与查图值如果接近 (差值小于查图值的20%) 则取二者的均值作为正式预报结果;否则应进一步综合分析修正预报结果。
(3) 文章采用相关图法和预报模型法分别从经验相关和数理统计两个方面对示例流域的洪水预报方法进行了分析。由于有效预报因子较少, 加之对流域冻土特性规律掌握不够全面等原因, 预报方程历史拟合精度较低, 今后还有待对冻土影响流域的洪水预报方法进行深入探讨研究。
参考文献
[1]邓先俊.陆地水文学[M].北京:水利电力出版社, 1984:61-114.
[2]李慧珑.水文预报[M].北京:水利电力出版社, 1979:60-122.
实时洪水预报 篇8
(1) 水利工程、农田蓄放水误差。流域中, 常有许多中小型水利工程, 遇洪水, 先拦蓄洪水, 若长期洪水拦蓄不下, 又大量放水泄洪, 遇干旱、农业需水季节, 放水灌溉, 泄空库容。如此常给洪水预报带来大的误差, 其大小取决于流域内中小型水利工程的多少。
(2) 设备故障, 导致不合理的观测数据或资料缺测。 有许多雨量站和水位站, 在水文遥测系统的运行过程中, 常遇到各种故障, 给实时洪水预报带来误差。
(3) 水文规律简化误差, 即模型结构误差。 如农业活动作用的忽略、降雪作为降雨处理及产流机理简化为蓄满产流、超渗产流等, 都属于模型结构误差, 当与实际出入大时, 就会带来大的误差。
(4) 流域水文规律的变化。流域水文规律受气候条件和下垫面条件的改变而改变。如北方高寒地区融雪径流形成的洪水与暴雨型洪水的差异, 雷暴雨、台风雨引起的洪水特征与锋面雨引起的洪水特征差异等, 及人类活动如森林的大面积砍伐、开挖人工河渠、修建大型水库、水土保持影响等, 都会给水文规律带来大的影响, 也会给实时洪水预报带来一定的误差。
(5) 水文资料观测误差和雨量资料代表性误差。
2 预报模型的误差
清河水库洪水预报包括产流预报和汇流预报两部分。清河水库的产流预报现有4 种模型:即水箱模型、新安江模型、大伙房模型、蓄满产流模型进行洪水预报, 经多年实践经验, 蓄满产流模型预报精度比较高, 达到84%左右, 新安江预报精度为83%, 接近蓄满产流模型, 大伙房、水箱模型预报精度相对较低, 为75%左右;蓄满产流模型适合本流域的特点, 是清河水库产流预报一直采用的模型, 其他模型预报起到参考、对比的作用。汇流预报主要应用单位线法, 包括雨强单位线和检验单位线两部分[1]。
3 蓄满产流模型原理
清河水库流域蓄水量每年从5 月1 日开始计算, 5 月1日起土壤上下层初始含水量分别为Pa上=0 mm, Pa下=50 mm, 上、下层最大土壤含水量分别为Pam上=45 mm、Pam下=85 mm, 最大土壤含水总量130 mm。
P′=P+Pa上-45-E
R=P′- (85-Pa下) × (1-e-0.012×P)
其中, R─净雨 (mm) ;P─流域平均降雨 (mm) ;Pam上─上层土壤最大含水量 (mm) ;Pam下─下层土壤最大含水量 (mm) ;Pa上─ 降雨起始日上层土壤含水量 (mm) ;Pa下─降雨起始日下层土壤含水量 (mm) ;E─雨期蒸发量 (mm) 。
4 蓄满产流模型误差的解决方法
流域平均降雨量的误差, 直接影响洪水预报的精度, 清河水库面雨量是一种“以点代面”特殊方法计算得出的经推算值, 并不是真正意义的实测值。采用泰森多边形法 (垂直平分法) 来计算流域平均降雨量, 各雨量站点多少、控制面积的权重, 也就是雨量站点布设密度, 直接影响流域平均降雨量的精度, 雨量站点少, 必然会降低计算面雨量的精度;雨量站点多, 计算出的面雨量值代表性好, 精度高, 使流域平均降雨量更接近真值, 提高洪水预报的精度[2]。
净雨的误差, 清河水库产流预报采用蓄满产流模型, 在清河水库洪水预报中, 土壤含水量 (Pa) 分为上层土壤含水量 (Pa上) 和下层土壤含水量 (Pa下) 。 每年4 月末有大的降雨, 则从这场降雨开始计算Pa, 否则从每年5 月1 日开始连续计算, 降雨开始日Pa上起始值设为0 mm, Pa下起始值设为50 mm, 汛后就不再计算, Pa计算的不连续性, 可带来一定的误差, Pa初始值的设定直接影响净雨的准确性, 而每年的春汛来水量不同, 选用固定的Pa初始值显然不合理, 尤其在春汛偏枯、久旱无雨年份的情况下, 计算出本流域产生净雨的误差较大, 适当降低Pa的初始值, 使产生净雨趋近于合理值, 提高洪水预报的精度[3]。
清河水库汇流预报主要应用单位线法, 包括雨强单位线和检验单位线两部分, 雨强单位线只有“降雨均匀和中下游”与“降雨在上游”2 种分类。当暴雨中心在下游、库区附近时, 预报精度低, 尤其是暴雨强度大时, 预报误差很大, 不能满足要求, 如2010 年7 月21 日洪水, 针对降雨不均匀、库区降雨大的特点, 采用以下方法提高预报精度:一是尝试了用单站产流做汇流计算的方法, 效果很好, 使汇流预报精度有很大提高;二是汇流预报不单纯依靠预报软件, 灵活应用电子表格做汇流过程, 解决了产流预报净雨的实时修正, 从而达到每时段的汇流过程可灵活削减, 更贴近实际过程;三是有效利用历史洪水, 在汇流预报过程中, 对多个预报峰值、峰现时间的选判, 参考类似的历史洪水数据进行对比分析[4]。
经验单位线只有“降雨在中下游”“先小后大, 先均匀后下游”“降雨在中游均匀”3 类。一是经验单位线数量少, 二是极端天气下, 非典型降雨过程越来越多, 水库已有的单位线不可能满足汇流预报的要求, 影响预报的整体精度。经验单位线数量、类型的多少, 直接影响洪水预报的精度, 本流域内选择典型的历史洪水, 分析推求出各种类型的经验单位线, 汇流预报时有可供选择的相似单位线;若经验单位线不能满足雨型的要求, 可选取已有类似单位线, 在不改变单位线整体形状的前提下, 按比例将峰提前或滞后、 扩大或缩小, 该方法能使峰现时间、 洪峰流量的预报结果更接近真值, 进行汇流预报, 从而提高洪水预报的精度[5]。
5 结语
实时洪水预报中产生误差的原因极其复杂, 为了能够提供足够准确的洪水预报方案, 必须对实时洪水预报中可能出现的误差加以考虑。针对清河水库流域的实时洪水预报问题, 提出上述保证预报精度的方法和措施, 将误差尽量减到最小, 从而提高清河水库洪水预报的精度。
参考文献
[1]姜凤海.浅析如何提高清河水库洪水预报精度[J].内蒙古水利, 2015 (3) :133-134.
[2]胡环.清河水库洪水预报误差评定及解决办法[J].内蒙古水利, 2015 (1) :129-130.
[3]姜玉婷.误差相似性综合修正技术在洪水预报中的应用[J].东北水利水电, 2013 (4) :30-31.
[4]张艳平, 王国利, 彭勇, 等.考虑洪水预报误差的水库汛限水位动态控制风险分析[J].中国科学:技术科学, 2011, 41 (9) :1256-1261.