Mine earthquake master control horizon determination method

Abstract

The invention provides a mine earthquake master control horizon determination method, which comprises the following steps of: establishing a mine area stratum model, and determining a plurality of key layers which are positioned above a fissure zone in the mine area stratum model and have mine earthquake risks as mine earthquake control layers according to a key layer theory; main control layer analysis is carried out on all the mine earthquake control layers, a step-by-step determination and governance analysis principle is formulated for the mine earthquake control layers, and multi-stage mine earthquake main control layers are determined from all the mine earthquake control layers; and performing logging correction analysis, and dynamically correcting the position of the target rock stratum of the multi-stage mine earthquake main control layer. According to the method, different mine earthquake disaster-causing conditions are analyzed, mine earthquake main control layers under the scene that mine earthquakes generate disaster-causing influence on the ground and under the scene that mine earthquakes generate disaster-causing influence on an underground mining working face are respectively determined, and the target rock layer position of the mine earthquake main control layers is dynamically corrected through logging analysis; and it is guaranteed that the target layer of the mine earthquake main control layer can be subjected to construction treatment.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of coal mine disaster control, and particularly relates to a mine shock main control horizon determination method. BACKGROUND

[0002] As an important basic energy in China, in recent years, although the proportion of coal consumption in the total primary energy consumption in China has been continuously decreasing, due to the continuous increase in the total energy consumption, the absolute amount of coal consumption has been continuously increasing. For the area where the main mining coal seam has a relatively large buried depth, the main mining coal seam has a buried depth of about 600-750m, and the high static load caused by the low thick hard roof and the strong dynamic load caused by the breakage are easy to induce the rock burst, and the high position thick hard roof breakage releases a large amount of elastic energy, which is easy to cause the frequent occurrence of ground mine shock. According to the actual monitoring, a large number of mine shock events are distributed above 200m of the roof, which is beyond the existing technical control range in the underground mine, and the rock burst and mine shock have become the main problem restricting the safe mining of coal in the area.

[0003] Moreover, the mine shock prevention and control has an inaccurate control target horizon, and the mine shock problem is often controlled after the problem occurs. Due to the unclear control target horizon, the main control horizon of the mine shock cannot be reasonably determined, the control effect is not obvious, and the control range has to be as large as possible to cover the control target horizon, and the control is blind. SUMMARY

[0004] The present application aims at at least one of the problems in the related art.

[0005] To achieve the above-mentioned purpose, the present application provides a mine shock main control horizon determination method, which comprises the following steps:

[0006] S1. According to the three-dimensional seismic exploration, the drilling core and the microseismic positioning analysis, a mine area stratum model is established, a working face is selected, and the key layer theory is used to determine the key layer distributed above the fracture zone, a plurality of key layers located above the fracture zone in the mine area stratum model and having a mine shock risk are determined as mine shock control layers;

[0007] S2. The main control layer analysis is performed on all the mine shock control layers above the fracture zone, the step-by-step determination and control analysis principle is formulated for the mine shock control layer, the mine shock control layer which cannot cause disaster after being controlled is removed, and a plurality of mine shock main control layers are determined from all the mine shock control layers;

[0008] S3. The well logging correction analysis is performed, and the horizon of the target rock layer of the plurality of mine shock main control layers is dynamically corrected.

[0009] The application determines the mine earthquake main control layer under the scenario that the mine earthquake has a disaster-causing impact on the ground and the scenario that the mine earthquake has a disaster-causing impact on the underground mining working face respectively by analyzing different mine earthquake disaster-causing situations, and dynamically corrects the target rock layer position of the mine earthquake main control layer through well logging analysis, so that the target layer position of the mine earthquake main control layer can be constructed and managed.

[0010] Optionally, in the S1, when the mine earthquake control layer is determined, the following steps are included:

[0011] S11, determining a key layer having a mine earthquake disaster-causing impact on the underground mining working face;

[0012] calculating the energy released by the breaking of the current key layer, determining the vibration energy attenuation law by analyzing the measured energy of each sensor in the historical mine earthquake cases of the current area, calculating the residual energy transferred to the working face according to the energy released by the breaking of the key layer and the vibration energy attenuation law, and evaluating the mine earthquake induced shock risk according to the residual energy, and determining the rock layer satisfying the induced shock condition as the mine earthquake control layer;

[0013] S12, determining a mine earthquake main control layer having a disaster-causing impact on the ground;

[0014] monitoring the velocity transferred to the ground when the mine earthquake occurs by arranging vibration velocity monitoring sensors on the ground, obtaining the maximum vibration velocity of different rock layers by theoretical analysis, determining the mine earthquake intensity level corresponding to different layer positions according to the maximum vibration velocity of different layer positions, and determining the mine earthquake control layer with a mine earthquake intensity level of II and above as the mine earthquake control layer.

