露天矿边坡稳定性并非静态条件，而是随开挖推进和地质力学条件变化而持续演化的动态状态。在本研究中，研究人员将地基雷达监测、卫星InSAR时序分析与数值稳定性建模相结合的集成方法应用于评估一处具有复杂地质和构造特征的大规模露天铜矿的边坡行为。雷达数据揭示了集中于特定边坡区的渐进式和间歇性变形，而InSAR观测表明在主变形阶段之后变形以较低速率持续，提供了边坡响应的长期视角。采用极限平衡法（Limit Equilibrium Method, LEM）和有限元法（Finite Element Method, FEM）的稳定性分析表明，边坡处于接近极限平衡的状态，特别是在饱和条件下，安全系数（Factor of Safety, FoS）接近临界水平，且应变局部化（strain localization）更为明显。研究结果显示观测变形模式与计算稳定性条件之间存在明确关联，结构不连续面（structural discontinuities）和地下水在控制边坡行为中起主导作用。基于上述发现，研究人员提出了一种将监测数据与稳定性评估相链接的集成工作流程，实现了临界区域的识别，并支持在 ongoing mining operations 中对边坡条件进行评估。该方法有助于提高决策可靠性，并支持更安全、更可持续的露天采矿实践。
**研究背景与问题**
全球经济发展高度依赖矿产资源，其在基础设施建设和能源转型技术中至关重要。露天开采作为金属矿产最主要的 surface mining 技术，凭借高效的作业能力和生产速率占据主导地位。然而，露天矿边坡稳定性是一个受多重因素影响的动态地质力学过程：随着系统性开挖，围岩约束减小，剪应力 tends to localize in pre-existing weakness planes，可能引发局部破坏乃至大规模垮塌。边坡性能不仅影响作业安全，还涉及设备暴露度、剥采比及长期经济边界，因此边坡稳定性评估已成为设计与实时矿山管理的核心环节，而非一次性分析工作。
传统边坡稳定性评估包含工程分析（量化潜在破坏）与现场测量（反映实际墙坡行为）双重轨道。典型方法包括提供安全系数的极限平衡分析，以及可研究应力重分布模式、变形响应和强度折减型稳定性指标的数值方法。随着开挖阶段和地下水条件的变化，监测成为关键反馈机制。地面雷达系统提供实时变形测量，卫星InSAR则提供更大范围、更长时间尺度的信息。尽管近年研究逐渐从纯几何安全评估转向基于机制的变形响应解释，但现有研究多将变形检测、预测、监测可行性或稳定性估计视为独立任务，缺乏系统整合。
**研究开展与核心结论**
本研究的核心贡献不在于单独使用雷达、InSAR、极限平衡分析或有限元建模，也无意提出新的InSAR处理算法或数值稳定方法；其创新在于将这些组件组织于一个回顾性的、基于监测的地面控制框架内，在同一工程解释流程中链接短期雷达位移行为、长期InSAR变形趋势以及基于FoS/SRF的稳定性指标。该工作流程旨在作为运营露天矿边坡管理的决策支持程序，使临界区域能够依据监测证据和计算稳定性条件进行识别、解释和再评估。
研究背景设定于土耳其北部一处年产约2 Mt矿石的露天铜矿，其高墙高度近700 m，具有复杂边坡几何形态、非均质岩性单元和构造控制岩体条件。矿区位于 Central Pontides 构造带，主要由绿泥石-绿帘石片岩、千枚岩单元和变质玄武岩层组成，NW向倾伏的叶理面和构造接触面在更大尺度上控制边坡力学响应。
**关键技术方法**
研究采用IBIS-Rover干涉雷达系统进行地基雷达监测（Ku波段17.05–17.35 GHz，2022年4月至11月），获取视线方向（Line-of-Sight, LOS）位移，精度±0.1 mm，范围分辨率0.75 m；采用Sentinel-1升降轨数据经ATLAS InSAR处理流程进行时序分析（2022年11月至2023年6月），通过Goldstein滤波和相干像元选取，将LOS位移分解为垂直与东西向分量；变形行为依据位移-时间关系和速度趋势进行分类（退缩型、渐进型、过渡型）。稳定性分析采用Bishop简化法和Janbu简化法进行极限平衡计算，并采用Rocscience RS2中二维平面应变模型的强度折减技术（Shear Strength Reduction, SRF）进行有限元分析，基于Hoek-Brown准则转化等效Mohr-Coulomb参数，分别评估干燥和饱和地下水条件。
**研究结果**
**雷达监测结果**：累积位移图显示变形集中于三个主要区域。位移时间序列揭示渐进变形模式：早期位移较低，4月下旬起明显加速，5月至9月持续累积，呈阶段性变形特征。速度记录显示多次短暂加速-减速事件，与局部应变调整而非全球性破坏相关。降雨与短期位移增量存在时间关联，但仅为近表面剪切强度的临时降低而非主要触发机制。10月后位移速率和速度值显著下降，进入渐进稳定阶段。据此识别出初始响应、渐进 mobilization、间歇性变形和最终稳定四个变形阶段。
**InSAR结果**：雷达监测期后的卫星观测显示变形持续，累积东西向位移约36.8 mm，垂直位移约21.7 mm，水平分量占优，水平/垂直位移比1.70，表观倾伏角约30°，反映以侧向位移为主的变形机制，与NW向倾伏构造叶理控制的滑动几何学兼容。
**稳定性分析结果**：极限平衡分析显示，干燥条件下Bishop法和Janbu法最小FoS分别为1.15和1.12；饱和条件下分别降至1.06和1.00，接近临界状态。有限元分析在饱和条件下获得临界SRF值0.94，最大剪应变显示变形局部化于上部边坡区并沿构造控制的软弱带扩展，总位移模式与之一致。虽然SRF<1不应解读为现场已发生全球崩溃，但表明在二维模型假设下，该剖面具有极低的稳定性裕度，易发生局部渐进变形。
**监测与稳定性结果的定量比较**：雷达监测的三个变形集中区（Area 1–3）与FEM最大剪应变集中区位于同一关键边坡 sector，形心偏移小于40 m，空间吻合良好。InSAR位移矢量与模型化的滑动几何学运动学兼容。饱和条件下FoS降低8–10%与监测记录的地下水敏感性一致。FoS描述预设滑面上力/力矩平衡，SRF反映临界应变局部化的发生水平，两者为低稳定性裕度的互补指标。
**内部一致性检验**：第一，雷达与InSAR两种基于不同传感器、几何、处理链和观测周期的独立监测数据集在变形 sector 位置和位移持续性上达成一致；第二，极限平衡和有限元两种基于不同假设的稳定分析方法在同一关键区域识别出低稳定性裕度；第三，监测变形区与FEM应变局部化模式空间重合，InSAR位移矢量与解释滑动几何学运动学兼容，地下水敏感性解释与FoS降低及SRF<1一致。
**讨论**
雷达与InSAR的联合使用提供了更完整的边坡变形图景：雷达有效检测短期位移变化和活动期的速度波动，InSAR则揭示长期变形演化。两者作为连续监测记录表明，变形在雷达监测活动期后并未终止，而是以降低速率持续。这种独立监测技术的一致性支持位移模式为真实边坡响应而非单一技术伪影的解释。
变形行为与稳定性条件的关联在将监测结果与分析结果对比时尤为清晰。观测到的位移和速度模式表明为阶段性渐进变形，与饱和条件下接近 unity 的安全系数一致；有限元分析显示局部应变发展而未形成明确全球破坏面，表明边坡在极限平衡状态附近运行，变形由渐进应力重分布而非突然崩溃控制。
近 unity 的FoS和 sub-unity 的SRF反映饱和条件下模型剖面的低稳定性裕度。需注意，FoS为静态比值，描述沿假设破坏面抗滑力与下滑力的平衡关系，但不描述变形的速率、样式或时机。临界SRF值在二维数值框架内解释，表明在给定模型假设下边坡接近临界状态，通过局部构造控制变形而非快速连续破坏面发展来调节。这与监测记录中渐进、间歇性位移而无突然全球破坏证据的特征一致。
地下水是该敏感性的主控因素之一。饱和分析虽为保守边界而非标定地下水模型，但揭示了孔隙水沿不连续面、破碎带或软弱岩性接触进入时，有效正应力降低和抗剪强度减弱的影响。监测记录中降雨与间歇性变形的时序关联支持了地下水敏感边坡行为的解释。
从运营角度，该结果凸显了将监测系统与稳定性分析整合于露天矿管理中的重要性。所提工作流程不依赖每次监测事件后对强度参数的重新标定，而是将计算FoS和SRF值与雷达、InSAR观测共同作为运营期间的低稳定性裕度指标。当此类值与监测变形重合时，应触发对受影响 sector 的更密切检查和设计或运营措施的重新评估（如台阶几何、排水条件、开挖顺序或 access control）。该方法提供了观测变形与工程稳定性指标之间的直接联系，支持运营期间边坡条件的重新评估，有助于更安全、更可持续的地面控制实践。
研究结果具有明确的现场和事件特定性：位移量级、FoS值、SRF值、地下水敏感性和变形阈值均源于单一矿山和单一记录的渐进变形事件，不应直接转移。本研究的-plan transferable contribution 在于解释序列本身：从监测数据检测局部变形，检验其与稳定性分析的力学一致性，并利用综合证据支持运营期间的边坡再评估。未来需在更多矿山 site、不同变形机制和具备历史破坏记录的 retrospective back-analysis 中检验该工作流程的广泛适用性。
**研究结论部分翻译**
本研究通过结合地基雷达监测、卫星InSAR观测和数值稳定性分析，考察了露天采矿环境中的边坡行为。结果表明，变形随时间渐进发展而非以单一破坏事件形式出现，位移集中于边坡特定区域。雷达数据捕获了短期变化和间歇性加速阶段，而InSAR观测表明主变形期后低速率变形持续，累积水平和垂直位移分别达到约36.8 mm和21.7 mm。
数值分析与监测结果一致。极限平衡分析显示边坡接近临界条件，干燥条件下安全系数从1.15–1.12降至饱和条件下的1.06–1.00。有限元结果进一步支持该行为，饱和条件下临界SRF值为0.94，变形沿构造控制带局部化。这些发现表明边坡响应 strongly influenced by groundwater conditions and structural discontinuities。饱和条件下，边坡在极限平衡附近运行，降低的静态稳定性裕度表现为渐进、构造控制的变形——而非突然全球破坏——可通过集成雷达和InSAR监测进行检测和追踪。
本研究的主要成果是变形监测与稳定性分析的综合解释。结果表明，雷达和InSAR导出的变形模式可直接与稳定性条件关联，从而对运营条件下的边坡性能进行更一致的评估。这使得临界区域识别和采矿过程中边坡行为演变的更好理解成为可能。
所提工作流程旨在作为将监测数据与露天矿运营中稳定性评估相集成的实用框架。在此框架内，工作流程不依赖对每个监测记录的强度参数常规重新标定。相反，低FoS或SRF值被作为定量警示指标，应导致受影响边坡 sector 的强化监测以及对设计、排水和运营措施的审查。虽然本研究未进行明确的边坡重新设计，但该方法通过将观测变形与工程稳定性指标关联，为此类决策提供了基础。以此方式，它支持运营期间边坡条件的重新评估，并有助于更安全、更可持续的地面控制实践。本研究以一处露天矿和一处记录的渐进变形事件展示了所提工作流程的应用。案例研究中的数值结果具有现场和事件特定性；然而，解释序列可适应于其他具备监测数据和稳定性分析资料的露天矿运营。未来工作应将该工作流程应用于更多矿山 site、不同变形机制，以及在存在 post-failure 数据集的情况下进行历史边坡破坏的回顾性反分析。
打赏
