CN116432296A - Method and device for calculating minimum safe thickness of crossing movable fracture waterproof rock mass - Google Patents
Method and device for calculating minimum safe thickness of crossing movable fracture waterproof rock mass Download PDFInfo
- Publication number
- CN116432296A CN116432296A CN202310667132.9A CN202310667132A CN116432296A CN 116432296 A CN116432296 A CN 116432296A CN 202310667132 A CN202310667132 A CN 202310667132A CN 116432296 A CN116432296 A CN 116432296A
- Authority
- CN
- China
- Prior art keywords
- fracture
- fractal
- rock mass
- type
- expansion
- 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
- 239000011435 rock Substances 0.000 title claims abstract description 201
- 238000000034 method Methods 0.000 title claims abstract description 40
- 230000003068 static effect Effects 0.000 claims abstract description 51
- 230000000977 initiatory effect Effects 0.000 claims abstract description 45
- 238000004364 calculation method Methods 0.000 claims abstract description 27
- 150000001875 compounds Chemical class 0.000 claims abstract description 23
- 230000000149 penetrating effect Effects 0.000 claims abstract description 12
- 239000004973 liquid crystal related substance Substances 0.000 claims description 64
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 49
- 238000005336 cracking Methods 0.000 claims description 25
- 230000006870 function Effects 0.000 claims description 17
- 230000002968 anti-fracture Effects 0.000 claims description 16
- 230000006835 compression Effects 0.000 claims description 14
- 238000007906 compression Methods 0.000 claims description 14
- 239000002131 composite material Substances 0.000 claims description 9
- 238000012937 correction Methods 0.000 claims description 9
- 238000005326 angular distribution function Methods 0.000 claims description 7
- 238000009933 burial Methods 0.000 claims description 7
- 238000004590 computer program Methods 0.000 claims description 7
- 238000012512 characterization method Methods 0.000 claims description 6
- 238000003860 storage Methods 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 5
- 230000005484 gravity Effects 0.000 claims description 4
- 238000005516 engineering process Methods 0.000 abstract description 7
- 206010017076 Fracture Diseases 0.000 description 267
- 208000010392 Bone Fractures Diseases 0.000 description 262
- 238000004891 communication Methods 0.000 description 7
- 230000008859 change Effects 0.000 description 5
- 238000010008 shearing Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 238000009412 basement excavation Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000003825 pressing Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000003491 array Methods 0.000 description 2
- 238000003776 cleavage reaction Methods 0.000 description 2
- 230000001808 coupling effect Effects 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 241001391944 Commicarpus scandens Species 0.000 description 1
- 208000002565 Open Fractures Diseases 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000000802 evaporation-induced self-assembly Methods 0.000 description 1
- 239000003673 groundwater Substances 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 239000002352 surface water Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/08—Construction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
- G06Q50/265—Personal security, identity or safety
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Tourism & Hospitality (AREA)
- Computer Hardware Design (AREA)
- Health & Medical Sciences (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Evolutionary Computation (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Physics (AREA)
- General Business, Economics & Management (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Data Mining & Analysis (AREA)
- Civil Engineering (AREA)
- Algebra (AREA)
- Architecture (AREA)
- Structural Engineering (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Computer Security & Cryptography (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The application discloses a minimum safe thickness calculation method and device for penetrating through a movable fracture waterproof rock mass, wherein the method comprises the following steps: the fractal fracture trace of the target water-proof rock mass is extracted, the fracture fractal dimension of the target water-proof rock mass is calculated, so that static stress intensity factors of the crack initiation and expansion of the type I and type II fractal fractures are corrected according to the fractal fracture tip stress field, the magnitude relation between the static stress intensity factors and the fracture toughness of the rock is combined, the compound crack initiation and expansion criterion of the type I-type II fractal fracture tensile shear and the compound crack initiation and expansion criterion of the type II fractal fracture compressive shear are deduced, the lateral pressure coefficient expression of the position of the minimum safe thickness of the water-proof rock mass against the fracture is deduced based on the relation between the fracture thickness and the lateral pressure coefficient, and the minimum safe thickness of the tunnel crossing active fracture water-proof rock mass is calculated. Therefore, the problem that the true form of the rock mass fracture cannot be completely reflected by assuming that the rock mass fracture is straight and smooth in the related technology is solved, and the accuracy of the rock mass fracture is reduced.
Description
Technical Field
The application relates to the technical field of tunnels, in particular to a method and a device for calculating minimum safe thickness of a through-movable fracture waterproof rock mass.
Background
Because of large creep rate and high earthquake frequency, the rock-soil body cracks develop, loosen and break in the active fracture, which is beneficial to the enrichment of groundwater to form a water-containing structure, or the water-containing structure is connected with an underground aquifer and an earth surface water body as a water guide channel, when the tunnel excavation approaches the active fracture, the water-proof rock fracture is extremely easy to occur under the coupling effect of earthquake and hydrodynamic force to induce water-bursting and mud-bursting disasters.
In the related technology, as the water-proof rock body is extremely easy to break under the coupling action of earthquake and hydrodynamic force to induce water-bursting and mud-bursting disasters, the water-proof rock body fracture is assumed to be a straight smooth fracture, and the minimum safety thickness of the tunnel water-proof rock body is calculated so as to reduce the threat to the life safety of tunnel staff.
However, in the related art, the marine rock mass fracture is assumed to be a straight and smooth fracture, the real form of the rock mass fracture cannot be completely reflected, the accuracy of the marine rock mass fracture is reduced, the potential safety hazard is increased, and the actual requirements of tunnel staff cannot be met, so that the problem is to be solved.
Disclosure of Invention
The application provides a minimum safe thickness calculation method and device for a through-movable fracture waterproof rock body, which are used for solving the problems that in the related art, a waterproof rock body crack is assumed to be a straight smooth crack, the real form of the rock body crack cannot be completely reflected, the accuracy of the waterproof rock body crack is reduced, the potential safety hazard is increased, and the actual demands of tunnel staff cannot be met.
Embodiments of a first aspect of the present application provide a minimum safe thickness calculation across a mobile fractured water-resistant rock mass comprising the steps of: extracting a fractal fracture trace of a target water-proof rock mass, and calculating a fracture fractal dimension of the target water-proof rock mass according to the fractal fracture trace; correcting static stress intensity factors of crack initiation and expansion of type I and type II fractal cracks according to a fractal crack tip stress field based on the crack fractal dimension; based on critical conditions of rock mass fracture initiation and expansion in fracture mechanics, deducing a compound initiation and expansion criterion of the type I-II fractal fracture tensile shear and a compound initiation and expansion criterion of the type II fractal fracture compressive shear by combining the magnitude relation of the static stress intensity factor and the rock fracture toughness; based on a relation between the anti-splitting thickness and the lateral pressure coefficient, deducing a lateral pressure coefficient expression of the anti-splitting minimum safe thickness position of the waterproof rock body by comparing the relation between the critical water pressure and the movable fracture water pressure; and calculating the minimum safe thickness of the tunnel penetrating through the movable fracture waterproof rock according to the lateral pressure coefficient expression.
Optionally, in one embodiment of the present application, the calculating the fracture fractal dimension of the target marine rock mass from the fractal fracture trace includes: covering the fractal crack trace for multiple times by adopting grids with different side lengths, and counting the number of grids occupied by the crack trace under each side length scale; and performing linear fitting according to the number of the square grids to obtain a linear slope of the correlation, and obtaining the fractal dimension of the fracture according to the linear slope.
Optionally, in one embodiment of the present application, the correcting the static stress intensity factor of the type I and type II fractal fracture initiation and propagation according to the fractal fracture tip stress field includes: characterizing the fractal fracture tip stress field through fracture average stress intensity factors, fractal dimension related parameters, fracture tip radial polar coordinates, annular coordinates and angular distribution functions; and obtaining a static stress intensity factor of the cracking expansion of the type I fractal crack and a static stress intensity factor of the cracking expansion of the type II fractal crack based on the fractal crack tip stress field and a preset stress function.
Optionally, in one embodiment of the present application, the type i-ii fractal fracture pull-shear composite fracture expansion criterion is:
Wherein, the liquid crystal display device comprises a liquid crystal display device,fractal dimension for fissure->Is the maximum principal stress +.>Is the minimum principal stress->Is the long axis of the fracture andincluded angle (I)>Is the pressure of crevice water>For seismic shear stress->Is seismic wave circle frequency +.>Fracture mean stress intensity factor->Rock fracture toughness for type I fractal fracture, < +.>Is->,/>For the moment of->1/2 of the parting line length, < >>And->Dividing the dynamic stress intensity factors by their valuesωCorresponding static value when=0, +.>And->Is the phase angle;
the II type fractal fracture compression shear cracking expansion criterion is as follows:
wherein, the liquid crystal display device comprises a liquid crystal display device,rock fracture toughness for type II fractal fracture, +.>Is the friction angle in the fracture surface.
Optionally, in one embodiment of the present application, the calculating the minimum safe thickness of the tunnel through the active fracture rock-insulating body according to the lateral pressure coefficient expression includes: substituting parameters of a tunnel section radius, a tunnel burial depth, a tunnel ground stress, a rock mass gravity, a rock mass original side pressure coefficient, a rock mass fracture toughness, a fracture zone thickness of the waterproof rock mass, a movable fracture water pressure, a fracture fractal dimension of the waterproof rock mass, a fracture dip angle and a SV seismic wave dynamic stress intensity factor, and calculating a specific numerical value of an anti-fracture thickness region; and adding the specific value to the fracture zone thickness to obtain the minimum safe thickness of the tunnel penetrating movable fracture water-proof rock mass.
Embodiments of a second aspect of the present application provide a minimum safe thickness calculation apparatus for traversing a mobile fractured water-resistant rock mass, comprising: the extraction module is used for extracting a fractal fracture trace of the target water-proof rock mass and calculating the fracture fractal dimension of the target water-proof rock mass according to the fractal fracture trace; the correction module is used for correcting static stress intensity factors of crack initiation and expansion of the type I and type II fractal cracks according to the fractal crack tip stress field based on the crack fractal dimension; the first deduction module is used for deducting a compound cracking expansion criterion of the type I-II fractal fracture tensile shear and a cracking expansion criterion of the type II fractal fracture compressive shear according to the critical condition of the cracking expansion of the rock mass fracture in fracture mechanics and combining the magnitude relation of the static stress intensity factor and the rock fracture toughness; the second deducing module is used for deducing a lateral pressure coefficient expression of the minimum anti-splitting safe thickness position of the waterproof rock body by comparing the relation between critical water pressure and movable fracture water pressure based on the relation between the anti-splitting thickness and the lateral pressure coefficient; and the calculation module is used for calculating the minimum safe thickness of the tunnel penetrating through the movable fracture waterproof rock body according to the lateral pressure coefficient expression.
Optionally, in one embodiment of the present application, the extracting module includes: the statistics unit is used for covering the fractal crack trace for multiple times by adopting grids with different side lengths, and counting the number of grids occupied by the crack trace under each side length scale; the first acquisition unit is used for carrying out linear fitting according to the number of the square grids to obtain a linear slope of a correlation, and obtaining the fracture fractal dimension according to the linear slope.
Optionally, in one embodiment of the present application, the correction module includes: the characterization unit is used for characterizing the fractal fracture tip stress field through fracture average stress intensity factors, fractal dimension related parameters, fracture tip radial polar coordinates, annular coordinates and angular distribution functions; the second acquisition unit is used for acquiring a static stress intensity factor of the initiation and expansion of the type I fractal fracture and a static stress intensity factor of the expansion of the type II fractal fracture based on the fractal fracture tip stress field and a preset stress function.
Optionally, in one embodiment of the present application, the type i-ii fractal fracture pull-shear composite fracture expansion criterion is:
wherein, the liquid crystal display device comprises a liquid crystal display device,fractal dimension for fissure->Is the maximum principal stress +.>Is the minimum principal stress- >Is the long axis of the fracture andincluded angle (I)>Is the pressure of crevice water>For seismic shear stress->Is seismic wave circle frequency +.>Fracture mean stress intensity factor->Rock fracture toughness for type I fractal fracture, < +.>Is->,/>For the moment of->1/2 of the parting line length, < >>And->Dividing the dynamic stress intensity factors by their valuesωCorresponding static value when=0, +.>And->Is the phase angle;
the II type fractal fracture compression shear cracking expansion criterion is as follows:
wherein, the liquid crystal display device comprises a liquid crystal display device,rock fracture toughness for type II fractal fracture, +.>Is the friction angle in the fracture surface.
Optionally, in one embodiment of the present application, the computing module includes: the calculation unit is used for substituting parameters of a tunnel section radius, a tunnel burial depth, a tunnel ground stress, a rock weight, a rock original side pressure coefficient, a rock fracture toughness, a fracture zone thickness of the waterproof rock, a movable fracture water pressure, a fracture fractal dimension of the waterproof rock, a fracture dip angle and an SV seismic wave dynamic stress intensity factor to calculate a specific numerical value of an anti-fracture thickness zone; and the third acquisition unit is used for adding the specific value and the fracture zone thickness to obtain the minimum safe thickness of the tunnel-crossing movable fracture waterproof rock mass.
An embodiment of a third aspect of the present application provides an electronic device, including: the system comprises a memory, a processor and a computer program stored on the memory and capable of running on the processor, wherein the processor executes the program to realize the minimum safe thickness calculation method for crossing the active fracture waterproof rock mass.
A fourth aspect of the present application provides a computer readable storage medium storing a computer program which when executed by a processor implements a method of calculating a minimum safe thickness across an active fracture riser rock as above.
According to the method and the device, the fractal fracture trace of the target water-proof rock mass can be extracted, the fracture fractal dimension of the target water-proof rock mass is calculated, so that static stress intensity factors of the initiation and expansion of the type I and type II fractal fractures are corrected according to the fractal fracture tip stress field, the magnitude relation between the static stress intensity factors and the fracture toughness of the rock is combined, the I-type II fractal fracture tensile shear compound initiation and expansion criterion and the II type fractal fracture compression shear initiation and expansion criterion are deduced, the lateral pressure coefficient expression of the position of the minimum anti-fracture safety thickness of the water-proof rock mass is deduced based on the relation between the anti-fracture thickness and the lateral pressure coefficient, the minimum safety thickness of the tunnel crossing movable fracture water-proof rock mass is calculated, the real form of the rock mass fracture can be completely reflected, the accuracy of the water-proof rock mass fracture is improved, the potential safety hazard is reduced, and the actual demands of tunnel workers are met. Therefore, the problems that in the related technology, the rock mass fracture is assumed to be straight and smooth, the real form of the rock mass fracture cannot be completely reflected, the accuracy of the rock mass fracture is reduced, the potential safety hazard is increased, and the actual demands of tunnel staff cannot be met are solved.
Additional aspects and advantages of the application will be set forth in part in the description which follows and, in part, will be obvious from the description, or may be learned by practice of the application.
Drawings
The foregoing and/or additional aspects and advantages of the present application will become apparent and readily appreciated from the following description of the embodiments, taken in conjunction with the accompanying drawings, in which:
FIG. 1 is a flow chart of a method of minimum safe thickness calculation across a mobile fractured marine rock mass according to an embodiment of the present application;
FIG. 2 is a schematic structural view of a minimum safe thickness calculation device for traversing a mobile fractured marine rock mass according to an embodiment of the present application;
fig. 3 is a schematic structural diagram of an electronic device according to an embodiment of the present application.
Detailed Description
Embodiments of the present application are described in detail below, examples of which are illustrated in the accompanying drawings, wherein the same or similar reference numerals refer to the same or similar elements or elements having the same or similar functions throughout. The embodiments described below by referring to the drawings are exemplary and intended for the purpose of explaining the present application and are not to be construed as limiting the present application.
The method and the device for calculating the minimum safe thickness of the crossing movable fracture waterproof rock body according to the embodiment of the application are described below with reference to the accompanying drawings. Aiming at the problems that in the related technology mentioned in the background technology center, a water-proof rock crack is assumed to be a straight smooth crack, the real form of the rock crack cannot be completely reflected, the accuracy of the water-proof rock crack is reduced, the potential safety hazard is increased, and the actual demand of tunnel workers cannot be met, the application provides a minimum safe thickness calculation method for traversing the movable fracture water-proof rock, in the method, a fractal crack trace of the target water-proof rock can be extracted, the crack fractal dimension of the target water-proof rock is calculated, and therefore the static stress intensity factors of the crack initiation expansion of the type I and type II fractal cracks are corrected according to the fractal crack tip stress field, the static stress intensity factors and the size relation of the rock fracture toughness are combined, the type I-II fractal crack pulling-shearing compound initiation expansion criterion and the type II fractal crack pressing initiation expansion criterion are deduced, and the lateral pressure coefficient expression of the position of the minimum safe thickness of the movable fracture water-proof rock is deduced based on the relation of the anti-splitting thickness and the lateral pressure coefficient, and the minimum safe thickness of the movable fracture water-proof rock is calculated. Therefore, the problems that in the related technology, the rock mass fracture is assumed to be straight and smooth, the real form of the rock mass fracture cannot be completely reflected, the accuracy of the rock mass fracture is reduced, the potential safety hazard is increased, and the actual demands of tunnel staff cannot be met are solved.
Specifically, fig. 1 is a schematic flow chart of a method for calculating a minimum safe thickness of a through-movable fracture water-proof rock according to an embodiment of the present application.
As shown in fig. 1, the method for calculating the minimum safe thickness of the crossing movable fracture waterproof rock mass comprises the following steps:
in step S101, a fractal fracture trace of the target marine riser body is extracted, and a fracture fractal dimension of the target marine riser body is calculated from the fractal fracture trace.
It can be understood that the fractal fracture trace of the target water-proof rock mass can be extracted, and the fracture fractal dimension of the target water-proof rock mass can be calculated according to the fractal fracture trace in the following steps, so that the executable performance of the minimum safe thickness of the water-proof rock mass is effectively improved.
Wherein in one embodiment of the present application, calculating a fracture fractal dimension of the target marine riser body from the fractal fracture trace comprises: covering fractal crack traces for multiple times by using grids with different side lengths, and counting the number of grids occupied by the crack traces under each side length scale; and performing linear fitting according to the number of the square grids to obtain a linear slope of the correlation, and obtaining the fractal dimension of the fracture according to the linear slope.
In the actual implementation process, the embodiment of the application can use a square grid coverage method to adopt different side lengths LThe grids cover the crack trace repeatedly, and the number of grids occupied by the crack trace under each side length scale is countedNFor a plurality of ln #N) And ln%L) Linear fitting of the data points can obtain a linear slope of the correlation, wherein the slope counter-value is the fractal dimension of the fractureDThereby effectively improving the feasibility of the minimum safe thickness of the waterproof rock mass.
In step S102, based on the fracture fractal dimension, correcting the static stress intensity factors of the fracture initiation and expansion of the type I and type II fractal fractures according to the fractal fracture tip stress field.
It can be appreciated that the embodiment of the application can correct the static stress intensity factors of the type I and type II fractal fracture initiation and expansion according to the fractal fracture tip stress field in the following steps based on the fracture fractal dimension, for example, the fractal fracture tip stress field distribution expression is obtained based on fracture mechanics theory, and the static stress intensity factors of the type I and type II fractal fracture initiation and expansion are corrected through Westergaard stress function based on the fracture fractal dimension, wherein the dynamic stress intensity factors of the type I and type II fracture of the water-proof rock under the action of SV seismic waves can be obtained according to the active fracture dynamics theory, so that the reality of the rock fracture morphology is effectively increased.
Wherein in one embodiment of the present application, modifying the static stress intensity factors of the type I and type II fractal fracture initiation and propagation according to the fractal fracture tip stress field comprises: characterizing a fractal fracture tip stress field by fracture average stress intensity factors, fractal dimension related parameters, fracture tip radial polar coordinates, circumferential coordinates and an angular distribution function; and obtaining a static stress intensity factor of the initiation and expansion of the type I fractal fracture and a static stress intensity factor of the expansion of the type II fractal fracture based on the fractal fracture tip stress field and a preset stress function.
As one possible implementation, the fractal fracture tip stress field in embodiments of the present applicationσ ij Can pass through fracture average stress intensity factorK I Related parameters of fractal dimensionα,α=(2-D) 2) radial polar coordinates of fracture tiprCircumferential coordinatesθAnd an angular distribution functionf ij (θ) Characterization, namely:
σ ij =K I r -α f ij (θ) ,
wherein, the liquid crystal display device comprises a liquid crystal display device,σ ij is the fractal fracture tip stress field,K I is the average stress intensity factor of the fracture,ris the radial polar coordinate of the fracture tip,θin the form of a circular coordinate, the circular coordinate,f ij (θ) As a function of the angular distribution,αas the fractal dimension related parameters,α=(2-D)/2。
wherein the Westergaard stress function of the type I fractal fracture can pass the fracture surface load at the coordinate xP(x) Complex variables of stress function zAnd 1/2 of the apparent length of the fractal fracture) Characterization, namely:
wherein, the liquid crystal display device comprises a liquid crystal display device,P(x) For the fracture surface load,zas a complex variable of the stress function,is 1/2 of the apparent length of the fractal fracture.
The expression of the static stress intensity factor of the cracking expansion of the type I fractal crack is as follows:
in addition, the upper formula is in the interval [ sleeve [ ],/>]Upper integration, the expression:
wherein, the liquid crystal display device comprises a liquid crystal display device,is a dimensionless coordinate parameter.
Wherein, the liquid crystal display device comprises a liquid crystal display device,s=x/ ,xand is on the abscissa.
The correction method of the II type fractal crack expansion static stress intensity factor is the same, namely:
wherein, the liquid crystal display device comprises a liquid crystal display device,is shear stress.
In addition, the fracture tip dynamic stress intensity factor under the action of SV seismic waves is as follows:
wherein, the liquid crystal display device comprises a liquid crystal display device,shear stress of earthquakeK (2) I is the dynamic stress intensity factor divided by the stress intensity factor in the followingωThe corresponding static value when=0,ωfor the circular frequency of the seismic wave,δ 2 is the phase angle.
Wherein, the liquid crystal display device comprises a liquid crystal display device,the expression is:
wherein, the liquid crystal display device comprises a liquid crystal display device,νin order for the lame constant to be set,α 2 the wave number is the wave number of the SV seismic wave,polarization amplitude for seismic waves.
Wherein, the liquid crystal display device comprises a liquid crystal display device,α 2 =ω/c 2 ,c 2 is the wave velocity of the elastic transverse wave of the rock,=/> c 2 /ω 3 ,/>is the seismic peak acceleration.
In step S103, based on the critical conditions of the crack initiation and expansion of the rock mass in fracture mechanics, the magnitude relation of static stress intensity factors and the rock fracture toughness is combined, and the I-II type fractal crack pulling and shearing compound initiation and expansion criterion and the II type fractal crack pressing and shearing initiation and expansion criterion are deduced.
It can be understood that according to the embodiment of the application, based on the critical conditions of initiation and expansion of rock mass cracks in fracture mechanics in the following steps, the magnitude relation of static stress intensity factors and rock fracture toughness is combined, and the I-II type fractal crack pulling and shearing composite initiation and expansion criterion and the II type fractal crack pressing and shearing initiation and expansion criterion in the following steps are deduced, so that initiation, expansion and penetration processes of the rock mass cracks can be completely and truly reflected, and the accuracy of the waterproof rock mass cracks is improved.
In some embodiments, when the static stress intensity factor is greater than or equal to the rock fracture toughness, the fracture of the marine rock will initiate crack propagation, i.e.:
wherein, the liquid crystal display device comprises a liquid crystal display device,and->The stress intensity factors of the fractal fracture compression shear expansion are respectively I type and II type, and the stress intensity factors are +.>And->Rock fracture toughness of type i and type ii fractal fractures, respectively.
At this time, the fracture tip surface stress state expression is:
wherein, the liquid crystal display device comprises a liquid crystal display device,normal stress of fracture surface->Is the shear stress of fracture surface->Is the maximum principal stress +.>Is the minimum principal stress->Is the long axis and->Included angle (I)>Is the pressure of the fracture water.
Secondly, when the normal fracture stress is tensile stress, the compound expansion mode of the I-II type fractal fracture tensile shear is shown, and at the moment, the static stress intensity factor of the I type fractal fracture expansion is as follows:
Wherein, the liquid crystal display device comprises a liquid crystal display device,1/2 of the fracture length.
The type II fractal crack expansion static stress intensity factor is:
the compound cracking extension criterion of the I-II type fractal fracture tensile shear is as follows:
wherein, the liquid crystal display device comprises a liquid crystal display device,fractal dimension for fissure->Is the maximum principal stress +.>Is the minimum principal stress->Is the long axis of the fracture andincluded angle (I)>Is the pressure of crevice water>For seismic shear stress->Is seismic wave circle frequency +.>Fracture mean stress intensity factor->Rock fracture toughness for type I fractal fracture, < +.>Is->,/>For the moment of->1/2 of the parting line length, < >>And->Dividing the dynamic stress intensity factors by their valuesωCorresponding static value when=0, +.>And->Is the phase angle.
In addition, when the normal fracture stress is compressive stress, the fracture surface is in a type II fractal fracture compression shear expansion mode, the fracture surface is effectively shear stress due to fracture expansion caused by shear stressThe method comprises the following steps:
the II-type fractal fracture compression shear expansion stress intensity factor is as follows:
deducing a compound cracking extension criterion of the I-II type fractal fracture tension shear, namely:
wherein, the liquid crystal display device comprises a liquid crystal display device,rock fracture toughness for type II fractal fracture, +.>Is the friction angle in the fracture surface.
In step S104, based on the relation between the anti-fracture thickness and the lateral pressure coefficient, the lateral pressure coefficient expression of the anti-fracture minimum safe thickness position of the waterproof rock mass is deduced by comparing the magnitude relation between the critical water pressure and the movable fracture water pressure.
It can be understood that the waterproof rock mass in the embodiment of the application may include an anti-fracture thickness zone and a fracture zone, the rock mass lateral pressure coefficient of the anti-fracture thickness zone is affected by excavation unloading disturbance, a certain change rule is shown along with the change of the distance between the rock mass and the face of the anti-fracture thickness zone, the fracture initiation and expansion critical water pressure of the waterproof rock mass fracture is reduced along with the decrease of the lateral pressure coefficient, if the fracture initiation and expansion critical water pressure of the boundary position of the anti-fracture thickness zone and the fracture zone is equal to the movable fracture water pressure, the anti-fracture thickness zone reaches the minimum safe thickness, and accordingly the relation between the anti-fracture thickness and the lateral pressure coefficient is deduced and established.
Further, the lateral pressure coefficient of the rock mass can change to cause the opposite directions of the maximum main stress and the minimum main stress, so that the crack expansion critical water pressure is changed, the calculation formula of the lateral pressure coefficient at the front side of the tunnel face is represented by the radius of the tunnel section and the original lateral pressure coefficient, meanwhile, the critical water pressure of the crack initiation expansion of the rock mass at the position can be obtained, and the lateral pressure coefficient expression of the fracture-resistant minimum safe thickness position of the waterproof rock mass can be deduced by comparing the relation between the critical water pressure and the movable fracture water pressure, so that the accuracy of calculating the minimum safe thickness of the waterproof rock mass is effectively improved, and the safety of tunnel workers is improved.
In practice, the marine rock mass of the embodiments of the present application includes a fracture-resistant thickness zoneS r And fracture zoneS f Is influenced by the disturbance of excavation unloading and is resistant to the lateral pressure coefficient of rock mass in a split thickness zoneA certain change rule is presented along with the change of the distance between the device and the face, namely:
wherein, the liquid crystal display device comprises a liquid crystal display device,is the original side pressure coefficient of rock mass, < > and->Is the distance between a certain position in the anti-cleavage thickness region and the face surface>Is the radius of the tunnel section.
The compound expansion critical water pressure of the I-II type fracture tension shear is as follows:
the expansion critical water pressure of the II type fractal fracture compression shear is as follows:
when the fracture crack initiation expansion critical water pressure and the movable fracture water pressure are generated at the junction position of the anti-fracture thickness region and the fracture zone regionpEqual, the anti-cleavage thickness region reaches the minimum safe thicknessThe method comprises the following steps:
wherein, the liquid crystal display device comprises a liquid crystal display device,is->Side pressure coefficient of the position.
Normally the vertical dead weight stress of the tunnel is the minimum principal stressσ 3 The method comprises the following steps:
σ 3 =γH,
wherein, the liquid crystal display device comprises a liquid crystal display device,γis the rock weight, and H is the tunnel burial depth.
The horizontal maximum principal stress is:
σ 1 =λγH,
when the original lateral pressure coefficient is larger than 1 and the face front lateral pressure coefficient is smaller than 1, two main stresses are in opposite directionsλThe position of the section where 1 is:
S 1 the I-II fractal fracture tensile shear compound expansion critical water pressure at the section position is as follows:
The expansion critical water pressure of the II type fractal fracture compression shear is as follows:
when (when)p>In the time-course of which the first and second contact surfaces,S r ´>S 1 ,σ 1 =λγH,σ 3 =γh, I-II type fractal fracture tension shear compound expansion modeλ 1 The method comprises the following steps:
wherein, the liquid crystal display device comprises a liquid crystal display device,pfor the pressure of the water of the movable fracture,and->The dynamic stress intensity factors of the tip of the type I and type II fissures are respectively shown.
Expansion mode of II-type fractal fracture compression shearλ 1 The method comprises the following steps:
when (when)p<In the time-course of which the first and second contact surfaces,S r ´<S 1 ,σ 1 and (3) withσ 3 In the different direction, the I-II fractal fracture tensile shear compound expansion modeλ 1 The method comprises the following steps:
expansion mode of II-type fractal fracture compression shearλ 1 The method comprises the following steps:
in step S105, the minimum safe thickness of the tunnel crossing active fracture riser is calculated from the lateral pressure coefficient expression.
It can be understood that the embodiment of the application can calculate the minimum safe thickness of the tunnel crossing movable fracture water-proof rock mass according to the lateral pressure coefficient expression in the following steps, so that the life safety of tunnel workers is effectively improved, the economic loss is reduced, the construction period delay is reduced, and an important theoretical basis is provided for the safety prevention and control of the water and mud burst of the tunnel crossing movable fracture.
Optionally, in one embodiment of the present application, calculating the minimum safe thickness of the tunnel through the active fracture riser from the lateral pressure coefficient expression includes: substituting parameters of a tunnel section radius, a tunnel burial depth, a tunnel ground stress, a rock mass gravity, a rock mass original side pressure coefficient, a rock mass fracture toughness, a fracture zone thickness of the waterproof rock mass, a movable fracture water pressure, a fracture fractal dimension of the waterproof rock mass, a fracture dip angle and a SV seismic wave dynamic stress intensity factor, and calculating a specific numerical value of an anti-fracture thickness region; and adding the specific numerical value and the fracture zone thickness to obtain the minimum safe thickness of the tunnel penetrating movable fracture water-proof rock mass.
In some embodiments, tunnel section radius, tunnel burial depth, tunnel ground stress, rock mass gravity, rock mass original side pressure coefficient and active fracture water pressure parameters in the embodiments of the application are obtained through on-site actual measurement, fracture dip angle can be obtained through simulation results of tunnel crossing active fracture numerical values, a grid coverage method is adopted to calculate fracture fractal dimension of the waterproof rock mass, and the fracture zone thickness of the waterproof rock mass is generally 1-3 m.
If the fracture is in an I-II type tensile shear composite expansion mode, the rock mass fracture toughness calculation formula is as follows:
wherein, the liquid crystal display device comprises a liquid crystal display device,is rock uniaxial compressive strength->Is an empirical parameter of the Hoek-Brown intensity criteria.
Wherein, the liquid crystal display device comprises a liquid crystal display device,
wherein, the liquid crystal display device comprises a liquid crystal display device,is the tensile strength of the rock.
If the fracture is in a type II compression-shear composite expansion mode, the rock mass fracture toughness calculation formula is as follows:
wherein, the liquid crystal display device comprises a liquid crystal display device,expansion of rock fracture toughness for type II compression shear composite +.>Is the minimum horizontal principal stress.
wherein, the liquid crystal display device comprises a liquid crystal display device,the fracture toughness of the rock is expanded for the I-II type pull-shear composite.
The fracture tip dynamic stress intensity factor under the action of SV seismic waves is as follows:
wherein, the liquid crystal display device comprises a liquid crystal display device,and->Representing a stress intensity factor of dimension 1 by dynamic Intensity factor graph queries specific values, < >>For fractal dimension related parameters +.>Is the rock elastic transverse wave velocity.
In the embodiment of the application, the critical water pressure of the compound expansion of the I-II type fractal fracture tension shear and the expansion of the II type fractal fracture tension shear of the marine-proof rock mass is calculated respectively, the crack initiation expansion mode of the marine-proof rock mass fracture can be determined by comparing the critical water pressure and the critical water pressure, and then the critical water pressure is calculatedThe minimum safe thickness of the anti-splitting thickness area of the waterproof rock mass is finally obtained, the life safety of tunnel workers is effectively improved, the economic loss is reduced, the construction period delay is reduced, and the national important tunnel engineering construction is effectively ensured.
According to the method for calculating the minimum safe thickness of the crossing active fracture water-proof rock, the fractal fracture trace of the target water-proof rock can be extracted, the fracture fractal dimension of the target water-proof rock is calculated, and therefore the static stress intensity factors of the cracking expansion of the type I and type II fractal fractures are corrected according to the fractal fracture tip stress field, the magnitude relation between the static stress intensity factors and the fracture toughness of the rock is combined, the compound cracking expansion criterion of the type I-type II fractal fracture pull-shear and the compound cracking expansion criterion of the type II fractal fracture pull-shear are deduced, the lateral pressure coefficient expression of the position of the minimum safe thickness of the water-proof rock is deduced based on the relation between the anti-cracking thickness and the lateral pressure coefficient, the minimum safe thickness of the tunnel crossing active fracture water-proof rock is calculated, the true form of the rock fracture can be completely reflected, the accuracy of the fracture of the water-proof rock is improved, the potential safety hazard is reduced, and the actual requirements of tunnel workers are met. Therefore, the problems that in the related technology, the rock mass fracture is assumed to be straight and smooth, the real form of the rock mass fracture cannot be completely reflected, the accuracy of the rock mass fracture is reduced, the potential safety hazard is increased, and the actual demands of tunnel staff cannot be met are solved.
A minimum safe thickness calculation device for traversing a mobile fracture marine rock according to an embodiment of the present application will be described next with reference to the accompanying drawings.
Fig. 2 is a block schematic diagram of a minimum safe thickness calculation device traversing a mobile fractured marine rock mass according to an embodiment of the present application.
As shown in fig. 2, the minimum safe thickness calculating apparatus 10 for crossing a movable fracture marine rock mass includes: the system comprises an extraction module 100, a correction module 200, a first derivation module 300, a second derivation module 400 and a calculation module 500.
Specifically, the extraction module 100 is configured to extract a fractal fracture trace of the target marine rock mass, and calculate a fracture fractal dimension of the target marine rock mass according to the fractal fracture trace.
The correction module 200 is used for correcting the static stress intensity factors of the crack initiation and expansion of the type I and type II fractal cracks according to the fractal crack tip stress field based on the crack fractal dimension.
The first deriving module 300 is used for deriving a type I-II fractal fracture tensile shear compound fracture initiation and expansion criterion and a type II fractal fracture compression shear fracture initiation and expansion criterion based on critical conditions of rock fracture initiation and expansion in fracture mechanics and in combination with the magnitude relation of static stress intensity factors and rock fracture toughness.
The second deriving module 400 is configured to derive a lateral pressure coefficient expression of the fracture-resistant minimum safe thickness position of the waterproof rock body by comparing the magnitude relation between the critical water pressure and the movable fracture water pressure based on the relation between the fracture-resistant thickness and the lateral pressure coefficient.
The calculation module 500 is used for calculating the minimum safe thickness of the tunnel passing through the movable fracture waterproof rock mass according to the lateral pressure coefficient expression.
Optionally, in one embodiment of the present application, the extracting module includes: the device comprises a statistics unit and a first acquisition unit.
The statistics unit is used for adopting different side length grids to cover fractal crack traces for multiple times and counting the number of grids occupied by the crack traces under each side length scale.
The first acquisition unit is used for carrying out linear fitting according to the number of the square grids to obtain a linear slope of the correlation, and obtaining the fractal dimension of the fracture according to the linear slope.
Optionally, in one embodiment of the present application, the correction module includes: a characterization unit and a second acquisition unit.
The characterization unit characterizes the fractal fracture tip stress field through fracture average stress intensity factors, fractal dimension related parameters, fracture tip radial polar coordinates, circumferential coordinates and angular distribution functions.
The second acquisition unit is used for acquiring a static stress intensity factor of the initiation and expansion of the type I fractal fracture and a static stress intensity factor of the expansion of the type II fractal fracture based on the fractal fracture tip stress field and a preset stress function.
Optionally, in one embodiment of the present application, the type i-ii fractal fracture pull shear composite fracture expansion criteria is:
wherein, the liquid crystal display device comprises a liquid crystal display device,fractal dimension for fissure->Is the maximum principal stress +.>Is the minimum principal stress->Is the long axis of the fracture andincluded angle (I)>Is the pressure of crevice water>For the shear stress of an earthquake,/>is seismic wave circle frequency +.>Fracture mean stress intensity factor->Rock fracture toughness for type I fractal fracture, < +.>Is->,/>For the moment of->1/2 of the parting line length, < >>And->Dividing the dynamic stress intensity factors by their valuesωCorresponding static value when=0, +.>And->Is the phase angle.
The cracking extension criterion of the II type fractal fracture compression shear is as follows:
wherein, the liquid crystal display device comprises a liquid crystal display device,rock fracture toughness for type II fractal fracture, +.>Is the friction angle in the fracture surface.
Optionally, in one embodiment of the present application, the computing module includes: a calculation unit and a third acquisition unit.
The calculation unit is used for substituting parameters of a tunnel section radius, a tunnel burial depth, a tunnel ground stress, a rock mass weight, a rock mass original side pressure coefficient, a rock mass fracture toughness, a water-proof rock mass fracture zone thickness, a movable fracture water pressure, a water-proof rock mass fracture fractal dimension, a fracture dip angle and an SV seismic wave dynamic stress intensity factor to calculate specific values of an anti-fracture thickness region.
And the third acquisition unit is used for adding the specific numerical value and the fracture zone thickness to obtain the minimum safe thickness of the tunnel crossing movable fracture water-proof rock mass.
It should be noted that the foregoing explanation of the embodiment of the method for calculating the minimum safe thickness of the penetrating movable fracture water-proof rock is also applicable to the device for calculating the minimum safe thickness of the penetrating movable fracture water-proof rock of the embodiment, and will not be repeated herein.
According to the minimum safe thickness calculation device for penetrating through the movable fracture waterproof rock body, the fractal fracture trace of the target waterproof rock body can be extracted, the fracture fractal dimension of the target waterproof rock body is calculated, and therefore the static stress intensity factors of the cracking expansion of the type I and type II fractal fractures are corrected according to the fractal fracture tip stress field, the magnitude relation between the static stress intensity factors and the fracture toughness of the rock is combined, the compound cracking expansion criterion of the type I-type II fractal fracture pull-shear and the compound cracking expansion criterion of the type II fractal fracture pull-shear are deduced, the lateral pressure coefficient expression of the minimum safe thickness position of the waterproof rock body is deduced based on the relation between the anti-cracking thickness and the lateral pressure coefficient, the minimum safe thickness of the tunnel penetrating through the movable fracture waterproof rock body is calculated, the true form of the rock fracture can be completely reflected, the accuracy of the fracture of the waterproof rock body is improved, the potential safety hazard is reduced, and the practical requirements of tunnel workers are met. Therefore, the problems that in the related technology, the rock mass fracture is assumed to be straight and smooth, the real form of the rock mass fracture cannot be completely reflected, the accuracy of the rock mass fracture is reduced, the potential safety hazard is increased, and the actual demands of tunnel staff cannot be met are solved.
Fig. 3 is a schematic structural diagram of an electronic device according to an embodiment of the present application. The electronic device may include:
The processor 302, when executing the program, implements the minimum safe thickness calculation method for traversing a mobile fractured marine rock mass provided in the above-described embodiments.
Further, the electronic device further includes:
a communication interface 303 for communication between the memory 301 and the processor 302.
A memory 301 for storing a computer program executable on the processor 302.
The memory 301 may comprise a high-speed RAM memory or may further comprise a non-volatile memory (non-volatile memory), such as at least one disk memory.
If the memory 301, the processor 302, and the communication interface 303 are implemented independently, the communication interface 303, the memory 301, and the processor 302 may be connected to each other through a bus and perform communication with each other. The bus may be an industry standard architecture (Industry Standard Architecture, abbreviated ISA) bus, an external device interconnect (Peripheral Component, abbreviated PCI) bus, or an extended industry standard architecture (Extended Industry Standard Architecture, abbreviated EISA) bus, among others. The buses may be divided into address buses, data buses, control buses, etc. For ease of illustration, only one thick line is shown in fig. 3, but not only one bus or one type of bus.
Alternatively, in a specific implementation, if the memory 301, the processor 302, and the communication interface 303 are integrated on a chip, the memory 301, the processor 302, and the communication interface 303 may communicate with each other through internal interfaces.
The processor 302 may be a central processing unit (Central Processing Unit, abbreviated as CPU), or an application specific integrated circuit (Application Specific Integrated Circuit, abbreviated as ASIC), or one or more integrated circuits configured to implement embodiments of the present application.
Embodiments of the present application also provide a computer readable storage medium having stored thereon a computer program which, when executed by a processor, implements a method of minimum safe thickness calculation for traversing an active fractured marine rock mass as described above.
In the description of the present specification, a description referring to terms "one embodiment," "some embodiments," "examples," "specific examples," or "some examples," etc., means that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the present application. In this specification, schematic representations of the above terms are not necessarily directed to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or N embodiments or examples. Furthermore, the different embodiments or examples described in this specification and the features of the different embodiments or examples may be combined and combined by those skilled in the art without contradiction.
Furthermore, the terms "first," "second," and the like, are used for descriptive purposes only and are not to be construed as indicating or implying a relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defining "a first" or "a second" may explicitly or implicitly include at least one such feature. In the description of the present application, the meaning of "N" is at least two, such as two, three, etc., unless explicitly defined otherwise.
Any process or method descriptions in flow charts or otherwise described herein may be understood as representing modules, segments, or portions of code which include one or more executable instructions for implementing specific logical functions or steps of the process, and additional implementations are included within the scope of the preferred embodiment of the present application in which functions may be executed out of order from that shown or discussed, including substantially concurrently or in reverse order from that shown or discussed, depending on the functionality involved, as would be understood by those reasonably skilled in the art of the embodiments of the present application.
Logic and/or steps represented in the flowcharts or otherwise described herein, e.g., a ordered listing of executable instructions for implementing logical functions, can be embodied in any computer-readable medium for use by or in connection with an instruction execution system, apparatus, or device, such as a computer-based system, processor-containing system, or other system that can fetch the instructions from the instruction execution system, apparatus, or device and execute the instructions. For the purposes of this description, a "computer-readable medium" can be any means that can contain, store, communicate, propagate, or transport the program for use by or in connection with the instruction execution system, apparatus, or device. More specific examples (a non-exhaustive list) of the computer-readable medium would include the following: an electrical connection (electronic device) having one or N wires, a portable computer cartridge (magnetic device), a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber device, and a portable compact disc read-only memory (CDROM). In addition, the computer readable medium may even be paper or other suitable medium on which the program is printed, as the program may be electronically captured, via optical scanning of the paper or other medium, then compiled, interpreted or otherwise processed in a suitable manner if necessary, and then stored in a computer memory.
It is to be understood that portions of the present application may be implemented in hardware, software, firmware, or a combination thereof. In the above-described embodiments, the N steps or methods may be implemented in software or firmware stored in a memory and executed by a suitable instruction execution system. If implemented in hardware, as in another embodiment, may be implemented using any one or combination of the following techniques, as known in the art: discrete logic circuits having logic gates for implementing logic functions on data signals, application specific integrated circuits having suitable combinational logic gates, programmable Gate Arrays (PGAs), field Programmable Gate Arrays (FPGAs), and the like.
Those of ordinary skill in the art will appreciate that all or a portion of the steps carried out in the method of the above-described embodiments may be implemented by a program to instruct related hardware, where the program may be stored in a computer readable storage medium, and where the program, when executed, includes one or a combination of the steps of the method embodiments.
In addition, each functional unit in each embodiment of the present application may be integrated in one processing module, or each unit may exist alone physically, or two or more units may be integrated in one module. The integrated modules may be implemented in hardware or in software functional modules. The integrated modules may also be stored in a computer readable storage medium if implemented in the form of software functional modules and sold or used as a stand-alone product.
The above-mentioned storage medium may be a read-only memory, a magnetic disk or an optical disk, or the like. Although embodiments of the present application have been shown and described above, it will be understood that the above embodiments are illustrative and not to be construed as limiting the application, and that variations, modifications, alternatives, and variations may be made to the above embodiments by one of ordinary skill in the art within the scope of the application.
Claims (10)
1. The method for calculating the minimum safe thickness of the crossing movable fracture waterproof rock mass is characterized by comprising the following steps of:
extracting a fractal fracture trace of a target water-proof rock mass, and calculating a fracture fractal dimension of the target water-proof rock mass according to the fractal fracture trace;
correcting static stress intensity factors of crack initiation and expansion of type I and type II fractal cracks according to a fractal crack tip stress field based on the crack fractal dimension;
based on critical conditions of rock mass fracture initiation and expansion in fracture mechanics, deducing a compound initiation and expansion criterion of the type I-II fractal fracture tensile shear and a compound initiation and expansion criterion of the type II fractal fracture compressive shear by combining the magnitude relation of the static stress intensity factor and the rock fracture toughness;
based on a relation between the anti-splitting thickness and the lateral pressure coefficient, deducing a lateral pressure coefficient expression of the anti-splitting minimum safe thickness position of the waterproof rock body by comparing the relation between the critical water pressure and the movable fracture water pressure; and
And calculating the minimum safe thickness of the tunnel penetrating through the movable fracture waterproof rock according to the lateral pressure coefficient expression.
2. The method of calculating a minimum safe thickness across an active fracture riser of claim 1, wherein calculating a fracture fractal dimension of the target riser from the fractal fracture trace comprises:
covering the fractal crack trace for multiple times by adopting grids with different side lengths, and counting the number of grids occupied by the crack trace under each side length scale;
and performing linear fitting according to the number of the square grids to obtain a linear slope of the correlation, and obtaining the fractal dimension of the fracture according to the linear slope.
3. The method of calculating a minimum safe thickness across a mobile fracture riser of claim 1, wherein correcting the static stress intensity factors of type I and type II fractal fracture initiation and propagation based on the fractal fracture tip stress field comprises:
characterizing the fractal fracture tip stress field through fracture average stress intensity factors, fractal dimension related parameters, fracture tip radial polar coordinates, annular coordinates and angular distribution functions;
and obtaining a static stress intensity factor of the cracking expansion of the type I fractal crack and a static stress intensity factor of the cracking expansion of the type II fractal crack based on the fractal crack tip stress field and a preset stress function.
4. The method of claim 1, wherein the type i-ii fractal fracture pull shear composite fracture expansion criteria is:
wherein, the liquid crystal display device comprises a liquid crystal display device,fractal dimension for fissure->Is the maximum principal stress +.>Is the minimum principal stress->Is the long axis and->Included angle (I)>Is the pressure of crevice water>For seismic shear stress->Is seismic wave circle frequency +.>Fracture mean stress intensity factor->Rock fracture toughness for type I fractal fracture, < +.>Is->,/>For the moment of->1/2 of the parting line length, < >>Anddivided by dynamic stress intensity factorωCorresponding static value when=0, +.>And->Is the phase angle;
the II type fractal fracture compression shear cracking expansion criterion is as follows:
5. The method for calculating the minimum safe thickness across a mobile fracture riser as recited in claim 1, wherein calculating the minimum safe thickness across a mobile fracture riser for a tunnel from the lateral pressure coefficient expression comprises:
substituting parameters of a tunnel section radius, a tunnel burial depth, a tunnel ground stress, a rock mass gravity, a rock mass original side pressure coefficient, a rock mass fracture toughness, a fracture zone thickness of the waterproof rock mass, a movable fracture water pressure, a fracture fractal dimension of the waterproof rock mass, a fracture dip angle and a SV seismic wave dynamic stress intensity factor, and calculating a specific numerical value of an anti-fracture thickness region;
And adding the specific value to the fracture zone thickness to obtain the minimum safe thickness of the tunnel penetrating movable fracture water-proof rock mass.
6. A minimum safe thickness calculation device for traversing a mobile fractured water-resistant rock mass, comprising:
the extraction module is used for extracting a fractal fracture trace of the target water-proof rock mass and calculating the fracture fractal dimension of the target water-proof rock mass according to the fractal fracture trace;
the correction module is used for correcting static stress intensity factors of crack initiation and expansion of the type I and type II fractal cracks according to the fractal crack tip stress field based on the crack fractal dimension;
the first deduction module is used for deducting a compound cracking expansion criterion of the type I-II fractal fracture tensile shear and a cracking expansion criterion of the type II fractal fracture compressive shear according to the critical condition of the cracking expansion of the rock mass fracture in fracture mechanics and combining the magnitude relation of the static stress intensity factor and the rock fracture toughness;
the second deducing module is used for deducing a lateral pressure coefficient expression of the minimum anti-splitting safe thickness position of the waterproof rock body by comparing the relation between critical water pressure and movable fracture water pressure based on the relation between the anti-splitting thickness and the lateral pressure coefficient; and
And the calculation module is used for calculating the minimum safe thickness of the tunnel penetrating through the movable fracture waterproof rock body according to the lateral pressure coefficient expression.
7. The minimum safe thickness computing device for traversing a mobile fracture riser as recited in claim 6, wherein the extraction module comprises:
the statistics unit is used for covering the fractal crack trace for multiple times by adopting grids with different side lengths, and counting the number of grids occupied by the crack trace under each side length scale;
the first acquisition unit is used for carrying out linear fitting according to the number of the square grids to obtain a linear slope of a correlation, and obtaining the fracture fractal dimension according to the linear slope.
8. The minimum safe thickness calculation device for traversing a mobile fracture riser as recited in claim 6, wherein the correction module comprises:
the characterization unit is used for characterizing the fractal fracture tip stress field through fracture average stress intensity factors, fractal dimension related parameters, fracture tip radial polar coordinates, annular coordinates and angular distribution functions;
the second acquisition unit is used for acquiring a static stress intensity factor of the initiation and expansion of the type I fractal fracture and a static stress intensity factor of the expansion of the type II fractal fracture based on the fractal fracture tip stress field and a preset stress function.
9. An electronic device, comprising: a memory, a processor and a computer program stored on the memory and executable on the processor, the processor executing the program to implement the method of minimum safe thickness calculation across an active fracture riser as claimed in any one of claims 1 to 5.
10. A computer readable storage medium having stored thereon a computer program, wherein the program is executed by a processor for implementing a method of minimum safe thickness calculation across an active fracture riser as claimed in any one of claims 1 to 5.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310667132.9A CN116432296A (en) | 2023-06-07 | 2023-06-07 | Method and device for calculating minimum safe thickness of crossing movable fracture waterproof rock mass |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310667132.9A CN116432296A (en) | 2023-06-07 | 2023-06-07 | Method and device for calculating minimum safe thickness of crossing movable fracture waterproof rock mass |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116432296A true CN116432296A (en) | 2023-07-14 |
Family
ID=87089334
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310667132.9A Pending CN116432296A (en) | 2023-06-07 | 2023-06-07 | Method and device for calculating minimum safe thickness of crossing movable fracture waterproof rock mass |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116432296A (en) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112131628A (en) * | 2020-07-29 | 2020-12-25 | 北京工业大学 | Method for calculating minimum water inrush prevention safety thickness by comprehensively considering inclination angle and trend of tunnel crossing water-rich fault |
CN113534291A (en) * | 2021-07-20 | 2021-10-22 | 中国石油大学(华东) | Quantitative prediction method for different-scale fractures of low-permeability reservoir under constraint of rock mechanical layer |
US20210334434A1 (en) * | 2020-04-24 | 2021-10-28 | Southwest Petroleum University | Optimization method for dense cutting, temporary plugging and fracturing in shale horizontal well stage |
-
2023
- 2023-06-07 CN CN202310667132.9A patent/CN116432296A/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20210334434A1 (en) * | 2020-04-24 | 2021-10-28 | Southwest Petroleum University | Optimization method for dense cutting, temporary plugging and fracturing in shale horizontal well stage |
CN112131628A (en) * | 2020-07-29 | 2020-12-25 | 北京工业大学 | Method for calculating minimum water inrush prevention safety thickness by comprehensively considering inclination angle and trend of tunnel crossing water-rich fault |
CN113534291A (en) * | 2021-07-20 | 2021-10-22 | 中国石油大学(华东) | Quantitative prediction method for different-scale fractures of low-permeability reservoir under constraint of rock mechanical layer |
Non-Patent Citations (1)
Title |
---|
孔凡猛: "川藏铁路深埋隧道穿越活动断裂动力突水突泥致灾机理研究", 中国优秀博士学位论文库, pages 102 - 116 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Li et al. | Assessment of underground tunnel stability to adjacent tunnel explosion | |
Trivino et al. | Assessment of crack initiation and propagation in rock from explosion-induced stress waves and gas expansion by cross-hole seismometry and FEM–DEM method | |
Li et al. | Influence of unloading disturbance on adjacent tunnels | |
CN106767515A (en) | A kind of tunnel defect quick diagnosis prevention and controls | |
Chen et al. | Blasting excavation induced damage of surrounding rock masses in deep-buried tunnels | |
CN105160138B (en) | Tunnel flow consolidates analysis of coupled system's stability method | |
Xiao et al. | Large deformation characteristics and reinforcement measures for a rock pillar in the Houziyan underground powerhouse | |
Wang et al. | Experimental study on failure process of arch dam based on acoustic emission technique | |
CN111553101A (en) | Method for forecasting cracking of overlying rock layer in tunnel excavation and surrounding rock supporting method | |
US11966672B2 (en) | Method and system for simulating contact and interaction between support member and chamber surrounding rock mass | |
Lin et al. | Development of a real-time scour monitoring system for bridge safety evaluation | |
CN110888161A (en) | Method and system for calculating b value in micro-seismic magnitude-frequency relation | |
Wu et al. | Reliability analysis and prediction on tunnel roof under blasting disturbance | |
Zareifard et al. | A simplified solution for stresses around lined pressure tunnels considering non-radial symmetrical seepage flow | |
Liu et al. | Nondestructive testing on cumulative damage of watery fractured rock mass under multiple cycle blasting | |
Shi et al. | Failure mechanism analysis for tunnel construction crossing the water-rich dense fracture zones: A case study | |
CN116432296A (en) | Method and device for calculating minimum safe thickness of crossing movable fracture waterproof rock mass | |
Zhou et al. | Safety monitoring and stability analysis of left bank high slope at Jinping-I hydropower station | |
Kositsyn et al. | Geotechnical projection of the influence of the construction of the designed metropolitene tunnel by the method of shield passage on the sedimentation of the earth’s surface | |
CN111379567B (en) | Reliability analysis method for initial support stage of expansive rock tunnel | |
CN112948944A (en) | Pipeline strain calculation method under action of general continuous surface displacement | |
Ishanxodjaev et al. | Influence of the inertia force of underground pipeline systems under seismic loads | |
Oreste | Evaluation of the tunnel face stability through a ground stress analysis with a hemispherical geometry approximation | |
Curran et al. | A two-dimensional approach for designing tunnel support in weak rock | |
Gao et al. | Limit analysis of tunnel collapse according to the Hoek–Brown criterion and bolt parameter research |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20230714 |