有限数据条件下边坡可靠度分析的Bootstrap方法

基于Bootstrap抽样技术提出了有限数据条件下边坡可靠度分析方法。简要介绍了传统的边坡可靠度分析方法。采用Bootstrap方法模拟了抗剪强度参数概率分布函数的统计不确定性。以无限边坡为例研究了抗剪强度分布参数和分布类型不确定性对边坡可靠度的影响规律。结果表明：基于有限数据估计的样本均值、样本标准差和AIC值具有较大的变异性,这种变异性进一步导致了抗剪强度参数概率分布函数存在明显的统计不确定性。在考虑抗剪强度参数概率分布函数的统计不确定性时,边坡可靠度指标应为具有一定置信度水平的置信区间,而不是传统可靠度分析中的固定值。边坡可靠度指标的置信区间变化范围随安全系数的增加而增大,同时考虑分布参数和分布类型不确定性计算的可靠度指标具有更大的变异性和更宽的置信区间变化范围。Bootstrap方法为有限数据条件下抗剪强度参数概率分布函数统计不确定性的模拟以及边坡可靠度的评估提供了一条有效的途径。     

基于BUS方法的无限长边坡可靠度更新

某一特定岩土场地的试验数据、监测资料和观测信息等通常十分有限,然而贝叶斯方法却可充分利用有限的场地信息克服试验数据样本量较小的不足。为有效估计有限样本条件下参数统计特征,提出了基于结构可靠度方法和贝叶斯更新(BUS)的边坡可靠度更新方法,通过融入直剪试验数据更新无限长边坡可靠度验证了提出方法的有效性,并系统探讨了岩土体参数先验信息如试验样本量、概率分布和似然函数模型对边坡可靠度更新的影响规律。结果表明:BUS方法能够考虑岩土体参数概率分布和似然函数模型的影响,融入有限的场地信息准确地估计参数统计特征和更新边坡可靠度,为解决有限样本条件下边坡可靠度更新问题提供了一条有效的途径。土体参数概率分布对边坡可靠度更新结果(参数后验均值、标准差以及更新的失效概率)具有重要的影响,基于常用的正态和对数正态分布的边坡可靠度更新结果偏于保守,相比之下,似然函数模型对边坡可靠度更新结果的影响相对较小。此外,岩土体参数不确定性和更新的边坡失效概率均随着试验样本量的增大而减小,但当样本量增大到一定程度时它们的变化不大。     

基于多源试验数据空间变异土体参数概率反演及边坡可靠度更新

岩土工程现场勘察试验通常只能获得有限的试验数据,据此难以真实地量化土体参数的空间变异性。提出了考虑土体参数空间变异性的概率反演和边坡可靠度更新方法,基于室内和现场两种不同来源的试验数据概率反演空间变异参数统计特征和更新边坡可靠度水平,并给出了计算流程。此外为合理地描述土体参数先验信息,发展了不排水抗剪强度非平稳随机场模型。最后通过不排水饱和黏土边坡算例验证了提出方法的有效性,并探讨了试验数据和钻孔位置对边坡后验失效概率的影响。结果表明：提出方法实现了空间变异土体参数概率反演与边坡可靠度更新的一体化,基于有限的多源试验数据概率反演得到的土体参数均值与试验数据非常吻合,明显降低了对参数不确定性的估计,更新的边坡可靠度水平显著增加。受土体参数空间自相关性的影响,试验数据对钻孔取样点附近区域土体参数统计特征更新的影响明显大于距离取样点较远区域。     

考虑参数空间变异结构的结构化交叉约束随机场模拟方法研究

斜坡岩土体抗剪强度参数的空间变异性具有一定的结构性。为研究岩土体参数空间变异结构对边坡失效概率的影响,依据变异函数的内涵推导出变程与相关距离的数学变换关系,并在此基础上提出了结构化交叉约束随机场模拟方法,用以模拟具有互相关性的参数随机场。建立了结构化交叉约束随机场计算模型,研究不同空间变异结构的抗剪强度参数对边坡失效概率的影响。研究结果表明:结构化交叉约束随机场可用于生成模拟具有复杂各向异性空间变异结构的参数随机场,由于考虑了随机偏差、条件数据和空间变异结构,能较为真实地反映地层实际参数,数据波动较条件参数插值场小。可靠性分析结果表明:不考虑抗剪强度参数空间结构分析易高估边坡的失效概率;考虑c′和φ′互相关性时,失效概率随着相关系数的增加而增加,当参数间呈负相关性时更容易高估边坡的失效概率。     

基于Copula函数的岩土结构物系统可靠度分析

