CN112907130A - Construction risk grade determination method, device and medium based on shield construction parameters - Google Patents

Construction risk grade determination method, device and medium based on shield construction parameters Download PDF

Info

Publication number
CN112907130A
CN112907130A CN202110313415.4A CN202110313415A CN112907130A CN 112907130 A CN112907130 A CN 112907130A CN 202110313415 A CN202110313415 A CN 202110313415A CN 112907130 A CN112907130 A CN 112907130A
Authority
CN
China
Prior art keywords
risk
evaluation
shield
index
construction
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
Application number
CN202110313415.4A
Other languages
Chinese (zh)
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.)
Shanghai Jiaotong University
Shantou University
Original Assignee
Shanghai Jiaotong University
Shantou 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 Shanghai Jiaotong University, Shantou University filed Critical Shanghai Jiaotong University
Priority to CN202110313415.4A priority Critical patent/CN112907130A/en
Publication of CN112907130A publication Critical patent/CN112907130A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • 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/0635Risk analysis of enterprise or organisation activities
    • 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
    • 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
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/08Construction

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Entrepreneurship & Innovation (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Educational Administration (AREA)
  • Marketing (AREA)
  • Development Economics (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • Quality & Reliability (AREA)
  • Operations Research (AREA)
  • Game Theory and Decision Science (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Excavating Of Shafts Or Tunnels (AREA)
  • Lining And Supports For Tunnels (AREA)

Abstract

The invention provides a construction risk grade determination method based on shield construction parameters, which comprises the following steps: s1, identifying risks of the actual shield tunnel engineering, and establishing a risk index system consisting of shield tunnel engineering risk evaluation indexes; s2, determining the evaluation level of the shield tunnel engineering risk and the evaluation standard of the evaluation index according to the actual situation; s3, determining the weight of the shield tunnel engineering risk evaluation index in the step S1 by using an analytic hierarchy process; s4, based on the evaluation standard of S2, calculating the contact degree of each evaluation index in S1 by using an improved set pair analysis method; and S5, revising the contact degree of each evaluation index obtained in the S4 based on the weight obtained in the S3 to obtain a comprehensive contact degree, and determining the risk level of each ring of the shield tunnel according to the comprehensive contact degree and the risk evaluation level of the S2. The invention also provides a construction risk grade determination device and medium based on the shield construction parameters, and the accuracy of risk analysis is obviously improved.

Description

Construction risk grade determination method, device and medium based on shield construction parameters
Technical Field
The invention relates to the field of water resource environments, in particular to a construction risk grade determination method, a construction risk grade determination device and a construction risk grade determination medium based on shield construction parameters.
Background
With the continuous development of urban rail transit construction in China, the shield tunnel construction technology is more and more widely applied due to economy and high efficiency. The shield construction may face a lot of unpredictable dangers, and reasonable shield construction parameters play a crucial role in shield tunneling; improper setting of shield parameters may cause ground collapse or uplift, instability of an excavation surface, mud outburst and water inrush of a tunnel and the like, thereby causing construction period delay, cost overbooking and other losses. Therefore, reasonable analysis and research on shield tunneling parameters, risk evaluation and control research and engineering risk level determination in the construction stage are necessary, which has important influence on the safety and cost of the tunnel project. The importance of operation parameters in the shield construction process is emphasized in shield short-distance downward-penetrating existing subway shield tunnel construction parameter control published in the journal of southwest university of transportation by Marvinghui equal to 2018, and the research of the method has great engineering significance for ensuring smooth tunneling of tunnels. Therefore, in order to ensure the safety and economy of shield tunnel tunneling, a determination method for directly evaluating the construction risk level based on shield construction parameters is provided.
At present, the research on the shield parameter construction risk is generally divided into three methods, namely a deterministic method, a numerical simulation method and a data statistical analysis method. Deterministic methods, such as Analytic Hierarchy Process (AHP), are qualitative methods based on expert experience, but they are subjective and ignore uncertainty of risk; the numerical simulation method is a typical risk assessment quantitative method, but the establishment of a model needs to consider more relevant parameters; the data statistical analysis method is a method for performing statistical analysis on the shield tunnel construction parameters, and compared with the former two methods, the method for evaluating the shield tunnel construction risk is simpler and more reliable.
The search of documents in the prior art shows that a patent with the patent publication number of CN110059963A discloses a tunnel risk evaluation method based on a fuzzy polymorphic Bayesian network, a risk accident tree is established according to a historical case of a tunnel accident, and then the fuzzy polymorphic Bayesian network is established by the accident tree, however, the evaluation method has the defects of complicated steps, complex calculation and large workload, and the randomness of risk evaluation and the fuzziness of human thinking are not considered. At present, a technology for analyzing and researching tunnel construction risks based on shield tunneling parameters does not exist, so that a simple, efficient and reliable method for determining the shield tunnel construction risks is urgently needed to be provided for the shield tunneling parameters.
Disclosure of Invention
Aiming at the defects of the existing set pair analysis method in urban water quality risk level determination, the invention provides a determination method, a device and a medium for directly evaluating the construction risk level based on shield construction parameters.
According to a first aspect of the present invention, a determination method for directly evaluating a construction risk level based on shield construction parameters is provided, which includes:
s1, identifying risks of the actual shield tunnel engineering, selecting risk evaluation indexes, and establishing a risk index system consisting of the shield tunnel engineering risk evaluation indexes;
s2, determining the evaluation grade of the shield tunnel engineering risk and the evaluation standard of the evaluation index according to the actual situation based on the risk index system of S1;
s3, determining the weight of the shield tunnel engineering risk evaluation index in the step S1 by using an analytic hierarchy process;
s4, based on the evaluation standard of S2, calculating the contact degree of each evaluation index in S1 by using an improved set pair analysis method;
and S5, revising the contact degree of each evaluation index obtained in the S4 based on the weight obtained in the S3 to obtain a comprehensive contact degree, and determining the risk level of each ring of the shield tunnel according to the comprehensive contact degree and the evaluation level of S2.
Optionally, in S1, the risk identification is to investigate, classify and summarize information such as the cause and result of a potential and objective risk event of a project, so as to discover and identify risks during the shield tunnel construction.
The risk index system is an evaluation index system which is established by only considering objective factors and is used for comprehensive evaluation of shield tunnel construction risks.
Optionally, the determining the weight of the shield tunnel engineering risk evaluation index by using the analytic hierarchy process includes:
comparing all factors of the risk index system, and constructing pairwise judgment matrixes;
and calculating the maximum characteristic value and the corresponding characteristic vector by adopting a root method, wherein the component of the characteristic vector is the weight of the factor.
Optionally, the factors of the risk indicator system are compared with each other, wherein a scale of 1-9 is used to represent the relative importance among the factors.
Optionally, the computing the maximum eigenvalue and the corresponding eigenvector by using a root method includes:
multiplying elements in the judgment matrix by rows and opening the multiplied elements to the power of n, and normalizing the obtained vector to obtain a sequencing weight vector, namely a characteristic vector W, wherein n is the number of factors;
computing the maximum feature root λmaxDetermined by the following equation:
AW=λmaxW
in the formula, A is a judgment matrix.
Optionally, the calculating the degree of association of each evaluation index by using an improved set-pair analysis method includes: calculating an absolute value of a distance between an index measured value and a corresponding value in an evaluation standard, and improving the existing set pair analysis method based on an interval fuzzy number to obtain an improved contact membership degree, wherein the interval fuzzy number is a closed interval on an actual number axis, and the interval fuzzy number is adopted to evaluate a risk index; and generating interval fuzzy numbers by using two adjacent values in the index evaluation standard, and calculating a middle value between the two intervals.
Optionally, the calculating the degree of association of each evaluation index by using an improved set-to-analysis method specifically includes:
Figure BDA0002990859030000031
in the formula, mukIs a multi-element degree of relation, mu is more than or equal to 0k≤1;mk,mk-1,mk+1Are respectively (b)k-1,bk),(bk,bk+1),(bk+1,bk+2) Is a median value of
Figure BDA0002990859030000032
Optionally, the comprehensive degree of association is determined by the following formula:
Figure BDA0002990859030000033
in the formula, mupkFor evaluating the comprehensive degree of relation corresponding to the index p, mu is more than or equal to 0pk≤1;wjA weight vector as an evaluation index; mu.sjkThe degree of relation of the evaluation index j relative to the k level is obtained; n is the number of evaluation indexes; k is the number of evaluation grades. When mu ispkThe index p is close to k level when the value of (d) is close to 1.
Optionally, determining a risk level of each ring of the shield tunnel according to the comprehensive contact degree, wherein: and determining the risk level of each ring of the shield tunnel according to the comprehensive contact degree and the maximum membership principle.
According to a second aspect of the present invention, there is provided a determination apparatus for directly assessing a construction risk level based on shield construction parameters, comprising a memory, a processor and a computer program stored in the memory and executable on the processor, wherein the processor when executing the program is operable to execute the above determination method for directly assessing a construction risk level based on shield construction parameters.
According to a third aspect of the present invention, there is provided a computer readable storage medium having stored thereon a computer program, which when executed by a processor is operable to perform the above-mentioned determination method for directly assessing a construction risk level based on shield construction parameters.
Compared with the prior art, the invention has at least one of the following beneficial effects:
the method, the device and the medium of the invention determine the risk grade of each ring of the shield tunnel based on the shield tunneling parameters, comprehensively adopt the improved set-pair analysis method and the analytic hierarchy process, realize the subjective and objective combination, give consideration to the fuzziness and uncertainty of risk evaluation, effectively solve the problem of complex uncertainty and improve the accuracy of risk evaluation.
The method, the device and the medium improve the set-pair analysis method, do not need to consider index types, simplify the original set-pair analysis method, improve the safety and the stability of decision, and provide a more simple, convenient, reasonable and efficient new method for determining the shield tunnel construction risk.
Drawings
Other features, objects and advantages of the invention will become more apparent upon reading of the detailed description of non-limiting embodiments with reference to the following drawings:
FIG. 1 is a flow chart of a shield tunnel construction risk level determination method of the present invention;
FIG. 2 is a schematic diagram of a risk indicator system in accordance with an embodiment of the present invention;
FIG. 3 is a schematic diagram illustrating a relationship membership comparison between an improved set pair analysis and a conventional method according to an embodiment of the present invention;
fig. 4 is a schematic diagram of results of shield tunnel connectivity and construction risk level in the embodiment of the present invention.
Detailed Description
The present invention will be described in detail with reference to specific examples. The following examples will assist those skilled in the art in further understanding the invention, but are not intended to limit the invention in any way. It should be noted that variations and modifications can be made by persons skilled in the art without departing from the spirit of the invention. All falling within the scope of the present invention.
The embodiment of the invention provides a determination method for directly evaluating construction risk grade based on shield construction parameters. The method comprises the steps of firstly carrying out risk identification on the actual shield tunnel engineering, establishing a risk index system, simultaneously determining the evaluation level of the construction risk and the evaluation standard of the index according to the actual condition, then determining the weight of the risk evaluation index of the shield tunnel engineering by utilizing an analytic hierarchy process, calculating a connection coefficient by utilizing an improved set pair analytic process, and finally revising the connection degree of each evaluation index based on the weight to obtain a comprehensive connection degree, thereby determining the risk level of each ring of the shield tunnel.
In a preferred embodiment of the invention, the determination method for directly evaluating the construction risk level based on the shield construction parameters comprehensively adopts an improved Set Pair Analysis (SPA) method and an Analytic Hierarchy Process (AHP) by analyzing and researching the shield tunneling parameters, combines the principle and the objective, considers the fuzziness and the uncertainty of the risk evaluation, and obviously improves the accuracy of the risk analysis, thereby providing a new method which is simpler, more reasonable and more efficient for determining the shield tunnel construction risk. Specifically, referring to fig. 1, a flowchart of a method according to a preferred embodiment of the present invention is shown, in which the determination method for directly evaluating the construction risk level based on the shield construction parameters includes the following steps:
the first step is as follows: and collecting shield construction parameter data.
In this step, the shield construction parameters are specific parameters of the shield machine during tunneling, assembling, grouting and other operations at each ring of segments, and specific numerical values can be obtained by statistics from a shield construction record table (ring report). The shield construction actual measurement data collected in the step can be used for analyzing the subsequent steps, and is the basis for evaluating the risk grade of each ring of the shield by the method.
Specifically, the shield construction record table is a table recorded by a shield recorder during shield construction, and the recorder records one table (namely, a ring report) at each ring pipe piece. The shield construction record table records specific information and construction parameter data corresponding to the number of rings.
The second step is that: and carrying out risk identification on the shield tunnel engineering and establishing a risk index system.
The risk identification is to investigate, classify and summarize the information such as the cause and the result of the potential and objective risk events of the engineering, and can discover and identify the risk during the construction of the shield tunnel. The risk index system is an evaluation index system which is established by only considering objective factors and is used for comprehensive evaluation of shield tunnel construction risks.
Specifically, the objective factors refer to parameters in the shield tunneling process, and the evaluation index system is a single-layer index system. For example, the shield construction parameters include shield total thrust, tunneling speed, cutter head rotating speed, cutter head torque, screw machine rotating speed, screw machine torque, soil output, penetration degree, soil bin pressure, grouting amount, foam dosage, elevation deviation, horizontal deviation and the like, and the selected factor indexes are determined by the actual condition of the shield tunnel.
The third step: and determining the evaluation grade of the construction risk and the evaluation standard of the index.
And (4) dividing the construction risk evaluation grade according to the actual engineering condition by a risk assessment worker.
The evaluation standard of the index is determined by the hydrogeology and the surrounding environment of the actual engineering, namely the grading standard of each index parameter of the engineering example.
Specifically, the grade of the index evaluation standard should be consistent with the evaluation grade of the construction risk.
The fourth step: determining the weight of the shield tunnel engineering risk evaluation index by using an analytic hierarchy process, specifically:
1) and comparing all factors of the risk index system, and constructing a pairwise judgment matrix A.
Specifically, a scale of 1-9 may be used to represent the relative importance between the various factors.
2) And calculating the maximum characteristic value and the corresponding characteristic vector by adopting a root method, wherein the component of the characteristic vector is the weight of the factor.
Firstly multiplying elements in a judgment matrix according to rows and opening the multiplied elements to the power of n, and then normalizing the obtained vector to obtain a sequencing weight vector, namely a characteristic vector W, wherein n is the number of factors.
Computing the maximum feature root λmaxDetermined by the following equation:
AW=λmaxW
in the formula, A is a judgment matrix.
3) And (4) calculating a consistency ratio CR (consistency ratio) and carrying out consistency check on the judgment matrix A.
Specifically, when CR is less than or equal to 0.1, the consistency of the judgment matrix is considered to be acceptable, and the component of the obtained feature vector W is the weight of the index; when CR is more than or equal to 0.1, the judgment matrix is corrected properly, namely, each element in the judgment matrix A is checked, if the judgment matrix A has a judgment error, the element with the error needs to be corrected, and if the judgment matrix A has the error accumulation, the element in the matrix needs to be finely adjusted. Due to the complexity of shield construction and the diversity of human recognition, certain differences may be caused when the judgment matrix A is constructed by comparing two factors in the step 1), and the step can ensure that the inconsistency degree is within an allowable range.
The fifth step: and constructing a set pair based on the index measured value set and the index evaluation standard set.
The index measured value set is a set formed by actual operation parameters of the shield machine at each ring pipe piece in the construction and tunneling process, namely a set formed by shield construction parameter data collected in the first step.
The index evaluation criterion set is a set formed by evaluation criteria of the indexes, that is, a set formed by the grading criteria of each index parameter determined in the third step.
The set of the measured values of the indexes has a certain relation with the set of the evaluation standards.
A set pair is a pair of two sets that have some relationship.
And a sixth step: and calculating the degree of multivariate connection on the basis of the set pair obtained in the fifth step. The method comprises the following steps of improving the existing set pair analysis method based on interval fuzzy numbers, and establishing a mathematical system for calculating the multivariate contact degree by using the improved set pair analysis method, specifically:
1) an existing set pair analysis method is improved based on interval fuzzy numbers.
The interval fuzzy number is a closed interval on the real number axis, and the interval fuzzy number is adopted to evaluate the risk index. And generating interval fuzzy numbers by using two adjacent values in the index evaluation standard, and calculating a middle value between the two intervals.
The set pair analysis theory provides concepts of the degree of contact and the joint coefficient, and the concepts can express the uncertain relation between two sets in a set pair, namely an index measured value set and an index evaluation standard set, wherein the joint coefficient is the mathematical expression of the degree of contact.
The key point of the set-pair analysis method is the calculation of the degree of association, and the essence of the method is the distance between the measured value of the calculated index and the corresponding value in the evaluation standard. The index types of the existing set pair analysis method are divided into two types, namely a smaller optimal type and a larger optimal type, the mathematical expressions of the two types are different, and when the contact degree is calculated, the specific index types need to be judged in advance and then the corresponding expressions need to be selected. In the embodiment of the invention, the absolute value of the distance between the measured index value and the corresponding value in the evaluation standard is calculated, and the existing set pair analysis method is improved based on the interval fuzzy number, so that the improved contact membership degree can be obtained. Because the calculated contact degree is an absolute value, whether the index is smaller, more optimal or larger, more optimal is not needed to be considered, the original set pair analysis method is simplified, and the safety and the stability of decision making are improved.
2) Establishing a mathematical system for calculating the multivariate contact degree by adopting an improved set-pair analysis method, specifically:
Figure BDA0002990859030000071
in the formula, mukIs a multi-element degree of relation, mu is more than or equal to 0k≤1;mk,mk-1,mk+1Are respectively (b)k-1,bk),(bk,bk+1),(bk+1,bk+2) Is a median value of
Figure BDA0002990859030000072
The relation of each evaluation index (construction parameter) can be calculated through the formula.
The seventh step: and performing weighted synthesis on the contact degrees of the evaluation indexes obtained in the sixth step to obtain the comprehensive contact degree of the evaluation indexes, so as to evaluate the risk level.
In the step, the contact degrees of all the evaluation indexes are weighted and synthesized, wherein the weight of each evaluation index is calculated by utilizing an analytic hierarchy process in the fourth step.
Specifically, the comprehensive degree of association of each evaluation index is determined by the following formula:
Figure BDA0002990859030000073
in the formula, mupkFor evaluating the comprehensive degree of relation corresponding to the index p, mu is more than or equal to 0pk≤1;wjA weight vector as an evaluation index; mu.sjkThe degree of relation of the evaluation index j relative to the k level is obtained; n is the number of evaluation indexes, and K is the number of evaluation grades. When mu ispkThe index p is close to k level when the value of (d) is close to 1.
In a specific embodiment, the risk rating may be a risk rating of each ring of the shield tunnel determined according to the comprehensive degree of association and the maximum membership principle.
In another embodiment of the present invention, a construction risk level determining apparatus based on shield construction parameters is further provided, which includes a memory, a processor, and a computer program stored in the memory and executable on the processor, and the processor, when executing the program, may be configured to execute the above construction risk level determining method based on shield construction parameters.
In another embodiment of the present invention, a computer-readable storage medium is further provided, on which a computer program is stored, which when executed by a processor is operable to execute the above-mentioned construction risk level determination method based on shield construction parameters. Specifically, in order to better illustrate the above technical solution of the present invention, specific application examples are provided below:
an intercity railway tunnel is subjected to continuous shield construction work with the length of about 3.296km in an airport area, wherein the distance of 1.7km needs to pass through an airport flight area, namely, pass through airport taxiways and parking ramps in sections from ZDK81+150 to ZDK82+ 930. The tunnel is constructed by adopting an earth pressure balance shield, the buried depth is 11-22 meters, the inner diameter and the outer diameter of the tunnel are 7.7 meters and 8.5 meters respectively, and the buried depth of underground water is 0-14 meters. Complex geological conditions such as a weak composite stratum, a muddy stratum and the like exist below an airport flight area, and the main geological structures are that silty clay and sedimentary rock pass through a tunnel. The engineering adopts a determination method for directly evaluating the construction risk grade based on the shield construction parameters to carry out risk analysis on the right line of the tunnel, and the determination method comprises seven steps, as shown in figure 1.
The first step is as follows: and collecting shield construction parameter data.
In this embodiment, 2126 ring data is counted from the shield construction record table (ring report), wherein 2098 ring complete data is counted. The data content comprises construction data such as total thrust, tunneling speed, cutter torque, cutter rotating speed, top soil pressure, bottom soil pressure, shield machine posture, soil output, surrounding rock conditions, propulsion oil cylinder stroke, shield tail grease consumption, foaming agent consumption, adjacent segment deviation, splicing time, grouting hole number, grouting amount, shield tail grouting pressure and the like.
The second step is that: and carrying out risk identification on the shield tunnel engineering and establishing a risk index system.
In this embodiment, the evaluation index in the evaluation index system is the shield total thrust a1And a tunneling speed A2Rotating speed A of cutter head3Penetration A4Grouting pressure A5Earth pressure A6And shield machine attitude A7The total number of the shield construction parameters is 7, and the specific reference is shown in FIG. 2.
The third step: and determining the evaluation grade of the construction risk and the evaluation standard of the index.
In the embodiment, the construction risk evaluation grade is divided into five grades according to the actual situation of the project, namely safety (I), relatively safety (II), general safety (III), relatively danger (IV) and danger (V).
In this embodiment, the classification standards of the index parameters are divided according to the hydrogeology and the surrounding environment of the actual project, as shown in table 1 below.
TABLE 1 grading evaluation criteria for index parameters
Figure BDA0002990859030000081
Figure BDA0002990859030000091
The fourth step: and determining the weight of the shield tunnel engineering risk evaluation index by using an analytic hierarchy process.
In the embodiment, a is obtained by comparing 7 factors of a risk index system with each other by adopting a scale of 1-9 proportionsijAnd constructing a pairwise judgment matrix A. Then, the root method is adopted to calculate the maximum eigenvalue lambda of the rootmaxAnd corresponding eigenvector W, multiplying elements in the judgment matrix A by rows and opening the multiplied elements to the power of 7, and normalizing the obtained vector to obtain an ordering weight vector, namely the eigenvector W, wherein the component of the eigenvector W is the weight W of the factor.
The feature vector W is determined by the following formula:
Figure BDA0002990859030000092
maximum characteristic root λmaxDetermined by the following equation:
AW=λmaxW
in this embodiment, the weight W of the component derivation factor of the feature vector W is calculated, specifically the weight Wi(0.2584,0.2480,0.1575,0.1292,0.0891,0.0635,0.0543), i.e. shield total thrust a1And a tunneling speed A2Rotating speed A of cutter head3Penetration A4Grouting pressure A5Soil and its preparation methodPressure in the cabin A6And shield machine attitude A7The weights of these 7 evaluation indices are 0.2584,0.2480,0.1575,0.1292,0.0891,0.0635 and 0.0543, respectively.
In this embodiment, the judgment matrix is subjected to consistency check, and the consistency ratio CR is 0.0155 < 0.1, which indicates that the consistency of the judgment matrix is considered acceptable.
The fifth step: and constructing a set pair based on the index measured value set and the index evaluation standard set.
In this embodiment, the index measured value set a is composed of 2126 shield construction parameter data collected in the first step, the index evaluation criterion set B is composed of the index parameter grading criteria divided in the third step and shown in table 1, and a set pair H is (a, B).
And a sixth step: an existing set pair analysis method is improved based on interval fuzzy numbers, and a mathematical system for calculating the quinary relation degree is established by using the improved set pair analysis method.
In this embodiment, the risk level has five levels, and the interval fuzzy number is generated by using two adjacent values in the index evaluation criteria in table 1, and each index is respectively in the form of I e (b)1,-∞),II∈(b1,b2),III∈(b2,b3),IV∈(b3,b4) And V e (b)4, + ∞), wherein the specific value of each index parameter is b1,b2,b3,b4See table 1, while calculating the median m between the two intervals1,m2,m3,m4,m5Specifically, there are
Figure BDA0002990859030000101
In this embodiment, the relationship degree of each evaluation factor is defined as an interval distance, an absolute value of a distance between an index measured value and a corresponding value in an evaluation standard is calculated, and when K is 5, the comparison between the relationship membership degree of the improved set pair analysis method and the conventional method is shown in fig. 3, so that the calculated relationship degree is an absolute value. The degree of relation among the first level, the second level, the third level, the fourth level and the fifth level can be determined by mu1(x),μ2(x),μ3(x),μ4(x) And mu5(x) To obtain, specifically:
Figure BDA0002990859030000102
Figure BDA0002990859030000103
the seventh step: and performing weighted synthesis on the index contact degrees to obtain the comprehensive contact degrees of all the evaluation indexes, thereby determining the risk level.
In this embodiment, the association degree of each evaluation index is calculated in the sixth step, and the weight w of each evaluation index is calculated in combination with the fourth stepi(0.2584,0.2480,0.1575,0.1292,0.0891,0.0635,0.0543), the overall degree of connectivity can be determined by the following equation:
Figure BDA0002990859030000104
in this embodiment, according to the maximum membership principle, the risk level of each ring of the shield tunnel can be determined by using the obtained comprehensive contact degree. The comprehensive degree of contact of each ring of each evaluation index obtained by risk analysis and the evaluated risk level of each ring of the shield by the method are shown in fig. 4. As shown in fig. 4, in the 2126-ring shield tunnel under study, most of the shield tunnel construction sections are at risks of class II and class III, and there is no risk of class V, and for the tunnel ring determined as class IV, there is a relatively high construction risk, specifically: 10.4% is safe (class I), 52.7% is safer (class II), 32.1% is generally safe (class III), and 5.9% is more dangerous (class IV).
According to the specific application example, the method realizes subjective and objective combination, considers the fuzziness and uncertainty of risk evaluation, can effectively solve the problem of complex uncertainty, improves the accuracy of risk evaluation, does not need to consider index types, simplifies the original set-pair analysis method, improves the safety and stability of decision, and provides a new method which is simpler, more reasonable and more efficient for determining the shield tunnel construction risk.
The foregoing description of specific embodiments of the present invention has been presented. It is to be understood that the present invention is not limited to the specific embodiments described above, and that various changes and modifications may be made by one skilled in the art within the scope of the appended claims without departing from the spirit of the invention. The above-described preferred features may be used in any combination without conflict with each other.