[0015] Further, in the process of determining the multiple mine earthquake main control layers, the mine earthquake control layers causing the mine earthquake disaster of the underground mining working face are preferentially determined and analyzed level by level, and N-level mine earthquake main control layers are determined;

[0016] The mine earthquake control layers causing the mine earthquake disaster of the ground are further added to the level-by-level determination and analysis, and n-level mine earthquake main control layers are determined;

[0017] N+n-level mine earthquake main control layers are formed in a unified sequence, and the N+n-level mine earthquake main control layers are arranged level by level according to the determination order.

[0018] Further, in the S2, when the multiple mine earthquake control layers causing the mine earthquake of the underground mining working face are determined level by level, the level-by-level determination and management analysis principles include:

[0019] The residual energy of the multiple mine earthquake control layers is compared, the mine earthquake control layer with the maximum residual energy transferred to the working face is taken as the first-level mine earthquake main control layer, and the first-level mine earthquake control layer is managed and analyzed;

[0020] Repeat S1 step, other control layer of the adjacent main control layer of the governance to mine earthquake risk judgment, eliminate mine earthquake risk mine earthquake control layer, the mine earthquake control layer adjacent to the main control layer of the governance and still exist mine earthquake risk as the next level of disaster-causing main control layer;

[0021] Circulate this step until no other key layer exists mine earthquake disaster risk, form a multi-level mine earthquake main control layer.

[0022] Further, in the determination analysis of the mine earthquake control layer causing ground mine earthquake disaster, including:

[0023] The last level of mine earthquake main control layer in the adjacent N level mine earthquake main control layer, and the mine earthquake control layer which can cause surface mine earthquake disaster is analyzed. The mine earthquake main control layer analysis is carried out according to S12 step, and this step is repeated until the n level mine earthquake main control layer which can cause ground mine earthquake is determined.

[0024] Further, in the S2, the calculation of the energy released by the key layer breaking includes the calculation of the elastic energy released by the initial breaking of the key layer and the elastic energy released by the periodic breaking, and the energy released by the key layer breaking at different heights from the working face is obtained.

[0025] Further, the elastic energy released by the initial breaking of the key layer is calculated, and the elastic energy released by the initial breaking of the key layer is calculated as:

[0026]

[0027] In the formula, U dc The elastic energy released by the initial breaking of the key layer is calculated as: dc The unit length conversion load of the self weight of the basic roof and the additional load on the overlying strata; E d The elastic modulus of the basic roof beam; R t The tensile strength of the basic roof rock; h d The thickness of the basic roof; b dc The suspended span of the basic roof at the initial breaking.

[0028] Further, the elastic energy released by the periodic breaking of the key layer is calculated, and the elastic energy released by the periodic breaking of the thick basic roof is calculated as:

[0029]

[0030] In the formula, U dz The elastic energy released by the initial breaking of the key layer is calculated as: dz The unit length conversion load of the self weight of the basic roof and the additional load on the overlying strata; b dz The suspended span of the basic roof at the initial breaking; E d The elastic modulus of the basic roof beam; Rt is the tensile strength of the basic roof rock; h d is the thickness of the basic roof.

[0031] Further, in S12, when analyzing the measured energy of each sensor in the current area historical mine earthquake case, the following is included:

[0032] At multiple measuring points, the vibration wave of the mine earthquake occurring in the working face is measured to form a velocity time history curve, which reflects the maximum vibration velocity of each sensor;

[0033] The velocity time history curve is integrated to obtain a displacement time history curve, which reflects the position change of the medium particle in the vibration wave propagation process;

[0034] When the energy transmitted to the working face is calculated, the energy of the vibration wave at a measuring point position is calculated according to the fact that the energy of the vibration wave at a certain place is proportional to the square of the maximum particle motion velocity at that place, and an energy and maximum particle motion velocity relationship equation at that place is established:

[0035]

[0036] Where A0 is the maximum particle motion velocity at a certain place; C is the fitting coefficient; E is the radiation energy of the vibration wave;

[0037] The energy calculation results of the vibration wave measuring points are sorted to obtain a maximum vibration velocity-distance curve, where the distance refers to the distance between each sensor and the seismic source, and the vibration energy attenuation law is obtained from the maximum vibration velocity-distance curve.

[0038] Further, in S11, when evaluating the mine earthquake induced shock risk based on the residual energy, a residual energy induced shock risk determination system is established based on the energy accumulation, transfer, release and dissipation after mining in the mine, the induced shock conditions are determined, and the rock stratum with residual energy satisfying the induced shock conditions is determined as the main control layer of the mine earthquake that can cause rock burst;

