CN109359393B - Method for Determination of Reasonable Mining Period of Four Quantities in Mine Under Uncertain Information Conditions - Google Patents

Method for Determination of Reasonable Mining Period of Four Quantities in Mine Under Uncertain Information Conditions Download PDF

Info

Publication number
CN109359393B
CN109359393B CN201811227594.4A CN201811227594A CN109359393B CN 109359393 B CN109359393 B CN 109359393B CN 201811227594 A CN201811227594 A CN 201811227594A CN 109359393 B CN109359393 B CN 109359393B
Authority
CN
China
Prior art keywords
reasonable
quantities
probability
mining
model
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.)
Expired - Fee Related
Application number
CN201811227594.4A
Other languages
Chinese (zh)
Other versions
CN109359393A (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.)
Liaoning Technical University
Shanxi Coking Coal Group Co Ltd
Original Assignee
Liaoning Technical University
Shanxi Coking Coal 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 Liaoning Technical University, Shanxi Coking Coal Group Co Ltd filed Critical Liaoning Technical University
Priority to CN201811227594.4A priority Critical patent/CN109359393B/en
Publication of CN109359393A publication Critical patent/CN109359393A/en
Application granted granted Critical
Publication of CN109359393B publication Critical patent/CN109359393B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0635Risk analysis of enterprise or organisation activities

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Economics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • General Physics & Mathematics (AREA)
  • Strategic Management (AREA)
  • Educational Administration (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Tourism & Hospitality (AREA)
  • Marketing (AREA)
  • General Business, Economics & Management (AREA)
  • Game Theory and Decision Science (AREA)
  • Development Economics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention provides a method for determining four reasonable mining periods of a mine under uncertain information conditions, and relates to the technical field of mine safety mining. The method comprises the following steps: determining a four-quantity reasonable mining period model; measuring the influence of uncertain information conditions on the reasonable sampling period of four quantities; analyzing uncertain information conditions which cause model failure, and quantitatively calculating the risk of the model; and (5) performing reliability analysis on the model. The method for determining the four reasonable mining periods of the mine under the uncertain information condition provided by the invention researches the change rule of the four reasonable mining period models under various uncertain information conditions, perfects the traditional four reasonable mining period calculation models and theoretical basis, establishes a corresponding mathematical model by analyzing the risk and reliability of the models under the uncertain information condition, provides guidance opinion for the establishment of a mine safety production plan, and provides theoretical basis for normal production connection of a working face.

Description

不确定信息条件下的矿井四量合理可采期确定方法Method for determining the reasonable mining period of four quantities in a mine under uncertain information conditions

技术领域Technical Field

本发明涉及矿山安全开采技术领域,具体涉及一种不确定信息条件下的矿井四量合理可采期确定方法。The invention relates to the technical field of safe mining of mines, and in particular to a method for determining a reasonable mining period of four quantities in a mine under uncertain information conditions.

背景技术Background Art

在煤炭生产过程中,保持采掘平衡是确保矿井正常、连续、稳定生产的前提条件,也是降低掘进工程维护成本的重要手段。随着我国煤炭行业的迅猛发展、煤炭生产科技技术进步以及开采条件复杂程度等的不断变化,国家制定的为保证矿井顺利接替的《“三量”规定》已逐渐出现采掘失调、不能适应现阶段矿井接替的需要。“四量”即开拓煤量、准备煤量期、回采煤量,及“四量”合理可采期确定方法的提出在很大程度上改进了传统“三量”理论在瓦斯预抽方面的不足,使得在计算高瓦斯和瓦斯突出矿井的合理可采期时有了可参考的理论基础。In the process of coal production, maintaining the balance between mining and excavation is a prerequisite for ensuring the normal, continuous and stable production of the mine, and is also an important means to reduce the maintenance cost of the tunneling project. With the rapid development of my country's coal industry, the progress of coal production science and technology, and the continuous changes in the complexity of mining conditions, the "Three Quantities" regulations formulated by the state to ensure the smooth succession of mines have gradually become unbalanced in mining and excavation and cannot meet the needs of mine succession at this stage. The "four quantities" are the development coal quantity, the preparation coal quantity period, and the recovery coal quantity, and the method for determining the reasonable mining period of the "four quantities" has greatly improved the shortcomings of the traditional "three quantities" theory in gas pre-extraction, providing a theoretical basis for reference when calculating the reasonable mining period of high-gas and gas outburst mines.

现阶段,矿井的生产实践活动越来越受到多种不确定信息条件的影响,如:不可预见的复杂地质突变、突发灾害、外部经济与市场需求、合理配采、开拓开采方法等,这些无论是生产过程中可能遇到的内在因素还是周围环境所带来的外部因素,都是不能够预先确定的,但它们又对矿井“四量”合理可采期都产生着不可忽视的影响。然而,目前人们在矿井不确定信息条件下“四量”合理可采期变化规律的研究甚少,在“四量”可采期的风险性与可靠性方面也尚无可借鉴的成果。At present, the production practice of mines is increasingly affected by a variety of uncertain information conditions, such as: unforeseen complex geological mutations, sudden disasters, external economic and market demand, reasonable allocation and mining, and pioneering mining methods. These factors, whether internal factors that may be encountered in the production process or external factors brought by the surrounding environment, cannot be determined in advance, but they have an impact on the reasonable mining period of the "four quantities" of the mine that cannot be ignored. However, there is currently little research on the changing laws of the reasonable mining period of the "four quantities" under the conditions of uncertain information in mines, and there is no reference to the risk and reliability of the mining period of the "four quantities".

发明内容Summary of the invention

针对现有技术存在的问题,本发明提供一种不确定信息条件下的矿井四量合理可采期确定方法,通过探究不确定信息条件下“四量”合理可采期的变化规律,对“四量”合理可采期模型进行风险性分析以及可靠性分析,从而得到“四量”合理可采期模型的风险函数,为工作面的正常生产接续提供指导作用。In view of the problems existing in the prior art, the present invention provides a method for determining the reasonable mining period of four quantities in a mine under uncertain information conditions. By exploring the changing law of the reasonable mining period of the "four quantities" under uncertain information conditions, risk analysis and reliability analysis are performed on the reasonable mining period model of the "four quantities", so as to obtain the risk function of the reasonable mining period model of the "four quantities", which provides guidance for the normal production continuation of the working face.

为了实现上述目的,一种不确定信息条件下的矿井四量合理可采期确定方法,包括以下步骤:In order to achieve the above purpose, a method for determining the reasonable mining period of four quantities in a mine under uncertain information conditions includes the following steps:

步骤1:确定“四量”合理可采期模型;Step 1: Determine the reasonable mining period model of the "four quantities";

步骤2:确定不确定信息条件下“四量”合理可采期变化规律,并以掘进速度作为随机变量来度量不确定信息条件对“四量”合理可采期的影响;Step 2: Determine the changing rules of the reasonable recoverable period of the "four quantities" under uncertain information conditions, and use the tunneling speed as a random variable to measure the impact of uncertain information conditions on the reasonable recoverable period of the "four quantities";

步骤3:通过主逻辑图分析导致“四量”合理可采期模型失效的不确定信息条件,并以故障树与事件树的概率计算为基础对模型的风险性进行定量计算;Step 3: Analyze the uncertain information conditions that lead to the failure of the "four quantities" reasonable recoverable period model through the main logic diagram, and quantitatively calculate the risk of the model based on the probability calculation of the fault tree and event tree;

步骤4:对“四量”合理可采期模型进行可靠性分析。Step 4: Conduct reliability analysis on the reasonable mining period model of the “four quantities”.

进一步地,所述步骤3中以故障树与事件树的概率计算为基础对模型的风险性进行定量计算的方法如下:Furthermore, the method for quantitatively calculating the risk of the model based on the probability calculation of the fault tree and the event tree in step 3 is as follows:

(1)故障树的概率计算:(1) Probability calculation of fault tree:

令导致采掘失衡的n个因素为H1,H2,H3,...,Hn,且因素间有交集,则导致采掘失衡的总概率的计算公式如下:Let n factors that lead to mining imbalance be H 1 , H 2 , H 3 , ..., H n , and there is an intersection between the factors, then the calculation formula for the total probability of leading to mining imbalance is as follows:

Figure BDA0001836360820000021
Figure BDA0001836360820000021

(2)事件树的概率计算:(2) Probability calculation of event tree:

令“四量”合理可采期模型中事件树的m个事故序列组为ASG1,ASG2,...,ASGm,且第j个事故序列组中包含的l个环节事件为Tj1,Tj2,...,Tjl,则第j个事故序列组条件概率的计算公式如下:Let the m accident sequence groups in the event tree of the “four-quantity” reasonable mining period model be ASG 1 , ASG 2 , ... , ASG m , and the l link events contained in the jth accident sequence group be T j1 , T j2 , ... , T jl , then the calculation formula for the conditional probability of the jth accident sequence group is as follows:

Figure BDA0001836360820000022
Figure BDA0001836360820000022

式中,P(ASGj|IE)为不确定事件IE发生时第j个连锁事故事件发生的概率,p(Tkl)为事故序列组中第l个环节事件发生的概率;Where, P(ASG j |IE) is the probability of the jth chain accident event occurring when the uncertain event IE occurs, and p(T kl ) is the probability of the lth link event occurring in the accident sequence group;

(3)“四量”合理可采期模型的风险概率计算:(3) Risk probability calculation of the “four quantities” reasonable mining period model:

根据导致采掘失衡的总概率和每个事故序列组的条件概率得到“四量”合理可采期模型的风险概率计算公式如下:According to the total probability of causing mining imbalance and the conditional probability of each accident sequence group, the risk probability calculation formula of the "four quantities" reasonable mining period model is as follows:

P(A)=∑ijP(IEi)P(ASGj|IEi)P(A|ASGj);P(A)=∑ ij P(IE i )P(ASG j |IE i )P(A|ASG j );

式中,P(A)为“四量”合理可采期模型失效概率,P(IEi)为发生第i个可能引发采掘失衡的不确定事件概率,P(ASGj|IEi)为不确定事件i发生时第j个连锁事故事件发生的概率,P(A|ASGj)为第j个连锁事故事件是否导致模型失效的概率。Where P(A) is the failure probability of the reasonable mining period model of the “four quantities”, P(IE i ) is the probability of the occurrence of the i-th uncertain event that may cause mining imbalance, P(ASG j |IE i ) is the probability of the occurrence of the j-th chain accident event when the uncertain event i occurs, and P(A|ASG j ) is the probability of whether the j-th chain accident event will lead to model failure.

进一步地,所述步骤4中对“四量”合理可采期模型进行可靠性分析的方法如下:Furthermore, the method for performing reliability analysis on the reasonable mining period model of the "four quantities" in step 4 is as follows:

令连续型随机变量X服从参数为(α,λ)的伽马分布,记作X~Ga(α,λ),其中α>0且λ>0,则变量X的数学期望和方差分别为

Figure BDA0001836360820000023
Figure BDA0001836360820000024
故变量X的概率密度函数为:Let the continuous random variable X obey the gamma distribution with parameters (α, λ), denoted as X~Ga(α, λ), where α>0 and λ>0, then the mathematical expectation and variance of the variable X are
Figure BDA0001836360820000023
and
Figure BDA0001836360820000024
Therefore, the probability density function of variable X is:

Figure BDA0001836360820000025
Figure BDA0001836360820000025

通过MATLAB计算得到伽马函数的累计分布函数为:The cumulative distribution function of the gamma function calculated by MATLAB is:

Figure BDA0001836360820000026
Figure BDA0001836360820000026

其中,X为“四量”合理可采期中任一合理可采期,x为合理可采期函数下的变量值,igamma为倒伽马分布;Among them, X is any reasonable mining period among the “four quantities”, x is the variable value under the reasonable mining period function, and igamma is the inverse gamma distribution;

则在给定可靠度的前提下,计算满足方程Fx(x)为给定可靠度下的x值,即为x值对应下的“四量”合理可采期模型可靠度。Under the premise of a given reliability, the value of x that satisfies the equation F x (x) is calculated to be the value of x under the given reliability, that is, the reliability of the reasonable mining period model of the "four quantities" corresponding to the value of x.

本发明的有益效果:Beneficial effects of the present invention:

本发明提出一种考虑不确定信息条件影响的矿井四量合理可采期确定方法,研究了“四量”合理可采期模型在多种不确定信息条件下的变化规律,完善了传统“四量”合理可采期计算模型和理论基础,并通过对模型在不确定信息条件下的风险性和可靠性进行分析,建立了相应的数学模型,为矿井安全生产计划的制定提供了指导意见,为工作面的正常生产接续提供了理论依据。The present invention proposes a method for determining the reasonable mining period of four quantities in a mine taking into account the influence of uncertain information conditions, studies the changing law of the reasonable mining period model of the "four quantities" under various uncertain information conditions, improves the traditional "four quantities" reasonable mining period calculation model and theoretical basis, and establishes a corresponding mathematical model by analyzing the risk and reliability of the model under uncertain information conditions, which provides guidance for the formulation of mine safety production plans and provides a theoretical basis for the normal production continuation of the working face.

附图说明BRIEF DESCRIPTION OF THE DRAWINGS

图1为本发明实施例中考虑不确定信息条件影响的矿井四量合理可采期确定方法流程图;FIG1 is a flow chart of a method for determining a reasonable mining period of four quantities in a mine taking into account the influence of uncertain information conditions in an embodiment of the present invention;

图2为本发明实施例中理论上的伽马分布示意图;FIG2 is a schematic diagram of a theoretical gamma distribution in an embodiment of the present invention;

图3为本发明实施例中实际结合矿井“四量”合理可采期应用的伽马分布示意图。3 is a schematic diagram of the gamma distribution actually applied in combination with the reasonable mining period of the “four quantities” of a mine in an embodiment of the present invention.

具体实施方式DETAILED DESCRIPTION

为了使本发明的目的、技术方案及优势更加清晰,下面结合附图和具体实施例对本发明做进一步详细说明。此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。In order to make the purpose, technical solution and advantages of the present invention clearer, the present invention is further described in detail below in conjunction with the accompanying drawings and specific embodiments. The specific embodiments described here are only used to explain the present invention and are not used to limit the present invention.

一种考虑不确定信息条件影响的矿井四量合理可采期确定方法,流程如图1所示,具体方法如下所述:A method for determining the reasonable mining period of four quantities in a mine considering the influence of uncertain information conditions, the process is shown in Figure 1, and the specific method is as follows:

步骤1:确定“四量”合理可采期模型;Step 1: Determine the reasonable mining period model of the "four quantities";

本实施例中,“四量”合理可采期模型如下:In this embodiment, the reasonable mining period model of the “four quantities” is as follows:

1)开拓煤量合理可采期模型如下:1) The model of reasonable mining period of coal development is as follows:

Figure BDA0001836360820000031
Figure BDA0001836360820000031

式中,Tk为开拓煤量合理可采期,Lk为掘进工程总长度,p为开拓过程的平均掘进队个数,Vk为开拓过程的平均掘进速度,Tk0为形成开拓煤量后到形成准备煤量前的平均准备接替时间,Tz2为石门揭煤平均预抽时间,Tz3为井筒揭煤平均预抽时间,Vd为掘进过程中受复杂地质构造影响到所降低的最大平均掘进速度,

Figure BDA0001836360820000032
为减缓速度的影响系数,且
Figure BDA0001836360820000033
为区间[0,1]内的常数;Where, Tk is the reasonable mining period of the development coal volume, Lk is the total length of the tunneling project, p is the average number of tunneling teams in the development process, Vk is the average tunneling speed in the development process, Tk0 is the average preparation replacement time from the formation of the development coal volume to the formation of the preparation coal volume, Tz2 is the average pre-pumping time for coal uncovering in the stone gate, Tz3 is the average pre-pumping time for coal uncovering in the shaft, Vd is the maximum average tunneling speed reduced by the influence of complex geological structures during the tunneling process,
Figure BDA0001836360820000032
is the influence coefficient of slowing down speed, and
Figure BDA0001836360820000033
is a constant in the interval [0, 1];

