CN107526920B - Method for analyzing corresponding relation between new surface subsidence dynamic process and overlying strata abscission layer development - Google Patents

Method for analyzing corresponding relation between new surface subsidence dynamic process and overlying strata abscission layer development Download PDF

Info

Publication number
CN107526920B
CN107526920B CN201710590471.6A CN201710590471A CN107526920B CN 107526920 B CN107526920 B CN 107526920B CN 201710590471 A CN201710590471 A CN 201710590471A CN 107526920 B CN107526920 B CN 107526920B
Authority
CN
China
Prior art keywords
subsidence
overburden
development
value
layer
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
CN201710590471.6A
Other languages
Chinese (zh)
Other versions
CN107526920A (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.)
Qingdao University of Technology
Original Assignee
Qingdao University of Technology
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 Qingdao University of Technology filed Critical Qingdao University of Technology
Priority to CN201710590471.6A priority Critical patent/CN107526920B/en
Publication of CN107526920A publication Critical patent/CN107526920A/en
Application granted granted Critical
Publication of CN107526920B publication Critical patent/CN107526920B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16ZINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00Subject matter not provided for in other main groups of this subclass

Landscapes

  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

The invention discloses a method for analyzing the corresponding relation between a new surface subsidence dynamic process and overburden abscission layer development, which comprises the steps of measuring a sticking measuring point on a model overburden section to obtain the time process of abscission layer from cracking to expanding to closing and quantitative development indexes of the length and the width of the abscission layer; arranging a subsidence measuring point on the model ground surface, and measuring the whole process data of the ground surface subsidence gradually developing from zero to the maximum value at the same time with the separation layer; and obtaining a time schedule by a model test, displaying the dynamic growth data of the surface subsidence and the dynamic development data of the overburden bed separation, statistically analyzing the corresponding relation of the dynamic growth data and the dynamic development data of the overburden bed separation, finding the corresponding relation of the initial period, the active period and the decline period of the surface subsidence and the initiation, development and closing time points of the bed separation, and newly proposing and dividing bed separation time nodes and each development stage. The characteristic of breaking the overburden separation layer can be more clearly reflected by indirectly obtaining data when the physical model and the digital model are used.

Description