[0039] When establishing the residual energy induced shock risk determination system, the system equation is established as follows:

[0040] E e + ΔE = E P + E r ;

[0041] Where E e is the elastic deformation energy of the coal and rock mass before excavation;

[0042] ΔE is the deformation energy increment caused by mining disturbance;

[0043] E P is the energy dissipated by plastic deformation and destruction of the coal and rock mass after excavation;

[0044] Er is the residual energy.

[0045] Further, the determination of the bumping condition comprises: when the energy E P <E e + ΔE, then E r > 0, at this time E r is the energy source of the bumping ground pressure, when the residual energy E r is completely converted into the kinetic energy of the coal and rock thrown, that is, the bumping ground pressure occurs.

[0046] Further, when the well logging analysis in S3 is performed, the well logging analysis method comprises resistivity logging, acoustic logging and natural gamma logging, wherein:

[0047] The resistivity logging judges the lithology and fluid properties by measuring the resistivity of the formation;

[0048] The acoustic logging measures the propagation velocity and amplitude attenuation of the acoustic wave in the formation, and calculates the porosity and lithology;

[0049] The natural gamma logging detects the natural radioactivity of the formation, and is used for analyzing the lithology composition, formation boundary and mineral content.

[0050] Further, when the target rock layer corresponding to the main control layer appears a mudstone layer or the argillaceous content of the rock layer is greater than the argillaceous content requirement of the fine sandstone, the fracturing test is further performed to judge the fracturing property of the rock layer and whether the fracturing treatment is performed;

[0051] According to the rock layer condition obtained by the well logging, in the order from the current target layer to the farthest, the target rock layer of the main control layer is modified and selected to the nearest rock layer with the sandstone content reaching the fine sandstone content standard and above.

[0052] The additional aspects and advantages of the present application will be partially given in the following description, partially will become obvious from the following description, or will be understood by the practice of the present application. BRIEF DESCRIPTION OF DRAWINGS

[0053] The above and / or additional aspects and advantages of the present application will become apparent and be readily understood from the following description, taken in conjunction with the accompanying drawings, in which:

[0054] Figure 1 is a schematic diagram of the overall method steps of the mine shock main control layer determination method according to the present application;

[0055] Figure 2 is a schematic diagram of the detailed steps of S1 of the mine shock main control layer determination method according to the present application

[0056] Figure 3 is a schematic diagram of the low-position thick and hard roof initial breaking bumping of the mine shock main control layer determination method according to the present application;

[0057] Figure 4 It is a low position thick hard roof breaking period and inducing shock schematic view according to the mine earthquake main control layer position determination method of the application;

[0058] Figure 5 It is a management analysis schematic view of a certain level mine earthquake main control layer according to the mine earthquake main control layer position determination method of the application;

[0059] Figure 6 It is a determination and management analysis schematic view of a lower level mine earthquake main control layer after a certain level mine earthquake main control layer is managed according to the mine earthquake main control layer position determination method of the application;

[0060] Figure 7 It is a mine earthquake corresponding maximum vibration velocity-distance curve schematic view according to the mine earthquake main control layer position determination method of the application. DETAILED DESCRIPTION

[0061] Embodiments of the application are described below in detail, examples of which are shown in the accompanying drawings, wherein the same or similar reference signs represent the same or similar elements or elements having the same or similar functions throughout. The embodiments described below by referring to the accompanying drawings are exemplary and are intended to explain the application, and cannot be understood as a limitation of the application.

[0062] The application proposes a mine earthquake main control layer position determination method, which is described in detail below with reference to Figures 1 to 7 .

[0063] A mine earthquake main control layer position determination method comprises the following steps:

[0064] S1, according to three-dimensional seismic exploration, drilling core, and microseismic positioning analysis, a mine area stratum model is established, a working face is selected, and a key layer distributed above a fracture zone is determined through a key layer theory, a plurality of key layers located above the fracture zone in the mine area stratum model and having mine earthquake risks are determined as mine earthquake control layers;

[0065] S2, a main control layer analysis is performed on all mine earthquake control layers above the fracture zone, a step-by-step determination and management analysis principle is formulated for the mine earthquake control layers, mine earthquake control layers that cannot produce disaster-causing effects after the current mine earthquake control layers are managed are eliminated, and a plurality of mine earthquake main control layers are determined from all mine earthquake control layers;

[0066] S3, a well logging correction analysis is performed, and a layer position of a target rock layer of the plurality of mine earthquake main control layers is dynamically corrected.

[0067] The application determines the mine earthquake main control layer under the scenario that the mine earthquake has a disaster-causing impact on the ground and the scenario that the mine earthquake has a disaster-causing impact on the underground mining working face by analyzing different mine earthquake disaster-causing situations, and dynamically corrects the target rock layer position of the mine earthquake main control layer through well logging analysis, so that the target layer position of the mine earthquake main control layer can be constructed and managed.

[0068] In some embodiments, in the S1, when the mine earthquake control layer is determined, the following steps are included:

[0069] S11, determining a mine earthquake control layer having a mine earthquake disaster-causing impact on the underground mining working face;

[0070] The energy released by the breaking of the current key layer is calculated, the vibration energy attenuation law is determined by analyzing the measured energy of each sensor in the historical mine earthquake cases of the current area, the residual energy transmitted to the working face is calculated according to the breaking release energy of the key layer and the vibration energy attenuation law, the mine earthquake triggering risk is evaluated by the residual energy, and the rock layer satisfying the triggering condition is determined as the mine earthquake control layer;

[0071] S12, determining a mine earthquake main control layer having a disaster-causing impact on the ground;

[0072] The vibration velocity monitoring sensor is arranged on the ground to monitor the velocity transmitted to the ground when the mine earthquake occurs, the maximum vibration velocity of the breaking of different rock layers is obtained through theoretical analysis, the mine earthquake intensity level corresponding to different layer positions is determined according to the maximum vibration velocity of different layer positions, and the mine earthquake control layer with a mine earthquake intensity level of II and above is determined as the mine earthquake control layer.

[0073] When the mine earthquake control layer having a mine earthquake disaster-causing impact on the underground mining working face and the mine earthquake main control layer having a disaster-causing impact on the ground appear at the same time, the mine earthquake control layer having a mine earthquake disaster-causing impact on the underground mining working face is determined first.

[0074] In some embodiments, in the S11, when the energy released by the breaking of the key layer is calculated, the elastic energy released by the initial breaking of the key layer and the elastic energy released by the periodic breaking of the key layer are calculated to obtain the energy of the key layer at different heights from the working face. Refer to Figures 2 to 3 Taking the low-position thick hard roof as an example, the impact load analysis of the initial breaking triggering of the key layer and the periodic breaking triggering is shown.

[0075] In some embodiments, the elastic energy released by the initial breaking of the key layer is calculated as follows according to the fixed-end beam:

[0076]

[0077] In the formula, U dc is the elastic energy released by the initial breaking of the basic roof;

[0078] q dc unit length conversion load of self-weight of basic roof and additional load of overlying strata;

[0079] E d beam elastic modulus of basic roof;

[0080] R t tensile strength of basic roof strata;

[0081] h d thickness of basic roof;

[0082] b dc unsupported span of basic roof at initial fracture.

[0083] In some embodiments, the elastic energy released by periodic fracture of key strata is calculated, and the elastic energy released by periodic fracture of thick basic roof is calculated as a cantilever beam:

[0084]

[0085] In the formula, U dz elastic energy released by periodic fracture of basic roof;

[0086] q dz unit length conversion load of self-weight of basic roof cantilever beam model and additional load of overlying strata;

[0087] b dz unsupported span of basic roof at periodic fracture;

[0088] E d beam elastic modulus of basic roof;

[0089] R t tensile strength of basic roof strata;

[0090] h d thickness of basic roof.

[0091] In some embodiments, when analyzing the measured energy of each sensor in the current area historical mine earthquake case, it includes:

[0092] The vibration wave of the mine earthquake occurring in the working face is measured at multiple measuring points to obtain the measured waveform. The velocity time history curve is formed according to the measured waveform, and the velocity time history curve is used to reflect the maximum vibration velocity of each sensor;

[0093] The displacement time history curve is obtained by integrating the velocity time history curve, and the displacement time history curve is used to reflect the position change of the medium particle in the vibration wave propagation process;

[0094] To calculate the energy of a vibration wave at a specific measuring point, based on the principle that the energy of a vibration wave at a given point is proportional to the square of the maximum velocity of the particle at that point, an equation relating the energy to the maximum velocity of the particle at that point is established:

[0095]

[0096] In the formula, A0 is the maximum velocity of a particle at a certain point;

[0097] C is the fitting coefficient;

[0098] E represents the energy radiated by the vibration wave;

[0099] By analyzing the energy calculation results from the seismic wave measurement points, a maximum vibration velocity-distance curve was obtained. Here, distance refers to the distance between each sensor and the seismic source. Figure 7 As shown, the vibration energy attenuation law is obtained from the maximum vibration velocity-distance curve.

[0100] In some embodiments, in S11, when assessing the risk of mine-induced shock using residual energy, a residual energy-induced shock risk determination system is established based on the analysis of energy accumulation, transfer, release and dissipation after mine mining, to determine the induced shock conditions, and to identify the rock strata whose residual energy meets the induced shock conditions as the main control layer of mine-induced shock that can cause rockburst.

