CN105427376A - Three-dimensional dynamic visualization method of water inrush of coal seam roof - Google Patents
Three-dimensional dynamic visualization method of water inrush of coal seam roof Download PDFInfo
- Publication number
- CN105427376A CN105427376A CN201510712900.3A CN201510712900A CN105427376A CN 105427376 A CN105427376 A CN 105427376A CN 201510712900 A CN201510712900 A CN 201510712900A CN 105427376 A CN105427376 A CN 105427376A
- Authority
- CN
- China
- Prior art keywords
- water
- roof
- zone
- dimensional dynamic
- bearing zone
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 259
- 239000003245 coal Substances 0.000 title claims abstract description 36
- 238000007794 visualization technique Methods 0.000 title abstract 2
- 238000000034 method Methods 0.000 claims abstract description 67
- 238000005065 mining Methods 0.000 claims abstract description 52
- 238000004088 simulation Methods 0.000 claims abstract description 39
- 239000003673 groundwater Substances 0.000 claims abstract description 34
- 230000008569 process Effects 0.000 claims abstract description 25
- 238000011156 evaluation Methods 0.000 claims abstract description 18
- 210000004027 cell Anatomy 0.000 claims description 23
- 238000009826 distribution Methods 0.000 claims description 22
- 238000005192 partition Methods 0.000 claims description 22
- 238000004519 manufacturing process Methods 0.000 claims description 17
- 230000015572 biosynthetic process Effects 0.000 claims description 10
- 238000000638 solvent extraction Methods 0.000 claims description 9
- 238000012360 testing method Methods 0.000 claims description 9
- 238000004422 calculation algorithm Methods 0.000 claims description 7
- 238000010276 construction Methods 0.000 claims description 6
- 238000001514 detection method Methods 0.000 claims description 6
- 239000002360 explosive Substances 0.000 claims description 6
- 238000009877 rendering Methods 0.000 claims description 6
- 238000005520 cutting process Methods 0.000 claims description 5
- 230000002452 interceptive effect Effects 0.000 claims description 5
- 238000003325 tomography Methods 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 4
- 230000008676 import Effects 0.000 claims description 4
- 230000008595 infiltration Effects 0.000 claims description 4
- 238000001764 infiltration Methods 0.000 claims description 4
- 239000007788 liquid Substances 0.000 claims description 4
- 238000011084 recovery Methods 0.000 claims description 4
- 238000013459 approach Methods 0.000 claims description 3
- 210000005056 cell body Anatomy 0.000 claims description 3
- 230000005465 channeling Effects 0.000 claims description 3
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 230000004927 fusion Effects 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 230000001788 irregular Effects 0.000 claims description 3
- 238000007726 management method Methods 0.000 claims description 3
- 230000002265 prevention Effects 0.000 claims description 3
- 230000002441 reversible effect Effects 0.000 claims description 3
- 230000006641 stabilisation Effects 0.000 claims description 3
- 238000011105 stabilization Methods 0.000 claims description 3
- 238000003860 storage Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 2
- 238000012800 visualization Methods 0.000 abstract description 6
- 238000005516 engineering process Methods 0.000 abstract description 5
- 238000013316 zoning Methods 0.000 abstract 2
- 230000000875 corresponding effect Effects 0.000 description 7
- 238000013461 design Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000000007 visual effect Effects 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 230000009172 bursting Effects 0.000 description 4
- 239000011435 rock Substances 0.000 description 4
- 241000209094 Oryza Species 0.000 description 3
- 235000007164 Oryza sativa Nutrition 0.000 description 3
- 238000005094 computer simulation Methods 0.000 description 3
- 235000009566 rice Nutrition 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- IYLGZMTXKJYONK-ACLXAEORSA-N (12s,15r)-15-hydroxy-11,16-dioxo-15,20-dihydrosenecionan-12-yl acetate Chemical compound O1C(=O)[C@](CC)(O)C[C@@H](C)[C@](C)(OC(C)=O)C(=O)OCC2=CCN3[C@H]2[C@H]1CC3 IYLGZMTXKJYONK-ACLXAEORSA-N 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000009412 basement excavation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- IYLGZMTXKJYONK-UHFFFAOYSA-N ruwenine Natural products O1C(=O)C(CC)(O)CC(C)C(C)(OC(C)=O)C(=O)OCC2=CCN3C2C1CC3 IYLGZMTXKJYONK-UHFFFAOYSA-N 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000012732 spatial analysis Methods 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 239000011800 void material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Remote Sensing (AREA)
- Computer Graphics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a three-dimensional dynamic visualization method of water inrush of a coal seam roof. The method includes: building a data set of a mining area; according to the data set, building a 3D geologic model of the mining area; building a roof water inrush condition comprehensive evaluation zoning graph of the mining area; based on the roof water inrush condition comprehensive evaluation zoning graph, determining prepared mining areas, main tunnels, working surfaces, and water inrush point positions of the coal seam; calculating the water level of each node of a water-bearing layer in each typical period before and after water inrush, and obtaining 3D streamlines of corresponding periods; performing three-dimensional dynamic visualization simulation of stratum fracture; and performing three-dimensional dynamic visualization simulation of a groundwater flow field. The method is firstly based on the "three-map method", groundwater simulation and scientific visualization technology are employed, three-dimensional dynamic temporal-spatial expression of the water inrush of the coal seam roof is realized, and a brand-new platform for understanding and analyzing the process of water inrush of the coal seam roof is provided.
Description
Technical field
The present invention relates to three-dimensional geological modeling technical field, refer to a kind of three dimensional dynamic FEM method based on the roof gushing water of " three figure methods " especially.
Background technology
Along with the mining degree of depth strengthens and lower group coal mining gradually, roof collapse is linked up overlying aquifer and causes top board to gush (dash forward) water disaster occurring or to worsen the example of working surface production environment increasing, and roof gushing water problem is perplex the main flood of North China type coalfield coal industry sustainable development always.Mine water bursting disaster is one of restriction coal in China industrial security the economic exploitation main disasters kind for a long time, be known as colliery " the second killer ", the whole nation 65% mine main, 40% raw coal reserves are threatened by water damage all in various degree, cause casualties and property loss, and impact mining is produced.
Because roof water bursting disaster has sudden, ambiguity and uncertain feature, art methods lacks the expressed intact to water bursting disaster generating process, accurately cannot reflect the space-time dynamic reality that seam mining process and roof gushing water occur, thus become the bottleneck realizing roof water-bursting predicting and visual analyzing thereof.
Summary of the invention
In view of this, the object of the invention is to a kind of three dimensional dynamic FEM method proposing roof gushing water, to improve, the integrality of roof water bursting disaster generating process is expressed.
Based on above-mentioned purpose, the three dimensional dynamic FEM method of roof gushing water provided by the invention comprises the following steps:
Build the data set in mining area, described data set comprises boring, section, tomography, fold, DTM/DEM data, remotely-sensed data, cloud data, water-bearing zone watery, water-resisting layer space structure, Fault water channeling situation, hydrology hole and relevant geologic map, topomap, hydrogeologic map, the comprehensive result map of bailing test, work arrangement planimetric map;
According to described data set, build the 3D geologic model in this mining area;
Build the roof water inrush condition comprehensive evaluation block plan in this mining area;
Based on described roof water inrush condition comprehensive evaluation block plan, determine the plan production zone in coal seam, main tunnel, workplace and projective water point position;
On the basis analyzing this groundwater in mining area Field Characteristics point, adopt GMS groundwater simulation system, parameter subregion is carried out to water-bearing zone, calculate the water level of each node in water-bearing zone each typical period of time before gushing water and after generation gushing water, and obtain the 3D streamline of relevant time period;
Formation is emitted to split and is carried out three dimensional dynamic FEM simulation;
Three dimensional dynamic FEM simulation is carried out to ground water field.
In some embodiments of the invention, the step of the roof water inrush condition comprehensive evaluation block plan in described this mining area of structure comprises:
1) thickness in water-bearing zone, specific capacity, infiltration coefficient, core recovery and consumption of rinsing liquid is calculated, adopt analytical hierarchy process to determine the weighted value that each polynary Geo-informatic Tupu is corresponding in the watery of reflection water-bearing zone, thus form the watery block plan in top board water-bearing zone;
2) according to constructed 3D geologic model, the height H of stratum, water-bearing zone to roof is determined, according to the THICKNESS CALCULATION caving zone height H of working seam
m, water guide Breaking belt height H
1i, thus form top board and emit and split security partition figure
3) top board water-bearing zone watery block plan and top board are emitted split security partition figure and superpose, thus obtain roof water inrush condition comprehensive evaluation block plan.
In some embodiments of the invention, described step 2) comprising:
According to described 3D geologic model and data set, calculate and generate the spatial distribution model of roof, coal seam height, caving zone height, water guide Breaking belt height;
According to constructed 3D geologic model, determine the height H of stratum, water-bearing zone to roof, calculate caving zone height according to the thickness of working seam in conjunction with the experimental formula of water prevention code
Water guide Breaking belt height
Wherein, M is that effectively adopting of coal seam is thick;
If H
1i>H, so this region is emit the explosive area splitting security partitioning figure; If H
1i<H, so this region is emit the comparatively safe district splitting security partitioning figure, thus forms top board and emit and split security partition figure.
In some embodiments of the invention, described step 3) comprising:
Emit at top board in the result splitting security partition figure the region being defined as emitting the comparatively safe district splitting security partitioning figure, watery result according to top board water-bearing zone watery block plan is divided into several ranks, emits the explosive area splitting security partition figure to superpose form roof water inrush condition comprehensive evaluation block plan with top board.
In some embodiments of the invention, described plan production zone meets the following conditions: A) gushing water, B do not occur) have water to enter goaf but without gushing water, C) there is gushing water.
In some embodiments of the invention, also comprise the method for longitudinal profile, construction work face and workplace transverse section, be specially:
Utilizing interactive approach, through intending production zone, with the method based on Z depth value, gathering a series of point set p1, p2 ... pm, and form an open curved surface; Make all stratum in itself and 3D geologic model carry out cutting operation, by reconstructed mesh, form the vertical/horizontal section model of workplace.
In some embodiments of the invention, also comprise:
Adopt the MFC method basis of formation soft environment that .net platform provides, in conjunction with OpenGL shape library, by definition class classCOpenGL:publicCGraBase{}, coordinate system is set, form, rendering parameter, pickup, auxiliary queue and rendering parameter, and COpenGL*p_Graphic is defined in homophony function, complete the structure of roof gushing water three dimensional dynamic FEM scene, to hole again model, stratigraphic model, " three figure ", main tunnel, intend production zone, workplace, and water level, in the scene that the streamline model importing * p_Graphic creates, and complete fusion treatment,
Tree control is adopted to carry out Classification Management to described model;
Reconstruct the Breaking belt set up, cutting obtains its emitting on workplace and splits distribution situation, in conjunction with the stratum branch in tree control, successively with * L [i] (i=1 in geologic model, ..., 6) the stratum grid pointed by carries out boolean and asks friendship, and crossing grid cell occurs Partial Reconstruction;
Corpus--based Method method, interpolation discrete point carry out subdivision respectively in the range of influence of Breaking belt, caving zone: the dot density of Approaching Coal Seam is slightly less than the dot density away from coal seam, and the tendency of inbreak presents arch and upper little lower wide, form a series of irregular broken Block Model, be designated as V, to realize broken simulation in real time; Fissure zone: close to the dot density of caving zone slightly larger than the dot density away from caving zone, form a series of convex hull fracture body Model, be designated as F, and reconstruct correspondingly layer grid cell, the degree that width and the coal seam top covering rockmass in crack are caving is associated;
In .net, define timer, started 1/multiple assignment of nIDEvent by OnTimer () event, in conjunction with workplace strike length L, carry out step-length setting, make step delta t=L/FN, FN is totalframes; Import reconstruct dynamic frame, realize dynamic demonstration function.
In some embodiments of the invention, described formation is emitted and is split the step of carrying out three dimensional dynamic FEM simulation and comprise:
Step a, the grid cell arranged in model V are active unit, and other stratigraphic boundary adjacent with V unit is stabilization element, build tree construction and carry out these unit of storage administration;
Step b, exploitation along with coal seam, according to set step delta t, select the affected relevant a series of cell cube v of t caving zone
i∈ V, i≤n, n are the unit number in V, and inbreak tendency presents arch, to be caving unit as current;
Step c, to utilize based on dynamic (dynamical) Rigid Body Collision detection algorithm, carry out Fast Collision Detection during One-male unit, and analyzed by gravity and cell cube intermolecular forces, acquiring unit speed, can calculate the displacement x of node then
i, Δ y
i, Δ z
i, thus realize dynamically falling in real time;
Affected associative cell body f in steps d, selection fissure zone
j∈ F, j≤m, m are the unit number in F, can fracture constraint by setting up correlation unit body, to form crack;
Step e, current all selected cells
when meeting end condition, if exploitation process length is less than L or not yet reaches projective water point position, then proceed to step b, make t+=Δ t, continue the operation of next frame; Otherwise, emit the process of splitting to terminate.
In some embodiments of the invention, described the step that ground water field carries out three dimensional dynamic FEM simulation to be comprised:
Seam mining starts to be the 1st stage to gushing water, for AB workplace longitudinal profile, and by arranging pRP->alpha=0.6 and the glBlendFunc (GL_SRC_ALPHA in CRunPara*pRP, GL_ONE_MINUS_SRC_ALPHA), water-bearing zone is set to transparent mode;
In simulation mining process, test current by selection cell cube v
i,
If
region, inbreak does not have influence on water-bearing zone, and streamline does not change the flow direction, is the 1st kind of situation, gushing water does not occur;
If v
i∈ Z2 region, water-bearing zone is arrived in inbreak, and water-bearing zone watery is poor, only has a small amount of groundwater flow to enter goaf, is the 2nd kind of situation, though goaf has current to enter, but gushing water does not occur;
If v
i∈ Z2 region, water-bearing zone is arrived in inbreak, and water-bearing zone watery is comparatively strong, then gushing water occurs, a large amount of groundwater flow enters goaf, and water level forms funnel, is the 3rd kind of situation, and gushing water occurs;
Gushing water was the 2nd stage after occurring, and by the water level of typical period of time after groundwater simulation acquisition gushing water, water level is reconstructed into grid surface, and carries out Contours Fill, utilized the differing heights of different color reflection water levels; Afterwards, layout at projective water point said units, adopt reverse Partical trace algorithm, obtain the spatial distribution data of a series of streamline, finally realize roof gushing water three dimensional dynamic FEM function;
The space distribution of the hydrogeologic data flow line in this mining area is utilized to test.
Inventors herein propose and solve " the three figure-bis-predicted method " that roof gushes the quantitative evaluation of (dashing forward) water disaster, namely the watery block plan of top board direct water filling aquifer, roof collapse security partition figure, top board gushes (dashing forward) water condition comprehensive sub-areas figure and coal face is overall and segmentation prediction of engineering inflow, top board direct water filling aquifer adopt front pre-unrestrained program prediction, wherein gush (dashing forward) water condition comprehensive sub-areas figure and formed by watery and inbreak security partition figure complex superposition.Along with the development of scientific visualization technology, Chinese scholars is devoted to groundwater simulation and visual research always, have developed MODFLOW simulation system, utilizes ε-machines method for visualizing to describe fluid, visual and Unsteady Groundwater is intuitively visual based on the stream of Surface.
The invention has the beneficial effects as follows:
1) the present invention is first based on " three figure methods ", in conjunction with groundwater simulation and scientific visualization technology, devises the overall realization flow of roof gushing water three dimensional dynamic FEM.Obtain real data collection from mining area, build through 3D geologic model and " three figure ", seam mining simulated domain and projective water point position are determined in design, by groundwater simulation, three dimensional dynamic FEM pre-service, realize emitting and split and the dynamic similation of ground water field.
2) based on constructed " three figure ", by main tunnel, the various candidate scheme of planning and design intending production zone, workplace and projective water point position, to contingent various actual complex situation in seam mining process truly can be reflected, and show that it is interrelated at 3d space, carry out seam mining and gushing water sunykatuib analysis and prediction, reduce man power and material's cost of enterprise, reduce the risk of actual mining engineering.
3) by boolean's cap of the foundation of roof gushing water three dimensional dynamic FEM scene, Breaking belt and stratum grid and reconstruct breaking of 3D range of influence, three dimensional dynamic FEM pre-service is realized.Simplifying grid operations, reduce the processing time in simulation process, strengthen real-time, adapt to enterprise's hardware environment configuration requirement, emitting the real time dynamic simulation split to lay the foundation for realizing in seam mining simulation process.
4) proposed by the inventionly to split and the Dynamic Simulation Method of ground water field about emitting, provide not only the brand-new platform of a understanding and analyses and prediction roof gushing water process, and can be generalized in the Dynamic Simulation Analysis predicted application of underground engineering disasters, surface collapse, underground water and association area.
Therefore, the present invention based on " three figure methods ", utilizes groundwater simulation and scientific visualization technology first, realizes roof gushing water Three-Dimensional Dynamic spatial and temporal expression, provides a brand-new platform understood and analyze roof gushing water process.
Accompanying drawing explanation
Fig. 1 is the schematic flow sheet of the three dimensional dynamic FEM method of the roof gushing water of the embodiment of the present invention;
Fig. 2 is the mining area boring of the embodiment of the present invention, stratum, the partial data of tomography and relation thereof;
Fig. 3 a, 3b are respectively water-bearing zone, the mining area model of the embodiment of the present invention and comprise part FAULT MODEL;
Fig. 4 a, 4b, 4c and 4d are respectively the spatial distribution model of the roof of the embodiment of the present invention, coal seam height, caving zone height, water guide Breaking belt height;
Fig. 5 a, 5b are respectively the surface chart at caving zone and interface, fissure zone top, Fig. 5 c, 5d are respectively the distribution relation figure of caving zone and K1 aquifer floor elevation, Z2 water-bearing zone top board, and Fig. 5 e, 5f are respectively the distribution relation figure of fissure zone and K1 aquifer floor elevation, Z2 water-bearing zone top board;
Fig. 6 a, 6d are respectively the top board water-bearing zone watery block plan in K1, Z2 water-bearing zone, and the top board that Fig. 6 b, 6e are respectively K1, Z2 water-bearing zone emits and splits security partition figure, and Fig. 6 c, 6f are respectively the roof water inrush condition comprehensive evaluation block plan in K1, Z2 water-bearing zone;
Fig. 7 a, Fig. 7 b are respectively a design proposal of seam mining process simulation and the incidence relation at 3d space thereof;
Fig. 8 is certain period 3D streamline distribution schematic diagram of the embodiment of the present invention;
Fig. 9 a be the embodiment of the present invention workplace on emit and split distribution situation schematic diagram, Fig. 9 b is the schematic diagram that stratum grid carries out that boolean asks friendship, and Fig. 9 c is that the grid cell that stratum grid Partial Reconstruction occurs to intersect hands over schematic diagram;
Figure 10 a is that inbreak does not have influence on water-bearing zone, streamline does not change the ground water field simulation drawing of the flow direction, and Figure 10 b is that water-bearing zone is arrived in inbreak, but the ground water field simulation drawing that water-bearing zone watery is poor, Figure 10 c is that water-bearing zone is arrived in inbreak, but the ground water field simulation drawing that water-bearing zone watery is stronger;
Figure 11 is the three dimensional dynamic FEM schematic diagram of embodiment of the present invention roof gushing water.
Embodiment
For making the object, technical solutions and advantages of the present invention clearly understand, below in conjunction with specific embodiment, and with reference to accompanying drawing, the present invention is described in more detail.
The environment that the present invention realizes comprises hardware environment and software environment:
Hardware environment comprises 3-d seismic exploration equipment, hydrology boring, conventional measurement devices, electronic metering equipment, GPS, laser scanner, power PC, server.
Software environment comprises groundwater simulation system, Geographic Information System .net platform and OpenGL shape library.
With reference to accompanying drawing 1, it is the schematic flow sheet of the three dimensional dynamic FEM method of roof gushing water provided by the invention.As one embodiment of the present of invention, the three dimensional dynamic FEM method of described roof gushing water can comprise:
Step 101: the data set building mining area, described data set comprises boring, section, tomography, fold, DTM/DEM data, remotely-sensed data, cloud data, water-bearing zone watery, water-resisting layer space structure, Fault water channeling situation, hydrology hole and relevant geologic map, topomap, hydrogeologic map, the comprehensive result map of bailing test, work arrangement planimetric map.The relative Repeat in mining area is (minX:1648, minY:659) and (maxX:13227, maxY:12441).
To be positioned at a midwestern mining area, Inner Mongolia Autonomous Region, this mining area is located in the Western Margin in Ordos Basin, belong to Ordos Plateau western edge underground water subregion, the forming and distribution of underground water controls by physical geography and geologic condition, presents the hydrogeological characteristics of the distinctive arid in the Northwest, semiarid region.According to the void formation of water-bearing media, the occurrence condition of underground water and hydraulic property, loose rock class pore water and petroclastic rock pores gap crevice water can be divided into.Be chalk water-bearing zone (K1), chalk aquitard, straight sieve water-bearing zone (Z2), straight sieve aquitard, totally 5 layers, coal seam by mining area stratigraphic division.From this mining area image data, then by digitized processing, build corresponding data set, and stored in a database.Fig. 2 shows this district's boring, stratum, the partial data of tomography and relation thereof.
Step 102: the data set obtained according to step 101, builds the 3D geologic model in this mining area.Preferably, described 3D geologic model mainly comprises water-bearing zone model (as shown in Figure 3 a) and part FAULT MODEL (as shown in Figure 3 b).It should be noted that, Fig. 3 a is followed successively by chalk water-bearing zone (K1), chalk aquitard, straight sieve water-bearing zone (Z2), straight sieve aquitard, coal seam from top to bottom.
Preferably, the way of thinking of multi-source data integration can be adopted to build the 3D geologic model in mining area.The way of thinking of multi-source data integration is mainly passed through data normalization, stress release treatment, inconsistency and uncertainty, in order to keep the border consistance of water-bearing zone, water level, streamline, when water-bearing zone is set up with parameter subregion for border, make the data from multi-source reliable and available.
Step 103: build the top board water-bearing zone watery block plan in this mining area, top board emits and split security partition figure and roof water inrush condition comprehensive evaluation block plan.Namely " three figure " is built, for the selection in later stage simulation mining region lays the foundation.
Particularly, in a preferred embodiment of the present invention, this step 103 mainly comprises the following steps:
Step 301: calculate the thickness in water-bearing zone, specific capacity, infiltration coefficient, core recovery and consumption of rinsing liquid, analytical hierarchy process is adopted to determine the weighted value that each polynary Geo-informatic Tupu is corresponding in the watery of reflection water-bearing zone, form the watery block plan in top board water-bearing zone, as shown in Fig. 6 a and Fig. 6 d.Wherein, Fig. 6 a is the top board water-bearing zone watery block plan in chalk water-bearing zone (K1), and Fig. 6 d is the top board water-bearing zone watery block plan of straight sieve water-bearing zone (Z2).
In an embodiment of the present invention, the watery subregion in top board water-bearing zone have employed the complex superposition of polynary Geo-informatic Tupu, determine the influential geologic agent of top board water-bearing zone watery, setting up the thickness of filled water bearing strata, specific capacity, infiltration coefficient, core recovery, on the basis of consumption of rinsing liquid five Geo-informatic Tupu special topic storehouses and figure, the spatial analysis functions of combining geographic information system, analytical hierarchy process is used to determine the weighted value that each polynary Geo-informatic Tupu is corresponding in the watery of reflection water-bearing zone, superpose after being multiplied with corresponding weighted value according to the numerical value after thematic map normalization and form top board water-bearing zone watery block plan.
Step 302: according to described 3D geologic model and data set, calculate and generate roof (Fig. 4 a), the spatial distribution model of coal seam height (Fig. 4 b), caving zone height (Fig. 4 c), water guide Breaking belt height (Fig. 4 d), the unit in Fig. 4 a-4d is rice.
According to constructed 3D geologic model, determine the height H of stratum, water-bearing zone to roof.Caving zone height H is calculated in conjunction with the experimental formula of water prevention code according to the thickness of working seam
m, water guide Breaking belt height H
1i.In conjunction with the situation in exploitation mining area, be middle hardness rock for rock hardness, then have:
The computing formula of caving zone height is:
The computing formula of water guide Breaking belt height is:
Wherein, M is that effectively adopting of coal seam is thick, and unit is rice.
If H
1i>H, so this region is emit the explosive area splitting security partitioning figure; If H
1i<H, so this region is emit the comparatively safe district splitting security partitioning figure, and formation top board emits and splits security partition figure.
Utilize formula (1), (2), node in 3D geologic model is calculated, and reconstruct this mining area caving zone (Fig. 5 a) and fissure zone (Fig. 5 b) push up the curved surface at interface, and then portray the spatial relationship before they and water-bearing zone.Fig. 5 c is the distribution relation figure of caving zone and K1 aquifer floor elevation, Fig. 5 d is the distribution relation figure of caving zone and Z2 water-bearing zone top board, Fig. 5 e is the distribution relation figure of fissure zone and K1 aquifer floor elevation, Fig. 5 f is the distribution relation figure of fissure zone and Z2 aquifer floor elevation, thus form top board and emit and split security partition figure, as shown in Fig. 6 b and Fig. 6 e.Wherein, Fig. 6 b is that the top board in K1 water-bearing zone emits and splits security partition figure, and Fig. 6 e is that the top board in Z2 water-bearing zone emits and splits security partition figure.
Step 303: respectively by the top board water-bearing zone watery block plan in chalk water-bearing zone (Fig. 6 a) and the top board in chalk water-bearing zone emit and split security partition figure (Fig. 6 b) superposition, the top board water-bearing zone watery block plan (Fig. 6 d) in straight sieve water-bearing zone and the top board in straight sieve water-bearing zone are emitted and splits security partition figure (Fig. 6 e) and superpose, form the roof water inrush condition comprehensive evaluation block plan in chalk water-bearing zone and straight sieve water-bearing zone, as shown in Fig. 6 c and Fig. 6 f.
Particularly, emit at top board in the result splitting security partition figure the region being defined as emitting the comparatively safe district splitting security partitioning figure, watery result according to top board water-bearing zone watery block plan is divided into several ranks, emits the explosive area splitting security partition figure to superpose form roof water inrush condition comprehensive evaluation block plan with top board.
Step 104: based on described roof water inrush condition comprehensive evaluation block plan, determine the plan production zone in coal seam, main tunnel, workplace and projective water point position, described plan production zone meets the following conditions: A) gushing water, B do not occur) have water to enter goaf but without gushing water, C) there is gushing water.
Particularly, as shown in Figure 7a, according to roof water inrush condition comprehensive evaluation block plan, a main tunnel, a plan production zone and a workplace and a projective water point position of predicting can be designed.In the present embodiment, projective water point position is X:5535, Y:3360, Z:1020, AB is the workplace longitudinal profile model set up through intending production zone, and shows its interrelated relation at 3d space, as shown in Figure 7b.Truly can reflect contingent various actual complex situation in seam mining process, and show that it is interrelated at 3d space.
Alternatively, by importing related data based on roof Interactive Design main tunnel circuit or database, then for the stratum be correlated with in coal seam, roadway excavation operation can be carried out, to form the main Roadway model of 3D.By importing shp file or Interactive Design, polygonal region can be formed as plan production zone.
Seam mining and gushing water generating process is observed for the ease of geologist, the invention provides the method for longitudinal profile, construction work face and workplace transverse section, main method is: utilize interactive approach, through intending production zone, with the method based on Z depth value, gather a series of point set p1, p2, ... pm, and form an open curved surface; Make all stratum in itself and 3D geologic model carry out cutting operation, by reconstructed mesh, form the vertical/horizontal section model of workplace.
Step 105: on the basis analyzing this groundwater in mining area Field Characteristics point, adopt GMS groundwater simulation system, parameter subregion is carried out to water-bearing zone, calculates the water level of each node in water-bearing zone each typical period of time before gushing water and after generation gushing water, and obtain the 3D streamline of relevant time period.
Particularly, for design proposal shown in Fig. 7, parameter subregion is carried out to water-bearing zone, and calculate the water level of each node in water-bearing zone each typical period of time before gushing water and after generation gushing water, see table 1 (unit is rice).On this basis, adopt Partical trace algorithm, obtain the 3D streamline of relevant time period, the part streamline data that table 2 is the period described in corresponding table 1, the 3D Streamline traced result of this period as shown in Figure 8.Table 3 essential record period T
itime corresponding water level and streamline data, be convenient to carry out the analysis of roof water-bursting predicting.
Step 106: before realizing three dimensional dynamic FEM, has needed pretreatment operation.Concrete steps are as follows:
Step 601: adopt the MFC method basis of formation soft environment that .net platform provides, in conjunction with OpenGL shape library, by definition class classCOpenGL:publicCGraBase{}, coordinate system, form, rendering parameter, pickup, auxiliary queue and rendering parameter are set, and COpenGL*p_Graphic is defined in homophony function, complete the structure of roof gushing water three dimensional dynamic FEM scene.To hole model, stratigraphic model, " three figure ", main tunnel, intend production zone, workplace and water level, the streamline model imports in the scene that * p_Graphic creates, and completes fusion treatment.
Adopt tree control to carry out Classification Management to these models, table 4 is this mining area part classifying information, and wherein, water level can be set to isoline, and Transparence Display.
Table 4. tree control
Step 602: utilize by calculating and reconstruct the Breaking belt (Fig. 5) set up, can cut and obtain its emitting on workplace AB and split distribution situation (Fig. 9 a).In conjunction with the stratum branch in tree control, successively with * L [i] (i=1 in geologic model, ..., 6) the stratum grid pointed by (L [6] is basic unit) carries out boolean and asks friendship (Fig. 9 b), and there is crossing grid cell (Fig. 9 c) in Partial Reconstruction, wherein, namely the bottom grid cell in Fig. 9 b and 9c constitutes the some effects region of caving zone and stratum L [3].
Step 603: Corpus--based Method method, interpolation discrete point carry out subdivision respectively in the range of influence of Breaking belt.Caving zone: the dot density of Approaching Coal Seam is slightly less than the dot density away from coal seam, and the tendency of inbreak presents arch and upper little lower wide, forms a series of irregular broken Block Model (being designated as V), to realize broken simulation in real time; Fissure zone: close to the dot density of caving zone slightly larger than the dot density away from caving zone, form a series of convex hull fracture body Model (being designated as F), and reconstruct correspondingly layer grid cell, the degree that width and the coal seam top covering rockmass in crack are caving is associated.
Step 604: define timer in .net, starts 1/multiple assignment of nIDEvent by OnTimer () event.In conjunction with workplace strike length L, carry out step-length setting, if step delta t=L/FN, FN are totalframes.Here, step-length is less, and Dynamic Announce effect is more true to nature, but needs to take into account the hardware and software environment supporting Visual Scene.Import reconstruct dynamic frame, realize dynamic demonstration function.
Step 107: formation is emitted to split and carried out three dimensional dynamic FEM simulation, and performing step is as follows:
Step 701: the grid cell arranged in model V is active unit, and other stratigraphic boundary adjacent with V unit is stabilization element, builds tree construction and carrys out these unit of storage administration.Simulation mining coal seam from the North (Fig. 7) in Ω region, if time t=0 is initial value.
Step 702: along with the exploitation in coal seam, according to set step delta t, selects the affected relevant a series of cell cube v of t caving zone
i∈ V (i≤n, n are the unit number in V), inbreak tendency presents arch, to be caving unit as current.
Step 703: utilize based on dynamic (dynamical) Rigid Body Collision detection algorithm, carry out Fast Collision Detection during One-male unit, and analyzed by gravity and cell cube intermolecular forces, acquiring unit speed, can calculate displacement (the Δ x of node then
i, Δ y
i, Δ z
i), thus realize dynamically falling in real time.
Step 704: select affected associative cell body f in fissure zone
j∈ F (j≤m, m are the unit number in F), can fracture constraint by setting up correlation unit body, to form crack.
Step 705: current all selected cells
when meeting end condition, if exploitation process length is less than L or not yet reaches projective water point position (Figure 10 (b)), then proceed to step 702, make t+=Δ t, continue the operation of next frame; Otherwise, emit the process of splitting to terminate (Figure 10 (c)).
Step 108. pair ground water field carries out three dimensional dynamic FEM simulation.Groundwater flow field stimulation is divided into 2 stages, 3 kinds of situations, and concrete steps are as follows:
Step 801. seam mining starts to be the 1st stage to gushing water.For the ease of observing, for AB workplace longitudinal profile, and by arranging pRP->alpha=0.6 and the glBlendFunc (GL_SRC_ALPHA in CRunPara*pRP, GL_ONE_MINUS_SRC_ALPHA), water-bearing zone K1, Z2 are set to transparent mode.
Step 802: in simulation mining process, tests current by selection cell cube v
i,
If
region, inbreak does not have influence on water-bearing zone, and streamline does not change the flow direction, and (Figure 10 a), is the 1st kind of situation, gushing water does not occur;
If v
i∈ Z2 region, water-bearing zone is arrived in inbreak, and water-bearing zone watery is poor, only has a small amount of underground water (streamline) to flow into goaf (Figure 10 b), is the 2nd kind of situation, though goaf has current to enter, but gushing water does not occur;
If v
i∈ Z2 region, water-bearing zone is arrived in inbreak, and water-bearing zone watery is comparatively strong, then gushing water occurs, a large amount of underground water (streamline) flows into goaf, and water level forms funnel (Figure 10 c), is the 3rd kind of situation, and gushing water occurs.
Step 803: gushing water was the 2nd stage after occurring.By the water level of typical period of time after groundwater simulation acquisition gushing water, water level is reconstructed into grid surface, and carries out Contours Fill, utilize the differing heights of different color reflection water levels; Afterwards, layout at projective water point said units, adopt reverse Partical trace algorithm, obtain the spatial distribution data of a series of streamline.In environment set by step 604, finally realize roof gushing water three dimensional dynamic FEM function, as shown in figure 11.
Step 804: utilize the space distribution of local area hydrogeologic data flow line to test, the streamline of acquisition can clearly describe ground water movement rule and effectively disclose geologic feature point to the control action of groundwater flow.
Those of ordinary skill in the field are to be understood that: the discussion of above any embodiment is only exemplary, and not intended to be implies that the scope of the present disclosure (comprising claim) is limited to these examples; Under thinking of the present invention, can combine between the technical characteristic in above embodiment or different embodiment yet, and there are other changes many of different aspect of the present invention as above, they do not provide in details for the sake of simplicity.Therefore, within the spirit and principles in the present invention all, any omission made, amendment, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.
Claims (9)
1. a three dimensional dynamic FEM method for roof gushing water, is characterized in that, comprise the following steps:
Build the data set in mining area, described data set comprises boring, section, tomography, fold, DTM/DEM data, remotely-sensed data, cloud data, water-bearing zone watery, water-resisting layer space structure, Fault water channeling situation, hydrology hole and relevant geologic map, topomap, hydrogeologic map, the comprehensive result map of bailing test, work arrangement planimetric map;
According to described data set, build the 3D geologic model in this mining area;
Build the roof water inrush condition comprehensive evaluation block plan in this mining area;
Based on described roof water inrush condition comprehensive evaluation block plan, determine the plan production zone in coal seam, main tunnel, workplace and projective water point position;
On the basis analyzing this groundwater in mining area Field Characteristics point, adopt GMS groundwater simulation system, parameter subregion is carried out to water-bearing zone, calculate the water level of each node in water-bearing zone each typical period of time before gushing water and after generation gushing water, and obtain the 3D streamline of relevant time period;
Formation is emitted to split and is carried out three dimensional dynamic FEM simulation;
Three dimensional dynamic FEM simulation is carried out to ground water field.
2. the three dimensional dynamic FEM method of roof gushing water according to claim 1, is characterized in that, the step of the roof water inrush condition comprehensive evaluation block plan in described this mining area of structure comprises:
1) thickness in water-bearing zone, specific capacity, infiltration coefficient, core recovery and consumption of rinsing liquid is calculated, adopt analytical hierarchy process to determine the weighted value that each polynary Geo-informatic Tupu is corresponding in the watery of reflection water-bearing zone, thus form the watery block plan in top board water-bearing zone;
2) according to constructed 3D geologic model, the height H of stratum, water-bearing zone to roof is determined, according to the THICKNESS CALCULATION caving zone height H of working seam
m, water guide Breaking belt height H
1i, thus form top board and emit and split security partition figure
3) top board water-bearing zone watery block plan and top board are emitted split security partition figure and superpose, thus obtain roof water inrush condition comprehensive evaluation block plan.
3. the three dimensional dynamic FEM method of roof gushing water according to claim 2, is characterized in that, described step 2) comprising:
According to described 3D geologic model and data set, calculate and generate the spatial distribution model of roof, coal seam height, caving zone height, water guide Breaking belt height;
According to constructed 3D geologic model, determine the height H of stratum, water-bearing zone to roof, calculate caving zone height according to the thickness of working seam in conjunction with the experimental formula of water prevention code
Water guide Breaking belt height
Wherein, M is that effectively adopting of coal seam is thick;
If H
1i>H, so this region is emit the explosive area splitting security partitioning figure; If H
1i<H, so this region is emit the comparatively safe district splitting security partitioning figure, thus forms top board and emit and split security partition figure.
4. the three dimensional dynamic FEM method of roof gushing water according to claim 2, is characterized in that, described step 3) comprising:
Emit at top board in the result splitting security partition figure the region being defined as emitting the comparatively safe district splitting security partitioning figure, watery result according to top board water-bearing zone watery block plan is divided into several ranks, emits the explosive area splitting security partition figure to superpose form roof water inrush condition comprehensive evaluation block plan with top board.
5. the three dimensional dynamic FEM method of roof gushing water according to claim 1, is characterized in that, described plan production zone meets the following conditions: A) gushing water, B do not occur) have water to enter goaf but without gushing water, C) there is gushing water.
6. the three dimensional dynamic FEM method of roof gushing water according to claim 1, is characterized in that, also comprises the method for longitudinal profile, construction work face and workplace transverse section, is specially:
Utilizing interactive approach, through intending production zone, with the method based on Z depth value, gathering a series of point set p1, p2 ... pm, and form an open curved surface; Make all stratum in itself and 3D geologic model carry out cutting operation, by reconstructed mesh, form the vertical/horizontal section model of workplace.
7. the three dimensional dynamic FEM method of roof gushing water according to claim 1, is characterized in that, also comprise:
Adopt the MFC method basis of formation soft environment that .net platform provides, in conjunction with OpenGL shape library, by definition class classCOpenGL:publicCGraBase{}, coordinate system is set, form, rendering parameter, pickup, auxiliary queue and rendering parameter, and COpenGL*p_Graphic is defined in homophony function, complete the structure of roof gushing water three dimensional dynamic FEM scene, to hole again model, stratigraphic model, " three figure ", main tunnel, intend production zone, workplace, and water level, in the scene that the streamline model importing * p_Graphic creates, and complete fusion treatment,
Tree control is adopted to carry out Classification Management to described model;
Reconstruct the Breaking belt set up, cutting obtains its emitting on workplace and splits distribution situation, in conjunction with the stratum branch in tree control, successively with * L [i] (i=1 in geologic model, ..., 6) the stratum grid pointed by carries out boolean and asks friendship, and crossing grid cell occurs Partial Reconstruction;
Corpus--based Method method, interpolation discrete point carry out subdivision respectively in the range of influence of Breaking belt, caving zone: the dot density of Approaching Coal Seam is slightly less than the dot density away from coal seam, and the tendency of inbreak presents arch and upper little lower wide, form a series of irregular broken Block Model, be designated as V, to realize broken simulation in real time; Fissure zone: close to the dot density of caving zone slightly larger than the dot density away from caving zone, form a series of convex hull fracture body Model, be designated as F, and reconstruct correspondingly layer grid cell, the degree that width and the coal seam top covering rockmass in crack are caving is associated;
In .net, define timer, started 1/multiple assignment of nIDEvent by OnTimer () event, in conjunction with workplace strike length L, carry out step-length setting, make step delta t=L/FN, FN is totalframes; Import reconstruct dynamic frame, realize dynamic demonstration function.
8. the three dimensional dynamic FEM method of roof gushing water according to claim 7, is characterized in that, described formation is emitted and split the step of carrying out three dimensional dynamic FEM simulation and comprise:
Step a, the grid cell arranged in model V are active unit, and other stratigraphic boundary adjacent with V unit is stabilization element, build tree construction and carry out these unit of storage administration;
Step b, exploitation along with coal seam, according to set step delta t, select the affected relevant a series of cell cube v of t caving zone
i∈ V, i≤n, n are the unit number in V, and inbreak tendency presents arch, to be caving unit as current;
Step c, to utilize based on dynamic (dynamical) Rigid Body Collision detection algorithm, carry out Fast Collision Detection during One-male unit, and analyzed by gravity and cell cube intermolecular forces, acquiring unit speed, can calculate the displacement x of node then
i, Δ y
i, Δ z
i, thus realize dynamically falling in real time;
Affected associative cell body f in steps d, selection fissure zone
j∈ F, j≤m, m are the unit number in F, can fracture constraint by setting up correlation unit body, to form crack;
Step e, current all selected cells
when meeting end condition, if exploitation process length is less than L or not yet reaches projective water point position, then proceed to step b, make t+=Δ t, continue the operation of next frame; Otherwise, emit the process of splitting to terminate.
9. the three dimensional dynamic FEM method of roof gushing water according to claim 7, is characterized in that, describedly comprises the step that ground water field carries out three dimensional dynamic FEM simulation:
Seam mining starts to be the 1st stage to gushing water, for AB workplace longitudinal profile, and by arranging pRP->alpha=0.6 and the glBlendFunc (GL_SRC_ALPHA in CRunPara*pRP, GL_ONE_MINUS_SRC_ALPHA), water-bearing zone is set to transparent mode;
In simulation mining process, test current by selection cell cube v
i,
If
region, inbreak does not have influence on water-bearing zone, and streamline does not change the flow direction, is the 1st kind of situation, gushing water does not occur;
If v
i∈ Z2 region, water-bearing zone is arrived in inbreak, and water-bearing zone watery is poor, only has a small amount of groundwater flow to enter goaf, is the 2nd kind of situation, though goaf has current to enter, but gushing water does not occur;
If v
i∈ Z2 region, water-bearing zone is arrived in inbreak, and water-bearing zone watery is comparatively strong, then gushing water occurs, a large amount of groundwater flow enters goaf, and water level forms funnel, is the 3rd kind of situation, and gushing water occurs;
Gushing water was the 2nd stage after occurring, and by the water level of typical period of time after groundwater simulation acquisition gushing water, water level is reconstructed into grid surface, and carries out Contours Fill, utilized the differing heights of different color reflection water levels; Afterwards, layout at projective water point said units, adopt reverse Partical trace algorithm, obtain the spatial distribution data of a series of streamline, finally realize roof gushing water three dimensional dynamic FEM function;
The space distribution of the hydrogeologic data flow line in this mining area is utilized to test.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510712900.3A CN105427376B (en) | 2015-10-28 | 2015-10-28 | The three dimensional dynamic FEM method of roof gushing water |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510712900.3A CN105427376B (en) | 2015-10-28 | 2015-10-28 | The three dimensional dynamic FEM method of roof gushing water |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105427376A true CN105427376A (en) | 2016-03-23 |
| CN105427376B CN105427376B (en) | 2017-12-15 |
Family
ID=55505555
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510712900.3A Expired - Fee Related CN105427376B (en) | 2015-10-28 | 2015-10-28 | The three dimensional dynamic FEM method of roof gushing water |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105427376B (en) |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106846474A (en) * | 2016-12-29 | 2017-06-13 | 中国科学院电子学研究所苏州研究院 | WebGIS spatio temporal process simulation methods based on temporal aspect and particIe system |
| CN106875471A (en) * | 2017-01-13 | 2017-06-20 | 山东科技大学 | Coal measures contains or water barrier Visualization Modeling method |
| CN106910139A (en) * | 2017-02-22 | 2017-06-30 | 北京石油化工学院 | A kind of prominent flooding disaster emergency evacuation analogy method in colliery |
| CN107564101A (en) * | 2017-10-30 | 2018-01-09 | 中国煤炭地质总局水文地质局 | A kind of method that visualization of 3 d geological model is established based on AutoCAD |
| CN108086977A (en) * | 2017-12-14 | 2018-05-29 | 中国矿业大学 | A kind of coal seam prevents water-resisting floor coal-mining method for water-bearing layer |
| CN109063230A (en) * | 2018-06-11 | 2018-12-21 | 中国矿业大学(北京) | Different gushing waters or permeable condition go down into a mine the evaluation method of the process of flooding |
| CN109767136A (en) * | 2019-01-17 | 2019-05-17 | 内蒙古上海庙矿业有限责任公司 | Roof water-bearing layer water fathering hazard assessment and the duplexing practice of prediction four |
| CN109871648A (en) * | 2019-03-11 | 2019-06-11 | 山东科技大学 | The construction method of groundwater resources three-dimensional visualization dynamic monitoring structural model |
| CN110163966A (en) * | 2019-05-16 | 2019-08-23 | 中煤科工集团北京华宇工程有限公司 | Establish the method for coal mine three-dimensional geological information model automatically using point cloud data |
| CN110491267A (en) * | 2019-07-10 | 2019-11-22 | 中国矿业大学(北京) | A kind of simulator of roof water damage Disaster mechanism and risk |
| CN110851991A (en) * | 2019-11-18 | 2020-02-28 | 核工业二〇八大队 | Underground water flow numerical simulation method |
| CN112233239A (en) * | 2020-10-29 | 2021-01-15 | 中国科学院地质与地球物理研究所 | Visual evaluation method and system based on dynamic three-dimensional modeling feedback |
| CN117218258A (en) * | 2023-11-08 | 2023-12-12 | 山东大学 | Shield geological and tunnel visualization method, system, medium and equipment |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| SU1617145A1 (en) * | 1989-01-20 | 1990-12-30 | Всесоюзный научно-исследовательский институт горной геомеханики и маркшейдерского дела | Method of reducing water inflow into mine workings |
| CN101660411A (en) * | 2009-05-08 | 2010-03-03 | 中国矿业大学(北京) | Quantitative evaluation method for coal-bed top-plate water burst (inrush) conditions |
| CN101660427A (en) * | 2009-05-18 | 2010-03-03 | 中国矿业大学(北京) | Three-chart two-prediction method for quantitative evaluation and prediction of coal-bed top-plate water burst (inrush) conditions |
| CN101894189A (en) * | 2010-07-14 | 2010-11-24 | 中国矿业大学(北京) | New method for evaluating coal seam bottom water bursting |
| US20140003672A1 (en) * | 2012-06-29 | 2014-01-02 | Weyerhaeuser Nr Company | Method and apparatus for deriving stream network vectors from digital elecation map data |
-
2015
- 2015-10-28 CN CN201510712900.3A patent/CN105427376B/en not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| SU1617145A1 (en) * | 1989-01-20 | 1990-12-30 | Всесоюзный научно-исследовательский институт горной геомеханики и маркшейдерского дела | Method of reducing water inflow into mine workings |
| CN101660411A (en) * | 2009-05-08 | 2010-03-03 | 中国矿业大学(北京) | Quantitative evaluation method for coal-bed top-plate water burst (inrush) conditions |
| CN101660427A (en) * | 2009-05-18 | 2010-03-03 | 中国矿业大学(北京) | Three-chart two-prediction method for quantitative evaluation and prediction of coal-bed top-plate water burst (inrush) conditions |
| CN101894189A (en) * | 2010-07-14 | 2010-11-24 | 中国矿业大学(北京) | New method for evaluating coal seam bottom water bursting |
| US20140003672A1 (en) * | 2012-06-29 | 2014-01-02 | Weyerhaeuser Nr Company | Method and apparatus for deriving stream network vectors from digital elecation map data |
Cited By (21)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106846474A (en) * | 2016-12-29 | 2017-06-13 | 中国科学院电子学研究所苏州研究院 | WebGIS spatio temporal process simulation methods based on temporal aspect and particIe system |
| CN106846474B (en) * | 2016-12-29 | 2020-04-03 | 中国科学院电子学研究所苏州研究院 | WebGIS (Web geographic information System) time-space process simulation method based on time sequence characteristics and particle systems |
| CN106875471B (en) * | 2017-01-13 | 2019-12-06 | 山东科技大学 | Three-dimensional visual modeling method for coal-series water-containing or water-resisting layer |
| CN106875471A (en) * | 2017-01-13 | 2017-06-20 | 山东科技大学 | Coal measures contains or water barrier Visualization Modeling method |
| CN106910139A (en) * | 2017-02-22 | 2017-06-30 | 北京石油化工学院 | A kind of prominent flooding disaster emergency evacuation analogy method in colliery |
| CN107564101B (en) * | 2017-10-30 | 2021-04-30 | 中煤地华盛水文地质勘察有限公司 | Method for establishing visual three-dimensional geological model based on AutoCAD |
| CN107564101A (en) * | 2017-10-30 | 2018-01-09 | 中国煤炭地质总局水文地质局 | A kind of method that visualization of 3 d geological model is established based on AutoCAD |
| CN108086977A (en) * | 2017-12-14 | 2018-05-29 | 中国矿业大学 | A kind of coal seam prevents water-resisting floor coal-mining method for water-bearing layer |
| CN109063230B (en) * | 2018-06-11 | 2023-04-07 | 中国矿业大学(北京) | Method for evaluating mine flooding process under different water inrush or water permeability conditions |
| CN109063230A (en) * | 2018-06-11 | 2018-12-21 | 中国矿业大学(北京) | Different gushing waters or permeable condition go down into a mine the evaluation method of the process of flooding |
| CN109767136A (en) * | 2019-01-17 | 2019-05-17 | 内蒙古上海庙矿业有限责任公司 | Roof water-bearing layer water fathering hazard assessment and the duplexing practice of prediction four |
| CN109871648A (en) * | 2019-03-11 | 2019-06-11 | 山东科技大学 | The construction method of groundwater resources three-dimensional visualization dynamic monitoring structural model |
| CN110163966A (en) * | 2019-05-16 | 2019-08-23 | 中煤科工集团北京华宇工程有限公司 | Establish the method for coal mine three-dimensional geological information model automatically using point cloud data |
| CN110163966B (en) * | 2019-05-16 | 2022-12-09 | 中煤科工集团北京华宇工程有限公司 | Method for automatically establishing three-dimensional geological information model of coal mine by using point cloud data |
| CN110491267B (en) * | 2019-07-10 | 2021-02-26 | 中国矿业大学(北京) | Simulation device for disaster mechanism and danger of water disaster of coal seam roof |
| CN110491267A (en) * | 2019-07-10 | 2019-11-22 | 中国矿业大学(北京) | A kind of simulator of roof water damage Disaster mechanism and risk |
| CN110851991A (en) * | 2019-11-18 | 2020-02-28 | 核工业二〇八大队 | Underground water flow numerical simulation method |
| CN110851991B (en) * | 2019-11-18 | 2023-07-14 | 核工业二〇八大队 | Underground water flow numerical simulation method |
| CN112233239A (en) * | 2020-10-29 | 2021-01-15 | 中国科学院地质与地球物理研究所 | Visual evaluation method and system based on dynamic three-dimensional modeling feedback |
| CN117218258A (en) * | 2023-11-08 | 2023-12-12 | 山东大学 | Shield geological and tunnel visualization method, system, medium and equipment |
| CN117218258B (en) * | 2023-11-08 | 2024-03-22 | 山东大学 | Shield geological and tunnel visualization method, system, medium and equipment |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105427376B (en) | 2017-12-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105427376B (en) | The three dimensional dynamic FEM method of roof gushing water | |
| Xiong et al. | A 3D multi-scale geology modeling method for tunnel engineering risk assessment | |
| CN107503797B (en) | Mine Methane tendency early warning system based on 3D emulation platform | |
| Kaufmann et al. | 3D geological modelling from boreholes, cross-sections and geological maps, application over former natural gas storages in coal mines | |
| EP3526628B1 (en) | Geologic structural model generation | |
| CN105469443A (en) | Method for generating three-dimensional geological map based on geological route (PRB) process double modeling | |
| Guo et al. | Explicit-implicit-integrated 3-D geological modelling approach: A case study of the Xianyan Demolition Volcano (Fujian, China) | |
| Wu et al. | A 3D modeling approach to complex faults with multi-source data | |
| Zhang et al. | A case study on integrated modeling of spatial information of a complex geological body | |
| CN106709988A (en) | Construction method of engineering geological section map | |
| CN106709987B (en) | Dynamic construction method of three-dimensional geological profile model | |
| CN105389451B (en) | A kind of in good time accurate modeling method for unmanned getting working face | |
| Park et al. | Determination of the hydraulic conductivity components using a three-dimensional fracture network model in volcanic rock | |
| CN105184862B (en) | A kind of three dimension strata model dynamic fixing method | |
| CN117197377A (en) | Three-dimensional dynamic visual modeling method for mining area | |
| Zhou et al. | Lithology-based 3d modeling of urban geological attributes and their engineering application: a case study of Guang’an city, SW China | |
| Manevich et al. | Geoecological aspects of stress-strain state modeling results of Leninsky coal deposit (Kuzbass, Russia) | |
| Kaufmann et al. | Reprint of “3D geological modelling from boreholes, cross-sections and geological maps, application over former natural gas storages in coal mines”[Comput. Geosci. 34 (2008) 278–290] | |
| CN111859687A (en) | Mixed geological modeling method and system for depicting geological structure of uranium-bearing sand layer | |
| CN102221707B (en) | Seismic attribute extracting method for geologic structure containing over-thrust fault | |
| Lønøy et al. | Geocellular rendering of cave surveys in paleokarst reservoir models | |
| Liu et al. | From the area to the point-study on the key technology of 3D geological hazard modeling in Three Gorges Reservoir area | |
| Wang et al. | 3D geological model of a tunnel for improving construction measures: a case study | |
| Vincenzi et al. | Parametric and numerical modeling tools to forecast hydrogeological impacts of a tunnel | |
| CN105205863A (en) | Three-dimensional (3D) digital graded basin molding 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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171215 Termination date: 20181028 |