本文档属于类型a：一篇关于单一原创性研究的学术论文。以下是根据文档内容撰写的详细学术报告。

基于各向异性空间变异性的边坡失效概率评估研究学术报告

一、 主要作者、机构、发表信息

本研究由龙龙陈 (Longlong Chen)、张文昌 (Wengang Zhang) 等学者合作完成。第一作者及部分团队成员来自重庆大学 (Chongqing University) 土木工程学院、山地城镇建设与新技术教育部重点实验室，以及库区环境地质灾害防治国家地方联合工程研究中心。其中，第一作者龙龙陈在论文研究期间亦在意大利米兰理工大学 (Politecnico di Milano) 土木与环境工程系进行相关工作。通讯作者为重庆大学的张文昌教授（邮箱：zhangwg@cqu.edu.cn）。

本研究以题为《Probabilistic assessment of slope failure considering anisotropic spatial variability of soil properties》（考虑土体性质各向异性空间变异性的边坡失效概率评估）的论文形式，发表于爱思唯尔 (Elsevier) 旗下的期刊《Geoscience Frontiers》2022年第13卷。文章于2021年7月19日收稿，2021年12月30日修订，2022年2月12日接受，并于2022年2月16日在线发表。该文章遵循CC BY-NC-ND 4.0知识共享许可协议开放获取。

二、 研究背景与目的

科学领域与研究背景： 本研究属于岩土工程可靠性分析与边坡稳定性评估领域。在传统确定性边坡稳定性分析中，土体的抗剪强度参数（如黏聚力c和内摩擦角φ）通常被视为确定值。然而，由于沉积、风化、填充等复杂自然过程，即使在同一土层内，土体性质也表现出显著的空间变异性 (Spatial Variability)。这种变异性本质上是空间相关的，即相邻土体性质更为相似。忽略这种空间相关性，而简单地将参数视为完全独立的随机变量，可能会导致对边坡失效概率的高估，从而无法真实反映工程风险。

现有研究局限与问题： 尽管随机场理论已被引入来刻画土体的空间变异性，但现有研究大多集中于各向同性 (Isotropic) 或横向各向异性 (Transverse Anisotropic) 的空间变异模式。这两种模式假设土体性质波动的主方向（即波动尺度 (Scale of Fluctuation) 最大的方向）要么处处相同（各向同性），要么恰好与水平或垂直坐标轴对齐（横向各向异性）。然而，实际地质环境（如受断层影响的沉积、崩积过程、倾斜地层中的风化土等）常常导致土体性质呈现出更复杂的一般各向异性空间变异性 (General Anisotropic Spatial Variability)，其波动尺度的两个主方向既不相等，也不一定与坐标轴平行。此外，现有研究主要关注失效概率本身，很少系统量化空间变异性对边坡失效模式（浅层、中层、深层或多滑面破坏）及其滑体体积的影响，而这些信息对于边坡加固设计和风险控制至关重要。

研究目的： 因此，本研究旨在：1）建立一个综合概率评估框架，以量化一般各向异性空间变异性对边坡失效概率、失效模式和滑体体积的影响；2）揭示一般各向异性空间变异性模式相较于传统各向同性及横向各向异性模式对边坡可靠性的独特影响机制；3）通过参数分析，探究土体参数间互相关性 (Cross-correlation)、波动尺度等关键因素在一般各向异性空间变异条件下对边坡可靠性的敏感性。

三、 详细研究流程与方法

本研究采用了一个集成的数值模拟框架——随机有限差分法 (Random Finite Difference Method， RFDM)，并结合蒙特卡洛模拟 (Monte Carlo Simulation， MCS) 进行概率评估。整个工作流程包含以下核心步骤：

1. 建立确定性边坡模型与随机场参数设定：

   • 研究对象： 研究选取了一个经典的有黏性-摩擦性土质边坡作为算例。边坡高10米，坡比为1:1（垂直:水平），基础厚度5米，总长30米。模型通过FLAC3D软件建立，共离散为4820个单元（主要为0.25米尺寸的六面体单元）。
   • 土体参数： 将黏聚力c和内摩擦角φ的正切值 (tan φ) 建模为服从对数正态分布 (Lognormal Distribution) 的互相关随机变量。均值分别为10 kPa和tan(30°)，变异系数分别为0.3和0.2。其他参数（如弹性模量、泊松比、重度）视为常数。
   • 空间变异性模式与表征： 研究重点对比了三种空间变异模式：各向同性、横向各向异性和一般各向异性。对于一般各向异性，其空间结构由两个主方向上的波动尺度（最大波动尺度d_max和最小波动尺度d_min）以及这两个主方向之间的夹角β来定义。研究采用了平方指数自相关函数 (Squared Exponential Auto-correlation Function) 来刻画空间相关性，因其能更好地满足各向同性条件下各方向波动尺度相等的数学要求。波动尺度d_h随方向角h的变化关系通过公式严格定义。

2. 生成各向异性随机场：

   • 方法： 采用矩阵分解法 (Matrix Decomposition Method)，或称中点法，来生成对数正态随机场的实现。该方法首先基于选定的自相关函数和模型各单元的中心坐标，构建n×n维的自相关矩阵Q_x（n为单元数）。然后通过Cholesky分解将Q_x分解为下三角矩阵L及其转置的乘积。接着，生成一个n×1维的独立标准正态随机向量Z，并通过U_d = L * Z得到相关标准正态随机场U_d。最后，通过等概率变换U = exp(μ_ln + σ_ln * U_d)将其转换为具有指定均值和标准差的对数正态随机场U。
   • 验证： 研究生成了一典型一般各向异性随机场（d_max = 4 m， d_min = 10 m， β = 45°），并验证了其统计特性（均值、标准差）以及波动尺度随方向角变化的曲线与理论公式高度吻合，证明了矩阵分解法的有效性。

3. 执行随机有限差分法 (RFDM) 与蒙特卡洛模拟 (MCS)：

   • 流程： 在FLAC3D中，将每个随机场实现（即一组在空间上变化的c和φ值）分配到对应的有限差分单元中。对每一个实现，采用强度折减法 (Shear Strength Reduction， SSR) 计算边坡的安全系数 (Factor of Safety， FoS)。强度折减法通过逐步折减土体的c和tan φ值，直至边坡达到临界破坏状态，此时的折减系数即为FoS。该方法无需预先假定滑面形状。
   • 概率计算： 重复进行大量（例如1000次）MCS，得到一系列FoS值。边坡的失效概率 (Probability of Failure, Pf) 计算为FoS < 1的实现次数占总模拟次数的比例。为确保计算精度与效率的平衡，研究设置了收敛准则：若模拟次数N<1000时，Pf的变异系数 (Coefficient of Variation， COV) 小于或等于0.1，则停止模拟；否则，模拟至1000次，并根据FoS的累积分布函数估算Pf。

4. 确定滑体体积与失效模式分类：

   • 滑体体积计算（创新方法）： 研究提出了一种基于剪切应变增量 (Shear Strain Increment， SSI) 聚类的算法来自动识别滑面并计算滑体体积。步骤包括：a) 从FLAC3D中提取临界破坏状态下的SSI分布；b) 利用K-means聚类算法自动识别高SSI区域，即潜在的剪切带 (Shear Band)；c) 构造穿过剪切带的虚拟竖直线，选取每条线上SSI最大的单元作为临界单元；d) 连接所有临界单元的中心点，形成临界滑面；e) 若剪切带不连续，则可能存在多个独立滑面（多滑面破坏），需分别识别；f) 计算位于最宽剪切带内、临界滑面之上的所有单元的体积之和，即为滑体体积。
   • 失效模式定义： 基于临界滑面出口点与坡趾点的水平相对位置，定义了四类失效模式：1) 浅层破坏：出口点在坡趾点后方（坡面上）；2) 中层破坏：出口点在坡趾点处；3) 深层破坏：出口点在坡趾点前方（坡脚地基出露）；4) 多滑面破坏：存在多个独立的剪切带和滑面。当中层破坏的出口点位于坡趾附近（0.05倍坡高范围内）时，也归类为中层破坏。

5. 参数化敏感性分析：

   • 分析场景： 针对一般各向异性空间变异性，研究细分为两种子场景进行深入分析：“水平场景”（d_max方向固定为水平，变化最小波动尺度的方向角h_{ga}^{d_min}）和“非水平场景”（d_min方向固定为水平，变化最大波动尺度的方向角h_{ga}^{d_max}）。
   • 分析参数： 在基础参数上，系统地变化以下参数，研究其对边坡失效概率、失效模式和滑体体积的影响：
     • 最小/最大波动尺度的方向角（h_{ga}^{d_min} 或 h_{ga}^{d_max}，从15°到165°）。
     • 黏聚力c与内摩擦角φ之间的互相关系数 (ρ_c,φ)（从-0.7到0.5）。
     • 最大波动尺度 (d_max)（从10米到60米）。
     • 最小波动尺度 (d_min)（从1米到6米）。

四、 主要研究结果

1. 空间变异性模式及方向角的显著影响：

   • 水平场景： 当d_max方向水平时，失效概率随最小波动尺度方向角h_{ga}^{d_min}变化呈现以h_{ga}^{d_min}=90°（即横向各向异性）为对称轴的近似对称分布，并在h_{ga}^{d_min}=90°处达到峰值（Pf≈8.1%）。这表明在水平层状边坡中，传统的横向各向异性假设可能导致最保守（失效概率最高）的评估结果。当h_{ga}^{d_min}偏离90°时，土壤强度的空间变异性增强，但不同随机场实现间整体平均强度的方差减小，导致FoS的变异性降低，从而失效概率下降。
   • 非水平场景： 当d_max方向不固定于水平时，其方向角h_{ga}^{d_max}对失效概率影响更为剧烈。当h_{ga}^{d_max}=150°（即土体强度层理以约30°倾角向坡外倾斜）时，失效概率达到全局最大值（Pf显著高于水平场景峰值）。这是因为这种“临空倾角 (Daylighting Orientation)”的弱层理最有利于滑体沿层面滑出，是滑坡发生的有利条件。相反，当层理反倾（如h_{ga}^{d_max}=15°）时，边坡可靠性更高。

2. 一般各向异性对参数敏感性的放大效应：

   • 与横向各向异性的对比： 研究将h_{ga}^{d_max}=150°的一般各向异性案例（最危险情况）与文献中仅考虑横向各向异性的结果进行了对比。
   • 互相关性 (ρ_c,φ) 的影响： ρ_c,φ从-0.7增至0.5时，一般各向异性条件下的Pf增长幅度（从约0.013增至0.283）远大于横向各向异性条件下的增长幅度（从约0.005增至0.12），放大倍数约为2-3倍。这是因为正相关性增加了土体总强度的方差，而在一般各向异性下，这种方差放大效应更为显著。
   • 波动尺度 (d_max, d_min) 的影响： 对于d_max，一般各向异性下的Pf是横向各向异性下的约1.3-2倍。对于d_min，放大效应更为惊人：当d_min从1米增至6米，一般各向异性下的Pf从0.016增至0.072，而横向各向异性下的Pf仅从约0.0003增至0.035，放大倍数可达2-50倍。这主要是因为低估d_min（或采用各向同性/横向各向异性模型）会高估土壤的空间变异性，从而可能得出过于保守或偏离实际的可靠性评估。

3. 失效模式与滑体体积的统计规律：

   • 主导失效模式： 在绝大多数考虑空间变异性的情况下（无论是哪种各向异性模式），边坡的主导失效模式均为中层破坏。这与均质状态下（确定性分析）发生的浅层破坏模式有显著区别。空间变异性促使滑面倾向于寻找最薄弱路径，导致破坏深度增加。
   • 滑体体积排序： 不同失效模式导致的平均滑体体积大小顺序为：深层破坏 > 多滑面破坏 ≈ 中层破坏 > 浅层破坏。
   • 关键参数对后果的影响： 在一般各向异性条件下，高估ρ_c,φ (正值方向)、d_max或d_min，会显著增加深层破坏和多滑面破坏模式的平均滑体体积。例如，ρ_c,φ从-0.7升至0.5，深层破坏的平均滑体体积从80.7 m³/m增至90.3 m³/m。d_max从10米增至60米，也会导致这两种破坏模式的滑体体积明显上升。这揭示了在不利的空间变异结构下，虽然深层和多滑面破坏的发生概率可能不高，但一旦发生，其后果（滑体体积）可能极其严重，为边坡风险管控提供了重要洞见。

