基于Copula函数的边坡可靠度高效分析方法

方旭; 冯晓波; Xu FANG; Xiao-Bo FENG

中国农村水利水电 ›› 2021 ›› (4) : 201-206,211.

水电建设

基于Copula函数的边坡可靠度高效分析方法

作者信息 +

Efficient Ways to Analyze Slope Reliability Based on Copula Function

Author information +

稿件信息 +

摘要

将Copula函数与高效蒙特卡洛方法结合，提出了含相关随机变量的边坡可靠度高效分析方法。以两个岩质边坡稳定性问题为例验证了所提方法的有效性。结果表明，该方法相比于直接蒙特卡洛方法在保证边坡失效概率的准确性的同时计算效率更高。Copula函数可以构造出具有不同相关结构的岩土体参数的联合概率分布，与高效蒙特卡洛方法（即蒙特卡洛重要抽样方法和子集模拟）结合能高效地处理含多种相关随机变量的边坡可靠度计算问题，相比于现行的Nataf变换方法结果更能体现岩土边坡真实稳定性。此外，该方法也能高效地计算含有复杂的隐式功能函数的边坡可靠度，研究成果拓展了高效蒙特卡洛方法在边坡可靠度分析中的应用。

Abstract

In this paper， the copula function is combined with efficient Monte Carlo Method， methods which can efficiently analyze the slope reliability problem involving correlated random variables are proposed. Two rock slope stabilities are studied to verify the effectiveness of the proposed methods as examples. The results show that these methods are more efficient than the direct Monte Carlo Method and the accuracy of the results is guaranteed. Based on the copula function which can be used to construct the joint distribution of different related structural parameters of rock and soil mass， efficient Monte Carlo Method （Monte Carlo important sampling method and subset simulation） can efficiently deal with the slope reliability problem with correlated random variables. Comparing with the existing Nataf transformation method， the results can reflect real geotechnical slope stability more. In addition， the methods can efficiently analyze the reliability of slope with complex implicit functions. The results expand the application of Monte Carlo important sampling method and subset simulation in slope reliability analysis.

关键词

Copula函数 / 边坡可靠度 / 相关随机变量 / 高效蒙特卡洛方法

Key words

copula function / slope stability / correlation random variables / efficient Monte Carlo Method

基金

国家自然科学基金面上项目(51879204)

引用本文

EndNote

Ris (Procite)

Bibtex

导出引用

方旭 , 冯晓波. 基于Copula函数的边坡可靠度高效分析方法[J].中国农村水利水电, 2021(4): 201-206,211

Xu FANG , Xiao-Bo FENG. Efficient Ways to Analyze Slope Reliability Based on Copula Function[J].China Rural Water and Hydropower, 2021(4): 201-206,211

0 引言

近年来，基于可靠度理论的分析和设计方法在岩土边坡工程的安全性评价中越来越受到重视。在英美等一些国家，岩土设计规范正渐渐向可靠度设计转变^［1，2］。陈祖煜等^［3］认为，我国岩土工程分析与设计也应往可靠度分析方向发展。作为一种常见的岩土结构物，边坡的稳定性分析是岩土工程中一个重要的工作。边坡岩土性质多变，其各部分岩土参数（比如岩土抗剪强度参数）存在不确定性。传统的确定性分析将土体参数视为唯一确定值^［4］，以有限元、极限平衡法等方法计算得到单一的边坡安全系数为依据进行边坡稳定性的评估，忽略了测量误差以及土体不均匀性等不确定性因素，导致分析结果不合理。因此，应对岩土边坡进行边坡稳定性可靠度分析^［5，6］。

对边坡稳定进行可靠度分析，往往需要事先获得岩土体物理力学参数的概率分布函数。同时，已有研究表明，岩土体物理力学参数间具有明显的相关关系，且有着不同的相关结构。如岩土体抗剪强度参数间具有很强的统计负相关性^［7-9］，土体剪切应力参数间也有负相关性。为了计算简便，很多情况下涉及相关变量的可靠度计算只依靠参数的边缘分布进行可靠度分析，忽略了参数间的相关性。文献［7］-［9］研究发现岩土体抗剪强度参数间负相关性对边坡稳定性有很大影响，忽略这种负相关性将会低估边坡可靠度，导致边坡设计方案不经济。也有研究考虑了相关性，但是将相关结构默认为高斯相关结构，运用Nataf变换解决相关非正态变量的独立化问题，与实际上参数间的相关结构可能并不符合，这样取得的分析结果不准确。

