Optimal selection method for setting scheme of auxiliary tunnel of extra-long railway tunnel and computer equipment

Abstract

The invention discloses an extra-long railway tunnel auxiliary tunnel setting scheme optimization method and computer equipment. The method comprises the following steps: establishing an auxiliary tunnel scheme evaluation index system which comprises a plurality of first-level indexes, wherein each first-level index comprises a plurality of second-level indexes; according to geological conditions, construction requirements and technical feasibility of a project site, setting alternative schemes for the auxiliary tunnel of the extra-long railway tunnel; the evaluation index system is used for collecting data related to each secondary index of each auxiliary tunnel setting scheme to be evaluated, calculating corresponding numerical values and forming an evaluation matrix X; on the basis of the matrix X, using a CRITIC method to calculate the weight of each first-level index and the weight of each second-level index to obtain a comprehensive weight matrix W; and on the basis of the comprehensive weight matrix W, a comprehensive evaluation value of each scheme is obtained through calculation by using a VIKOR method, and optimization is carried out. According to the method, comprehensive evaluation is carried out through the CRITIC method in combination with the VIKOR method, so that the auxiliary tunnel arrangement scheme of the extra-long railway tunnel is scientifically and effectively optimized.

Description

Technical Field

[0001] The present invention relates to the technical field of tunnel construction, and specifically to an optimal selection method for the auxiliary adit setting scheme of extra-long railway tunnels and computer equipment. Background Art

[0002] In mountain railway construction, in order to cope with complex terrains and difficult construction conditions and ensure the smooth progress of extra-long tunnel projects as planned, the construction strategy of "shortening long tunnels by building auxiliary adits" is usually adopted. These auxiliary adits not only provide conditions for multi-point parallel excavation of the main tunnel, significantly accelerating the construction progress, but also greatly facilitate the transportation of construction materials and the efficient external transportation of tunneling muck. Auxiliary adits are also crucial for ensuring construction safety and future operation safety. They provide additional guarantees for the ventilation and drainage systems of the tunnel, ensuring the health and safety of the construction environment; after the tunnel is put into operation, the auxiliary adits can also be used as emergency evacuation channels, improving the safety and emergency response capabilities of the tunnel. Building auxiliary adits requires comprehensive consideration of various factors such as construction period, cost, safety, and environmental protection. In terms of the construction period, the construction time of the auxiliary adit directly affects the construction progress of the main tunnel; in terms of cost, the construction of the auxiliary adit involves a large amount of resource investment, and cost control is a major challenge; in terms of safety, the design and construction of the auxiliary adit must meet strict safety standards to prevent accidents; in terms of environmental protection, the impact on the surrounding environment also needs to be considered when building the auxiliary adit, and measures should be taken to minimize ecological damage. These complex requirements make the planning and design of auxiliary adits extremely challenging, and there is an urgent need for a comprehensive and scientific optimal selection method to support decision-making. Summary of the Invention

[0003] The technical problem to be solved by the present invention is to provide a scientific and comprehensive optimal selection method for the auxiliary adit setting scheme of extra-long railway tunnels and computer equipment.

[0004] To solve the above technical problem, the present invention adopts the following technical solutions: An optimal selection method for the auxiliary adit setting scheme of extra-long railway tunnels, comprising the following steps:

[0005] S1. Establish an evaluation index system for the auxiliary adit scheme, which includes multiple first-level indicators, and each first-level indicator includes multiple second-level indicators;

[0006] S2. According to the geological conditions, construction requirements, and technical feasibility of the project site, formulate alternative schemes for the setting of auxiliary adits for extra-long railway tunnels;

[0007] S3. Based on the evaluation index system established in S1, collect data related to each second-level indicator of each auxiliary adit setting scheme to be evaluated, calculate the corresponding values, and form an evaluation matrix X;

[0008] S4. Based on matrix X, use the CRITIC method to calculate the weights of each first-level index and second-level index, and obtain the comprehensive weight matrix W;

[0009] S5. Based on the comprehensive weight matrix W, use the VIKOR method to calculate the comprehensive evaluation value of each scheme and conduct optimization.

[0010] Furthermore, the first-level indicators specifically include the construction period indicator A 1 , the cost indicator A 2 , the safety indicator A 3 and the environmental protection indicator A 4 ;

[0011] The construction period indicator A 1 includes the total project construction period A 11 , the construction period of auxiliary adits A 12 , the construction time of construction access roads A 13 , the construction time of auxiliary facilities A 14 , the slag removal time from the work area A 15 These are the five second-level indicators;

[0012] The cost indicator A 2 includes the cost of opening a working face A 21 , the construction cost of auxiliary adits A 22 , the construction cost of construction access roads A 23 , the construction cost of auxiliary facilities A 24 , the land acquisition cost A 25 , the land reclamation cost A 26 , the ventilation cost of the work area A 27 , the slag removal cost from the work area A 28 , the drainage cost of the work area A 29 These are the nine second-level indicators;

[0013] The safety indicator A 3 includes the safety of the entrance slope A 31 , the risk distribution in the auxiliary adit A 32 , the safety level of slag removal A 33 , the emergency rescue ability A 34 These are the four second-level indicators;

[0014] The environmental protection indicator A 4 includes the occupied area of cultivated land A 41 , the occupied area of forest land A 42 , the energy consumption of slag removal from the work area A 43 , the energy consumption of ventilation in the work area A 44 , the energy consumption of drainage in the work area A 45 These are the five second-level indicators.

[0015] Furthermore, the collection and calculation of each second-level indicator in S3 are as follows:

[0016] Total project duration A 11 :

[0017] Its calculation formula is:

[0018] max T i (1)

[0019] Where T i represents the construction duration of the i-th working face, and its calculation method is:

[0020]

[0021] Where Q ri represents the engineering quantity of the r-class surrounding rock of the i-th working face, V ri represents the construction rate of the unit engineering quantity of the i-th working face under the r-class surrounding rock, and ATi represents the construction time of the auxiliary adit of the i-th working face;

[0022] Auxiliary adit construction duration A 12 :

[0023] Its calculation method is:

[0024]

[0025] Where, AQ ri represents the engineering quantity of the r-class surrounding rock of the auxiliary adit of the i-th working face, AV ri represents the construction rate of the unit engineering quantity of the auxiliary adit of the i-th working face under the r-class surrounding rock;

[0026] Construction access road construction time A 13 :

[0027] Its calculation method is:

[0028]

[0029] Where RT i represents the construction time of the construction access road of the i-th working face, RQ i represents the engineering quantity of the construction access road of the i-th working face, RV i represents the construction rate of the unit engineering quantity of the construction access road of the i-th working face;

[0030] Ancillary facilities construction time A 14 :

[0031] Its calculation formula is:

[0032]

[0033] AF iis the construction time of the auxiliary facilities for the $i$-th work area, $W$ i is the total project volume of the $i$-th work area, $C$ i is the construction capacity of the auxiliary facilities construction, $\alpha$ is the proportionality coefficient related to the type of auxiliary facilities, which is used to reflect the input ratio required for different types of auxiliary facilities;

[0034] The slag transportation time out of the work area $A$ 15 :

[0035] Its calculation formula is:

