[发明专利]一种微震或声发射异常到时的确定方法

权利要求书

1.一种微震或声发射异常到时的确定方法，其特征在于，包含以下步骤：

步骤一：在微震或声发射试验中，分情况对获得的n个到时数据进行初步处理，得到m组到时数据：

步骤二：利用m组到时数据，求得m组微震源或声发射源坐标；

步骤三：用Logistic概率分布函数拟合m组微震或声发射定位结果，分别求得关于微震源或声发射源坐标x、y、z的Logistic概率分布函数F_x、F_y、F_z；

步骤四：求出概率分布函数F_x<0.005、F_y<0.005、F_z<0.005时的解，并在m组定位结果中找出满足解的定位结果，称为异常定位结果，这些异常定位结果所对应的到时数据中即包含异常到时数据。

2.根据权利要求1所述的微震或声发射异常到时的确定方法，其特征在于，步骤一的具体实现过程为：

(1)若波速已知，则对于触发传感器数量n超过4的微震或声发射事件，选取4个到时数据作为一组，获得组到时数据；

(2)若波速未知，则对于触发传感器数量n超过5的微震或声发射事件，选取5个到时数据作为一组，获得组到时数据。

3.根据权利要求1所述的微震或声发射异常到时的确定方法，其特征在于，步骤二中，利用共轭梯度法或麦夸特法中的一种求得m组微震源或声发射源坐标。

4.根据权利要求1所述的微震或声发射异常到时的确定方法，其特征在于，步骤三中，Logistic概率分布函数的基本形式为：

$F_x=\frac{1}{1+e^{-\frac{x-\mathrm{\&mu;}}{s}}};F_y=\frac{1}{1+e^{-\frac{y-\mathrm{\&mu;}1}{s1}}}:F_z=\frac{1}{1+e^{-\frac{z-\mathrm{\&mu;}2}{s2}}};$

其中，μ、μ1、μ2分别为变量x、y、z的平均值，s、s1、s2分别为与变量x、y、z的标准差σ、σ1、σ2相关的尺度参数，且x，y，z分别为微震源或声发射源坐标的x值、y值、z值。