近年来，表征参数间相关性的Copula函数在岩土工程可靠度研究领域被广泛用于构建含多种相关结构的随机变量间的联合概率分布函数。如唐小松^［10］等学者采用Copula函数表征了抗剪强度参数间的相关性以及构建了基桩荷载位移双曲线参数间的联合概率模型。在这一情形下，需要发展基于Copula函数构建了参数间的联合概率分布模型的边坡可靠度分析方法。目前的研究中，极少有学者研究含多种相关结构随机变量的边坡可靠度分析方法。仅唐小松^［10］等在应用Copula函数表征参数间相关性后，运用直接积分法或者直接蒙特卡洛方法来计算边坡可靠度。直接积分法在遇到高维且具有复杂表达式的可靠度问题时，失效概率计算将不能实现^［10］。直接蒙特卡洛方法效率低，尤其在小失效概率水平时，计算非常困难。因此，亟须发展考虑了变量间的多种相关结构的边坡可靠度高效分析方法。

常用的蒙特卡洛模拟方法除了直接蒙特卡洛方法（MCS）外还有蒙特卡洛重要抽样法（ISM）以及子集模拟方法（SS）。蒙特卡洛重要抽样法是在直接蒙特卡洛方法的基础上通过改变抽样中心令抽样样本更大概率地落入失效域中，从而减少抽样数量，提高计算效率；子集模拟是在直接蒙特卡洛方法的基础上通过引入合理的中间失效域，将小失效概率转换为一系列容易实现的较大条件失效概率的乘积，逐层筛选，大大减少了抽样数目。而子集模拟方法在岩土边坡工程的应用中很少^［11］，仅曹子君^［12］做了基于子集模拟方法的边坡可靠度分析，且文献^［12］中没有考虑岩土体参数间的相关性和相关结构；文献^［12］中用Nataf解决非正态相关变量的变换问题，将变量间的相关结构默认为高斯相关结构，并没有考虑参数间不同相关结构对结果的影响。本文在采用Copula函数表征参数间相关性的基础上，将重要抽样法和子集模拟方法用于含多种相关结构相关随机变量的边坡可靠分析，提出了高效地实施含有相关结构随机变量的边坡可靠度计算的方法。

1 基于Copula函数的边坡可靠度直接蒙特卡洛方法和高效分析方法

1.1 基于Copula函数的边坡岩土体参数联合概率分布模型

Copula理论最早由Sklar于1959年提出。Copula函数是将参数的联合分布函数与其边缘函数联结在一起的函数，不受边缘分布类型的影响，可以构造出具有任意边缘分布类型和任意相关结构的联合分布函数。岩土边坡工程参数常常是二维参数，如抗剪强度参数。所以本文以二维参数为例运用Copula函数构造岩土体参数的联合概率分布模型。根据Copula理论，参数X ₁和X ₂的联合概率分布函数F（x ₁ ，x ₂）及联合概率密度函数 f（x ₁ ，x ₂）分别为：

F（x₁， x₂） = C［F₁（x₁），F₂（x₂）；θ］ = C（u₁， u₂； θ）（1）

f（x₁， x₂） = D（u₁， u₂； θ） f₁（x₁）f₂（x₂）（2）

式中：u ₁ =F ₁（x ₁）和u ₂ =F₂ （x ₂）分别是参数X ₁和X ₂的边缘分布函数；f ₁（x ₁）和f ₂（x ₂）分别是参数X ₁和X ₂的边缘密度函数；C（u ₁，u ₂； θ）和D（u ₁， u ₂； θ）分别是Copula函数的分布函数和密度函数；θ为Copula函数的相关参数。

由式（1）和（2），已知参数X ₁和X ₂的边缘分布和参数间的相关性就可以构造参数X ₁和X ₂的联合概率分布模型。由于岩土体参数间有较强的统计负相关性，而且相关结构基本是对称的，本文依据这种特性选取了4种比较适合描述岩土体参数间相关关系的4种Copula函数：Gaussian Copula、Frank Copula、No.16 Copula和Plackett Copula。其中Gaussian Copula所表征的相关结构即为常用的二维正态或者Nataf变换中隐含的高斯相关结构。有关这4种Copula函数的概率分布函数和参数取值范围参见文献^［10］。

由上可知，采用Copula函数可以表征边坡岩土体参数间的任意相关结构（包括采用Gaussian Copula构建的高斯相关结构和采用其他类型Copula函数构建的非高斯相关结构），从而建立起岩土体参数的联合概率模型。接下来进行可靠度分析的关键，就是选择合适的可靠度分析方法。本文重点提出和上述概率模型结合的蒙特卡洛重要抽样方法（第1.3节）和子集模拟方法（第1.4节）。同时介绍了已有的直接蒙特卡洛方法（第1.2节），以供和其余两种方法对比分析。

