一种基于蒙特卡洛模拟的边坡稳定可靠度敏感性分析方法
一种基于蒙特卡洛模拟的边坡稳定可靠度敏感性分析方法
【技术领域】
[0001] 本发明涉及一种边坡参数敏感性分析方法，尤其涉及一种基于蒙特卡洛模拟的边 坡稳定可靠度敏感性分析方法。
【背景技术】
[0002] 基于可靠度理论的边坡敏感性分析能有效地表征输入参数的不确定性及其统计 特征（均值、方差等）对边坡稳定可靠度的影响，对边坡设计、加固等具有重要指导意义，近 年来受到广泛关注。蒙特卡洛模拟方法作为一种求解可靠度的随机模拟方法，因其概念简 单、适用性强在岩土工程可靠度领域得到广泛应用。然而蒙特卡洛模拟方法存在不确定性 传播机理不明确的缺点，无法直接获得失效概率随着输入参数不确定性变化的响应规律。 因此，在利用蒙特卡洛模拟进行不确定性参数的敏感性分析时，通常需要重复模拟，计算输 入随机变量统计特征的一系列尚散点对应的失效概率，然后分析统计特征与失效概率之间 的变化规律。在工程实践中，一般边坡设计方案对应的失效概率较小（比如，低于0.001)， 而对于小概率问题，蒙特卡洛模拟方法计算效率较低，重复模拟需要耗费大量的时间和计 算机资源，所需的计算时间和资源随着边坡稳定性分析模型的复杂程度的增加而增加。
[0003]Wang[1]等提出了一种基于蒙特卡洛模拟的敏感性分析方法，该方法能够有效地 定量分析不确定性参数对于失效概率的相对贡献以及计算不同不确定性参数值对应的条 件失效概率，但是未对失效概率与不确定性参数的统计特征间的响应规律进行研宄。随后 Wang[2]提出了一种结合贝叶斯理论的敏感性分析方法，该方法不需要重新执行蒙特卡洛模 拟，就可以计算单个不确定性参数的统计特征变化后边坡的失效概率，获得失效概率与不 确定性参数的统计特征之间的响应规律。但是，该方法需要制作频率直方图，计算过程较为 复杂，且无法解决多个变量统计特征同时变化的敏感性分析问题。
[0004] 文中涉及如下参考文献：
[0005] [l]ffangY,CaoZJ,AuSK.EfficientMonteCarloSimulationofparameter sensitivityinprobabilisticslopestabilityanalysis[J].Computersand Geotechnics, 2010, 37(7 - 8):1015 - 1022.
[0006] [2]WangY.UncertainparametersensitivityinMonteCarloSimulationby samplereassembling[J].ComputersandGeotechnics, 2012, 46:39 - 47.
[0007] [3]KiureghianAD,LiuPL.Structuralreliabilityunderincomplete probabilityinformation[J].JournalofEngineeringMechanics, 1986, 112:86- 114.
【发明内容】
[0008] 针对现有技术存在的不足，本发明基于蒙特卡洛模拟提出了一种计算过程简便、 计算效率高的边坡稳定可靠度敏感性分析方法。
[0009] 为解决上述技术问题，本发明采用如下的技术方案：
[0010] 基于蒙特卡洛模拟的边坡稳定可靠度敏感性分析方法，包括步骤：
[0011] 步骤1，确定不确定性参数的统计特征，并构造不确定性参数的联合概率密度函数 f(x)；
[0012]步骤2,基于联合概率密度函数f(x)，利用蒙特卡洛模拟方法获得边坡失效概率 Pf，并获得失效样本；
[0013] 步骤3,设计多种敏感性分析方案，根据敏感性分析方案对应的统计特征，分别构 造各敏感性分析方案下不确定性参数的联合概率密度函数fk(x);
[0014] 所述的敏感性分析方案采用如下方法获得：基于步骤1确定的统计特征，考虑不 确定参数的实际变化范围，改变一个或多个不确定性参数的统计特征值，保持其它不确定 性参数的统计特征概率分布不变，即获得一种敏感性分析方案；
[0015] 步骤4,获得各敏感性分析方案下的边坡失效概率，本步骤进一步包括：
[0016] 4. 1根据联合概率密度函数fk(x)和失效样本，获得失效样本的权重指标
f(\)和fk(\)分别为第j个失效样本的权重指标、不确定性参数联合概 率密度函数和第k个敏感性分析方案下的联合概率密度函数；
[0017] 4. 2根据权重指标获得敏感性分析方案的边坡失效概率
，P/为第k 个敏感性分析方案的边坡失效概率，^为失效样本数，N为步骤2蒙特卡洛模拟方法中产生 的随机向量样本数；
[0018] 步骤5,根据各敏感性分析方案的边坡失效概率，获得边坡失效概率随不确定性参 数统计特征的变化趋势，从而识别影响边坡稳定可靠度的关键不确定性参数。
[0019] 步骤1进一步包括子步骤：
[0020] 1.1根据不确定性参数的试验数据，或现有文献中不确定性参数的统计特征取值， 确定不确定性参数的统计特征以及各不确定性参数间的相关系数；
[0021] 1. 2根据不确定性参数的统计特征及各不确定性参数间的相关系数，构造不确定 性参数的联合概率密度函数f(x)。
[0022] 子步骤1. 2具体为：
[0023] (1)若不确定性参数全部相互独立，则不确定性参数的联合概率密度函数为各不 确定性参数概率密度函数的乘积；
[0024] (2)若不确定性参数全部相关，则基于Copula理论或Nataf变换法构造不确定性 参数的联合概率密度函数；
[0025] (3)若不确定性参数部分独立、部分相关，首先，采用步骤（1)中方法获得相互独 立的不确定性参数的联合概率密度函数，记为
，其中，q为相互独立的不确定性参 数个数；然后，采用步骤（2)中方法获得相关的不确定参数的联合概率密度函数f(xq+1，一 ，xn);则不确定性参数的联合概率密度函数
[0026] 步骤2进一步包括子步骤：
[0027] 2. 1产生服从于联合概率密度函数f(x)的随机向量样本；
[0028] 2. 2对随机向量样本进行边坡稳定性分析，判断各随机向量样本对应的边坡是否 失效；
[0029] 2. 3统计失效样本数，获得边坡失效概率
其中，^为失效样本数，N为随机 向量样本数；所述的失效样本即失效边坡对应的随机向量样本；
[0030] 2. 4输出并保存失效样本。
[0031] 与现有技术相比，本发明具有以下优点和有益效果：
[0032] 1、适用范围广：
[0033] 本发明基于蒙特卡洛模拟，概念清晰、简单易懂。此外，本发明能够有效地解决单 一不确定性参数和多种不确定性参数的边坡稳定可靠度敏感性分析问题，同时也能较好地 适用于功能函数表达式为隐式的复杂边坡稳定可靠度敏感性分析问题，适用范围广。
[0034] 2、计算过程简单，计算效率高：
[0035] 传统蒙特卡洛模拟方法在计算敏感性分析方案所对应的边坡失效概率时，需重复 执行蒙特卡洛模拟，计算效率较低，尤其对于复杂边坡的稳定可靠度敏感性分析而言，蒙特 卡洛模拟中边坡确定性分析所需时间长，重复执行蒙特卡洛模拟进行敏感性分析计算时间 长，对计算资源要求高。本发明充分利用不确定性参数原概率分布蒙特卡洛模拟的计算信 息，在求解敏感性分析方案失效概率过程中，只需计算失效样本的权重指标，避免了蒙特卡 洛模拟的重复执行，简化了计算过程，提高了计算效率，是一种简单高效的边坡稳定可靠度 敏感性分析方法。尤其对于上述复杂边坡的稳定可靠度敏感性分析而言，本发明对计算效 率的提高显得尤为突出，有利于推动蒙特卡洛模拟方法在边坡可靠度工程的应用。
[0036] 3、揭示不确定性传播机理，为边坡工程设计提供重要的参考依据：
[0037] 本发明能有效地揭示边坡可靠度水平与不确定性参数的统计特征间的响应关系， 从而明确了蒙特卡洛模拟过程中不确定性的传播机理，对边坡风险控制、设计优化等具有 一定的指导意义。
[0038] 综上，本发明可简化计 算过程，提高计算效率，揭示蒙特卡洛模拟过程中不确定性 的传播机理，对推动岩土工程可靠度与风险理论在工程实践中应用有着重要的实用价值。
【附图说明】
[0039] 图1为本发明的具体流程示意图；
[0040] 图2为基于蒙特卡洛模拟的边坡可靠度分析流程图；
[0041] 图3为实施例中边坡的剖面图；
[0042] 图4为实施例中不同敏感性方案的不确定性参数的概率分布图，以硬质粘土厚度 T"为例；
[0043] 图5~7为实施例中敏感性分析结果。
【具体实施方式】
[0044] 下面将结合附图和【具体实施方式】进一步说明本发明。
[0045] 本发明的具体流程见图1，具体步骤如下：
[0046] 步骤1，确定不确定性参数的统计特征，并构造不确定性参数的联合概率密度函数 f(X)，即不确定性参数原分布的联合概率密度函数f(X)，X为表征不确定性参数的随机向 量。
[0047] 岩土工程可靠度分析中，土体参数可分为两类参数：确定性参数和不确定性参数， 不确定性参数包括土体粘聚力、内摩擦角、土体容重、土层厚度等。对于不确定性参数，其不 确定性的量化可根据有无试验数据分为两种情况：一，若边坡设计中有现场或室内试验结 果，可根据不确定性参数的试验数据，采取统计方法量化其不确定性。比如对于不确定性参 数土体重度而言，在现场取多个测点试样，根据各测点试样的土体重度，利用统计方法确定 其统计特征，例如均值、方差等，然后根据获得的统计特征量化土体重度的不确定性。二，若 没有关于不确定性参数的试验资料，可根据现有文献不确定性参数的统计特征取值，对所 研宄边坡中相应不确定性参数的统计特征进行合理的设定。
[0048] 本步骤进一步包括以下子步骤：
[0049] 1. 1根据不确定性参数的试验数据，或现有文献不确定性参数的统计特征取值，确 定不确定性参数的统计特征以及各不确定性参数间的相关系数P。
[0050] 1. 2根据不确定性参数的统计特征以及各不确定性参数间的相关系数，构造不确 定性参数的联合概率密度函数f(x)。
[0051] 联合概率密度函数f(X)的构造可分为三种情况：（1)不确定性参数全部相互独立 情况下联合概率密度函数的构造；(2)不确定性参数全部相关情况下联合概率密度函数的 构造；（3)不确定性参数部分独立、部分相关的情况下联合概率密度函数的构造。
[0052] 对于相互独立的不确定性参数而言，n维不确定性参数的联合概率密度函数f(x) 为各不确定性参数的概率密度函数的乘积，因此f(x)可以表示为：
[0054] 式⑴中：f(Xi)为第i个不确定性参数\的概率密度函数，其中取值， n为不确定性参数个数，n表示连乘运算。
[0055] 对于相关的不确定性参数而言，联合概率密度函数的构造可采用Copula理论、Nataf?变换法等。本具体实施中采用较常用的Nataf?变换法构造相关的不确定性参数的联 合概率密度函数。现以二维相关不确定性参数XJPX2为例，介绍Nataf变换法的基本原 理。
[0056] 将相关的不确定性参数&和X2的概率密度函数分别记为f(xJ和f(x2)，其中，Xl 和&分别表示X:和X2的取值，其累计分布函数分别记为Fi(Xl)和F2 (x2)，不确定性参数& 和&间的皮尔逊相关系数为P。对不确定性参数X:和X2进行如下等概率变换：
[0058] 式（2)中(yi)为标准正态变量^的累计分布函数，〇<( ?)为标准正态变量累 计分布函数〇 ( ?)的反函数，1为X,等概率变换后对应的标准正态分布变量，x,和y1表 示XJPY^取值。
[0059] 根据Nataf变换法，可推导不确定性参数XJPX2的联合概率密度函数f(x)为：
[0061] 式⑶中：P。为变量YjpY2的皮尔逊相关系数；<i> ( ?)为标准正态分布的概率 密度函数；巾2 (?)为二维标准正态分布联合概率密度函数，其表达式为：
[0063] 根据相关系数定义，P和口^具有关系：
[0065] 式（5)中：yjpy2分别为不确定性参数XjPX^均值；〇 〇 2分别为不确 定性参数XJPX2的标准差。
[0066]变量YjPY2的相关系数P。可通过式（5)迭代求解，也可以根据DerKiureghian 和Liu[3]提供的针对不同概率统计特征的经验公式求解，还可以根据高斯-埃尔米特积分 求解。但是一般情况下，pjPP相差不大，在近似求解中可以直接取PP。
[0067] 需要指出的是，岩土体不确定性参数通常服从正态分布。若不确定性参数XJPX^ 服从正态分布，此时根据式（3)和（4)，不确定性参数XJPX2的联合概率密度函数f(x)可 以表达为：
[0069] 对于部分独立、部分相关的不确定性参数而言，根据公式（1)获得独立的不确定 性参数的联合概率密度函数，记为1
?，根据上述Nataf变换方法或Copula理论可获 得相关的不确定参数的联合概率密度函数f(xq+1，…，xn)，故部分独立部分相关情况下不确 定性参数的联合概率密度函数：
，其中，f(Xi)表示独立不确定 性参数Xi,(i= 1，. . .q)的概率密度函数，Xq+pX#、…Xn表示相关的不确定性参数。
[0070] 步骤2,根据步骤1构建的联合概率密度函数f(x)，利用蒙特卡洛模拟方法计算边 坡失效概率Pf，并输出和保存失效样本j= 1，2, . . .，ns。
[0071] 本步骤采用传统蒙特卡洛模拟方法实现，具体流程见图2,主要包括如下子步骤：
[0072] 2. 1根据步骤1构建的联合概率密度函数f(x)，产生服从于联合概率密度函数 f(x)的随机向量样本Xpi= 1，2,...N，i表示随机向量样本编号。
[0073] 2. 2构建边坡稳定性分析模型，计算随机向量样本xj#应的边坡功能函数值g(xJ =FS(Xi)-l，并根据功能函数值g(Xi)判断边坡是否失效。
[0074] 目前常用的边坡稳定性分析法有极限平衡法和有限元法，其中极限平衡法包括 瑞典条分法、简化毕肖普法、摩根斯坦-普赖斯法等。用于边坡稳定性的商业软件主要有 GeoStudio(极限平衡法）、ANSYS(有限元法），ABAQUS(有限元法）等。本发明适用任何一 种边坡稳定性分析法。
[0075] 对于具体的边坡工程，根据要采用的边坡稳定性分析方法，运用相应的商业软件 建立边坡稳定性分析模型。随后，代入随机向量Xi，计算其对应的边坡功能函数值g(xj= FS(xj-1,当边坡功能函数值g(xj>0,即边坡安全系数FS(xD小于1时，表示边坡处于稳定 状态；当g(xD= 0时，表示边坡处于临界状态；当g(x)〈0时，表示边坡失稳，S卩边坡失效。
[0076] 循环执行子步骤2. 1和子步骤2. 2,直至i=N。
[0077] 2. 3统计失效样本数目ns，计算边坡失效概率
[0078] 满足g(Xi)〈0的随机向量样本Xi为失效样本。统计失效样本数量计为ns，失效样 本记为j= 1，2,. . .，ns。根据蒙特卡洛模拟方法原理，边坡失效概率Pf可以表示为：
[0080] 式（7)中，巧表示边坡失效概率估值，I[g(x)]为指示函数，当g(x)〈0时， I[g (X)]取 1;当g (x)彡 0 时，I[g (x)]取为 0。
[0081] 2. 4输出并保存失效样本Xj，j= 1，2, ? ? ?，ns。
[0082] 步骤3,基于步骤1确定的不确定性参数统计特征，通过增加和减小不确定性参数 统计特征获得一系列敏感性分析方案，根据各敏感性分析方案对应的不确定性参数的统计 特征构造联合概率密度函数。
[0083] 针对所研宄的不确定性参数，通过改变其统计特征值，在确保其不超出相应不确 定性参数的实际变化范围前提下，设置M种敏感性分析方案；根据各敏感性分析方案对应 的统计特征，分别构造相应的不确定性参数联合概率密度函数fk (x)，k= 1，2,...M，M为敏 感性分析方案数量，k表示 敏感性分析方案编号，fk(x)表示第k个敏感性分析方案相应的 不确定性参数的联合概率密度函数。
[0084] 利用蒙特卡洛模拟方法分析不确定性参数敏感性时，需要通过计算输入不确定性 参数统计特征（如均值，方差）的一系列离散点的边坡失效概率，以探宄边坡失效概率和不 确定性参数的统计特征间的响应规律。
[0085] 如对土体粘聚力参数进行敏感性分析时，可通过增加和减小粘聚力的均值y或 标准差〇，设置M种粘聚力的不同分布作为敏感性分析方案。然后，根据各敏感性分析方案 下的统计特征，构造各敏感性分析方案下的不确定性参数联合概率密度函数fk(x)。联合概 率密度函数fk(x)的具体构造方法同子步骤1.2。
[0086] 这里需要指出的是，本发明能够对多个不确定性参数同时进行敏感性分析，如可 以设置一系列同时改变多个不确定性参数的均值y和标准差〇的敏感性分析方案。
[0087] 步骤4,利用步骤3获得的联合概率密度函数fk(x)和步骤2获得的失效样本Xj， j= 1，2,. . .，ns，求解各敏感性分析方案下的边坡失效概率Pfk，根据边坡失效概率Pfk随不 确定性参数统计特征的变化规律对不确定性参数的统计特征和变异系数进行敏感性分析。
[0088] 根据边坡失效概率通用公式，第k个敏感性分析方案下的边坡失效概率Pfk可表示 为：
[0090] 式⑶中为在联合概率密度函数fk(x)下的数学期望算子。
[0091] 结合不确定性参数原分布的联合概率密度函数f(x)，即步骤1构造的联合概率密 度函数f(x)，第k个敏感性分析方案下的边坡失效概率Pfk可进一步表示为：
[0093] 式（9)中，《k为本发明自定义的权重指标，其计算表达式如下：
[0095] 通过权重指标《k可以把求解敏感性分析方案条件下的边坡失效概率Pfk转化为 求解I[g(x) ] ?15在f(x)分布下的数学期望。因此，可以直接利用步骤2. 1中蒙特卡洛模拟 的随机样本Xi，i=l，2, ? ? *，N，计算边坡失效概率Pfk。故基于联合概率密度函数f(x) 产生的N组随机样本求解第k个敏感性分析方案的边坡失效概率Pfk可表示为：
[0097] 考虑到指示函数I[g(x)]的取值特征，仅失效样本Xj,j= 1，2，? ? ?，ns会对失 效概率Pfk产生贡献，公式（11)可进一步表达为：
[0099] Pfk估计值的方差]可表示为：
[0101] 的变异系数COV[$]可表示为：
[0103] #的变异系数cov[#]:能够有效地表征计算结果的准确性，为本发明方法计算 结果的可靠性提供依据。
[0104] 由式（12)可以看出，在求解敏感性分析方案的边坡失效概率时，本发明只需要计 算原分布产生的失效样本处的权重指标并对其求和，无需重新进行蒙特卡洛模拟，计算过 程较为简单、计算效率高。
[0105] 需要说明的是，本步骤中计算失效概率所用联合概率密度函数为失效样本范围内 的联合概率密度函数，有别于步骤1到3中全体样本范围内的联合概率密度函数。
[0106] 本步骤的具体子步骤如下：
[0107] 本步骤的具体子步骤如下：
[0108] 4. 1根据失效样本的联合概率密度函数f(Xj)和fk (Xj)，获得第j个失效样本Xj的 权重指标
，勿、f(xp和fk(xp分别为第j个失效样本的权重指标、不确定性 参数联合概率密度函数和第k个敏感性分析方案下的联合概率密度函数。
[0109] 4. 2根据权重指标 < 计算敏感性分析方案的边坡失效概率
[0110] 步骤5,根据敏感性分析方案的边坡失效概率，绘制边坡失效概率随不确定性参数 统计特征变化的趋势曲线，识别影响边坡稳定可靠度的关键不确定性参数。
[0111] 实施例1
[0112] 一、工程简介
[0113] 詹姆斯湾堤坝位于加拿大魁北克地区的一项水利工程，堤坝长50km，其中一种设 计方案剖面如图3所示，堤坝高12m，中间有一个56m宽平台，边坡坡角为18.4°，坡比为 3:1。坝体以下土层从上至下依次为硬质粘土、海成粘土、湖积粘土和冰碛层。
[0114] 参考文献[1]，本实施例考虑6个随机变量，即待研宄的不确定性参数，分别为硬 质粘土厚度、海成粘土不排水抗剪强度、湖积粘土的不排水抗剪强度S&、坝基到冰碛 层的距离Dmi、筑坝填土的摩擦角<i>Fill和重度YFill。不确定性参数均服从正态分布且相 互独立，相应的统计特征见表1。其它确定性参数如下：海成粘土的厚度为8m，硬质粘土 不排水抗剪强度为41kPa，硬质粘土、海成粘土、湖积粘土的重度分别为19kN/m3、19kN/m3、 20. 5kN/m3。此外，由几何关系可知湖积粘土厚度IY=Dmi-Tra-8。
[0115] 为简化计算，文献[1]对最危险滑动面增加两个限制：（1)最危险滑动面与冰碛层 相切并过点（X，y)，其中，x= 4. 9m，y= 36m; (2)最危险滑动面圆心横坐标x= 85. 9m。这 种处理方法使得最危险滑动面的位置仅与坝基到冰碛层的距离Dmi有关。
[0116] 表1不确定性参数的统计特征
[0117]
[0118] 二、具体实施流程
[0119] 步骤1，根据表1中各不确定性参数的统计特征，构造不确定性参数原分布的联合 概率密度函数f(ci)Fill，丫Fill, Tcr，SuL，Dml)〇
[0120] 由于本实施例中各不确定性参数服从正态分布且相互独立，因此联合概率密度函 数MFill，TFin?TcrJS^,SuLJDTill )可以表示为：
[0121] f( ^Fiii)TFill)Tcr,SuM,SuL,DTill) -f( 4) Fill)f(yFill)f(Tcr)f(S^)f(SuL)f(DTill) (15)
[0122] 其中：
[0124] 式中，yp y2、y3、y4、y5、y6和o工、o 2、o 3、o 4、o 5、o6分别为参数（J)Fm、 YFill、T。,、SfSupDm^均值和标准差，其具体数值见表1。
[0125] 步骤2,根据不确定性参数原分布的联合概率密度函数f(<i)Fill，yFill，Tra，S^，SuU Dmi)，利用蒙特卡洛模拟法获得詹姆斯湾堤坝的边坡失效概率Pf，并输出和保存失效样本 Xj，j l，2，...，ris〇
[0126]2. 1利用VBA程序产生服从于联合概率密度函数f(d>Fill， 的随机向量样本Xi，i= 1，2,. . .N，i表示随机向量样本编号，N= 2X105。
[0127] 2. 2构建边坡稳定性分析模型，计算随机向量样本xj#应的边坡功能函数值g(xJ =FS(Xi)_l，并根据边坡功能函数值g(Xi)判断詹姆斯湾堤坝是否失效。
[0128] 首先，采用简化毕肖普法，利用EXCEL表单构建詹姆斯湾堤坝稳定性分析模型。然 后，将随机向量样本作为输入参数带入詹姆斯堤湾坝的稳定性分析模型中。最后，根据稳定 性分析模型计算随机向量变量Xi对应的边坡安全系数FS(xJ和边坡功能函数值g(Xi)。
[0129] 循环执行子步骤2. 1和子步骤2. 2,直至i=N。
[0130] 2. 3统计失效样本数\，共有失效样本\= 452个。根据式（7)计算詹姆斯湾堤 坝的边坡失效概率Pf= 452/200000 = 2. 26X10 '
[0131] 2. 4将失效样本&储存在EXCEL表单内，供敏感性分析使用。
[0132] 以上基于不确定性参数原分布的蒙特卡洛模拟法（20万次）计算时间约10小时， 所采用计算机配置为内存4GB、CPUIntelCorei3和主频3. 3GHz。因此，如果采用传统的 蒙特卡洛模拟法求解詹姆斯湾堤坝的敏感性分析问题，对于每种敏感性分析方案，在相同 的模拟次数下所需要的计算时间都约为10个小时。由此可见，采用传统蒙特卡洛模拟法求 解边坡稳定可靠度敏感性分析问题，计算效率非常低。
[0133] 步骤3,通过改变所研宄 不确定性参数的统计特征值，设置一系列敏感性分析方 案。根据各敏感性分析方案对应的统计特征，构造各敏感性分析方案下不确定性参数的联 合概率浩度函数f(巾Fiii，yFiii，Tcr,S-SuL,DTill)，f(<}> Fill,yFill,Tcr,S-Si,DTill)表不弟k 种敏感性分析方案下的不确定性参数联合概率密度函数。
[0134] 本步骤通过改变某不确定性参数的均值、方差或其他统计特征，来获得多种敏感 性分析方案，并通过分析各敏感性分析方案的边坡失效概率，获得该不确定性参数的统计 特征对边坡稳定可靠度的影响程度，从而得到影响边坡稳定可靠度的关键不确定性参数。
[0135] 由于本实施例中不确定性参数较多，此处仅选取硬质粘土厚度、冰碛层厚度 Dmi和湖积粘土的不排水抗剪强度S&，研宄的均值变化对詹姆斯湾堤坝边坡失 效概率的影响。如表2所示，本实施例中通过改变不确定性参数UmpSi的均值，共设置 18种敏感性分析方案：方案1~6中，硬质粘土厚度的均值分别设为u3-3 〇 3、u3-2 〇 3、 u3- 〇 3、u3+ 〇 3、u3+2 〇 3、u3+3 〇 3,u3、〇 3分别为T原分布的均值和标准差，其概率分布曲线 见图4;其他不确定性参数的统计特征值均保持原分布的统计特征值。类似的，分别对参数 Sul和Dmi做同样的处理，可以得到方案7~12和方案13~18。
[0136] 各敏感性分析方案的联合概率密度函数fk((i)mi，YFill，Tra，s^，SuUDmi)的构造 与案例分析步骤1类似，在此不再赘述。
[0137] 表2敏感性分析方案
[0138]
[0139] 步骤4,根据求得的联合概率密度函数fk(<i>Fill，YFill，T。,，S^，DTill)和步骤2中 输出的失效样本\，j=l，2, ? ? ?，&，求解各敏感性分析方案的边坡失效概率Pfk。
[0140] 根据式（10)，第k个敏感性分析方案对应的权重指标《k可表示为：
[0142] 如当k取1时，权重指标o1可以表示为
[0144] 失效样本的权重指标值带入式（12)，则可以计算出各敏感性分析方案所对应 的边坡失效概率。
[0145] 在此需要指出的是，权重指标的计算可以结合EXCEL内置的概率密度函数（如正 态分布概率密度函数NORM.DIST)，计算十分简便。同时，由步骤3~4可以看出，在求解敏 感性分析方案的过程中，由于充分利用了原概率分布蒙特卡洛模拟的失效样本，直接利用 失效样本的权重指标求解失效概率，无需重新执行蒙特卡洛模拟，计算过程非常简单高效， 计算一种敏感性分析方案对应的边坡失效概率所需时间仅为5. 4X1(T3s。
[0146] 步骤5,根据敏感性分析方案的边坡失效概率，绘制失效概率随不确定性参数的统 计特征变化的趋势曲线，识别影响詹姆斯湾堤坝稳定可靠度的关键不确定性参数。
[0147] 图5~7分别为詹姆斯湾堤坝边坡失效概率随不确定性参数Tra、SjPDTill均值 变化的趋势曲线。如图5所示，随着硬质粘土厚度均值从2. 56m逐步增加到5. 44m、失 效概率从8X1(T3缓慢下降到1X1(T3左右。如图6所示，当湖积粘土粘聚力S&均值的从 12. 27kN/m2增加50. 13kN/m2时，詹姆斯湾堤坝失效概率则从KT1左右急剧下降到KT7,说 明湖积粘土粘聚力S&的均值的增加会导致詹姆斯湾堤坝可靠度水平急剧降低。如图7所 示，随着冰碛层厚度Dmi均值从15. 5m逐步增加21. 5m，詹姆斯湾堤坝的失效概率从1(T5增 加到1〇_2左右，可见詹姆斯湾堤坝的失效概率会随着不确定性参数冰碛层厚度DTill均值的 增加而增加。
[0148] 通过上述分析可以得出，在本实施例所研宄三个参数之中，詹姆斯湾堤坝的失效 概率对于冰碛层厚度Dmi、湖积粘土粘聚力S&较敏感，其中，湖积粘土粘聚力S&均值的变 化对詹姆斯湾堤坝稳定可靠度影响最大。故在工程上，应该采取相应的工程技术获得更多 关于湖积粘土粘聚力S&的信息，从而更加准确的评估边坡的安全性。
[0149] 上述实施例表明，本发明所提出的蒙特卡洛模拟的边坡稳定可靠度敏感性分析方 法，能够高效地求解敏感性分析方案的失效概率，有效地揭示边坡可靠度水平与不确定参 数的统计特征之间的响应关系，并识别出影响边坡稳定可靠度的关键不确定性参数，为边 坡风险控制、稳定加固，优化设计等提供了参考依据。
【主权项】
1. 基于蒙特卡洛模拟的边坡稳定可靠度敏感性分析方法，其特征是，包括步骤： 步骤1，确定不确定性参数的统计特征，并构造不确定性参数的联合概率密度函数 f(x)； 步骤2,基于联合概率密度函数f (X)，利用蒙特卡洛模拟方法获得边坡失效概率Pf，并 获得失效样本； 步骤3,设计多种敏感性分析方案，根据敏感性分析方案对应的统计特征，分别构造各 敏感性分析方案下不确定性参数的联合概率密度函数fk(x); 所述的敏感性分析方案采用如下方法获得：基于步骤1确定的统计特征，考虑不确定 参数的实际变化范围，改变一个或多个不确定性参数的统计特征值，保持其它不确定性参 数的统计特征概率分布不变，即获得一种敏感性分析方案； 步骤4,获得各敏感性分析方案下的边坡失效概率，本步骤进一步包括： 4. 1根据联合概率密度函数fk(x)和失效样本，获得失效样本的权重指标<、f(Xp和fk(Xp分别为第j个失效样本的权重指标、不确定性参数联合概率密度函数 和第k个敏感性分析方案下的联合概率密度函数； 4. 2根据权重指标获得敏感性分析方案的边坡失效概率· β/为第k个敏 感性分析方案的边坡失效概率，\为失效样本数，N为步骤2蒙特卡洛模拟方法中产生的随 机向量样本数； 步骤5,根据各敏感性分析方案的边坡失效概率，获得边坡失效概率随不确定性参数统 计特征的变化趋势，从而识别影响边坡稳定可靠度的关键不确定性参数。2. 如权利要求1所述的基于蒙特卡洛模拟的边坡稳定可靠度敏感性分析方法，其特征 是： 步骤1进一步包括子步骤： 1. 1根据不确定性参数的试验数据，或现有文献中不确定性参数的统计特征取值，确定 不确定性参数的统计特征以及各不确定性参数间的相关系数； 1. 2根据不确定性参数的统计特征及各不确定性参数间的相关系数，构造不确定性参 数的联合概率密度函数f(x)。3. 如权利要求2所述的基于蒙特卡洛模拟的边坡稳定可靠度敏感性分析方法，其特征 是： 子步骤1. 2具体为： (1) 若不确定性参数全部相互独立，则不确定性参数的联合概率密度函数为各不确定 性参数概率密度函数的乘积； (2) 若不确定性参数全部相关，则基于Copula理论或Nataf变换法构造不确定性参数 的联合概率密度函数； (3) 若不确定性参数部分独立、部分相关，首先，采用步骤（1)中方法获得相互独立的 不确定性参数的联合概率密度函数，记为其中，q为相互独立的不确定性参数个 数；然后，采用步骤（2)中方法获得相关的不确定参数的联合概率密度函数f(Xq+1，…，xn); 则不确定性参数的联合概率密度函数.
【专利摘要】本发明提供了一种基于蒙特卡洛模拟的边坡稳定可靠度敏感性分析方法，包括：步骤1，构造不确定性参数的联合概率密度函数；步骤2，利用蒙特卡洛模拟方法获得边坡失效概率，并获得失效样本；步骤3，设计多种敏感性分析方案，并分别构造各敏感性分析方案下不确定性参数的联合概率密度函数；步骤4，获得各敏感性分析方案下的边坡失效概率；步骤5，根据各敏感性分析方案的边坡失效概率，获得边坡失效概率随不确定性参数的统计特征的变化趋势。本发明适用范围广、计算过程简单、计算效率高，可有效揭示边坡可靠度水平与不确定性参数的统计特征间的响应规律，对边坡风险控制、设计优化等具有一定的指导意义。
【IPC分类】G06F17/50
【公开号】CN104899380
【申请号】CN201510319305
【发明人】李典庆, 吕天健, 张浮平, 曹子君, 唐小松, 周创兵
【申请人】武汉大学
【公开日】2015年9月9日
【申请日】2015年6月11日
【技术领域】
[0001] 本发明涉及一种边坡参数敏感性分析方法，尤其涉及一种基于蒙特卡洛模拟的边 坡稳定可靠度敏感性分析方法。
【背景技术】
[0002] 基于可靠度理论的边坡敏感性分析能有效地表征输入参数的不确定性及其统计 特征（均值、方差等）对边坡稳定可靠度的影响，对边坡设计、加固等具有重要指导意义，近 年来受到广泛关注。蒙特卡洛模拟方法作为一种求解可靠度的随机模拟方法，因其概念简 单、适用性强在岩土工程可靠度领域得到广泛应用。然而蒙特卡洛模拟方法存在不确定性 传播机理不明确的缺点，无法直接获得失效概率随着输入参数不确定性变化的响应规律。 因此，在利用蒙特卡洛模拟进行不确定性参数的敏感性分析时，通常需要重复模拟，计算输 入随机变量统计特征的一系列尚散点对应的失效概率，然后分析统计特征与失效概率之间 的变化规律。在工程实践中，一般边坡设计方案对应的失效概率较小（比如，低于0.001)， 而对于小概率问题，蒙特卡洛模拟方法计算效率较低，重复模拟需要耗费大量的时间和计 算机资源，所需的计算时间和资源随着边坡稳定性分析模型的复杂程度的增加而增加。
[0003]Wang[1]等提出了一种基于蒙特卡洛模拟的敏感性分析方法，该方法能够有效地 定量分析不确定性参数对于失效概率的相对贡献以及计算不同不确定性参数值对应的条 件失效概率，但是未对失效概率与不确定性参数的统计特征间的响应规律进行研宄。随后 Wang[2]提出了一种结合贝叶斯理论的敏感性分析方法，该方法不需要重新执行蒙特卡洛模 拟，就可以计算单个不确定性参数的统计特征变化后边坡的失效概率，获得失效概率与不 确定性参数的统计特征之间的响应规律。但是，该方法需要制作频率直方图，计算过程较为 复杂，且无法解决多个变量统计特征同时变化的敏感性分析问题。
[0004] 文中涉及如下参考文献：
[0005] [l]ffangY,CaoZJ,AuSK.EfficientMonteCarloSimulationofparameter sensitivityinprobabilisticslopestabilityanalysis[J].Computersand Geotechnics, 2010, 37(7 - 8):1015 - 1022.
[0006] [2]WangY.UncertainparametersensitivityinMonteCarloSimulationby samplereassembling[J].ComputersandGeotechnics, 2012, 46:39 - 47.
[0007] [3]KiureghianAD,LiuPL.Structuralreliabilityunderincomplete probabilityinformation[J].JournalofEngineeringMechanics, 1986, 112:86- 114.
【发明内容】
[0008] 针对现有技术存在的不足，本发明基于蒙特卡洛模拟提出了一种计算过程简便、 计算效率高的边坡稳定可靠度敏感性分析方法。
[0009] 为解决上述技术问题，本发明采用如下的技术方案：
[0010] 基于蒙特卡洛模拟的边坡稳定可靠度敏感性分析方法，包括步骤：
[0011] 步骤1，确定不确定性参数的统计特征，并构造不确定性参数的联合概率密度函数 f(x)；
[0012]步骤2,基于联合概率密度函数f(x)，利用蒙特卡洛模拟方法获得边坡失效概率 Pf，并获得失效样本；
[0013] 步骤3,设计多种敏感性分析方案，根据敏感性分析方案对应的统计特征，分别构 造各敏感性分析方案下不确定性参数的联合概率密度函数fk(x);
[0014] 所述的敏感性分析方案采用如下方法获得：基于步骤1确定的统计特征，考虑不 确定参数的实际变化范围，改变一个或多个不确定性参数的统计特征值，保持其它不确定 性参数的统计特征概率分布不变，即获得一种敏感性分析方案；
[0015] 步骤4,获得各敏感性分析方案下的边坡失效概率，本步骤进一步包括：
[0016] 4. 1根据联合概率密度函数fk(x)和失效样本，获得失效样本的权重指标
f(\)和fk(\)分别为第j个失效样本的权重指标、不确定性参数联合概 率密度函数和第k个敏感性分析方案下的联合概率密度函数；
[0017] 4. 2根据权重指标获得敏感性分析方案的边坡失效概率
，P/为第k 个敏感性分析方案的边坡失效概率，^为失效样本数，N为步骤2蒙特卡洛模拟方法中产生 的随机向量样本数；
[0018] 步骤5,根据各敏感性分析方案的边坡失效概率，获得边坡失效概率随不确定性参 数统计特征的变化趋势，从而识别影响边坡稳定可靠度的关键不确定性参数。
[0019] 步骤1进一步包括子步骤：
[0020] 1.1根据不确定性参数的试验数据，或现有文献中不确定性参数的统计特征取值， 确定不确定性参数的统计特征以及各不确定性参数间的相关系数；
[0021] 1. 2根据不确定性参数的统计特征及各不确定性参数间的相关系数，构造不确定 性参数的联合概率密度函数f(x)。
[0022] 子步骤1. 2具体为：
[0023] (1)若不确定性参数全部相互独立，则不确定性参数的联合概率密度函数为各不 确定性参数概率密度函数的乘积；
[0024] (2)若不确定性参数全部相关，则基于Copula理论或Nataf变换法构造不确定性 参数的联合概率密度函数；
[0025] (3)若不确定性参数部分独立、部分相关，首先，采用步骤（1)中方法获得相互独 立的不确定性参数的联合概率密度函数，记为
，其中，q为相互独立的不确定性参 数个数；然后，采用步骤（2)中方法获得相关的不确定参数的联合概率密度函数f(xq+1，一 ，xn);则不确定性参数的联合概率密度函数
[0026] 步骤2进一步包括子步骤：
[0027] 2. 1产生服从于联合概率密度函数f(x)的随机向量样本；
[0028] 2. 2对随机向量样本进行边坡稳定性分析，判断各随机向量样本对应的边坡是否 失效；
[0029] 2. 3统计失效样本数，获得边坡失效概率
其中，^为失效样本数，N为随机 向量样本数；所述的失效样本即失效边坡对应的随机向量样本；
[0030] 2. 4输出并保存失效样本。
[0031] 与现有技术相比，本发明具有以下优点和有益效果：
[0032] 1、适用范围广：
[0033] 本发明基于蒙特卡洛模拟，概念清晰、简单易懂。此外，本发明能够有效地解决单 一不确定性参数和多种不确定性参数的边坡稳定可靠度敏感性分析问题，同时也能较好地 适用于功能函数表达式为隐式的复杂边坡稳定可靠度敏感性分析问题，适用范围广。
[0034] 2、计算过程简单，计算效率高：
[0035] 传统蒙特卡洛模拟方法在计算敏感性分析方案所对应的边坡失效概率时，需重复 执行蒙特卡洛模拟，计算效率较低，尤其对于复杂边坡的稳定可靠度敏感性分析而言，蒙特 卡洛模拟中边坡确定性分析所需时间长，重复执行蒙特卡洛模拟进行敏感性分析计算时间 长，对计算资源要求高。本发明充分利用不确定性参数原概率分布蒙特卡洛模拟的计算信 息，在求解敏感性分析方案失效概率过程中，只需计算失效样本的权重指标，避免了蒙特卡 洛模拟的重复执行，简化了计算过程，提高了计算效率，是一种简单高效的边坡稳定可靠度 敏感性分析方法。尤其对于上述复杂边坡的稳定可靠度敏感性分析而言，本发明对计算效 率的提高显得尤为突出，有利于推动蒙特卡洛模拟方法在边坡可靠度工程的应用。
[0036] 3、揭示不确定性传播机理，为边坡工程设计提供重要的参考依据：
[0037] 本发明能有效地揭示边坡可靠度水平与不确定性参数的统计特征间的响应关系， 从而明确了蒙特卡洛模拟过程中不确定性的传播机理，对边坡风险控制、设计优化等具有 一定的指导意义。
[0038] 综上，本发明可简化计 算过程，提高计算效率，揭示蒙特卡洛模拟过程中不确定性 的传播机理，对推动岩土工程可靠度与风险理论在工程实践中应用有着重要的实用价值。
【附图说明】
[0039] 图1为本发明的具体流程示意图；
[0040] 图2为基于蒙特卡洛模拟的边坡可靠度分析流程图；
[0041] 图3为实施例中边坡的剖面图；
[0042] 图4为实施例中不同敏感性方案的不确定性参数的概率分布图，以硬质粘土厚度 T"为例；
[0043] 图5~7为实施例中敏感性分析结果。
【具体实施方式】
[0044] 下面将结合附图和【具体实施方式】进一步说明本发明。
[0045] 本发明的具体流程见图1，具体步骤如下：
[0046] 步骤1，确定不确定性参数的统计特征，并构造不确定性参数的联合概率密度函数 f(X)，即不确定性参数原分布的联合概率密度函数f(X)，X为表征不确定性参数的随机向 量。
[0047] 岩土工程可靠度分析中，土体参数可分为两类参数：确定性参数和不确定性参数， 不确定性参数包括土体粘聚力、内摩擦角、土体容重、土层厚度等。对于不确定性参数，其不 确定性的量化可根据有无试验数据分为两种情况：一，若边坡设计中有现场或室内试验结 果，可根据不确定性参数的试验数据，采取统计方法量化其不确定性。比如对于不确定性参 数土体重度而言，在现场取多个测点试样，根据各测点试样的土体重度，利用统计方法确定 其统计特征，例如均值、方差等，然后根据获得的统计特征量化土体重度的不确定性。二，若 没有关于不确定性参数的试验资料，可根据现有文献不确定性参数的统计特征取值，对所 研宄边坡中相应不确定性参数的统计特征进行合理的设定。
[0048] 本步骤进一步包括以下子步骤：
[0049] 1. 1根据不确定性参数的试验数据，或现有文献不确定性参数的统计特征取值，确 定不确定性参数的统计特征以及各不确定性参数间的相关系数P。
[0050] 1. 2根据不确定性参数的统计特征以及各不确定性参数间的相关系数，构造不确 定性参数的联合概率密度函数f(x)。
[0051] 联合概率密度函数f(X)的构造可分为三种情况：（1)不确定性参数全部相互独立 情况下联合概率密度函数的构造；(2)不确定性参数全部相关情况下联合概率密度函数的 构造；（3)不确定性参数部分独立、部分相关的情况下联合概率密度函数的构造。
[0052] 对于相互独立的不确定性参数而言，n维不确定性参数的联合概率密度函数f(x) 为各不确定性参数的概率密度函数的乘积，因此f(x)可以表示为：
[0054] 式⑴中：f(Xi)为第i个不确定性参数\的概率密度函数，其中取值， n为不确定性参数个数，n表示连乘运算。
[0055] 对于相关的不确定性参数而言，联合概率密度函数的构造可采用Copula理论、Nataf?变换法等。本具体实施中采用较常用的Nataf?变换法构造相关的不确定性参数的联 合概率密度函数。现以二维相关不确定性参数XJPX2为例，介绍Nataf变换法的基本原 理。
[0056] 将相关的不确定性参数&和X2的概率密度函数分别记为f(xJ和f(x2)，其中，Xl 和&分别表示X:和X2的取值，其累计分布函数分别记为Fi(Xl)和F2 (x2)，不确定性参数& 和&间的皮尔逊相关系数为P。对不确定性参数X:和X2进行如下等概率变换：
[0058] 式（2)中(yi)为标准正态变量^的累计分布函数，〇<( ?)为标准正态变量累 计分布函数〇 ( ?)的反函数，1为X,等概率变换后对应的标准正态分布变量，x,和y1表 示XJPY^取值。
[0059] 根据Nataf变换法，可推导不确定性参数XJPX2的联合概率密度函数f(x)为：
[0061] 式⑶中：P。为变量YjpY2的皮尔逊相关系数；<i> ( ?)为标准正态分布的概率 密度函数；巾2 (?)为二维标准正态分布联合概率密度函数，其表达式为：
[0063] 根据相关系数定义，P和口^具有关系：
[0065] 式（5)中：yjpy2分别为不确定性参数XjPX^均值；〇 〇 2分别为不确 定性参数XJPX2的标准差。
[0066]变量YjPY2的相关系数P。可通过式（5)迭代求解，也可以根据DerKiureghian 和Liu[3]提供的针对不同概率统计特征的经验公式求解，还可以根据高斯-埃尔米特积分 求解。但是一般情况下，pjPP相差不大，在近似求解中可以直接取PP。
[0067] 需要指出的是，岩土体不确定性参数通常服从正态分布。若不确定性参数XJPX^ 服从正态分布，此时根据式（3)和（4)，不确定性参数XJPX2的联合概率密度函数f(x)可 以表达为：
[0069] 对于部分独立、部分相关的不确定性参数而言，根据公式（1)获得独立的不确定 性参数的联合概率密度函数，记为1
?，根据上述Nataf变换方法或Copula理论可获 得相关的不确定参数的联合概率密度函数f(xq+1，…，xn)，故部分独立部分相关情况下不确 定性参数的联合概率密度函数：
，其中，f(Xi)表示独立不确定 性参数Xi,(i= 1，. . .q)的概率密度函数，Xq+pX#、…Xn表示相关的不确定性参数。
[0070] 步骤2,根据步骤1构建的联合概率密度函数f(x)，利用蒙特卡洛模拟方法计算边 坡失效概率Pf，并输出和保存失效样本j= 1，2, . . .，ns。
[0071] 本步骤采用传统蒙特卡洛模拟方法实现，具体流程见图2,主要包括如下子步骤：
[0072] 2. 1根据步骤1构建的联合概率密度函数f(x)，产生服从于联合概率密度函数 f(x)的随机向量样本Xpi= 1，2,...N，i表示随机向量样本编号。
[0073] 2. 2构建边坡稳定性分析模型，计算随机向量样本xj#应的边坡功能函数值g(xJ =FS(Xi)-l，并根据功能函数值g(Xi)判断边坡是否失效。
[0074] 目前常用的边坡稳定性分析法有极限平衡法和有限元法，其中极限平衡法包括 瑞典条分法、简化毕肖普法、摩根斯坦-普赖斯法等。用于边坡稳定性的商业软件主要有 GeoStudio(极限平衡法）、ANSYS(有限元法），ABAQUS(有限元法）等。本发明适用任何一 种边坡稳定性分析法。
[0075] 对于具体的边坡工程，根据要采用的边坡稳定性分析方法，运用相应的商业软件 建立边坡稳定性分析模型。随后，代入随机向量Xi，计算其对应的边坡功能函数值g(xj= FS(xj-1,当边坡功能函数值g(xj>0,即边坡安全系数FS(xD小于1时，表示边坡处于稳定 状态；当g(xD= 0时，表示边坡处于临界状态；当g(x)〈0时，表示边坡失稳，S卩边坡失效。
[0076] 循环执行子步骤2. 1和子步骤2. 2,直至i=N。
[0077] 2. 3统计失效样本数目ns，计算边坡失效概率
[0078] 满足g(Xi)〈0的随机向量样本Xi为失效样本。统计失效样本数量计为ns，失效样 本记为j= 1，2,. . .，ns。根据蒙特卡洛模拟方法原理，边坡失效概率Pf可以表示为：
[0080] 式（7)中，巧表示边坡失效概率估值，I[g(x)]为指示函数，当g(x)〈0时， I[g (X)]取 1;当g (x)彡 0 时，I[g (x)]取为 0。
[0081] 2. 4输出并保存失效样本Xj，j= 1，2, ? ? ?，ns。
[0082] 步骤3,基于步骤1确定的不确定性参数统计特征，通过增加和减小不确定性参数 统计特征获得一系列敏感性分析方案，根据各敏感性分析方案对应的不确定性参数的统计 特征构造联合概率密度函数。
[0083] 针对所研宄的不确定性参数，通过改变其统计特征值，在确保其不超出相应不确 定性参数的实际变化范围前提下，设置M种敏感性分析方案；根据各敏感性分析方案对应 的统计特征，分别构造相应的不确定性参数联合概率密度函数fk (x)，k= 1，2,...M，M为敏 感性分析方案数量，k表示 敏感性分析方案编号，fk(x)表示第k个敏感性分析方案相应的 不确定性参数的联合概率密度函数。
[0084] 利用蒙特卡洛模拟方法分析不确定性参数敏感性时，需要通过计算输入不确定性 参数统计特征（如均值，方差）的一系列离散点的边坡失效概率，以探宄边坡失效概率和不 确定性参数的统计特征间的响应规律。
[0085] 如对土体粘聚力参数进行敏感性分析时，可通过增加和减小粘聚力的均值y或 标准差〇，设置M种粘聚力的不同分布作为敏感性分析方案。然后，根据各敏感性分析方案 下的统计特征，构造各敏感性分析方案下的不确定性参数联合概率密度函数fk(x)。联合概 率密度函数fk(x)的具体构造方法同子步骤1.2。
[0086] 这里需要指出的是，本发明能够对多个不确定性参数同时进行敏感性分析，如可 以设置一系列同时改变多个不确定性参数的均值y和标准差〇的敏感性分析方案。
[0087] 步骤4,利用步骤3获得的联合概率密度函数fk(x)和步骤2获得的失效样本Xj， j= 1，2,. . .，ns，求解各敏感性分析方案下的边坡失效概率Pfk，根据边坡失效概率Pfk随不 确定性参数统计特征的变化规律对不确定性参数的统计特征和变异系数进行敏感性分析。
[0088] 根据边坡失效概率通用公式，第k个敏感性分析方案下的边坡失效概率Pfk可表示 为：
[0090] 式⑶中为在联合概率密度函数fk(x)下的数学期望算子。
[0091] 结合不确定性参数原分布的联合概率密度函数f(x)，即步骤1构造的联合概率密 度函数f(x)，第k个敏感性分析方案下的边坡失效概率Pfk可进一步表示为：
[0093] 式（9)中，《k为本发明自定义的权重指标，其计算表达式如下：
[0095] 通过权重指标《k可以把求解敏感性分析方案条件下的边坡失效概率Pfk转化为 求解I[g(x) ] ?15在f(x)分布下的数学期望。因此，可以直接利用步骤2. 1中蒙特卡洛模拟 的随机样本Xi，i=l，2, ? ? *，N，计算边坡失效概率Pfk。故基于联合概率密度函数f(x) 产生的N组随机样本求解第k个敏感性分析方案的边坡失效概率Pfk可表示为：
[0097] 考虑到指示函数I[g(x)]的取值特征，仅失效样本Xj,j= 1，2，? ? ?，ns会对失 效概率Pfk产生贡献，公式（11)可进一步表达为：
[0099] Pfk估计值的方差]可表示为：
[0101] 的变异系数COV[$]可表示为：
[0103] #的变异系数cov[#]:能够有效地表征计算结果的准确性，为本发明方法计算 结果的可靠性提供依据。
[0104] 由式（12)可以看出，在求解敏感性分析方案的边坡失效概率时，本发明只需要计 算原分布产生的失效样本处的权重指标并对其求和，无需重新进行蒙特卡洛模拟，计算过 程较为简单、计算效率高。
[0105] 需要说明的是，本步骤中计算失效概率所用联合概率密度函数为失效样本范围内 的联合概率密度函数，有别于步骤1到3中全体样本范围内的联合概率密度函数。
[0106] 本步骤的具体子步骤如下：
[0107] 本步骤的具体子步骤如下：
[0108] 4. 1根据失效样本的联合概率密度函数f(Xj)和fk (Xj)，获得第j个失效样本Xj的 权重指标
，勿、f(xp和fk(xp分别为第j个失效样本的权重指标、不确定性 参数联合概率密度函数和第k个敏感性分析方案下的联合概率密度函数。
[0109] 4. 2根据权重指标 < 计算敏感性分析方案的边坡失效概率
[0110] 步骤5,根据敏感性分析方案的边坡失效概率，绘制边坡失效概率随不确定性参数 统计特征变化的趋势曲线，识别影响边坡稳定可靠度的关键不确定性参数。
[0111] 实施例1
[0112] 一、工程简介
[0113] 詹姆斯湾堤坝位于加拿大魁北克地区的一项水利工程，堤坝长50km，其中一种设 计方案剖面如图3所示，堤坝高12m，中间有一个56m宽平台，边坡坡角为18.4°，坡比为 3:1。坝体以下土层从上至下依次为硬质粘土、海成粘土、湖积粘土和冰碛层。
[0114] 参考文献[1]，本实施例考虑6个随机变量，即待研宄的不确定性参数，分别为硬 质粘土厚度、海成粘土不排水抗剪强度、湖积粘土的不排水抗剪强度S&、坝基到冰碛 层的距离Dmi、筑坝填土的摩擦角<i>Fill和重度YFill。不确定性参数均服从正态分布且相 互独立，相应的统计特征见表1。其它确定性参数如下：海成粘土的厚度为8m，硬质粘土 不排水抗剪强度为41kPa，硬质粘土、海成粘土、湖积粘土的重度分别为19kN/m3、19kN/m3、 20. 5kN/m3。此外，由几何关系可知湖积粘土厚度IY=Dmi-Tra-8。
[0115] 为简化计算，文献[1]对最危险滑动面增加两个限制：（1)最危险滑动面与冰碛层 相切并过点（X，y)，其中，x= 4. 9m，y= 36m; (2)最危险滑动面圆心横坐标x= 85. 9m。这 种处理方法使得最危险滑动面的位置仅与坝基到冰碛层的距离Dmi有关。
[0116] 表1不确定性参数的统计特征
[0117]
[0118] 二、具体实施流程
[0119] 步骤1，根据表1中各不确定性参数的统计特征，构造不确定性参数原分布的联合 概率密度函数f(ci)Fill，丫Fill, Tcr，SuL，Dml)〇
[0120] 由于本实施例中各不确定性参数服从正态分布且相互独立，因此联合概率密度函 数MFill，TFin?TcrJS^,SuLJDTill )可以表示为：
[0121] f( ^Fiii)TFill)Tcr,SuM,SuL,DTill) -f( 4) Fill)f(yFill)f(Tcr)f(S^)f(SuL)f(DTill) (15)
[0122] 其中：
[0124] 式中，yp y2、y3、y4、y5、y6和o工、o 2、o 3、o 4、o 5、o6分别为参数（J)Fm、 YFill、T。,、SfSupDm^均值和标准差，其具体数值见表1。
[0125] 步骤2,根据不确定性参数原分布的联合概率密度函数f(<i)Fill，yFill，Tra，S^，SuU Dmi)，利用蒙特卡洛模拟法获得詹姆斯湾堤坝的边坡失效概率Pf，并输出和保存失效样本 Xj，j l，2，...，ris〇
[0126]2. 1利用VBA程序产生服从于联合概率密度函数f(d>Fill， 的随机向量样本Xi，i= 1，2,. . .N，i表示随机向量样本编号，N= 2X105。
[0127] 2. 2构建边坡稳定性分析模型，计算随机向量样本xj#应的边坡功能函数值g(xJ =FS(Xi)_l，并根据边坡功能函数值g(Xi)判断詹姆斯湾堤坝是否失效。
[0128] 首先，采用简化毕肖普法，利用EXCEL表单构建詹姆斯湾堤坝稳定性分析模型。然 后，将随机向量样本作为输入参数带入詹姆斯堤湾坝的稳定性分析模型中。最后，根据稳定 性分析模型计算随机向量变量Xi对应的边坡安全系数FS(xJ和边坡功能函数值g(Xi)。
[0129] 循环执行子步骤2. 1和子步骤2. 2,直至i=N。
[0130] 2. 3统计失效样本数\，共有失效样本\= 452个。根据式（7)计算詹姆斯湾堤 坝的边坡失效概率Pf= 452/200000 = 2. 26X10 '
[0131] 2. 4将失效样本&储存在EXCEL表单内，供敏感性分析使用。
[0132] 以上基于不确定性参数原分布的蒙特卡洛模拟法（20万次）计算时间约10小时， 所采用计算机配置为内存4GB、CPUIntelCorei3和主频3. 3GHz。因此，如果采用传统的 蒙特卡洛模拟法求解詹姆斯湾堤坝的敏感性分析问题，对于每种敏感性分析方案，在相同 的模拟次数下所需要的计算时间都约为10个小时。由此可见，采用传统蒙特卡洛模拟法求 解边坡稳定可靠度敏感性分析问题，计算效率非常低。
[0133] 步骤3,通过改变所研宄 不确定性参数的统计特征值，设置一系列敏感性分析方 案。根据各敏感性分析方案对应的统计特征，构造各敏感性分析方案下不确定性参数的联 合概率浩度函数f(巾Fiii，yFiii，Tcr,S-SuL,DTill)，f(<}> Fill,yFill,Tcr,S-Si,DTill)表不弟k 种敏感性分析方案下的不确定性参数联合概率密度函数。
[0134] 本步骤通过改变某不确定性参数的均值、方差或其他统计特征，来获得多种敏感 性分析方案，并通过分析各敏感性分析方案的边坡失效概率，获得该不确定性参数的统计 特征对边坡稳定可靠度的影响程度，从而得到影响边坡稳定可靠度的关键不确定性参数。
[0135] 由于本实施例中不确定性参数较多，此处仅选取硬质粘土厚度、冰碛层厚度 Dmi和湖积粘土的不排水抗剪强度S&，研宄的均值变化对詹姆斯湾堤坝边坡失 效概率的影响。如表2所示，本实施例中通过改变不确定性参数UmpSi的均值，共设置 18种敏感性分析方案：方案1~6中，硬质粘土厚度的均值分别设为u3-3 〇 3、u3-2 〇 3、 u3- 〇 3、u3+ 〇 3、u3+2 〇 3、u3+3 〇 3,u3、〇 3分别为T原分布的均值和标准差，其概率分布曲线 见图4;其他不确定性参数的统计特征值均保持原分布的统计特征值。类似的，分别对参数 Sul和Dmi做同样的处理，可以得到方案7~12和方案13~18。
[0136] 各敏感性分析方案的联合概率密度函数fk((i)mi，YFill，Tra，s^，SuUDmi)的构造 与案例分析步骤1类似，在此不再赘述。
[0137] 表2敏感性分析方案
[0138]
[0139] 步骤4,根据求得的联合概率密度函数fk(<i>Fill，YFill，T。,，S^，DTill)和步骤2中 输出的失效样本\，j=l，2, ? ? ?，&，求解各敏感性分析方案的边坡失效概率Pfk。
[0140] 根据式（10)，第k个敏感性分析方案对应的权重指标《k可表示为：
[0142] 如当k取1时，权重指标o1可以表示为
[0144] 失效样本的权重指标值带入式（12)，则可以计算出各敏感性分析方案所对应 的边坡失效概率。
[0145] 在此需要指出的是，权重指标的计算可以结合EXCEL内置的概率密度函数（如正 态分布概率密度函数NORM.DIST)，计算十分简便。同时，由步骤3~4可以看出，在求解敏 感性分析方案的过程中，由于充分利用了原概率分布蒙特卡洛模拟的失效样本，直接利用 失效样本的权重指标求解失效概率，无需重新执行蒙特卡洛模拟，计算过程非常简单高效， 计算一种敏感性分析方案对应的边坡失效概率所需时间仅为5. 4X1(T3s。
[0146] 步骤5,根据敏感性分析方案的边坡失效概率，绘制失效概率随不确定性参数的统 计特征变化的趋势曲线，识别影响詹姆斯湾堤坝稳定可靠度的关键不确定性参数。
[0147] 图5~7分别为詹姆斯湾堤坝边坡失效概率随不确定性参数Tra、SjPDTill均值 变化的趋势曲线。如图5所示，随着硬质粘土厚度均值从2. 56m逐步增加到5. 44m、失 效概率从8X1(T3缓慢下降到1X1(T3左右。如图6所示，当湖积粘土粘聚力S&均值的从 12. 27kN/m2增加50. 13kN/m2时，詹姆斯湾堤坝失效概率则从KT1左右急剧下降到KT7,说 明湖积粘土粘聚力S&的均值的增加会导致詹姆斯湾堤坝可靠度水平急剧降低。如图7所 示，随着冰碛层厚度Dmi均值从15. 5m逐步增加21. 5m，詹姆斯湾堤坝的失效概率从1(T5增 加到1〇_2左右，可见詹姆斯湾堤坝的失效概率会随着不确定性参数冰碛层厚度DTill均值的 增加而增加。
[0148] 通过上述分析可以得出，在本实施例所研宄三个参数之中，詹姆斯湾堤坝的失效 概率对于冰碛层厚度Dmi、湖积粘土粘聚力S&较敏感，其中，湖积粘土粘聚力S&均值的变 化对詹姆斯湾堤坝稳定可靠度影响最大。故在工程上，应该采取相应的工程技术获得更多 关于湖积粘土粘聚力S&的信息，从而更加准确的评估边坡的安全性。
[0149] 上述实施例表明，本发明所提出的蒙特卡洛模拟的边坡稳定可靠度敏感性分析方 法，能够高效地求解敏感性分析方案的失效概率，有效地揭示边坡可靠度水平与不确定参 数的统计特征之间的响应关系，并识别出影响边坡稳定可靠度的关键不确定性参数，为边 坡风险控制、稳定加固，优化设计等提供了参考依据。
【主权项】
1. 基于蒙特卡洛模拟的边坡稳定可靠度敏感性分析方法，其特征是，包括步骤： 步骤1，确定不确定性参数的统计特征，并构造不确定性参数的联合概率密度函数 f(x)； 步骤2,基于联合概率密度函数f (X)，利用蒙特卡洛模拟方法获得边坡失效概率Pf，并 获得失效样本； 步骤3,设计多种敏感性分析方案，根据敏感性分析方案对应的统计特征，分别构造各 敏感性分析方案下不确定性参数的联合概率密度函数fk(x); 所述的敏感性分析方案采用如下方法获得：基于步骤1确定的统计特征，考虑不确定 参数的实际变化范围，改变一个或多个不确定性参数的统计特征值，保持其它不确定性参 数的统计特征概率分布不变，即获得一种敏感性分析方案； 步骤4,获得各敏感性分析方案下的边坡失效概率，本步骤进一步包括： 4. 1根据联合概率密度函数fk(x)和失效样本，获得失效样本的权重指标<、f(Xp和fk(Xp分别为第j个失效样本的权重指标、不确定性参数联合概率密度函数 和第k个敏感性分析方案下的联合概率密度函数； 4. 2根据权重指标获得敏感性分析方案的边坡失效概率· β/为第k个敏 感性分析方案的边坡失效概率，\为失效样本数，N为步骤2蒙特卡洛模拟方法中产生的随 机向量样本数； 步骤5,根据各敏感性分析方案的边坡失效概率，获得边坡失效概率随不确定性参数统 计特征的变化趋势，从而识别影响边坡稳定可靠度的关键不确定性参数。2. 如权利要求1所述的基于蒙特卡洛模拟的边坡稳定可靠度敏感性分析方法，其特征 是： 步骤1进一步包括子步骤： 1. 1根据不确定性参数的试验数据，或现有文献中不确定性参数的统计特征取值，确定 不确定性参数的统计特征以及各不确定性参数间的相关系数； 1. 2根据不确定性参数的统计特征及各不确定性参数间的相关系数，构造不确定性参 数的联合概率密度函数f(x)。3. 如权利要求2所述的基于蒙特卡洛模拟的边坡稳定可靠度敏感性分析方法，其特征 是： 子步骤1. 2具体为： (1) 若不确定性参数全部相互独立，则不确定性参数的联合概率密度函数为各不确定 性参数概率密度函数的乘积； (2) 若不确定性参数全部相关，则基于Copula理论或Nataf变换法构造不确定性参数 的联合概率密度函数； (3) 若不确定性参数部分独立、部分相关，首先，采用步骤（1)中方法获得相互独立的 不确定性参数的联合概率密度函数，记为其中，q为相互独立的不确定性参数个 数；然后，采用步骤（2)中方法获得相关的不确定参数的联合概率密度函数f(Xq+1，…，xn); 则不确定性参数的联合概率密度函数.
【专利摘要】本发明提供了一种基于蒙特卡洛模拟的边坡稳定可靠度敏感性分析方法，包括：步骤1，构造不确定性参数的联合概率密度函数；步骤2，利用蒙特卡洛模拟方法获得边坡失效概率，并获得失效样本；步骤3，设计多种敏感性分析方案，并分别构造各敏感性分析方案下不确定性参数的联合概率密度函数；步骤4，获得各敏感性分析方案下的边坡失效概率；步骤5，根据各敏感性分析方案的边坡失效概率，获得边坡失效概率随不确定性参数的统计特征的变化趋势。本发明适用范围广、计算过程简单、计算效率高，可有效揭示边坡可靠度水平与不确定性参数的统计特征间的响应规律，对边坡风险控制、设计优化等具有一定的指导意义。
【IPC分类】G06F17/50
【公开号】CN104899380
【申请号】CN201510319305
【发明人】李典庆, 吕天健, 张浮平, 曹子君, 唐小松, 周创兵
【申请人】武汉大学
【公开日】2015年9月9日
【申请日】2015年6月11日
最新回复(0)