[0036]

[0037] where $OT$ i is the slag transportation time out of the $i$-th work area, $OQ$ i is the total amount of muck in the $i$-th work area, $OV$ i is the muck transportation rate of the $i$-th work area, $L$ i is the length of the $i$-th work area, $\beta_O$ is the coefficient of the slag transportation time increasing with the length;

[0038] The cost of opening a working face $A$ 21 :

[0039] Its calculation formula is:

[0040] $C$ WFi $=$ $N$ i $\times$ $W$ i+ $M$ im $\times$ $MC$ im (7)

[0041] where, $C$ WFi is the total cost of the $i$-th working face, $N$ i is the total number of people required for the $i$-th working face, $W$ i is the per capita labor cost of the $i$-th working face, $M$ im is the total number of m types of mechanical equipment required for the $i$-th working face, $MC$ im is the unit mechanical cost of the m types of machinery for the $i$-th working face;

[0042] The cost of auxiliary adit construction $A$ 22 :

[0043] Its calculation formula is:

[0044]

[0045] where, $C$ APi is the total cost of the $i$-th auxiliary adit, $AQ$ ri is the project volume of the $i$-th auxiliary adit under the r-class surrounding rock grade, $AU$ riThe comprehensive unit price per running meter of the i-th auxiliary adit project volume under the r-class surrounding rock grade;

[0046] Construction access road construction cost A 23 :

[0047] Its calculation formula is:

[0048] C RDi =RQ i ×RU i (9)

[0049] Wherein, C RDi is the total construction access road construction cost of the i-th work area, RQi is the construction access road project volume of the i-th work area, and RUi is the comprehensive unit price of the unit project volume of the construction access road in the i-th work area;

[0050] Ancillary facilities construction cost A 24 :

[0051] Its calculation formula is:

[0052]

[0053] Wherein, AC j is the unit construction cost of the j-th ancillary facility, and Q j is the project volume of the corresponding ancillary facility;

[0054] Land acquisition cost A 25 :

[0055] Its calculation formula is:

[0056] LC i =A i ×PC i (11)

[0057] Wherein, A i is the area of the expropriated land, and PC i is the land acquisition cost standard per unit area;

[0058] Land reclamation cost A 26 :

[0059] Its calculation formula is:

[0060] RC i =A i ×β×R i (12)

[0061] Wherein, β represents the reclamation rate, and R i is the reclamation cost standard per unit area;

[0062] Ventilation cost of the work area A 27:

[0063] The calculation formula is:

[0064]

[0065] Among them, CV i is the total ventilation cost of the \(i\)-th work area, E Vik is the quantity of the \(k\)-th ventilation equipment in the \(i\)-th work area, U Vik is the unit time cost of the \(k\)-th ventilation equipment in the \(i\)-th work area, L i is the length of the \(i\)-th work area, \(n\) is the number of types of ventilation equipment used in this work area, β V is the coefficient of the ventilation cost increasing with the length;

[0066] The slag transportation cost \(A\) of the work area 28 :

[0067] Its calculation formula is:

[0068]

[0069] Among them, C Oi is the total slag transportation cost of the \(i\)-th work area, OQ i is the total amount of muck in the \(i\)-th work area, UC Oi is the unit muck transportation cost. L i is the length of the \(i\)-th work area, β O is the coefficient of the muck transportation cost increasing with the length;

[0070] The drainage cost \(A\) of the work area 29 :

[0071] Its calculation formula is:

[0072]

[0073] Among them, C Di is the total drainage cost of the \(i\)-th work area, E Dik is the quantity of the \(k\)-th type of drainage equipment in the \(i\)-th work area, U Dik is the unit time cost of the \(k\)-th drainage equipment in the \(i\)-th work area, L i is the length of the \(i\)-th work area, \(m\) is the number of types of drainage equipment used in this work area, β D is the coefficient of the drainage cost increasing with the length;

[0074] The safety level of the portal slope \(A\) 31 :

[0075] Its calculation formula is:

[0076] S i =α 1 ×SFi +α 2 ×G i +α 3 ×M i (16)

[0077] Among them, S i is the comprehensive safety score of the portal slope in the i-th work area, SF i is the stability coefficient of the slope, G i is the geological condition score, M i is the design measure score, α 1 ,α 2 ,α 3 , are the weight coefficients of each factor;

[0078] Risk distribution A in the auxiliary adit 32 :

[0079] Its calculation formula is:

[0080]

[0081] Among them, W j is the weight of the j-th risk, P i is the probability of the j-th risk occurring in the i-th adit, D i is the corresponding loss;

[0082] Safety degree of transporting muck out A 33 :

[0083] Its calculation formula is:

[0084]

[0085] Among them, S OTi is the comprehensive safety score of transporting muck out in the i-th work area, RL i is the safety score of the transportation route, RE i is the reliability score of the transportation equipment, θ i is the slope angle of the auxiliary adit;

[0086] Emergency rescue capacity A 34 :

[0087] Its calculation formula is:

[0088]

[0089] Among them, ER i is the emergency rescue capacity score of the i-th work area, R i is the risk distribution score of the i-th work area, L i is the length of the i-th work area, β Lis the emergency rescue coefficient related to the work area length;

[0090] Cultivated land occupation area A 41 : directly obtained from actual data;

[0091] Forest land occupation area A 42 : directly obtained from actual data;

[0092] Energy consumption of slag transported out of the work area A 43 :

[0093] Its calculation formula is:

[0094]

[0095] where E Oi is the total energy consumption of slag transported out of the i-th work area, OQ i is the total amount of muck in the i-th work area, UC E is the basic transportation energy consumption per unit of muck, θ i is the slope angle of the auxiliary tunnel, tan(θ i ) represents the incremental impact of slope on energy consumption, L i is the length of the i-th work area, β E is the coefficient of energy consumption increasing with length;

[0096] Energy consumption of ventilation in the work area A 44 :

[0097] The calculation formula is:

[0098]

[0099] where E Vi is the total ventilation energy consumption of the i-th work area, E Vik is the number of the k-th type of ventilation equipment in the i-th work area, UE Vik is the unit energy consumption of the k-th ventilation equipment in the i-th work area, L i is the length of the i-th work area, β V is the coefficient of energy consumption increasing with length;

[0100] Energy consumption of drainage in the work area A 45 :

[0101] Its calculation formula is:

[0102]

[0103] where E Di is the total drainage energy consumption of the i-th work area. E Dik is the number of the k-th type of drainage equipment in the i-th work area, UE Dik is the unit energy consumption of the k-th type of drainage equipment in the i-th work area, Li is the length of the \(i\)-th work area, and \(\beta\) D is the coefficient of energy consumption increasing with length.

[0104] Furthermore, S4 specifically includes the following steps:

[0105] S4.1. Standardize the evaluation matrix \(X\) to obtain the standardized matrix \(R\):

[0106] The specific formula is as follows:

[0107]

[0108] where \(r\) ij is the value of the \(i\)-th scheme on the \(j\)-th secondary index after standardization, \(x\) ij is the value of the original \(i\)-th scheme on the \(j\)-th secondary index, \(\min(x\) j ) is the minimum value of the \(j\)-th secondary index, \(\max(x\) j ) is the maximum value of the \(j\)-th secondary index;

