工程爆破中震源能量计算新方法.docxVIP

工程爆破中震源能量计算新方法 微地震监测技术以时间、远程、动态和三维为特点，在矿山动力灾害的监测和报警中发挥着越来越重要的作用。作为描述岩石破裂释放的灵活性能力的定量，源能量是源“时空强”的三个参数之一。合理计算源能量是分析岩石破裂场和岩体动态气候变化的基础，因此，源能量的计算具有重要意义。 微震震源能量计算主要借鉴地震学的原理及方法,如Gutenberg-Richter法、能量密度法、震动持续时间法、地震图积分法和里氏震级法.一方面这些方法比较复杂,需要参数众多;另一方面,矿山微震在震源机制、震级、震相、地震波频率、检波器类型等方面又不同于天然地震,不能直接照搬天然地震的方法及参数.目前,适用于矿山微震的经验公式还没有建立起来.震源地震波能量研究最大的难点在于震源能量的无法直接量测性和能量衰减规律的复杂性,相较于震源定位研究,震源能量计算是微震监测领域的难点和薄弱环节. 本文的主要思路是借鉴爆破工程中有关地震波能量传播规律的研究成果,依据测点峰值振动速度、震源距及震源能量之间的统计规律反演震源能量,为微地震震源能量计算提供一种可行的方案. 1 微震监测技术 峰峰集团梧桐庄煤矿182106工作面存在突水危险性,为了避免突水灾害的发生,峰峰集团与北京科技大学合作,采用北京科技大学研制的微地震监测系统BMS对开采过程中的岩层破裂规律、构造活化过程等进行监测,开展微震监测技术用于底板突水的监测和预警研究. 在检波器布设优化研究的基础上,结合现场施工条件,沿工作面上下顺槽,在顶、底板中布设检波器,检波器的平面投影如图1所示. 为进行微震监测系统的校验和相关参数的计算,进行了放炮标定.采用煤矿许用乳化炸药16.6 kg,齐发爆破方式,爆破后各检波器获得了清晰的地震波形图,如图2所示,各测点的主要数据如表1所示. 2 计算源能量的方法 2.1 振幅平方和和能量的线性回归 要精确地计算出地震波传播过程中,各个检波器位置处地震波能量密度是非常困难的,工程中通常采用某一波形参数来表示测点能量大小,如最大峰值振动速度、峰值振动速度平方、地震波持续时间内各采样点振幅的平方和等,一般认为地震波持续时间内各采样点振幅平方和表示能量更准确.利用表1中数据,计算得测点峰值振动速度与能量(地震波持续时间内采样点振幅平方和)之间的线性相关系数r=0.951,大于0.7,说明测点峰值速度与能量之间具有高度线性相关性.进行峰值速度–能量的线性回归,回归曲线及散点图如图3所示. 如此高的线性相关性,说明在地震波能量传播规律研究中,完全可以用质点振动峰值速度作为测点能量的标志,而无需计算微震事件起止时间段内的振幅平方和.这样一方面可以减少计算工作量;另一方面因为峰值速度的计算不需要事件起止时间的拾取,可以实现自动计算,有利于提高微震监测系统整体的效能. 2.2 地震波能量衰减规律 2.2.1 地震地表震波的吸收衰减 影响地震波传播过程中能量衰减的主要因素为波前扩散和吸收衰减.在均匀介质中,点震源的波前为球面,随着传播距离的增大,球面逐渐扩展,但总能量仍保持不变,而使单位面积上的能量减小,振动的最大振幅与震源距呈反比关系,表示为 式中,A为测点最大振幅,A0为震源振幅,r为震源距. 由于实际岩层并非理想弹性介质,地震波在传播过程中会因摩擦、液体流动、黏性张弛等原因造成衰减,称之为吸收衰减.吸收衰减的程度可用单位距离内振幅的衰减率表示,称为吸收系数β.理论上已经证明,吸收系数是频率的函数,对某一种岩石介质,在相同的振幅下,频率越高,吸收系数越大.因此,地震波随传播距离的增大,高频部分很快被吸收,而只保留较低的频率成分,因此吸收系数也是震源距的函数.地震波的吸收衰减表示为 地震波受波前扩散和吸收衰减的综合影响,可表示成 不同震源距位置处β值不同,在天然地震的震动衰减规律计算中通常采用分段函数的形式,而在微震监测中,由于震源距范围有限,为简化计算,在一定的震源距范围内采用固定β值. 2.2.2 衰减系数q 爆破工程中质点振动速度衰减规律常用萨道夫斯基经验公式表示,即 式中:V为质点振动速度,cm·s-1;r为震源距,m;K为场地系数;α为衰减系数;Q为最大段药量,kg;V0为震源初始振动速度,cm·s-1. 这里V0只是回归曲线统计意义的震源振动速度,或称为虚震源振动速度,因为在爆破近区质点振动速度呈现完全不同的规律特性. 对于表1中数据,回归得K=111.45,α=1.3. 工程爆破中经常采用萨氏公式预测质点振动速度或计算最大允许炸药量,也是我国《爆破安全规程》(GB6722—2003)推荐使用的回归方式. 2.2.3 小尺度下地震波能量衰减规律 可用于地震波能量衰减规律的回归模型主要有两种:第一种幂函数形式,如式(4);第二种指数函数形式,如式(3).采用表1中数据,将