Method for analyzing corresponding relation between new surface subsidence dynamic process and overlying strata abscission layer development
Technical Field
The invention relates to the technical field of surface subsidence dynamic processes and overburden abscission layers, in particular to a method for analyzing a corresponding relation between a new surface subsidence dynamic process and overburden abscission layer development.
Background
In the process of underground mining of the coal mine, the spontaneous combustion stress balance state of the overlying strata is destroyed along with the pushing of the underground mining working face, and the stress in the overlying strata of the goaf is redistributed under the action of the redistributed secondary stress field so as to achieve new balance. The overburden rock is broken, caving, bending and sinking, the movement change gradually develops from bottom to top to the ground surface to form ground surface subsidence, and the movement of the overburden rock is not completely synchronous and consistent due to different lithology, distance from a goaf position, thickness and the like, so that the displacement difference of overburden rock movement in the vertical direction occurs, and a separation space is formed.
The production conditions of the overburden separation layer are as follows: the delamination of rock stratum means that normal displacement is generated between two adjacent rock strata, the generation condition of the delamination mainly depends on the bending rigidity of the upper rock stratum and the lower rock stratum, when the bending rigidity of the upper rock stratum is larger than that of the lower rock stratum, the delamination is generated, and the bending rigidity of the rock stratum depends on the elastic modulus and the thickness of the rock stratum.
Disclosure of Invention
The invention aims to provide a method for analyzing the corresponding relation between a new surface subsidence dynamic process and overburden separation layer development.
In order to achieve the purpose, the invention provides the following technical scheme: measuring a sticking measuring point on a model overburden section to obtain a time process from cracking to expanding to closing of a separation layer and quantitative development indexes of the length and the width of the separation layer;
arranging a subsidence measuring point on the model ground surface, and measuring the whole process data of the ground surface subsidence gradually developing from zero to the maximum value at the same time with the separation layer;
and obtaining a time schedule by a model test, displaying the dynamic growth data of the surface subsidence and the dynamic development data of the overburden bed separation, statistically analyzing the corresponding relation of the dynamic growth data and the dynamic development data of the overburden bed separation, finding the corresponding relation of the initial period, the active period and the decline period of the surface subsidence and the initiation, development and closing time points of the bed separation, and newly proposing and dividing bed separation time nodes and each development stage.
The analysis method for the corresponding relation between the new surface subsidence dynamic process and the overburden abscission layer development is characterized in that the analysis is carried out by measuring a model overburden profile pasting measuring point, when the mining depth of a main overburden key layer on a working surface is 30-40m, the surface subsidence speed is the fastest, the subsidence speed reaches 7.5mm/d, therefore, the maximum subsidence speed of the surface is synchronous with the maximum subsidence speed of the key layer, the rapid subsidence of the surface is caused by the primary breakage of the key layer,
the calculation formula is as follows:
dW(t)/dt=c[W(t)/(1-1-e-ct)-Wk(1-1-e-ct)]
w (t) is a dynamic convergence value at a certain time
WkAs the final subsidence of the earth's surface
1-e-ctTo obtain the variation in a certain time
dW (t)/dt is the ground surface sinking speed, when W (t)/(1-1-e)-ct)-Wk(1-1-e-ct) When the value is greater than 0, the sinking speed of the time table is increased, and when W: (t)/(1-1-e-ct)-Wk(1-1-e-ct) The table descent speed decreases when the value is less than 0, but the resulting value deviates from the actual value to some extent due to the condition of the slant length and the length in the model.
Preferably, the analysis method for the corresponding relation between the new surface subsidence dynamic process and the overburden bed separation development simulates the movement of the key layer rock mass before breaking by setting the parameters of the joints among the key layer blocks to be large so that the blocks do not move, when the mining width reaches the breaking distance of the key layer, the parameters of the joints of the blocks are reduced, the movement of the key layer after breaking is simulated, the breaking of the overburden rock key layer influences the dynamic process of the surface subsidence,
the calculation formula is as follows:
W(t)=W0ω(1-ae-kt)
w (t) is a dynamic convergence value at a certain time
W0Is the maximum subsidence of the earth's surface at the end of the decline period
ΔωIs the residual subsidence value of the earth surface after the end of the decay period
Due to W0The amount of W (t) sinking is related to a, k and t when 1-ae is a fixed value-ktWhen the value is more than 0, the value of W (t) is increased, and when the value is 1-ae-ktWhen the value is less than 0, the W (t) sinking value is reduced, and the breaking of the key layer of the overburden rock causes the obvious increase of the surface sinking speed and does not change periodically.
Preferably, the analysis method for the corresponding relation between the new surface subsidence dynamic process and the overburden bed separation development is used for simulating excavation in a coal bed according to irregular depths, starting next excavation after model operation is stable after each excavation, when the depth is advanced to 150m, a primary breaking, a first period breaking and a second period breaking occur on a key layer respectively, the surface subsidence increment is smaller before the key layer is broken, the surface subsidence increment is increased after the key layer is broken, and a surface movement influence angle and a movement influence boundary are not constant,
the calculation formula is as follows:
W(x,y,t)=Wm(1-x2/a2-y2/b2)(1-e-ct)k
w (x, y, t) is the amount of sinking of any point of the sunken basin at any time
1-e-ctTo obtain the variation in a certain time
(1-x2/a2-y2/b2)(1-e-ct)kWhen the product is positive, the amount of W (x, y, t) sinking increases, (1-x)2/a2-y2/b2)(1-e-ct)kWhen the product is negative, the subsidence of W (x, y, t) is reduced, and the influence boundary of the surface subsidence is obviously enlarged along with the increase of a and b.
Compared with the prior art, the invention has the beneficial effects that: the method is characterized in that reasonable subsidence reduction mining technology and parameters are determined according to the characteristics of the overlying strata structure and the key layer under specific conditions.
Drawings
FIG. 1 is a schematic view of a broken-line structure of sinking speed and ground distance according to a new method for analyzing the corresponding relationship between the dynamic process of surface sinking and the development of overburden separation;
FIG. 2 is a table diagram of rock mechanical parameters of a new surface subsidence dynamic process and overburden separation development corresponding relation analysis method of the present invention;
fig. 3 is a table diagram illustrating the influence angle between mining width and surface movement in the method for analyzing the correspondence between the new surface subsidence dynamic process and overburden separation development according to the present invention.
Detailed Description
The invention provides a technical scheme that: measuring a sticking measuring point on a model overburden section to obtain a time process from cracking to expanding to closing of a separation layer and quantitative development indexes of the length and the width of the separation layer;
arranging a subsidence measuring point on the model ground surface, and measuring the whole process data of the ground surface subsidence gradually developing from zero to the maximum value at the same time with the separation layer;
and obtaining a time schedule by a model test, displaying the dynamic growth data of the surface subsidence and the dynamic development data of the overburden bed separation, statistically analyzing the corresponding relation of the dynamic growth data and the dynamic development data of the overburden bed separation, finding the corresponding relation of the initial period, the active period and the decline period of the surface subsidence and the initiation, development and closing time points of the bed separation, and newly proposing and dividing bed separation time nodes and each development stage.
A new surface subsidence dynamic process and overburden absciss layer development corresponding relation analysis method is characterized in that measurement of a model overburden profile pasting measuring point is used for analysis, when the mining depth of a main overburden key layer on a working surface is 30-40m, the surface subsidence speed is the fastest, the subsidence speed reaches 7.5mm/d, therefore, the maximum subsidence speed of the surface is synchronous with the maximum subsidence speed of the key layer, the initial breakage of the key layer leads to the rapid subsidence of the surface,
the calculation formula is as follows:
dW(t)/dt=c[W(t)/(1-1-e-ct)-Wk(1-1-e-ct)]
w (t) is a dynamic convergence value at a certain time
WkAs the final subsidence of the earth's surface
1-e-ctTo obtain the variation in a certain time
dW (t)/dt is the ground surface sinking speed, when W (t)/(1-1-e)-ct)-Wk(1-1-e-ct) When the value is larger than 0, the sinking speed of the time table is increased, and when W (t)/(1-1-e)-ct)-Wk(1-1-e-ct) The table descent speed decreases when the value is less than 0, but the resulting value deviates from the actual value to some extent due to the condition of the slant length and the length in the model.
A method for analyzing the relation between the dynamic sinking process of earth surface and the development of overlying strata separation layer includes such steps as setting the joint parameters between blocks in key layer to be very big for preventing the blocks from moving, simulating the movement before breaking the key layer, decreasing the joint parameters when the mining width reaches the breaking distance of key layer, simulating the movement after breaking the key layer, and breaking the key layer to influence the dynamic process of earth surface,
the calculation formula is as follows:
W(t)=W0ω(1-ae-kt)
w (t) is a dynamic convergence value at a certain time
W0Is the maximum subsidence of the earth's surface at the end of the decline period
ΔωIs the residual subsidence value of the earth surface after the end of the decay period
Due to W0The amount of W (t) sinking is related to a, k and t when 1-ae is a fixed value-ktWhen the value is more than 0, the value of W (t) is increased, and when the value is 1-ae-ktWhen the value is less than 0, the W (t) sinking value is reduced, and the breaking of the key layer of the overburden rock causes the obvious increase of the surface sinking speed and does not change periodically.
A new surface subsidence dynamic process and overlying strata abscission layer development corresponding relation analysis method is characterized in that excavation is carried out in a simulated coal bed according to irregular depths, model operation is carried out to be stable after each excavation, then next excavation is carried out, when the depth is advanced to 150m, a key layer is subjected to primary breakage, primary period breakage and secondary period breakage respectively, the surface subsidence increment is smaller before the key layer is broken, the surface subsidence increment is increased after the key layer is broken, and the surface movement influence angle and the movement influence boundary are not constant,
the calculation formula is as follows:
W(x,y,t)=Wm(1-x2/a2-y2/b2)(1-e-ct)k
w (x, y, t) is the amount of sinking of any point of the sunken basin at any time
1-e-ctTo obtain the variation in a certain time
(1-x2/a2-y2/b2)(1-e-ct)kWhen the product is positive, the amount of W (x, y, t) sinking increases, (1-x)2/a2-y2/b2)(1-e-ct)kWhen the product is negative, the amount of W (x, y, t) subsidence is reduced, and as a and b are increased, the influence boundary of the surface subsidence is increasedIs obviously enlarged.
On a similar simulation test model, obtaining the time process of the separation layer from cracking, expanding to closing and the quantitative development index of the length and width of the separation layer through the measurement of a model overlying rock section sticking measuring point; meanwhile, a subsidence measuring point is arranged on the model ground surface, and the whole process data of the ground surface subsidence gradually developing from zero to the maximum value is measured at the same time with the separation layer. Displaying dynamic growth data of surface subsidence and dynamic development data of overburden strata according to a time schedule, carrying out statistical analysis on the corresponding relation between the dynamic growth data and the dynamic development data of overburden strata, finding out the corresponding relation between the initial period, the active period and the decline period of the surface subsidence and the initiation, development and closing time points of the strata, carrying out statistics on the contribution rate of the strata in each period to the surface subsidence, and newly proposing and dividing stratum separation time nodes and each development stage; and converting the corresponding relation between the ground surface subsidence on the spot and the development of the overburden bed separation according to the experimental similarity criterion, and giving the optimal time standard for controlling the ground surface subsidence by the overburden bed separation grouting.
Although the present invention has been described in detail with reference to the foregoing embodiments, it will be apparent to those skilled in the art that changes may be made in the embodiments and/or equivalents thereof without departing from the spirit and scope of the invention. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims (2)