2)准备煤量合理可采期模型如下:2) Prepare a reasonable mining period model for coal quantity as follows:

Figure BDA0001836360820000034
Figure BDA0001836360820000034

式中,Tz为准备煤量合理可采期,Lz为准备工作所需的掘进总长度,q为准备过程的平均掘进队个数,Vz为准备过程的平均掘进速度,Tz0为形成准备煤量后到形成回采煤量前的平均准备接替时间,Tz1为准备巷道揭露煤层前的平均预抽时间;Where, Tz is the reasonable mining period of the prepared coal volume, Lz is the total excavation length required for the preparation work, q is the average number of excavation teams in the preparation process, Vz is the average excavation speed in the preparation process, Tz0 is the average preparation replacement time from the formation of the prepared coal volume to the formation of the recovered coal volume, and Tz1 is the average pre-extraction time before the preparation roadway exposes the coal seam;

3)回采煤量合理可采期模型如下:3) The model of reasonable recoverable period of coal recovery volume is as follows:

Th=Lm/(VJ-φVd)+max{To,Tc};T h =L m /(V J -φV d )+max{T o , T c };

式中,Th为工作面回采煤量合理可采期,To为工作面平均设备安装调试时间,Tc为工作面平均瓦斯抽采达标时间,Lm为回采工作所需的掘进总长度,VJ为工作面平均掘进速度;Where, Th is the reasonable period of coal recovery at the working face, To is the average equipment installation and commissioning time at the working face, Tc is the average gas extraction standard time at the working face, Lm is the total excavation length required for the recovery work, and VJ is the average excavation speed at the working face;