1.2 基于Copula函数的直接蒙特卡洛方法

在进行简单失效模式的边坡可靠度计算时，直接蒙特卡罗模拟方法能得到精度较高的结果。直接蒙特卡洛方法的计算过程如下面所叙述的3个步骤：首先，模拟出服从给定Copula函数构造的参数联合分布模型的岩土体参数的N个样本点；然后，将上述模拟的样本点代入边坡失效模式中计算相应的功能函数响应值；最后，统计失效样本数目L从而得到失效概率的无偏估计值P_f = L/N。为了满足计算的失效概率变异系数小于10%，一般样本数量需大于100/P_f，因此当失效概率P_f 很小时，需要的抽样样本数量N很大，计算效率很低。

1.3 基于Copula函数的蒙特卡洛重要抽样方法

直接蒙特卡洛方法编程简单，概念清晰，容易实现，但对于低失效概率边坡可靠度分析问题，其计算效率非常低^［13］。很多研究人员致力于提高蒙特卡洛方法的抽样效率，提出了重要抽样法。其基本思想是提高随机样本出现在贡献率较大的“重要区域”的频率，从而提高抽样效率。基于这个原理我们可以将随机抽样密度中心取在设计点上，因为设计点是失效边界上对失效概率贡献最大的点^［12］。而设计点可由一次可靠度方法求得。本文借鉴这一方法，将其应用于采用Copula函数表征了参数间相关结构的边坡可靠分析中。取正态分布密度函数px为重要抽样密度函数，以设计点为抽样中心，对参数X抽样，得到样本X_i = （x_i ₁， x_i ₂， …， x_in ）（i = 1， 2，…，N），则P_f 的无偏估计值为：

{\hat{P}}_{f} = \frac{1}{N} \sum_{i = 1}^{N} \frac{I [g_{X} (X_{i})] f_{X} (X_{i})}{p x (X_{i})}

（3）

式中：

I (x)

为示性函数，当x＞0时，I（x） = 1；当x≤0时，I（x） = 0；f_X （X）为参数X的联合分布密度函数，结合Copula函数，可由基于Copula函数构造的参数间联合密度函数［即式（2）］代替，将其带入式（3）中即可求得失效概率。

1.4 基于Copula函数的子集模拟方法

子集模拟方法最早由Au和Beck^［14］提出，通过引入中间失效事件，将小失效概率用一系列较大的条件失效概率的乘积表示。子集模拟方法可以一层一层快速逼近功能函数失效域，快速寻找最终失效域，大大降低随机抽样数目，从而提高了可靠度分析效率，抽样过程不涉及设计点的寻找，相比于重要抽样法，进一步提高了计算效率，且不会产生维数灾难，可以更高效率地分析小失效概率岩土边坡稳定性。考虑到边坡岩土体参数间含有相关性，且相关结构类型不定，本文提出基于Copula函数的子集模拟方法来进行边坡可靠度分析。子集模拟概念及实现过程如下：

设岩土边坡失效临界响应值为b，失效域

F

定义为F = ｛g ≤ b｝。使中间事件满足嵌套关系：

F_{1} \supset F_{2} \supset \dots \supset F_{m} = F

，那么目标失效概率可以表示为：

P_{F} = P (F) = P (F_{m}) = P (F_{m} | F_{m - 1}) P (F_{m - 1}) = \dots = P (F_{1}) \prod_{i = 2}^{m} P (F_{i} | F_{i - 1})

（4）

选取中间事件F_i =｛g ≤ b_i， i =1，2…m｝，

m

为中间事件总数，中间事件条件概率

P (F_{i} | F_{i - 1})

足够大，即可以保证高效地模拟P _F。设条件失效概率为P ₀ = P（F_i |F_i _-1），取P ₀ = 0.1，那么由式（4）可得：

P _F = P（F₁）

\prod_{i = 2}^{m} P

（F _i |F _i _-1） = P（F _m ）P₀ ^（ ^m ^-1）（5）

可以看出，当P_F 为10^-4量级的时候，只需要四轮抽样，每轮抽样100个即可达到目标概率，而直接蒙特卡洛方法需要10⁵个抽样样本才能达到目标失效概率，子集模拟极大提高了抽样模拟效率。与Copula函数结合，该方法的具体实现过程如下：

（1）通过蒙特卡洛随机抽样得到n组相互独立的标准正态分布随机变量V = ｛（V ₁ _i，V ₂ _i ）， i=1，2，…，n｝，通过反变换，利用Copula函数表征的参数概率联合分布模型，将V变换为n组服从给定Copula函数相关结构的样本X = ｛（X ₁ _i，X ₂ _i ）， i=1，2，…，n｝，将X带入功能函数g，得到一系列功能函数响应值，并将响应值按由小到大排序，取第nP ₀个功能函数相应值为中间事件域值b_i。前nP ₀组样本是落入第一个失效域F ₁ =｛g ≤ b ₁ ｝的样本点。失效概率为P（F_i |F_i _-1） = P ₀ = 0.1。

（2）分别利用Metropolis-Hastings算法以落在失效域F_i _-1中的样本点V_i _-1为“种子”，产生1/P ₀组新样本V_i。

（3）通过步骤（2）产生的nP ₀ /P ₀ =n组样本

V

，通过Copula函数构造的参数X ₁和X ₂的联合概率密度函数将V变换为n组符合样本相关结构和相关关系的样本X = ｛（X ₁ _i，X ₂ _i ）， i=1，2，…，n｝带入功能函数g计算功能函数响应值，将功能函数响应值按从小到大排序，以第nP ₀个功能函数值为中间事件阈值b_i。

（4）重复（2）、（3）步，直到第m层有超出nP ₀个功能函数值小于零时，终止重复。

（5）统计第m层落入失效域F_m 的样本组数N _F，则可以计算出失效概率为P_F = （N_F /n）P ₀ ^（ ^m ^-1）。

2 功能函数为显式表达式的边坡可靠度分析

本节将所提的考虑了参数间相关结构的高效可靠度分析方法（即重要抽样法和子集模拟方法）用在功能函数为显式表达式的边坡。

2.1 边坡分析模型

以文献^［15］中的一个功能函数为显式的算例进行分析，基本资料：边坡岩体的岩质成分为红色黏土层（SC）、紫红色矽质岩（SQ）以及震旦纪矽质的云岩（SLS），如图1。该边坡岩体发育有一组顺坡产出缓倾角节理，平均产状为N70^oE/SE∠30^o，贯通层面程度不高，岩体结构为顺坡层状结构。边坡岩体节理与黏土岩层面可能构成双滑面剪切破坏模式，前者为被动滑面，后者为主动滑面。

图1 边坡岩体稳定分析简图

Fig.1 Diagram of slope rock mass analysis

Full size|PPT slide

采用刚体极限平衡法计算边坡安全系数FS，据此可得边坡稳定可靠度分析功能函数为：

g = F S - 1 = \frac{f_{1} [G_{1} c o s α_{1} + P s i n (α_{1} + α_{2})] + c_{1} l_{1}}{G_{1} s i n α_{1} + P c o s (α_{2} - α_{1})} - 1

（6）

其中：

P = G_{2} (s i n α_{2} - f_{2} c o s α_{2}) - c_{2} l_{2}

式中：

G_{1}

和

G_{2}

为块体重量；

c_{1}

和

c_{2}

为滑动面上的黏聚力；f ₁和f ₂为滑动面上的摩擦系数；l ₁和l ₂为块体基面长度。

将f ₁、f ₂、c ₁、c ₂、α ₁、α ₂视为随机变量，其中f ₁、f ₂、c ₁、c ₂根据室内试验数据确定，服从对数正态分布，α ₁和α ₂根据现场试验数据确定，服从正态分布，统计参数见表1。确定性参数取值为：块体重量G ₁=46 669 kN，G ₁=34 274 kN；块体基面长度l ₁ = 67.2 m，l ₂=96.8 m；岩体重度γ=27 kN/m³。运用本文用所提的方法：基于Copula函数的蒙特卡洛重要抽样法、基于Copula函数的子集模拟方法计算该显式岩质边坡的失效概率。为了验证所提方法计算得到的失效概率的准确性，也采用直接蒙特卡洛方法计算了该边坡的失效概率，供对比分析。用4种Copula函数（Gaussian、Frank、No.16和Plackett）表征相关性，目的是为了对比分析参数间相关结构对可靠度计算结果的影响。为了简化计算过程，设f ₁和c ₁及f ₂和c ₂两对参数的相关系数相同，α ₁和α ₂相互独立。设定参数间的相关系数为-0.5。计算结果如表2所示。

表1 边坡滑动面强度参数的统计特征

Tab.1 Statistical characteristics of strength parameters of sliding surface of slope

项目	f ₁	f ₂	c ₁/kPa	c ₂/kPa	α ₁	α ₂
均值	0.939	0.277	249.09	39.22	30°	59°08′
标准差	0.123	0.010	84.34	4.90	7°15′	6°20′
分布类型	对数正态	对数正态	对数正态	对数正态	正态	正态

