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 PDFInfo
- 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
Links
- 238000005065 mining Methods 0.000 title claims abstract description 95
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000004364 calculation method Methods 0.000 claims abstract description 29
- 238000004458 analytical method Methods 0.000 claims abstract description 12
- 238000009826 distribution Methods 0.000 claims description 11
- 230000005641 tunneling Effects 0.000 claims description 9
- 238000010586 diagram Methods 0.000 claims description 7
- 230000001186 cumulative effect Effects 0.000 claims description 3
- 238000005315 distribution function Methods 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 abstract description 17
- 238000013178 mathematical model Methods 0.000 abstract description 2
- 238000011160 research Methods 0.000 abstract description 2
- 238000005070 sampling Methods 0.000 abstract 1
- 239000003245 coal Substances 0.000 description 28
- 238000009412 basement excavation Methods 0.000 description 14
- 238000011161 development Methods 0.000 description 9
- 238000002360 preparation method Methods 0.000 description 8
- 238000000605 extraction Methods 0.000 description 7
- 238000011084 recovery Methods 0.000 description 5
- 230000015572 biosynthetic process Effects 0.000 description 4
- 101100012902 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) FIG2 gene Proteins 0.000 description 2
- 238000005086 pumping Methods 0.000 description 2
- 101001121408 Homo sapiens L-amino-acid oxidase Proteins 0.000 description 1
- 102100026388 L-amino-acid oxidase Human genes 0.000 description 1
- 101100233916 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) KAR5 gene Proteins 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 238000012502 risk assessment Methods 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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
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
Description
技术领域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:
(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:
式中,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)=∑i∑jP(IEi)P(ASGj|IEi)P(A|ASGj);P(A)=∑ i ∑ j 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
令连续型随机变量X服从参数为(α,λ)的伽马分布,记作X~Ga(α,λ),其中α>0且λ>0,则变量X的数学期望和方差分别为和故变量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 and Therefore, the probability density function of variable X is:
通过MATLAB计算得到伽马函数的累计分布函数为:The cumulative distribution function of the gamma function calculated by MATLAB is:
其中,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:
式中,Tk为开拓煤量合理可采期,Lk为掘进工程总长度,p为开拓过程的平均掘进队个数,Vk为开拓过程的平均掘进速度,Tk0为形成开拓煤量后到形成准备煤量前的平均准备接替时间,Tz2为石门揭煤平均预抽时间,Tz3为井筒揭煤平均预抽时间,Vd为掘进过程中受复杂地质构造影响到所降低的最大平均掘进速度,为减缓速度的影响系数,且为区间[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, is the influence coefficient of slowing down speed, and is a constant in the interval [0, 1];
2)准备煤量合理可采期模型如下:2) Prepare a reasonable mining period model for coal quantity as follows:
式中,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为瓦斯预抽安全系数,且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 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:
(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:
式中,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)=∑i∑jP(IEi)P(ASGj|IEi)P(A|ASGj);P(A)=∑ i ∑ j 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的数学期望和方差分别为和故变量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 and Therefore, the probability density function of variable X is:
通过MATLAB计算得到伽马函数的累计分布函数为:The cumulative distribution function of the gamma function calculated by MATLAB is:
其中,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,则问题即为计算满足方程下的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 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)=∑i∑jP(IEi)P(ASGj|IEi)P(A|ASGj)=33.6%;P(A)=∑ i ∑ j 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:
通过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)
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)
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 |
-
2018
- 2018-10-22 CN CN201811227594.4A patent/CN109359393B/en not_active Expired - Fee Related
Patent Citations (8)
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)
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 |