[0101] When establishing the residual energy-induced impulse risk assessment system, the system equations are as follows:

[0102] E e +ΔE=E P +E r ;

[0103] In the formula, E e This refers to the elastic deformation energy of the coal and rock mass before excavation.

[0104] ΔE represents the increase in deformation energy caused by mining disturbance;

[0105] E P The energy dissipated during the plastic deformation and destruction of the excavated coal and rock mass;

[0106] E r This represents the remaining energy.

[0107] Induced impulse conditions include: when E in the system P <E e When +ΔE, then E r >0, at this time E r This is the energy source of the rockburst, when the remaining energy E r All of this energy is converted into the kinetic energy of the ejected coal and rock, resulting in rockburst.

[0108] The vibration energy attenuation law of mine earthquake is determined by analyzing the measured energy of each sensor in the current regional historical mine earthquake cases, and the main control rock stratum of mine earthquake is determined combined with the energy release of key stratum breakage. The possibility of inconsistency between pure theoretical calculation and actual situation is reduced by combining theoretical calculation with field monitoring.

[0109] In some embodiments, with reference to Figures 5 to 6 In the process of determining the multi-level mine earthquake main control stratum in S2 step, the mine earthquake control stratum causing the mine earthquake disaster of underground mining working face is preferentially determined and analyzed level by level, and N-level mine earthquake main control stratum is determined.

[0110] Then, the mine earthquake control stratum causing the mine earthquake disaster of ground surface is added to the level-by-level determination and analysis, and n-level mine earthquake main control stratum is determined.

[0111] The N+n-level mine earthquake main control stratum is formed in a unified order, and the N+n-level mine earthquake main control stratum is arranged level by level according to the determination order.

[0112] In some embodiments, when the multiple mine earthquake control strata causing the mine earthquake of underground mining working face are determined level by level, the level-by-level determination and management analysis principles include:

[0113] The residual energy of the multiple mine earthquake control strata is compared, the mine earthquake control stratum with the maximum residual energy transmitted to the working face is taken as the first-level mine earthquake main control stratum, and the first-level mine earthquake control stratum is managed and analyzed;

[0114] S1 step is repeated to judge the mine earthquake risk of other control strata near the managed main control stratum, the mine earthquake control stratum with no mine earthquake risk is removed, and the mine earthquake control stratum near the managed main control stratum and still with mine earthquake risk is taken as the disaster-causing main control stratum to be managed in the next level;

[0115] This step is cycled until there is no other key stratum with mine earthquake disaster risk, and the multi-level mine earthquake main control stratum is formed.

[0116] In some embodiments, when the mine earthquake control stratum causing the mine earthquake disaster of ground surface is added to the level-by-level determination and analysis, it includes:

[0117] The mine earthquake control stratum near the last level of N-level mine earthquake main control stratum and capable of causing ground surface mine earthquake disaster is analyzed as a mine earthquake main control stratum, and this mine earthquake main control stratum analysis is performed according to S12 step. This step is cycled until n-level mine earthquake main control stratum capable of causing ground surface mine earthquake is determined.

[0118] More specifically, when the multi-level mine earthquake main control stratum is determined level by level:

[0119] The control layer causing the mining shock risk of the underground mining working face is analyzed and determined in stages, the control layer of the mining shock causing the maximum residual energy transferred to the working face is selected as the primary mining shock control layer in the stratum model of the mining area, and the primary mining shock control layer is simulated and analyzed for treatment;

[0120] After the primary mining shock control layer is treated, the residual mining shock control layers on the upper and lower sides of the primary mining shock control layer are determined for mining shock risk, the mining shock control layer still causing the mining shock risk is determined as the secondary mining shock control layer, and the step is repeated until the N-level mining shock control layer causing the underground mining shock is determined;

[0121] The mining shock control layer causing the ground mining shock risk is considered, the N-level mining shock control layer causing the underground mining shock is analyzed and treated in stages, the last level mining shock control layer in the N-level mining shock control layer and the mining shock control layer causing the ground mining shock disaster are analyzed, the mining shock control layer analysis is performed according to the step S122, and the step is repeated until the n-level mining shock control layer causing the ground mining shock is determined.

[0122] The N+n-level mining shock control layer is determined in stages, and the N+n-level mining shock control layer is arranged in the determined order in stages, and after the N+n-level mining shock control layer is treated and analyzed, the mining shock control layer outside the two zones no longer has the mining shock risk.

[0123] In some embodiments, when the well logging analysis in S3 is performed, the target rock layer corresponding to the determined mining shock control rock layer is observed and measured by well logging, the lithology, rock strength, sand content, and mud content factors are obtained, the fracturing property of the target rock layer is judged, and the target rock layer position of the mining shock control rock layer in the treatment analysis stage is dynamically corrected.

