CN116184500A - Real-time inversion method and device for ground stress of tunnel based on microseismic information - Google Patents
Real-time inversion method and device for ground stress of tunnel based on microseismic information Download PDFInfo
- Publication number
- CN116184500A CN116184500A CN202310325332.6A CN202310325332A CN116184500A CN 116184500 A CN116184500 A CN 116184500A CN 202310325332 A CN202310325332 A CN 202310325332A CN 116184500 A CN116184500 A CN 116184500A
- Authority
- CN
- China
- Prior art keywords
- ground stress
- stress
- microseismic
- fault plane
- tensor
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000012544 monitoring process Methods 0.000 claims abstract description 6
- 238000013507 mapping Methods 0.000 claims description 9
- 238000012216 screening Methods 0.000 claims description 7
- 238000004590 computer program Methods 0.000 claims description 6
- 238000001514 detection method Methods 0.000 abstract description 5
- 238000010586 diagram Methods 0.000 description 18
- 238000005259 measurement Methods 0.000 description 4
- 239000011435 rock Substances 0.000 description 4
- 230000006835 compression Effects 0.000 description 3
- 238000007906 compression Methods 0.000 description 3
- 241000220317 Rosa Species 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 208000010392 Bone Fractures Diseases 0.000 description 1
- 206010017076 Fracture Diseases 0.000 description 1
- 208000013201 Stress fracture Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a real-time inversion method and a device for ground stress of a tunnel based on microseismic information, belonging to the technical field of tunnel safety monitoring, wherein the method comprises the following steps: constructing a spatial distribution model of each microseismic event in a tunnel region to be detected; determining a seismic source mechanism solution of each elastic wave according to the microseism event space distribution model; according to the seismic source mechanism solution of each elastic wave in the tunnel region to be detected, solving a first ground stress tensor of the fault plane; performing P/T axis to ground stress linear inversion according to the azimuth angle and the inclination angle of the first ground stress tensor to determine a second ground stress tensor; and determining the principal stress of each direction of the fault plane according to the shape ratio of the second ground stress tensor. The method solves the problems of how to improve the detection precision of the ground stress of the tunnel area and how to feed back the front ground stress condition in real time.
Description
Technical Field
The invention belongs to the technical field of tunnel safety monitoring, and particularly relates to a real-time inversion method and device for ground stress of a tunnel based on microseismic information.
Background
With the development of tunnel engineering, the ground stress is taken as data for judging the stability of tunnels and providing important support for the construction and support of tunnels, and the measurement requirement thereof is also rapidly increasing.
The traditional ground stress measuring method such as hydraulic fracturing, stress relieving and the like can only measure the stress magnitude and direction of one point location, but cannot obtain the distribution condition of all ground stress fields of a certain area, so that the ground stress detection precision of a tunnel area is not high, and meanwhile, the traditional ground stress measuring method cannot feed back the front ground stress condition in real time along with the progress of construction.
Disclosure of Invention
The invention aims to provide a real-time inversion method and device for ground stress of a tunnel based on microseismic information, which are used for solving the problems of how to improve the detection precision of the ground stress of a tunnel area and feeding back the condition of the front ground stress in real time.
The invention adopts the following technical scheme:
the invention provides a real-time inversion method of ground stress of a tunnel based on microseismic information, which comprises the following steps:
s101, constructing a spatial distribution model of each microseismic event in a tunnel region to be detected;
s102, determining a seismic source mechanism solution of each elastic wave according to a microseism event space distribution model;
s103, according to the seismic source mechanism solution of each elastic wave in the tunnel region to be detected, solving a first ground stress tensor of the fault plane;
s104, performing P/T axis to ground stress linear inversion according to the azimuth angle and the inclination angle of the first ground stress tensor to determine a second ground stress tensor;
s105, determining the principal stress of each direction of the fault plane according to the shape ratio of the second ground stress tensor.
Optionally, constructing the spatial distribution model of each microseismic event in the tunnel region to be detected includes:
acquiring the position information and the pressure intensity of each microseismic event in a tunnel region to be detected through a microseismic monitoring system;
and constructing a microseismic event space distribution model of each microseismic event according to the position information and the pressure intensity of each microseismic event.
Optionally, determining a source mechanism solution for each elastic wave according to the microseismic event spatial distribution model includes:
and analyzing the microseismic event space distribution model to determine a seismic source mechanism solution of each fault plane, wherein the seismic source mechanism solution comprises the trend, the sliding direction and the inclination angle of the fault plane.
Optionally, determining the type of the fault plane according to the trend, the sliding direction and the inclination angle of the fault plane by a seismic source mechanism.
Optionally, according to a source mechanism solution of each elastic wave in the tunnel region to be detected, the solving the first ground stress tensor of the fault plane includes:
according to the seismic source mechanism solution of each elastic wave in the tunnel region to be detected, determining the normal vector, the tangential stress direction vector and the sliding vector of the fault plane;
and then solving a first ground stress tensor of the fault plane according to the fault plane normal vector, the tangential stress direction vector and the sliding vector.
Optionally, determining the second ground stress tensor by performing the P/T axis-to-ground stress linear inversion according to the azimuth angle and the inclination angle of the first ground stress tensor includes:
according to the azimuth angle and the inclination angle of each component of the first ground stress tensor, mapping each component of the first ground stress tensor onto a beach ball to form each corresponding mapping point;
then determining the effective component of the first ground stress tensor according to the position relation between each mapping point and the known observation P point and the known observation T point on the beach ball;
and finally forming a second ground stress tensor according to the effective component.
Optionally, determining the principal stress of each direction of the fault plane from the shape ratio of the second ground stress tensor includes:
firstly, calculating a first difference value between the maximum principal stress and the intermediate principal stress in the second ground stress tensor;
then calculating a second difference value between the maximum principal stress and the minimum principal stress;
then calculating the ratio of the first difference value to the second difference value to obtain a shape ratio;
and finally, screening out a group of stresses with the shape ratio not being 1, and taking the stress as the stress corresponding to each direction of the fault plane.
The invention also provides a real-time inversion device of the ground stress of the tunnel based on the microseismic information, which comprises a memory, a processor and a computer program which is stored in the memory and can run on the processor, wherein the real-time inversion method of the ground stress of the tunnel based on the microseismic information is realized when the processor executes the computer program.
The beneficial effects of the invention are as follows: according to the method, all microseismic events of the whole tunnel region to be detected are detected and analyzed, a plurality of elastic wave equations related to each microseismic event are established through P-wave inversion, and all seismic source mechanism solutions with unknown quantity less than that of an equation set are solved; further, accurate ground stress tensors are determined through P/T analysis and shape ratio analysis according to the seismic source mechanism solution, the problem of how to improve the detection accuracy of the ground stress is solved through the method provided by the invention, and the technical effect of feeding back the front ground stress condition in real time is achieved.
Drawings
FIG. 1 is a flow chart of a real-time inversion method of the ground stress of a tunnel based on microseismic information;
FIG. 2 is a schematic diagram of a spatial distribution model of a tunnel microseismic event according to the present invention;
fig. 3 is a schematic diagram of a mechanism solution of a microseism focus according to the present invention, in which:
(a) forward off-slip, (b) forward fault, (c) forward off-slip, (d) reverse fault, (e) reverse fault, (f) forward fault, (g) forward fault, and (h) forward slip fault;
fig. 4 is a rose diagram of a source mechanism according to the present invention, wherein:
(a), (b) and (c) are rose figures of the fault plane 1 with respect to azimuth, inclination and sliding angle, respectively, (d), (e) and (f) are rose figures of the fault plane 2 with respect to azimuth, inclination and sliding angle, respectively;
FIG. 5 is a schematic illustration of a principal stress multi-solution provided by the present invention;
FIG. 6 is a schematic diagram of the relationship between the P/T axis and the principal stress provided by the present invention;
fig. 7 is a linear ground stress inversion rose diagram provided by the invention, wherein:
(a) a rose diagram of the maximum principal stress of the fault plane 1, (b) a rose diagram of the intermediate principal stress of the fault plane 1 and (c) a rose diagram of the minimum principal stress of the fault plane 1;
fig. 8 is a schematic diagram of a real-time inversion device of the ground stress of a tunnel based on microseismic information.
Detailed Description
The invention will be described in detail below with reference to the drawings and the detailed description.
The invention provides a real-time inversion method of the ground stress of a tunnel based on microseismic information, as shown in figure 1, comprising the following steps:
and step 101, constructing a spatial distribution model of each microseismic event in the tunnel region to be detected.
In one embodiment, by selecting a tunnel region to be detected at the tunnel surrounding rock, the region may select a tunnel region of a cube about 15 meters from the tunnel entrance, for example: 5*5. The geological formation data for each tunnel zone is shown in table 1, where the microseismic monitoring system(s) is/are installed, and any rock mass will typically produce many tiny micro-fractures prior to macroscopic damage. The micro-cracks generate elastic waves in an elastic wave performance release mode, a micro-vibration system arranged on the rock mass is monitored in real time, and the moment, the position and the property of the occurrence of the micro-cracks of the rock mass are obtained through an inversion method. And the method can capture multiple elastic waves emitted by all the seismic sources in the area, and screen out effective elastic waves by comparing the elastic waves with microseismic waveforms in a database. And constructing a microseismic event spatial distribution model according to the screened effective elastic waves, wherein the microseismic event spatial distribution model can reflect the position information of each seismic source and the distribution condition of each seismic source, as shown in fig. 2, each ball in the diagram represents one seismic source, the distribution is more dense, and the correlation between the seismic sources with a longer distance is smaller.
Table 1 geological formation data for tunnel region
Step 102, determining a seismic source mechanism solution of each elastic wave according to the micro-seismic event space distribution model.
In one embodiment, the physical process or source physical process of the occurrence of the earthquake, referred to as the source mechanism, may be determined from a seismic record of a plurality of seismic stations or microseismic systems, and some characteristic quantities at the source or some physical quantities of the source physical process at the time of the occurrence of the earthquake, referred to as the source parameters. The parameters of the seismic source comprise the trend, the tendency and the inclination angle of a fault surface of the seismic source, the direction and the amplitude of the dislocation of two plates of the fault of the seismic source, the length and the width of the fault of the seismic source, the expansion speed of fault fracture, the main stress state of the seismic source and the like, and can be obtained through the approaches of combining macroscopic seismic measurement with the analysis of a microseism event space distribution model and the like through a seismic source mechanism. The method comprises the steps of determining a seismic source mechanism solution of a fault plane by analyzing macroscopic seismic survey, a microseismic event space distribution model and elastic waves, wherein the seismic source mechanism solution comprises a seismic source mechanism and seismic source parameters, and judging the type of the fault plane by the trend, the tendency and the dip angle of the fault plane in the seismic source parameters in the seismic source mechanism solution when the description is needed.
And step 103, according to the seismic source mechanism solution of each elastic wave in the tunnel region to be detected, solving a first ground stress tensor of the fault plane.
In one embodiment, the fault plane related geologic structure information analyzed from macroscopic seismic measurements, microseismic event spatial distribution models, and elastic waves, such as fault plane type, may be forward or reverse sliding, forward fault, etc. Taking a normal fault as an example, by judging the fault plane where the microseismic event is located, measuring the fault plane normal vector, the fault plane tangential stress direction vector, the fault plane sliding vector and the fault plane stress tensor of the fault plane, establishing equations about the fault plane normal vector, the fault plane tangential stress direction vector, the fault plane sliding vector and the normal fault plane stress tensor, and solving the first ground stress tensor of the fault plane. The specifically established equation is shown in formulas (1) - (3), and the first ground stress tensor of the fault plane is solved.
The first ground stress tensor solved here is the resultant force of the principal stresses in each direction of the fault plane. The solved seismic source mechanism solutions have a plurality of groups of solutions, the solutions of different groups correspond to one ground stress tensor, the equations establish a joint equation through a plurality of seismic sources of the tunnel area to be monitored, and an equation set with the number of equations larger than the unknown number is established, so that the plurality of groups of seismic source mechanism solutions can be solved. On the one hand, when the more seismic sources of the tunnel region to be detected are considered, the more accurate the solved seismic source mechanism solution is, and the more accurate the first ground stress tensor of the solved fault plane is; on the other hand, by solving the source mechanism solution of all the microseismic events in the tunnel region, compared with the traditional single source mechanism solution for solving one microseismic event at a certain point, the efficiency of judging each fault plane of the region is improved.
And 104, performing P/T axis to ground stress linear inversion according to the azimuth angle and the inclination angle of the first ground stress tensor to determine a second ground stress tensor.
In one embodiment, according to the first ground stress tensor of the solved fault plane, in order to further detect whether the solved first ground stress meets the azimuth requirement, a P-wave analysis is performed on each microseismic source, and the P-wave analysis is described on a P-wave initial motion quadrant distribution diagram, wherein the distribution diagram comprises a P-axis (a maximum main stress axis) at a position of a bisecting expansion quadrant and a position T-axis (a minimum main stress axis) at a position of a compression quadrant, the ground stress is linearly inverted by performing a P/T axis on the first ground stress, the expansion and compression quadrants are calculated and analyzed on the P-wave initial motion quadrant distribution diagram through a point description, and whether a P point corresponding to the first ground stress tensor is located between two actually measured P observation points and a T point corresponding to the first ground stress tensor is located between two actually measured T observation points. If the P point of the first ground stress tensor used for description is positioned between two observation P points and the T point of the first ground stress tensor is positioned between two observation T points, the first ground stress tensor accords with the azimuth requirement, so that the first ground stress which does not accord with the azimuth requirement in the seismic source mechanism solution is removed; and the first ground stress tensor is selected to be in accordance with the P/T axis distribution.
And 105, determining the principal stress of each direction of the fault plane according to the shape ratio of the second ground stress tensor.
In one embodiment, it is to be noted that each component of the second ground stress tensor is a main stress in each direction on the fault plane, and the maximum main stress, the intermediate main stress and the minimum main stress of the second ground stress are determined according to the order of the magnitudes of the stresses, respectively, and the shape ratio of the second ground stress tensor is obtained by calculating a first difference value of the maximum main stress and the intermediate main stress, a second difference value of the maximum main stress and the minimum main stress, and calculating a ratio of the first difference value and the second difference value. And screening a group of maximum principal stress, intermediate principal stress and minimum principal stress which meet the shape ratio not equal to 1. And screening out the maximum principal stress, the intermediate principal stress and the minimum principal stress as principal stresses in each direction with higher fault plane precision.
Optionally, constructing the spatial distribution model of each microseismic event in the tunnel region to be detected includes:
acquiring the position information and the pressure intensity of each microseismic event in a tunnel region to be detected through a microseismic monitoring system;
and constructing a microseismic event space distribution model of each microseismic event according to the position information and the pressure intensity of each microseismic event.
In one embodiment, by establishing a 3D model of the tunnel region to be detected, the position information of each microseismic event is converted into coordinate information of each point in the 3D model according to the actually detected position information, and each point contains pressure information of the microseismic event, and a microseismic event spatial distribution model of each microseismic event through containing the position information and the pressure information is established, as shown in fig. 2.
Optionally, determining a source mechanism solution for each elastic wave according to the microseismic event spatial distribution model includes:
and analyzing the microseismic event space distribution model to determine a seismic source mechanism solution of each fault plane, wherein the seismic source mechanism solution comprises a trend, a sliding direction and an inclination angle.
Optionally, determining the type of the fault plane according to the trend, the sliding direction and the inclination angle of the fault plane by a seismic source mechanism.
In one embodiment, geologic structure analysis is performed on each microseismic event through a microseismic event space distribution model, the parameters of the seismic source such as the trend, the sliding direction and the inclination angle of the seismic source fault plane in the seismic source mechanism solution are solved, and beach balls of the corresponding type of fault plane are drawn. As shown in fig. 3, the types of the resolved fault planes generally include (a 1) forward sliding, (b 1) forward sliding, (c 1) forward sliding, (d 1) reverse sliding, (e 1) reverse sliding, (f 1) forward sliding, (g 1) forward sliding, and (h 1) forward sliding.
It should be noted that, in order to more clearly study the characteristics of the elastic wave of each microseismic event in each seismic source mechanism, a rose diagram is used to perform statistical analysis on the azimuth angle, the inclination angle and the sliding angle in the seismic source mechanism, as shown in fig. 4, fig. 4 illustrates two rose diagrams with a turn-off layer 1 and a fault layer 2 determined according to the analyzed seismic source mechanism respectively, the azimuth angle, the inclination angle and the sliding angle of each fault layer 1 are respectively described by (a 2), (b 2) and (c 2), corresponding angle degrees and amplitude values can be read in each diagram, and the rose diagrams of the fault layers 2 and the rose diagrams of the fault layers 1 are respectively consistent in description.
Optionally, according to a source mechanism solution of each elastic wave in the tunnel region to be detected, the solving the first ground stress tensor of the fault plane includes:
according to the seismic source mechanism solution of each elastic wave in the tunnel region to be detected, determining the normal vector, the tangential stress direction vector and the sliding vector of the fault plane;
and then solving a first ground stress tensor of the fault plane according to the fault plane normal vector, the tangential stress direction vector and the sliding vector.
In one embodiment, the corresponding type of ground stress tensor equation of the fault plane is established according to the resolved focus mechanism, and each component of the first ground stress tensor is respectively represented by delta i I=1, 2,3 …, i representing the component number.
Assuming that the principal stress on the fault plane is as shown in formula (1):
δ i ＝τ ij n i n j (1)
wherein delta i Is the principal stress of the first stress tensor, n i Is the normal vector of the ith fault plane, n j Is the normal vector of the j-th fault plane.
τN i ＝T i -δ i (2)
Wherein T is i As the first stress tensor of the fault plane, N i The vector of the shear stress on the fault plane is τ, which is the shear stress of the fault plane.
Assuming that the sliding vector s of the fault plane is identical to the tangential stress direction of the fault plane, and the following formula is obtained after normalization processing:
At＝S (3)
wherein t is a vector of a first ground stress tensor, andeach component satisfying sigma 1 +σ 2 +σ 3 =0, the matrix with coefficient matrix a of 3X5 is:
s is a unit vector of the sliding vector of the fault plane. When there are M sets of source mechanism solutions, equation (3) becomes 3M equations for 5 unknowns of ground stress, and the desired first ground stress solution can be obtained from the generalized linear inversion of the L2 norm.
Based on the above formula (1) -formula (3), the first ground stress tensor vector is solved to be t= [ sigma ] 11 σ 12 σ 13 σ 22 σ 23 ]。
At this time, t is a plurality of sets of solutions, and the specific first ground stress inversion multi-solution is shown in table 2:
TABLE 2 stress inversion results
Wherein sigma 1 Is the maximum principal stress of the first ground stress, sigma 2 Is the intermediate principal stress of the first ground stress, sigma 3 As shown in fig. 5, since each component of the first ground stress is analyzed from different orientations, the first ground stress has multiple resolvability and sigma due to internal angle redundancy 1 Either 30 or 70. According to table 2, the polynomials of the first ground stress are expressed as: a set of solutions of sigma 1 Azimuth angle N89.87 ° E, inclination 3.01 °; sigma (sigma) 2 Azimuth angle N4.77 ° W, inclination angle 7.59 °; sigma (sigma) 3 Azimuth angle N14.21 ° W, inclination angle 6.93 °; another set of solutions is sigma 1 Azimuth angle N89.87 ° E, inclination 85.86 °; sigma (sigma) 2 Azimuth angle N4.77 ° W, inclination angle 7.59 °; sigma (sigma) 3 Azimuth angle N14.21 ° W, inclination angle 46.11 °, etc.
Optionally, determining the second ground stress tensor by performing the P/T axis-to-ground stress linear inversion according to the azimuth angle and the inclination angle of the first ground stress tensor includes:
according to the azimuth angle and the inclination angle of each component of the first ground stress tensor, mapping each component of the first ground stress tensor onto a beach ball to form each corresponding mapping point;
then determining the effective component of the first ground stress tensor according to the position relation between each mapping point and the known observation P point and the known observation T point on the beach ball;
and finally forming a second ground stress tensor according to the effective component.
In one embodiment, based on the source mechanism solution, the first ground stress has the property of multiple solutions, in order to determine each relatively accurate component in the first ground stress, the P/T axis to ground stress linear inversion is performed on the first ground stress, the multiple solutions map each component of the first ground stress on a beach ball in the form of a point according to azimuth angle and inclination angle, and the position relation between the known observation P point and the observation T point on the beach solution is analyzed. As shown in FIG. 6, on a beach ball, one component of the first principal stress is the maximum principal stress σ 1 Is located between the positions of two preset observation points P, indicating sigma 1 The observation point P is a point corresponding to macroscopic seismic measurement; similarly minimum sigma 3 Between the positions of two preset observation points T, indicating sigma 3 In the compressional wave quadrant, the observation point T is also the point corresponding to the macroscopic seismic survey.
Determining the maximum principal stress sigma by screening out the first ground stress component meeting the expansion wave quadrant and the compression wave quadrant 1 Intermediate principal stress sigma 2 And a minimum principal stress sigma 3 And forming a second ground stress. And screening out the effective components of the first ground stress meeting the azimuth requirement, removing the data which do not meet the requirement, qualitatively measuring the geological features of the fault plane, and improving the calculation accuracy of the first ground stress.
Optionally, determining the principal stress of each direction of the fault plane from the shape ratio of the second ground stress tensor includes:
firstly, calculating a first difference value between the maximum principal stress and the intermediate principal stress in the second ground stress tensor;
then calculating a second difference value between the maximum principal stress and the minimum principal stress;
then calculating the ratio of the first difference value to the second difference value to obtain a shape ratio;
and finally, screening out a group of stresses with the shape ratio not being 1, and taking the stress as the stress corresponding to each direction of the fault plane.
In one embodiment, the maximum principal stress σ is determined first 1 Intermediate principal stress sigma 2 And minimum principal stress sigma 3 A plurality of parameters, determining a shape ratio R, the expression of R being as shown in formula (4):
R＝(σ 1 -σ 2 )/(σ 1 -σ 3 ) (4)
wherein (sigma) 1 -σ 2 ) Characterizing a first difference between the calculated maximum principal stress and the intermediate principal stress, (σ) 1 -σ 3 ) And characterizing and calculating a second difference value of the maximum principal stress and the minimum principal stress, wherein R is the shape ratio.
The closer the main stress and the minimum main stress are, the limit r=1 is taken, and the tension states of the minimum main stress and the main stress are consistent; when R is relatively small, the maximum principal stress and the intermediate principal stress approach each other, and when r=0 is taken as the limit, the compressive state of the intermediate principal stress matches the maximum principal stress. When R is not equal to 1, it is stated that the region is faulted and that according to the determined maximum principal stress sigma 1 Intermediate principal stress sigma 2 And minimum principal stress sigma 3 To determine important parameters of geologic features of a fault plane.
In fig. 7, the maximum principal stresses σ of segment level 1 are depicted by rosettes, respectively 1 As shown in fig. (a 3), the section plane 1 has an inter-principal stress σ 2 As shown in the graph (b 3), and the minimum principal stress sigma of the layer 1 3 As shown in fig. (c 3), each graph depicts the distribution of azimuth angles and inclination angles corresponding to the principal stress.
Finally determining the more accurate maximum principal stress sigma 1 Intermediate principal stress sigma 2 And minimum principal stress sigma 3 As shown in table 3:
TABLE 3 principal stress inversion results Table
The geological features of the fault plane can be quantitatively measured through the shape ratio, and the detection accuracy of the fault plane is improved.
The present invention also provides a real-time inversion apparatus 800 of the ground stress of the tunnel based on the microseismic information, as shown in fig. 8, the apparatus 800 includes a memory 810, a processor 820, and a computer program 830 stored in the memory and executable on the processor, and the processor 820 implements a real-time inversion method of the ground stress of the tunnel based on the microseismic information according to any one of the above embodiments when executing the computer program 830.
It should be noted that, the real-time inversion method of the ground stress of the tunnel based on the microseismic information, which can be implemented by the device 800, is consistent with the above method embodiments, and will not be described herein.
Claims (8)
1. The real-time inversion method for the ground stress of the tunnel based on the microseismic information is characterized by comprising the following steps of:
s101, constructing a spatial distribution model of each microseismic event in a tunnel region to be detected;
s102, determining a seismic source mechanism solution of each elastic wave according to the microseism event space distribution model;
s103, according to the seismic source mechanism solution of each elastic wave in the tunnel region to be detected, solving a first ground stress tensor of the fault plane;
s104, performing P/T axis to ground stress linear inversion according to the azimuth angle and the inclination angle of the first ground stress tensor to determine a second ground stress tensor;
s105, determining the principal stress of each direction of the fault plane according to the shape ratio of the second ground stress tensor.
2. The method for real-time inversion of the ground stress of a tunnel based on microseismic information according to claim 1, wherein the constructing a spatial distribution model of each microseismic event in the tunnel region to be detected comprises:
acquiring the position information and the pressure intensity of each microseismic event in a tunnel region to be detected through a microseismic monitoring system;
and constructing a microseismic event space distribution model of each microseismic event according to the position information and the pressure intensity of each microseismic event.
3. The method of real-time inversion of the ground stress of a tunnel based on microseismic information of claim 1, wherein determining a source mechanism solution for each elastic wave from the microseismic event spatial distribution model comprises:
and analyzing the microseismic event space distribution model to determine a seismic source mechanism solution of each fault plane, wherein the seismic source mechanism solution comprises the trend, the sliding direction and the inclination angle of the fault plane.
4. A method of real-time inversion of the ground stress of a tunnel based on microseismic information as claimed in claim 3, wherein the type of fault plane is determined from the direction of the fault plane, the sliding direction and the dip angle of the source mechanism solution.
5. The method of real-time inversion of the ground stress of a tunnel based on microseismic information according to claim 1, wherein solving the first ground stress tensor of the fault plane according to the source mechanism solution of each elastic wave in the tunnel region to be detected comprises:
according to the seismic source mechanism solution of each elastic wave in the tunnel region to be detected, determining a normal vector, a tangential stress direction vector and a sliding vector of the fault plane;
and then solving a first ground stress tensor of the fault plane according to the fault plane normal vector, the shear stress direction vector and the sliding vector.
6. The method of real-time inversion of the ground stress of a tunnel based on microseismic information according to claim 1, wherein the step of performing the linear inversion of the P/T axis to the ground stress according to the azimuth angle and the inclination angle of the first ground stress tensor to determine the second ground stress tensor comprises:
mapping each component of the first ground stress tensor onto a beach ball according to the azimuth angle and the inclination angle of each component of the first ground stress tensor to form each corresponding mapping point;
then determining the effective component of the first ground stress tensor according to the position relation between each mapping point and the known observation P point and the known observation T point on the beach ball;
and finally, forming the second ground stress tensor according to the effective component.
7. The method of real-time inversion of the ground stress of a tunnel based on microseismic information of claim 1 wherein determining the principal stress of each direction of the fault plane from the shape ratio of the second ground stress tensor comprises:
firstly, calculating a first difference value between the maximum principal stress and the intermediate principal stress in the second ground stress tensor;
then calculating a second difference value between the maximum principal stress and the minimum principal stress;
then calculating the ratio of the first difference value to the second difference value to obtain the shape ratio;
and finally screening out a group of stresses with the shape ratio not being 1, and taking the group of stresses as the main stresses corresponding to each direction of the fault plane.
8. A real-time inversion device of the ground stress of a tunnel based on microseismic information, which is characterized by comprising a memory, a processor and a computer program stored in the memory and capable of running on the processor, wherein the processor realizes the real-time inversion method of the ground stress of the tunnel based on microseismic information according to any one of claims 1 to 7 when executing the computer program.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310325332.6A CN116184500A (en) | 2023-03-30 | 2023-03-30 | Real-time inversion method and device for ground stress of tunnel based on microseismic information |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310325332.6A CN116184500A (en) | 2023-03-30 | 2023-03-30 | Real-time inversion method and device for ground stress of tunnel based on microseismic information |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116184500A true CN116184500A (en) | 2023-05-30 |
Family
ID=86449008
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310325332.6A Pending CN116184500A (en) | 2023-03-30 | 2023-03-30 | Real-time inversion method and device for ground stress of tunnel based on microseismic information |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116184500A (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104536047A (en) * | 2015-01-04 | 2015-04-22 | 奥菲（北京）石油技术有限公司 | Technology for applying microseism to rock burst monitoring |
| CN105093274A (en) * | 2014-05-07 | 2015-11-25 | 中国石油化工股份有限公司 | Hydrofracture focal mechanism inversion method and system |
| US20170254909A1 (en) * | 2016-03-02 | 2017-09-07 | Microseismic, Inc. | Method for determining maximum horizontal stress magnitude and direction using microseismic derived fracture attributes and its application to evaluating hydraulic fracture stimulation induced stress changes |
| CN107807385A (en) * | 2017-09-18 | 2018-03-16 | 成都理工大学 | A kind of regional stress field change in time and space inversion method |
| WO2022053945A1 (en) * | 2020-09-08 | 2022-03-17 | Matrix Jvco Ltd | Artificial intelligence tool for geostress and geostrain analysis |
| CN114779330A (en) * | 2022-04-26 | 2022-07-22 | 中国矿业大学 | Mining working face main fracture azimuth analysis and prediction method based on microseismic monitoring |
-
2023
- 2023-03-30 CN CN202310325332.6A patent/CN116184500A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105093274A (en) * | 2014-05-07 | 2015-11-25 | 中国石油化工股份有限公司 | Hydrofracture focal mechanism inversion method and system |
| CN104536047A (en) * | 2015-01-04 | 2015-04-22 | 奥菲（北京）石油技术有限公司 | Technology for applying microseism to rock burst monitoring |
| US20170254909A1 (en) * | 2016-03-02 | 2017-09-07 | Microseismic, Inc. | Method for determining maximum horizontal stress magnitude and direction using microseismic derived fracture attributes and its application to evaluating hydraulic fracture stimulation induced stress changes |
| CN107807385A (en) * | 2017-09-18 | 2018-03-16 | 成都理工大学 | A kind of regional stress field change in time and space inversion method |
| WO2022053945A1 (en) * | 2020-09-08 | 2022-03-17 | Matrix Jvco Ltd | Artificial intelligence tool for geostress and geostrain analysis |
| CN114779330A (en) * | 2022-04-26 | 2022-07-22 | 中国矿业大学 | Mining working face main fracture azimuth analysis and prediction method based on microseismic monitoring |
Non-Patent Citations (2)
| Title |
|---|
| 任青山;艾德春;刘鸿;: "震源机制解反演与实测相结合研究平顶山矿区地应力分布特征", 煤炭技术, vol. 34, no. 10, 31 October 2015 (2015-10-31) * |
| 唐超等: "基于微震监测的应力场反演方法及其工程应用研究", 岩土工程学报, vol. 43, no. 9, 30 September 2021 (2021-09-30), pages 1730 - 1736 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US5996726A (en) | System and method for determining the distribution and orientation of natural fractures | |
| DONG et al. | Three-dimensional analytical solution of acoustic emission or microseismic source location under cube monitoring network | |
| RU2494418C1 (en) | Method of measuring coordinates of microseismic sources and parameters of mechanisms of centres thereof in conditions of strong seismic interference (versions) | |
| MXPA05006835A (en) | Methods for determining formation and borehole parameters using fresnel volume tomography. | |
| MX2011003851A (en) | Time reverse imaging operators for source location. | |
| Peters et al. | Application of prediction analysis to hypocenter determination using a local array | |
| Shang et al. | Acoustic emission source location from P-wave arrival time corrected data and virtual field optimization method | |
| US20160161620A1 (en) | System and method for estimating repeatability using base data | |
| Wang et al. | Micromechanisms of coal fracture: Insights from quantitative AE technique | |
| CN105549068A (en) | 3D anisotropic micro seismic interference inverse positioning method and 3D anisotropic micro seismic interference inverse positioning system | |
| WO2014164970A2 (en) | Methods and systems for locating seismic events | |
| WO2017107862A1 (en) | Method and device for determining overall measurement of seismic observation system repeatability | |
| Dai et al. | Focal mechanism determination for microseismic events and its application to the left bank slope of the Baihetan hydropower station in China | |
| US6980482B2 (en) | Seismic survey system | |
| Li et al. | Dynamic Time Warping Distance Method for Similarity Test of Multipoint Ground Motion Field. | |
| CN116184500A (en) | Real-time inversion method and device for ground stress of tunnel based on microseismic information | |
| Abel et al. | An analytic model for microseismic event location estimate accuracy | |
| CN113608262B (en) | Seismic data processing method and device for calculating rotation component by using translation component | |
| Favorskaya | A novel method for wave phenomena investigation | |
| Basu et al. | On the design of a dense array to extract rotational components of earthquake ground motion | |
| Leśniak et al. | Combined mine tremors source location and error evaluation in the Lubin Copper Mine (Poland) | |
| CN104406681B (en) | A kind of method of testing determining microseism velocity of wave in real time | |
| Huo et al. | A study of impact point detecting method based on seismic signal | |
| Yang et al. | A combined technique of implantable sensors and probabilistic localization method for monitoring acoustic events on concrete slab | |
| Rui et al. | A robust triaxial localization method of AE source using refraction path |
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 |