研究了基于Copula函数的抗剪强度参数联合分布模型对岩土结构物系统可靠度的影响规律,揭示了抗剪强度参数间相关结构对岩土结构物系统可靠度影响的过程。介绍了不完备概率信息条件下基于Copula函数的抗剪强度参数联合分布模型构造方法,提出了Copula理论框架下岩土结构物系统失效概率计算的蒙特卡罗模拟方法,以典型岩土结构物如挡土墙和岩质边坡为例研究了Copula函数对系统可靠度的影响。结果表明:Copula函数为不完备概率信息条件下抗剪强度参数联合分布模型的构建提供了一种有效的工具。不完备概率信息条件下岩土结构物系统可靠度不唯一,不同Copula函数计算的系统失效概率差别非常大,这种差别在实际工程可靠度设计中应该引起重视。Copula函数对系统可靠度的影响分为两个机制不同的阶段,首先是参数间相关结构对系统中单一失效模式可靠度的影响,其次是每一个失效模式对系统可靠度的影响。     

水利水电工程岩基抗剪强度参数二维分布模型构造的Copula方法

基于国内103个水利水电工程1 174组岩基抗剪强度试验数据,采用Copula函数研究岩基抗剪强度参数联合分布模型,探讨水利水电工程中岩基抗剪强度参数联合分布模型构建方法。利用最小二乘法求出岩基抗剪强度参数试验数据的相关统计参数,基于AIC准则识别出岩基抗剪强度参数边缘分布。选择4种Copula函数构造岩基抗剪强度参数二维分布模型,探讨了基于Copula函数的岩基抗剪强度参数二维分布模型的优越性。结果表明：水利水电工程岩基抗剪强度参数存在明显的统计负相关性。Copula方法能够构造具有任意边缘分布和任意相关结构的岩基抗剪强度参数联合分布模型,它为构造抗剪强度参数联合分布模型提供了一种简便的工具。已知岩基抗剪强度参数的边缘分布函数和相关系数不能唯一确定岩基抗剪强度参数的联合概率分布模型,在抗剪强度参数边缘分布函数和相关系数完全相同的前提下,不同Copula函数建立的抗剪强度参数联合概率分布模型差异显著。与常用的抗剪强度参数二维正态分布模型相比,基于Copula函数的抗剪强度参数二维分布模型具有较强的灵活性,它能更好地拟合原始观测数据。水利水电工程中惯用小值平均法确定标准值,当摩擦系数取较小值时,不同Copula函数构造的黏聚力的条件累积分布函数差异显著,这将对抗剪强度参数标准值的选取以及相应的设计方案具有明显的影响。     

基于贝叶斯理论的抗剪强度参数最优Copula函数识别

提出基于贝叶斯理论的抗剪强度参数最优Copula函数识别方法,首先简要介绍了基于Copula函数的岩土体抗剪强度参数相关结构表征方法,给出常用的识别最优Copula函数的最小平方欧氏距离法和AIC(akaike information criterion)准则。其次,采用蒙特卡洛模拟方法验证了贝叶斯理论识别最优Copula函数的有效性,比较了3种方法的最优Copula函数识别能力,并分析了影响贝叶斯理论识别精度的主要因素。最后,收集了实际工程共23组抗剪强度参数试验数据,研究了贝叶斯理论在抗剪强度参数最优Copula函数识别中的应用。结果表明,贝叶斯理论能够有效地识别表征抗剪强度参数间相关结构的最优Copula函数,且能有效考虑先验信息对识别结果的影响;与传统的最小平方欧氏距离法和AIC准则相比,贝叶斯理论的识别能力和识别精度都更高;抗剪强度参数的样本数目、相关性大小、真实Copula函数类型以及先验信息都对贝叶斯理论的识别精度具有重要的影响。此外,常用的Gaussian Copula函数并不总是表征抗剪强度参数间相关结构的最优Copula函数。     

先验概率分布及似然函数模型的选择对边坡可靠度评价影响的定量评估

