CN104179148B - 一种基于不可控洪水特性的洪水调度安全度评价方法 - Google Patents
一种基于不可控洪水特性的洪水调度安全度评价方法 Download PDFInfo
- Publication number
- CN104179148B CN104179148B CN201410177185.3A CN201410177185A CN104179148B CN 104179148 B CN104179148 B CN 104179148B CN 201410177185 A CN201410177185 A CN 201410177185A CN 104179148 B CN104179148 B CN 104179148B
- Authority
- CN
- China
- Prior art keywords
- flood
- starting
- reservoir
- point detection
- time
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Alarm Systems (AREA)
Abstract
本发明涉及一种基于不可控洪水特性的洪水调度安全度评价方法,其特征在于步骤包括:(1)根据预报入库洪水洪峰和水库实时水位判断入库洪水类型;(2)根据预报入库洪水过程和水库面临水位,按敞泄方式调度,计算滞蓄洪水的最快下泄时间;(3)定义洪水调度起调水位取值范围的区间,该区间分为相对区间和绝对区间;(4)分别计算洪水调度方案的动态相对安全度或动态绝对安全度,判断洪水调度方案对下游防洪或水库自身安全的影响;(5)选择是否采取预泄调度的方式,降低起调水位提高安全度;(6)若有修正信息,重复步骤(1)至(5)方案。本发明解决稀遇洪水及中小可控洪水在调度中的安全度评价,为解决洪水资源安全利用提供了支撑条件。
Description
技术领域
本发明属于水利水电工程技术领域,特别是涉及一种基于不可控洪水特性的洪水调度安全度评价方法。
背景技术
汛期水库洪水调度主要解决防洪安全问题,尤其是遭遇重现期较长的洪水,无法实现预见期长、精度高的洪水预报,工程管理者往往被动应对入库洪水,甚至呈现失控、误判状态。其根本原因是对洪水调度的基本原理、工程设计参数和稀遇洪水的特性等缺乏更深入的分析和研究,特别是不利因素组合条件下洪水调度安全水平评估分析方法,更是国内外研究和实践的空白。如对过高的不适当起调水位、工程泄洪设施突发故障等情况,再遭遇稀遇洪水如果不在事前进行全面、系统和有针对性的安全度量化评估,一旦发生这些情况,虽有可能侥幸涉险度汛,但更可能的结果是工程失事引发灾难。如何克服现有技术的不足已成为当今水利水电工程技术领域中亟待解决的重大难题之一。
发明内容
本发明的目的是为克服现有技术的不足而提供一种基于不可控洪水特性的洪水调度安全度评价方法,本发明建立以最快下泄时间为参数的不可控洪水调度安全度评价技术系统,既可解决稀遇洪水在调度过程中的安全度评价,又可解决中小可控洪水调度的安全度评价,为实现洪水资源安全利用提供了强有力的支撑条件。
根据本发明提出的一种基于不可控洪水特性的洪水调度安全度评价方法,其特征在于包括如下具体步骤:
(1)根据预报入库洪水洪峰和水库实时水位判断入库洪水类型;
(2)根据预报入库洪水过程和水库面临水位,按敞泄方式调度,计算滞蓄洪水的最快下泄时间;
(3)定义洪水调度起调水位取值范围的区间,该区间分为相对区间和绝对区间,其中:相对区间定义为起调水位在汛限水位与防洪高水位之间;绝对区间定义为起调水位在汛限水位与校核洪水位之间;
(4)分别计算洪水调度方案的动态相对安全度或动态绝对安全度,其中:计算洪水调度方案的动态相对安全度,并判断洪水调度方案对下游防洪安全的影响;计算洪水调度方案的动态绝对安全度,并判断洪水调度方案对水库自身安全的影响;
(5)选择是否采取预泄调度的方式,降低起调水位以提高安全度;
(6)若预报洪水有修正信息,再次重复上述步骤(1)至(5)的方案。
本发明进一步的优选方案是:
本发明步骤(1)所述的入库洪水类型包括可控洪水、临界洪水或不可控洪水,各入库洪水类型的判别条件包括:
a.当洪水起调水位对应的最大下泄能力比面临洪水洪峰流量大时,该洪水定义为可控洪水;
b.当起调水位对应的最大下泄能力恰好等于面临洪水洪峰流量时,该洪水定义为临界洪水;
c.当面临洪水的洪峰流量大于起调水位对应的最大下泄能力时,该洪水定义为不可控洪水。
本发明步骤(2)所述的最快下泄时间,是指在水库调洪至最高水位,执行敞泄调度方式下泄滞蓄洪水,使水位恢复至起调水位的时间,最快下泄时间的影响因素包括:
a.洪水量级影响:在同一起调水位,洪水量级愈大,则最快下泄时间愈长;
b.起调水位影响:同一洪水,起调水位愈高,其最快下泄时间愈长;
c.退水缓急影响:同样量级洪峰的洪水,退水缓的洪水其最快下泄时间比退水急的洪水的最快下泄时间要长。
本发明步骤(4)所述的动态相对安全度,是指对应步骤(3)所述的起调水位在相对区间时,确定洪水调度方案的安全程度,该动态相对安全度计算公式:
式中Tmax、Tmin、Ts分别为该洪水在Z防高、Z汛限、Zs水位起调对应的最快下泄时间,其中Z防高>Zs>Z汛限。
本发明步骤(4)所述的动态绝对安全度,是指对应步骤(3)所述的起调水位在绝对区间时,确定洪水调度方案的安全程度,该动态绝对安全度计算公式与述的动态相对安全度计算公式相同,其中Tmax为该洪水在校核水位起调的最快下泄时间。
本发明步骤(5)所述的选择是否采取预泄调度的方式,是指预泄调度对不可控洪水影响表现在两个方面:
a.预泄不改变洪水的不可控性;
b.预泄可缩短最快下泄时间。
本发明与现有技术相比其显著优点在于:第一,本发明通过对洪水及起调水位相互耦合关系的分析,结合水库库容特性和最大下泄能力特性,针对同一起调水位将洪水量化为不可控、临界、可控三个级别,使调度人员能明确地判断本次洪水调度的边界,以避免盲目和误判。第二,本发明提炼最快下泄时间(Tmin)这个核心参数,并用来表达水库防洪调度的安全水平,既能准确地反应洪水调度由已知条件决定的确定性特征(如洪水大小、退水缓急、下泄能力、库容曲线、起调水位等),又能表达后续洪水未知对水库防洪安全构成的不确定性特征(如后续洪水大小、发生时间等)。第三,本发明构建的防洪调度安全度概念及其计算模式涵盖不可控、临界、可控等各类洪水和任意水位起调的各种工况,极大地方便了调度人员对水库各种组合方案实施离线调度演练,为防洪调度决策提供系列化的成果。第四,本发明构建的理论概念及其计算模式还可推广应用到非正常状态下的调度模式如泄水建筑物故障导致的泄流能力变化、洪水预报失误产生的违规偏差以及下游例外的防洪要求等。本发明能够从根本上结束汛期洪水调度中大量弃水而汛后无水可用的状态,将我国清洁可再生的水电能源开发利用提高到新水平。
附图说明
图1为本发明提出的一种基于不可控洪水特性的洪水调度安全度评价方法的实施步骤方框示意图。
图2为由图2-1和图2-2组成的洪水起调水位与洪水可控性关系示意图;其中:图2-1为洪水起调水位确定,由预报入库洪水洪峰判别洪水可控性,即为可控、临界和不可控洪水定义示意图;图2-2为入库洪水确定,由起调水位判别洪水可控性,即为起调水位对洪水可控性影响示意图。
图3为由图3-1、图3-2、图3-3、图3-4和图3-5组成的最快下泄时间Ts示意图;其中:图3-1为最快下泄时间Ts的物理概念示意图;图3-2为洪水大小和Ts关系,即洪水大小对Tmin的影响示意图;图3-3为起调水位和Ts关系,即起调水位对洪水Tmin的影响示意图;图3-4为预泄对Ts影响关系,即预泄和洪水可控性关系示意图;图3-5为退水曲线对Ts影响关系,即退水缓急对Tmin的影响示意图。
具体实施方式
下面结合附图和实施例对本发明的具体实施方式作进一步的详细说明。
结合图1,本发明提出的一种基于不可控洪水特性的洪水调度安全度评价方法,包括如下具体步骤:
(1)根据预报入库洪水洪峰和水库实时水位判断入库洪水类型;所述的入库洪水类型包括可控洪水、临界洪水或不可控洪水,各入库洪水类型的判别条件包括:
a.当洪水起调水位对应的最大下泄能力比面临洪水洪峰流量大时,该洪水定义为可控洪水;
b.当起调水位对应的最大下泄能力恰好等于面临洪水洪峰流量时,该洪水定义为临界洪水;
c.当面临洪水的洪峰流量大于起调水位对应的最大下泄能力时,该洪水定义为不可控洪水。
图2为由图2-1和图2-2组成的洪水起调水位与洪水可控性关系示意图,其中包括:图2-1所示的洪水起调水位确定,由预报入库洪水洪峰判别洪水可控性,即可控、临界和不可控洪水定义,其中:Z0为洪水起调水位;洪水曲线①为可控洪水(q0m>Qm1)、洪水曲线②为临界洪水(q0m=Qm2)、洪水曲线③为不可控洪水(q0m<Qm3);图2-2所示的入库洪水确定,由起调水位判别洪水可控性,即起调水位对洪水可控性影响;当在洪峰Qm的洪水给定条件下:在Z″0水位起调,该洪水为不可控洪水(Qm>q″0m);在Z'0水位起调,该洪水为临界洪水(Qm=q'0m);在Z0水位起调,该洪水为可控洪水(Qm<q0m)。
(2)根据预报入库洪水过程和水库面临水位,按敞泄方式调度,计算滞蓄洪水的最快下泄时间;所述的最快下泄时间,是指在水库调洪至最高水位,执行敞泄调度方式下泄滞蓄洪水,使水位恢复至起调水位的时间,最快下泄时间的影响因素包括:
a.洪水量级影响:在同一起调水位,洪水量级愈大,则最快下泄时间愈长;
b.起调水位影响:同一洪水,起调水位愈高,其最快下泄时间愈长;
c.退水缓急影响:同样量级洪峰的洪水,退水缓的洪水其最快下泄时间比退水急的洪水的最快下泄时间要长。
图3为由图3-1、图3-2、图3-3、图3-4和图3-5组成的最快下泄时间Ts示意图,其中包括:图3-1所示的最快下泄时间Ts的物理概念;图3-2所示的洪水大小和Ts关系,即洪水大小对Tmin的影响关系,其中:Z0为确定的起调水位,当入库洪水Qm2>Qm1时,对应T2>T1;图3-3所示的起调水位和Ts关系,即起调水位对洪水Tmin的影响关系,其中:Qm为确定的入库洪水洪峰,当起调水位Z'0>Z0,对应T2>T1;图3-4所示的预泄对Ts影响关系,即预泄和洪水可控性关系,其中:Qm为确定的入库洪水洪峰,通过预泄,水库起调水位由Z0降至Z'0,对应T0'<T0;图3-5所示的退水曲线对Ts影响关系,即退水缓急对Tmin的影响关系,其中:Qm为确定的入库洪水洪峰,①洪水过程线退水急、②洪水过程线退水缓,对应T2>T1。
(3)定义洪水调度起调水位取值范围的区间,该区间分为相对区间和绝对区间,其中:相对区间定义为起调水位在汛限水位与防洪高水位之间;绝对区间定义为起调水位在汛限水位与校核洪水位之间。
(4)分别计算洪水调度方案的动态相对安全度或动态绝对安全度,其中:计算洪水调度方案的动态相对安全度,并判断洪水调度方案对下游防洪安全的影响;计算洪水调度方案的动态绝对安全度,并判断洪水调度方案对水库自身安全的影响;所述的动态相对安全度,是指对应步骤(3)所述的起调水位在相对区间时,确定洪水调度方案的安全程度,该动态相对安全度计算公式:
式中Tmax、Tmin、Ts分别为该洪水在Z防高、Z汛限、Zs水位起调对应的最快下泄时间,其中Z防高>Zs>Z汛限。
所述的动态绝对安全度,是指对应步骤(3)所述的起调水位在绝对区间时,确定洪水调度方案的安全程度,该动态绝对安全度计算公式与上述的动态相对安全度计算公式相同,其中Tmax为该洪水在校核水位起调的最快下泄时间。
(5)选择是否采取预泄调度的方式,降低起调水位以提高安全度;所述的选择是否采取预泄调度的方式是指预泄调度对不可控洪水影响,表现在预泄不改变洪水的不可控性和预泄可缩短最快下泄时间两个方面,其中:
a.预泄不改变洪水的不可控性;需要说明的是:不可控洪水一般都是较大的洪水,通常对大洪水存在预报难、预见期短或预报精度差,预泄对洪水可控性的影响是消极的,如图3-4所示,在Z0水位起调时洪水不可控,在T预代表的预见期内预泄,使水位降至Z0′,该洪水仍然不可控,即预泄无法将不可控洪水变为可控洪水。其根本原因是随着起调水位降低,水库最大下泄能力也降低了,因而导致水库在更低水位面临同一场不可控洪水。
b.预泄可缩短最快下泄时间:需要说明的是:如图3-4所示,若在Z0水位起调并遭遇不可控洪水,水库自b点开始,洪水进入不可控阶段,执行敞泄后仍被动滞蓄洪水,直至达到本次调洪最高水位Z1。在此水位仍执行敞泄,以实现库水位尽快恢复至起调水位Z0。其最快下泄时间如图示为T0。
而执行预泄,库水位降至Z0′。水库从c点开始面临洪水不可控阶段且遵循cd线敞泄,直至达到最高水位Z1′,再次敞泄恢复至水位Z0所需最快下泄时间为T0′。数值模拟计算结果显示,T0′<T0。即预泄可以缩短水库自调洪高水位降至起调水位时间,因此也相应提高了水库防洪安全水平。而图中abcd所围面积的洪水量恰好等于Z0Z0′之间预泄的库容,即预泄以后,自Z0′起调至d点,水库将达到Z0水位。预泄调度在Z0水位仅面临ad线以右的洪水;而不预泄在Z0水位起调,水库面临b点以右洪水,即不预泄比预泄多了一段(ab段)不可控洪水要滞蓄。
在退水段同样由于两者最高水位差异,预泄后最高水位Z1′低于不预泄的最高水位Z1。因此预泄比不预泄要少泄ef段洪水,结果使得T0′<T0。
(6)若预报洪水有修正信息,再次重复上述步骤(1)至(5)的方案。
现以本发明应用于长江流域东方水库水利工程为例,进一步说明本发明的具体实施方式。
一、基本资料:
a.库容水位关系(V~Z),详见附表1;
b.最大下泄能力水位关系(qm~Z),详见附表1;
c.典型洪水资料(p=0.01%),详见附表2;
d.初始条件:
汛限水位:Z0=1138m;起调水位:Zs=1138.6m;防洪高水位(Z防高)为1141m;校核洪水位(Z校核)为1145m;洪水取典型洪水0.8倍;下游无防洪要求。
二、应用实施例的实现过程:
步骤1,根据预报入库洪水洪峰和水库实时水位判断入库洪水类型。
由上述基本资料qm~Z关系表2中查得Zs=1138.6m、qm=6870m3/s;0.8倍典型洪水的洪峰流量为Qm=0.8×9020=7216m3/s;Qm>qm;得出洪水在(ZS)1138.6m起调水位时为不可控洪水。
步骤2,按敞泄方式调度,计算滞蓄洪水的最快下泄时间(Ts),分析最快下泄时间(Ts)的影响因素。其中:
计算最快下泄时间Ts=15.63h;
退水缓急对最快下泄时间(Ts)的影响:
设0.8倍校核洪水为入库洪水,即Qm=7216m3/s,洪水历时为167h,峰现时间设置三种退水缓急类型,起调水位(ZS)为1138.6m,Qm>qm,仍属不可控洪水,执行敞泄调度方式;对洪峰相同,历时相同,峰现时间不同的三个洪水分别得出其最快下泄时间(Ts),详见表3,由表1结果可见退水缓(峰现早)的洪水,其最快下泄时间(Ts)最长,退水愈缓,Ts愈长。
步骤3,定义洪水调度起调水位取值范围的区间分为相对区间和绝对区间。
相对区间定义为起调水位在汛限水位与防洪高水位之间,东方水库汛限水位(Z汛限)为1138m,防洪高水位(Z防高)为1141m;
绝对区间定义为起调水位在汛限水位与校核洪水位之间,东方水库汛限水位(Z汛限)为1138m,校核洪水位(Z校核)为1145m;
步骤4,分别计算洪水调度方案的动态相对安全度或动态绝对安全度。
一是计算动态相对安全度,得出:
Tmin=6.17h;Tmax=30.75h;α相对=0.608;
二是计算动态绝对安全度,得出:
Tmin=6.17h;Tmax=41.65h;Ts=15.63h;α绝对=0.734;
验证在1138.6m水位起调,遭遇0.8倍校核量级洪水,对水库自身安全的影响。经计算在汛限水位起调,遭遇校核洪水且采用敞泄调度方式,其最快下泄时间为20.72h,本调度方案最快下泄时间为15.63h,显然本方案对水库自身安全无影响。
步骤5,选择是否采取预泄调度的方式,降低起调水位以提高安全度。
若调度人员对预报洪水信息信心不足,采用更为保守偏安全的预泄调度方式,在预见期内将水库预泄至水位1138.2m,即Zs=1138.2m,查得该水位最大下泄能力qm=6650m3/s,入库洪水洪峰流量(Qm)为7216m3/s,Qm>qm,洪水仍然不可控;计算最快下泄时间(Ts)为9.86h,此时间比1138.6m水位起调对应的最快下泄时间15.63h为短、α绝对=0.896。
步骤6,预报洪水信息修正(修正洪峰)信息。
假设预报洪水洪峰(Qm)进行10%放大修正,即上述计算中0.8倍校核洪水同倍增大10%,起调水位(ZS)仍为1138.6m,得出Ts=19.26h,即洪水增大,同一起调水位最快下泄时间延长,安全度下降。
表1:东方水库水位——库容——最大下泄流量关系表
水位(m) | 1055 | 1076 | 1080 | 1085 | 1090 | 1095 | 1100 | 1105 | 1110 | 1112 |
库容(106m3) | 612 | 1136 | 1260 | 1430 | 1620 | 1820 | 2030 | 2270 | 2520 | 2628 |
下泄流量(m3/s) | 153 | 651 | 747 | 851 | 944 | 1028 | 1106 | 1179 | 1247 | 1273 |
水位(m) | 1114 | 1116 | 1118 | 1120 | 1122 | 1124 | 1126 | 1128 | 1130 | 1131 |
库容(106m3) | 2738 | 2848 | 2964 | 3080 | 3208 | 3336 | 3468 | 3604 | 3740 | 3814 |
下泄流量(m3/s) | 1299 | 1324 | 1349 | 1374 | 1398 | 1422 | 1445 | 1468 | 3190 | 3529 |
水位(m) | 1132 | 1133 | 1134 | 1135 | 1136 | 1137 | 1138 | 1139 | 1140 | 1141 |
库容(106m3) | 3888 | 3962 | 4036 | 4110 | 4187 | 4265 | 4332 | 4420 | 4497 | 4580 |
下泄流量(m3/s) | 3888 | 4263 | 4656 | 5065 | 5540 | 6032 | 6547 | 7075 | 7626 | 8184 |
水位(m) | 1142 | 1143 | 1144 | 1145 | 1145 | 1138 | 1139 | 1140 | 1141 | |
库容(106m3) | 4662 | 4745 | 4827 | 4910 | 4956 | 4332 | 4420 | 4497 | 4580 | |
下泄流量(m3/s) | 8765 | 9463 | 9965 | 10585 | 10633 | 6547 | 7075 | 7626 | 8184 |
表2:东方水库0.01%洪水过程线表
时间(h) | 0 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 |
流量(m3/s) | 1260 | 1350 | 1420 | 1490 | 1560 | 1740 | 1920 | 2110 | 2250 | 2500 |
时间(h) | 10 | 11 | 12 | 13 | 14 | 15 | 16 | 17 | 18 | 19 |
流量(m3/s) | 2830 | 3030 | 3230 | 3430 | 3480 | 3540 | 3590 | 3660 | 3740 | 3810 |
时间(h) | 20 | 21 | 22 | 23 | 24 | 25 | 26 | 27 | 28 | 29 |
流量(m3/s) | 3770 | 3740 | 3700 | 3810 | 3920 | 4030 | 4180 | 4330 | 4480 | 4640 |
时间(h) | 30 | 31 | 32 | 33 | 34 | 35 | 36 | 37 | 38 | 39 |
流量(m3/s) | 4800 | 4960 | 5000 | 5030 | 5060 | 5060 | 5060 | 5060 | 4940 | 4810 |
时间(h) | 40 | 41 | 42 | 43 | 44 | 45 | 46 | 47 | 48 | 49 |
流量(m3/s) | 4690 | 5210 | 5490 | 5760 | 6280 | 6800 | 7320 | 7420 | 7520 | 7680 |
时间(h) | 50 | 51 | 52 | 53 | 54 | 55 | 56 | 57 | 58 | 59 |
流量(m3/s) | 7960 | 8280 | 8610 | 8750 | 8880 | 9020 | 8830 | 8640 | 8450 | 8230 |
时间(h) | 60 | 61 | 62 | 63 | 64 | 65 | 66 | 67 | 68 | 69 |
流量(m3/s) | 8010 | 7800 | 7550 | 7430 | 7320 | 7190 | 7070 | 6940 | 6820 | 6690 |
时间(h) | 70 | 71 | 72 | 73 | 74 | 75 | 76 | 77 | 78 | 79 |
流量(m3/s) | 6570 | 6450 | 6330 | 6210 | 6090 | 5960 | 5840 | 5690 | 5540 | 5390 |
时间(h) | 80 | 81 | 82 | 83 | 84 | 85 | 86 | 87 | 88 | 89 |
流量(m3/s) | 5300 | 5200 | 5060 | 4870 | 4670 | 4470 | 4340 | 4210 | 4080 | 4030 |
时间(h) | 90 | 91 | 92 | 93 | 94 | 95 | 96 | 97 | 98 | 99 |
流量(m3/s) | 3980 | 3920 | 3900 | 3880 | 3850 | 3830 | 3810 | 3780 | 3760 | 3730 |
时间(h) | 100 | 101 | 102 | 103 | 104 | 105 | 106 | 107 | 108 | 109 |
流量(m3/s) | 3700 | 3680 | 3650 | 3620 | 3590 | 3570 | 3540 | 3510 | 3480 | 3400 |
时间(h) | 110 | 111 | 112 | 113 | 114 | 115 | 116 | 117 | 118 | 119 |
流量(m3/s) | 3320 | 3240 | 3160 | 3070 | 2990 | 2950 | 2920 | 2880 | 2850 | 2810 |
时间(h) | 120 | 121 | 122 | 123 | 124 | 125 | 126 | 127 | 128 | 129 |
流量(m3/s) | 2770 | 2730 | 2680 | 2640 | 2590 | 2540 | 2500 | 2440 | 2390 | 2340 |
时间(h) | 130 | 131 | 132 | 133 | 134 | 135 | 136 | 137 | 138 | 139 |
流量(m3/s) | 2300 | 2260 | 2220 | 2170 | 2140 | 2110 | 2090 | 2060 | 2030 | 1990 |
时间(h) | 140 | 141 | 142 | 143 | 144 | 145 | 146 | 147 | 148 | 149 |
流量(m3/s) | 1960 | 1930 | 1890 | 1860 | 1830 | 1800 | 1760 | 1730 | 1700 | 1670 |
时间(h) | 150 | 151 | 152 | 153 | 154 | 155 | 156 | 157 | 158 | 159 |
流量(m3/s) | 1640 | 1620 | 1600 | 1580 | 1550 | 1530 | 1510 | 1490 | 1470 | 1450 |
时间(h) | 160 | 161 | 162 | 163 | 164 | 165 | 166 | 167 | ||
流量(m3/s) | 1430 | 1410 | 1380 | 1370 | 1350 | 1340 | 1320 | 1310 |
表3:退水缓急对最快下泄时间影响结果情况表
起调水位(m) | 汛限水位(m) | 洪峰流量(m3/s) | 峰现时间(h) | 下泄时间(h) |
1138.6 | 1138 | 7216 | 94 | 18.79 |
1138.6 | 1138 | 7216 | 100 | 17.95 |
1138.6 | 1138 | 7216 | 106 | 17.08 |
注:表3中数据为持续时间相同峰现时间不同的洪水比较。
本发明的具体实施方式中凡未涉到的说明属于本领域的公知技术,可参考公知技术加以实施。
以上具体实施方式及实施例是对本发明提出的一种基于不可控洪水特性的洪水调度安全度评价方法技术思想的具体支持,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在本技术方案基础上所做的任何等同变化或等效的改动,均仍属于本发明技术方案保护的范围。
Claims (2)
1.一种基于不可控洪水特性的洪水调度安全度评价方法,其特征在于包括如下具体步骤:
(1)根据预报入库洪水洪峰和水库实时水位判断入库洪水类型;
(2)根据预报入库洪水过程和水库面临水位,按敞泄方式调度,计算滞蓄洪水的最快下泄时间;
(3)定义洪水调度起调水位取值范围的区间,该区间分为相对区间和绝对区间,其中:相对区间定义为起调水位在汛限水位与防洪高水位之间;绝对区间定义为起调水位在汛限水位与校核洪水位之间;
(4)分别计算洪水调度方案的动态相对安全度或动态绝对安全度,其中:计算洪水调度方案的动态相对安全度,并判断洪水调度方案对下游防洪安全的影响;计算洪水调度方案的动态绝对安全度,并判断洪水调度方案对水库自身安全的影响;
(5)选择是否采取预泄调度的方式,降低起调水位以提高安全度;
(6)若预报洪水有修正信息,再次重复上述步骤(1)至(5)的方案;
其中:步骤(1)所述的入库洪水类型包括可控洪水、临界洪水和不可控洪水,各入库洪水类型的判别条件包括:
a.当洪水起调水位对应的最大下泄能力比面临洪水洪峰流量大时,该洪水定义为可控洪水;
b.当起调水位对应的最大下泄能力恰好等于面临洪水洪峰流量时,该洪水定义为临界洪水;
c.当面临洪水洪峰流量大于起调水位对应的最大下泄能力时,该洪水定义为不可控洪水;
步骤(2)所述的最快下泄时间是指水库调洪至最高水位,执行敞泄调度方式下泄滞蓄洪水,使水位恢复至起调水位的时间,最快下泄时间的影响因素包括:
a.洪水量级影响:在同一起调水位,洪水量级愈大,则最快下泄时间愈长;
b.起调水位影响:同一洪水,起调水位愈高,其最快下泄时间愈长;
c.退水缓急影响:同样量级洪峰的洪水,退水缓的洪水的最快下泄时间比退水急的洪水的最快下泄时间要长;
步骤(4)所述的动态相对安全度是指对应步骤(3)所述的起调水位在相对区间时,确定洪水调度方案的安全程度,该动态相对安全度计算公式为:
式中Tmax、Tmin、Ts分别为该洪水在Z防高、Z汛限、Zs水位起调对应的最快下泄时间,其中Z防高>Zs>Z汛限;
步骤(4)所述的动态绝对安全度是指对应步骤(3)所述的起调水位在绝对区间时,确定洪水调度方案的安全程度,该动态绝对安全度计算公式与上述的动态相对安全度计算公式相同,其中Tmax为该洪水在校核洪水位起调对应的最快下泄时间。
2.根据权利要求1所述的一种基于不可控洪水特性的洪水调度安全度评价方法,其特征在于步骤(5)所述的选择是否采取预泄调度的方式是指预泄调度对不可控洪水影响表现在以下两个方面:
a.预泄不改变洪水的不可控性;
b.预泄可缩短最快下泄时间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410177185.3A CN104179148B (zh) | 2014-04-29 | 2014-04-29 | 一种基于不可控洪水特性的洪水调度安全度评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410177185.3A CN104179148B (zh) | 2014-04-29 | 2014-04-29 | 一种基于不可控洪水特性的洪水调度安全度评价方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104179148A CN104179148A (zh) | 2014-12-03 |
CN104179148B true CN104179148B (zh) | 2016-02-03 |
Family
ID=51960468
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410177185.3A Expired - Fee Related CN104179148B (zh) | 2014-04-29 | 2014-04-29 | 一种基于不可控洪水特性的洪水调度安全度评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104179148B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106320259B (zh) * | 2016-09-06 | 2018-07-17 | 长江水利委员会长江科学院 | 一种兼顾排沙的三峡水库汛期沙峰调度方法 |
CN106320258B (zh) * | 2016-09-06 | 2018-07-17 | 长江水利委员会长江科学院 | 一种三峡水库汛期“蓄清排浑”动态使用的泥沙调度方法 |
CN106592504B (zh) * | 2016-12-21 | 2018-07-20 | 国网福建省电力有限公司 | 一种根据短期降水预报计算水电站入库洪水量的方法 |
CN110852525B (zh) * | 2019-11-19 | 2023-04-18 | 长安大学 | 基于预报误差越级设防的水库防洪限制水位动态调节方法 |
CN112287528B (zh) * | 2020-10-16 | 2021-07-20 | 黄河勘测规划设计研究院有限公司 | 一种基于高滩高槽的多沙河流水库防洪高水位确定方法 |
CN113177695A (zh) * | 2021-04-07 | 2021-07-27 | 武汉艾信科技有限公司 | 一种洪水调度方案评价与选取方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102776869A (zh) * | 2012-07-12 | 2012-11-14 | 贵州乌江水电开发有限责任公司 | 一种梯级水电站群的无预报洪水调度方法 |
CN103205950A (zh) * | 2013-04-09 | 2013-07-17 | 河海大学 | 一种保障下游通江湖泊供水安全的水库调控方法 |
CN103646294A (zh) * | 2013-11-12 | 2014-03-19 | 国网电力科学研究院 | 一种考虑闸门操作规则的水电站洪水优化调度方法 |
CN103729550A (zh) * | 2013-12-18 | 2014-04-16 | 河海大学 | 基于传播时间聚类分析的多模型集成洪水预报方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001098522A (ja) * | 1999-09-28 | 2001-04-10 | Tokyo Metropolis | 河川氾濫防止システム |
-
2014
- 2014-04-29 CN CN201410177185.3A patent/CN104179148B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102776869A (zh) * | 2012-07-12 | 2012-11-14 | 贵州乌江水电开发有限责任公司 | 一种梯级水电站群的无预报洪水调度方法 |
CN103205950A (zh) * | 2013-04-09 | 2013-07-17 | 河海大学 | 一种保障下游通江湖泊供水安全的水库调控方法 |
CN103646294A (zh) * | 2013-11-12 | 2014-03-19 | 国网电力科学研究院 | 一种考虑闸门操作规则的水电站洪水优化调度方法 |
CN103729550A (zh) * | 2013-12-18 | 2014-04-16 | 河海大学 | 基于传播时间聚类分析的多模型集成洪水预报方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104179148A (zh) | 2014-12-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104179148B (zh) | 一种基于不可控洪水特性的洪水调度安全度评价方法 | |
CN101794495B (zh) | 洪水预报系统中实时校正模型的优选方法 | |
US11236867B2 (en) | Method for detecting pipe burst in water distribution systems based on pressure disturbance extraction | |
CN103198235B (zh) | 基于风功率纵向时刻概率分布的风功率预测值预评估方法 | |
CN101625732B (zh) | 江河潮水水位的预测方法 | |
CN104901306B (zh) | 一种考虑连锁故障的电网运行安全裕度计算方法 | |
US11550976B2 (en) | Node flow optimization distribution method for improving accuracy of transient hydraulic simulation of water supply in-series pipeline | |
CN107067127B (zh) | 一种电力系统高风险连锁故障事故链搜索方法 | |
CN105426971A (zh) | 基于混沌优化bp神经网络模型的江河涌潮短期预报方法 | |
CN103207948A (zh) | 基于风速相关性的风电场测风仪风速缺失数据插补方法 | |
CN104615845A (zh) | 一种超临界洪水风险分析方法 | |
Dai et al. | Flood‐risk analysis based on a stochastic differential equation method | |
Raso et al. | Optimal and centralized reservoir management for drought and flood protection on the Upper Seine–Aube river system using stochastic dual dynamic programming | |
CN105449667A (zh) | 一种发输电系统可靠性预测方法 | |
CN103473453B (zh) | 梯级枢纽联合调度中水流滞时问题处理方法 | |
CN103088783B (zh) | 一种面向生态的水库调度方案生成方法 | |
CN103366099A (zh) | 一种水文模型参数调试方法 | |
CN105067935B (zh) | 基于无功不平衡量的变压器分接头参数估计方法 | |
CN104578050A (zh) | 一种电网暂态稳定强相关输电断面的识别方法 | |
Valamanesh et al. | Probabilistic modeling of joint hurricane-induced wind and wave hazards to offshore wind farms on the Atlantic coast | |
CN104037803A (zh) | 区域风电场群时序电量的统计与分析方法 | |
CN111028487B (zh) | 一种水处理监测方法及系统 | |
Ma et al. | Study on flood control risk of flood control engineering system based on the clustering of measured data | |
CN104866974A (zh) | 一种考虑预安排停电的配电网可靠性评估设备 | |
Hua et al. | Preliminary study on dynamic control of flood limited level of reservoir |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160203 |