CN108829636B - Prediction method for dynamic load stress intensity in lateral coal seam of goaf - Google Patents

Prediction method for dynamic load stress intensity in lateral coal seam of goaf Download PDF

Info

Publication number
CN108829636B
CN108829636B CN201810360498.0A CN201810360498A CN108829636B CN 108829636 B CN108829636 B CN 108829636B CN 201810360498 A CN201810360498 A CN 201810360498A CN 108829636 B CN108829636 B CN 108829636B
Authority
CN
China
Prior art keywords
dynamic load
load stress
follows
unit
passes
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.)
Active
Application number
CN201810360498.0A
Other languages
Chinese (zh)
Other versions
CN108829636A (en
Inventor
神文龙
王襄禹
陈淼
王春
肖同强
王猛
许磊
柏建彪
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Henan University of Technology
Original Assignee
Henan University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Henan University of Technology filed Critical Henan University of Technology
Priority to CN201810360498.0A priority Critical patent/CN108829636B/en
Publication of CN108829636A publication Critical patent/CN108829636A/en
Application granted granted Critical
Publication of CN108829636B publication Critical patent/CN108829636B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • G06F17/13—Differential equations
    • 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/15—Correlation function computation including computation of convolution operations
    • 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
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06—Power analysis or power optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • Operations Research (AREA)
  • Computer Hardware Design (AREA)
  • Computing Systems (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention discloses a method for predicting dynamic load stress intensity in a lateral coal seam of a goaf, which comprises the following steps: 1. establishing a mechanical model of propagation attenuation of dynamic load stress in an interlayer structure surface, and analyzing to obtain an intensity attenuation rule and a fluctuation direction change rule when the dynamic load stress passes through the interlayer structure surface; 2. establishing a mathematical model of propagation attenuation of dynamic load stress in a rock stratum to obtain an intensity attenuation rule along with the increase of a propagation distance; 3. based on the above mechanical model and mathematical model, establishing an analysis model for analyzing the bearing and force transmission of the roof in the mining influence area, and revealing the transmission attenuation rule when the dynamic load passes through the multilayer laminated roof; 4. and (4) bringing the known model parameters into the transmission attenuation rule, and calculating the distribution rule of the dynamic load stress intensity in the coal seam on one side of the goaf. The method determines the analytic solution of the dynamic load stress intensity in the coal seam on one side of the goaf, considers the influences of the bearing stress, the coal rock property, the structural plane property and the dynamic load stress wave parameter, and can efficiently and accurately predict the dynamic load stress intensity distribution.

Description

Prediction method for dynamic load stress intensity in lateral coal seam of goaf
Technical Field
The invention relates to the technical field of prediction of underground dynamic load stress intensity, in particular to a prediction method of dynamic load stress intensity in a lateral coal seam of a mining area.
Background
Stope roof activity is the root of mine pressure display, when the strength and rigidity of underground engineering structures (roadway surrounding rocks, stope faces, supporting bodies, goaf roofs and the like) are not enough to bear loads, the structures evolve from one mechanical equilibrium state to another equilibrium state, the essence of evolution is that the changes of stress fields, displacement fields and fracture fields are caused by the changes of the structures, the evolution results are mine pressure displays (roof caving, rib caving, heaving, sedimentation, instantaneous large deformation, impact mine pressure, failure of the supporting bodies and the like) of different grades and different types, and the effect of evolution is that the life safety of underground workers is endangered, and the safe and efficient production of mines is restricted.
When dynamic load acts on a coal-series stratum in a pressure-bearing state, coal rock masses in an action area vibrate away from a balance position and move relative to adjacent coal rock masses, acting force given by the adjacent coal rock masses is applied, and meanwhile, the adjacent coal rock masses are also vibrated away from the balance position by reacting force, the vibration is transmitted to surrounding coal rock strata in the form of stress waves, when the dynamic load stress strength is large enough, underground excavation spaces such as adjacent roadways and the like are threatened safely, and coal rock dynamic disasters such as impact mine pressure and the like can occur seriously.
Under the influence of medium properties, interlayer structural planes, incident angles, propagation distances and other factors, the intensity of the stress wave is in an attenuation trend, the spatial and temporal evolution rule of the stress wave in the coal rock is mastered, and the method has important guiding significance for safe and efficient exploitation of underground solid mineral resources.
Disclosure of Invention
Aiming at the defects in the prior art, the invention aims to provide the prediction method of the dynamic load stress intensity in the lateral coal seam of the goaf, which is economical and practical, has high operation efficiency and higher precision.
The invention discloses a method for predicting dynamic load stress intensity in a lateral coal seam of a goaf, which comprises the following steps:
the first step is as follows: based on the transfer theory of stress waves in the layered rock stratum, establishing a mechanical model of propagation attenuation of dynamic load stress in the surface of the interlayer structure, and analyzing to obtain an intensity attenuation rule and a fluctuation direction change rule when the dynamic load stress passes through the surface of the interlayer structure;
the second step is that: establishing a balance differential equation of propagation attenuation of dynamic load stress in a rock stratum based on an absorption theory of stress wave propagation in a solid medium, solving a mathematical model of the balance differential equation, and obtaining an intensity attenuation rule of the dynamic load stress along with the increase of a propagation distance;
the third step: based on the above mechanical model and mathematical model, establishing an analysis model for analyzing the bearing and force transmission of the roof in the mining influence area, and revealing the transmission attenuation rule when the dynamic load passes through the multilayer laminated roof;
the fourth step: and (4) bringing the known model parameters into the transmission attenuation rule, and calculating the distribution rule of the dynamic load stress intensity in the coal seam on one side of the goaf.
Furthermore, in the first step, the intensity attenuation law T when the dynamic load stress passes through the interlayer structural surfaceiComprises the following steps:
Figure BDA0001635837060000021
in the formula:
Tiand the attenuation coefficient when the dynamic load stress passes through the interlayer structural surface of the ith layer, the unit is: 1;
i. the ith layer of the interlayer structural surface and the rock stratum, unit: a layer;
vTPithe vibration speed of the transmitted longitudinal wave when the dynamic load stress passes through the laminated structural surface of the ith layer is as follows: m/s;
vIPithe vibration speed of the incident longitudinal wave when the dynamic load stress passes through the interlayer structure surface of the ith layer is as follows: m/s;
solving for vTPiAnd vIPiThe iterative formula of (a) is:
Figure BDA0001635837060000022
Figure BDA0001635837060000023
in the formula:
vRPithe vibration speed of the reflected longitudinal wave when the dynamic load stress passes through the interlayer structure surface of the ith layer is as follows: m/s;
vRSithe vibration speed of the reflected transverse wave when the dynamic load stress passes through the interlayer structural surface of the ith layer is as follows, unit: m/s;
vTSithe vibration speed of the transmission transverse wave when the dynamic load stress passes through the interlayer structural surface of the ith layer is as follows, unit: m/s;
(iteration step number, unit: 1;
gi1、gi2、gi3、gi4、gi5、gi6、fi1、fi2、fi3、fi4、fi5、fi6、fi7、fi8、fi9、fi10the coefficient of the iterative equation is related to the physical and mechanical parameters of the rock strata on two sides of the structural plane, and the specific operation function is as follows:
Figure BDA0001635837060000031
Figure BDA0001635837060000032
in the formula:
ρidensity of the ith formation, unit: kg/m3；
μiPoisson's ratio of the i-th formation, unit: 1;
knithe normal stiffness of the ith interlayer structural surface is as follows: MPa;
ksithe tangential rigidity of the structural surface between the ith layer is as follows: MPa;
Δ t, iteration time step, unit: s;
βPithe incident angle and the reflection angle of the longitudinal wave when the dynamic load stress passes through the interlayer structure surface of the ith layer, the unit is as follows: (iv) DEG;
βSithe reflection angle of the transverse wave when the dynamic load stress passes through the interlayer structure surface of the ith layer is as follows: (iv) DEG;
βPi+1the transmission angle of the longitudinal wave when the dynamic load stress passes through the interlayer structure surface of the ith layer is as follows: (iv) DEG;
βsi+1the transmission angle of the transverse wave when the dynamic load stress passes through the interlayer structure surface of the ith layer is as follows: (iv) DEG;
CPiand CSiThe propagation speed of longitudinal wave and the propagation speed of transverse wave in the ith stratum are respectively as follows: m/s, can be calculated by the following formula:
Figure BDA0001635837060000043
in the formula:
Eielastic modulus of the i-th formation layer, unit: MPa;
k. stress concentration coefficient, unit: 1;
γ, average volume weight of rock formation, unit: n/m3；
H. Average buried depth of stratum, unit: and m is selected.
Further, in the first step, the fluctuation direction change rule when the dynamic load stress passes through the interlayer structural surface is as follows:
Figure BDA0001635837060000041
furthermore, in the second step, the intensity attenuation law Λ of the dynamic load stress in the rock stratumiComprises the following steps:
Figure BDA0001635837060000042
in the formula:
Λithe intensity attenuation coefficient of the dynamic load stress in the ith stratum is as follows, unit: 1;
ζiand the intensity attenuation factor of the dynamic load stress in the ith stratum, unit: 1;
hithickness of the ith layer, unit: and m is selected.
Further, in the third step, the transmission attenuation law T when the dynamic load passes through the multilayer layered top plate is as follows:
Figure BDA0001635837060000051
in the formula:
t, the transmission attenuation coefficient when the dynamic load stress passes through the n layers of laminated top plates, the unit is: 1;
n, the number of layers of the rock stratum and the interlayer structural surface which are penetrated by the dynamic load stress, unit: and (3) a layer.
Compared with the prior art, the invention has the following advantages and beneficial effects:
according to the method, based on the intensity attenuation rule and the fluctuation direction change rule of the dynamic load stress in the interlayer structural surface and the intensity attenuation rule of the dynamic load stress in the rock stratum, the functional relation between the dynamic load source dynamic load intensity of the top plate and the dynamic load response intensity of any point of the lower coal rock stratum is established, so that the analytical solution of the dynamic load stress intensity in the coal seam on one side of the goaf is determined, the influence of the bearing stress, the coal rock attribute, the structural surface attribute and the dynamic load stress wave parameter is considered in the analytical solution, and the distribution rule of the dynamic load stress intensity in the coal seam on one side of the goaf can be economically, efficiently and accurately predicted.
Drawings
FIG. 1 is a schematic diagram of a physical model of the propagation of a dynamic carrier stress wave in a formation according to the present invention.
FIG. 2 is a schematic diagram of a dynamic carrier stress wave propagation attenuation model in an i-th stratum and an i-th stratum interlayer structural plane in the invention.
FIG. 3 is a cross-sectional view of a bottom layer near a working surface in an embodiment of the invention.
FIG. 4 is a graph showing the distance-attenuation coefficient relationship of a dynamic load stress wave in an embodiment of the present invention.
FIG. 5 is a diagram illustrating the distribution of dynamic stress intensity in an embodiment of the present invention.
In the figure: 1. the dynamic carrier wave propagation method comprises the following steps of 1 st rock stratum, 2 nd rock stratum, 3 rd rock stratum, i-1 th rock stratum, 4 th rock stratum, i th rock stratum, 5 th rock stratum, n-1 th rock stratum, 6 th rock stratum, 7 th 1 st interlaminar structural surface, 8 th 2 nd interlaminar structural surface, 9 th interlaminar structural surface, i-1 th interlaminar structural surface, 10 th interlaminar structural surface, 11 th interlaminar structural surface, i +1 th interlaminar structural surface, 12 th interlaminar structural surface, 13 th interlaminar structural surface, a dynamic carrier wave source, 14 dynamic carrier stress wave disturbance area, 15 th dynamic carrier stress wave disturbance range, 16 th rock stratum of the i-1 th rock stratum, density, Poisson ratio, elastic modulus, longitudinal wave propagation speed and transverse wave propagation speed, 17 th rock stratum density, Poisson ratio, elastic modulus, longitudinal wave propagation speed and transverse wave propagation speed, and 18 th layer dynamic carrier wave propagation distance in the i th rock stratum.
Detailed Description
The present invention will be described in further detail with reference to examples, but the embodiments of the present invention are not limited thereto.
Examples
As shown in fig. 1-2, a method for predicting dynamic load stress intensity in a lateral coal seam of a goaf comprises the following steps:
the first step is as follows: based on the transfer theory of stress waves in the layered rock stratum, establishing a mechanical model of propagation attenuation of dynamic load stress in the surface of the interlayer structure, and analyzing to obtain an intensity attenuation rule and a fluctuation direction change rule when the dynamic load stress passes through the surface of the interlayer structure;
the second step is that: establishing a balance differential equation of propagation attenuation of dynamic load stress in a rock stratum based on an absorption theory of stress wave propagation in a solid medium, solving a mathematical model of the balance differential equation, and obtaining an intensity attenuation rule of the dynamic load stress along with the increase of a propagation distance;
the third step: based on the above mechanical model and mathematical model, establishing an analysis model for analyzing the bearing and force transmission of the roof in the mining influence area, and revealing the transmission attenuation rule when the dynamic load passes through the multilayer laminated roof;
the fourth step: and (4) bringing the known model parameters into the transmission attenuation rule, and calculating the distribution rule of the dynamic load stress intensity in the coal seam on one side of the goaf.
Furthermore, in the first step, the intensity attenuation law T when the dynamic load stress passes through the interlayer structural surfaceiComprises the following steps:
Figure BDA0001635837060000061
in the formula:
Tiand the attenuation coefficient when the dynamic load stress passes through the interlayer structural surface of the ith layer, the unit is: 1;
i. the ith layer of the interlayer structural surface and the rock stratum, unit: a layer;
vTPithe vibration speed of the transmitted longitudinal wave when the dynamic load stress passes through the laminated structural surface of the ith layer is as follows: m/s;
vIPithe vibration speed of the incident longitudinal wave when the dynamic load stress passes through the interlayer structure surface of the ith layer is as follows: m/s;
solving for vTPiAnd vIPiThe iterative formula of (a) is:
Figure BDA0001635837060000071
Figure BDA0001635837060000072
in the formula:
vRPithe vibration speed of the reflected longitudinal wave when the dynamic load stress passes through the interlayer structure surface of the ith layer is as follows: m/s;
vRSithe vibration speed of the reflected transverse wave when the dynamic load stress passes through the interlayer structural surface of the ith layer is as follows, unit: m/s;
vTSithe vibration speed of the transmission transverse wave when the dynamic load stress passes through the interlayer structural surface of the ith layer is as follows, unit: m/s;
(iteration step number, unit: 1;
gi1、gi2、gi3、gi4、gi5、gi6、fi1、fi2、fi3、fi4、fi5、fi6、fi7、fi8、fi9、fi10the coefficient of the iterative equation is related to the physical and mechanical parameters of the rock strata on two sides of the structural plane, and the specific operation function is as follows:
Figure BDA0001635837060000073
Figure BDA0001635837060000081
in the formula:
ρidensity of the ith formation, unit: kg/m3；
μiPoisson's ratio of the i-th formation, unit: 1;
knithe normal stiffness of the ith interlayer structural surface is as follows: MPa;
ksithe tangential rigidity of the structural surface between the ith layer is as follows: MPa;
Δ t, iteration time step, unit: s;
βPithe incident angle and the reflection angle of the longitudinal wave when the dynamic load stress passes through the interlayer structure surface of the ith layer, the unit is as follows: (iv) DEG;
βSidynamic loadThe reflection angle of the transverse wave when the stress passes through the interlayer structure surface of the ith layer is as follows: (iv) DEG;
βPi+1the transmission angle of the longitudinal wave when the dynamic load stress passes through the interlayer structure surface of the ith layer is as follows: (iv) DEG;
βSi+1the transmission angle of the transverse wave when the dynamic load stress passes through the interlayer structure surface of the ith layer is as follows: (iv) DEG;
CPiand CSiThe propagation speed of longitudinal wave and the propagation speed of transverse wave in the ith stratum are respectively as follows: m/s, can be calculated by the following formula:
Figure BDA0001635837060000091
in the formula:
Eielastic modulus of the i-th formation layer, unit: MPa;
k. stress concentration coefficient, unit: 1;
γ, average volume weight of rock formation, unit: n/m3；
H. Average buried depth of stratum, unit: and m is selected.
Further, in the first step, the fluctuation direction change rule when the dynamic load stress passes through the interlayer structural surface is as follows:
Figure BDA0001635837060000092
furthermore, in the second step, the intensity attenuation law Λ of the dynamic load stress in the rock stratumiComprises the following steps:
Figure BDA0001635837060000093
in the formula:
Λithe intensity attenuation coefficient of the dynamic load stress in the ith stratum is as follows, unit: 1;
ζiand the intensity attenuation factor of the dynamic load stress in the ith stratum, unit: 1;
hithickness of the ith layer, unit: and m is selected.
Further, in the third step, the transmission attenuation law T when the dynamic load passes through the multilayer layered top plate is as follows:
Figure BDA0001635837060000101
in the formula:
t, the transmission attenuation coefficient when the dynamic load stress passes through the n layers of laminated top plates, the unit is: 1;
n, the number of layers of the rock stratum and the interlayer structural surface which are penetrated by the dynamic load stress, unit: and (3) a layer.
Based on an engineering example, the method comprises the following specific steps:
as shown in fig. 3, for a stratigraphic section near a certain working face of a certain mine, in order to ensure stability of a surrounding rock for tunneling in a coal body on one side of a goaf, dynamic load stress intensity in the coal layer on one side of the goaf needs to be determined, it is known that a dynamic load stress wave source is located at a thin sandstone layer above the coal layer and close to the edge of the goaf, the dynamic load stress wave is generated by fracture of a cantilever structure of the rock stratum, the intensity of the dynamic load stress is 26.7MPa, the peak vibration speed is 5m/s, the vibration frequency is 50Hz, and the dynamic load stress is propagated to the periphery of the surrounding rock, and the distribution rule when the dynamic load stress is propagated to the middle layer of the coal layer on the side of a lower goaf is solved by applying the models shown in fig. 1 and fig. 2, and the concrete solution can be decomposed into the following steps:
firstly, according to a stratum profile diagram shown in fig. 3, determining that a dynamic load stress wave in the fine sandstone should pass through an interlayer structural surface between the fine sandstone and a thick mudstone, the thick mudstone layer, an interlayer structural surface formed by the thick mudstone and the limestone, the limestone rock layer, an interlayer structural surface formed by the limestone rock and the coal bed, and part of the coal bed rock layer can reach the coal bed on one side of the goaf only after three interlayer structural surfaces and three rock layers are formed.
Secondly, determining the normal stiffness of the three-layer interlayer structural plane to be k respectively according to the geological conditions of the mining engineeringn1＝2500MPa、kn2＝2000MPa、kn31500MPa, and a tangential stiffness ks1＝2500MPa、ks2＝2000MPa、ks31500MPa, the layer thickness of each of the three strata is h1＝33m、h2＝13.5m、h3At a density of p, 3.25m1＝2200kg/m3、ρ2＝2600kg/m3、ρ3＝1400kg/m3The modulus of elasticity is respectively E1＝20000MPa、E2＝25000MPa、E38000MPa, Poisson's ratio P1＝0.26、P2＝0.18、μ30.35, 25000kN/m stratum average volume weight gamma3The buried depth H is 600m, the stress concentration coefficient k is 1, and the vibration speed v of the dynamic load stress wave sourceIP15sin (100 pi t), iteration time step Δ t 10-5s, angle of incidence taken as betaP1The attenuation indexes of dynamic load stress waves in the three layers of rock layers are zeta respectively1＝0.031、ζ2＝0.019、ζ3＝0.038；
Thirdly, the density (rho) of the rock stratum in the second step is measured1、ρ2、ρ3) Elastic modulus (E)1、E2、E3) Poisson's ratio (P)1、P2、P3) And the stress concentration coefficient k, the average volume weight gamma of the stratum and the buried depth H of the stratum are brought into the formula (6) to calculate that the propagation speeds of longitudinal waves in the three rock layers are respectively CP1、CP2、CP3The propagation velocity of the transverse wave is CS1、CS2、CS3；
Fourthly, the propagation velocity C of the longitudinal wave in each rock stratum solved in the third stepP1、CP2、CP3Transverse wave propagation velocity CS1、CS2、CS3And angle of incidence betaP1The incidence angle beta of the longitudinal wave upon incidence on the three-layered interlayer structure can be obtained by substituting the formula (7)P1、βP2、βP3Angle of reflection of transverse wave betaS1、βS2、βS3Angle of transmission of longitudinal wave betaP2、βP3、βP4Transverse wave transmission angle betaS2、βS3、βS4；
The fifth step, the calculation result of the fourth step and the calculation of the third step are carried outAs a result, structural plane normal stiffness (k)n1、kn2、kn3) Tangential height of structural plane (k)s1、ks2、ks3) Iteration time step Δ t, formation density (ρ)1、ρ2、ρ3) Poisson's ratio (P)1、μ2、P3) The iterative coefficient (g) of the attenuation coefficient of the dynamic load stress wave passing through three interlayer structural surfaces can be solved by substituting the equations (4) and (5)11、g12、g13、g14、g15、g16)、
(g21、g22、g23、g24、g25、g26)、(g31、g32、g33、g34、g35、g36)、
(f11、f12、f13、f14、f15、f16、f17、f18、f19、f110)、
(f21、f22、f23、f24、f25、f26、f27、f28、f29、f210)、
(f31、f32、f33、f34、f35、f36、f37、f38、f39、f310)。
Sixthly, calculating the result of the fifth step and the vibration speed v of the dynamic load stress wave source1P1The transmission longitudinal wave vibration speed v at two sides of the three-layer interlayer structural plane can be solved by substituting the formula (2) and the formula (3)TP1、vTP2、vTP3And the incident longitudinal wave vibration velocity v1P1、v1P2、v1P3；
And seventhly, substituting the calculation result of the sixth step into the formula (1) to solve the attenuation coefficient that the dynamic load stress wave passes through the three-layer structural surface as T1、T2、T3；
Eighthly, attenuating indexes zeta of dynamic load stress waves in the three layers of rock layers1、ζ2、ζ3Thickness of rock formation h1、h2、h3And substituting the correlation calculation result of the fourth step into the formula (8) to solve the problem that the attenuation coefficients of the dynamic carrier stress wave in the three-layer rock stratum are respectively lambda1、Λ2、Λ3；
A ninth step of substituting the calculation results of the seventh step and the eighth step into a formula (9) to solve the attenuation coefficient T of the dynamic load stress wave reaching the coal seam on one side of the goaf through the three rock strata and the three interlayer structural surface, wherein the result is shown in FIG. 4;
the distribution law of the dynamic load stress intensity transmitted into the coal seam is obtained by multiplying the intensity of the dynamic load stress wave seismic source by the attenuation coefficient, and is shown in fig. 5.
In summary, through the description of the embodiment, a person skilled in the art can better implement the present solution.

