CN114856706A - Method for comprehensively evaluating disturbance influence of main key layer based on subsidence - Google Patents
Method for comprehensively evaluating disturbance influence of main key layer based on subsidence Download PDFInfo
- Publication number
- CN114856706A CN114856706A CN202210530296.2A CN202210530296A CN114856706A CN 114856706 A CN114856706 A CN 114856706A CN 202210530296 A CN202210530296 A CN 202210530296A CN 114856706 A CN114856706 A CN 114856706A
- Authority
- CN
- China
- Prior art keywords
- main key
- key layer
- disturbance
- layer
- subsidence
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 21
- 239000003245 coal Substances 0.000 claims abstract description 100
- 238000005065 mining Methods 0.000 claims abstract description 27
- 238000010606 normalization Methods 0.000 claims abstract description 18
- 230000008878 coupling Effects 0.000 claims abstract description 12
- 238000010168 coupling process Methods 0.000 claims abstract description 12
- 238000005859 coupling reaction Methods 0.000 claims abstract description 12
- 238000005553 drilling Methods 0.000 claims abstract description 11
- 238000011156 evaluation Methods 0.000 claims abstract description 11
- 238000009826 distribution Methods 0.000 claims abstract description 8
- 238000002474 experimental method Methods 0.000 claims abstract description 8
- 230000002452 interceptive effect Effects 0.000 claims abstract description 8
- 239000010410 layer Substances 0.000 claims description 222
- 239000011435 rock Substances 0.000 claims description 38
- 238000004364 calculation method Methods 0.000 claims description 26
- 230000015572 biosynthetic process Effects 0.000 claims description 22
- 238000009933 burial Methods 0.000 claims description 21
- 230000004044 response Effects 0.000 claims description 13
- 239000011229 interlayer Substances 0.000 claims description 10
- 238000004088 simulation Methods 0.000 claims description 5
- 238000012216 screening Methods 0.000 claims description 3
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims 4
- 238000012545 processing Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000008961 swelling Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000003449 preventive effect Effects 0.000 description 2
- 238000004062 sedimentation Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21F—SAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
- E21F17/00—Methods or devices for use in mines or tunnels, not covered elsewhere
- E21F17/18—Special adaptations of signalling or alarm devices
- E21F17/185—Rock-pressure control devices with or without alarm devices; Alarm devices in case of roof subsidence
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21F—SAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
- E21F17/00—Methods or devices for use in mines or tunnels, not covered elsewhere
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21F—SAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
- E21F17/00—Methods or devices for use in mines or tunnels, not covered elsewhere
- E21F17/18—Special adaptations of signalling or alarm devices
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Engineering & Computer Science (AREA)
- Mining & Mineral Resources (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geology (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
The invention discloses a subsidence-based main key layer disturbance influence comprehensive evaluation method, and belongs to the field of coal mining disturbance influence evaluation. The method mainly comprises the following steps: (1) collecting, sorting and analyzing geological drilling data; (2) designing an interactive experiment scheme; (3) calculating the disturbance degree of a main key layer of each experimental scheme; (4) calculating the subsidence of the main key layer of each experimental scheme; (5) carrying out normalization processing on the data, and determining the coupling relation between the disturbance degree and the subsidence of the main key layer; (6) based on the coupling relation, classifying the disturbance influence level of the main key layer; (7) calculating the disturbance degree normalization values of all drilling main key layers in the region; (8) and matching the coordinates of the drilled holes with the disturbance degree normalized value to finally obtain the grade distribution of the disturbance degree normalized value of the main key layer in the region. The method provides a theoretical basis for qualitatively describing the relative size of the main key layer in the area influenced by the disturbance of coal seam mining and the variation trend of the subsidence of the main key layer.
Description
Technical Field
The invention relates to a comprehensive evaluation method for the disturbance influence degree of coal seam mining on an overlying main key layer, in particular to a comprehensive evaluation method for the disturbance influence degree of main key layer subsidence on the coal seam mining.
Background
Since the beginning of coal mining activities, the law of overburden movement has been one of the core problems in mining science research. With the extraction of coal resources, different degrees of influence on overburden rock can be inevitably caused, expert and scholars analyze the influence law of coal seam mining on overburden rock movement under different conditions from different angles or by different methods respectively, and a series of relatively complete theoretical technical systems and surrounding rock control methods are established.
Because of different diagenesis time and composition, multiple layers of rock strata with different thicknesses and strengths are formed above the coal seam, and the rock strata generally control the activities of partial or all rock strata of overlying rock strata of a stope, and are one of important research objects in coal mining pressure and rock stratum control.
Expert scholars put forward a key layer theory of rock stratum control according to years of research and practice on roof rock stratum control, the rock stratum which plays a role in determining all or part of rock mass activity in the stratum is called a key layer, the former is called a main key layer, the latter is called a sub-key layer, and a method for judging the position of the key layer in overlying strata is established.
The existing research mainly focuses on the application of the movement of a key layer to surrounding rock mining fractures, surface movement, stope mine pressure and the rock stratum control theory, the disturbance influence degree of the key layer on the main key layer when factors such as coal thickness and buried depth change are seldom concerned, the coal mining geological conditions are complex and changeable, the mining of coal under different geological conditions can cause the disturbance influence of different degrees on an overlying key layer, and the integrity degree, the breaking position, the stability and the like of the key layer can affect the stress distribution of the overlying rock, the mining difficulty of a coal bed and the safety of stope operation.
Disclosure of Invention
The invention aims to provide a comprehensive evaluation method for disturbance influence of a main key layer, which has strong universality and high reliability of evaluation results, predicts the relative size of the disturbance influence degree of coal seam mining in an area on overlying strata in advance, reflects the relative size of subsidence of the main key layer in the area, takes precautionary measures in advance for subsequent coal seam mining, and provides reference for ensuring safe production of a working face.
In order to solve the technical problems, the invention adopts the technical scheme that: a method for comprehensively evaluating disturbance influence of a main key layer based on subsidence comprises the following steps:
s1, collecting, sorting and analyzing geological drilling data: the method comprises the steps of numbering drill holes by collecting geological drill hole data of the whole mining area or a local mining area, calculating a main key layer position of each drill hole, counting information such as coordinates of each drill hole, the thickness of a first layer of coalable layer, the burial depth of a coal layer, the distance between the coal layer and the main key layer, physical and mechanical parameters (including tensile strength, crushing and swelling coefficients and the like) of rock layers between the coal layer and the main key layer, the thickness of each rock layer between the coal layer and the main key layer and the like, finally counting and analyzing the maximum value and the minimum value of relevant factors such as the thickness of the coalable layer in the area, the burial depth of the coal layer, the distance between the coal layer and the main key layer and the like, and continuing to perform the step S2;
s2, designing an interactive experiment scheme: according to the related drilling data arranged in the step S1, designing at least 100 groups of interactive experimental schemes by taking the maximum and minimum burial depth, the coal thickness and the distance between the coal bed and the main key layer as two end limits;
s3, respectively calculating the disturbance degree of the main key layer of each experimental scheme: numbering rock stratums from the first layer of the drilled holes to the overlying main key layer from bottom to top by a formula
Calculating the disturbance degree of a main key layer of each drill hole, wherein lambda represents a main key layer disturbance degree influence factor, C represents a geological constant,
representing the disturbance degree of the main key layer, M representing the thickness of the coal bed, H representing the thickness of a rock stratum between the main key layer and the coal bed, D representing the buried depth of the coal bed, H 1 ，H 2 ···H i Indicating the thickness, R, of the i-th rock above the seam of coal being mined t1 ，R t2 ···R ti Indicating tensile strength, K, of the i-th rock above the seam being mined p1 ，K p2 ···K pi Expressing the crushing and swelling coefficient of the ith layer of rock above the mined coal bed, expressing lambda as a main key layer disturbance degree influence factor, expressing C as a geological constant, and entering the step S5;
s4, calculating the subsidence of the main key layer of each experimental scheme: screening and sorting out the maximum and minimum coal thickness, burial depth and interlayer spacing according to the interactive experimental scheme designed in the step S2, designing a response surface experimental scheme by taking the maximum coal thickness, burial depth and interlayer spacing as horizontal '1', taking the minimum coal thickness, burial depth and interlayer spacing as horizontal '-1', taking the mean value between the maximum value and the minimum value as horizontal '0', and designing a response surface experimental scheme by using the FLAC 3D Performing numerical simulation on the designed response surface scheme by using numerical simulation software, extracting the maximum subsidence of the main key layer as a response value to obtain a main key layer subsidence response surface calculation formula, performing main key layer subsidence calculation on the rest experiment schemes, and continuing to perform the step S5;
s5, determining the coupling relation between the disturbance degree and the subsidence of the main key layer: respectively normalizing the disturbance degree and the subsidence of the main key layer in the experimental data, obtaining a scatter diagram by taking the disturbance degree of the normalized main key layer as an abscissa and the subsidence of the normalized main key layer as an ordinate, performing linear fitting on scatter points to obtain a coupling relation between the disturbance degree and the subsidence normalized value of the main key layer, and continuing to perform the step S6;
s6, disturbance grade division: based on the coupling relation between the normalized main key layer disturbance degree and the sinkage, performing main key layer disturbance influence grade division, and continuing to perform step S7;
s7, calculating the disturbance degree normalization values of all the main key layers of the drill holes in the region: calculating the disturbance degrees of the main key layers of all the drill holes according to the calculation formula of the disturbance degrees of the main key layers in the step S3, and normalizing the disturbance degrees of the main key layers of all the drill holes in the area according to the calculation formula of the normalized value of the disturbance degrees of the main key layers in the step S5;
s8, guiding engineering practice: the method comprises the steps of matching the coordinates of a drill hole with the disturbance degree values of the main key layers one by one, establishing a corresponding database through ArcMap, obtaining a distribution map of the disturbance degree normalized values of the main key layers in an area by adopting a spline function interpolation method in software, finally analyzing the relative size of the main key layers affected by coal mining disturbance and the variation trend of the subsidence of the main key layers in the area according to the grade distribution of the disturbance degree normalized values of the main key layers in the area, and taking advance preventive measures for the strong disturbance area.
In step S3, the calculation formula of the disturbance factor of the main key layer in the area is:
in the formula: h max Represents the maximum distance, M, of the in-zone mined coal seam from the main key seam min Indicating the minimum thickness of the seam of coal mined in the zone, D min Indicating the minimum depth of burial, K, of the coal seam in the zone Pimax Representing the maximum coefficient of crushing and expansion, R, in a single formation within a zone timax Indicates the maximum tensile strength, H, in a single formation in a zone imax Represents the maximum thickness, K, of a single formation within a zone Pimin Representing the minimum coefficient of crushing and expansion, R, in a single formation within a zone timin Indicates the minimum tensile strength, H, in a single formation within a zone imin The minimum thickness of a single formation in the zone is indicated, and H represents the distance between the coal seam and the main critical zone in the borehole when the single formation in the zone has the maximum or minimum thickness.
In step S3, the calculation formula of the regional intrinsic constants is:
the symbols in the formula represent the influence factors of the disturbance degree of the same main key layer.
In step S5, the main key layer disturbance degree normalization calculation formula is:
wherein w represents a disturbance degree normalization value of the main key layer;
representing the maximum value of the disturbance degree of the main key layer in the region;
representing the minimum value of the disturbance degree of the main key layer in the region;
representing the primary key layer perturbation value of the borehole i.
In step S5, the main key layer subsidence normalization calculation formula is:
wherein eta represents the sinkage normalized value of the main key layer; y is max Representing the maximum value of the subsidence of the main key layer in the region; y is min Representing the minimum value of the subsidence of the main key layer in the region; y is i Representing the main critical layer subsidence of the borehole i.
In step S6, the disturbance influence level of the main key layer may be divided into five levels according to the obtained normalized value of the disturbance degree of the main key layer, and the division criteria of the five levels are as follows: w is 0.8-1.0 (large subsidence in strong disturbance area), w is 0.6-0.8 (large subsidence in strong disturbance area), w is 0.4-0.6 (medium subsidence in medium disturbance area, w is 0.2-0.4 (small subsidence in weak disturbance area), w is 0.0-0.2 (small subsidence in weak disturbance area).
In step S8, after the prediction and evaluation result guides the field production, the coupling relationship and the disturbance influence level partitioning criterion may be optimized according to the field practice data.
Compared with the prior art, the invention has the following beneficial effects:
(1) the comprehensive evaluation method for disturbance influence of the main key layer based on the subsidence can qualitatively describe the degree of disturbance influence of the main key layer on coal seam mining.
(2) The practical application of the invention can effectively predict the variation trend of the subsidence of the main key layer in the area, evaluate the intensity of disturbance influence on the main key layer in the area, and provide reliable reference for adjusting and reducing the disturbance influence on the main key layer, improving the stability of the main key layer, strengthening rock stratum control and ensuring the safe production of a working face.
Drawings
FIG. 1 is a flow chart of an implementation of a method for comprehensively evaluating disturbance influence of a main key layer.
FIG. 2 is a schematic diagram of implementation steps of the main key layer disturbance influence comprehensive evaluation method.
FIG. 3 is a schematic illustration of a borehole formation numbering according to an embodiment of the invention.
FIG. 4 is a diagram illustrating a distribution of perturbation degree normalization values of the main key layer according to an embodiment of the present invention.
In fig. 3: a, drilling hole numbering, B, main key layer and E, coal layer; in fig. 4: s-working face number, G-working face direction of advance, Z-bore, J 1 Regions of strong disturbance, J 2 Regions of greater disturbance, J 3 Moderate disturbance zone, J 4 Zone of weaker perturbation, J 5 -a weakly disturbed region.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below, and it is obvious that the described embodiments are some embodiments of the present invention, but not all embodiments; all other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
As shown in fig. 1, the invention provides a method for comprehensively evaluating disturbance influence of a main key layer based on subsidence, which comprises the following steps:
s1, collecting, sorting and analyzing the geological drilling data on site: the method comprises the steps of numbering drill holes by collecting geological drill hole data of a whole mining area or a local mining area, calculating a main key layer position of each drill hole, counting information such as coordinates of each drill hole, the thickness of a first layer of coalable layer, the burial depth of a coal layer, the distance between the coal layer and the main key layer, physical and mechanical parameters (including tensile strength, crushing and swelling coefficients and the like) of a stratum between the coal layer and the main key layer, the thickness of each stratum between layers and the like, and finally counting and analyzing the maximum value and the minimum value of relevant factors such as the thickness of the coalable layer, the burial depth of the coal layer, the distance between the coal layer and the main key layer in the area;
s2, respectively taking the maximum and minimum burial depth, the coal thickness and the distance between the coal bed and the main key layer as two end limits, and designing at least 100 groups of interactive experimental schemes;
s3, numbering rock stratums from the first layer of the coal seam to the overlying main key layer of each drill hole from bottom to top, and calculating the disturbance degree of the main key layer of each drill hole, wherein the calculation formula of the disturbance degree of the main key layer is as follows:
in the formula (1), lambda represents a main key layer disturbance degree influence factor, C represents a geological constant,
representing the disturbance degree of the main key layer, M representing the thickness of the coal bed, H representing the thickness of a rock stratum between the main key layer and the coal bed, D representing the buried depth of the coal bed, H 1 ，H 2 ···H i Indicating the thickness, R, of the i-th rock above the seam of coal being mined t1 ，R t2 ···R ti Denotes the tensile strength, K, of the i-th rock above the mined coal seam p1 ，K p2 ···K pi The coefficient of crushing and expansion of the ith layer of rock above the mined coal seam is shown,lambda represents a main key layer disturbance degree influence factor, C represents a geological constant, and the step S5 is carried out;
the calculation formula of the disturbance degree influence factor of the main key layer in the area is as follows:
in formula (2): h max Represents the maximum distance, M, of the in-zone mined coal seam from the main key seam min Indicating the minimum thickness of the seam of coal mined in the zone, D min Indicating the minimum depth of burial, K, of the coal seam in the zone Pimax Representing the maximum coefficient of crushing and expansion, R, in a single formation within a zone timax Indicates the maximum tensile strength, H, in a single formation in a zone imax Represents the maximum thickness, K, of a single formation within a zone Pimin Representing the minimum coefficient of crushing and expansion, R, in a single formation within a zone timin Indicates the minimum tensile strength, H, in a single formation within a zone imin The minimum thickness of a single formation in the zone is indicated, and H represents the distance between the coal seam and the main critical zone in the borehole when the single formation in the zone has the maximum or minimum thickness.
The calculation formula of the geological constant in the region is as follows:
in the formula (3), the symbols represent the same formula (2).
S4, according to the interactive experiment scheme designed in the step S2, the maximum and minimum coal thickness, the burial depth and the interlayer spacing are screened and sorted out, a response surface experiment scheme is designed with the maximum coal thickness, the burial depth and the interlayer spacing as horizontal '1', the minimum coal thickness, the burial depth and the interlayer spacing as horizontal '-1', the mean value between the maximum value and the minimum value as horizontal '0', and the FLAC is utilized 3D The numerical simulation software carries out numerical simulation on the designed response surface scheme, extracts the maximum sinking amount of the main key layer as a response value to obtain a calculation formula of the sinking amount of the main key layer, and carries out calculation of the sinking amount of the main key layer on the rest experiment schemesContinuing with step S5;
s5, respectively carrying out normalization processing on the disturbance degree and the subsidence of the main key layer in the experimental data, taking the disturbance degree of the normalized main key layer as an abscissa and the subsidence of the normalized main key layer as an ordinate to obtain a scatter diagram, carrying out linear fitting on scatter points to obtain a coupling relation between the disturbance degree and the subsidence normalization value of the main key layer, and continuing to carry out the step S6;
the normalized calculation formula of the disturbance degree of the main key layer is as follows:
wherein w represents a disturbance degree normalization value of the main key layer;
representing the maximum value of the disturbance degree of the main key layer in the region;
representing the minimum value of the disturbance degree of the main key layer in the region;
representing the primary key layer perturbation value of the borehole i.
The main key layer subsidence quantity normalization calculation formula is as follows:
wherein eta represents the sinkage normalized value of the main key layer; y is max Representing the maximum value of the subsidence of the main key layer in the region; y is min Representing the minimum value of the subsidence of the main key layer in the region; y is i Representing the main critical layer subsidence of the borehole i.
S6, based on the coupling relation between the normalized main key layer disturbance degree and the sinkage, carrying out main key layer disturbance influence grade division, and continuing to carry out the step S7;
s7, calculating the main key layer disturbance degree of all the drill holes according to the main key layer disturbance degree calculation formula in the step S3, and normalizing the main key layer disturbance degree of all the drill holes in the area according to the main key layer disturbance degree normalization value calculation formula in the step S5;
s8, matching the coordinates of the drilled holes with the disturbance degree values of the main key layers one by one, establishing a corresponding database through ArcMap, obtaining a distribution diagram of the disturbance degree normalized values of the main key layers in the area by adopting a spline function interpolation method in software, and finally analyzing the relative size of the main key layers in the area influenced by the coal seam mining disturbance and the change trend of the subsidence of the main key layers according to the level distribution of the disturbance degree normalized values of the main key layers in the area to take advance preventive measures on the strong disturbance area.
The comprehensive evaluation result can be divided into five grades of a strong disturbance area, a medium disturbance area, a weak disturbance area and a weak disturbance area, and the five grades are divided according to the following standard:
a strong disturbance area: the sedimentation amount in the area is large, and w is 0.8-1.0;
a strong disturbance area: the settlement in the area is large, and w is 0.6-0.8;
medium disturbance area: the settlement amount in the area is medium, and w is 0.4-0.6;
weaker perturbation zone: the settlement in the area is small, and w is 0.2-0.4;
and (3) a weak disturbance area: the sedimentation amount in the area is small, and w is 0.0-0.2.
The following description will explain embodiments of the present invention by taking a coal mine as an example.
The coal seam occurrence condition difference of a certain mining area 23 is large, the coal seam thickness change is large, the coal seam thickness is 1-5 m, the whole coal seam is a slowly inclined coal seam, the burial depth is 300-700 m, the coal seam is 30-70 m away from a main key layer, the lithology of the main key layer is mainly sandstone and is influenced by a water-bearing layer of a top direct-compass group, the roof is subjected to dripping, sprinkling and water inrush during the stoping of a working face, and a plurality of problems seriously influence the economic benefit of a coal mine enterprise, so the disturbance influence of the coal seam mining on the main key layer in the area with the large coal seam occurrence change is explored, the site production is guided, and specific implementation steps are as follows aiming at the conditions:
1. according to geological report data, 12 drill holes on 3 exploration lines of a 23 mining area of the mine are selected, main key layer positions of all the drill holes are calculated respectively, information such as coordinates of all the drill holes, the thickness of a coal layer to be mined, the burial depth of the coal layer, the distance between the coal layer and the main key layer, physical and mechanical parameters (including tensile strength, crushing and expansion coefficients and the like) of rock layers between the coal layer and the main key layer, the thickness of each rock layer between the layers and the like is calculated, and finally the maximum value and the minimum value of relevant factors such as the thickness of the coal layer to be mined, the burial depth of the coal layer, the distance between the coal layer and the main key layer in the area are screened out.
2. According to geological drilling data, five groups of coal beds with the thickness of 1m, 2m, 3m, 4m and 5m, five groups of buried depths of 300m, 400m, 500m, 600m and 700m and five groups of coal beds and main key layer intervals of 30m, 40m, 50m, 60m and 70m are selected, and 125 groups of three-factor five-level interaction experiments are obtained.
3. According to the step S3, numbering is carried out on rock strata of each drilling hole from bottom to top, and related data of each drilling hole are substituted, so that the main key layer disturbance degree value of each experimental scheme can be obtained through a formula (1-3).
4. According to the maximum coal thickness of 5m and the minimum coal thickness of 1 m; the maximum buried depth is 700m, and the minimum buried depth is 300 m; the maximum interlayer spacing is 70m and the minimum interlayer spacing is 30m, and a three-factor three-level response surface experimental scheme is designed and passes through FLAC 3D Performing simulation calculation, screening out the maximum vertical displacement of the main key layer of each model, and obtaining a response plane regression model calculation formula of three factors and the subsidence of the main key layer:
in the formula: y represents the subsidence of the main key layer, M represents the coal thickness, D represents the buried depth, and H represents the distance between the coal layer and the main key layer.
And respectively calculating the main key layer subsidence in the rest schemes by the formula (6).
5. And (4) carrying out normalization processing on the disturbance degree and the subsidence obtained by calculating the interactive experimental scheme according to the formula (4-5), and analyzing and fitting the coupling relation between the disturbance degree and the subsidence of the main key layer.
6. According to the coupling relation between the disturbance degree and the subsidence of the main key layer, the disturbance degrees of the main key layer of 12 drill holes are calculated by the formula (1) and are respectively as follows: 4.124, 4.028, 3.429, 2.816, 2.180, 4.499, 2.887, 4.046, 2.056, 1.767, 1.452 and 1.731, and the perturbation degree of each main key layer of the drill hole is normalized by the formula (4) to obtain the corresponding: 0.88, 0.85, 0.65, 0.45, 0.24, 1.00, 0.47, 0.85, 0.20, 0.10, 0.00, 0.09.
7. And matching the coordinates of the drilled holes with the disturbance degree values of the main key layer one by one, establishing a corresponding database through ArcMap, and obtaining the graph 4 by adopting a spline function interpolation method in software.
In conclusion, the disturbance degree of the main key layer in the southeast area of the mining area 23 is more severely changed compared with the northwest area, the disturbance influence on the main key layer is increased along the advancing direction of the working face, the sinking increment of the main key layer is relatively increased, and precautionary measures should be taken in advance for the strong disturbance area.
Finally, it should be noted that: the above embodiments are only used to illustrate the technical solution of the present invention, and not to limit the same; while the invention has been described in detail and with reference to the foregoing embodiments, it will be understood by those skilled in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some or all of the technical features may be equivalently replaced; and the modifications or the substitutions do not make the essence of the corresponding technical solutions depart from the scope of the technical solutions of the embodiments of the present invention.
Claims (5)
1. A settlement-based main key layer disturbance influence comprehensive evaluation method is characterized by comprising the following steps:
s1, numbering the drill holes by collecting geological drill hole data of the whole mining area or a local mining area, calculating a main key layer position of each drill hole, counting information such as coordinates of each drill hole, the thickness of a first layer of coalbed methane layer, the coal seam burial depth, the interval between the coal seam and the main key layer, physical and mechanical parameters (including tensile strength, crushing and expansion coefficients and the like) of a coal seam and a main key interlayer rock layer, the thickness of each interlayer rock layer and the like, finally counting and analyzing the maximum value and the minimum value of relevant factors such as the thickness of the coalbed methane layer, the coal seam burial depth, the interval between the coal seam and the main key layer and the like in the area, and performing the step S2;
s2, designing at least 100 groups of interactive experimental schemes by taking the maximum and minimum burial depth, the coal thickness and the distance between the coal bed and the main key layer as two end limits according to the related drilling data arranged in the step S1;
s3, numbering rock stratums from the first layer of the coal seam to the overlying main key layer of each drill hole from bottom to top, respectively calculating the disturbance degree of the main key layer of each experimental scheme, and entering the step S5;
s4, the screening part scheme utilizes the response surface to combine numerical simulation to obtain a main key layer subsidence response surface calculation formula, main key layer subsidence calculation is carried out on the rest experiment schemes, and the step S5 is carried out;
s5, normalizing the disturbance degree and the subsidence of the main key layer in the experimental data to obtain the coupling relation of the disturbance degree and the subsidence normalization value of the main key layer;
s6, disturbance grade division: based on the coupling relation between the normalized main key layer disturbance degree and the sinkage, performing main key layer disturbance influence grade division, and continuing to perform step S7;
s7, calculating the disturbance degree normalization values of all drilling main key layers in the region;
s8, matching the coordinates of the drilled holes with the disturbance degree values of the main key layer one by one, establishing a corresponding database through ArcMap, and obtaining a main key layer disturbance degree normalized value distribution map in the area by adopting a spline function interpolation method in software.
2. The method for comprehensively evaluating the disturbance influence of the main key layer based on the subsidence according to claim 1, wherein in the step S3, the calculation formula of the disturbance degree of the main key layer is as follows:
in the formula, lambda represents a main key layer disturbance degree influence factor, C represents a geological constant,
representing the disturbance degree of the main key layer, M representing the thickness of the coal bed, H representing the thickness of a rock stratum between the main key layer and the coal bed, D representing the buried depth of the coal bed, H 1 ，H 2 ···H i Indicating the thickness, R, of the i-th rock above the seam of coal being mined t1 ，R t2 ···R ti Indicating tensile strength, K, of the i-th rock above the seam being mined p1 ，K p2 ···K pi And the coefficient of crushing and expansion of the ith layer of rock above the mined coal bed is represented, lambda represents a main key layer disturbance degree influence factor, and C represents a geological constant.
The calculation formula of the disturbance degree influence factor of the main key layer in the area is as follows:
in the formula: h max Represents the maximum distance, M, of the in-zone mined coal seam from the main key seam min Indicating the minimum thickness of the seam of coal mined in the zone, D min Represents the minimum burial depth of the coal seam in the region,
representing the maximum coefficient of fracture and expansion in a single formation within a zone,
denotes the maximum tensile strength, H, in a single formation within a zone imax Representing the maximum thickness of a single formation within the region,
indicating the minimum coefficient of fracture and expansion in a single formation within the zone,
indicates the minimum tensile strength, H, in a single formation within a zone imin The minimum thickness of a single formation in the zone is indicated, and H represents the distance between the coal seam and the main critical zone in the borehole when the single formation in the zone has the maximum or minimum thickness.
The calculation formula of the geological constant in the region is as follows:
the symbols in the formula represent the influence factors of the disturbance degree of the same main key layer.
3. The method for comprehensively evaluating the disturbance influence of the main key layer based on the subsidence according to claim 1, wherein in the step S4, the calculation formula of the subsidence of the main key layer is obtained by combining a response surface with a numerical simulation mode.
4. The method for comprehensively evaluating the disturbance influence of the main key layer based on the subsidence according to claim 1, wherein in the step S5, the main key layer disturbance degree normalization calculation formula is:
wherein w represents a disturbance degree normalization value of the main key layer;
representing the maximum value of the disturbance degree of the main key layer in the region;
representing the minimum value of the disturbance degree of the main key layer in the region;
representing the primary key layer perturbation value of the borehole i.
The main key layer subsidence quantity normalization calculation formula is as follows:
wherein eta represents the sinkage normalized value of the main key layer; y is max Representing the maximum value of the subsidence of the main key layer in the region; y is min Representing the minimum value of the subsidence of the main key layer in the region; y is i Representing the main critical layer subsidence of the borehole i.
5. The method for comprehensively evaluating the disturbance influence of the main key layer based on the subsidence according to claim 1, wherein in the step S6, the disturbance influence level of the main key layer can be divided into five levels according to the obtained normalized value of the disturbance degree of the main key layer, and the division criteria of the five levels are as follows: w is 0.8-1.0 (large subsidence in strong disturbance area), w is 0.6-0.8 (large subsidence in strong disturbance area), w is 0.4-0.6 (medium subsidence in medium disturbance area, w is 0.2-0.4 (small subsidence in weak disturbance area), w is 0.0-0.2 (small subsidence in weak disturbance area).
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210530296.2A CN114856706B (en) | 2022-05-16 | 2022-05-16 | Method for comprehensively evaluating disturbance influence of main key layer based on subsidence |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210530296.2A CN114856706B (en) | 2022-05-16 | 2022-05-16 | Method for comprehensively evaluating disturbance influence of main key layer based on subsidence |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114856706A true CN114856706A (en) | 2022-08-05 |
| CN114856706B CN114856706B (en) | 2023-04-11 |
Family
ID=82637135
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210530296.2A Active CN114856706B (en) | 2022-05-16 | 2022-05-16 | Method for comprehensively evaluating disturbance influence of main key layer based on subsidence |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114856706B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115600401A (en) * | 2022-10-17 | 2023-01-13 | 中国矿业大学（北京）（Cn） | Lower key layer stability evaluation method for maintaining safety of coal mine underground reservoir |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030036851A1 (en) * | 2001-08-07 | 2003-02-20 | Stankus John C. | Method of roof instability rating |
| CN112711847A (en) * | 2020-12-28 | 2021-04-27 | 西安科技大学 | Method for determining surface subsidence coefficients of key layer at different positions of overlying strata |
| CN113339073A (en) * | 2021-07-06 | 2021-09-03 | 中国矿业大学 | Impact risk evaluation method based on roof rock stratum structure |
| CN113435014A (en) * | 2021-06-04 | 2021-09-24 | 华北水利水电大学 | Dynamic prediction method for mining overburden rock movement deformation |
| CN113486517A (en) * | 2021-07-07 | 2021-10-08 | 安徽理工大学 | Ground control method and device for mining disasters in coal mining area |
| CN114109506A (en) * | 2021-11-26 | 2022-03-01 | 中国矿业大学 | Coal mine earthquake risk assessment method |
-
2022
- 2022-05-16 CN CN202210530296.2A patent/CN114856706B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030036851A1 (en) * | 2001-08-07 | 2003-02-20 | Stankus John C. | Method of roof instability rating |
| CN112711847A (en) * | 2020-12-28 | 2021-04-27 | 西安科技大学 | Method for determining surface subsidence coefficients of key layer at different positions of overlying strata |
| CN113435014A (en) * | 2021-06-04 | 2021-09-24 | 华北水利水电大学 | Dynamic prediction method for mining overburden rock movement deformation |
| CN113339073A (en) * | 2021-07-06 | 2021-09-03 | 中国矿业大学 | Impact risk evaluation method based on roof rock stratum structure |
| CN113486517A (en) * | 2021-07-07 | 2021-10-08 | 安徽理工大学 | Ground control method and device for mining disasters in coal mining area |
| CN114109506A (en) * | 2021-11-26 | 2022-03-01 | 中国矿业大学 | Coal mine earthquake risk assessment method |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115600401A (en) * | 2022-10-17 | 2023-01-13 | 中国矿业大学（北京）（Cn） | Lower key layer stability evaluation method for maintaining safety of coal mine underground reservoir |
| CN115600401B (en) * | 2022-10-17 | 2023-09-08 | 中国矿业大学(北京) | Lower key layer stability evaluation method for maintaining safety of underground coal mine reservoir |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114856706B (en) | 2023-04-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111691872B (en) | Dynamic and accurate detection method for development range and degree of coal seam overlying strata fracture | |
| CN105986816A (en) | Method for recognizing sweet spots in shale stratum | |
| CN103472484A (en) | Horizontal well track optimization method based on RS three-dimensional sensitivity seismic attribution analysis | |
| CN110359895B (en) | Water exploration and drainage method for staged fracturing of horizontal well of heterogeneous huge and thick sandstone | |
| WO2023000257A1 (en) | Geological-seismic three-dimensional prediction method for favorable metallogenic site of sandstone-type uranium deposit | |
| Yang et al. | Analysis of abnormally high water production in coalbed methane vertical wells: A case study of the Shizhuangnan block in the southern Qinshui Basin, China | |
| CN114856706B (en) | Method for comprehensively evaluating disturbance influence of main key layer based on subsidence | |
| CN111897002B (en) | Roof pre-splitting measure effect evaluation method based on microseismic monitoring | |
| CN110851991B (en) | Underground water flow numerical simulation method | |
| CN115422740A (en) | Method for predicting height of water flowing fractured zone of layered fully-mechanized caving mining of huge thick coal seam | |
| Stiles et al. | Reservoir management in the means San Andres Unit | |
| Chen et al. | Geology-engineering integration to improve drilling speed and safety in ultra-deep clastic reservoirs of the Qiulitage structural belt | |
| CN113622952A (en) | Comprehensive treatment method for water damage of mountain green coal in structurally complex area | |
| CN113435014A (en) | Dynamic prediction method for mining overburden rock movement deformation | |
| Male et al. | Statistical analysis of fractures from the hydraulic fracture test site 1 | |
| Casabianca et al. | The Machar Oil Field: waterflooding a fractured chalk reservoir | |
| Abolfazlzadeh | Application of seismic monitoring in caving mines | |
| CN115374611A (en) | Method for judging tunnel soft rock large deformation disaster risk level under complex geological condition | |
| CN115236740A (en) | Water-flowing fractured zone development height prediction method based on seismic interpretation achievement | |
| Liu et al. | Study on stability of underlying room and pillar old goaf in close coal seam and mining of the upper coal seam | |
| Karacaer et al. | Integrated numerical modeling of a Huff-N-Puff pilot in Eagle Ford | |
| Jiang et al. | Geomechanical approach lead to successful Implementation of the first highly deviated well at kuqa depression in Tarim Basin | |
| CN116910886B (en) | Impact energy early warning method and system for karst region | |
| Gao | Research on infill adjustment method for XX narrow and thin sand body oilfield | |
| CN116611567B (en) | Mining area mining overlying strata roof water inrush composite disaster risk partition prediction method |
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 | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Li Yang Inventor after: Wang Nan Inventor after: Lei Xinghai Inventor after: Li Tiezheng Inventor after: Ren Yuqi Inventor after: Jin Xiangyang Inventor before: Li Yang Inventor before: Wang Nan Inventor before: Lei Xinghai Inventor before: Li Tiezheng Inventor before: Ren Yuqi Inventor before: Jin Xiangyang |