受工程勘察成本及试验场地限制,可获得的试验数据通常有限,基于有限的试验数据难以准确估计岩土参数统计特征和边坡可靠度。贝叶斯方法可以融合有限的场地信息降低对岩土参数不确定性的估计进而提高边坡可靠度水平。但是,目前的贝叶斯更新研究大多假定参数先验概率分布为正态、对数正态和均匀分布,似然函数为多维正态分布,这种做法的合理性有待进一步验证。总结了岩土工程贝叶斯分析常用的参数先验概率分布及似然函数模型,以一个不排水黏土边坡为例,采用自适应贝叶斯更新方法系统探讨了参数先验概率分布和似然函数对空间变异边坡参数后验概率分布推断及可靠度更新的影响。计算结果表明：参数先验概率分布对空间变异边坡参数后验概率分布推断及可靠度更新均有一定的影响,选用对数正态和极值I型分布作为先验概率分布推断的参数后验概率分布离散性较小。选用Beta分布和极值I型分布获得的边坡可靠度计算结果分别偏于保守和危险,选用对数正态分布获得的边坡可靠度计算结果居中。相比之下,似然函数的影响更加显著。与其他类型似然函数相比,由多维联合正态分布构建的似然函数可在降低对岩土参数不确定性估计的同时,获得与场地信息更为吻合的计算结果。另外,构建似然函数时不同位置处测量误差之间的自相关性对边坡后验失效概率也具有一定的影响。     

双滑块边坡锚固系统时变可靠性分析

考虑锚杆锚固段从岩体中拔出、拉杆拉断、拉杆从注浆体中拔出等失效模式,利用系统可靠性原理和极限平衡分析方法,并基于Monte-Carlo抽样原理提出双滑块岩质边坡锚固系统破坏概率的直接求解方法。同时考虑锚杆钢筋的腐蚀与软弱滑动面抗剪强度c、? 的时变性,建立了考虑锚杆多失效模式双滑块岩质边坡锚固系统的时变可靠性模型。算例计算结果表明：软弱滑动面上的强度参数c、? 的时变性和注浆体与围岩之间的抗力时变性对锚固系统的破坏概率的影响较大,而锚杆的腐蚀对锚固系统破坏概率的影响不明显。     

岩石抗剪强度参数的理论概率分布形态研究

岩石抗剪强度参数的概率分布形态是岩石工程可靠度分析和设计的基础。在考虑完整岩石压缩强度为服从正态分布随机变量的条件下,针对Mohr-Coulomb和Drucker-Prager屈服准则的线性系数形式,基于随机变量函数分布理论推导出岩石抗剪强度参数内摩擦角 和黏聚力 的概率密度函数。 和 概率密度函数显示：不仅岩石压缩强度和抗剪强度参数概率分布具有非一致性,而且根据不同屈服准则计算得出的岩石抗剪强度参数概率分布也具有非一致性。在进一步分析屈服准则系数具有不同变异系数和相关系数时的抗剪强度参数概率密度函数特征的基础上,提出根据概率分布的偏度和峰度确定一般情况下抗剪强度参数概率分布形态的方法,从理论上解决岩石压缩强度与抗剪强度参数分布的协调性问题。最后,对大理岩常规三轴压缩试验得出的抗剪强度参数进行大样本统计分析,验证了其概率密度函数理论推导的正确性以及概率分布形态确定方法的合理性。该研究为实际岩石强度概率分析时选择抗剪强度参数合理概率分布形式提供了理论指导。     

Modeling uncertainty in stability analysis for design of embankment dams on difficult foundations

A probabilistic 3-D slope stability analysis model (PTDSSAM) is developed to evaluate the stability of embankment dams and their foundations under conditions of staged construction taking into consideration uncertainty, spatial variabilities and correlations of shear strength parameters, as well as the uncertainties in pore water pressure. The model has the following capabilities: (1) conducting undrained shear strength analysis (USA) and effective stress analysis (ESA) slope stability analysis of staged construction, (2) incorporation of field monitored data of pore water pressure, and (3) incorporation of increase of undrained shear strength with depth, effective stress, and pore water pressure dissipation. The PTDSSAM model is incorporated in a computer program that can analyze slopes located in multilayered deposits, considering the total slope width.

The main outputs of the program are the geometric parameters of the most critical sliding surface (i.e., center of rotation/radius of rotation and critical width of failure), mean 2-D safety factor, mean 3-D safety factor, squared coefficient of variation of resisting moment, and the probability of slope failure. The program is applied to a case study, Karameh dam in Jordan. Monitored data of induced pore water pressure in the dam embankment and soft foundation were gathered during dam construction.

The stability of Karameh dam embankment and foundation was evaluated during staged construction using deterministic and probabilistic analysis. Foundation stability was evaluated based on the monitored data of pore water pressure.

The study showed that the mean values of the corrective factors which account for the discrepancies between the in situ and laboratory-measured values of soil properties and for the modeling errors have significant influence on the 2-D safety factor, 3-D safety factor, slope probability of failure, and on the expected failure width.