Claims (4)

1. A prediction method for dynamic load stress intensity in a lateral coal seam of a goaf is characterized by comprising the following steps:
the first step is as follows: based on the transfer theory of stress waves in the layered rock stratum, establishing a mechanical model of propagation attenuation of dynamic load stress in the surface of the interlayer structure, and analyzing to obtain an intensity attenuation rule and a fluctuation direction change rule when the dynamic load stress passes through the surface of the interlayer structure;
the second step is that: establishing a balance differential equation of propagation attenuation of dynamic load stress in a rock stratum based on an absorption theory of stress wave propagation in a solid medium, solving a mathematical model of the balance differential equation, and obtaining an intensity attenuation rule of the dynamic load stress along with the increase of a propagation distance;
the third step: based on the above mechanical model and mathematical model, establishing an analysis model for analyzing the bearing and force transmission of the roof in the mining influence area, and revealing the transmission attenuation rule when the dynamic load passes through the multilayer laminated roof;
the fourth step: the known model parameters are brought into the transmission attenuation rule, and the distribution rule of the dynamic load stress intensity in the coal seam on one side of the goaf is calculated;
in the first step, the intensity attenuation law T when the dynamic load stress passes through the interlayer structure surfaceiComprises the following steps:
Figure FDA0003382255690000011
in the formula:
Tiand the attenuation coefficient when the dynamic load stress passes through the interlayer structural surface of the ith layer, the unit is: 1;
i. the ith layer of the interlayer structural surface and the rock stratum, unit: a layer;
vTPithe vibration speed of the transmitted longitudinal wave when the dynamic load stress passes through the laminated structural surface of the ith layer is as follows: m/s;
vIPithe vibration speed of the incident longitudinal wave when the dynamic load stress passes through the interlayer structure surface of the ith layer is as follows: m/s;
solving for vTPiAnd vIPiThe iterative formula of (a) is:
Figure FDA0003382255690000012
Figure FDA0003382255690000013
in the formula:
vRPithe vibration speed of the reflected longitudinal wave when the dynamic load stress passes through the interlayer structure surface of the ith layer is as follows: m/s;
vRSithe vibration speed of the reflected transverse wave when the dynamic load stress passes through the interlayer structural surface of the ith layer is as follows, unit: m/s;
vTSithe vibration speed of the transmission transverse wave when the dynamic load stress passes through the interlayer structural surface of the ith layer is as follows, unit: m/s;
j. iteration step number, unit: 1;
gi1、gi2、gi3、gi4、gi5、gi6、fi1、fi2、fi3、fi4、fi5、fi6、fi7、fi8、fi9、fi10the coefficient of the iterative equation is related to the physical and mechanical parameters of the rock strata on two sides of the structural plane, and the specific operation function is as follows:
Figure FDA0003382255690000021
Figure FDA0003382255690000022
in the formula:
ρidensity of the ith formation, unit: kg/m3；
μiPoisson's ratio of the i-th formation, unit: 1;
knithe normal stiffness of the ith interlayer structural surface is as follows: MPa;
ksithe tangential rigidity of the structural surface between the ith layer is as follows: MPa;
Δ t, iteration time step, unit: s;
βPithe incident angle and the reflection angle of the longitudinal wave when the dynamic load stress passes through the interlayer structure surface of the ith layer, the unit is as follows: (iv) DEG;
βSithe reflection angle of the transverse wave when the dynamic load stress passes through the interlayer structure surface of the ith layer is as follows: (iv) DEG;
βPi+1the transmission angle of the longitudinal wave when the dynamic load stress passes through the interlayer structure surface of the ith layer is as follows: (iv) DEG;
βSi+1the transmission angle of the transverse wave when the dynamic load stress passes through the interlayer structure surface of the ith layer is as follows: (iv) DEG;
CPiand CSiThe propagation speed of longitudinal wave and the propagation speed of transverse wave in the ith stratum are respectively as follows: m/s, can be calculated by the following formula:
Figure FDA0003382255690000031
in the formula:
Eielastic modulus of the i-th formation layer, unit: MPa;
k. stress concentration coefficient, unit: 1;
γ, average volume weight of rock formation, unit: n/m3；
H. Average buried depth of stratum, unit: and m is selected.
2. The method for predicting the intensity of the dynamic load stress in the lateral coal seam of the goaf according to claim 1, wherein in the first step, the fluctuation direction change rule of the dynamic load stress when the dynamic load stress passes through the interlayer structural surface is as follows:
Figure FDA0003382255690000032
in the formula:
βSithe reflection angle of the transverse wave when the dynamic load stress passes through the interlayer structure surface of the ith layer is as follows: (iv) DEG;
CPiand CSiThe propagation speed of longitudinal wave and the propagation speed of transverse wave in the ith stratum are respectively as follows: m/s;
βPithe incident angle and the reflection angle of the longitudinal wave when the dynamic load stress passes through the interlayer structure surface of the ith layer, the unit is as follows: (iv) DEG;
βPi+1the transmission angle of the longitudinal wave when the dynamic load stress passes through the interlayer structure surface of the ith layer is as follows: (iv) DEG;
βSi+1the transmission angle of the transverse wave when the dynamic load stress passes through the interlayer structure surface of the ith layer is as follows: (iv) DEG;
CPi+1and CSi+1The propagation speed of longitudinal wave and the propagation speed of transverse wave in the i +1 th layer of rock stratum respectively are as follows: m/s.
3. The method for predicting the intensity of the dynamic stress in the lateral coal seam of the goaf according to claim 1, wherein in the second step, the intensity attenuation law Λ of the dynamic stress in the coal seam is determinediComprises the following steps:
Figure FDA0003382255690000041
in the formula:
Λithe intensity attenuation coefficient of the dynamic load stress in the ith stratum is as follows, unit: 1;
ζiand the intensity attenuation factor of the dynamic load stress in the ith stratum, unit: 1;
hithickness of the ith layer, unit: m;
βPi+1the transmission angle of the longitudinal wave when the dynamic load stress passes through the interlayer structure surface of the ith layer is as follows: degree.
4. The method for predicting the dynamic load stress intensity in the lateral coal seam of the goaf according to claim 1, wherein in the third step, the transfer attenuation law T when the dynamic load passes through the multilayer layered roof is as follows:
Figure FDA0003382255690000042
in the formula:
t, the transmission attenuation coefficient when the dynamic load stress passes through the n layers of laminated top plates, the unit is: 1;
Tiand the attenuation coefficient when the dynamic load stress passes through the interlayer structural surface of the ith layer, the unit is: 1;
Λithe intensity attenuation coefficient of the dynamic load stress in the ith stratum is as follows, unit: 1;
n, the number of layers of the rock stratum and the interlayer structural surface which are penetrated by the dynamic load stress, unit: and (3) a layer.
CN201810360498.0A 2018-04-20 2018-04-20 Prediction method for dynamic load stress intensity in lateral coal seam of goaf Active CN108829636B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810360498.0A CN108829636B (en) 2018-04-20 2018-04-20 Prediction method for dynamic load stress intensity in lateral coal seam of goaf

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810360498.0A CN108829636B (en) 2018-04-20 2018-04-20 Prediction method for dynamic load stress intensity in lateral coal seam of goaf

