CN115048821B - Prediction method of lagging rock burst - Google Patents

Prediction method of lagging rock burst Download PDF

Info

Publication number
CN115048821B
CN115048821B CN202210973548.9A CN202210973548A CN115048821B CN 115048821 B CN115048821 B CN 115048821B CN 202210973548 A CN202210973548 A CN 202210973548A CN 115048821 B CN115048821 B CN 115048821B
Authority
CN
China
Prior art keywords
rock burst
cracks
value
acoustic emission
crack
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202210973548.9A
Other languages
Chinese (zh)
Other versions
CN115048821A (en
Inventor
刘冬桥
孙杰
杨园园
张梓谦
李娜
胡义
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Mining and Technology Beijing CUMTB
Original Assignee
China University of Mining and Technology Beijing CUMTB
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China University of Mining and Technology Beijing CUMTB filed Critical China University of Mining and Technology Beijing CUMTB
Priority to CN202210973548.9A priority Critical patent/CN115048821B/en
Publication of CN115048821A publication Critical patent/CN115048821A/en
Application granted granted Critical
Publication of CN115048821B publication Critical patent/CN115048821B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/14Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object using acoustic emission techniques
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/44Processing the detected response signal, e.g. electronic circuits specially adapted therefor
    • G01N29/4409Processing the detected response signal, e.g. electronic circuits specially adapted therefor by comparison
    • G01N29/4418Processing the detected response signal, e.g. electronic circuits specially adapted therefor by comparison with a model, e.g. best-fit, regression analysis
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00Indexing codes associated with group G01N29/00
    • G01N2291/02Indexing codes associated with the analysed material
    • G01N2291/023Solids
    • G01N2291/0232Glass, ceramics, concrete or stone
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Pathology (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Immunology (AREA)
  • Health & Medical Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Signal Processing (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

The invention provides a prediction method of lagging rock burst, and belongs to the technical field of rock burst prediction. The invention discloses a prediction method of lagging rock burst, which comprises the following steps: carrying out an acoustic emission experiment on the rock sample to construct a classification model of the crack; performing an indoor hysteretic rockburst experiment on the rock sample to obtain an RA value and an AF value of acoustic emission of hysteretic rockburst; substituting the RA value and the AF value of the acoustic emission of the lag rock burst into a classification model to obtain the crack type distribution of the whole process of the lag rock burst; calculating the accumulated proportion of the cracks in the whole process of the lag rock burst to obtain the corresponding relation between the accumulated proportion of the cracks and time; and monitoring the corresponding relation in real time to judge whether the lagging rock burst occurs. The invention provides a prediction method of a lag rock burst, which is used for realizing accurate prediction of the lag rock burst.

Description

Prediction method of lagging rock burst
Technical Field
The invention relates to the technical field of rockburst prediction, in particular to a prediction method of lagging rockburst.
Background
The lag rock burst is a common rock burst type, is usually generated in a deep space after excavation is finished, and is different from hours to months, the hazard of the lag rock burst is strong, the generation time is unknown, and the generation mechanism and the influence factors are complex.
In the prior art, the grade of the rock burst which is possible to occur is mostly predicted based on lithology, engineering geological conditions and the stress state of surrounding rocks, the method is mostly based on judgment of experience, and has great subjective factors, for example, in on-site engineering application, the micro-seismic monitoring is often adopted to predict the rock burst, but a gap exists between the performance and the actual requirement, the gap can be influenced by a plurality of factors such as exploitation factors, geological factors, station measuring factors, monitoring factors and algorithm factors, and therefore the precision of the micro-seismic monitoring has errors; when rock burst occurs, the micro-seismic energy or event number increases sharply, but rock burst cannot occur when the micro-seismic energy or event number increases sharply, so the micro-seismic monitoring method is only a necessary condition but not a sufficient condition for rock burst occurrence. In addition, due to the complexity of the actual engineering, the stress state of the surrounding rock is often in an unknown state, so the rock burst cannot be predicted accurately in real time.
Therefore, how to accurately predict the hysteresis rock burst is a technical problem which needs to be solved urgently.
Disclosure of Invention
The invention aims to provide a method for predicting a lagging rock burst, which is used for improving the accuracy of predicting the lagging rock burst.
The above object of the present invention can be achieved by the following technical solutions:
the invention provides a prediction method of lagging rock burst, which comprises the following steps:
carrying out an acoustic emission experiment on the rock sample to construct a classification model of the cracks;
performing an indoor hysteretic rockburst experiment on the rock sample to obtain an RA value and an AF value of acoustic emission of hysteretic rockburst;
substituting the RA value and the AF value of the acoustic emission of the lag rock burst into a classification model to mark out the crack type to which the acoustic emission signal generated at each time point in the whole process of the lag rock burst belongs;
calculating the accumulated proportion of the cracks in the whole process of the lag rock burst to obtain the corresponding relation between the accumulated proportion of the cracks and time;
and monitoring the corresponding relation in real time to judge whether the lagging rock burst occurs.
In an embodiment of the invention, the crack is a shear crack and/or a tension crack.
In an embodiment of the invention, a direct shear acoustic emission test is performed on a rock sample to obtain an RA value and an AF value of the direct shear acoustic emission test;
drawing a first density cloud picture according to the obtained RA value and AF value of the direct shear acoustic emission test;
defining an area with a slope of 0-40 in the first density cloud picture as a shear crack aggregation area S1, wherein a label of a data point in the shear crack aggregation area S1 is {1};
and training data points in the shear crack aggregation zone S1 through a support vector machine to construct a classification model of the shear cracks.
In an embodiment of the invention, a Brazilian acoustic emission test is performed on a rock sample to obtain an RA value and an AF value of the Brazilian acoustic emission test;
drawing a second density cloud picture according to the obtained RA value and AF value of the Brazilian splitting acoustic emission test;
defining the region with the slope larger than 400 in the second density cloud map as a tension crack aggregation region S2, wherein the label of the data point is {2};
and training data points in the tension crack aggregation area S2 through a support vector machine to construct a classification model of the tension cracks.
In an embodiment of the present invention, the classification model isfIt is represented as:
Figure 140602DEST_PATH_IMAGE001
(ii) a Wherein the content of the first and second substances,
an input set I { (RA, AF) | RA ∈ S1 { (RA, AF) | S2, AF ∈ S1 { (RA, AF) | S2};
output set O = { "1", "2" }.
In the embodiment of the invention, in the RA-AF density cloud picture corresponding to the classification model, the boundary line of the tension crack and the shear crack
Figure 746770DEST_PATH_IMAGE002
The acoustic emission signal in the area below the boundary is a shear crack, and the acoustic emission signal in the area above the boundary is a tension crack.
In the embodiment of the invention, the loading mode of the hysteresis rock burst test is three-way six-sided loading-single-sided sudden unloading-load retention-axial loading.
In an embodiment of the present invention, the substrate is,
the initial stress level of the loading was: sigma v =26MPa,σ h =33MPa,σ H =21MPa, the whole loading process is divided into five stages; wherein the content of the first and second substances,
the first stage is loading three-way stress to initial ground stress, the second stage is loading holding stage after loading to initial ground stress state, and the third stage is sigma h Unloading to 0MPa, and the fourth stage is sigma h After unloading to 0MPa, the fifth stage is loading sigma v And finally, a rock burst stage.
In the embodiment of the present invention, the correspondence relationship is a crack cumulative percentage map obtained based on the cumulative percentage of cracks and the time corresponding to the cumulative percentage of cracks.
In the embodiment of the invention, the crack cumulative proportion graph is a line graph, the abscissa of the line graph is time, and the ordinate of the line graph is the cumulative proportion;
and monitoring the trend of the accumulated proportion in the line graph changing along with time in real time to judge whether the lag rock burst occurs.
The invention has the characteristics and advantages that:
the prediction method of the lag rock burst disclosed by the invention starts from the angle of crack evolution, quantitatively describes the change of the accumulated occupation ratios of different types of cracks in the whole process of the lag rock burst, can accurately realize the prediction of the lag rock burst through the change trend of the accumulated occupation ratios of the cracks, and does not predict the possible grade of the rock burst based on lithology, engineering geological conditions and the stress state of surrounding rocks in the prior art, namely, the prediction method is not determined by experience and is not easily influenced by factors such as mining factors, geological factors, station measuring factors, monitoring factors and algorithm factors, so that the prediction method has the advantages of small error and convenience in application, and can accurately realize the prediction of the lag rock burst.
Drawings
The above and other features and advantages of the present invention will become more apparent by describing in detail exemplary embodiments thereof with reference to the attached drawings.
FIG. 1 is a flow chart of a method of predicting a lag rock burst in accordance with the present invention;
FIG. 2 is another flow chart of a method of predicting a lag rock burst in accordance with the present invention;
FIG. 3 is a stress path diagram of a lag rock burst in accordance with the present invention;
FIG. 4 shows the results of the classification of cracks of a lag rock burst according to the present invention;
FIG. 5 is a cumulative event number distribution characteristic of different types of cracks in the present invention;
FIG. 6 is a cumulative percentage of cracks curve according to the present invention;
fig. 7 is an enlarged view of the cumulative percentage of tensile cracks curve in stage v of fig. 6.
Detailed Description
Example embodiments will now be described more fully with reference to the accompanying drawings. Example embodiments may, however, be embodied in many different forms and should not be construed as limited to the embodiments set forth herein; rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the concept of example embodiments to those skilled in the art. The same reference numerals in the drawings denote the same or similar structures, and thus their detailed description will be omitted.
The terms "a," "an," "the," "said" are used to indicate the presence of one or more elements/components/etc.; the terms "comprising" and "having" are intended to be inclusive and mean that there may be additional elements/components/etc. other than the listed elements/components/etc.
As shown in fig. 1 to 7, the present invention provides a method for predicting a lag rock burst, including: carrying out an acoustic emission experiment on the rock sample to construct a classification model of the crack, wherein the crack can be a shear crack and/or a tension crack; performing an indoor lag rockburst experiment on the rock sample to obtain an RA value and an AF value of acoustic emission of lag rockburst; substituting the RA value and the AF value of the acoustic emission of the lag rock burst into a classification model (or a crack signal identification library) to mark out the crack type to which the acoustic emission signal generated at each time point in the whole process of the lag rock burst belongs; calculating the accumulated proportion of the cracks in the whole process of the lag rock burst to obtain the corresponding relation between the accumulated proportion of the cracks and time; and monitoring the corresponding relation in real time, for example, observing the change trend of the cumulative percentage of the cracks to judge whether the lagging rock burst occurs.
The prediction method of the lag rock burst can quantitatively describe the change of the accumulated occupation ratios of different types of cracks in the whole process of the lag rock burst from the perspective of crack evolution, can accurately predict the lag rock burst through the change trend of the accumulated occupation ratios of the cracks, and does not predict the possible grade of the rock burst based on lithology, engineering geological conditions and the stress state of surrounding rocks in the prior art, namely, the prediction method is not determined by experience and is not easily influenced by factors such as mining factors, geological factors, station measuring factors, monitoring factors, algorithm factors and the like, so that the prediction method has the advantages of small error and convenience in application, and can accurately predict the lag rock burst.
In one embodiment of the present invention, performing an acoustic emission experiment on a rock sample to construct a classification model of cracks includes: performing a direct shear acoustic emission test on the rock sample to obtain an RA value and an AF value of the direct shear acoustic emission test; drawing a first density cloud picture according to the obtained RA value and AF value of the direct shear acoustic emission test; defining an area with a slope ranging from 0 to 40 in the first density cloud chart as a shear crack accumulation zone S1, wherein a label of a data point in the shear crack accumulation zone S1 is {1}; most acoustic emission signals of the direct shear acoustic emission test are shear signals so as to construct a classification model of the shear cracks; the support vector machine can train data points in the shear crack gathering area S1 to construct a classification model of the shear cracks, has good advantages in dividing nonlinear data as a machine learning algorithm with supervised learning, and can accurately divide the shear cracks.
Further, carrying out Brazilian splitting acoustic emission test on the rock sample to obtain an RA value and an AF value of the Brazilian splitting acoustic emission test; drawing a second density cloud picture according to the obtained RA value and AF value of the Brazilian splitting acoustic emission test, wherein most of signals obtained from the Brazilian splitting test are tension crack signals so as to construct a classification model of tension cracks; defining a region with a slope larger than 400 in the second density cloud image as a tension crack aggregation region S2, wherein the label of a data point is {2}; and training data points in the tension crack aggregation area S2 through a support vector machine to construct a classification model of the tension cracks.
In the present invention, the classification model isfIt is expressed as:
Figure 610821DEST_PATH_IMAGE003
(ii) a Wherein, the input set I = { (RA, AF) | RA ∈ S1 { (RA, AF) | S2}, AF ∈ S1 { (RA, AF) }; output set O = { "1", "2" }.
Specifically, in the RA-AF density cloud chart corresponding to the classification model, the boundary line of the tension crack and the shear crack
Figure 18669DEST_PATH_IMAGE004
The acoustic emission signal in the area below the boundary is a shear crack, and the acoustic emission signal in the area above the boundary is a tension crack.
Specifically, as shown in fig. 3, the indoor hysteresis rockburst experiment performed on the rock sample includes: the red sandstone can be selected for carrying out a lag rock burst test, and the loading mode of the lag rock burst test is three-dimensional six-surface loading, single-surface sudden unloading, load retention and axial loading.
More specifically, the initial stress level of loading may be: sigma v =26MPa,σ h =33MPa,σ H =21MPa, the whole loading process is divided into five stages; wherein the content of the first and second substances,
the first stage is loading three-way stress to initial ground stress, the second stage is loading holding stage after loading to initial ground stress state, and the third stage is sigma h And unloading to 0MPa, wherein the fourth stage is a load-holding stage after single-side unloading, and the fifth stage is a rock burst stage.
In specific application, an acoustic emission test can be carried out according to the hysteresis rockburst test to obtain acoustic emission data, an RA value and an AF value of an acoustic emission parameter are calculated, and the calculated RA value and AF value of the acoustic emission parameter are substituted into the classification modelfThe crack type distribution of the whole process of the lagging rock burst can be obtained.
Further, as shown in fig. 4 and 5, the crack evolution characteristic of the entire hysteresis rock burst process can be obtained by calculating the cumulative crack occupation ratio of the entire hysteresis rock burst process according to the following formulas (1) and (2) based on the crack species distribution (or the crack distribution characteristic).
Figure 575552DEST_PATH_IMAGE005
(1)
Figure 311427DEST_PATH_IMAGE006
(2)
Wherein p is t And p s Respectively the cumulative percentage of the tensile cracks and the cumulative percentage of the shear cracks from the beginning of loading to the time t; m is a group of t The cumulative number of the tension cracks from the beginning of loading to the time t; m is a group of s The accumulated number of the shear cracks from the beginning of loading to the moment t is shown; m is the total number of acoustic emission events from the start of loading to time t.
Monitoring the corresponding relation in real time to judge whether the hysteresis rock burst occurs comprises the following steps: the corresponding relation is a crack accumulative ratio graph obtained based on the accumulative ratio of the cracks and the time corresponding to the accumulative ratio of the cracks, the relation of the accumulative ratio of the cracks changing along with the time can be observed according to the crack accumulative ratio graph, and the relation of the accumulative ratio of the cracks changing along with the time is monitored in real time to judge whether the lag rock burst occurs.
In a possible embodiment, the crack cumulative proportion map is a line graph, the abscissa of the line graph is time, and the ordinate of the line graph is the cumulative proportion; and monitoring the trend of the accumulated proportion in the line graph along with the change of time in real time to judge whether the hysteresis rock burst occurs, and in specific application, observing the trend (such as slope) of the line graph along with the change of time to judge whether the hysteresis rock burst occurs.
In the present invention, as can be seen from fig. 6, the cumulative percentage of cracks curve can show different variation trends in different stages, taking the cumulative percentage of tensile cracks as an example, stage i is a stage of loading to initial stress, and the cumulative percentage of tensile cracks is fasterLifting; with the application of stress, the red sandstone gradually enters an elastic stage, and the cumulative percentage of tensile cracks is slowed down and kept stable; the stage II is a load-holding stage, the stress state of the red sandstone is not changed, so that an acoustic emission signal is not generated, and a crack accumulated proportion curve is kept unchanged; stage III is σ h When the crack is quickly unloaded to the 0MPa stage, the accumulated crack occupation ratio curve is kept stable firstly, and the crack occupation ratio curve is slightly reduced along with the falling of the cushion block and is kept stable along with the stretching of the empty surface; the stage IV is a load-holding stage, the stress state of the red sandstone is not changed, so that an acoustic emission signal is not generated, and a crack accumulated proportion curve is kept unchanged; stage V is loading σ v And at the rock burst stage, the cumulative occupation ratio curve of the tension cracks shows three processes of keeping stability, slowly descending and violently and continuously descending, and reaches the lowest point when the rock burst occurs.
The key of rock burst prediction is to facilitate the distinction of the catastrophe points, as shown in fig. 6 and 7, the cumulative occupancy curve of the tension cracks and the shear crack curve have distinct changing trends in different stages, taking the cumulative occupancy curve of the tension cracks as an example, in the stages i to iv, the cumulative occupancy curve of the tension cracks increases with the application of stress or keeps unchanged with the increase of stress, and in the stage v, the cumulative occupancy curve of the tension cracks decreases with the application of stress, which is different from the first four stages and is in strong contrast with the phenomenon in the first four stages, so that once the cumulative occupancy curve of the tension cracks is severely and continuously decreased, the occurrence of the lagging rock burst can be represented.
In embodiments of the present invention, the term "plurality" means two or more unless explicitly defined otherwise. The terms "mounted," "connected," "fixed," and the like are to be construed broadly and include, for example, fixed connections, removable connections, or integral connections. Specific meanings of the above terms in the embodiments of the present invention can be understood by those of ordinary skill in the art according to specific situations.
In the description of the embodiments of the present invention, it should be understood that the terms "upper", "lower", and the like indicate orientations or positional relationships based on those shown in the drawings, and are only for convenience in describing the embodiments of the present invention and simplifying the description, but do not indicate or imply that the referred devices or units must have a specific direction, be configured in a specific orientation, and operate, and thus, should not be construed as limiting the embodiments of the present invention.
In the description herein, the appearances of the phrase "one embodiment," "a preferred embodiment," or the like, are intended to mean that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the embodiments of the invention. In this specification, the schematic representations of the terms used above do not necessarily refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.
The above description is only a preferred embodiment of the present invention, and is not intended to limit the present invention, and various modifications and changes may be made to the present embodiment by those skilled in the art. Any modification, equivalent replacement, improvement and the like made within the spirit and principle of the embodiments of the present invention should be included in the protection scope of the embodiments of the present invention.

Claims (4)

1. A method for predicting a lag rock burst, comprising:
carrying out an acoustic emission experiment on the rock sample to construct a classification model of the cracks;
performing an indoor hysteretic rockburst experiment on the rock sample to obtain an RA value and an AF value of acoustic emission of hysteretic rockburst;
substituting the RA value and the AF value of the acoustic emission of the lag rock burst into the classification model to mark out the crack type to which the acoustic emission signal generated at each time point in the whole process of the lag rock burst belongs;
calculating the accumulated proportion of the cracks in the whole process of the lag rock burst so as to obtain the corresponding relation between the accumulated proportion of the cracks and time;
monitoring the corresponding relation in real time to judge whether hysteresis rock burst occurs or not;
the cracks are shear cracks and/or tension cracks;
carrying out Brazilian splitting acoustic emission test on the rock sample to obtain an RA value and an AF value of the Brazilian splitting acoustic emission test;
drawing a first density cloud picture according to the obtained RA value and AF value of the Brazilian acoustic emission test;
defining an area with a slope of 0-40 in the first density cloud picture as a tension crack aggregation area S1, wherein a label of a data point in the tension crack aggregation area S1 is {1};
training data points in the tension crack gathering area S1 through a support vector machine to construct a classification model of the tension cracks;
performing a direct shear acoustic emission test on the rock sample to obtain an RA value and an AF value of the direct shear acoustic emission test;
drawing a second density cloud picture according to the obtained RA value and AF value of the direct shear acoustic emission test;
defining a region with a slope larger than 400 in the second density cloud image as a shear crack aggregation zone S2, wherein the label of a data point is {2};
training data points in the shear crack aggregation zone S2 through the support vector machine to construct a classification model of the shear crack;
in the RA-AF density cloud picture corresponding to the classification model, the boundary of the tension crack and the shear crack
Figure 113400DEST_PATH_IMAGE002
Acoustic emission signals in the area above the boundary are shear cracks, and acoustic emission signals in the area below the boundary are tension cracks;
the loading mode of the hysteresis rock burst test is three-dimensional six-surface loading, single-surface sudden unloading, load retention and axial loading;
the initial ground stress level of loading was: sigma v =26MPa,σ h =33MPa,σ H =21MPa, the whole loading process is divided into five stages; wherein, the first and the second end of the pipe are connected with each other,
first orderThe section is a stage of loading three-dimensional stress to initial ground stress, the second stage is a load-holding stage after loading to the initial ground stress state, and the third stage is sigma h Unloading to 0MPa, and the fourth stage is sigma h After unloading to 0MPa, the fifth stage is loading sigma v And finally, a rock burst stage.
2. The method of predicting a lag rock burst of claim 1,
the classification model is f, which is expressed as:
Figure 499382DEST_PATH_IMAGE004
(ii) a Wherein, the first and the second end of the pipe are connected with each other,
the input set I = { (RA, AF) | RA ∈ S1 { (RA, AF) | S2}, AF ∈ S1 { (RA, AF) | S2};
output set O = { "1", "2" }.
3. The method of predicting a lag rock burst of claim 1,
the correspondence relationship is a crack cumulative proportion map obtained based on the cumulative proportion of cracks and the time corresponding to the cumulative proportion of cracks.
4. The method of predicting a lag rock burst of claim 3,
the crack accumulated proportion graph is a line graph, the abscissa of the line graph is time, and the ordinate of the line graph is the accumulated proportion;
and monitoring the trend of the accumulated proportion in the line graph changing along with time in real time to judge whether the lagging rock burst occurs.
CN202210973548.9A 2022-08-15 2022-08-15 Prediction method of lagging rock burst Active CN115048821B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210973548.9A CN115048821B (en) 2022-08-15 2022-08-15 Prediction method of lagging rock burst

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210973548.9A CN115048821B (en) 2022-08-15 2022-08-15 Prediction method of lagging rock burst

Publications (2)

Publication Number Publication Date
CN115048821A CN115048821A (en) 2022-09-13
CN115048821B true CN115048821B (en) 2022-11-18

Family

ID=83168053

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210973548.9A Active CN115048821B (en) 2022-08-15 2022-08-15 Prediction method of lagging rock burst

Country Status (1)

Country Link
CN (1) CN115048821B (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111272883A (en) * 2020-03-04 2020-06-12 成都理工大学 Rock fracture mode intelligent detection and identification method based on acoustic emission model
CN111812211A (en) * 2020-07-09 2020-10-23 武汉理工大学 RA-AF-E rock material brittle fracture crack classification method based on acoustic emission parameters
CN113063857A (en) * 2020-06-15 2021-07-02 中国科学院武汉岩土力学研究所 Acoustic emission identification method for rock structural surface tension-shear failure in direct shear test
CN113702509A (en) * 2021-10-29 2021-11-26 中国矿业大学(北京) Method for stretching and shearing cracks in divided rock burst process

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100596328C (en) * 2007-05-16 2010-03-31 中国矿业大学(北京) Test method for deep rock explosion process model
CN109827846A (en) * 2019-02-02 2019-05-31 中南大学 It is a kind of based on add unloading response lag than index Rock burst proneness Grade Judgment
CN111044367B (en) * 2019-12-31 2020-12-18 成都理工大学 Rock crack subcritical propagation rate experimental test method based on triaxial stress-strain curve

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111272883A (en) * 2020-03-04 2020-06-12 成都理工大学 Rock fracture mode intelligent detection and identification method based on acoustic emission model
CN113063857A (en) * 2020-06-15 2021-07-02 中国科学院武汉岩土力学研究所 Acoustic emission identification method for rock structural surface tension-shear failure in direct shear test
CN111812211A (en) * 2020-07-09 2020-10-23 武汉理工大学 RA-AF-E rock material brittle fracture crack classification method based on acoustic emission parameters
CN113702509A (en) * 2021-10-29 2021-11-26 中国矿业大学(北京) Method for stretching and shearing cracks in divided rock burst process

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于蠕变试验的时滞型岩爆孕育机理研究;李栋等;《地下空间与工程学报》;20200831;第71-78页 *

Also Published As

Publication number Publication date
CN115048821A (en) 2022-09-13

Similar Documents

Publication Publication Date Title
CN110210144B (en) Optimization design method for promoting uniform expansion of horizontal well fracturing fracture by temporary plugging agent
CN104406623B (en) Slope dynamic stability coefficient determination method based on level of ground water and displacement monitoring
US8694292B2 (en) Method and system for estimating building performance
CN110532694B (en) Arch dam safety evaluation method based on geomechanical model comprehensive method test
CN110909802A (en) Improved PSO (particle swarm optimization) based fault classification method for optimizing PNN (portable network) smoothing factor
CN109635325A (en) Reservoir landslide stability prediction method based on compound hydrodynamic force and displacement monitoring
Nguyen et al. Forecasting Damage Mechanics By Deep Learning.
CN114169254A (en) Abnormal energy consumption diagnosis method and system based on short-term building energy consumption prediction model
CN113376005B (en) Method for predicting final rheological strain of rock-fill material based on long-term strength level of particles
Olea A practical primer on geostatistics
CN104570109B (en) A kind of method of reservoir hydrocarbons prediction
CN113378967A (en) Structural health monitoring multivariate data anomaly diagnosis method based on convolutional neural network and transfer learning
Jia et al. Bridge seismic damage assessment model applying artificial neural networks and the random forest algorithm
CN115048821B (en) Prediction method of lagging rock burst
Wu et al. Assessment of environmental and nondestructive earthquake effects on modal parameters of an office building based on long-term vibration measurements
Jagerman et al. Vessel arrival process and queueing in marine ports handling bulk materials
Hendrick et al. Automated prediction of wet-snow avalanche activity in the Swiss Alps
Bagherkhani et al. Structure-MR damper reliability analysis using weighted uniform simulation method
CN113836625A (en) Power transmission tower anti-seismic performance grade division method and system based on power spectrum method
CN103413022B (en) Oil recovery factor calibration auxiliary device and householder method thereof
Oberguggenberger et al. From probability to fuzzy sets: the struggle for meaning in geotechnical risk assessment
CN111753399A (en) Method for predicting filling slurry ring pipe pressure drop by machine learning
CN110044717A (en) Determine the method that Loading ami deloading response ratio plays height in Rock Under Uniaxial Compression compression classification plus unloading test
CN116150667A (en) XGBoost-based flight delay reason pre-judging method
CN111479235B (en) Formalization method for detecting wireless sensor network node

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