表2 不同可靠度计算方法的边坡可靠度结果

Tab.2 Results using different slope reliability methods

可靠度计算方法	N	Copula函数
		Guassian		Frank		Plackett		No.16
		P_f	误差/%	P_f	误差/%	P_f	误差/%	P_f	误差/%
直接蒙特卡洛(MCS)	1×10⁶	0.027 4	0	0.028 3	0	0.028 7	0	0.032 6	0
重要性抽样(ISM)	1×10⁴	0.022 0	19.71	0.028 6	1.09	0.024 3	15.32	0.030 0	7.98
子集模拟(SS)	2×10²	0.027 2	0.73	0.029 5	4.24	0.029 5	2.79	0.034 4	5.52

注：相关系数为-0.5；“误差”指的是ISM、SS方法计算的P _f与MCS方法计算的P_f 之间的绝对误差。

2.2 边坡可靠度高效性分析

表2中，以N=1×10⁶次直接蒙特卡洛方法计算的可靠度指标为准确值，同时比较了4种Copula函数相关结构下的3种可靠度方法。从表2可以看出，同种相关结构，蒙特卡洛重要抽样方法误差小于20%，子集模拟方法误差小于6%，说明了本文所提方法的有效性和准确性；在模拟效率上，基于Copula函数的重要抽样法只要1×10⁴次就能得到准确的结果，效率是直接蒙特卡洛的100倍，是高效的边坡可靠度分析方法；基于Copula函数的子集模拟只用了200次就可以得到相对准确的结果，这说明在计算精确度较好时，子集模拟效率相比直接蒙特卡洛法和蒙特卡洛重要抽样法有很大提高。

2.3 不同相关结构结果对比分析

对比不同相关结构的结果的差异时，相关系数从-0.1到 -0.9变化，采用最高效的子集模拟方法计算的结果作比较，如图2。可以看出，本算例中采用不同Copula函数计算的失效概率存在差异，这是由于考虑了随机变量间不同相关结构的结果。No.16 Copula函数计算的失效概率高于其他3种Copula函数计算的失效概率，其他3种Copula函数的计算结果十分接近。若以No.16 Copula函数计算的结果为准确值，Gaussian Copula、Plackett Copula和Frank Copula函数计算出来的失效概率偏低，采取这3种Copula函数的计算结果进行可靠度设计会使设计不经济；若以Gaussian Copula函数计算的失效概率为准确值，运用No.16 Copula函数计算出的失效概率结果偏高，采取No.16 Copula函数方法进行可靠度设计会导致设计不安全。因此，运用Nataf变换（即认为边坡岩土体参数间隐含高斯相关结构）的做法可能使得计算结果不准确，而基于Copula函数，可以考虑边坡岩土体参数间任意相关结构，得出更符合实际的边坡失效概率。同时，我们可以看出，该岩土边坡的失效概率随着相关系数的绝对值增大而减小，说明岩土参数间的负相关性对岩土边坡的稳定性有积极影响，忽略这种相关性会导致岩土边坡可靠度设计不经济。

图2 考虑不同相关结构时边坡可靠度的比较(子集模拟计算结果)

Fig.2 Results of subset simulation considering different related structures (SS results)

Full size|PPT slide

3 功能函数为显式表达式的边坡可靠度分析

上一节中边坡的功能函数为显式函数关系。下面以四川省大岗山水电站左岸边坡工程为例^［10］说明本文所提方法在功能函数为隐式函数的边坡可靠度分析中的有效性。

3.1 边坡分析模型

大岗山水电站位于四川省西部大渡河中游石棉县境内^［10］，坝址处控制流域面积达6.27 万km²，占全流域的81%，多年平均流量1 010 m³/s，是大渡河干流近期开发的大型水电工程之一，电站装机容量2 600 MW，最大坝高210 m，电站正常蓄水位 1 130 m，总库容7.42 亿m³.该电站地理位置上处于川滇南北向构造带北段，为北东向、南北向与北西向等多组构造的交汇结合部位。坝区位于由磨西断裂、大渡河断裂和金坪断裂所切割的黄草山断块西侧边缘，大地构造部位属扬子淮台地西部二级构造单元康滇地轴范畴。坝区河谷深切，呈

Ω

