CN113642865A - 一种矿山风险等级的确定方法、装置及服务器 - Google Patents

一种矿山风险等级的确定方法、装置及服务器 Download PDF

Info

Publication number
CN113642865A
CN113642865A CN202110870685.5A CN202110870685A CN113642865A CN 113642865 A CN113642865 A CN 113642865A CN 202110870685 A CN202110870685 A CN 202110870685A CN 113642865 A CN113642865 A CN 113642865A
Authority
CN
China
Prior art keywords
microseismic
determining
risk
target
magnitude
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.)
Pending
Application number
CN202110870685.5A
Other languages
English (en)
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.)
Fankou Lead Zinc Mine of Shenzhen Zhongjin Lingnan Nonfemet Co Ltd
Original Assignee
Fankou Lead Zinc Mine of Shenzhen Zhongjin Lingnan Nonfemet 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 Fankou Lead Zinc Mine of Shenzhen Zhongjin Lingnan Nonfemet Co Ltd filed Critical Fankou Lead Zinc Mine of Shenzhen Zhongjin Lingnan Nonfemet Co Ltd
Priority to CN202110870685.5A priority Critical patent/CN113642865A/zh
Publication of CN113642865A publication Critical patent/CN113642865A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0635—Risk analysis of enterprise or organisation activities
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/02—Agriculture; Fishing; Mining

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Educational Administration (AREA)
  • Development Economics (AREA)
  • General Physics & Mathematics (AREA)
  • Marketing (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Animal Husbandry (AREA)
  • Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Mining & Mineral Resources (AREA)
  • Agronomy & Crop Science (AREA)
  • Alarm Systems (AREA)

Abstract

本申请适用于矿山安全技术领域，提供了一种矿山风险等级的确定方法、装置及服务器，所述方法包括：获取监测区域在预设时间段内的多次微震的震级和微震指标信息；根据每次微震的震级和微震指标信息确定所述监测区域的微震特征信息；将微震指标信息和微震特征信息导入预设的风险预测模型，计算得到监测区域的目标风险系数；根据目标风险系数以及预设风险转换关系确定监测区域的目标风险等级。本申请提供的矿山风险等级的确定方法可以准确地确定矿山的风险等级，提高了对矿山风险的监测准确率。

Description

一种矿山风险等级的确定方法、装置及服务器
技术领域
本申请属于矿山安全技术领域，尤其涉及一种矿山风险等级的确定方法、矿山风险等级的确定装置、服务器及计算机可读存储介质。
背景技术
目前，矿山开采过程中形成的巷道和采空区，破坏了矿山区域原岩应力平衡。采场、围岩、矿体内应力重新分布，导致矿柱、工作面顶板和围岩发生位移和变形。同时，由于开采扰动会大范围改变岩体的应力分布，使得应力集中区域的岩体破裂、岩层急剧弹性能瞬间释放形成矿震(矿山微地震)，从而导致顶板冒落、矿柱压裂和围岩开裂等地压灾害，危及人身安全。因此，需要对矿山的矿震风险等级进行监测。
现有技术通常是通过微震监测技术对矿山地压的矿震风险等级进行预测。然而，现有技术通常只是通过是否发生大事件、过去一段时间的b值(某段时间段内小震数与大震数的比值)2个方面的评分来评估矿震风险，无法准确地确定矿山的矿震风险等级。
发明内容
本申请实施例提供了一种矿山风险等级的确定方法、矿山风险等级的确定装置、服务器及计算机可读存储介质，可以解决现有技术无法准确地确定矿山的风险等级的问题。
第一方面，本申请实施例提供了一种矿山风险等级的确定方法，包括：
获取监测区域在预设时间段内的多次微震的震级和微震指标信息；
根据每次微震的震级和微震指标信息确定所述监测区域的微震特征信息；
将所述微震指标信息和所述微震特征信息导入预设的风险预测模型，计算得到所述监测区域的目标风险系数；
根据所述目标风险系数以及预设风险转换关系确定所述监测区域的目标风险等级。
可选的，所述获取监测区域在预设时间段内的多次微震的震级和微震指标信息，包括：
获取所述监测区域在所述预设时间段内的多种微震参数；所述多种微震参数包括视体积、能量指数、能量释放量和/或施密特数；
针对任一微震参数，生成所述微震参数随时间变化的变化曲线；
根据所述多种微震参数的变化曲线确定所述微震指标信息。
可选的，所述微震特征信息包括第一微震特征信息和第二微震特征信息，所述根据每次微震的震级和微震指标信息确定所述监测区域的微震特征信息，包括：
获取所述监测区域在所述预设时间段内每天发生的微震的次数；
根据每次所述微震的震级确定所述预设时间段内大小震级的比值；
根据所述大小震级的比值和所述每天发生的微震的次数确定所述第一微震特征信息；
根据目标次数和所述微震指标信息确定所述第二微震特征信息；所述目标次数指发生的震级大于第一预设等级的微震事件的次数。
可选的，所述述目标风险系数根据以下公式计算得到：
F＝w1*f1+w2*f2+w3*f3；
其中，F表示所述目标风险系数，f1表示所述微震指标信息对应的第一风险因子，f2表示所述第一微震特征信息对应的第二风险因子，f3表示所述第二微震特征信息对应的第三风险因子，w1表示所述第一风险因子对应的权重系数，w2表示所述第二风险因子对应的权重系数，w3表示所述第三风险因子对应的权重系数。
可选的，所述将所述微震指标信息和所述微震特征信息导入预设的风险预测模型，计算得到所述监测区域的目标风险系数之前，还包括：
根据预设权重计算方法确定所述第一风险因子、第二风险因子及第三风险因子分别对应的权重系数。
可选的，所述目标风险等级根据以下公式计算得到：
nr＝A*nf；
其中，nr表示所述目标风险等级，A表示转换因子，表示所述目标风险系数nf。
可选的，所述根据所述目标风险系数确定所述监测区域的目标风险等级之后，还包括：
根据所述目标风险等级以及预设风险等级与预设策略之间的对应关系，确定所述目标风险等级对应的目标策略，并执行所述目标策略。
第二方面，本申请实施例提供了一种矿山风险等级的确定装置，包括：
第一获取单元，用于获取监测区域在预设时间段内的多次微震的震级和微震指标信息；
第一确定单元，用于根据每次微震的震级和微震指标信息确定所述监测区域的微震特征信息；
计算单元，用于将所述微震指标信息和所述微震特征信息导入预设的风险预测模型，计算得到所述监测区域的目标风险系数；
第二确定单元，用于根据所述目标风险系数以及预设风险转换关系确定所述监测区域的目标风险等级。
第三方面，本申请实施例提供了一种服务器，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，上述处理器执行所述计算机程序时实现如上述第一方面中任一项所述的矿山风险等级的确定方法的步骤。
第四方面，本申请实施例提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现如上述第一方面中任一项所述的矿山风险等级的确定方法的步骤。
第五方面，本申请实施例提供了一种计算机程序产品，当计算机程序产品在服务器上运行时，使得服务器可执行上述第一方面中任一项所述的矿山风险等级的确定方法。
本申请实施例与现有技术相比存在的有益效果是：
本申请实施例提供的一种矿山风险等级的确定方法，通过获取到的监测区域在预设时间段内的所有微震的震级和微震指标信息，可以确定监测区域的第一微震特征信息和第二微震特征信息，并将上述微震指标信息、第一微震特征信息及第二微震特征信息导入预设的风险预测模型，可以计算得到监测区域的目标风险系数；最后，根据目标风险系数以及预设风险转换关系确定监测区域的目标风险等级。上述风险等级确定方法可以准确地确定矿山的风险等级，提高了对矿山风险的监测准确率。
附图说明
为了更清楚地说明本申请实施例中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本申请一实施例提供的矿山风险等级的确定方法的实现流程图；
图2是本申请一实施例提供的矿山风险等级的确定方法中S101的具体实现流程图；
图3是本申请一实施例提供的矿山风险等级的确定方法中S102的具体实现流程图；
图4是本申请一实施例提供的矿山风险等级的确定装置的结构示意图；
图5是本申请一实施例提供的服务器的结构示意图。
具体实施方式
以下描述中，为了说明而不是为了限定，提出了诸如特定系统结构、技术之类的具体细节，以便透彻理解本申请实施例。然而，本领域的技术人员应当清楚，在没有这些具体细节的其它实施例中也可以实现本申请。在其它情况中，省略对众所周知的系统、装置、电路以及方法的详细说明，以免不必要的细节妨碍本申请的描述。
应当理解，当在本申请说明书和所附权利要求书中使用时，术语“包括”指示所描述特征、整体、步骤、操作、元素和/或组件的存在，但并不排除一个或多个其它特征、整体、步骤、操作、元素、组件和/或其集合的存在或添加。
还应当理解，在本申请说明书和所附权利要求书中使用的术语“和/或”是指相关联列出的项中的一个或多个的任何组合以及所有可能组合，并且包括这些组合。
如在本申请说明书和所附权利要求书中所使用的那样，术语“如果”可以依据上下文被解释为“当...时”或“一旦”或“响应于确定”或“响应于检测到”。类似地，短语“如果确定”或“如果检测到[所描述条件或事件]”可以依据上下文被解释为意指“一旦确定”或“响应于确定”或“一旦检测到[所描述条件或事件]”或“响应于检测到[所描述条件或事件]”。
另外，在本申请说明书和所附权利要求书的描述中，术语“第一”、“第二”、“第三”等仅用于区分描述，而不能理解为指示或暗示相对重要性。
在本申请说明书中描述的参考“一个实施例”或“一些实施例”等意味着在本申请的一个或多个实施例中包括结合该实施例描述的特定特征、结构或特点。由此，在本说明书中的不同之处出现的语句“在一个实施例中”、“在一些实施例中”、“在其他一些实施例中”、“在另外一些实施例中”等不是必然都参考相同的实施例，而是意味着“一个或多个但不是所有的实施例”，除非是以其他方式另外特别强调。术语“包括”、“包含”、“具有”及它们的变形都意味着“包括但不限于”，除非是以其他方式另外特别强调。
请参阅图1，图1是本申请一实施例提供的矿山风险等级的确定方法的实现流程图。本申请实施例中，该矿山风险等级的确定方法的执行主体为服务器。如图1所示，本申请一实施例提供的矿山风险等级的确定方法可以包括S101～S104，详述如下：
在S101中，获取监测区域在预设时间段内的多次微震的震级和微震指标信息。
在本申请实施例中，微震指标信息包括但不限于能量指数、视体积、施密特数及能量释放量等微震参数。
预设时间段可以根据实际需要设置，此处不作限制，示例性的，为了提高对确定矿山风险等级的准确率，预设时间段可以设置为3个月。
在本申请实施例的一种实现方式中，服务器可以从监测区域内部署的微震监测系统获取到预设时间段内的微震指标信息和多次微震的震级。
在本申请实施例的另一种实现方式中，服务器还可以实时从监测区域内部署的微震监测系统获取到每次微震的微震指标信息和震级并进行存储。在服务器需要对监测区域的风险等级进行确定时，直接从其内部获取到预设时间段内的微震指标信息和多次微震的震级。
在本申请的一个实施例中，S101可以具体包括如图2所示的步骤S201～S203，详述如下：
在S201中，获取所述监测区域在所述预设时间段内的多种微震参数；所述多种微震参数包括视体积、能量指数、能量释放量和/或施密特数。
在本实施例的一种实现方式中，服务器可以从监测区域内部署的微震监测系统获取到预设时间段内的多种微震参数。
在本实施例的另一种实现方式中，服务器还可以实时从监测区域内部署的微震监测系统获取到每次微震的微震参数并进行存储。在服务器需要获取微震指标信息时，直接从其内部获取到预设时间段内的多种微震参数。
能量指数指一次微震所产生的实测辐射微震能与监测区域内所有微震的平均微震能之比。服务器可以根据以下公式计算得到平均微震能：
Figure BDA0003188791540000071
其中，
Figure BDA0003188791540000072
表示平均微震能，P表示微震体变势，d表示比例因子，c表示常数系数。
视体积用于描述震源非弹性变形区岩体的体积，可以根据震源体变势和应变改变量计算得到。具体的，服务器可以根据以下公式计算得到视体积：
V＝P/Δε；
其中，V表示视体积，P表示微震体变势，Δε表示应变改变量。
应变改变量可以根据微震体变势和角频率计算得到。具体的，服务器可以根据以下公式计算得到应变改变量：
Δε＝cPf0 3；
其中，Δε表示应变改变量，c表示常数系数，P表示微震体变势，f0表示角频率。
施密特数是一个无量纲的标量，定义为运动黏性系数和扩散系数的比值，用来描述同时有动量扩散及质量扩散的流体。
能量释放量指一次微震释放的能量。
在S202中，针对任一微震参数，生成所述微震参数随时间变化的变化曲线。
本实施例中，变化曲线用于描述微震参数在预设时间段内各个时刻对应的数值。可选的，服务器可以以预设时间段内的各个时刻为X轴，以给个时刻对应的微震参数的数值为Y轴，生成该微震参数的变化曲线。
在S203中，根据所述多种微震参数的变化曲线确定所述微震指标信息。
本实施例中，服务器可以根据生成的微震参数的变化曲线，计算该变化曲线的曲线斜率，并将多种微震参数的变化曲线对应的曲线斜率确定为微震指标信息。
在S102中，根据每次微震的震级和微震指标信息确定所述监测区域的微震特征信息。
本申请实施例中，微震特征信息包括但不限于第一微震特征信息和第二微震特征信息。其中，第一微震特征信息包括但不限于预设时间段内每天发生的微震的次数和预设时间段内大小震级的比值。第二微震特征信息包括但不限于微震指标信息和目标次数。其中，目标次数指发生的震级大于第一预设等级的微震事件的次数。优选的，第一预设等级可以设置为1.0。
基于此，S102可以具体包括如图3所示的步骤S301～S304，详述如下：
在S301中，获取所述监测区域在所述预设时间段内每天发生的微震的次数。
在本申请的一个实施例中，由于监测区域内部署的微震监测系统可以实时监测到微震事件的发生，因此，服务器可以通过微震监测系统在预设时间段内，每天统计发生微震的次数。
在本申请另一个实施例中，服务器可以预先从监测区域内部署的微震监测系统获取到预设时间段内的每次微震发生的时间并进行存储。当服务器需要确定预设时间段内每天发生的微震次数时，服务器可以根据每次微震发生的时间确定每天发生的微震次数。
在S302中，根据每次所述微震的震级确定所述预设时间段内大小震级的比值。
本实施例中，服务器在获取到预设时间段内每次微震的震级时，可以将每次微震的震级与第二预设等级进行比较。
服务器在检测到微震的震级大于或等于第二预设等级时，确定该微震的震级为大震级。服务器在检测到微震的震级小于第二预设等级时，确定该微震的震级为小震级。
服务器在确定大小震级后，可以确定预设时间段内大小震级的比值。其中，大小震级的比值特指小震级与大震级之间的比值。
在S303中，根据所述大小震级的比值和所述每天发生的微震的次数确定所述第一微震特征信息。
本实施例中，服务器可以根据预设时间段内，每天发生的微震的次数生成每天发生的微震的次数随时间变化的变化曲线，并通过上述变化曲线确定曲线斜率。
基于此，服务器可以根据预设时间段内的大小震级的比值和上述曲线斜率确定第一微震特征信息。
在本申请的一个实施例中，第一微震特征信息还可以包括视应力。其中，视应力用于描述震源非弹性变形区岩体的应力。
因此，服务器可以生成预设时间段内视应力随时间变化的变化曲线，并确定上述变化曲线的曲线斜率。
本实施例中，服务器可以根据预设时间段内的大小震级的比值、每天发生的微震的次数随时间变化的变化曲线的曲线斜率及视应力随时间变化的变化曲线的曲线斜率确定第一微震特征信息。
在S304中，根据目标次数和所述微震指标信息确定所述第二微震特征信息；所述目标次数指发生的震级大于第一预设等级的微震事件的次数。
本实施例中，结合S201，微震指标信息包括能量指数、视体积、能量释放量和/或施密特数，因此，服务器可以根据预设时间段内能量指数的变化曲线、视体积的变化曲线、能量释放量的变化曲线和/或施密特数的变化曲线，以及大于或等于第三预设等级的微震发生的具体时刻(具体到时)，服务器可以确定在大于或等于第三预设等级的微震发生的时刻之前，能量指数、视体积、能量释放量和/或施密特数的变化规律，并以此确定前兆特征信息。服务器还可以在大于或等于第三预设等级的微震发生的具体时刻时，获取该具体时刻时能量指数、视体积、能量释放量和/或施密特数的具体数值，并以此确定主震特征信息。其中，主震指大于或等于第三预设等级的微震。第三预设等级可以根据实际需要设置，此处不作限制。
基于此，服务器可以根据预设时间段内的目标次数、前兆特征信息及主震特征信息确定第二微震特征信息。
在S103中，将所述微震指标信息和所述微震特征信息导入预设的风险预测模型，计算得到所述监测区域的目标风险系数。
本申请实施例中，由于微震特征信息包括第一微震特征信息和第二微震特征信息，因此，服务器可以根据以下公式计算得到目标风险系数：
F＝w1*f1+w2*f2+w3*f3；
其中，F表示目标风险系数，f1表示微震指标信息对应的第一风险因子，f2表示第一微震特征信息对应的第二风险因子，f3表示第二微震特征信息对应的第三风险因子，w1表示第一风险因子对应的权重系数，w2表示第二风险因子对应的权重系数，w3表示第三风险因子对应的权重系数。
在本申请的一个实施例中，由于第一风险因子、第二风险因子及第三风险因子在变化时，对目标风险系数的影响是不同的，因此，服务器可以根据预设权重计算方法计算得到第一风险因子对应的权重系数、第二风险因子对应的权重系数及第三风险因子对应的权重系数。其中，预设权重计算方法可以根据实际需要确定，此处不作限制，优选的，预设权重系数计算方法可以是层次分析法。
在本申请的一个实施例中，由于微震指标信息可以包括能量指数、视体积、能量释放量和/或施密特数，第一微震特征信息可以包括大小震级的比值、每天发生微震的次数和/或视应力，第二微震特征信息可以包括目标次数和/或前兆特征信息、主震特征信息，因此，服务器可以根据以下公式计算目标风险系数：
Figure BDA0003188791540000101
其中，F表示目标风险系数，fi表示第i个参数对应的风险因子，wi表示第i个风险因子对应的权重系数。其中，参数可以包括能量指数、视体积、能量释放量、施密特数、大小震级的比值、每天发生微震的次数、视应力、包括目标次数、前兆特征信息及主震特征信息。
基于此，服务器可以根据层次分析法确定上述参数的风险因子各自对应的权重系数。优选的，能量指数、视体积、能量释放量、施密特数、大小震级的比值、每天发生微震的次数、视应力、包括目标次数、前兆特征信息及主震特征信息的风险因子对应的标准权重系数分别为：0.2276、0.1404、0.073、0.0524、0.0609、0.1533、0.0966、0.0979、0.049及0.049。
需要说明的是，在本申请的一个实施例中，当视体积和能量释放量的变化曲线的曲线斜率处于(0.5,0.75]时，视体积和能量释放量的风险因子对应的实际权重系数为其标准权重系数的一半，当视体积和能量释放量的变化曲线的曲线斜率处于(0.75,1]时，视体积和能量释放量的风险因子对应的实际权重系数等于其标准权重系数。其中，视体积的风险因子即为视体积的变化曲线的曲线斜率，能量释放量的风险因子即为能量释放量的变化曲线的曲线斜率。
当能量指数、施密特数及视应力的变化曲线的曲线斜率处于(-0.5,-0.75]时，能量指数、施密特数及视应力的风险因子对应的实际权重系数为其标准权重系数的一半，当能量指数、施密特数及视应力的变化曲线的曲线斜率处于(-0.75,-1]时，能量指数、施密特数及视应力风险因子对应的实际权重系数等于其标准权重系数。其中，能量指数的风险因子即为能量指数的变化曲线的曲线斜率，施密特数的风险因子即为施密特数的变化曲线的曲线斜率，视应力的风险因子即为视应力的变化曲线的曲线斜率。
当每天发生的微震的次数处于(5,10)时，每天发生的微震的次数的风险因子对应的实际权重系数为其标准权重系数的一半，当每天发生的微震的次数大于10时，每天发生的微震的次数的风险因子对应的实际权重系数等于其标准权重系数。其中，每天发生的微震的次数的风险因子即为每天发生的微震的次数。
在S104中，根据所述目标风险系数以及预设风险转换关系确定所述监测区域的目标风险等级。
本申请实施例中，预设风险转换关系可以根据实际需要确定，此处不作限制。
在本申请的一个实施例中，服务器可以根据以下公式计算目标风险等级：
nr＝A*nf；
其中，nr表示所述目标风险等级，A表示转换因子，表示所述目标风险系数nf。
在本申请实施例的一种实现方式中，风险系数可以设置为百分制，风险等级包括预设数目个等级，因此，服务器可以根据将满分风险系数与预设数目之间的比值确定为转换因子。其中，预设数目可以根据实际需要设置，此处不作限制，优选的，预设数目为4。风险等级可以用于描述矿山岩体状态和地压活动强度。
基于此，当目标风险等级为第一等级时，说明矿山岩体稳定，地压活动明显；当目标风险等级为第二等级时，说明矿山岩体较稳定，地压活动较明显；当目标风险等级为第三等级时，说明矿山岩体不稳定，地压活动明显。当目标风险等级为第四等级时，说明矿山岩体极不稳定，地压活动剧烈。
在本申请的一个实施例中，服务器在确定目标风险等级后，还可以根据目标风险等级以及预设风险等级与预设策略之间的对应关系，确定目标风险等级对应的目标策略，并执行目标策略。
本实施例中，不同的预设风险等级可以对应不同的预设策略。其中，预设策略包括输出提示信息和输出报警信息。
提示信息用于提示相关人员加强监测区域内的安全排查，特别注意监测区域内的浮石冒落和小范围片帮风险。基于此，服务器可以通过网络连接的显示屏显示上述提示信息，也可以通过麦克风播报上述提示信息。
报警信息用于描述矿山风险大，需要工作人员撤出监测区域，并搬走监测区域内的重要设备。基于此，服务器可以通过警报器输出报警信息。报警信息可以是警铃声。
由于当预设风险等级大于第三等级时，矿山岩体极不稳定，地压活动剧烈，容易造成人员伤害和设备损坏。因此，当预设风险等级为第四等级时，对应的预设策略可以为输出报警信息。当预设风险等级小于或等于第三等级时，对应的预设策略可以为输出提示信息。
基于此，服务器在检测到目标风险等级大于第三等级时，可以确定其对应的目标策略为输出报警信息。
服务器在检测到目标风险等级小于或等于第三等级时，可以确定其对应的目标策略为输出提示信息。
以上可以看出，本申请一实施例提供的矿山风险等级的确定方法，通过获取到的监测区域在预设时间段内的所有微震的震级和微震指标信息，可以确定监测区域的第一微震特征信息和第二微震特征信息，并将上述微震指标信息、第一微震特征信息及第二微震特征信息导入预设的风险预测模型，可以计算得到监测区域的目标风险系数；最后，根据目标风险系数以及预设风险转换关系确定监测区域的目标风险等级。上述风险等级确定方法可以准确地确定矿山的风险等级，提高了对矿山风险的监测准确率。
应理解，上述实施例中各步骤的序号的大小并不意味着执行顺序的先后，各过程的执行顺序应以其功能和内在逻辑确定，而不应对本申请实施例的实施过程构成任何限定。
对应于上文实施例所述的一种矿山风险等级的确定方法，图4示出了本申请一实施例提供的矿山风险等级的确定装置的结构框图，为了便于说明，仅示出了与本申请实施例相关的部分。参照图4，该矿山风险等级的确定装置400包括：第一获取单元41、第一确定单元42、计算单元43及第二确定单元44。其中：
第一获取单元41用于获取监测区域在预设时间段内的多次微震的震级和微震指标信息。
第一确定单元42用于根据每次微震的震级和微震指标信息确定所述监测区域的微震特征信息。
计算单元43用于将所述微震指标信息和所述微震特征信息导入预设的风险预测模型，计算得到所述监测区域的目标风险系数。
第二确定单元44用于根据所述目标风险系数以及预设风险转换关系确定所述监测区域的目标风险等级。
在本申请的一个实施例中，第一获取单元41具体包括：第三获取单元、生成单元及第三确定单元。其中：
第三获取单元用于获取所述监测区域在所述预设时间段内的多种微震参数；所述多种微震参数包括视体积、能量指数、能量释放量和/或施密特数。
生成单元用于针对任一微震参数，生成所述微震参数随时间变化的变化曲线。
第三确定单元用于根据所述多种微震参数的变化曲线确定所述微震指标信息。
在本申请的一个实施例中，所述微震特征信息包括第一微震特征信息和第二微震特征信息，第一确定单元42具体包括：第四获取单元、第四确定单元、第五确定单元及第六确定单元。其中：
第四获取单元用于获取所述监测区域在所述预设时间段内每天发生的微震的次数。
第四确定单元用于根据每次所述微震的震级确定所述预设时间段内大小震级的比值。
第五确定单元用于根据所述大小震级的比值和所述每天发生的微震的次数确定所述第一微震特征信息。
第六确定单元用于根据目标次数和所述微震指标信息确定所述第二微震特征信息；所述目标次数指发生的震级大于第一预设等级的微震事件的次数。
在本申请的一个实施例中，所述所述目标风险系数根据以下公式计算得到：
F＝w1*f1+w2*f2+w3*f3；
其中，F表示所述目标风险系数，f1表示所述微震指标信息对应的第一风险因子，f2表示所述第一微震特征信息对应的第二风险因子，f3表示所述第二微震特征信息对应的第三风险因子，w1表示所述第一风险因子对应的权重系数，w2表示所述第二风险因子对应的权重系数，w3表示所述第三风险因子对应的权重系数。
在本申请的一个实施例中，矿山风险等级的确定装置400还包括：第七确定单元。
第七确定单元用于根据预设权重计算方法确定所述第一风险因子、第二风险因子及第三风险因子分别对应的权重系数。
在本申请的一个实施例中，所述目标风险等级根据以下公式计算得到：
nr＝A*nf；
其中，nr表示所述目标风险等级，A表示转换因子，表示所述目标风险系数nf。
在本申请的一个实施例中，矿山风险等级的确定装置400还包括：执行单元。
执行单元用于根据所述目标风险等级以及预设风险等级与预设策略之间的对应关系，确定所述目标风险等级对应的目标策略，并执行所述目标策略。
以上可以看出，本申请一实施例提供的矿山风险等级的确定装置，通过获取到的监测区域在预设时间段内的所有微震的震级和微震指标信息，可以确定监测区域的第一微震特征信息和第二微震特征信息，并将上述微震指标信息、第一微震特征信息及第二微震特征信息导入预设的风险预测模型，可以计算得到监测区域的目标风险系数；最后，根据目标风险系数以及预设风险转换关系确定监测区域的目标风险等级。上述风险等级确定方法可以准确地确定矿山的风险等级，提高了对矿山风险的监测准确率。
需要说明的是，上述装置/单元之间的信息交互、执行过程等内容，由于与本申请方法实施例基于同一构思，其具体功能及带来的技术效果，具体可参见方法实施例部分，此处不再赘述。
所属领域的技术人员可以清楚地了解到，为了描述的方便和简洁，仅以上述各功能单元、模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能单元、模块完成，即将所述装置的内部结构划分成不同的功能单元或模块，以完成以上描述的全部或者部分功能。实施例中的各功能单元、模块可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中，上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。另外，各功能单元、模块的具体名称也只是为了便于相互区分，并不用于限制本申请的保护范围。上述系统中单元、模块的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
图5为本申请一实施例提供的服务器的结构示意图。如图5所示，该实施例的服务器5包括：至少一个处理器50(图5中仅示出一个)处理器、存储器51以及存储在所述存储器51中并可在所述至少一个处理器50上运行的计算机程序52，所述处理器50执行所述计算机程序52时实现上述任意各个矿山风险等级的确定方法实施例中的步骤。
该服务器可包括，但不仅限于，处理器50、存储器51。本领域技术人员可以理解，图5仅仅是服务器5的举例，并不构成对服务器5的限定，可以包括比图示更多或更少的部件，或者组合某些部件，或者不同的部件，例如还可以包括输入输出设备、网络接入设备等。
所称处理器50可以是中央处理单元(Central Processing Unit，CPU)，该处理器50还可以是其他通用处理器、数字信号处理器(Digital Signal Processor，DSP)、专用集成电路(Application Specific Integrated Circuit，ASIC)、现成可编程门阵列(Field-Programmable Gate Array，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。
所述存储器51在一些实施例中可以是所述服务器5的内部存储单元，例如服务器5的内存。所述存储器51在另一些实施例中也可以是所述服务器5的外部存储设备，例如所述服务器5上配备的插接式硬盘，智能存储卡(Smart Media Card,SMC)，安全数字(SecureDigital,SD)卡，闪存卡(Flash Card)等。进一步地，所述存储器51还可以既包括所述服务器5的内部存储单元也包括外部存储设备。所述存储器51用于存储操作系统、应用程序、引导装载程序(BootLoader)、数据以及其他程序等，例如所述计算机程序的程序代码等。所述存储器51还可以用于暂时地存储已经输出或者将要输出的数据。
本申请实施例还提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现可实现上述各个方法实施例中的步骤。
本申请实施例提供了一种计算机程序产品，当计算机程序产品在服务器上运行时，使得服务器执行时实现可实现上述各个方法实施例中的步骤。
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本申请实现上述实施例方法中的全部或部分流程，可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于计算机可读存储介质中，该计算机程序在被处理器执行时，可实现上述各个方法实施例的步骤。其中，所述计算机程序包括计算机程序代码，所述计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。所述计算机可读介质至少可以包括：能够将计算机程序代码携带到终端设备的任何实体或装置、记录介质、计算机存储器、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、电载波信号、电信信号以及软件分发介质。例如U盘、移动硬盘、磁碟或者光盘等。在某些司法管辖区，根据立法和专利实践，计算机可读介质不可以是电载波信号和电信信号。
在上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详述或记载的部分，可以参见其它实施例的相关描述。
本领域普通技术人员可以意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本申请的范围。
在本申请所提供的实施例中，应该理解到，所揭露的矿山风险等级的确定装置和方法，可以通过其它的方式实现。例如，以上所描述的装置/服务器实施例仅仅是示意性的，例如，所述模块或单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通讯连接可以是通过一些接口，装置或单元的间接耦合或通讯连接，可以是电性，机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
以上所述实施例仅用以说明本申请的技术方案，而非对其限制；尽管参照前述实施例对本申请进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本申请各实施例技术方案的精神和范围，均应包含在本申请的保护范围之内。

Claims (10)

1.一种矿山风险等级的确定方法，其特征在于，包括：
获取监测区域在预设时间段内的多次微震的震级和微震指标信息；
根据每次微震的震级和微震指标信息确定所述监测区域的微震特征信息；
将所述微震指标信息和所述微震特征信息导入预设的风险预测模型，计算得到所述监测区域的目标风险系数；
根据所述目标风险系数以及预设风险转换关系确定所述监测区域的目标风险等级。
2.如权利要求1所述的矿山风险等级的确定方法，其特征在于，所述获取监测区域在预设时间段内的多次微震的震级和微震指标信息，包括：
获取所述监测区域在所述预设时间段内的多种微震参数；所述多种微震参数包括视体积、能量指数、能量释放量和/或施密特数；
针对任一微震参数，生成所述微震参数随时间变化的变化曲线；
根据所述多种微震参数的变化曲线确定所述微震指标信息。
3.如权利要求1或2所述的矿山风险等级的确定方法，其特征在于，所述微震特征信息包括第一微震特征信息和第二微震特征信息，所述根据每次微震的震级和微震指标信息确定所述监测区域的微震特征信息，包括：
获取所述监测区域在所述预设时间段内每天发生的微震的次数；
根据每次所述微震的震级确定所述预设时间段内大小震级的比值；
根据所述大小震级的比值和所述每天发生的微震的次数确定所述第一微震特征信息；
根据目标次数和所述微震指标信息确定所述第二微震特征信息；所述目标次数指发生的震级大于第一预设等级的微震事件的次数。
4.如权利要求3所述的矿山风险等级的确定方法，其特征在于，所述目标风险系数根据以下公式计算得到：
F＝W1*f1+W2*f2+W3*f3；
其中，F表示所述目标风险系数，f1表示所述微震指标信息对应的第一风险因子，f2表示所述第一微震特征信息对应的第二风险因子，f3表示所述第二微震特征信息对应的第三风险因子，w1表示所述第一风险因子对应的权重系数，w2表示所述第二风险因子对应的权重系数，w3表示所述第三风险因子对应的权重系数。
5.如权利要求4所述的矿山风险等级的确定方法，其特征在于，所述将所述微震指标信息和所述微震特征信息导入预设的风险预测模型，计算得到所述监测区域的目标风险系数之前，还包括：
根据预设权重计算方法确定所述第一风险因子、第二风险因子及第三风险因子分别对应的权重系数。
6.如权利要求1或2或4或5任一项所述的矿山风险等级的确定方法，其特征在于，所述目标风险等级根据以下公式计算得到：
nr＝A*nf；
其中，nr表示所述目标风险等级，A表示转换因子，表示所述目标风险系数nf。
7.如权利要求6所述的矿山风险等级的确定方法，其特征在于，所述根据所述目标风险系数确定所述监测区域的目标风险等级之后，还包括：
根据所述目标风险等级以及预设风险等级与预设策略之间的对应关系，确定所述目标风险等级对应的目标策略，并执行所述目标策略。
8.一种矿山风险等级的确定装置，其特征在于，包括：
第一获取单元，用于获取监测区域在预设时间段内的多次微震的震级和微震指标信息；
第一确定单元，用于根据每次微震的震级和微震指标信息确定所述监测区域的微震特征信息；
计算单元，用于将所述微震指标信息和所述微震特征信息导入预设的风险预测模型，计算得到所述监测区域的目标风险系数；
第二确定单元，用于根据所述目标风险系数以及预设风险转换关系确定所述监测区域的目标风险等级。
9.一种服务器，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现如权利要求1至7任一项所述的矿山风险等级的确定方法的步骤。
10.一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1至7任一项所述的矿山风险等级的确定方法的步骤。
CN202110870685.5A 2021-07-30 2021-07-30 一种矿山风险等级的确定方法、装置及服务器 Pending CN113642865A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110870685.5A CN113642865A (zh) 2021-07-30 2021-07-30 一种矿山风险等级的确定方法、装置及服务器

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110870685.5A CN113642865A (zh) 2021-07-30 2021-07-30 一种矿山风险等级的确定方法、装置及服务器

Publications (1)

Publication Number Publication Date
CN113642865A true CN113642865A (zh) 2021-11-12

Family

ID=78419051

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110870685.5A Pending CN113642865A (zh) 2021-07-30 2021-07-30 一种矿山风险等级的确定方法、装置及服务器

Country Status (1)

Country Link
CN (1) CN113642865A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114330991A (zh) * 2021-11-17 2022-04-12 华能核能技术研究院有限公司 核电生产工作风险等级的确定方法及装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103291364A (zh) * 2013-06-21 2013-09-11 中国矿业大学 一种冲击矿压的微震多维信息综合时序预警方法
CN105257339A (zh) * 2015-11-16 2016-01-20 北京安科兴业科技股份有限公司 掘进工作面多参量综合监测预警方法
CN112379425A (zh) * 2020-10-21 2021-02-19 中国矿业大学 基于震源参数反演的区域冲击危险等级及畛域的判识方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103291364A (zh) * 2013-06-21 2013-09-11 中国矿业大学 一种冲击矿压的微震多维信息综合时序预警方法
CN105257339A (zh) * 2015-11-16 2016-01-20 北京安科兴业科技股份有限公司 掘进工作面多参量综合监测预警方法
CN112379425A (zh) * 2020-10-21 2021-02-19 中国矿业大学 基于震源参数反演的区域冲击危险等级及畛域的判识方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张楚旋等: "基于微震活动性参数的矿山灾害风险性评价", 《安全与环境学报》, vol. 17, no. 4, pages 1276 - 1280 *
陈祥祥等: "某钨矿地压微震监测系统技术应用研究", 《中国钨业》, vol. 31, no. 4, pages 26 - 30 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114330991A (zh) * 2021-11-17 2022-04-12 华能核能技术研究院有限公司 核电生产工作风险等级的确定方法及装置

Similar Documents

Publication Publication Date Title
Jordan et al. Operational earthquake forecasting can enhance earthquake preparedness
US20210173107A1 (en) Microseismic Wireless Monitoring, Receiving and Early Warning System of Rock
Kong et al. Numerical analysis of roadway rock-burst hazard under superposed dynamic and static loads
CN111290924A (zh) 监控方法、装置及电子设备
Taborda et al. Earthquake ground‐motion simulation including nonlinear soil effects under idealized conditions with application to two case studies
CN113642865A (zh) 一种矿山风险等级的确定方法、装置及服务器
CN107370766A (zh) 一种网络流量异常检测方法及系统
CN106289371A (zh) 水下爆破安全测定装置与方法
Dang et al. A novel method based on maximum likelihood estimation for the construction of seismic fragility curves using numerical simulations
CN114415116A (zh) 一种煤矿开采监测方法、装置及计算机设备
CN109578075A (zh) 冲击地压危险的微震监测预警方法、装置和系统
Xu et al. Monitoring of train-induced vibrations on rock slopes
Chen et al. Rockburst precursors and the dynamic failure mechanism of the deep tunnel: A review
CN115620851B (zh) 数据校验方法、装置、电子设备及可读存储介质
Li et al. Research on Engineering Practice and Effect Evaluation Method of Pressure Relief in Deep Rock Burst Danger Area of Coal Mine
CN115688428A (zh) 爆破近区支护结构爆炸药量临界值计算方法、装置和设备
Ma et al. Numerical investigation of safety strategy for gas disaster prevention in successive panels using upper protective layer mining: a case study
CN112699559A (zh) 一种基于沉降速率信息熵安全值的安全评价方法及系统
Zhou et al. Study on Dominant Frequency Attenuation of Blasting Vibration for Ultra-Small-Spacing Tunnel
Yan et al. Research on the method and model for calculating impact load in the rockburst tunnel
Ding et al. Surrounding rock pressure calculation based on time functions and stress release rate determination of deep soft rock tunnel: taking Zhonghe Tunnel as an example
Wang et al. Optimization of BP neural network model for rockburst prediction under multiple influence factors
Yang et al. Predicting the accuracy and applicability of micro-seismic monitoring of rock burst in TBM tunneling using the data from two case studies in China
Li et al. Stability analysis and fracture mechanism of surrounding rock in a hydraulic tunnel based on microseismic monitoring
Rodríguez et al. Representative-Area Approach to Define Blast-Induced Ground Vibrations—Damage Prevention Criterion Abacus

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