基于Monte-Carlo方法的土石围堰挡水导流风险分析.PDFVIP

基于Monte-Carlo方法的土石围堰挡水导流风险分析 1 1 1 1 1 胡志根 ，刘 全 ，贺昌海 ，肖焕雄 ，周宜红 ， 2 2 2 傅 峥 ，李定葵 ，郑家祥 (1 武汉大学水利水电学院,湖北武汉 430072; 2 国家电力公司成都勘测设计研究院,四川成都 610072) 摘要：运用 Monte-Carlo 方法模拟施工洪水入库过程和导流建筑物泄流的随机性,通过仿真 分析确定上游围堰堰前水位变化过程及其分布函数。在此基础上,根据导流建筑物的设计规 模确定导流系统的风险。通过实例验证分析说明,研究方法和计算模型是可靠的、适用的。 关键词：土石围堰；施工导流；风险分析； Monte-Carlo 方法 施工导流中的不确定性因素很多,研究施工导流风险率的方法与途径也多种多样。但由 [1～4] 于风险分析方法引入到水电施工界的时间较短,成果还有一定的局限性 。 Borgman and Shanne 、Lynn 等将 Poisson 过程引入风险分析模型,取部分历时洪峰序列, 假定在(0,ｔ)内洪峰个数服从齐次 Poisson 过程,洪峰大小服从指数分布,研究了最大超定量 洪峰分布问题。之后,在这方面还有不少学者都做了工作,给出了一些风险率模型。1989 年 以后,邓永录、徐宗学、叶守泽、肖焕雄、韩采燕等[5～9]以随机点过程理论为依据,提出了一 些风险率计算模型,但由于只考虑了水文不确定性,而未考虑系统泄流能力的水力不确定性, 很难直接应用于施工导流工程的实际。姜树海[10,11]基于调洪演算和水库坝前水位变化的随机 微分方程,以 Ito 随机微分方程的理论为基础,建立了与水库坝前水位变化直接联系的风险率 模型,但水库坝前水位的分布推求和失效概率的计算困难。应用调洪演算的随机微分方程方 法求解施工导流调洪演算和堰前水位变化时,要求推求其微分方程组。对于不同工程项目施 工洪水的特性,微分方程的求解可能存在困难。为此,可利用 Monte-Carlo 方法模拟施工洪水 过程和导流建筑物泄流工况,通过系统仿真进行施工洪水调洪演算,用统计分析模型确定上 游围堰堰前水位分布和在不同的导流标准条件下围堰运行的动态风险。 1. 围堰挡水时导流系统风险因素 1 1 施工洪水入库的随机性 由于入库施工洪水过程的随机性,水文资料的收集、整理和设计洪水过程线推求结果与 实际洪水过程间的偏差,使施工设计洪水具有随机性。根据坝址的实测水文资料,按放大典型 [13] 洪水过程线方法确定计算洪水过程线,最大洪峰流量为Ｐ-Ⅲ随机抽样均值,其密度函数为 βα α−1 −β(x −a0 ) f (x ) (X −a ) e （1） 0 Γ(α) 式中： α,β,a0 分别为 P-Ⅲ型分布的形状、刻度和位置参数；Γ(α) 为α的伽玛参数。 2 2C 且 α 4 / Cs ，β 2 ，α μ (1− V ) ；C 为 P-Ⅲ型分布的离差系数；C μ C C 0 Q C s