CN112818439A - Soft rock tunnel surrounding rock sub-grade grading method - Google Patents

Soft rock tunnel surrounding rock sub-grade grading method Download PDF

Info

Publication number
CN112818439A
CN112818439A CN202011632957.XA CN202011632957A CN112818439A CN 112818439 A CN112818439 A CN 112818439A CN 202011632957 A CN202011632957 A CN 202011632957A CN 112818439 A CN112818439 A CN 112818439A
Authority
CN
China
Prior art keywords
rock
grading
grade
sub
test
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011632957.XA
Other languages
Chinese (zh)
Other versions
CN112818439B (en
Inventor
王永东
万善通
彭浩
刘晓
燕新
梁辉如
刘俊锋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangdong Road And Bridge Construction Development Co ltd
Changan University
Original Assignee
Guangdong Road And Bridge Construction Development Co ltd
Changan University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Guangdong Road And Bridge Construction Development Co ltd, Changan University filed Critical Guangdong Road And Bridge Construction Development Co ltd
Priority to CN202011632957.XA priority Critical patent/CN112818439B/en
Publication of CN112818439A publication Critical patent/CN112818439A/en
Application granted granted Critical
Publication of CN112818439B publication Critical patent/CN112818439B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/08Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/24Investigating strength properties of solid materials by application of mechanical stress by applying steady shearing forces
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/24Earth materials
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0639Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06393Score-carding, benchmarking or key performance indicator [KPI] analysis
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0014Type of force applied
    • G01N2203/0016Tensile or compressive
    • G01N2203/0019Compressive
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0014Type of force applied
    • G01N2203/0025Shearing

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Immunology (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • Geometry (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Theoretical Computer Science (AREA)
  • Strategic Management (AREA)
  • Computer Hardware Design (AREA)
  • Development Economics (AREA)
  • Educational Administration (AREA)
  • Economics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Mathematical Analysis (AREA)
  • Structural Engineering (AREA)
  • Marketing (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Game Theory and Decision Science (AREA)
  • Civil Engineering (AREA)
  • Architecture (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Remote Sensing (AREA)

Abstract

A soft rock tunnel surrounding rock sub-grade grading method comprises the following steps: determining grading indexes and grading; calculating the numerical characteristic value of each grading index; MATLAB generates a corresponding cloud model diagram, and the test value certainty is read; acquiring four test parameter values of the sub-level section to be divided through an indoor test; determining each grading index weight; calculating comprehensive certainty and judging the grade; according to the invention, four parameters which do not pass through a rock sample are obtained through an indoor test, the soft rock tunnel surrounding rock sub-level classification is carried out by combining a targeted classification method, the IV-level surrounding rock and the V-level surrounding rock are divided into five sub-levels of IV-a, IV-b, IV-c, V-a and V-b, and corresponding support measures are formulated according to different sub-levels, so that the subdivided surrounding rock brings great convenience for subsequent tunnel construction, plays a role in saving a construction period and avoiding construction safety accidents, and can be popularized and applied to the fields of soft rock tunnel surrounding rock classification and tunnel support design.

Description

Soft rock tunnel surrounding rock sub-grade grading method
Technical Field
The invention belongs to the technical field of supporting equipment or devices for graded design of tunnel surrounding rocks, and particularly relates to a soft rock tunnel surrounding rock sub-grade grading method.
Background
The surrounding rock grading refers to selecting some indexes related to the tunnel surrounding rock, such as rock hardness degree, rock integrity degree, rock structure surface state and the like, and dividing the tunnel surrounding rock into a plurality of grades according to specific standards or methods according to the requirements of tunnel engineering design and construction. In the current design process of the highway tunnel, after the surrounding rock grade is divided, the excavation method, the supporting measure and the supporting parameter corresponding to the surrounding rock grade can be further determined according to the specification and the standard.
For IV-grade and V-grade soft surrounding rocks, the grading system has large span and is very easy to cause adverse effects on the aspects of construction progress, support measures, support parameters and the like of the soft rock tunnel, and the tunnel needs more detailed surrounding rock grade division because the soft rock has strong complexity and uncertainty.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provide the soft rock tunnel surrounding rock sub-grade classification method which is reasonable in design, brings great convenience to subsequent tunnel construction, saves the construction period and avoids construction safety accidents.
The technical scheme for solving the technical problems comprises the following steps:
A. determining grading indexes and grading;
(A1) determining the value range of the uniaxial compressive strength Rc in IV and V-grade surrounding rock to be 0-60 MPa;
(A2) determining the value range of the rock integrity coefficient Kv at IV and V level surrounding rocks to be 0-0.75, and determining the value range of the cohesive force c at IV and V level surrounding rocks to be 0-0.7 MPa;
(A3) softening coefficient K of surrounding rockfThe value range of the IV-grade and V-grade surrounding rocks is determined to be 0-0.75;
B. calculating the numerical characteristic value of each grading index;
determining expected E of normal cloud with surrounding rock sub-grade grading index corresponding to certain grade standardxEntropy of normal cloud EnNormal cloud super entropy He
Figure BDA0002880515440000021
Figure BDA0002880515440000022
He=K
In the formula, CmaxAnd CminThe variable fuzzy threshold value is used for adjusting the variable fuzzy threshold value according to the constant K;
C. MATLAB generates a corresponding cloud model diagram, and the test value certainty is read;
(C1) producing the expected value is ExEntropy is EnThe normal random number x of (a); calculating the expected value is ExEntropy is EnIs normal membership cloud, the degree of membership at the desired curve x
Figure BDA0002880515440000023
Respectively calculating the entropy E of the changes along the expected curves MA and MB according to a reduced semi-normal rulenl、Enr:
Figure BDA0002880515440000024
(c2) Calculating points
Figure BDA0002880515440000025
Variance σ of (c)x
Figure BDA0002880515440000026
Figure BDA0002880515440000027
(C3) Generating the expected value is
Figure BDA0002880515440000028
Variance is σxNormal random number μ of (a):
Figure BDA0002880515440000029
(C4) finally obtaining the cloud droplet xi (x, mu);
(C5) after the cloud model digital characteristic values of all indexes under different sub-level grades are obtained, the cloud model graph corresponding to the tunnel surrounding rock sub-level grading indexes is obtained after the cloud model digital characteristic values of all surrounding rock sub-level grading indexes are input by using a forward cloud generator;
D. acquiring four test parameter values of the sub-level section to be divided through an indoor test;
(D1) adopting an indoor uniaxial compression test to obtain uniaxial compression strength and softening coefficient of the surrounding rock, and dividing the uniaxial compression test into a uniaxial compression test in a drying state, a uniaxial compression test in a water saturation state and a uniaxial compression test in a natural state according to different water containing states of the sample, wherein the ratio of the uniaxial compression strength in a rock saturation state to the uniaxial compression strength in a drying state is the softening coefficient of the rock;
(D2) after the longitudinal wave velocity of the rock mass and the elastic longitudinal wave velocity of the rock are respectively obtained through tests, the ratio of the longitudinal wave velocity to the elastic longitudinal wave velocity of the rock is the rock integrity coefficient;
(D3) obtaining corresponding shear strength indexes, namely cohesive force c and internal friction angle through direct shear test of rock
Figure BDA0002880515440000035
E. Determining each grading index weight; d, after obtaining a grading index test value of the required sub-grade grading section through an indoor test, determining the index weight through an entropy weight method by MATLAB software;
(E1) carrying out standardized processing on the data, and selecting j tunnel surrounding rock grading indexes, namely X1,X2,X3,···,XjIf i evaluation objects are set, the evaluation index value of the ith object is xi1,xi2,xi3,···,xijFormula y of normalization processijThe following were used:
Figure BDA0002880515440000031
(E2) according to the information entropy definition, calculating the information entropy corresponding to each index value of the n evaluation objects:
Figure BDA0002880515440000032
Figure BDA0002880515440000033
(E3) obtaining weight values omega of all indexes of i groups of surrounding rock grading objectsj
Figure BDA0002880515440000034
F. Calculating comprehensive certainty and judging the grade;
(F1) firstly, calculating the test value certainty, obtaining the grading index test value of the grading section of the soft rock tunnel according to an indoor test, reading the test value certainty mu (x) of each grading index test value corresponding to each grade according to a generated table by using MATLAB software, and marking the test value certainty mu (x) as 0 if the test value does not belong to a certain grade;
(F2) calculating comprehensive certainty factor, determining degree mu (x) by the calculated test value, and calculating each grading index weight omega by combining entropy weight methodjAnd calculating the comprehensive certainty U of each section:
Figure BDA0002880515440000041
(F3) and judging the sub-grade level of the surrounding rock of each section, and judging the final sub-grade level attribution of the surrounding rock of each section according to the obtained numerical value of the comprehensive certainty factor U of each section.
In step B of the present invention, K is 0.01[95 ].
According to the invention, four parameters including uniaxial compressive strength, rock integrity coefficient, softening coefficient and cohesive force which do not pass through a rock sample are obtained through an indoor test, the soft rock tunnel surrounding rock sub-grade classification is carried out by combining a targeted classification method, the IV-grade surrounding rock and the V-grade surrounding rock are divided into five sub-grades of IV-a, IV-b, IV-c, V-a and V-b, and corresponding support measures are formulated according to different sub-grades, so that the subdivided surrounding rock brings great convenience for subsequent tunnel construction, plays a role in saving construction period and avoiding construction safety accidents, and can be popularized and applied to the field of soft rock tunnel surrounding rock classification and tunnel support design.
Drawings
FIG. 1 is a flow chart of the present invention.
FIG. 2 is a schematic diagram of a cloud model of the present invention, wherein the origin is defined as point A.
FIG. 3 is a diagram of a uniaxial compressive strength cloud model of the present invention.
FIG. 4 is a diagram of an integrity coefficient cloud model of the present invention.
FIG. 5 is a cloud model of cohesion according to the present invention.
FIG. 6 is a diagram of a softening coefficient cloud model of the present invention.
Detailed Description
The present invention will be described in further detail with reference to the following drawings and examples, but the present invention is not limited to these examples.
Example 1
In fig. 1, the soft rock tunnel surrounding rock sub-grade grading method related by the invention comprises the following steps:
A. determining grading indexes and grading; through statistical arrangement of rock mass grading and tunnel surrounding rock grading index systems at home and abroad and combing of engineering grading and grade evaluation research results aiming at soft rocks or weak surrounding rocks, four parameters including uniaxial compressive strength, rock integrity coefficient, softening coefficient and cohesive force of the rocks are used as the grading index system for soft rock tunnel surrounding rock sub-grade grading, and the IV-grade surrounding rocks and the V-grade surrounding rocks are divided into five sub-grades of IV-a, IV-b, IV-c, V-a and V-b.
(A1) When surrounding rock stability analysis is carried out, the hardness degrees of IV-level surrounding rocks and V-level surrounding rocks mainly comprise three types of harder rocks, softer rocks and soft rocks, and in addition, extreme conditions, namely the conditions of completely weathered soft rocks and completely disintegrated soaked water are considered, so when index numerical level classification of uniaxial compressive strength is carried out, the value ranges of the uniaxial compressive strength Rc in the IV-level and V-level surrounding rocks are determined to be 0-60 MPa according to the relationship of qualitative classification of the uniaxial compressive strength and the hardness degrees of the rocks in the Highway tunnel design Specification (JTG 3370.1-2018);
(A2) determining the value range of the rock integrity coefficient Kv at IV and V level surrounding rocks to be 0-0.75, and determining the value range of the cohesive force c at IV and V level surrounding rocks to be 0-0.7 MPa;
(A3) the rock softening tendency can be defined as that the value is less than 0.75 by the specification of geotechnical engineering investigation Specification (GB 50021-2001), so that the surrounding rock softening coefficient K is definedfThe value range of the IV-grade and V-grade surrounding rocks is determined to be 0-0.75;
for the value ranges of the indexes among different surrounding rock grades, the related specifications and documents are mainly referred, and finally determined surrounding rock sub-grade grading index grades are classified and shown in table 1.
TABLE 1 carbonaceous mudstone tunnel surrounding rock grading index grading
Figure BDA0002880515440000051
B. Calculating the numerical characteristic value of each grading index;
determining expected E of normal cloud with surrounding rock sub-grade grading index corresponding to certain grade standardxEntropy of normal cloud EnNormal cloud super entropy He
Figure BDA0002880515440000061
Figure BDA0002880515440000062
He=K
In the formula, CmaxAnd CminThe variable fuzzy threshold value is used for adjusting the variable fuzzy threshold value according to the constant K; in the invention, K is uniformly 0.01[95]]. The numerical characteristic value calculation results of each grade index are shown in table 2.
TABLE 2 numerical characteristics of the respective grading indices
Figure BDA0002880515440000063
In the above steps, the cloud definition and the cloud model digital features are as follows:
and x, Y and C are set, wherein x belongs to Y, Y is a common quantitative set and is called a domain of discourse, and C represents a qualitative concept on Y. If any object x in Y has a random number mu (x) epsilon [0,1 ] with a stable tendency]It is called the membership of x to C, also calledAnd (5) determining the degree. The distribution of the certainty factor in Y is called a cloud model. If x satisfies: x to N (E)x,E'n) Wherein, E'n~N(En,H2 e) μ (x) and satisfies the following formula:
Figure BDA0002880515440000071
the distribution of the certainty μ on the domain of discourse Y is called a normal cloud, as shown in FIG. 1, in which a cubic normal distribution rule is implied, i.e. three numerical characteristics of a normal cloud model, denoted as N3 (E)x,E'n,H2 e) Wherein E isx、En、HeExpectation, entropy and hyper-entropy, respectively, referred to as normal clouds.
(1) Expectation of normal cloud Ex: representing the expectation of the spatial distribution of cloud droplets in the universe of discourse, i.e. the centroid of the area under the coverage of the cloud model
Figure BDA0002880515440000072
Is the information center value which can reflect the fuzzy concept most;
(2) entropy E of Normal cloudn: the uncertainty measure represents a qualitative concept and represents the value range size which can be accepted by the qualitative concept in the discourse space, namely the ambiguity;
(3) hyper-entropy of normal cloud He: i.e. points on the expected curve of the cloud model
Figure BDA0002880515440000073
The corresponding variance of the deterministic random distribution, which describes a measure of uncertainty in entropy, reflects the randomness of the appearance of the samples representing qualitative conceptual values, revealing a correlation of ambiguity and randomness.
C. MATLAB generates a corresponding cloud model diagram, and the test value certainty is read;
(C1) producing the expected value is ExEntropy is EnThe normal random number x of (a); calculating the expected value is ExEntropy is EnIs normal membership cloud, the degree of membership at the desired curve x
Figure BDA0002880515440000074
Respectively calculating the entropy E of the changes along the expected curves MA and MB according to a reduced semi-normal rulenl、Enr:
Figure BDA0002880515440000075
(c2) Calculating points
Figure BDA0002880515440000076
Variance σ of (c)x
Figure BDA0002880515440000077
Figure BDA0002880515440000081
(C3) Generating the expected value is
Figure BDA0002880515440000082
Variance is σxNormal random number μ of (a):
Figure BDA0002880515440000083
(C4) finally obtaining the cloud droplet xi (x, mu);
(C5) after the cloud model digital characteristic values of all indexes under different sub-level grades are obtained, the cloud model graph corresponding to the tunnel surrounding rock sub-level grading indexes is obtained after the cloud model digital characteristic values of all surrounding rock sub-level grading indexes are input by using a forward cloud generator;
after all the index numerical characteristic values calculated in table 2 are imported into MATLAB for calculation, cloud model diagrams of all five sub-level certainty degrees of iv-av-b corresponding to uniaxial compressive strength, integrity coefficient, cohesive force and softening coefficient in the sub-level hierarchical index system can be obtained, as shown in fig. 3 to fig. 6.
In the above steps, a normal cloud Generator (MCG) is an important implementation tool for cloud model theory, and the formula is as follows:
Figure BDA0002880515440000084
forward cloud generator is an expectation of a known cloud model ExEntropy EnAnd entropy HeAnd generating a two-dimensional point xi (x, mu) meeting the normal cloud distribution rule by using the three digital characteristic values, namely realizing the conversion from a qualitative concept to a quantitative value, wherein the two-dimensional point xi (x, mu) is also called cloud drop.
The normal cloud model is formed by thousands of cloud droplets generated by the forward cloud generator, and the mathematical model of the normal cloud generator is as follows without loss of generality:
(1) the expected curve of the normal cloud satisfies N (E)x,En 2) The normal distribution of (1) is shown in the above formula.
(2) Points on the desired curve
Figure BDA0002880515440000085
The random number is a normal random number with the variance of the point as the convergence center in the direction of certainty degree, and the normal distribution N (E) isx,En 2) Is generated from a normal distribution
Figure BDA0002880515440000091
And (4) generating.
(3) Varying along the expected curve of a normal cloud, reaching a maximum value σ at point Mx=HeAt points A and B σAσ B0. On both sides of point M, the curve is changed according to two decreasing semi-normal rules, and conforms to the "3 b rule".
D. Acquiring four test parameter values of the sub-level section to be divided through an indoor test;
(D1) adopting an indoor uniaxial compression test to obtain uniaxial compression strength and softening coefficient of the surrounding rock, and dividing the uniaxial compression test into a uniaxial compression test in a drying state, a uniaxial compression test in a water saturation state and a uniaxial compression test in a natural state according to different water containing states of the sample, wherein the ratio of the uniaxial compression strength in a rock saturation state to the uniaxial compression strength in a drying state is the softening coefficient of the rock;
(D2) after the longitudinal wave velocity of the rock mass and the elastic longitudinal wave velocity of the rock are respectively obtained through tests, the ratio of the longitudinal wave velocity to the elastic longitudinal wave velocity of the rock is the rock integrity coefficient;
(D3) obtaining corresponding shear strength indexes, namely cohesive force c and internal friction angle through direct shear test of rock
Figure BDA0002880515440000095
E. Determining each grading index weight; d, after obtaining a grading index test value of the required sub-grade grading section through an indoor test, determining the index weight through an entropy weight method by MATLAB software;
(E1) carrying out standardized processing on the data, and selecting j tunnel surrounding rock grading indexes, namely X1,X2,X3,···,XjIf i evaluation objects are set, the evaluation index value of the ith object is xi1,xi2,xi3,···,xijFormula y of normalization processijThe following were used:
Figure BDA0002880515440000092
(E2) according to the information entropy definition, calculating the information entropy corresponding to each index value of the n evaluation objects:
Figure BDA0002880515440000093
Figure BDA0002880515440000094
(E3) obtaining i groups of surrounding rocksWeight value omega of each index of grading objectj
Figure BDA0002880515440000101
F. Calculating comprehensive certainty and judging the grade;
(F1) firstly, calculating the test value certainty, obtaining the grading index test value of the grading section of the soft rock tunnel according to an indoor test, reading the test value certainty mu (x) of each grading index test value corresponding to each grade according to a generated table by using MATLAB software, and marking the test value certainty mu (x) as 0 if the test value does not belong to a certain grade;
(F2) calculating comprehensive certainty factor, determining degree mu (x) by the calculated test value, and calculating each grading index weight omega by combining entropy weight methodjAnd calculating the comprehensive certainty U of each section:
Figure BDA0002880515440000102
(F3) and judging the sub-grade level of the surrounding rock of each section, and judging the final sub-grade level attribution of the surrounding rock of each section according to the obtained numerical value of the comprehensive certainty factor U of each section.
The classification method of the invention is adopted to perform the sub-grade classification application of the surrounding rock on five sections selected from K9+525, K9+660, K9+910, K10+405 and K11+968 of the ancient top tunnel as follows:
and step A to step C are as described above, and in the step D, the rock samples of five sections are subjected to indoor tests, and the obtained grading index test values of the sections are arranged.
And E, according to the efficient operation of the entropy weight method realized by using MATLAB software in the step E, calculating the weights of the uniaxial compressive strength, the rock integrity coefficient, the cohesion and the softening coefficient, and finally obtaining the result shown in Table 3.
TABLE 3 test values and weights of ranking indices
Figure BDA0002880515440000103
And D, calculating the certainty factor according to the step F, judging the grade, and acquiring the certainty factor of the test value of each section corresponding to each grading index according to the research content in the step c, wherein the certainty factor is shown in the table 4, and the test value of each section corresponding to each grading index corresponds to each grading index.
TABLE 4 Each index certainty of five sections
Figure BDA0002880515440000111
Figure BDA0002880515440000121
From the calculated certainty μ (x), the comprehensive certainty U of each section was calculated in combination with the classification index weight calculated by the entropy weight method, and the result is shown in table 5.
TABLE 5 comprehensive determination of each section and quality grading results of surrounding rocks
Figure BDA0002880515440000122
From the results in Table 5, it was found that in the K9+525 section, U (V-a) > U (IV-c) > U (V-b) > U (IV-a) ═ U (IV-b), and therefore the grade of the surrounding rock was determined to be V-a in the K9+525 section. Similarly, the grade of the surrounding rock of the K9+660 section is V-b, the grade of the surrounding rock of the K9+910 section is V-a, the grade of the surrounding rock of the K10+405 section is IV-b, and the grade of the surrounding rock of the K11+968 section is IV-c.
Finally, according to the tunnel geological survey report, in the five sections selected at this time, the sections K9+525, K9+660 and K9+910 all belong to V-level surrounding rock, and the sections K10+405 and K11+968 belong to IV-level surrounding rock, so that the results of the five sections after a new grading index system and a cloud model grading method are adopted are consistent with the results of BQ-method surrounding rock grading in the geological survey report on the main level, and on the basis, the five sections are subjected to sub-level division by combining with the soft rock tunnel surrounding rock sub-level grading standard, so that the soft rock tunnel surrounding rock sub-level grading system has scientificity and practicability.

Claims (2)

1. A soft rock tunnel surrounding rock sub-grade grading method is characterized by comprising the following steps:
A. determining grading indexes and grading;
(A1) determining the value range of the uniaxial compressive strength Rc in IV and V-grade surrounding rock to be 0-60 MPa;
(A2) determining the value range of the rock integrity coefficient Kv at IV and V level surrounding rocks to be 0-0.75, and determining the value range of the cohesive force c at IV and V level surrounding rocks to be 0-0.7 MPa;
(A3) softening coefficient K of surrounding rockfThe value range of the IV-grade and V-grade surrounding rocks is determined to be 0-0.75;
B. calculating the numerical characteristic value of each grading index;
determining expected E of normal cloud with surrounding rock sub-grade grading index corresponding to certain grade standardxEntropy of normal cloud EnNormal cloud super entropy He
Figure FDA0002880515430000011
Figure FDA0002880515430000012
He=K
In the formula, CmaxAnd CminThe variable fuzzy threshold value is used for adjusting the variable fuzzy threshold value according to the constant K;
C. MATLAB generates a corresponding cloud model diagram, and the test value certainty is read;
(C1) producing the expected value is ExEntropy is EnThe normal random number x of (a); calculating the expected value is ExEntropy is EnNormal affiliation ofMembership at the desired curve x of the genus cloud
Figure FDA0002880515430000013
Respectively calculating the entropy E of the changes along the expected curves MA and MB according to a reduced semi-normal rulenl、Enr:
Figure FDA0002880515430000014
(c2) Calculating points
Figure FDA0002880515430000015
Variance σ of (c)x
Figure FDA0002880515430000021
Figure FDA0002880515430000022
(C3) Generating the expected value is
Figure FDA0002880515430000023
Variance is σxNormal random number μ of (a):
Figure FDA0002880515430000024
(C4) finally obtaining the cloud droplet xi (x, mu);
(C5) after the cloud model digital characteristic values of all indexes under different sub-level grades are obtained, the cloud model graph corresponding to the tunnel surrounding rock sub-level grading indexes is obtained after the cloud model digital characteristic values of all surrounding rock sub-level grading indexes are input by using a forward cloud generator;
D. acquiring four test parameter values of the sub-level section to be divided through an indoor test;
(D1) adopting an indoor uniaxial compression test to obtain uniaxial compression strength and softening coefficient of the surrounding rock, and dividing the uniaxial compression test into a uniaxial compression test in a drying state, a uniaxial compression test in a water saturation state and a uniaxial compression test in a natural state according to different water containing states of the sample, wherein the ratio of the uniaxial compression strength in a rock saturation state to the uniaxial compression strength in a drying state is the softening coefficient of the rock;
(D2) after the longitudinal wave velocity of the rock mass and the elastic longitudinal wave velocity of the rock are respectively obtained through tests, the ratio of the longitudinal wave velocity to the elastic longitudinal wave velocity of the rock is the rock integrity coefficient;
(D3) obtaining corresponding shear strength indexes, namely cohesive force c and internal friction angle through direct shear test of rock
Figure FDA0002880515430000025
E. Determining each grading index weight; d, after obtaining a grading index test value of the required sub-grade grading section through an indoor test, determining the index weight through an entropy weight method by MATLAB software;
(E1) carrying out standardized processing on the data, and selecting j tunnel surrounding rock grading indexes, namely X1,X2,X3,···,XjIf i evaluation objects are set, the evaluation index value of the ith object is xi1,xi2,xi3,···,xijFormula y of normalization processijThe following were used:
Figure FDA0002880515430000031
(E2) according to the information entropy definition, calculating the information entropy corresponding to each index value of the n evaluation objects:
Figure FDA0002880515430000032
Figure FDA0002880515430000033
(E3) obtaining weight values omega of all indexes of i groups of surrounding rock grading objectsj
Figure FDA0002880515430000034
F. Calculating comprehensive certainty and judging the grade;
(F1) firstly, calculating the test value certainty, obtaining the grading index test value of the grading section of the soft rock tunnel according to an indoor test, reading the test value certainty mu (x) of each grading index test value corresponding to each grade according to a generated table by using MATLAB software, and marking the test value certainty mu (x) as 0 if the test value does not belong to a certain grade;
(F2) calculating comprehensive certainty factor, determining degree mu (x) by the calculated test value, and calculating each grading index weight omega by combining entropy weight methodjAnd calculating the comprehensive certainty U of each section:
Figure FDA0002880515430000035
(F3) and judging the sub-grade level of the surrounding rock of each section, and judging the final sub-grade level attribution of the surrounding rock of each section according to the obtained numerical value of the comprehensive certainty factor U of each section.
2. The soft rock tunnel surrounding rock sub-grade grading method according to claim 1, characterized in that: and in the step B, K is uniformly 0.01[95 ].
CN202011632957.XA 2020-12-31 2020-12-31 Soft rock tunnel surrounding rock sub-level grading method Active CN112818439B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011632957.XA CN112818439B (en) 2020-12-31 2020-12-31 Soft rock tunnel surrounding rock sub-level grading method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011632957.XA CN112818439B (en) 2020-12-31 2020-12-31 Soft rock tunnel surrounding rock sub-level grading method

Publications (2)

Publication Number Publication Date
CN112818439A true CN112818439A (en) 2021-05-18
CN112818439B CN112818439B (en) 2023-07-28

Family

ID=75856585

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011632957.XA Active CN112818439B (en) 2020-12-31 2020-12-31 Soft rock tunnel surrounding rock sub-level grading method

Country Status (1)

Country Link
CN (1) CN112818439B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117077027A (en) * 2023-07-24 2023-11-17 西南交通大学 Surrounding rock sub-level grading method and device based on intelligent grading model grading probability

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019042483A2 (en) * 2017-08-30 2019-03-07 中铁工程装备集团有限公司 Tbm boring rock state real-time sensing system and method
CN110717689A (en) * 2019-10-18 2020-01-21 山西中煤平朔爆破器材有限责任公司 Method for evaluating explosibility of bench rock mass of strip mine rock by grades
WO2020125668A1 (en) * 2018-12-18 2020-06-25 中国铁建重工集团股份有限公司 Method and system for automatically identifying surrounding rock level by applying while-drilling parameters
US20200240268A1 (en) * 2019-01-24 2020-07-30 Huaneng Tibet Yarlungzangbo River Hydropower Development Investment Co., Ltd. Tunnel boring robot and remote mobile terminal command system

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019042483A2 (en) * 2017-08-30 2019-03-07 中铁工程装备集团有限公司 Tbm boring rock state real-time sensing system and method
WO2020125668A1 (en) * 2018-12-18 2020-06-25 中国铁建重工集团股份有限公司 Method and system for automatically identifying surrounding rock level by applying while-drilling parameters
US20200240268A1 (en) * 2019-01-24 2020-07-30 Huaneng Tibet Yarlungzangbo River Hydropower Development Investment Co., Ltd. Tunnel boring robot and remote mobile terminal command system
CN110717689A (en) * 2019-10-18 2020-01-21 山西中煤平朔爆破器材有限责任公司 Method for evaluating explosibility of bench rock mass of strip mine rock by grades

Non-Patent Citations (15)

* Cited by examiner, † Cited by third party
Title
于丽等: "岩质围岩施工阶段亚级分级的数量化理论研究", 《岩土力学》 *
于丽等: "岩质围岩施工阶段亚级分级的数量化理论研究", 《岩土力学》, no. 12, 10 December 2009 (2009-12-10) *
卢宇明: "区间数型隧道围岩综合分级云优化理论模型", 《科学技术与工程》 *
卢宇明: "区间数型隧道围岩综合分级云优化理论模型", 《科学技术与工程》, no. 09, 28 March 2020 (2020-03-28) *
卢宇明;: "区间数型隧道围岩综合分级云优化理论模型", 科学技术与工程, no. 09 *
吴友银;左清军;闫天玺;: "富水条件下板溪群板岩隧道围岩分级研究" *
吴友银;左清军;闫天玺;: "富水条件下板溪群板岩隧道围岩分级研究", 现代隧道技术, no. 03 *
李科: "基于熵权-云模型的隧道围岩分级方法研究", 《现代隧道技术》 *
李科: "基于熵权-云模型的隧道围岩分级方法研究", 《现代隧道技术》, no. 04, 15 August 2018 (2018-08-15), pages 75 - 81 *
王明年;刘大刚;刘彪;李海军;: "公路隧道岩质围岩亚级分级方法研究" *
王明年;刘大刚;刘彪;李海军;: "公路隧道岩质围岩亚级分级方法研究", 岩土工程学报, no. 10 *
王明年等: "公路隧道围岩亚级物理力学参数研究", 《岩石力学与工程学报》 *
王明年等: "公路隧道围岩亚级物理力学参数研究", 《岩石力学与工程学报》, no. 11, 15 November 2008 (2008-11-15) *
王明年等: "公路隧道岩质围岩亚级分级方法研究", 《岩土工程学报》 *
王明年等: "公路隧道岩质围岩亚级分级方法研究", 《岩土工程学报》, no. 10, 15 October 2009 (2009-10-15), pages 1590 - 1594 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117077027A (en) * 2023-07-24 2023-11-17 西南交通大学 Surrounding rock sub-level grading method and device based on intelligent grading model grading probability
CN117077027B (en) * 2023-07-24 2024-03-15 西南交通大学 Surrounding rock sub-level grading method and device based on intelligent grading model grading probability

Also Published As

Publication number Publication date
CN112818439B (en) 2023-07-28

Similar Documents

Publication Publication Date Title
Weiss Fracture and fragmentation of ice: a fractal analysis of scale invariance
CN109577972B (en) Glutenite reservoir rock mechanical parameter logging evaluation method based on lithology classification
Li et al. Rainfall and earthquake-induced landslide susceptibility assessment using GIS and Artificial Neural Network
CN110472363B (en) Surrounding rock deformation grade prediction method and system suitable for high-speed railway tunnel
CN109598015B (en) Grading evaluation method for rock mass fragmentation degree of fragmentation structure
CN110717689A (en) Method for evaluating explosibility of bench rock mass of strip mine rock by grades
Ou et al. Ubiquitiform in applied mechanics
Aghajani et al. Application of artificial neural network for calculating anisotropic friction angle of sands and effect on slope stability
CN112365054A (en) Comprehensive grading prediction method for deep well roadway surrounding rock
CN112818439A (en) Soft rock tunnel surrounding rock sub-grade grading method
Xu et al. Risk evaluation of debris flow hazard based on asymmetric connection cloud model
CN106872579B (en) Normal distribution is fitted the method that rock mass velocity divides rock-mass quality classification
Jamshidi Chenari et al. New method for estimation of the scale of fluctuation of geotechnical properties in natural deposits
CN109583003B (en) Face-crack polygon-based method for quantifying cracking degree of cracked structure rock mass
Bera et al. A multi-attribute decision making approach of mix design based on experimental soil characterization
Qu et al. Multi-axle moving train loads identification on simply supported bridge by using simulated annealing genetic algorithm
Rassoul et al. Predicting maximum dry density, optimum moisture content and California bearing ratio (CBR) through soil index using ordinary least squares (OLS) and artificial neural networks (ANNS)
CN105718668A (en) Open pit mine cast blasting effect analysis method
Ding et al. Evaluation of Landslide Susceptibility in Mountainous Areas of Changji City at the Northern Foot of Tianshan Mountain based on Coupled Model of Weight of Evidence and Shanon’s Entropy
Yu et al. Rock burst classification prediction method based on weight inverse analysis cloud model
Yin et al. Effect of fabric factors on the mechanical behavior of foliated rocks: A particle flow approach
CN117078106B (en) Comprehensive evaluation and index weight sensitivity analysis method for blasting rock mass quality
Liu et al. Cloud Model Membership Degree of Rock Slope Stability Evaluation: Method and a Case Study.
Peng Study on seismic stability of loess landslide based on fuzzy comprehensive evaluation
Yunxia et al. A grey fuzzy comprehensive model for evaluation of geological structure complexity

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant