《1 概述》
1 概述
随着水电站规模的不断扩大, 工程设计标准也相应提高, 泄洪设施规模也随之提高, 以确保大坝安全。水电站所选择的防洪设计标准较高, 这是必要的, 但稀遇洪水出现的机遇很少, 如为一工程设计的运行寿命期内都可能不遇的洪水而统一高标准的修建泄洪设施, 势必会使部分泄洪设施闲置或得不到充分利用, 造成经济损失。此时, 在确保大坝安全的前提下, 把泄洪设施分为正常泄洪设施和非正常泄洪设施, 以不同的设计标准进行修建是十分必要的。在常遇洪水的情况下, 启用正常泄洪设施即可确保大坝安全, 而在稀遇洪水时, 启用正常泄洪设施和非常泄洪设施, 并允许非常泄洪建筑物的自身部分破坏, 来确保大坝的安全, 这样可降低非常建筑物的设计标准, 降低工程造价。
目前, 国际上对大坝泄洪建筑物布置选型有简化的趋势, 同时尽量考虑采用非常溢洪设施, 以减少工程投资。采用既能保持一定泄洪可靠性, 又能减少常设泄洪设施的方案, 是高坝泄洪枢纽优化布置的方向。
采用非常泄洪设施对大坝泄洪风险会带来一定影响, 科学地评估这一影响是水电站采用非常泄洪设施的重要依据。通常水电站枢纽工程是根据在使用年限内可能遇到的自然事件来设计的, 其中包括洪水的估算和设计标准及其重现期的确定。通过水库调洪演算, 可确定水电站枢纽工程泄洪建筑物的规模。传统的调洪演算是按照设计洪水过程和水库调度方案, 根据水量平衡原理, 推算出库水位变化过程, 以确定水库的防洪特征水位及相应的泄洪建筑物型式、尺寸、高程等。然而由于在确定泄洪建筑物规模时, 存在着许多不确定性因素, 如入库洪水过程的水文条件不确定性、大坝泄洪能力的水力条件不确定性、库容与水位关系的不确定性和防洪起调水位的初始条件不确定性等, 计算的许多物理量均为随机量, 从而导致了库水位变化过程和泄水建筑物泄洪能力的随机性, 这些随机因素的存在增加了按照传统方法确定的泄洪设施规模泄洪能力风险, 降低了大坝的可靠度。因此有必要考虑这些随机量对大坝安全可靠度的影响。
本文引入随机微分方程, 采用泄洪风险分析方法评估水电站采用非常泄洪设施对大坝防洪安全所带来的影响, 并结合某一工程设置两条非常泄洪洞的可能性进行论证, 分析调洪方案, 确定两条非常泄洪洞启用的条件及相关的泄洪风险, 为导流洞改建方案的确定提供依据。
《2 泄洪风险分析方法》
2 泄洪风险分析方法
《2.1调洪过程随机微分方程》
2.1调洪过程随机微分方程
传统的调洪计算过程是确定性的, 直接利用水库的水量平衡原理建立计算方程, 考虑初始条件, 则为
式中Z (t0) 为库水位, Z0为初始库水位, Q (t) 为t时刻入库洪水流量, q (Z, m) 为相应时刻出库泄流量。
(1) 式为确定性方程, 实际上不同水位的库水面面积F (Z) 、洪水过程Q (t) 以及泄水建筑物过流流量系数m的测量均存在一定的随机性, 是个随机过程, 可认为在过程的不同时刻它们的概率密度符合正态分析
式中
比较 (1) 式和 (3) 式, 方程 (3) 右边增加了
简化 (3) 式, 令
(4) 式是一个典型的Ito方程, 它带有一个随机输入项, 并可具有随机初始条件, 其解过程为Markov过程。
B (t) 是一Wiener过程, dB (t) /dt是一正态白噪声, B (t) 的一维概率密度函数f (B) 为:
此式表明, B (t) 的均值E[B (t) ]=0, 方差D[B (t) ]=σ2t, 其中σ2为常数, 称为过程强度, 它取决于W (t) 的离散程度, 也即取决于入库洪水、出库泄量和库容自身的变异性。通过实际调查和资料分析, 可给出以上三个随机过程的标准差σQ (t) 、σq (t) 和σw (t) , 并考虑到三个随机输入过程的相互不相关性, 则可得输出过程的方差D[B (t) ]为:
由此可计算出确定均值的σ2作为 (5) 式的过程强度参数, σ2t反映了Q、q和W对B (t) 的综合作用。
[对于初始条件Z0, 可以是确定的, 也可以是随机的, 由起调时刻库水位的控制条件所决定的, 表明洪水过程开始时刻的影响。
《2.2随机微分方程的求解》
2.2随机微分方程的求解
对于典型的Ito方程 (4) 式, 由于非线性性, 很难直接求出其具有随机性质的解, 考虑到其解具有Markov性质, 可以利用随机过程的动力学方程来求解。就考察泄洪风险率而言, 关注的是解过程的概率密度, 即在调洪过程的各个时刻, 库水位的随机分布状况。
库水位解过程Z (t) 的一阶概率密度函数f (h, t) 应能满足随机过程的动力学方程, 并由此根据B (t) 的Wiener过程的性质, 导出实用的向前方程 (Fokker-Planck方程)
(7) 式为一确定性的偏微分方程, 适用于非线性的Ito方程问题。采用有限差分法求解。
为此, 对于时间导数, 采用后差分格式, 坐标导数采用中心差分格式, 由于方程中的φ (t, Z) 和g (t, Z) 与Z有关, 给求解方程带来困难, 假定φ (t, Z) 和g (t, Z) 在 (tn, tn+Δt) 时间内变化不大, 可采用tn时刻的φ (t, Z) 和g (t, Z) 值来替代, 则方程离散形式为:
其中, [A]为一三对角矩阵, {f}为密度函数矩阵, {F}为源项矩阵, 其元素为:
利用初始条件和边界条件, 可以迭代求解出f (t, z) 分布, 从而可以确定泄流量及库水位随机变化的均值过程。
《2.3调洪过程的泄洪风险计算方法》
2.3调洪过程的泄洪风险计算方法
对水库调洪过程中的泄洪风险极限状况, 应直接以库水位Z (t) 不超过坝顶高程Zt为标准, 即满足极限状态方程:
由此定义泄洪风险率Pf为:在一定的洪水重现期和泄洪建筑物设计规模以及各种可能的水库自然、工程和运行条件下, 发生洪水漫顶事件的概率:
这里Zt为坝顶高度, 而Z为一随机过程, 可由前述随机方程计算确定。由此可以Pf值来具体度量水库的泄洪风险率, 其物理概念明确、清晰, 更能反映问题的本质。
在求出库水位的概率密度函数f (h, t) 后, 水库漫顶的概率可由下式计算:
采用Pf值, 可以比较不同运行方案的泄洪风险。
《3 计算条件对泄洪风险灵敏度影响的分析》
3 计算条件对泄洪风险灵敏度影响的分析
某工程在设计中考虑了七种风险较高的洪水过程频率, 来分析设置非常泄洪设施方案的可行性及其风险率。
考虑随机过程输入对Z (t) 的影响, 选取几组变异系数 (δQ (t) , δq (t) , δw (z) ) 的组合, 分别选取三个不同时刻 (泄洪后150 h洪峰时刻、泄洪后350 h) 的值, 用 (6) 式求得B (t) 的σ值如表1所示。
Table 1 σ values for B (t)
《表1》
P/% | δQ | 0.1 | 0.1 | 0.1 | 0.1 |
δq | 0.01 | 0.01 | 0.01 | 0.01 | |
δw | 0.0 | 0.01 | 0.02 | 0.05 | |
0.01 | 62 847 | 137 616 | 252 738 | 615 168 | |
0.02 | 59 749 | 133 542 | 246 163 | 599 938 | |
0.05 | 55 938 | 128 072 | 237 049 | 578 540 | |
0.1 | 52 397 | 123 264 | 229 145 | 560 082 | |
0.2 | 49 117 | 118 080 | 220 232 | 538 894 | |
0.5 | 44 800 | 11 114 | 208 177 | 510 165 | |
1 | 41 354 | 104 943 | 197 239 | 483 870 |
计算条件的选择往往会影响泄洪风险率, 为定性地了解这些选择所存在的误差对泄洪风险率的影响程度, 这里选择P=0.01%的洪水过程为基准, 分析三个变异系数δQ (t) 、δq (t) 和δw (z) 的选择、起调初始水位随机性对泄洪风险灵敏度的影响。
《3.1变异系数的灵敏度分析》
3.1变异系数的灵敏度分析
以δQ (t) =10%、δq (t) =1%和δw (z) =5%为基准, 分别改变δQ (t) 、δq (t) 和δw (z) , 分析它们对泄洪风险的影响, 计算结果如表2所示。
由表2可见, δQ (t) 和δq (t) 的变化对于泄洪过程最大风险的影响成缩小趋势, δQ (t) 、δq (t) 的变化对泄洪库水位随机过程的均方差和泄洪风险的影响最大为0.79%, 而δw (z) 对泄洪风险影响十分明显, 其值越大, 风险增加越快。三个变异系数中, 以δw (z) 对泄洪风险影响最大, δQ (t) 和δq (t) 对泄洪风险影响不敏感。
《3.2 δw (z) 对风险影响分析》
3.2 δw (z) 对风险影响分析
为了更细致地分析δw (z) 对泄洪风险的影响, δQ (t) 和δq (t) 分别选为10%和1%, δw (z) 的选择考虑到目前地形测量的精度, 选择四种情况作为计算依据, 分别为0%、1%、2%和5%, 相应的B (t) 的σ值如表1所示, 计算不同δw (z) 值的水库泄洪风险。不同δw (z) 值的计算结果如表2、图1和图2所示。
分析可见, δw (z) 越大, 泄洪风险越大, 而且δw (z) 与最大泄洪风险的关系也表明这一点, 说明库容的随机性较大程度地影响泄洪过程的随机性大小。
由于目前对于入库洪水过程、出库泄洪过程以及水位库容的随机性研究尚不多见, 无法较为精确地确定三个变异系数, 为了相对比较不同泄洪运行方案的泄洪风险率, 这里选取δQ (t) 、δq (t) 和δw (z) 分别为10%、1%和5%。
Table 2 The results of sensitivity analysis
《表2》
序号 | δQ (t) | δq (t) | δw (z) | σB | σZ | Pf[Z≥Zt]/% | 误差/% |
1 | 0.10 | 0.01 | 0.05 | 615 168 | 5.74 | 32.73 | 0 |
2 | 0.05 | 0.01 | 0.05 | 612 770 | 5.72 | 32.68 | -0.15 |
3 | 0.20 | 0.01 | 0.05 | 624 668 | 5.83 | 32.99 | 0.79 |
4 | 0.10 | 0.05 | 0.05 | 615 855 | 5.74 | 32.76 | 0.09 |
5 | 0.10 | 0.10 | 0.05 | 617 997 | 5.76 | 32.81 | 0.24 |
6 | 0.10 | 0.01 | 0.00 | 62 847 | 0.59 | 1.16 | -96.46 |
7 | 0.10 | 0.01 | 0.01 | 137 616 | 1.28 | 5.25 | -83.96 |
8 | 0.10 | 0.01 | 0.02 | 252 738 | 2.36 | 14.50 | -55.70 |
9 | 0.10 | 0.01 | 0.10 | 1 225 496 | 11.39 | 41.13 | 25.67 |
注:误差栏为与δQ (t) =0.1, δq (t) =0.01和δw (z) =0.05情况比较的误差百分数。
Table 3 The analyses for the effect of random of initiative reservoir level on the discharge risk
《表3》
σz0/m | 0.0 | 0.5 | 1.0 | 2.0 | 3.0 | 4.0 | 5.0 |
μz/m | 602.46 | 602.46 | 602.46 | 602.46 | 602.46 | 602.46 | 602.46 |
σz/m | 5.74 | 5.74 | 5.81 | 6.07 | 6.47 | 6.99 | 7.61 |
Pfmax/% | 32.74 | 32.74 | 32.95 | 33.68 | 34.67 | 35.78 | 36.91 |
《3.3起调水位随机性对泄洪风险影响分析》
3.3起调水位随机性对泄洪风险影响分析
在前面的分析计算中, 没有考虑起调水位Z0的随机性, 即认为起调水位Z0是确定的, 然而实际上并不能保证起调水位Z0是一个确定值。如果考虑Z0的随机性, 则其对泄洪风险的影响有多大, 需作一分析。
为简化起见, 假定起调水位的随机性符合正态分布, 即
式中,
分析不同σz0值对库水位随机分布的均方差σz以及泄洪风险的影响, 如表3所示。在本算例中, 当σz0值小于1.0 m时, 起调水位随机性对于最大泄洪风险影响微弱, 可以忽略不计;当σz0>1.0 m时, 其影响会超过最大泄洪风险发生时间, 但增加的有限, 并会随时间增长, 这种影响会逐渐削弱。为简便起见, 在下面的计算中, 不考虑起调水位的随机性。
《4 采用非常泄洪设施方案对泄洪风险影响的分析》
4 采用非常泄洪设施方案对泄洪风险影响的分析
某水电站枢纽由大坝、泄水建筑物、地下厂房等组成, 大坝为混凝土双曲拱坝, 坝顶高程605 m, 最大坝高273 m。为了解决高水头、大流量的泄洪消能问题, 设计中采用分散泄洪、分区消能的原则, 泄洪设施由坝身7个表孔、8个深孔和5条泄洪洞组成, 共同宣泄洪水。该电站枢纽处于高山峡谷, 没有布置开敞式非常溢洪道的地形条件, 两岸各有一条16 m×22 m的漂木导流隧洞未予利用。设计上曾设想将这两条导流隧洞后期改造为永久泄洪隧洞, 但因水头高 (200 m左右) 、泄量大 (单洞泄量3 500 m3/s) , 国内外均无成熟的技术可资借鉴, 故未予采用。如将这两条导流洞改建为非常泄洪洞, 降低其设计标准, 减少正常泄洪设施规模 (原五条永久泄洪洞改由三条永久泄洪洞和两条非常泄洪洞组成) , 则可大大减少泄洪设施的工程造价, 并优化枢纽布置。为此利用上述风险分析方法, 论证该水电站采用非常泄洪设施方案的可行性, 并分析非常泄洪设施启用条件。
《4.1设计方案的泄洪风险分析》
4.1设计方案的泄洪风险分析
图3为不同入库洪水过程下水库泄洪风险计算结果。分析可见, 泄洪风险率与入库洪水过程的频率P有直接的关系, P值越小, 泄洪风险率越大, 当发生校核洪水 (P=0.01%) 时, 泄洪风险率最大为32.73%, 而且泄洪风险随泄洪过程而变, 在库水位达到最高水位时, 泄洪风险率最大, 当P≫0.5% (二百年一遇以下) 洪水过程时, 入库洪水、库容及泄洪流量的随机性基本上不会导致发生漫顶的事件, 风险率趋于0。
《4.2五条泄洪洞改为三条泄洪洞的泄洪风险分析》
4.2五条泄洪洞改为三条泄洪洞的泄洪风险分析
五条泄洪洞改为三条泄洪洞后, 泄洪风险分析计算结果如图4和表4所示。分析可见, 取消两条泄洪洞后, 由于泄洪能力下降, 各种频率洪水过程的泄洪风险都有不同程度的增加, 主要体现在低频率洪水过程, 尤其是P=0.01%和P=0.05%的洪水过程, 三条泄洪洞的泄洪风险超过了50%, 说明其库水位均值已超过坝顶。三条泄洪洞方案在P<0.10% (千年一遇以上) 的洪水过程中泄洪风险偏高, 存在着较高机率的洪水漫顶可能性, 因此, 五个泄洪洞改为三个泄洪洞方案不可取。
Table 4 The computational results of discharge risk under difference cases
《表4》
P/% | 0.01 | 0.02 | 0.05 | 0.10 | 0.20 | 0.50 | 1.00 | |
μz/m | 设计方案 | 602.47 | 600.73 | 598.16 | 595.96 | 593.40 | 589.29 | 584.96 |
三条洞方案 | 607.64 | 606.07 | 603.82 | 601.97 | 599.93 | 596.95 | 594.17 | |
584 m | 602.54 | 600.84 | 598.38 | 596.32 | 594.04 | 590.66 | 587.68 | |
586 m | 502.55 | 600.86 | 598.40 | 596.37 | 594.12 | 590.86 | 588.22 | |
588 m | 602.56 | 600.87 | 598.44 | 596.42 | 594.21 | 591.14 | 588.77 | |
590 m | 602.58 | 600.89 | 598.47 | 596.48 | 594.32 | 591.53 | 590.04 | |
592 m | 602.59 | 600.92 | 598.51 | 596.56 | 594.50 | 592.07 | 592.03 | |
594 m | 602.61 | 600.95 | 598.57 | 596.68 | 594.82 | 594.01 | 594.00 | |
δz/m | 设计方案 | 5.74 | 5.61 | 5.45 | 5.30 | 5.12 | 4.79 | 4.46 |
三条洞方案 | 5.85 | 5.66 | 5.43 | 5.25 | 5.09 | 4.88 | 4.67 | |
584 m | 5.87 | 5.70 | 5.47 | 5.30 | 5.13 | 4.92 | 4.63 | |
586 m | 5.87 | 5.69 | 5.46 | 5.29 | 5.12 | 4.88 | 4.60 | |
588 m | 5.87 | 5.68 | 5.46 | 5.29 | 5.11 | 4.85 | 4.57 | |
590 m | 5.87 | 5.68 | 5.45 | 5.28 | 5.10 | 4.83 | 4.51 | |
592 m | 5.87 | 5.68 | 5.45 | 5.27 | 5.09 | 4.80 | 4.28 | |
594 m | 5.87 | 5.68 | 5.45 | 5.26 | 5.07 | 4.74 | 4.52 | |
Pfmax/% | 设计方案 | 32.73 | 22.17 | 10.52 | 4.57 | 1.33 | 0.09 | 0.00 |
三条洞方案 | 67.34 | 57.36 | 41.17 | 28.03 | 15.90 | 5.13 | 1.22 | |
584 m | 33.59 | 23.12 | 11.33 | 5.25 | 1.83 | 0.26 | 0.03 | |
586 m | 33.64 | 23.18 | 11.38 | 5.31 | 1.87 | 0.28 | 0.04 | |
588 m | 33.71 | 23.26 | 11.48 | 5.39 | 1.94 | 0.31 | 0.05 | |
590 m | 33.79 | 23.37 | 11.59 | 5.50 | 2.02 | 0.37 | 0.07 | |
592 m | 33.88 | 23.49 | 11.73 | 5.64 | 2.16 | 0.46 | 0.23 | |
594 m | 34.00 | 23.63 | 11.91 | 5.87 | 2.44 | 0.82 | 0.98 |
《4.3非常泄洪洞方案的风险分析》
4.3非常泄洪洞方案的风险分析
正如以上分析, 三条泄洪洞方案在较低频率 (千年一遇以上) 的洪水过程情况下, 低水位时泄洪风险并没有明显的增加, 只是高水位时风险增加较大, 这为设置两条非常泄洪洞在高水位时启用来降低泄洪风险提供了机会。下面分析两条非常泄洪洞在不同水位启用时的泄洪风险。
表4为泄洪风险计算结果, 图5为各频率洪水情况下泄洪风险变化过程。由表4可见, 在一定库水位下启用两条非常泄洪洞后, 其泄洪风险基本上维持原设计方案的风险水平, 启用水位对于泄洪风险的影响也较小, 计算的几种启用水位下泄洪风险差异在1.5%以内, 表明非常泄洪洞的启用对于降低泄洪风险具有较大的作用。图5表明了泄洪风险变化过程, 可见, 不同启用水位对于泄洪风险的影响主要体现在启用水位以下的时间段内, 对于稀遇洪水最终的结果影响不大, 对于常遇洪水 (五百年一遇以下) , 启用水位过高时风险有所增加, 但因其最高水位较低, 整个泄洪风险也较低。
从上面的分析可见, 把两条导流洞改建为非常泄洪洞后, 只要在来水超过五百年一遇的洪水情况下才需启用, 而非常泄洪设施启用的初始水位对大坝泄洪设施防洪安全的影响并不十分明显, 并可以维持原设计方案的等同泄洪风险水平, 表明这一方案是可行的。
值得说明的是, 上述计算的泄洪风险率是计算洪水过程的频率下出现的概率, 即是条件概率, 例如P=0.01%洪水过程的泄洪风险率为33.59%, 这时发生漫顶的总概率为0.01%×33.59%=3.359×10-5, 是相当低的。
《图5》
Fig.5 The computational results of discharge risk under difference reservoir level using unusual tunnels
《5 结论》
5 结论
非常泄洪设施是优化水电站枢纽布置、降低工程成本的有效途径, 针对采用非常泄洪设施方案, 本文提出了一种基于随机微分方程的泄洪风险分析方法。在泄洪风险分析中, 引入随机微分方程, 可以全面地反映各种随机因素的影响, 是一种物理概念清楚、合理可靠的方法。
某工程非常泄洪设施方案泄洪风险分析结果表明, 非常泄洪洞方案通常只要在五百年一遇以上的稀遇洪水条件下起用。在不同库水位情况下启用, 均能达到基本维持原设计方案的泄流能力及泄洪风险水平, 库最高水位也基本与设计方案一致, 因此该工程在泄洪枢纽布置上采用非常泄洪设施方案是可行的, 基本上保持原设计方案的泄洪风险水平。