[0124] In some embodiments, the well logging analysis method includes resistivity logging, acoustic logging, and natural gamma logging, wherein,

[0125] The resistivity logging principle is to judge the lithology and fluid properties by measuring the resistivity of the stratum;

[0126] The specific techniques include:

[0127] The lateral logging technique transmits focused current to the stratum, reduces the influence of the wellbore and surrounding rock, and is suitable for high-resistivity strata;

[0128] The induction logging technique measures the stratum conductivity using the electromagnetic induction principle and is suitable for low-resistivity or wellbores filled with conductive mud;

[0129] Acoustic logging, principle: measure the propagation velocity (time difference) and amplitude attenuation of acoustic waves in the formation, calculate the porosity and lithology; acoustic waves propagate fast in dense rock (such as limestone, hard sandstone), and slow in loose or porous rock; the ratio of P-wave and S-wave velocity can also be used to evaluate the mechanical properties of the formation (such as fracture pressure);

[0130] Natural gamma logging, principle: detect the natural radioactivity of the formation (mainly from potassium, uranium, thorium isotopes), among which the radioactivity of mudstone is high, and that of sandstone, limestone, etc. is low, which is used to analyze the lithology composition, formation boundary and mineral content;

[0131] In some embodiments, when the target rock layer corresponding to the main control layer position appears a mudstone layer or the rock layer has a shale content greater than the shale content requirement of fine sandstone, then the fracturing test is performed to determine the fracturing property of the rock layer and whether to perform fracturing treatment.

[0132] According to the rock layer condition obtained by logging, the target rock layer of the main control layer position is selected and modified to the nearest rock layer with sandstone content reaching or above the fine sandstone content standard in order from the current target layer position.

[0133] In some embodiments, in order to further accurately determine the disaster-causing layer position of mine shock, box plot + event scatter normal distribution diagram are used in stages to statistically analyze the distribution law of microseismic events of different energy levels in the vertical direction, and to further verify whether the determined layer position of the mine shock main control layer is correct.

[0134] In the description of the present application, it should be understood that the orientations or positional relationships indicated by the terms "center", "longitudinal", "transverse", "length", "width", "thickness", "upper", "lower", "front", "rear", "left", "right", "vertical", "horizontal", "top", "bottom", "inner", "outer", "clockwise", "counterclockwise", "axial", "radial", "circumferential" are based on the orientations or positional relationships shown in the drawings, and are only for the convenience of describing the present application and simplifying the description, and therefore cannot be understood as indicating or implying that the devices or elements indicated must have a particular orientation, be constructed and operated in a particular orientation, and therefore cannot be understood as limiting the present application.

[0135] In addition, the terms "first" and "second" are only for descriptive purposes, and cannot be understood as indicating or implying relative importance or implicitly indicating the number of the indicated technical features. Therefore, the features defined with "first" and "second" can explicitly or implicitly include at least one of the features. In the description of the present application, the meaning of "a plurality of" is at least two, such as two, three, etc., unless otherwise specifically limited.

[0136] In the present application, unless specifically defined otherwise, the terms "mounting", "connected", "connecting", "fixed", "fixedly connected", "connection", "fixedly connected", "connected", "fixed", and the like should be construed broadly, for example, can be fixed connection, can be detachable connection, or integral; can be mechanical connection, can be electrical connection or communication with each other; can be directly connected, or indirectly connected through an intermediate medium, can be the internal communication of two elements or the interaction relationship between two elements, unless otherwise specifically defined. For those skilled in the art, the specific meaning of the above terms in the present application can be understood according to the specific circumstances.

[0137] In the present application, unless specifically defined otherwise, the first feature is "on" or "under" the second feature can be that the first and second features are in direct contact, or the first and second features are in indirect contact through an intermediate medium. Moreover, the first feature "on", "above" and "on" the second feature can be that the first feature is directly above or obliquely above the second feature, or only indicates that the first feature is higher than the second feature in horizontal height. The first feature "under", "below" and "below" the second feature can be that the first feature is directly below or obliquely below the second feature, or only indicates that the first feature is lower than the second feature in horizontal height.

[0138] In the present application, the terms "one embodiment", "some embodiments", "example", "specific example", or "some examples" mean that the specific features, structures, materials or characteristics described in connection with the embodiment or example are included in at least one embodiment or example of the present application. In the present specification, the illustrative expressions of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described can be combined in any appropriate manner in any one or more embodiments or examples. In addition, those skilled in the art can combine and combine different embodiments or examples described in the specification and the features of different embodiments or examples without contradiction.

