CN105926569B - A kind of old goaf site stability quantitative evaluation method in colliery based on settlement monitoring data - Google Patents
A kind of old goaf site stability quantitative evaluation method in colliery based on settlement monitoring data Download PDFInfo
- Publication number
- CN105926569B CN105926569B CN201610272795.0A CN201610272795A CN105926569B CN 105926569 B CN105926569 B CN 105926569B CN 201610272795 A CN201610272795 A CN 201610272795A CN 105926569 B CN105926569 B CN 105926569B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msup
- mfrac
- msub
- mtr
- 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.)
- Expired - Fee Related
Links
Classifications
-
- E—FIXED CONSTRUCTIONS
- E02—HYDRAULIC ENGINEERING; FOUNDATIONS; SOIL SHIFTING
- E02D—FOUNDATIONS; EXCAVATIONS; EMBANKMENTS; UNDERGROUND OR UNDERWATER STRUCTURES
- E02D1/00—Investigation of foundation soil in situ
- E02D1/02—Investigation of foundation soil in situ before construction work
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Soil Sciences (AREA)
- Analytical Chemistry (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Paleontology (AREA)
- Civil Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Structural Engineering (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
Abstract
The invention discloses a kind of old goaf site stability quantitative evaluation method in colliery based on settlement monitoring data, and old goaf place is exploited suitable for the main caving method being over that settles.Its step is 1)In the old goaf earth's surface arrangement monitoring line in colliery, surface subsidence monitoring data are obtained;2)Earth's surface deformation behaviour initial analysis based on Monitoring Data;3)The old goaf earth's surface deformation in colliery calculates;4)Old goaf site stability evaluation.This method can utilize the probability integration process for considering time effect, when parameter determines, this relational expression of the decrement of residual settlement of the incrementss=be calculated settled according to actual monitoring, appraisal thus can be carried out using the Monitoring Data of more than 2 years, it is easy to operate, be easy to carry out, and Appreciation gist professional standard, code, evaluation effect is good, the development of the old goaf site stability quantitative assessment in colliery is greatly facilitated.
Description
Technical field
The present invention relates to a kind of goaf site stability evaluation technical method, more particularly to it is a kind of based on Monitoring Data
Old goaf stability quantitative evaluation method, old goaf place is exploited suitable for the main caving method being over that settles.
Background technology
With the development of the social economy, the utilization of land resource and the underground space is increasingly taken seriously, adopted always in the past
Land reclamation is carried out on dead zone, builds industrial premises and residential building necessitates, but the building yard stability in old goaf
Evaluation is a urgent problem to be solved.Old goaf is covered, in addition due to buried underground by unconsolidated sediment, broken sillar
Coal mining is unordered, and the spatial distribution in goaf and duration are failed to understand, it is difficult to effectively be commented goaf site stability
Valency.At present, the new technology of surface movement monitoring has obtained large development, including total powerstation, GPS global positioning systems with new instrument
System, ADAS, InSar real-time monitoring system etc., by the analysis to measured data, surface deformation law is found, and then to mined out
It is a kind of effective method that area's subsidence, which carries out evaluation,.
Due to the importance of measured data, the accumulation and analysis work of field data, China are all paid much attention in countries in the world
In the whole nation, major mineral bureau has carried out observation work since the fifties, according to incompletely statistics, the existing line of observation more than 1000 in China
Bar, the method for establishing data processing and analysis, such as Coal Mining Subsidence Estimation System.The problem of being primarily present at present has:(1) open
Adopting settlement prediction model can be predicted using Monitoring Data to the stability of rock-soil body on new working seam periphery, but can not be right
The old goaf site stability for having exploited the several years is predicted and evaluated.For example, Yin Yanbo (2008) is by mined out
Area's peripheral rock body is monitored in real time, while collecting Monitoring Data, using nonlinear science theoretical prediction model method pair
The stability of country rock body carries out prediction research.Li Leis etc. (2012) establish Wavelet Denoising Method model, Settlement Prediction Model and mined out
Area's unstability distinguishes model, and goaf unstability early warning is inquired into reference to Xi Hao village iron ore VI-4 ore bodies goaf situation.(2) utilize
InSar real-time monitoring systems, can be to regionality in large area although the history settlement monitoring data of long period can be obtained
Table settling characteristics and trend are predicted, but are not so good as traditional Monitoring method of the subsidence to the less goaf place of area, its precision
Height, and do not establish corresponding stability quantitative evaluation method.For example, Deng's noise made in coughing or vomiting medium (2015) utilizes In SAR technologies to Shan
Certain western scape Terra SAR image of mining area 19 is handled and analyzed, establish earth's surface remnants subsidence velocities circulation peak value and adopt it is thick,
Subsidence velocity cycle period sinks with deep ratio and face propulsion speed, working face accumulation and stops adopting the empirical relation of time
Formula.Liu Xiaofei etc. (2014) thinks that D-In SAR technologies can apply to the old goaf residual deformation of monitor large-area, for
The small old goaf of residual deformation can't obtain preferable result after settlement stability.(3) traditional Monitoring method of the subsidence, although
Precision is high, but is difficult to obtain prolonged settlement monitoring data in practical operation, and Monitoring Data often only has 2 years or 3 years, such as
What carries out stability quantitative assessment according to the Monitoring Data of 2 years or 3 years to old goaf, is that be currently needed for solving one is important
Problem.
In view of the above-mentioned problems, utilize the Monitoring Data of obtained short period (2 years or 3 years), with reference to being carried in invention
The probability integration process of duration, proposes that one kind is easy to operate, can carry out the old goaf in colliery after the consideration seam mining gone out
The method of site stability quantitative assessment is increasingly necessary.
The content of the invention
The technical problem to be solved in the present invention is to provide a kind of easily operated old goaf based on settlement monitoring data
Site stability quantitative evaluation method.
In order to solve the above technical problems, the technical solution used in the present invention is:A kind of coal based on settlement monitoring data
The old goaf site stability quantitative evaluation method of ore deposit, its key technology are:Comprise the steps
1) in the old goaf earth's surface arrangement monitoring line in colliery, surface subsidence monitoring data are obtained
Method for arranging is:For exploiting working face, at least 1 monitoring line is arranged by parallel bearing, vertical coal seam is walked
To at least 3 monitoring lines of arrangement;
2) the earth's surface deformation behaviour initial analysis based on Monitoring Data
Analyzed according to Monitoring Data, carry out the analysis of two aspects:
(1) the maximum lower sinker position of ground settlement is determined
The subsidence curve figure of each settlement point and settling amount is established according to the measured data of acquisition, is then analyzed parallel and vertical
The monitoring point settling amount of bearing, find out the position of sinker and numbering under maximum;
(2) judge whether main sedimentation terminates
The deflection and time curve and subsidence velocity curve of maximum lower sinker according to determined by step (1), with
Figure comprising general goaf law of subsidence curve compares, if the subsidence curve pattern that actual monitoring obtains enters
Decline phase pattern, it is possible to determine that the main sedimentation in goaf is over, if the main sedimentation in old goaf is over, you can think
The surplus value of subsidence in old goaf place is equal to residual settlement amount;Earth's surface deformation can be carried out according to this data to calculate;
3) the old goaf earth's surface deformation in colliery calculates
For gently inclined seam, using probability integration process calculation formula, sinking:
Tilt:
Curvature:
Move horizontally:
Horizontal distortion:
In formula:W (x, y) is the sinking of earth's surface arbitrfary point (coordinate is x and y), mm;
ix(x,y)、KX(x,y)、UX(x,y)、εx(x, y) is respectively the inclination in earth's surface arbitrfary point (coordinate is x and y) x directions
It is worth (mm/m), curvature value (10-3/ m), horizontal movement value (mm) and horizontal distortion value (mm/m);
iy(x,y)、Ky(x,y)、Uy(x,y)、εy(x, y) is respectively the inclination in earth's surface arbitrfary point (coordinate is x and y) y directions
It is worth (mm/m), curvature value (10-3/ m), horizontal movement value (mm) and horizontal distortion value (mm/m);W0、U0Respectively earth's surface is fully adopted
Dynamic maximum sinking value and maximum horizontal movement value (mm);
Wherein,
Limit of integration (the m that D is working face tilt length and strike length determines2);
R is the main radius of influence (m);
η is subsidence factor;M is working seam thickness (m);
α is seam inclination (°);B is displacement factor;
Wherein, remaining subsidence factor is with the expression formula of sinking time:
ηt'=Ae-Bt
(1-8),
Wherein, A and B is the parameter of determination;
By remaining subsidence factor ηtSubsidence factor η in ' alternate form (3-6), it is possible to calculate the maximum residual of different times
Remaining settling amount:
We=ηt′M cos(α) (1-9)
Again by WeW in alternate form (3-2)~(3-5)0, you can using Matlab software programmings program calculate residual tilt,
The size of curvature and horizontal distortion, then draw remaining sinking isopleth, tilting value isopleth, horizontal distortion isopleth and song
Rate value isoline figure;
4) old goaf site stability evaluation
(1) according to obtained above-mentioned data and figure, the stability in old goaf place can be carried out according to evaluation criterion
Evaluation;Site stability evaluation standard foundation《Goaf Road Design and construction technology detailed rules and regulations》(JTG/T-D31-03-2011)
Carry out;
(2) according to above-mentioned standard, by 4 the remaining sinking being calculated, tilting value, horizontal distortion value and curvature value fingers
Mark carries out global alignment, and draws comparison diagram;
(3) zone boundary determined by stability criterion is drawn according to estimation of stability standard and comparison diagram, extracts coordinate,
It is projected on CAD diagram, has just obtained old goaf site stability evaluation block plan.
Land subsidence monitoring network is made up of work benchmark and each monitoring point in the step 1), is seen with second grade leveling
Survey.
The parameter A is issued according to exploitation working face lithologic character and 2011 year traffic Department of Transportation《Goaf highway
Technique of design and construction detailed rules and regulations》After Appendix D .0.1-1 determines parameter η, determined by formula A=1- η;Parameter B is according to being calculated
Sinker is contrasted in the sinking in a certain year and the actual monitoring settling amount of the position under certain maximum, is repeated several times after calculating just
It can determine that the size of B values.
It is using beneficial effect caused by above-mentioned technical proposal:This method can utilize the probability for considering time effect
Integration method, when parameter determines, the decrement of the residual settlement of the incrementss=be calculated settled according to actual monitoring this
Relational expression, thus appraisal can be carried out using the Monitoring Data of more than 2 years, it is easy to operate, be easy to carry out, and evaluation according to
According to professional standard, code, evaluation effect is good, greatly facilitates the development of the old goaf site stability quantitative assessment in colliery.
Brief description of the drawings
A certain working face surface subsidence monitoring line arrangement schematic diagram in the old goaf places of Fig. 1;
The settling amount of each monitoring point of parallel (vertical) bearings of Fig. 2;
Fig. 3 is the determination method of surface movement perdurabgility;
Fig. 4-1-a are sinking isopleth;4-1-b is tilts isogram, and 4-1-c is horizontal distortion isogram, 4-
1-d is curvature value isogram;
Fig. 4-2 index comprehensive comparison diagrams;
Fig. 5 is the workflow diagram of the present invention;
Embodiment
The present invention is further detailed below by specific the drawings and specific embodiments.
1st, the arrangement of goaf earth's surface monitoring line
The larger old goaf of one scope is often made up of multiple working faces exploited, the recovery time of working face
Differ, thus the settling phase of different operating face overlying Rock And Soil also differs.If adjacent relatively near between working face, the later stage
The working face of exploitation also can produce material impact to the settling amount in previous work face.These influence factors are considered, for opening
Mining face, 1 monitoring line generally is arranged by parallel bearing, vertical bearing arranges 3 monitoring lines, as shown in Figure 1.
Land subsidence monitoring network is made up of work benchmark and each monitoring point, is observed with second grade leveling.
2nd, the earth's surface deformation behaviour initial analysis based on Monitoring Data
, it is necessary to make the analysis of 2 aspects after surface subsidence monitoring data are obtained:(1) the maximum lower sinker of ground settlement
The determination of position;(2) it is main to settle the judgement whether terminated.
(1) determination of the maximum lower sinker position of ground settlement
After measured data is obtained, by curve map as shown in Figure 2, the monitoring of parallel and vertical bearing is analyzed
Point settling amount, can find out the position of sinker and numbering under maximum.
(2) it is main to settle the judgement whether terminated
Referring to accompanying drawing 3, specific decision method is:The deflection of maximum lower sinker and time according to determined by above-mentioned (1)
Relation curve and subsidence velocity curve (subsidence velocity curve foundation《Goaf Road Design and construction technology detailed rules and regulations》(JTG/T-
D31-03-2011) page 68), (come from the 3D.0.3-1 comprising general goaf law of subsidence curve《Goaf is public
Road technique of design and construction detailed rules and regulations》(JTG/T-D31-03-2011)) compare, or with reference to technical manual Appendix D .0.3 other
Method, if the subsidence curve pattern that actual monitoring obtains enters decline phase pattern, it is possible to determine that the main sedimentation in goaf is
Terminate.If the main sedimentation in old goaf is over, you can thinks that the surplus value of subsidence in old goaf place is equal to remaining sink
Drop amount, suitable for old goaf site stability quantitative evaluation method provided by the present invention.
3rd, the old goaf earth's surface Method for Calculating Deformation based on Monitoring Data
In numerous prediction of mining subsidence methods, probability integration process be current China be widely used, ripe method.
What State Administration of Coal Industry in 2000 promulgated《Building, water body, railway and main roadway coal column are stayed located at mining pressed coal code》
In, probability integration process is classified as main prediction of mining subsidence method.The issue of 2011 Department of Transportation of year traffic portion《Goaf highway
Technique of design and construction detailed rules and regulations》Recommendation is used cooperatively surface movement and deformation observation method and probability integration process, to calculate the surplus of place
Remaining deformation values.
3.1 consider the probability integration process expression formula of time effects
Probability integration process predicts surface movement and deformation such as formula (1-1)~formula (1-5):(quoted from State Administration of Coal Industry《Building
Thing, water body, railway and main roadway coal column, which stay, to be set and mining pressed coal code》, Coal Industry Press, 2000.)
Sinking:
Tilt:
Curvature:
Move horizontally:
Horizontal distortion:
In formula:W (x, y) is the sinking of earth's surface arbitrfary point (coordinate is x and y), mm;
ix(x,y)、KX(x,y)、UX(x,y)、εx(x, y) is respectively the inclination in earth's surface arbitrfary point (coordinate is x and y) x directions
It is worth (mm/m), curvature value (10-3/ m), horizontal movement value (mm) and horizontal distortion value (mm/m);iy(x,y)、Ky(x,y)、Uy(x,
y)、εy(x, y) is respectively tilting value (mm/m), the curvature value (10 in earth's surface arbitrfary point (coordinate is x and y) y directions-3/ m), it is horizontal
Movement value (mm) and horizontal distortion value (mm/m);W0、U0The respectively maximum sinking value of earth's surface sufficient mining and maximum horizontal movement
It is worth (mm);
Limit of integration (the m that D is working face tilt length and strike length determines2);R is the main radius of influence (m);η is
Subsidence factor;M is working seam thickness (m);α is seam inclination (°);B is displacement factor.
When sinking time influences after considering to exploit, due to formed after exploitation underground cavity, crack etc. with sedimentation progressively
Densification, original lithologic character change, and subsidence factor should change with time and stress state, thus it needs to be determined that sinking
Coefficient changes with time relation by for predicting that the residual settlement amount in old goaf is most important.Residual settlement amount refers under earth's surface
Sink after terminating substantially, crush the deflection of slight void earth's surface after the very long time is further closely knit between sillar.
It is assumed that under load action or under slacking, space, crack and absciss layer gap close completely existing for broken rock
Close, subsidence factor reaches 1, then the calculation expression (1-7) of maximum remaining subsidence factor is:
η '=1- η (1-7)
Hard rock in being belonged to due to mining area lithology, adopting thick larger and bottom has more tomography, according to《Goaf Road Design
With construction technology detailed rules and regulations》, to adopting the goaf that depth is 200~300m, the time limit of sinking is about 1 year, and main sedimentation is basically completed.At any time
Between increase stress equilibrium destroy, rock stratum activation, continue to sink, subsidence factor now progressively subtracts with time and densification degree
Small, Guo Guangli etc. combines Xuzhou mine of western region data and analyzed, under what different times goaf was taken has completed
Heavy coefficient is shown in Table 1, and model flood winter etc. (2009) gives remaining subsidence factor and is with the expression formula (1-8) of sinking time:
ηt'=Ae-Bt (1-8)
In formula, A, B are respectively parameter;ηt' it is remaining subsidence factor;T is the time (a).The formula is negative exponential function, when residual
During remaining heavy just beginning, sinking is the remaining subsidence factor of maximum, and sinking when being t=0 is A=1- η.When parameter B is true
After fixed, it is possible to which the residual settlement coefficient of different times is calculated.
The subsidence factor reference value of the Site in Xuzhou Mining Area region different times of table 1
By remaining subsidence factor ηtSubsidence factor η in ' alternate form (1-6), it is possible to calculate the maximum residual of different times
Remaining settling amount, expression formula (1-9):
We=ηt′Mcos(α) (1-9)
Again by WeW in alternate form (1-2)~(1-5)0, you can the size of residual tilt, curvature and horizontal distortion is calculated,
And then the stability in old goaf place can be evaluated according to evaluation criterion.3.2 earth's surfaces deform calculating parameter and condition
It is determined that
(1) evaluation work region feature and numbering
Full basin is mobile and deformation calculation formula is only applicable to the mined out region of computation rule rectangle, it is contemplated that mining area exists more
Individual goaf, and gob edge shape is more complicated, and goaf can be simplified to multiple rectangle evaluation work faces.Each evaluation work
The essential characteristic indicator-specific statistics in face, filled in by the statistics of table 2.
The working face essential characteristic of table 2 and indicator-specific statistics table
(2) remaining subsidence factor ηt' determination
By formula (1-8) as long as determining parameter A and B, remaining subsidence factor η just can determine thatt′。
Parameter A determination:Parameter A is according to exploitation working faceLithologic characterIssued with 2011 year traffic Department of Transportation《Adopt
Dead zone Road Design and construction technology detailed rules and regulations》Appendix D .0.1-1After determining parameter η, determined by formula A=1- η (from formula 1- above
Discussion in 8), wherein η is the subsidence factor q in Appendix D .0.1-1.Appendix D .0.1-1 is shown in Table lattice 3.
Parameter B determination:By research and application data, determine that setting basin has completed main sedimentation, in residual settlement
In the stage, then the maximum point of settling amount in selective analysis Monitoring Data, finds out position of this in setting basin, provides one
B values carry out tentative calculation, and calculating utilizes formula (1-1) and (1-9), and by WeW in alternate form (1-1)0, you can it is maximum to calculate this
Sinking of the lower sinker in a certain year.Then result of calculation and the actual monitoring settling amount of the position are contrasted, it is repeatedly several
The size of B values is just can determine that after secondary calculating.Institute according to principle be:Obtain incrementss=calculating residual settlement of monitoring sedimentation
Decrement, the time typically calculated in units of year.Such as
W11 years-W10 years=S10 years-S11 years (1-10)
In formula, W10 years、W11 yearsRespectively the monitoring deflection of the 10th year and the monitoring deflection of the 11st year, S10 years、S11 yearsRespectively
Remaining deflection and the remaining deflection of the 11st year for the 10th year.
The Appendix D .0.1-1 of table 3
It is the parameter A and parameter B that can determine that each working face in goaf by the above method, may be referred to table 4 and list it
It is worth size.
Table 4 parameter A and B value
(3) main moving influence angle beta
The main angle tangent value that influences is the ratio between mining depth and the main radius of influence.According to the lithologic character in goaf, and
Reference《Goaf Road Design and construction technology detailed rules and regulations》Appendix D .0.1-1 carries out value.
(4) displacement factor b
Displacement factor reflects the proportionate relationship between earth's surface maximum horizontal movement value and maximum sinking value.Can be according to adopting
The lithologic character of dead zone, and reference《Goaf Road Design and construction technology detailed rules and regulations》Appendix D .0.1-1 carries out value, generally takes
B=0.25.
4th, old goaf site stability evaluation
4.1 residual deformation quantitative calculation methods
On Matlab software platforms, by working out simple small routine, you can solve formula (1-1)~(1~9) and count
Calculate, and then obtain remaining sinking of any location point after a certain year in each working face in goaf, residual tilt value, residual
Remaining curvature value, residual level deformation values.For the larger neighbouring exploitation working face that influences each other, above-mentioned each value needs are folded
Add.Acquired results can be counted by the form of table 5.
The residual deformation result of calculation of the ××× working face of table 5
Then the data obtained using surfer softwares and table 5, are depicted as remaining sinking in the range of value of mining area
Isopleth, tilting value isopleth, horizontal distortion isopleth and curvature value isoline figure.With reference to example such as Fig. 4-1-a, 4-1-b, 4-
Shown in 1-c, 4-1-d.
4.2 old goaf site stability grade evaluation criterion, principle and methods
Site stability evaluation standard foundation《Goaf Road Design and construction technology detailed rules and regulations》(JTG/T-D31-03-
2011), quantitative assessing index is as shown in table 6, and wherein surface movement and deformation value is the remaining shifting of earth's surface after building levelling of the land
Dynamic deformation value.
Table 6 presses the site stability grade evaluation criterion that surface movement and deformation value determines
The principle of evaluation is:A bit in area is evaluated, as long as having 1 to reach unstable in 4 indexs, that is, is thought at this not
Stable, according to the worst one-level of stability, division determines stability grade.
According to mentioned above principle, by 4 the remaining sinking being calculated, tilting value, horizontal distortion value and curvature value indexs
Global alignment is carried out, example is as shown in the Fig. 4-2.Solid line is remaining sinking in figure, real thick line is horizontal distortion value, empty thick line is
Curvature value.Because tilting value is both less than Stability Assessment standard 3, do not contrasted on the diagram.
With reference to the stability grade evaluation criterion in table 6, Fig. 4-2 is analyzed, draws and is determined by stability criterion
Zone boundary, extract coordinate, it is projected on CAD diagram, has just obtained old goaf site stability evaluation block plan.
The evaluation method is successfully applied to Xingtai Mining Area by us.Xingtai Mining Area is located at Hebei South Xingtai urban district and city
Area is western, and landform Xi Gaodong is low, and western part is Taihang Mountain mountains low-to-middle in height landforms;It is hills landforms to the west of the Beijing-Guangzhou Railway of east, Yi Dongwei
Plain landforms.Upper group of coal include 1,2,4 under, No. 5 coal seams.The group 1,2 coal seams top, bottom plate are essentially sandstone aquifer, under 4,5#
Coal seam top, bottom plate are wild blue or green, the steel gray rock of volt.Recovery method is moves towards longwell, retrogressing in area, descending across falling coal-mining method.Xingtai coal
Ore deposit only adopts No. 2 coals inside planning region and nearby.To grasp coal mine gob deformation characteristicses, deformation cycle and deflection, for prison
Survey area's Evaluation of The Ground Stability and foundation is provided.Carry out monitoring to mining area, the monitoring phase is 2 years, is monitored five times altogether, from
In October, 2012 in October, 2014, it is in October, 2012, in April, 2013, in October, 2013, in April, 2014 and 2014 respectively
October in year.There are 8 working faces that have of Monitoring Data covering at present, institute's mining coal seam is No. 2 coal seams, and working face numbering is respectively
7816、7802、7820、7804、7806、7808、7810、7812。
Using such a evaluation method, obtained the remaining sinking in Xingtai Mining Area goaf place, tilting value, curvature value and
The distribution of contours figure of horizontal distortion value, and then establishing criteria has carried out estimation of stability to old goaf, meets industry standard
It is required that obtain the accreditation of expert.Moreover, observed in the recent period the place and surrounding building, the stable region that is divided,
Basically stable area and understable area are consistent with observation result.
Claims (3)
- A kind of 1. old goaf site stability quantitative evaluation method in colliery based on settlement monitoring data, it is characterised in that:Bag Include following step1) in the old goaf earth's surface arrangement monitoring line in colliery, surface subsidence monitoring data are obtainedMethod for arranging is:For exploiting working face, at least 1 monitoring line, vertical bearing cloth are arranged by parallel bearing Put at least 3 monitoring lines;2) the earth's surface deformation behaviour initial analysis based on Monitoring DataAnalyzed according to Monitoring Data, carry out the analysis of two aspects:(1) the maximum lower sinker position of ground settlement is determinedThe subsidence curve figure of each settlement point and settling amount is established according to the measured data of acquisition, then analyzes parallel and vertical coal seam The monitoring point settling amount of trend, find out the position of sinker and numbering under maximum;(2) judge whether main sedimentation terminatesThe deflection of maximum lower sinker and time curve and subsidence velocity curve according to determined by step (1), with comprising The figure of general goaf law of subsidence curve compares, if the subsidence curve pattern that actual monitoring obtains enters decline Phase pattern, it is possible to determine that the main sedimentation in goaf is over, if the main sedimentation in old goaf is over, you can think to adopt always The surplus value of subsidence in dead zone place is equal to residual settlement amount;Earth's surface deformation can be carried out according to this data to calculate;3) the old goaf earth's surface deformation in colliery calculatesFor gently inclined seam, using probability integration process calculation formula, sinking:<mrow> <mi>W</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>W</mi> <mn>0</mn> </msub> <munder> <mrow> <mo>&Integral;</mo> <mo>&Integral;</mo> </mrow> <mi>D</mi> </munder> <mfrac> <mn>1</mn> <msup> <mi>r</mi> <mn>2</mn> </msup> </mfrac> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>&pi;</mi> <mfrac> <mrow> <msup> <mrow> <mo>(</mo> <mi>s</mi> <mo>-</mo> <mi>x</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mi>y</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> <msup> <mi>r</mi> <mn>2</mn> </msup> </mfrac> <mo>)</mo> </mrow> <mi>d</mi> <mi>s</mi> <mi>d</mi> <mi>t</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>Tilt:<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>i</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>W</mi> <mn>0</mn> </msub> <munder> <mrow> <mo>&Integral;</mo> <mo>&Integral;</mo> </mrow> <mi>D</mi> </munder> <mfrac> <mrow> <mn>2</mn> <mi>&pi;</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>-</mo> <mi>x</mi> <mo>)</mo> </mrow> </mrow> <msup> <mi>r</mi> <mn>4</mn> </msup> </mfrac> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>&pi;</mi> <mfrac> <mrow> <msup> <mrow> <mo>(</mo> <mi>s</mi> <mo>-</mo> <mi>x</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mi>y</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> <msup> <mi>r</mi> <mn>2</mn> </msup> </mfrac> <mo>)</mo> </mrow> <mi>d</mi> <mi>s</mi> <mi>d</mi> <mi>t</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>i</mi> <mi>y</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>W</mi> <mn>0</mn> </msub> <munder> <mrow> <mo>&Integral;</mo> <mo>&Integral;</mo> </mrow> <mi>D</mi> </munder> <mfrac> <mrow> <mn>2</mn> <mi>&pi;</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mi>y</mi> <mo>)</mo> </mrow> </mrow> <msup> <mi>r</mi> <mn>4</mn> </msup> </mfrac> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>&pi;</mi> <mfrac> <mrow> <msup> <mrow> <mo>(</mo> <mi>s</mi> <mo>-</mo> <mi>x</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mi>y</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> <msup> <mi>r</mi> <mn>2</mn> </msup> </mfrac> <mo>)</mo> </mrow> <mi>d</mi> <mi>s</mi> <mi>d</mi> <mi>t</mi> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>Curvature:<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>K</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>W</mi> <mn>0</mn> </msub> <munder> <mrow> <mo>&Integral;</mo> <mo>&Integral;</mo> </mrow> <mi>D</mi> </munder> <mfrac> <mrow> <mn>2</mn> <mi>&pi;</mi> </mrow> <msup> <mi>r</mi> <mn>4</mn> </msup> </mfrac> <mrow> <mo>(</mo> <mfrac> <mrow> <mn>2</mn> <mi>&pi;</mi> <msup> <mrow> <mo>(</mo> <mi>s</mi> <mo>-</mo> <mi>x</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> <msup> <mi>r</mi> <mn>2</mn> </msup> </mfrac> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>&pi;</mi> <mfrac> <mrow> <msup> <mrow> <mo>(</mo> <mi>s</mi> <mo>-</mo> <mi>x</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mi>y</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> <msup> <mi>r</mi> <mn>2</mn> </msup> </mfrac> <mo>)</mo> </mrow> <mi>d</mi> <mi>s</mi> <mi>d</mi> <mi>t</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>K</mi> <mi>y</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>W</mi> <mn>0</mn> </msub> <munder> <mrow> <mo>&Integral;</mo> <mo>&Integral;</mo> </mrow> <mi>D</mi> </munder> <mfrac> <mrow> <mn>2</mn> <mi>&pi;</mi> </mrow> <msup> <mi>r</mi> <mn>4</mn> </msup> </mfrac> <mrow> <mo>(</mo> <mfrac> <mrow> <mn>2</mn> <mi>&pi;</mi> <msup> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mi>y</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> <msup> <mi>r</mi> <mn>2</mn> </msup> </mfrac> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>&pi;</mi> <mfrac> <mrow> <msup> <mrow> <mo>(</mo> <mi>s</mi> <mo>-</mo> <mi>x</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mi>y</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> <msup> <mi>r</mi> <mn>2</mn> </msup> </mfrac> <mo>)</mo> </mrow> <mi>d</mi> <mi>s</mi> <mi>d</mi> <mi>t</mi> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>Move horizontally:<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>U</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>U</mi> <mn>0</mn> </msub> <munder> <mrow> <mo>&Integral;</mo> <mo>&Integral;</mo> </mrow> <mi>D</mi> </munder> <mfrac> <mrow> <mn>2</mn> <mi>&pi;</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>-</mo> <mi>x</mi> <mo>)</mo> </mrow> </mrow> <msup> <mi>r</mi> <mn>3</mn> </msup> </mfrac> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>&pi;</mi> <mfrac> <mrow> <msup> <mrow> <mo>(</mo> <mi>s</mi> <mo>-</mo> <mi>x</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mi>y</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> <msup> <mi>r</mi> <mn>2</mn> </msup> </mfrac> <mo>)</mo> </mrow> <mi>d</mi> <mi>s</mi> <mi>d</mi> <mi>t</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>U</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>U</mi> <mn>0</mn> </msub> <munder> <mrow> <mo>&Integral;</mo> <mo>&Integral;</mo> </mrow> <mi>D</mi> </munder> <mfrac> <mrow> <mn>2</mn> <mi>&pi;</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>-</mo> <mi>x</mi> <mo>)</mo> </mrow> </mrow> <msup> <mi>r</mi> <mn>3</mn> </msup> </mfrac> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>&pi;</mi> <mfrac> <mrow> <msup> <mrow> <mo>(</mo> <mi>s</mi> <mo>-</mo> <mi>x</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mi>y</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> <msup> <mi>r</mi> <mn>2</mn> </msup> </mfrac> <mo>)</mo> </mrow> <mi>d</mi> <mi>s</mi> <mi>d</mi> <mi>t</mi> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>Horizontal distortion:<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>&epsiv;</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>U</mi> <mn>0</mn> </msub> <munder> <mrow> <mo>&Integral;</mo> <mo>&Integral;</mo> </mrow> <mi>D</mi> </munder> <mfrac> <mrow> <mn>2</mn> <mi>&pi;</mi> </mrow> <msup> <mi>r</mi> <mn>3</mn> </msup> </mfrac> <mrow> <mo>(</mo> <mfrac> <mrow> <mn>2</mn> <mi>&pi;</mi> <msup> <mrow> <mo>(</mo> <mi>s</mi> <mo>-</mo> <mi>x</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> <msup> <mi>r</mi> <mn>2</mn> </msup> </mfrac> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>&pi;</mi> <mfrac> <mrow> <msup> <mrow> <mo>(</mo> <mi>s</mi> <mo>-</mo> <mi>x</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mi>y</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> <msup> <mi>r</mi> <mn>2</mn> </msup> </mfrac> <mo>)</mo> </mrow> <mi>d</mi> <mi>s</mi> <mi>d</mi> <mi>t</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>&epsiv;</mi> <mi>y</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>U</mi> <mn>0</mn> </msub> <munder> <mrow> <mo>&Integral;</mo> <mo>&Integral;</mo> </mrow> <mi>D</mi> </munder> <mfrac> <mrow> <mn>2</mn> <mi>&pi;</mi> </mrow> <msup> <mi>r</mi> <mn>3</mn> </msup> </mfrac> <mrow> <mo>(</mo> <mfrac> <mrow> <mn>2</mn> <mi>&pi;</mi> <msup> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mi>y</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> <msup> <mi>r</mi> <mn>2</mn> </msup> </mfrac> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>&pi;</mi> <mfrac> <mrow> <msup> <mrow> <mo>(</mo> <mi>s</mi> <mo>-</mo> <mi>x</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mi>y</mi> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> <msup> <mi>r</mi> <mn>2</mn> </msup> </mfrac> <mo>)</mo> </mrow> <mi>d</mi> <mi>s</mi> <mi>d</mi> <mi>t</mi> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>In formula:W (x, y) is the sinking of earth's surface arbitrfary point, and coordinate is x and y, unit mm;ix(x,y)、KX(x,y)、UX(x,y)、εx(x, y) is respectively the tilting value in earth's surface arbitrfary point x directions, wherein earth's surface arbitrfary point Coordinate is x and y, unit mm/m;Curvature value, unit 10-3/m;Horizontal movement value, unit mm;With horizontal distortion value, unit mm/m;iy(x,y)、Ky(x,y)、Uy(x,y)、εy(x, y) is respectively the tilting value in earth's surface arbitrfary point y directions, the seat of earth's surface arbitrfary point It is designated as x and y, unit mm/m;Curvature value, unit 10-3/m;Horizontal movement value, unit mm;With horizontal distortion value, unit mm/m; W0、U0The respectively maximum sinking value of earth's surface sufficient mining and maximum horizontal movement value, unit mm;Wherein,<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>W</mi> <mn>0</mn> </msub> <mo>=</mo> <mi>&eta;</mi> <mi>M</mi> <mi> </mi> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mrow> <mo>(</mo> <mi>&alpha;</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>U</mi> <mn>0</mn> </msub> <mo>=</mo> <msub> <mi>bW</mi> <mn>0</mn> </msub> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>The limit of integration that D is working face tilt length and strike length determines, m2;R is the main radius of influence, m;η is subsidence factor;M is working seam thickness, m;α is seam inclination, unit degree;B is displacement factor;It is assumed that under load action or under slacking, space, crack and absciss layer gap close completely existing for broken rock, under Heavy coefficient reaches 1, then the calculation expression of maximum remaining subsidence factor is:η '=1- η (1-7)Remaining subsidence factor is with the expression formula of sinking time:η′t=Ae-Bt(1-8),Wherein, A and B is the parameter of determination;By remaining subsidence factor η 'tSubsidence factor η in alternate form (1-6), it is possible to calculate the maximum residual settlement of different times Amount:We=η 'tMcos(α) (1-9)Again by WeW in alternate form (1-2)~(1-3)0, use WeB replaces the U in (1-4)~(1-5)0, you can it is soft using Matlab Part, which programs, calculates the size of residual tilt, curvature and horizontal distortion, then draws remaining sinking isopleth, tilting value etc. It is worth line, horizontal distortion isopleth and curvature value isogram;4) old goaf site stability evaluation(1) according to obtained above-mentioned data and figure, the stability in old goaf place can be evaluated according to evaluation criterion; Site stability evaluation standard foundation《Goaf Road Design and construction technology detailed rules and regulations》Carry out, JTG/T-D31-03-2011;(2) according to above-mentioned standard, 4 the remaining sinking being calculated, tilting value, horizontal distortion value and curvature value indexs are entered Row global alignment, and draw comparison diagram;(3) zone boundary determined by stability criterion is drawn according to estimation of stability standard and comparison diagram, coordinate is extracted, by it Project on CAD diagram, just obtained old goaf site stability evaluation block plan.
- 2. the old goaf site stability quantitative assessment side in the colliery according to claim 1 based on settlement monitoring data Method, it is characterised in that:Land subsidence monitoring network is made up of work benchmark and each monitoring point in the step 1), with second-class water Quasi- observation.
- 3. the old goaf site stability quantitative assessment side in the colliery according to claim 1 based on settlement monitoring data Method, it is characterised in that:The parameter A is issued according to exploitation working face lithologic character and 2011 year traffic Department of Transportation《It is mined out Area's Road Design and construction technology detailed rules and regulations》After Appendix D .0.1-1 determines parameter η, determined by formula A=1- η;Parameter B is according to calculating Sinker is contrasted in the sinking in a certain year and the actual monitoring settling amount of the position under certain obtained maximum, is repeated several times meter The size of B values is just can determine that after calculation.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610272795.0A CN105926569B (en) | 2016-04-28 | 2016-04-28 | A kind of old goaf site stability quantitative evaluation method in colliery based on settlement monitoring data |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610272795.0A CN105926569B (en) | 2016-04-28 | 2016-04-28 | A kind of old goaf site stability quantitative evaluation method in colliery based on settlement monitoring data |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105926569A CN105926569A (en) | 2016-09-07 |
CN105926569B true CN105926569B (en) | 2018-02-16 |
Family
ID=56836543
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610272795.0A Expired - Fee Related CN105926569B (en) | 2016-04-28 | 2016-04-28 | A kind of old goaf site stability quantitative evaluation method in colliery based on settlement monitoring data |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105926569B (en) |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106407569B (en) * | 2016-09-22 | 2019-10-08 | 河南理工大学 | Surface subsidence value calculating method under the conditions of a kind of Moderate and Thick Unconsolidated Formation thin bedrock |
CN107067135A (en) * | 2016-12-22 | 2017-08-18 | 煤炭科学技术研究院有限公司 | Sinking land in coalmining areas Land Surface Stability appraisal procedure for building photovoltaic plant |
CN107092740B (en) * | 2017-04-12 | 2020-08-04 | 乌鲁木齐城市轨道集团有限公司 | Prediction method for residual deformation of old goaf of steep-dip extra-thick coal seam |
CN108399497B (en) * | 2018-02-11 | 2021-06-22 | 通用技术集团工程设计有限公司 | Strip mining goaf stability evaluation method |
CN108595834B (en) * | 2018-04-24 | 2022-02-01 | 山东科技大学 | Coal seam top and bottom plate power partition evaluation method based on multiple geological factors |
CN108647181B (en) * | 2018-05-10 | 2021-07-02 | 中国矿业大学(北京) | Long-wall mining earth surface residual subsidence coefficient calculation method |
CN110306522B (en) * | 2019-06-14 | 2021-03-09 | 华侨大学 | Method for predicting subsidence depth of explosive compaction filling stone |
CN111366129A (en) * | 2020-02-25 | 2020-07-03 | 神华神东煤炭集团有限责任公司 | Key layer breaking area determining method |
CN111323776B (en) * | 2020-02-27 | 2021-04-13 | 长沙理工大学 | Method for monitoring deformation of mining area |
CN111602546B (en) * | 2020-05-22 | 2022-04-15 | 河北地质大学 | Desertification soil vegetation protector |
CN111680261A (en) * | 2020-06-10 | 2020-09-18 | 淮南矿业(集团)有限责任公司 | Setting method and setting system for protecting coal pillar |
CN112085248A (en) * | 2020-07-22 | 2020-12-15 | 何满潮 | Site planning design method influenced by mining subsidence |
CN112749477B (en) * | 2020-11-30 | 2021-09-24 | 安徽理工大学 | Prediction method for deflection amount of vertical shaft of coal mine with thick surface soil and thin bedrock |
CN114755681A (en) * | 2021-01-08 | 2022-07-15 | 中国地质大学(北京) | Goaf surface residual deformation prediction method based on InSAR multi-point joint solution |
CN113338929B (en) * | 2021-07-07 | 2022-08-30 | 中国矿业大学 | Prediction method for surface residual deformation of single narrow working face extremely-insufficient mining |
CN114202143A (en) * | 2021-08-25 | 2022-03-18 | 中国建筑股份有限公司 | Goaf safety evaluation method and device and storage medium |
CN116485218B (en) * | 2022-12-23 | 2024-06-25 | 中煤科工开采研究院有限公司 | Stability area dividing method and treatment time planning method for subsidence area |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7966165B2 (en) * | 2005-10-17 | 2011-06-21 | National University Corporation Nagoya University | Soil-water coupled analyzer and soil-water coupled analysis method |
CN103606019B (en) * | 2013-12-04 | 2016-09-28 | 江西理工大学 | Mine goaf overlying stratum sedimentation dynamic prediction method based on time-space relationship |
CN104408277B (en) * | 2014-09-28 | 2017-05-03 | 沈阳大学 | Method for predicting, preventing and controlling earth surface residual movement and deformation caused by newly-built building in mine lot |
CN104504461A (en) * | 2014-12-09 | 2015-04-08 | 鞍钢集团矿业公司 | Method for predicating mine pit slope deformation destroy induced by conversion from surface mining to underground mining |
-
2016
- 2016-04-28 CN CN201610272795.0A patent/CN105926569B/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN105926569A (en) | 2016-09-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105926569B (en) | A kind of old goaf site stability quantitative evaluation method in colliery based on settlement monitoring data | |
CN111750822B (en) | Coal mining induced overlying strata and surface subsidence collaborative dynamic prediction method | |
CN106326528B (en) | A kind of opencut end side underground mining induces Distribution Law of Ground Crack prediction technique | |
CN109577982A (en) | Wall type Lian Cailian fills water protection mining and water resource migration monitoring, water damage method for early warning | |
CN103902780B (en) | Solid filling coal mining earth's surface Deformation prediction method | |
CN106339528A (en) | Method for forecasting range of surface movement induced by underground mining of open-pit iron mine end slope | |
CN106150546B (en) | A kind of method that the filling of goaf part is carried out according to earth's surface subsidence | |
CN1873674B (en) | Estimating mining subsidence under turnout cluster of dedicated railroad, and governing expert system | |
CN106894817A (en) | A kind of mechanization Upward slicing wedging mixed filling mining codes | |
CN106055892A (en) | Method for determining exploitation feasibility of middle-layer abandoned coal in composite pillar type ascending mining area | |
WO2019149146A1 (en) | Grading method of mines/mining areas for water-preserved mining | |
Li et al. | Fuzzy probability measures (FPM) based non-symmetric membership function: Engineering examples of ground subsidence due to underground mining | |
CN109063285B (en) | Design method of earth slope drilling arrangement scheme | |
Kılıç et al. | Microzonation of Zeytinburnu region with respect to soil amplification: a case study | |
CN110555559A (en) | Ground settlement prediction method under insufficient mining condition of thick unconsolidated formation | |
CN106251053A (en) | The suitability assessment method of high voltage power transmission transmission tower is built on earth's surface, a kind of sinking land in coalmining areas | |
Guo et al. | Optimization of Land Saving and Loss Reducing and Slope Stability Variation Patterns in Open‐Pit Mine | |
Meshcheryakov et al. | Reclamation technology of land destroyed by mining and logistics monitoring criteria | |
CN102678124A (en) | Method for back-filling surface subsidence area based on underground metal mineral deposit mining | |
CN111400659A (en) | Method for calculating farmland protection rate of mining-selecting-charging integrated mining technology | |
Singh et al. | Slope assessment systems: A review and evaluation of current techniques used for cut slopes in the mountainous terrain of West Malaysia | |
Issabek et al. | TO THE PROBLEM OF SELECTING VERTICAL SHAFTS LOCATION WITH COMBINED GEOTECHNOLOGY OF DEVELOPING DEPOSITS. | |
CN101561359B (en) | Grading method for gas tunnel | |
CN114140027A (en) | Evaluation method for closing coal mine residual coal resource reutilization | |
Tizzano et al. | Characterization of failure parameters and preliminary slope stability analysis of the Cedar Canyon Landslide, Iron County, Utah |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180216 Termination date: 20210428 |
|
CF01 | Termination of patent right due to non-payment of annual fee |