CN116539413A - Rock characteristic stress determining method based on acoustic emission counting nonlinear evolution - Google Patents
Rock characteristic stress determining method based on acoustic emission counting nonlinear evolution Download PDFInfo
- Publication number
- CN116539413A CN116539413A CN202310760691.4A CN202310760691A CN116539413A CN 116539413 A CN116539413 A CN 116539413A CN 202310760691 A CN202310760691 A CN 202310760691A CN 116539413 A CN116539413 A CN 116539413A
- Authority
- CN
- China
- Prior art keywords
- curve
- acoustic emission
- strain
- stress
- sub
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000011435 rock Substances 0.000 title claims abstract description 145
- 238000000034 method Methods 0.000 title claims abstract description 66
- 230000001186 cumulative effect Effects 0.000 claims abstract description 149
- 230000008859 change Effects 0.000 claims abstract description 129
- 238000012417 linear regression Methods 0.000 claims abstract description 28
- 238000007906 compression Methods 0.000 claims abstract description 22
- 230000006835 compression Effects 0.000 claims abstract description 18
- 238000002474 experimental method Methods 0.000 claims abstract description 15
- 238000009825 accumulation Methods 0.000 claims description 53
- 238000005336 cracking Methods 0.000 claims description 34
- 238000004590 computer program Methods 0.000 claims description 11
- 238000004422 calculation algorithm Methods 0.000 claims description 7
- 230000001174 ascending effect Effects 0.000 claims description 5
- 230000006870 function Effects 0.000 claims 1
- 230000008569 process Effects 0.000 description 11
- 230000000977 initiatory effect Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000009471 action Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000003247 decreasing effect Effects 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 230000000630 rising effect Effects 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 206010010214 Compression fracture Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/02—Details
- G01N3/06—Special adaptations of indicating or recording means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L1/00—Measuring force or stress, in general
- G01L1/25—Measuring force or stress, in general using wave or particle radiation, e.g. X-rays, microwaves, neutrons
- G01L1/255—Measuring force or stress, in general using wave or particle radiation, e.g. X-rays, microwaves, neutrons using acoustic waves, or acoustic emission
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating 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/14—Investigating 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating 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/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/23—Dune restoration or creation; Cliff stabilisation
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Immunology (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Pathology (AREA)
- Acoustics & Sound (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Toxicology (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
The application relates to a rock characteristic stress determining method based on acoustic emission counting nonlinear evolution, and relates to the technical field of underground engineering. The method comprises the following steps: and carrying out an indoor rock uniaxial compression acoustic emission experiment on the rock sample to be tested, obtaining a stress-strain curve and an acoustic emission cumulative count curve, dividing the stress-strain curve and the acoustic emission cumulative count curve into a plurality of acoustic emission cumulative count sub-curves, and carrying out linear regression on each acoustic emission cumulative count sub-curve to obtain a linear curve corresponding to each acoustic emission cumulative count sub-curve. And generating a fitting degree curve between the acoustic emission accumulated count sub-curve and the corresponding linear curve according to each acoustic emission accumulated count sub-curve. And generating a counting change rate curve according to the linear curve corresponding to each acoustic emission accumulated counting sub-curve. And determining target characteristic strain according to the counting change rate curve and the fitting degree curve. And determining the rock characteristic stress of the rock sample to be tested according to the target characteristic strain and the stress-strain curve. Rock characteristic stresses can be determined using the present application.
Description
Technical Field
The application relates to the technical field of underground engineering, in particular to a rock characteristic stress determining method based on acoustic emission counting nonlinear evolution.
Background
The rock contains a large number of cracks inside, which are affected by complex geological conditions. When the rock is compressed, the state of cracks such as closure, expansion and penetration affect the macroscopic deformation of the rock. Therefore, the research on the relation between the crack state in the rock and the rock characteristic stress can provide scientific guidance for the stability evaluation of the rock engineering.
The conventional rock characteristic stress determination method is to divide rock characteristic stress through a compression broken rock volume strain curve and a crack volume strain curve. The rock characteristic stress includes the dilatation stress sigma d Crack closure stress sigma c And crack initiation stress sigma i . Wherein the expansion stress sigma d The crack closure stress sigma is the stress corresponding to the highest point of the rock volume strain curve c The crack initiation stress sigma is the stress corresponding to the starting point of the horizontal segment of the crack volume strain curve i Is the stress corresponding to the end point of the horizontal segment of the crack volume strain curve.
However, the conventional determination of rock characteristic stresses relies mainly on experience, and since rock is a heterogeneous and anisotropic material, the strain of the rock is not uniform throughout the compression fracture process, resulting in inaccurate classification and calculation of rock characteristic stresses.
Thus, a method for determining rock characteristic stress is needed.
Disclosure of Invention
Based on the above, it is necessary to provide a rock characteristic stress determining method based on acoustic emission count nonlinear evolution, aiming at the technical problems.
In a first aspect, there is provided a rock characteristic stress determination method based on acoustic emission count nonlinear evolution, the method comprising:
carrying out an indoor rock uniaxial compression acoustic emission experiment on a rock sample to be tested to obtain a stress-strain curve and an acoustic emission cumulative count curve of which acoustic emission count increases along with strain of the rock sample to be tested;
dividing the acoustic emission cumulative count curve into a plurality of acoustic emission cumulative count sub-curves, and performing linear regression on each acoustic emission cumulative count sub-curve to obtain a linear curve corresponding to each acoustic emission cumulative count sub-curve;
generating a fitting degree curve between each acoustic emission cumulative count sub-curve and the corresponding linear curve according to each acoustic emission cumulative count sub-curve and the corresponding linear curve;
generating a counting change rate curve of the acoustic emission counting along with the change of the strain according to the linear curve corresponding to each acoustic emission accumulated counting sub-curve;
Determining target characteristic strain of the rock sample to be tested according to the counting change rate curve and the fitness curve;
and determining the rock characteristic stress of the rock sample to be tested according to the target characteristic strain and the stress-strain curve.
As an optional implementation manner, the acoustic emission cumulative count curve is a nonlinear curve, the dividing the acoustic emission cumulative count curve into a plurality of acoustic emission cumulative count sub-curves, and performing linear regression on each acoustic emission cumulative count sub-curve to obtain a linear curve corresponding to each acoustic emission cumulative count sub-curve, including:
dividing the acoustic emission cumulative count curve into a plurality of acoustic emission cumulative count sub-curves according to a preset strain sliding window and a preset sliding step length;
and carrying out linear regression on the acoustic emission accumulation counting sub-curves aiming at each acoustic emission accumulation counting sub-curve to obtain a linear curve corresponding to the acoustic emission accumulation counting sub-curve.
As an optional implementation manner, the generating a fitness curve between each acoustic emission accumulation count sub-curve and the corresponding linear curve according to each acoustic emission accumulation count sub-curve and the corresponding linear curve includes:
Determining the fitting degree between the acoustic emission cumulative count sub-curve and the corresponding linear curve based on a preset curve fitting degree algorithm aiming at each acoustic emission cumulative count sub-curve and the corresponding linear curve;
and generating a fitting degree curve of which the fitting degree changes along with the strain according to the fitting degree between each acoustic emission cumulative count sub-curve and the corresponding linear curve.
As an optional implementation manner, the generating a count change rate curve of acoustic emission count along with strain according to the linear curve corresponding to each acoustic emission accumulation count sub-curve includes:
acquiring the slope of a linear curve corresponding to each acoustic emission accumulation counting sub-curve, and determining the slope as the counting change rate of the acoustic emission accumulation counting sub-curve corresponding to the linear curve;
and generating a counting change rate curve of the acoustic emission counting along with the change of the strain according to the counting change rate of each acoustic emission accumulated counting sub-curve.
As an alternative embodiment, the target characteristic strain includes a crack closure strain, a cracking strain, and a dilatation strain, and the determining the target characteristic strain of the rock sample to be measured according to the count change rate curve and the fitness curve includes:
Acquiring a first strain when both the counting change rate in the counting change rate curve and the fitting degree in the fitting degree curve reach a preset stable state, and determining the first strain as the crack closure strain;
acquiring a second strain when the fitting degree in the fitting degree curve is in the stable state and the counting change rate in the counting change rate curve enters the ascending state, and determining the second strain as the cracking strain;
and obtaining a third strain when the fitting degree in the fitting degree curve is changed from the stable state to the descending state, and determining the third strain as the dilatation strain.
As an alternative embodiment, the preset steady state is that the count change rate in the count change rate curve is within a preset count change rate threshold; and/or the number of the groups of groups,
and the fitting degree in the fitting degree curve is in a preset fitting degree threshold range.
As an alternative embodiment, the target characteristic strain includes a crack closure strain, a cracking strain, and a dilatation strain, the rock characteristic stress includes a crack closure stress, a cracking stress, and a dilatation stress, and determining the rock characteristic stress of the rock sample to be tested according to the target characteristic strain and the stress-strain curve includes:
Determining a first stress corresponding to the crack closure strain as the crack closure stress in the stress-strain curve;
determining a second stress corresponding to the cracking stress as the cracking stress in the stress-strain curve;
and in the stress-strain curve, determining a third stress corresponding to the dilatation strain as the dilatation stress.
In a second aspect, there is provided a rock characteristic stress determination device based on acoustic emission count nonlinear evolution, the device comprising:
the acquisition module is used for carrying out an indoor rock uniaxial compression acoustic emission experiment on a rock sample to be tested, and acquiring a stress-strain curve and an acoustic emission cumulative count curve of acoustic emission count increasing along with strain of the rock sample to be tested;
the first determining module is used for dividing the acoustic emission accumulation counting curve into a plurality of acoustic emission accumulation counting sub-curves, and carrying out linear regression on each acoustic emission accumulation counting sub-curve to obtain a linear curve corresponding to each acoustic emission accumulation counting sub-curve;
the first generation module is used for generating a fitting degree curve between each acoustic emission accumulation counting sub-curve and the corresponding linear curve according to each acoustic emission accumulation counting sub-curve and the corresponding linear curve;
The second generation module is used for generating a counting change rate curve of the acoustic emission count along with the change of the strain according to the linear curve corresponding to each acoustic emission accumulated counting sub-curve;
the second determining module is used for determining target characteristic strain of the rock sample to be tested according to the counting change rate curve and the fitting degree curve;
and the third determining module is used for determining the rock characteristic stress of the rock sample to be tested according to the target characteristic strain and the stress-strain curve.
As an alternative embodiment, the first determining module is specifically configured to:
dividing the acoustic emission cumulative count curve into a plurality of acoustic emission cumulative count sub-curves according to a preset strain sliding window and a preset sliding step length;
and carrying out linear regression on the acoustic emission accumulation counting sub-curves aiming at each acoustic emission accumulation counting sub-curve to obtain a linear curve corresponding to the acoustic emission accumulation counting sub-curve.
As an alternative embodiment, the first generating module is specifically configured to:
determining the fitting degree between the acoustic emission cumulative count sub-curve and the corresponding linear curve based on a preset curve fitting degree algorithm aiming at each acoustic emission cumulative count sub-curve and the corresponding linear curve;
And generating a fitting degree curve of which the fitting degree changes along with the strain according to the fitting degree between each acoustic emission cumulative count sub-curve and the corresponding linear curve.
As an alternative embodiment, the second generating module is specifically configured to:
acquiring the slope of a linear curve corresponding to each acoustic emission accumulation counting sub-curve, and determining the slope as the counting change rate of the acoustic emission accumulation counting sub-curve corresponding to the linear curve;
and generating a counting change rate curve of the acoustic emission counting along with the change of the strain according to the counting change rate of each acoustic emission accumulated counting sub-curve.
As an alternative embodiment, the target characteristic strain includes a crack closure strain, a cracking strain, and a dilatation strain, and the second determining module is specifically configured to:
acquiring a first strain when both the counting change rate in the counting change rate curve and the fitting degree in the fitting degree curve reach a preset stable state, and determining the first strain as the crack closure strain;
acquiring a second strain when the fitting degree in the fitting degree curve is in the stable state and the counting change rate in the counting change rate curve enters the ascending state, and determining the second strain as the cracking strain;
And obtaining a third strain when the fitting degree in the fitting degree curve is changed from the stable state to the descending state, and determining the third strain as the dilatation strain.
As an alternative embodiment, the preset steady state is that the count change rate in the count change rate curve is within a preset count change rate threshold; and/or the number of the groups of groups,
and the fitting degree in the fitting degree curve is in a preset fitting degree threshold range.
As an alternative embodiment, the target characteristic strain includes a crack closure strain, a cracking strain, and a dilatation strain, the rock characteristic stress includes a crack closure stress, a cracking stress, and a dilatation stress, and the third determining module is specifically configured to:
determining a first stress corresponding to the crack closure strain as the crack closure stress in the stress-strain curve;
determining a second stress corresponding to the cracking stress as the cracking stress in the stress-strain curve;
and in the stress-strain curve, determining a third stress corresponding to the dilatation strain as the dilatation stress.
In a third aspect, there is provided a computer device comprising a memory and a processor, the memory having stored thereon a computer program executable on the processor, the processor implementing the method steps according to any of the first aspects when the computer program is executed.
In a fourth aspect, there is provided a computer readable storage medium having stored thereon a computer program which, when executed by a processor, implements the method steps of any of the first aspects.
The application provides a rock characteristic stress determining method based on acoustic emission counting nonlinear evolution, and the technical scheme provided by the embodiment of the application at least brings the following beneficial effects: and carrying out an indoor rock uniaxial compression acoustic emission experiment on the rock sample to be tested, and obtaining a stress-strain curve and an acoustic emission cumulative count curve of which the acoustic emission count increases along with the strain of the rock sample to be tested. Dividing the acoustic emission accumulation counting curve into a plurality of acoustic emission accumulation counting sub-curves, and carrying out linear regression on each acoustic emission accumulation counting sub-curve to obtain a linear curve corresponding to each acoustic emission accumulation counting sub-curve. And generating a fitting degree curve between each acoustic emission cumulative count sub-curve and the corresponding linear curve according to each acoustic emission cumulative count sub-curve and the corresponding linear curve. And generating a counting change rate curve of the acoustic emission counting along with the change of the strain according to the linear curve corresponding to each acoustic emission accumulated counting sub-curve. And determining the target characteristic strain of the rock sample to be measured according to the counting change rate curve and the fitness curve. And determining the rock characteristic stress of the rock sample to be tested according to the target characteristic strain and the stress-strain curve. The rock corresponds to different characteristic stresses in a crack closing stage, a cracking stage and a capacity expansion stage respectively, and the method and the device utilize an indoor rock uniaxial compression acoustic emission experiment to acquire acoustic emission accumulated counts and reflect the crack activity degree of a rock sample to be tested according to the acoustic emission accumulated counts. The evolution process based on the acoustic emission cumulative count is a theoretical basis of a process of converting from linear to nonlinear, and the acoustic emission cumulative count curve is subjected to piecewise linear regression, so that the fitting degree of the linear curve after the linear regression and the acoustic emission cumulative count curve is judged in a piecewise manner. In addition, the slope of the linear curve after the linear return can represent the count change rate of the acoustic emission count along with the change of the strain, and therefore, the speed of the count change rate (namely the magnitude of the slope) can be determined. According to the method and the device, according to the counting change rate and the fitting degree of the acoustic emission cumulative count, the time for converting the linearity to the nonlinearity in the acoustic emission cumulative count curve can be accurately judged, namely the acoustic emission nonlinearity count evolution characteristic of the rock sample to be measured in the compression process is quantitatively calculated, and further, the rock characteristic stress is accurately determined.
It is to be understood that both the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the application.
Drawings
In order to more clearly illustrate the embodiments of the invention or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described, it being obvious that the drawings in the following description are only some embodiments of the invention, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of a rock characteristic stress determination method based on acoustic emission count nonlinear evolution provided by an embodiment of the present application;
FIG. 2 is a flow chart of another rock characteristic stress determination method based on acoustic emission count nonlinear evolution provided by an embodiment of the present application;
FIG. 3 is a flow chart of yet another rock characteristic stress determination method based on acoustic emission count nonlinear evolution provided by an embodiment of the present application;
FIG. 4 is a flow chart of yet another rock characteristic stress determination method based on acoustic emission count nonlinear evolution provided by an embodiment of the present application;
FIG. 5 is a flow chart of yet another rock characteristic stress determination method based on acoustic emission count nonlinear evolution provided by an embodiment of the present application;
FIG. 6 is a flow chart of yet another rock characteristic stress determination method based on acoustic emission count nonlinear evolution provided by an embodiment of the present application;
FIG. 7 is a schematic diagram of a stress path combined with acoustic emission cumulative count provided in an embodiment of the present application;
FIG. 8 is a graph of a fitness, count rate of change versus number and stress path profile provided in an embodiment of the present application;
fig. 9 is a schematic structural diagram of a rock characteristic stress determining device based on acoustic emission count nonlinear evolution according to an embodiment of the present application;
fig. 10 is a schematic structural diagram of a computer device according to an embodiment of the present application.
Detailed Description
In order to make the objects, technical solutions and advantages of the present application more apparent, the present application will be further described in detail with reference to the accompanying drawings and examples. It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the present application.
The following will describe a rock feature stress determining method based on acoustic emission count nonlinear evolution provided in the embodiment of the present application in detail with reference to a specific implementation manner, and fig. 1 is a flowchart of a rock feature stress determining method based on acoustic emission count nonlinear evolution provided in the embodiment of the present application, as shown in fig. 1, and specific steps are as follows:
And 101, performing an indoor rock uniaxial compression acoustic emission experiment on a rock sample to be tested to obtain a stress-strain curve and an acoustic emission cumulative count curve of which the acoustic emission count increases along with the strain of the rock sample to be tested.
In the implementation, a technician selects typical hard rock (such as red sandstone) as a rock sample to be tested, then develops an indoor rock uniaxial compression acoustic emission experiment, can adopt a displacement loading rate of 0.004mm/s to enable the rock sample to be tested to be damaged statically, and acquires acoustic emission counts of the rock sample to be tested, which change along with the strain in the compression process, in real time through acoustic emission monitoring equipment. After the experiment is finished, the computer generates a stress-strain curve of the rock sample to be tested. Fig. 7 is a schematic diagram of a combination of a stress path and acoustic emission cumulative count according to an embodiment of the present application, where, as shown in fig. 7, a solid line is a stress path, a horizontal axis is strain, a left vertical axis is stress, a stress-strain curve initially shows a linear rising trend after a downward convex growth, and finally shows an upward convex growth. The technician counts the collected acoustic emission counts through computer equipment, and generates a whole accumulated overall process diagram of the acoustic emission counts through a computer, namely an acoustic emission accumulated count curve of the acoustic emission counts increasing along with the strain, wherein a broken line in fig. 7 is the acoustic emission accumulated count curve, a horizontal axis is the strain, and a right vertical axis is a logarithmic value of the acoustic emission accumulated count. The acoustic emission cumulative count curve increases convexly, then linearly, and then the concave increasing trend finally increases suddenly.
And 102, dividing the acoustic emission cumulative count curve into a plurality of acoustic emission cumulative count sub-curves, and performing linear regression on each acoustic emission cumulative count sub-curve to obtain a linear curve corresponding to each acoustic emission cumulative count sub-curve.
In practice, since the acoustic emission cumulative count curve is a curve that evolves from linear to nonlinear, the embodiments of the present application divide the acoustic emission cumulative count curve into several acoustic emission cumulative count sub-curves. And (3) carrying out linear regression on each acoustic emission cumulative count sub-curve based on a preset linear regression method (such as least square regression, theil's regression or Siegel regression) to obtain a linear curve corresponding to each acoustic emission cumulative count sub-curve.
As an optional implementation manner, the acoustic emission cumulative count curve is a nonlinear curve, fig. 2 is a flowchart of another rock characteristic stress determining method based on acoustic emission cumulative count nonlinear evolution provided in the embodiment of the present application, as shown in fig. 2, in step 102, the acoustic emission cumulative count curve is divided into a plurality of acoustic emission cumulative count sub-curves, and linear regression is performed on each acoustic emission cumulative count sub-curve, so that the specific steps for obtaining the linear curve corresponding to each acoustic emission cumulative count sub-curve are as follows:
Step 201, dividing the acoustic emission cumulative count curve into a plurality of acoustic emission cumulative count sub-curves according to a preset strain sliding window and a preset sliding step length.
In practice, the technician may divide the acoustic emission cumulative count curve into several acoustic emission cumulative count sub-curves by a computer. Specifically, sliding interception can be performed in a sliding window mode according to a preset strain sliding window and a preset sliding step length. For example, the strain sliding window is 0.1 and the sliding step is 0.05. As shown in FIG. 7, the computer can determine an acoustic emission cumulative count curve with a strain value between 0 and 0.1 as a first acoustic emission cumulative count sub-curve, the computer can determine an acoustic emission cumulative count curve with a strain value between 0.05 and 0.15 as a second acoustic emission cumulative count sub-curve, the computer can determine an acoustic emission cumulative count curve with a strain value between 0.1 and 0.2 as a third acoustic emission cumulative count sub-curve, and so on, and finally divide the acoustic emission cumulative count curve into a plurality of acoustic emission cumulative count sub-curves.
Step 202, for each acoustic emission accumulation count sub-curve, performing linear regression on the acoustic emission accumulation count sub-curve to obtain a linear curve corresponding to the acoustic emission accumulation count sub-curve.
In implementation, for each acoustic emission cumulative count sub-curve, a technician performs linear regression on the acoustic emission cumulative count sub-curve through a computer and a preset linear regression algorithm to obtain a linear curve corresponding to the acoustic emission cumulative count sub-curve.
And 103, generating a fitting degree curve between each acoustic emission cumulative count sub-curve and the corresponding linear curve according to each acoustic emission cumulative count sub-curve and the corresponding linear curve.
In practice, a technician can generate a fitness curve between each acoustic emission cumulative count sub-curve and a corresponding linear curve from each acoustic emission cumulative count sub-curve and a corresponding linear curve through a computer. FIG. 8 is a graph of a fitness, count rate of change versus number and stress path profile provided in an embodiment of the present application. As shown in fig. 8, the short-dashed line in fig. 8 represents a fitness curve, the vertical axis (left side) of the coordinate system in which the fitness curve is located represents the fitness, and the horizontal axis represents the strain.
As an alternative implementation manner, fig. 3 is a flowchart of still another rock feature stress determining method based on the nonlinear evolution of acoustic emission counts, as shown in fig. 3, and in step 103, a specific step of generating, according to each acoustic emission cumulative count sub-curve and a corresponding linear curve, a fitness curve between each acoustic emission cumulative count sub-curve and a corresponding linear curve is as follows:
Step 301, determining, for each acoustic emission cumulative count sub-curve and a corresponding linear curve, a fitting degree between the acoustic emission cumulative count sub-curve and the corresponding linear curve based on a preset curve fitting degree algorithm.
In practice, for each acoustic emission cumulative count sub-curve and corresponding linear curve, a technician, through a computer, may determine the fit between the acoustic emission cumulative count sub-curve and the corresponding linear curve based on a preset curve fit algorithm.
Step 302, generating a fitting degree curve with the fitting degree changing along with the strain according to the fitting degree between each acoustic emission cumulative count sub-curve and the corresponding linear curve.
In the implementation, a technician fits the fitting degree between each acoustic emission cumulative count sub-curve and the corresponding linear curve into a fitting degree curve with the fitting degree changing along with the strain through a computer, and the technician can also carry out smoothing treatment on the fitting degree curve through the computer. As indicated by the short dashed line in fig. 8.
And 104, generating a counting change rate curve of the acoustic emission count along with the change of the strain according to the linear curve corresponding to each acoustic emission accumulated counting sub-curve.
In the implementation, a technician generates a counting change rate curve of the acoustic emission count along with the change of the strain according to a linear curve corresponding to each acoustic emission accumulated counting sub-curve through a computer, as shown in fig. 8, a long dashed line in fig. 8 is the counting change rate curve, a vertical axis (right side) of a coordinate system where the counting change rate curve is located is a counting change rate logarithmic value, and a horizontal axis is the strain.
As an optional implementation manner, fig. 4 is a flowchart of still another rock feature stress determining method based on the nonlinear evolution of the acoustic emission count, as shown in fig. 4, and in step 104, according to a linear curve corresponding to each acoustic emission cumulative count sub-curve, a specific step of generating a count change rate curve of the acoustic emission count changing with strain is as follows:
step 401, for each acoustic emission accumulation count sub-curve corresponding linear curve, acquiring a slope of the linear curve, and determining the slope as a count change rate of the acoustic emission accumulation count sub-curve corresponding to the linear curve.
In an implementation, the slope of the linear curve corresponding to the acoustic emission cumulative count sub-curve is the count change rate of the acoustic emission cumulative count. The technician can obtain the slope of the linear curve corresponding to each acoustic emission cumulative count sub-curve through the computer, and determine the slope as the count change rate of the acoustic emission cumulative count sub-curve corresponding to the linear curve.
Step 402, generating a counting change rate curve of acoustic emission counting along with the change of the strain according to the counting change rate of each acoustic emission accumulation counting sub-curve.
In implementation, a technician can fit the count change rate of each acoustic emission cumulative count sub-curve to a count change rate curve of acoustic emission count change with strain through a computer, and the technician can also perform smoothing on the count change rate curve through the computer, as shown by a long-dashed line in fig. 8.
And 105, determining the target characteristic strain of the rock sample to be tested according to the counting change rate curve and the fitting degree curve.
In practice, the computer determines the target characteristic strain of the rock sample to be measured according to the count change rate curve and the fitness curve.
As an alternative implementation manner, the target characteristic strain includes crack closure strain, cracking strain and dilatation strain, fig. 5 is a flowchart of still another rock characteristic stress determining method based on acoustic emission counting nonlinear evolution provided in the embodiment of the present application, and as shown in fig. 5, the specific steps of determining the target characteristic strain of the rock sample to be measured according to the counting change rate curve and the fitness curve in step 105 are as follows:
Step 501, a first strain when both the count change rate in the count change rate curve and the fitting degree in the fitting degree curve reach a preset stable state is obtained, and the first strain is determined to be a crack closure strain.
In an implementation, the computer may obtain a first strain when both the count rate of change in the count rate of change curve and the fitness in the fitness curve reach a preset steady state, and determine the first strain as a crack closure strain. When the fitting degree curve reaches a stable state, the acoustic emission cumulative count sub-curve and the corresponding linear curve reach higher fitting degree, and further the acoustic emission cumulative count sub-curve at the moment is in a linear state. When the counting change rate curve reaches a stable state, the slope of the acoustic emission accumulation counting curve is almost unchanged. The starting point of the acoustic emission cumulative count curve entering the linear stage can be determined through two indexes of the count change rate and the fitting degree, and the strain corresponding to the starting point is the crack closure strain. As shown in fig. 8, in the early stage of compression of the rock sample to be measured, the indexes such as stress and acoustic emission count have large fluctuation, and the strain value starts to be gradually stable from 0.3351, so that the strain value 0.3351 can be determined as crack closure strain.
Step 502, obtaining a second strain when the fitting degree in the fitting degree curve is in a stable state and the counting change rate in the counting change rate curve enters an ascending state, and determining the second strain as the cracking strain.
In an implementation, the computer may obtain a second strain when the fitness in the fitness curve is in a steady state and the count rate of change in the count rate of change curve enters a rising state, and determine the second strain as the cracking strain. The fitting degree curve keeps a stable state, which shows that the fitting degree of the acoustic emission accumulation counting sub-curve and the corresponding linear curve is kept higher, the counting change rate curve starts to rise, and the slope of the acoustic emission accumulation counting curve increases, so that the acoustic emission accumulation counting increases faster. The strain at this time is the cracking strain. As shown in fig. 8, the strain value is 1.8955 to 2.2300, the fitting degree curve is kept stable, and the count change rate curve is slightly increased, so that the strain value 1.8955 can be determined as the crack initiation strain.
In step 503, a third strain when the fitness in the fitness curve changes from the steady state to the decreasing state is obtained, and the third strain is determined as the dilatation strain.
In practice, the computer may obtain a third strain in the fitness curve when the fitness changes from a steady state to a decreasing state and determine the third strain as the dilatation strain. The fitting degree curve is changed from a stable state to a descending state, which indicates that the fitting degree of the acoustic emission accumulation counting sub-curve and the corresponding linear curve is reduced, further indicates that the acoustic emission accumulation counting curve is converted from linear to nonlinear, the counting change rate curve is continuously increased, indicates that the slope of the acoustic emission accumulation counting curve is further increased, the acoustic emission accumulation counting is further increased, and the strain at the moment is the dilatation strain. As shown in fig. 8, the strain value is from 2.2300 to 2.9188, and the fitting degree in the fitting degree curve is changed from a steady state to a sharp decrease, and thus, the strain value 2.2300 can be determined as the dilatation strain.
As an alternative embodiment, the preset steady state is that the count change rate in the count change rate curve is within a preset count change rate threshold; and/or the number of the groups of groups,
the fitting degree in the fitting degree curve is in a preset fitting degree threshold range.
In practice, the technician may set the count change rate threshold range and the fitness threshold range in advance in the computer. When the counting change rate in the counting change rate curve is within a preset fitness threshold range, the computer judges that the counting change rate in the counting change rate curve is in a stable state. Similarly, when the fitting degree in the fitting degree curve is within a preset threshold range of the fitting degree, the computer judges that the fitting degree in the fitting degree curve is in a stable state.
And 106, determining the rock characteristic stress of the rock sample to be tested according to the target characteristic strain and the stress-strain curve.
In practice, the computer may determine rock characteristic stresses corresponding to the target characteristic strain in a stress-strain curve based on the target characteristic strain.
As an alternative implementation manner, the target characteristic strain includes a crack closure strain, a cracking strain and a dilatation strain, the rock characteristic stress includes a crack closure stress, a cracking stress and a dilatation stress, fig. 6 is a flowchart of yet another rock characteristic stress determining method based on acoustic emission count nonlinear evolution provided in the embodiment of the present application, as shown in fig. 6, in step 106, a specific step of determining a rock characteristic stress of a rock sample to be tested according to the target characteristic strain and a stress-strain curve is as follows:
In step 601, a first stress corresponding to the crack closure strain is determined as a crack closure stress in a stress-strain curve.
In implementations, the computer may determine a first stress corresponding to the crack closure strain as the crack closure stress in the stress-strain curve. As shown in fig. 8, the solid line in the figure represents the stress-strain curve, the horizontal axis represents strain, and the vertical axis represents stress. The computer may determine the crack closure stress as 11.817Mpa in the stress-strain curve based on the crack closure strain value 0.3351.
In step 602, in the stress-strain curve, the second stress corresponding to the determined cracking strain is determined as the cracking stress.
In implementations, the computer may determine a first stress corresponding to the crack closure strain as the crack closure stress in the stress-strain curve. As shown in fig. 8, the computer may determine that the cracking stress is 75.78Mpa in the stress-strain curve based on the cracking strain 1.8955.
In step 603, in the stress-strain curve, the third stress corresponding to the determined dilatation strain is determined as the dilatation stress.
In implementations, the computer may determine a first stress corresponding to the crack closure strain as the crack closure stress in the stress-strain curve. As shown in fig. 8, the computer may determine that the dilatation stress is 89.58MPa in the stress-strain curve based on the dilatation strain 2.2300.
The embodiment of the application provides a rock characteristic stress determining method based on acoustic emission counting nonlinear evolution, which comprises the following steps: and carrying out an indoor rock uniaxial compression acoustic emission experiment on the rock sample to be tested to obtain a stress-strain curve and an acoustic emission cumulative count curve of which the acoustic emission count increases along with the strain of the rock sample to be tested. Dividing the acoustic emission cumulative count curve into a plurality of acoustic emission cumulative count sub-curves, and performing linear regression on each acoustic emission cumulative count sub-curve to obtain a linear curve corresponding to each acoustic emission cumulative count sub-curve. And generating a fitting degree curve between each acoustic emission accumulation counting sub-curve and the corresponding linear curve according to each acoustic emission accumulation counting sub-curve and the corresponding linear curve. And generating a counting change rate curve of the acoustic emission counting along with the change of the strain according to the linear curve corresponding to each acoustic emission accumulated counting sub-curve. And determining the target characteristic strain of the rock sample to be measured according to the counting change rate curve and the fitting degree curve. And determining the rock characteristic stress of the rock sample to be tested according to the target characteristic strain and the stress-strain curve. The rock corresponds to different characteristic stresses in a crack closing stage, a cracking stage and a capacity expansion stage respectively, and the method and the device utilize an indoor rock uniaxial compression acoustic emission experiment to acquire acoustic emission accumulated counts and reflect the crack activity degree of a rock sample to be tested according to the acoustic emission accumulated counts. The evolution process based on the acoustic emission cumulative count is a theoretical basis of a process of converting from linear to nonlinear, and the acoustic emission cumulative count curve is subjected to piecewise linear regression, so that the fitting degree of the linear curve after the linear regression and the acoustic emission cumulative count curve is judged in a piecewise manner. In addition, the slope of the linear curve after the linear return can represent the count change rate of the acoustic emission count along with the change of the strain, and therefore, the speed of the count change rate (namely the magnitude of the slope) can be determined. According to the method and the device, according to the counting change rate and the fitting degree of the acoustic emission cumulative count, the time for converting the linearity to the nonlinearity in the acoustic emission cumulative count curve can be accurately judged, namely the acoustic emission nonlinearity count evolution characteristic of the rock sample to be measured in the compression process is quantitatively calculated, and further, the rock characteristic stress is accurately determined.
It should be understood that, although the steps in the flowcharts of fig. 1 to 6 are sequentially shown as indicated by arrows, these steps are not necessarily sequentially performed in the order indicated by the arrows. The steps are not strictly limited to the order of execution unless explicitly recited herein, and the steps may be executed in other orders. Moreover, at least a portion of the steps of fig. 1-6 may include steps or stages that are not necessarily performed at the same time, but may be performed at different times, nor does the order in which the steps or stages are performed necessarily occur sequentially, but may be performed alternately or alternately with other steps or at least a portion of the steps or stages in other steps.
It should be understood that the same/similar parts of the embodiments of the method described above in this specification may be referred to each other, and each embodiment focuses on differences from other embodiments, and references to descriptions of other method embodiments are only needed.
The embodiment of the application also provides a rock characteristic stress determining device based on acoustic emission counting nonlinear evolution, as shown in fig. 9, the device comprises:
The acquisition module 910 is configured to perform an indoor rock uniaxial compression acoustic emission experiment on a rock sample to be tested, and acquire a stress-strain curve and an acoustic emission cumulative count curve of the rock sample to be tested, where the acoustic emission count increases along with the strain;
the first determining module 920 is configured to divide the acoustic emission cumulative count curve into a plurality of acoustic emission cumulative count sub-curves, and perform linear regression on each acoustic emission cumulative count sub-curve to obtain a linear curve corresponding to each acoustic emission cumulative count sub-curve;
a first generating module 930, configured to generate a fitness curve between each acoustic emission cumulative count sub-curve and the corresponding linear curve according to each acoustic emission cumulative count sub-curve and the corresponding linear curve;
the second generating module 940 is configured to generate a count change rate curve of the acoustic emission count according to the linear curve corresponding to each acoustic emission accumulated count sub-curve;
a second determining module 950, configured to determine a target characteristic strain of the rock sample to be measured according to the count rate of change curve and the fitness curve;
a third determining module 960 is configured to determine a rock characteristic stress of the rock sample to be tested according to the target characteristic strain and the stress-strain curve.
As an alternative embodiment, the first determining module 920 is specifically configured to:
dividing the acoustic emission cumulative count curve into a plurality of acoustic emission cumulative count sub-curves according to a preset strain sliding window and a preset sliding step length;
and carrying out linear regression on the acoustic emission cumulative count sub-curve aiming at each acoustic emission cumulative count sub-curve to obtain a linear curve corresponding to the acoustic emission cumulative count sub-curve.
As an alternative embodiment, the first generating module 930 is specifically configured to:
determining the fitting degree between each acoustic emission cumulative count sub-curve and the corresponding linear curve based on a preset curve fitting degree algorithm aiming at each acoustic emission cumulative count sub-curve and the corresponding linear curve;
and generating a fitting degree curve of which the fitting degree changes along with the strain according to the fitting degree between each acoustic emission accumulation counting sub-curve and the corresponding linear curve.
As an alternative embodiment, the second generating module 940 is specifically configured to:
acquiring the slope of a linear curve corresponding to each acoustic emission accumulation counting sub-curve, and determining the slope as the counting change rate of the acoustic emission accumulation counting sub-curve corresponding to the linear curve;
And generating a counting change rate curve of the acoustic emission counting along with the change of the strain according to the counting change rate of each acoustic emission accumulated counting sub-curve.
As an alternative embodiment, the target characteristic strain includes a crack closure strain, a crack initiation strain, and a dilatation strain, and the second determining module 950 is specifically configured to:
acquiring a first strain when the counting change rate in the counting change rate curve and the fitting degree in the fitting degree curve reach a preset stable state, and determining the first strain as a crack closure strain;
acquiring a second strain when the fitting degree in the fitting degree curve is in a stable state and the counting change rate in the counting change rate curve enters an ascending state, and determining the second strain as a cracking strain;
and acquiring a third strain when the fitting degree in the fitting degree curve is changed from the stable state to the descending state, and determining the third strain as the dilatation strain.
As an alternative embodiment, the preset steady state is that the count change rate in the count change rate curve is within a preset count change rate threshold; and/or the number of the groups of groups,
the fitting degree in the fitting degree curve is in a preset fitting degree threshold range.
As an alternative embodiment, the target characteristic strain includes a crack closure strain, a crack initiation strain, and a dilatation strain, the rock characteristic stress includes a crack closure stress, a crack initiation stress, and a dilatation stress, and the third determining module 960 is specifically configured to:
in the stress-strain curve, determining a first stress corresponding to the crack closure strain as a crack closure stress;
in the stress-strain curve, determining a second stress corresponding to the cracking stress as the cracking stress;
in the stress-strain curve, the third stress corresponding to the determined dilatation strain is determined as the dilatation stress.
The embodiment of the application provides a rock characteristic stress determining device based on acoustic emission counting nonlinear evolution, which comprises: the acquisition module 910 is configured to perform an indoor rock uniaxial compression acoustic emission experiment on a rock sample to be tested, and acquire a stress-strain curve and an acoustic emission cumulative count curve of the rock sample to be tested, where the acoustic emission count increases with strain. The first determining module 920 is configured to divide the acoustic emission cumulative count curve into a plurality of acoustic emission cumulative count sub-curves, and perform linear regression on each acoustic emission cumulative count sub-curve to obtain a linear curve corresponding to each acoustic emission cumulative count sub-curve. The first generating module 930 is configured to generate a fitness curve between each acoustic emission cumulative count sub-curve and the corresponding linear curve according to each acoustic emission cumulative count sub-curve and the corresponding linear curve. The second generating module 940 is configured to generate a count change rate curve of the acoustic emission count according to the linear curve corresponding to each acoustic emission cumulative count sub-curve. A second determining module 950, configured to determine a target characteristic strain of the rock sample to be measured according to the count change rate curve and the fitness curve. A third determining module 960 is configured to determine a rock characteristic stress of the rock sample to be tested according to the target characteristic strain and the stress-strain curve. The rock corresponds to different characteristic stresses in a crack closing stage, a cracking stage and a capacity expansion stage respectively, and the method and the device utilize an indoor rock uniaxial compression acoustic emission experiment to acquire acoustic emission accumulated counts and reflect the crack activity degree of a rock sample to be tested according to the acoustic emission accumulated counts. The evolution process based on the acoustic emission cumulative count is a theoretical basis of a process of converting from linear to nonlinear, and the acoustic emission cumulative count curve is subjected to piecewise linear regression, so that the fitting degree of the linear curve after the linear regression and the acoustic emission cumulative count curve is judged in a piecewise manner. In addition, the slope of the linear curve after the linear return can represent the count change rate of the acoustic emission count along with the change of the strain, and therefore, the speed of the count change rate (namely the magnitude of the slope) can be determined. According to the method and the device, according to the counting change rate and the fitting degree of the acoustic emission cumulative count, the time for converting the linearity to the nonlinearity in the acoustic emission cumulative count curve can be accurately judged, namely the acoustic emission nonlinearity count evolution characteristic of the rock sample to be measured in the compression process is quantitatively calculated, and further, the rock characteristic stress is accurately determined.
For specific limitations of the rock characteristic stress determination device based on the acoustic emission count nonlinear evolution, reference may be made to the above limitation of the rock characteristic stress determination method based on the acoustic emission count nonlinear evolution, and the description thereof will not be repeated here. The above-mentioned rock characteristic stress determining device based on acoustic emission count nonlinear evolution can be implemented by all or part of software, hardware and their combination. The above modules may be embedded in hardware or may be independent of a processor in the computer device, or may be stored in software in a memory in the computer device, so that the processor may call and execute operations corresponding to the above modules.
In one embodiment, a computer device is provided, as shown in fig. 10, including a memory and a processor, where the memory stores a computer program that can be run on the processor, and the processor implements the method steps of rock compression deformation crack characteristic stress determination based on acoustic emission count nonlinear evolution when executing the computer program.
In one embodiment, a computer readable storage medium has stored thereon a computer program which, when executed by a processor, implements the steps of the above method of rock compression deformation crack signature determination based on acoustic emission count nonlinear evolution.
Those skilled in the art will appreciate that implementing all or part of the above described methods may be accomplished by way of a computer program stored on a non-transitory computer readable storage medium, which when executed, may comprise the steps of the embodiments of the methods described above. Any reference to memory, storage, database, or other medium used in the various embodiments provided herein may include non-volatile and/or volatile memory. The nonvolatile memory can include Read Only Memory (ROM), programmable ROM (PROM), electrically Programmable ROM (EPROM), electrically Erasable Programmable ROM (EEPROM), or flash memory. Volatile memory can include Random Access Memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in a variety of forms such as Static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double Data Rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous Link DRAM (SLDRAM), memory bus direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM), among others.
It is noted that relational terms such as first and second, and the like are used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. Moreover, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising one … …" does not exclude the presence of other like elements in a process, method, article, or apparatus that comprises the element.
It should be noted that, user information (including but not limited to user equipment information, user personal information, etc.) and data (including but not limited to data for presentation, analyzed data, etc.) referred to in the present application are information and data authorized by the user or sufficiently authorized by each party.
In this specification, each embodiment is described in a related manner, and identical and similar parts of each embodiment are all referred to each other, and each embodiment mainly describes differences from other embodiments. In particular, for system embodiments, since they are substantially similar to method embodiments, the description is relatively simple, as relevant to see a section of the description of method embodiments.
The technical features of the above embodiments may be arbitrarily combined, and all possible combinations of the technical features in the above embodiments are not described for brevity of description, however, as long as there is no contradiction between the combinations of the technical features, they should be considered as the scope of the description.
The above examples merely represent a few embodiments of the present application, which are described in more detail and are not to be construed as limiting the scope of the invention. It should be noted that it would be apparent to those skilled in the art that various modifications and improvements could be made without departing from the spirit of the present application, which would be within the scope of the present application. Accordingly, the scope of protection of the present application is to be determined by the claims appended hereto.
Claims (10)
1. A rock characteristic stress determination method based on acoustic emission count nonlinear evolution, the method comprising:
carrying out an indoor rock uniaxial compression acoustic emission experiment on a rock sample to be tested to obtain a stress-strain curve and an acoustic emission cumulative count curve of which acoustic emission count increases along with strain of the rock sample to be tested;
dividing the acoustic emission cumulative count curve into a plurality of acoustic emission cumulative count sub-curves, and performing linear regression on each acoustic emission cumulative count sub-curve to obtain a linear curve corresponding to each acoustic emission cumulative count sub-curve;
generating a fitting degree curve between each acoustic emission cumulative count sub-curve and the corresponding linear curve according to each acoustic emission cumulative count sub-curve and the corresponding linear curve;
generating a counting change rate curve of the acoustic emission counting along with the change of the strain according to the linear curve corresponding to each acoustic emission accumulated counting sub-curve;
determining target characteristic strain of the rock sample to be tested according to the counting change rate curve and the fitness curve;
and determining the rock characteristic stress of the rock sample to be tested according to the target characteristic strain and the stress-strain curve.
2. The method of claim 1, wherein the acoustic emission cumulative count curve is a nonlinear curve, the dividing the acoustic emission cumulative count curve into a plurality of acoustic emission cumulative count sub-curves, and performing linear regression on each of the acoustic emission cumulative count sub-curves to obtain a linear curve corresponding to each of the acoustic emission cumulative count sub-curves, comprising:
dividing the acoustic emission cumulative count curve into a plurality of acoustic emission cumulative count sub-curves according to a preset strain sliding window and a preset sliding step length;
and carrying out linear regression on the acoustic emission accumulation counting sub-curves aiming at each acoustic emission accumulation counting sub-curve to obtain a linear curve corresponding to the acoustic emission accumulation counting sub-curve.
3. The method of claim 1, wherein said generating a fitness curve between each of said acoustic emission cumulative count sub-curves and a corresponding linear curve from each of said acoustic emission cumulative count sub-curves and a corresponding linear curve comprises:
determining the fitting degree between the acoustic emission cumulative count sub-curve and the corresponding linear curve based on a preset curve fitting degree algorithm aiming at each acoustic emission cumulative count sub-curve and the corresponding linear curve;
And generating a fitting degree curve of which the fitting degree changes along with the strain according to the fitting degree between each acoustic emission cumulative count sub-curve and the corresponding linear curve.
4. The method of claim 1, wherein generating a count rate of change curve of acoustic emission count as a function of strain from the linear curve corresponding to each of the acoustic emission cumulative count sub-curves comprises:
acquiring the slope of a linear curve corresponding to each acoustic emission accumulation counting sub-curve, and determining the slope as the counting change rate of the acoustic emission accumulation counting sub-curve corresponding to the linear curve;
and generating a counting change rate curve of the acoustic emission counting along with the change of the strain according to the counting change rate of each acoustic emission accumulated counting sub-curve.
5. The method of claim 1, wherein the target characteristic strain comprises a crack closure strain, a cracking strain, and a dilatation strain, the determining the target characteristic strain of the rock sample under test from the count rate of change curve and the fitness curve comprising:
acquiring a first strain when both the counting change rate in the counting change rate curve and the fitting degree in the fitting degree curve reach a preset stable state, and determining the first strain as the crack closure strain;
Acquiring a second strain when the fitting degree in the fitting degree curve is in the stable state and the counting change rate in the counting change rate curve enters the ascending state, and determining the second strain as the cracking strain;
and obtaining a third strain when the fitting degree in the fitting degree curve is changed from the stable state to the descending state, and determining the third strain as the dilatation strain.
6. The method of claim 5, wherein the predetermined steady state is a count rate of change in the count rate of change curve being within a predetermined count rate of change threshold; and/or the number of the groups of groups,
and the fitting degree in the fitting degree curve is in a preset fitting degree threshold range.
7. The method of claim 1, wherein the target characteristic strain comprises a crack closure strain, a cracking strain, and a dilatation strain, the rock characteristic stress comprises a crack closure stress, a cracking stress, and a dilatation stress, and the determining the rock characteristic stress of the rock sample to be tested from the target characteristic strain and the stress-strain curve comprises:
determining a first stress corresponding to the crack closure strain as the crack closure stress in the stress-strain curve;
Determining a second stress corresponding to the cracking stress as the cracking stress in the stress-strain curve;
and in the stress-strain curve, determining a third stress corresponding to the dilatation strain as the dilatation stress.
8. Rock characteristic stress determining device based on acoustic emission count nonlinear evolution, characterized in that the device comprises:
the acquisition module is used for carrying out an indoor rock uniaxial compression acoustic emission experiment on a rock sample to be tested, and acquiring a stress-strain curve and an acoustic emission cumulative count curve of acoustic emission count increasing along with strain of the rock sample to be tested;
the first determining module is used for dividing the acoustic emission accumulation counting curve into a plurality of acoustic emission accumulation counting sub-curves, and carrying out linear regression on each acoustic emission accumulation counting sub-curve to obtain a linear curve corresponding to each acoustic emission accumulation counting sub-curve;
the first generation module is used for generating a fitting degree curve between each acoustic emission accumulation counting sub-curve and the corresponding linear curve according to each acoustic emission accumulation counting sub-curve and the corresponding linear curve;
the second generation module is used for generating a counting change rate curve of the acoustic emission count along with the change of the strain according to the linear curve corresponding to each acoustic emission accumulated counting sub-curve;
The second determining module is used for determining target characteristic strain of the rock sample to be tested according to the counting change rate curve and the fitting degree curve;
and the third determining module is used for determining the rock characteristic stress of the rock sample to be tested according to the target characteristic strain and the stress-strain curve.
9. A computer device comprising a memory and a processor, the memory having stored thereon a computer program executable on the processor, characterized in that the processor implements the steps of the method according to any of claims 1 to 7 when the computer program is executed.
10. A computer readable storage medium, on which a computer program is stored, characterized in that the computer program, when being executed by a processor, implements the steps of the method of any of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310760691.4A CN116539413B (en) | 2023-06-27 | 2023-06-27 | Rock characteristic stress determining method based on acoustic emission counting nonlinear evolution |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310760691.4A CN116539413B (en) | 2023-06-27 | 2023-06-27 | Rock characteristic stress determining method based on acoustic emission counting nonlinear evolution |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116539413A true CN116539413A (en) | 2023-08-04 |
| CN116539413B CN116539413B (en) | 2023-11-07 |
Family
ID=87458054
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310760691.4A Active CN116539413B (en) | 2023-06-27 | 2023-06-27 | Rock characteristic stress determining method based on acoustic emission counting nonlinear evolution |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116539413B (en) |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101179134B1 (en) * | 2012-03-23 | 2012-09-07 | 국방과학연구소 | System for Measuring High Velocity Impact Acoustic Emissions and the Method Therefor |
| CN103852377A (en) * | 2013-09-24 | 2014-06-11 | 核工业北京地质研究院 | Method for recognizing uniaxial compression initiation crack strength of rock based on accumulative sound emission impacting number |
| CN106405675A (en) * | 2016-08-25 | 2017-02-15 | 山东科技大学 | Dynamic monitoring system and method for early warning against slope slide of tailing pond of strip mining pit |
| CN109374682A (en) * | 2018-11-26 | 2019-02-22 | 中国工程物理研究院化工材料研究所 | A kind of monitoring device of fragile material crack initiation time |
| CN110987674A (en) * | 2019-12-23 | 2020-04-10 | 神华神东煤炭集团有限责任公司 | Crustal stress testing method based on core Kaiser effect |
| US20210088429A1 (en) * | 2019-09-24 | 2021-03-25 | Southwest Petroleum University | Device and Method for Evaluating Fracture Initiation and Propagation, and Stress Sensitivity of Propped Fracture |
| CN113108732A (en) * | 2021-04-13 | 2021-07-13 | 江西理工大学 | Guided wave monitoring method for slope landslide early warning |
| CN113218766A (en) * | 2021-05-27 | 2021-08-06 | 江西理工大学 | Rock crack initiation stress and damage stress identification method based on moment tensor analysis |
| CN113776943A (en) * | 2021-11-05 | 2021-12-10 | 中国矿业大学(北京) | Rock compressive strength prediction method |
| CN114279842A (en) * | 2021-12-27 | 2022-04-05 | 中国科学院武汉岩土力学研究所 | Method and system for determining rock crack initiation stress and damage stress |
| CN115856092A (en) * | 2023-01-30 | 2023-03-28 | 北京科技大学 | Method for determining rock crack initiation stress based on acoustic emission data and stress data |
-
2023
- 2023-06-27 CN CN202310760691.4A patent/CN116539413B/en active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101179134B1 (en) * | 2012-03-23 | 2012-09-07 | 국방과학연구소 | System for Measuring High Velocity Impact Acoustic Emissions and the Method Therefor |
| CN103852377A (en) * | 2013-09-24 | 2014-06-11 | 核工业北京地质研究院 | Method for recognizing uniaxial compression initiation crack strength of rock based on accumulative sound emission impacting number |
| CN106405675A (en) * | 2016-08-25 | 2017-02-15 | 山东科技大学 | Dynamic monitoring system and method for early warning against slope slide of tailing pond of strip mining pit |
| CN109374682A (en) * | 2018-11-26 | 2019-02-22 | 中国工程物理研究院化工材料研究所 | A kind of monitoring device of fragile material crack initiation time |
| US20210088429A1 (en) * | 2019-09-24 | 2021-03-25 | Southwest Petroleum University | Device and Method for Evaluating Fracture Initiation and Propagation, and Stress Sensitivity of Propped Fracture |
| CN110987674A (en) * | 2019-12-23 | 2020-04-10 | 神华神东煤炭集团有限责任公司 | Crustal stress testing method based on core Kaiser effect |
| CN113108732A (en) * | 2021-04-13 | 2021-07-13 | 江西理工大学 | Guided wave monitoring method for slope landslide early warning |
| CN113218766A (en) * | 2021-05-27 | 2021-08-06 | 江西理工大学 | Rock crack initiation stress and damage stress identification method based on moment tensor analysis |
| CN113776943A (en) * | 2021-11-05 | 2021-12-10 | 中国矿业大学(北京) | Rock compressive strength prediction method |
| CN114279842A (en) * | 2021-12-27 | 2022-04-05 | 中国科学院武汉岩土力学研究所 | Method and system for determining rock crack initiation stress and damage stress |
| CN115856092A (en) * | 2023-01-30 | 2023-03-28 | 北京科技大学 | Method for determining rock crack initiation stress based on acoustic emission data and stress data |
Non-Patent Citations (3)
| Title |
|---|
| QINGSONG WANG ET AL.: "Acoustic emission characteristics and energy mechanism in karst limestone failure under uniaxial and triaxial compression", 《 BULLETIN OF ENGINEERING GEOLOGY AND THE ENVIRONMENT》, pages 1427 - 1442 * |
| 张昕 等: "单轴压缩下砂岩声发射及分形特征研究", 《中国测试》, vol. 43, no. 02, pages 13 - 19 * |
| 郭清露 等: "大理岩热损伤声发射力学特性试验研究", 《岩石力学与工程学报》, vol. 34, no. 12, pages 2388 - 2400 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116539413B (en) | 2023-11-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5582211B1 (en) | Stress-strain relationship simulation method, springback amount prediction method, and springback analysis device | |
| JP5866892B2 (en) | Stress-strain relationship evaluation method and springback amount prediction method | |
| López-Querol et al. | Drained cyclic behaviour of loose Dogs Bay sand | |
| Behnia et al. | Damage detection of SFRC concrete beams subjected to pure torsion by integrating acoustic emission and Weibull damage function | |
| Katsanos et al. | Prediction of inelastic response periods of buildings based on intensity measures and analytical model parameters | |
| CN113486570A (en) | Method for predicting tunnel seismic vulnerability based on random IDA and machine learning | |
| CN116539413B (en) | Rock characteristic stress determining method based on acoustic emission counting nonlinear evolution | |
| CN115081221A (en) | Method for establishing rock nonlinear creep model based on fractional derivative | |
| KR101851075B1 (en) | Reliability improvement method of finite element analysis of the stress relaxation prediction through the metal material heat treatment | |
| CN110907296A (en) | Method for identifying flow times of dynamic creep test of asphalt mixture | |
| CN116227008B (en) | Prestress anchor bolt support analysis method and device based on surrounding rock deformation control | |
| JP2006337343A (en) | Estimation system of true stress-logarithmic strain curve of structure comprising crystalline polymer | |
| CN112859192A (en) | Rock core saturation prediction model construction method and rock core saturation prediction method | |
| CN116499881A (en) | Method for establishing rock theoretical damage evolution model | |
| CN112364517B (en) | Method for evaluating service life of cylinder cover under low cycle fatigue and related equipment | |
| CN108616315B (en) | Power amplifier output power detection method and device, computer equipment and storage medium | |
| CN116106125B (en) | Method for evaluating brittleness degree of rock | |
| Hinnerichs et al. | A hybrid experimental-numerical procedure for determining creep crack growth rates | |
| CN117686309B (en) | Rock mass property-based rock stratum maximum horizontal principal stress prediction method | |
| CN110274836B (en) | Creep strength testing method based on acoustic emission activity analysis | |
| Weathers et al. | Automated Determination of Felicity Ratio for Composite Overwrapped Pressure Vessels | |
| CN114781278B (en) | Electromechanical equipment service life prediction method and system based on data driving | |
| CN116297868B (en) | Method for calculating rock burst acoustic emission precursor characteristics | |
| CN114386525B (en) | Adaptive adjustment method, device, equipment and storage medium for evaluation standard | |
| Levin et al. | Making Shear Simple–Validation of the Shear Compression Specimen 0 (SCS0) for Shear Testing |
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 |