土坝边坡风险计算模型研究.docxVIP

土坝边坡风险计算模型研究 水库是水库的重要组成部分。建国以来我国共修建各类水库8万多座,其中土石坝占90%以上。这些水库在我国国民经济基础产业中具有重要地位,可在防洪、灌溉、发电、城乡供水等方面创造巨大的社会财富和经济效益。然而,大坝一旦失事将对国家和人民的生命财产带来不可估量的损失。据统计,在国内外溃坝失事事件中,土石坝所占的比例高达70.5%,滑坡是仅次于洪水漫顶的溃坝重要原因。目前,人们对土坝漫坝风险已经做了比较深入的研究,对土坝滑坡失稳风险研究则较少。因此,研究土坝滑坡失稳风险具有重要的现实意义。 1 风险定义 1.1 风险或后果 风险是某一事件产生我们所不希望的后果的可能性。风险是对危险的一种度量。具体来说,风险是指系统在规定的工作条件和规定的时间内,其不能完成预定功能的概率及由此所产生的后果。即风险可表示为失事概率及其后果的函数 ˉR=f(Ρ,C)(1)Rˉˉˉ=f(P,C)(1) 式中:ˉRRˉˉˉ——风险度; P——失事概率或风险概率; C——对应的失事后果或风险损失,常用相应于失事概率的损失期望值表示。 规定的工作条件是指系统(或结构)所处的环境即作用因素的状态和系统内部特性(抗力、荷载)。规定的时间则是指进行风险估算的基准期。 1.2 生命损失、经济损失和环境恶化 失事后果是指对系统失事进行划分后,相互独立的、某种形式和特定程度的失事影响。一般将失事后果分为三个方面,即:生命损失、经济损失和环境恶化。 C=C(s1,s2,s3)=3Σi=1si(2)C=C(s1,s2,s3)=Σi=13si(2) 式中:s1、s2和s3分别表示生命损失、经济损失和环境恶化的损失或治理费用的期望值。 由于风险分析的重点是计算风险概率,因此,本文主要研究土坝滑坡风险概率的计算。 1.3 坝坡失稳抗力 土坝滑坡失稳风险可定义为:当洪水位超过某一极限时,作用于土坝的荷载超过其抗力,导致坝坡失稳的概率。 ˉR=Ρ(R＜S)(3)Rˉˉˉ=P(R＜S)(3) 这里荷载R可以是土坝边坡的抗滑力矩、坝顶高程或设计水位,而S则可以是土坝边坡的滑动力矩,也可以是土坝上游水位。 2 设计不出土错 影响土坝滑坡失稳的因素很多,如水位、土的物理力学性质指标、结构尺寸和施工质量等。如果土坝在运用过程中,这些因素的值和设计值相吻合,那么,只要设计不出差错,土坝在边坡滑动方面就不会产生风险。但实际情况上述因素存在不确定性,由此产生了土坝滑坡失稳的风险。 2.1 水位的不确定性 水对土坝的作用包括静水压力、动水压力、浪压力以及渗透力等。这些作用均与水位有关,包括上游库水位和下游坝址水位以及上下游的水位差。水位是一个随机变量,具有较大的不确定性。上游水位与天然来水量、上一时段末的库水位及下泄流量有关,而天然来水量和下泄流量都是不确定的随机变量,具有不确定性。但每一个水位的发生均服从一定的概率分布,即有一定的发生概率。下游水位一般与下泄流量有关。下泄流量也有一定的不确定性。所以,下游水位也是不确定的。据文献研究成果,库水位的不确定性可以用库水位概率曲线来表示,而库水位概率曲线可由入流频率曲线和水库状态转移方程进行概率演算得到。 2.2 土坝边坡风险计算中重力密度的确定 影响土坝滑坡风险的主要是以下几种土质的物理力学指标,即凝聚力c、内磨擦角φ和重力密度γ。文献的研究结果表明,在影响土坝边坡风险计算的土质物理力学指标中,c、φ变异性影响程度较大,而重力密度γ的变异性则影响很少。因此,本文只考虑c和φ的变异影响,而把重力密度当作定值处理。上述随机变量的概率分布形式根据经验:凝聚力c的变异性较大,采用极值Ⅰ型分布,摩擦系数tgφ采用对数正态分布。 3 水库滑动风险模型及其求解 3.1 土坝坝坡滑动力的概率密度函数法 根据土力学和水工建筑物的有关理论,当坝上游水位达到一定高度时,作用于坝体的滑动力矩S将大于其抗滑力矩R,从而导致土坝滑坡失稳。据此,土坝坝坡失稳风险模型为: 设功能函数 G(R,S)=R-S=G(.) (4) 则ˉR=Ρ(S＞R)=∫∫G(.)＜0fR,S(r,s)drds(5)Rˉˉˉ=P(S＞R)=∫∫G(.)＜0fR,S(r,s)drds(5) 式中:fR,S(r,s)是土坝坝坡滑动力矩S和抗滑力矩R的联合概率密度函数。 显然,直接用(5)式计算风险概率是极为困难的。考虑到S主要由上游水位的不确定性所决定,而R主要由大坝土质的物理力学特性所决定,所以可以认为两者是相互独立的随机变量。于是(5)式变为: 利用概率组合的方法可以间接估算ˉRRˉˉˉ。设滑动力矩S与坝上游水位H的联合概率密度函数为 fH,S(s,h)=f(s|h)f0(h) (7) 式中:f(s|h)——给定某一水位h情况下滑动力矩s的条件概率密度函数; f0(h)