形嵌入河曲形态，其新构造运动总体表现为以整体间歇性强烈抬升为主，有作用方向表现为NWW-SEE向挤压的区域构造应力场。已有测试结果显示，大岗山坝区应力场是自重力和构造应力叠加的应力场，且构造应力是坝址左岸应力场的主要组成部分。该边坡边坡稳定性问题突出，其中、陡倾角的顺坡向节理裂隙较为发育，具有卸荷风化强烈、高地应力以及高地震烈度的特点。在不加固情况下，一点小的外来荷载，就很可能产生规模巨大的滑坡和山体变形。大岗山水电站坝址50年超越概率10%基岩水平向加速度为0.251 g，相应的地震基本烈度为VIII度。通过分析坝区VII~VII剖面左岸天然边坡（见图3，图中符号V2、V1、IV、III1、II为岩体质量分级），发现岩脉β ₂₈和断层f ₆₈切割边坡岩体组合成潜在滑动体。将岩脉β ₂₈的抗剪强度f ₁和c ₁，断层f ₆₈的抗剪断强度f ₂和c ₂视为对数边坡岩体重度为γ = 25 kN/m³，左岸边坡滑动面强度参数统计特征值见表3。

图3 坝址区VII~VII剖面左岸天然边坡

Fig.3 Left natural slope of dam site section VII~VII

Full size|PPT slide

表3 左岸边坡滑动面强度参数的统计特征

Tab.3 Statistical characteristics of the strength parameters of the sliding surface of the left bank slope

项目	f ₁	f ₂	c ₁/kPa	c ₂/kPa
均值	0.400 0	0.450 0	30.0	50
标准差	0.600 0	0.076 5	5.4	10
分布类型	对数正态	对数正态	对数正态	对数正态

计算该边坡稳定性时考虑3种工况：天然工况（不考虑降雨和地震因素）、降雨工况（地下水位高度取为滑体高度的10%）、VIII度地震工况（计算中水平向地震加速度为0.251 g）。边坡稳定安全系数采用刚体极限平衡法中的剩余推力法计算，需要反复迭代才能求得，极限状态方程g=FS-1=0，即没有显式表达式。运用本文所提的基于Copula函数的蒙特卡洛重要抽样法和基于Copula函数的子集模拟方法计算该隐式岩质边坡的失效概率。为了验证所提方法计算得到的失效概率的准确性，也采用直接蒙特卡洛方法计算了该边坡的失效概率，供对比分析。此外，为了简化计算，设f ₁和c ₁及f ₂和c ₂两对参数的相关系数相同，在对比不同可靠度计算方法时设定相关系数为-0.5不变，计算结果如表4、表5和表6。

表4 不同方法计算的可靠度结果（天然工况）

Table.4 Results of different slope reliability methods (Natural condition)

可靠度计算方法	N	Copula函数
		Guassian		Frank		Plackett		No.16
		P_f	误差/%	P_f	误差/%	P_f	误差/%	P_f	误差/%
直接蒙特卡洛(MCS)	1×10⁷	0.000 220	0	0.000 389	0	0.000 397	0	0.000 476	0
重要性抽样(ISM)	1×10⁵	0.000 292	24.65	0.000 355	8.74	0.000 345	13.10	0.000 389	18.27
子集模拟(SS)	2×10³	0.000 232	5.45	0.000 358	7.97	0.000 343	13.60	0.000 488	2.52

注：相关系数为-0.5；“误差”指的是ISM、SS方法计算的P_f 与MCS方法计算的P_f 之间的绝对误差。

表5 不同方法计算的可靠度结果（降雨工况）

Table.5 Results of different slope reliability methods (Rainfall condition)

可靠度计算方法	N	Copula函数
		Guassian		Frank		Plackett		No.16
		P_f	误差/%	P_f	误差/%	P_f	误差/%	P_f	误差/%
直接蒙特卡洛(MCS)	1×10⁶	0.001 4	0	0.001 4	0	0.001 3	0	0.003 5	0
重要性抽样(ISM)	1×10⁴	0.001 2	14.3	0.001 9	35.71	0.001 8	38.46	0.002 5	28.57
子集模拟(SS)	8×10²	0.001 2	14.3	0.001 7	21.43	0.001 3	0	0.002 9	17.14

注：相关系数为-0.5；“误差”指的是ISM、SS方法计算的P_f 与MCS方法计算的P_f 之间的绝对误差。

表6 不同方法计算的可靠度结果（地震工况）

Tab.6 Results of different slope reliability methods (Earthquake condition)

可靠度计算方法	N	Copula函数
		Guassian		Frank		Plackett		No.16
		P_f	误差/%	P_f	误差/%	P_f	误差/%	P_f	误差/%
直接蒙特卡洛(MCS)	1×10⁷	0.271 7	0	0.268 5	0	0.268 5	0	0.264 6	0
重要性抽样(ISM)	1×10⁵	0.299 5	10.23	0.296 7	10.50	0.296 7	10.50	0.291 0	9.97
子集模拟(SS)	2×10³	0.305 6	12.48	0.308 3	14.82	0.308 3	14.82	0.305 3	22.41

