CN105868472A - Gob air leakage flow field dynamic numerical simulation method based on deformation geometry - Google Patents
Gob air leakage flow field dynamic numerical simulation method based on deformation geometry Download PDFInfo
- Publication number
- CN105868472A CN105868472A CN201610187845.5A CN201610187845A CN105868472A CN 105868472 A CN105868472 A CN 105868472A CN 201610187845 A CN201610187845 A CN 201610187845A CN 105868472 A CN105868472 A CN 105868472A
- Authority
- CN
- China
- Prior art keywords
- goaf
- coal
- gob
- face
- flow field
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 70
- 238000004088 simulation Methods 0.000 title claims abstract description 40
- 239000003245 coal Substances 0.000 claims abstract description 128
- 230000008859 change Effects 0.000 claims abstract description 27
- 239000012530 fluid Substances 0.000 claims abstract description 15
- 238000004364 calculation method Methods 0.000 claims abstract description 14
- 230000035699 permeability Effects 0.000 claims description 20
- 238000004519 manufacturing process Methods 0.000 claims description 10
- 238000006073 displacement reaction Methods 0.000 claims description 8
- 239000011435 rock Substances 0.000 claims description 4
- 238000004886 process control Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 abstract description 21
- 238000009826 distribution Methods 0.000 abstract description 19
- 238000005065 mining Methods 0.000 abstract description 7
- 230000001737 promoting effect Effects 0.000 abstract 2
- 239000000463 material Substances 0.000 description 17
- 238000011084 recovery Methods 0.000 description 14
- 238000007796 conventional method Methods 0.000 description 9
- 230000009931 harmful effect Effects 0.000 description 5
- 230000002269 spontaneous effect Effects 0.000 description 5
- 231100000331 toxic Toxicity 0.000 description 5
- 230000002588 toxic effect Effects 0.000 description 5
- 230000015572 biosynthetic process Effects 0.000 description 4
- 238000000205 computational method Methods 0.000 description 4
- 230000009467 reduction Effects 0.000 description 4
- 230000001052 transient effect Effects 0.000 description 4
- 238000013459 approach Methods 0.000 description 3
- 238000005056 compaction Methods 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 238000011105 stabilization Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000002485 combustion reaction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 230000002085 persistent effect Effects 0.000 description 2
- 239000000047 product Substances 0.000 description 2
- 230000006641 stabilisation Effects 0.000 description 2
- 208000032843 Hemorrhage Diseases 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000012467 final product Substances 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000008595 infiltration Effects 0.000 description 1
- 238000001764 infiltration Methods 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000002045 lasting effect Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
- 230000008733 trauma Effects 0.000 description 1
- 238000009423 ventilation Methods 0.000 description 1
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 provides a gob air leakage flow field dynamic numerical simulation method based on deformation geometry. The method comprises the steps that an initial gob geometric model for mining on a coal face is established according to a gob formed after mining is conducted for one day at the daily average promoting speed; a free deformation area, a fixed boundary and a moving boundary of the gob are arranged; a gob air leakage flow field is arranged; the initial gob geometric model is divided into boxes, and gob air leakage flow field numerical simulation calculation is carried out; the gob geometric model is promoted step by step at the daily average promoting speed and is divided into boxes again, and gob air leakage flow field numerical simulation calculation is carried out till the gob geometric model is promoted to a stop line or the gob geometric model does not arrive at the stop line but the set calculation stop time is up. According to the method, deformation geometry is adopted for controlling the dynamic evolution process of boundary movement and area changes of the coal face and the gob, and a free and porous medium fluid flow control working face and the flow field of the gob are coupled to calculate and describe the dynamic change and distribution law of the gob air leakage flow field.
Description
Technical field
The present invention relates to goaf air leaking analogue technique field, be specifically related to a kind of goaf air leaking based on deformation geometry
Flow field dynamic numeric simulation method.
Background technology
Goaf air leaking is the key factor affecting work surface safe working, leaks out and causes spontaneous combustion in goaf and induce
Goaf toxic and harmful is gushed out to work surface in a large number, therefore, it is necessary to clear and definite goaf air leaking rule is Safety of Coal Mine Production
Theoretical foundation and technical support are provided.Conventional stope wind flow field numerical value resolves mainly has 20th century 80 with emulation mode
～the self-compiling program stope Flow Field Calculation method based on formula translation of the nineties, the G3 software emulation side since 21 century
Method, bi-coordinate system method, moving coordinate system, method, also have business software PHOENICS, fluent, comsol etc. of commonly using square
Method, the level that any of the above method improves stope flow field Theoretical Calculation to a certain extent to the simulation calculation of flow field in goaf is
Work surface safety in production provides guidance, but during these method parts can not realize the advance of the face, the dynamic of goaf becomes
Change, is based partially on self-compiling program exploitation and the technology of obtaining and theoretical difficulty have a foregone conclusion relatively greatly in the promotion and application of the industry
Limit.
Summary of the invention
It is an object of the invention to provide a kind of goaf air leaking flow field based on deformation geometry dynamic numeric simulation method.
The technical scheme is that
A kind of goaf air leaking flow field based on deformation geometry dynamic numeric simulation method, including:
According to the goaf formed after exploiting a day with every per day fltting speed, set up at the beginning of coal-face exploitation
Beginning goaf geometric model: advance the speed for goaf length with every per day fltting speed, controls mined out based on deformation geometry
District's geometric model extends with the back production pusher of coal-face, arranges air intake lane and the return airway of coal-face；
Set Free Transform region, goaf, goaf fixed boundary, goaf moving boundary；
Goaf air leaking flow field is set: goaf is set to porous media flows region, is set to by coal-face
Free fluid flow region, is with coal-face fltting speed by goaf bulking factor, advances time, y-coordinate control
The variable of system, is with the change of goaf bulking factor by the porosity of erecting dense individual props coal petrography, permeability, something lost coal thickness
Variable；The air intake lane of coal-face is blast velocity control border, and the return airway of coal-face is that blast controls border, uses wind
Pressure reduction pointwise is distributed along the dividing potential drop of process control goaf external boundary；
By initial goaf geometric model grid division and carry out goaf air leaking Numerical Simulation calculate, including goaf
Bulking factor calculating, the porosity of erecting dense individual props coal petrography, permeability, something lost coal thickness, blast, speed calculation of leaking out；
With every per day fltting speed iterative method goaf geometric model and again subdivision grid and carry out goaf leakage
Numerical Simulation of Wind calculates, and stops adopting line until being advanced into or does not reaches to complete when setting is held the clock mined out to stopping adopting line
Leak out Flow Field Numerical Calculation in district.
Initial goaf geometric model is set as Free Transform region, described goaf；Open-off cut is set as described solid
Deckle circle, as global coordinate system.
Described goaf moving boundary includes x direction is moved, y direction is fixing single-way moving border and y direction is fixed, x
The Moving and free boundary that direction is freely pushed ahead by described single-way moving border.
The x direction translational speed on described single-way moving border is coal-face drawing speed, described single-way moving border
X direction displacement use grid displacement according to fltting speed and to advance time product to specify.
Described single-way moving border includes the side up and down in coal face and goaf.
Beneficial effect:
Gob edge dynamic mobile and goaf regional dynamics that this method advances continuously based on coal-face increase
Situ extraction actual, use the built-in deformation geometry method in comsol software to control coal-face and gob edge
Mobile and the Dynamic Evolution of regional change, employing freedom and fluid flow in porous medium control work surface and the stream in goaf
, the two coupling calculates and dynamically changes with the goaf air leaking flow field in factoid stope Dynamic Evolution and be distributed rule
Rule.The Dynamic profiling in goaf air leaking flow field during the advance of the face can be realized by this method, both meet engineering site
Reality also can improve Theoretical Calculation level, can for the prediction of the flow field dynamic changing process that leaks out in the Dynamic Evolution of goaf and
Describe and more convenient and method accurately is provided, moreover it is possible to pre-for problems such as the gushing out of goaf toxic and harmful, spontaneous fires
Survey the proposition with further investigation and administering method and base support is provided.
Accompanying drawing explanation
Fig. 1 is the initial goaf geometric model of the specific embodiment of the invention；
Erecting dense individual props coal petrography broken swollen system when Fig. 2 is the coal-face propelling different time of the specific embodiment of the invention
Number, (a) be coal-face advance 20 days time the broken swollen coefficient of erecting dense individual props coal petrography, (b) be coal-face advance 40 days time
The broken swollen coefficient of erecting dense individual props coal petrography, (c) be coal-face advance 60 days time the broken swollen coefficient of erecting dense individual props coal petrography, (d) is
Erecting dense individual props coal petrography broken swollen coefficient when coal-face advances 90 days；
Fig. 3 is the change continued with recovery time away from the open-off cut broken expansion factor of 30m across strike of the specific embodiment of the invention
Law；
Fig. 4 is that open-off cut y=85m, 45m of the specific embodiment of the invention continue along the broken swollen coefficient of journey with recovery time
Changing Pattern；
Fig. 5 is the Changing Pattern that the Residual coal in goaf thickness of the specific embodiment of the invention continues with recovery time, and (a) adopts
Residual coal in goaf thickness during the coal advance of the face 20 days, Residual coal in goaf thickness when (b) coal-face advances 40 days, (c) adopts
Residual coal in goaf thickness during the coal advance of the face 60 days, Residual coal in goaf thickness when (d) coal-face advances 90 days；
Fig. 6 be the specific embodiment of the invention away from open-off cut 30m along inclined direction lose coal thickness with recovery time continue
Changing Pattern, (a) is along the journey overall situation regularity of distribution, and (b) is the something lost coal thickness distribution rule near twice, and (c) is goaf
The something lost coal thickness distribution rule at middle part；
Fig. 7 is that at the open-off cut of the specific embodiment of the invention, y=85m, 45m lose coal thickness with recovery time continuity along journey
Changing Pattern；
Fig. 8 be the specific embodiment of the invention exploitation fltting speed 3m/d working condition under different times erecting dense individual props
The porosity distribution rule of coal petrography, (a) be coal-face advance 20 days time erecting dense individual props coal petrography porosity distribution rule,
(b) be coal-face advance 40 days time erecting dense individual props coal petrography porosity distribution rule, (c) be coal-face advance 60
It time erecting dense individual props coal petrography porosity distribution rule, (d) be coal-face advance 90 days time erecting dense individual props coal petrography
Porosity distribution rule；
Fig. 9 is that at the distance open-off cut 30m of the specific embodiment of the invention, the porosity of inbreak coal petrography changes over rule
Rule；
Figure 10 is the distance open-off cut 10m of the specific embodiment of the invention, and at y=45,85m, the porosity of inbreak coal petrography becomes
Law；
Figure 11 is the porosity attenuation law of the inbreak coal petrography of the specific embodiment of the invention；
Figure 12 is the permeability attenuation law of the specific embodiment of the invention；
Figure 13 be the specific embodiment of the invention open tangent line y=5m, at y=45m, y=85m, 3 points are with coal work
With coal-getter when velocity variations rule that face advances and time lengthening is leaked out does not changes with region, goaf under conventional method and border
Make face to advance and time lengthening is leaked out velocity variations rule；Wherein (a), (c), (e) are respectively described a position for opening tangent line y=
At 5m, y=45m, y=85m, 3 points advance with coal-face and time lengthening is leaked out velocity variations rule, (b), (d), (f)
Opening tangent line y=5m respectively, at y=45m, y=85m, 3 somes region, goaf and borders under conventional approaches do not change；
Figure 14 be the specific embodiment of the invention different fltting speeds under goaf air leaking rule, (a) is wind speed 0.8m/
Goaf air leaking rule under s, (b) is goaf air leaking rule under wind speed 1m/s, and (c) is goaf air leaking rule under wind speed 1.4m/s
Rule, (d) is goaf air leaking rule under wind speed 1.8m/s.
Detailed description of the invention
Below in conjunction with the accompanying drawings the detailed description of the invention of the present invention is elaborated.
A kind of goaf air leaking flow field based on deformation geometry dynamic numeric simulation method, including:
The goaf that step 1, basis are formed after exploiting one day with every per day fltting speed, sets up coal-face and opens
The initial goaf geometric model adopted: advance the speed with every per day fltting speed for goaf length, based on deformation geometry control
Goaf processed geometric model extends with the back production pusher of coal-face, arranges air intake lane and the return airway of coal-face；
Initial goaf geometric model in present embodiment is as it is shown in figure 1, coal work face length 90m, air entering and returning tunnel width
As a example by 3.5m, high 4m, back production area width about 7m, i.e. rib and goaf ultimate range 7m, as it is shown in figure 1, L1=3.5m, L2=
7m.Use deformation geometry to carry out goaf control, i.e. with 1 day as time step, be the mined out head of district with every per day fltting speed
Degree is advanced the speed and is extended goaf with the retrogressing of coal-face, first day mined out as shown in the territory A in Fig. 1 during initial setting up
District is formed, L3=vt 1 [d] (vt is every per day fltting speed, and unit is m/d), the mined out head of district in follow-up recovery process
Degree is continuously increased, and its a length of L3=vt t [d], in formula, d represents sky.Figure center line 3-3 is the location of starting cut beginning to adopt, simulation
In use rollback stoping method due to coal-face, be in coal-face move the horizontal seat of geometric coordinate system so arranging this line
Target zero point, i.e. x1=0m, the like.X2=-L3=-vt 1 [d], x3=x2-7m, in figure, the distance between x3, x4 is
Air entering and returning tunnel length in simulation process, in order to simplified model and saving calculate this segment length L4 in resource present embodiment
It is set to 3m, then x4=x3-3m.Territory C1, territory C2 are respectively return aircourse and intake part, and territory B is that coal-face is empty
Between, wherein omit the equipment such as coal-winning machine, cable trough, two row's hydraulic supports, by rib, roof and floor and various equipment for ventilation
Impact considers to get final product relative reduction in windage owing to omitting the error that above equipment is brought.Coal-face moves geometry and sits
The y0=0m of y-axis, y1=3.5m, y2=86.5m, y3=90m in mark system.
The territory of geometry deformation selects and controls: deformation geometry master control region is territory C1, territory C2, territory B, territory A；Set geometry
Deformation joint time be 2 i.e. geometry deformation control function be second order derived function, Free Transform sets Mesh Smoothing control method and uses
Winslow Mesh Smoothing method.
Step 2, setting Free Transform region, goaf, goaf fixed boundary, goaf moving boundary；
Initial goaf geometric model is set as Free Transform region, described goaf；Open-off cut is set as fixed edge
Boundary, as global coordinate system.
Goaf moving boundary includes x direction is moved, y direction is fixing single-way moving border and y direction is fixed, x direction
The Moving and free boundary freely pushed ahead by described single-way moving border.
The x direction translational speed on single-way moving border is coal-face drawing speed, the x side on described single-way moving border
Grid displacement is used according to fltting speed and to advance time product to specify to displacement.Single-way moving border include coal face and
The side up and down in goaf.
Present embodiment specifically sets as follows:
Set Free Transform region, goaf: due to during coal-face back production pusher this coal-face and
Goaf all moves forward, and therefore by territory C1, C2, territory B, territory A are disposed as Free Transform region.Free Transform region sets
Determine grid initial displacement dx0=0m, dy0=0m.
Setting goaf fixed boundary: the open-off cut in Fig. 1 arranges 3-3 and is set as fixed boundary, this border is not because mining
Work surface and the movement in goaf and change, the coordinate system of this fixed boundary is global coordinate system, i.e. the material on this border and several
What coordinate system is all fixing does not moves and deforms, it is intended that its x, y direction grid displacement is 0m always.
Setting goaf moving boundary one moving boundary 1 and moving boundary 2:
(1) moving boundary 1 is set: retreating with coal-face and goaf length is stepped up, taking time step is 1
My god, during coal-face pusher, border 2-1,3-1 are the single-way moving border that x direction is moved, y direction is fixing；X direction is moved
Dynamic speed is coal-face drawing speed.
(2) moving boundary 2 is set: coal-face advances and there is also another kind of y direction during the prolongation of goaf
Fixing do not move, moving boundary that x direction is freely pushed ahead by above-mentioned moving boundary 1, such border is concentrated mainly on coal mining
Work surface and the side up and down in goaf, such as 1-2,1-3,1-4,1-5,1-6,2-2,2-3,3-2,3-4 in Fig. 1, this kind of limit
The x direction on boundary is moved with work surface and is not provided with speed, and y does not moves in direction.
Step 3, goaf air leaking flow field is set: goaf is set to porous media flows region, by coal-face
It is set to free fluid flow region, goaf bulking factor is with coal-face fltting speed, advances time, y
The variable of coordinate control, is with goaf bulking factor by the porosity of erecting dense individual props coal petrography, permeability, something lost coal thickness
The variable of change；The air intake lane of coal-face is blast velocity control border, and the return airway of coal-face is that blast controls border,
Leeway pointwise is used to be distributed along the dividing potential drop of process control goaf external boundary；
Arranging freedom and fluid flow in porous medium region: free fluid flow region is territory C1, territory C2, territory B, porous is situated between
Matter fluid flow region is territory A；
(1) free fluid flow region
The control territory of free fluid flowing is territory C1, territory C2, territory B, and numerical computations uses equation (1)～(2) to control this district
The fluids free flow in territory moves:
In formula: F is by the airflow temperature of coal-face, humidity, density, coal-face in the free fluid flowing space
The buoyancy that the road discrepancy in elevation, gravity etc. cause up and down, unit is N/m3；ρ is thermal current density, and unit is kg/m3；P2 is coal-getter
Make the absolute blast in face, the i.e. air pressure of free turbulence, Pa；U is the vector flow velocity of air-flow, m/s；U is the dynamic viscosity of gas, Pa/
s2；For current gradient.
(2) fluid flow in porous medium region
The control territory of fluid flow in porous medium is territory A, and the equation of its fluid flow inside uses (3)～(4) to control
System:
In formula: εpPorosity for inbreak coal petrography；κ is the permeability of porous media；QbrFor flow quality in porous media
Flow, unit is kg/m3/s。
Porosity ε of erecting dense individual props coal petrographypComputational methods such as formula (1), permeability κ computational methods such as formula (2).
εp=1-1/kp (1)
κ=εp 3dp2/150(1-εp)2 (2)
In formula: kp is broken swollen coefficient, its computational methods such as formula (3)；Dp is the mean diameter of Residual coal in goaf, and unit is m.
Kp=kp0+(kpmax-kp0)exp(a1(vt·t+L3+x)(1-exp(a2·a3Min (y, L-y)))) (3)
In formula, kp0For initial broken swollen coefficient；kpmaxFor maximum broken swollen coefficient；a1、a2、a3It is respectively according to field condition true
Fixed constant term；Vt is the every per day fltting speed of coal-face, and unit is m/d.
Losing coal/float coal thickness h1 in goaf uses (3) to control.
H1=m1 (1-f1) kp (3)
In formula: h1 is for losing coal/float coal thickness, m；M1 is minable thickness, and present embodiment is according to Ruian company 014N1-1Combine
Put that work surface is actual takes it for 10.4m；Being cut to broken swollen coefficient, computational methods such as formula (3), f1 is the rate of extraction, present embodiment root
According to Ruian company 014N1-1Fully mechanized coal face is actual peeks according to formula (4).
Ne-leakage wall (i.e. without border of leaking out) is border 1-3,2-2,3-2,3-3,3-4,2-3,1-4,1-6,2-1,1-5.
Blast velocity control border is set and controls border with blast: border 1-2 is set to work surface air intake lane entrance, this entrance
Use normal direction inflow velocity to control, i.e. the air intake lane of coal-face is blast velocity control border；Border 1-1 is set to work surface
Return air exports, and uses outlet pressure to control, i.e. the return airway of coal-face is that blast controls border；By coal-face and adopting
The common boundary 3-1 of dead zone uses leeway pointwise to control border blast.
Step 4, by initial goaf geometric model grid division and carry out goaf air leaking Numerical Simulation calculate, including
Goaf bulking factor calculating, the porosity of erecting dense individual props coal petrography, permeability, something lost coal thickness, blast, leak out velometer
Calculate；
Grid division divides employing General Physics field to be controlled, and size of mesh opening uses special micronization processes, largest unit chi
Very little 0.5m, minimum cell size 0.02m, largest unit production rate 1.13, tortuosity 0.3.
Step 5, with every per day fltting speed iterative method goaf geometric model and again subdivision grid carrying out adopt
Dead zone leak out Numerical Simulation calculate, stop adopting line or not reaching complete when setting is held the clock to stopping adopting line until being advanced into
Become goaf air leaking Flow Field Numerical Calculation.
By the goaf air leaking Field Flow Numerical Simulation method based on deformation geometry of present embodiment, can analyze and adopt
Dead zone material parameter and flow field change rule:
(1) the broken swollen coefficient in goaf and something lost coal thickness Changing Pattern
(1.1) broken swollen coefficient
Do not consider only to adopt near open-off cut do not put coal region time, Fig. 2 (a)～(d) are respectively coal-face and advance 20
My god, 40 days, 60 days time, 90 days time the broken swollen coefficient of erecting dense individual props coal petrography, near coal-face and up and down inbreak near twice
The broken swollen coefficient of coal petrography is relatively big, similar with near work surface near open-off cut, and the porosity of broken swollen coefficient big inbreak coal petrography is inevitable
Greatly, therefore " twice ", " two lines " are easy spontaneous fire region under the conditions of mining production practice.
In Fig. 2, the broken swollen coefficient at deep, goaf and middle part is relatively low, and continuing along with recovery time, goaf length by
Gradually strengthening, some region of broken swollen coefficient also will change, such as Fig. 2～Fig. 3.
Fig. 3 is the Changing Pattern continued with recovery time away from the broken swollen coefficient of open-off cut 30m across strike, inclines at coal-face
On tilted direction, the broken swollen big middle part of coefficient twice is low about center line symmetrical (nearly horizontal coal-face can be considered symmetrical)；This region with
Open-off cut keeps constant relative to position, and continuing of coal-face recovery time, this region moves away from coal-face, broken swollen
Coefficient is gradually lowered, and initial stage of production broken swollen coefficient relatively low comparatively fast, exploitation the persistent period longer time as 70～90 days time, mined out
The broken swollen coefficient of district's central region is basically identical, and goaf 190～250m has been imbedded in this section of this region of time, existing according to ore deposit pressure
As the most basicly stable compacting in this region of rule.To sum up, coal-face advances backward, the process that goaf length constantly extends
In, each region in goaf experiencings the descending gradually decay of broken swollen coefficient and is finally fully compacted it in this region
Rear broken swollen coefficient maintains minimum process, so the goaf numerical simulation carried out under coal-face propulsioning condition can not be ignored
This principle key rule.
Fig. 4 is to describe at open-off cut, the Changing Pattern of y=45m and y=85m 2 broken swollen coefficient time at any time, D in figure
Representing and use the inventive method, C represents the conventional method not using deformation geometry method, occurs having same meaning in subsequent figure
Justice, shows in figure that the lasting broken swollen coefficient used under conditions of the inventive method in time is descending and progressively decays, goaf
Central region decay is faster than twice, the minimum under broken swollen coefficient stabilization of longer this point afterwards of recovery time to abundant compaction state
Value, at y=85m, 45m, minimum broken swollen coefficient is respectively 1.081,1.05；And show the conventional method not using deformation geometry method
Under conditions of, the size in goaf pre-sets fixing, so the broken swollen coefficient of any point is the most relevant in its coordinate, and normal
In rule method, material coordinate system is immovable, thus causes broken swollen coefficient to change the most in time, unrelated with the time, broken swollen system
Number always remains as 1.081,1.05 respectively with recovery time continuity, and therefore in coal-face progradation, be correlated with in goaf
Discounting for dynamic process during problem, the common movement that just can not realize gob material framework and geometric Framework can not be as
The real reaction real Changing Pattern of each physical field.
(1.2) coal thickness Changing Pattern is lost
Fig. 5 (a)～(d) are that coal-face advances Residual coal in goaf thickness distribution rule after different time, show in figure
It is thicker that twice lose coal, and will be above minable coal seam thickness near the region something lost coal thickness that twice wall is helped due to broken swollen impact；With
Time Residual coal in goaf distribution be also discrete, the something lost coal thickness lost in coal thickness and coverage beyond twice influence area
Gap is bigger.Residual coal in goaf thickness time Changing Pattern such as Fig. 6～Fig. 7.
Fig. 6 is the goaf distance open-off cut 30m something lost coal regularity of distribution along inclined direction, and (a) is along journey overall situation distribution rule
Rule, (b), (c) are respectively the something lost coal thickness distribution rule near twice and in the middle part of goaf, show and lose coal relatively near twice in figure
Thickness is 10.92m to the maximum more than the minable thickness of 10.4m, and this is not put at the bottom of the rate of extraction and broken swollen coefficient owing to only adopting near twice
Bigger reason.As shown in Fig. 6 (a), it is heterogeneous that fully mechanized coal face Residual coal in goaf is along inclined direction distributed, especially two
Adnexa exploitation in road carries out 90d and loses coal thickness period in the range of 8-10.92m, and central region is in the range of 1.1～1.15m, and
Owing in the middle part of goaf, compaction is preferable, lose coal thickness in the middle part of goaf minimum.Due to roof pressure, during with exploitation
Between continue goaf and lose coal thickness everywhere and be gradually lowered, as shown in Fig. 6 (b), (c).
Fig. 7 describes the float coal thickness Changing Pattern taking at open-off cut at y=45,85m 2, shows and conscientiously examine in figure
Consider goaf dynamic formation process time, this 2 lose the continuity in time of coal thickness be gradually lowered, its main cause be with
Time duration erecting dense individual props coal petrography under roof pressure effect is compacted repeatedly, and broken swollen coefficient diminishes and causes something lost coal thickness relatively low,
When 90 days minimum at y=85,45m lose coal thickness be respectively 8.664,1.092m；The something lost coal thickness of twice reduces in time
Speed is less than in the middle part of goaf；In conjunction with Fig. 6 (c), after the exploitation persistent period is longer, top board progressively steady pressure also tends to stable,
The trend that something lost coal thickness is the most gradually stablized and extended in time and reduces disappears substantially.Fig. 7 does not consider the dynamic change in goaf
Change and be directly the most at any time by the Residual coal in goaf thickness under the shape rigid condition in goaf in numerical simulation modeling process
Between change, the something lost coal thickness at y=85,45m remains 8.664 respectively, 1.092m, and reason does not becomes with broken swollen coefficient
The reason changed is similar to；The most do not consider the dynamic propagation process in goaf and when simulation directly by certain time under certain fltting speed
The direct fixed boundary in region, goaf of interior formation and area size, the Residual coal in goaf thickness obtained is above-mentioned setting exploitation
The something lost coal thickness of the last moment of time, this numerical computation method is ignored something lost coal thickness and is increased in the advance of the face and goaf
During change, belong to and carry out under the assumed condition that goaf moment is formed, be unscientific.
(2) porosity of erecting dense individual props coal petrography, permeability variation rule
Fig. 8 (a)～(d) are respectively under exploitation fltting speed 3m/d working condition, and coal-face advances 20,40,60,90
It time erecting dense individual props coal petrography porosity distribution, in figure show, the porosity of erecting dense individual props coal petrography is higher mined out in twice
In the middle part of district relatively low；Under the different mining period of contrast, each cloud atlas display arbitrary region is imbedded mined out after advance of the face certain time
Deep, district, the porosity of the inbreak coal petrography in this region extends in time and diminishes.Gob edge moves, region propagation process
The detailed rule of porosity change of inbreak coal petrography is as shown in figs. 9 to 12.
When Fig. 9 is the advance of the face 10,30,50,70,90 days, distance open-off cut 30m emits along coal-face incline direction
The porosity of coal breakage rock, along journey distribution and situation of change, shows the less twice in porosity deep of erecting dense individual props coal petrography relatively in figure
Greatly, the porosity reduction of this line region inbreak coal petrography and is continued with recovery time；Coal-face emits after advancing 30 days
The trend that the porosity of coal breakage rock reduces in time dies down, and illustrates that the porosity of inbreak coal petrography of arbitrary region in goaf is with adopting
Propelling and the growth in region, goaf of coal work surface and reduce, tend towards stability, such as Figure 10 institute after the certain depth of embedment goaf
Show.
The Changing Pattern of the porosity of erecting dense individual props coal petrography is more clearly described by Figure 10, for distance open-off cut
The porosity change rule of inbreak coal petrography at 10m, y=45,85m.Showing in figure, under the numerical simulation of the present invention, goaf is any
Point (with distance open-off cut 10m, at y=45,85m as a example by) inbreak coal petrography porosity in time be decay, in goaf
Between the porosity rate of decay of inbreak coal petrography be faster than twice, under the conditions of the fltting speed 3m/d that this legend is selected, in goaf
Heart region porosity of inbreak coal petrography after exploitation 60d starts stable about 0.05, at advance of the face 90d near upper and lower road
Time still there is larger porosity and still at continuous decrement.And conventional method (gob edge does not moves, region does not increases) in figure
The porosity of gained inbreak coal petrography is the porosity of the inbreak coal petrography under the last moment steady statue of this method, is not with work
Make face to advance and change, respectively 0.048,0.16, because initially setting up 3m/d*90d=when using conventional traditional method
The goaf of 270m length, the coordinate of arbitrfary point is all fixing, and the porosity of distance coal-face its inbreak coal petrography the most remote is more
Little, but can not change in time with the distance of coal-face, and in this method, grow out of nothing in goaf, from short to long
Progressively increase according to fltting speed, and the work surface border in goaf is also progressively moved according to fltting speed and time.
Under conventional method, the porosity of the inbreak coal petrography of certain point is a fixed constant, only with its residing for the position in goaf
Being equipped with pass, unrelated with the time, during dynamic formation, this porosity changes in time, because of the propelling goaf with work surface
During being formed and extending, deep regional is gradually compacted therefore its porosity and is gradually lowered.Theoretical according to round section joint ring, emit coal breakage
The porosity of rock position relative with goaf point relevant (x, y), and for goaf has been formed and exist arbitrfary point with
The distance of open-off cut is fixing, but be gradually increased with the distance of work surface, must experience progressively and repeatedly be compacted
Process, when it enters deeper region, goaf, roof pressure and compaction could be basicly stable, therefore, in difference exploitation speed
Under degree, different from the distance of coal-face in the identical time, fltting speed is the slowest the nearest with work surface and rate of increase is relative
Relatively low, so the time that embedment deeper region, goaf needs is relatively long, and when fltting speed is big, arbitrfary point or region with adopt
Distance between coal work surface quickly increases, and the speed required time in embedment goaf is shorter, and therefore fltting speed is bigger
Time the decay of every material parameter very fast, as Figure 11～Figure 12 is respectively away from porosity and the infiltration cutting eye tendency midpoint inbreak coal petrography
Rate attenuation law.
According to the porosity attenuation law of the inbreak coal petrography shown in Figure 11, fltting speed when coal-face advances 90 days
Under the conditions of 2m/d, the porosity of open-off cut tendency midpoint inbreak coal petrography is reduced to 0.04843, in the case of fltting speed 3～12m/d
When 90 days, the porosity of inbreak coal petrography has been reduced to 0.04766, and permeability reduction is 3.5077e-8m2, and before 90 days certain
One time stablized in this value.Under conditions of median is certain, the porosity of erecting dense individual props coal petrography is by broken swollen coefficient control
System, permeability is controlled by the porosity of inbreak coal petrography, therefore the hole of erecting dense individual props coal petrography on the premise of broken swollen coefficient stabilization
Gap rate and permeability should also be as being stable, fltting speed 3,5,6,8,10,12m/d time, porosity and permeability reach stable
Required time is 89,53,45,34,27,23d.
Gob edge and regional extent are fixed by conventional method, although use algorithm, moving coordinate system method windward
Etc. method, but gas concentration field, chemical reaction field, seepage field, atmosphere transmission can only be solved with the field such as diffusion field, thermal field in the time
On priority problem, it is impossible to solve or ignore such as critical material parameters such as broken swollen coefficient, the porosity of inbreak coal petrography, permeabilities
Transient changing problem, necessarily result in numerical simulation result owing to these basic critical material parameters can not realize transient changing
Mistake.
(3) goaf air leaking rule
The blast wind distributional class quasi-stability distribution in goaf in conventional method numerical simulation, it does not constitute aobvious with the time
The dependency relation write, each length calculating time step goaf of dynamic approach increases, so in each time step
The physical parameter redistributions such as the porosity of inbreak coal petrography in goaf, permeability, its blast etc. also redistributes, due to blast
Redistribution cause the redistribution leaked out, so dynamic approach is different from the air leakage distribution of conventional method, this is also
More meet goaf reality.
As a example by face propulsion speed 3m/d, during working face mining 90d, a length of 270m in goaf, this 270m are 90d
Inside gradually forming, therefore internal hemorrhage due to trauma wind velocity in goaf also will change with this process.Figure 13 (a), (c), (e) retouch respectively
Stating a position for opening tangent line y=5m, at y=45m, y=85m, 3 points advance with coal-face and time lengthening is leaked out speed
Changing Pattern, (b), (d), (f) open tangent line y=5m respectively, and at y=45m, y=85m, 3 points are the most mined out trivial
Territory and border do not change, and mined out head of district 270m is constant.Figure 13 (a), (c), (e) display form rule according to goaf and are configured
Numerical simulation result reflect the result of coincidence theory and practice, i.e. advance the progressively growth in goaf with coal-face,
The speed of leaking out of arbitrfary point, goaf is gradually lowered, and when this some embedment deeper region, goaf, hourglass wind velocity is stable at one relatively
Little value.
Figure 13 (b), (d), (f) are conventional method acquired results, show, the visitor not considering that goaf dynamically changes in figure
Seeing actual time, the wind speed that leaks out of arbitrfary point, goaf does not changes over time, and similar with each material parameter under the method is all placed in
This coordinate in goaf is relevant, and Figure 13 (b), (d), air speed value and the left column of (f) leak out velocity-stabilization period simultaneously
Air speed value is identical, absolutely proves the transient problem not considering that the numerical simulation etc. of goaf dynamic process is not the most standard, can only
Steady result is described.
Figure 14 (a)～(d) be wind speed 0.8,1,1.4,1.8m/s time difference fltting speed under open-off cut midpoint leak out speed
Variation tendency, shows in figure with recovery time continuity speed of leaking out by decaying greatly and stable a lower value, and fltting speed
The velocity attenuation that leaks out the most greatly is the fastest, i.e. face propulsion speed more large scale mined area time needed for reaching stable state of leaking out is the shortest.Open
Time of adopting is identical, and with wind speed one regularly goaf air leaking speed increases with fltting speed and reduces, and this is also that fltting speed is accelerated
Be conducive to control to leak out, be conducive to the major reason that control spontaneous combustion in goaf and toxic and harmful are gushed out to work surface.
Coal-face constantly advances during exploiting backward, grow out of nothing and progressively extend in goaf, i.e. goaf
Dynamically change is the process that border is progressively moved, regional extent is cumulative, and being therefore analyzed goaf relevant issues can not be by
The border in goaf and regional extent are existing fixing, it is necessary to actual exploitation rate and time for foundation objective description goaf
Process is changed stepwise.
Present embodiment is on the basis of the goaf dynamic evolution numerical simulation set up, in the Dynamic Evolution of goaf
The material parameter broken swollen coefficient of such as inbreak coal petrography, something lost coal thickness, the porosity of inbreak coal petrography, permeability etc. and goaf blast
It is analyzed with air leakage distribution and Changing Pattern, draws following Main Conclusions:
(1) material parameter Changing Pattern in the Dynamic Evolution of goaf
In goaf dynamic evolution and expansion process, broken swollen coefficient, something lost coal thickness, the porosity of inbreak coal petrography, permeability
Descending with the decay of the time of propelling Deng critical material parameter, behind embedment deeper region, goaf, above material parameter could be stablized
To lower value, reach lower state；Coal-face fltting speed is the biggest, and material parameter decay is the fastest；Dynamic Evolution Model energy
Enough present dynamic changing process until reaching stable state comprehensively.Do not consider that the numerical simulation of goaf Dynamic Evolution only describes
It is steady result, it is impossible to describe the dynamic changing process reached before stable state strictly according to the facts.
(2) leak out in the Dynamic Evolution of goaf Changing Pattern
Speed of leaking out advances dynamically change with coal-face, continues gradually decay in time and to stable state, coal work
When face recovery time is identical, coal-face fltting speed leak out the most greatly speed the least and decay faster, beneficially goaf prevent
The toxic and harmfuls such as fire and control gas are gushed out to coal-face.
(3) assume based on goaf moment that the numerical simulation carried out is fixed from advance due to gob edge and region
Realize real goaf dynamic numeric simulation, thus cause based on carried out analog result assumed above at little fltting speed
Lower error is big, and under big fltting speed, error is relatively low；Point out to form, based on goaf transient state, the method assumed the suitableeest simultaneously
By compacting zone under the conditions of coal-face stops to adopt (stop the time of adopting must be considerably larger than simulated domain form required time) for a long time
Outer end is to the relevant issues in this relatively small scope of work surface, and as stopped adopting spontaneous fire during removing frame, gas is gushed out, leaked
Wind rule etc..
Result above shows, stope dynamic evolution numerical simulation can describe gob material parameter (broken swollen coefficient, something lost strictly according to the facts
Coal thickness, the porosity of inbreak coal petrography, permeability) and blast and the dynamic changing process of speed of leaking out, and general existing will adopt
The method for numerical simulation that border, dead zone and region are fixed is all based on the formation of goaf moment, thus can not be the most real
Describe the problems referred to above, steady result can only be presented, be not the goaf dynamic numeric simulation method of true meaning.
Stope dynamic evolution method for numerical simulation and model can truly present gob edge and move, in the propagation process of region
Each variable and the Changing Pattern of parameter.The most broken swollen coefficient of gob material, something lost coal thickness, the porosity of inbreak coal petrography, permeability
Steady statue, the biggest decay of coal-face fltting speed is reached Deng after the most progressively decaying in the Dynamic Evolution of goaf
The fastest；The decay of material parameter and change cause goaf air leaking speed dynamic attenuation in time, and fltting speed gets over large scale mined area
Decay of leaking out is the fastest, and inleakage is the least, is more beneficial to goaf fire prevention and control toxic and harmful is gushed out.Above gob material
Parameter, the dynamic rule of speed of leaking out show goaf carries out the vacation that numerical simulation can not be formed based on goaf moment
If, it is impossible to gob edge and region are fixed in advance, it is necessary to the dynamic mobile and the regional dynamics that realize border increase.
Claims (5)
1. goaf air leaking flow field based on a deformation geometry dynamic numeric simulation method, it is characterised in that including:
According to the goaf formed after exploiting a day with every per day fltting speed, set up initially adopting of coal-face exploitation
Dead zone geometric model: advance the speed for goaf length with every per day fltting speed, controls goaf based on deformation geometry several
What model extends with the back production pusher of coal-face, arranges air intake lane and the return airway of coal-face；
Set Free Transform region, goaf, goaf fixed boundary, goaf moving boundary；
Goaf air leaking flow field is set: goaf is set to porous media flows region, coal-face is set to freedom
Fluid flow region, is to control with coal-face fltting speed, propelling time, y-coordinate by goaf bulking factor
Variable, is the change with the change of goaf bulking factor by the porosity of erecting dense individual props coal petrography, permeability, something lost coal thickness
Amount；The air intake lane of coal-face is blast velocity control border, and the return airway of coal-face is that blast controls border, uses blast
Difference pointwise is distributed along the dividing potential drop of process control goaf external boundary；
By initial goaf geometric model grid division and carry out goaf air leaking Numerical Simulation calculate, including goaf rock
Broken swollen coefficient calculations, the porosity of erecting dense individual props coal petrography, permeability, something lost coal thickness, blast, speed calculation of leaking out；
With every per day fltting speed iterative method goaf geometric model subdivision grid carry out goaf air leaking field again
Numerical simulation calculation, stops adopting line or not reaching to complete goaf when setting is held the clock to stopping adopting line until being advanced into
Leak out Flow Field Numerical Calculation.
Goaf air leaking flow field based on deformation geometry the most according to claim 1 dynamic numeric simulation method, its feature
It is, initial goaf geometric model is set as Free Transform region, described goaf；Open-off cut is set as described fixing
Border, as global coordinate system.
Goaf air leaking flow field based on deformation geometry the most according to claim 1 dynamic numeric simulation method, its feature
Being, described goaf moving boundary includes x direction is moved, y direction is fixing single-way moving border and y direction is fixed, x direction
The Moving and free boundary freely pushed ahead by described single-way moving border.
Goaf air leaking flow field based on deformation geometry the most according to claim 3 dynamic numeric simulation method, its feature
Being, the x direction translational speed on described single-way moving border is coal-face drawing speed, the x on described single-way moving border
Direction displacement uses grid displacement according to fltting speed and to advance time product to specify.
5., according to goaf air leaking flow field based on the deformation geometry dynamic numeric simulation method described in claim 3 or 4, it is special
Levying and be, described single-way moving border includes the side up and down in coal face and goaf.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610187845.5A CN105868472B (en) | 2016-03-29 | 2016-03-29 | A kind of goaf air leaking flow field dynamic numeric simulation method based on deformation geometry |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610187845.5A CN105868472B (en) | 2016-03-29 | 2016-03-29 | A kind of goaf air leaking flow field dynamic numeric simulation method based on deformation geometry |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105868472A true CN105868472A (en) | 2016-08-17 |
| CN105868472B CN105868472B (en) | 2019-03-12 |
Family
ID=56625232
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610187845.5A Expired - Fee Related CN105868472B (en) | 2016-03-29 | 2016-03-29 | A kind of goaf air leaking flow field dynamic numeric simulation method based on deformation geometry |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105868472B (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107862123A (en) * | 2017-11-01 | 2018-03-30 | 重庆工程职业技术学院 | A kind of analysis of goaf air movement resistance and computational methods |
| CN108228952A (en) * | 2017-04-26 | 2018-06-29 | 中国矿业大学 | A kind of electrical anisotropy method for numerical simulation of electrical survey (-ing) goaf fissure zone |
| CN108303158A (en) * | 2018-02-08 | 2018-07-20 | 中国矿业大学 | A kind of computational methods of shallow buried coal seam work face goaf earth's surface inleakage |
| CN109281698A (en) * | 2018-09-29 | 2019-01-29 | 天地（常州）自动化股份有限公司 | Mensuration of Mine Ventilation Resistance data processing method based on relative pressure |
| CN109840370A (en) * | 2019-01-22 | 2019-06-04 | 中国矿业大学(北京) | Stop spontaneous combustion in goaf analogy method and application under the state of adopting |
| CN113700472A (en) * | 2021-08-20 | 2021-11-26 | 中国矿业大学(北京) | Method for determining air leakage direction and measuring air leakage amount of goaf |
| CN115828704A (en) * | 2022-12-22 | 2023-03-21 | 山东省地质矿产勘查开发局八〇一水文地质工程地质大队(山东省地矿工程勘察院) | Rapid prediction method for underground water pollution |
| CN118568804A (en) * | 2024-08-01 | 2024-08-30 | 中国矿业大学 | Three-dimensional dynamic numerical simulation method for spontaneous combustion of goaf coal during working face propulsion |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102828767A (en) * | 2012-07-19 | 2012-12-19 | 大同煤矿集团有限责任公司 | Gob natural fire control method |
| CN103806943A (en) * | 2013-12-13 | 2014-05-21 | 中国神华能源股份有限公司 | Device and method for monitoring air leakage of goaf of coal seam |
| CN105223337A (en) * | 2015-10-22 | 2016-01-06 | 重庆大学 | The three-dimensional simulation method of goaf cracked district form |
-
2016
- 2016-03-29 CN CN201610187845.5A patent/CN105868472B/en not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102828767A (en) * | 2012-07-19 | 2012-12-19 | 大同煤矿集团有限责任公司 | Gob natural fire control method |
| CN103806943A (en) * | 2013-12-13 | 2014-05-21 | 中国神华能源股份有限公司 | Device and method for monitoring air leakage of goaf of coal seam |
| CN105223337A (en) * | 2015-10-22 | 2016-01-06 | 重庆大学 | The three-dimensional simulation method of goaf cracked district form |
Non-Patent Citations (1)
| Title |
|---|
| 王美巧 等: "基于CFD对采空区"三场"的数值模拟", 《煤矿安全》 * |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108228952A (en) * | 2017-04-26 | 2018-06-29 | 中国矿业大学 | A kind of electrical anisotropy method for numerical simulation of electrical survey (-ing) goaf fissure zone |
| CN108228952B (en) * | 2017-04-26 | 2020-08-14 | 中国矿业大学 | Electrical anisotropy numerical simulation method for detecting goaf fracture zone by electrical method |
| CN107862123A (en) * | 2017-11-01 | 2018-03-30 | 重庆工程职业技术学院 | A kind of analysis of goaf air movement resistance and computational methods |
| CN108303158A (en) * | 2018-02-08 | 2018-07-20 | 中国矿业大学 | A kind of computational methods of shallow buried coal seam work face goaf earth's surface inleakage |
| CN108303158B (en) * | 2018-02-08 | 2019-11-29 | 中国矿业大学 | A kind of calculation method of shallow buried coal seam work face goaf earth's surface inleakage |
| CN109281698A (en) * | 2018-09-29 | 2019-01-29 | 天地（常州）自动化股份有限公司 | Mensuration of Mine Ventilation Resistance data processing method based on relative pressure |
| CN109840370A (en) * | 2019-01-22 | 2019-06-04 | 中国矿业大学(北京) | Stop spontaneous combustion in goaf analogy method and application under the state of adopting |
| CN113700472A (en) * | 2021-08-20 | 2021-11-26 | 中国矿业大学(北京) | Method for determining air leakage direction and measuring air leakage amount of goaf |
| CN113700472B (en) * | 2021-08-20 | 2022-12-06 | 中国矿业大学(北京) | Method for determining air leakage direction and measuring air leakage amount of goaf |
| CN115828704A (en) * | 2022-12-22 | 2023-03-21 | 山东省地质矿产勘查开发局八〇一水文地质工程地质大队(山东省地矿工程勘察院) | Rapid prediction method for underground water pollution |
| CN118568804A (en) * | 2024-08-01 | 2024-08-30 | 中国矿业大学 | Three-dimensional dynamic numerical simulation method for spontaneous combustion of goaf coal during working face propulsion |
| CN118568804B (en) * | 2024-08-01 | 2024-10-29 | 中国矿业大学 | Three-dimensional dynamic numerical simulation method for spontaneous combustion of goaf coal during working face propulsion |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105868472B (en) | 2019-03-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105868472A (en) | Gob air leakage flow field dynamic numerical simulation method based on deformation geometry | |
| CN107330194B (en) | Numerical simulation analysis method for influence of shield tunnel excavation on adjacent pile group foundation | |
| CN112943364A (en) | Dynamo and Revit based long tunnel parametric modeling method | |
| CN108194081A (en) | Equipment system of roadway-free coal-pillar-free self-retained roadway mining method | |
| CN108131122A (en) | Improve the CO2 amounts of sealing up for safekeeping and the method for oil recovery factor | |
| CN107133381B (en) | Method for determining length of anchor rod for roadway support | |
| CN103821554B (en) | Based on the boring method for arranging without coal pillar mining Y type ventilation goaf | |
| CN103628894A (en) | Water-rich fault fracture zone collapsed roadway repairing method based on angle-variable umbrella-shaped advancing grouting | |
| CN104632219A (en) | Large mining height gob-side entry retaining method for thick coal seams | |
| CN104632221A (en) | Liquid carbon dioxide blasting induced caving mining method | |
| CN103306713A (en) | Networked reinforcement method for cross-layer drilling of outburst coal seam | |
| CN109505606B (en) | Pre-roof-control mechanized sublevel open stoping subsequent filling mining method | |
| CN102926759A (en) | Long-distance tunneling construction method for big-diameter slurry shield underwater shallow earthing | |
| CN110030013B (en) | Gob-side entry retaining method for three-seam periodic roof cutting self-entry side of transition support area | |
| CN104806243B (en) | A kind of wedge shape switching transition method of open air trestle exploitation | |
| CN103510958A (en) | Gentle dip super-high seam rock cross-cut coal uncovering method | |
| CN104846848A (en) | Method for controlling earth surface water and soil loss of underground layered mineral product mining subsidence area | |
| CN104265294A (en) | Coal pillar-free mining technology for blasting mining face of steeply dipping seam | |
| CN109446602A (en) | A kind of Numerical Experimental Method of surface vertical borehole extraction super high seam gas | |
| 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 | |
| CN115680653A (en) | Hydraulic fracturing and permeability increasing extraction method for coal seam roof of coal roadway driving face | |
| CN102587982A (en) | Method for constructing underground salt cavern gas storage groups | |
| CN116401816A (en) | Horizontal well injection/vertical well in-situ leaching uranium extraction streamline simulation and visualization method and system | |
| CN106401591A (en) | Construction method for small-angle downward-penetrating high-speed grid tunnel | |
| CN106815423A (en) | A kind of W types ventilation goaf spontaneous combustion of remaining coal analogy method |
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: 20190312 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |