干支流洪水作用下的水利工程洪水漫顶风险率计算研究.docVIP

干支流洪水作用下的水利工程洪水漫顶风险率计算研究 摘 要：洪水漫顶是影响水利工程安全的一个重要因素。干支流交汇处支流上水利工程不仅受到支流洪水的直接影响还会受到干流洪水顶托的影响，采用 Copula-Monte Carlo 模拟方法构建干支流洪水作用下大坝的洪水漫顶率计算模型， 并分析不同 Copula 函数、不同干支流洪水相关性对漫顶风险率的影响。以桂平航运枢纽为例，研究结果表明：最优的Clayton Copula 计算的洪水漫顶率为 1.3810-4;由于干支流洪水间相关性较弱，不同 Copula 函数对洪水漫顶的影响并不明显，GH、AMH 、CClayton、Frank Copula 都会一定程度地低估洪水漫顶风险;随着干支流洪水流量相关性的增强，洪水漫顶率也会相应的增加。 关键词：洪水漫顶;Copula-Monte Carlo 模拟;干支流洪水;桂平航运枢纽 水库的防洪标准[1]是与大坝的防洪安全密切相关的。大坝防洪安全最值得关注的是洪水漫顶[2]，洪水漫顶是影响水利工程安全的一个重要因素。因此，洪水漫顶风险率的计算分析对于大坝的防洪安全至关重要。大坝防洪风险分析的方法主要有直接积分方法、一次二阶矩方法和蒙特卡罗模拟方法。直接积分法[3] 可求得精确解，但由于涉及多重积分的问题，只适用于低维且调洪函数为简单表达式的情况。一次二阶矩方法要求计算函数对变量的偏导数，只能给出近似结果，如莫崇勋、董增川等[4]采用一次二阶矩法计算澄碧河水库漫坝风险。蒙特卡罗模拟方法是通过对实测资料统计分析，确定合适的随机模型，然后应用随机模型，生成足够多的、各种可能组合形式的洪水过程，供水库调洪计算所用，如梅亚东和谈广鸣[5]采用蒙特卡罗模拟法计算大坝防洪安全综合风险率;吕弯弯等[6]运用蒙特卡罗法计算土石坝洪水漫顶风险率的模型，分析洪水、风浪、水库运行调度方式和水库特性对洪水漫顶风险率的影响。 传统的大坝洪水漫顶风险率计算方法通常只对洪峰流量进行单变量频率分析经调洪演算计算漫坝风险率。后来，越来越多学者认为洪水是一个多变量的随机事件，研究了同时考虑洪峰洪量的两变量洪水频率分析的大坝洪水漫顶风险率[7]，表明两变量洪水频率分析更合理，以上都是仅基于单一洪水作用下的洪水漫坝风险率研究。然而随着我国经济的快速发展，越来越多沿江城市在穿越河流的支流入汇处修建了大量的水利工程，这些工程同时起着抵御外江洪水倒灌和保证城市排涝及生态需水等多重作用。对于这些干支流洪水作用下的水利工程而言，漫顶风险率的计算需要考虑干支流洪水的共同作用。目前国内外研究较多的是入库洪水、泄流能力、水库调度规则等对漫顶风险率计算的影响，还未见对于同时考虑干支流洪水对漫顶风险影响的研究。为此，本文通过引入 Copula 函数，开展基于 Copula 函数的干支流洪水作用下的水利工程洪水漫顶风险率的研究，分析不同 Copula 函数以及不同干支流洪水相关性对漫顶风险率计算的影响。 1 干支流洪水作用下的水利工程洪水漫顶风险率计算模型 1.1 大坝洪水漫顶风险分析模型 洪水漫顶风险率是坝前最高水位 Z 超过坝顶高程 Zd 发生的概率，可表示为Pf = P(Z Zd ) (1)一般而言，坝顶高程 Zd 随机变动幅度变动不大，可将它视为确定值。坝前最高水位 Z 则存在着诸多不确定的因素，如入库洪水过程的不确定性，泄流能力的不确定性等，因此坝前最高水位 Z 属于随机变量，漫顶率的计算可转化为坝前最高水位 Z 的分布研究，一旦累积分布 F(Z) 确定， 漫顶率就可以通过Pf = P(Z Zd )=1- F(Z) 计算求得。由于干支流洪水作用下水利工程的设计洪水位不仅会受支流洪水的影响还会受到干流洪水顶托影响，利用坝址洪水进行调洪演算反推坝前水位 Z 时需考虑干支流洪水流量间的相关性。因此，本文将采用Copula-Monte Carlo 模拟法求解干支流洪水作用下水利工程的洪水漫顶风险率。 1.2 基于 Copula 的干支流洪水联合分布 Copula 函数可用来描述多变量水文要素间的相关性，可构造出边缘分布服从任意分布的联合分布。假设 X 干和 Y 支分别表示干流和支流的洪水流量，则二维 Copula 函数为= F(y 支 )分别为边缘分布函数。 本文选取 5 种常用的阿基米德 Copula 函数，即 GH、Frank 、AMH、Clayton 和 CClayton Copula 函数来描述干支流洪水间相关结构，其中，GH 和 CClayton Copula 均具有上尾相关性，相关系数分别为 2-2-1/ 和2-1/，Clayton Copula 具有下尾相关性，Frank 和 AMH Copula 函