The degree of spatial correlation associated with shear strength parameters within a soil deposit also influences the probability of slope failure and the expected failure width. This correlation is quantified by scale of fluctuation. It is found that a larger scale of fluctuation gives an increase in the probability of slope failure and a reduction in the critical failure width.     

Impact of copula selection on geotechnical reliability under incomplete probability information

This paper aims to investigate the impact of copula selection on geotechnical reliability under incomplete probability information. The copula theory is introduced briefly. Thereafter, four copulas, namely Gaussian, Plackett, Frank, and No. 16 copulas, are selected to model the dependence structure between cohesion and friction angle. A copula-based approach is used to construct the joint probability density function of cohesion and friction angle with given marginal distributions and correlation coefficient. The reliability of an infinite slope and a retaining wall is presented to demonstrate the impact of copula selection on reliability. The results indicate that the probabilities of failure of geotechnical structures with given marginal distributions and correlation coefficient of shear strength parameters cannot be determined uniquely. The resulting probabilities of failure associated with different copulas can differ considerably. Such a difference increases with decreasing probability of failure. Significant difference in probabilities of failure could be observed for relatively small coefficients of variation of the shear strength parameters or a strong negative correlation between cohesion and friction angle. The Gaussian copula, often adopted out of expedience without proper validation, may not capture the dependence structure between cohesion and friction angle properly. Furthermore, the Gaussian copula may greatly underestimate the probability of failure for geotechnical structures.     

基于极限应变判据-动态局部强度折减的边坡破坏演化数值模拟

岩土材料具有极限应变特征,其值可以通过室内试验或数值模拟方法求得。边坡岩土体主要为剪切破坏,可采用极限剪应变作为材料破坏的判据。事实上,边坡岩土体强度参数劣化并非整体性的,而是一个渐进性局部损伤至整体失稳的过程。本文基于极限应变判据,建立了边坡破坏的动态局部强度折减方法,此方法对边坡中超过极限剪应变值的单元进行强度折减。通过不断折减计算过程中产生的超过极限剪应变值单元,直到坡内超过极限剪应变值的单元贯通,认为边坡整体发生破坏。将该方法应用到实际边坡工程中,计算结果和边坡破坏模式及变形监测数据基本吻合。此方法在分析边坡渐进性破坏及稳定性评价方面具有较好应用前景。     

Probabilistic slope stability analysis by a copula-based sampling method

In probabilistic slope stability analysis, the influence of cross correlation of the soil strength parameters, cohesion and internal friction angle, on the reliability index has not been investigated fully. In this paper, an expedient technique is presented for probabilistic slope stability analysis that involves sampling a series of combinations of soil strength parameters through a copula as input to an existing conventional deterministic slope stability program. The approach organises the individual marginal probability density distributions of componential shear strength as a bivariate joint distribution by the copula function to characterise the dependence between shear strengths. The technique can be used to generate an arbitrarily large sample of soil strength parameters. Examples are provided to illustrate the use of the copula-based sampling method to estimate the reliability index of given slopes, and the computed results are compared with the first-order reliability method, considering the correlated random variables. A sensitivity study was conducted to assess the influence of correlational measurements on the reliability index. The approach is simple and can be applied in practice with little effort beyond what is necessary in a conventional analysis.     

Influence of the spatial variability of shear strength parameters on rainfall induced landslides: a case study of sandstone slope in Japan

Rainfall-induced landslides frequently occur in humid temperate regions worldwide. Research activity in understanding the mechanism of rainfall-induced landslides has recently focused on the probability of slope failure involving non-homogeneous soil profiles. This paper presents probabilistic analyses to assess the stability of unsaturated soil slope under rainfall. The influence of the spatial variability of shear strength parameters on the probability of rainfall-induced slope failure is conducted by means of a series of seepage and stability analyses of an infinite slope based on random fields. A case study of shallow failure located on sandstone slopes in Japan is used to verify the analysis framework. The results confirm that a probabilistic analysis can be efficiently used to qualify various locations of failure surface caused by spatial variability of soil shear strength for a shallow infinite slope failure due to rainfall.     

黄河小浪底水库风雨沟西侧边坡稳定性研究

运用结构分析法、关键块体分析法和工程地质类比法定性评价边坡稳定性,抽取边坡岩体破坏模式。根据自制简易纵剖面仪现场绘制了结构面表面轮廓曲线,运用修正直边法和JRC—JCS模型定向统计估算结构面抗剪强度。通过对比验证定向统计估算结果与试验结果,确定计算参数取值标准。在以上两方面研究的基础上,建立边坡极限平衡计算模型,进行边坡破坏概率计算。最后,对边坡稳定性进行了综合评价。