CN115345372B - Surface subsidence prediction method for deformation area control of coal-grain composite area - Google Patents

Surface subsidence prediction method for deformation area control of coal-grain composite area Download PDF

Info

Publication number
CN115345372B
CN115345372B CN202211002816.9A CN202211002816A CN115345372B CN 115345372 B CN115345372 B CN 115345372B CN 202211002816 A CN202211002816 A CN 202211002816A CN 115345372 B CN115345372 B CN 115345372B
Authority
CN
China
Prior art keywords
mining
range
subsidence
filling
coal
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
Application number
CN202211002816.9A
Other languages
Chinese (zh)
Other versions
CN115345372A (en
Inventor
郭广礼
李怀展
查剑锋
宫亚强
蒋乾
刘文锴
胡青峰
殷和建
汪铁宁
张学伟
王跃宗
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Mining and Technology CUMT
Original Assignee
China University of Mining and Technology CUMT
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China University of Mining and Technology CUMT filed Critical China University of Mining and Technology CUMT
Priority to CN202211002816.9A priority Critical patent/CN115345372B/en
Publication of CN115345372A publication Critical patent/CN115345372A/en
Application granted granted Critical
Publication of CN115345372B publication Critical patent/CN115345372B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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

Abstract

The invention discloses a ground subsidence prediction method for controlling a deformation area of a coal-grain composite area, which is suitable for ground subsidence prediction. Firstly, determining the address information of a coal bed, then determining the information of a collapse working surface, the information of a filling mining working surface and the information of backfill filling according to the nonlinear relation between mining degree and sinking coefficient, calculating the equivalent mining height of the filling working surface based on the equivalent mining height theory, and calculating the earth surface subsidence, earth surface horizontal movement, earth surface inclination, earth surface horizontal deformation and earth surface curvature caused by the mining range by utilizing a probability integral model, thereby realizing the earth surface subsidence prediction of the control of the deformation area facing the tillage protection of the coal-grain composite area. The method has the advantages of simple calculation steps and accurate calculation results, is more in line with the actual situation, and has important practical value for underground coal mining and ground tillage precision protection.

Description

Surface subsidence prediction method for deformation area control of coal-grain composite area
Technical Field
The invention relates to a surface subsidence prediction method for controlling a deformation area, which is applicable to a surface subsidence prediction method for controlling the deformation area of a coal-grain composite area, and belongs to the field of surface subsidence prediction.
Technical Field
However, due to the fact that underground coal mining damages earth surface building, cultivated land and the like, the conventional filling mining, strip mining and other methods commonly used for reducing ground subsidence have the problems of high cost, large resource loss and the like, and considering the beneficial factors of high deformation tolerance, large allowable deformation value and the like of the cultivated land, a new method for controlling coal mining in a deformation area facing to the cultivation protection of a coal-grain composite area is gradually developed, and prediction of earth surface subsidence is a problem to be solved.
Prediction of surface subsidence is related to a number of factors including the production thickness, the dip angle of the coal seam, the production depth, the prediction method, and the prediction parameters, among others, where a scientific and reasonable prediction method relates to whether the final predicted result is accurate. Although many scholars develop researches on the ground subsidence prediction of the conventional caving method exploitation, filling exploitation and strip exploitation at present, no ground subsidence prediction method facing to the deformation area control of the cultivated land protection of the coal-grain composite area exists at present.
Disclosure of Invention
Aiming at the technical problems, the ground surface subsidence prediction method for controlling the deformation area of the coal-grain composite area, which is simple in calculation steps, high in accuracy and easy to obtain calculation parameters, is provided.
In order to achieve the technical aim, the invention provides a ground surface subsidence prediction method for controlling a deformation area of a coal-grain composite area, which comprises the following steps:
firstly, determining parameters of an underground coal bed in a coal-grain composite area through geological mining data;
then according to the nonlinear relation of the mining range degree and the sinking coefficient, a mining method of filling while mining is used, mining of the caving working face is carried out after filling is completed on the filling working faces at two sides, the mining width b of the caving working face is determined, and then according to the width a and the filling rate of the filling mining working faceThe number of the caving working surfaces is n, and the number of the filling working surfaces is n-1;
then calculating the equivalent mining height m of the filling working surface according to the equivalent mining height theory 1 ；
Dividing the coal seam mining range R into an upper part and a lower part by the section of the filling working face, wherein the upper part of the section is equivalent to the sinking space of the top plate of the whole filling working face, and the mining range is R 1 The lower half of the section is the caving space under the sinking roof, and the mining range is R 2 Wherein the exploitation range is R 1 All are in the empty state, and the mining range is R 2 The unfilled space behind the filling body formed by the filling material is reserved for the working face interval, and the width of the filling body is equal to the width a of the filling mining working face; wherein the lower half of the section is the sinking height of the sinking space of the sinking top plateFilling rate of filling body formed by filled working face +.>Influence;
calculating the exploitation range R by using a probability integration model 1 And R is 2 Induced ground surfaceSink W 1 And W is 2 Horizontal movement 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 correspondingly superposing the two groups of deformation values and finally budgeting parameters for controlling the surface deformation area: ground subsidence W, horizontal movement U, inclination I, horizontal deformation epsilon, and curvature K.
Further, the section exploitation range r=r of the coal seam 1 +nR 2 Wherein R is 1 ＝[nb+(n-1)a]m 1 ，m 1 For the roof sinking height/equivalent height, +.>Is 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-grain composite area which needs to be accurately protected in ground cultivated land, collecting and researching information such as geological mining data, a drilling histogram and the like of the area to be mined, and determining mining depth H, mining thickness m and coal seam inclination angle alpha information of the coal seam to be mined;
step 2: according to the ground subsidence monitoring data of underground working surface mining in the mining area or similar geological mining conditions, a nonlinear relation between the mining degree and the subsidence coefficient is established, the mining width b of the caving working surface is determined based on the nonlinear relation, and then the filling mining working surface width a and the filling rate are designed according to the range of the area to be mined, the ground plowing protection condition, the filling method and the filling material informationThe method comprises the steps of designing a caving working face and a filling and mining working face at intervals, wherein the number of the caving working face is n, and the number of the filling and mining working face is n-1;
step 3: according to the mining thickness m and the filling rateCalculating the equivalent mining height of the filling mining working face according to the equivalent mining height theory>The section exploitation range R of the surface deformation area control facing the coal grain composite area cultivated land protection is regarded as the length nb+ (n-1) a and the height m 1 Is set in the range R 1 Plus n lengths of b and heights of +.>Is set in the range R 2 ；
Step 4: reference to existing mass surface subsidence measured predicted parameters for a mine: the sinking coefficient q, the horizontal movement coefficient B, the main influence tangent tan beta, the inflection point offset distance S and the exploitation influence propagation angle theta are respectively determined to determine the exploitation range R 1 And a mining range R 2 Probability integral prediction parameters of (a): dip coefficient q 1 、q 2 Coefficient of horizontal movement B 1 、B 2 Mainly affects tan beta 1 、tanβ 2 Inflection point offset S 1 、S 2 Mining affects propagation angle θ 1 、θ 2 ；
Step 5: calculating the mining range R by using a probability integration model according to the subsidence prediction parameters 1 Induced subsurface subsidence W 1 Horizontal movement U of ground surface 1 Inclination of ground surface I 1 Horizontal deformation epsilon of ground surface 1 And surface curvature K 1 The exploitation range R can also be calculated 2 Induced subsurface subsidence W 2 Horizontal movement U of ground surface 2 Inclination of ground surface I 2 Horizontal deformation epsilon of ground surface 2 And surface curvature K 2 Ground subsidence w=w controlled by ground deformation zone 1 +W 2 Ground level shift u=u 1 +U 2 Surface inclination i=i 1 +I 2 Horizontal deformation of the earth's surface epsilon=epsilon 1 +ε 2 And surface curvature k=k 1 +K 2 And further, the ground surface subsidence prediction for controlling the deformation area of the cultivated land protection of the coal-grain composite area is realized.
Further, the exploitation range R 1 The dip coefficient of (2) is: referring to the actual measurement parameters under the mining conditions of the mining area thin coal seam collapse method by adopting the equivalent mining height theory, wherein the sinking coefficient is increased by 5% -10%, and the method specifically comprises the following steps:
q 1 ＝q(1+x)
wherein: q 1 Is the exploitation range R based on the equivalent exploitation height theory 1 A dip coefficient; q is the subsidence coefficient of mining by the mining area collapse method; the value range of x is between 0.05 and 0.1;
production range R 1 Horizontal movement coefficient b of (2) 1 =b, mainly affecting the tangent tan β 1 =tan β, inflection point offset S 1 ＝S。
Further, the exploitation range R 2 Dip coefficient of (1) is equal to q 2 And/q is expressed as:
wherein: q 2 For the range R 2 Is a dip coefficient of (a); q is a mining subsidence coefficient of a mining area by a caving method; m is the mining range R 2 Mining thicknessA unit m; b is the mining range R 2 Is a unit m; a is the remaining width, namely the width of the filling working surface, and the unit is m; h is the sampling depth, and the unit is m.
Further, the exploitation range R 2 Horizontal movement coefficient B of (2) 2 The calculation method of (1) is as follows:
wherein: b (B) 2 For the range R 2 Is a horizontal movement coefficient of (2); b is a mining area collapse method mining horizontal movement coefficient; b is the mining range R 2 Is a unit m; a is the remaining width, namely the width of the filling working surface, and the unit is m; h is the sampling depth, and the unit is m.
Further, the exploitation range R 2 Is the main effect of (a) on tan beta 2 Is calculated according to the following steps:
tanβ 2 /tanβ＝0.7847e -0.0012PH
wherein: tan beta 2 For the range R 2 Is a major impact on the tangent; tan beta is the tangent of the mining area which is mainly affected by mining through a mining area caving method; p is the comprehensive evaluation coefficient of the overburden, and each mine area has own parameter value; h is the sampling depth, and the unit is m.
Further, inflection point offset S 2 Is calculated according to the following steps:
wherein: s is S 2 For the range R 2 Is a inflection point offset distance; b is the mining range R 2 Is a unit m; a is the remaining width, namely the width of the filling working surface, and the unit is m; h is the sampling depth, and the unit is m.
Further, the exploitation influence propagation angle θ takes the exploitation influence propagation angle when exploiting by all collapse methods under similar geological conditions.
The beneficial effects are that:
the method considers the compression characteristic of the filling body controlled by the deformation area of the coal-grain composite area, gives consideration to the earth surface subsidence control and the cultivated land protection, scientifically and reasonably provides an earth surface subsidence prediction method for the deformation area control of the cultivated land protection of the coal-grain composite area, fills the blank of the earth surface subsidence prediction method for controlling the green mining of coal in the regional subsidence of the coal-grain composite area, provides technical support for the collaborative development of the green mining of coal and the cultivated land protection of the coal-grain composite area, and has important practical significance and application value for realizing the regional earth surface subsidence control, the cultivated land protection and the like.
Based on the prior art, the method refers to the equivalent mining height m of the equivalent mining height theory through scientific and reasonable analysis 1 And a filler compaction thicknessInnovative ofDividing underground coal mining space suitable for controlling coal mining in farmland protection deformation area of coal-grain composite area into R 1 And R is 2 R is then determined based on the existing subsidence prediction parameters 1 、R 2 And (3) independently calculating by adopting a probability integral model, and finally superposing corresponding calculation results to obtain the ground surface movement and deformation caused by the whole mining space, thereby being particularly suitable for ground surface subsidence calculation when the deformation of the coal-grain composite region is controlled to be used for farmland protection.
Drawings
FIG. 1 is a schematic calculation diagram of a surface subsidence prediction method facing to control of deformation areas of a coal-grain composite area.
Detailed Description
The invention will be described in further detail with reference to the drawings and the detailed description below:
as shown in FIG. 1, the method for predicting the subsidence of the earth surface facing the control of the deformation area of the coal-grain composite area comprises the following steps:
firstly, determining coal bed parameters through geological mining data; then determining the mining width b of the caving working face according to the nonlinear relation between the mining range and the sinking coefficient; then designing the width a and the filling rate of the filling mining working faceAnd the number n of the caving working surfaces, then the number n-1 of the filling working surfaces can be determined; then calculating the equivalent mining height m of the filling working surface according to the equivalent mining height theory 1 The method comprises the steps of carrying out a first treatment on the surface of the Finally, the profile mining range R can be regarded as having a length nb+ (n-1) a and a height m 1 Is set in the range R 1 Plus n lengths of b and heights of +.>Is set in the range R 2 . And then, the mining range and the caused surface subsidence, horizontal movement, inclination, horizontal deformation and curvature can be calculated by utilizing a probability integral model, and the surface subsidence W, the horizontal movement U, the inclination I, the horizontal deformation epsilon and the curvature K controlled by the surface deformation area can be obtained through superposition calculation.
The method comprises the following specific steps:
step 1: aiming at a coal-grain composite area which needs to be accurately protected in ground cultivated land, collecting and researching information such as geological mining data, a drilling histogram and the like of the area to be mined, and determining information such as mining depth H, mining thickness m, coal seam inclination angle alpha and the like of the coal seam to be mined.
Step 2: according to the ground subsidence monitoring data of underground working surface mining in the mining area or similar geological mining conditions, a nonlinear relation between the mining degree and the subsidence coefficient is established, the mining width b of the caving working surface is determined based on the nonlinear relation, and then the filling mining working surface width a and the filling rate are designed according to the information such as the area range to be mined, the ground tillage protection condition, the filling method, the filling material and the likeAnd simultaneously determining the number n of the caving working faces, and then determining the number n-1 of the filling and mining working faces.
Step 3: according to the mining thickness m determined in the step 1 and the filling rate determined in the step 2According to the equivalent mining height theory, the equivalent mining height of the filling mining working face can be calculated>The section exploitation range R of the surface deformation area control facing the coal grain composite area cultivated land protection is regarded as the length nb+ (n-1) a and the height m 1 Is set in the range R 1 Plus n lengths of b and heights of +.>Is set in the range R 2 。
Step 4: determining the exploitation range R by referring to the actual measured predicted parameters of a large number of surface subsidence existing in the mining area 1 And a mining range R 2 Probability integral prediction parameters of (a): the dip coefficient q, the horizontal movement coefficient B, the principal influencing tangent tan beta, the inflection point offset S, the mining influencing propagation angle theta,the method comprises the following steps:
(1) Production range R 1 Is of length nb+ (n-1) a and heightThe predicted parameters of the method can refer to the actual measurement parameters under the mining conditions of the mining area thin coal seam caving method by the equivalent mining height theory, wherein the sinking coefficient is increased by 5% -10%, and the method specifically comprises the following steps:
q 1 ＝q(1+x)
wherein: q 1 Is the exploitation range R based on the equivalent exploitation height theory 1 A dip coefficient; q is the subsidence coefficient of mining by the mining area collapse method; the value range of x is between 0.05 and 0.1.
(2) Production range R 2 From n lengths b and heightsIs used for the production of the working face, the production range R is obtained because the filling body of the filling face is compacted 2 The predicted parameters of (a) can be determined by referring to the mining width b and the reserved width a, and the specific calculation formula is as follows:
1) Calculating a sinking coefficient:
wherein: q 2 For the range R 2 Is a dip coefficient of (a); q is a mining subsidence coefficient of a mining area by a caving method; m is the mining range R 2 Mining thicknessA unit m; b is the mining range R 2 Is a unit m; a is the remaining width, namely the width of the filling working surface, and the unit is m; h is the sampling depth, and the unit is m.
2) Calculating a horizontal movement coefficient:
wherein: b (B) 2 For the range R 2 Is a horizontal movement coefficient of (2); b is a mining area collapse method mining horizontal movement coefficient; b is the mining range R 2 Is a unit m; a is the remaining width, namely the width of the filling working surface, and the unit is m; h is the sampling depth, and the unit is m.
3) Calculation of dominant impact tangent
tanβ 2 /tanβ＝0.7847e -0.0012PH
Wherein: tan beta 2 For the range R 2 Is a major impact on the tangent; tan beta is the tangent of the mining area which is mainly affected by mining through a mining area caving method; p is the comprehensive evaluation coefficient of the overburden, and each mine area has own parameter value; h is the sampling depth, and the unit is m.
4) Calculation of inflection point offset
Wherein: s is S 2 For the range R 2 Is a inflection point offset distance; b is the mining range R 2 Is a unit m; a is the remaining width, namely the width of the filling working surface, and the unit is m; h is the sampling depth, and the unit is m.
5) Calculation of mining influence propagation angle
The exploitation influence propagation angle can be similar to the exploitation influence propagation angle when exploiting all the collapse methods under geological conditions.
Step 5: calculating the mining range S by using a probability integration model according to the subsidence prediction parameters calculated in the step 4 1 Induced subsurface subsidence W 1 Horizontal movement U of ground surface 1 Inclination of ground surface I 1 Horizontal deformation epsilon of ground surface 1 And surface curvature K 1 The exploitation range S can also be calculated 2 Induced subsurface subsidence W 2 Horizontal movement U of ground surface 2 Inclination of ground surface I 2 Horizontal deformation epsilon of ground surface 2 And surface curvature K 2 Thus, the surface subsidence W=W of the surface deformation zone control can be calculated 1 +W 2 Ground horizontal movement U =U 1 +U 2 Surface inclination i=i 1 +I 2 Horizontal deformation of the earth's surface epsilon=epsilon 1 +ε 2 And surface curvature k=k 1 +K 2 And further, the ground surface subsidence prediction for controlling the deformation area of the cultivated land protection of the coal-grain composite area is realized.
Embodiment 1,
Step 1: and determining the coal bed parameters of the mining area according to the geological mining data, the drilling histogram and other information of the mining area, and obtaining the coal bed thickness m=3m, the coal bed inclination angle alpha=0° and the mining depth H=500m.
Step 2: determining the width b=180m of the caving working 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 rateThe number of caving faces n=3, and the number of filling faces n-1=2.
Step 3: calculating the equivalent mining height of the filling working surface based on the equivalent mining height theory The profile extraction range R can be considered as having a length nb+ (n-1) a=840 m and a height m 1 Mining range r=0.6m 1 Plus n=3 lengths b=180m and heights +.>Is set in the range R 2 。
Step 4: obtaining the estimated subsidence parameters of the mining area by referring to the estimated parameters of the mining area in actual measurement: dip coefficient q=0.5, horizontal movement coefficient b=0.27; mainly affecting the tangent: tan β=2.3; inflection point offset s=0.05h=25m, mining affects propagation angle θ=90°; respectively calculating the exploitation range R according to a calculation formula 1 Is expected to be parameter for dip: q 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 Is expected to be parameter for dip: q 2 ＝0.366，B 2 ＝0.188，tanβ 2 ＝1.6，S 2 ＝24.6m，θ 2 ＝90°。
Step 5: substituting the subsidence predicted parameter calculated in the step 4 into a probability integration model to calculate the mining range S 1 Induced subsurface subsidence W 1 Horizontal movement U of ground surface 1 Inclination of ground surface I 1 Horizontal deformation epsilon of ground surface 1 And surface curvature K 1 The exploitation range S can also be calculated 2 Induced subsurface subsidence W 2 Horizontal movement U of ground surface 2 Inclination of ground surface I 2 Horizontal deformation epsilon of ground surface 2 And surface curvature K 2 Thus, the surface subsidence W=W of the surface deformation zone control can be calculated 1 +W 2 Ground level shift u=u 1 +U 2 Surface inclination i=i 1 +I 2 Horizontal deformation of the earth's surface epsilon=epsilon 1 +ε 2 And surface curvature k=k 1 +K 2 And further, the ground surface subsidence prediction for controlling the deformation area of the cultivated land protection of the coal-grain composite area is realized. In the embodiment, the ground surface movement and deformation extremum controlled by the ground surface deformation area calculated by adopting the probability integral model is as follows: w (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 ground surface subsidence prediction method for controlling deformation areas of coal-grain composite areas is characterized by comprising the following steps:
firstly, determining parameters of a coal bed under a coal-grain composite area through geological mining data, wherein the parameters comprise mining depth H, mining thickness m and coal bed inclination angle alpha information;
then according to the nonlinear relation between the mining range and the sinking coefficient, a mining method of filling while mining is used, mining of the caving working face is carried out after filling is completed on the filling working faces at two sides, and the caving working face is determinedFace mining width b, then face width a and fill rateDetermining the number n of the caving working faces, wherein the number n of the filling working faces is n-1;
then calculating the equivalent mining height m of the filling working surface according to the equivalent mining height theory 1 ；
Dividing the coal seam mining range R into an upper part and a lower part by the section of the filling working face, wherein the upper part of the section is equivalent to the sinking space of the top plate of the whole filling working face, and the mining range is R 1 The lower half of the section is the caving space under the sinking roof, and is marked as the mining range R 2 Wherein the exploitation range is R 1 All are in the empty state, and the mining range is R 2 The unfilled space behind the filling body formed by the filling material is reserved for the working face interval, and the width of the filling body is equal to the width a of the filling mining working face; wherein the lower half of the section is the sinking height of the sinking space of the sinking top plateFilling rate of filling body formed by filled working face +.>Influence;
reference to existing mass surface subsidence measured predicted parameters for a mine: the sinking coefficient q, the horizontal movement coefficient B, the main influence tangent tan beta, the inflection point offset distance S and the exploitation influence propagation angle theta are calculated to obtain the exploitation range R 1 Dip predicted parameter dip coefficient q 1 Coefficient of horizontal movement B 1 Mainly affects tan beta 1 Inflection point offset S 1 Mining affects propagation angle θ 1， Production range R 2 Is expected to be parameter for dip: dip coefficient q 2 Coefficient of horizontal movement B 2 Mainly affects tan beta 2 Inflection point offset S 2 Mining affects propagation angle θ 2 ；
Substituting the calculated estimated parameter of subsidence into probability integralModel, calculate the exploitation range R 1 Induced subsurface subsidence W 1 Horizontal movement U of ground surface 1 Inclination of ground surface I 1 Horizontal deformation epsilon of ground surface 1 And surface curvature K 1 Range of exploitation R 2 Induced subsurface subsidence W 2 Horizontal movement U of ground surface 2 Inclination of ground surface I 2 Horizontal deformation epsilon of ground surface 2 And surface curvature K 2 Will extract the range R 1 Induced subsurface subsidence W 1 Horizontal movement U of ground surface 1 Inclination of ground surface I 1 Horizontal deformation epsilon of ground surface 1 And surface curvature K 1 And the exploitation range R 2 Induced subsurface subsidence W 2 Horizontal movement U of ground surface 2 Inclination of ground surface I 2 Horizontal deformation epsilon of ground surface 2 And surface curvature K 2 And finally budgeting out parameters for controlling the surface deformation area after superposition: ground subsidence W, horizontal movement U, inclination I, horizontal deformation epsilon, and curvature K.
2. The method for predicting the subsidence of the earth for controlling the deformation area of the composite region of coal and grain according to claim 1, wherein the method comprises the following steps: profile mining range r=r of coal seam 1 +nR 2 Wherein R is 1 ＝[nb+(n-1)a]m 1 ，m 1 For the sinking height of the top plate, namely equivalent height, the method comprises the steps of ++>Is the height of the filling body after being compacted by the top plate.
3. The method for predicting the subsidence of the ground surface facing the control of the deformation area of the coal-grain composite area according to claim 1, which is characterized by comprising the following specific steps:
step 1: aiming at a coal-grain composite area which needs to be accurately protected in ground cultivated land, collecting and researching geological mining data and drilling histogram information of the area to be mined, and determining mining depth H, mining thickness m and coal seam inclination angle alpha information of the coal seam to be mined;
step 2: according to the ground subsidence monitoring data of underground working surface mining in the mining area or similar geological mining conditions, a nonlinear relation between the mining degree and the subsidence coefficient is established, the mining width b of the caving working surface is determined based on the nonlinear relation, and then the filling mining working surface width a and the filling rate are designed according to the range of the area to be mined, the ground plowing protection condition, the filling method and the filling material informationThe method comprises the steps of designing a caving working face and a filling and mining working face at intervals, wherein the number of the caving working face is n, and the number of the filling and mining working face is n-1;
step 3: according to the mining thickness m and the filling rateCalculating the equivalent mining height of the filling mining working face according to the equivalent mining height theory>The section exploitation range R of the surface deformation area control facing the coal grain composite area cultivated land protection is regarded as the length nb+ (n-1) a and the height m 1 Is set in the range R 1 Plus n lengths of b and heights of +.>Is set in the range R 2 ；
Step 4: reference to existing mass surface subsidence measured predicted parameters for a mine: the sinking coefficient q, the horizontal movement coefficient B, the main influence tangent tan beta, the inflection point offset distance S and the exploitation influence propagation angle theta are respectively determined to determine the exploitation range R 1 And a mining range R 2 Is expected to be parameter for dip: dip coefficient q 1 、q 2 Coefficient of horizontal movement B 1 、B 2 Mainly affects tan beta 1 、tanβ 2 Inflection point offset S 1 、S 2 Mining affects propagation angle θ 1 、θ 2 ；
Step 5: based on the expected parameter of the subsidence,calculating the mining range R by using a probability integration model 1 Induced subsurface subsidence W 1 Horizontal movement U of ground surface 1 Inclination of ground surface I 1 Horizontal deformation epsilon of ground surface 1 And surface curvature K 1 The exploitation range R can also be calculated 2 Induced subsurface subsidence W 2 Horizontal movement U of ground surface 2 Inclination of ground surface I 2 Horizontal deformation epsilon of ground surface 2 And surface curvature K 2 Ground subsidence w=w controlled by ground deformation zone 1 +W 2 Ground level shift u=u 1 +U 2 Surface inclination i=i 1 +I 2 Horizontal deformation of the earth's surface epsilon=epsilon 1 +ε 2 And surface curvature k=k 1 +K 2 And further, the ground surface subsidence prediction for controlling the deformation area of the cultivated land protection of the coal-grain composite area is realized.
4. The method for predicting subsidence of earth for controlling deformation zone of composite region of coal and grain according to claim 2, wherein the exploitation range R 1 The dip coefficient of (2) is: referring to the actual measurement parameters under the mining conditions of the mining area thin coal seam collapse method by adopting the equivalent mining height theory, wherein the sinking coefficient is increased by 5% -10%, and the method specifically comprises the following steps:
q 1 ＝q(1+x)
wherein: q 1 Is the exploitation range R based on the equivalent exploitation height theory 1 A dip coefficient; q is the subsidence coefficient of mining by the mining area collapse method; the value range of x is between 0.05 and 0.1;
production range R 1 Horizontal movement coefficient B of (2) 1 =b, mainly affecting the tangent tan β 1 =tan β, inflection point offset S 1 ＝S。
5. The method for predicting subsidence of earth for controlling deformation zone of composite region of coal and grain according to claim 2, wherein the exploitation range R 2 Dip coefficient of (1) is equal to q 2 And/q is expressed as:
wherein: q 2 For the range R 2 Is a dip coefficient of (a); q is a mining subsidence coefficient of a mining area by a caving method; m is the mining range R 2 Is used for the production of the thickness of the steel,a unit m; b is the mining range R 2 Is a unit m; a is the remaining width, namely the width of the filling working surface, and the unit is m; h is the sampling depth, and the unit is m.
6. The method for predicting subsidence of earth for controlling deformation zone of composite region of coal and grain according to claim 2, wherein the exploitation range R 2 Horizontal movement coefficient B of (2) 2 The calculation method of (1) is as follows:
wherein: b (B) 2 For the range R 2 Is a horizontal movement coefficient of (2); b is a mining area collapse method mining horizontal movement coefficient; b is the mining range R 2 Is a unit m; a is the remaining width, namely the width of the filling working surface, and the unit is m; h is the sampling depth, and the unit is m.
7. The method for predicting subsidence of earth for controlling deformation zone of composite region of coal and grain according to claim 2, wherein the exploitation range R 2 Is the main effect of (a) on tan beta 2 Is calculated according to the following steps:
tanβ 2 /tanβ＝0.7847e -0.0012PH
wherein: tan beta 2 For the range R 2 Is a major impact on the tangent; tan beta is the tangent of the mining area which is mainly affected by mining through a mining area caving method; p is the comprehensive evaluation coefficient of the overburden, and each mine area has own parameter value; h is the sampling depth, and the unit is m.
8. The earth's surface subsidence facing the deformation zone control of the coal-grain composite zone according to claim 2The notch prediction method is characterized in that the inflection point offset distance S 2 Is calculated according to the following steps:
wherein: s is S 2 For the range R 2 Is a inflection point offset distance; b is the mining range R 2 Is a unit m; a is the remaining width, namely the width of the filling working surface, and the unit is m; h is the sampling depth, and the unit is m.
9. The method for predicting the subsidence of the earth surface facing the deformation area control of the composite region of the coal and the grain according to claim 2, wherein the exploitation influence propagation angle theta is taken as the exploitation influence propagation angle when the whole caving method is exploited under the similar geological condition.
CN202211002816.9A 2022-08-19 2022-08-19 Surface subsidence prediction method for deformation area control of coal-grain composite area Active CN115345372B (en)

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 CN115345372A (en) 2022-11-15
CN115345372B true 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)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115829133B (en) * 2022-12-12 2023-07-18 通用技术集团工程设计有限公司 Surface subsidence prediction method for filling process recovery strip coal pillar

