深埋隧道开挖过程微震事件空间分形维数计算方法
【专利摘要】本发明涉及一种深埋隧道开挖过程微震事件空间分形维数计算方法。首先根据深埋隧道开挖过程中微震事件沿洞轴、洞径方向距离掌子面的产生范围确定总计算区域沿洞轴方向总长度L、总计算区域及部分计算区域沿洞径方向断面长度H，然后统计总计算区域内所有微震事件数、部分计算区域内的微震事件偶对的数目，并根据所选取的部分计算区域范围长度及其所对应的相关指数建立直角坐标系进行线性拟合，如果两者之间具有很好的线性相关性，其斜率为隧道开挖过程中微震事件空间分布的分形维数值。采用以上方法可以定量分析隧道开挖过程中掌子面应力型灾害孕育及发生过程中微震事件空间分布的演化机制，对灾害的预测与防治也具有重要的指导意义。
【专利说明】深埋隧道开挖过程微震事件空间分形维数计算方法
【技术领域】
[0001] 本发明涉及隧道开挖过程中的微震监测方法。尤其是一种深埋隧道开挖过程微震事件空间分形维数计算方法。适用于深埋隧道开挖灾害孕育及发生过程中破裂空间分布特征及演化规律的定量描述，属岩土力学【技术领域】。
【背景技术】
[0002]随着隧道工程岩体深度的不断增加，工程灾害日趋增多，如岩爆、塌方、巷道围岩开裂、片帮等，对工程的安全性造成了巨大威胁，因此深部隧道工程中的岩石力学问题也是目前国内外岩石力学界研究的焦点。针对深部岩体隧道工程展开研究的方法有很多，其中以微震监测法应用的最为广泛。在深部岩体隧道开挖过程中，围岩当受到应力作用时就会在产生微破裂的同时释放应变能并产生应力波，这些应力波由微震仪器接收并定位，运用微震监测的手段可以测量这些微破裂的发生及发展情况，在灾害发生之前对这些物理机制点进行测量从而达到对灾害的孕育和发生过程进行分析的目的。
[0003]Fractal character and mechanism of rock bursts (International Journalof Rock Mechanics and Mining Sciences, 1993, 30(4): 343 — 350) 一文中针对矿山开采过程中采矿区微震事件的分布特征提出了球形半径r的分形测量方法，并运用分形维数成功的对深部矿山开采过程中岩爆孕育及发生过程中的微震活动性进行了定量描述，取得了显著成果。但是上述计算方法无法根据深埋隧道的线性特点有效的针对深埋隧道开挖过程中随着掌子面的不断推进呈线性展开的微震事件进行测量与计算，并且容易将相邻洞段及掌子面开挖过程中所产生的微震事件计算其中使结果产生较大偏差。
【发明内容】
[0004]针对上述存在问题，本发明的目的在于提供一种深埋隧道开挖过程微震事件空间分形维数计算方法。
[0005]为了实现上述目的，本发明的技术解决方案为:
深埋隧道开挖过程微震事件空间分形维数计算方法，所述计算方法的具体步骤包括: a)根据深埋隧道开挖过程中微震事件在掌子面后方沿洞轴方向距离掌子面分布范围 ￡
长度确定总计算区域沿洞轴方向总长度为L ;微震事件在掌子面附近沿洞径方向距离 2
掌子面分布范围长度确定总计算区域及部分计算区域沿洞径方向正方形断面边长为
2
H0
[0006]b)以灾害发生位置为中心，统计总微震事件总计算区域范围内的所有微震事件
￡ 2L
数，一，_，…，z范围长度下部分计算区域内的微震事件偶对数目，其中η为常数。 μ η
[0007]c)根据相关指数计算公式:
【权利要求】
1.深埋隧道开挖过程微震事件空间分形维数计算方法，其特征在于:所述计算方法的具体步骤包括: a)根据深埋隧道(I)开挖过程中微震事件在掌子面后方沿洞轴方向距离掌子面分布范 T围长度一，确定总计算区域(2)沿洞轴方向总长度为L ;微震事件在掌子面附近沿洞径方向 2
H距离掌子面分布范围长度确定总计算区域(2)及部分计算区域(3)沿洞径方向正方形
2断面(5)边长为H; b)以灾害发生位置(4)为中心，统计总微震事件总计算区域(2)内的所有微震事件数，τ 2L一，_，…，z范围长度下部分计算区域(3)内的微震事件偶对数目，其中η为常数；M /I c)根据相关指数计算公式: eg )=..........公式一` “ N(N-1)Λ、
一 xL式中:χ_| 2 为总计算区域范围内的所有微震事件数)表不—范围长度
χ?.下部分计算区域内的微震事件偶对数目，计算出；长度范围内的相关指数; CO分别以1gl1, 1gl2,…，1014?横坐标，kgC(4)，，…，logC(4)为纵坐标，建立直角坐标系并进行线性拟合: D = Iimlogc^公式一
M iog(4) 当拟合直线的自相似系数大于0.9，表明岩石微裂隙的产生在空间上是具有分形分布关系的，求出的斜率D值为隧道开挖过程中微震事件空间分布的分形维数值。
【文档编号】G01V1/28GK103605160SQ201310658333
【公开日】2014年2月26日 申请日期:2013年12月9日 优先权日:2013年12月9日
【发明者】冯夏庭, 于洋, 丰光亮, 肖亚勋, 郝宪杰 申请人:中国科学院武汉岩土力学研究所