[0139] Although the embodiments of the present application have been shown and described above, it is understood that the above embodiments are exemplary and should not be construed as limiting the present application, and those skilled in the art can make changes, modifications, replacements and variations to the above embodiments within the scope of the present application.

Claims

1. A method for determining a main control horizon of a mine shock, characterized in that, The method comprises the following steps: S1, according to three-dimensional seismic exploration, drilling core, microseismic positioning analysis, a stratum model of the mining area is established, a working face is selected, and a key layer distributed above a fissure zone is determined by a key layer theory, a plurality of key layers above the fissure zone in the stratum model of the mining area and having a mine earthquake risk are determined as mine earthquake control layers; S2, main control layer analysis is performed on all mine earthquake control layers above the fissure zone, a step-by-step determination and management analysis principle is formulated for the mine earthquake control layers, mine earthquake control layers that cannot cause disaster after being managed are removed, and a plurality of levels of mine earthquake main control layers are determined from all mine earthquake control layers; S3, logging correction analysis is performed on the layer position of the target rock layer of the plurality of levels of mine earthquake main control layers.

2. The method of claim 1, wherein the master horizon is determined by the steps of: determining a first master horizon; determining a second master horizon; and determining a third master horizon. In S1, when the mine earthquake control layers are determined, the following steps are included: S11, mine earthquake control layers having a mine earthquake disaster-causing effect on the underground mining working face are determined; The energy released by the breaking of the current key layer is calculated, the vibration energy attenuation law is determined by analyzing the measured energy of each sensor in the historical mine earthquake cases of the current area, the residual energy transmitted to the working face is calculated according to the breaking energy of the key layer and the vibration energy attenuation law, the mine earthquake triggering risk is evaluated by the residual energy, and the rock layer satisfying the triggering condition is determined as the mine earthquake control layer; S12, mine earthquake control layers having a disaster-causing effect on the ground are determined; The vibration velocity monitoring sensor is arranged on the ground to monitor the velocity transmitted to the ground when the mine earthquake occurs, the maximum vibration velocity of the breaking of different rock layers is obtained through theoretical analysis, the mine earthquake intensity grade corresponding to different layer positions is determined according to the maximum vibration velocity of different layer positions, and the mine earthquake control layer with a mine earthquake intensity grade of II and above is determined as the mine earthquake control layer.

3. The method of claim 2, wherein the master horizon is determined by: In the process of determining the plurality of levels of mine earthquake main control layers, the mine earthquake control layers causing the mine earthquake disaster of the underground mining working face are preferentially determined and analyzed step by step, and N levels of mine earthquake main control layers are determined; The step-by-step determination and analysis of the mine earthquake control layers causing the mine earthquake disaster of the ground are added, and n levels of mine earthquake main control layers are determined; N+n levels of mine earthquake main control layers are formed in a unified sequence, and the N+n levels of mine earthquake main control layers are arranged step by step according to the determination sequence.

4. The method of claim 3, wherein the master horizon is determined by the steps of: determining a first master horizon; determining a second master horizon; and determining a third master horizon. In S2, when the plurality of mine earthquake control layers causing the mine earthquake of the underground mining working face are determined step by step, the step-by-step determination and management analysis principle includes: The residual energy of the plurality of mine earthquake control layers is compared, the mine earthquake control layer with the maximum residual energy transmitted to the working face is taken as the first level of mine earthquake main control layer, and the first level of mine earthquake control layer is managed and analyzed; S1 is repeated, the other control layers near the managed main control layer are judged for the mine earthquake risk, the mine earthquake control layer without the mine earthquake risk is removed, the mine earthquake control layer near the managed main control layer and still having the mine earthquake risk is taken as the disaster-causing main control layer to be managed in the next level; The step is repeated until there is no other key layer having the mine earthquake disaster risk, and the plurality of levels of mine earthquake main control layers are formed.

5. The method for determining the main seismic stratigraphic position as described in claim 4, characterized in that, When the step-by-step determination and analysis of the mine earthquake control layers causing the mine earthquake disaster of the ground are added, the following steps are included: The mine earthquake main control layer analysis is performed on the mine earthquake control layer which is the last mine earthquake main control layer of the mine earthquake main control layer close to the Nth mine earthquake main control layer and can cause surface mine earthquake disaster. The mine earthquake main control layer analysis is performed according to the S12 step. The step is repeated until the nth mine earthquake main control layer which can cause surface mine earthquake is determined.

6. The method of claim 2, wherein, In the S11, the energy released by the key layer breaking is calculated, including the elastic energy released by the key layer initial breaking and the elastic energy released by the key layer periodic breaking, to obtain the energy released by the key layer breaking at different heights from the working face.