Citations (6)

* Cited by examiner, † Cited by third party
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

Patent Citations (6)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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》.2022,全文. *
概率积分法预计参数求取及应用研究;韩永斌;《矿山测量》;第49卷(第4期);全文 *
煤岩开采与地表沉陷的预计程序;朱广轶, 王玮, 于莉;沈阳大学学报(第02期);全文 *
煤矿区土地损毁程度预测分析――以李粮店矿区破坏程度为例;李艳;邱志伟;马文明;;山东国土资源(第06期);全文 *
近水平煤层开采地表沉陷预测系统研究;郝长胜;马占一;孙宝雷;李日东;;矿业安全与环保(第01期);全文 *
重庆市矿区开采沉陷预计参数类比分析研究;王启春;李建;郭广礼;;煤炭科学技术(第10期);全文 *

Also Published As

Publication number Publication date
CN115345372A (en) 2022-11-15

Similar Documents

Publication Publication Date Title
Jachens et al. Differential compaction mechanism for earth fissures near Casa Grande, Arizona
CN107526873B (en) Shallow tunnel surrounding rock collapse mode identification and supporting structure calculation method
CN107489424B (en) A kind of predictor method that shield subway work induces stratum deformation and influences on ancient building
CN111750822B (en) Coal mining induced overlying strata and surface subsidence collaborative dynamic prediction method
CN104965994A (en) Determining and estimating method for surface subsidence characteristic parameters caused by subway tunnel construction
CN110043312B (en) Control method for grouting filling ground surface settlement range
CN115345372B (en) Surface subsidence prediction method for deformation area control of coal-grain composite area
CN109555556A (en) A kind of solid filling coal mining Full Ratio design method towards cultivated land protection
CN106555607A (en) A kind of old goaf grouting filling control earth's surface residual settlement method
Li et al. Experimental study on instability mechanism and critical intensity of rainfall of high-steep rock slopes under unsaturated conditions
Asadi et al. Profiling function for surface subsidence prediction in mining inclined coal seams
CN106405678A (en) Mining overlying strata water flowing fracture zone height detection method based on stress monitoring
CN115344928B (en) Surface subsidence superposition calculation method for deformation control of cultivated underground coal mining area
CN104123454A (en) Method for predicating ground surface settlement caused by underground excavation of single-joint hard rock
Saltus Upper-crustal structure beneath the Columbia River Basalt Group, Washington: Gravity interpretation controlled by borehole and seismic studies
CN112948924A (en) Near unconsolidated formation mining water flowing fractured zone height determination method based on overlying strata structure
CN111737895B (en) Method for dynamically evaluating stability of top plate of underground goaf of strip mine
CN108549780A (en) Shallow embedding deep space stays the load evaluation method of tunnel filling body
CN113375753A (en) Method for monitoring and analyzing influence of coal mine fully mechanized mining face mining on underground water
Zhang et al. Evaluation of surface subsidence due to inclined coal seam mining: A case study in the 1930 Coal Mine, China
CN104196569A (en) Prediction method for judging whether water-resisting layer on filling coal face is stable or not
CN109184705B (en) Method for controlling stratum deformation in shield tunnel complex tunneling environment
CN110851991A (en) Underground water flow numerical simulation method
CN111159643B (en) Method for influencing ground subsidence by Weibull time function based on measured data
CN115758067A (en) Pile foundation impact drilling geological layering 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