[0109] S4.2. Calculate the standard deviation \(\sigma\) of each secondary index j :

[0110] The specific formula is as follows:

[0111]

[0112] where \(\sigma\) j is the standard deviation of the \(j\)-th secondary index, \(n\) is the number of schemes, is the average value of the \(j\)-th secondary index;

[0113] S4.3. Calculate the correlation coefficient \(\rho\) between secondary indices jk :

[0114] The specific formula is as follows:

[0115]

[0116] where \(\rho\) jk is the correlation coefficient between the \(j\)-th and \(k\)-th secondary indices, \(r\) ij and \(r\) ik are the values after standardization respectively;

[0117] S4.4. According to the CRITIC method, calculate the information content \(C\) of each secondary index j and determine its weight \(w\) ij :

[0118] The specific formula is as follows:

[0119]

[0120] Among them, C j is the information amount of the j-th secondary index, w ij is the weight of the j-th secondary index under the i-th primary index, and m is the number of secondary indexes;

[0121] S4.5. Aggregate the weights w ij of the secondary indexes under each primary index, and calculate the comprehensive information amount C i of the primary index:

[0122] The specific formula is as follows:

[0123]

[0124] Among them, C i is the comprehensive information amount of the i-th primary index; m i is the number of secondary indexes under the i-th primary index; w ij is the weight of the j-th secondary index under the i-th primary index;

[0125] S4.5. According to the comprehensive information amount C i , calculate the weight W i of each primary index:

[0126] The specific formula is as follows:

[0127]

[0128] Among them, W i is the weight of the i-th primary index; n is the number of primary indexes;

[0129] S4.6. Combine the weight W i of the primary index with the weight w ij of the secondary index to generate the final comprehensive weight matrix W:

[0130] The specific formula is as follows:

[0131]

[0132] Among them, W is the comprehensive weight matrix, W i is the weight of the i-th primary index, and w ij is the weight of the j-th secondary index under the i-th primary index.

[0133] Furthermore, S5 specifically includes the following steps:

[0134] S5.1. For each secondary index, determine the ideal solution and the negative ideal solution

[0135] The specific formula is as follows:

[0136]

[0137] Among them, is the ideal solution of the j-th secondary index, is the negative ideal solution of the j-th secondary index;

[0138] S5.2. Using the comprehensive weight matrix W, calculate the comprehensive score S of each scheme i and the maximum regret value R i :

[0139] The specific formula is as follows:

[0140]

[0141]

[0142] Among them, S i is the comprehensive score of the i-th scheme, and R i is the maximum regret value of the i-th scheme;

[0143] S5.3. Calculate the Q i value according to the VIKOR method, and sort the schemes according to the Q i value to obtain the optimal scheme:

[0144] The specific formula is as follows:

[0145]

[0146] Among them, Q i is the comprehensive evaluation value of the i-th scheme, v is the decision maker's preference coefficient, S * = min(S i ); S - = max(S i ), R * = min(R i ), R - = max(R i );

[0147] S5.4. Sort each alternative scheme according to the Q i value.

[0148] A computer device includes a memory and a processor. The memory stores a computer program. When the computer program is executed by the processor, the processor executes the steps of the method.

[0149] The beneficial effects of the present invention are reflected in:

[0150] The method for optimizing the setting scheme of auxiliary adits in extra-long railway tunnels of the present invention comprehensively evaluates and optimizes different auxiliary adit setting schemes from key dimensions such as construction period, cost, safety, and environmental protection through a method combining CRITIC and VIKOR multi-index evaluation;

[0151] Through a systematic evaluation index system, the present invention can accurately evaluate the advantages and disadvantages of each scheme, and based on quantitative data and scientific analysis, select the optimal auxiliary adit setting scheme, effectively overcoming the limitations of traditional methods in scheme selection. It not only improves construction efficiency and optimizes resource utilization, but also significantly enhances the safety and environmental protection of the project, providing strong technical support and decision-making basis for the setting of auxiliary adits in extra-long railway tunnels. Brief Description of the Drawings

[0152] Figure 1 It is a schematic diagram of the steps of the method for optimizing the setting scheme of auxiliary adits in extra-long railway tunnels of the present invention. Detailed Embodiments

[0153] The following will clearly and completely describe the technical solutions in the embodiments of the present invention with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only a part of the embodiments of the present invention, rather than all the embodiments. Without conflict, the embodiments in the present application and the features in the embodiments can be combined with each other. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts shall fall within the protection scope of the present invention.

[0154] See Figure 1 。

[0155] The method for optimizing the setting scheme of auxiliary adits in extra-long railway tunnels of the present invention includes the following steps:

[0156] S1. Establish an evaluation index system for auxiliary adit schemes that includes multiple primary indicators, and each primary indicator includes multiple secondary indicators;

[0157] S2. According to the geological conditions, construction requirements, and technical feasibility of the project site, formulate alternative setting schemes for auxiliary adits in extra-long railway tunnels;

[0158] S3. Based on the evaluation index system established in S1, collect data related to each secondary indicator of each auxiliary adit setting scheme to be evaluated, calculate the corresponding values, and form an evaluation matrix X;

[0159] S4. Based on matrix X, and using the CRITIC method to calculate the weights of each primary indicator and secondary indicator, and obtain a comprehensive weight matrix W;

[0160] S5. Based on the comprehensive weight matrix W, the comprehensive evaluation value of each plan is calculated using the VIKOR method and optimized.

[0161] In one embodiment, the primary indicators specifically include the construction period indicator A 1 , the cost indicator A 2 , the safety indicator A 3 and the environmental protection indicator A 4 ;

[0162] The construction period indicator A 1 includes five secondary indicators: the total project construction period A 11 , the construction period of the auxiliary adit A 12 , the construction time of the construction access road A 13 , the construction time of the auxiliary facilities A 14 , and the slag transportation time out of the work area A 15 ;

[0163] The cost indicator A 2 includes nine secondary indicators: the cost of opening a working face A 21 , the construction cost of the auxiliary adit A 22 , the construction cost of the construction access road A 23 , the construction cost of the auxiliary facilities A 24 , the cost of land acquisition A 25 , the cost of land reclamation A 26 , the ventilation cost of the work area A 27 , the slag transportation cost out of the work area A 28 , and the drainage cost of the work area A 29 ;

[0164] The safety indicator A 3 includes four secondary indicators: the safety of the portal slope A 31 , the risk distribution in the auxiliary adit A 32 , the safety level of slag transportation A 33 , and the emergency rescue ability A 34 ;

[0165] The environmental protection indicator A 4 includes five secondary indicators: the occupied area of cultivated land A 41 , the occupied area of forest land A 42 , the energy consumption of slag transportation in the work area A 43 , the energy consumption of ventilation in the work area A 44 , and the energy consumption of drainage in the work area A 45 ;

[0166] In one embodiment, the collection and calculation of each secondary indicator in S3 are as follows:

[0167] The total project construction period A 11:Through the auxiliary adit, several working faces can be opened in the main tunnel. The total project duration is the duration of the working face with the longest duration among all working faces;