7. The method for determining the main seismic stratigraphic position as described in claim 6, characterized in that, The elastic energy released by the key layer initial breaking is calculated. The elastic energy released by the key layer initial breaking is calculated according to the fixed support beam as follows: In the formula, U dc is the initial release energy of the basic roof; q dc is the unit length load of the dead weight of the basic roof and the additional load of the overlying rock; E d is the beam elastic modulus of the basic roof; R t is the tensile strength of the basic roof; h d is the thickness of the basic roof; b dc is the suspended span of the basic roof at the initial break.

8. The method for determining the main seismic stratigraphic position as described in claim 7, characterized in that, The elastic energy released by the key layer periodic breaking is calculated. The elastic energy released by the key layer periodic breaking is calculated according to the cantilever beam as follows: In the formula, U dz is the basic roof periodic fracture release elastic energy; q dz is the unit length conversion load of the dead weight of the basic roof cantilever beam model and the overburden load thereon; b dz is the suspended span when the basic roof periodically breaks; E d is the beam elastic modulus of the basic roof; R t is the tensile strength of the basic roof rock; h d is the thickness of the basic roof.

9. The method for determining the main seismic stratigraphic position as described in claim 2, characterized in that, In the S11, the measured energy of each sensor in the current area historical mine earthquake case is analyzed, including: The vibration wave of the working face mine earthquake is measured at multiple measuring points to form a velocity time history curve to reflect the maximum vibration velocity of each sensor; The velocity time history curve is integrated to obtain a displacement time history curve to reflect the position change of the medium particle in the vibration wave propagation process; The energy transmitted to the working face is calculated. The energy of the vibration wave at a measuring point is calculated according to the fact that the energy of the vibration wave at a certain place is proportional to the square of the maximum particle motion velocity at the place. The energy and the maximum particle motion velocity at a certain place are related as follows: Wherein, A0 is the maximum particle motion velocity at a certain place; C is the fitting coefficient; E is the radiation energy of the vibration wave; The energy calculation results of the vibration wave measuring points are sorted to obtain a maximum vibration velocity-distance curve, wherein the distance refers to the distance between each sensor and the seismic source. The vibration energy attenuation law is obtained from the maximum vibration velocity-distance curve.

10. The method for determining the main seismic stratigraphic position as described in claim 2, characterized in that, In the S11, the mine earthquake inducing shock risk is evaluated according to the residual energy. A residual energy inducing shock risk judgment system is established according to the energy accumulation, transfer, release and dissipation after mining to determine the inducing shock conditions. The rock stratum with residual energy satisfying the inducing shock conditions is determined as the mine earthquake main control layer which can cause rock burst. When the residual energy inducing shock risk judgment system is established, the system equation is established as follows: E e + ΔE = E P + E r ; In the formula, E e is the elastic deformation energy of the coal rock mass before excavation; ΔE is the deformation energy increment caused by mining disturbance. E P To produce plastic deformation and destroy the energy dissipation of the coal and rock mass after excavation; E r For the remaining energy.

11. The method for determining the main seismic stratigraphic position as described in claim 10, characterized in that, The determination of the inducing condition includes: when E P <E e +ΔE, E r >0, E r is the energy source of the rock burst, when the residual energy E r is all converted into the kinetic energy of the coal and rock, the rock burst occurs.

12. The method of claim 1, wherein, In the well logging analysis in S3, the target rock stratum corresponding to the determined mine earthquake main control rock stratum is observed and measured by well logging to obtain the lithology, rock stratum strength, sandstone content, and clay content factors. The fracturing property of the target rock stratum is judged. The target rock stratum horizon of the mine earthquake main control rock stratum in the treatment analysis stage is dynamically corrected.

13. The method for determining the main seismic stratigraphic position as described in claim 12, characterized in that, In the well logging analysis in S3, the well logging analysis method includes resistivity logging, acoustic logging, and natural gamma logging, wherein: The resistivity logging is used to judge the lithology and fluid properties by measuring the resistivity of the stratum; The acoustic logging is used to measure the propagation velocity and amplitude attenuation of the acoustic wave in the stratum to calculate the porosity and lithology; The natural gamma logging is used to detect the natural radioactivity of the stratum to analyze the lithology composition, stratum boundary, and mineral content.

14. The method of claim 13, wherein, When the target rock stratum corresponding to the main control horizon has a mudstone layer or the mud content of the rock stratum is greater than the fine sandstone mud content requirement, the fracturing property test is performed again to judge the fracturing property of the rock stratum and whether the fracturing treatment is performed. According to the stratum condition obtained by logging, the target stratum of the main control horizon is selected and corrected to the nearest stratum with sand content reaching the fine sand content standard or above in the order from the current target horizon to the farthest.

Images