一种基于最小距离的微震或声发射异常到时数据识别方法与流程

本发明涉一种基于最小距离的微震或声发射异常到时数据识别方法。技术背景微震及声发射技术广泛应用于采矿过程及地下施工过程中的监测领域，为地下工程的进行提供安全信息。在微震及声发射定位过程中，传感器会接收到一些异常的到时信息(即传感器感应到微震或声发射信号的时刻)，这些信息如果直接用来进行微震及声发射定位，将会对定位结果产生极大的误差。因此，如何识别这些异常的到时数据并进行排除，在微震及声发射定位技术中意义重大。有必要设计一种异常到时数据识别方法。技术实现要素：本发明所解决的技术问题是，针对现有技术的不足，提出一种基于最小距离的微震或声发射异常到时数据识别方法，易于实现、准确度高。本发明的技术解决方案如下：一种基于最小距离的微震或声发射异常到时数据识别方法，包括以下步骤：步骤1：到时数据初步定位：将微震或声发射事件中，n个传感器的n个到时数据依次记为t1、t2、t3…tn；分别排除t1、t2、t3…tn中的一个，得到n组数据，每一组数据中包括n-1个到时数据；对n组数据分别利用定位方法进行定位，得到n组初步定位结果；步骤2：到时数据再定位：在n个到时数据中任意选取m个数据，利用定位方法进行再次定位，获得共组定位结果；步骤3：求取频率分布：分别求出步骤2的组定位结果中，各个x、y、z坐标值的出现频率，进而选取出现频率最大的x、y、z坐标值，作为最优的定位结果O*(x*,y*,z*)；步骤4：距离计算：分别求出步骤1所得n组初步定位结果与最优定位结果的距离d；步骤5：异常到时数据识别：比较得到n组初步定位结果与最优定位结果的距离d的最小值dmin，dmin对应的初步定位结果相应的一组数据所排除的到时数据即为异常到时数据。所述定位为利用解析定位方法或迭代定位方法其中的一种进行定位。所述m＝6。有益效果：本发明对微震或声发射所得n个到时数据选用6个到时数据并利用解析定位算法或数值定位方法的其中一种对其进行定位处理，以获得个定位结果。创新点在于利用个定位结果求其频率分布，并选取出现频率最高的x、y、z值作为最优定位结果，并以初步定位结果与最优定位结果的距离作为判据，当距离最小时，表明该组数据所得到的定位结果更接近于最优结果，进一步表明该组数据所排除的数据为异常到时数据。本发明很好地利用频率分布的优势，从一系列定位结果中能够准确地确定最优解，从数学的角度保证了选取的最优定位结果的可靠性和准确性。以此为基础利用距离判定各组数据的可靠性也突破了传统的判定方法，具有易于实现、可靠度高等优点，在微震及声发射信号识别领域能够广泛运用。附图说明图1是实施例1算法流程图；图2是实施例1传感器布置示意图。图3是实施例1中x的频率分布图。图4是实施例1中y的频率分布图。图5是实施例1中z的频率分布图。具体实施方式以下将结合附图和具体实施例对本发明做进一步详细说明：实施例1：以下基于本实施例对本发明基于最小距离的微震或声发射异常到时数据识别方法进行详细说明。如图2所示，假定某次声发射试验系统有8个传感器位于正方体的8个顶点上，取各传感器坐标A(0，0，0)，B(800，0，0)，C(800，800，0)，D(0，800，0)，E(0，0，800)，F(800，0，800)，G(800，800，800)，H(0，800，800)，单位均为m。试验中测得某点O处产生的声发射信号到达各传感器的时刻分别为139.2561ms、164.9302ms、175.1945ms、151.4098ms、103.2400ms、131.4610ms、144.3272ms、113.6219ms。具体实施步骤如下：(1)将8个传感器获得的8个到时数据由小到进行排列，可得：tE＝103.2400ms；tH＝113.6219ms；tF＝131.4610mstA＝139.2561ms；tG＝144.3272ms；tD＝151.4098mstB＝164.9302ms；tC＝175.1945ms(2)依次排除tE,tH,tF,tA,tG,tD,tB,tC，利用余下的7个到时数据进行定位，可以得到8个初步定位结果：OE(258,336,580)；OH(345.15,375.14,470.34)OF(329.73,371.85,487.95)；OA(339.97,372.73,485.23)OG(354.31,381.48,460.65)；OD(346.86,379.00,465.82)OB(340.57,375.73,475.74)；OC(350.65,379.12,462.82)(3)对所有的8个到时数据，选取6个利用解析定位法或数值定位方法进行定位，共可获得个定位结果，如表1所示。(4)利用28个定位结果，分别求出x、y、z的频率分布，进而选取出现频率最大的x、y、z值，作为最优的定位结果(x*,y*,z*)。由图3、图4及图5可以看出，最优定位结果为O*(258,336,580)。(5)计算8个初步结果与最优定位结果之间的距离，如表2所示，可以得出，OE组定位结果距最优结果距离最短，该组数据中排除的到时数据为tE，可以推知传感器E的到时数据为异常到时。表1算例的28组定位结果(单位：m)表2初步定位结果与最优定位结果的距离(单位：m)定位编号OEOHOFOAOGODOBOCO*2144.25120.92129.33159.65149.77137.62154.30以上所述仅为本发明的实施例而已，并不用以限制本发明，凡在本发明精神和原则之内，所作任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。当前第1页1 2 3