4)瓦斯抽采达标煤量合理可采期模型如下:4) The model of reasonable recoverable period of coal volume that meets gas extraction standards is as follows:

Tw=kTh Tw = kTh ;

式中,Tw为瓦斯抽采达标煤量可采期,k为瓦斯预抽安全系数,且

Figure BDA0001836360820000041
S1为回采煤量面积,S2为瓦斯预抽煤量面积。Where Tw is the period of coal mining that meets the gas extraction standard, k is the gas pre-extraction safety factor, and
Figure BDA0001836360820000041
S1 is the area of recovered coal volume, and S2 is the area of gas pre-extraction coal volume.

步骤2:确定不确定信息条件下“四量”合理可采期变化规律,并以掘进速度作为随机变量来度量不确定信息条件对“四量”合理可采期的影响;Step 2: Determine the changing rules of the reasonable recoverable period of the "four quantities" under uncertain information conditions, and use the tunneling speed as a random variable to measure the impact of uncertain information conditions on the reasonable recoverable period of the "four quantities";

本实施例中,以不可预见的复杂地质构造为例,通过分析,这些构造可能导致原计划可采的煤量掘巷后发现构造太多无法回采,之前掘进的巷道也可能报废,造成无效巷道的增多,等价于巷道采掘速度的降低,而回采工作又取决于掘进工作能否顺利完成,因此可认为不可预见的地质构造是通过造成掘进速度降低引发的“四量”合理可采期增长。同理,对于其他的不确定信息条件,如:突发灾害、外部经济与市场需求、合理配采、开拓开采方法等都可直接或间接转化为掘进速度的降低引起“四量”合理可采期的增长。In this embodiment, taking unforeseen complex geological structures as an example, through analysis, these structures may lead to the discovery that there are too many structures that cannot be mined after the originally planned mineable coal is excavated, and the previously excavated tunnels may also be scrapped, resulting in an increase in invalid tunnels, which is equivalent to a reduction in the tunnel mining speed. The mining work depends on whether the excavation work can be completed smoothly. Therefore, it can be considered that unforeseen geological structures cause the increase of the reasonable mineable period of the "four quantities" by causing a reduction in the excavation speed. Similarly, for other uncertain information conditions, such as sudden disasters, external economic and market demand, reasonable allocation and mining, and development mining methods, they can be directly or indirectly converted into a reduction in the excavation speed, causing an increase in the reasonable mineable period of the "four quantities".

步骤3:通过主逻辑图分析导致“四量”合理可采期模型失效的不确定信息条件,并以故障树与事件树的概率计算为基础对模型的风险性进行定量计算;Step 3: Analyze the uncertain information conditions that lead to the failure of the "four quantities" reasonable recoverable period model through the main logic diagram, and quantitatively calculate the risk of the model based on the probability calculation of the fault tree and event tree;

所述以故障树与事件树的概率计算为基础对模型的风险性进行定量计算的方法如下:The method for quantitatively calculating the risk of the model based on the probability calculation of the fault tree and the event tree is as follows:

(1)故障树的概率计算:(1) Probability calculation of fault tree:

令导致采掘失衡的n个因素为H1,H2,H3,...,Hn,且因素间有交集,则导致采掘失衡的总概率的计算公式如下:Let n factors that lead to mining imbalance be H 1 , H 2 , H 3 , ..., H n , and there is an intersection between the factors, then the calculation formula for the total probability of leading to mining imbalance is as follows:

Figure BDA0001836360820000042
Figure BDA0001836360820000042

(2)事件树的概率计算:(2) Probability calculation of event tree:

令“四量”合理可采期模型中事件树的m个事故序列组为ASG1,ASG2,...,ASGm,且第j个事故序列组中包含的l个环节事件为Tj1,Tj2,...,Tjl,则第j个事故序列组条件概率的计算公式如下:Let the m accident sequence groups in the event tree of the “four-quantity” reasonable mining period model be ASG 1 , ASG 2 , ... , ASG m , and the l link events contained in the jth accident sequence group be T j1 , T j2 , ... , T jl , then the calculation formula for the conditional probability of the jth accident sequence group is as follows:

Figure BDA0001836360820000051
Figure BDA0001836360820000051

式中,P(ASGj|IE)为不确定事件IE发生时第j个连锁事故事件发生的概率,p(Tkl)为事故序列组中第l个环节事件发生的概率;Where, P(ASG j |IE) is the probability of the jth chain accident event occurring when the uncertain event IE occurs, and p(T kl ) is the probability of the lth link event occurring in the accident sequence group;

