CN107103143A - The Forecasting Methodology of working face top covering rockmass permeability variation under the conditions of mining influence - Google Patents
The Forecasting Methodology of working face top covering rockmass permeability variation under the conditions of mining influence Download PDFInfo
- Publication number
- CN107103143A CN107103143A CN201710307809.2A CN201710307809A CN107103143A CN 107103143 A CN107103143 A CN 107103143A CN 201710307809 A CN201710307809 A CN 201710307809A CN 107103143 A CN107103143 A CN 107103143A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msup
- mfrac
- rock
- 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.)
- Withdrawn
Links
- 230000035699 permeability Effects 0.000 title claims abstract description 54
- 238000005065 mining Methods 0.000 title claims abstract description 44
- 238000000034 method Methods 0.000 title claims abstract description 35
- 239000011435 rock Substances 0.000 claims abstract description 97
- 239000003245 coal Substances 0.000 claims description 50
- 238000004364 calculation method Methods 0.000 claims description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 230000033001 locomotion Effects 0.000 abstract description 15
- 230000007547 defect Effects 0.000 abstract description 2
- 239000008398 formation water Substances 0.000 abstract description 2
- 238000005086 pumping Methods 0.000 abstract 1
- 239000010410 layer Substances 0.000 description 28
- 230000008569 process Effects 0.000 description 7
- 238000011160 research Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 238000006073 displacement reaction Methods 0.000 description 5
- 238000004519 manufacturing process Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 239000002245 particle Substances 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000009412 basement excavation Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000004836 empirical method Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 235000007926 Craterellus fallax Nutrition 0.000 description 1
- 240000007175 Datura inoxia Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005489 elastic deformation Effects 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 230000002427 irreversible effect Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 230000036244 malformation Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 239000012466 permeate Substances 0.000 description 1
- 230000002085 persistent effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 238000012797 qualification Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000002344 surface layer Substances 0.000 description 1
- 238000002834 transmittance 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/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Business, Economics & Management (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Databases & Information Systems (AREA)
- Computer Hardware Design (AREA)
- Algebra (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Software Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Agronomy & Crop Science (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Operations Research (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses the Forecasting Methodology of working face top covering rockmass permeability variation under the conditions of mining influence.Using based on influence function strata movement model, calculate movement of each beded rock mass after seam mining, then strain of the reaction to rock stratum both horizontally and vertically is obtained whereby, by setting up the relation between strain and permeability, " strain porosity permeability " relational model of complete set is obtained, the Permeability Distribution rule of mining overburden is studied.Instant invention overcomes traditional rock mass permeability variation Forecasting Methodology exist defect, and can be to such as gas pumping, flow of formation water the problems such as be further analyzed.
Description
Technical field
The invention belongs to mining technical field, more particularly to working face top covering rockmass permeates under the conditions of mining influence
The Forecasting Methodology of rate change.
Background technology
Mining causes overlying rock and earth's surface to produce depression movement, causes huge to activities such as environment, building, productions
Influence, especially large area, rock stratum caused by high tension coal mining and surface movement is to cause mining collapse disaster and region to become
The root of shape.Therefore, quantification and qualification is carried out to earth's surface and the depression situation of overlying strata not only to protection ground surface environment with building
Build thing safety etc. it is significant, and in advance precognition rock stratum inside displacement feature to mining period water damage, gas
Reach that pre-control also has important value in advance Deng dangerous downhole safety factor of production.Both at home and abroad rock stratum is caused for mining
The research of depression, researcher has been obtained for the achievement compared with horn of plenty, and it is various heavy to be established according to different theoretical foundation
Fall into forecast model, such as theoretical analysis forecast model, empirical method forecast model, and physical analogy model etc..
Influence function method is a kind of very effective ways being predicted to subsidence, compared with other method, the party
Method is a kind of method from empirical method to theoretical model transition, and its theoretical foundation is distribution function, therefore also known as influence function
Method.This method thinks, any one unit spot is often excavated among underground ore bed, can result in certain limit around it
Interior earth's surface produces depression, and further away from the point, subsidence amount is smaller.The final depression of earth's surface certain point is then institute in ore bed
There is the superposition for excavating unit spot to the influence caused by the Ground Point.
To cause after seam mining strata movement with it is disrumpent feelings, and form in overlying strata absciss layer and crack.Long-term theory is ground
Study carefully and show with production practices, overlying rock movement is a very complicated motion process.Pass through field observation research and theory
The transmittance process of depression is analyzed, finds there is a kind of inevitable contact, specific manifestation between subsidence and overlying rock movement
For the two correspondence spatially and in time persistent, during influence function method can be for expression rock depression
Correspondence and continuity.
Influence function method is used for the principle for predicting surface movement and deformation, is equally applicable to grinding to the movement and deformation on stratum
Study carefully.Certain rock stratum layer position under earth's surface as shown in Figure 1 on working face, the excavation of any one cell cube of working seam can cause
Predicted point in the rock stratum is to the movement excavated at cell cube.In figure, the origin of O-X global coordinate systems is arranged on the left of working face
Edge and roof bottom intersection.In the overlying strata of arbitrary height (h be coal seam height), excavate cell cube with it is expected that
Horizontal range between point is nearer, excavates bigger to the influence caused by predicted point.When in working face all cell cubes by one by one
During excavation, the final mean annual increment movement P (x, h) of predicted point is these production units in the folded of the influence caused by this predicted point in overlying strata
Plus.
Traditional rock mass permeability variation Forecasting Methodology, mainly passes through the sampling to live coal (rock) sample, Ran Hou
The experiment of correlation is carried out on special experimental facilities, permeability is obtained.Although this kind of method precision is preferably, there is following lack
Point：1) experiment condition of laboratory loading may be inconsistent with live truth, and the permeability variation reliability of acquisition is low；2)
Poor in timeliness, it is necessary to first drill on a large scale, and coring obtains sample, take longer.
The content of the invention
In order to solve the technical problem that above-mentioned background technology is proposed, it is contemplated that working face overlying under the conditions of mining influence
The Forecasting Methodology of rock mass permeability variation, the defect for overcoming traditional rock mass permeability variation Forecasting Methodology to exist improves permeability
The reliability and real-time of variation prediction.
In order to realize above-mentioned technical purpose, the technical scheme is that：
The Forecasting Methodology of working face top covering rockmass permeability variation, comprises the following steps under the conditions of mining influence：
(1) overlying rock is divided into n-layer by the lithology situation on stratum, each layering entered successively from coal seam to ground apparent bearing
Row 1,2 ..., n level number；
(2) by the 1st layer of rock mass, successively calculate the influence function of the sinking prediction of each layering and move horizontally prediction
Influence function；
(3) by the 1st layer of rock mass, successively influence function is integrated in selected integral domain, each point is obtained
The deflection and the amount of moving horizontally at any point in layer；
(4) calculate each point vertical direction according to the deflection of each point to strain, each point is calculated according to the amount of moving horizontally of each point
Horizontal direction is strained, and is strained further according to each point vertical direction and the face of horizontal direction strain calculation overlying strata is strained；
(5) relation of the face strain of coal and rock porosity and overlying strata is set up, coal and rock permeability and coal and rock hole is set up
The relation of rate, so as to obtain the relation of the face strain of permeability variation and overlying strata, characterizing overlying strata body permeability using the relation becomes
That changes predicts the outcome.
Further, in step (2), the influence function for the sinking prediction being respectively layered and the influence letter for moving horizontally prediction
Number is as follows：
In above formula, fs(x',zi) and fu(x',zi) it is respectively sinking predicted impact function and the influence letter for moving horizontally prediction
Number；X' represents the horizontal range between i-th layer of rock mass future position and underground mining unit spot, S (x', zi-1) represent the i-th -1 layer rock
The sinking that body occurs at x', aiRepresent the subsidence factor of i-th layer of rock mass, RiRepresent the mining influence radius of i-th layer of rock mass, h
Represent working face buried depth, ziRepresent the difference in height between i-th layer of rock mass future position and working face of coal seam；As i=1, S (x',
zi-1)=m, m represent seam mining height.
Further, in step (3), the deflection and the amount of moving horizontally are as follows：
In above formula, S (x, i) and U (x, i) are respectively the deflection and the amount of moving horizontally, d of i-th layer of certain point1iAnd d2iRespectively
For the flex point offset distance of i-th layer of rock mass left and right sides, the origin of O-X global coordinate systems is arranged at working face of coal seam left side edge, x
The coordinate of i-th layer of rock mass future position in O-X global coordinate systems is represented, W represents left side edge and the right side in O-X global coordinate systems
Level interval between lateral edges, selected integral domain is Wc=W-d1i-d2i。
Further, in step (4), the vertical direction strain and horizontal direction strain are as follows：
In above formula, εx(x,z)、εz(x, z) is respectively vertical direction strain and horizontal direction strain, and S, U are respectively deflection
With the amount of moving horizontally, x, z are represented horizontally and vertically respectively.
Further, in step (4), the face strain of the overlying strata is as follows：
εt=εx(x,z)+εz(x,z)+εx(x,y)εz(x,z)。
Further, in step (5), the relation of the face strain of coal and rock porosity and overlying strata is as follows：
In above formula, φ is the coal and rock porosity under mining influence, φ0For the original porosity of coal and rock.
Further, in step (5), the relation of coal and rock permeability and coal and rock porosity is as follows：
In above formula, K is coal and rock permeability, KcFor constant, ∑ s is the porous media mesopore surface area of unit volume.
Further, in step (5), the relation of the face strain of permeability variation and overlying strata is as follows：
In above formula, K0For the original permeability of coal and rock, using coal and rock permeability and original permeability under mining influence
RatioCharacterize permeability variation.
The beneficial effect brought using above-mentioned technical proposal：
The present invention calculates each beded rock mass after seam mining using the strata movement model based on influence function
It is mobile, obtain the strain of rock stratum both horizontally and vertically by being differentiated to displacement whereby, by set up strain and permeability it
Between relation, obtain " strain-porosity-permeability " relational model of complete set, obtain the permeability point of mining overburden
Cloth, reliability is higher, and real-time is good, predicts the outcome and is conducive to analysis gas migration, the rule of flow of formation water.
Brief description of the drawings
Fig. 1 is for influence function method integral and calculating schematic diagram；
Fig. 2 is the basic flow sheet of the present invention；
Fig. 3 is overlying rock layering schematic diagram in the present invention.
Embodiment
Below with reference to accompanying drawing, technical scheme is described in detail.
The Forecasting Methodology of working face top covering rockmass permeability variation under the conditions of mining influence, as shown in Fig. 2 detailed process is such as
Under.
First, by overlying rock carry out it is limited be evenly dividing, while to overlying rock until surface layer carry out 1,2 ...,
N level number, as shown in Figure 3.
According to the order that earth's surface n-th layer is upwardly into since the 1st layer of rock stratum (directly top) to the moving influence of rock stratum
Function is defined.The deformation of rock stratum is still described using the influence function provided in background technology, but will be to therein
Calculating parameter is modified, and the revised influence function formula for calculating rock depression is shown below：
In above formula, x' represents the horizontal range between i-th layer of rock mass future position and underground mining unit spot, S (x', zi-1)
Represent the sinking that the i-th -1 layer rock mass occurs at x', aiRepresent the subsidence factor of i-th layer of rock mass, RiI-th layer of rock mass of expression
Mining influence radius, h represents working face buried depth, ziRepresent the difference in height between i-th layer of rock mass future position and working face of coal seam；When
During i=1, S (x', zi-1)=m, m represent seam mining height.
Above formula expresses a kind of thought of iterative calculation in fact, when calculating i-th layer of rock depression amount, it is necessary to first
The rock depression result of calculation of the i-th -1 layer of acquisition, then and after the characteristic parameter with reference to this layer could be carried out.It should be noted
It is, it is more special when calculating the 1st layer of (i.e. during i=1) rock stratum, need to be by S (x', the z in formulai-1) substituted with m, to ensure
The continuity of calculating.
Similarly calculate the influence function for moving horizontally prediction of the i-th rock stratum：
According to the order that earth's surface n-th layer is upwardly into since the 1st layer of rock stratum successively to moving influence function suitable
Scope be integrated, obtain the deformation result of calculation of each rock stratum.
Selected integral domain is still according to shown in Fig. 1, but it is noted that any rock stratum i deviation of inflection point is mutual not
Identical, therefore, the deviation of inflection point for remembering the rock stratum i of the left and right sides is d1iAnd d2i, then influence function is integrated, obtained
The sinking at any point and the amount of moving horizontally in stratum：
In above formula, the origin of O-X global coordinate systems is arranged at working face of coal seam left side edge, and x represents to sit in the O-X overall situations
The coordinate of i-th layer of rock mass future position in mark system, W represents the water between left side edge and right side edge in O-X global coordinate systems
Flat spacing, selected integral domain is Wc=W-d1i-d2i。
Because the sinking of earth's surface each point under mining influence and the displacement of horizontal direction are different so that between points
In the presence of relative movement, so as to cause earth's surface to deform.In traditional subsidence research, earth's surface deformation be divided into inclination,
Curvature, horizontal strain.But, for the research of overlying strata depression, the distribution situation of overlying strata horizontal strain, vertical strain and face strain
It is more meaningful for engineering problem caused by research processing overlying strata depression strain.The horizontal strain of earth's surface refers to that earth's surface is adjacent at 2 points
The amount of moving horizontally difference and the ratio of horizontal range between this 2 points, reflect adjacent 2 points in unit length level move
Dynamic difference.According to definition, horizontal distortion can be regarded as the first derivative moved horizontally.It can be appreciated that the original of overburden deformation
The principle that reason is deformed with earth's surface is the same, i.e., the deflection of each point and the displacement of horizontal direction be not in overlying strata under mining influence
Together so that between points in the presence of relative movement, so as to cause overlying strata to deform.Therefore, the level of each point should in overlying strata
Change is also that overlying strata move horizontally first derivative to x：
Similarly, overlying strata vertical strain refers to the difference of adjacent 2 points of sinking in overlying strata vertical direction and at this 2 points
Between vertical range ratio, reflect it is adjacent up and down in point-to-point transmission unit height sinking difference.Therefore, according to definition,
The vertical strain of each point can be regarded as the p- h of overlying strata sinking first derivative in overlying strata, therefore the vertical strain of overlying strata each point
For：
According to the feature of subsidence basin principal section：Horizontal displacement is zero on the direction perpendicular to principal section, i.e., vertical
In horizontal strain is not present on principal section direction, therefore under two-dimensional case, the face strain of mining overburden is considered as two-dimensional surface
Strain.Relation with horizontal strain and vertical strain is strained by face, derives that overlying strata face strains formula under two-dimensional case：
εt=εx(x,z)+εz(x,z)+εx(x,y)εz(x,z)
Coal and rock is considered as porous media, and numerous studies show that its deformation is the summation of following two parts deformation.The
A part be medium body distortion, i.e., caused by the deformation of coal and rock skeleton particle deformation, the process be it is reversible,
For elastic deformation process；Part II is malformation, i.e., the causing with respect to the changing of the relative positions on space structure due to skeleton particle
Deformation, the process is typically irreversible.
In the relation that the coal and rock porosity of set-up below is strained with face, the faint body distortion of coal and rock is have ignored,
Then it is analyzed as follows on lower pertinent literature Research foundation：
Use VsRepresent the solid skeletal volume of porous media, Δ VsRepresent that it changes；Use VbRepresent the outward appearance body of porous media
Product, Δ VbRepresent that it changes；Use VpRepresent the pore volume of porous media, Δ VpRepresent that it changes.
Defined by porosity, it is assumed that coal and rock porosity is under initial transition：
When coal and rock by original state to a certain deformation state when, its porosity is：
Stereomutation under three-dimensional case for coal and rock unit volume in bulk strain, i.e. deformation process, and in two dimension
Take the face of coal and rock to strain under situation to strain for face, ignore the faint body distortion of coal and rock, make Δ Vs=0, then above formula be changed into：
As can be seen from the above equation, the function of face strain can be considered as by adopting the porosity of overlying coal and rock.
Coal and rock permeability can change with the change of porosity, so influence coal and rock in gas flowing,
Therefore pass of the permeability with face strain variation can be obtained by the relation of permeability and porosity in Kozeny-Carman equations
System.The Kozeny-Carman equations of relation are between coal and rock permeability and porosity：
In above formula, KcFor constant, 5 are about taken；SpFor the surface area of unit pore volume in porous media, ∑ s is unit volume
Porous media mesopore surface area.
Assuming that in an initial condition, coal and rock permeability is：
New permeability is with the ratio between permeability under reset condition under mining influence：
In view of for the coal and rock of specific structure, during its ess-strain, the total surface area of media particle
It is approximately considered constant, then has ∑ s=∑s s0, i.e.,：
Arrangement obtains the relation formula between permeability variation and face strain：
Above formula is the sign formula of coal and rock permeability variation under mining influence, is the function of coal petrography dignity strain.
The technological thought of above example only to illustrate the invention, it is impossible to which protection scope of the present invention is limited with this, it is every
According to technological thought proposed by the present invention, any change done on the basis of technical scheme each falls within the scope of the present invention
Within.
Claims (8)
1. the Forecasting Methodology of working face top covering rockmass permeability variation under the conditions of mining influence, it is characterised in that including following step
Suddenly：
(1) overlying rock is divided into n-layer by the lithology situation on stratum, 1 is carried out successively to each layering from coal seam to ground apparent bearing,
2 ..., n level number；
(2) by the 1st layer of rock mass, successively calculate the influence function of the sinking prediction of each layering and move horizontally the influence of prediction
Function；
(3) by the 1st layer of rock mass, successively influence function is integrated in selected integral domain, obtained in each layering
The deflection and the amount of moving horizontally at any point；
(4) calculate each point vertical direction according to the deflection of each point to strain, each point level is calculated according to the amount of moving horizontally of each point
Direction is strained, and is strained further according to each point vertical direction and the face of horizontal direction strain calculation overlying strata is strained；
(5) relation of the face strain of coal and rock porosity and overlying strata is set up, coal and rock permeability and coal and rock porosity is set up
Relation, so as to obtain the relation of the face strain of permeability variation and overlying strata, overlying strata body permeability variation is characterized using the relation
Predict the outcome.
2. according to claim 1 under the conditions of mining influence working face top covering rockmass permeability variation Forecasting Methodology, it is special
Levy and be, in step (2), the influence function and the influence function for moving horizontally prediction for the sinking prediction being respectively layered are as follows：
<mrow>
<msub>
<mi>f</mi>
<mi>s</mi>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>x</mi>
<mo>&prime;</mo>
</msup>
<mo>,</mo>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mi>S</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>x</mi>
<mo>&prime;</mo>
</msup>
<mo>,</mo>
<msub>
<mi>z</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
</mrow>
<msub>
<mi>R</mi>
<mi>i</mi>
</msub>
</mfrac>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>&pi;</mi>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msup>
<mi>x</mi>
<mo>&prime;</mo>
</msup>
<msub>
<mi>R</mi>
<mi>i</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msup>
</mrow>
<mrow>
<msub>
<mi>f</mi>
<mi>u</mi>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>x</mi>
<mo>&prime;</mo>
</msup>
<mo>,</mo>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>-</mo>
<mn>2</mn>
<mi>&pi;</mi>
<mfrac>
<mrow>
<mi>S</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>x</mi>
<mo>&prime;</mo>
</msup>
<mo>,</mo>
<msub>
<mi>z</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mi>n</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>R</mi>
<mi>i</mi>
</msub>
<mo>&CenterDot;</mo>
<mi>h</mi>
</mrow>
</mfrac>
<msup>
<mi>x</mi>
<mo>&prime;</mo>
</msup>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>&pi;</mi>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msup>
<mi>x</mi>
<mo>&prime;</mo>
</msup>
<msub>
<mi>R</mi>
<mi>i</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msup>
</mrow>
In above formula, fs(x',zi) and fu(x',zi) it is respectively sinking predicted impact function and the influence function for moving horizontally prediction；
X' represents the horizontal range between i-th layer of rock mass future position and underground mining unit spot, S (x', zi-1) represent the i-th -1 layer rock mass
The sinking occurred at x', aiRepresent the subsidence factor of i-th layer of rock mass, RiRepresent the mining influence radius of i-th layer of rock mass, h tables
Show working face buried depth, ziRepresent the difference in height between i-th layer of rock mass future position and working face of coal seam；As i=1, S (x', zi-1)
=m, m represent seam mining height.
3. according to claim 2 under the conditions of mining influence working face top covering rockmass permeability variation Forecasting Methodology, it is special
Levy and be, in step (3), the deflection and the amount of moving horizontally are as follows：
<mrow>
<mi>S</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mo>&Integral;</mo>
<mrow>
<msub>
<mi>d</mi>
<mrow>
<mn>1</mn>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>x</mi>
</mrow>
<mrow>
<mi>W</mi>
<mo>-</mo>
<msub>
<mi>d</mi>
<mrow>
<mn>2</mn>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>x</mi>
</mrow>
</msubsup>
<msub>
<mi>f</mi>
<mi>s</mi>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>x</mi>
<mo>&prime;</mo>
</msup>
<mo>,</mo>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<msup>
<mi>dx</mi>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<mfrac>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<msub>
<mi>R</mi>
<mi>i</mi>
</msub>
</mfrac>
<msubsup>
<mo>&Integral;</mo>
<mrow>
<msub>
<mi>d</mi>
<mrow>
<mn>1</mn>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>x</mi>
</mrow>
<mrow>
<mi>W</mi>
<mo>-</mo>
<msub>
<mi>d</mi>
<mrow>
<mn>2</mn>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>x</mi>
</mrow>
</msubsup>
<mi>S</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>x</mi>
<mo>&prime;</mo>
</msup>
<mo>,</mo>
<msub>
<mi>z</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>&pi;</mi>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msup>
<mi>x</mi>
<mo>&prime;</mo>
</msup>
<msub>
<mi>R</mi>
<mi>i</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msup>
<msup>
<mi>dx</mi>
<mo>&prime;</mo>
</msup>
</mrow>
<mrow>
<mi>U</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mo>&Integral;</mo>
<mrow>
<msub>
<mi>d</mi>
<mrow>
<mn>1</mn>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>x</mi>
</mrow>
<mrow>
<mi>W</mi>
<mo>-</mo>
<msub>
<mi>d</mi>
<mrow>
<mn>2</mn>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>x</mi>
</mrow>
</msubsup>
<msub>
<mi>f</mi>
<mi>u</mi>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>x</mi>
<mo>&prime;</mo>
</msup>
<mo>,</mo>
<msub>
<mi>z</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<msup>
<mi>dx</mi>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<mo>-</mo>
<mn>2</mn>
<mi>&pi;</mi>
<mfrac>
<mrow>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>&CenterDot;</mo>
<mi>n</mi>
</mrow>
<mrow>
<msub>
<mi>R</mi>
<mi>i</mi>
</msub>
<mo>&CenterDot;</mo>
<mi>h</mi>
</mrow>
</mfrac>
<msubsup>
<mo>&Integral;</mo>
<mrow>
<msub>
<mi>d</mi>
<mrow>
<mn>1</mn>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>x</mi>
</mrow>
<mrow>
<mi>W</mi>
<mo>-</mo>
<msub>
<mi>d</mi>
<mrow>
<mn>2</mn>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>x</mi>
</mrow>
</msubsup>
<mi>S</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>x</mi>
<mo>&prime;</mo>
</msup>
<mo>,</mo>
<msub>
<mi>z</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msup>
<mi>x</mi>
<mo>&prime;</mo>
</msup>
<mo>&CenterDot;</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>&pi;</mi>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msup>
<mi>x</mi>
<mo>&prime;</mo>
</msup>
<msub>
<mi>R</mi>
<mi>i</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msup>
<msup>
<mi>dx</mi>
<mo>&prime;</mo>
</msup>
</mrow>
In above formula, S (x, i) and U (x, i) are respectively the deflection and the amount of moving horizontally, d of i-th layer of certain point1iAnd d2iRespectively i-th
The flex point offset distance of the layer rock mass left and right sides, the origin of O-X global coordinate systems is arranged at working face of coal seam left side edge, and x is represented
The coordinate of i-th layer of rock mass future position in O-X global coordinate systems, W represents left side edge and right edge in O-X global coordinate systems
Level interval between edge, selected integral domain is Wc=W-d1i-d2i。
4. according to claim 1 under the conditions of mining influence working face top covering rockmass permeability variation Forecasting Methodology, it is special
Levy and be, in step (4), the vertical direction strain and horizontal direction strain are as follows：
<mrow>
<msub>
<mi>&epsiv;</mi>
<mi>x</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mi>d</mi>
<mi>S</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>d</mi>
<mi>x</mi>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>&epsiv;</mi>
<mi>z</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mi>d</mi>
<mi>U</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>d</mi>
<mi>z</mi>
</mrow>
</mfrac>
</mrow>
In above formula, εx(x,z)、εz(x, z) is respectively vertical direction strain and horizontal direction strain, and S, U are respectively deflection and water
Momentum is translated, x, z are represented horizontally and vertically respectively.
5. according to claim 4 under the conditions of mining influence working face top covering rockmass permeability variation Forecasting Methodology, it is special
Levy and be, in step (4), the face strain of the overlying strata is as follows：
εt=εx(x,z)+εz(x,z)+εx(x,y)εz(x,z)。
6. according to claim 5 under the conditions of mining influence working face top covering rockmass permeability variation Forecasting Methodology, it is special
Levy and be, in step (5), the relation of the face strain of coal and rock porosity and overlying strata is as follows：
<mrow>
<mi>&phi;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&phi;</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&epsiv;</mi>
<mi>t</mi>
</msub>
</mrow>
<mrow>
<mn>1</mn>
<mo>+</mo>
<msub>
<mi>&epsiv;</mi>
<mi>t</mi>
</msub>
</mrow>
</mfrac>
</mrow>
In above formula, φ is the coal and rock porosity under mining influence, φ0For the original porosity of coal and rock.
7. according to claim 6 under the conditions of mining influence working face top covering rockmass permeability variation Forecasting Methodology, it is special
Levy and be, in step (5), the relation of coal and rock permeability and coal and rock porosity is as follows：
<mrow>
<mi>K</mi>
<mo>=</mo>
<mfrac>
<msup>
<mi>&phi;</mi>
<mn>3</mn>
</msup>
<mrow>
<msub>
<mi>K</mi>
<mi>c</mi>
</msub>
<msup>
<mi>&Sigma;s</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
In above formula, K is coal and rock permeability, KcFor constant, ∑ s is the porous media mesopore surface area of unit volume.
8. according to claim 7 under the conditions of mining influence working face top covering rockmass permeability variation Forecasting Methodology, it is special
Levy and be, in step (5), the relation of the face strain of permeability variation and overlying strata is as follows：
<mrow>
<mfrac>
<mi>K</mi>
<msub>
<mi>K</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<msub>
<mi>&epsiv;</mi>
<mi>t</mi>
</msub>
<msub>
<mi>&phi;</mi>
<mn>0</mn>
</msub>
</mfrac>
</mrow>
<mrow>
<mn>1</mn>
<mo>+</mo>
<msub>
<mi>&epsiv;</mi>
<mi>t</mi>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
</mrow>
In above formula, K0For the original permeability of coal and rock, using the ratio of coal and rock permeability and original permeability under mining influenceCharacterize permeability variation.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710307809.2A CN107103143A (en) | 2017-05-04 | 2017-05-04 | The Forecasting Methodology of working face top covering rockmass permeability variation under the conditions of mining influence |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710307809.2A CN107103143A (en) | 2017-05-04 | 2017-05-04 | The Forecasting Methodology of working face top covering rockmass permeability variation under the conditions of mining influence |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107103143A true CN107103143A (en) | 2017-08-29 |
Family
ID=59656664
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710307809.2A Withdrawn CN107103143A (en) | 2017-05-04 | 2017-05-04 | The Forecasting Methodology of working face top covering rockmass permeability variation under the conditions of mining influence |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107103143A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108844818A (en) * | 2018-04-28 | 2018-11-20 | 辽宁工程技术大学 | A kind of method of determining rock in load any time porosity |
| CN111898300A (en) * | 2020-07-31 | 2020-11-06 | 重庆大学 | Slope rock deformation prediction method and system under mining action and storage medium |
| CN113435014A (en) * | 2021-06-04 | 2021-09-24 | 华北水利水电大学 | Dynamic prediction method for mining overburden rock movement deformation |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102708278A (en) * | 2012-04-09 | 2012-10-03 | 北方工业大学 | Method for predicating surface deformation under combined mining effects |
| CN104462654A (en) * | 2014-11-11 | 2015-03-25 | 中国矿业大学 | Shallow burial coal mining earth surface interpenetrated crack distribution and air leakage characteristic judgment method |
-
2017
- 2017-05-04 CN CN201710307809.2A patent/CN107103143A/en not_active Withdrawn
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102708278A (en) * | 2012-04-09 | 2012-10-03 | 北方工业大学 | Method for predicating surface deformation under combined mining effects |
| CN104462654A (en) * | 2014-11-11 | 2015-03-25 | 中国矿业大学 | Shallow burial coal mining earth surface interpenetrated crack distribution and air leakage characteristic judgment method |
Non-Patent Citations (2)
| Title |
|---|
| LUO YI等: "Enhanced subsurface subsidence prediction model that considers overburden stratification", 《INTERNATIONAL JOURNAL OF MINNING ENGINEERING》 * |
| 陈海栋: "保护层开采过程中卸载煤体损伤及渗透性演化特征研究", 《中国博士学位论文全文数据库 工程科技I辑》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108844818A (en) * | 2018-04-28 | 2018-11-20 | 辽宁工程技术大学 | A kind of method of determining rock in load any time porosity |
| CN111898300A (en) * | 2020-07-31 | 2020-11-06 | 重庆大学 | Slope rock deformation prediction method and system under mining action and storage medium |
| CN111898300B (en) * | 2020-07-31 | 2022-08-02 | 重庆大学 | Slope rock deformation prediction method and system under mining action and storage medium |
| CN113435014A (en) * | 2021-06-04 | 2021-09-24 | 华北水利水电大学 | Dynamic prediction method for mining overburden rock movement deformation |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Huang et al. | Influence of filling-drawdown cycles of the Three Gorges reservoir on deformation and failure behaviors of anaclinal rock slopes in the Wu Gorge | |
| Xu et al. | Discrete element modeling of strata and surface movement induced by mining under open-pit final slope | |
| Lin et al. | Large deformation analysis of a high steep slope relating to the Laxiwa Reservoir, China | |
| Zhang et al. | Time-dependent subsidence prediction model and influence factor analysis for underground gas storages in bedded salt formations | |
| CN103791884B (en) | A kind of research method of the integral law of coal mining overlying strata and surface movement and deformation | |
| CN103982179A (en) | Paleopressure quantitative inversion detection method of oil reservoir | |
| CN104713738B (en) | Intelligent two-way multi-angle overturn model test device and test method | |
| Do et al. | Verifying discontinuous deformation analysis simulations of the jointed rock mass behavior of shallow twin mountain tunnels | |
| Li et al. | Impact of longwall mining on groundwater above the longwall panel in shallow coal seams | |
| Aksoy et al. | Estimation of limit angle using laminated displacement discontinuity analysis in the Soma coal field, Western Turkey | |
| CN107103143A (en) | The Forecasting Methodology of working face top covering rockmass permeability variation under the conditions of mining influence | |
| CN103425877A (en) | Geological disaster prediction system and method based on porous media fluid-structure interaction model | |
| CN101942992A (en) | Method for predicting pore pressure of regional high-pressure saltwater layer by utilizing curvature of face of geologic structure | |
| Griffith et al. | 3D topographic stress perturbations and implications for ground control in underground coal mines | |
| CN107013209A (en) | The Forecasting Methodology of working face overlying stratiform Rock Displacement Movement under the conditions of mining influence | |
| Odonne et al. | Abnormal reverse faulting above a depleting reservoir | |
| Su et al. | Stability prediction and optimal angle of high slope in open-pit mine based on two-dimension limit equilibrium method and three-dimension numerical simulation | |
| YongYao et al. | Comprehensive critical mechanical model of covered karst collapse under the effects of positive and negative pressure | |
| CN106372295B (en) | Method for determining safe horizontal distance between shield and karst cave in sandy karst stratum | |
| CN112733412B (en) | Speed equivalent characterization method for hydrodynamic force action landslide motion mechanism research | |
| Wu et al. | Development characteristics of the fault system and its control on basin structure, Bodong Sag, East China | |
| ASLAN et al. | Wide-area ground deformation monitoring in geothermal fields in western Turkey | |
| Chen et al. | Beneficial use of deep ordovician limestone water from mine safety dewatering at the Xinglongzhuang coal mine, north China | |
| Ding et al. | Numerical simulations of how staged dewatering and mining influence surface subsidence | |
| Kong et al. | A review of the mass-loss-induced non-Darcy seepage and seepage suffosion in the fractured zone: the concepts discrimination and connection |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB02 | Change of applicant information |
Address after: 221116 No. 1 University Road, copper mountain, Jiangsu, Xuzhou Applicant after: China University of Mining & Technology Address before: 221116 Xuzhou University Road, Jiangsu, No. 1 Applicant before: China University of Mining & Technology |
|
| CB02 | Change of applicant information | ||
| WW01 | Invention patent application withdrawn after publication |
Application publication date: 20170829 |
|
| WW01 | Invention patent application withdrawn after publication |