一种新的矩阵传播方法，用于理论分析具有平行冰填充裂缝的低温梯度层状岩体中的波传播现象

《International Journal of Rock Mechanics and Mining Sciences》：A new matrix propagation method for theoretical analysis of wave propagation in low-temperature gradient layered rock mass with parallel ice-filled joints

高海拔寒冷地区露天矿爆破中，冰裂隙层状岩体中应力波传播机制研究。提出融合矩阵传播法（MPM）与非线位移连续模型（DDM）的新理论模型，揭示温度梯度与多次反射对波传播的影响规律，发现低温梯度抑制波传播且冰裂隙显著加剧能量耗散，间接传播效应随入射频率和裂隙初始模量变化，为优化爆破参数提供理论支撑。

　　
高海拔寒冷地区露天矿爆破作业中应力波传播机制研究

摘要部分揭示了研究的核心问题：在低纬度高原寒区，由于地表与岩层间热交换效率差异，形成独特的分层温度场分布，这显著改变了爆破应力波在含冰层理岩体中的传播特性。研究创新性地将矩阵传播法（MPM）与非线性的位移间断模型（DDM）相结合，构建了适用于多组平行冰岩界面反射叠加效应的理论模型。通过对比传统单界面模型与现有数值模拟方法，验证了该模型在预测应力波多次反射导致的能量传输增强效应方面具有显著优势。研究证实低温环境与冰层界面共同作用，使应力波在传播过程中产生显著的能量衰减现象，其衰减程度与入射波频率、界面初始刚度及冰层填充数量存在非线性关联。

引言部分系统梳理了该领域的研究进展。现有理论模型存在三大局限性：其一，单界面模型无法处理多组平行界面间的复杂反射叠加效应；其二，传统矩阵传播法假设界面完全粘合，无法表征冰岩界面间的动态滑移与脱粘行为；其三，现有方法未充分考虑温度梯度导致的波阻抗空间异质性。针对这些问题，本研究建立了包含以下创新要素的理论体系：
1. 非线性位移间断模型（DDM2.0）：通过引入动态刚度修正系数，准确描述冰岩界面在冲击荷载下的应变硬化与脆性破碎行为。该模型特别考虑了低温环境下冰体粘弹性特性随时间的变化规律。
2. 矩阵传播法改进算法：开发适用于多界面反射叠加的快速矩阵合成技术，将传统需逐次计算的反射过程转化为矩阵相乘操作，计算效率提升两个数量级。
3. 温度场耦合分析模块：建立温度梯度-波阻抗-应力波传播的耦合数学关系，揭示低温环境如何通过改变岩体波阻抗分布影响应力波传播路径。

模型验证部分采用SHPB试验进行多维度验证，测试结果显示理论预测值与实测数据的吻合度达92%以上。通过设置不同参数组合的对比试验，证实了模型对以下关键问题的表征能力：
- 含冰层理界面间的非对称反射效应
- 多次反射波叠加产生的相干增强现象
- 低温环境导致的波阻抗梯度变化规律

案例研究部分选取典型高原矿区进行工程验证，分析发现：
1. 当入射波频率低于15Hz时，多界面反射效应可使应力波能量传输系数提升40%-60%
2. 界面初始刚度下降至岩石基体值的30%以下时，能量衰减速率降低约25%
3. 含冰界面数量超过5组时，间接能量传输效应开始呈现非线性衰减特征

工程应用建议包括：
- 优化装药结构：在层理密集区域采用预裂减能技术
- 调整起爆时序：根据温度梯度分布设置不同区段的微差间隔
- 优化钻孔布局：在冰层密集区采用交错式布孔方案
- 建立动态监测系统：实时监测应力波在层理界面的反射系数变化

该研究突破了传统单界面模型的理论框架，建立了考虑温度梯度耦合作用的多界面反射叠加模型。在理论层面，揭示了低温环境下应力波传播的"三重耦合效应"：温度梯度-波阻抗分布-界面损伤演化的动态耦合关系。在工程应用方面，提出了基于界面反射特性的爆破能量调控方法，为寒区露天矿爆破安全高效作业提供了新的理论依据。研究结果已成功应用于玉龙铜矿、甘蒙金矿等高原矿区，使爆破效率提升18%-22%，周边岩体损伤降低35%-45%。

后续研究建议：
1. 开展多场耦合试验：集成温度-应力-渗流耦合监测系统
2. 开发智能预测模型：基于机器学习构建动态参数反演算法
3. 研制新型冻融防护材料：提升含冰层理岩体的抗爆性能
4. 建立寒区爆破数据库：积累不同地质条件下的反射系数实测数据

该研究在基础理论层面实现了三个突破：首次系统揭示温度梯度对多界面应力波传播的调控机制；创新性地建立冰岩界面动态力学参数辨识方法；开发了适用于寒区复杂地质条件的爆破效应预测系统。这些理论成果和实践经验为高寒地区深部开采工程中的爆破技术优化提供了重要支撑，对保障极端环境下的采矿作业安全具有重要参考价值。