1. A method for analyzing the corresponding relation between a new surface subsidence dynamic process and overburden separation layer development is characterized by comprising the following steps:
measuring a sticking measuring point on a model overburden section to obtain a time process from cracking to expanding to closing of a separation layer and quantitative development indexes of the length and the width of the separation layer;
arranging a subsidence measuring point on the model ground surface, and measuring the whole process data of the ground surface subsidence gradually developing from zero to the maximum value at the same time with the separation layer;
obtaining a time schedule by a model test, displaying the dynamic growth data of the surface subsidence and the dynamic development data of the overburden bed separation, statistically analyzing the corresponding relation of the two, finding the corresponding relation of the initial period, the active period and the decline period of the surface subsidence and the initiation, development and closing time points of the bed separation, and re-proposing and dividing bed separation time nodes and each development stage;
analyzing by measuring the model overburden profile pasting measuring point, when the mining depth of the main overburden key layer on the working surface is 30-40m, the ground surface sinking speed is the fastest, the sinking speed reaches 7.5mm/d, the maximum ground surface sinking speed is synchronous with the maximum sinking speed of the key layer, the initial breakage of the key layer leads to the rapid sinking of the ground surface,
the calculation formula is as follows:
dW(t)/dt=c[W(t)/(1-1-e-ct)-Wk(1-1-e-ct)]
w (t) is a dynamic convergence value at a certain time
WkAs the final subsidence of the earth's surface
1-e-ctIs the amount of change in a certain time
dW (t)/dt is the ground surface sinking speed, when W (t)/(1-1-e)-ct)-Wk(1-1-e-ct) When the value is larger than 0, the sinking speed of the time table is increased, and when W (t)/(1-1-e)-ct)-Wk(1-1-e-ct) When the value is less than 0, the descending speed of the schedule is reduced, but the value is limited by the conditions of the slant length and the length in the model, so that the obtained value has certain deviation from the actual value;
by setting parameters of joints among the blocks of the key layer, the blocks do not move so as to simulate the movement before breaking of the rock of the key layer, when the mining width reaches the breaking distance of the key layer, the joint parameters of the blocks are reduced, the movement after breaking of the key layer is simulated, the breaking of the key layer of overlying strata influences the dynamic process of ground subsidence,
W(t)=W0+Δω(1-ae-kt)
W0is the maximum subsidence of the earth's surface at the end of the decline period
Delta omega is the residual subsidence value of the earth surface after the end of the decay period
Since W0 is a fixed number, the amount of W (t) sinking is related to a, k and t when 1-ae-ktWhen the value is larger than 0, the W (t) sinking value is increasedWhen 1-ae-ktWhen the value is less than 0, the W (t) sinking value is reduced, and the breaking of the key layer of the overburden rock causes the obvious increase of the surface sinking speed and does not change periodically.
2. The method for analyzing the correspondence between the new surface subsidence dynamic process and the overburden separation layer development according to claim 1, wherein the method comprises the following steps:
excavating according to irregular depth in the simulated coal bed, starting next excavating after model operation after each excavating is stable, when the depth is advanced to 150m, the key layer is subjected to primary breaking, primary period breaking and secondary period breaking respectively, the ground surface subsidence increment is smaller before the key layer is broken, and after the key layer is broken, the ground surface subsidence increment is increased, and the ground surface movement influence angle and the movement influence boundary are not constant,
W(x,y,t)=Wm(1-x2/a2-y2/b2)(1-e-ct)k
w (x, y, t) is the amount of sinking of any point of the sunken basin at any time
1-e-ctIs the amount of change in a certain time
(1-x2/a2-y2/b2)(1-e-ct)kWhen the product is positive, the amount of W (x, y, t) sinking increases, (1-x)2/a2-y2/b2)(1-e-ct)kWhen the product is negative, the subsidence of W (x, y, t) is reduced, and the influence boundary of the surface subsidence is obviously enlarged along with the increase of a and b.
CN201710590471.6A 2017-07-19 2017-07-19 Method for analyzing corresponding relation between new surface subsidence dynamic process and overlying strata abscission layer development Active CN107526920B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710590471.6A CN107526920B (en) 2017-07-19 2017-07-19 Method for analyzing corresponding relation between new surface subsidence dynamic process and overlying strata abscission layer development

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710590471.6A CN107526920B (en) 2017-07-19 2017-07-19 Method for analyzing corresponding relation between new surface subsidence dynamic process and overlying strata abscission layer development