五、 结论与研究意义

结论： 本研究成功建立并应用了考虑土体性质各向异性空间变异性的边坡失效概率RFDM评估框架。主要结论包括：1）土体空间变异性的模式（特别是最大波动尺度的方向）对边坡可靠性有决定性影响，其中最大波动尺度以约30°向坡外倾斜的一般各向异性模式是最危险的情况。2）与广泛采用的横向各向异性模型相比，一般各向异性模型会显著放大互相关性ρ_c,φ和波动尺度d_max、d_min对边坡失效概率的影响（放大倍数最高可达50倍）。3）空间变异性改变了边坡的破坏机制，使中层破坏成为主导模式。4）对ρ_c,φ、d_max、d_min的高估会急剧增加深层和多滑面破坏的潜在滑体体积，加重灾害后果。

科学价值与应用价值： * 科学价值： 本研究深化了对复杂各向异性空间变异性影响边坡失稳机理的认识，将可靠性分析从单一的失效概率计算拓展到对失效模式和灾害后果（滑体体积）的联合概率预测，为岩土工程概率风险评估理论提供了更完善的框架。 * 应用价值： 研究结果对工程实践具有重要指导意义：1）勘察与建模建议： 工程师应基于现场地质条件和勘探试验数据，审慎确定土体的空间变异模式及相关参数（如d_max、d_min的方向和大小，ρ_c,φ），避免简单地采用各向同性或横向各向异性假设，以防导致非保守或过度保守的设计。2）设计与加固启示： 分析表明，中层破坏是主要风险，且深层/多滑面破坏后果严重。这提示加固措施（如抗滑桩）应有足够的锚固深度以抵御中层滑面，并对潜在的深层软弱区给予特别关注。3）风险评估： 提出的自动化滑面识别与滑体体积计算方法，可为滑坡风险评估中的灾害范围与规模预测提供实用工具。

六、 研究亮点

1. 研究对象新颖： 首次系统、深入地研究了一般各向异性空间变异性这一在实际地质环境中常见但既往研究忽视的复杂模式对边坡可靠性的全面影响。
2. 方法集成创新： 将随机场理论、有限差分法（FLAC3D）、蒙特卡洛模拟和创新的基于剪切应变增量的滑面自动识别算法相结合，形成了一个自动化、高效的随机有限差分法 (RFDM) 概率评估框架。该框架不仅能计算失效概率，还能自动判别失效模式和精确计算滑体体积。
3. 研究维度拓展： 超越了传统可靠性分析仅关注失效概率的局限，开创性地对空间变异性导致的不同失效模式的发生概率及其对应的平均滑体体积进行了系统量化统计，实现了从“会不会破坏”到“以何种方式破坏、破坏后果多大”的更全面风险评估。
4. 揭示放大效应： 明确量化了在一般各向异性条件下，关键参数（ρ_c,φ， d_max， d_min）对可靠性的影响被显著放大的现象，这对正确理解和应用概率分析结果至关重要。

七、 其他有价值的讨论

研究还对样本选择策略进行了简要讨论。为了获得足够稳定的失效模式和滑体体积统计数据，作者使用了全部1000次MCS样本（包括稳定和不稳定样本）进行统计分析，而非仅使用失效样本。结果显示，使用全样本得到的统计结果（如失效模式比例）比仅用失效样本更为稳定。作者指出，采用子集模拟 (Subset Simulation) 等更高效的随机模拟方法，有望用更少的计算量获得更多失效样本，从而进一步提升对低概率破坏模式（如深层破坏）统计特性的估计精度。这一讨论为后续相关研究的方法选择提供了参考。