[0168] Its calculation formula is:

[0169] max T i (1)

[0170] Where T i represents the duration of the i-th working face, and its calculation method is:

[0171]

[0172] Where Q ri represents the engineering quantity of the r-class surrounding rock of the i-th working face, V ri represents the construction rate per unit engineering quantity of the i-th working face under the r-class surrounding rock, and ATi represents the construction time of the auxiliary adit of the i-th working face;

[0173] Auxiliary adit construction duration A 12 : The auxiliary adit construction duration represents the construction time for building the auxiliary adit; its calculation method is:

[0174]

[0175] Where, AQ ri represents the engineering quantity of the r-class surrounding rock of the auxiliary adit of the i-th working face, AV ri represents the construction rate per unit engineering quantity of the auxiliary adit of the i-th working face under the r-class surrounding rock;

[0176] Construction time of the construction access road A 13 : Since most tunnels are located in remote mountainous areas with inconvenient transportation, it is mostly necessary to build temporary construction access roads in advance for building auxiliary adits;

[0177] Its calculation method is:

[0178]

[0179] Where RT i represents the construction time of the construction access road of the i-th working face, RQ i represents the engineering quantity of the construction access road of the i-th working face, RV i represents the construction rate per unit engineering quantity of the construction access road of the i-th working face;

[0180] Construction time of auxiliary facilities A 14 : It refers to the time required for various auxiliary facilities built to ensure the smooth progress of construction during the construction of the auxiliary adit. These facilities may include temporary water supply, power supply, ventilation, lighting systems, etc.;

[0181] Its calculation formula is:

[0182]

[0183] AF i is the construction time of the auxiliary facilities in the \(i\)-th work area, and \(W\) i is the total engineering quantity (such as earthwork volume, concrete volume, etc.) in the \(i\)-th work area, and \(C\) i is the construction capacity of the auxiliary facilities (i.e., the engineering quantity that can be completed per unit time), and \(\alpha\) is the proportionality coefficient related to the type of auxiliary facilities, which is used to reflect the input ratio required for different types of auxiliary facilities;

[0184] The slag transportation time \(A\) of the work area 15 : During the tunnel construction process, the excavated muck needs to be transported out of the work area in a timely manner to ensure the construction progress and the cleanliness and safety of the work area. The slag transportation time refers to the time required from the generation of muck to the complete removal of muck from the work area;

[0185] Its calculation formula is:

[0186]

[0187] where \(OT\) i is the slag transportation time of the \(i\)-th work area, \(OQ\) i is the total amount of muck in the \(i\)-th work area, \(OV\) i is the muck transportation rate of the \(i\)-th work area, \(L\) i is the length of the \(i\)-th work area, and \(\beta_O\) is the coefficient of the increase in slag transportation time with the increase in length; usually \(\beta\) O > 1 indicates the effect of the increase in slag transportation time with the increase in the length of the work area. As the length of the work area increases, the slag transportation time shows an increasing characteristic;

[0188] The cost \(A\) of opening a working face 21 : The cost of opening a working face refers to the comprehensive cost of labor and mechanical equipment invested at one time for opening the \(i\)-th working face during the tunnel construction process. This cost includes all personnel salaries and the rental or purchase costs of mechanical equipment for completing the initial operations;

[0189] Its calculation formula is:

[0190] C WFi = N i ×W i + M im ×MC im (7)

[0191] where, \(C\) WFi is the total cost of the \(i\)-th working face, \(N\) i is the total number of people required for the \(i\)-th working face, \(W\) iis the per capita labor cost of the i-th working face, M im is the total number of m-type mechanical equipment required for the i-th working face, MC im is the unit mechanical cost of the m-type machinery for the i-th working face;

[0192] Auxiliary adit construction cost A 22 : The auxiliary adit construction cost refers to all the expenses required for the construction of the auxiliary adit. These expenses usually include civil engineering expenses, material expenses, equipment expenses, labor expenses, and related management expenses, etc.;

[0193] Its calculation formula is:

[0194]

[0195] Among them, C APi is the total cost of the i-th auxiliary adit, AQ ri is the engineering quantity of the i-th auxiliary adit under the r-type surrounding rock grade, AU ri is the comprehensive per-meter unit price of the i-th auxiliary adit per unit engineering quantity under the r-type surrounding rock grade;

[0196] Construction access road construction cost A 23 : The construction access road construction cost refers to the comprehensive cost required for building the construction access road in the i-th work area. Calculated based on the comprehensive unit price per unit engineering quantity, it includes all expenses such as earth and rock excavation, road surface paving, and drainage facilities of the construction access road;

[0197] Its calculation formula is:

[0198] C RDi =RQ i ×RU i (9)

[0199] Among them, C RDi is the total construction cost of the construction access road in the i-th work area, RQi is the engineering quantity of the construction access road in the i-th work area, and RUi is the comprehensive unit price of the construction access road per unit engineering quantity in the i-th work area;

[0200] Ancillary facilities construction cost A 24 : Includes all costs for building the above ancillary facilities, such as material costs, labor costs, equipment costs, etc.;

[0201] Its calculation formula is:

[0202]

[0203] Among them, AC j is the unit construction cost of the j-th ancillary facility, Q j is the engineering quantity of the corresponding ancillary facility;

[0204] Land acquisition cost A 25 :

[0205] Its calculation formula is:

[0206] LC i = A i ×PC i (11)

[0207] Among them, A i is the area of the expropriated land, and PC i is the land acquisition cost standard per unit area;

[0208] Land reclamation cost A 26 : It refers to the compensation cost required for land expropriation during tunnel construction. The land acquisition cost is usually determined according to factors such as the type, use, and geographical location of the land;

[0209] Its calculation formula is:

[0210] RC i = A i ×β×R i (12)

[0211] Among them, β represents the reclamation rate R i is the reclamation cost standard per unit area;

[0212] Ventilation cost of the work area A 27 : The ventilation cost of the work area refers to the comprehensive cost of the ventilation equipment used to keep the air in the work area circulating during tunnel construction; these equipment may include multiple ventilators and related supporting facilities;

[0213] The calculation formula is:

[0214]

[0215] Among them, CV i is the total ventilation cost of the i-th work area, E Vik is the number of the k-th ventilation equipment in the i-th work area, U Vik is the unit time cost of the k-th ventilation equipment in the i-th work area, L i is the length of the i-th work area, n is the number of types of ventilation equipment used in this work area, β V is the coefficient of the ventilation cost increasing with the length; usually β V > 1 reflects the progressive influence of the work area length on the ventilation cost;

[0216] Cost of transporting slag out of the work area A 28:The cost of transporting muck out of the work area refers to the total cost required to transport the muck out of the work area during tunnel construction. As the length of the work area increases, the transportation distance of the muck also increases. Therefore, the cost of transporting muck out will increase with the increase in the length of the work area;

[0217] Its calculation formula is:

[0218]

[0219] Where C Oi is the total cost of transporting muck out of the i-th work area, OQ i is the total amount of muck in the i-th work area, UC Oi is the unit muck transportation cost. L i is the length of the i-th work area, β O is the coefficient of the increase in the muck transportation cost with the length; β O > 1 represents the increasing effect of the transportation cost with the increase in the length of the work area;