Publications (2)

Publication Number Publication Date
CN107526920A CN107526920A (en) 2017-12-29
CN107526920B true CN107526920B (en) 2022-01-14

Family

ID=60748967

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710590471.6A Active CN107526920B (en) 2017-07-19 2017-07-19 Method for analyzing corresponding relation between new surface subsidence dynamic process and overlying strata abscission layer development

Country Status (1)

Country Link
CN (1) CN107526920B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103216236A (en) * 2013-03-13 2013-07-24 中国矿业大学(北京) Method for reducing subsidence of stope overburden bed separation by grouting
CN103558646A (en) * 2013-10-31 2014-02-05 中国矿业大学 Method of detecting dynamic developmental situation of overburden mining-induce fissure in radon gas earth surface
CN105067294A (en) * 2015-07-17 2015-11-18 青岛理工大学 Mining overburden bed separation different-speed growth test monitoring system and analysis method
CN106772678A (en) * 2016-12-16 2017-05-31 安徽理工大学 A kind of wellhole many reference amounts method for surveying of rock deformation destructive characteristics
CN106884657A (en) * 2017-03-31 2017-06-23 青岛理工大学 Method for displaying surface subsidence dynamic process based on overburden bed separation dynamic development

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103216236A (en) * 2013-03-13 2013-07-24 中国矿业大学(北京) Method for reducing subsidence of stope overburden bed separation by grouting
CN103558646A (en) * 2013-10-31 2014-02-05 中国矿业大学 Method of detecting dynamic developmental situation of overburden mining-induce fissure in radon gas earth surface
CN105067294A (en) * 2015-07-17 2015-11-18 青岛理工大学 Mining overburden bed separation different-speed growth test monitoring system and analysis method
CN106772678A (en) * 2016-12-16 2017-05-31 安徽理工大学 A kind of wellhole many reference amounts method for surveying of rock deformation destructive characteristics
CN106884657A (en) * 2017-03-31 2017-06-23 青岛理工大学 Method for displaying surface subsidence dynamic process based on overburden bed separation dynamic development

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于异速生长理论的采动覆岩离层时空分布规律与沉陷模型研究;胡晓阳;《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》;20160615;论文第2-3章 *