(3)“四量”合理可采期模型的风险概率计算:(3) Risk probability calculation of the “four quantities” reasonable mining period model:

根据导致采掘失衡的总概率和每个事故序列组的条件概率得到“四量”合理可采期模型的风险概率计算公式如下:According to the total probability of causing mining imbalance and the conditional probability of each accident sequence group, the risk probability calculation formula of the "four quantities" reasonable mining period model is as follows:

P(A)=∑ijP(IEi)P(ASGj|IEi)P(A|ASGj);P(A)=∑ ij P(IE i )P(ASG j |IE i )P(A|ASG j );

式中,P(A)为“四量”合理可采期模型失效概率,P(IEi)为发生第i个可能引发采掘失衡的不确定事件概率,P(ASGj|IEi)为不确定事件i发生时第j个连锁事故事件发生的概率,P(A|ASGj)为第j个连锁事故事件是否导致模型失效的概率。Where P(A) is the failure probability of the reasonable mining period model of the “four quantities”, P(IE i ) is the probability of the occurrence of the i-th uncertain event that may cause mining imbalance, P(ASG j |IE i ) is the probability of the occurrence of the j-th chain accident event when the uncertain event i occurs, and P(A|ASG j ) is the probability of whether the j-th chain accident event will lead to model failure.

本实施例中,在“四量”合理可采期模型中,主逻辑图(MLD)对风险事故发生的必要条件进行分级描述,从一个不期望发生的工作面采掘失衡出发,自上而下对模型采取分而治之的策略进行展开,逐步将事故分解成由初始事件所描述的子块,直到所有可能的初始事件集合为止。事件树分析(ETA)则是通过事故序列审查与事故发展分析,得到引起采掘失衡因素所导致的事故序列组ASG。事件树中的事件节点表示由上一层事件所可能导致的后果事件。对事件树的每个不确定因素要求建立发生的联合概率,以确定该点的发生概率。故障树(FTA)是一个演绎推理的过程,这种推理通常为事件树的失效环节事件提供更为详实的细节。通过FTA可以清楚地了解模型是由于何种不确定因素引起失效的,从而找出导致采掘失衡的基本原因。对于“四量”合理可采期在实际应用中,通过统计相关的数据,利用定量计算的方法可以确切的得到“四量”合理可采期模型失效的概率,或生产过程中不能正常接续的概率,从而能够根据实际情况预先对生产计划作出相应的调整,把风险带来的损失降到最低。In this embodiment, in the "four-quantity" reasonable mining period model, the main logic diagram (MLD) describes the necessary conditions for the occurrence of risk accidents in a hierarchical manner. Starting from an undesirable mining imbalance at a working face, the model is developed from top to bottom with a divide-and-conquer strategy, and the accident is gradually decomposed into sub-blocks described by the initial event until all possible initial events are collected. Event tree analysis (ETA) obtains the accident sequence group ASG caused by the factors causing mining imbalance through accident sequence review and accident development analysis. The event node in the event tree represents the consequence event that may be caused by the event of the previous layer. For each uncertain factor in the event tree, it is required to establish a joint probability of occurrence to determine the probability of occurrence at that point. Fault tree (FTA) is a process of deductive reasoning, which usually provides more detailed details for the failure link events of the event tree. Through FTA, it is possible to clearly understand which uncertain factors cause the model to fail, so as to find out the basic cause of mining imbalance. In practical applications of the reasonable mining period of the "four quantities", by statistically analyzing relevant data and using quantitative calculation methods, we can accurately obtain the probability of failure of the reasonable mining period model of the "four quantities" or the probability of failure to continue normally during the production process. This will enable us to make corresponding adjustments to the production plan in advance based on actual conditions and minimize the losses caused by risks.

步骤4:对“四量”合理可采期模型进行可靠性分析。Step 4: Conduct reliability analysis on the reasonable mining period model of the “four quantities”.

所述对“四量”合理可采期模型进行可靠性分析的方法如下:The method for performing reliability analysis on the reasonable mining period model of the "four quantities" is as follows:

令连续型随机变量X服从参数为(α,λ)的伽马分布,记作X~Ga(α,λ),其中α>0且λ>0,则变量X的数学期望和方差分别为

Figure BDA0001836360820000052
Figure BDA0001836360820000053
故变量X的概率密度函数为:Let the continuous random variable X obey the gamma distribution with parameters (α, λ), denoted as X~Ga(α, λ), where α>0 and λ>0, then the mathematical expectation and variance of the variable X are
Figure BDA0001836360820000052
and
Figure BDA0001836360820000053
Therefore, the probability density function of variable X is:

Figure BDA0001836360820000054
Figure BDA0001836360820000054

通过MATLAB计算得到伽马函数的累计分布函数为:The cumulative distribution function of the gamma function calculated by MATLAB is:

Figure BDA0001836360820000061
Figure BDA0001836360820000061

其中,X为“四量”合理可采期中任一合理可采期,x为合理可采期函数下的变量值,igamma为倒伽马分布;Among them, X is any reasonable mining period among the “four quantities”, x is the variable value under the reasonable mining period function, and igamma is the inverse gamma distribution;

则在给定可靠度的前提下,计算满足方程F(x)为给定可靠度下的x值,即为x值对应下的“四量”合理可采期模型可靠度。Under the premise of a given reliability, the x value that satisfies the equation F(x) is calculated to be the value of x under the given reliability, that is, the reliability of the "four quantities" reasonable mining period model corresponding to the x value.