Claims (10)

1. A construction risk grade determining method based on shield construction parameters is characterized by comprising the following steps:
s1, identifying risks of the actual shield tunnel engineering, selecting risk evaluation indexes, and establishing a risk index system consisting of the shield tunnel engineering risk evaluation indexes;
s2, determining the evaluation grade of the shield tunnel engineering risk and the evaluation standard of the evaluation index according to the actual situation based on the risk index system of S1;
s3, determining the weight of the shield tunnel engineering risk evaluation index in the step S1 by using an analytic hierarchy process;
s4, based on the evaluation standard of S2, calculating the contact degree of each evaluation index in S1 by using an improved set pair analysis method;
and S5, revising the contact degree of each evaluation index obtained in the S4 based on the weight obtained in the S3 to obtain a comprehensive contact degree, and determining the risk level of each ring of the shield tunnel according to the comprehensive contact degree and the evaluation level of S2.
2. The method for determining the construction risk level based on the shield construction parameters according to claim 1, wherein the determining the weight of the risk evaluation index of the shield tunnel engineering by using an analytic hierarchy process comprises:
comparing all factors of the risk index system, and constructing pairwise judgment matrixes;
and calculating the maximum characteristic value and the corresponding characteristic vector by adopting a root method, wherein the component of the characteristic vector is the weight of the factor.
3. The shield construction parameter-based construction risk level determination method according to claim 2, wherein the factors of the risk index system are compared with each other, wherein a scale of 1-9 is adopted to represent relative importance among the factors.
4. The method for determining the construction risk level based on the shield construction parameters according to claim 2, wherein the calculating the maximum eigenvalue and the corresponding eigenvector by using a root method comprises:
multiplying elements in the judgment matrix by rows and opening the multiplied elements to the power of n, and normalizing the obtained vector to obtain a sequencing weight vector, namely a characteristic vector W, wherein n is the number of factors;
computing the maximum feature root λmaxDetermined by the following equation:
AW=λmaxW
in the formula, A is a judgment matrix.
5. The method for determining the construction risk level based on the shield construction parameters according to claim 1, wherein the calculating the degree of association of each evaluation index by using an improved set pair analysis method comprises:
calculating the absolute value of the distance between the measured index value and the corresponding value in the evaluation standard, and improving the existing set pair analysis method based on the interval fuzzy number to obtain the improved contact membership degree, wherein,
the interval fuzzy number is a closed interval on the real number axis, and the interval fuzzy number is adopted to evaluate the risk index; and generating interval fuzzy numbers by using two adjacent values in the index evaluation standard, and calculating a middle value between the two intervals.
6. The method for determining the construction risk level based on the shield construction parameters according to claim 1, wherein the degree of association of each evaluation index is calculated by using an improved set-pair analysis method, specifically:
Figure FDA0002990859020000021
in the formula, mukIs a multi-element degree of relation, mu is more than or equal to 0k≤1;mk,mk-1,mk+1Are respectively (b)k-1,bk),(bk,bk+1),(bk+1,bk+2) Intermediate values of, i.e.
Figure FDA0002990859020000022
7. The method for determining the construction risk level based on the shield construction parameters according to claim 1, wherein the comprehensive degree of association is determined by the following formula:
Figure FDA0002990859020000023
in the formula, mupkFor evaluating the comprehensive degree of relation corresponding to the index p, mu is more than or equal to 0pk≤1;wjA weight vector as an evaluation index; mu.sjkThe degree of relation of the evaluation index j relative to the k level is obtained; n is the number of evaluation indexes, and K is the number of evaluation grades. When mu ispkThe index p is close to k level when the value of (d) is close to 1.
8. The method for determining the construction risk level based on the shield construction parameters according to claim 1, wherein the risk level of each ring of the shield tunnel is determined according to the comprehensive contact degree, wherein: and determining the risk level of each ring of the shield tunnel according to the comprehensive contact degree and the maximum membership principle.
9. A construction risk level determination apparatus based on shield construction parameters, comprising a memory, a processor and a computer program stored in the memory and operable on the processor, wherein the processor, when executing the program, is operable to perform the method of any one of claims 1 to 8.
10. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, is adapted to carry out the method of determining a construction risk level based on shield construction parameters according to any one of claims 1 to 8.
CN202110313415.4A 2021-03-24 2021-03-24 Construction risk grade determination method, device and medium based on shield construction parameters Pending CN112907130A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110313415.4A CN112907130A (en) 2021-03-24 2021-03-24 Construction risk grade determination method, device and medium based on shield construction parameters

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110313415.4A CN112907130A (en) 2021-03-24 2021-03-24 Construction risk grade determination method, device and medium based on shield construction parameters