注：相关系数为-0.5；“误差”指的是ISM、SS方法计算的P_f 与MCS方法计算的P_f 之间的绝对误差。

3.2 边坡可靠度高效性分析

由表4、表5和表6所示的任意一种工况的计算结果可以看出，以直接蒙特卡洛方法计算的失效概率为准确值，同一种Copula函数相关结构条件下，蒙特卡洛重要抽样方法和子集模拟方法都能得到相对准确的结果，绝大部分情况失效概率计算结果误差小于20%，精度最差的为降雨工况时蒙特卡洛重要抽样法计算的结果，误差为38.46%，说明了本文所提方法处理含隐式功能函数表达式的边坡可靠度问题的有效性；在模拟效率上，蒙特卡洛重要抽样方法相比于直接蒙特卡洛方法，抽样次数是其1%，效率很高；而子集模拟方法的抽样模拟次数仅需800~2 000 次，相比于蒙特卡洛重要抽样方法的计算效率又进一步提高。此外，还可以看出，大岗山水电站坝址处左岸边坡在天然工况下失效概率在考虑几种相关结构的情况下最大为4.8×10^-4左右，失效风险非常小，在天然状态下具有足够的安全度；而在地震工况下失效概率达到了0.27，该左岸边坡非常不安全，有必要对该左岸边坡增加锚固措施，防止失效破坏，危及人民生命和财产安全。

3.3 不同相关结构结果对比分析

同时为了对比不同相关结构对可靠度分析结果的影响，在天然工况下，设定抗剪断强度参数f ₁和c ₁及f ₂和c ₂的相关系数从-0.1~-0.9变化，选取最高效的子集模拟方法的结果，结果如图4。由图4可以看出，随着相关系数变小，即参数间的负相关性增加，失效概率下降，说明岩土体参数间的负相关性对岩土体稳定性有积极作用。此外，不同Copula函数计算的失效概率有着显著差异。No.16 Copula函数计算的失效概率相比其他3种Copula函数的计算结果较高，Gaussian Copula函数计算的失效概率偏低。若以Plackett Copula函数计算的失效概率为准确值，运用No.16 Copula函数的计算结果对该进行边坡可靠度设计，会造成设计方案偏保守，而如果运用Gaussian Copula函数的计算结果作该边坡可靠度设计，会造成设计方案不安全的结果。因此，运用Nataf变换（即默认为Gaussian Copula表征的相关性）认为边坡岩土体参数间隐含高斯相关结构的做法会使得计算结果不准确，而基于Copula函数，可以考虑边坡岩土体参数间任意相关结构，能得出更合理的失效概率。此外，由于NO.16 Copula函数计算的失效概率偏高，应用该函数的计算结果进行边坡可靠度设计会偏于保守，设计更安全，所以对于该类工程边坡，在进行边坡可靠度分析时建议应用NO.16 Copula函数。

图4 考虑不同相关结构时边坡可靠度的比较（天然工况，子集模拟计算结果）

Fig.4 Results of subset simulation considering different related structures (Natural condition; SS results)

Full size|PPT slide

4 结论

（1）引入Copula函数表征岩土体参数间的相关性并进行边坡可靠度分析十分必要。Copula函数可以表征参数间包括高斯相关结构在内的任意相关结构。不同的相关结构类型情况下，边坡失效概率差别较大。在岩土体参数间具体相关结构未知时，需要谨慎选取Copula函数类型。有试验数据时，建议根据数据选择最适合的Copula函数类型。盲目地选取任意Copula函数都会导致设计结果不合理。尤其需要指出的是，如果为了简便而运用Gaussian Copula函数的计算结果做边坡可靠度设计，会造成设计方案不安全或不经济。

（2）所提的含相关结构随机变量的边坡可靠度高效分析方法可以克服直接蒙特卡洛方法效率低的缺陷，且计算精度有保障。尤其是子集模拟算法在基于Copula函数表征了变量间的相关结构的边坡可靠度分析中效率极高，只需模拟少量样本即可获得失效概率值。同时，该方法能有效处理含隐式功能函数边坡可靠度分析问题，进一步拓展了高效蒙特卡洛方法在边坡可靠度分析中的应用。

（3）大岗山水电站坝址处左岸边坡在自然工况和降雨工况下失效概率很小，一般情况下有足够的安全度；但在地震工况下的失效概率很大，应对其进行锚固支护措施提高其稳定性。

参考文献

