CN116485218A - Stability area dividing method and treatment time planning method for subsidence area - Google Patents
Stability area dividing method and treatment time planning method for subsidence area Download PDFInfo
- Publication number
- CN116485218A CN116485218A CN202211667763.2A CN202211667763A CN116485218A CN 116485218 A CN116485218 A CN 116485218A CN 202211667763 A CN202211667763 A CN 202211667763A CN 116485218 A CN116485218 A CN 116485218A
- Authority
- CN
- China
- Prior art keywords
- time
- period
- mining
- subsidence area
- deformation
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000005065 mining Methods 0.000 claims abstract description 154
- 238000012544 monitoring process Methods 0.000 claims description 28
- 238000011084 recovery Methods 0.000 claims description 24
- 238000011160 research Methods 0.000 claims description 16
- 238000000638 solvent extraction Methods 0.000 claims description 3
- 239000003245 coal Substances 0.000 abstract description 48
- 238000012360 testing method Methods 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 230000007774 longterm Effects 0.000 description 4
- 239000002689 soil Substances 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000005192 partition Methods 0.000 description 3
- 238000005067 remediation Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 201000004569 Blindness Diseases 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000011449 brick Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000000979 retarding effect Effects 0.000 description 1
- 238000012502 risk assessment Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- 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/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06313—Resource planning in a project environment
-
- 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
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/12—Timing analysis or timing optimisation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Educational Administration (AREA)
- Tourism & Hospitality (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Life Sciences & Earth Sciences (AREA)
- Agronomy & Crop Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Biodiversity & Conservation Biology (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
The invention provides a stability area dividing method and a treatment time planning method for a subsidence area, wherein the method comprises the following steps: fitting the earth surface movement deformation of the earth surface movement observation station and the subsidence area sample monitored by the InSAR based on the Bertalanffy time function, and solving a functional relation equation of the earth surface movement deformation and the time; dividing the period of the earth surface movement deformation into a plurality of time phases based on a functional relation equation; determining the time corresponding relation between the mining stage and the time stage of the mine working face sample corresponding to the subsidence area sample according to the mining engineering plan of the subsidence area sample; and determining the stability type of the current subsidence area according to the time corresponding relation between the time phase and the mining phase of the current mine working face. The method and the device divide the effective stability area of the coal mining subsidence area, have high division accuracy, and provide references for the treatment time planning of the coal mining subsidence area.
Description
Technical Field
The invention relates to the technical field of subsidence area treatment, in particular to a subsidence area stability area dividing method and a treatment time planning method.
Background
After the underground coal mine is mined, the formed goaf causes the stress balance state of surrounding original rock to be destroyed, and the movement deformation of the overlying strata and the surface soil layer is propagated layer by layer and causes the ground surface to sink. Mining subsidence in mining areas can cause damage to surrounding underground and ground surface structures, so that loss of fixed property and damage to the environment are caused, and economic sustainable development of the mining areas is affected.
The rock-soil layer movement deformation is a complex and periodic movement process caused by underground exploitation, and the observation of surface movement and deformation by arranging a surface movement observation station (net) is the most direct method for knowing the complex surface movement process. Conventional earth surface movement observation station burying and observation techniques typically monitor earth surface movement deformation in a mainly concentrated deformation phase, i.e., typically only the initial, active and decay phases of earth surface movement deformation. For the earth surface movement deformation duration time which is large in mining depth and is influenced by repeated mining operations of multiple coal seam spaces, the conventional earth surface movement observation station monitoring technology needs to input long-time manpower and material resources, has poor economic benefit, is difficult to store and monitor long-term monitoring points (the common monitoring period is 1 to 2 years), and has high cost, long period and great difficulty.
Meanwhile, for the moving deformation of an overburden soil layer with large mining depth and repeated disturbance of multiple coal seam spaces, an effective coal mining subsidence area dividing method and a time analysis method do not exist, so that comprehensive treatment and planning of the coal mining subsidence area are scientifically and effectively guided, larger area dividing and treatment time blindness is brought to planning and utilization of reclamation, treatment and the like of the coal mining subsidence area, and even larger moving deformation values still exist in an unstable stage of the subsidence area, so that unnecessary secondary reclamation treatment is caused.
The standards of building, water, railway and main roadway coal pillar reserving and coal pressing exploitation recommend the levels of damage to brick and concrete building (construction) caused by different mining operations and the corresponding surface movement deformation values, and divide the surface movement deformation levels into I, II, III, IV four levels. The goaf is divided into a place which is not suitable for construction and a place which is relatively stable according to the characteristics and the deformation of the ground surface subsidence basin according to the exploitation condition. The geological disaster assessment under the influence of underground mining in geological disaster risk assessment norms is limited to judging the influence degree of mining on the ground surface according to the size of the movement deformation characteristic value of the mining ground surface, and is divided into three grades, namely strong, stronger and non-strong. Meanwhile, the division is carried out according to the classification method based on the basis of the stability division of the coal mining subsidence area.
For the treatment time planning of the coal mining subsidence area, the earth surface movement duration parameters of the mining area and the peripheral mining areas are often referred to firstly, when the earth surface can reach stability is judged, and the treatment time is arranged according to the earth surface movement cycle rule; when the earth surface movement period rule is not obtained, the earth surface is estimated to be stable according to experience in a fuzzy manner, and the reclamation and treatment work of the coal mining subsidence area is carried out after the earth surface movement deformation is stable.
Disclosure of Invention
The invention provides a stability area dividing method and a treatment time planning method for a subsidence area, which are used for solving the defects that in the prior art, an upper rock-covered soil layer repeatedly disturbed by a large mining depth and a multi-coal-seam space moves and deforms, and an effective coal mining subsidence area stability area dividing method is not available, and realizing the stability area dividing and the treatment time planning of the subsidence area based on InSAR and observation stations.
The invention provides a method for dividing stability areas of a subsidence area, which comprises the following steps:
fitting the earth surface movement deformation of the earth surface movement observation station and the subsidence area sample monitored by the InSAR based on the Bertalanffy time function, and solving a functional relation equation of the earth surface movement deformation and the time;
dividing the period of the earth surface movement deformation into a plurality of time phases based on a functional relation equation of the earth surface movement deformation and time, wherein the time phases comprise an initial phase, an active phase, a decay phase and a residual deformation phase;
determining the time corresponding relation between the mining stage of the mine working face sample corresponding to the subsidence area sample and the time stage according to the mining engineering plan of the subsidence area sample;
and determining the stability type of the current subsidence area according to the mining stage of the current mine working face corresponding to the current subsidence area and the time corresponding relation between the time stage and the mining stage of the mine working face sample.
According to the method for dividing the stability area of the subsidence area, disclosed by the invention, the earth surface movement deformation of the subsidence area sample monitored by the earth surface movement observation station and the InSAR is fitted based on the Bertalanffy time function, and the step of solving the function relation equation of the earth surface movement deformation and the time comprises the following steps:
monitoring the earth surface movement deformation of the subsidence area sample in the initial period and the active period through an earth surface movement observation station;
monitoring the earth surface movement deformation of the subsidence area sample in the decay period and the residual deformation period through InSAR;
splicing the ground surface movement observation station with the ground surface movement deformation monitored by the InSAR;
fitting the spliced earth surface movement deformation based on the Bertalaffy time function, and solving a functional relation equation of the earth surface movement deformation and time.
According to the method for dividing the stability area of the subsidence area provided by the invention, the step of dividing the period of the earth surface movement deformation into a plurality of time phases based on the function relation equation of the earth surface movement deformation and time comprises the following steps:
determining a function curve of the subsidence amount and time of the maximum subsidence point of the earth surface of the subsidence area sample and a function curve of the subsidence speed and time of the maximum subsidence point of the earth surface;
determining an initial period, an active period and a decay period of the ground surface movement deformation according to the function curve of the sinking amount and the time and the function curve of the sinking speed and the time;
and determining the ending time of the decay period according to the function relation equation of the earth surface movement deformation and time, and taking the ending time of the decay period as the starting time of the residual deformation period.
According to the method for dividing the stability area of the subsidence area, which is provided by the invention, the step of determining the ending time of the decay period according to the function relation equation of the ground surface movement deformation and time comprises the following steps:
and if the continuous preset duration of the earth surface movement deformation after entering the decay period is smaller than a preset threshold value according to the function relation equation of the earth surface movement deformation and time, taking the ending time corresponding to the first preset duration as the ending time of the decay period.
According to the method for dividing the stability area of the subsidence area provided by the invention, the step of determining the time corresponding relation between the mining stage of the mine working face sample corresponding to the subsidence area sample and the time stage according to the mining engineering plan of the subsidence area sample comprises the following steps:
determining that the initial period and the active period of the surface movement deformation correspond to a first time period from a ground research point at which the mine working surface sample starts to be subjected to stope disturbance propagation to the subsidence area sample to a point at which the mine working surface sample is pushed through the ground research point according to the mining engineering plan of the subsidence area sample;
determining that the decay period of the surface movement deformation corresponds to a period from a second period after the mining of the mine working face sample is finished to a third period after the mining of the mine working face sample is finished, wherein the second period is smaller than the third period;
and determining that the residual deformation period of the surface movement deformation corresponds to a third time period after the mining of the mine working face sample is finished.
According to the method for dividing the stability area of the subsidence area provided by the invention, the step of determining the stability type of the current subsidence area according to the mining stage of the current mine working face corresponding to the current subsidence area and the time corresponding relation between the time stage and the mining stage of the mine working face sample comprises the following steps:
if the mining stage of the current mine working face is from a ground research point at which mining disturbance is transmitted to the subsidence area sample from the beginning of the mine working face sample to a first time period at which the mine working face sample is pushed through the ground research point, determining that the stability type of the current subsidence area is an unstable area according to the initial period and the active period of the mining stage;
if the mining stage of the current mine working face is from the second time length after the mining of the mine working face sample is finished to the third time length after the mining of the mine working face sample is finished, the mining stage corresponds to the decay period, and the stability type of the current subsidence area is determined to be a basic stable area;
and if the mining stage of the current mine working face is a third time period after the mining of the mine working face sample is finished, determining that the stability type of the current subsidence area is a residual deformation area according to the residual deformation period corresponding to the mining stage.
The invention also provides a treatment time planning method based on the stability area dividing method of the subsidence area, which comprises the following steps:
when the stability type of the current subsidence area is an unstable area, performing no ground vegetation recovery and treatment work on the current subsidence area;
when the stability type of the current subsidence area is a basic stable area, selectively carrying out surface vegetation recovery and treatment work on the current subsidence area;
and when the stability type of the current subsidence area is a residual deformation area, carrying out surface vegetation recovery and treatment work on the current subsidence area.
According to the subsidence area stability region dividing method and the treatment time planning method, the surface movement observation station is used for acquiring high-precision large deformation and severe deformation monitoring characteristics of a short period and a medium period, and the InSAR is used for observing small deformation, long period and slow deformation characteristics of a border region of a coal mining subsidence area in a range of the subsidence area, so that a surface movement deformation and time function relation equation with large mining depth, long time span, full life cycle and high monitoring precision under the influence of space repeated mining is established, and the surface movement time stage division and working face mining stage division are carried out on the movement deformation period of the subsidence area based on the relation equation, so that the effective stability region division is carried out on the subsidence area based on a time corresponding relation between the surface movement time stage and the mining stage, the division accuracy is high, and a reference is provided for the treatment time planning of the coal mining subsidence area.
Drawings
In order to more clearly illustrate the invention or the technical solutions of the prior art, the following description will briefly explain the drawings used in the embodiments or the description of the prior art, and it is obvious that the drawings in the following description are some embodiments of the invention, and other drawings can be obtained according to the drawings without inventive effort for a person skilled in the art.
FIG. 1 is a schematic flow chart of a method for dividing stability areas of a subsidence area;
FIG. 2 is a graph showing the surface movement deformation and time function relation curve fitting of the full life cycle in the subsidence area stability zone division method provided by the invention;
FIG. 3 is a schematic diagram of the change curve and the phase division of the relationship between the ground surface movement deformation and the time in the stability area division method of the subsidence area;
FIG. 4 is a schematic diagram of stability zones for recovery and remediation of vegetation on the surface of a coal mining subsidence area in a coal mine in a subsidence area stability zone division method provided by the invention;
FIG. 5 is a schematic flow chart of a treatment time planning method based on a stability zone division method of a subsidence area.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the present invention more apparent, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings, and it is apparent that the described embodiments are some embodiments of the present invention, not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
A method for dividing stability area of a subsidence area according to the present invention is described below with reference to fig. 1, including:
step 101, fitting the earth surface movement deformation of an earth surface movement observation station and an InSAR monitored subsidence area sample based on a Bertalanffy time function, and solving a functional relation equation of the earth surface movement deformation and time;
and combining the upper and lower control charts of the coal mine, and taking a working surface sample of the subsidence area sample as a test working surface to be researched to design a ground surface movement observation station, wherein the design comprises the number of measurement points, the distance, the observation times and the like.
And burying earth surface movement observation stations with observation piles on the earth surface corresponding to the test working surface, and periodically observing plane coordinates and elevation coordinates of earth surface movement deformation monitoring points to obtain high-precision large deformation and severe deformation monitoring characteristics of the subsidence area samples in a short term, a medium term and a long term.
And determining an InSAR data acquisition area comprising a test working surface and a residual deformation working surface of the acquired area by combining satellite image orbit coordinates.
And observing the characteristics of small deformation and long-period slow deformation of the boundary region of the subsidence area in the sample range of the coal mining subsidence area through InSAR.
And analyzing the measured data of the earth surface movement observation station based on a probability integration method, and fitting to obtain earth surface movement deformation calculation parameters.
And calculating the earth surface movement angle value parameter according to the earth surface movement deformation calculation parameter, and acquiring earth surface subsidence rules, earth surface stabilization time characteristics and other earth surface movement deformation characteristics under the geological mining condition according to the earth surface movement angle value parameter.
The earth surface movement observation station is used for mainly monitoring the earth surface movement deformation excitation in the concentrated deformation stage, namely, only monitoring the initial period, the active period and the decay period of the earth surface movement deformation, wherein the monitoring period is 1-2 years. InSAR can retrospectively acquire historical data of 5 to 10 years and residual deformation observation with small quantity value in a later period. And combining the earth surface movement observation station with the earth surface movement deformation monitored by the InSAR to obtain the earth surface movement deformation of the whole life cycle of the mine.
Fitting the data monitored by the earth surface movement observation station and the InSAR, and establishing an earth surface movement deformation and time function relation equation with large mining depth, long time span under the influence of space repeated mining, full life cycle and high monitoring precision.
Taking the conditions of large mining depth, multiple coal seam spaces and space repeated mining operation of a certain coal mine as an example, a surface movement deformation and time function relation curve with longer time span, full life cycle and high monitoring precision under the influence of space repeated mining operation and larger mining depth is established based on MATLAB (matrix laboratory), and is shown in figure 2. Wherein the axis of abscissa is time in days and the axis of ordinate is the amount of subsidence in millimeters. The established full life cycle earth surface movement deformation and time function relation equation is as follows:
W＝-2667.4·(1-0.9940·e -0.0243·T ) 48 。
102, dividing the period of the earth surface movement deformation into a plurality of time phases based on a functional relation equation of the earth surface movement deformation and time, wherein the time phases comprise an initial phase, an active phase, a decay phase and a residual deformation phase;
based on a functional relation equation of the ground surface movement deformation and time, the time T1 when the ground surface movement deformation enters an initial period, the time T2 when the ground surface movement deformation enters an active period, the time T3 when the ground surface movement deformation enters a decay period and the time T4 when the ground surface movement deformation enters a residual deformation period are determined.
The period of the earth moving deformation is divided into a plurality of time phases, namely T0 to T1, T1 to T2, T2 to T3 and T3 to T4, corresponding to the initial phase, the active phase, the decay phase and the residual deformation phase of the earth moving deformation.
Step 103, determining a time corresponding relation between a mining stage of a mine working face sample corresponding to the subsidence area sample and the time stage according to a mining engineering plan of the subsidence area sample;
and determining how much time or how far ahead the mining of the mine working face sample ends when the mine working face sample enters the initial period, the active period, the decay period and the residual deformation period by combining the mining engineering plan of the subsidence area sample, so that the mining of the mine working face sample is divided into a plurality of mining stages. Such as the distance of the earth's surface from the retarding face to the advancing face for each production stage.
Each time phase of the earth's surface movement deformation corresponds to a respective production phase.
And 104, determining the stability type of the current subsidence area according to the mining stage of the current mine working face corresponding to the current subsidence area and the time corresponding relation between the time stage and the mining stage of the mine working face sample.
And determining the partition basis and partition category of the suitability of the surface vegetation restoration and treatment work of the coal mining subsidence area according to the technical requirements of reclamation and comprehensive treatment of the coal mining subsidence area and the control requirement of the future ground surface movement deformation, and partitioning according to the basis.
And determining the earth surface movement deformation time phase corresponding to the mining phase of the current mine working surface according to the time corresponding relation between the mining phase of the current mine working surface corresponding to the current subsidence area and the mining phase of the mine working surface sample.
And determining the stability type corresponding to the time phase of the surface movement deformation according to the time phase of the surface movement deformation corresponding to the mining phase of the current mine working face, wherein the time phase of the surface movement deformation and the stability type are associated in advance.
Optionally, in combination with the related regulation specification, if the mining stage of the current mine working face corresponds to the initial period and the active period of the surface movement deformation, the stability type of the current subsidence area is an unstable area, namely an unsuitable vegetation recovery and treatment area.
If the mining stage of the current mine working face corresponds to the ground surface movement deformation decay stage, the stability type of the current subsidence area is a basic stable area, namely a vegetation recovery and treatment area with proper selectivity. Whether vegetation recovery and treatment are carried out or not can be selected according to actual conditions and requirements.
If the current mining stage of the mine working face corresponds to the residual deformation area, the stability type of the current subsidence area is the residual deformation area, namely the vegetation recovery and treatment area is suitable.
According to the embodiment, the high-precision large deformation and severe deformation monitoring characteristics of a short period, a medium period and a long period are obtained through the earth surface movement observation station, the small deformation, long period and slow deformation characteristics of the boundary area of the coal mining subsidence area are observed through the InSAR, the earth surface movement deformation and time function relation equation with large mining depth, long time span under the influence of space repeated mining, full life period and high monitoring precision is established, the earth surface movement deformation period of the subsidence area is divided into earth surface movement time phases and working surface mining phases based on the relation equation, and therefore the effective stability area division is carried out on the coal mining subsidence area based on the time corresponding relation between the earth surface movement time phases and the mining phases, and the division accuracy is high, so that references are provided for the treatment time planning of the coal mining subsidence area.
On the basis of the above embodiment, in this embodiment, the step of fitting the earth surface movement deformation of the earth surface movement observation station and the subsidence area sample monitored by the InSAR based on the Bertalanffy time function, and obtaining the function relation equation of the earth surface movement deformation and time includes:
monitoring the earth surface movement deformation of the subsidence area sample in the initial period and the active period through an earth surface movement observation station;
the monitoring of the ground surface mobile observation station needs to be put into manpower and material resources for a long time, the economic benefit is poor, long-term monitoring point storage and long-term monitoring are difficult to achieve, while the InSAR can carry out long-period ground surface small deformation monitoring, and has the advantages of low cost, reliable data and convenience in processing. Therefore, in the present embodiment, the earth movement deformation in the initial period and the active period, that is, the case where the earth movement deformation amount is large is monitored by the earth movement observation station.
Monitoring the earth surface movement deformation of the subsidence area sample in the decay period and the residual deformation period through InSAR;
under the condition that the movement deformation of the ground surface is relatively low, the historical deformation data of the surface corresponding to the collected working surface on the periphery of the test working surface and the test working surface are analyzed through the monitoring data of the InSAR, and the deformation characteristics of the movement deformation recession period and the residual deformation stage of the ground surface are mastered.
Splicing the ground surface movement observation station with the ground surface movement deformation monitored by the InSAR;
and combining the earth surface movement observation station with the earth surface movement deformation monitored by the InSAR to construct the earth surface movement deformation of the whole life cycle of the mine.
Fitting the spliced earth surface movement deformation based on the Bertalaffy time function, and solving a functional relation equation of the earth surface movement deformation and time.
Fitting the data monitored by the earth surface movement observation station and the InSAR, and establishing an earth surface movement deformation and time function relation equation with large mining depth, long time span under the influence of space repeated mining, full life cycle and high monitoring precision.
On the basis of the foregoing embodiment, the step of dividing the period of the surface movement deformation into a plurality of time phases based on the functional relation equation of the surface movement deformation and time in this embodiment includes:
determining a function curve of the subsidence amount and time of the maximum subsidence point of the earth surface of the subsidence area sample and a function curve of the subsidence speed and time of the maximum subsidence point of the earth surface;
taking the conditions of large mining depth, multiple coal seam spaces and space repeated mining operation of a certain coal mine as an example, a surface movement deformation and time function relation change curve and deformation stage division with large mining depth, long time span under the influence of space repeated mining operation, full life cycle and high monitoring precision are established as shown in figure 3.
In fig. 3, W (t) represents a time function curve of the subsidence amount of the maximum subsidence point of the earth surface, V (t) represents a time function curve of the subsidence speed of the maximum subsidence point of the earth surface, L (t) represents a relative position of the maximum subsidence point of the earth surface and the working surface, and a time function curve, a negative value represents a distance that the working surface is located in front of the working surface when the working surface is not mined to pass through the working surface, and a positive value represents a distance that the working surface is pushed through the working surface.
Determining an initial period, an active period and a decay period of the ground surface movement deformation according to the function curve of the sinking amount and the time and the function curve of the sinking speed and the time;
the initial period, the active period and the decay period of the ground surface movement deformation are divided by combining the ' three-lower ' coal mining standard '.
And determining the ending time of the decay period according to the function relation equation of the earth surface movement deformation and time, and taking the ending time of the decay period as the starting time of the residual deformation period.
And determining the end time of the decay period based on the established earth surface movement deformation and time function relation equation with large mining depth, long time span under the influence of space repeated mining, full life cycle and high monitoring precision, and determining the subsequent residual deformation of the decay period as the residual deformation period.
On the basis of the foregoing embodiment, the step of determining the end time of the decay period according to the equation of the functional relation between the surface movement deformation and time in this embodiment includes:
and if the continuous preset duration of the earth surface movement deformation after entering the decay period is smaller than a preset threshold value according to the function relation equation of the earth surface movement deformation and time, taking the ending time corresponding to the first preset duration as the ending time of the decay period.
Based on a functional relation equation of the ground surface movement deformation and time, finding out the condition that the first duration preset time length of the ground surface movement deformation is smaller than a preset threshold value after the decay period of the ground surface movement deformation, and taking the ending time corresponding to the first preset time length under the condition as the ending time of the decay period. The subsequent residual deformation after the end of the decay period is referred to as the residual deformation period.
On the basis of the foregoing embodiments, in this embodiment, the step of determining, according to the plan view of the mining engineering of the subsidence area sample, a time correspondence between the mining stage and the time stage of the mine working face sample corresponding to the subsidence area sample includes:
determining that the initial period and the active period of the surface movement deformation correspond to a first time period from a ground research point at which the mine working surface sample starts to be subjected to stope disturbance propagation to the subsidence area sample to a point at which the mine working surface sample is pushed through the ground research point according to the mining engineering plan of the subsidence area sample;
the first time length is t 1 To t 2 And (3) determining the first time according to the initial period of the mobile deformation and the time corresponding to the active period.
Mine working face sample self-openingThe initial recovery disturbance propagates to the ground research point and the mine working face sample is pushed through the ground research point t 1 To t 2 Within a month, the range from the front L1 of the corresponding working face to the L1 of the working face after the working face is pushed is the initial period and the active period of the ground surface movement deformation.
Determining that the decay period of the surface movement deformation corresponds to a period from a second period after the mining of the mine working face sample is finished to a third period after the mining of the mine working face sample is finished, wherein the second period is smaller than the third period;
the second time length is t 3 To t 4 For one month. T after working face stoping is finished 3 To t 4 Month to end of recovery t 5 The annual period is the period of movement, deformation and recession of the earth surface.
And determining that the residual deformation period of the surface movement deformation corresponds to a third time period after the mining of the mine working face sample is finished.
Working face mining end t 5 After the year, the earth surface enters a residual deformation period with small deformation value and long duration.
On the basis of the foregoing embodiment, in this embodiment, the step of determining the stability type of the current subsidence area according to the mining phase of the current mine working face corresponding to the current subsidence area and the time correspondence between the time phase and the mining phase of the mine working face sample includes:
if the mining stage of the current mine working face is from a ground research point at which mining disturbance is transmitted to the subsidence area sample from the beginning of the mine working face sample to a first time period at which the mine working face sample is pushed through the ground research point, determining that the stability type of the current subsidence area is an unstable area according to the initial period and the active period of the mining stage;
if the mining stage of the current mine working face is from the second time length after the mining of the mine working face sample is finished to the third time length after the mining of the mine working face sample is finished, the mining stage corresponds to the decay period, and the stability type of the current subsidence area is determined to be a basic stable area;
and if the mining stage of the current mine working face is a third time period after the mining of the mine working face sample is finished, determining that the stability type of the current subsidence area is a residual deformation area according to the residual deformation period corresponding to the mining stage.
Taking the conditions of large mining depth, multiple coal seam space and space repeated mining operation of a certain coal mine as an example, and combining with the plan view of the coal mine mining engineering and the stability time analysis of the surface movement deformation, the stability partition basis of the recovery and treatment of the surface vegetation of the coal mine mining subsidence area can be obtained as follows:
(1) unsuitable vegetation recovery and remediation zone (unstable zone): the working face is in an initial period and an active period of surface movement deformation within 2.5 to 3 months from the beginning of stoping disturbance transmission to the ground research point to the pushing of the working face through the ground research point (corresponding to the range from 190m in front of the working face to 190m after the pushing of the working face), and the surface movement deformation value is large, the deformation is severe in the period, and the method is not suitable for the surface vegetation recovery and treatment work of a coal mining subsidence area.
(2) Selectively suitable vegetation recovery and remediation zone (substantially stable zone): the period from 2.5 to 3 months after the stoping of the working face is finished to 2 years after the stoping is finished is the ground surface movement deformation decay period, and the ground surface vegetation recovery and treatment work of the coal mining subsidence area can be selectively carried out in the period.
(3) Suitable vegetation recovery and treatment area (residual deformation area): after the working face is exploited for 2 years, the earth surface enters a residual deformation period with a tiny deformation value and a long duration, and the earth surface vegetation recovery and treatment work of the coal mining subsidence area can be carried out in the period and later.
And combining the mine area up-and-down comparison charts, calculating parameters, angle parameters and the like based on an earth surface movement deformation probability integration method, calculating earth surface movement deformation values of the coal mining subsidence areas based on corresponding earth surface subsidence calculation auxiliary software, defining earth surface vegetation restoration and treatment suitability areas of each coal mining subsidence area, and providing corresponding comprehensive treatment time schedule suggestions according to earth surface movement deformation characteristics.
Taking the repeated mining conditions of large mining depth and multiple coal seam space of a certain coal mine as an example, the surface vegetation recovery and treatment suitability areas of each coal mining subsidence area are defined as shown in fig. 4. Referring to fig. 4, the comprehensive treatment of the coal mining subsidence area is scientifically arranged in regions and steps, such as personnel, machinery, backfill and the like, are prepared, and the comprehensive treatment of the coal mining subsidence area is sequentially implemented.
The embodiment can scientifically and effectively guide comprehensive treatment and planning of the coal mining subsidence area, and scientific guidance and reasonable planning are provided for regional division and treatment time in the planning and utilization processes of reclamation, treatment and the like of the coal mining subsidence area.
As shown in fig. 5, the present embodiment provides a method for planning treatment time of a stability zone division method of a subsidence area, including:
step 501, when the stability type of the current subsidence area is an unstable area, performing no surface vegetation recovery and treatment work on the current subsidence area;
step 502, selectively performing surface vegetation recovery and treatment work on the current subsidence area when the stability type of the current subsidence area is a basic stable area;
and step 503, when the stability type of the current subsidence area is the residual deformation area, performing surface vegetation recovery and treatment work on the current subsidence area.
The embodiment can scientifically and effectively guide comprehensive treatment and planning of the coal mining subsidence area, and scientific guidance and reasonable planning are provided for regional division and treatment time in the planning and utilization processes of reclamation, treatment and the like of the coal mining subsidence area.
Finally, it should be noted that: the above embodiments are only for illustrating the technical solution of the present invention, and are not limiting; although the invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit and scope of the technical solutions of the embodiments of the present invention.
Claims (7)
1. A method for dividing stability areas of a subsidence area, comprising:
fitting the earth surface movement deformation of the earth surface movement observation station and the subsidence area sample monitored by the InSAR based on the Bertalanffy time function, and solving a functional relation equation of the earth surface movement deformation and the time;
dividing the period of the earth surface movement deformation into a plurality of time phases based on a functional relation equation of the earth surface movement deformation and time, wherein the time phases comprise an initial phase, an active phase, a decay phase and a residual deformation phase;
determining the time corresponding relation between the mining stage of the mine working face sample corresponding to the subsidence area sample and the time stage according to the mining engineering plan of the subsidence area sample;
and determining the stability type of the current subsidence area according to the mining stage of the current mine working face corresponding to the current subsidence area and the time corresponding relation between the time stage and the mining stage of the mine working face sample.
2. The method for dividing stability areas of a subsidence area according to claim 1, wherein the step of fitting the earth movement deformation of the subsidence area sample monitored by the earth movement observation station and the InSAR based on the Bertalanffy time function, and solving a functional relation equation of the earth movement deformation and time comprises:
monitoring the earth surface movement deformation of the subsidence area sample in the initial period and the active period through an earth surface movement observation station;
monitoring the earth surface movement deformation of the subsidence area sample in the decay period and the residual deformation period through InSAR;
splicing the ground surface movement observation station with the ground surface movement deformation monitored by the InSAR;
fitting the spliced earth surface movement deformation based on the Bertalaffy time function, and solving a functional relation equation of the earth surface movement deformation and time.
3. The method of claim 1, wherein the step of dividing the period of the surface movement deformation into a plurality of time phases based on the surface movement deformation and a time functional relation equation comprises:
determining a function curve of the subsidence amount and time of the maximum subsidence point of the earth surface of the subsidence area sample and a function curve of the subsidence speed and time of the maximum subsidence point of the earth surface;
determining an initial period, an active period and a decay period of the ground surface movement deformation according to the function curve of the sinking amount and the time and the function curve of the sinking speed and the time;
and determining the ending time of the decay period according to the function relation equation of the earth surface movement deformation and time, and taking the ending time of the decay period as the starting time of the residual deformation period.
4. A method of zonal stability partitioning of a subsidence area as set forth in claim 3 wherein said step of determining the end time of said decay period from a functional relationship equation of said surface movement deformation and time comprises:
and if the continuous preset duration of the earth surface movement deformation after entering the decay period is smaller than a preset threshold value according to the function relation equation of the earth surface movement deformation and time, taking the ending time corresponding to the first preset duration as the ending time of the decay period.
5. The method for dividing stability areas of a subsidence area according to any one of claims 1 to 4, wherein the step of determining a time correspondence between a mining stage of a mine working face sample corresponding to the subsidence area sample and the time stage according to a mining engineering plan of the subsidence area sample comprises:
determining that the initial period and the active period of the surface movement deformation correspond to a first time period from a ground research point at which the mine working surface sample starts to be subjected to stope disturbance propagation to the subsidence area sample to a point at which the mine working surface sample is pushed through the ground research point according to the mining engineering plan of the subsidence area sample;
determining that the decay period of the surface movement deformation corresponds to a period from a second period after the mining of the mine working face sample is finished to a third period after the mining of the mine working face sample is finished, wherein the second period is smaller than the third period;
and determining that the residual deformation period of the surface movement deformation corresponds to a third time period after the mining of the mine working face sample is finished.
6. The method of claim 5, wherein the step of determining the stability type of the current subsidence area according to the mining phase of the current mine face corresponding to the current subsidence area and the time correspondence between the time phase and the mining phase of the mine face sample comprises:
if the mining stage of the current mine working face is from a ground research point at which mining disturbance is transmitted to the subsidence area sample from the beginning of the mine working face sample to a first time period at which the mine working face sample is pushed through the ground research point, determining that the stability type of the current subsidence area is an unstable area according to the initial period and the active period of the mining stage;
if the mining stage of the current mine working face is from the second time length after the mining of the mine working face sample is finished to the third time length after the mining of the mine working face sample is finished, the mining stage corresponds to the decay period, and the stability type of the current subsidence area is determined to be a basic stable area;
and if the mining stage of the current mine working face is a third time period after the mining of the mine working face sample is finished, determining that the stability type of the current subsidence area is a residual deformation area according to the residual deformation period corresponding to the mining stage.
7. A method of governance time planning based on the method of stability zone partitioning of a subsidence area as claimed in any one of claims 1 to 6, comprising:
when the stability type of the current subsidence area is an unstable area, performing no ground vegetation recovery and treatment work on the current subsidence area;
when the stability type of the current subsidence area is a basic stable area, selectively carrying out surface vegetation recovery and treatment work on the current subsidence area;
and when the stability type of the current subsidence area is a residual deformation area, carrying out surface vegetation recovery and treatment work on the current subsidence area.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211667763.2A CN116485218A (en) | 2022-12-23 | 2022-12-23 | Stability area dividing method and treatment time planning method for subsidence area |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211667763.2A CN116485218A (en) | 2022-12-23 | 2022-12-23 | Stability area dividing method and treatment time planning method for subsidence area |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116485218A true CN116485218A (en) | 2023-07-25 |
Family
ID=87214346
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211667763.2A Pending CN116485218A (en) | 2022-12-23 | 2022-12-23 | Stability area dividing method and treatment time planning method for subsidence area |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116485218A (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103091676A (en) * | 2013-01-22 | 2013-05-08 | 中国矿业大学 | Mining area surface subsidence synthetic aperture radar interferometry monitoring and calculating method |
| US20130185040A1 (en) * | 2012-01-13 | 2013-07-18 | Westerngeco L.L.C. | Parameterizing a geological subsurface feature |
| CN105926569A (en) * | 2016-04-28 | 2016-09-07 | 河北地质大学 | Method for quantitatively evaluating site stability of old goaf in coalmine based on settlement monitoring data |
| CN111750822A (en) * | 2019-10-31 | 2020-10-09 | 河南理工大学 | Coal mining induced overlying strata and surface subsidence collaborative dynamic prediction method |
| CN114674277A (en) * | 2022-03-17 | 2022-06-28 | 辽宁大学 | Method for monitoring surface subsidence of deep mining and evaluating land reclamation suitability by combining whole-field survey lines |
| CN115096257A (en) * | 2022-06-22 | 2022-09-23 | 平安煤炭开采工程技术研究院有限责任公司 | Mining area mining subsidence monitoring method and device |
-
2022
- 2022-12-23 CN CN202211667763.2A patent/CN116485218A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130185040A1 (en) * | 2012-01-13 | 2013-07-18 | Westerngeco L.L.C. | Parameterizing a geological subsurface feature |
| CN103091676A (en) * | 2013-01-22 | 2013-05-08 | 中国矿业大学 | Mining area surface subsidence synthetic aperture radar interferometry monitoring and calculating method |
| CN105926569A (en) * | 2016-04-28 | 2016-09-07 | 河北地质大学 | Method for quantitatively evaluating site stability of old goaf in coalmine based on settlement monitoring data |
| CN111750822A (en) * | 2019-10-31 | 2020-10-09 | 河南理工大学 | Coal mining induced overlying strata and surface subsidence collaborative dynamic prediction method |
| CN114674277A (en) * | 2022-03-17 | 2022-06-28 | 辽宁大学 | Method for monitoring surface subsidence of deep mining and evaluating land reclamation suitability by combining whole-field survey lines |
| CN115096257A (en) * | 2022-06-22 | 2022-09-23 | 平安煤炭开采工程技术研究院有限责任公司 | Mining area mining subsidence monitoring method and device |
Non-Patent Citations (2)
| Title |
|---|
| 陈炳乾: ""面向矿区沉降监测的InSAR技术及应用研究"", 《中国博士学位论文全文数据库 工程科技Ⅰ辑》, no. 4, pages 021 - 3 * |
| 高超等: ""基于Bertalanffy时间函数的地表动态沉陷预测模型"", 《煤炭学报》, vol. 45, no. 8, pages 2740 - 2748 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Marschall et al. | Hydro-mechanical evolution of the EDZ as transport path for radionuclides and gas: insights from the Mont Terri rock laboratory (Switzerland) | |
| CN110043312B (en) | Control method for grouting filling ground surface settlement range | |
| Masset et al. | Quantitative hydraulic analysis of pre-drillings and inflows to the Gotthard Base Tunnel (Sedrun Lot, Switzerland) | |
| Swift | Relationship between joint movement and mining subsidence | |
| Guney et al. | Analysis of surface subsidence due to longwall mining under weak geological conditions: Turgut basin of Yatağan-Muğla (Turkey) case study | |
| CN110145309B (en) | Shallow coal seam grouting filling green coal mining method and settlement observation device thereof | |
| Yang et al. | Deformation patterns and failure mechanism of high and steep stratified rock slopes with upper steep and lower gentle style induced by step-by-step excavations | |
| CN113742995A (en) | Mine water inflow prediction method and system based on coal mine big data | |
| Konyukhov | Analysis of mechanized tunneling parameters to determine the overcutting characteristics | |
| CN110851991B (en) | Underground water flow numerical simulation method | |
| CN116485218A (en) | Stability area dividing method and treatment time planning method for subsidence area | |
| CN116756807A (en) | Linear underground engineering water inflow numerical prediction method in coupling construction process | |
| Vervoort | Impact of the closure of a coal district on the environmental issue of long-term surface movements | |
| RU2067166C1 (en) | Method for development of oil deposit in tectonically complicated sedimentary strata | |
| CN116910886B (en) | Impact energy early warning method and system for karst region | |
| Krasilnikov | GEOINFORMATIONAL APPROACH TO SOLVING FORECASTING ENGINEERING AND GEOLOGICAL TASKS IN THE DEVELOPMENT OF MINERAL DEPOSITS | |
| CN108645552B (en) | Ground stress measurement method based on double-stress system | |
| Kannan et al. | Geological assessment of underground storage cavern sites | |
| Rapisarda | Landslide analysis in Apennine chain areas | |
| Eftekhari et al. | Complexity of the ground conditions and non compliance with basic assumptions in the trench stability analysis: a Case Study in Iran | |
| Contreras | Reviewing Laubscher’s empirical method to estimate subsidence limits | |
| Majcherczyk et al. | Numerical modeling used for designing of coal mine roadway support | |
| Leroux et al. | Tunnel Database: An Information System Useful for Underground Construction in Montreal | |
| Evin et al. | Pit slope optimization based on hydrogeologic inputs | |
| Zeng et al. | Exploration and Analysis of the Law of Overburden and Surface Ecological Damage Caused by Mining |
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 |