CN109446659A - A kind of method that integrated forecasting heading emits heights of roofs and hidden danger classification - Google Patents
A kind of method that integrated forecasting heading emits heights of roofs and hidden danger classification Download PDFInfo
- Publication number
- CN109446659A CN109446659A CN201811285772.9A CN201811285772A CN109446659A CN 109446659 A CN109446659 A CN 109446659A CN 201811285772 A CN201811285772 A CN 201811285772A CN 109446659 A CN109446659 A CN 109446659A
- Authority
- CN
- China
- Prior art keywords
- roof
- fall
- hidden danger
- roof fall
- rock
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Abstract
The present invention provides a kind of methods that integrated forecasting heading emits heights of roofs and hidden danger classification, it is characterized in that, include: step 1, is four kinds by fall of ground type i points: crack country rock body back part block roof fall, the break wall rock body back roof fall of butterfly plastic zone, weak bonding stratiform country rock body tunnel Combined Roof Cave and loose wall rock body back roof fall；Step 2, for every kind of roof fall type founding mathematical models, calculate separately on mine roof above-mentioned four kinds of roof fall types of k-th of coordinate points emits heights of roofs hidden danger value Hi；Step 3 obtains the weighted average coefficients of every kind of roof fall probability of happening, emits heights of roofs hidden danger value H according to four kinds of roof fall types according to the empirical probability that every kind of roof fall occursiWith the weighted average coefficients of every kind of roof fall probability of happening, the integrated forecasting value H for emitting heights of roofs is obtainedk.To solve the technical problem that back strata structure complexity existing in the prior art emits heights of roofs and hidden danger classification is difficult to comprehensive determination.
Description
Technical field
The present invention relates to fall of ground technical fields, emit heights of roofs and hidden danger more particularly, to integrated forecasting heading
The method of classification.
Background technique
The energy is the basis of Chinese national economy and social development, and one of China's main energy sources are still coal at present,
Ratio shared by coal is respectively 76% and 69% in all non-renewable energy production and consumption structures, and therefore, coal is quite long
It will be the main energy sources in China in period.
China's coal-mine is mainly pit mining, and mine top plate disaster threat is serious.The research of mine top plate disaster is related to ground
The comprehensive crossover project of the subjects such as matter, mechanics, mining, mathematics, computer, mine top plate disaster are typically no bright when occurring
Aobvious sign, it has sudden, instantaneity and concealed feature, once generation is destructive huge, should find in time in time
Prevention.During top plate disaster is mainly mining work activities, roof rock mass in process of deformation and failure stress field and displacement field transfer,
The mechanical process of accumulation and release is its external rock mass environment, internal lithologic structure, geological structure and its physico-mechanical properties
Concentrated expression, have apparent temporal-spatial evolution characteristics, therefore, mine top plate disaster has close contact with rock mass environment.
Gradually popularization with roadway bolting support technology in China, roadway supported by bolt is using more and more extensive, accordingly
Design of bolt support and support theory application also gradually from preparatory working, entity tunnel, small section tunnel and simple condition
The tunnel of stope drift active workings, gob side entry, large section roadway (open-off cut) and various complex conditions is developed in tunnel etc., different
Roadway bolting support Parameters design also occurs therewith, and some then needs by many and diverse calculating and obtains more parameter
Design can be provided, some then must accurately grasp mine construction geological conditions, however underground complex geologic conditions, coal mine
Technical staff is more to carry out design of its support to tunnel using the method for engineering analogy, same in order to save time and cost
One tunnel, the different tunnels in same exploiting field or even entire mine coal seam tunnel all use same roadway support form and supporting
Parameter has that supporting intensity is superfluous, wastes support material for the good tunnel of surrounding rock stability, for adjoining rock stability
Property difference there are the hidden danger of roof fall in tunnel.Therefore, only pay much attention to surrounding rock mass stability classification or Immediate roof stability classification
Importance, could fundamentally solve the problems, such as one, fall of ground accident is preferably minimized.
The information disclosed in the background technology section is intended only to deepen understanding of the general background technology to the present invention,
And it is not construed as recognizing or implying in any form that the information constitutes existing skill known to those skilled in the art
Art.
Summary of the invention
The purpose of the present invention is to provide a kind of integrated forecasting headings to emit the method that heights of roofs and hidden danger are classified, with
Solve the technology that back strata structure complexity existing in the prior art emits heights of roofs and hidden danger classification is difficult to comprehensive determination
Problem.
To achieve the goals above, the invention adopts the following technical scheme:
A kind of method that integrated forecasting heading emits heights of roofs and hidden danger classification, comprising:
Fall of ground type i points are four kinds by step 1: crack country rock body back part block roof fall, butterfly are moulded
Property area's break wall rock body back roof fall, weak bonding stratiform country rock body tunnel Combined Roof Cave and loose wall rock body tunnel
Top plate roof fall；
Step 2 calculates separately the kth coordinate points on mine roof for every kind of roof fall type founding mathematical models
Above-mentioned four kinds of roof fall types emit heights of roofs hidden danger value Hi；
Step 3 obtains the weighted average system of every kind of roof fall probability of happening according to the empirical probability that every kind of roof fall occurs
Number, emits heights of roofs hidden danger value H according to four kinds of roof fall typesiWith the weighted average coefficients of every kind of roof fall probability of happening, obtain
The integrated forecasting value H of heights of roofs must be emittedk；
Wherein, i=1,2,3,4；K=1,2 ... ..., n.
Further, crack country rock body back part block roof fall height hidden danger value H1, obtain as follows
: it establishes crack country rock body back Key Blocks and falls mechanical model, it is critical to obtain back hexahedron sillar
Quantitative relation formula between Metacentre Height and rock mass environment, it may be assumed that
Wherein, block free face length l, block structure face inclination alpha, angle of frictionSurvey pressure coefficient lambda.
Further, the butterfly plastic zone break wall rock body back emits heights of roofs hidden danger value H2, as follows
It obtains:
The implicit function expression formula on circular tunnel plastic zone boundary under non-isopiestic state:
According to tunnel buried depth H, country rock bulk density γ, bidirection press ratio λ, country rock cohesive force C and internal friction angleUsing upper
Plastic zone radial depth r when formula obtains roadway surrounding rock difference θ is stated, and then according to formula:
Obtain butterfly blade profile roof fall hidden danger height H2；
Further, the weak bonding stratiform country rock body tunnel Combined Roof Cave height hidden danger value H3, as follows
It obtains:
First according to roof strata lift height hn, distribution load qx, delamination strength σtx, mining influence coefficient fc, buried depth
Influence coefficient fm, ground stress anomally influence coefficient fd, rock stratum integrality influence coefficient fw, the stabilization of rock stratum is determined from bottom to up
Property, until encountering stable rock stratum, distinguishing rule and criterion are as follows:
When the limited span [L] of the rock stratum is less than the practical span L in tunnel, then the rock stratum unstability is caving；The pole of the rock stratum
When limiting the practical span L of span [L] greater than tunnel, and two kinds of forms: form are subdivided into according to the length for forming girder construction
One, if the mistake girder construction formed tends to clamped beam, the rock stability；Form two, if the mistake girder construction formed is between clamped beam
When between simply supported beam, judge whether the rock stratum meets the sliding instability condition and hinged sillar structure of hinged sillar structure
Deformation instability condition, meeting any instability condition, then the rock stratum unstability is caving, otherwise the rock stability；
The unstability formula of stratiform rock beam limited span:
The sliding instability formula of hinged sillar structure:
The deformation instability formula of hinged sillar structure:
Thus it obtains,Obtain H3。
Further, the loose wall rock body back emits heights of roofs hidden danger value H4, it obtains as follows:
Wherein, tunnel radius α, head-room H, country rock internal friction anglePu Shi
Coefficient fTop。
Further, hkIt (i) is the i-th seed type of k-th of coordinate points of mine fall of ground height hidden danger value；
wkIt (i) is the corresponding weighted average coefficients of k-th of coordinate points of mine, i-th kind of fall of ground type；
HkFor the fall of ground high integrity predicted value of k-th of coordinate points of mine；
Wherein, i=1,2,3,4；K=1,2 ..., n.
Further, emitting heights of roofs and determine its roof fall rank according to each coordinate points, so that back be divided
For the region with different roof fall ranks.
By adopting the above technical scheme, the invention has the following beneficial effects:
The method that integrated forecasting heading provided by the invention emits heights of roofs and hidden danger classification, by the complete of fall of ground
Portion's situation specifically sums up are as follows: blocky roof fall, stratiform roof fall, butterfly-shaped roof fall and loose sandbox roof fall these four roof fall forms, it is first
First obtain every kind of roof fall form emits heights of roofs, heights of roofs is emitted using fuzzy mathematics method integrated forecasting heading, by tunnel
Roof fall hidden danger has according to the different demarcation rank for emitting heights of roofs convenient for drawing out roof roof fall hidden danger rank zoning map
Help propose targeted design of its support, for the accident hair for instructing mine roadway support design, reducing high tunnel of taking risks
Raw probability guarantees that mine safety production has important theory significance and practical value.
Detailed description of the invention
It, below will be to tool in order to illustrate more clearly of the specific embodiment of the invention or technical solution in the prior art
Body embodiment or attached drawing needed to be used in the description of the prior art be briefly described, it should be apparent that, it is described below
Attached drawing be some embodiments of the present invention, for those of ordinary skill in the art, what is do not made the creative labor
Under the premise of, it is also possible to obtain other drawings based on these drawings.
Fig. 1 is top plate hexahedron stability of block analysis model provided in an embodiment of the present invention；
Fig. 2 is roadway surrounding rock butterfly provided in an embodiment of the present invention plastic zone stability analysis model；
Fig. 3 is layered roof rock beam mechanical model in tunnel provided in an embodiment of the present invention；
Fig. 4 is loose wall rock body back mechanical model provided in an embodiment of the present invention；
Fig. 5 is country rock body back in crack provided in an embodiment of the present invention part block neutrality height H1It is related
Parameter typing and calculating interface；
Fig. 6 is weak bonding stratiform country rock body tunnel compound roof neutrality height H provided in an embodiment of the present invention2Phase
It closes parameter typing and calculates interface；
Fig. 7 is butterfly plastic zone break wall rock body back neutrality height H provided in an embodiment of the present invention3Phase
It closes parameter typing and calculates interface；
Fig. 8 is loose wall rock body back neutrality height H provided in an embodiment of the present invention4Relevant parameter typing
With calculating interface；
Fig. 9 is heading roof fall hidden danger classification results provided in an embodiment of the present invention interface；
Figure 10 is one borehole roof fall hidden danger rank distribution schematic diagram provided in an embodiment of the present invention；
Figure 11 is full mine roof roof fall hidden danger rank graduation schematic diagram provided in an embodiment of the present invention.
Specific embodiment
Technical solution of the present invention is clearly and completely described below in conjunction with attached drawing, it is clear that described reality
Applying example is a part of the embodiment of the present invention, instead of all the embodiments.Based on the embodiments of the present invention, this field is common
Technical staff's every other embodiment obtained without making creative work belongs to what the present invention protected
Range.
In the description of the present invention, it should be noted that term " center ", "upper", "lower", "left", "right", "vertical",
The orientation or positional relationship of the instructions such as "horizontal", "inner", "outside" be based on the orientation or positional relationship shown in the drawings, be only for
Convenient for the description present invention and simplify description, rather than the device or element of indication or suggestion meaning there must be specific side
Position is constructed and operated in a specific orientation, therefore is not considered as limiting the invention.In addition, term " first ", " the
Two ", " third " is used for descriptive purposes only and cannot be understood as indicating or suggesting relative importance.
In the description of the present invention, it should be noted that unless otherwise clearly defined and limited, term " installation ",
" connected ", " connection " shall be understood in a broad sense, for example, it may be being fixedly connected, may be a detachable connection, or integrally connect
It connects；It can be mechanical connection, be also possible to be electrically connected；Can be directly connected, can also indirectly connected through an intermediary,
It can be the connection inside two elements.For the ordinary skill in the art, can understand as the case may be
The concrete meaning of above-mentioned term in the present invention.
Below in conjunction with attached drawing, detailed description of the preferred embodiments.It should be understood that this place is retouched
The specific embodiment stated is merely to illustrate and explain the present invention, and is not intended to restrict the invention.
As shown in Figure 1, Figure 2, Fig. 3, Fig. 4, Fig. 5, Fig. 6, Fig. 7, Fig. 8, Fig. 9, Figure 10 and Figure 11 show integrated forecasting of the present invention
Heading emits the specific embodiment of the method for heights of roofs and hidden danger classification.Following steps have been specifically included in the present embodiment:
Whole fall of ground type i points are four kinds: crack country rock body back part block roof fall by step 1
(i.e. blocky roof fall), butterfly plastic zone break wall rock body back roof fall (butterfly-shaped roof fall), weak bonding stratiform country rock body tunnel
Combined Roof Cave (stratiform roof fall) and loose wall rock body back roof fall (sandbox roof fall)；
Step 2 calculates separately the kth coordinate points on mine roof for every kind of roof fall type founding mathematical models
Above-mentioned four kinds of roof fall types emit heights of roofs hidden danger value Hi, i=1,2,3,4；
1) crack country rock body back part block roof fall height hidden danger value H is being calculated1When, with reference to Fig. 1, by as follows
Mode obtains:
It establishes crack country rock body back Key Blocks to fall mechanical model, obtains back hexahedron rock
Quantitative relation formula between block neutrality height and rock mass environment, it may be assumed that
Wherein, block free face length l, block structure face inclination alpha, angle of frictionSurvey pressure coefficient lambda.
2) heights of roofs hidden danger value H is emitted in calculating butterfly plastic zone break wall rock body back2When, with reference to Fig. 2, by such as
Under type obtains:
The implicit function expression formula on circular tunnel plastic zone boundary under non-isopiestic state:
According to tunnel buried depth H, country rock bulk density γ, bidirection press ratio λ, country rock cohesive force C and internal friction angleUsing upper
Plastic zone radial depth r when formula obtains roadway surrounding rock difference θ is stated, and then according to formula:
Obtain butterfly blade profile roof fall hidden danger height H2；
3) weak bonding stratiform country rock body tunnel Combined Roof Cave height hidden danger value H is being calculated3When, with reference to Fig. 3, by such as
Under type obtains:
First according to roof strata lift height hn, distribution load qx, delamination strength σtx, mining influence coefficient fc, buried depth
Influence coefficient fm, ground stress anomally influence coefficient fd, rock stratum integrality influence coefficient fw, the stabilization of rock stratum is determined from bottom to up
Property, until encountering stable rock stratum, distinguishing rule and criterion are as follows:
When the limited span [L] of the rock stratum is less than the practical span L in tunnel, then the rock stratum unstability is caving；The pole of the rock stratum
When limiting the practical span L of span [L] greater than tunnel, and two kinds of forms: form are subdivided into according to the length for forming girder construction
One, if the mistake girder construction formed tends to clamped beam, the rock stability；Form two, if the mistake girder construction formed is between clamped beam
When between simply supported beam, judge whether the rock stratum meets the sliding instability condition and hinged sillar structure of hinged sillar structure
Deformation instability condition, meeting any instability condition, then the rock stratum unstability is caving, otherwise the rock stability；
The unstability formula of stratiform rock beam limited span:
The sliding instability formula of hinged sillar structure:
The deformation instability formula of hinged sillar structure:
Thus it obtains,Obtain H3。
4) heights of roofs hidden danger value H is emitted in calculating loose wall rock body back4When, with reference to Fig. 4, obtain as follows:
Wherein, tunnel radius α, head-room H, country rock internal friction angleIt is general
Family name's coefficient fTop。
Step 3 obtains the weighted average system of every kind of roof fall probability of happening according to the empirical probability that every kind of roof fall occurs
Number, emits heights of roofs hidden danger value H according to four kinds of roof fall typesiWith the weighted average coefficients of every kind of roof fall probability of happening, obtain
The integrated forecasting value H of heights of roofs must be emittedk, wherein i=1,2,3,4；K=1,2 ... ..., n.
Specific description utilizes fuzzy variable weight step analysis integrated forecasting fall of ground height method
1) prediction technique introduction
K-th of coordinate points of mine share i kind fall of ground type, wherein i=1, and 2,3,4；pkIt (i) is k-th of mine seat
The empirical probability that punctuate the i-th seed type fall of ground occurs, and
hkIt (i) is the i-th seed type of k-th of coordinate points of mine fall of ground Height Prediction value, wherein i=1,2,3,4；K=
1,2 ..., n；
Fk=(fij) it is the complementary type fuzzy judgment matrix that is generated by p of k-th of coordinate points of mine, wherein i=1,2,3,4；
J=1,2,3,4；K=1,2 ..., n；rkIt (i) is the fuzzy judgment matrix F of k-th of coordinate points of minekThe sum of i-th row, wherein i
=1,2,3,4；K=1,2 ..., n；
Rk=(rij) it is k-th of coordinate points of mine by the Judgement Matrix with Fuzzy Consistency gone with r construction, wherein i=1,2,
3,4；J=1,2,3,4；K=1,2 ..., n；
wkIt (i) is the corresponding weighted average coefficients of k-th of coordinate points of mine, i-th kind of fall of ground type, wherein i=1,
2,3,4；
HkFor the fall of ground high integrity predicted value of k-th of coordinate points of mine, wherein k=1,2 ..., n；
2) fuzzy variable weight index calculates step
(1) assume that the empirical probability that the i seed type fall of ground of k-th of coordinate points of mine occurs is Pk(i), wherein i=
1,2,3,4；AndIf the i-th seed type fall of ground can not be sent out
Raw, i.e. the empirical probability p that the i-th seed type fall of ground occurs is 0, does not consider this type when establishing fuzzy judgment matrix
Fall of ground.
(2) fuzzy judgment matrix is established
Wherein
(3) it to determine the corresponding weighted average coefficients of k-th of coordinate points of mine, i-th kind of fall of ground type, first to construct
Fuzzy consistent matrix, so-called Fuzzy consistent matrix R=(rij)4×4Refer to the matrix for meeting following relationship:
i)0≤rij≤ 1 (i, j=1,2,3,4)；
ii)rij+rji=1 (i, j=1,2,3,4)；
iii)
Ask row andAnd utilize conversion formulaBy fuzzy judgment matrix F=(fij)4×4
It transform Judgement Matrix with Fuzzy Consistency R=(r asij)4×4。
(4) if R is fuzzy consistent judgment matrix, weighted average coefficients can use formulaTo count
It calculates.It should be noted that working asIt needs to modify, until Wk>0.Then weighted average coefficients vector is
3) fall of ground high integrity predicted value
Acquire WkAfterwards, h is recycledk(i), so that it may the fall of ground height of k-th of coordinate points of mine be made comprehensive pre-
Measured value, value are calculated by following formula:
Heights of roofs is emitted using fuzzy variable weight analytic hierarchy process (AHP) integrated forecasting heading, by each surface drilling and
After the roof fall hidden danger rank that top plate pries through drilling determines, according to the corresponding coordinate that drills, roof fall hidden danger rank is corresponded into the coal
On mine mining engineering plan.Different roof fall ranks region is distinguished by different isopleth, the model that identical isopleth includes
Enclosing is exactly the identical top plate type of roof fall hidden danger rank, draws out heading top plate hidden danger rank block plan, with reference to Figure 10,
Figure 11.
Inventor is based on the above method and is programmed using java applet language, exploitation heading roof fall hidden danger classification
System software forms visualization interface, convenient for operating the judgement with top plate roof fall hidden danger rank.
Into after system interface, occur first interface first as shown in figure 5, listing [office on interface in toolbar
Portion's block falls], [compound roof is caving], [inbreak of butterfly plastic zone], [loose top plate leakage emits] and [classification results] is each
A module will input corresponding design parameter value according to user in different calculating parameter modules, utilize the public affairs of program editing
Formula and algorithm carry out COMPREHENSIVE CALCULATING to parameter.When calculating first module, first to input mine title, one borehole title,
Drill X the and Y actual coordinate on mining engineering plan.
Local block is clicked in system interface falls button into this information interface.It falls module for local block
It calculates, parameter free face length need to be inputted, survey and press coefficient, behind structural plane inclination angle and structural plane angle of friction, click computed push-buttom,
Back part block roof fall height H can be obtained1Value, click save button, if in input process, discovery
Mistake can click removing button, from new input.
Compound roof is clicked in system face is caving button into this information interface, as shown in Figure 6.For compound roof
Be caving module calculating, it is corresponding buy a surface drilling and pry through drilling back rock stratum sequentially input from bottom to up, single monolayer thick
Degree, lithology, lithology code, compression strength, tensile strength, elasticity modulus, body force, tunnel span, tunnel buried depth, actual measurement ground
Stress clicks computed push-buttom, and rock stratum effective thickness can be obtained, and crustal stress influences coefficient, tunnel buried depth coefficient and each rock
The limited span of layer, finally obtains tunnel lamellar composite top plate and emits heights of roofs H2Value, click save button.
Inbreak button in butterfly plastic zone is clicked in system interface into this information interface, as shown in Figure 7.For butterfly
Inbreak module in plastic zone calculates, and need to input parameter country rock internal friction angle, country rock cohesive force, tunnel buried depth, tunnel radius, volume
Power surveys pressure coefficient and butterfly leaf angle, clicks computed push-buttom, obtains back butterfly plastic zone and emit heights of roofs H3Value, click
Save button.
Loose top plate leakage is clicked in system interface emits button into this information interface, as shown in Figure 8.For loosely pushing up
Plate leakage emits module calculating, need to input parameter country rock internal friction angle, head-room, Protodyakonov coefficient, tunnel radius and safety coefficient,
Computed push-buttom is clicked, the loose top plate in tunnel is obtained and emits heights of roofs H4Value clicks save button.
Classification results button is clicked in main interface and enters this information interface, as shown in Figure 9.Computed push-buttom is clicked, is obtained
The roof fall of block top plate, Combined Roof Cave, the roof fall of butterfly leaf top plate, loose top plate roof fall fuzzy variable weight coefficient, heading emit
Heights of roofs integrated forecasting value H and classification results are clicked and generate documents button, and txt Text Mode is stored as, and click to close and press
Button carries out the typing of next drill hole information.During each column inputs parameter and calculates, saves, if carelessness omission is defeated
Enter to drill and relevant information or forget a save button, system can automatic prompt, return to modification.
Finally, it should be noted that the above embodiments are only used to illustrate the technical solution of the present invention., rather than its limitations；
Although present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that: its
It is still possible to modify the technical solutions described in the foregoing embodiments, or special to some or all of technologies
Sign is equivalently replaced；And these are modified or replaceed, the present invention that it does not separate the essence of the corresponding technical solution is each to be implemented
The range of example technical solution.
Claims (7)
1. a kind of method that integrated forecasting heading emits heights of roofs and hidden danger classification characterized by comprising
Fall of ground type i points are four kinds: crack country rock body back part block roof fall, butterfly plastic zone by step 1
The roof fall of break wall rock body back, weak bonding stratiform country rock body tunnel Combined Roof Cave and loose wall rock body back emit
Top；
Step 2 calculates separately on mine roof the above-mentioned of k-th of coordinate points for every kind of roof fall type founding mathematical models
Four kinds of roof fall types emit heights of roofs hidden danger value Hi；
Step 3 obtains the weighted average coefficients of every kind of roof fall probability of happening according to the empirical probability that every kind of roof fall occurs, according to
Four kinds of roof fall types emit heights of roofs hidden danger value HiWith the weighted average coefficients of every kind of roof fall probability of happening, acquisition, which emits, is risen
The integrated forecasting value H of degreek；
Wherein, i=1,2,3,4；K=1,2 ... ..., n.
2. the method as described in claim 1, which is characterized in that crack country rock body back part block roof fall height
Hidden danger value H1, obtain as follows: establishing crack country rock body back Key Blocks and fall mechanical model, obtain lane
Quantitative relation formula between road top plate hexahedron sillar neutrality height and rock mass environment, it may be assumed that
Wherein, block free face length l, block structure face inclination alpha, angle of frictionSurvey pressure coefficient lambda.
3. the method as described in claim 1, which is characterized in that the butterfly plastic zone break wall rock body back, which emits, to be risen
Spend hidden danger value H2, it obtains as follows:
The implicit function expression formula on circular tunnel plastic zone boundary under non-isopiestic state:
According to tunnel buried depth H, country rock bulk density γ, bidirection press ratio λ, country rock cohesive force C and internal friction angleUtilize above-mentioned public affairs
Formula obtains plastic zone radial depth r when roadway surrounding rock difference θ, and then according to formula:
Obtain butterfly blade profile roof fall hidden danger height H2。
4. the method as described in claim 1, which is characterized in that the weak bonding stratiform country rock body tunnel Combined Roof Cave is high
Spend hidden danger value H3, it obtains as follows:
First according to roof strata lift height hn, distribution load qx, delamination strength σtx, mining influence coefficient fc, buried depth influence system
Number fm, ground stress anomally influence coefficient fd, rock stratum integrality influence coefficient fw, the stability of rock stratum is determined from bottom to up, until meeting
To stablizing rock stratum, distinguishing rule and criterion are as follows:
When the limited span [L] of the rock stratum is less than the practical span L in tunnel, then the rock stratum unstability is caving；The limit of the rock stratum across
When being greater than the practical span L in tunnel away from [L], and two kinds of forms: form one are subdivided into according to the length for forming girder construction, if shape
At mistake girder construction tend to clamped beam when, the rock stability；Form two, if the mistake girder construction formed is between clamped beam and simply supported beam
Between when, judge whether the rock stratum meets the sliding instability condition of hinged sillar structure and the deformation instability item of hinged sillar structure
Part, meeting any instability condition, then the rock stratum unstability is caving, otherwise the rock stability；
The unstability formula of stratiform rock beam limited span:
The sliding instability formula of hinged sillar structure:
The deformation instability formula of hinged sillar structure:
Thus it obtains,Obtain H3。
5. the method as described in claim 1, which is characterized in that the loose wall rock body back emits heights of roofs hidden danger value H4,
It obtains as follows:
Wherein, tunnel radius α, head-room H, country rock internal friction angleProtodyakonov coefficient
fTop。
6. the method as described in claim 1, which is characterized in that hkIt (i) is k-th of coordinate points of mine the i-th seed type fall of ground
Height hidden danger value；
wkIt (i) is the corresponding weighted average coefficients of k-th of coordinate points of mine, i-th kind of fall of ground type；
HkFor the fall of ground high integrity predicted value of k-th of coordinate points of mine；
Wherein, i=1,2,3,4；K=1,2 ..., n.
7. the method as described in claim 1, which is characterized in that emitting heights of roofs and determine that it is top that it is emitted according to each coordinate points
Not, so that back is divided into the region with different roof fall ranks.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811285772.9A CN109446659B (en) | 2018-10-31 | 2018-10-31 | Method for comprehensively predicting roof fall height and hidden danger classification of coal seam roadway |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811285772.9A CN109446659B (en) | 2018-10-31 | 2018-10-31 | Method for comprehensively predicting roof fall height and hidden danger classification of coal seam roadway |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109446659A true CN109446659A (en) | 2019-03-08 |
| CN109446659B CN109446659B (en) | 2023-04-07 |
Family
ID=65550115
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811285772.9A Active CN109446659B (en) | 2018-10-31 | 2018-10-31 | Method for comprehensively predicting roof fall height and hidden danger classification of coal seam roadway |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109446659B (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110276162A (en) * | 2019-07-02 | 2019-09-24 | 内蒙古科技大学 | Mining area repeated mining fall of ground hidden danger Height Prediction method and device |
| CN111537699A (en) * | 2020-05-29 | 2020-08-14 | 煤炭科学技术研究院有限公司 | Method for distinguishing hard rock stratum of coal mine through multiple indexes |
| CN111578853A (en) * | 2020-05-27 | 2020-08-25 | 呼伦贝尔学院 | Method for predicting height of roof fall hidden danger of repeated mining roadway in mining area |
| CN111734492A (en) * | 2020-07-02 | 2020-10-02 | 中国矿业大学 | Intelligent working surface roof fall early warning method based on U-shaped area characteristic parameter extraction |
| CN113553646A (en) * | 2021-07-13 | 2021-10-26 | 中煤科工集团北京土地整治与生态修复科技研究院有限公司 | Method for judging stability of roadway in open pit mining treatment process |
| CN113553704A (en) * | 2021-07-13 | 2021-10-26 | 中煤科工集团北京土地整治与生态修复科技研究院有限公司 | Method for judging stability of roadway after completion of open pit control |
| CN115795773A (en) * | 2022-01-20 | 2023-03-14 | 山东科技大学 | Analysis method for influence factors of roof elastic energy contribution rate during coal disaster |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4437520A (en) * | 1981-06-15 | 1984-03-20 | In Situ Technology, Inc. | Method for minimizing subsidence effects during production of coal in situ |
| CN101777152A (en) * | 2010-01-28 | 2010-07-14 | 同济大学 | Life-cycle risk analysis-based cutting and tunnel scheme decision model |
| CN103924967A (en) * | 2014-03-04 | 2014-07-16 | 中国矿业大学(北京) | Roadway roof collapse hidden danger detection method |
| CN104008290A (en) * | 2014-05-27 | 2014-08-27 | 中国矿业大学(北京) | Roadway roof fall hidden danger grading method taking bidirectional anisobaric stress fields into consideration |
| CN104500090A (en) * | 2014-11-25 | 2015-04-08 | 兖州煤业股份有限公司 | Method for treating coal-mine roadway roof caving by roof filling through foam-material filler |
| CN105298523A (en) * | 2015-11-13 | 2016-02-03 | 河南理工大学 | Mining roadway roof different-level control method |
| CN106121717A (en) * | 2016-07-04 | 2016-11-16 | 河南理工大学 | A kind of driving face in coal mine caving place injecting treatment method |
| CN108062439A (en) * | 2017-12-08 | 2018-05-22 | 西安科技大学 | Roadway support quantifying design method based on plastic zone of surrounding rock size |
| CN108678807A (en) * | 2018-05-08 | 2018-10-19 | 河南理工大学 | A kind of tunnel mine pressure monitoring method for early warning |
-
2018
- 2018-10-31 CN CN201811285772.9A patent/CN109446659B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4437520A (en) * | 1981-06-15 | 1984-03-20 | In Situ Technology, Inc. | Method for minimizing subsidence effects during production of coal in situ |
| CN101777152A (en) * | 2010-01-28 | 2010-07-14 | 同济大学 | Life-cycle risk analysis-based cutting and tunnel scheme decision model |
| CN103924967A (en) * | 2014-03-04 | 2014-07-16 | 中国矿业大学(北京) | Roadway roof collapse hidden danger detection method |
| CN104008290A (en) * | 2014-05-27 | 2014-08-27 | 中国矿业大学(北京) | Roadway roof fall hidden danger grading method taking bidirectional anisobaric stress fields into consideration |
| CN104500090A (en) * | 2014-11-25 | 2015-04-08 | 兖州煤业股份有限公司 | Method for treating coal-mine roadway roof caving by roof filling through foam-material filler |
| CN105298523A (en) * | 2015-11-13 | 2016-02-03 | 河南理工大学 | Mining roadway roof different-level control method |
| CN106121717A (en) * | 2016-07-04 | 2016-11-16 | 河南理工大学 | A kind of driving face in coal mine caving place injecting treatment method |
| CN108062439A (en) * | 2017-12-08 | 2018-05-22 | 西安科技大学 | Roadway support quantifying design method based on plastic zone of surrounding rock size |
| CN108678807A (en) * | 2018-05-08 | 2018-10-19 | 河南理工大学 | A kind of tunnel mine pressure monitoring method for early warning |
Non-Patent Citations (5)
| Title |
|---|
| 吕坤等: "弱黏结复合层状顶板冒顶隐患分级方法与应用", 《矿业科学学报》 * |
| 吴颖龙等: "采空区大面积悬顶危险性预警模型研究——以陕西省榆林市榆阳区二墩煤矿为例", 《中国煤炭》 * |
| 张润来等: "基于组合赋权与有限区间云模型的巷道冒顶预测研究", 《矿冶工程》 * |
| 王永军: "自适应变权重巷道冒顶高度综合预测技术的应用", 《煤炭科技》 * |
| 赵志强等: "煤层巷道冒顶机理与预警方法", 《煤炭学报》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110276162A (en) * | 2019-07-02 | 2019-09-24 | 内蒙古科技大学 | Mining area repeated mining fall of ground hidden danger Height Prediction method and device |
| CN111578853A (en) * | 2020-05-27 | 2020-08-25 | 呼伦贝尔学院 | Method for predicting height of roof fall hidden danger of repeated mining roadway in mining area |
| CN111537699A (en) * | 2020-05-29 | 2020-08-14 | 煤炭科学技术研究院有限公司 | Method for distinguishing hard rock stratum of coal mine through multiple indexes |
| CN111734492A (en) * | 2020-07-02 | 2020-10-02 | 中国矿业大学 | Intelligent working surface roof fall early warning method based on U-shaped area characteristic parameter extraction |
| CN111734492B (en) * | 2020-07-02 | 2021-08-24 | 中国矿业大学 | Intelligent working surface roof fall early warning method based on U-shaped area characteristic parameter extraction |
| CN113553646A (en) * | 2021-07-13 | 2021-10-26 | 中煤科工集团北京土地整治与生态修复科技研究院有限公司 | Method for judging stability of roadway in open pit mining treatment process |
| CN113553704A (en) * | 2021-07-13 | 2021-10-26 | 中煤科工集团北京土地整治与生态修复科技研究院有限公司 | Method for judging stability of roadway after completion of open pit control |
| CN115795773A (en) * | 2022-01-20 | 2023-03-14 | 山东科技大学 | Analysis method for influence factors of roof elastic energy contribution rate during coal disaster |
| CN115795773B (en) * | 2022-01-20 | 2023-06-23 | 山东科技大学 | Analysis method for influence factors of roof elastic energy contribution rate in coal disaster |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109446659B (en) | 2023-04-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109446659A (en) | A kind of method that integrated forecasting heading emits heights of roofs and hidden danger classification | |
| Yang et al. | Ground subsidence and surface cracks evolution from shallow-buried close-distance multi-seam mining: a case study in Bulianta coal mine | |
| CN105260575A (en) | Roadway surrounding rock deformation predicting method based on neural network | |
| CN102999873A (en) | Digital mine safety monitoring and inspecting system | |
| CN105320817B (en) | A kind of Shallow-buried Large-span city tunnel flat ratio determines method | |
| Maruvanchery et al. | Early construction cost and time risk assessment and evaluation of large-scale underground cavern construction projects in Singapore | |
| Wu et al. | A RBFNN-based method for the prediction of the developed height of a water-conductive fractured zone for fully mechanized mining with sublevel caving | |
| CN109736796A (en) | A kind of intensification blasthole forward probe mechanism and its forecasting procedure | |
| He et al. | Subsidence prediction of overburden strata and surface based on the voussoir beam structure theory | |
| Wang et al. | Reasonable coal pillar design and remote control mining technology for highwall residual coal resources | |
| CN115830829A (en) | General forecast early warning system of colliery water damage | |
| Liang et al. | Optimization of mining method in subsea deep gold mines: A case study | |
| Mohamed et al. | Preliminary rib support requirements for solid coal ribs using a coal pillar rib rating (CPRR) | |
| Coşar | Application of rock mass classification systems for future support design of the Dim Tunnel near Alanya | |
| CN106709988A (en) | Construction method of engineering geological section map | |
| Sheng et al. | Investigation of the optimization of unloading mining scheme in large deep deposit based on vague set theory and its application | |
| Klemetti et al. | Longwall gateroad yield pillar response and model verification–a case study | |
| Guo et al. | Analysis and prediction of nonuniform deformation in composite strata during tunnel excavation | |
| Guo et al. | A Three-Dimensional Supporting Technology, Optimization and Inspiration from a Deep Coal Mine in China | |
| Liu et al. | Stope structure evaluation based on the damage model driven by microseismic data and Mathews stability diagram method in Xiadian Gold Mine | |
| Zhang et al. | Disaster-causing mechanism of extremely thick igneous rock induced by mining and prevention method by backfill mining | |
| Yang et al. | Predicting the height of the water-conducting fractured zone based on a multiple regression model and information entropy in the Northern Ordos Basin, China | |
| Jia et al. | Stability analysis of shallow goaf based on field monitoring and numerical simulation: A case study at an Open-Pit Iron Mine, China | |
| Liu et al. | Mining method selection and optimization for hanging-wall ore-body at Yanqianshan Iron Mine, China | |
| CN115236740A (en) | Water-flowing fractured zone development height prediction method based on seismic interpretation achievement |
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 |