原文顺序 | 文献年度倒序 | 文中引用次数倒序

1	PHOON K K， KULHAWY F H， GRIGORIU M D. Development of a reliability-based design framework for transmission line structure foundations［J］. Journal of Geotechnical and Geoenvironmental Engineering， 2003，129（9）：798-806. 本文引用 [1]

2	FENTON G A， NAGHIBI F， DUNDAS D， et al. Reliabilitybased geotechnical design in 2014 Canadian Highway Bridge Design Code［J］. Canadian Geotechnical Journal， 2016，53 （2）：236-251. 本文引用 [1]

3	陈祖煜，陈立宏，王玉杰，等. 滑坡和建筑物抗滑稳定分析中的可靠度分析和分项系数设计方法［C］//高坝大库安全建设与风险管理高端论坛，2011. 本文引用 [1]

4	王禹玺，李典庆，唐小松.边坡系统可靠度区间估计方法［J］.武汉大学学报（工学版），2019，52（11）：941-950. 本文引用 [1]

5	SILVA F， LAMBE T W， MARR W A. Probability and risk of slopefailure［J］. Journal of Geotechnical & Geoenvironmental Engineering， 2008，134（12）：1 691-1 699. 本文引用 [1]

6	唐小松，李典庆，周创兵，等. 不完备概率信息条件下边坡可靠度分析方法［J］. 岩土工程学报，2013，35（6）：1 027-1 034. 本文引用 [1]

7	范明桥，盛金保. 土强度指标φ， c的互相关性［J］.岩土工程学报，1997（4）：100-104. 本文引用 [2]

8	LOW B K. Reliability analysis of rock slopes involving correlated nonnormals ［J］. International Journal of Rock Mechanics and Mining Sciences， 2007，44（6）.

9	TANG X S， LI D Q， CHEN Y F， ZHOU C B， ZHANG L M. Improved knowledge-based clustered partitioning approach and its application to slope reliability analysis［J］. Computers and Geotechnics， 2012，45. 本文引用 [1]

10	唐小松. 基于Copula理论的岩土体参数不确定性建模与可靠度分析［D］. 武汉：武汉大学，2014. 本文引用 [6]

11	张曼，唐小松，李典庆.含相关非正态变量边坡可靠度分析的子集模拟方法［J］.武汉大学学报（工学版），2012，45（1）：41-45，53. 本文引用 [1]

12	曹子君.子集模拟在边坡可靠性分析中的应用［D］. 成都：西南交通大学，2009. 本文引用 [4]

13	张明.结构可靠度分析：方法与程序= Structural Reliability Analysis： Methods and Procedures［M］. 北京：科学出版社，2009. 本文引用 [1]

14	李洪双，马远卓. 结构可靠性分析与随机优化设计的统一方法［M］. 北京：国防工业出版社，2015. 本文引用 [1]

15	祝玉学. 边坡可靠性分析［M］. 北京：冶金工业出版社，1993. 本文引用 [1]

PDF(1415 KB)

647

访问

引用

详细情况

段落导航

摘要
Abstract
关键词
Key words
基金
引用本文
0 引言
1 基于Copula函数的边坡可靠度直接蒙特卡洛方法和高效分析方法
1.1 基于Copula函数的边坡岩土体参数联合概率分布模型
1.2 基于Copula函数的直接蒙特卡洛方法
1.3 基于Copula函数的蒙特卡洛重要抽样方法
1.4 基于Copula函数的子集模拟方法
2 功能函数为显式表达式的边坡可靠度分析
2.1 边坡分析模型
图1 边坡岩体稳定分析简图
表1 边坡滑动面强度参数的统计特征
表2 不同可靠度计算方法的边坡可靠度结果
2.2 边坡可靠度高效性分析
2.3 不同相关结构结果对比分析
图2 考虑不同相关结构时边坡可靠度的比较(子集模拟计算结果)
3 功能函数为显式表达式的边坡可靠度分析
3.1 边坡分析模型
图3 坝址区VII~VII剖面左岸天然边坡
表3 左岸边坡滑动面强度参数的统计特征
表4 不同方法计算的可靠度结果（天然工况）
表5 不同方法计算的可靠度结果（降雨工况）
表6 不同方法计算的可靠度结果（地震工况）
3.2 边坡可靠度高效性分析
3.3 不同相关结构结果对比分析
图4 考虑不同相关结构时边坡可靠度的比较（天然工况，子集模拟计算结果）
4 结论
参考文献

收稿日期	出版日期
2020-04-09	2021-04-15
发布日期
2021-04-28

选择文件类型/文献管理软件名称

选择包含的内容