CN118607238B - Method, device and equipment for estimating grouting quantity of surrounding rock cracks of repeated mining roadway - Google Patents
Method, device and equipment for estimating grouting quantity of surrounding rock cracks of repeated mining roadway Download PDFInfo
- Publication number
- CN118607238B CN118607238B CN202410761924.7A CN202410761924A CN118607238B CN 118607238 B CN118607238 B CN 118607238B CN 202410761924 A CN202410761924 A CN 202410761924A CN 118607238 B CN118607238 B CN 118607238B
- Authority
- CN
- China
- Prior art keywords
- area
- damage
- roadway
- mining
- fracture
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000005065 mining Methods 0.000 title claims abstract description 130
- 239000011435 rock Substances 0.000 title claims abstract description 80
- 238000000034 method Methods 0.000 title claims abstract description 51
- 230000006378 damage Effects 0.000 claims abstract description 152
- 238000011161 development Methods 0.000 claims abstract description 94
- 239000003245 coal Substances 0.000 claims abstract description 65
- 238000012544 monitoring process Methods 0.000 claims abstract description 64
- 238000012937 correction Methods 0.000 claims abstract description 38
- 230000008859 change Effects 0.000 claims abstract description 34
- 230000005641 tunneling Effects 0.000 claims abstract description 29
- 238000005553 drilling Methods 0.000 claims description 27
- 238000004364 calculation method Methods 0.000 claims description 22
- 238000006073 displacement reaction Methods 0.000 claims description 12
- 230000006870 function Effects 0.000 claims description 10
- 238000004519 manufacturing process Methods 0.000 abstract description 11
- 238000010586 diagram Methods 0.000 description 20
- 208000008918 voyeurism Diseases 0.000 description 11
- 238000004088 simulation Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 8
- 238000004590 computer program Methods 0.000 description 7
- 238000003860 storage Methods 0.000 description 5
- JHGSLSLUFMZUMK-UHFFFAOYSA-N [2-hydroxy-2-(4-hydroxyphenyl)ethyl]-[4-(4-hydroxyphenyl)butan-2-yl]azanium;chloride Chemical compound Cl.C=1C=C(O)C=CC=1C(O)CNC(C)CCC1=CC=C(O)C=C1 JHGSLSLUFMZUMK-UHFFFAOYSA-N 0.000 description 4
- 238000012423 maintenance Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000002457 bidirectional effect Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000004873 anchoring Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 230000009172 bursting Effects 0.000 description 1
- 238000013499 data model Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000003801 milling Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000002787 reinforcement Effects 0.000 description 1
- 238000007789 sealing Methods 0.000 description 1
- 230000008961 swelling Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Civil Engineering (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Architecture (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
本发明提供了一种重复采动巷道围岩裂隙注浆量的估计方法，包括根据裂隙破坏范围的监测数据计算裂隙理论破坏范围面积；根据裂隙发育深度的监测数据计算平均发育深度，根据平均发育深度计算巷道内空间的总实际等效破坏面积；根据每次采动或掘进前后总实际等效破坏面积的变化量除以理论破坏范围面积的变化量得到修正系数，根据修正系数对理论破坏范围面积进行修正得到修正后的破坏范围面积，根据修正后的破坏面积得到每次采动或掘进前后修正后破坏面积变化量；根据修正后破坏面积变化量乘以巷道长度得到注浆量。本方案解决了现有技术中数值模的方法预测的裂隙变化量不够精确，从而导致巷道围岩支护工作效率低下影响了煤炭生产安全和生产效率的问题。
The present invention provides a method for estimating the amount of grouting in the surrounding rock fissures of a repeated mining tunnel, including calculating the theoretical fissure damage range area according to the monitoring data of the fissure damage range; calculating the average development depth according to the monitoring data of the fissure development depth, and calculating the total actual equivalent damage area of the tunnel space according to the average development depth; obtaining a correction coefficient by dividing the change in the total actual equivalent damage area before and after each mining or tunneling by the change in the theoretical damage range area, correcting the theoretical damage range area according to the correction coefficient to obtain a corrected damage range area, and obtaining a corrected damage area change before and after each mining or tunneling according to the corrected damage area; and obtaining the grouting amount by multiplying the corrected damage area change by the tunnel length. This solution solves the problem that the change in fissure predicted by the numerical model method in the prior art is not accurate enough, which leads to low efficiency of tunnel surrounding rock support work and affects coal production safety and production efficiency.
Description
Technical Field
The invention belongs to the technical field of roadway support safety, and particularly relates to a method, a device and equipment for estimating grouting quantity of surrounding rock cracks of a repeated mining roadway.
Background
Along with the wide application of the double-roadway tunneling technology in longwall mining, in order to prevent roof fall and water bursting accidents caused by unfavorable geological structures of coal mines, so that coal production safety and production efficiency are guaranteed, two stoping roadways are simultaneously dug when the stoping roadway of the working face is tunneling, one roadway serves the working face, the stoping roadway is gradually abandoned along with the working face, the other roadway is left to serve the next working face, and therefore the remaining roadways can be influenced by repeated mining, stress of rocks around the roadways is reissued, surrounding rocks of the roadways are damaged under the action of high stress, cracks are generated inside the surrounding rocks, and the expansion and penetration of the cracks are needed to be predicted.
In recent years, a great deal of research and discussion is carried out on repeated mining tunnel by a plurality of scholars, but the estimated research on the grouting amount of the repeated mining tunnel is still very few, at present, the grouting amount of the tunnel is predicted by theoretical calculation by a plurality of scholars, the method predicts the crack amount of surrounding rock of the repeated mining tunnel by combining theory and numerical simulation in the early period of mining, however, the crack variation predicted by the numerical simulation method is not accurate enough, and when an emergency is encountered in the actual exploitation process, the theoretical calculation amount is greatly different from the actual crack amount, so that the working efficiency of the roadway surrounding rock support is low, and the production safety and the production efficiency of coal are affected.
Disclosure of Invention
In order to solve the problems that in the prior art, the crack variation predicted by a numerical mode method is not accurate enough, and the theoretical calculation amount and the actual crack amount are larger when an emergency is encountered in the actual exploitation process, so that the working efficiency of roadway surrounding rock support is low, and the coal production safety and the production efficiency are affected, the invention provides a method, a device and equipment for estimating the crack grouting amount of repeated exploitation roadway surrounding rock.
In order to achieve the above object, the present invention provides the following technical solutions:
The method for estimating the grouting quantity of the surrounding rock cracks of the repeated mining roadway comprises the following steps:
Constructing a repeated mining roadway surrounding rock numerical model, acquiring monitoring data of crack damage ranges of roadway surrounding rock in each repeated mining stage according to the numerical model, and acquiring monitoring data of crack development depth according to on-site monitoring stations of the roadway;
Calculating average development depth of the cracks according to the monitoring data of the development depth of the cracks before and after each mining or tunneling, and calculating total actual equivalent damage area of the space in the roadway according to the average development depth;
obtaining a correction coefficient according to the ratio of the change amount of the total actual equivalent damage area before and after each mining or tunneling to the change amount of the theoretical damage area, correcting the theoretical damage area according to the correction coefficient to obtain the corrected damage area, and obtaining the change amount of the damage area before and after each mining or tunneling according to the difference between the corrected damage area and the corrected damage area before and after each mining or tunneling;
and multiplying the corrected damage area change by the roadway length to obtain grouting quantity.
Further, the step of calculating the fracture theoretical damage range area according to the monitoring data of the fracture damage range comprises the following steps:
And establishing a two-dimensional plane rectangular coordinate system according to the monitoring data of the fracture damage range, fitting a piecewise function curve of the boundary line according to the point data on the boundary line of the fracture damage range, and carrying out integral summation on the fracture damage range according to the piecewise function curve to obtain the area of the fracture theoretical damage range.
Further, the step of calculating the total actual equivalent destruction area of the space in the roadway according to the average development depth comprises the following steps:
The crack development depth comprises crack development depths on four surfaces of a top plate, a bottom plate, a coal pillar upper and a coal wall upper of a roadway;
the relation expression of the average development depth of the fissures is as follows:
Wherein j= { top, bottom, column and wall }, i is the index of a monitoring point drilling, L j′i is the maximum value of the development depth of the ith drilling crack on the top plate, the bottom plate, the coal pillar upper or the coal wall upper, and X j is the displacement, namely the deformation, of the inner surface of the top plate, the bottom plate, the coal pillar upper or the coal wall upper roadway;
the related formula for calculating the total actual equivalent destruction area of the space in the roadway according to the average development depth is as follows:
Wherein, Is the total actual equivalent failure area of the device,Is the actual equivalent damage area of the top plate, the bottom plate, the coal pillar side or the coal wall side, n is the number of the maximum allowed holes, R is the radius of the holes, and the smallest distance between the two holes with the same radius is half of the smallest mutual influence.
Further, the correction coefficient calculation formula is as follows:
Wherein, Is the total equivalent destruction area of the current mining influence stage,Is the total equivalent damage area of the last mining influence stage, S 3 is the theoretical damage range area of the current mining influence stage, and S 2 is the theoretical damage range area of the last mining influence stage.
Further, the step of correcting the theoretical damage range area according to the correction coefficient to obtain a corrected damage range area includes:
If the correction coefficient sigma is less than or equal to1, the fracture theoretical damage range area does not need to be corrected, and if the correction coefficient sigma is more than 1, the fracture theoretical damage range area is corrected according to the correction coefficient sigma, and a specific correction formula is as follows:
S'2＝S2*σ
wherein S 2 is the fracture theoretical damage range area of the current mining influence stage, and S' 2 is the damage range area corrected by the current mining influence stage.
Further, the grouting amount calculation formula is as follows:
δ1=ΔS*c
ΔS=S'2-S1
Wherein δ 1 is the grouting amount of the crack in the current mining influence stage, c represents the roadway length, Δs is the corrected damage area variation, S' 2 is the corrected damage area in the current mining influence stage, and S 1 is the corrected damage area in the last mining influence stage.
An estimation device of repeated mining roadway surrounding rock crack grouting amount comprises:
The data acquisition module is used for constructing a repeated mining roadway surrounding rock numerical model, acquiring monitoring data of crack damage ranges of roadway surrounding rock in each stage of repeated mining according to the numerical model, and acquiring monitoring data of crack development depth according to on-site monitoring stations of the roadway;
the area calculation module is used for calculating the area of the fracture theoretical damage range according to the monitoring data of the fracture damage range, calculating the average development depth of the fracture according to the monitoring data of the development depth of the fracture before and after each mining or tunneling, and calculating the total actual equivalent damage area of the space in the roadway according to the average development depth;
The correction module is used for obtaining a correction coefficient according to the ratio of the change amount of the total actual equivalent damage area before and after each mining or tunneling to the change amount of the theoretical damage area, correcting the theoretical damage area according to the correction coefficient to obtain a corrected damage area, and obtaining the change amount of the damage area after each mining or tunneling according to the difference between the corrected damage area and the corrected damage area before each mining or tunneling;
and the result acquisition module is used for multiplying the corrected damage area variation by the roadway length to obtain grouting quantity.
The computer equipment comprises a memory and a processor, wherein computer execution instructions are stored in the memory, and the processor executes the computer execution instructions stored in the memory so as to realize the estimation method of the repeated mining roadway surrounding rock crack grouting amount.
The method for estimating the grouting quantity of the surrounding rock cracks of the repeated mining roadway has the following beneficial effects:
The method comprises the steps of constructing a repeated mining roadway surrounding rock numerical model, acquiring monitoring data of crack damage ranges of roadway surrounding rock at each stage of repeated mining according to the numerical model, acquiring monitoring data of crack development depth according to on-site monitoring stations of the roadway, calculating the area of a crack theoretical damage range according to the monitoring data of the crack damage ranges, calculating the average development depth of the crack according to the monitoring data of the crack development depth before and after each mining or tunneling, calculating the total actual equivalent damage area of a space in the roadway according to the average development depth, obtaining a correction coefficient according to the change amount of the total actual equivalent damage area before and after each mining or tunneling divided by the change amount of the area of the theoretical damage range, correcting the area of the theoretical damage range according to the correction coefficient, obtaining the change amount of the damage area after each mining or tunneling according to the corrected damage area, and multiplying the change amount of the damage area after correction by the length of the roadway to obtain grouting amount.
According to the method, the actual data acquired by the roadway site monitoring station and the numerical model simulation theoretical data are combined to predict the change area of the crack quantity of the repeated mining roadway surrounding rock, the total actual equivalent damage area is used for correcting the crack theoretical damage range area, the grouting quantity is calculated, the change quantity of the corrected damage range area is more accurate, simple addition, subtraction, multiplication and division operation is only involved in the calculation process, complex theoretical calculation is not involved, the operation efficiency is improved, theoretical basis can be provided for predicting the grouting quantity in time, grouting is timely carried out on the crack of the surrounding rock, and then the stability of the surrounding rock is improved, so that the roadway surrounding rock supporting working efficiency and the coal production efficiency are improved, safety guarantee is provided for coal production safety, the anchoring performance of anchor rods and anchor cables is improved, an accurate range is provided for roadway surrounding rock supporting, and the method has important significance for roadway supporting safety.
Drawings
In order to more clearly illustrate the embodiments of the present invention and the design thereof, the drawings required for the embodiments will be briefly described below. The drawings in the following description are only some of the embodiments of the present invention and other drawings may be made by those skilled in the art without the exercise of inventive faculty.
FIG. 1 is a flow chart of a prediction method according to an embodiment of the method of the present invention;
FIG. 2 is a schematic diagram of a development curve of a roadway surrounding rock fracture at a mining stage according to an embodiment of the method of the present invention;
FIG. 3 is a schematic diagram of a development curve of a roadway surrounding rock fracture at the primary and secondary mining stages according to an embodiment of the method of the present invention;
FIG. 4 is a schematic diagram of a range of roadway surrounding rock fracture damage in a numerically simulated roadway driving stage according to an embodiment of the method of the present invention;
FIG. 5 is a schematic diagram of a range of fracture damage of roadway surrounding rock at a mining stage of numerical simulation according to an embodiment of the method of the present invention;
FIG. 6 is a schematic diagram of a range of fracture damage of roadway surrounding rock at a numerical simulation secondary mining stage according to an embodiment of the method of the present invention;
FIG. 7 is a schematic diagram of a two-dimensional rectangular coordinate system of roadway surrounding rock at a roadway driving stage according to an embodiment of the method of the invention;
FIG. 8 is a schematic diagram of a rectangular coordinate system of a roadway surrounding rock two-dimensional plane at the primary mining stage according to an embodiment of the method of the present invention;
FIG. 9 is a schematic diagram of a rectangular coordinate system of two-dimensional plane of roadway surrounding rock at the secondary mining stage according to an embodiment of the method of the present invention;
FIG. 10 is a schematic diagram of a casing deformation monitoring system according to an embodiment of the method of the present invention;
FIG. 11 is a schematic diagram of a system for monitoring deformation of a casing according to an embodiment of the present invention;
In the figure, 1, a sleeve, 2, a strain collector, 3, a data receiver, 4, a data processing center, 5, a monitoring and dispatching center and 6, a wire.
Detailed Description
The present invention will be described in detail below with reference to the drawings and the embodiments, so that those skilled in the art can better understand the technical scheme of the present invention and can implement the same. The following examples are only for more clearly illustrating the technical aspects of the present invention, and are not intended to limit the scope of the present invention.
Method embodiment
The invention provides an estimation method of repeated mining roadway surrounding rock crack grouting amount, which is specifically shown in fig. 1 and comprises the following steps:
S1, constructing a repeated mining roadway surrounding rock numerical model, acquiring monitoring data of crack damage ranges of roadway surrounding rock in each repeated mining stage according to the data model, and simulating the whole process of mining influence after the repeated mining roadway surrounding rock numerical model is constructed, recording the damage ranges of the roadway surrounding rock in each repeated mining stage, and establishing a mathematical formula to represent the damage ranges by areas.
And establishing a repeated mining stope roadway surrounding rock numerical model, simulating the whole process of primary mining and secondary mining, and recording the damage range of roadway surrounding rock in each stage of repeated mining.
And establishing a mathematical formula, defining the damage range of each stage of repeated mining by using the area, wherein the change amount of the area of each stage of repeated mining is the change amount of the crack.
And establishing a two-dimensional plane rectangular coordinate system according to the monitoring data of the fracture damage range, fitting a piecewise function curve of a boundary line according to the point data on the boundary line of the fracture damage range, carrying out integral summation on the fracture damage range according to the piecewise function curve to obtain fracture damage areas, and taking the fracture damage areas before and after each mining or tunneling as differences to obtain the variation of the fracture damage areas.
S2, establishing a roadway site monitoring station, and observing and recording the maximum value of the development depth of surrounding rock cracks affecting the whole process of mining.
Selecting a roadway surrounding rock crack monitoring station, peeping the monitoring station to the top and bottom plates and the two sides of the roadway, putting a high-strength sleeve into the peeping hole after the peeping hole is formed, protecting the shape of the drilling hole, and judging whether the sleeve needs to be taken out for maintenance according to a sleeve monitoring system.
And taking out the sleeve 1 in the monitoring point peeping borehole at each stage of the influence of primary mining and secondary mining, observing and recording the maximum value of the development depth of the surrounding rock fracture in each borehole, and after each observation is completed, putting the sleeve 1 into the borehole again to keep the borehole shape, so that the repeated observation is convenient, and finally recovering the sleeve 1.
Peeping holes are drilled at monitoring sites towards the top and bottom plates and the two sides of the roadway, the peeping holes are perpendicular to the surface of the roadway, the number of the holes is designed according to the geometric dimension of the roadway, mutual influence among the observation holes is avoided, and the influence range of the bidirectional isobaric holes is ensured:
And R is the radius of the drilled hole, the minimum distance between two drilled holes with the same radius, which are not affected by each other, is 2R, if the width of the roadway is a, the drilled holes are symmetrically distributed in the roadway, and the distance from the initial drilled hole to the boundary of the roadway is m.
The maximum number of allowed holesIt should be noted that, when the number of two-side holes is calculated, a represents the roadway height.
The average development depth of the fissures has the following relation expression:
Wherein j= { top, bottom, column and wall }, i is the index of the drilling hole at the monitoring point, L j′i is the maximum value of the development depth of the ith drilling crack on the top plate, the bottom plate, the coal pillar upper or the coal wall upper, and X j is the displacement, namely the deformation, of the inner surface of the top plate, the bottom plate, the coal pillar upper or the coal wall upper roadway.
And in the stage of primary mining influence, taking out the casing 1, observing and recording the maximum value of the development depth of the surrounding rock fracture in each drilling hole, and putting the casing 1 into the drilling hole again.
At the stage of primary mining influence, the minimum development depth S Top min＝L Top min+X Top , the maximum development depth S Top max＝L Top max+X Top and the average development depth of roof cracksThe calculation formula of the crack development depth of the bottom plate, the coal pillar upper and the coal wall upper is the same as the calculation formula of the average development depth of the top plate crack. The drawing of the roadway surrounding rock crack development diagram is shown in figure 2.
And in the secondary mining influence stage, taking out the casing, observing and recording the maximum value of the development depth of the surrounding rock fracture in each drilling hole, and putting the casing into the drilling hole again.
The secondary mining influence stage comprises a minimum development depth S ' Top min＝L′ Top min+X′ Top and a maximum development depth S' Top max＝L′ Top max+X′ Top of the roof fracture, and an average development depthThe calculation formula of the crack development depth of the bottom plate, the coal pillar upper and the coal wall upper is the same as the calculation formula of the average development depth of the top plate crack. Drawing the roadway surrounding rock crack development graphic illustration, as shown in fig. 3.
And S3, correcting the damage area in the numerical simulation according to the field actual measurement value, and defining the variation of the crack by using the variation of the area, wherein the variation of the crack is the grouting quantity.
Correcting the fracture damage areas of the three-stage roadway driving, primary mining and secondary mining according to the monitored development depth, wherein the change amount of the area defines the change amount of the fracture, and the change amount of the fracture is grouting amount.
The method for predicting the grouting amount of the surrounding rock fracture is further described below by taking specific engineering as an example.
In the embodiment, the influence of two mining operations of 22205 working face mining and 22204 working face mining on the 22205 working face air return tunnel is taken as an engineering background, the Boolean coal mine 22205 working face air return tunnel is located in a 2-2 coal two-disc area, the length is 4865m, the tunnel is tunneled along a 2-2 coal seam, the average burial depth is 300m, the tunnel section is rectangular, the width is multiplied by the height=5000 mm multiplied by 3500mm, the stoping tunnel adopts double-tunnel arrangement, and the coal pillar width between the double tunnels is 20m. When 22204 working face is mined, a retaining roadway (22205 working face return roadway) is used as a working face auxiliary transportation roadway, when 22205 working face is mined, the retaining roadway (22205 working face return roadway) is used as a return roadway to serve the working face auxiliary transportation roadway, the retaining roadway is influenced by two mining operations, modeling parameters are obtained, and a numerical simulation model is built.
And simulating the whole process of driving, primary mining and secondary mining, and recording the fracture damage range of the surrounding rock of the roadway in each stage of repeated mining, as shown in figures 4-6.
And (3) establishing a two-dimensional plane rectangular coordinate system according to the model, establishing a mathematical formula, and calculating fracture destruction areas of three stages of roadway driving, primary mining and secondary mining.
In the tunneling stage, the fracture destruction area is marked as S 1, points R, S, T, U, V, W, X, Y, Z are marked on a coordinate system respectively,The actual post-drilling correction is shown in fig. 7.
Establishing a mathematical formula to obtain:
In the primary mining stage, the fracture destruction area is marked as S 2, the points A, B, C, D, E, F, G, H are marked on the coordinate system respectively, The actual post-drilling correction is shown in fig. 8.
Establishing a mathematical formula to obtain:
in the secondary mining stage, the fracture destruction area is marked as S 3, the points I, J, K, L, M, N, Q, P are marked on the coordinate system respectively, The actual post-drilling correction is shown in fig. 9.
Establishing a mathematical formula to obtain:
And taking out the sleeve 1 in the monitoring point peeping borehole at each stage of the influence of primary mining and secondary mining, observing and recording the maximum value of the development depth of the surrounding rock fracture in each borehole, and after each observation is completed, putting the sleeve 1 into the borehole again to keep the borehole shape, so that the repeated observation is convenient, and finally recovering the sleeve 1.
Selecting a roadway monitoring station, drilling holes on the top and bottom plates and two sides of a roadway at the monitoring station, and observing the development depth of surrounding rock cracks of the roadway, wherein the method comprises the following steps:
Selecting a surrounding rock crack monitoring site to ensure that the stress in the deep part of the coal wall is close to the stress of the original rock, avoiding a stress reduction area and a stress increase area in front of a working face, wherein the surrounding rock of the top and bottom plates of the selected monitoring site roadway and the coal bodies of the two sides are relatively complete, so that drilling construction and maintenance and later crack observation are facilitated;
Peeping holes are drilled at monitoring sites towards the top and bottom plates and the two sides of the roadway, the peeping holes are perpendicular to the surface of the roadway, the number of the holes is designed according to the geometric dimension of the roadway, mutual influence among the observation holes is avoided, and the influence range of the bidirectional isobaric holes is ensured:
And R is the radius of the drilled hole, the minimum distance between two drilled holes with the same radius, which are not affected by each other, is 2R, if the width of the roadway is a, the drilled holes are symmetrically distributed in the roadway, and the distance from the initial drilled hole to the boundary of the roadway is m.
The maximum number of allowed holesIt should be noted that, when the number of two-side holes is calculated, a represents the roadway height.
After peeping the borehole, the high-strength sleeve 1 is put into the borehole to protect the borehole shape, and judging whether the sleeve 1 needs to be taken out for maintenance according to the sleeve monitoring and dispatching center 5 comprises the following steps:
The high-strength sleeve 1 can reduce or counteract the damage of mining stress to a drilled hole so as to keep the shape of the drilled hole, the top of the front end of the sleeve is spherical, surrounding rock pressure can be evenly distributed, the tail end of the sleeve is sealed by using a threaded cover, a monitoring lead 6 outlet hole is reserved on the sealing cover, and a rubber ring is arranged between the outlet hole and the lead 6 to seal so as to prevent foreign matters from entering the sleeve 1. The length of each section of sleeve 1 is 3.5m, construction is facilitated in a roadway, the number of the connections of each section of sleeve 1 is reduced, the risk that the strength of the connection is low and the connection is easy to be damaged by mining stress is reduced, the sleeve 1 and the sleeve 1 are connected through threads, a strain collector 2 which is a strain gauge is arranged at the connection as shown in fig. 10, the strain gauge is connected with a strain early warning device, the stress deformation of the sleeve 2 is monitored, the sleeve 2 is required to be taken out when the stress deformation exceeds a set threshold value, so that the sleeve is prevented from being clamped, the sleeve 2 is taken out for maintenance, and when a working face is not pushed, a drill bit is used for milling a hole again, and the shape of the drilled hole is maintained.
And checking whether the sleeve pipe has a pipe clamping phenomenon or not according to the distance between the stope face and the monitoring point, checking every other overhaul shift when the distance is far, and checking every overhaul shift when the distance is near.
When the sleeve 1 is put in, the strain collector 2 is arranged in the middle part and at the joint of each section of sleeve 1 and is connected with the external data receiver 3, the data receiver 3 transmits received data to the data processing center 4 and the ground monitoring and dispatching center 5, the stress condition of the sleeve 1 is automatically judged by the preset programs of the data processing center 4 and the ground monitoring and dispatching center 5, and finally, whether the sleeve needs to be taken out or maintained is output, so that workers can check the real-time state of the sleeve 1 at the underground monitoring station and the ground monitoring and dispatching center 5. The system is shown in fig. 10, and the judgment logic diagram is shown in fig. 11.
Taking out the sleeve 1 in the monitoring point peeping borehole at each stage of the influence of one mining, observing and recording the maximum value of the development depth of the surrounding rock crack in each borehole, and after each observation, putting the sleeve 1 into the borehole again to keep the borehole shape so as to facilitate repeated observation, and finally recovering the sleeve 1, wherein the measuring of the deformation of the top and bottom plates of the roadway and the deformation of two sides of the roadway by using the top and bottom plate displacement meter and the roadway surface three-dimensional laser measuring instrument comprises the following steps:
The maximum value of the crack development depth is divided into values on four surfaces of a top plate, a bottom plate, a coal pillar upper and a coal wall upper, the maximum value of the crack development depth of x holes on the top plate is marked as L Top 1、L Top 2、L Top 3…L Top x, the maximum value of the crack development depth of x holes on the bottom plate is marked as L Bottom 1、L Bottom 2、L Bottom 3…L Bottom x, the maximum value of the crack development depth of y holes on the coal pillar upper is marked as L Column 1、L Column 2、L Column 3…L Column y, and the maximum value of the crack development depth of y holes on the coal wall upper is marked as L Wall with a wall body 1、L Wall with a wall body 2、L Wall with a wall body 3…L Wall with a wall body y.
The displacement of the top plate, the bottom plate and the two sides in the one mining influence stage is measured by adopting a top plate displacement meter and a roadway surface three-dimensional laser measuring instrument, so that the average sinking amount of the top plate is X Top , the average bulging amount of the bottom plate is X Bottom , the displacement amount of the coal pillar side is X Column , and the displacement amount of the coal wall side is X Wall with a wall body .
And in the secondary mining influence stage, taking out the casing, observing and recording the maximum value of the development depth of the surrounding rock fracture in each drilling hole, and finally recovering the casing, wherein the method comprises the following steps:
The maximum value of the crack development depth is divided into values on four surfaces of a top plate, a bottom plate, a coal pillar upper and a coal wall upper, the maximum value of the crack development depth of x holes on the top plate is marked as L ' Top 1、L′ Top 2、L′ Top 3…L′ Top x, the maximum value of the crack development depth of x holes on the bottom plate is marked as L' Bottom 1、L′ Bottom 2、L′ Bottom 3…L′ Bottom x, the maximum value of the crack development depth of y holes on the coal pillar upper is marked as L ' Column 1、L′ Column 2、L′ Column 3…L′ Column y, and the maximum value of the crack development depth of y holes on the coal wall upper is marked as L' Wall with a wall body 1、L′ Wall with a wall body 2、L′ Wall with a wall body 3…L′ Wall with a wall body y.
The displacement of the top plate and the bottom plate and the two sides in the secondary mining influence stage is obtained, the average sinking amount of the top plate is X ' Top , the average swelling amount of the bottom plate is X' Bottom , the displacement amount of the coal pillar sides is X ' Column , and the displacement amount of the coal wall sides is X' Wall with a wall body .
And in the stage of primary mining influence, taking out the casing, observing and recording the maximum value of the development depth of the surrounding rock cracks in each drilling hole, calculating the stress and crack distribution range of the surrounding rock, and putting the casing into the drilling hole again.
The primary mining influence stage is to calculate the average development depth according to the maximum value of the development depth of x drill cracks of the top plateL Top i is the maximum value of the development depth of the ith drilling crack of the top plate, and then the calculation formula of the development depth of the crack of the bottom plate, the coal pillar upper and the coal wall upperAverage depth of development of roof fissuresThe calculation formulas are the same, and the drawing of the roadway surrounding rock crack development diagram is shown in figure 2.
Calculating the equivalent damage area of the roof fracture: Then the calculation formula of the equivalent fracture area of the bottom plate, the coal pillar upper and the coal wall upper Equivalent area of failure to roof fractureThe calculation formulas are the same, and the total equivalent destruction area is as follows
And in the secondary mining influence stage, taking out the casing, observing and recording the maximum value of the development depth of the surrounding rock fracture in each drilling hole, and putting the casing into the drilling hole again.
The secondary mining influence stage is to calculate the average development depth according to the development depth of x drilling cracks of the top plateL Top i is the maximum value of the development depth of the ith drilling crack of the top plate, and then the calculation formula of the development depth of the crack of the bottom plate, the coal pillar upper and the coal wall upperAverage depth of development of roof fissuresThe calculation formulas are the same, and the roadway surrounding rock crack development diagram is drawn, as shown in fig. 3.
Calculating the equivalent damage area of the roof fracture: Then the calculation formula of the equivalent fracture area of the bottom plate, the coal pillar upper and the coal wall upper Equivalent area of failure to roof fractureThe calculation formulas are the same, and the total equivalent destruction area is as follows
Calculating a correction coefficient according to the fracture equivalent area obtained by drilling and the fracture destruction area obtained in numerical simulation:
and correcting the fracture damage areas of the primary mining and the secondary mining at two stages in the numerical simulation according to the calculated correction coefficient sigma.
If the correction coefficient sigma is less than or equal to 1, the fracture damage area is not required to be corrected, and if the correction coefficient sigma is more than 1, the fracture damage area is corrected according to the correction coefficient sigma.
And correcting the primary mining influence stage S ' 2＝S2 sigma, and the secondary mining influence stage S ' 3＝S3 sigma, wherein the fracture destruction areas obtained after correction are respectively S ' 2、S'3.
The change amount of fracture destruction area after primary mining is Δs1=s '2-S1, and the change amount of fracture destruction area after secondary mining is Δs2=s' 3-S'2.
The grouting amount of the crack after primary mining is delta 1＝ΔS1*c*103, the grouting amount of the crack after secondary mining is delta 2＝ΔS2*c*103, and c represents the roadway length.
The grouting quantity can be obtained, the structure and the performance of broken surrounding rock can be effectively improved through grouting reinforcement, and the bearing capacity of the surrounding rock of a roadway is integrally improved.
Device embodiment
The invention provides an estimation device for grouting quantity of surrounding rock cracks of a repeated mining roadway, which comprises the following components:
The data acquisition module is used for constructing a repeated mining roadway surrounding rock numerical model, acquiring monitoring data of crack damage ranges of roadway surrounding rock in each stage of repeated mining according to the numerical model, and acquiring monitoring data of crack development depth according to on-site monitoring stations of the roadway.
The area calculation module is used for calculating the area of the fracture theoretical damage range according to the monitoring data of the fracture damage range, calculating the average development depth of the fracture according to the monitoring data of the development depth of the fracture before and after each mining or tunneling, and calculating the total actual equivalent damage area of the space in the roadway according to the average development depth.
The correction module is used for obtaining a correction coefficient according to the ratio of the change amount of the total actual equivalent damage area before and after each mining or tunneling to the change amount of the theoretical damage area, correcting the theoretical damage area according to the correction coefficient to obtain the corrected damage area, and obtaining the change amount of the damage area before and after each mining or tunneling according to the difference between the corrected damage area and the corrected damage area before each mining or tunneling.
And the result acquisition module is used for multiplying the corrected damage area variation by the roadway length to obtain grouting quantity.
Device embodiment
The invention provides a computer device, which comprises a memory and a processor, wherein the memory stores computer execution instructions, and the processor executes the computer execution instructions stored in the memory to realize the method for estimating the grouting amount of the surrounding rock cracks of the repeated mining roadway, which is described above, wherein the method is described in detail in the embodiment of the method and is not repeated here.
Storage medium embodiment
The present invention provides a computer readable storage medium for storing computer executable instructions, which when executed by a processor, are configured to implement a method for estimating grouting amount of surrounding rock cracks of a repeated mining roadway as described above, where the method is described in detail in the method embodiments, and is not described herein.
It will be appreciated by those skilled in the art that embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
It should be noted that the above-described embodiments will enable those skilled in the art to more fully understand the invention, but do not limit it in any way. Therefore, although the present invention has been described in detail with reference to the present specification and examples, it will be understood by those skilled in the art that the present invention may be modified or substituted for others, and all such modifications and improvements which do not depart from the spirit and scope of the present invention are intended to be included in the scope of the present invention. Any reference sign in a claim should not be construed as limiting the claim concerned.
Claims (5)
1. The method for estimating the grouting quantity of the surrounding rock cracks of the repeated mining roadway is characterized by comprising the following steps of:
Constructing a repeated mining roadway surrounding rock numerical model, acquiring monitoring data of crack damage ranges of roadway surrounding rock in each repeated mining stage according to the numerical model, and acquiring monitoring data of crack development depth according to on-site monitoring stations of the roadway;
Calculating average development depth of the cracks according to the monitoring data of the development depth of the cracks before and after each mining or tunneling, and calculating total actual equivalent damage area of the space in the roadway according to the average development depth;
obtaining a correction coefficient according to the ratio of the change amount of the total actual equivalent damage area before and after each mining or tunneling to the change amount of the theoretical damage area, correcting the theoretical damage area according to the correction coefficient to obtain the corrected damage area, and obtaining the change amount of the damage area before and after each mining or tunneling according to the difference between the corrected damage area and the corrected damage area before and after each mining or tunneling;
multiplying the corrected damage area variation by the roadway length to obtain grouting quantity;
The step of calculating the total actual equivalent destruction area of the space in the roadway according to the average development depth comprises the following steps:
The crack development depth comprises crack development depths on four surfaces of a top plate, a bottom plate, a coal pillar upper and a coal wall upper of a roadway;
the relation expression of the average development depth of the fissures is as follows:
Wherein j= { top, bottom, column and wall }, i is the index of a monitoring point drilling, L j′i is the maximum value of the development depth of the ith drilling crack on the top plate, the bottom plate, the coal pillar upper or the coal wall upper, and X j is the displacement, namely the deformation, of the inner surface of the top plate, the bottom plate, the coal pillar upper or the coal wall upper roadway;
the related formula for calculating the total actual equivalent destruction area of the space in the roadway according to the average development depth is as follows:
Wherein, Is the total actual equivalent failure area of the device,Is the actual equivalent damage area of the top plate, the bottom plate, the coal pillar side or the coal wall side, n is the number of the maximum allowed holes, R is the radius of the holes, and R is half of the minimum distance between two holes with the same radius and not affected by each other.
2. The method for estimating a fracture grouting amount of a surrounding rock of a repeated mining roadway according to claim 1, wherein the step of calculating a fracture theoretical damage range area according to the monitoring data of the fracture damage range comprises:
And establishing a two-dimensional plane rectangular coordinate system according to the monitoring data of the fracture damage range, fitting a piecewise function curve of the boundary line according to the point data on the boundary line of the fracture damage range, and carrying out integral summation on the fracture damage range according to the piecewise function curve to obtain the area of the fracture theoretical damage range.
3. The method for estimating a grouting amount of a surrounding rock fracture of a repeated mining roadway according to claim 2, wherein the step of correcting the theoretical damage range area according to the correction coefficient to obtain the corrected damage range area comprises the following steps:
If the correction coefficient sigma is less than or equal to1, the fracture theoretical damage range area does not need to be corrected, and if the correction coefficient sigma is more than 1, the fracture theoretical damage range area is corrected according to the correction coefficient sigma, and a specific correction formula is as follows:
S'2＝S2*σ
Wherein S 2 is the fracture theoretical damage range area of the current mining influence stage, and S' 2 is the fracture damage range area corrected by the current mining influence stage.
4. An estimation device of repeated mining roadway surrounding rock crack grouting quantity is characterized by comprising:
The data acquisition module is used for constructing a repeated mining roadway surrounding rock numerical model, acquiring monitoring data of crack damage ranges of roadway surrounding rock in each stage of repeated mining according to the numerical model, and acquiring monitoring data of crack development depth according to on-site monitoring stations of the roadway;
the area calculation module is used for calculating the area of the fracture theoretical damage range according to the monitoring data of the fracture damage range, calculating the average development depth of the fracture according to the monitoring data of the development depth of the fracture before and after each mining or tunneling, and calculating the total actual equivalent damage area of the space in the roadway according to the average development depth;
The correction module is used for obtaining a correction coefficient according to the ratio of the change amount of the total actual equivalent damage area before and after each mining or tunneling to the change amount of the theoretical damage area, correcting the theoretical damage area according to the correction coefficient to obtain a corrected damage area, and obtaining the change amount of the damage area after each mining or tunneling according to the difference between the corrected damage area and the corrected damage area before each mining or tunneling;
the result acquisition module is used for multiplying the corrected damage area variation by the roadway length to obtain grouting quantity;
The step of calculating the total actual equivalent destruction area of the space in the roadway according to the average development depth comprises the following steps:
The crack development depth comprises crack development depths on four surfaces of a top plate, a bottom plate, a coal pillar upper and a coal wall upper of a roadway;
the relation expression of the average development depth of the fissures is as follows:
Wherein j= { top, bottom, column and wall }, i is the index of a monitoring point drilling, L j′i is the maximum value of the development depth of the ith drilling crack on the top plate, the bottom plate, the coal pillar upper or the coal wall upper, and X j is the displacement, namely the deformation, of the inner surface of the top plate, the bottom plate, the coal pillar upper or the coal wall upper roadway;
the related formula for calculating the total actual equivalent destruction area of the space in the roadway according to the average development depth is as follows:
Wherein, Is the total actual equivalent failure area of the device,Is the actual equivalent damage area of the top plate, the bottom plate, the coal pillar side or the coal wall side, n is the number of the maximum allowed holes, R is the radius of the holes, and R is half of the minimum distance between two holes with the same radius and not affected by each other.
5. A computer device comprising a memory and a processor, wherein the memory has stored therein computer-executable instructions, and the processor executes the computer-executable instructions stored in the memory to implement a method of estimating a repeated mining roadway surrounding rock fracture grouting volume as claimed in any one of claims 1-3.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410761924.7A CN118607238B (en) | 2024-06-13 | 2024-06-13 | Method, device and equipment for estimating grouting quantity of surrounding rock cracks of repeated mining roadway |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410761924.7A CN118607238B (en) | 2024-06-13 | 2024-06-13 | Method, device and equipment for estimating grouting quantity of surrounding rock cracks of repeated mining roadway |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN118607238A CN118607238A (en) | 2024-09-06 |
| CN118607238B true CN118607238B (en) | 2025-01-03 |
Family
ID=92560653
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410761924.7A Active CN118607238B (en) | 2024-06-13 | 2024-06-13 | Method, device and equipment for estimating grouting quantity of surrounding rock cracks of repeated mining roadway |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN118607238B (en) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106246162A (en) * | 2016-09-21 | 2016-12-21 | 山东科技大学 | Floor undulation is across borescopic imaging device and slip casting effect monitoring method |
| CN111271129A (en) * | 2020-01-31 | 2020-06-12 | 天地科技股份有限公司 | Method for acquiring deformation and fracture expansion rule of stoping roadway surrounding rock |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112127909B (en) * | 2020-09-08 | 2021-12-07 | 河海大学 | Accurate grouting repairing and reinforcing method for broken surrounding rock of tunnel |
| CN113586091B (en) * | 2021-08-16 | 2023-10-17 | 许昌学院 | A grouting method to prevent and control water damage considering the impact of mining cracks |
| CN114233394B (en) * | 2021-11-26 | 2023-10-31 | 安徽理工大学 | Surrounding rock monitoring and supporting method for stoping roadway |
| CN114635697B (en) * | 2022-02-25 | 2023-01-24 | 中国矿业大学 | Pre-grouting treatment method for water damage of shaft well construction |
| CN116503511B (en) * | 2023-05-06 | 2023-09-19 | 中国矿业大学(北京) | Immediate mapping imaging method for morphology of surrounding rock damage area of tunneling roadway |
-
2024
- 2024-06-13 CN CN202410761924.7A patent/CN118607238B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106246162A (en) * | 2016-09-21 | 2016-12-21 | 山东科技大学 | Floor undulation is across borescopic imaging device and slip casting effect monitoring method |
| CN111271129A (en) * | 2020-01-31 | 2020-06-12 | 天地科技股份有限公司 | Method for acquiring deformation and fracture expansion rule of stoping roadway surrounding rock |
Also Published As
| Publication number | Publication date |
|---|---|
| CN118607238A (en) | 2024-09-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Hoek | Support for very weak rock associated with faults and shear zones | |
| Lin et al. | Experimental study on failure behaviour of deep tunnels under high in-situ stresses | |
| CN104763432B (en) | A kind of method that high stress tunnel country rock release controls large deformation | |
| CN103293560B (en) | A kind of method of testing adopting triple stress court | |
| CN104653226A (en) | Stress-gradient-based method for dividing coal impact ground pressure danger area | |
| CN110374636A (en) | A kind of construction method of tunnel rock burst security protection | |
| Broch et al. | The design of unlined hydropower tunnels and shafts: 100 years of Norwegian experience | |
| CN104035142A (en) | Early warning method for sudden gushing water danger of undersea tunnel | |
| Duddeck | Application of numerical analyses for tunnelling | |
| CN103216265B (en) | A kind of high stress tunnel surrouding rock deformation control method | |
| CN118607238B (en) | Method, device and equipment for estimating grouting quantity of surrounding rock cracks of repeated mining roadway | |
| CN111709065B (en) | Method for determining hole distance of ultra-front deep hole grouting drilling under complex surrounding environment | |
| CN106290016B (en) | A kind of underground air impact disaster experimental rig | |
| Jiang et al. | Automatic Modeling and Optimization Design Platform for Highway Tunnels with Application to Stress and Displacement Prediction of Surrounding Rocks during Tunnel Excavation | |
| CN117034794A (en) | Method, device, equipment and medium for simulating ortho-position advanced pre-grouting drilling | |
| CN114329935B (en) | Testing method and device for three-dimensional mining stress of surrounding rock in coal mine tunnel | |
| CN112685926A (en) | Construction period discrimination method for shallow coal seam air leakage prevention project | |
| CN119437543B (en) | Test platform and comprehensive evaluation method for water and soil load of shield tunneling machine | |
| Zhang et al. | Study on the evaluation of surrounding rock stability and limit displacement of high-pressure stratum tunnel | |
| CN114856587B (en) | Optimized construction method of composite stratum slurry shield based on BIM technology | |
| CN119914356B (en) | Accurate coal and gas outburst prediction method based on gas geological analysis | |
| CN120175353A (en) | Mechanized rapid construction method for underground chamber | |
| Huang et al. | OPTIMIZATION OF MECHANIZED RAPID CONSTRUCTION STEP DISTANCE FOR HARD ROCK BODY OF MA BAISHAN TUNNEL | |
| Niklasson et al. | The Sture LPG Storage Project | |
| Zhang | General report for modelling and testing of tunnels and deep excavations in soft ground |
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 |