Also Published As

Publication number Publication date
CN107526920A (en) 2017-12-29

Similar Documents

Publication Publication Date Title
US11209415B2 (en) Discrete element method for modelling a fracture evolution of a roadway surrounding rock
Peng et al. Physical simulation of ground fissures triggered by underground fault activity
Bizjak et al. Displacement analysis of tunnel support in soft rock around a shallow highway tunnel at Golovec
CN107092740B (en) Prediction method for residual deformation of old goaf of steep-dip extra-thick coal seam
CN110397470A (en) It is a kind of based on crack develop gob side entry driving thin coal pillar rational width determine method
CN108978676B (en) Construction method of foundation pit group and foundation pile group
Corkum et al. Numerical analysis of longwall mining layout for a Wyoming Trona mine
Jalbout et al. Rock mechanics tools for mining in high stress ground conditions at Nickel Rim South Mine
Protosenya et al. Prediction of the stress-strain state and stability of the front of tunnel face at the intersection of disturbed zones of the soil mass
CN112488455A (en) Slope stability analysis and evaluation method, system, medium, equipment and terminal
CN111695790A (en) Mining method for security pillar
CN115903078A (en) Hydraulic fracturing layer position, determination method thereof and determination method of hydraulic fracturing position
He et al. A method for predicting the water-flowing fractured zone height based on an improved key stratum theory
CN110259443A (en) A kind of coal seam borehole wall stability prediction method based on 3DEC discrete element
CN110851991B (en) Underground water flow numerical simulation method
CN107526920B (en) Method for analyzing corresponding relation between new surface subsidence dynamic process and overlying strata abscission layer development
CN115659706B (en) Design method for roof joint cutting in coal pillar-free self-entry mining
CN110043262A (en) A kind of coal mine tight roof fractured horizontal well crack well combines monitoring method up and down
Fuenzalida et al. REBOP–FLAC3D hybrid approach to cave modelling
Ansari et al. Investigating Effect of Induced Stresses due to Coal Panel Extraction on Next Panel Strata behavior during Mechanized Longwall Mining: a Case Study
Ghazvinian et al. Back analysis of cave propagation and subsidence at Cadia East Mine
CN115324556A (en) Comprehensive prediction method for fracture-induced deformation risk level of oil-gas casing
Garrido et al. Interaction Between Tunnels and Slopes. Analysis of Weak Rock Masses Considering the Hoek-Brown Failure Criterion
Lv et al. Breaking of key layers and surface subsidence in the loess mountainous areas
e Sousa et al. Three-dimensional analyses of the Paraíso tunnel, Brazil

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