基于力学的冰崩启动判据：连接断裂参数与冰川尺度失稳

《International Journal of Mining Science and Technology》：A mechanics-based initiation criterion for ice avalanches: Linking fracture parameters to glacier-scale instability

　　

冰川速度是冰川灾害防灾减灾的关键指标，但其与冰崩（IA）的内在联系仍未得到充分认识。为阐明连接冰川速度与IA的物理机制，本研究利用Sentinel-2B遥感影像分析了色东普谷冰崩前的冰川运动特征。据此，研究人员提出一种力学路径，即冰川加速提高了冰的变形速率，降低其断裂韧性，从而促进IA启动。低温四点弯曲（4 PB）实验定量揭示了冰的断裂韧性和裂纹尖端屈服强度（CTYS）随应变率和温度升高而逐渐降低的规律。无支撑冰体被理想化为含单条裂缝的悬臂梁。几何敏感性分析表明，悬伸长度和裂缝长度是控制拉应力集中的主要因素，为悬冰川IA启动的拉伸破坏准则奠定了基础。这些结果揭示了一种力学路径，即冰川加速降低冰的裂纹扩展阻力，从而促进IA启动。针对艾格（Eiger）悬冰川事件的验证进一步表明，冰川加速可作为应变率增加、断裂韧性退化及渐进损伤累积的外部标志。本研究为阐明悬冰川拉伸型IA的启动机制提供了新视角。

本文解读了《International Journal of Mining Science and Technology》上发表的题为“基于力学的冰崩启动判据：连接断裂参数与冰川尺度失稳”的研究。冰崩（Ice Avalanche, IA）是高海拔冰川区频发的破坏性块体运动过程，常诱发碎屑流、河道堵塞和冰川湖溃决洪水（Glacial Lake Outburst Flood, GLOF）等链生灾害。尽管已有大量研究关注其动力学过程与灾害链模拟，但IA的启动机制仍缺乏力学层面的统一框架。现有观测研究多集中于表面速度模式或经验阈值，较少将实验室尺度的断裂参数转化为可用的失稳判据；同时，观测手段在偏远高山区受恶劣天气、云层遮挡和成本限制，难以捕捉内部应力重分布与断裂抗力退化过程。因此，建立基于物理机制的IA启动判据，对冷山区工程风险管理和基础设施规划具有重要意义。

该研究以西藏色东普谷为研究区。该流域面积约66.89 km2，高差达4548 m，在2017—2018年间多次发生IA诱发的泥石流事件，尤其2017年11月18日林芝地震后事件频率显著增加。研究人员利用Sentinel-2B遥感影像，通过特征追踪计算了冰崩前的冰川运动速度。结果显示，IA事件主要发生在冰川速度较高或加速阶段，尤其是地震后的数月内。这一观测支持了研究假设：当地震扰动改变冰川运动特征、提高冰体变形速率时，冰的力学抗性被削弱，进而促进IA启动。

为量化这一力学退化路径，研究人员采用四点弯曲（4 PB）实验。实验在?2、?10、?20、?30 °C四种温度与0.001、0.01、0.1、1、10 mm/s五种加载速率下，对200 mm × 40 mm × 40 mm的预制裂纹柱状冰试样进行测试，共计60组有效试验，并辅以声发射（Acoustic Emission, AE）监测。基于峰值载荷和试样几何计算了临界I型应力强度因子（断裂韧性, K_IC）与裂纹尖端屈服强度（CTYS）。结果发现，K_IC和CTYS均随应变率增加呈对数下降，随温度降低而升高。AE事件率与损伤因子演化表明，高加载速率压缩稳定裂纹扩展阶段，使试样更趋于瞬时脆性失稳；低加载速率下损伤演化更加渐进，表现为累积损伤特征。

研究结果部分包含以下关键发现：

