CN115345372A - Surface subsidence prediction method for deformation area control of coal and grain composite area - Google Patents
Surface subsidence prediction method for deformation area control of coal and grain composite area Download PDFInfo
- Publication number
- CN115345372A CN115345372A CN202211002816.9A CN202211002816A CN115345372A CN 115345372 A CN115345372 A CN 115345372A CN 202211002816 A CN202211002816 A CN 202211002816A CN 115345372 A CN115345372 A CN 115345372A
- Authority
- CN
- China
- Prior art keywords
- mining
- area
- filling
- coal
- subsidence
- 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
- 238000000034 method Methods 0.000 title claims abstract description 67
- 239000003245 coal Substances 0.000 title claims abstract description 56
- 239000002131 composite material Substances 0.000 title claims abstract description 36
- 238000005065 mining Methods 0.000 claims abstract description 251
- 238000004364 calculation method Methods 0.000 claims abstract description 18
- 238000005259 measurement Methods 0.000 claims description 8
- 239000000463 material Substances 0.000 claims description 5
- 238000005553 drilling Methods 0.000 claims description 4
- 230000000694 effects Effects 0.000 claims description 3
- 238000011156 evaluation Methods 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 claims description 3
- 210000003127 knee Anatomy 0.000 claims description 2
- 238000004519 manufacturing process Methods 0.000 claims description 2
- 238000000605 extraction Methods 0.000 claims 1
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000002195 synergetic effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Tourism & Hospitality (AREA)
- Development Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- General Business, Economics & Management (AREA)
- Game Theory and Decision Science (AREA)
- Marketing (AREA)
- Quality & Reliability (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a method for predicting surface subsidence for deformation area control of a coal and grain composite area, which is suitable for surface subsidence prediction. The method comprises the steps of firstly determining address information of a coal seam, then determining caving face information, filling mining face information and backfill filling information according to a nonlinear relation between mining degree and a sinking coefficient, calculating equivalent mining height of a filling working face based on an equivalent mining height theory, and calculating surface subsidence, surface horizontal movement, surface inclination, surface horizontal deformation and surface curvature caused by a mining range by utilizing a probability integral model, so that surface subsidence prediction of deformation area control facing to cultivated land protection of a coal and grain composite area is realized. The method has the advantages of simple calculation steps, accurate calculation results, more accordance with actual conditions and important practical value for underground coal mining and accurate protection of ground cultivated land.
Description
Technical Field
The invention relates to a method for predicting surface subsidence controlled by a deformation area, which is suitable for a method for predicting surface subsidence controlled by a deformation area oriented to a coal and grain composite area and belongs to the field of surface subsidence prediction.
Technical Field
However, because underground coal mining often damages surface buildings (structures), cultivated land and the like, and the conventional methods for reducing the surface subsidence, such as filling mining and strip mining, have the problems of high cost, large resource loss and the like, and in consideration of the beneficial factors of high deformation tolerance of the cultivated land, large allowable deformation value and the like, new methods for controlling coal mining in deformation areas for cultivated land protection of coal-grain composite areas are gradually developed, and the prediction of the surface subsidence becomes an urgent problem to be solved.
The prediction of surface subsidence is related to many factors including mining thickness, coal seam dip, mining depth, prediction method and prediction parameters, etc., where a scientifically reasonable prediction method relates to whether the final prediction result is accurate. Although many scholars currently carry out research aiming at the prediction of surface subsidence in conventional caving mining, filling mining and strip mining, no surface subsidence prediction method for deformation area control of coal-grain composite area cultivated land protection exists at present.
Disclosure of Invention
Aiming at the technical problems, the earth surface subsidence prediction method for the deformation area control of the coal and grain composite area is simple in calculation steps, high in accuracy and easy to obtain calculation parameters.
In order to achieve the technical purpose, the invention provides a surface subsidence prediction method for deformation area control of a coal and grain composite area, which comprises the following steps:
firstly, determining underground coal seam parameters of a coal-grain composite area through geological mining data;
then according to the nonlinear relation between the mining degree and the sinking coefficient, a mining method of mining while mining and filling is used, mining of the caving face is carried out after filling of the filling face on the two sides is completed, the mining width b of the caving face is determined, and then according to the width a and the filling rate of the filling mining face
The number of the caving working faces is n, and the number of the filling working faces is n-1;
then calculating the equivalent mining height m of the filling working face according to the equivalent mining height theory 1 ；
The mining range R of the coal seam is divided into an upper part and a lower part by the section of the filling working face, the upper part of the section is equivalent to the sinking space of the top plate of the whole filling mining working face, and the mining range is recorded as R 1 The lower half part of the profile is a caving space below a sinking top plate, and the mining range is recorded as R 2 Wherein the production range is R 1 All are in a mining empty state, and the mining range is R 2 Then reserving unfilled space behind a filling body made of filling materials for the working face at intervals, wherein the width of the filling body is equal to the width a of the filling mining working face; wherein the lower half part of the profile is the sinking height of the caving space under the sinking top plate
Filling rate of filling body formed by filling working face
Influence;
respectively calculating the mining range R by utilizing a probability integral model 1 And R 2 Induced surface subsidence W 1 And W 2 Horizontally moving U 1 And U 2 Inclination I 1 And I 2 Horizontal deformation epsilon 1 And epsilon 2 And curvature K 1 And K 2 And after correspondingly superposing the two groups of deformation values, finally estimating the parameters for controlling the surface deformation area: surface subsidence W, horizontal movement U, inclination I, horizontal deformation epsilon, and curvature K.
Further, the section mining range R = R of the coal seam 1 +nR 2 Wherein R is 1 ＝[nb+(n-1)a]m 1 ，
m 1 Is the sinking height/equivalent mining height of the top plate,
the height of the filling body after being compacted by the top plate.
The method comprises the following specific steps:
step 1: aiming at a coal and grain composite region needing precise protection of ground cultivated land, collecting and researching geological mining data, a drilling histogram and other information of the region to be mined, and determining mining depth H, mining thickness m and coal seam inclination angle alpha information of a coal seam to be mined;
and 2, step: establishing a nonlinear relation between mining degree and subsidence coefficient according to surface subsidence monitoring data mined by an underground working face in the mining area or similar geological mining conditions, determining mining width b of the caving working face based on the nonlinear relation, and designing width a and filling rate of the filling mining working face according to the area range to be mined, ground farmland protection conditions, filling methods and filling material information
Designing caving working faces and filling mining working faces at intervals, wherein the number of the caving working faces is n, and the number of the filling mining working faces is n-1;
and step 3: according to mining thickness m and filling rate
Calculating the equivalent mining height of the filling mining working face according to the equivalent mining height theory
Thus, the section mining range R for controlling the surface deformation region for the farmland protection of the coal and grain composite region is regarded as the length nb + (n-1) a and the height m 1 Mining area R 1 Plus n lengths b and heights
Mining area R 2 ；
And 4, step 4: and (3) actual measurement predicted parameters of a large amount of existing surface subsidences of the reference mining area: a sinking coefficient q, a horizontal movement coefficient B, a main influence angle tangent tan beta, a knee offset S and a mining influence propagation angle theta are respectively determined to be mining ranges R 1 And a mining range R 2 The probability integral prediction parameter of (2): coefficient of subsidence q 1 、q 2 Coefficient of horizontal movement B 1 、B 2 The main influence of the tangent tan beta 1 、tanβ 2 Inflection offset S 1 、S 2 Angle of propagation of mining effect theta 1 、θ 2 ；
And 5: calculating the mining range R by using a probability integral model according to the subsidence prediction parameters 1 Induced surface subsidence W 1 Horizontal movement of the earth's surface U 1 Surface inclination I 1 Horizontal deformation of earth's surface epsilon 1 And surface curvature K 1 The mining range R can also be calculated 2 Induced subsidence of the earth's surface W 2 Horizontal movement of the earth's surface U 2 Surface inclination I 2 Horizontal deformation of earth's surface epsilon 2 And surface curvature K 2 Surface subsidence controlled by surface deformation area W = W 1 +W 2 Horizontal ground surface movement U = U 1 +U 2 Surface inclination I = I 1 +I 2 Surface horizontal deformation epsilon = epsilon 1 +ε 2 And surface curvature K = K 1 +K 2 And further realize the coal and grain oriented composite areaSurface subsidence prediction for deformation zone control for farmland protection.
Further, the mining area R 1 The sinking coefficient of (A) is: the actual measurement parameters under the mining condition of the mining area thin coal seam caving method are adopted according to the equivalent mining height theory, wherein the sinking coefficient is increased by 5-10 percent, and the method specifically comprises the following steps:
q 1 ＝q(1+x)
in the formula: q. q.s 1 For the exploitation range R based on the equivalent mining height theory 1 A sinking coefficient; q is the mining subsidence coefficient of the caving method of the mining area; the value range of x is between 0.05 and 0.1;
mining area R 1 Horizontal movement coefficient b of 1 = b, main influence angle tangent tan β 1 = tan β, inflection offset S 1 ＝S。
Further, the mining area R 2 Through the sinking coefficient of q 2 The/q is expressed as:
in the formula: q. q.s 2 To the mining area R 2 The sinking coefficient of (a); q is mining subsidence coefficient by mining area caving method; m is the mining area R 2 Thickness of mining
The unit m; b is the mining area R 2 The width of mining of (1) in units of m; a is the width of the filling working face, namely the width of the remaining width, and the unit m; h is the depth of cut in m.
Further, the mining area R 2 Horizontal movement coefficient B of 2 The calculation method comprises the following steps:
in the formula: b 2 To the mining area R 2 A horizontal movement coefficient of (a); b is mining horizontal movement coefficient of mining area caving method; b is the mining area R 2 The width of mining (1) is in unit m; a is left wideFilling face width, unit m; h is the depth of cut in m.
Further, the mining area R 2 Has a major influence on the tangent tan beta 2 The calculation process of (2):
tanβ 2 /tanβ＝0.7847e -0.0012PH
in the formula: tan beta 2 To the mining area R 2 The main influence of (1) on the tangent; tan beta is the angle tangent mainly influenced by mining by a mining area caving method; p is an overburden comprehensive evaluation coefficient, and each mine area has a parameter value thereof; h is the depth of cut in m.
Further, inflection offset S 2 The calculation process of (2):
in the formula: s 2 To the mining area R 2 (iv) inflection offset of; b is the mining area R 2 The width of mining (1) is in unit m; a is the width of the filling working face, namely the width of the remaining width, and the unit m; h is the depth of cut in m.
Further, the mining influence propagation angle theta is the mining influence propagation angle in the mining process of all caving methods under the similar geological conditions.
Has the beneficial effects that:
the method considers the filling body compression characteristic of deformation area control of the coal and grain composite area, gives consideration to surface subsidence control and farmland protection, scientifically and reasonably provides a surface subsidence prediction method for deformation area control of farmland protection of the coal and grain composite area, fills the blank of the surface subsidence prediction method defect of green mining of regional subsidence control coal of the coal and grain composite area, provides technical support for realizing the synergistic development of green mining and farmland protection of coal of the coal and grain composite area, and has important practical significance and application value for realizing regional surface subsidence control, farmland protection and the like.
The method is based on the prior art, and by scientific and reasonable analysis, the equivalent mining height m of the equivalent mining height theory is referred 1 And compacted thickness of the filling body
Innovative underground coal mining space suitable for controlling coal mining in farmland protection deformation area of coal and grain composite area is divided into R 1 And R 2 Then, R is respectively determined according to the existing subsidence prediction parameters 1 、R 2 The prediction parameters are independently calculated by adopting a probability integral model, and finally, corresponding calculation results are superposed to obtain the movement and deformation of the ground surface caused by the whole mining space, so that the method is particularly suitable for calculating the subsidence of the ground surface when the deformation of the coal and grain composite area region controls the farmland protection.
Drawings
FIG. 1 is a schematic diagram of calculation of the surface subsidence prediction method for deformation area control of a coal and grain composite area.
Detailed Description
The invention will be further described in detail with reference to the figures and the specific implementation process:
as shown in fig. 1, the invention relates to a surface subsidence prediction method for deformation area control of a coal-grain composite area, which comprises the following steps:
firstly, determining coal seam parameters through geological mining data; then determining the mining width b of the caving face according to the nonlinear relation between the mining degree and the sinking coefficient; then designing the width a and the filling rate of the filling mining working face
And the number n of the caving working faces, and then the number n-1 of the filling working faces can be determined; then calculating the equivalent mining height m of the filling working face according to the equivalent mining height theory 1 (ii) a Finally, the profile mining area R can be regarded as the length nb + (n-1) a and the height m 1 Mining area R 1 Plus n lengths b and heights
Mining area R 2 . Then, the probability integral model can be used for respectively calculating the mining range and the caused surface subsidence, horizontal movement, inclination, horizontal deformation and curvature, and the surface deformation area control can be obtained through superposition calculationSurface subsidence W, horizontal movement U, inclination I, horizontal deformation epsilon, and curvature K.
The method comprises the following specific steps:
step 1: aiming at a coal and grain composite region needing precise protection of a ground farmland, collecting and researching geological mining data, a drilling histogram and other information of the region to be mined, and determining mining depth H, mining thickness m, coal bed inclination angle alpha and other information of a coal bed to be mined.
Step 2: establishing a nonlinear relation between mining degree and subsidence coefficient according to surface subsidence monitoring data mined by an underground working face in the mining area or similar geological mining conditions, determining mining width b of a caving working face based on the nonlinear relation, and designing filling mining working face width a and filling rate according to information such as the range of a region to be mined, ground cultivated land protection conditions, filling methods and filling materials
And simultaneously determining the number n of the caving working faces, and then determining the number of the filling mining working faces as n-1.
And step 3: according to the mining thickness m determined in the step 1 and the filling rate determined in the step 2
The equivalent mining height of the filling mining working face can be calculated according to the equivalent mining height theory
Thus, the section mining range R for controlling the surface deformation region for the farmland protection of the coal and grain composite region is regarded as the length nb + (n-1) a and the height m 1 Mining area R 1 Plus n lengths b and heights
Mining area R 2 。
And 4, step 4: respectively determining mining ranges R by referring to actual measurement predicted parameters of a large number of existing surface subsidences in a mining area 1 And a mining range R 2 Probability integral prediction parameter of (2): coefficient of subsidence q, horizontal movement systemNumber B, main influence tangent tan β, inflection offset S, mining influence propagation angle θ, as follows:
(1) Mining area R 1 Is a length nb + (n-1) a and a height
The estimated parameters can refer to the equivalent mining height theory, and the actual measurement parameters under the mining conditions of the mining area thin coal seam caving method are adopted, wherein the sinking coefficient should be increased by 5-10 percent, and the method specifically comprises the following steps:
q 1 ＝q(1+x)
in the formula: q. q of 1 For the exploitation range R based on the equivalent mining height theory 1 A sinking coefficient; q is the mining subsidence coefficient of the caving method of the mining area; the value range of x is between 0.05 and 0.1.
(2) Mining area R 2 N pieces of length b and height
The mining area R is formed because the filling body of the filling face is compacted 2 The predicted parameters can be determined by referring to the mining of a strip with the mining width b and the reserved width a, and the specific calculation formula is as follows:
1) And (3) calculating a sinking coefficient:
in the formula: q. q.s 2 To the mining area R 2 The sinking coefficient of (a); q is mining subsidence coefficient by mining area caving method; m is the mining area R 2 Thickness of mining
The unit m; b is the mining area R 2 The width of mining of (1) in units of m; a is the width of the filling working face, namely the width of the filling working face, and the unit is m; h is the depth of cut in m.
2) Calculation of horizontal movement coefficient:
in the formula: b 2 To the mining area R 2 A horizontal movement coefficient of; b is mining horizontal movement coefficient of mining area caving method; b is the mining area R 2 The width of mining of (1) in units of m; a is the width of the filling working face, namely the width of the filling working face, and the unit is m; h is the depth of cut in m.
3) Calculation of the dominant influence tangent
tanβ 2 /tanβ＝0.7847e -0.0012PH
In the formula: tan beta 2 To the mining area R 2 The main influence of (c) on the tangent; tan beta is the angle tangent mainly influenced by mining by a mining area caving method; p is an overburden comprehensive evaluation coefficient, and each mine area has a parameter value thereof; h is the depth of cut in m.
4) Calculation of inflection offset
In the formula: s. the 2 To the mining area R 2 Inflection offset of (2); b is the mining area R 2 The width of mining of (1) in units of m; a is the width of the filling working face, namely the width of the filling working face, and the unit is m; h is the depth of cut in m.
5) Calculation of mining impact propagation angle
The mining influence propagation angle can be similar to that of mining by a total caving method under geological conditions.
And 5: calculating a mining range S by using a probability integral model according to the subsidence prediction parameters calculated in the step 4 1 Induced surface subsidence W 1 Horizontal movement of the earth's surface U 1 Surface inclination I 1 Horizontal deformation of earth's surface epsilon 1 And surface curvature K 1 The mining range S can also be calculated 2 Induced subsidence of the earth's surface W 2 The earth's surface moves horizontally U 2 Surface inclination I 2 Horizontal deformation of earth's surface epsilon 2 And surface curvature K 2 So that the deformation zone of the earth surface can be calculatedDomain controlled surface subsidence W = W 1 +W 2 Horizontal ground surface movement U = U 1 +U 2 Surface inclination I = I 1 +I 2 Horizontal deformation of earth surface epsilon = epsilon 1 +ε 2 And surface curvature K = K 1 +K 2 And further, the surface subsidence prediction of deformation area control for the farmland protection of the coal-grain composite area is realized.
The first embodiment,
Step 1: determining coal seam parameters of the mining area according to geological mining data, a drilling histogram and other information of the area to be mined, and acquiring the thickness m =3m of the coal seam, the inclination angle alpha =0 ° of the coal seam, and the mining depth H =500m.
And 2, step: determining the width b =180m of the caving face according to the nonlinear relation between the mining degree and the sinking coefficient established by the actual measurement data of the mining area; design filling face width a =150m, filling rate
And if the number of the collapsed working faces is n =3, the number of the filled working faces is n-1=2.
And step 3: calculating equivalent mining height of filling working face based on equivalent mining height theory
The profile mining range R can then be regarded as a length nb + (n-1) a =840m and a height m 1 Mining range R of =0.6m 1 Plus n =3 with a length of b =180m and a height of
Mining area R 2 。
And 4, step 4: and (3) obtaining the expected subsidence parameters of the whole caving method of the mining area by referring to the actual measurement expected parameters of the mining area: the sinking coefficient q =0.5, and the horizontal movement coefficient B =0.27; the main influence angle tangent: tan β =2.3; inflection point offset S =0.05h =25m, mining impact propagation angle θ =90 °; respectively calculating according to a calculation formulaOut of the range of exploitation R 1 Sinkage prediction parameters of (1): q. q.s 1 ＝q(1+0.05)＝0.525，B 1 ＝0.27，tanβ 1 ＝2.3，S 1 ＝25m，θ 1 =90 ° and mining range R 2 Sinkage prediction parameters of (1): q. q.s 2 ＝0.366，B 2 ＝0.188，tanβ 2 ＝1.6，S 2 ＝24.6m，θ 2 ＝90°。
And 5: substituting the subsidence prediction parameters calculated in the step 4 into a probability integral model to calculate the mining range S 1 Induced subsidence of the earth's surface W 1 Horizontal movement of the earth's surface U 1 Surface inclination I 1 Horizontal deformation of earth's surface epsilon 1 And surface curvature K 1 The mining range S can also be calculated 2 Induced subsidence of the earth's surface W 2 The earth's surface moves horizontally U 2 Surface inclination I 2 Horizontal deformation of earth's surface epsilon 2 And surface curvature K 2 Thus, the surface subsidence W = W controlled by the surface deformation area can be calculated 1 +W 2 Horizontal ground surface movement U = U 1 +U 2 Surface inclination I = I 1 +I 2 Surface horizontal deformation epsilon = epsilon 1 +ε 2 And surface curvature K = K 1 +K 2 And further, the surface subsidence prediction of deformation area control for the farmland protection of the coal-grain composite area is realized. In this embodiment, the extreme values of the earth surface movement and deformation controlled by the earth surface deformation region calculated by the probability integral model are as follows: w max ＝1189mm，U max ＝249mm，I max ＝4.2mm/m，ε max ＝1.3mm/m，K max ＝0.02mm/m 2 。
Claims (9)
1. A surface subsidence prediction method for deformation area control of a coal-grain composite area is characterized by comprising the following steps:
firstly, determining underground coal seam parameters of a coal-grain composite area through geological mining data;
then, according to the nonlinear relation between the mining degree and the sinking coefficient, a mining method of mining while filling is used, the mining of the caving face is carried out after the filling of the filling face at the two sides is completed,determining the mining width b of the caving face, and then according to the width a and the filling rate of the filling face
The number of the collapsed working faces is n, and the number of the filled working faces is n-1;
then calculating the equivalent mining height m of the filling working face according to the equivalent mining height theory 1 ；
Dividing the mining range R of the coal seam into an upper part and a lower part by the profile of the filling working face, wherein the upper part of the profile is equivalent to the sinking space of the top plate of the whole filling mining working face and is recorded as the mining range R 1 The lower half part of the profile is a caving space below a sinking top plate, and the mining range is recorded as R 2 Wherein the production range is R 1 All are in a mining empty state, and the mining range is R 2 Then reserving unfilled space behind a filling body made of filling materials for the working face at intervals, wherein the width of the filling body is equal to the width a of the filling mining working face; wherein the lower half part of the profile is the sinking height of the caving space under the sinking top plate
Filling rate of filling body formed by filling working face
(ii) an effect;
respectively calculating mining ranges R by utilizing a probability integral model 1 And R 2 Induced surface subsidence W 1 And W 2 Horizontally moving U 1 And U 2 Inclination I 1 And I 2 Horizontal deformation epsilon 1 And epsilon 2 And curvature K 1 And K 2 And after the two groups of deformation values are correspondingly superposed, finally, estimating parameters for controlling the surface deformation area: surface subsidence W, horizontal movement U, inclination I, horizontal deformation epsilon, and curvature K.
2. The method for predicting surface subsidence oriented to deformation region control of coal and grain composite area according to claim 1, wherein the method is characterized in thatThe method comprises the following steps: section mining range R = R of coal bed 1 +nR 2 Wherein R is 1 ＝[nb+(n-1)a]m 1 ，
m 1 Is the roof sinking height (i.e. equivalent mining height),
the height of the filling body after being compacted by the top plate.
3. The method for predicting the surface subsidence for the deformation area control of the coal and grain composite area according to claim 1, which comprises the following specific steps:
step 1: aiming at a coal and grain composite region needing precise protection of ground cultivated land, collecting and researching geological mining data, a drilling histogram and other information of the region to be mined, and determining mining depth H, mining thickness m and coal seam inclination angle alpha information of a coal seam to be mined;
step 2: establishing a nonlinear relation between mining degree and subsidence coefficient according to surface subsidence monitoring data mined by an underground working face in the mining area or similar geological mining conditions, determining mining width b of a caving working face based on the nonlinear relation, and designing filling mining working face width a and filling rate according to the range of a region to be mined, ground cultivated land protection conditions, filling methods and filling material information
Designing caving working faces and filling mining working faces at intervals, wherein the number of the caving working faces is n, and the number of the filling mining working faces is n-1;
and 3, step 3: according to mining thickness m and filling rate
Calculating the equivalent mining height of the filling mining working face according to the equivalent mining height theory
Thus, the section mining range R for controlling the surface deformation region for the farmland protection of the coal and grain composite region is regarded as the length nb + (n-1) a and the height m 1 Mining area R 1 Plus n lengths b and heights
Mining area R 2 ；
And 4, step 4: and (3) actual measurement predicted parameters of a large amount of existing surface subsidences of the reference mining area: a sinking coefficient q, a horizontal movement coefficient B, a main influence angle tangent tan beta, a knee offset S and a mining influence propagation angle theta are respectively determined to be mining ranges R 1 And a mining range R 2 The probability integral prediction parameter of (2): coefficient of subsidence q 1 、q 2 Coefficient of horizontal movement B 1 、B 2 The main influence of the tangent tan beta 1 、tanβ 2 Inflection offset S 1 、S 2 Angle of propagation of mining effect theta 1 、θ 2 ；
And 5: calculating the mining range R by utilizing a probability integral model according to the subsidence prediction parameters 1 Induced subsidence of the earth's surface W 1 The earth's surface moves horizontally U 1 Surface inclination I 1 Horizontal deformation of earth's surface epsilon 1 And surface curvature K 1 The mining range R can also be calculated 2 Induced subsidence of the earth's surface W 2 Horizontal movement of the earth's surface U 2 Surface inclination I 2 Horizontal deformation of earth's surface epsilon 2 And surface curvature K 2 Surface subsidence controlled by surface deformation area W = W 1 +W 2 The earth surface horizontally moves U = U 1 +U 2 Surface inclination I = I 1 +I 2 Horizontal deformation of earth surface epsilon = epsilon 1 +ε 2 And surface curvature K = K 1 +K 2 And further, the surface subsidence prediction of deformation area control for the farmland protection of the coal-grain complex area is realized.
4. The method for predicting subsidence of ground surface of deformation area control for arable land protection of coal and grain composite area according to claim 2, wherein the method comprises the steps ofIn the mining area R 1 The sinking coefficient is: the actual measurement parameters under the mining condition of the mining area thin coal seam caving method are adopted according to the equivalent mining height theory, wherein the sinking coefficient is increased by 5-10 percent, and the method specifically comprises the following steps:
q 1 ＝q(1+x)
in the formula: q. q of 1 For the exploitation range R based on the equivalent mining height theory 1 A sinking coefficient; q is the mining subsidence coefficient of the caving method of the mining area; the value range of x is between 0.05 and 0.1;
mining area R 1 Horizontal movement coefficient B of 1 N, mainly influences the tangent tan β 1 = tan β, inflection offset S 1 ＝S。
5. The method for predicting subsidence of ground surface facing deformation area control of arable land protection in coal and grain composite area according to claim 2, wherein the mining range R 2 Through the sinking coefficient of q 2 The/q is expressed as:
in the formula: q. q of 2 To the mining area R 2 The sinking coefficient of (a); q is mining subsidence coefficient by mining area caving method; m is the mining area R 2 Thickness of mining
The unit m; b is the mining area R 2 The width of mining (1) is in unit m; a is the width of the filling working face, namely the width of the filling working face, and the unit is m; h is the depth of cut in m.
6. The method for predicting subsidence of ground surface facing deformation area control of arable land protection in coal and grain composite area according to claim 2, wherein the mining range R 2 Horizontal movement coefficient B of 2 The calculation method comprises the following steps:
in the formula: b is 2 To the mining area R 2 A horizontal movement coefficient of; b is mining horizontal movement coefficient of mining area caving method; b is the mining area R 2 The width of mining (1) is in unit m; a is the width of the filling working face, namely the width of the remaining width, and the unit m; h is the depth of cut in m.
7. The method for predicting subsidence of ground surface for deformation area control of arable land protection in coal and grain complex area as claimed in claim 2, wherein the extraction range R 2 Is mainly influencing the tangent tan beta 2 The calculation process of (2):
tanβ 2 /tanβ＝0.7847e -0.0012PH
in the formula: tan beta 2 To the mining area R 2 The main influence of (1) on the tangent; tan beta is the angle tangent mainly influenced by mining by a mining area caving method; p is the comprehensive evaluation coefficient of the overburden stratum, and each mine area has own parameter value; h is the depth of cut in m.
8. The method for predicting surface subsidence of deformation region control for arable land protection of coal-grain composite areas as claimed in claim 2, wherein the inflection offset S 2 The calculation process of (2):
in the formula: s 2 To the mining area R 2 (iv) inflection offset of; b is the mining area R 2 The width of mining of (1) in units of m; a is the width of the filling working face, namely the width of the filling working face, and the unit is m; h is the depth of cut in m.
9. The method for predicting subsidence of ground surface facing deformation area control of arable land protection of coal and grain composite area of claim 2, wherein the mining influence propagation angle θ is the mining influence propagation angle of the whole caving method mining under similar geological conditions.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211002816.9A CN115345372B (en) | 2022-08-19 | 2022-08-19 | Surface subsidence prediction method for deformation area control of coal-grain composite area |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211002816.9A CN115345372B (en) | 2022-08-19 | 2022-08-19 | Surface subsidence prediction method for deformation area control of coal-grain composite area |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115345372A true CN115345372A (en) | 2022-11-15 |
| CN115345372B CN115345372B (en) | 2024-02-09 |
Family
ID=83953618
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211002816.9A Active CN115345372B (en) | 2022-08-19 | 2022-08-19 | Surface subsidence prediction method for deformation area control of coal-grain composite area |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115345372B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115829133A (en) * | 2022-12-12 | 2023-03-21 | 通用技术集团工程设计有限公司 | Surface subsidence prediction method for strip coal pillar recovered by filling process |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103902780A (en) * | 2014-04-08 | 2014-07-02 | 中国矿业大学 | Method for predicting deformation of solid-filled coal mining surfaces |
| CN109555556A (en) * | 2019-01-21 | 2019-04-02 | 中国矿业大学 | A kind of solid filling coal mining Full Ratio design method towards cultivated land protection |
| CN109992837A (en) * | 2019-03-06 | 2019-07-09 | 太原理工大学 | A method of using skewness influence function Mountainous Area coal-mining subsidence |
| CN111255453A (en) * | 2020-03-10 | 2020-06-09 | 安徽理工大学 | Method for reducing ground surface subsidence range of underground mining area |
| CN211692557U (en) * | 2019-12-29 | 2020-10-16 | 寇雅威 | A fill device for collecting space area |
| CN112184902A (en) * | 2020-09-21 | 2021-01-05 | 东华理工大学 | Underground mining face inversion method for boundary crossing mining identification |
-
2022
- 2022-08-19 CN CN202211002816.9A patent/CN115345372B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103902780A (en) * | 2014-04-08 | 2014-07-02 | 中国矿业大学 | Method for predicting deformation of solid-filled coal mining surfaces |
| CN109555556A (en) * | 2019-01-21 | 2019-04-02 | 中国矿业大学 | A kind of solid filling coal mining Full Ratio design method towards cultivated land protection |
| CN109992837A (en) * | 2019-03-06 | 2019-07-09 | 太原理工大学 | A method of using skewness influence function Mountainous Area coal-mining subsidence |
| CN211692557U (en) * | 2019-12-29 | 2020-10-16 | 寇雅威 | A fill device for collecting space area |
| CN111255453A (en) * | 2020-03-10 | 2020-06-09 | 安徽理工大学 | Method for reducing ground surface subsidence range of underground mining area |
| CN112184902A (en) * | 2020-09-21 | 2021-01-05 | 东华理工大学 | Underground mining face inversion method for boundary crossing mining identification |
Non-Patent Citations (6)
| Title |
|---|
| GUOJIAN ZHANG, GUANGLI GUO, WEI WEI, SIFENG ZHANG, GUOBIAO YAO: "A NEW METHOD TO CONTROL THE REGIONAL STRATA MOVEMENT OF SUPER-THICK WEAK CEMENTATION OVERBURDEN IN DEEP MINING", 《THE CIVIL ENGINEERING JOURNAL》 * |
| 朱广轶, 王玮, 于莉: "煤岩开采与地表沉陷的预计程序", 沈阳大学学报, no. 02 * |
| 李艳;邱志伟;马文明;: "煤矿区土地损毁程度预测分析――以李粮店矿区破坏程度为例", 山东国土资源, no. 06 * |
| 王启春;李建;郭广礼;: "重庆市矿区开采沉陷预计参数类比分析研究", 煤炭科学技术, no. 10 * |
| 郝长胜;马占一;孙宝雷;李日东;: "近水平煤层开采地表沉陷预测系统研究", 矿业安全与环保, no. 01 * |
| 韩永斌: "概率积分法预计参数求取及应用研究", 《矿山测量》, vol. 49, no. 4 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115829133A (en) * | 2022-12-12 | 2023-03-21 | 通用技术集团工程设计有限公司 | Surface subsidence prediction method for strip coal pillar recovered by filling process |
| CN115829133B (en) * | 2022-12-12 | 2023-07-18 | 通用技术集团工程设计有限公司 | Surface subsidence prediction method for filling process recovery strip coal pillar |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115345372B (en) | 2024-02-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107526873B (en) | Shallow tunnel surrounding rock collapse mode identification and supporting structure calculation method | |
| CN111750822B (en) | Coal mining induced overlying strata and surface subsidence collaborative dynamic prediction method | |
| Jachens et al. | Differential compaction mechanism for earth fissures near Casa Grande, Arizona | |
| CN109577982A (en) | Wall type Lian Cailian fills water protection mining and water resource migration monitoring, water damage method for early warning | |
| CN106326528B (en) | A kind of opencut end side underground mining induces Distribution Law of Ground Crack prediction technique | |
| CN112446143B (en) | Strip filling mining parameter design method under thick unconsolidated formation and thin bedrock | |
| CN107489424A (en) | A kind of shield subway work induces stratum deformation and the predictor method influenceed on ancient building | |
| CN107784150B (en) | Crack sealing performance prediction method | |
| CN109798149A (en) | A kind of tight roof high seam stope ground pressure strength grade predicting method | |
| CN110765642A (en) | Zonal evaluation method for roof rock stratum structure of coal seam, mining area or working face | |
| CN110727988A (en) | Deep soft foundation consolidation settlement layering summation algorithm based on soft soil parameter space anisotropy | |
| Ren et al. | The overburden strata caving characteristics and height determination of water conducting fracture zone in fully mechanized caving mining of extra thick coal seam | |
| CN110555559A (en) | Ground settlement prediction method under insufficient mining condition of thick unconsolidated formation | |
| CN115345372A (en) | Surface subsidence prediction method for deformation area control of coal and grain composite area | |
| CN113420457B (en) | End slope stability analysis method for steep slope mining of open pit coal mine | |
| CN113339073B (en) | Impact risk evaluation method based on roof rock stratum structure | |
| CN115344928B (en) | Surface subsidence superposition calculation method for deformation control of cultivated underground coal mining area | |
| CN113010993B (en) | Method for predicting height of water-guiding crack zone of thick-punching laminated mining area | |
| CN111737895B (en) | Method for dynamically evaluating stability of top plate of underground goaf of strip mine | |
| CN112948924A (en) | Near unconsolidated formation mining water flowing fractured zone height determination method based on overlying strata structure | |
| CN108549780A (en) | Shallow embedding deep space stays the load evaluation method of tunnel filling body | |
| CN110851991A (en) | Underground water flow numerical simulation method | |
| CN110630316A (en) | Method for judging development characteristics of waste working face gas guide crack zone | |
| CN111159643B (en) | Method for influencing ground subsidence by Weibull time function based on measured data | |
| CN114817826A (en) | Equivalent subsidence-based mining overburden bed separation calculation method |
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 |