摘要
We propose a new method for inverting source function of microseismic event induced in mining. The observed data from microseismic monitoring during mining are represented by a wave equation in a spherical coordinate system and then the data are transformed from the time-space domain to the time-slowness domain based on tomographic principle, from whichwe can obtain the signals related to the source in the time-slowness domain. Through analyzing the relationship between the signal located at the maximum energy and the source function, we derive the tomographic equations to compute the source function from the signals and to calculate the effective radiated energy based on the source function. Moreover, we fit the real amplitude spectrum of the source function computed from the observed data into the co-2 model based on the least squares principle and determine the zero-frequency level spectrum and the corner frequency, finally, the source rupture radius of the event is calculated and The synthetic and field examples demonstrate that the proposed tomographic inversion methods are reliable and efficient
本文提出了一种矿山微地震震源函数反演的新方法。以球坐标系波动方程为基础,对时空域中的微地震观测信号进行层析成像投影到慢度时间域,我们可以获得微地震波成像能量最大值及其对应的震源位置、发震时间和传播速度以及慢度时间信号。通过分析能量最大值处的慢度时间信号与震源函数之间的关系,推导出了计算矿山微地震震源函数的层析成像公式和利用震源函数计算微地震有效辐射能量的公式。进一步利用最小二乘法,将震源函数的振幅谱拟合成模型震源谱,确定了微地震波的零频极限值和拐角频率,最终计算出该微震事件的震源破裂半径。利用这一方法对理论模型数据和实际资料处理,结果表明:本方法具有既能求出任意一个地震事件的震源位置、发震时间和传播速度,又能同时求出震源函数及其频率特征值和震源破裂半径等震源参数的优点,因此,本方法是一种实时快速计算方法。
基金
supported jointly by projects of the National Natural Science Fund Project(No.51174016)
the National Key Basic Research and Development Plan 973(No.2010CB226803)