理论上的伽马函数分布如图2所示,本实施例中,以回采煤量合理可采期Th进行可靠性分析,在计算掘进时间的问题上,给定可靠度(概率)为0.9,则问题即为计算满足方程

Figure BDA0001836360820000063
下的x值,即在90%概率下可实现采掘顺利对接的回采煤量可采期Th=x,且在伽马函数分布图中x表现为在k点处的取值。The theoretical distribution of the gamma function is shown in FIG2. In this embodiment, the reliability analysis is performed based on the reasonable recoverable period of the recovered coal volume T h. In the problem of calculating the excavation time, the reliability (probability) is given as 0.9, and the problem is to calculate the equation satisfying
Figure BDA0001836360820000063
The x value under this condition, that is, the recoverable period of coal volume that can achieve smooth mining and excavation connection under 90% probability is Th = x, and in the gamma function distribution diagram, x is represented as the value at point k.

本实施例中,以某瓦斯突出矿井为例,通过工作人员已统计数据结合实际经验得知,在采掘过程事先没有预知且遇到较大断层的情况下,受断层影响必须重新开掘工作面的概率为80%;受该断层影响可能发生周围环境其他突发灾害的概率为60%;因后续可能产生的突发灾害影响导致模型失效的概率为70%,最终求得“四量”合理可采期模型失效概率如下:In this embodiment, a gas outburst mine is taken as an example. Through the statistical data collected by the staff and the actual experience, it is known that when the mining process is not predicted in advance and a large fault is encountered, the probability that the working face must be re-excavated due to the influence of the fault is 80%; the probability that other sudden disasters in the surrounding environment may occur due to the influence of the fault is 60%; the probability of model failure due to the influence of subsequent sudden disasters is 70%. Finally, the failure probability of the "four quantities" reasonable mining period model is obtained as follows:

P(A)=∑ijP(IEi)P(ASGj|IEi)P(A|ASGj)=33.6%;P(A)=∑ ij P(IE i )P(ASG j |IE i )P(A|ASG j )=33.6%;

通过计算结果得知,由于受到不可预知的断层影响,使得“四量”合理可采期失效的概率达到33.6%,显然给生产过程的正常接续带来了较大的风险。The calculation results show that due to the influence of unpredictable faults, the probability of failure of the reasonable mining period of the "four quantities" reaches 33.6%, which obviously brings great risks to the normal continuation of the production process.

本实施例中,实际应用时的伽马函数分布如图3所示,为了降低不确定信息条件给生产接续带来的风险,从而得到最佳的回采煤量合理可采期Th,计算模型在达到90%可靠度下的掘进时间,具体如下:In this embodiment, the gamma function distribution in actual application is shown in FIG3 . In order to reduce the risk of uncertain information conditions to production continuity, and thus obtain the optimal reasonable recoverable period Th of the recovered coal volume, the excavation time of the calculation model under a reliability of 90% is calculated as follows:

经统计该矿井某采区相邻几条巷道掘进时间后,计算求得期望E(X)和方差Var(X)分别为8.95和2.01,参数α和λ分别为40.051和4.475,由于本实施例中x为时间,需取正值,因此将参数α和λ代入如下公式求解x:After counting the excavation time of several adjacent tunnels in a mining area of the mine, the expected E(X) and variance Var(X) are calculated to be 8.95 and 2.01 respectively, and the parameters α and λ are 40.051 and 4.475 respectively. Since x is time in this embodiment and needs to be a positive value, the parameters α and λ are substituted into the following formula to solve x:

Figure BDA0001836360820000062
Figure BDA0001836360820000062

通过MATLAB求得时间x为7.192,即在图3中k点的取值。The time x is obtained by MATLAB as 7.192, which is the value of point k in Figure 3.

由图3可以看出,在不确定信息条件存在的前提下,回采煤量合理可采期为7.192个月能够使矿井生产接续正常达到90%以上的概率,可以为矿井的生产计划提供指导意见,从而保证生产的正常接续。As can be seen from Figure 3, under the premise of uncertain information conditions, the reasonable recoverable period of coal recovery is 7.192 months, which can make the probability of normal mine production continuation reach more than 90%, and can provide guidance for the mine's production plan, thereby ensuring the normal continuation of production.

最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解;其依然可以对前述实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;因而这些修改或者替换,并不使相应技术方案的本质脱离本发明权利要求所限定的范围。Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that the technical solutions described in the aforementioned embodiments may still be modified, or some or all of the technical features therein may be replaced by equivalents. Therefore, these modifications or replacements do not cause the essence of the corresponding technical solutions to deviate from the scope defined by the claims of the present invention.

Claims (2)

