CN114814940A - 微震监测台网评价方法、装置、电子设备及介质 - Google Patents
微震监测台网评价方法、装置、电子设备及介质 Download PDFInfo
- 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
Links
- 238000012544 monitoring process Methods 0.000 title claims abstract description 202
- 238000011156 evaluation Methods 0.000 title claims description 23
- 239000011159 matrix material Substances 0.000 claims abstract description 48
- 230000001960 triggered effect Effects 0.000 claims abstract description 38
- 238000000034 method Methods 0.000 claims abstract description 36
- 239000011435 rock Substances 0.000 claims abstract description 14
- 238000004364 calculation method Methods 0.000 claims description 13
- 239000000126 substance Substances 0.000 claims description 8
- 238000004590 computer program Methods 0.000 claims description 6
- 238000009434 installation Methods 0.000 abstract description 7
- 238000010586 diagram Methods 0.000 description 7
- 238000004891 communication Methods 0.000 description 4
- 241000209094 Oryza Species 0.000 description 3
- 235000007164 Oryza sativa Nutrition 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 235000009566 rice Nutrition 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000004807 localization Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000009412 basement excavation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000000802 evaporation-induced self-assembly Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000000528 statistical test Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/288—Event detection in seismic signals, e.g. microseismics
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
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
技术领域
本发明涉及微震定位技术领域,尤其是涉及一种微震监测台网评价方法、装置、电子设备及介质。
背景技术
微震监测系统是一种覆盖面广、可实时有效监测深井矿山、隧道掘进、边坡治理等工程的技术手段。技术人员会在台网设计过程中,对微震监测台网的性能进行评估以确定最优传感器安装布置方案,但是,受施工条件限制和工作状态影响,在实际工程应用中,传感器安装点会偏离微震监测系统的最优传感器安装布置方案,使得微震监测台网的监测结果产生偏差。因此,需要对目标区域的微震监测台网的性能进行再次评估。
发明内容
有鉴于此,本发明的目的在于提供一种微震监测台网评价方法、装置、电子设备及介质,能够基于实际传感器的安装位置对微震监测台网的性能进行评估,从而提高微震监测台网监测结果准确性的判断力。
为了实现上述目的,本发明实施例采用的技术方案如下:
第一方面,本发明实施例提供了一种微震监测台网评价方法,包括:获取目标微震监测台网在目标监测区域内安装的传感器的传感器坐标和目标监测点坐标;基于传感器坐标和目标监测点坐标确定目标监测点监测到的微震事件的最小震级,并基于最小震级的微震事件确定第一误差;基于矿山岩体结构和目标监测区域结构确定第二误差;基于最小震级的微震事件触发的传感器对最小震级的微震事件进行源定位,得到协方差矩阵;基于第一误差、第二误差和协方差矩阵,得到目标微震监测台网的总误差。
在一种实施方式中,基于传感器坐标和目标监测点坐标确定目标监测点监测到的微震事件的最小震级,并基于最小震级的微震事件确定第一误差的步骤,包括:基于传感器坐标和目标监测点坐标计算震源距;基于震源距计算目标监测点监测到的微震事件的最小震级;确定最小震级的微震事件触发的传感器与目标监测点之间的传感器距离;基于最小震级的微震事件的信噪比和传感器距离,计算第一误差。
在一种实施方式中,基于震源距计算目标监测点监测到的微震事件的最小震级的步骤,包括:按照以下公式计算目标监测点监测到的微震事件的最小震级:
其中,M(x 0,y 0,z 0)为目标监测点(x 0,y 0,z 0)监测到的微震事件的最小震级,r为震源距,A、B、α为与目标监测区域的介质属性相关的常数。
在一种实施方式中,基于最小震级的微震事件的信噪比和传感器距离,计算第一误差的步骤,包括:基于目标微震监测台网采集到的最小震级的微震事件的波形计算最小震级的微震事件的信噪比;基于信噪比和传感器距离,计算最小震级的微震事件触发的传感器的到时拾取误差;将到时拾取误差的标准差确定为第一误差。
在一种实施方式中,基于最小震级的微震事件触发的传感器对最小震级的微震事件进行源定位,得到协方差矩阵的步骤,包括:基于最小震级的微震事件触发的传感器,计算每个被触发的传感器的地震波走时残差;基于每个被触发的传感器的地震波走时残差的系数矩阵,确定协方差矩阵。
在一种实施方式中,基于第一误差、第二误差和协方差矩阵,得到目标微震监测台网的总误差的步骤,包括:按照以下公式确定目标微震监测台网的总误差:
在一种实施方式中,上述方法还包括:对监测区域进行网格划分,得到多个目标监测区域。
第二方面,本发明实施例提供了一种微震监测台网评价装置,包括:数据获取模块,用于获取目标微震监测台网在目标监测区域内安装的传感器的传感器坐标和目标监测点坐标;第一计算模块,用于基于传感器坐标和目标监测点坐标确定目标监测点监测到的微震事件的最小震级,并基于最小震级的微震事件确定第一误差;第二计算模块,用于基于矿山岩体结构和目标监测区域结构确定第二误差;第三计算模块,用于基于最小震级的微震事件触发的传感器对最小震级的微震事件进行源定位,得到协方差矩阵;第四计算模块,用于基于第一误差、第二误差和协方差矩阵,得到目标微震监测台网的总误差。
第三方面,本发明实施例提供了一种电子设备,包括处理器和存储器,存储器存储有能够被处理器执行的计算机可执行指令,处理器执行计算机可执行指令以实现上述第一方面提供的任一项的方法的步骤。
第四方面,本发明实施例提供了一种计算机可读存储介质,计算机可读存储介质上存储有计算机程序,计算机程序被处理器运行时执行上述第一方面提供的任一项的方法的步骤。
本发明实施例带来了以下有益效果:
本发明实施例提供的上述微震监测台网评价方法、装置、电子设备及介质,首先,获取目标微震监测台网在目标监测区域内安装的传感器的传感器坐标和目标监测点坐标;然后,基于传感器坐标和目标监测点坐标确定目标监测点监测到的微震事件的最小震级,并基于最小震级的微震事件确定第一误差;接着,基于矿山岩体结构和目标监测区域结构确定第二误差;之后,基于最小震级的微震事件触发的传感器对最小震级的微震事件进行源定位,得到协方差矩阵;最后,基于第一误差、第二误差和协方差矩阵,得到目标微震监测台网的总误差。上述方法结合目标微震监测台网的到时拾取能力产生的第一误差和由矿山岩体结构和目标监测区域结构产生的第二误差,对目标微震监测台网的性能进行评价,全面衡量了影响微震监测台网性能的误差来源,从而提高了微震监测台网监测结果的准确性。
本发明的其他特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
为使本发明的上述目的、特征和优点能更明显易懂,下文特举较佳实施例,并配合所附附图,作详细说明如下。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种微震监测台网评价方法的流程图;
图2为本发明实施例提供的一种目标监测区域可监测到的最小震级示意图;
图3为本发明实施例提供的一种目标微震监测台网的定位误差示意图;
图4为本发明实施例提供的另一种微震监测台网评价方法的流程图;
图5为本发明实施例提供的一种微震监测台网评价装置的结构示意图;
图6为本发明实施例提供的一种电子设备的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
目前,受施工条件限制和工作状态影响,在实际工程应用中,传感器安装点会偏离微震监测系统的最优传感器安装布置方案,使得微震监测台网的监测结果产生偏差。因此,需要对目标区域的微震监测台网的性能进行再次评估。
基于此,本发明实施例提供的一种微震监测台网评价方法、装置、电子设备及介质,能够基于实际传感器的安装位置对微震监测台网的性能进行评估,从而提高微震监测台网监测结果的准确性。
为便于对本实施例进行理解,首先对本发明实施例所公开的一种微震监测台网评价方法进行详细介绍,该方法可以由电子设备执行,诸如:电脑、智能手机、平板电脑等。参见图1所示的一种微震监测台网评价方法的流程图,示意出该方法主要包括以下步骤S101至步骤S105:
步骤S101:获取目标微震监测台网在目标监测区域内安装的传感器的传感器坐标和目标监测点坐标。
在一种实施方式中,可以先对监测区域进行网格划分,得到多个目标监测区域,然后对每个目标监测区域内的微震监测台网的性能进行评价。在具体应用中,由于本发明实施例提供的方法是对已经安装完成的微震监测台网的性能进行评价,因此,可以获取现场微震监测台网工程安装数据,从而得到目标微震监测台网在目标监测区域内安装的传感器的传感器坐标,同时根据实际需求确定目标监测点坐标。
步骤S102:基于传感器坐标和目标监测点坐标确定目标监测点监测到的微震事件的最小震级,并基于最小震级的微震事件确定第一误差。
在一种实施方式中,微震监测台网监测到的有效触发的微震事件需要触发四个以上的传感器,基于此,本发明实施例中可以根据传感器与目标监测点之间的震源距,结合地震波震级与传播距离经验关系计算目标监测点可以监测到的最小震级的微震事件,并根据监测到的最小震级以及目标微震监测台网记录的波形确定微震事件的到时拾取误差,即第一误差。
步骤S103:基于矿山岩体结构和目标监测区域结构确定第二误差。
第二误差即为介质波速模型引起的误差。具体的,可以根据矿山岩体结构和目标监测区域结构,进行非均质速度结构等效分析。
步骤S104:基于最小震级的微震事件触发的传感器对最小震级的微震事件进行源定位,得到协方差矩阵。
在一种实施方式中,可以根据最小震级的微震事件触发的传感器序列建立求解该微震事件震源的二范数走时方程,通过目标函数的矩阵表达式,确定最小震级的微震事件定位的协方差矩阵表达式。
步骤S105:基于第一误差、第二误差和协方差矩阵,得到目标微震监测台网的总误差。
具体可以按照以下公式确定目标微震监测台网的总误差:
本发明实施例提供的上述微震监测台网评价方法,结合目标微震监测台网的到时拾取能力产生的第一误差以及由矿山岩体结构和目标监测区域结构产生的二误差,基于目标微震监测台网实际布设的传感器位置定位,对目标微震监测台网的性能进行评价,全面衡量了影响微震监测台网性能的误差来源,从而提高了微震监测台网监测结果准确性的判断力。
为了便于理解,本发明实施例还提供了一种基于传感器坐标和目标监测点坐标确定目标监测点监测到的微震事件的最小震级,并基于最小震级的微震事件确定第一误差的具体实现方式,即对于前述步骤S102,可以采用包括但不限于以下步骤1至步骤4实现:
步骤1:基于传感器坐标和目标监测点坐标计算震源距。
具体的,可以按照以下公式计算震源距:
式中,r为震源距,(x 0,y 0,z 0)为目标监测点坐标,(x i ,y i ,z i )为第i个传感器的传感器坐标。
步骤2:基于震源距计算目标监测点监测到的微震事件的最小震级;
在具体应用中,可以按照以下公式计算目标监测点监测到的微震事件的最小震级:
其中,M(x 0,y 0,z 0)为目标监测点(x 0,y 0,z 0)监测到的微震事件的最小震级,r为震源距,A、B、α为与目标监测区域的介质属性相关的常数,可以根据目标监测区域的特征进行校正,对于破碎岩体,地震波衰减速度较快,α取值应相应减小。在一种实施方式中,目标监测区域可监测到的最小震级的结果可以参见图2所示,图中横坐标表示网格点沿x轴位置,单位米,纵坐标表示网格点沿y轴位置,单位米,颜色越深表示能够观测到的最小震级越小。
步骤3:确定最小震级的微震事件触发的传感器与目标监测点之间的传感器距离。
步骤4:基于最小震级的微震事件的信噪比和传感器距离,计算第一误差。
在一种实施方式中,对于第一误差的计算可以采用包括但不限于以下方式:
首先,基于目标微震监测台网采集到的最小震级的微震事件的波形计算最小震级的微震事件的信噪比。具体的,最小震级的微震事件的信噪比S=信号能量/噪声能量;其中,信号能量即为最小震级的微震事件的能量,噪声能量即为微震监测台网记录的波形的背景噪声。
然后,基于信噪比和传感器距离,计算最小震级的微震事件触发的传感器的到时拾取误差。
最后,将到时拾取误差的标准差确定为第一误差。
进一步,本发明实施例还提供了一种基于最小震级的微震事件触发的传感器对最小震级的微震事件进行源定位,得到协方差矩阵的具体实现方式,即对于前述步骤S104,可以采用包括但不限于以下方式实现:
首先,基于最小震级的微震事件触发的传感器,计算每个被触发的传感器的地震波走时残差。
在一种实施方式中,根据最小震级的微震事件触发的传感器序列建立求解微震震源的二范数走时方程,其中,二范数走时方程的目标函数为,为第i个传感器的地震波走时残差;接着,对地震波走时残差的未知数震源位置和震源发生时间t 0进行一阶泰勒近似得到:,将其写为矩阵形式即为:,其中,A为系数矩阵,,T i 为第i个传感器的地震波走时,为第i个传感器的地震波走时对未知数的偏导,=(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:根据最小震级的微震事件触发的传感器序列确定定位协方差矩阵。
具体的,根据最小震级的微震事件触发的传感器序列建立求解微震震源的二范数走时方程,对地震波走时残差的未知数震源位置和震源发生时间t 0进行一阶泰勒近似得到走时方程的矩阵形式,根据矩阵形式的系数矩阵确定定位协方差矩阵。
步骤S407:根据到时拾取误差、定位协方差矩阵、介质波速模型误差和统计检验水平误差,计算总误差。
具体的,总误差为:
步骤S408:重复步骤S403至步骤S407,得到每个网格化的目标监测区域的最小震级和总误差。
需要说明的是,本发明实施例所提供的方法,其实现原理及产生的技术效果和前述实施例相同,为简要描述,本实施例部分未提及之处,可参考前述实施例中相应内容。
对于前述实施例提供的微震监测台网评价方法,本发明实施例还提供了一种微震监测台网评价装置,参见图5所示的一种微震监测台网评价装置的结构示意图,该装置可以包括以下部分:
数据获取模块501,用于获取目标微震监测台网在目标监测区域内安装的传感器的传感器坐标和目标监测点坐标;
第一计算模块502,用于基于传感器坐标和目标监测点坐标确定目标监测点监测到的微震事件的最小震级,并基于最小震级的微震事件确定第一误差;
第二计算模块503,用于基于矿山岩体结构和目标监测区域结构确定第二误差;
第三计算模块504,用于基于最小震级的微震事件触发的传感器对最小震级的微震事件进行源定位,得到协方差矩阵;
第四计算模块505,用于基于第一误差、第二误差和协方差矩阵,得到目标微震监测台网的总误差。
本发明实施例提供的上述微震监测台网评价装置,结合目标微震监测台网实际布设的传感器位置的到时拾取能力产生的第一误差以及由矿山岩体结构和目标监测区域结构产生的第二误差,对目标微震监测台网的性能进行评价,全面衡量了影响微震监测台网性能的误差来源,从而提高了微震监测台网监测结果的准确性。
在一种实施方式中,上述第一计算模块502进一步还用于:基于传感器坐标和目标监测点坐标计算震源距;基于震源距计算目标监测点监测到的微震事件的最小震级;确定最小震级的微震事件触发的传感器与目标监测点之间的传感器距离;基于最小震级的微震事件的信噪比和传感器距离,计算第一误差。
在一种实施方式中,上述第一计算模块502进一步还用于:按照以下公式计算目标监测点监测到的微震事件的最小震级:
其中,M(x 0,y 0,z 0)为目标监测点(x 0,y 0,z 0)监测到的微震事件的最小震级,r为震源距,A、B、α为与目标监测区域的介质属性相关的常数。
在一种实施方式中,上述第一计算模块502进一步还用于:基于目标微震监测台网采集到的最小震级的微震事件的波形计算最小震级的微震事件的信噪比;基于信噪比和传感器距离,计算最小震级的微震事件触发的传感器的到时拾取误差;将到时拾取误差的标准差确定为第一误差。
在一种实施方式中,上述第三计算模块504进一步还用于:基于最小震级的微震事件触发的传感器,计算每个被触发的传感器的地震波走时残差;基于每个被触发的传感器的地震波走时残差的系数矩阵,确定协方差矩阵。
在一种实施方式中,上述第四计算模块505进一步还用于:按照以下公式确定目标微震监测台网的总误差:
在一种实施方式中,上述装置还包括区域划分模块,用于对监测区域进行网格划分,得到多个目标监测区域。
本发明实施例所提供的装置,其实现原理及产生的技术效果和前述方法实施例相同,为简要描述,装置实施例部分未提及之处,可参考前述方法实施例中相应内容。
本发明实施例还提供了一种电子设备,具体的,该电子设备包括处理器和存储装置;存储装置上存储有计算机程序,计算机程序在被处理器运行时执行如上实施方式的任一项的方法。
图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所述的方法,其特征在于,基于所述传感器坐标和所述目标监测点坐标确定所述目标监测点监测到的微震事件的最小震级,并基于所述最小震级的微震事件确定第一误差的步骤,包括:
基于所述传感器坐标和所述目标监测点坐标计算震源距;
基于所述震源距计算所述目标监测点监测到的微震事件的最小震级;
确定所述最小震级的微震事件触发的传感器与所述目标监测点之间的传感器距离;
基于所述最小震级的微震事件的信噪比和所述传感器距离,计算第一误差。
4.根据权利要求2所述的方法,其特征在于,基于所述最小震级的微震事件的信噪比和所述传感器距离,计算第一误差的步骤,包括:
基于所述目标微震监测台网采集到的所述最小震级的微震事件的波形计算所述最小震级的微震事件的信噪比;
基于所述信噪比和所述传感器距离,计算所述最小震级的微震事件触发的传感器的到时拾取误差;
将所述到时拾取误差的标准差确定为第一误差。
5.根据权利要求1所述的方法,其特征在于,基于所述最小震级的微震事件触发的传感器对所述最小震级的微震事件进行源定位,得到协方差矩阵的步骤,包括:
基于最小震级的微震事件触发的传感器,计算每个被触发的传感器的地震波走时残差;
基于每个被触发的传感器的地震波走时残差的系数矩阵,确定协方差矩阵。
7.根据权利要求1所述的方法,其特征在于,所述方法还包括:
对监测区域进行网格划分,得到多个目标监测区域。
8.一种微震监测台网评价装置,其特征在于,包括:
数据获取模块,用于获取目标微震监测台网在目标监测区域内安装的传感器的传感器坐标和目标监测点坐标;
第一计算模块,用于基于所述传感器坐标和所述目标监测点坐标确定所述目标监测点监测到的微震事件的最小震级,并基于所述最小震级的微震事件确定第一误差;
第二计算模块,用于基于矿山岩体结构和目标监测区域结构确定第二误差;
第三计算模块,用于基于所述最小震级的微震事件触发的传感器对所述最小震级的微震事件进行源定位,得到协方差矩阵;
第四计算模块,用于基于所述第一误差、所述第二误差和所述协方差矩阵,得到所述目标微震监测台网的总误差。
9.一种电子设备,其特征在于,包括处理器和存储器,所述存储器存储有能够被所述处理器执行的计算机可执行指令,所述处理器执行所述计算机可执行指令以实现权利要求1至7任一项所述的方法的步骤。
10.一种计算机可读存储介质,所述计算机可读存储介质上存储有计算机程序,其特征在于,所述计算机程序被处理器运行时执行上述权利要求1至7任一项所述的方法的步骤。
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)
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)
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 | 中煤科工开采研究院有限公司 | 微震监测能量校核方法 |
-
2022
- 2022-06-28 CN CN202210737968.7A patent/CN114814940B/zh active Active
Patent Citations (12)
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)
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)
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 |