Calculation method for adaptability of small turning radius TBM tunneling in pumped storage power stations

Abstract

The present invention discloses a calculation method for the tunneling adaptability of a small turning radius TBM in a pumped storage power station, including determining the tunneling adaptability index of the small turning radius TBM in the pumped storage power station and constructing a hierarchical diagram; respectively establishing paired comparison tables between the base layer and the target layer and the index layer; converting the paired comparison tables into paired comparison matrices; respectively calculating the total weights of each element in the index layer and calculating the membership degrees of each index F <subgt;1< / subgt> to F <subgt;20< / subgt>; and finally calculating the tunneling adaptability of the small turning radius TBM in the pumped storage power station S . The present invention can scientifically calculate and evaluate the adaptability of the small turning radius TBM in the tunneling construction of the pumped storage power station project, and provide a scientific basis for the evaluation and optimization improvement of the TBM tunneling construction state.

Description

Technical Field

[0001] The present invention relates to the technical field of small turning radius TBM tunneling construction in pumped storage power stations, and particularly relates to a calculation method for the tunneling fitness of a small turning radius TBM in a pumped storage power station. Background Art

[0002] A full-face tunnel boring machine (TBM) integrates technologies such as hydraulics, control, machinery, and lasers, and can realize highly mechanized and automated excavation and support of tunnels. Due to its advantages such as fast tunneling speed, good construction environmental protection, high construction accuracy, and considerable comprehensive benefits, it has been widely used in engineering construction such as water conservancy and hydropower, highway engineering, railway engineering, urban rail transit, and oil and gas pipelines in recent years. As the length of tunnel construction becomes longer and longer, the advantages of TBM construction become more and more obvious, and the use of TBM construction has become an inevitable trend in long tunnel construction. Further, as an important part of the development of new energy technologies, pumped storage power stations can effectively improve the coordination of power grid operation and the safety and stability of power grid operation; the project of a pumped storage power station is a new field for the application of small turning radius TBM construction. Compared with traditional diversion tunnels, the project of a pumped storage power station has characteristics such as small line curves, frequent passing through stations, and limited site environment. The previous TBM design and construction experience is no longer generally applicable. Therefore, it is urgent to carry out relevant calculations on the tunneling fitness of small turning radius TBM in tunneling projects such as diversion and pumped storage to facilitate a scientific evaluation of the tunneling adaptability of TBM. Summary of the Invention

[0003] In view of this, the purpose of the present invention is to provide a calculation method for the tunneling fitness of a small turning radius TBM in a pumped storage power station to calculate and evaluate the tunneling adaptability of a small turning radius TBM in the operation of a pumped storage power station, aiming at the deficiencies of the prior art.

[0004] To achieve the above purpose, the present invention adopts the following technical solutions:

[0005] A calculation method for the adaptability of a small turning radius TBM tunneling in a pumped - storage power station, comprising the following steps: (1) Determine the adaptability indexes of the small turning radius TBM tunneling in the pumped - storage power station and construct a hierarchical diagram; the hierarchical diagram includes an objective layer, a criterion layer, and an index layer; the objective layer is to calculate the adaptability degree (S) of the small turning radius TBM tunneling in the pumped - storage power station, the criterion layer includes TBM tunneling construction parameters Q1, engineering geological conditions Q2, complex geological problems Q3, and construction technology and management level Q4, and the index layer respectively includes the elements belonging to each item of the criterion layer, that is, Q1 includes TBM tunneling speed F1, TBM tunneling thrust F2, TBM tunneling cutterhead rotation speed F3, and TBM tunneling cutterhead torque F4, Q2 includes uniaxial compressive strength of surrounding rock F5, integrity of surrounding rock F6, abrasiveness of surrounding rock F7, quartz content of surrounding rock F8, in - situ stress of stratum F9, and tunnel depth F 10 , Q3 includes stratum fault fracture zone F 11 , large deformation of surrounding rock F 12 , rock burst of stratum F 13 , composite stratum F 14 , sudden water inrush of tunnel F 15 , water permeability of surrounding rock F 16 , ambient temperature F 17 , and harmful gases F 18 , Q4 includes construction technology level F 19 and construction management level F 20 .

[0006] (2) Establish pairwise comparison tables between the criterion layer and the objective layer and the index layer respectively;

[0007] (3) Convert the pairwise comparison tables in step (2) into pairwise comparison matrices;

[0008] (4) Calculate the total weights of the elements in the index layer;

[0009] (5) Calculate the membership degrees of each index F1~F 20 ;

[0010] (6) Calculate the adaptability degree S of the small turning radius TBM tunneling in the pumped - storage power station,

[0011] P = ω·F T

[0012] where, ω=(ν1′ω1′…ν1′ω4′ ν2′ω5′…ν2′ω 10 ′ ν3′ω 11 ′…ν3′ω 18 ′ ν4′ω 19 ′ν4′ω 20 ′), F=(F1(x1) F2(x2)…F 20 (x 20))。