**峰值载荷与损伤演化**：4 PB试验中，峰值载荷随加载速率增加从1786.05 N单调降至590.21 N，降幅达66.95%，表明冰体存在显著的应变率弱化效应。AE监测将破坏过程划分为弹性变形、裂纹萌生、稳定扩展和破坏汇合四个阶段。高加载速率或低温下，稳定扩展阶段缩短，AE事件集中释放，损伤因子突增，呈明显脆性失稳特征。

**断裂韧性与CTYS**：冰的K_IC与应变率呈对数负相关，各温度下拟合决定系数R2≥0.93；温度从?2 °C降至?30 °C时，K_IC显著增大。CTYS同样随应变率增加而降低，且数值远高于宏观抗拉强度，反映裂纹尖端强烈的应力集中和塑性约束。CTYS与AE监测耦合表明，CTYS较低的试样对应更密集的AE爆发和更短的前兆阶段，能灵敏指示微观失稳机制。

**力学模型与启动判据**：研究人员将无支撑冰体理想化为含单条表面裂缝的悬臂梁，利用Airy应力函数推导出冰体内应力场解析解，提出双重启动判据：一是应力驱动失稳，即局部应力超过CTYS；二是断裂抗力退化，即冰川加速降低CTYS，使有效抗力低于既有应力。几何敏感性分析显示，悬伸长度和裂缝长度是控制拉应力集中的主导因素，冰厚度和床坡角影响相对较弱。

**模型验证**：以2016年瑞士艾格峰（Eiger）悬冰川IA事件为例，该事件约15000 m3冰体于8月25日脱离。结合冰震记录与冰川速度数据，分别基于冰震事件数和冰川速度推算损伤因子演化。两种方法均显示IA前损伤持续累积，并于24日出现明显加速，与现场观测的裂缝C扩展一致，验证了模型合理性。

在讨论部分，研究人员对比了巴西劈裂、三点弯曲、楔形劈裂和四点弯曲四种方法测得的冰断裂韧性，指出4 PB法在趋势一致性、数值适中性和分散度方面表现最优。进一步分析表明，应变率对K_IC的影响强于温度；K_IC随应变率呈三阶段演化：先准线性下降、再快速下降、最后进入平台，类似蠕变曲线，暗示断裂退化与时间相关变形过程存在潜在类比。基于此，提出应变率驱动的IA启动概念模型：外部扰动（降雨、升温、地震）加速冰川流动，提高应变率，降低断裂韧性，当驱动力超过CTYS时，裂缝进入失稳扩展并触发IA。同时，讨论也指出研究局限性：实验室冰与天然冰在尺度、杂质和结构上存在差异；悬臂梁模型简化了基底接触和侧向约束；尚未考虑热-水-力耦合及多裂缝相互作用；未来需更多案例验证以提升模型普适性。

研究结论翻译如下：
(1) Sentinel-2B遥感影像识别出色东普谷IA前的异常冰川运动特征。结果表明，IA事件主要发生在冰川高速或加速阶段。2017年林芝地震后，色东普谷冰川运动显著增强，随后数月内发生多次IA诱发的泥石流事件。
(2) 低温条件下含预制裂纹冰试样的四点弯曲实验结合AE监测，揭示了冰断裂参数对应变率和温度的响应。当加载速率从0.001 mm/s增至10 mm/s时，平均峰值载荷从1786.05 N降至590.21 N，降幅达66.95%。这表明冰川加速可降低冰的裂纹扩展阻力，从而促进IA启动。
(3) 建立了悬冰川拉伸型IA启动的力学判据，识别出两条启动路径：一是无支撑冰体的几何演化增强局部拉应力集中，当应力超过CTYS时触发IA；二是冰川加速提高应变率并降低断裂韧性和CTYS，当退化的抗裂阻力低于既有应力水平时发生破坏。
(4) 几何敏感性分析和艾格悬冰川案例验证进一步确定了IA启动的关键控制因素和模型适用性。结果表明，悬伸长度和裂缝长度是控制拉应力集中的主导因素。这些结果提示，冰川加速可被解释为应变率增加、断裂韧性退化和渐进损伤累积的外部表现，为理解悬冰川拉伸型IA启动提供了力学基础。