[0220] The drainage cost A of the work area 29 : The drainage cost of the work area refers to the comprehensive cost of the drainage equipment used to ensure that there is no water accumulation in the work area during tunnel construction;

[0221] Its calculation formula is:

[0222]

[0223] Where C Di is the total drainage cost of the i-th work area, E Dik is the number of the k-th type of drainage equipment in the i-th work area, U Dik is the unit time cost of the k-th drainage equipment in the i-th work area, L i is the length of the i-th work area, m is the number of types of drainage equipment used in this work area, β D is the coefficient of the increase in the drainage cost with the length; usually β D > 1 represents the progressive influence of the length of the work area on the drainage cost;

[0224] The safety level of the portal slope A 31 : Evaluate the stability and safety of the portal slope of the auxiliary adit. Usually, a comprehensive evaluation is carried out through the stability coefficient of the slope, geological conditions, support measures, etc.;

[0225] Its calculation formula is:

[0226] S i =α 1 ×SF i +α 2 ×G i +α 3 ×M i (16)

[0227] Among them, S i is the comprehensive safety score of the portal slope in the i-th work area, and SF i is the stability coefficient of the slope, G i is the geological condition score, M i is the design measure score, α 1 , α 2 , α 3 , are the weight coefficients of each factor; these weights can be set according to the actual situation and expert opinions, and satisfy α 1 +α 2 +α 3 +α 4 = 1;

[0228] Risk distribution A in the auxiliary adit 32 : Evaluate the possible risks in the auxiliary adit during the construction process, including risks such as collapse, water inrush, gas, etc. The distribution can be described by a multi-factor evaluation model.

[0229] Its calculation formula is:

[0230]

[0231] Among them, W j is the weight of the j-th risk, P i is the probability of the j-th risk occurring in the i-th adit, and D i is the corresponding loss;

[0232] Safety degree of mucking out A 33 : The safety degree of mucking out refers to the safety when the muck is transported out of the work area during tunnel construction. The factors affecting the safety of muck transportation include the slope of the transportation route, the reliability of the transportation equipment, and the slope of the auxiliary adit, etc.;

[0233] Its calculation formula is:

[0234]

[0235] Among them, S OTi is the comprehensive safety score of mucking out in the i-th work area, RL i is the safety score of the transportation route (such as the smoothness and curvature of the route, scored by experts), RE i is the reliability score of the transportation equipment (such as the stability and failure rate of the equipment, scored by experts), θ i is the slope angle (in radians) of the auxiliary adit, indicating the impact of the slope on safety. The greater the slope, the lower the safety score;

[0236] Emergency rescue capacity A 34Emergency rescue capacity refers to the emergency response ability during tunnel construction in case of emergencies. Emergency rescue capacity is not only affected by the length of the work area but also associated with the risk distribution in the auxiliary adit. The more extensive the risk distribution and the longer the work area, the stronger the emergency rescue capacity required to ensure safety during construction;

[0237] Its calculation formula is:

[0238]

[0239] Where ER i is the emergency rescue capacity score of the i-th work area (the larger the value, the stronger the emergency rescue capacity), R i is the risk distribution score of the i-th work area (calculated from A 32 , the larger the value, the greater the risk), L i is the length of the i-th work area, β L is the emergency rescue coefficient related to the work area length, β L > 0, S i is the effectiveness score of emergency equipment and measures (such as the availability of rescue equipment and rescue teams, scored by experts);

[0240] Cultivated land occupation area A 41 Cultivated land occupation area refers to the cultivated land area temporarily or permanently occupied during construction, directly obtained from actual data;

[0241] Forest land occupation area A 42 Forest land occupation area refers to the forest land area temporarily or permanently occupied during construction, directly obtained from actual data;

[0242] Energy consumption for transporting muck out of the work area A 43 Energy consumption for transporting muck out of the work area refers to the energy consumed for transporting muck out of the work area during tunnel construction. The energy consumption for muck transportation is not only related to the length of the work area but also to the slope of the auxiliary adit. The greater the slope, the higher the transportation difficulty and energy consumption;

[0243] Its calculation formula is:

[0244]

[0245] Where E Oi is the total energy consumption for transporting muck out of the i-th work area, OQ i is the total amount of muck in the i-th work area, UC E is the basic transportation energy consumption per unit of muck, θ i is the slope angle of the auxiliary adit, tan(θ i ) represents the incremental impact of the slope on energy consumption, L i is the length of the i-th work area, βE is the coefficient of energy consumption increasing with the length; usually β E > 1 indicates the increasing effect of energy consumption with the increase of the work area length, reflecting the cumulative impact of longer transportation distance on energy consumption;

[0246] The ventilation energy consumption A of the work area 44 : The ventilation energy consumption of the work area refers to the energy consumption of the ventilation equipment used to keep the air flowing in the work area during tunnel construction. As the length of the work area increases, the required ventilation equipment and operation time increase;

[0247] The calculation formula is:

[0248]

[0249] Among them, E Vi is the total ventilation energy consumption of the i-th work area, E Vik is the number of the k-th type of ventilation equipment in the i-th work area, UE Vik is the unit energy consumption of the k-th ventilation equipment in the i-th work area, L i is the length of the i-th work area, β V is the coefficient of energy consumption increasing with the length; usually β V > 1, reflecting the progressive impact of the work area length on ventilation energy consumption;

[0250] The drainage energy consumption A of the work area 45 : The drainage energy consumption of the work area refers to the energy consumption of the drainage equipment used to ensure that there is no water accumulation in the work area during tunnel construction. As the length of the work area increases, the required drainage equipment and operation time increase;

[0251] Its calculation formula is:

[0252]

[0253] Among them, E Di is the total drainage energy consumption of the i-th work area. E Dik is the number of the k-th type of drainage equipment in the i-th work area, UE Dik is the unit energy consumption of the k-th type of drainage equipment in the i-th work area, L i is the length of the i-th work area, β D is the coefficient of energy consumption increasing with the length.; usually β D > 1, reflecting the progressive impact of the work area length on drainage energy consumption.

[0254] In one embodiment, S4 specifically includes the following steps:

[0255] S4.1. Standardize the evaluation matrix X to eliminate the influence between different index dimensions and obtain the standardized matrix R:

[0256] The specific formula is as follows:

[0257]

[0258] Among them, r ij is the value of the i-th solution on the j-th secondary index after standardization, and x ij is the value of the original i-th solution on the j-th secondary index. min(x j ) is the minimum value of the j-th secondary index, and max(x j ) is the maximum value of the j-th secondary index;

[0259] S4.2. Calculate the standard deviation σ j of each secondary index, reflecting the differences of each index in different solutions:

[0260] The specific formula is as follows:

[0261]

[0262] Among them, σ j is the standard deviation of the j-th secondary index, n is the number of solutions, is the average value of the j-th secondary index;

[0263] S4.3. Calculate the correlation coefficient ρ jk between secondary indexes, measuring the similarity degree between different secondary indexes:

[0264] The specific formula is as follows:

[0265]

[0266] Among them, ρ jk is the correlation coefficient between the j-th and k-th secondary indexes, and r ij and r ik are the values after standardization respectively;

