CN106021762A - Numerical simulation method for coal mine abscission layer water periodical water inrush disaster - Google Patents
Numerical simulation method for coal mine abscission layer water periodical water inrush disaster Download PDFInfo
- Publication number
- CN106021762A CN106021762A CN201610368278.3A CN201610368278A CN106021762A CN 106021762 A CN106021762 A CN 106021762A CN 201610368278 A CN201610368278 A CN 201610368278A CN 106021762 A CN106021762 A CN 106021762A
- Authority
- CN
- China
- Prior art keywords
- water
- layer
- height
- model
- excavation
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a numerical simulation method for a coal mine abscission layer water periodical water inrush disaster. The method comprises the following steps that a numerical calculation model is established, the whole model is composed of 12 coal rock layers, the fourth layer above a second coal bed is a main sandstone aquifer, and the third mudstone layer is a key waterproof layer; when the key waterproof layer is 20 m, the length of the model in the moving direction is 600 m, the height is 300 m, the model is divided into 65,750 250*263 units, and rock mass only bears gravity stress ad water pressure; according to boundary conditions, the two ends are horizontally restrained, the bottom end is fixed, and the periphery is arranged to be a waterproof boundary. The development process of a water guide fracture is simulated through step-by-step excavation; excavation is started from the position 100 m away from the left side of the moving direction of the model, 400 m is excavated in total with the height of 5 m, 10 m is excavated in each step, and totally 40 steps are included. A computer numerical simulation means is adopted to simulate an abscission layer formation process, and periodicity of the secondary roof abscission layer water disaster is summarized.
Description
Technical field
The present invention relates to roof absciss layer water damage prevention and control field, be specifically related to a kind of colliery absciss layer water cycle gushing water and cause the method for numerical simulation of calamity.
Background technology
China's coal-mine water damage takes place frequently, especially under conditions of there is no relatively strong aquifer occur top board burst water cause calamity, if Shenhua Ning Mei group red building colliery 2 coal mining water source of inundation in mine is essentially from gritstone crack Porosity Aquifer bottom its top board Jurassic system straight sieve group, this water-bearing layer watery is weak.1121 work surfaces start back production to 2010 in JIUYUE, 2009 and advance 186m altogether March, but experienced by 4 times and concentrate retention of excessive fluid, and maximum flooding quantity is up to 3000m3/h.Tongchuan, Shaanxi Province Jiaoping coalfield jade China colliery 4-2Coal seam water filling source is essentially from Luo He group hypomere sandstone crack phreatic aquifer of uniting under top board Cretaceous System, and this water-bearing layer is aquitard.1412 work surfaces occur 6 water inrush accidents, instantaneous gushing water maximum to reach 2000m in March, 2012 in May, 20133/ h, has had a strong impact on working face extraction speed.There is roof water inrush on May 21st, 2005 in Huaibei mining industry group Hai Zi colliery, maximum stream flow reaches 3887m3/ h, instantaneous floods work surface, machine lane and air way, causes 5 people dead.Similar accident is refined long flint nozzle and ore deposit, Giant Buddha Temple, forever colliery, Gansu Province Cui Mu, Guo Jia river in Shaanxi；Hua Feng colliery, Shandong；Nan Tong colliery, Sichuan, Yutianbao coalmine；An ore deposit, Xin Ji colliery, Anhui, Laohutai colliery, Fushun etc. are got through in Chongqing all generation.This type of disaster has that harm is big, omen is inconspicuous, the gushing water moment water yield is big but total amount is less, gushing water sign is inconspicuous and there is periodic roof water inrush feature, attracting wide attention and pay attention to, Wang Jingming, Li Wen equality is defined as roof secondary absciss layer water damage through further investigation.Build congruence and Forming Mechanism and the law of development of seam mining absciss layer have been done correlational study.Qiao Wei decile is isolated the primary condition that layer water is formed, and carries out differentiating to the position that may produce absciss layer and calculates, and have studied the relatedness of mining height and the relation of absciss layer Shuitu water and face propulsion speed and the absciss layer space product water yield.Zhu soldier[9]Deng the technique study Hai Zi colliery 745 water gush face mechanism using engineering sounding and theory analysis, it is believed that this gushing water is to destroy key strata of covering rock structure under the load transmission effect of absciss layer district hydrops and cause.By there is the water charge condition analysis of especially big roof water inrush work surface in Cheng Xinming, Zhao Tuanzhi etc., it is believed that work surface adopts overlaying influence, and the separated strata seeper formed in top board thick-layer sandstone, is the main water filling source of gushing water to Hai Zi ore deposit.Li Zhongkai etc. are by hydrogeological engineering geological condition analysis, it is believed that the water source of Hai Zi ore deposit 5.21 water inrush accident, from roof sandstone separated strata seeper, is under complicated hydrogeological and engineering geological conditions, a kind of dynamically water bursting sources that coal-face top board produces.Wang Zhengming etc. utilize similar material model and numerical model to simulate Xin Ji mining area 1307 work surface seam mining overlying strata separation layer phenomenon, it is determined that top board water damage is caused by secondary absciss layer water bags burst.At present, achievement in research focuses primarily upon the aspects such as the engineer applied of absciss layer formation mechenism, the regularity of distribution and absciss layer separation in reducing subsidence by grouting, but less for researchs such as absciss layer water formation condition and water fathering mechanism thereof, and the particularly periodic rule of absciss layer hydrodynamic(al) state is rarely reported.
Summary of the invention
For solving the problems referred to above, the invention provides a kind of colliery absciss layer water cycle gushing water and cause the method for numerical simulation of calamity, simulate absciss layer forming process, summarize the periodicity of top board secondary absciss layer water damage.
For achieving the above object, the technical scheme that the present invention takes is:
Absciss layer water cycle gushing water in colliery causes the method for numerical simulation of calamity, comprises the steps:
Setting up mathematical calculation model, whole model is made up of 12 layers of coal rock layer, and wherein above No. 2 coal seams, the 4th floor is main sandstone aquifer, and the 3rd mfs layer is water-resisting key strata；When this layer of water-resisting key strata is 20m, it is 600m that model prolongs strike length, a height of 300m, is divided into 250 × 263 totally 65750 unit, and rock mass is solely subjected to weight stress and water pressure；Boundary condition is both ends horizontal constraint, and bottom is fixed, and sets periphery as water proof border.The process edge that water producing fractures is grown is simulated by step excavation；Model is walked 100m to the left and is started excavation, advances 400m, mining height 5m altogether, and often step excavation 10m, is divided into 40 steps
The method have the advantages that
Use Computer Numerical Simulation means, simulate absciss layer forming process, summarize the periodicity of top board secondary absciss layer water damage.Research shows: absciss layer is owing to the overlying strata differences of mechanical properties of main mining coal seam is relatively big, causes overlying strata to produce inharmonious consistent bending subsidence when seam mining and produces.Along with the propelling of work surface, absciss layer forms rear continued development to extreme value, when there is water filling source top, then forms absciss layer water.Now, if the convex development that jumps of water flowing fractured zone development height turns on absciss layer water body, this absciss layer water causes absciss layer water damage by bursting into down-hole.Then absciss layer is gradually compacted atrophy, and water body runs off totally.Along with roof is the most disrumpent feelings, absciss layer water damage periodically occurs.The scale that absciss layer cycle dynamics occurs is different, and wherein with the absciss layer Space Scale that occurs for the first time as maximum, the corresponding secondary water damage of absciss layer is the most maximum.Study area work surface often advances 100m, absciss layer water damage dynamically to occur 1 time.
Accompanying drawing explanation
Fig. 1 is the RFPA mathematical calculation model in the embodiment of the present invention.
Fig. 2 is springform spirogram during model excavation during advance of the face 10m in the embodiment of the present invention
Fig. 3 is the springform spirogram during model excavation of advance of the face 30m in the embodiment of the present invention.
Fig. 4 is the springform spirogram during model excavation of advance of the face 60m in the embodiment of the present invention.
Fig. 5 is the springform spirogram during model excavation of advance of the face 90m in the embodiment of the present invention.
Fig. 6 is the springform spirogram during model excavation of advance of the face 120m in the embodiment of the present invention.
Fig. 7 is the springform spirogram during model excavation of embodiment of the present invention advance of the face 150m.
Fig. 8 is the springform spirogram during model excavation of advance of the face 190m in the embodiment of the present invention.
Fig. 9 is the springform spirogram during model excavation of advance of the face 210m in the embodiment of the present invention.
Figure 10 is the springform spirogram during model excavation of advance of the face 240m in the embodiment of the present invention.
Figure 11 is the springform spirogram during model excavation of advance of the face 240m in the embodiment of the present invention.
Figure 12 is the springform spirogram during model excavation of embodiment of the present invention advance of the face 270m.
Figure 13 is the springform spirogram during model excavation of advance of the face 300m in the embodiment of the present invention.
Figure 14 is the springform spirogram during model excavation of embodiment of the present invention advance of the face 330m.
Figure 15 is the springform spirogram during model excavation of embodiment of the present invention advance of the face 360m.
Figure 16 is the springform spirogram during model excavation of advance of the face 400m in the embodiment of the present invention.
Figure 17 is model I-14 stress vector figure in the embodiment of the present invention.
In figure, (a) advance of the face 30m；(b) advance of the face 60m；(c) advance of the face 90m；(d) advance of the face 120m；(e) advance of the face 210m；(f) advance of the face 240m
Figure 18 is the height of water flowing fractured zone growth course figure of I-14 model in the embodiment of the present invention.
Figure 19 is surface subsidence value and the face advanced distance relation of I-14 model in the embodiment of the present invention.
Detailed description of the invention
In order to make objects and advantages of the present invention clearer, below in conjunction with embodiment, the present invention is further elaborated.Should be appreciated that specific embodiment described herein, only in order to explain the present invention, is not intended to limit the present invention.
Embodiment
There are boring H1504, H1604 and H1704 in the range of colliery 1121 fully-mechanized mining working in following example, in conjunction with Mine Integrated block diagram, using H1604 boring as main, set up experimental model.See Tables 1 and 2.The present embodiment uses the true rupture process analysis system RFPA seepage flow version of power soft science and technology (Dalian) limited company exploitation to carry out numerical simulation.The main mining coal seam overlying strata used produce the method for discrimination of absciss layer, comprises the steps:
Rock stratum different for which floor lithology is formed combination key stratum, and a single rock stratum can be regarded as, by several single rock stratum combination some groups of compoboards of composition, in each group of compoboard, the load of 1st layer rock stratum nethermost to this group, n-th layer rock stratum:
In formula, E1, E2..., EnFor the elastic modelling quantity of each layer, n is this group rock stratum number (from bottom to top arrangement)；h1, h2..., hnThickness for each layer；ρ1, ρ2..., ρnDensity for each layer.If calculating qn,1>qn+1,1, just the 1st layer of rock stratum is applied less than load by explanation n+1 layer rock stratum, then can conclude that product absciss layer between n-th layer and (n+1)th layer；
The generation of absciss layer must meet inequality, qn,1>qn+1,1, i.e.
After arrangement absciss layer occur discriminant:
If meeting above-mentioned discriminant, then judge main mining coal seam overlying strata have produced absciss layer.
No. 2 Seam Roof And Floor lithology abridged tables of table 1 fully-mechanized mining working
2 No. 2 coal seam overlying strata main physical parameters of table
Whole model is made up of 12 layers of coal rock layer, and wherein above No. 2 coal seams, the 4th floor is main sandstone aquifer, and the 3rd mfs layer is water-resisting key strata.When this layer of water-resisting key strata is 20m, it is 600m that model prolongs strike length, a height of 300m, is divided into 250 × 263 totally 65750 unit, and rock mass is solely subjected to weight stress and water pressure.Boundary condition is both ends horizontal constraint, and bottom is fixed, and sets periphery as water proof border.The process edge that water producing fractures is grown is simulated by step excavation.Model is walked 100m to the left and is started excavation, advances 400m, mining height 5m altogether, and often step excavation 10m, is divided into 40 steps.Mathematical calculation model such as Fig. 1.
The water pressure in water-bearing layer is often ignored in research before, and this causes experimental result to have greater difference with actual.In order to overcome the impact of this factor, this experiment sets up numerical model according to practical situation, consider the impact that water flowing fractured zone is grown in recovery process by the hydraulic pressure in water-bearing layer in a model, compose the head of 50m, research water producing fractures law of development and the seam mining impact on overlying strata sandstone aquifer to water-bearing layer during modeling.
(1) absciss layer generation evolution
In coal mine underground recovery process, along with the passage of underground mining work surface, destroying the natural stress equilibrium state of overlying strata body, goaf overlying rock is due to the effect of the secondary stress field of redistribution, and the stress within rock stratum also redistributes to reach new balance.Overlying rock to rupture, inbreak, bend and sink, and this motion change the most progressively develops, and forms subsidence to earth's surface.Overlying rock due to lithology, distance position, goaf, the difference of thickness etc., and causing the motion of rock stratum is not that Complete Synchronization is consistent, thus occur in that overlying strata motion displacement difference on vertical, form absciss layer space.
Due to the 1121 work surface average 20m of water-resisting key strata thickness, therefore this section by the water-resisting key strata (mud stone) with 20m as main study subject, simulation working face mining during, overlying strata separation layer development characteristics and water flowing fractured zone development height rule.
In order to eliminate boundary effect, each side stay 100m coal column at distance model.When excavation the 1st step, i.e. during advance of the face 10m, work surface overlying strata are not caving phenomenon (Fig. 2)；
When excavation the 3rd step, when the advance of the face to 30m, the outstanding dew of overlying rock, after bending rock beam hangs the span that the span of dew reaches certain under gravity, bending sedimentation develops into critical intensity, at the cracking ends of rock beam, there is absciss layer phenomenon in siltstone during coal seam is directly pushed up, and be caving for the first time along with part sandstone, there is the trend upwards developed in the crack in overlying strata, and now caving zone height is about 8.5m (Fig. 3).
When excavating 6 step, when i.e. excavation length reaches 60m, roof weighting is more fierce, occurs that second time is caving on a large scale.Absciss layer the most upwards develops, and in overlying strata, uppermost bedseparated fissures distance roof reaches 18m, absciss layer crack Breadth Maximum 5m.Being caving angle is 54 °, and caving zone height is 12.5m.Now in overlying strata, the subsoil water in the 2nd layer of rock stratum will import goaf, and i.e. IIUnderWater-bearing layer is switched on, and IIOnWater-bearing layer is not turned on.As shown in Figure 4.
When excavation to the 9th step, i.e. excavation length is 90m, between the 3rd layer of the most crucial water barrier (mud stone) in rock stratum and the 4th layer of rock stratum i.e. basal water (lithology is gritstone), for the first time absciss layer occurs.Now absciss layer Breadth Maximum is 10.2m, and height is 0.5m.This absciss layer is likely to become water space.Being caving and highly continue up extension growth, reach 25.8m, being caving angle is 60 °, as shown in Figure 5.
When excavation to the 12nd step, i.e. excavation length be 120m time, between the 3rd layer of rock stratum and the 4th layer before
Absciss layer continued development occurs.Now absciss layer Breadth Maximum is 20.5m, and height is 1.5m.This absciss layer storage capacity is increased dramatically.Being caving and highly continue up extension growth, reach 42.0m, being caving angle is 58 °.Now caving zone height has developed to the 3rd layer of rock stratum i.e. bottom of water-resisting key strata.3rd layer of inside occurs in that tension fracture near advance of the face direction, but this crack is not linked up with caving zone, as shown in Figure 6.
When excavation to the 15th step, i.e. excavation length is 150m, between the 3rd layer of rock stratum and the 4th layer, occur that abscission layer development arrives maximum before.Now absciss layer Breadth Maximum is 76.8m, and height is 2.5m.This absciss layer storage capacity increases to maximum.Being caving to the 15th step from excavation the 12nd step and the most can not extend up growth, keep 42.0m, being caving angle is 58 °.Now caving zone height has developed to the 3rd layer of rock stratum i.e. bottom of water-resisting key strata.The tension fracture continuous enlargement that 3rd layer of inside occurs near advance of the face direction, has the sign that same caving zone is through, as shown in Figure 7.
When excavation to the 19th step, i.e. excavation length is 190m, occurring that absciss layer occurs in second time on advance of the face direction between the 3rd layer of rock stratum and the 4th layer, absciss layer occurs the phenomenon that Guan Bi compression reduces on advance of the face direction for the first time.Now second time absciss layer Breadth Maximum is 5.5m, and height is 0.4m.Goaf, work surface rear part is upcoated rock stratum and is caving rear filling.The 2nd layer in roof there is again absciss layer.Be caving grown complete, keep 42.0m, being caving angle is 58 °.Now caving zone height has developed to the 3rd layer of rock stratum i.e. bottom of water-resisting key strata.As shown in Figure 8.
When excavation to the 21st step, i.e. excavation length is 210m, occurring that the absciss layer occurred for the second time continues extension on advance of the face direction between the 3rd layer of rock stratum and the 4th layer, absciss layer continues up compression reduction near advance of the face side for the first time.Now second time absciss layer Breadth Maximum is 17.3m, and height is 0.8m.Goaf, work surface rear part is upcoated rock stratum and is caving rear filling.In coal seam overlying strata the 2nd layer, abscission layer development is violent.Height of water flowing fractured zone is developed to the 3rd layer of middle part, about 50.0m.Now the 4th layer of water-bearing layer is in critical conduction mode.As shown in Figure 9.
When excavation to the 24th step, i.e. excavation length is 240m, advance of the face direction occurs between the 3rd layer of rock stratum and the 4th layer the absciss layer that second time occurs continues extension.Absciss layer Breadth Maximum is 33.3m, and height is 1.3m.In coal seam overlying strata the 2nd layer, abscission layer development is violent.5th layer has 140m, the absciss layer of maximum high about 0.3m with growth on the 4th layer of rock stratum contact surface.Height of water flowing fractured zone is developed to the 4th layer of middle part, about 60m.Now the 4th layer of water-bearing layer is in the conduction state.But finding in experimentation, owing to the 3rd layer of water barrier is mud stone, along with the advance of the face, overlying confined pressure increases, and water flowing fractured zone is crushed body by mud stone and fills, and has made water producing fractures up, has prevented the subsoil water in water-bearing layer to burst into work surface.There is large-scale subsidence basin in earth's surface.As shown in Figure 10 and Figure 11.
When excavation to the 27th step, i.e. excavation length is 270m, advance of the face direction occurs between the 3rd layer of rock stratum and the 4th layer third time absciss layer.Absciss layer Breadth Maximum reaches 33.6m for the third time, and height is 0.9m.Along with work surface persistently advances, overlying strata the 2nd layer of front abscission layer development in coal seam is violent, and rear absciss layer is compacted Guan Bi successively.In overlying strata, absciss layer continues up development, and the 6th layer absciss layer occurs with the 5th layer of rock stratum contact surface.Water flowing fractured zone development height is still 60m.As shown in figure 12.
When excavation to the 30th step, i.e. excavation length is 300m, the third time absciss layer occurred between the 3rd layer of rock stratum and the 4th layer on advance of the face direction continues to grow.Absciss layer no longer increases for the second time, and absciss layer is gradually compacted for the first time, and length and height are all greatly reduced.6th layer increases suddenly with the 5th layer of rock stratum contact surface absciss layer, and length is up to 152m, and water flowing fractured zone development height is still 60m.As shown in figure 13.
When excavation to the 33rd step, i.e. excavation length is 330m, the third time absciss layer occurred between the 3rd layer of rock stratum and the 4th layer on advance of the face direction continues to grow.This absciss layer width reaches 15.9m, and maximum height reaches 1.1m, and absciss layer is gradually compacted for the second time, and length and height are all greatly reduced.6th layer is gradually compacted with the 5th layer of rock stratum contact surface absciss layer.Water flowing fractured zone development height is still 60m.As shown in figure 14.
When excavation to the 36th step, i.e. excavation length is 360m, the third time absciss layer occurred between the 3rd layer of rock stratum and the 4th layer on advance of the face direction continues to grow.This absciss layer width reaches 57.3m, and maximum height reaches 2.1m, second time and the absciss layer formed for the first time and is the most gradually compacted, and length and height are all greatly reduced.6th layer is gradually compacted with the 5th layer of rock stratum contact surface absciss layer.Now overlying strata all have abscission layer development from the 2nd layer to the 8th layer rock stratum contact surface and inside.Water flowing fractured zone development height is still 60m, as shown in figure 15.
When excavation to the 40th step, i.e. excavation length is 400m (adopting eventually), the 4th absciss layer occurred between the 3rd layer of rock stratum and the 4th layer on advance of the face direction continues to grow.This absciss layer width reaches 9.3m, and maximum height reaches 0.8m, and the absciss layer formed for the third time, for the second time and for the first time is the most gradually compacted, and length and height are all greatly reduced.6th layer is gradually compacted with the 5th layer of rock stratum contact surface absciss layer.Now overlying strata all have abscission layer development, rear absciss layer to be gradually compacted from the 2nd layer to the 8th layer rock stratum contact surface and inside.Water flowing fractured zone development height is still 60m, as shown in figure 16.
(2) water flowing fractured zone growth course
Stress vector figure can reflect in recovery process, the situation of change of cover stress and the developmental state of water producing fractures.The seepage flow trace development maximum height of Figure 17 Green is water flowing fractured zone and grows maximum height.Make along with work surface pushes ahead the stress vector figure (see Figure 18) when height of water flowing fractured zone changes, then it is analyzed, it can be deduced that the feature that water flowing fractured zone is grown.
From Figure 17 (a), as advance of the face 30m, major developmental tensile stress in the middle part of goaf, and two ends, goaf are mainly extrusion stress, the direction of stress is being closed at goaf based on horizontal direction, and distance place farther out, goaf is the most tiltedly pointed in the middle part of goaf.In conjunction with Fig. 2, now rock stratum there occurs rupture, absciss layer phenomenon, but have water guide meaning crack only goaf top plate exist, height be 8.5m.When the advance of the face to 60m, the size generation large change of stress distribution, but stress direction change is inconspicuous, and height of water flowing fractured zone now increases, for 12.5m (Figure 17 (b)).Along with carrying out further of exploitation, the characteristic distributions in extrusion stress direction is still so that as horizontal direction, both sides are for V-arrangement in the middle part of goaf, and numerical value increased；Meanwhile, the tensile stress distribution in both sides, goaf gradually expands, and numerical value increases.When the advance of the face to 90m, based on horizontal direction, then there is tensile stress in the middle part of goaf in both sides in extrusion stress, and now height of water flowing fractured zone increases to 25.8m (Figure 17 (c)).When the advance of the face to 120m, in lateral symmetry, extrusion stress based on horizontal direction, and occurs that tensile stress, height of water flowing fractured zone occur in the surface in goaf in the middle part of goaf, highly increase to 42m (Figure 17 (d)).Hereafter there is relative stabilization, the i.e. advance of the face to 140m, 160m, 180 in height of water flowing fractured zone, and during 200m, height of water flowing fractured zone is the most a small amount of to be increased；When work surface continues to be advanced to 210m, the extruding of above goaf basement rock rock stratum, tensile stress the most upwards develop, and the extruding of near surface ravelly ground, tensile stress develop downwards, and stress distribution feature is like before.Now, height of water flowing fractured zone shows the feature that saltatory increases, and increases to 50m (Figure 17 (e)).When work surface continues to be advanced to 240m, above goaf basement rock rock stratum tensile stress the most upwards develops, and the extruding of earth's surface ravelly ground, tensile stress develop downwards, and height of water flowing fractured zone increases to 620m (Figure 17 (f)).Hereafter, height of water flowing fractured zone the most no longer continues along with the lasting propelling of work surface to increase.That is this work surface maximum height of water flowing fractured zone is 60m.
In order to preferably study water flowing fractured zone law of development, in simulation being calculated, the water flowing fractured zone developmental state of each step is extracted data and is made water flowing fractured zone development height figure, as shown in figure 18.As can be seen from Figure 18, the advance distance of work surface directly affects the development height in crack, water flowing fractured zone has from occurring, growing (rising), maximum height, stable growth course along with the exploitation in coal seam, the sedimentation of overlying strata, absciss layer, the formation of destruction, and this meets the law of development of water flowing fractured zone.This time experiment is it can be seen that twice jump evolution occurs in height of water flowing fractured zone.During advance of the face 30m, fissure zone gradually manifests, and height is 8.5m；The advance of the face is to before 60m, although Overburden Rock Failure scope gradually expands, but crack height is still grown slowly；Between 60m to 120m, fissure zone height steeply rises, 12.5m saltatory rise to 42m, the development height of this stage fissure zone almost with face advanced distance proportional, for first time jump evolution.Hereafter water flowing fractured zone temporary stabilization is at this height, absciss layer phenomenon occurs, is typical cantilever beam structure, prevent the upwards development in crack above work surface.When work surface increases to 240m from 190m, fissure zone height steeply rises, 42m saltatory rise to 60m, the development height of this stage fissure zone almost with face advanced distance proportional, for second time jump evolution.
(3) surface subsidence process
As seen from Figure 19, along with advance of the face mining subsidence starting distance is about 85m, when the advance of the face to 240m, surface subsidence is increased dramatically to 1560mm.When the advance of the face to 350m, surface subsidence, to maximum, reaches sufficient mining, and maximum surface subsidence is 3500mm.
Understanding in sum, the development height of water flowing fractured zone determines gathering and the loss of absciss layer water, and cover stress distribution is the risk factor that water flowing fractured zone is grown, and surface movement and deformation is its performance on earth's surface.Twice saltatory of the height of water flowing fractured zone in Figure 18 increases corresponding earth's surface and violent sinking performance all occurs, and absciss layer water has the phenomenon of seepage flow.Therefore, abscission layer development feature and water flowing fractured zone are grown, cover stress is distributed and surface movement sinks to having close relatedness.
Numerical simulation result is analyzed
Due to seam mining, the water in lower leaf water-bearing layer is imported into work surface, and because watery is weak, norm al discharge rate is 15m3/h.Along with the advance of the face, between water barrier and higher slice water-bearing layer, absciss layer occurs, higher slice aquifer recharge absciss layer forms absciss layer water, absciss layer is become one metastable " water storage tank ", and the most whole stratum is in poised state, along with the continuation of work surface advances, roof destruction intensity strengthens, original balance is broken, and absciss layer instant water is burst into mine, causes gushing water.Absciss layer water runs off totally in the short period of time, reduce with absciss layer space, inflow rate of mine water recovers normal, along with periodic weighting, main roof caving, fracture, re-form absciss layer water, propelling along with work surface, new balance is broken again, causes again bursting into down-hole of absciss layer water, forms new water inrush accident.Gushing water goes round and begins again, the generation of cycle dynamics.
In sum,
(1) the main mining coal seam in study area overlying strata can form absciss layer water between the 3rd layer of thick for 20m mud stone and the 4th layer of 40m thickness sandstone.Through numerical simulation calculation, when the advance of the face to 100m, large-scale water storage absciss layer occurs for the first time.It is 76.8m that this absciss layer finally develops into width, and height is the absciss layer water space of 2.5m.
(2) numerical simulation shows, the work surface of exploitation 400m, absciss layer water damage periodically occurs 4 times, the most often advances 100m, forms 1 large-scale absciss layer water damage.Therefore, the cycle dynamics that absciss layer water damage is formed is about 100m.Along with the propelling of work surface, absciss layer forms rear continued development to extreme value, the most gradually atrophy.The scale that absciss layer cycle dynamics occurs is different, wherein with the absciss layer Space Scale that occurs for the first time as maximum.The periodicity that absciss layer water damage occurs is the most disrumpent feelings with roof close relationship.
(3) absciss layer water damage is periodically formed 4 times, and corresponding water flowing fractured zone is grown and presented 2 convex jumping property growths, and therefore the cycle dynamics feature of absciss layer water damage is grown with height of water flowing fractured zone and presented positive correlation.The cycle absciss layer water damage that absciss layer is formed is formed, and is owing to along with the lasting propelling of work surface, the convex development that jumps of water flowing fractured zone development height turns on the result of absciss layer water body.Study area height of water flowing fractured zone is finally 60m.
The above is only the preferred embodiment of the present invention; it should be pointed out that, for those skilled in the art, under the premise without departing from the principles of the invention; can also make some improvements and modifications, these improvements and modifications also should be regarded as protection scope of the present invention.
Claims (1)
1. colliery absciss layer water cycle gushing water causes the method for numerical simulation of calamity, it is characterised in that comprise the steps:
Setting up mathematical calculation model, whole model is made up of 12 layers of coal rock layer, and wherein above No. 2 coal seams, the 4th floor is
Main sandstone aquifer, the 3rd mfs layer is water-resisting key strata；When this layer of water-resisting key strata is 20m,
It is 600m that model prolongs strike length, a height of 300m, is divided into 250 × 263 totally 65750 unit, rock mass
It is solely subjected to weight stress and water pressure；Boundary condition is both ends horizontal constraint, and bottom fixes, set periphery as
Water proof border.The process edge that water producing fractures is grown is simulated by step excavation；Model is walked 100m to the left and is opened
Begin excavation, advances 400m, mining height 5m altogether, and often step excavation 10m, is divided into 40 steps.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610368278.3A CN106021762A (en) | 2016-05-30 | 2016-05-30 | Numerical simulation method for coal mine abscission layer water periodical water inrush disaster |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610368278.3A CN106021762A (en) | 2016-05-30 | 2016-05-30 | Numerical simulation method for coal mine abscission layer water periodical water inrush disaster |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106021762A true CN106021762A (en) | 2016-10-12 |
Family
ID=57091457
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610368278.3A Pending CN106021762A (en) | 2016-05-30 | 2016-05-30 | Numerical simulation method for coal mine abscission layer water periodical water inrush disaster |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106021762A (en) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106405678A (en) * | 2016-11-22 | 2017-02-15 | 临沂大学 | Mining overlying strata water flowing fracture zone height detection method based on stress monitoring |
| CN108038282A (en) * | 2017-11-30 | 2018-05-15 | 安徽理工大学 | Load transmission method for numerical simulation is exploited under a kind of unconsolidated confined aquifer |
| CN109470624A (en) * | 2018-12-19 | 2019-03-15 | 山东科技大学 | A kind of interior overlying rock absciss layer water-filling and simulation of infiltration system |
| CN109611084A (en) * | 2018-11-20 | 2019-04-12 | 中国矿业大学 | A kind of mining overburden can ponding separating location method of discrimination |
| CN109667577A (en) * | 2019-02-14 | 2019-04-23 | 中国矿业大学 | Absciss layer hydrodynamic(al) water pressure change procedure analogy method when a kind of upper rock stratum fracture of absciss layer |
| CN109813486A (en) * | 2019-02-14 | 2019-05-28 | 中国矿业大学 | Absciss layer hydrodynamic(al) water pressure monitoring device when a kind of upper rock stratum fracture of absciss layer |
| CN109933954A (en) * | 2019-04-12 | 2019-06-25 | 宿州学院 | A kind of identification and diaster prevention and control method of remnants absciss layer |
| CN116644614A (en) * | 2023-06-30 | 2023-08-25 | 中国矿业大学 | Mining overburden water burst separation layer position identification method |
| CN116957325A (en) * | 2023-07-05 | 2023-10-27 | 中国矿业大学 | Roof separation layer water inrush prediction method based on multiple nonlinear indexes |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105332738A (en) * | 2015-08-24 | 2016-02-17 | 西安科技大学 | Early warning method of coal mining subsidence |
| CN105548483A (en) * | 2015-12-30 | 2016-05-04 | 山东科技大学 | System and method for simulating water inrush from coal mining floor above artesian water aquifers |
-
2016
- 2016-05-30 CN CN201610368278.3A patent/CN106021762A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105332738A (en) * | 2015-08-24 | 2016-02-17 | 西安科技大学 | Early warning method of coal mining subsidence |
| CN105548483A (en) * | 2015-12-30 | 2016-05-04 | 山东科技大学 | System and method for simulating water inrush from coal mining floor above artesian water aquifers |
Non-Patent Citations (3)
| Title |
|---|
| 于庆磊: "煤层开采引起覆岩三带分布规律的数值模拟", 《第九届全国岩石力学与工程学术大会论文集》 * |
| 朱卫兵等: "覆岩离层区积水引发的采场突水机制", 《岩石力学与工程学报》 * |
| 王长俊等: "关于岩层移动过程中的离层位置判别", 《煤炭技术》 * |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106405678A (en) * | 2016-11-22 | 2017-02-15 | 临沂大学 | Mining overlying strata water flowing fracture zone height detection method based on stress monitoring |
| CN108038282A (en) * | 2017-11-30 | 2018-05-15 | 安徽理工大学 | Load transmission method for numerical simulation is exploited under a kind of unconsolidated confined aquifer |
| CN109611084B (en) * | 2018-11-20 | 2021-10-22 | 中国矿业大学 | Method for judging position of water-accumulatible separation layer of mining overburden rock |
| CN109611084A (en) * | 2018-11-20 | 2019-04-12 | 中国矿业大学 | A kind of mining overburden can ponding separating location method of discrimination |
| CN109470624A (en) * | 2018-12-19 | 2019-03-15 | 山东科技大学 | A kind of interior overlying rock absciss layer water-filling and simulation of infiltration system |
| CN109667577A (en) * | 2019-02-14 | 2019-04-23 | 中国矿业大学 | Absciss layer hydrodynamic(al) water pressure change procedure analogy method when a kind of upper rock stratum fracture of absciss layer |
| CN109813486A (en) * | 2019-02-14 | 2019-05-28 | 中国矿业大学 | Absciss layer hydrodynamic(al) water pressure monitoring device when a kind of upper rock stratum fracture of absciss layer |
| CN109933954A (en) * | 2019-04-12 | 2019-06-25 | 宿州学院 | A kind of identification and diaster prevention and control method of remnants absciss layer |
| CN109933954B (en) * | 2019-04-12 | 2022-06-24 | 宿州学院 | Residual separation layer identification and disaster prevention method |
| CN116644614A (en) * | 2023-06-30 | 2023-08-25 | 中国矿业大学 | Mining overburden water burst separation layer position identification method |
| CN116644614B (en) * | 2023-06-30 | 2023-11-14 | 中国矿业大学 | Mining overburden water burst separation layer position identification method |
| CN116957325A (en) * | 2023-07-05 | 2023-10-27 | 中国矿业大学 | Roof separation layer water inrush prediction method based on multiple nonlinear indexes |
| CN116957325B (en) * | 2023-07-05 | 2024-02-13 | 中国矿业大学 | Roof separation layer water inrush prediction method based on multiple nonlinear indexes |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106021762A (en) | Numerical simulation method for coal mine abscission layer water periodical water inrush disaster | |
| Miao et al. | The height of fractured water-conducting zone in undermined rock strata | |
| CN104314570B (en) | Long thick ore body Wall ore reclaims and the transitional exploitation method of open air trestle | |
| Tang | Applications of rock failure process analysis (RFPA) method | |
| CN103899352B (en) | The design of solid filling Full Ratio and control method in coal mining | |
| CN104775816B (en) | Isolated grouting-filling mining method for overburden rock without village migration under cut-side local pressed-coal conditions | |
| CN113175325B (en) | Coal and intergrown sandstone type uranium ore coordinated mining method based on key layer protection | |
| CN107939399A (en) | A kind of more stage casings level pillar recovery method and into line structure under exploitation pattern at the same time | |
| CN105332738A (en) | Early warning method of coal mining subsidence | |
| CN109162713B (en) | Coal-water dual-resource mine mining method without changing overlying strata hydrogeological conditions | |
| CN103758519A (en) | Thick and large mineral deposit stage subsequent filling method mining piecewise optimization design and implementation method | |
| CN115030722A (en) | Efficient water-retaining coal mining method for goaf lag filling | |
| CN102678119A (en) | High-efficiency filling mining method for controlling rock mass movement | |
| Zhai et al. | A reinforcement method of floor grouting in high-water pressure working face of coal mines: a case study in Luxi coal mine, North China | |
| Liu et al. | Characteristics analysis of roof overburden fracture in thick coal seam in deep mining and engineering application of super high water material in backfill mining | |
| Hu et al. | Research on the reasonable width of the waterproof coal pillar during the mining of a shallow coal seam located close to a reservoir | |
| CN104196569B (en) | The Forecasting Methodology that on filling coal mining work surface, whether water barrier is stable | |
| Shen et al. | Research on novel method of gob-side entry retaining under the synergistic effect of roof cutting and roadside filling in thick coal seams | |
| Dong et al. | Study on overburden stability and development height of water flowing fractured zone in roadway mining with cemented backfill | |
| He | Numerical simulation analysis of karst tunnel water bursting movement | |
| Li et al. | Dynamic evolution characteristics and mechanism of surrounding rock fractures during the repeated mining of closed distance deep coal seam | |
| Lai et al. | Research of the Backfill Body Compaction Ratio Based on Upward Backfill Safety Mining of the Close‐Distance Coal Seam Group | |
| CN108921350A (en) | Metal mine goaf top covering rockmass Mobile Space-time law forecasting method | |
| Guo et al. | Ground control by L-shaped cemented paste backfilling technology in underground coal seam mining: a case study | |
| Fan et al. | Evolution law of mining-induced overburden stratum stress and fracture fields in inclined coal seam |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20161012 |