[0013] Preferably, the pairwise comparison table of the reference layer and the target layer in step (2) is as follows:

[0014] Table 1 Pairwise comparison table of the reference layer and the target layer

[0015]

[0016]

[0017] In the table, a ij are respectively the importance evaluation scales between the elements Q of the reference layer i Among them, a ij ≥0, a ij = 1 / a ji , a ii = 1.

[0018] The importance evaluation scale is determined according to the following table after statistical analysis based on social surveys, expert opinions, etc.

[0019] Table 2 Importance evaluation scale table

[0020]

[0021]

[0022] Preferably, the pairwise comparison table of the reference layer and the index layer in step (2) Q i -F are respectively:

[0023] Table 3 Pairwise comparison table of the index layer and the reference layer (Q1-F)

[0024] <![CDATA[Q1-F]]> <![CDATA[F1]]> <![CDATA[F2]]> <![CDATA[F3]]> <![CDATA[F4]]> <![CDATA[F1]]> <![CDATA[b 11 > <![CDATA[b 12 > <![CDATA[b 13 > <![CDATA[b 14 > <![CDATA[F2]]> <![CDATA[b 21 > <![CDATA[b 22 > <![CDATA[b 23 > <![CDATA[b 24 > <![CDATA[F3]]> <![CDATA[b 31 > <![CDATA[b 32 > <![CDATA[b 33 > <![CDATA[b 34 > <![CDATA[F4]]> <![CDATA[b 41 > <![CDATA[b 42 > <![CDATA[b 43 > <![CDATA[b 44 >

[0025] Table 4 Pairwise comparison table of the index layer and the reference layer (Q2-F)

[0026] <![CDATA[Q2-F]]> <![CDATA[F5]]> <![CDATA[F6]]> <![CDATA[F7]]> <![CDATA[F8]]> <![CDATA[F9]]> <![CDATA[F 10 <!-- 3 -->]]> <![CDATA[F5]]> <![CDATA[c 11 > <![CDATA[c 12 > <![CDATA[c 13 > <![CDATA[c 14 > <![CDATA[c 15 > <![CDATA[c 16 > <![CDATA[F6]]> <![CDATA[c 21 > <![CDATA[c 22 > <![CDATA[c 23 > <![CDATA[c 24 > <![CDATA[c 25 > <![CDATA[c 26 > <![CDATA[F7]]> <![CDATA[c 31 > <![CDATA[c 32 > <![CDATA[c 33 > <![CDATA[c 34 > <![CDATA[c 35 > <![CDATA[c 36 > <![CDATA[F8]]> <![CDATA[c 41 > <![CDATA[c 42 > <![CDATA[c 43 > <![CDATA[c 44 > <![CDATA[c 45 > <![CDATA[c 46 > <![CDATA[F9]]> <![CDATA[c 51 > <![CDATA[c 52 > <![CDATA[c 53 > <![CDATA[c 54 > <![CDATA[c 55 > <![CDATA[c 56 > <![CDATA[F 10 > <![CDATA[c 61 > <![CDATA[c 62 > <![CDATA[c 63 > <![CDATA[c 64 > <![CDATA[c 65 > <![CDATA[c 66 >

[0027] Table 5 Pairwise comparison table of the index layer and the reference layer (Q3-F)

[0028] <![CDATA[Q3-F]]> <![CDATA[F 11 > <![CDATA[F 12 > <![CDATA[F 13 > <![CDATA[F 14 > <![CDATA[F 15 > <![CDATA[F 16 > <![CDATA[F 17 > <![CDATA[F 18 > <![CDATA[F 11 > <![CDATA[d 11 > <![CDATA[d 12 > <![CDATA[d 13 > <![CDATA[d 14 > <![CDATA[d 15 > <![CDATA[d 16 > <![CDATA[d 17 > <![CDATA[d 18 > <![CDATA[F 12 > <![CDATA[d 21 > <![CDATA[d 22 > <![CDATA[d 23 > <![CDATA[d 24 > <![CDATA[d 25 > <![CDATA[d 26 > <![CDATA[d 27 > <![CDATA[d 28 > <![CDATA[F 13 > <![CDATA[d 31 > <![CDATA[d 32 > <![CDATA[d 33 > <![CDATA[d 34 > <![CDATA[d 35 > <![CDATA[d 36 > <![CDATA[d 37 > <![CDATA[d 38 > <![CDATA[F 14 > <![CDATA[d 41 > <![CDATA[d 42 > <![CDATA[d 43 > <![CDATA[d 44 > <![CDATA[d 45 > <![CDATA[d 46 > <![CDATA[d 47 > <![CDATA[d 48 > <![CDATA[F 15 > <![CDATA[d 51 > <![CDATA[d 52 > <![CDATA[d 53 > <![CDATA[d 54 > <![CDATA[d 55 > <![CDATA[d 56 > <![CDATA[d 57 > <![CDATA[d 58 > <![CDATA[F 16 > <![CDATA[d 61 > <![CDATA[d 62 > <![CDATA[d 63 > <![CDATA[d 64 > <![CDATA[d 65 > <![CDATA[d 66 > <![CDATA[d 67 > <![CDATA[d 68 > <![CDATA[F 17 > <![CDATA[d 71 > <![CDATA[d 72 > <![CDATA[d 73 > <![CDATA[d 74 > <![CDATA[d 75 > <![CDATA[d 76 > <![CDATA[d 77 > <![CDATA[d 78 > <![CDATA[F 18 > <![CDATA[d 81 > <![CDATA[d 82 > <![CDATA[d 83 > <![CDATA[d 84 > <![CDATA[d 85 > <![CDATA[d 86 > <![CDATA[d 87 > <![CDATA[d 88 >

[0029] Table 6 Pairwise comparison table of the index layer and the reference layer (Q4-F)

[0030] <![CDATA[Q4-F]]> <![CDATA[F 19 > <![CDATA[F 20 > <![CDATA[F 19 > <![CDATA[e 11 > <![CDATA[e 12 > <![CDATA[F 20 > <![CDATA[e 21 > <![CDATA[e 22 >

[0031] In the table, each element b ij , c ij , d ij , e ij and aij The definition is similar.

[0032] Preferably, the pairwise comparison matrix in step (3) is as follows:

[0033]

[0034]

[0035]

[0036]

[0037]

[0038] Preferably, the calculation method for the total weights of the elements in the index layer in step (4) is as follows:

[0039] ① Calculate the maximum eigenvalue λ ij ), 4×4 of matrices A = (a ij ), 4×4 B = (b ij ), 6×6 C = (C ij ), 8×8 D = (d ij ), 2×2 E = (e A,max ), B,max λ C,max ), D,max λ E,max ), A and their corresponding eigenvectors ν B = (ν1, ν2, ν3, ν4), ω C = (ω1, ω2, ω3, ω4), ω 10 ), ω D = (ω 11 , ω 12 ), 18 ω E = (ω 19 , ω 20 );

[0040] ② Normalize the eigenvectors ν A , ω B , ω C , ω D , ω D respectively to obtain

[0041] ν A ′ = (ν1′, ν2′, ν3′, ν4′), ω Bω′=(ω1′,ω2′,ω3′,ω4′), ω C ω′=(ω5′,ω6′,…,ω 10 ′),

[0042] ω D ω′=(ω 11 ′,ω 12 ′,…,ω 18 ′), ω E ω′=(ω 19 ′,ω 20 ′), then ν1′~ν4′ are the weights of each element Q1~

[0043] Q4 of the reference layer with respect to the target layer S, ω1′~ω 20 ′ are the weights of each element F1~F 20 of the index layer with respect to the reference layer Q1~Q4, ν1′ω1′~ν1′ω4′, ν2′ω5′~ν2′ω 10 ′, ν3′ω 11 ′~ν3′ω 18 ′, ν4′ω 19 ′, ν4′ω 20 ′

[0044] are the total weights of each element F1~F 20 of the index layer with respect to the target layer S. And denote

[0045] ω=(ν1′ω1′…ν1′ω4′ν2′ω5′…ν2′ω 10 ′ν3′ω 11 ′…ν3′ω 18 ′ν4′ω 19 ′ν4′ω 20 ′).

[0046] Preferably, the calculation of the membership degree in step (5) is as follows:

[0047] 1) Calculate the membership degree of the TBM tunneling speed F1,

[0048]

[0049] In the formula, the input value of the independent variable x1 is the TBM penetration index in kN / (mm / r), and σ c is the uniaxial compressive strength of the surrounding rock in MPa;

[0050] 2) Calculate the membership degree of the TBM tunneling thrust F2,

[0051]

[0052] In the formula, the input value of the independent variable x2 is the TBM tunneling thrust, with the unit of kN. Among them, F0 is the rated thrust of the TBM, with the unit of kN; F max is the maximum design thrust of the TBM, with the unit of kN.

[0053] 3) Calculate the membership degree of the TBM tunneling cutterhead rotation speed F3

[0054]

[0055] In the formula, the input value of the independent variable x3 is the TBM tunneling cutterhead rotation speed, with the unit of rad / min,

[0056] where n1 = 22.266D -0.7453 , n2 = 47.7 / D, n3 = 42.6D -0.5 , and D is the TBM diameter, with the unit of m;

[0057] 4) Calculate the membership degree of the TBM tunneling cutterhead torque F4,

[0058]

[0059] In the formula, the input value of the independent variable x4 is the TBM tunneling cutterhead torque, with the unit of kN·m;

[0060] where T min = 45D 2 , T max = 60D 2 , and D is the TBM diameter, with the unit of m;

[0061] 5) Calculate the membership degree of the uniaxial compressive strength of the surrounding rock F5,

[0062]

[0063] In the formula, the input value of the independent variable x5 is the uniaxial compressive strength of the surrounding rock, with the unit of MPa;

[0064] 6) Calculate the membership degree of the surrounding rock integrity F6,

[0065]

[0066] In the formula, the input value of the independent variable x6 is the surrounding rock integrity coefficient, dimensionless;

[0067] 7) Calculate the membership degree of the surrounding rock abrasiveness F7,

[0068]

[0069] In the formula, the input value of the independent variable x7 is the surrounding rock abrasiveness coefficient, dimensionless;

[0070] 8) Calculate the membership degree of the quartz content F8 of the surrounding rock,

[0071]

[0072] where the input value of the independent variable x8 is the percentage content of quartz in the surrounding rock;

[0073] 9) Calculate the membership degree of the in-situ stress F9 of the strata,

[0074]

[0075] where the input value of the independent variable x9 is the in-situ stress ratio T s = σ max / σ c where σ max is the maximum principal stress of the tunnel section, in MPa, and σ c is the uniaxial compressive strength of the rock, in MPa;

[0076] 10) Calculate the membership degree of the tunnel depth F 10 ,

[0077]

[0078] where the input value of the independent variable x 10 is the tunnel depth, in m;

[0079] where h = 0.45×2 s-1 ·ω, ω = 1 + i(D - 5), s is the surrounding rock grade, D is the tunnel diameter, in m, and when D < 5m, take i = 0.2, when D > 5m, take i = 0.1; in addition, when the surrounding rock grade is Ⅰ - Ⅲ, take H = 2h, when the surrounding rock grade is Ⅳ, Ⅴ, take H = 2.5h;

[0080] 11) Calculate the membership degree of the formation fault fracture zone F 11 ,

[0081]

[0082] where the input value of the independent variable x 11 is the width of the formation fault fracture zone, in m;

[0083] 12) Calculate the membership degree of the large deformation of the surrounding rock F 12 ,

[0084]

[0085] where the input value of the independent variable x 12The input value is determined by the surrounding rock grade. When the surrounding rock grade is Class I, take x = 1; when the surrounding rock grade is Class II, take x = 0.75; when the surrounding rock grade is Class III, take x = 0.5; when the surrounding rock grades are Class IV and V, take x = 0.25;

[0086] 13) Calculate the membership degree of the formation rockburst F 13 ;

[0087]

[0088] In the formula, the independent variable x 13 has an input value of the surrounding rock strength-stress ratio σ c / σ max , where σ max is the maximum principal stress of the tunnel cross-section, in MPa, and σ

[0089] is the uniaxial compressive strength of the surrounding rock, in MPa; c

[0090] 14) Calculate the membership degree of the composite formation F 14 ;

[0091]

[0092] In the formula, the independent variable x 14 has an input value of the formation composite ratio, dimensionless;

[0093] 15) Calculate the membership degree of the tunnel water inrush and breakthrough F 15 ;

[0094]

[0095] In the formula, the independent variable x 15 has an input value of the tunnel water inflow, in m 3 / d;

[0096] 16) Calculate the membership degree of the surrounding rock permeability F 16 ;

[0097]

[0098] In the formula, the independent variable x 16 has an input value of the surrounding rock permeability, in Lu;

[0099] 17) Calculate the membership degree of the ambient temperature F 17 ;

[0100]

[0101] In the formula, the independent variable x 17 has an input value of the ambient temperature inside the tunnel, in °C; ​

[0102] 18) Calculate the membership degree of harmful gas F 18 ,

[0103]

[0104] where the input value of the independent variable x 18 is the gas concentration in the tunnel (%);

[0105] 19) Calculate the membership degree of construction technology level F 19 ,

[0106]

[0107] where the input value of the independent variable x 19 is determined according to experience, and the specific values can be referred to the following table;

[0108] Table 7 Fitness values of construction technology level

[0109]

[0110] 20) Calculate the membership degree of construction management level F 20 ,

[0111]

[0112] where the input value of the independent variable x 20 is determined according to experience and can be referred to the following table.

[0113] Table 8 Fitness values of construction management level

[0114]

[0115] After calculating the fitness of the small turning radius TBM tunneling in the pumped-storage power station, evaluate the adaptability of the TBM tunneling according to Table 9.

[0116] Table 9 Evaluation comparison table of the adaptability of the small turning radius TBM tunneling in the pumped-storage power station

[0117] Level Ⅰ Ⅱ Ⅲ Ⅳ Ⅴ Adaptive evaluation Fully adapted Adapted Weakly adapted Not adapted Highly not adapted Overall fitness ≥0.9 0.6～0.9 0.4～0.6 0.1～0.4 ≤0.1

[0118] The beneficial effects of the present invention are as follows:

[0119] The present invention can scientifically calculate and evaluate the adaptability of the small turning radius TBM tunneling in the pumped-storage power station project. The calculation process comprehensively considers the TBM tunneling parameters, engineering geological conditions, complex geological problems, and the technical and management levels of the construction unit. The calculation results are scientific and reliable, and can provide a scientific basis for the evaluation of the TBM tunneling construction status and the optimization design. Description of the Drawings

[0120] Figure 1 This is the hierarchical diagram of the present invention. Detailed implementation manners

[0121] The present invention will be further described below in conjunction with the accompanying drawings and embodiments.

[0122] A calculation method for the adaptability of a small turning radius TBM tunneling in a pumped storage power station includes the following steps:

[0123] (1) Determine the adaptability indexes of the small turning radius TBM tunneling in the pumped storage power station, and construct a hierarchical diagram as Figure 1 shown; the hierarchical diagram includes an objective layer, a reference layer, and an index layer; the objective layer is to calculate the adaptability degree (S) of the small turning radius TBM tunneling in the pumped storage power station, the reference layer includes TBM tunneling construction parameters Q1, engineering geological conditions Q2, complex geological problems Q3, and construction technology and management level Q4, and the index layer respectively includes the elements belonging to each item of the reference layer, that is, Q1 includes TBM tunneling speed F1, TBM tunneling thrust F2, TBM tunneling cutterhead rotation speed F3, and TBM tunneling cutterhead torque F4, Q2 includes uniaxial compressive strength of surrounding rock F5, integrity of surrounding rock F6, abrasiveness of surrounding rock F7, quartz content of surrounding rock F8, in-situ stress of strata F9, and tunnel depth F 10 , Q3 includes formation fault fracture zone F 11 , large deformation of surrounding rock F 12 , rock burst of strata F 13 , composite stratum F 14 , sudden water inrush in tunnel F 15 , water permeability of surrounding rock F 16 , environmental temperature F 17 , and harmful gases F 18 , Q4 includes construction technology level F 19 and construction management level F 20 .

[0124] The relevant parameters for calculating the adaptability degree of this project are given in Table 10 as follows.

[0125] Table 10 Values of relevant indexes for adaptability degree calculation

[0126]

[0127]

[0128] (2) Respectively establish paired comparison tables between the reference layer and the objective layer and the index layer, as shown in Tables 11 - 15;

[0129] Table 11 Paired comparison table between the reference layer and the objective layer

[0130] S-Q <![CDATA[Q1]]> <![CDATA[Q2]]> <![CDATA[Q3]]> <![CDATA[Q4]]> <![CDATA[Q1]]> 1 1 / 2 1 / 2 4 <![CDATA[Q2]]> 2 1 1 / 2 5 <![CDATA[Q3]]> 2 2 1 7 <![CDATA[Q4]]> 1 / 4 1 / 5 1 / 7 1

[0131] Pairwise Comparison Table between Index Layer and Benchmark Layer (Q1-F)

[0132] Q1-F <![CDATA[F1]]> <![CDATA[F2]]> <![CDATA[F3]]> <![CDATA[F4]]> F1 1 1 4 5 F2 1 1 4 5 F3 1 / 4 1 / 4 1 2 F4 1 / 5 1 / 5 1 / 2 1

[0133] Table 13 Pairwise Comparison Table between Index Layer and Benchmark Layer (Q2-F)

[0134] <![CDATA[Q2-F]]> <![CDATA[F5]]> <![CDATA[F6]]> <![CDATA[F7]]> <![CDATA[F8]]> <![CDATA[F9]]> <![CDATA[F 10 > <![CDATA[F5]]> 1 4 7 7 6 6 <![CDATA[F6]]> 1 / 4 1 4 4 3 3 <![CDATA[F7]]> 1 / 7 1 / 4 1 1 1 / 2 1 / 2 <![CDATA[F8]]> 1 / 7 1 / 4 1 1 1 / 2 1 / 2 <![CDATA[F9]]> 1 / 6 1 / 3 2 2 1 1 <![CDATA[F 10 > 1 / 6 1 / 3 2 2 1 1

[0135] Table 14 Pairwise Comparison Table between Index Layer and Benchmark Layer (Q3-F)

[0136] <![CDATA[Q3-F]]> <![CDATA[F 11 > <![CDATA[F 12 > <![CDATA[F 13 > <![CDATA[F 14 > <![CDATA[F 15 > <![CDATA[F 16 > <![CDATA[F 17 > <![CDATA[F 18 > <![CDATA[F 11 > 1 1 1 / 2 2 3 4 7 8 <![CDATA[F 12 > 1 1 1 / 2 2 3 4 7 8 <![CDATA[F 13 > 2 2 1 3 4 6 7 8 <![CDATA[F 14 > 1 / 2 1 / 2 1 / 3 1 2 3 6 7 <![CDATA[F 15 > 1 / 3 1 / 3 1 / 4 1 / 2 1 2 5 6 <![CDATA[F 16 > 1 / 4 1 / 4 1 / 6 1 / 3 1 / 2 1 4 5 <![CDATA[F 17 > 1 / 7 1 / 7 1 / 7 1 / 6 1 / 5 1 / 4 1 2 <![CDATA[F 18 > 1 / 8 1 / 8 1 / 8 1 / 7 1 / 6 1 / 5 1 / 2 1

[0137] Table 15 Pairwise Comparison Table between Index Layer and Benchmark Layer (Q4-F)

[0138] <![CDATA[Q4-F]]> <![CDATA[F 19 > <![CDATA[F 20 > <![CDATA[F 19 > 1 4 <![CDATA[F 20 > 1 / 4 1

[0139] (3) Convert the pairwise comparison table in step (2) into a pairwise comparison matrix;

[0140]

[0141]

[0142]

[0143]

[0144]

[0145] Calculate the eigenvector corresponding to the maximum eigenvalue of each matrix, and perform normalization to obtain the weights of each element in the benchmark layer and the index layer,

[0146]

[0147] (4) Calculate the total weights of each element in the index layer;

[0148] Table 16 Total Weights of Each Element in the Index Layer for the Target Layer

[0149]

[0150] (5) Calculate the membership degrees of each index F1 to F20 according to formulas (6) to (25), as shown in Table 17;

[0151] Table 17 Membership Degrees of Each Element in the Index Layer

[0152] Index Membership degree Index Membership degree <![CDATA[TBM tunneling speed (F1(x1))]]> 0.385 <![CDATA[Fault fracture zone (F 11 (x 11 ))]]> 0 <![CDATA[TBM driving thrust (F2(x2))]]> 0.963 <![CDATA[Large deformation of surrounding rock (F 12 (x 12 ))]]> 0.667 <![CDATA[TBM cutterhead rotation speed (F3(x3))]]> 1 <![CDATA[Rock burst in strata (F 13 (x 13 ))]]> 0.85 <![CDATA[TBM cutterhead torque (F4(x4))]]> 1 <![CDATA[Composite stratum (F 14 (x 14 ))]]> 1 <![CDATA[Uniaxial compressive strength of surrounding rock (F5(x5))]]> 1 <![CDATA[Sudden water breakthrough in tunnel (F 15 (x 15 ))]]> 0.558 <![CDATA[Integrity of surrounding rock (F6(x6))]]> 1 <![CDATA[Permeability rate of surrounding rock (F 16 (x 16 ))]]> 0.964 <![CDATA[Abrasion of surrounding rock (F7(x7))]]> 0 <![CDATA[Ambient temperature (F 17 (x 17 ))]]> 0.915 <![CDATA[Quartz content of surrounding rock (F8(x8))]]> 0.636 <![CDATA[Harmful gas (F 18 (x 18 ))]]> 1 <![CDATA[In-situ stress of formation (F9(x9))]]> 0.473 <![CDATA[Construction technology level (F 19 (x 19 ))]]> 0.8 <![CDATA[Tunnel burial depth (F 10 (x 10 ))]]> 0.153 <![CDATA[Construction management level (F 20 (x 20 ))]]> 0.8

[0153] (6) Calculate the adaptability P of the small turning radius TBM tunneling in the pumped storage power station,

[0154] P = ω·F T = 0.722

[0155] where ω = (ν1′ω1′…ν1′ω4′ ν2′ω5′…ν2′ω 10 ′ ν3′ω 11 ′…ν3′ω 18 ′ ν4′ω 19 ′ν4′ω 20 ′), F = (F1(x1) F2(x2)…F 20 (x 20 ))。

[0156] The fitness of TBM tunneling for this engineering example is calculated and can be obtained according to Table 9. The adaptability of TBM tunneling for this engineering example is "adaptable", indicating that the comprehensive tunneling construction ability of the TBM is relatively strong. At the same time, it can also be seen from the fitness of each index that there are certain deficiencies in the coping ability of the TBM for this project in terms of surrounding rock abrasiveness and fault fracture zones. In addition, the tunneling construction fitness of similar small-turn TBM projects of pumped-storage power stations in the future can be calculated and determined according to this method, providing a scientific basis for the evaluation and optimization of TBM tunneling construction status.

[0157] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and are not intended to limit them. Any other modifications or equivalent replacements made by those of ordinary skill in the art to the technical solutions of the present invention shall fall within the scope of the claims of the present invention as long as they do not depart from the spirit and scope of the technical solutions of the present invention.

Claims

1. A calculation method for the adaptability of a small turning radius TBM tunneling in a pumped storage power station, characterized in that, It includes the following steps: (1) Determine the adaptability index of TBM tunneling with a small turning radius in a pumped-storage power station and construct a hierarchical diagram; the hierarchical diagram includes an objective layer, a reference layer, and an index layer; the objective layer is to calculate the adaptability degree S of TBM tunneling with a small turning radius in a pumped-storage power station, the reference layer includes TBM tunneling construction parameters Q1, engineering geological conditions Q2, complex geological problems Q3, and construction technology and management level Q4, and the index layer respectively includes the elements belonging to each item of the reference layer, that is, Q1 includes TBM tunneling speed F1, TBM tunneling thrust F2, TBM tunneling cutterhead rotation speed F3, and TBM tunneling cutterhead torque F4, Q2 includes uniaxial compressive strength of surrounding rock F5, integrity of surrounding rock F6, abrasiveness of surrounding rock F7, quartz content of surrounding rock F8, in-situ stress of stratum F9, and tunnel depth F 10 , Q3 includes stratum fault fracture zone F 11 , large deformation of surrounding rock F 12 , rock burst of stratum F 13 , composite stratum F 14 , sudden water inrush in tunnel F 15 , permeability of surrounding rock F 16 , ambient temperature F 17 , and harmful gases F 18 , Q4 includes construction technology level F 19 , and construction management level F 20 ; (2) Establish pairwise comparison tables between the reference layer and the target layer and the index layer respectively; (3) Convert the pairwise comparison tables in step (2) into pairwise comparison matrices; (4) Calculate the total weights of the elements in the index layer. The calculation method of the total weights of the elements in the index layer is as follows: ① Calculate the pairwise comparison matrix A=(a ij ) 4×4 , B=(b ij ) 4×4 、C=(C ij ) 6×6 、D=(d ij ) 8×8 、E=(e ij ) 2×2 The largest characteristic root λ A,max , B,max , C,max , D,max , E,max and its corresponding eigenvector ν A =(ν1,ν2,ν3,ν4),ω B =(ω1,ω2,ω3,ω4),ω C =(ω5,ω6,…,ω 10 ),ω D =(ω 11 ,ω 12 ,…,ω 18 ),ω E =(ω 19 ,ω 20 ); ②Normalize the eigenvectors ν A , ω B , ω C , ω D , ω D respectively to obtain ν A ′ = (ν1′, ν2′, ν3′, ν4′), ω B ′ = (ω1′, ω2′, ω3′, ω4′), ω C ′ = (ω5′, ω6′, …, ω 10 ′), ω D ′ = (ω 11 ′, ω 12 ′, …, ω 18 ′), ω E ′ = (ω 19 ′, ω 20 ′). Then ν1′ to ν4′ are the weights of each element Q1 to Q4 in the benchmark layer for the target layer S in turn, and ω1′ to ω 20 ′ are the weights of each element F1 to F 20 in the index layer for the benchmark layer Q1 to Q4 in turn. ν1′ω1′ to ν1′ω4′, ν2′ω5′ to ν2′ω 10 ′, ν3′ω 11 ′ to ν3′ω 18 ′, ν4′ω 19 ′, ν4′ω 20 ′ are the total weights of each element F1 to F 20 in the index layer for the target layer S in turn. Denote ω = (ν1′ω1′…ν1′ω4′ν2′ω5′…ν2′ω 10 ′ν3′ω 11 ′…ν3′ω 18 ′ν4′ω 19 ′ν4′ω 20 ′); (5) Calculate the membership degrees of each index F1 to F 20 ; (6) Calculate the tunneling fitness P of the pumped-storage power station with a small turning radius TBM; P = ω·F T where ω = (ν1′ω1′…ν1′ω4′ν2′ω5′…ν2′ω 10 ′ν3′ω 11 ′…ν3′ω 18 ′ν4′ω 19 ′ν4′ω 20 ′), F = (F1(x1)F2(x2)…F 20 (x 20 ))), ν1′ω1′ is the TBM tunneling speed, ν1′ω2′ is the TBM tunneling thrust, ν1′ω3′ is the TBM cutterhead rotation speed, ν1′ω4′ is the TBM cutterhead torque, ν2′ω5′ is the uniaxial compressive strength of the surrounding rock, ν2′ω6′ is the integrity coefficient of the surrounding rock, ν2′ω7′ is the abrasion coefficient of the surrounding rock, ν2′ω8′ is the quartz content of the surrounding rock, ν2′ω9′ is the in-situ stress level of the stratum, ν2′ω 10 ′ is the tunnel depth, ν3′ω 11 ′ is the width of the fault fracture zone of the stratum, ν3′ω 12 ′ is the large deformation condition of the surrounding rock, ν3′ω 13 ′ is the rock burst of the stratum, ν3′ω 14 ′ is the proportion of composite strata, ν3′ω 15 ′ is the sudden water inflow of the tunnel, ν3′ω 16 ′ is the permeability of the surrounding rock, ν3′ω 17 ′ is the environmental temperature inside the tunnel, ν3′ω 18 ′ is the emission amount of harmful gases inside the tunnel, ν4′ω 19 ′ is the construction technology level, ν4′ω 20 ′ is the construction management level.

2. The calculation method for the adaptability of a pumped-storage power station's TBM tunneling with a small turning radius according to claim 1, wherein The pairwise comparison table between the reference layer and the target layer in step (2) is as follows: In the table, a ij are respectively the importance evaluation scales between each element Q i in the reference layer. Among them, a ij ≥0, a ij = 1 / a ji , a ii = 1.

3. The calculation method for the adaptability of a pumped-storage power station's TBM tunneling with a small turning radius according to claim 2, characterized in that, Pairwise comparison table Q of the reference layer and the index layer in step (2) i -F are respectively: Each element b in the table ij , c ij , d ij , e ij is similar to the definition of a ij .

4. The calculation method for the adaptability of a pumped-storage power station's small turning radius TBM tunneling according to claim 1, wherein The pairwise comparison matrix in step (3) is as follows:

5. The calculation method for the adaptability of a pumped-storage power station's small turning radius TBM tunneling according to claim 1, wherein The calculation of the membership degree in step (5) is as follows: 1) Calculate the membership degree of the TBM tunneling speed F1; In the formula, the input value of the independent variable x1 is the TBM penetration index with the unit of kN / (mm / r), and σ c is the uniaxial compressive strength of the surrounding rock, with the unit of MPa; 2) Calculate the membership degree of the TBM tunneling thrust F2; In the formula, the input value of the independent variable x2 is the TBM tunneling thrust, with the unit of kN; where F0 is the rated thrust of the TBM, with the unit of kN; F max is the maximum design thrust of the TBM, with the unit of kN; 3) Calculate the membership degree of the TBM cutterhead rotation speed F3 where the input value of the independent variable x3 is the TBM cutterhead rotation speed, with the unit of rad / min; Among them, n1 = 22.266D -0.7453 , n2 = 47.7 / D, n3 = 42.6D -0.5 , where D is the diameter of the TBM, with the unit of m; 4) Calculate the membership degree of the TBM cutterhead torque F4; where the input value of the independent variable x4 is the TBM cutterhead torque, with the unit of kN·m; Among them, T min = 45D 2 , T max = 60D 2 , where D is the diameter of the TBM, in m; 5) Calculate the membership degree of the uniaxial compressive strength F5 of the surrounding rock; where the input value of the independent variable x5 is the uniaxial compressive strength of the surrounding rock, with the unit of MPa; 6) Calculate the membership degree of the surrounding rock integrity F6; where the input value of the independent variable x6 is the surrounding rock integrity coefficient, dimensionless; 7) Calculate the membership degree of the surrounding rock abrasiveness F7; where the input value of the independent variable x7 is the surrounding rock abrasiveness coefficient, dimensionless; 8) Calculate the membership degree of the quartz content F8 of the surrounding rock; where the input value of the independent variable x8 is the percentage of quartz in the surrounding rock; 9) Calculate the membership degree of the in-situ stress F9 of the stratum In the formula, the input value of the independent variable x9 is the formation stress ratio T s = σ max / σ c , where σ max is the maximum principal stress of the tunnel section, with the unit of MPa, and σ c is the uniaxial compressive strength of the rock, with the unit of MPa; 10) Calculate the membership degree of the buried depth F of the tunnel 10 of where the independent variable x 10 has an input value of the tunnel burial depth, in m; where h = 0.45 × 2 s-1 ·ω, ω = 1 + i(D - 5), s is the surrounding rock grade, D is the tunnel diameter in m, and when D < 5m, i = 0.2 is taken, when D > 5m, i = 0.1 is taken; in addition, when the surrounding rock grade is from grade Ⅰ to Ⅲ, H = 2h is taken, when the surrounding rock grade is Ⅳ or Ⅴ, H = 2.5h is taken; 11) Calculate the membership degree of the formation fault fracture zone F 11 of, In the formula, the independent variable x 11 has an input value of the width of the formation fault fracture zone, with the unit of m; 12) Calculate the membership degree of the large deformation F of the surrounding rock 12 of In the formula, the independent variable x 12 's input value is determined by the surrounding rock grade. When the surrounding rock grade is grade I, take x 12 = 1. When the surrounding rock grade is grade II, take x 12 = 0.

75. When the surrounding rock grade is grade III, take x 12 = 0.

5. When the surrounding rock grade is grades IV and V, take x 12 = 0.25; 13) Calculate the membership degree of the formation rockburst F 13 and In the formula, the independent variable x 13 has an input value of the surrounding rock strength-stress ratio σ c / σ max , where σ max is the maximum principal stress of the tunnel section, with the unit of MPa, and σ c is the uniaxial compressive strength of the surrounding rock, with the unit of MPa; 14) Calculate the membership degree of the composite stratum F 14 and where the independent variable x 14 has an input value of formation composite ratio, dimensionless; 15) Calculate the membership degree of sudden water inflow F in the tunnel 15 of In the formula, the independent variable x 15 has an input value of the water inflow of the tunnel, with the unit of m 3 / d; 16) Calculate the membership degree of the water permeability rate F of the surrounding rock 16 ​ where the independent variable x 16 has an input value of the surrounding rock permeability in Lu 17) Calculate the membership degree of the ambient temperature F 17 and where the independent variable x 17 has an input value of the environmental temperature in the tunnel, with the unit of °C; 18) Calculate the membership degree of harmful gas F 18 and where the independent variable x 18 has an input value of the gas concentration percentage in the tunnel; 19) Calculate the membership degree of the construction technology level F 19 of In the formula, the input value of the independent variable x 19 is determined according to experience; 20) Calculate the membership degree of the construction management level F 20 of, In the formula, the input value of the independent variable x 20 is determined according to experience.

Images