[0267] S4.4. According to the CRITIC method, calculate the information amount C j of each secondary index and determine its weight w ij :

[0268] The specific formula is as follows:

[0269]

[0270] Among them, C j is the information amount of the j-th secondary index, and w ij is the weight of the j-th secondary index under the i-th primary index, and m is the number of secondary indexes;

[0271] S4.5. For the secondary index weights w under each primary indexij Aggregate and calculate the comprehensive information quantity C of the first-level indicators i :

[0272] The specific formula is as follows:

[0273]

[0274] Among them, C i is the comprehensive information quantity of the i-th first-level indicator; m i is the number of second-level indicators under the i-th first-level indicator; w ij is the weight of the j-th second-level indicator under the i-th first-level indicator;

[0275] S4.5. Calculate the weight W of each first-level indicator according to the comprehensive information quantity C i : i :

[0276] The specific formula is as follows:

[0277]

[0278] Among them, W i is the weight of the i-th first-level indicator; n is the number of first-level indicators;

[0279] S4.6. Combine the weight W of the first-level indicator with the weight w of the second-level indicator i to generate the final comprehensive weight matrix W: ij The specific formula is as follows:

[0280] The specific formula is as follows:

[0281]

[0282] Among them, W is the comprehensive weight matrix, W i is the weight of the i-th first-level indicator, w ij is the weight of the j-th second-level indicator under the i-th first-level indicator.

[0283] In one embodiment, S5 specifically includes the following steps:

[0284] S5.1. For each second-level indicator, determine the ideal solution and the negative ideal solution

[0285] The specific formula is as follows:

[0286]

[0287] Among them, is the ideal solution of the j-th second-level indicator, is the negative ideal solution of the j-th second-level indicator;

[0288] S5.2. Calculate the comprehensive score S of each scheme by using the comprehensive weight matrix W i and the maximum regret value R i :

[0289] The specific formula is as follows:

[0290]

[0291] where S i is the comprehensive score of the i-th scheme, and R i is the maximum regret value of the i-th scheme;

[0292] S5.3. Calculate the Q i value according to the VIKOR method, and sort the schemes according to the Q i value to obtain the optimal scheme:

[0293] The specific formula is as follows:

[0294]

[0295] where Q i is the comprehensive evaluation value of the i-th scheme, v is the decision maker's preference coefficient, usually taking 0.5; S * = min(S i ); S - = max(S i ), R * = min(R i ), R - = max(R i );

[0296] S5.4. Sort the alternative schemes according to the Q i value; the lower the Q i value, the higher the superiority of the spoil ground setting scheme for the extra-long railway tunnel.

[0297] A computer device includes a memory and a processor. The memory stores a computer program. When the computer program is executed by the processor, the processor executes the steps of the method.

[0298] The present invention will be described below with reference to application examples.

[0299] The starting and ending mileage of a certain tunnel is DK70+026.1~DK86+300, with a total length of 16,273.9m. The first 415m at the tunnel entrance is a double-track tunnel, and the rest are single-track tunnels. The tunnel has a one-way downhill slope of -5‰~-16.7‰. Except for the section from DK83+559.22 to DK85+665.65, which is located on a left-hand curve with a radius of R = 2000m, the rest are on straight lines. Both the tunnel entrance and exit are close to local roads, with relatively convenient transportation. There are currently three auxiliary adit setting schemes as shown in Table 1.

[0300] Table 1 Auxiliary Adit Setting Schemes for a Certain Extra-long Railway Tunnel

[0301]

[0302]

[0303] For Scheme 1, the entrance of the 1# inclined shaft is located in a mountain gully. There are many newly built access roads, and most of the land occupied is forest land, making land acquisition difficult and affecting the construction of the inclined shaft entrance.

[0304] For Scheme 2, the entrance of the 1# inclined shaft is in a densely wooded area with a deep gully at the entrance and a large height difference.

[0305] For Scheme 3, the entrance of the 1# inclined shaft has a wide field of vision, a flat site, and convenient construction conditions.

[0306] For Scheme 1, the entrance of the 2# inclined shaft is located in a mountain gully. The construction road needs to pass through a village. The original construction access road cannot meet the traffic requirements and 800m of new road needs to be built and concretely hardened, resulting in great obstruction to later construction.

[0307] For Scheme 2, the entrance of the 2# inclined shaft is close to the existing rural road, and the entrance site is open. After on-site investigation, the existing concrete road between the national road and the entrance cannot meet the traffic requirements, and the access road is built along the gully, making it difficult to widen the road.

[0308] For Scheme 3, the 2# inclined shaft is located beside the national road, with convenient transportation and an open site.

[0309] The implementation steps of a comprehensive optimization method for the auxiliary adit setting scheme of an extra-long railway tunnel are as follows:

[0310] 1. Data collection and index classification

[0311] In an example of the present invention, for the extra-long railway tunnel project, the following four first-level indicators are first determined, and each second-level indicator is classified under the corresponding first-level indicator.

[0312] 2. In the embodiment of the present invention, through the calculations of formulas (1)-(22) for the collected data, all index values are obtained, as shown in Table 2.

[0313] Table 2 Scheme Index Values

[0314]

[0315]

[0316] 3. Standardize the data to eliminate the dimensional differences between different indicators. After standardization, a standardized matrix R is formed,

[0317] where each element r ij represents the standardized value of the i-th solution on the j-th secondary indicator.

[0318] Next, calculate the weights of each secondary indicator using the CRITIC method. The specific steps are as follows:

[0319] (1) Calculate the standard deviation σ j : Reflects the differences in secondary indicators among different solutions.

[0320] (2) Calculate the correlation coefficient ρ jk : Measures the correlation between secondary indicators.

[0321] (3) Calculate the information content C j and determine the weight w of the secondary indicator ij， as shown in Table 3.

[0322] Table 3 Weights of secondary indicators

[0323]

[0324]

[0325] Subsequently, calculate the weights of the primary indicators. Aggregate the information content of the secondary indicators under each primary indicator to calculate the comprehensive information content C i , and thereby determine the weight W of each primary indicator i , as shown in Table 4.

[0326] Table 4 Weights of primary indicators

[0327]

[0328] Finally, combine the weights of the primary and secondary indicators to form the final comprehensive weight matrix W.

[0329] 4. Comprehensive evaluation and optimization using the VIKOR method

[0330] In the embodiments of the present invention, the VIKOR method is used to conduct a comprehensive evaluation of each solution. The specific implementation steps are as follows:

[0331] (1) Determine the ideal solution and the negative ideal solution: For each secondary index, determine the optimal value (ideal solution) and the worst value (negative ideal solution).

[0332] (2) Calculate the comprehensive score and the maximum regret value:

[0333] The comprehensive score S i : It is obtained by weighted calculation of all secondary indexes and reflects the overall performance of each scheme.

[0334] The maximum regret value R i : It reflects the weighted score of the worst performance of each scheme in all secondary indexes.

[0335] (4) Calculate the Q i value and sort: Combining the comprehensive score and the maximum regret value, calculate the Q i value of each scheme by the VIKOR method, and sort the schemes according to the Q i value. The scheme with the lowest Q i value is the optimal scheme and is recommended for the auxiliary adit setting in actual projects. The results are shown in Table 5.

[0336] Table 5 Evaluation results

[0337]

[0338] Finally, it can be seen from the results of the embodiments that Scheme Three is the optimal auxiliary adit setting scheme.

[0339] In summary, the present invention comprehensively and scientifically determines the weights of the primary and secondary indexes through the CRITIC method, and combines the VIKOR method for comprehensive evaluation, so as to scientifically and effectively optimize the auxiliary adit setting scheme of extra-long railway tunnels. This method can comprehensively consider factors such as construction period, cost, safety and environmental protection during the construction process, and provide solid decision-making support for tunnel construction.

[0340] It should be noted that if there are directional indications such as up, down, left, right, front, back... in the embodiments of the present invention, the directional indications are only used to explain the relative positional relationship and movement conditions between components in a certain specific posture as shown in the drawings. If the specific posture changes, the directional indications will also change accordingly.

[0341] In addition, if there are descriptions involving "first", "second", etc. in the embodiments of the present invention, such descriptions of "first", "second", etc. are for descriptive purposes only and should not be construed as indicating or implying their relative importance or implicitly specifying the quantity of the indicated technical features. Thus, features defined with "first" and "second" may explicitly or implicitly include at least one such feature. In addition, the meaning of "and / or" appearing throughout the text includes three parallel scenarios. Taking "A and / or B" as an example, it includes Scenario A, or Scenario B, or the scenario where both A and B are satisfied simultaneously. In addition, "a plurality of" means two or more. In addition, the technical solutions between various embodiments may be combined with each other, but it must be based on the ability of those of ordinary skill in the art to implement. When the combination of technical solutions results in contradictions or cannot be implemented, it should be considered that such a combination of technical solutions does not exist and is not within the scope of protection required by the present invention.

[0342] It should be understood that the examples and embodiments described herein are for illustrative purposes only and are not intended to limit the present invention. Those skilled in the art can make various modifications or changes based on it. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principles of the present invention shall be included within the scope of protection of the present invention.

Claims

1. A method for optimizing the auxiliary tunnel setting scheme for an extra-long railway tunnel, characterized in that: The following steps are involved: S1. Establish an auxiliary tunnel scheme evaluation index system that includes multiple first-level indicators, and each first-level indicator includes multiple second-level indicators; S2. Formulate alternative plans for the setting up of auxiliary tunnels for extra-long railway tunnels based on the geological conditions, construction requirements and technical feasibility of the project site; S3. Based on the evaluation index system established in S1, collect data related to each secondary index of each auxiliary tunnel setting scheme to be evaluated, calculate the corresponding values, and form an evaluation matrix X; S4. Based on the matrix X, the CRITIC method is used to calculate the weights of the first-level indicators and the second-level indicators to obtain the comprehensive weight matrix W; S5. Based on the comprehensive weight matrix W, the VIKOR method is used to calculate the comprehensive evaluation value of each scheme and optimize it.

2. The method for optimizing the auxiliary tunnel setting scheme for a super-long railway tunnel according to claim 1 is characterized in that: The first-level indicators specifically include construction period indicator A1, cost indicator A2, safety indicator A3 and environmental protection indicator A4; The construction period indicator A1 includes the total construction period of the project A 11 , Auxiliary tunnel construction period A 12 , Construction access road construction time A 13 、Construction time of auxiliary facilities A 14 , Slag removal time in the work area A 15 Five secondary indicators; Cost indicator A2 includes the cost of opening a working face A 21 , Auxiliary tunnel construction cost A 22 , Construction access road construction cost A 23 、Cost of construction of auxiliary facilities A 24 , Land acquisition cost A 25 , Land reclamation cost A 26 、Work area ventilation cost A 27 、Slag transportation cost in the work area A 28 、Work area drainage cost A 29 Nine secondary indicators; Safety index A3 includes the safety index of the tunnel entrance slope A 31 、Risk distribution in auxiliary tunnel A 32 、Safety level of slag transport A 33 、Emergency rescue capability A 34 Four secondary indicators; Environmental protection index A4 includes cultivated land area A 41 、Forest land area A 42 、Energy consumption of slag transport in the work area A 43 , Work area ventilation energy consumption A 44 、Energy consumption of drainage in work area A 45 Five secondary indicators.

3. The method for optimizing the auxiliary tunnel arrangement scheme for a super-long railway tunnel according to claim 2 is characterized in that: The collection and calculation of various secondary indicators in S3 are as follows: Total project duration 11 : The calculation formula is: max T i (1) Where T i represents the duration of the i-th working surface, which is calculated as follows: Where Q ri V represents the engineering volume of the r-type surrounding rock at the i-th working face, ri represents the unit construction rate of the i-th working face under the r-type surrounding rock, and ATi represents the construction time of the auxiliary tunnel of the i-th working face; Auxiliary tunnel construction period A 12 : The calculation method is: Among them, AQ ri represents the engineering volume of the auxiliary tunnel of the i-th working face with type r surrounding rock, AV ri It represents the construction rate per unit engineering quantity of the auxiliary tunnel of the i-th working face under the surrounding rock of type r; Construction time of construction access road 13 : The calculation method is: RT i represents the construction time of the construction access road of the i-th working face, RQ i represents the construction access road quantity of the i-th working face, RV i It represents the construction rate per unit of construction work for the construction access road of the i-th working face; Construction time of auxiliary facilities 14 : The calculation formula is: AF i is the construction time of the auxiliary facilities in the i-th work area, W i is the total engineering quantity of the i-th work area, C i is the construction capacity of ancillary facilities, α is the proportional coefficient related to the type of ancillary facilities, which is used to reflect the investment ratio required for different types of ancillary facilities; Slag removal time in the work area A 15 : The calculation formula is: Among them, OT i is the slag removal time of the i-th work area, OQ i is the total amount of slag in the i-th work area, OV i is the soil transportation rate of the i-th construction area, L i is the length of the i-th work area, βO is the coefficient of the time for transporting the slag increasing with the length; Cost of opening working face A 21 : The calculation formula is: C WFi =N i ×W i+ M im ×MC im (7) Among them, C WFi is the total cost of the i-th working surface, N i is the total number of people required for the i-th work surface, W i is the average labor cost per person at the i-th work surface, M im is the total number of m-type mechanical equipment required for the i-th working surface, MC im is the unit machinery cost of type m machinery at the i-th working surface; Auxiliary tunnel construction cost A 22 : The calculation formula is: Among them, C APi is the total cost of the ith auxiliary tunnel, AQ ri is the engineering quantity of the ith auxiliary tunnel under the surrounding rock grade of type r, AU ri is the comprehensive unit price per meter of the i-th auxiliary tunnel unit engineering quantity under the r-type surrounding rock grade; Construction access road construction cost A 23 : The calculation formula is: C RDi =RQ i ×RU i (9) Among them, C RDi is the total cost of construction access road construction in the ith work area, RQi is the construction access road quantity in the ith work area, and RUi is the comprehensive unit price of the construction access road unit quantity in the ith work area; Cost of construction of auxiliary facilities A 24 : The calculation formula is: Among them, AC j is the unit construction cost of the jth ancillary facility, Q j The engineering volume of the corresponding ancillary facilities; Land acquisition cost A 25 : The calculation formula is: LC i =A i ×PC i (11) Among them, A i The land area to be acquired, PC i The standard of land acquisition fee per unit area; Land reclamation cost A 26 : The calculation formula is: RC i =A i ×β×R i (12) Among them, β represents the reclamation rate R i is the reclamation cost standard per unit area; Ventilation cost of work area A 27 : The calculation formula is: Among them, CV i is the total ventilation cost of the i-th work area, E Vik is the number of the kth type of ventilation equipment in the ith work area, U Vik is the unit time cost of the kth ventilation equipment in the ith work area, L i is the length of the i-th work area, n is the number of types of ventilation equipment used in the work area, β V is the coefficient of ventilation cost increasing with length; Cost of slag transport from the work area A 28 : The calculation formula is: Among them, C Oi is the total cost of slag transportation in the i-th work area, OQ i is the total amount of slag in the i-th work area, UC Oi is the unit slag transportation cost. i is the length of the i-th work area, β O is the coefficient of the cost of transporting the muck increasing with the length; Work area drainage cost A 29 : The calculation formula is: Among them, C Di is the total drainage cost of the i-th work area, E Dik is the number of the kth type of drainage equipment in the ith work area, U Dik is the unit time cost of the kth drainage equipment in the ith work area, L i is the length of the i-th work area, m is the number of types of drainage equipment used in the work area, β D is the coefficient of drainage cost increasing with length; Safety level of tunnel entrance slope A 31 : The calculation formula is: S i =α1×SF i +α2×G i +α3×M i (16) Among them, S i is the comprehensive safety score of the tunnel slope in the i-th work area, SF i is the slope stability factor, G i Score the geological conditions, M i is the design measure score, α1, α2, α3, are the weight coefficients of each factor; Risk distribution in auxiliary tunnel A 32 : The calculation formula is: Among them, W j is the weight of the j-th risk, P i is the probability of the jth risk occurring in the i-th tunnel, D i For the corresponding loss; Safety level of slag transport A 33 : The calculation formula is: Among them, S OTi is the comprehensive safety score of slag transport in the i-th work area, RL i Score the safety of the transport route, RE i is the reliability score of the transport equipment, θ i is the slope angle of the auxiliary tunnel; Emergency rescue capability A 34 : The calculation formula is: Among them, ER i is the emergency rescue capability score of the i-th work area, R i is the risk distribution score of the i-th work area, L i is the length of the i-th work area, β L is the emergency rescue coefficient related to the length of the work area; Area of ​​cultivated land A 41 : Derived directly from actual data; Forest land area A 42 : Derived directly from actual data; Energy consumption of slag transport in the work area A 43 : The calculation formula is: Among them, E Oi is the total energy consumption of slag transport in the i-th work area, OQ i is the total amount of slag in the i-th work area, UC E is the basic transportation energy consumption per unit of slag, θ i is the slope angle of the auxiliary tunnel, tan(θ i ) represents the incremental effect of slope on energy consumption, L i is the length of the i-th work area, β E is the coefficient of energy consumption increasing with length; Energy consumption of ventilation in work area A 44 : The calculation formula is: Among them, E Vi is the total ventilation energy consumption of the i-th work area, E Vik is the number of the kth type of ventilation equipment in the ith work area, UE Vik is the unit energy consumption of the kth ventilation equipment in the ith work area, L i is the length of the i-th work area, β V is the coefficient of energy consumption increasing with length; Energy consumption of drainage in work area A 45 : The calculation formula is: Among them, E Di is the total energy consumption of drainage in the i-th work area. Dik is the number of the kth type of drainage equipment in the ith work area, UE Dik is the unit energy consumption of the kth type of drainage equipment in the ith work area, L i is the length of the i-th work area, β D is the coefficient by which energy consumption increases with length.

4. The method for optimizing the auxiliary tunnel arrangement scheme for a super-long railway tunnel according to claim 1 or 2, characterized in that: S4 specifically includes the following steps: S4.

1. Standardize the evaluation matrix X to obtain the standardized matrix R: The specific formula is as follows: Among them, r ij is the value of the i-th solution on the j-th secondary indicator after standardization, x ij is the value of the original i-th solution on the j-th secondary index, min(x j ) is the minimum value of the j-th secondary index, max(x j ) is the maximum value of the j-th secondary index; S4.

2. Calculate the standard deviation σ of each secondary indicator j : The specific formula is as follows: Among them, σ j is the standard deviation of the jth secondary index, n is the number of solutions, is the average value of the jth secondary index; S4.

3. Calculate the correlation coefficient ρ between the secondary indicators jk : The specific formula is as follows: Among them, ρ jk is the correlation coefficient between the jth and kth secondary indicators, r ij and r ik are the standardized values ​​respectively; S4.

4. Calculate the information content C of each secondary indicator according to the CRITIC method j And determine its weight w ij : The specific formula is as follows: Among them, C j is the information content of the jth secondary index, w ij is the weight of the jth secondary indicator under the ith first-level indicator, and m is the number of secondary indicators; S4.

5. The weight of the secondary indicator under each primary indicator w ij Aggregation, calculate the comprehensive information volume C of the first-level index i : The specific formula is as follows: Among them, C i is the comprehensive information of the i-th first-level indicator; m i is the number of secondary indicators under the i-th primary indicator; w ij is the weight of the jth secondary indicator under the i-th primary indicator; S4.

5. Based on the comprehensive information C i , calculate the weight W of each first-level indicator i : The specific formula is as follows: Among them, W i is the weight of the i-th first-level indicator; n is the number of first-level indicators; S4.

6. Set the first-level indicator weight W i and the secondary index weight w ij Combined, the final comprehensive weight matrix W is generated: The specific formula is as follows: Among them, W is the comprehensive weight matrix, W i is the weight of the i-th first-level indicator, w ij is the weight of the jth secondary indicator under the i-th primary indicator.

5. The method for optimizing the auxiliary tunnel arrangement scheme for a super-long railway tunnel according to claim 4 is characterized in that: S5 specifically includes the following steps: S5.

1. For each secondary indicator, determine the ideal solution and negative ideal solution The specific formula is as follows: Among them, f j * is the ideal solution for the jth secondary index, is the negative ideal solution of the j-th secondary index; S5.

2. Use the comprehensive weight matrix W to calculate the comprehensive score S of each solution i and the maximum regret value R i : The specific formula is as follows: Among them, S i is the comprehensive score of the ith solution, R i is the maximum regret value of the i-th solution; S5.

3. Calculation of Q according to the VIKOR method i value, and according to Q i Sort the solutions by value and get the optimal solution: The specific formula is as follows: Among them, Q i is the comprehensive evaluation value of the i-th option, v is the decision maker’s preference coefficient, S * =min(S i );S - =max(S i ), R * =min(R i ), R - =max(R i ); S5.

4. According to Q i value, and rank the alternatives.

6. A computer device comprising a memory and a processor, wherein the memory stores a computer program, and when the computer program is executed by the processor, the processor executes the steps of the method according to any one of claims 1 to 5.

Images