计算水库运用设计洪水必须首先分析水库运用面临径流年内随机过程概率规律,并考虑水库调度决策特点、要求。本文为水库运用设计洪水预报概率模型提供了理论框架。并以三峡水库为例予以验证。结果表明,概率模型给出的三峡水库动态设计洪水预报是合理的,可靠的。<br/>To calculate the design flood for reservoir operation, we first need to analyze the regularity of annual runoff process’s stochastic probability of the reservoir, and take the characteristics and the requirements of reservoir operation decision into account. A theoretical framework for develop-ment of design flood forecast probability model has been provided in this paper and validated with the Three Gorges Reservoir. The results indicate that the dynamic design flood forecast of the Three Gorges Reservoir based on this probability model is reasonable and reliable.
王善序
长江水利委员会水文局,湖北 武汉
Email: wangsxabc@sina.com
收稿日期:2015年2月5日;录用日期:2015年2月12日;发布日期:2015年2月16日
计算水库运用设计洪水必须首先分析水库运用面临径流年内随机过程概率规律,并考虑水库调度决策特点、要求。本文为水库运用设计洪水预报概率模型提供了理论框架。并以三峡水库为例予以验证。结果表明,概率模型给出的三峡水库动态设计洪水预报是合理的,可靠的。
关键词 :水库运用设计洪水,径流年内不平穏随机过程,动态设计洪水
我国已建成数百座大型水库。但是,由于各种原因,我国的水库调节程度较低,而来水大多集中于汛期。水库汛期运用至为主要。传统设计规则的局限和相关技术研究的滞后,“我国的水库在工程设计洪水和调度运用中都存在一些值得进一步研究的问题”[见:水库设计运用专题研究招标文件,国家防汛抗旱总指挥部办公室、水利部水利水电规划设计总院,2002年4月]。其中一个重要问题是:水库在汛期经常控制水位过低,弃水较多、蓄水不多,汛后又常不能及时蓄满,严重制约了水库综合利用效益的发挥。显然,造成这一现象的直接原因是水库汛限水位设置过低。而后者与水库运用仍采用年最大设计洪水有关。于是,水库运用设计洪水成了优化水库运用的技术瓶颈。“水库运用采用年最大设计洪水是不合理的”逐渐成为业内人士的共识。但为什么是不合理的?怎样的设计洪水才是合理?仍众说纷纭。
水库是一种高风险产业。在汛期,它随时可能遭遇特大洪水袭击。水库运用需要水文学提供未来一段时间内的水情信息。现代水文预报技术发展为改善水库运用调度提供了重要信息。遗憾的是预见期过短,通常无法覆盖水库调洪时段。水库运用调度面临决策时段的水情仍是不确定的。为保证防洪安全,仍需要设计洪水。按[
现在的问题是:被人们使用了一个多世纪,作为水库防洪设计重要依据的年最大设计洪水为什么如今在水库运用时会引起人们的不满?水库运用究竟应采用什么样的设计洪水才既能保证水库防洪安全又能提高兴利效益呢?但有一点是很显然的,只有能反映水库运用面临的实时径流演变概率规律的设计洪水才可能是合格的水库运用设计洪水。本文的目的是为水库运用寻求合格的设计洪水。
早在上世纪中叶,我国水文界就关注水库运用设计洪水问题。有人提出以汛期分期设计洪水取代年最大设计洪水,以考虑汛期洪水“季节性”演变规律。并被少数水库采用。进入本世纪以来,由于官方的重视,汛期分期设计洪水得到了广泛的研究和应用。汛期洪水演变非常复杂。“季节性”演变仅是汛期洪水演变的一个侧面,以它来反映汛期设计洪水变化显然是不全面的。何况水库运用调度也非仅按季节调度。尤其是它是否能反映汛期实时径流演变概率规律?回答显然是否定的。
本文将从分析水库运用面临实际径流随机过程概率规律出发,并考虑水库实时运用特点和要求,研究水库运用设计洪水预报概率模型。并以长江三峡水库为例予以验证
河流径流过程受众多因素影响、非常复杂,通常可视作为以年为周期的随机过程。大量资料分析表明,它的以年为时间间隔的各种时间系列,如年、月径流量,年最大洪水(如最大洪峰流量,最大时段洪量)等是相互独立的;据此估计的统计特征值,只要系列足够长,是穏定的。说明径流多年(年间)过程是平穏的。因此,可以视径流年系列为简单样本。只要采用单变量分布或频率就可以完整地表述径流多年(或年间)变化过程的概率规律。
设以X表示径流,则其频率可写成
其中x为一正实数。它们的统计特征,如均值
和方差
都与时间无关。式中
设与防洪标准对应的频率为pd,则设计洪水xd可由
求得。
这些看似简单的表达式构成了传统水文设计中估计年最大设计洪水频率分析的理论基础。世界各国水文工程师在此基础上,尽管采用的具体计算方法,如抽样、频率曲线线型及其参数估计方法等可能不同,只要系列较长,估计的年最大设计洪水却不会有显著差别。
一个多世纪来,洪水频率分析虽屡遭非议,在我国甚至多次险被“打倒”。但随着观测资料累积,越来越多的人接受了它。洪水频率分析已成为不可取代的水库设计洪水主要计算方法。其根本原因是:频率分析成果能反映水库设计所涉及径流多年过程的概率规律。
然而,与年间过程不同,径流年内过程却要复杂得多。地球上大多数河流在一个水文年内有枯期和汛期之分。即使在汛期也不是天天有汛、有洪水期和非洪水期(平水期)之别,在洪水期,涨水段和退水段也不一样。径流年内过程变化多端。它是一种前后时段径流相关、统计特征随时间而变的不平穏随机过程[
设以
来描述[
及其参数,如均值函数
和
显然,一维频率不足以描述汛期径流过程。如为了研究它任意两个时间径流变量之间的关系,需要引入二维频率
与此有关还要引入协方差函数
显然有
上述诸式为分析径流年内不平穏随机过程概率规律提供了基本数学工具。
由此可见,径流随机过程具有多重特性。径流年间过程和年内过程是不同性质的随机过程。与它们相关联的,估计设计洪水的概率模型也应是不同的。
在传统水文设计中,只论及或研究了在水库设计阶段要求的、与径流多年(即年间)过程相关联的年最大设计洪水。而与径流年内过程相关联的、不平穏过程的设计洪水计算问题迄今仍是空白。
迄今,水库设计洪水通常指的是在水库规划设计阶段估计的年最大设计洪水。在水库设计阶段,设计者需要保证他们设计的水库在建成后长期运行中的防洪安全[
然而,与设计阶段不同,水库运用实时防洪调度的特点是在既成工程条件下,通过优化调度决策发挥工程效益。而且,“每次只要作出一个时段(面临时段)的决策”[
设在水库运用决策时间t,这时水库当时及前期入库流量是已知的。未来决策时段t,t + dt径流量的条件频率可写成
其中,x0表示当前和前期入库流量状态,笔者称之为前期影响流量或水量;为决策时段长。则决策时段设计洪水xd可由条件频率
求得。式中pd为与防洪标准、决策时间和时段长有关的频率。
以上各式为径流年内不平穏随机过程的频率分析和水库运用设计洪水计算概率模型奠定了基础。
近20多年来,水文学中多变量频率分析已有了长足进展。现在,已有可能通过多变量频率分析估计水库运用决策时段设计洪水。
由式(13)可知,水库运用决策时段设计洪水的基本特点是:a) 随时间推移而变,即它是时变的;b)与当时和前期入库水情有关,是动态的。笔者称之谓“动态设计洪水”。
在这一小节,笔者将应用上节提出的水库运用动态设计洪水概率模型,研制、分析三峡水库运用动态设计洪水。
三峡水库控制长江上游、集水面积100万平方公里,是一座以防洪为主,兼有发电、航运等综合利用水库。
长江上游暴雨洪水集中在每年6~9月。雨季开始时间一般从南往北,从东向西逐步推进。南岸乌江6月中下旬可能出现较大洪水,到7月嘉陵江进入汛期,川西岷、沱江及金沙江要到7月下旬、8月才先后入汛。年最大洪水大多出现在7~8月份。入秋后,长江上游也常出现较大暴雨洪水。一般认为,7月上旬至9月上旬是长江上游洪水最盛时期
根据多年分析,长江上游洪水一般由多次暴雨形成。宜昌站汛期流量过程线形状大体可分单峰型,双峰型和多峰型三种类型。单峰型一般为大水年。在实测系列中,最著名的是1954年和1998年,由多次连续洪水先后迭加而成。其最大流量多出现在7月下旬至8月中旬。双峰型即所谓前后期洪水。前峰多出现在7月上中旬,后峰多出现在9月上中旬。多峰型一般为小水年,每次洪水都很小。
根据文献[
根据宜昌站日平均流量系列,计算了宜昌站汛期主要统计特征,以分析长江上游汛期径流过程的特点。图1~图4分别为宜昌站汛期7日水量相关系数,均值,变差系数(cv)和偏差系数(cs)图。它们反映了宜昌站汛期径流过程的一般统计特征。
图1. 宜昌站汛期前后7日水量相关系数图
图2. 宜昌站汛期7日水量平均值 (水量单位:亿立方米)
图3. 宜昌站汛期7日水量变差系数
图4. 宜昌站汛期7日水量偏差系数
这些图显示:宜昌汛期径流过程明显相关。7日水量相关系数普遍大于0.5,而且,入汛后有逐步增大趋势,至9月底竟接近0.8。在三个最重要的统计参数中,均值随时间变化最为明显。它入汛后迅速增大,至7月中下旬达到峰值,然后又缓慢递减。它们显示,宜昌站汛期径流统计特征随时间而变。
由此可见,三峡水库运用所面临的长江上游径流年内过程是前后时段水量相关、统计特征随时间而变的不平穏随机过程。
根据上述有关水库运用设计洪水概率模型,采用宜昌站1877~2000年日平均流量系列,笔者制作了三峡水库动态设计洪水预报概率模型。现应用此模型,对三峡水库动态设计洪水作一简短分析。
三峡水库动态设计洪水的时变规律可用它的时变设计洪水表达。时变设计洪水是笔者于本世纪初提出的一种水库运用设计洪水[
应该指出的是,上述三峡水库设计洪水汛期时变性反映的是一种“气候性”的时变规律。在任何一个实际年汛期,由于实际入库水情的影响,当年实际设计洪水变化与之仍会有相当的差异。
三峡水库运用动态设计洪水还与入库水情状态有关。
图6和图7分别是1942年笔1998年三峡水库汛期1000年一遇相应7日设计水量及据此推算的动态汛限水位预报成果(假定下泄流量为50,000立方米/秒)。1942年和1998年分别是宜昌站观测系列中有名的枯水年和大水年。这二年汛期水情相差特大,动态设计洪水–动态汛限水位也大相径庭。反映了入库水情对动态洪水–动态汛限水位的影响。
图7还显示,当入库流量接近、达到高水(约等于、大于50,000立方米/秒)时,动态设计洪水逼近、等于年最大7日设计洪量。说明,动态设计洪水能保证防洪安全。在此过程中,动态汛限水位则快速下行靠近设计阶段确定的汛限水位。只是动态7日设计水量与年最大设计洪水过程线是有差别的,而且用的调洪演算的方法也不同,此处动态汛限水位高于145米。
这些图中还可看到:在一次洪水从起涨、涨至洪峰,再退水消落的过程中,动态设计洪水与流量过
图5. 三峡水库同频率7日时变设计水量(水量单位:亿立方米)
图6. 三峡水库1942年汛期1000年一遇相应7日设计水量及汛限水位
程演变相向而行;而动态汛限水位则背向而行。
2012年7月长江上游曾出现三峡水库建成后最大一次洪峰,最大流量近70,000立方米/秒。在图8中给出了这次洪水过程中、逐日计算的三峡水库未来7日洪灾风险率,或频率演变过程。其间,入库洪水出现二次洪峰,与之对应,洪灾风险率(7日水量大于480亿立方米的概率)也出现二次峰。说明水库运用面临的实际径流过程中,洪灾风险率也是动态的。这也符合人们的经验。
在水文学中,传统的设计观念认为,对应于同一防洪标准,设计洪水是唯一的。显然,这种观点仅当设计涉及的是径流多年演变过程时才正确。然而,对于水库运用,它涉及径流年内不平穏随机过程,对应同一防洪标准的设计洪水,汛限水位和洪灾风险率都随时间和入库水情而变。当一次洪水起涨,预示洪灾风险来临、并随洪水增大而增加,设计洪水也随之增大;反之,洪峰过后退水开始,洪水渐消退,洪灾风险、设计洪水随之渐减。这些说明,动态设计洪水是符合客观实际径流演变规律,因而能为水库科学处理防洪与兴利矛盾提供依据。
在分析设计洪水合理性和可靠性时,人们首先会关注洪水经验点据与洪水频率曲线间的拟合情况。
图7. 三峡水库1998年汛期1000年一遇相应7日设计水量及汛限水位
图8. 三峡水库2012年7日水量超过480亿立方米风险率图
类似地,在考察动态设计洪水的合理性和可靠性时,自然也会关心条件频率曲线(如式(13)所示)与经验点据的拟合情况。不过,应该指出的是:首先,与一般洪水频率分析不同,这里的条件频率曲线不是通过与经验点据拟合而得的,而是根据模型计算的;其次,也与一般洪水频率分析不同,这里的经验点据是近似的。这是因为,在历史系列中,前期水情状态与当前的前期水情相同的几乎是没有的。只能在一定的允许误差范围内选取经验点据;所以最后,经验点据的个数一般少于观测系列长。而不同时间的经验点据的个数一般也是不相同的。
图9~图14反映了2012年从6月中下旬入库流量起涨、至7月下旬出现最大洪峰(期间出现二次洪峰),此后整个8月持续退水(详见图16中入库流量过程)。上述过程中各代表性时间条件频率曲线与相应近似经验点据拟合情况如下。图8~图11四图代表从汛初起涨后经二次洪峰至7月下旬达年最大洪峰阶段。期间条件频率曲线是上行的,虽二峰间曾略有回落如图10。图12、图13则代表退水后,7日水量条件频率曲线逐渐下行。
由图可见:概率模型计算的三峡水库7日水量条件频率曲线与近似经验点据的拟合是良好的。这意味三峡水库动态设计洪水概率模型能较好模拟长江上游汛期径流年内过程概率规律。由此估计的三峡水库动态设计洪水是合理、可靠的。
图9. 三峡水库2012年07月01日~07月07日7日水量条件频率曲线(经验点对应的前期影响流量误差≤±10%),制作日期:06月30日
图10. 三峡水库2012年07月08日~07月14日7日水量条件频率曲线(经验点对应的前期影响流量误差≤±10%),制作日期:07月07日
图11. 三峡水库2012年07月20日~07月26日7日水量条件频率曲线(经验点对应的前期影响流量误差≤±10%),制作日期:07月19日
图12. 三峡水库2012年07月26日~08月01日7日水量条件频率曲线(经验点对应的前期影响流量误差≤±10%),制作日期:07月25日
图13. 三峡水库2012年08月11日~08月17日7日水量条件频率曲线(经验点对应的前期影响流量误差≤±10%),制作日期:08月10日
在估计了动态设计水量后,水库调度决策者就可按当时下游容许安全泄量得到决策时段所需防洪库容。由于防洪高水位是确定的,利用库容曲线即得决策时段相应的动态汛限水位。
自2010年以来,笔者利用三峡公司网页公布的逐日入库流量资料,进行三峡水库逐日1000年一遇相应7日动态设计洪水–动态汛限水位试验预报。图15~图19是最近五年的试验结果(其中2010年资料不完整)。由这些图可见,三峡水库现行预报调度是比较保守的,有待优化的空间是比较大的。
这些试验结果表明,本文上述有关水库运用设计洪水–汛限水位的论述是合理的。
通过以上各节叙述可以得出:
水库汛期运用的基本问题是处理防洪与兴利矛盾。为保证水库防洪安全水库必须随时预留有足够的
图14. 三峡水库2012年08月20日~08月26日7日水量条件频率曲线(经验点对应的前期影响流量误差≤±10%),制作日期:08月19日
图15. 三峡水库2010年汛期逐日1000年一遇相应动态汛限水位预报(当下泄流量 = 50,000立方米/秒时),制作日期:08月11日
图16. 三峡水库2011年汛期逐日1000年一遇相应动态汛限水位预报(当下泄流量 = 50,000立方米/秒时),制作日期:10月10日
图17. 三峡水库2012年汛期逐日1000年一遇相应动态汛限水位预报(当下泄流量 = 50,000立方米/秒时),制作日期:10月10日
图18. 三峡水库2013年汛期逐日1000年一遇相应动态汛限水位预报(当下泄流量 = 50,000立方米/秒时),制作日期:10月09日
图19. 三峡水库2014年汛期逐日1000年一遇相应动态汛限水位预报(当下泄流量 = 50,000立方米/秒时),制作日期:09月30日
防洪库容以调蓄可能发生的设计洪水。所以,水库运用欲合理处理防洪与兴利矛盾首先应合理采用设计洪水,使之反映实际径流过程概率规律。
河流径流随机过程具有多重概率规律或特性。径流多年过程是平穏随机过程。它的统计特征不随时间而变。但径流年内过程却是不平穏随机过程,它的统计特征随时间而变。研究水库设计洪水应该首先认清设计所涉及径流随机过程的概率规律或特性。
水库设计和运用涉及不同的径流随机过程。在水库设计阶段为保证水库长期运用防洪安全、需要确定水库防洪能力。设计洪水计算需分析径流多年演变概率规律。水库运用是在既定工程规模条件下,通过调度决策提高水库兴利。它面临实时径流过程,且只对未来一个时段制定运行决策。这时,就需要研究径流年内不平穏随机过程概率规律、计算设计洪水,以便在保证防洪安全前提下提高兴利效益。显然,年最大设计洪水,汛期分期设计洪水都不是合格的水库运用设计洪水。
本文根据随机过程理论,为构建径流年内不平稳随过程频率分析概率模型提供了理论框架。以构建相应的水库运用动态设计洪水预报概率模型。传统的水文设计认为,一个防洪标准对应一个设计洪水。这只有当设计涉及的是多年径流过程时才正确。在水库运用面临的不平穏过程,洪灾风险率、动态设计洪水都是时变的并和入库水情有关。这时,设计洪水也一定是时变的,并和入库水情有关。水库在汛期只有按动态设计洪水预留防洪库容(设置汛限水位)才能使水库顺应自然径流演变规律、合理处理防洪与兴利矛盾、既保证防洪安全,又能提高兴利。
总之,年最大设计洪水是为水库防洪能力设计服务的。动态设计洪水是为运用决策设定汛限水位而作。动态设计洪水模型能在任何决策时刻提供未来决策时段内动态设计洪水-汛限水位预报。作为水库设计洪水组成部分,两者皆不可或缺、不可替代。
三峡水库年内径流过程分析和动态设计洪水-汛限水位试验预报成果验证了上述论述。并证明三峡水库动态设计洪水是合理可靠的。
笔者还制作了长江流域大渡河瀑布沟,清江长阳站,汉江白河等水库(站)的动态设计洪水概率模型,都获得了良好成果。
随着大批水库建成运行,用好管好已建成水库已成为重大的科技、经济和社会问题。相当长时期来,人们重视观测,预报和优化技术应用,但对水库运用设计洪水未必予以足够重视。近20年来,笔者从研究汛期分期设计洪水到提出动态设计洪水,才逐步认识到传统的,以平穏随机过程为分析对象的频率分析已不能满足水库运用设计洪水计算的要求。设计洪水已成为优化水库运用的一大技术瓶颈。发展或开发年内径流不平穏随机过程的频率分析-动态设计洪水预报技术已成水文学当务之急。
水库设计洪水估计是水文行业中最具争议、甚为艰难的领域。多年来也少有深入研究、缺乏创新。动态设计洪水为笔者在数年前才提出,迄今了解者仍少。鲁迅先生说过:路本来是没有的,走的人多了也就成了路。希望有更多的水文学人来研究水库设计洪水新领域——动态设计洪水。
不当之处谨请批评指正!
王善序, (2015) 论水库运用设计洪水Design Flood for Reservoir Operation. 水资源研究,01,9-22. doi: 10.12677/JWRR.2015.41002