Publications (1)

Publication Number Publication Date
CN112907130A true CN112907130A (en) 2021-06-04

Family

ID=76106171

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110313415.4A Pending CN112907130A (en) 2021-03-24 2021-03-24 Construction risk grade determination method, device and medium based on shield construction parameters

Country Status (1)

Country Link
CN (1) CN112907130A (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113919186A (en) * 2021-12-14 2022-01-11 中国民航大学 Event tree-based method for calculating severity of synthetic consequence of primary overrun event
CN114066183A (en) * 2021-11-02 2022-02-18 山东康威大数据科技有限公司 Tunnel monitoring system based on big data analysis
CN115115240A (en) * 2022-07-04 2022-09-27 广西大学 Urban subway shield tunnel construction risk evaluation method based on two-dimensional cloud model

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112100851A (en) * 2020-09-16 2020-12-18 西南科技大学 Method for evaluating tunnel water inrush disaster risk based on set pair analysis

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112100851A (en) * 2020-09-16 2020-12-18 西南科技大学 Method for evaluating tunnel water inrush disaster risk based on set pair analysis

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
任红岗等: "基于集对分析―灰色关联理论的采空区处理多属性方案综合评价", 《有色金属(矿山部分)》 *
刘沐宇等: "大跨径悬索桥施工风险评估的集对分析方法", 《土木工程与管理学报》 *
吕海敏: "城市地铁系统沉涝灾害风险评估方法与防灾对策", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 *
赵金先等: "基于AHP-SPA的地铁建设项目施工风险评价", 《土木工程与管理学报》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114066183A (en) * 2021-11-02 2022-02-18 山东康威大数据科技有限公司 Tunnel monitoring system based on big data analysis
CN113919186A (en) * 2021-12-14 2022-01-11 中国民航大学 Event tree-based method for calculating severity of synthetic consequence of primary overrun event
CN115115240A (en) * 2022-07-04 2022-09-27 广西大学 Urban subway shield tunnel construction risk evaluation method based on two-dimensional cloud model

Similar Documents

Publication Publication Date Title
CN109034582B (en) Tunnel crossing fault water inrush and mud inrush risk evaluation method based on cloud model and combined weighting
CN112907130A (en) Construction risk grade determination method, device and medium based on shield construction parameters
Li et al. Risk assessment of water inrush in karst tunnels and software development
CN103268420B (en) A kind of method for evaluating hazard of high rock slope
CN110610285A (en) Underground metal mine goaf risk grading evaluation method
Leng et al. A hybrid data mining method for tunnel engineering based on real-time monitoring data from tunnel boring machines
Wang et al. Risk Assessment of Water Inrush in Karst Tunnels Based on the Efficacy Coefficient Method.
CN106407493A (en) Multi-dimensional Gaussian cloud model-based rock burst grade evaluation method
Yang et al. The fuzzy comprehensive evaluation of water and sand inrush risk during underground mining
CN113570226A (en) Method for evaluating occurrence probability grade of tunnel water inrush disaster in fault fracture zone
CN116975576B (en) Mountain road geological disaster risk evaluation method based on key information statistics
Suh et al. Subsidence hazard assessment at the Samcheok coalfield, South Korea: a case study using GIS
CN114881396A (en) Tunnel collapse risk assessment method based on AHP and TOPSIS
Shi et al. Data-driven sequential development of geological cross-sections along tunnel trajectory
Du et al. Risk evaluation of colluvial cutting slope based on fuzzy analytic hierarchy process and multilevel fuzzy comprehensive evaluation
Han et al. Topological mapping of complex networks from high slope deformation time series for landslide risk assessment
Wang et al. Research on rock mass quality classification based on an improved rough set–cloud model
CN117409541A (en) Strip mine slope monitoring and early warning method and system
CN116109033B (en) Sea area sedimentary basin carbon dioxide geological sequestration suitability evaluation method
CN116663891A (en) Method, equipment and storage medium for determining safety of building structure in tunnel construction
Zhang et al. Risk assessment of coal and gas outburst in driving face based on finite interval cloud model
CN114493295B (en) Method for determining tunnel monitoring project under salt rock geology and monitoring system
Zhou et al. Developing a cloud evidence method for dynamic early warning of tunnel construction safety risk in undersea environment
Song et al. Using Bayesian network and Intuitionistic fuzzy Analytic Hierarchy Process to assess the risk of water inrush from fault in subsea tunnel
Jamshidi et al. Evaluating digital soil mapping approaches for 3D mapping of soil organic carbon

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

Application publication date: 20210604

RJ01 Rejection of invention patent application after publication