Publications (2)

Publication Number Publication Date
CN108829636A CN108829636A (en) 2018-11-16
CN108829636B true CN108829636B (en) 2022-04-12

Family

ID=64155136

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810360498.0A Active CN108829636B (en) 2018-04-20 2018-04-20 Prediction method for dynamic load stress intensity in lateral coal seam of goaf

Country Status (1)

Country Link
CN (1) CN108829636B (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8719736B1 (en) * 2013-03-15 2014-05-06 Synopsys, Inc. Compact and accurate wafer topography proximity effect modeling for full chip simulation
CN105808818A (en) * 2016-01-28 2016-07-27 中煤科工集团唐山研究院有限公司 Method for evaluating foundation stability of coal mining subsidence area
CN107657084A (en) * 2017-09-05 2018-02-02 太原理工大学 A kind of quantitative judgement method of compound residual exploiting field Upward mining feasibility
CN107665285A (en) * 2017-10-31 2018-02-06 河南理工大学 One kind determines to leave coal road rational position analysis method under coal column based on sensitive factor

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103953390B (en) * 2014-05-08 2016-02-10 中国矿业大学 The water-retaining production Tight filling method of control overlying mining rock crack and subsidence

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8719736B1 (en) * 2013-03-15 2014-05-06 Synopsys, Inc. Compact and accurate wafer topography proximity effect modeling for full chip simulation
CN105808818A (en) * 2016-01-28 2016-07-27 中煤科工集团唐山研究院有限公司 Method for evaluating foundation stability of coal mining subsidence area
CN107657084A (en) * 2017-09-05 2018-02-02 太原理工大学 A kind of quantitative judgement method of compound residual exploiting field Upward mining feasibility
CN107665285A (en) * 2017-10-31 2018-02-06 河南理工大学 One kind determines to leave coal road rational position analysis method under coal column based on sensitive factor

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
In Situ Test Study of Characteristics of Coal Mining Dynamic Load;Jiang He等;《Shock and Vibration》;20150218;全文 *
Surrounding Rock Deformation Mechanism and Control Technology for Gob-Side Entry Retaining with Fully Mechanized Gangue Backfilling Mining:A Case Study;Peng Gong等;《Shock and Vibration》;20170514;全文 *
应力波在层状岩体中的传播特性和结构面弹性模量测试;刘单权;《中国优秀博硕士学位论文全文数据库(硕士)》;20170515;全文 *
应力波在层状节理岩体中的透反射与衰减规律;孙冰等;《防灾减灾工程学报》;20151231;第35卷(第6期);全文 *
爆炸应力波在层状节理岩体中传播规律及数值模拟;刘其兵;《中国优秀博硕士学位论文全文数据库(硕士)》;20160415;全文 *
采空区下近距离动压巷道矿压特征及支护模拟研究;黄博;《中国优秀博硕士学位论文全文数据库(博士)》;20170215;全文 *

Also Published As

Publication number Publication date
CN108829636A (en) 2018-11-16

Similar Documents

Publication Publication Date Title
Li et al. Study on three-dimensional dynamic stability of open-pit high slope under blasting vibration
Peng Topical areas of research needs in ground control–a state of the art review on coal mine ground control
Zhu et al. Stress wave propagation across jointed rock mass under dynamic extension and its effect on dynamic response and supporting of underground opening
Simangunsong et al. Effect of bedding plane on prediction blast-induced ground vibration in open pit coal mines
Liang et al. Theoretical models for simulating the failure range and stability of inclined floor strata induced by mining and hydraulic pressure
Hong et al. Experimental and numerical studies of the blast-induced overbreak and underbreak in underground roadways
Deng et al. Dynamic buckling mechanism of pillar rockbursts induced by stress waves
Feng et al. 3D modeling of the influence of a splay fault on controlling the propagation of nonlinear stress waves induced by blast loading
Zhu et al. Overburden movement characteristics of top-coal caving mining in multi-seam areas
Song et al. A dynamic ejection coal burst model for coalmine roadway collapse
Chen et al. Numerical simulation of abnormal roof water-inrush mechanism in mining under unconsolidated aquifer based on overburden dynamic damage
Han et al. Surface subsidence prediction method for coal mines with ultrathick and hard stratum
Chen et al. Numerical investigation of 3D distribution of mining-induced fractures in response to longwall mining
Li et al. Deformation and fracture at floor area and the correlation with main roof breakage in deep longwall mining
Zhao et al. Influence of excavation damaged zone on the dynamic response of circular cavity subjected to transient stress wave
Ma et al. An early warning method for water inrush in Dongjiahe coal mine based on microseismic moment tensor
Liu et al. Rock burst risk in an island longwall coal face by stress field
Wang et al. Process of overburden failure in steeply inclined multi-seam mining: insights from physical modelling
Pivarčiová et al. Design of the construction and research of vibrations and heat transfer of mine workings.
Tan et al. Researches on floor stratum fracturing induced by antiprocedure mining underneath close-distance goaf
CN108829636B (en) Prediction method for dynamic load stress intensity in lateral coal seam of goaf
Wang et al. Study on mechanism of rock burst in horizontal section mining of a steeply inclined extra-thick coal seam
Wang et al. Mechanisms influencing the lateral roof roadway deformation by mining-induced fault population activation: a case study
Luo et al. Cispm–ms: a Tool to predict surface subsidence and to study interactions associated with multi–seam mining Operations
Xue et al. Mechanism and prevention practice of coal burst in deep gradual residual coal pillar area: a case study

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
GR01 Patent grant
GR01 Patent grant