CN103670511B - A kind of analysis calculation method of roof strata level of breakage - Google Patents
A kind of analysis calculation method of roof strata level of breakage Download PDFInfo
- Publication number
- CN103670511B CN103670511B CN201310490614.8A CN201310490614A CN103670511B CN 103670511 B CN103670511 B CN 103670511B CN 201310490614 A CN201310490614 A CN 201310490614A CN 103670511 B CN103670511 B CN 103670511B
- Authority
- CN
- China
- Prior art keywords
- plane
- area
- height
- breakage
- roof strata
- 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.)
- Active
Links
Landscapes
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a kind of analysis calculation method of roof strata level of breakage, analytical procedure is as follows: after choosing seam mining, roof strata is caving height <i>h</iGreatT.Gr eaT.GT
collapse, the expansion development height <i>h</iGreatT.Gr eaT.GT of crack in rock stratum
split, caving zone area of plane <i>m</iGreatT.Gr eaT.GT
collapsewith the area of plane <i>m</iGreatT.Gr eaT.GT of top board shelly formation
splitfor evaluation index; Build the membership function of each evaluation index, after comprehensive analysis tentative calculation, the weight of selected each index, calculates roof strata level of breakage total degree of membership <i> η </i>; The expansion developmental state of rock stratum internal fissure is evaluated with <i> η </i> value, <i> η </i> value is larger, rock breakdown degree is larger, rock stratum internal fissure is just educated all the more, otherwise less; Set up the analysis and calculation model of roof strata level of breakage.The invention process is simple, according to clear and definite, and the level of breakage of the roof strata after accurate evaluation seam mining.
Description
(1) technical field
The present invention relates to mining engineering field, particularly a kind of analysis calculation method of roof strata level of breakage.
(2) background technology
Coal resources in China is abundant especially, and not only in Plain and mountain region etc., regional underground is composed and deposited a large amount of coal resources, and also composes under many water bodys or aquifer and deposit a large amount of coal resources.The safety of seam mining under water body or aquifer, key depends on the level of breakage of roof strata.For a long time, China mines under carrying out water body or aquifer always and have accumulated certain experience.But about the design formulas of safety pillar under wa ter bodies in Main Basis " building, water body, railway and main roadway coal column stay and establish and mining pressed coal code " of mining under China's water body or aquifer at present, and only considered this single factor of height of leaking crevice belt, not from the level of breakage of assay roof strata the scope of breaking, seam mining whether safety can not to be reflected under water body or aquifer accurately, bring danger to the exploitation in coal seam under water body or aquifer, time serious, even cause water leak accident.
(3) summary of the invention
The present invention, in order to make up the deficiencies in the prior art, provides a kind of analysis rationally, the analysis calculation method of the roof strata level of breakage of convenience of calculation.
The present invention is achieved through the following technical solutions:
An analysis calculation method for roof strata level of breakage, is characterized in that: analytical procedure is as follows:
1) after, choosing seam mining, roof strata is caving height, the area of plane of the expansion development height of crack in rock stratum, caving zone and the area of plane four influence factors of top board shelly formation evaluation index as roof strata level of breakage;
2) field measurement, according to after seam mining, probing, physical prospecting or micro-seismic monitoring data, determine to be caving height h
collapse, the expansion development height h of crack in rock stratum
split, be caving highly for h
collapsetime the area of plane m of caving zone
collapsebe h with the rupture height of roof strata
split 'time the area of plane m of top board shelly formation
split;
3), structure is caving height membership function α, the expansion development height membership function β of crack in rock stratum, the area of plane membership function γ of caving zone and the area of plane membership function θ of top board shelly formation respectively;
Be caving height membership function:
In formula: α is for being caving height degree of membership; h
collapsefor being caving height, m; H is roof to the minimum bedrock thickness of aquifer or water body, m.
The expansion development height membership function of crack in rock stratum:
In formula: β is the expansion development height degree of membership of crack in rock stratum; h
splitfor the expansion development height of crack in rock stratum, m; H is roof to the minimum bedrock thickness of aquifer or water body, m.
The area of plane membership function of caving zone:
In formula: γ is the area of plane degree of membership of caving zone; m
collapsefor being caving highly for h
collapsetime the area of plane, m
2; M is the area of plane in goaf, m
2.
The area of plane membership function of top board shelly formation:
In formula: θ is the area of plane degree of membership of top board shelly formation; m
splitrupture height for roof strata is h
split 'time the area of plane, m
2; M is the area of plane in goaf, m
2.
For the level of breakage of overall merit roof strata, superposition calculation is carried out to α, β, γ and θ:
In formula: η is the total degree of membership of roof strata level of breakage; B is the weight of α; C is the weight of β, and d is the weight of γ; E is the weight of θ.B, c, d and e directly choose according to the actual geological condition of achievement in research in the past and mine, carry out multiple stacking tentative calculation, determine after comprehensive analysis;
4), after, comprehensively analyzing each weight selected, the total degree of membership η of roof strata level of breakage is calculated;
5), evaluate the expansion developmental state of overlying rock internal fissure with the size of the total degree of membership η of roof strata level of breakage, η value is larger, and roof strata level of breakage is larger, and overlying rock internal fissure is just educated all the more, and the height in crack is also higher, otherwise less;
6) analysis and calculation model of roof strata level of breakage, is set up.
The invention has the beneficial effects as follows: the analysis calculation method of roof strata level of breakage, from the level of breakage that spatially have rated roof strata, can the situation that be caving of top board overlying rock and the expansion developmental state in crack after accurate analysis seam mining, for under water body or aquifer, the safe working in coal seam provides the foundation of science, be suitable for wide popularization and application.
Detailed description of the invention
For more clearly describing technical scheme of the present invention, advantage and practicality, complete description is carried out to the specific embodiment of the present invention below:
The analysis calculation method of roof strata level of breakage, is characterized in that: analytical procedure is as follows:
1) after, choosing seam mining, roof strata is caving height, the area of plane of the expansion development height of crack in rock stratum, caving zone and the area of plane four influence factors of top board shelly formation evaluation index as roof strata level of breakage;
2) field measurement, according to after seam mining, probing, physical prospecting or micro-seismic monitoring data, determine to be caving height h
collapse, the expansion development height h of crack in rock stratum
split, be caving highly for h
collapsetime the area of plane m of caving zone
collapsebe h with the rupture height of roof strata
split 'time the area of plane m of top board shelly formation
split;
3), structure is caving height membership function α, the expansion development height membership function β of crack in rock stratum, the area of plane membership function γ of caving zone and the area of plane membership function θ of top board shelly formation respectively;
Be caving height membership function:
In formula: α is for being caving height degree of membership; h
collapsefor being caving height, m; H is roof to the minimum bedrock thickness of aquifer or water body, m.
The expansion development height membership function of crack in rock stratum:
In formula: β is the expansion development height degree of membership of crack in rock stratum; h
splitfor the expansion development height of crack in rock stratum, m; H is roof to the minimum bedrock thickness of aquifer or water body, m.
The area of plane membership function of caving zone:
In formula: γ is the area of plane degree of membership of caving zone; m
collapsefor being caving highly for h
collapsetime the area of plane, m
2; M is the area of plane in goaf, m
2.
The area of plane membership function of top board shelly formation:
In formula: θ is the area of plane degree of membership of top board shelly formation; m
splitrupture height for roof strata is h
split 'time the area of plane, m
2; M is the area of plane in goaf, m
2.
For the level of breakage of overall merit roof strata, superposition calculation is carried out to α, β, γ and θ:
In formula: η is the total degree of membership of roof strata level of breakage; B is the weight of α; C is the weight of β, and d is the weight of γ; E is the weight of θ.B, c, d and e directly choose according to the actual geological condition of achievement in research in the past and mine, carry out multiple stacking tentative calculation, determine after comprehensive analysis;
4), after, comprehensively analyzing each weight selected, the total degree of membership η of roof strata level of breakage is calculated;
5), evaluate the expansion developmental state of overlying rock internal fissure with the size of the total degree of membership η of roof strata level of breakage, η value is larger, and roof strata level of breakage is larger, and overlying rock internal fissure is just educated all the more, and the height in crack is also higher, otherwise less;
6) analysis and calculation model of roof strata level of breakage, is set up.
Adopt the analysis calculation method of roof strata level of breakage of the present invention, reflect the expansion developmental state of roof strata internal fissure more accurately.According to the concrete evaluation result of this patent, the safety of exploiting under can accurately evaluating water body or aquifer.Compared with former achievement in research, this patent considers more influence factor, makes evaluation result have more science.
Claims (1)
1. an analysis calculation method for roof strata level of breakage, is characterized in that: analytical procedure is as follows:
1) after, choosing seam mining, roof strata is caving height, the area of plane of the expansion development height of crack in rock stratum, caving zone and the area of plane four influence factors of top board shelly formation evaluation index as roof strata level of breakage;
2) field measurement, according to after seam mining, probing, physical prospecting or micro-seismic monitoring data, determine to be caving height h
collapse, the expansion development height h of crack in rock stratum
split, be caving highly for h
collapsetime the area of plane m of caving zone
collapsebe h with the rupture height of roof strata
split 'time the area of plane m of top board shelly formation
split;
3), structure is caving height membership function α, the expansion development height membership function β of crack in rock stratum, the area of plane membership function γ of caving zone and the area of plane membership function θ of top board shelly formation respectively;
Be caving height membership function:
In formula: α is for being caving height degree of membership; h
collapsefor being caving height, m; H is roof to the minimum bedrock thickness of aquifer or water body, m;
The expansion development height membership function of crack in rock stratum:
In formula: β is the expansion development height degree of membership of crack in rock stratum; h
splitfor the expansion development height of crack in rock stratum, m; H is roof to the minimum bedrock thickness of aquifer or water body, m;
The area of plane membership function of caving zone:
In formula: γ is the area of plane degree of membership of caving zone; m
collapsefor being caving highly for h
collapsetime the area of plane, m
2; M is the area of plane in goaf, m
2;
The area of plane membership function of top board shelly formation:
In formula: θ is the area of plane degree of membership of top board shelly formation; m
splitrupture height for roof strata is h
split 'time the area of plane, m
2; M is the area of plane in goaf, m
2;
For the level of breakage of overall merit roof strata, superposition calculation is carried out to α, β, γ and θ:
In formula: η is the total degree of membership of roof strata level of breakage; B is the weight of α; C is the weight of β, and d is the weight of γ; E is the weight of θ;
B, c, d and e directly choose according to the actual geological condition of achievement in research in the past and mine, carry out multiple stacking tentative calculation, determine after comprehensive analysis;
4), after, comprehensively analyzing each weight selected, the total degree of membership η of roof strata level of breakage is calculated;
5), evaluate the expansion developmental state of overlying rock internal fissure with the size of the total degree of membership η of roof strata level of breakage, η value is larger, and roof strata level of breakage is larger, and overlying rock internal fissure is just educated all the more, and the height in crack is also higher, otherwise less;
6) analysis and calculation model of roof strata level of breakage, is set up.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310490614.8A CN103670511B (en) | 2013-10-19 | 2013-10-19 | A kind of analysis calculation method of roof strata level of breakage |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310490614.8A CN103670511B (en) | 2013-10-19 | 2013-10-19 | A kind of analysis calculation method of roof strata level of breakage |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103670511A CN103670511A (en) | 2014-03-26 |
CN103670511B true CN103670511B (en) | 2015-12-09 |
Family
ID=50309200
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310490614.8A Active CN103670511B (en) | 2013-10-19 | 2013-10-19 | A kind of analysis calculation method of roof strata level of breakage |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103670511B (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105804736B (en) * | 2016-05-16 | 2019-03-15 | 西安科技大学 | A kind of calculation method of flat seam base object model rock stratum limited span |
CN108446450B (en) * | 2018-02-23 | 2021-07-02 | 通用技术集团工程设计有限公司 | Method for analyzing and calculating damage degree of building influenced by mining |
CN108665136A (en) * | 2018-02-24 | 2018-10-16 | 煤炭工业济南设计研究院有限公司 | A kind of analysis calculation method of filling mining effect assessment |
CN113008307B (en) * | 2021-03-09 | 2021-11-30 | 华侨大学 | Method for determining fracture intervals of equal-interval fractured rock stratum |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101881182B (en) * | 2010-07-14 | 2011-10-26 | 中国矿业大学(北京) | Simulation and forecast evaluation method for delaying water inrush by mine fracture structure |
CN103089305A (en) * | 2013-01-14 | 2013-05-08 | 中国矿业大学(北京) | Risk assessment method for coal and gas outburst of coal beds |
-
2013
- 2013-10-19 CN CN201310490614.8A patent/CN103670511B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN103670511A (en) | 2014-03-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103591922B (en) | Short range seam mining upper Seam Floor Failure degree of depth radon gas detection method | |
CN104481587A (en) | Large-mining depth and long-span fully-mechanized top-coal caving face roof sandstone fracture water detecting and preventing method | |
CN102799955B (en) | Water bursting coefficient is less than 0.06MPa/m district Water Inrush and evaluates three figure methods | |
CN103993878A (en) | Water diversion fracture zone height estimating method | |
CN103161434A (en) | Mining method for low permeability reservoir of shale gas and the like | |
CN103670511B (en) | A kind of analysis calculation method of roof strata level of breakage | |
CN114046178A (en) | Coal mine underground water prevention and control early warning system, method and application | |
Meng et al. | In situ investigation and numerical simulation of the failure depth of an inclined coal seam floor: a case study | |
CN111379562B (en) | Water-controlled coal mining method and device under composite water body | |
CN103606019A (en) | Mine goaf overlying stratum sedimentation dynamic prediction method based on time-space relationship | |
Qiu et al. | Characteristics of vertical karst development and grouting reinforcement engineering practice of the Ordovician top in the Feicheng coalfield, China | |
Zhang et al. | Disaster-causing mechanism of extremely thick igneous rock induced by mining and prevention method by backfill mining | |
CN111680896B (en) | Coal mine underground reservoir safety distance determining method | |
Xu et al. | Water disaster types and water control measures of Hanxing coal mine area | |
Wei et al. | A case study of damage energy analysis and an early warning by microseismic monitoring for large area roof caving in shallow depth seams | |
CN101561359B (en) | Grading method for gas tunnel | |
Luan et al. | A Case Study on Development of Water-flow Fractured Zone During Close-Distance Repeated Coal Mining | |
Wang | Analysis of the karst development law based on multiple exploration technologies of Cambrian limestone | |
Zhang et al. | Detection and evaluation of crack development near the fault zone under the influence of coal mining | |
JINLONG et al. | 3D comprehensive advance exploration and prevention technology of mine water-inrush. | |
Chang et al. | Control of water inrush from longwall floor aquifers using a division paste backfilling method | |
Wang et al. | Evaluation of Hydraulic Characteristics of the Normal Fault Based on Key Data | |
Shi et al. | Study on the Water Source for Filling the Weak Rich (Conductive) Water Aquifer in the Changxing Formation Limestone of Coal Seam Roof | |
Cheng et al. | Microseismic monitoring and numerical simulation research of floor failure depth in extra-thick coal seam | |
Geng et al. | Risk assessment and cause analysis of water inrush in thin bedrock areas with shallow coal seams for safe, water-conserving mining |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CP01 | Change in the name or title of a patent holder |
Address after: 250031 Shandong province Ji'nan City flyover District Dikou Road No. 141 Patentee after: General Technology Group Engineering Design Co., Ltd. Address before: 250031 Shandong province Ji'nan City flyover District Dikou Road No. 141 Patentee before: Coal Industry Jinan Design & Research Co., Ltd. |
|
CP01 | Change in the name or title of a patent holder |