CN114814940A - 微震监测台网评价方法、装置、电子设备及介质 - Google Patents

微震监测台网评价方法、装置、电子设备及介质 Download PDF

Info

Publication number
CN114814940A
CN114814940A CN202210737968.7A CN202210737968A CN114814940A CN 114814940 A CN114814940 A CN 114814940A CN 202210737968 A CN202210737968 A CN 202210737968A CN 114814940 A CN114814940 A CN 114814940A
Authority
CN
China
Prior art keywords
error
microseismic
target
minimum magnitude
sensor
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202210737968.7A
Other languages
English (en)
Other versions
CN114814940B (zh
Inventor
常莹
张达
杨小聪
冀虎
戴锐
石雅倩
吴彦博
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
BGRIMM Technology Group Co Ltd
Original Assignee
BGRIMM Technology Group Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by BGRIMM Technology Group Co Ltd filed Critical BGRIMM Technology Group Co Ltd
Priority to CN202210737968.7A priority Critical patent/CN114814940B/zh
Publication of CN114814940A publication Critical patent/CN114814940A/zh
Application granted granted Critical
Publication of CN114814940B publication Critical patent/CN114814940B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/288Event detection in seismic signals, e.g. microseismics
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Acoustics & Sound (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Business, Economics & Management (AREA)
  • Emergency Management (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提供了一种微震监测台网评价方法、装置、电子设备及介质,包括:获取目标微震监测台网在目标监测区域内安装的传感器的传感器坐标和目标监测点坐标;基于传感器坐标和目标监测点坐标确定目标监测点监测到的微震事件的最小震级,并基于最小震级的微震事件确定第一误差;基于矿山岩体结构和目标监测区域结构确定第二误差;基于最小震级的微震事件触发的传感器对最小震级的微震事件进行源定位,得到协方差矩阵;基于第一误差、第二误差和协方差矩阵,得到目标微震监测台网的总误差。本发明能够基于实际传感器的安装位置对微震监测台网的性能进行评估,从而提高微震监测台网监测结果的准确性。

Description

微震监测台网评价方法、装置、电子设备及介质
技术领域
本发明涉及微震定位技术领域,尤其是涉及一种微震监测台网评价方法、装置、电子设备及介质。
背景技术
微震监测系统是一种覆盖面广、可实时有效监测深井矿山、隧道掘进、边坡治理等工程的技术手段。技术人员会在台网设计过程中,对微震监测台网的性能进行评估以确定最优传感器安装布置方案,但是,受施工条件限制和工作状态影响,在实际工程应用中,传感器安装点会偏离微震监测系统的最优传感器安装布置方案,使得微震监测台网的监测结果产生偏差。因此,需要对目标区域的微震监测台网的性能进行再次评估。
发明内容
有鉴于此,本发明的目的在于提供一种微震监测台网评价方法、装置、电子设备及介质,能够基于实际传感器的安装位置对微震监测台网的性能进行评估,从而提高微震监测台网监测结果准确性的判断力。
为了实现上述目的,本发明实施例采用的技术方案如下:
第一方面,本发明实施例提供了一种微震监测台网评价方法,包括:获取目标微震监测台网在目标监测区域内安装的传感器的传感器坐标和目标监测点坐标;基于传感器坐标和目标监测点坐标确定目标监测点监测到的微震事件的最小震级,并基于最小震级的微震事件确定第一误差;基于矿山岩体结构和目标监测区域结构确定第二误差;基于最小震级的微震事件触发的传感器对最小震级的微震事件进行源定位,得到协方差矩阵;基于第一误差、第二误差和协方差矩阵,得到目标微震监测台网的总误差。
在一种实施方式中,基于传感器坐标和目标监测点坐标确定目标监测点监测到的微震事件的最小震级,并基于最小震级的微震事件确定第一误差的步骤,包括:基于传感器坐标和目标监测点坐标计算震源距;基于震源距计算目标监测点监测到的微震事件的最小震级;确定最小震级的微震事件触发的传感器与目标监测点之间的传感器距离;基于最小震级的微震事件的信噪比和传感器距离,计算第一误差。
在一种实施方式中,基于震源距计算目标监测点监测到的微震事件的最小震级的步骤,包括:按照以下公式计算目标监测点监测到的微震事件的最小震级:
Figure P_220624093844469_469359001
其中,Mx 0,y 0,z 0)为目标监测点(x 0,y 0,z 0)监测到的微震事件的最小震级,r为震源距,A、B、α为与目标监测区域的介质属性相关的常数。
在一种实施方式中,基于最小震级的微震事件的信噪比和传感器距离,计算第一误差的步骤,包括:基于目标微震监测台网采集到的最小震级的微震事件的波形计算最小震级的微震事件的信噪比;基于信噪比和传感器距离,计算最小震级的微震事件触发的传感器的到时拾取误差;将到时拾取误差的标准差确定为第一误差。
在一种实施方式中,基于最小震级的微震事件触发的传感器对最小震级的微震事件进行源定位,得到协方差矩阵的步骤,包括:基于最小震级的微震事件触发的传感器,计算每个被触发的传感器的地震波走时残差;基于每个被触发的传感器的地震波走时残差的系数矩阵,确定协方差矩阵。
在一种实施方式中,基于第一误差、第二误差和协方差矩阵,得到目标微震监测台网的总误差的步骤,包括:按照以下公式确定目标微震监测台网的总误差:
Figure M_220624093844500_500648001
其中,
Figure M_220624093844565_565052001
为第一误差,C为协方差矩阵,
Figure M_220624093844596_596366002
为第二误差,
Figure M_220624093844611_611948003
Figure M_220624093844643_643212004
分布
Figure M_220624093844674_674452005
分位数。
在一种实施方式中,上述方法还包括:对监测区域进行网格划分,得到多个目标监测区域。
第二方面,本发明实施例提供了一种微震监测台网评价装置,包括:数据获取模块,用于获取目标微震监测台网在目标监测区域内安装的传感器的传感器坐标和目标监测点坐标;第一计算模块,用于基于传感器坐标和目标监测点坐标确定目标监测点监测到的微震事件的最小震级,并基于最小震级的微震事件确定第一误差;第二计算模块,用于基于矿山岩体结构和目标监测区域结构确定第二误差;第三计算模块,用于基于最小震级的微震事件触发的传感器对最小震级的微震事件进行源定位,得到协方差矩阵;第四计算模块,用于基于第一误差、第二误差和协方差矩阵,得到目标微震监测台网的总误差。
第三方面,本发明实施例提供了一种电子设备,包括处理器和存储器,存储器存储有能够被处理器执行的计算机可执行指令,处理器执行计算机可执行指令以实现上述第一方面提供的任一项的方法的步骤。
第四方面,本发明实施例提供了一种计算机可读存储介质,计算机可读存储介质上存储有计算机程序,计算机程序被处理器运行时执行上述第一方面提供的任一项的方法的步骤。
本发明实施例带来了以下有益效果:
本发明实施例提供的上述微震监测台网评价方法、装置、电子设备及介质,首先,获取目标微震监测台网在目标监测区域内安装的传感器的传感器坐标和目标监测点坐标;然后,基于传感器坐标和目标监测点坐标确定目标监测点监测到的微震事件的最小震级,并基于最小震级的微震事件确定第一误差;接着,基于矿山岩体结构和目标监测区域结构确定第二误差;之后,基于最小震级的微震事件触发的传感器对最小震级的微震事件进行源定位,得到协方差矩阵;最后,基于第一误差、第二误差和协方差矩阵,得到目标微震监测台网的总误差。上述方法结合目标微震监测台网的到时拾取能力产生的第一误差和由矿山岩体结构和目标监测区域结构产生的第二误差,对目标微震监测台网的性能进行评价,全面衡量了影响微震监测台网性能的误差来源,从而提高了微震监测台网监测结果的准确性。
本发明的其他特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
为使本发明的上述目的、特征和优点能更明显易懂,下文特举较佳实施例,并配合所附附图,作详细说明如下。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种微震监测台网评价方法的流程图;
图2为本发明实施例提供的一种目标监测区域可监测到的最小震级示意图;
图3为本发明实施例提供的一种目标微震监测台网的定位误差示意图;
图4为本发明实施例提供的另一种微震监测台网评价方法的流程图;
图5为本发明实施例提供的一种微震监测台网评价装置的结构示意图;
图6为本发明实施例提供的一种电子设备的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
目前,受施工条件限制和工作状态影响,在实际工程应用中,传感器安装点会偏离微震监测系统的最优传感器安装布置方案,使得微震监测台网的监测结果产生偏差。因此,需要对目标区域的微震监测台网的性能进行再次评估。
基于此,本发明实施例提供的一种微震监测台网评价方法、装置、电子设备及介质,能够基于实际传感器的安装位置对微震监测台网的性能进行评估,从而提高微震监测台网监测结果的准确性。
为便于对本实施例进行理解,首先对本发明实施例所公开的一种微震监测台网评价方法进行详细介绍,该方法可以由电子设备执行,诸如:电脑、智能手机、平板电脑等。参见图1所示的一种微震监测台网评价方法的流程图,示意出该方法主要包括以下步骤S101至步骤S105:
步骤S101:获取目标微震监测台网在目标监测区域内安装的传感器的传感器坐标和目标监测点坐标。
在一种实施方式中,可以先对监测区域进行网格划分,得到多个目标监测区域,然后对每个目标监测区域内的微震监测台网的性能进行评价。在具体应用中,由于本发明实施例提供的方法是对已经安装完成的微震监测台网的性能进行评价,因此,可以获取现场微震监测台网工程安装数据,从而得到目标微震监测台网在目标监测区域内安装的传感器的传感器坐标,同时根据实际需求确定目标监测点坐标。
步骤S102:基于传感器坐标和目标监测点坐标确定目标监测点监测到的微震事件的最小震级,并基于最小震级的微震事件确定第一误差。
在一种实施方式中,微震监测台网监测到的有效触发的微震事件需要触发四个以上的传感器,基于此,本发明实施例中可以根据传感器与目标监测点之间的震源距,结合地震波震级与传播距离经验关系计算目标监测点可以监测到的最小震级的微震事件,并根据监测到的最小震级以及目标微震监测台网记录的波形确定微震事件的到时拾取误差,即第一误差。
步骤S103:基于矿山岩体结构和目标监测区域结构确定第二误差。
第二误差即为介质波速模型引起的误差。具体的,可以根据矿山岩体结构和目标监测区域结构,进行非均质速度结构等效分析。
步骤S104:基于最小震级的微震事件触发的传感器对最小震级的微震事件进行源定位,得到协方差矩阵。
在一种实施方式中,可以根据最小震级的微震事件触发的传感器序列建立求解该微震事件震源的二范数走时方程,通过目标函数的矩阵表达式,确定最小震级的微震事件定位的协方差矩阵表达式。
步骤S105:基于第一误差、第二误差和协方差矩阵,得到目标微震监测台网的总误差。
具体可以按照以下公式确定目标微震监测台网的总误差:
Figure M_220624093844705_705727001
其中,
Figure M_220624093844753_753606001
为第一误差,C为协方差矩阵,
Figure M_220624093844770_770643002
为第二误差,
Figure M_220624093844801_801900003
Figure M_220624093844817_817519004
分布
Figure M_220624093844848_848784005
分位数。
本发明实施例提供的上述微震监测台网评价方法,结合目标微震监测台网的到时拾取能力产生的第一误差以及由矿山岩体结构和目标监测区域结构产生的二误差,基于目标微震监测台网实际布设的传感器位置定位,对目标微震监测台网的性能进行评价,全面衡量了影响微震监测台网性能的误差来源,从而提高了微震监测台网监测结果准确性的判断力。
为了便于理解,本发明实施例还提供了一种基于传感器坐标和目标监测点坐标确定目标监测点监测到的微震事件的最小震级,并基于最小震级的微震事件确定第一误差的具体实现方式,即对于前述步骤S102,可以采用包括但不限于以下步骤1至步骤4实现:
步骤1:基于传感器坐标和目标监测点坐标计算震源距。
具体的,可以按照以下公式计算震源距:
Figure P_220624093844864_864395001
式中,r为震源距,(x 0,y 0,z 0)为目标监测点坐标,(x i ,y i ,z i )为第i个传感器的传感器坐标。
步骤2:基于震源距计算目标监测点监测到的微震事件的最小震级;
在具体应用中,可以按照以下公式计算目标监测点监测到的微震事件的最小震级:
Figure P_220624093844895_895642001
其中,Mx 0,y 0,z 0)为目标监测点(x 0,y 0,z 0)监测到的微震事件的最小震级,r为震源距,A、B、α为与目标监测区域的介质属性相关的常数,可以根据目标监测区域的特征进行校正,对于破碎岩体,地震波衰减速度较快,α取值应相应减小。在一种实施方式中,目标监测区域可监测到的最小震级的结果可以参见图2所示,图中横坐标表示网格点沿x轴位置,单位米,纵坐标表示网格点沿y轴位置,单位米,颜色越深表示能够观测到的最小震级越小。
步骤3:确定最小震级的微震事件触发的传感器与目标监测点之间的传感器距离。
步骤4:基于最小震级的微震事件的信噪比和传感器距离,计算第一误差。
在一种实施方式中,对于第一误差的计算可以采用包括但不限于以下方式:
首先,基于目标微震监测台网采集到的最小震级的微震事件的波形计算最小震级的微震事件的信噪比。具体的,最小震级的微震事件的信噪比S=信号能量/噪声能量;其中,信号能量即为最小震级的微震事件的能量,噪声能量即为微震监测台网记录的波形的背景噪声。
然后,基于信噪比和传感器距离,计算最小震级的微震事件触发的传感器的到时拾取误差。
最后,将到时拾取误差的标准差确定为第一误差。
在一种实施方式中,可以由信噪比S和最小震级的微震事件触发的传感器距离R确定传感器到时拾取的误差,即
Figure M_220624093844926_926874001
,其中,S 0R 0SR的设定阈值。在具体应用中,可以将到时拾取误差的标准差
Figure M_220624093845084_084635002
确定为第一误差。
进一步,本发明实施例还提供了一种基于最小震级的微震事件触发的传感器对最小震级的微震事件进行源定位,得到协方差矩阵的具体实现方式,即对于前述步骤S104,可以采用包括但不限于以下方式实现:
首先,基于最小震级的微震事件触发的传感器,计算每个被触发的传感器的地震波走时残差。
在一种实施方式中,根据最小震级的微震事件触发的传感器序列建立求解微震震源的二范数走时方程,其中,二范数走时方程的目标函数为
Figure M_220624093845100_100719001
Figure M_220624093845131_131977002
为第i个传感器的地震波走时残差;接着,对地震波走时残差
Figure M_220624093845165_165173003
的未知数震源位置和震源发生时间t 0进行一阶泰勒近似得到:
Figure M_220624093845196_196508004
,将其写为矩阵形式即为:
Figure M_220624093845258_258926005
,其中,A为系数矩阵,
Figure M_220624093845290_290184006
T i 为第i个传感器的地震波走时,
Figure M_220624093845370_370249007
为第i个传感器的地震波走时对未知数
Figure M_220624093845385_385867008
的偏导,
Figure M_220624093845417_417157009
=(x,y,z)。
然后,基于每个被触发的传感器的地震波走时残差的系数矩阵,确定协方差矩阵。
在一种具体的实施方式中,可以根据前述方程的矩阵表达形式,得到微震定位的协方差矩阵表达式,即为:C=(A’*A)-1,即协方差矩阵。
进一步,可以根据协方差矩阵确定目标微震监测台网的定位误差,参见图3所示的目标微震监测台网的定位误差示意图,图中横坐标表示网格点沿x轴位置,单位米,纵坐标表示网格点沿y轴位置,单位米,颜色越深表示定位误差越大。
本发明实施例提供的上述微震监测台网评价方法,结合微震监测台网实际布设的传感器位置、微震监测台网的到时拾取能力、矿山岩石介质波速模型误差三方面,对实际工程应用中的微震监测台网进行可测最小震级和微震定位误差分析,实现对微震监测台网评价,对微震监测台网在矿山、隧道、边坡治理等岩土工程和安全监测工程应用领域提供实际评价价值;此外,还可以根据计算得到的微震监测台网的总误差对定位结果进行修正,从而提高微震监测台网监测结果的准确性。
本发明实施例还提供了一种具体的微震监测台网评价方法,参见图4所示的另一种微震监测台网评价方法的流程图,示意出该方法主要包括以下步骤S401至步骤S408:
步骤S401:获取微震监测台网的传感器坐标。
具体的,可参见表1所示的传感器的安装实际坐标。
表1 传感器安装实际坐标
传感器编号 <i>x </i>( 单位 : 米 ) <i>y </i>( 单位 : 米 ) <i>z </i>( 单位 : 米 )
1 21133 21240 440
2 21219 21246 440
3 21301 21133 320
4 21255 21155 320
5 21208 21172 320
6 21161 21189 320
步骤S402:确定目标监测区域以及目标监测点坐标。
步骤S403:计算各传感器坐标与目标检测点坐标之间的震源距。
具体的,也即计算各个传感器到目标监测点的距离。
步骤S404:根据震级与距离经验关系,确定可监测微震事件的最小震级。
步骤S405:根据微震监测台网的噪声水平和最小震级的微震事件触发的传感器距离计算到时拾取误差。
步骤S406:根据最小震级的微震事件触发的传感器序列确定定位协方差矩阵。
具体的,根据最小震级的微震事件触发的传感器序列建立求解微震震源的二范数走时方程,对地震波走时残差
Figure M_220624093845448_448365001
的未知数震源位置和震源发生时间t 0进行一阶泰勒近似得到走时方程的矩阵形式,根据矩阵形式的系数矩阵确定定位协方差矩阵。
步骤S407:根据到时拾取误差、定位协方差矩阵、介质波速模型误差和统计检验水平误差,计算总误差。
具体的,总误差为:
Figure M_220624093845479_479696001
其中,
Figure M_220624093845510_510873001
为第一误差,C为协方差矩阵,
Figure M_220624093845545_545519002
为第二误差,
Figure M_220624093845577_577334003
Figure M_220624093845608_608519004
分布
Figure M_220624093845624_624147005
分位数。
步骤S408:重复步骤S403至步骤S407,得到每个网格化的目标监测区域的最小震级和总误差。
需要说明的是,本发明实施例所提供的方法,其实现原理及产生的技术效果和前述实施例相同,为简要描述,本实施例部分未提及之处,可参考前述实施例中相应内容。
对于前述实施例提供的微震监测台网评价方法,本发明实施例还提供了一种微震监测台网评价装置,参见图5所示的一种微震监测台网评价装置的结构示意图,该装置可以包括以下部分:
数据获取模块501,用于获取目标微震监测台网在目标监测区域内安装的传感器的传感器坐标和目标监测点坐标;
第一计算模块502,用于基于传感器坐标和目标监测点坐标确定目标监测点监测到的微震事件的最小震级,并基于最小震级的微震事件确定第一误差;
第二计算模块503,用于基于矿山岩体结构和目标监测区域结构确定第二误差;
第三计算模块504,用于基于最小震级的微震事件触发的传感器对最小震级的微震事件进行源定位,得到协方差矩阵;
第四计算模块505,用于基于第一误差、第二误差和协方差矩阵,得到目标微震监测台网的总误差。
本发明实施例提供的上述微震监测台网评价装置,结合目标微震监测台网实际布设的传感器位置的到时拾取能力产生的第一误差以及由矿山岩体结构和目标监测区域结构产生的第二误差,对目标微震监测台网的性能进行评价,全面衡量了影响微震监测台网性能的误差来源,从而提高了微震监测台网监测结果的准确性。
在一种实施方式中,上述第一计算模块502进一步还用于:基于传感器坐标和目标监测点坐标计算震源距;基于震源距计算目标监测点监测到的微震事件的最小震级;确定最小震级的微震事件触发的传感器与目标监测点之间的传感器距离;基于最小震级的微震事件的信噪比和传感器距离,计算第一误差。
在一种实施方式中,上述第一计算模块502进一步还用于:按照以下公式计算目标监测点监测到的微震事件的最小震级:
Figure P_220624093845655_655402001
其中,Mx 0,y 0,z 0)为目标监测点(x 0,y 0,z 0)监测到的微震事件的最小震级,r为震源距,A、B、α为与目标监测区域的介质属性相关的常数。
在一种实施方式中,上述第一计算模块502进一步还用于:基于目标微震监测台网采集到的最小震级的微震事件的波形计算最小震级的微震事件的信噪比;基于信噪比和传感器距离,计算最小震级的微震事件触发的传感器的到时拾取误差;将到时拾取误差的标准差确定为第一误差。
在一种实施方式中,上述第三计算模块504进一步还用于:基于最小震级的微震事件触发的传感器,计算每个被触发的传感器的地震波走时残差;基于每个被触发的传感器的地震波走时残差的系数矩阵,确定协方差矩阵。
在一种实施方式中,上述第四计算模块505进一步还用于:按照以下公式确定目标微震监测台网的总误差:
Figure M_220624093845671_671038001
其中,
Figure M_220624093845717_717946001
为第一误差,C为协方差矩阵,
Figure M_220624093845750_750588002
为第二误差,
Figure M_220624093845766_766739003
Figure M_220624093845798_798003004
分布
Figure M_220624093845829_829237005
分位数。
在一种实施方式中,上述装置还包括区域划分模块,用于对监测区域进行网格划分,得到多个目标监测区域。
本发明实施例所提供的装置,其实现原理及产生的技术效果和前述方法实施例相同,为简要描述,装置实施例部分未提及之处,可参考前述方法实施例中相应内容。
本发明实施例还提供了一种电子设备,具体的,该电子设备包括处理器和存储装置;存储装置上存储有计算机程序,计算机程序在被处理器运行时执行如上实施方式的任一项的方法。
图6为本发明实施例提供的一种电子设备的结构示意图,该电子设备100包括:处理器60,存储器61,总线62和通信接口63,处理器60、通信接口63和存储器61通过总线62连接;处理器60用于执行存储器61中存储的可执行模块,例如计算机程序。
其中,存储器61可能包含高速随机存取存储器(RAM,Random Access Memory),也可能还包括非不稳定的存储器(non-volatile memory),例如至少一个磁盘存储器。通过至少一个通信接口63(可以是有线或者无线)实现该系统网元与至少一个其他网元之间的通信连接,可以使用互联网,广域网,本地网,城域网等。
总线62可以是ISA总线、PCI总线或EISA总线等。总线可以分为地址总线、数据总线、控制总线等。为便于表示,图6中仅用一个双向箭头表示,但并不表示仅有一根总线或一种类型的总线。
其中,存储器61用于存储程序,处理器60在接收到执行指令后,执行程序,前述本发明实施例任一实施例揭示的流过程定义的装置所执行的方法可以应用于处理器60中,或者由处理器60实现。
处理器60可能是一种集成电路芯片,具有信号的处理能力。在实现过程中,上述方法的各步骤可以通过处理器60中的硬件的集成逻辑电路或者软件形式的指令完成。上述的处理器60可以是通用处理器,包括中央处理器(Central Processing Unit,简称CPU)、网络处理器(Network Processor,简称NP)等;还可以是数字信号处理器(Digital SignalProcessing,简称DSP)、专用集成电路(Application Specific Integrated Circuit,简称ASIC)、现成可编程门阵列(Field-Programmable Gate Array,简称FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行本发明实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合本发明实施例所公开的方法的步骤可以直接体现为硬件译码处理器执行完成,或者用译码处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器,闪存、只读存储器,可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器61,处理器60读取存储器61中的信息,结合其硬件完成上述方法的步骤。
本发明实施例所提供的可读存储介质的计算机程序产品,包括存储了程序代码的计算机可读存储介质,程序代码包括的指令可用于执行前面方法实施例中的方法,具体实现可参见前述方法实施例,在此不再赘述。
功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
最后应说明的是:以上实施例,仅为本发明的具体实施方式,用以说明本发明的技术方案,而非对其限制,本发明的保护范围并不局限于此,尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化,或者对其中部分技术特征进行等同替换;而这些修改、变化或者替换,并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以权利要求的保护范围为准。

Claims (10)

1.一种微震监测台网评价方法,其特征在于,包括:
获取目标微震监测台网在目标监测区域内安装的传感器的传感器坐标和目标监测点坐标;
基于所述传感器坐标和所述目标监测点坐标确定所述目标监测点监测到的微震事件的最小震级,并基于所述最小震级的微震事件确定第一误差;
基于矿山岩体结构和目标监测区域结构确定第二误差;
基于所述最小震级的微震事件触发的传感器对所述最小震级的微震事件进行源定位,得到协方差矩阵;
基于所述第一误差、所述第二误差和所述协方差矩阵,得到所述目标微震监测台网的总误差。
2.根据权利要求1所述的方法,其特征在于,基于所述传感器坐标和所述目标监测点坐标确定所述目标监测点监测到的微震事件的最小震级,并基于所述最小震级的微震事件确定第一误差的步骤,包括:
基于所述传感器坐标和所述目标监测点坐标计算震源距;
基于所述震源距计算所述目标监测点监测到的微震事件的最小震级;
确定所述最小震级的微震事件触发的传感器与所述目标监测点之间的传感器距离;
基于所述最小震级的微震事件的信噪比和所述传感器距离,计算第一误差。
3.根据权利要求2所述的方法,其特征在于,基于所述震源距计算所述目标监测点监测到的微震事件的最小震级的步骤,包括:
按照以下公式计算所述目标监测点监测到的微震事件的最小震级:
Figure P_220624093841318_318014001
其中,Mx 0,y 0,z 0)为目标监测点(x 0,y 0,z 0)监测到的微震事件的最小震级,r为震源距,A、B、α为与目标监测区域的介质属性相关的常数。
4.根据权利要求2所述的方法,其特征在于,基于所述最小震级的微震事件的信噪比和所述传感器距离,计算第一误差的步骤,包括:
基于所述目标微震监测台网采集到的所述最小震级的微震事件的波形计算所述最小震级的微震事件的信噪比;
基于所述信噪比和所述传感器距离,计算所述最小震级的微震事件触发的传感器的到时拾取误差;
将所述到时拾取误差的标准差确定为第一误差。
5.根据权利要求1所述的方法,其特征在于,基于所述最小震级的微震事件触发的传感器对所述最小震级的微震事件进行源定位,得到协方差矩阵的步骤,包括:
基于最小震级的微震事件触发的传感器,计算每个被触发的传感器的地震波走时残差;
基于每个被触发的传感器的地震波走时残差的系数矩阵,确定协方差矩阵。
6.根据权利要求1所述的方法,其特征在于,基于所述第一误差、所述第二误差和所述协方差矩阵,得到所述目标微震监测台网的总误差的步骤,包括:
按照以下公式确定所述目标微震监测台网的总误差:
Figure M_220624093841478_478638001
其中,
Figure M_220624093841623_623659001
为第一误差,C为协方差矩阵,
Figure M_220624093841654_654930002
为第二误差,
Figure M_220624093841686_686182003
Figure M_220624093841717_717429004
分布
Figure M_220624093841750_750106005
分位数。
7.根据权利要求1所述的方法,其特征在于,所述方法还包括:
对监测区域进行网格划分,得到多个目标监测区域。
8.一种微震监测台网评价装置,其特征在于,包括:
数据获取模块,用于获取目标微震监测台网在目标监测区域内安装的传感器的传感器坐标和目标监测点坐标;
第一计算模块,用于基于所述传感器坐标和所述目标监测点坐标确定所述目标监测点监测到的微震事件的最小震级,并基于所述最小震级的微震事件确定第一误差;
第二计算模块,用于基于矿山岩体结构和目标监测区域结构确定第二误差;
第三计算模块,用于基于所述最小震级的微震事件触发的传感器对所述最小震级的微震事件进行源定位,得到协方差矩阵;
第四计算模块,用于基于所述第一误差、所述第二误差和所述协方差矩阵,得到所述目标微震监测台网的总误差。
9.一种电子设备,其特征在于,包括处理器和存储器,所述存储器存储有能够被所述处理器执行的计算机可执行指令,所述处理器执行所述计算机可执行指令以实现权利要求1至7任一项所述的方法的步骤。
10.一种计算机可读存储介质,所述计算机可读存储介质上存储有计算机程序,其特征在于,所述计算机程序被处理器运行时执行上述权利要求1至7任一项所述的方法的步骤。
CN202210737968.7A 2022-06-28 2022-06-28 微震监测台网评价方法、装置、电子设备及介质 Active CN114814940B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210737968.7A CN114814940B (zh) 2022-06-28 2022-06-28 微震监测台网评价方法、装置、电子设备及介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210737968.7A CN114814940B (zh) 2022-06-28 2022-06-28 微震监测台网评价方法、装置、电子设备及介质