1.一种不确定信息条件下的矿井四量合理可采期确定方法,其特征在于,包括以下步骤:1. A method for determining the reasonable mining period of four quantities in a mine under uncertain information conditions, characterized in that it comprises the following steps: 步骤1:确定“四量”合理可采期模型;Step 1: Determine the reasonable mining period model of the "four quantities"; 步骤2:确定不确定信息条件下“四量”合理可采期变化规律,并以掘进速度作为随机变量来度量不确定信息条件对“四量”合理可采期的影响;Step 2: Determine the changing rules of the reasonable recoverable period of the "four quantities" under uncertain information conditions, and use the tunneling speed as a random variable to measure the impact of uncertain information conditions on the reasonable recoverable period of the "four quantities"; 步骤3:通过主逻辑图分析导致“四量”合理可采期模型失效的不确定信息条件,并以故障树与事件树的概率计算为基础对模型的风险性进行定量计算;Step 3: Analyze the uncertain information conditions that lead to the failure of the "four quantities" reasonable recoverable period model through the main logic diagram, and quantitatively calculate the risk of the model based on the probability calculation of the fault tree and event tree; 步骤4:对“四量”合理可采期模型进行可靠性分析;Step 4: Conduct reliability analysis on the reasonable mining period model of the “four quantities”; 所述步骤3中以故障树与事件树的概率计算为基础对模型的风险性进行定量计算的方法如下:The method for quantitatively calculating the risk of the model based on the probability calculation of the fault tree and the event tree in step 3 is as follows: (1)故障树的概率计算:(1) Probability calculation of fault tree: 令导致采掘失衡的n个因素为H1,H2,H3,...,Hn,且因素间有交集,则导致采掘失衡的总概率的计算公式如下:Let n factors that lead to mining imbalance be H 1 , H 2 , H 3 , ..., H n , and there is an intersection between the factors, then the calculation formula for the total probability of leading to mining imbalance is as follows:
Figure FDA0004080168890000011
Figure FDA0004080168890000011
(2)事件树的概率计算:(2) Probability calculation of event tree: 令“四量”合理可采期模型中事件树的m个事故序列组为ASG1,ASG2,...,ASGm,且第j个事故序列组中包含的l个环节事件为Tj1,Tj2,...,Tjl,则第j个事故序列组条件概率的计算公式如下:Let the m accident sequence groups in the event tree of the “four-quantity” reasonable mining period model be ASG 1 , ASG 2 , ... , ASG m , and the l link events contained in the jth accident sequence group be T j1 , T j2 , ... , T jl , then the calculation formula for the conditional probability of the jth accident sequence group is as follows:
Figure FDA0004080168890000012
Figure FDA0004080168890000012
式中,P(ASGj|IE)为不确定事件IE发生时第j个连锁事故事件发生的概率,p(Tkl)为事故序列组中第l个环节事件发生的概率;Where, P(ASGj|IE) is the probability of the jth chain accident event occurring when the uncertain event IE occurs, and p(T kl ) is the probability of the lth link event occurring in the accident sequence group; (3)“四量”合理可采期模型的风险概率计算:(3) Risk probability calculation of the “four quantities” reasonable mining period model: 根据导致采掘失衡的总概率和每个事故序列组的条件概率得到“四量”合理可采期模型的风险概率计算公式如下:According to the total probability of causing mining imbalance and the conditional probability of each accident sequence group, the risk probability calculation formula of the "four quantities" reasonable mining period model is as follows: P(A)=∑ijP(IEi)P(ASGj|IEi)P(A|ASGj);P(A)=∑ ij P(IE i )P(ASG j |IE i )P(A|ASG j ); 式中,P(A)为“四量”合理可采期模型失效概率,P(IEi)为发生第i个可能引发采掘失衡的不确定事件概率,P(ASGj|IEi)为不确定事件i发生时第j个连锁事故事件发生的概率,P(A|ASGj)为第j个连锁事故事件是否导致模型失效的概率。Where P(A) is the failure probability of the reasonable mining period model of the “four quantities”, P(IE i ) is the probability of the occurrence of the i-th uncertain event that may cause mining imbalance, P(ASG j |IE i ) is the probability of the occurrence of the j-th chain accident event when the uncertain event i occurs, and P(A|ASG j ) is the probability of whether the j-th chain accident event will lead to model failure.
2.根据权利要求1所述的不确定信息条件下的矿井四量合理可采期确定方法,其特征在于,所述步骤4中对“四量”合理可采期模型进行可靠性分析的方法如下:2. The method for determining the reasonable mining period of four quantities in a mine under uncertain information conditions according to claim 1 is characterized in that the method for performing reliability analysis on the reasonable mining period model of the "four quantities" in step 4 is as follows: 令连续型随机变量X服从参数为(α,λ)的伽马分布,记作X~Ga(α,λ),其中α>0且λ>0,则变量X的数学期望和方差分别为
Figure FDA0004080168890000013
Figure FDA0004080168890000014
故变量X的概率密度函数为:
Let the continuous random variable X obey the gamma distribution with parameters (α, λ), denoted as X~Ga(α, λ), where α>0 and λ>0, then the mathematical expectation and variance of the variable X are
Figure FDA0004080168890000013
and
Figure FDA0004080168890000014
Therefore, the probability density function of variable X is:
Figure FDA0004080168890000015
Figure FDA0004080168890000015
通过MATLAB计算得到伽马函数的累计分布函数为:The cumulative distribution function of the gamma function calculated by MATLAB is:
Figure FDA0004080168890000021
Figure FDA0004080168890000021
其中,X为“四量”合理可采期中任一合理可采期,x为合理可采期函数下的变量值,igamma为倒伽马分布;Among them, X is any reasonable mining period among the “four quantities”, x is the variable value under the reasonable mining period function, and igamma is the inverse gamma distribution; 则在给定可靠度的前提下,计算满足方程F(x)为给定可靠度下的x值,即为x值对应下的“四量”合理可采期模型可靠度。Under the premise of a given reliability, the x value that satisfies the equation F(x) is calculated to be the value of x under the given reliability, that is, the reliability of the "four quantities" reasonable mining period model corresponding to the x value.
CN201811227594.4A 2018-10-22 2018-10-22 Method for Determination of Reasonable Mining Period of Four Quantities in Mine Under Uncertain Information Conditions Expired - Fee Related CN109359393B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811227594.4A CN109359393B (en) 2018-10-22 2018-10-22 Method for Determination of Reasonable Mining Period of Four Quantities in Mine Under Uncertain Information Conditions

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811227594.4A CN109359393B (en) 2018-10-22 2018-10-22 Method for Determination of Reasonable Mining Period of Four Quantities in Mine Under Uncertain Information Conditions

Publications (2)

Publication Number Publication Date
CN109359393A CN109359393A (en) 2019-02-19
CN109359393B true CN109359393B (en) 2023-04-25

Family

ID=65346095

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811227594.4A Expired - Fee Related CN109359393B (en) 2018-10-22 2018-10-22 Method for Determination of Reasonable Mining Period of Four Quantities in Mine Under Uncertain Information Conditions

Country Status (1)

Country Link
CN (1) CN109359393B (en)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101392648A (en) * 2008-09-28 2009-03-25 中国地质大学(武汉) A device for advanced detection and sampling of coalbed methane in exploration boreholes
CN102495935A (en) * 2011-12-15 2012-06-13 中国石油天然气集团公司 Determination method for risk of storage medium leakage of underground natural gas storage reservoir
CN103016057A (en) * 2012-12-20 2013-04-03 大同煤矿集团有限责任公司 Prevention method for dynamic pressure area of working face
CN103912310A (en) * 2014-04-15 2014-07-09 辽宁工程技术大学 Prediction method for mine dynamic disasters
CN104890523A (en) * 2014-03-04 2015-09-09 通用电气公司 System and method for controlling energy usage
CN104950904A (en) * 2015-06-11 2015-09-30 陕西西科美芯科技集团有限公司 Self-flying method for mine emergent rescue aircraft
CN105626151A (en) * 2016-02-28 2016-06-01 辽宁工程技术大学 Coalmine stoping roadway impact ground pressure pre-warning method
CN106199305A (en) * 2016-07-01 2016-12-07 太原理工大学 Underground coal mine electric power system dry-type transformer insulation health state evaluation method

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101392648A (en) * 2008-09-28 2009-03-25 中国地质大学(武汉) A device for advanced detection and sampling of coalbed methane in exploration boreholes
CN102495935A (en) * 2011-12-15 2012-06-13 中国石油天然气集团公司 Determination method for risk of storage medium leakage of underground natural gas storage reservoir
CN103016057A (en) * 2012-12-20 2013-04-03 大同煤矿集团有限责任公司 Prevention method for dynamic pressure area of working face
CN104890523A (en) * 2014-03-04 2015-09-09 通用电气公司 System and method for controlling energy usage
CN103912310A (en) * 2014-04-15 2014-07-09 辽宁工程技术大学 Prediction method for mine dynamic disasters
CN104950904A (en) * 2015-06-11 2015-09-30 陕西西科美芯科技集团有限公司 Self-flying method for mine emergent rescue aircraft
CN105626151A (en) * 2016-02-28 2016-06-01 辽宁工程技术大学 Coalmine stoping roadway impact ground pressure pre-warning method
CN106199305A (en) * 2016-07-01 2016-12-07 太原理工大学 Underground coal mine electric power system dry-type transformer insulation health state evaluation method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
三个煤量合理可采期的确定;吕子辉等;《煤炭科学技术》;20010625;第45-47页 *
基于主逻辑图的安全风险建模研究;曾士勇等;《中国安全科学学报》;20031130;第66-68页 *

Also Published As

Publication number Publication date
CN109359393A (en) 2019-02-19

Similar Documents

Publication Publication Date Title
CN114580752B (en) Engineering construction progress intelligent management system based on BIM technology
CN103513293B (en) Tunnel geology synthetical forecast expert system and its implementation
CN110059963A (en) A kind of tunnel risk evaluating method based on fuzzy polymorphism Bayesian network
Yu et al. Probabilistic risk analysis of diversion tunnel construction simulation
CN110619483A (en) Tunnel surrounding rock grade dynamic change and decision-making method based on multi-source data fusion analysis
CN117634308B (en) Decision system and analysis method for analyzing rail transit engineering data
CN112100727A (en) Early warning prevention and control method for water inrush of water-rich tunnel under influence of fault fracture zone
CN111324988B (en) Gas overrun early warning model construction method and early warning method based on machine learning
CN101634229A (en) Risk-based tunnel supporting structure design method
WO2019091134A1 (en) Method for predicting water-rich levels of sandstone in coal seam roof
CN116777085B (en) Coal mine water damage prediction system based on data analysis and machine learning technology
Bo et al. Prediction of tunnel deformation using PSO variant integrated with XGBoost and its TBM jamming application
CN107091085B (en) Multi-parameter discrimination method for stratum stability of shallow-buried and underground-excavated tunnel
Wu et al. A dynamic decision approach for risk analysis in complex projects
CN115620513A (en) Urban road network state monitoring method and visualization system based on toughness analysis
CN116756953A (en) Dynamic optimization method and device for anchor support design during tunnel construction period
CN103122772B (en) Method for judging burst accident in early burst happening period rapidly and predicting gas discharge scale
CN109359393B (en) Method for Determination of Reasonable Mining Period of Four Quantities in Mine Under Uncertain Information Conditions
CN115759351A (en) Slurry shield tunneling comprehensive early warning method and system and storage medium
CN105550412B (en) A kind of determination method and system of reality drilling well and planned well coincidence rate
CN108573328A (en) A Method for Stability Evaluation of Mined-out Area in Mine
CN117910336A (en) Prediction method and system for ground subsidence probability of shield construction in initial stage
CN108256191B (en) Tunnel excavation progress simulation method considering construction abnormity
CN116295230A (en) An Intelligent Monitoring System for Deformation of Weak Surrounding Rock in Shallow Buried Large Section Tunnel
CN109299886A (en) Subway tunnel structure safe condition numerical analysis method based on integrated data base

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20230425