Publications (2)

Publication Number Publication Date
CN114814940A true CN114814940A (zh) 2022-07-29
CN114814940B CN114814940B (zh) 2022-09-09

Family

ID=82522341

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210737968.7A Active CN114814940B (zh) 2022-06-28 2022-06-28 微震监测台网评价方法、装置、电子设备及介质

Country Status (1)

Country Link
CN (1) CN114814940B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115826037A (zh) * 2022-11-02 2023-03-21 北京国信安科技术有限公司 矿山采空区微震监测预警系统监测预警能力的评估方法
CN116842411A (zh) * 2023-06-12 2023-10-03 中国矿业大学 基于震动波探测概率的煤矿微震台网监测能力评价方法

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010276383A (ja) * 2009-05-26 2010-12-09 Sumitomo Electric System Solutions Co Ltd 緊急地震速報受信装置
CN105093297A (zh) * 2015-07-10 2015-11-25 中联煤层气有限责任公司 一种微地震定位精度的质量控制方法
CN105785436A (zh) * 2016-03-17 2016-07-20 北京矿冶研究总院 一种矿用微震监测方法
CN106199703A (zh) * 2016-08-26 2016-12-07 中国矿业大学 一种微震震源自动定位及可靠性综合评价方法
CN107861157A (zh) * 2017-10-27 2018-03-30 福州大学 一种地下水封洞库运营期微震监测方法
US9945970B1 (en) * 2011-08-29 2018-04-17 Seismic Innovations Method and apparatus for modeling microseismic event location estimate accuracy
CN108490485A (zh) * 2018-06-20 2018-09-04 成都理工大学 双线隧道微震监测系统及其定位精度评估方法
CN110987347A (zh) * 2019-12-09 2020-04-10 大连理工大学 开挖扰动作用下抽水蓄能电站厂房围岩稳定性判断方法
CN111308559A (zh) * 2020-03-06 2020-06-19 西南大学 基于p波入射方向的煤矿井下微震监测传感器布设方法
CN111458746A (zh) * 2020-05-29 2020-07-28 东北大学 一种基于U-Net神经网络的隧道微震波形到时拾取方法
CN113655543A (zh) * 2021-08-18 2021-11-16 西北核技术研究所 一种地震台网定位能力评估方法
CN113885074A (zh) * 2021-09-09 2022-01-04 中煤科工开采研究院有限公司 微震监测能量校核方法

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010276383A (ja) * 2009-05-26 2010-12-09 Sumitomo Electric System Solutions Co Ltd 緊急地震速報受信装置
US9945970B1 (en) * 2011-08-29 2018-04-17 Seismic Innovations Method and apparatus for modeling microseismic event location estimate accuracy
CN105093297A (zh) * 2015-07-10 2015-11-25 中联煤层气有限责任公司 一种微地震定位精度的质量控制方法
CN105785436A (zh) * 2016-03-17 2016-07-20 北京矿冶研究总院 一种矿用微震监测方法
CN106199703A (zh) * 2016-08-26 2016-12-07 中国矿业大学 一种微震震源自动定位及可靠性综合评价方法
CN107861157A (zh) * 2017-10-27 2018-03-30 福州大学 一种地下水封洞库运营期微震监测方法
CN108490485A (zh) * 2018-06-20 2018-09-04 成都理工大学 双线隧道微震监测系统及其定位精度评估方法
CN110987347A (zh) * 2019-12-09 2020-04-10 大连理工大学 开挖扰动作用下抽水蓄能电站厂房围岩稳定性判断方法
CN111308559A (zh) * 2020-03-06 2020-06-19 西南大学 基于p波入射方向的煤矿井下微震监测传感器布设方法
CN111458746A (zh) * 2020-05-29 2020-07-28 东北大学 一种基于U-Net神经网络的隧道微震波形到时拾取方法
CN113655543A (zh) * 2021-08-18 2021-11-16 西北核技术研究所 一种地震台网定位能力评估方法
CN113885074A (zh) * 2021-09-09 2022-01-04 中煤科工开采研究院有限公司 微震监测能量校核方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
LI NAN等: "The Influence Mechanism and Optimization of the Sensor Network on the MS/AE Source Location", 《SHOCK AND VIBRATION》 *
WUESTEFELD, A等: "Benchmarking earthquake location algorithms: A synthetic comparison", 《GEOPHYSICS》 *
刘晓明等: "基于综合评价法的矿山微震监测台网布设评价", 《中国安全生产科学技术》 *
唐礼忠等: "大规模深井开采微震监测系统站网布置优化", 《岩石力学与工程学报》 *
张达等: "BSN矿山微震监测技术及其应用", 《中国地震》 *
朱自强等: "基于阵列分析的传感器布设效果研究", 《矿业研究与开发》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115826037A (zh) * 2022-11-02 2023-03-21 北京国信安科技术有限公司 矿山采空区微震监测预警系统监测预警能力的评估方法
CN115826037B (zh) * 2022-11-02 2023-12-12 北京国信安科技术有限公司 矿山采空区微震监测预警系统监测预警能力的评估方法
CN116842411A (zh) * 2023-06-12 2023-10-03 中国矿业大学 基于震动波探测概率的煤矿微震台网监测能力评价方法

Also Published As

Publication number Publication date
CN114814940B (zh) 2022-09-09

Similar Documents

Publication Publication Date Title
CN114814940B (zh) 微震监测台网评价方法、装置、电子设备及介质
US10877084B2 (en) Nonlinear model transformation solving and optimization method for partial discharge positioning based on multi-ultrasonic sensor
Tsogka et al. The stretching method for vibration‐based structural health monitoring of civil structures
CN109828302B (zh) 一种基于多振动传感器的震源定位方法及装置
CN104537251A (zh) 一种风机叶片的冲击载荷识别方法
CN115239105A (zh) 一种在役风电场的风资源的评估方法以及装置
He et al. A fast acoustic emission beamforming localization method based on Hilbert curve
CN114089033A (zh) 一种基于频谱分析的异常信号检测方法及系统
CN110398775B (zh) 隧道突涌水灾害微震事件信号波动初至拾取方法及系统
JP2018179580A (ja) 中性子強度モニタリングシステムおよび方法
CN117590172A (zh) 应用于变压器的局部放电声电联合定位方法、装置和设备
CN115343757A (zh) 地压风险预警方法、装置、电子设备及介质
CN114518595B (zh) 微震震源定位方法、装置、电子设备及介质
CN115356767A (zh) 基于区域分割提高微震定位与波速反演精度的方法及装置
CN113919500B (zh) 微震事件时空分形维数计算方法、装置及电子设备
CN109682457A (zh) 基于角度检测的光纤事件识别方法
CN111929730B (zh) 小尺度地质异常体检测方法及装置
Hu et al. Damage detection in plate structures based on space-time autoregressive moving average processes
JP6107967B2 (ja) 分析装置、分析方法および分析プログラム
CN104406681B (zh) 一种实时确定微震波速的测试方法
Bramley et al. On the measurement of safe fault failure rates in high-performance compute processors
Haiden et al. Statistical evaluation of ECMWF extreme wind forecasts
CN114185084B (zh) 拉底作业影响范围的确定方法、装置、电子设备及介质
CN110849971B (zh) 基于双指数窗函数法的结构模态参数识别方法
CN108844616A (zh) 一种次声事件检测方法及系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant