# The Prediction of Metro Shield Construction Cost Based on a Backpropagation Neural Network Improved by Quantum Particle Swarm Optimization

Authors: Lanjun Liu; Denghui Liu; Han Wu; Xinyu Wang
Journal: Advances in Civil Engineering
 (2020)

---

## Abstract

The prediction of construction cost of metro shield engineering is of great significance to project management. In this study, we used the rough set theory, a backpropagation (BP) neural network, and quantum particle swarm optimization (QPSO) to establish a prediction model for predicting the metro shield construction costs. The model accounts for the complexity of metro shield construction and the nonlinear relationship between the construction cost factors. First, the factors affecting the construction cost were determined by referring to the Chinese National Standards and analysing the engineering practice of typical metro shield projects. The rough set theory was used to simplify the system of influencing factors to extract the dominant influencing factors and reduce the number of input variables in the BP neural network. Since the BP neural network easily falls into a local minimum and has a slow convergence speed, QPSO was used to optimize the weights and thresholds of the BP neural network. This method combined the strong nonlinear analysis capabilities of the BP and the global search capabilities of the QPSO. Finally, we selected 50 projects in China for a case analysis. The results showed the dominant factors affecting the construction cost of these projects included ten indicators, such as the type of tunnelling machine and the geological characteristics. The determination coefficient, mean absolute percentage error, root mean square error, and mean absolute error, which are frequently used error analysis tools, were used to analyse the calculation errors of different models (the proposed model, a multiple regression method, a traditional BP model, a BP model optimized by the genetic algorithm, and the BP model optimized by the particle swarm optimization). The results showed that the proposed method had the highest prediction accuracy and stability, demonstrating the effectiveness and excellent performance of this proposed method.


---

## Body

## 1. Introduction

With the characteristics of energy-saving, land-saving, large transportation volume, low pollution, and low-operational risk, the metro plays an increasingly important role in the urban public transportation system [1]. Nowadays, developing countries, represented by China, are at the peak of urban rail transit construction. According to the Thirteenth Five-Year Plan issued by the Chinese government, by 2020, the total length of the Chinese metro will reach 8,600 km, and the investment in metro construction will reach 300 billion dollars. The construction of metro projects includes the tunnel and metro stations. The cost of tunnel construction often accounts for more than 40% of the total construction costs [2]. Shield construction technology has the advantages of a small impact on the surrounding environment, fast construction speed, excellent engineering quality, construction safety, environmental protection, and wide application range. Shield construction is the most commonly used method of urban metro tunnel construction [3]. Improving the accuracy and reliability of the cost prediction of shield construction ensures the smooth execution of metro construction projects and protects the economic interests of the owners, construction companies, and the public.There are substantial differences in the definition of construction costs between different countries due to national regulations and different construction contracts. In this study, we selected the most authoritative definition of the contract unit price, which has been recommended by the International Federation of Consulting Engineers (FIDIC) and is accepted by the World Bank, to define the construction cost of the metro shield. The selection of the contract unit price also facilitates follow-up research. The total cost of metro shield construction comprises all costs paid by the owner during the metro shield construction, including direct costs, indirect costs, headquarters management fees, profits, taxes, tentative deposits, changes approved by engineers during the implementation process, and the amount claimed by the contractor [4]. The final cost of metro shield projects depends on the length of the project, but the unit price of the construction cost per kilometre (km) is comparable. Therefore, the construction cost analysed in this article is the unit price of metro shield construction per km.According to the above definition, the components of the metro shield construction cost are relatively complex, and the influencing factors include the construction details, construction environment, project management, capital structure, and many other factors. Combined with engineering practice, it is evident that the influence of these factors on the construction cost of the metro shield is relatively complicated. On the one hand, the construction process has high degrees of mechanization and standardization; thus, the construction costs are related linearly to the project size [5]. On the other hand, many nonlinear and qualitative factors, such as hydrological conditions, geological conditions, environmental characteristics, and project management capabilities, have a nonlinear influence on the metro construction costs [6].Many scholars have conducted in-depth research on construction cost prediction. Wang et al. [7] used an activity-based costing method to predict the supply chain cost of prefabricated building construction in an uncertain environment. Although the model could be used to determine the key factors to reduce the cost of the entire supply chain, the subjectivity and arbitrariness of the choice of the cost drivers led to poor cost forecasting performance. Chen [8] used multiple regression analysis to predict the cost of four construction projects. However, it was assumed that the relationship between the independent variable and the dependent variable was linear, affecting the accuracy of the prediction results. Tabei et al. [9] believed that the main factor leading to the large prediction error was the uncertainty of the factors, and a new method was proposed to predict the project cash flow in a fuzzy environment. However, this method requires a large amount of comprehensive statistical data to obtain consistent research results, adversely affecting the applicability of this method. Cheng et al. [10] used the golden section method, and the grey forecast method improved by dichotomy to establish a forecast model. The method provided good prediction accuracy for the case study, but the grey forecast method requires uniformity of the original data. In addition, the above methods either rely significantly on professional knowledge or are only suitable for specific situations; thus, they do not have modelling and analysis flexibility.The emergence of artificial intelligence algorithms, such as artificial neural networks (ANNs), provides a new approach for solving the complex problems of predicting the metro shield construction cost. ANNs can approach any nonlinear continuous function with arbitrary precision and are well suited for solving the highly nonlinear problem of cost prediction [11]. ANNs have strong self-learning and self-adaptive abilities. They can extract information from input and output data and store it in the weights of the network; thus, they possess good generalization ability [12]. In addition, ANNs are able to use a limited amount of data obtained in the field, rather than using large amounts of statistical data so that a data-driven method can be used to predict the construction costs.A backpropagation (BP) neural network is a multilayer feedforward network trained by an error BP algorithm, and it is the most widely used type of ANN [11]. However, traditional BP neural networks may have the following shortcomings when applied to the analysis of complex systems. (1) The traditional BP neural network algorithm is a local search optimization method, which easily falls into local extremes, and the weights converge to the local extremes, a small sample size resulting in network training failure [13]. (2) The BP neural network is very sensitive to the initial network weights. Initializing the network with different weights tends to converge to different local minima [14]. Therefore, a high-precision model has to be developed to improve the traditional BP neural network.Some scholars have used a genetic algorithm (GA) and particle swarm optimization (PSO) to improve the traditional BP model. The initial weights and thresholds of the BP neural network are optimized by these metaheuristic algorithms to achieve higher accuracy than other methods. Wen et al. [15] established a construction cost prediction model for a power transmission and transformation project using a BP model optimized by a GA. However, the GA algorithm had some problems, such as depending on the initial population selection, a slow convergence speed, and many parameter settings, resulting in relatively low accuracy of the model. Liu and Chi [16] used the PSO to optimize the weights and thresholds and establish a cost estimation model for a government investment project. Compared with classical metaheuristic algorithms, such as a GA, the PSO algorithm has the advantages of simple operation, fast convergence speed, and few parameter settings [17]. However, PSO has a limited search space and easily falls into a local extremum [18, 19]. The reason is that, in the classical PSO, the particles have a tendency to orbit rather than converging to zero, and the particles have a maximum speed limit during the search; thus, the search area is limited in each search. Sun et al. [18] proposed the novel quantum particle swarm optimization algorithm (QPSO), which is based on the PSO algorithm. During particle convergence, the particle keeps approaching the point until it falls inside. From the perspective of quantum mechanics [19], an attractive potential exists at the midpoint of the convergence process, attracting the particles to converge so that the whole group is aggregated. Particles in the quantum bound state can appear at any space point with a certain probability, which ensures good convergence performance and global search ability of the QPSO. Xu et al. [20] used the QPSO for the optimization of synthetic and actual hyperspectral datasets. The experimental results showed that the QPSO method obtained better results than vertex component analysis (VCA) and the standard PSO in remote sensing. Sun et al. [21] solved the problem of designing an optimal two-dimensional zero-phase infinite impulse response (IIR) digital filter using QPSO. Before analysing unmanned aerial vehicle (UAV) route planning, Fu et al. [22] tested various optimization algorithms, such as GA, differential evolution (DE), the standard PSO, the phase-coded PSO, the QPSO, and other optimization algorithms. The results showed that the computational power of QPSO was better than that of GA, DE, and classic PSO algorithms. Therefore, in this study, we used the QPSO to optimize the BP model.Many factors affect the construction cost of metro shields. If all factors are included in the input of the BP model, the computational burden for training the neural network will be high, and the prediction accuracy will be reduced. In the selection of input variables for the BP model, most studies used threshold analysis with an arbitrary threshold to select the input variables for the BP model [23]. This selection method is affected by subjective factors, causing information redundancy or loss. The accuracy of neural network prediction can be increased by using a suitable method for selecting the important factors affecting the construction cost of metro shields as model inputs.Rough set theory is a mathematical analysis tool to describe the completeness and uncertainty of data [24]. The method simplifies the data, obtains the minimum expression of knowledge while retaining key information, and identifies and evaluates the dependence between the data. Shiau et al. [25] successfully used rough set theory to screen factors affecting transportation sustainability in Taiwan. Cheng et al. [26] used the rough set theory to select and analyse 21 indicators affecting regional green competitiveness and formulate regional green competitiveness indicators. Ge et al. [27] extracted spatial relationship indicators using rough set theory to describe the complex relationship between geographical phenomena and the environment. Many factors affect the cost of subway shield construction. If we regard them as a knowledge expression system and use a fuzzy rough set method to analyse various factors that affect the cost of subway shield construction, we can reduce the number of variables and obtain many important attributes influencing the construction cost of a subway shield. These factors can be used as input to the neural network to optimize the input variables, and the importance of each attribute in the cost prediction of subway shield construction can be determined.Against this background, we established a prediction model based on a BP neural network improved by the QPSO. This paper provided the following contributions. (1) This paper determined the factors affecting the construction cost and provided details on the data acquisition or calculation method for each factor. We used the rough set theory to select the main influencing factors affecting the construction cost of 50 metro shield projects in a case analysis. This approach not only handled scientifically the complexity of metro shield construction and the numerous factors affecting construction costs but also substantially reduced the input variables of the BP neural network prediction model. (2) This paper used the BP algorithm optimized by the QPSO to establish a cost prediction model for metro shield construction. This method combined the global convergence ability of the QPSO and the speed of the local search of the BP algorithm, thus improving the stability of the training process of the BP network and the accuracy of the training results. (3) The error analysis of the case study showed that the BP neural network optimized by QPSO had higher prediction accuracy and stability than the multiple regression method, the conventional BP model, the BP model optimized by GA, and the BP model optimized by PSO.The remaining sections of this article are organized as follows. Section2 constructs a detailed construction cost impact index system of subway shield engineering and the cost prediction model based on the BP optimized by the QPSO. Section 3 is the case analysis, which analyses the construction costs of 50 subway shield projects in detail. Section 4 is the discussion. Section 4 details the error and stability analysis of different algorithms and also analyses the influence of the topological structure of the BP neural network. Section 5 is the conclusion. After a brief analysis of the limitations of the research results in this paper, Section 5 summarizes the research results of this article and points out further research directions.


## 2. Materials and Methods

### 2.1. The Indicator System of the Construction Cost of Metro Shield Engineering

#### 2.1.1. Analysis of the Influencing Factors of Metro Shield Construction Cost

We use a project breakdown structure (PBS) to decompose the construction process and to analyse accurately and comprehensively the factors affecting the construction cost of subway shield engineering [28]. Subsequently, we selected the factors affecting the construction cost of metro shield engineering. By referring to the Urban Rail Transit Project Investment Estimation Indicators (GCG 101-2008) promulgated by the Ministry of Housing and Urban-Rural Development of the People’s Republic of China, and the Regulations on Capital Construction Project Construction Cost Management promulgated by the Ministry of Finance of the People’s Republic of China (Caijian (2016) No. 504), we decomposed the construction aspects into four categories: project aspects, construction conditions, technical standards, and project management.(1) Project Aspects. Metro shield engineering refers to the installation of the shield and does not include other aspects of construction, such as connecting the passages and pump houses, reinforcement of risk sources, construction monitoring, and evacuation platforms. The metro shield project includes three parts: the selection of the shield machine, the shield construction at the start or exit stage, and the construction in the shield section. With reference to the engineering practice of the Chengdu Metro Line 11, the Wuhan Metro Line 8, the Tianjin Metro Line 7, and the Wuhan Metro Line 21, we analysed the engineering aspects of the metro shield project. The selection of the shield machine has not only a significant impact on the smooth completion of construction but also on the cost. A shield machine adapted to the engineering and geological characteristics can complete the task quickly with high quality and low risk. Thus, the shield machine, which costs more than 10 million dollars, substantially reduces the amortization cost of the project. At present, there are two types of shield machines: an Earth pressure balance shield machine and a mud-water balance shield machine, which are adapted to different geological characteristics. The permeability coefficient of the stratum is not only an important geological factor but also a vital parameter to determine the selection of the shield machine [29]. According to the national standard of the People’s Republic of China (Code for Construction and Acceptance of Shield Tunnelling Methods (GB 50446-2017)), when the ground seepage coefficient exceeds 10−4 m/s, the mud-water balance shield machine should be considered. When the seepage coefficient is less than 10−7 m/s, the Earth pressure balance shield machine is considered. When the ground seepage coefficient is between the two values, either model can be used, depending on the conditions. The minimum radius of the tunnel line, the maximum slope, and the clearance in the tunnel also affect the selection of the shield machine. The minimum radius of the line should match the size of the shield machine, and an appropriate shield size should be selected according to the clearance in the tunnel. It is worth mentioning that the geological characteristics affect the selection and the efficiency of the shield machine. Therefore, they also have a significant impact on the shield construction costs.In the initial stage of shield construction, the launch shaft is first established and reinforced, followed by the assembly, debugging, and trial excavation. After trial tunnelling has been performed successfully, normal tunnelling begins. During the trial tunnelling process, it is necessary to evaluate the tunnelling parameters of the shield machine to find the optimal parameters. Therefore, the construction cost of the trial tunnelling stage is often significantly higher than that of the formal construction stage. The construction cost of the trial tunnelling stage is proportional to the length of the trial tunnel. The cost of the shield machine’s entry and exit work is primarily determined by the number of times the machine enters and exits the field; the cost is generally calculated based on the number of times.In the construction of the shield section, in addition to hydrogeological conditions affecting the construction cost, segment engineering is the main factor. Different segment thicknesses entail different concrete requirements, and the concrete costs account for a large proportion of the total material cost. It is worth mentioning that since the construction cost analysed in this study is based on the unit price per km for metro shield construction, the length of the shield section has a negligible effect on the unit price.Based on the above analysis, the basic engineering features of the metro shield included the type of shield machine (X1), the permeability coefficient of the stratum (X2), the minimum radius of the tunnel line (X3), the maximum slope (X4), the clearance in the tunnel (X5), the geological characteristics (X6), hydrological conditions (X7), predrilling length (X8), the number of shield machine entrances and exits (X9), and the segment thickness (X10).(2) Construction Conditions. Similar to the construction cost, the construction period is also affected by many factors. Generally, a reduction in the construction period increases the direct costs of construction, whereas the indirect costs are reduced. In contrast, an increase in the construction period may lead to a reduction in direct costs but an increase in indirect costs. Therefore, the construction period is selected as an influencing factor. Due to differences in the area of metro projects, the area directly affects the wages and the price of materials and machinery required by the project. Since different regions have different fee collection standards, the location of the project affects the construction cost.Metro projects are often built in cities or suburbs, where there are already many buildings and municipal pipelines. Certain measures are required during shield construction to protect these buildings and municipal pipelines. Therefore, the more complicated the surrounding environment of a subway project, the more engineering measures are required and the higher the construction cost is. The seismic fortification intensity (seismic capacity of metro projects) is also an important factor affecting construction costs. Different project standards are used for different seismic requirements. The higher the seismic fortification intensity, the higher the engineering standards and quality standards are. Therefore, we regard the seismic fortification intensity as a basic engineering feature.We select the construction period (X11), environmental characteristics (X12), seismic fortification intensity (X13), and other related factors.(3) Technical Standards. According to the Code for the Design of Metro Systems (GB 50157-2013), which is the national standard of the People’s Republic of China, the current design speed of metro systems is divided into several grades, such as 60, 80, 100, 120, and 160. However, most of the lines are designed for a speed of 80 due to the station distance, train performance, and overall cost. If the designed speed is high, but the station distance is small, the train has to slow down to enter the station before. This approach does not represent an advanced design of the system and is a waste of resources. If the train runs fast, the higher the requirements for the train performance and the curve radius of the line, the higher the construction cost is. However, the accessibility is low, and the travel needs of residents in some areas are not considered, which is unreasonable. Therefore, the design speed (X14) is selected as a basic engineering feature that affects the construction cost of the metro shield project.At present, common structural forms of tunnel construction include double-line double tunnels and double-line single tunnels. In theory, double-line single tunnels are more economical because only one excavation is needed, and the cost is reduced [30]. However, the situation is the opposite in practical applications. Due to the need for high technology and a sound geological structure for large cross-section excavation, the construction of a double-line single tunnel is often more expensive. Therefore, a single-line single tunnel is typically used in many more cases. We select the tunnel structure (X15) as a basic engineering feature that affects the construction cost of metro shield engineering.(4) Project Management. Project management includes the management ability of the owner (X16), the management ability of the construction company (X17), and the financing structure (X18). The stronger the project management capabilities of the owners and construction companies are during the life cycle of the metro tunnel project, the better the allocation and coordination of resources are. Scientific decisions are made to facilitate project execution and lower construction costs. Due to extensive investments in metro projects, it requires a long time to recover the costs to benefit the public. When metros are constructed in developing countries, method such as Build-Transfer (BT), Build-Operate-Transfer (BOT), Public-Private Partnership (PPP), Build-Own-Operate (BOO), and other financing models are used [31]. These approaches save government interest expenses, thereby indirectly reducing construction costs. Different financing models have different cost inputs; therefore, the financing models affect project management.


#### 2.1.2. Data Collection

According to the analysis inSection 2.1.1, this section presents the data collection methods of the 19 primary and secondary indicators. X6, X7, X16, and X17 are qualitative indicators, whose scores are obtained from a questionnaire survey of experts. The other indicators are quantitative indicators, and the data are obtained by referring to design documents or field surveys. We focus on the practices used during the construction of Chengdu Metro Line 11, Wuhan Metro Line 8, Tianjin Metro Line 7, and Wuhan Metro Line 21. This section provides a detailed description of the indicators and how to obtain the data, which is convenient for subsequent research on the selection of the index system and case analysis based on rough set theory.Shield Machine Type (X1). The soil pressure balance shield machine has a score of 1, the mud-water balance shield machine has a score of 2, and the other types of shield machines have a score of 3. This indicator has no unit.The Permeability Coefficient of the Formation (X2). The score of this indicator can be obtained by consulting the geological survey or design documents, and the unit is m/s. The score can also be obtained from an on-site survey; the calculation formula [32] is(1)κ=kρgη,whereκ is the permeability coefficient of the formation, k is the permeability of the porous medium, η is the dynamic viscosity coefficient, ρ is the fluid density, and g is the acceleration of gravity. A field test is used to obtain the values of k, ρ, and η so that the permeability coefficient of the formation can be calculated.The Minimum Radius of the Tunnel Line (X3). The score of this indicator can be obtained by consulting the design documents or conducting an on-site investigation; the unit is m.Maximum Slope (X4). The score of this indicator can be obtained by consulting the design documents or conducting an on-site investigation. The unit of this indicator is ‰.Clearance in the Tunnel (X5). The score of this indicator can be obtained by consulting the design documents or conducting an on-site investigation; the unit is m.Geological Features (X6). According to the Technical Code for the Investigation of Geotechnical Engineering (YS 5202-2004), which is the industry standard of the People’s Republic of China, the geological factors that impact subway shield engineering and typical underground engineering construction include the topography, stratum lithology, geological structure, earthquake probability, hydrogeology, natural building materials, and undesirable physical geological phenomena, such as karst, high probability of landslide, collapse, sand liquefaction, and foundation deformation. This qualitative indicator is complex and comprehensive. We divide the impact of the geological features on the construction cost into five levels: (1) very unfavourable, (2) unfavourable, (3) general, (4) favourable, and (5) very favourable.Hydrological Conditions (X7). Hydrogeological conditions refer to the impact of various changes and movements of groundwater within the scope of the subway shield project on the construction cost of the shield project. We divide the impact of hydrological conditions on the construction cost into five levels: (1) very unfavourable, (2) unfavourable, (3) general, (4) favorable, and (5) very favourable.Predrilling Length (X8). The score of this indicator can be obtained by consulting design documents or project management documents, and the unit is m.The Number of Shield Machines Entering and Exiting the Field (X9). The score of this indicator can be obtained by consulting the design documents or conducting an on-site investigation. There is no unit for this indicator.Segment Thickness (X10). The score of this indicator can be obtained by consulting design documents or project management documents, and the unit is m.Construction Period (X11). The score of this indicator can be obtained by consulting the design documents or conducting an on-site investigation; the unit is day (d).Environmental Characteristics (X12). We divide the impact of the subway shield project on the construction cost into five levels: (1) very unfavourable, (2) unfavourable, (3) general, (4) favourable, and (5) very favourable.Seismic Fortification Intensity (X13). According to the National Standard of the People’s Republic of China Standard for the Classification of Seismic Protection of Building Construction (GB 50223-2008), the design basic seismic acceleration value is selected as the score of this index. The unit is gravitational acceleration (g).Design Speed (X14). The score of this indicator can be obtained by consulting the design documents, and the unit is km/h.Tunnel Structure (X15). At present, the common structural forms of tunnel construction are mainly double-line double tunnels and double-line single tunnels. The scores are 1 and 2, respectively.The Project Management Ability of the Owner (X16). Generally, the higher the management ability of the owner, the better the cost control is. The management ability of the contractor is determined according to the quality of the construction unit. This is a qualitative indicator. Since the indicator is dimensionless, the levels are (1) excellent, (2) good, (3) moderate, and (4) poor.The Project Management Ability of the Construction Companies (X17). Generally, the higher the management ability of the contractor, the better the cost control is. The management ability of the contractor is determined according to the quality of the construction unit. This is a qualitative indicator. Since the indicator is dimensionless, the levels are (1) excellent, (2) good, (3) moderate, and (4) poor.The Financing Structure of the Project (X18). There are currently four common financing structures for subway projects. The score of BT is 1, the score of BOT is 2, the score of PPP is 3, the score of BOO is 4, and that of the other modes is 5.


#### 2.1.3. Index Screening Based on the Rough Set Theory

Rough set theory is a mathematical theory proposed by the Polish mathematician Pawlak that can deal with uncertain, incomplete, and incompatible information and knowledge [24]. Attribute reduction based on the rough set is accomplished by deleting unimportant or irrelevant and redundant attributes while maintaining the classification or decision-making ability of the information system and deriving the classification or decision-making of the problem. An information system can be represented by a 4-tuple:(2)S=U,A,V,f,where U is the universe of discourse, A=C∪D=a1,a2,…,an is a limited set of nonempty attributes, C∩D=∅, C are conditional attribute sets, D is a decision attribute set, and A is a set of attribute values. Function A is used to assign a corresponding attribute value to each object in the universe of discourse. The knowledge expression system with conditional attributes and decision attributes is collectively called a decision table, in which the rows represent the records of the object and the columns represent the attributes of the object [25].Definition 1.
The degree of contribution of the conditional attributesC to the decision attributes D is described as(3)γCD=POSCDU,where POSCD is called the positive domain of D with respect to C, which describes the set of all elements in U that can be accurately classified into Class A based on knowledge of U|D. The γCD represents the proportion of objects that can be classified into Decision Class C under conditional attribute U|D. This defines the degree of contribution of the conditional attributes to the decision attributes.Definition 2.
Required attributes and unnecessary attributes.
Suppose there isq∈C,C⊆A in the information system S. If POSCD≠POSC−qD, it is said that q is a necessary attribute of C and must be retained. Otherwise, q is an unnecessary attribute of C and can be deleted. If both q∈C and C are necessary attributes of C, it can independently form a set of features expressing the classification of information systems. In this case, C is called independent; otherwise, it is dependent.Definition 3.
Reduction and core.
In the information systemS, E⊆C⊆A and E are independent and satisfy γED=γCD∧∀E⊂E⇒γE`D≠γCD. E is called a reduction of C, which is denoted as RedC. The set of necessary attributes of C forms the core of C, which is denoted as CoreC=∩RedC. An information system usually contains multiple reductions. If E is independent, then E is the smallest set that maintains the classification ability of the universe U [25, 26].
In this study, the Rosetta software is selected to screen the indicators of the metro shield construction cost. The software is a logical analysis and data processing tool based on rough set theory [33]. It was developed in cooperation with the Department of Computer and Information Science at the Norwegian University of Science and Technology and the Institute of Mathematics at the University of Warsaw, Poland. The software integrates a variety of data preprocessing methods, such as data complementation, discretization, and other functions. It can also perform operations such as rule acquisition and index reduction, as well as data mining and knowledge acquisition [34].
The Chengdu Metro Line 11 (16 shield sections), Wuhan Metro Line 8 (9 shield sections), Tianjin Metro Line 7 (13 shield sections), Wuhan Metro Line 21 (12 shield sections), and others were used in the case study. Quantitative index data were obtained by consulting design materials and field surveys and using other methods. In addition, 100 experts were invited to score each qualitative index. The criteria for selecting experts are work unit, working years, professional titles, and educational background. After two rounds of questionnaires, a total of 72 valid questionnaires were obtained. For simplification, 50 questionnaires were randomly selected. Cronbach’s alpha reliability coefficient of the questionnaires was greater than 0.7, so the questionnaire system was credible [34]. In fact, this paper randomly selected three times, and the values of Cronbach’s alpha of the three questionnaires were 0.7831, 0.831, and 0.847, which were all greater than 0.7. All index score data are shown in Table 1. Due to the length of the paper, Table 1 only lists some data.
The objective was to maintain the classification ability of the knowledge system. Rosetta’s built-in attribute reduction algorithm was used to delete the irrelevant or unimportant secondary indicators and to extract the sensitive predictive indicators of the metro shield construction costs. We extracted the type of shield machine (X1), the geological characteristics (X6), the number of entrances and exits of the shield machine (X9), and the segment thickness (X10) as the attributes of the project aspect. The construction period (X11) and environmental characteristics (X12) were used as the attributes of the construction conditions. The design speed (X14) and tunnel structure (X15) were used as the attributes of the technical standards. The project management ability of the owner (X16), the project management ability of the construction company (X17), and the financing structure of the project (X18) were used as the attributes of project management. In addition, since the design speed in the current metro engineering practice is 80 km/h, we deleted the design speed (X14). We discuss the reason for deleting this indicator in Section 3.3.Table 1
Scores of the indices and the actual construction cost.
No.Actual costX1X2X3X4⋯X15X16X17X18UnitMillions/km—m/sm‰⋯————1281.6217.0×10−82506⋯11232473.7819.2×10−52804⋯13133867.0324.3×10−43004⋯12124387.4911.0×10−52506⋯1222⋮⋮⋮⋮⋮⋮⋮⋮⋮⋮⋮46383.6913.7×10−52504⋯221347738.5014.0×10−52504⋯231348424.3112.8×10−73006⋯221349337.4716.4×10−82506⋯121350549.3418.5×10−62506⋯1113


#### 2.1.4. Development of the Indicator System

The primary and secondary indicators affecting the construction cost of metro tunnel engineering are shown in Table2.Table 2
The indicators affecting the metro tunnel construction cost.
Primary indicatorsSecondary indicatorsType of IndicatorProject aspectsType of shield machine (X1)QuantitativeGeological features (X6)QualitativeNumber of times of shield machine entering and leaving the site (X9)QuantitativeSegment thickness (X10)QuantitativeConstruction conditionsConstruction period (X11)QuantitativeEnvironmental characteristics (X12)QualitativeTechnical standardTunnel structure (X15)QuantitativeProject managementProject management ability of the owner (X16)QualitativeProject management ability of the construction company (X17)QualitativeFinancing structure of the project (X18)Quantitative


### 2.2. The Construction Cost Prediction Model Based on the QPSO-BP Algorithm

The QPSO algorithm has the advantages of speed and a strong global search, unlike the BP model. In addition, QPSO does not require specific characteristic information of the problem, such as derivative gradient information. Therefore, the QPSO and BP models can be combined to take advantage of the generalization ability of the BP model and the ability of the QPSU to optimize the initial weights and thresholds of the layers of the network. This integration improves the neural network’s convergence speed and learning ability.


#### 2.2.1. The BP Model

The BP network was proposed in 1986 by a team of scientists headed by Rumelhart and McClelland. The BP neural network has self-learning and adaptive capabilities. When the output value is different from the expected value, the error is calculated. In this case, the neuron runs the information transmission process opposite to the learning process [35]. Since the principle of the BP neural network learning algorithm is the gradient descent method, the output result, in which the error meets the requirements, can be obtained by repeated training. At present, the three-layer BP network is the most commonly used network. The model structure is shown in Figure 1, where X1, X2, ⋯, Xn are the input values and Y1, Y2, ⋯, Ym are the output results.Figure 1
The structure diagram of the BP neural network.


#### 2.2.2. The QPSO Model

PSO is an optimization method that emerged in the 1990s. It has the advantages of simple structure, easy programming, and small amount of calculation. The search space of the standard PSO is limited, and it is easy to fall into the local extremum [36]. In the standard PSO, the convergence of particles is realized in the form of orbit, and the maximum speed of particles is limited in the process of searching. So, the standard PSO cannot globally converge with probability [37]. In the process of particle convergence in QPSO, particle i keeps approaching Pi until it falls within Pi. There is an attractive potential at Pi during the convergence process, which attracts particles to converge to Pi and makes the whole population have aggregation [18]. Particles in the quantum bound state can appear at any space point with a certain probability [20]. Therefore, particles can search in the whole solution space and thus have better convergence. This is the theoretical reason why QPSO has better computational performance than PSO.Compared with the standard PSO, the QPSO has a different evolutionary search strategy and eliminates the unnecessary velocity vector in the evolutionary equation. In addition, the QPSO has a simpler evolutionary equation and fewer parameters and is easier to control; thus, QPSO has better global search performance.In a population withM particles, each particle moves in an D-dimensional space. The velocity and position vectors of the ith particle are Xi,n=Xi,n1,Xi,n2,…,Xi,nD and Vi,n=Vi,n1,Vi,n2,…,Vi,nD. In the classical PSO, the trajectory of particles [17] is updated using(4)Vi,n+1j=Vi,nj+c1ri,njPi,nj−Xi,nj+c2Ri,njGi,nj−Xi,nj,(5)Xi,n+1j=Xi,nj+Vi,n+1j,where i=1,2,…,M, j=1,2,…,D, c1 and c2 are acceleration factors, ri,nj and Ri,nj are random numbers in the range of 0,1, Vi,nj∈−Vmax,Vmax, Pi,n=Pi,n1,Pi,n2,…,Pi,nD is the individual optimal position of particle i, and Gn=Gn1,Gn2,…,GnD is the global optimal position.In the QPSO algorithm, the position and velocity vectors of each particle are replaced by wave functions [19]. The trajectory of the particles is updated using equations (6)–(8) [19–22]:(6)Xi,n+1j=Pi,nj±α⋅Xi,nj−Ci,nj⋅≤In1ui,n+1j,(7)Cnj=1M∑i=1Mpi,nj,(8)Pi,nj=φi,nj⋅Pi,nj+1−φi,nj⋅Gnj,where Cn=Cn1,Cn2,…,CnD is the average optimal position, which is determined by averaging the individual optimal positions of all particles, φi,nj and ui,n+1j are random numbers uniformly distributed in 0,1, and α is the scaling coefficient that controls the convergence speed. It is the only control parameter in the QPSO algorithm, except for the population size and iteration times.


#### 2.2.3. The Prediction Model Based on the QPSO-BP

The flowchart of the prediction model based on the QPSO-BP is shown in Figure2.Figure 2
The flowchart of the prediction model based on the QPSO-BP.The main steps of using QPSO to optimize the weights and thresholds of the BP neural network are as follows:Step 1. Data collection and normalization.The numerical dimensions differ for the influencing factors, causing a decrease in the convergence speed when using the BP and QPSO algorithms and affecting the accuracy of the model [38]. Equation (9) was used [38] to normalize the data to obtain a range from 0 to 1:(9)xi∗=xi−xminxmax−xmin−1,wherexi is the original data, xmax is the maximum value in the original data sequence, xmin is the minimum value in the original data sequence, and xi∗ is the normalized data.Step 2. Determining the input and output sample sets of the neural network and the topology of the neural network.In the QPSO algorithm, the first-generation population is initialized within the range of0,1, and the size of the initial particle swarm population M, the scaling factor α, and the maximum number of iterations N are set.Step 3. Determining the coding scheme of the neural network weights and thresholds.We use a threshold or a weight to correspond to one dimension in the particle vector. For a specific neural network, if the input hasn dimensions, the number of hidden layers is l, the number of nodes in each corresponding hidden layer is ki, and the output dimension is m; then, the dimension of the particles can be obtained using the following equation [13, 15]:(10)D=nk1+∑i=1l−1kiki+1+kl−1m+∑i=1lki+m.Step 4. Recoding the population individuals of the quantum particle swarm into the weights and thresholds of the neural network.Step 5. Selecting the appropriate fitness function according to the neural network principle and input and output sample set.The objective is to minimize the difference between the expected output of the network and the actual output. The fitness function is expressed as[13](11)fx=∑i=1syi−yi^2,wheres is the number of samples, yi is the actual output of sample i, and yi^ is the predicted output of sample i.Step 6. Updating the particles according to equations (6) to (8).Step 7. Evaluation of the particles and calculation of the fitness value using equation (11).Step 8. If the end conditions are met, the optimization ends and the weights and thresholds at this time are used as the initial weights and thresholds of the neural network. Otherwise, return to Step (4).


## 2.1. The Indicator System of the Construction Cost of Metro Shield Engineering

### 2.1.1. Analysis of the Influencing Factors of Metro Shield Construction Cost

We use a project breakdown structure (PBS) to decompose the construction process and to analyse accurately and comprehensively the factors affecting the construction cost of subway shield engineering [28]. Subsequently, we selected the factors affecting the construction cost of metro shield engineering. By referring to the Urban Rail Transit Project Investment Estimation Indicators (GCG 101-2008) promulgated by the Ministry of Housing and Urban-Rural Development of the People’s Republic of China, and the Regulations on Capital Construction Project Construction Cost Management promulgated by the Ministry of Finance of the People’s Republic of China (Caijian (2016) No. 504), we decomposed the construction aspects into four categories: project aspects, construction conditions, technical standards, and project management.(1) Project Aspects. Metro shield engineering refers to the installation of the shield and does not include other aspects of construction, such as connecting the passages and pump houses, reinforcement of risk sources, construction monitoring, and evacuation platforms. The metro shield project includes three parts: the selection of the shield machine, the shield construction at the start or exit stage, and the construction in the shield section. With reference to the engineering practice of the Chengdu Metro Line 11, the Wuhan Metro Line 8, the Tianjin Metro Line 7, and the Wuhan Metro Line 21, we analysed the engineering aspects of the metro shield project. The selection of the shield machine has not only a significant impact on the smooth completion of construction but also on the cost. A shield machine adapted to the engineering and geological characteristics can complete the task quickly with high quality and low risk. Thus, the shield machine, which costs more than 10 million dollars, substantially reduces the amortization cost of the project. At present, there are two types of shield machines: an Earth pressure balance shield machine and a mud-water balance shield machine, which are adapted to different geological characteristics. The permeability coefficient of the stratum is not only an important geological factor but also a vital parameter to determine the selection of the shield machine [29]. According to the national standard of the People’s Republic of China (Code for Construction and Acceptance of Shield Tunnelling Methods (GB 50446-2017)), when the ground seepage coefficient exceeds 10−4 m/s, the mud-water balance shield machine should be considered. When the seepage coefficient is less than 10−7 m/s, the Earth pressure balance shield machine is considered. When the ground seepage coefficient is between the two values, either model can be used, depending on the conditions. The minimum radius of the tunnel line, the maximum slope, and the clearance in the tunnel also affect the selection of the shield machine. The minimum radius of the line should match the size of the shield machine, and an appropriate shield size should be selected according to the clearance in the tunnel. It is worth mentioning that the geological characteristics affect the selection and the efficiency of the shield machine. Therefore, they also have a significant impact on the shield construction costs.In the initial stage of shield construction, the launch shaft is first established and reinforced, followed by the assembly, debugging, and trial excavation. After trial tunnelling has been performed successfully, normal tunnelling begins. During the trial tunnelling process, it is necessary to evaluate the tunnelling parameters of the shield machine to find the optimal parameters. Therefore, the construction cost of the trial tunnelling stage is often significantly higher than that of the formal construction stage. The construction cost of the trial tunnelling stage is proportional to the length of the trial tunnel. The cost of the shield machine’s entry and exit work is primarily determined by the number of times the machine enters and exits the field; the cost is generally calculated based on the number of times.In the construction of the shield section, in addition to hydrogeological conditions affecting the construction cost, segment engineering is the main factor. Different segment thicknesses entail different concrete requirements, and the concrete costs account for a large proportion of the total material cost. It is worth mentioning that since the construction cost analysed in this study is based on the unit price per km for metro shield construction, the length of the shield section has a negligible effect on the unit price.Based on the above analysis, the basic engineering features of the metro shield included the type of shield machine (X1), the permeability coefficient of the stratum (X2), the minimum radius of the tunnel line (X3), the maximum slope (X4), the clearance in the tunnel (X5), the geological characteristics (X6), hydrological conditions (X7), predrilling length (X8), the number of shield machine entrances and exits (X9), and the segment thickness (X10).(2) Construction Conditions. Similar to the construction cost, the construction period is also affected by many factors. Generally, a reduction in the construction period increases the direct costs of construction, whereas the indirect costs are reduced. In contrast, an increase in the construction period may lead to a reduction in direct costs but an increase in indirect costs. Therefore, the construction period is selected as an influencing factor. Due to differences in the area of metro projects, the area directly affects the wages and the price of materials and machinery required by the project. Since different regions have different fee collection standards, the location of the project affects the construction cost.Metro projects are often built in cities or suburbs, where there are already many buildings and municipal pipelines. Certain measures are required during shield construction to protect these buildings and municipal pipelines. Therefore, the more complicated the surrounding environment of a subway project, the more engineering measures are required and the higher the construction cost is. The seismic fortification intensity (seismic capacity of metro projects) is also an important factor affecting construction costs. Different project standards are used for different seismic requirements. The higher the seismic fortification intensity, the higher the engineering standards and quality standards are. Therefore, we regard the seismic fortification intensity as a basic engineering feature.We select the construction period (X11), environmental characteristics (X12), seismic fortification intensity (X13), and other related factors.(3) Technical Standards. According to the Code for the Design of Metro Systems (GB 50157-2013), which is the national standard of the People’s Republic of China, the current design speed of metro systems is divided into several grades, such as 60, 80, 100, 120, and 160. However, most of the lines are designed for a speed of 80 due to the station distance, train performance, and overall cost. If the designed speed is high, but the station distance is small, the train has to slow down to enter the station before. This approach does not represent an advanced design of the system and is a waste of resources. If the train runs fast, the higher the requirements for the train performance and the curve radius of the line, the higher the construction cost is. However, the accessibility is low, and the travel needs of residents in some areas are not considered, which is unreasonable. Therefore, the design speed (X14) is selected as a basic engineering feature that affects the construction cost of the metro shield project.At present, common structural forms of tunnel construction include double-line double tunnels and double-line single tunnels. In theory, double-line single tunnels are more economical because only one excavation is needed, and the cost is reduced [30]. However, the situation is the opposite in practical applications. Due to the need for high technology and a sound geological structure for large cross-section excavation, the construction of a double-line single tunnel is often more expensive. Therefore, a single-line single tunnel is typically used in many more cases. We select the tunnel structure (X15) as a basic engineering feature that affects the construction cost of metro shield engineering.(4) Project Management. Project management includes the management ability of the owner (X16), the management ability of the construction company (X17), and the financing structure (X18). The stronger the project management capabilities of the owners and construction companies are during the life cycle of the metro tunnel project, the better the allocation and coordination of resources are. Scientific decisions are made to facilitate project execution and lower construction costs. Due to extensive investments in metro projects, it requires a long time to recover the costs to benefit the public. When metros are constructed in developing countries, method such as Build-Transfer (BT), Build-Operate-Transfer (BOT), Public-Private Partnership (PPP), Build-Own-Operate (BOO), and other financing models are used [31]. These approaches save government interest expenses, thereby indirectly reducing construction costs. Different financing models have different cost inputs; therefore, the financing models affect project management.


### 2.1.2. Data Collection

According to the analysis inSection 2.1.1, this section presents the data collection methods of the 19 primary and secondary indicators. X6, X7, X16, and X17 are qualitative indicators, whose scores are obtained from a questionnaire survey of experts. The other indicators are quantitative indicators, and the data are obtained by referring to design documents or field surveys. We focus on the practices used during the construction of Chengdu Metro Line 11, Wuhan Metro Line 8, Tianjin Metro Line 7, and Wuhan Metro Line 21. This section provides a detailed description of the indicators and how to obtain the data, which is convenient for subsequent research on the selection of the index system and case analysis based on rough set theory.Shield Machine Type (X1). The soil pressure balance shield machine has a score of 1, the mud-water balance shield machine has a score of 2, and the other types of shield machines have a score of 3. This indicator has no unit.The Permeability Coefficient of the Formation (X2). The score of this indicator can be obtained by consulting the geological survey or design documents, and the unit is m/s. The score can also be obtained from an on-site survey; the calculation formula [32] is(1)κ=kρgη,whereκ is the permeability coefficient of the formation, k is the permeability of the porous medium, η is the dynamic viscosity coefficient, ρ is the fluid density, and g is the acceleration of gravity. A field test is used to obtain the values of k, ρ, and η so that the permeability coefficient of the formation can be calculated.The Minimum Radius of the Tunnel Line (X3). The score of this indicator can be obtained by consulting the design documents or conducting an on-site investigation; the unit is m.Maximum Slope (X4). The score of this indicator can be obtained by consulting the design documents or conducting an on-site investigation. The unit of this indicator is ‰.Clearance in the Tunnel (X5). The score of this indicator can be obtained by consulting the design documents or conducting an on-site investigation; the unit is m.Geological Features (X6). According to the Technical Code for the Investigation of Geotechnical Engineering (YS 5202-2004), which is the industry standard of the People’s Republic of China, the geological factors that impact subway shield engineering and typical underground engineering construction include the topography, stratum lithology, geological structure, earthquake probability, hydrogeology, natural building materials, and undesirable physical geological phenomena, such as karst, high probability of landslide, collapse, sand liquefaction, and foundation deformation. This qualitative indicator is complex and comprehensive. We divide the impact of the geological features on the construction cost into five levels: (1) very unfavourable, (2) unfavourable, (3) general, (4) favourable, and (5) very favourable.Hydrological Conditions (X7). Hydrogeological conditions refer to the impact of various changes and movements of groundwater within the scope of the subway shield project on the construction cost of the shield project. We divide the impact of hydrological conditions on the construction cost into five levels: (1) very unfavourable, (2) unfavourable, (3) general, (4) favorable, and (5) very favourable.Predrilling Length (X8). The score of this indicator can be obtained by consulting design documents or project management documents, and the unit is m.The Number of Shield Machines Entering and Exiting the Field (X9). The score of this indicator can be obtained by consulting the design documents or conducting an on-site investigation. There is no unit for this indicator.Segment Thickness (X10). The score of this indicator can be obtained by consulting design documents or project management documents, and the unit is m.Construction Period (X11). The score of this indicator can be obtained by consulting the design documents or conducting an on-site investigation; the unit is day (d).Environmental Characteristics (X12). We divide the impact of the subway shield project on the construction cost into five levels: (1) very unfavourable, (2) unfavourable, (3) general, (4) favourable, and (5) very favourable.Seismic Fortification Intensity (X13). According to the National Standard of the People’s Republic of China Standard for the Classification of Seismic Protection of Building Construction (GB 50223-2008), the design basic seismic acceleration value is selected as the score of this index. The unit is gravitational acceleration (g).Design Speed (X14). The score of this indicator can be obtained by consulting the design documents, and the unit is km/h.Tunnel Structure (X15). At present, the common structural forms of tunnel construction are mainly double-line double tunnels and double-line single tunnels. The scores are 1 and 2, respectively.The Project Management Ability of the Owner (X16). Generally, the higher the management ability of the owner, the better the cost control is. The management ability of the contractor is determined according to the quality of the construction unit. This is a qualitative indicator. Since the indicator is dimensionless, the levels are (1) excellent, (2) good, (3) moderate, and (4) poor.The Project Management Ability of the Construction Companies (X17). Generally, the higher the management ability of the contractor, the better the cost control is. The management ability of the contractor is determined according to the quality of the construction unit. This is a qualitative indicator. Since the indicator is dimensionless, the levels are (1) excellent, (2) good, (3) moderate, and (4) poor.The Financing Structure of the Project (X18). There are currently four common financing structures for subway projects. The score of BT is 1, the score of BOT is 2, the score of PPP is 3, the score of BOO is 4, and that of the other modes is 5.


### 2.1.3. Index Screening Based on the Rough Set Theory

Rough set theory is a mathematical theory proposed by the Polish mathematician Pawlak that can deal with uncertain, incomplete, and incompatible information and knowledge [24]. Attribute reduction based on the rough set is accomplished by deleting unimportant or irrelevant and redundant attributes while maintaining the classification or decision-making ability of the information system and deriving the classification or decision-making of the problem. An information system can be represented by a 4-tuple:(2)S=U,A,V,f,where U is the universe of discourse, A=C∪D=a1,a2,…,an is a limited set of nonempty attributes, C∩D=∅, C are conditional attribute sets, D is a decision attribute set, and A is a set of attribute values. Function A is used to assign a corresponding attribute value to each object in the universe of discourse. The knowledge expression system with conditional attributes and decision attributes is collectively called a decision table, in which the rows represent the records of the object and the columns represent the attributes of the object [25].Definition 1.
The degree of contribution of the conditional attributesC to the decision attributes D is described as(3)γCD=POSCDU,where POSCD is called the positive domain of D with respect to C, which describes the set of all elements in U that can be accurately classified into Class A based on knowledge of U|D. The γCD represents the proportion of objects that can be classified into Decision Class C under conditional attribute U|D. This defines the degree of contribution of the conditional attributes to the decision attributes.Definition 2.
Required attributes and unnecessary attributes.
Suppose there isq∈C,C⊆A in the information system S. If POSCD≠POSC−qD, it is said that q is a necessary attribute of C and must be retained. Otherwise, q is an unnecessary attribute of C and can be deleted. If both q∈C and C are necessary attributes of C, it can independently form a set of features expressing the classification of information systems. In this case, C is called independent; otherwise, it is dependent.Definition 3.
Reduction and core.
In the information systemS, E⊆C⊆A and E are independent and satisfy γED=γCD∧∀E⊂E⇒γE`D≠γCD. E is called a reduction of C, which is denoted as RedC. The set of necessary attributes of C forms the core of C, which is denoted as CoreC=∩RedC. An information system usually contains multiple reductions. If E is independent, then E is the smallest set that maintains the classification ability of the universe U [25, 26].
In this study, the Rosetta software is selected to screen the indicators of the metro shield construction cost. The software is a logical analysis and data processing tool based on rough set theory [33]. It was developed in cooperation with the Department of Computer and Information Science at the Norwegian University of Science and Technology and the Institute of Mathematics at the University of Warsaw, Poland. The software integrates a variety of data preprocessing methods, such as data complementation, discretization, and other functions. It can also perform operations such as rule acquisition and index reduction, as well as data mining and knowledge acquisition [34].
The Chengdu Metro Line 11 (16 shield sections), Wuhan Metro Line 8 (9 shield sections), Tianjin Metro Line 7 (13 shield sections), Wuhan Metro Line 21 (12 shield sections), and others were used in the case study. Quantitative index data were obtained by consulting design materials and field surveys and using other methods. In addition, 100 experts were invited to score each qualitative index. The criteria for selecting experts are work unit, working years, professional titles, and educational background. After two rounds of questionnaires, a total of 72 valid questionnaires were obtained. For simplification, 50 questionnaires were randomly selected. Cronbach’s alpha reliability coefficient of the questionnaires was greater than 0.7, so the questionnaire system was credible [34]. In fact, this paper randomly selected three times, and the values of Cronbach’s alpha of the three questionnaires were 0.7831, 0.831, and 0.847, which were all greater than 0.7. All index score data are shown in Table 1. Due to the length of the paper, Table 1 only lists some data.
The objective was to maintain the classification ability of the knowledge system. Rosetta’s built-in attribute reduction algorithm was used to delete the irrelevant or unimportant secondary indicators and to extract the sensitive predictive indicators of the metro shield construction costs. We extracted the type of shield machine (X1), the geological characteristics (X6), the number of entrances and exits of the shield machine (X9), and the segment thickness (X10) as the attributes of the project aspect. The construction period (X11) and environmental characteristics (X12) were used as the attributes of the construction conditions. The design speed (X14) and tunnel structure (X15) were used as the attributes of the technical standards. The project management ability of the owner (X16), the project management ability of the construction company (X17), and the financing structure of the project (X18) were used as the attributes of project management. In addition, since the design speed in the current metro engineering practice is 80 km/h, we deleted the design speed (X14). We discuss the reason for deleting this indicator in Section 3.3.Table 1
Scores of the indices and the actual construction cost.
No.Actual costX1X2X3X4⋯X15X16X17X18UnitMillions/km—m/sm‰⋯————1281.6217.0×10−82506⋯11232473.7819.2×10−52804⋯13133867.0324.3×10−43004⋯12124387.4911.0×10−52506⋯1222⋮⋮⋮⋮⋮⋮⋮⋮⋮⋮⋮46383.6913.7×10−52504⋯221347738.5014.0×10−52504⋯231348424.3112.8×10−73006⋯221349337.4716.4×10−82506⋯121350549.3418.5×10−62506⋯1113


### 2.1.4. Development of the Indicator System

The primary and secondary indicators affecting the construction cost of metro tunnel engineering are shown in Table2.Table 2
The indicators affecting the metro tunnel construction cost.
Primary indicatorsSecondary indicatorsType of IndicatorProject aspectsType of shield machine (X1)QuantitativeGeological features (X6)QualitativeNumber of times of shield machine entering and leaving the site (X9)QuantitativeSegment thickness (X10)QuantitativeConstruction conditionsConstruction period (X11)QuantitativeEnvironmental characteristics (X12)QualitativeTechnical standardTunnel structure (X15)QuantitativeProject managementProject management ability of the owner (X16)QualitativeProject management ability of the construction company (X17)QualitativeFinancing structure of the project (X18)Quantitative


## 2.1.1. Analysis of the Influencing Factors of Metro Shield Construction Cost

We use a project breakdown structure (PBS) to decompose the construction process and to analyse accurately and comprehensively the factors affecting the construction cost of subway shield engineering [28]. Subsequently, we selected the factors affecting the construction cost of metro shield engineering. By referring to the Urban Rail Transit Project Investment Estimation Indicators (GCG 101-2008) promulgated by the Ministry of Housing and Urban-Rural Development of the People’s Republic of China, and the Regulations on Capital Construction Project Construction Cost Management promulgated by the Ministry of Finance of the People’s Republic of China (Caijian (2016) No. 504), we decomposed the construction aspects into four categories: project aspects, construction conditions, technical standards, and project management.(1) Project Aspects. Metro shield engineering refers to the installation of the shield and does not include other aspects of construction, such as connecting the passages and pump houses, reinforcement of risk sources, construction monitoring, and evacuation platforms. The metro shield project includes three parts: the selection of the shield machine, the shield construction at the start or exit stage, and the construction in the shield section. With reference to the engineering practice of the Chengdu Metro Line 11, the Wuhan Metro Line 8, the Tianjin Metro Line 7, and the Wuhan Metro Line 21, we analysed the engineering aspects of the metro shield project. The selection of the shield machine has not only a significant impact on the smooth completion of construction but also on the cost. A shield machine adapted to the engineering and geological characteristics can complete the task quickly with high quality and low risk. Thus, the shield machine, which costs more than 10 million dollars, substantially reduces the amortization cost of the project. At present, there are two types of shield machines: an Earth pressure balance shield machine and a mud-water balance shield machine, which are adapted to different geological characteristics. The permeability coefficient of the stratum is not only an important geological factor but also a vital parameter to determine the selection of the shield machine [29]. According to the national standard of the People’s Republic of China (Code for Construction and Acceptance of Shield Tunnelling Methods (GB 50446-2017)), when the ground seepage coefficient exceeds 10−4 m/s, the mud-water balance shield machine should be considered. When the seepage coefficient is less than 10−7 m/s, the Earth pressure balance shield machine is considered. When the ground seepage coefficient is between the two values, either model can be used, depending on the conditions. The minimum radius of the tunnel line, the maximum slope, and the clearance in the tunnel also affect the selection of the shield machine. The minimum radius of the line should match the size of the shield machine, and an appropriate shield size should be selected according to the clearance in the tunnel. It is worth mentioning that the geological characteristics affect the selection and the efficiency of the shield machine. Therefore, they also have a significant impact on the shield construction costs.In the initial stage of shield construction, the launch shaft is first established and reinforced, followed by the assembly, debugging, and trial excavation. After trial tunnelling has been performed successfully, normal tunnelling begins. During the trial tunnelling process, it is necessary to evaluate the tunnelling parameters of the shield machine to find the optimal parameters. Therefore, the construction cost of the trial tunnelling stage is often significantly higher than that of the formal construction stage. The construction cost of the trial tunnelling stage is proportional to the length of the trial tunnel. The cost of the shield machine’s entry and exit work is primarily determined by the number of times the machine enters and exits the field; the cost is generally calculated based on the number of times.In the construction of the shield section, in addition to hydrogeological conditions affecting the construction cost, segment engineering is the main factor. Different segment thicknesses entail different concrete requirements, and the concrete costs account for a large proportion of the total material cost. It is worth mentioning that since the construction cost analysed in this study is based on the unit price per km for metro shield construction, the length of the shield section has a negligible effect on the unit price.Based on the above analysis, the basic engineering features of the metro shield included the type of shield machine (X1), the permeability coefficient of the stratum (X2), the minimum radius of the tunnel line (X3), the maximum slope (X4), the clearance in the tunnel (X5), the geological characteristics (X6), hydrological conditions (X7), predrilling length (X8), the number of shield machine entrances and exits (X9), and the segment thickness (X10).(2) Construction Conditions. Similar to the construction cost, the construction period is also affected by many factors. Generally, a reduction in the construction period increases the direct costs of construction, whereas the indirect costs are reduced. In contrast, an increase in the construction period may lead to a reduction in direct costs but an increase in indirect costs. Therefore, the construction period is selected as an influencing factor. Due to differences in the area of metro projects, the area directly affects the wages and the price of materials and machinery required by the project. Since different regions have different fee collection standards, the location of the project affects the construction cost.Metro projects are often built in cities or suburbs, where there are already many buildings and municipal pipelines. Certain measures are required during shield construction to protect these buildings and municipal pipelines. Therefore, the more complicated the surrounding environment of a subway project, the more engineering measures are required and the higher the construction cost is. The seismic fortification intensity (seismic capacity of metro projects) is also an important factor affecting construction costs. Different project standards are used for different seismic requirements. The higher the seismic fortification intensity, the higher the engineering standards and quality standards are. Therefore, we regard the seismic fortification intensity as a basic engineering feature.We select the construction period (X11), environmental characteristics (X12), seismic fortification intensity (X13), and other related factors.(3) Technical Standards. According to the Code for the Design of Metro Systems (GB 50157-2013), which is the national standard of the People’s Republic of China, the current design speed of metro systems is divided into several grades, such as 60, 80, 100, 120, and 160. However, most of the lines are designed for a speed of 80 due to the station distance, train performance, and overall cost. If the designed speed is high, but the station distance is small, the train has to slow down to enter the station before. This approach does not represent an advanced design of the system and is a waste of resources. If the train runs fast, the higher the requirements for the train performance and the curve radius of the line, the higher the construction cost is. However, the accessibility is low, and the travel needs of residents in some areas are not considered, which is unreasonable. Therefore, the design speed (X14) is selected as a basic engineering feature that affects the construction cost of the metro shield project.At present, common structural forms of tunnel construction include double-line double tunnels and double-line single tunnels. In theory, double-line single tunnels are more economical because only one excavation is needed, and the cost is reduced [30]. However, the situation is the opposite in practical applications. Due to the need for high technology and a sound geological structure for large cross-section excavation, the construction of a double-line single tunnel is often more expensive. Therefore, a single-line single tunnel is typically used in many more cases. We select the tunnel structure (X15) as a basic engineering feature that affects the construction cost of metro shield engineering.(4) Project Management. Project management includes the management ability of the owner (X16), the management ability of the construction company (X17), and the financing structure (X18). The stronger the project management capabilities of the owners and construction companies are during the life cycle of the metro tunnel project, the better the allocation and coordination of resources are. Scientific decisions are made to facilitate project execution and lower construction costs. Due to extensive investments in metro projects, it requires a long time to recover the costs to benefit the public. When metros are constructed in developing countries, method such as Build-Transfer (BT), Build-Operate-Transfer (BOT), Public-Private Partnership (PPP), Build-Own-Operate (BOO), and other financing models are used [31]. These approaches save government interest expenses, thereby indirectly reducing construction costs. Different financing models have different cost inputs; therefore, the financing models affect project management.


## 2.1.2. Data Collection

According to the analysis inSection 2.1.1, this section presents the data collection methods of the 19 primary and secondary indicators. X6, X7, X16, and X17 are qualitative indicators, whose scores are obtained from a questionnaire survey of experts. The other indicators are quantitative indicators, and the data are obtained by referring to design documents or field surveys. We focus on the practices used during the construction of Chengdu Metro Line 11, Wuhan Metro Line 8, Tianjin Metro Line 7, and Wuhan Metro Line 21. This section provides a detailed description of the indicators and how to obtain the data, which is convenient for subsequent research on the selection of the index system and case analysis based on rough set theory.Shield Machine Type (X1). The soil pressure balance shield machine has a score of 1, the mud-water balance shield machine has a score of 2, and the other types of shield machines have a score of 3. This indicator has no unit.The Permeability Coefficient of the Formation (X2). The score of this indicator can be obtained by consulting the geological survey or design documents, and the unit is m/s. The score can also be obtained from an on-site survey; the calculation formula [32] is(1)κ=kρgη,whereκ is the permeability coefficient of the formation, k is the permeability of the porous medium, η is the dynamic viscosity coefficient, ρ is the fluid density, and g is the acceleration of gravity. A field test is used to obtain the values of k, ρ, and η so that the permeability coefficient of the formation can be calculated.The Minimum Radius of the Tunnel Line (X3). The score of this indicator can be obtained by consulting the design documents or conducting an on-site investigation; the unit is m.Maximum Slope (X4). The score of this indicator can be obtained by consulting the design documents or conducting an on-site investigation. The unit of this indicator is ‰.Clearance in the Tunnel (X5). The score of this indicator can be obtained by consulting the design documents or conducting an on-site investigation; the unit is m.Geological Features (X6). According to the Technical Code for the Investigation of Geotechnical Engineering (YS 5202-2004), which is the industry standard of the People’s Republic of China, the geological factors that impact subway shield engineering and typical underground engineering construction include the topography, stratum lithology, geological structure, earthquake probability, hydrogeology, natural building materials, and undesirable physical geological phenomena, such as karst, high probability of landslide, collapse, sand liquefaction, and foundation deformation. This qualitative indicator is complex and comprehensive. We divide the impact of the geological features on the construction cost into five levels: (1) very unfavourable, (2) unfavourable, (3) general, (4) favourable, and (5) very favourable.Hydrological Conditions (X7). Hydrogeological conditions refer to the impact of various changes and movements of groundwater within the scope of the subway shield project on the construction cost of the shield project. We divide the impact of hydrological conditions on the construction cost into five levels: (1) very unfavourable, (2) unfavourable, (3) general, (4) favorable, and (5) very favourable.Predrilling Length (X8). The score of this indicator can be obtained by consulting design documents or project management documents, and the unit is m.The Number of Shield Machines Entering and Exiting the Field (X9). The score of this indicator can be obtained by consulting the design documents or conducting an on-site investigation. There is no unit for this indicator.Segment Thickness (X10). The score of this indicator can be obtained by consulting design documents or project management documents, and the unit is m.Construction Period (X11). The score of this indicator can be obtained by consulting the design documents or conducting an on-site investigation; the unit is day (d).Environmental Characteristics (X12). We divide the impact of the subway shield project on the construction cost into five levels: (1) very unfavourable, (2) unfavourable, (3) general, (4) favourable, and (5) very favourable.Seismic Fortification Intensity (X13). According to the National Standard of the People’s Republic of China Standard for the Classification of Seismic Protection of Building Construction (GB 50223-2008), the design basic seismic acceleration value is selected as the score of this index. The unit is gravitational acceleration (g).Design Speed (X14). The score of this indicator can be obtained by consulting the design documents, and the unit is km/h.Tunnel Structure (X15). At present, the common structural forms of tunnel construction are mainly double-line double tunnels and double-line single tunnels. The scores are 1 and 2, respectively.The Project Management Ability of the Owner (X16). Generally, the higher the management ability of the owner, the better the cost control is. The management ability of the contractor is determined according to the quality of the construction unit. This is a qualitative indicator. Since the indicator is dimensionless, the levels are (1) excellent, (2) good, (3) moderate, and (4) poor.The Project Management Ability of the Construction Companies (X17). Generally, the higher the management ability of the contractor, the better the cost control is. The management ability of the contractor is determined according to the quality of the construction unit. This is a qualitative indicator. Since the indicator is dimensionless, the levels are (1) excellent, (2) good, (3) moderate, and (4) poor.The Financing Structure of the Project (X18). There are currently four common financing structures for subway projects. The score of BT is 1, the score of BOT is 2, the score of PPP is 3, the score of BOO is 4, and that of the other modes is 5.


## 2.1.3. Index Screening Based on the Rough Set Theory

Rough set theory is a mathematical theory proposed by the Polish mathematician Pawlak that can deal with uncertain, incomplete, and incompatible information and knowledge [24]. Attribute reduction based on the rough set is accomplished by deleting unimportant or irrelevant and redundant attributes while maintaining the classification or decision-making ability of the information system and deriving the classification or decision-making of the problem. An information system can be represented by a 4-tuple:(2)S=U,A,V,f,where U is the universe of discourse, A=C∪D=a1,a2,…,an is a limited set of nonempty attributes, C∩D=∅, C are conditional attribute sets, D is a decision attribute set, and A is a set of attribute values. Function A is used to assign a corresponding attribute value to each object in the universe of discourse. The knowledge expression system with conditional attributes and decision attributes is collectively called a decision table, in which the rows represent the records of the object and the columns represent the attributes of the object [25].Definition 1.
The degree of contribution of the conditional attributesC to the decision attributes D is described as(3)γCD=POSCDU,where POSCD is called the positive domain of D with respect to C, which describes the set of all elements in U that can be accurately classified into Class A based on knowledge of U|D. The γCD represents the proportion of objects that can be classified into Decision Class C under conditional attribute U|D. This defines the degree of contribution of the conditional attributes to the decision attributes.Definition 2.
Required attributes and unnecessary attributes.
Suppose there isq∈C,C⊆A in the information system S. If POSCD≠POSC−qD, it is said that q is a necessary attribute of C and must be retained. Otherwise, q is an unnecessary attribute of C and can be deleted. If both q∈C and C are necessary attributes of C, it can independently form a set of features expressing the classification of information systems. In this case, C is called independent; otherwise, it is dependent.Definition 3.
Reduction and core.
In the information systemS, E⊆C⊆A and E are independent and satisfy γED=γCD∧∀E⊂E⇒γE`D≠γCD. E is called a reduction of C, which is denoted as RedC. The set of necessary attributes of C forms the core of C, which is denoted as CoreC=∩RedC. An information system usually contains multiple reductions. If E is independent, then E is the smallest set that maintains the classification ability of the universe U [25, 26].
In this study, the Rosetta software is selected to screen the indicators of the metro shield construction cost. The software is a logical analysis and data processing tool based on rough set theory [33]. It was developed in cooperation with the Department of Computer and Information Science at the Norwegian University of Science and Technology and the Institute of Mathematics at the University of Warsaw, Poland. The software integrates a variety of data preprocessing methods, such as data complementation, discretization, and other functions. It can also perform operations such as rule acquisition and index reduction, as well as data mining and knowledge acquisition [34].
The Chengdu Metro Line 11 (16 shield sections), Wuhan Metro Line 8 (9 shield sections), Tianjin Metro Line 7 (13 shield sections), Wuhan Metro Line 21 (12 shield sections), and others were used in the case study. Quantitative index data were obtained by consulting design materials and field surveys and using other methods. In addition, 100 experts were invited to score each qualitative index. The criteria for selecting experts are work unit, working years, professional titles, and educational background. After two rounds of questionnaires, a total of 72 valid questionnaires were obtained. For simplification, 50 questionnaires were randomly selected. Cronbach’s alpha reliability coefficient of the questionnaires was greater than 0.7, so the questionnaire system was credible [34]. In fact, this paper randomly selected three times, and the values of Cronbach’s alpha of the three questionnaires were 0.7831, 0.831, and 0.847, which were all greater than 0.7. All index score data are shown in Table 1. Due to the length of the paper, Table 1 only lists some data.
The objective was to maintain the classification ability of the knowledge system. Rosetta’s built-in attribute reduction algorithm was used to delete the irrelevant or unimportant secondary indicators and to extract the sensitive predictive indicators of the metro shield construction costs. We extracted the type of shield machine (X1), the geological characteristics (X6), the number of entrances and exits of the shield machine (X9), and the segment thickness (X10) as the attributes of the project aspect. The construction period (X11) and environmental characteristics (X12) were used as the attributes of the construction conditions. The design speed (X14) and tunnel structure (X15) were used as the attributes of the technical standards. The project management ability of the owner (X16), the project management ability of the construction company (X17), and the financing structure of the project (X18) were used as the attributes of project management. In addition, since the design speed in the current metro engineering practice is 80 km/h, we deleted the design speed (X14). We discuss the reason for deleting this indicator in Section 3.3.Table 1
Scores of the indices and the actual construction cost.
No.Actual costX1X2X3X4⋯X15X16X17X18UnitMillions/km—m/sm‰⋯————1281.6217.0×10−82506⋯11232473.7819.2×10−52804⋯13133867.0324.3×10−43004⋯12124387.4911.0×10−52506⋯1222⋮⋮⋮⋮⋮⋮⋮⋮⋮⋮⋮46383.6913.7×10−52504⋯221347738.5014.0×10−52504⋯231348424.3112.8×10−73006⋯221349337.4716.4×10−82506⋯121350549.3418.5×10−62506⋯1113


## 2.1.4. Development of the Indicator System

The primary and secondary indicators affecting the construction cost of metro tunnel engineering are shown in Table2.Table 2
The indicators affecting the metro tunnel construction cost.
Primary indicatorsSecondary indicatorsType of IndicatorProject aspectsType of shield machine (X1)QuantitativeGeological features (X6)QualitativeNumber of times of shield machine entering and leaving the site (X9)QuantitativeSegment thickness (X10)QuantitativeConstruction conditionsConstruction period (X11)QuantitativeEnvironmental characteristics (X12)QualitativeTechnical standardTunnel structure (X15)QuantitativeProject managementProject management ability of the owner (X16)QualitativeProject management ability of the construction company (X17)QualitativeFinancing structure of the project (X18)Quantitative


## 2.2. The Construction Cost Prediction Model Based on the QPSO-BP Algorithm

The QPSO algorithm has the advantages of speed and a strong global search, unlike the BP model. In addition, QPSO does not require specific characteristic information of the problem, such as derivative gradient information. Therefore, the QPSO and BP models can be combined to take advantage of the generalization ability of the BP model and the ability of the QPSU to optimize the initial weights and thresholds of the layers of the network. This integration improves the neural network’s convergence speed and learning ability.


### 2.2.1. The BP Model

The BP network was proposed in 1986 by a team of scientists headed by Rumelhart and McClelland. The BP neural network has self-learning and adaptive capabilities. When the output value is different from the expected value, the error is calculated. In this case, the neuron runs the information transmission process opposite to the learning process [35]. Since the principle of the BP neural network learning algorithm is the gradient descent method, the output result, in which the error meets the requirements, can be obtained by repeated training. At present, the three-layer BP network is the most commonly used network. The model structure is shown in Figure 1, where X1, X2, ⋯, Xn are the input values and Y1, Y2, ⋯, Ym are the output results.Figure 1
The structure diagram of the BP neural network.


### 2.2.2. The QPSO Model

PSO is an optimization method that emerged in the 1990s. It has the advantages of simple structure, easy programming, and small amount of calculation. The search space of the standard PSO is limited, and it is easy to fall into the local extremum [36]. In the standard PSO, the convergence of particles is realized in the form of orbit, and the maximum speed of particles is limited in the process of searching. So, the standard PSO cannot globally converge with probability [37]. In the process of particle convergence in QPSO, particle i keeps approaching Pi until it falls within Pi. There is an attractive potential at Pi during the convergence process, which attracts particles to converge to Pi and makes the whole population have aggregation [18]. Particles in the quantum bound state can appear at any space point with a certain probability [20]. Therefore, particles can search in the whole solution space and thus have better convergence. This is the theoretical reason why QPSO has better computational performance than PSO.Compared with the standard PSO, the QPSO has a different evolutionary search strategy and eliminates the unnecessary velocity vector in the evolutionary equation. In addition, the QPSO has a simpler evolutionary equation and fewer parameters and is easier to control; thus, QPSO has better global search performance.In a population withM particles, each particle moves in an D-dimensional space. The velocity and position vectors of the ith particle are Xi,n=Xi,n1,Xi,n2,…,Xi,nD and Vi,n=Vi,n1,Vi,n2,…,Vi,nD. In the classical PSO, the trajectory of particles [17] is updated using(4)Vi,n+1j=Vi,nj+c1ri,njPi,nj−Xi,nj+c2Ri,njGi,nj−Xi,nj,(5)Xi,n+1j=Xi,nj+Vi,n+1j,where i=1,2,…,M, j=1,2,…,D, c1 and c2 are acceleration factors, ri,nj and Ri,nj are random numbers in the range of 0,1, Vi,nj∈−Vmax,Vmax, Pi,n=Pi,n1,Pi,n2,…,Pi,nD is the individual optimal position of particle i, and Gn=Gn1,Gn2,…,GnD is the global optimal position.In the QPSO algorithm, the position and velocity vectors of each particle are replaced by wave functions [19]. The trajectory of the particles is updated using equations (6)–(8) [19–22]:(6)Xi,n+1j=Pi,nj±α⋅Xi,nj−Ci,nj⋅≤In1ui,n+1j,(7)Cnj=1M∑i=1Mpi,nj,(8)Pi,nj=φi,nj⋅Pi,nj+1−φi,nj⋅Gnj,where Cn=Cn1,Cn2,…,CnD is the average optimal position, which is determined by averaging the individual optimal positions of all particles, φi,nj and ui,n+1j are random numbers uniformly distributed in 0,1, and α is the scaling coefficient that controls the convergence speed. It is the only control parameter in the QPSO algorithm, except for the population size and iteration times.


### 2.2.3. The Prediction Model Based on the QPSO-BP

The flowchart of the prediction model based on the QPSO-BP is shown in Figure2.Figure 2
The flowchart of the prediction model based on the QPSO-BP.The main steps of using QPSO to optimize the weights and thresholds of the BP neural network are as follows:Step 1. Data collection and normalization.The numerical dimensions differ for the influencing factors, causing a decrease in the convergence speed when using the BP and QPSO algorithms and affecting the accuracy of the model [38]. Equation (9) was used [38] to normalize the data to obtain a range from 0 to 1:(9)xi∗=xi−xminxmax−xmin−1,wherexi is the original data, xmax is the maximum value in the original data sequence, xmin is the minimum value in the original data sequence, and xi∗ is the normalized data.Step 2. Determining the input and output sample sets of the neural network and the topology of the neural network.In the QPSO algorithm, the first-generation population is initialized within the range of0,1, and the size of the initial particle swarm population M, the scaling factor α, and the maximum number of iterations N are set.Step 3. Determining the coding scheme of the neural network weights and thresholds.We use a threshold or a weight to correspond to one dimension in the particle vector. For a specific neural network, if the input hasn dimensions, the number of hidden layers is l, the number of nodes in each corresponding hidden layer is ki, and the output dimension is m; then, the dimension of the particles can be obtained using the following equation [13, 15]:(10)D=nk1+∑i=1l−1kiki+1+kl−1m+∑i=1lki+m.Step 4. Recoding the population individuals of the quantum particle swarm into the weights and thresholds of the neural network.Step 5. Selecting the appropriate fitness function according to the neural network principle and input and output sample set.The objective is to minimize the difference between the expected output of the network and the actual output. The fitness function is expressed as[13](11)fx=∑i=1syi−yi^2,wheres is the number of samples, yi is the actual output of sample i, and yi^ is the predicted output of sample i.Step 6. Updating the particles according to equations (6) to (8).Step 7. Evaluation of the particles and calculation of the fitness value using equation (11).Step 8. If the end conditions are met, the optimization ends and the weights and thresholds at this time are used as the initial weights and thresholds of the neural network. Otherwise, return to Step (4).


## 2.2.1. The BP Model

The BP network was proposed in 1986 by a team of scientists headed by Rumelhart and McClelland. The BP neural network has self-learning and adaptive capabilities. When the output value is different from the expected value, the error is calculated. In this case, the neuron runs the information transmission process opposite to the learning process [35]. Since the principle of the BP neural network learning algorithm is the gradient descent method, the output result, in which the error meets the requirements, can be obtained by repeated training. At present, the three-layer BP network is the most commonly used network. The model structure is shown in Figure 1, where X1, X2, ⋯, Xn are the input values and Y1, Y2, ⋯, Ym are the output results.Figure 1
The structure diagram of the BP neural network.


## 2.2.2. The QPSO Model

PSO is an optimization method that emerged in the 1990s. It has the advantages of simple structure, easy programming, and small amount of calculation. The search space of the standard PSO is limited, and it is easy to fall into the local extremum [36]. In the standard PSO, the convergence of particles is realized in the form of orbit, and the maximum speed of particles is limited in the process of searching. So, the standard PSO cannot globally converge with probability [37]. In the process of particle convergence in QPSO, particle i keeps approaching Pi until it falls within Pi. There is an attractive potential at Pi during the convergence process, which attracts particles to converge to Pi and makes the whole population have aggregation [18]. Particles in the quantum bound state can appear at any space point with a certain probability [20]. Therefore, particles can search in the whole solution space and thus have better convergence. This is the theoretical reason why QPSO has better computational performance than PSO.Compared with the standard PSO, the QPSO has a different evolutionary search strategy and eliminates the unnecessary velocity vector in the evolutionary equation. In addition, the QPSO has a simpler evolutionary equation and fewer parameters and is easier to control; thus, QPSO has better global search performance.In a population withM particles, each particle moves in an D-dimensional space. The velocity and position vectors of the ith particle are Xi,n=Xi,n1,Xi,n2,…,Xi,nD and Vi,n=Vi,n1,Vi,n2,…,Vi,nD. In the classical PSO, the trajectory of particles [17] is updated using(4)Vi,n+1j=Vi,nj+c1ri,njPi,nj−Xi,nj+c2Ri,njGi,nj−Xi,nj,(5)Xi,n+1j=Xi,nj+Vi,n+1j,where i=1,2,…,M, j=1,2,…,D, c1 and c2 are acceleration factors, ri,nj and Ri,nj are random numbers in the range of 0,1, Vi,nj∈−Vmax,Vmax, Pi,n=Pi,n1,Pi,n2,…,Pi,nD is the individual optimal position of particle i, and Gn=Gn1,Gn2,…,GnD is the global optimal position.In the QPSO algorithm, the position and velocity vectors of each particle are replaced by wave functions [19]. The trajectory of the particles is updated using equations (6)–(8) [19–22]:(6)Xi,n+1j=Pi,nj±α⋅Xi,nj−Ci,nj⋅≤In1ui,n+1j,(7)Cnj=1M∑i=1Mpi,nj,(8)Pi,nj=φi,nj⋅Pi,nj+1−φi,nj⋅Gnj,where Cn=Cn1,Cn2,…,CnD is the average optimal position, which is determined by averaging the individual optimal positions of all particles, φi,nj and ui,n+1j are random numbers uniformly distributed in 0,1, and α is the scaling coefficient that controls the convergence speed. It is the only control parameter in the QPSO algorithm, except for the population size and iteration times.


## 2.2.3. The Prediction Model Based on the QPSO-BP

The flowchart of the prediction model based on the QPSO-BP is shown in Figure2.Figure 2
The flowchart of the prediction model based on the QPSO-BP.The main steps of using QPSO to optimize the weights and thresholds of the BP neural network are as follows:Step 1. Data collection and normalization.The numerical dimensions differ for the influencing factors, causing a decrease in the convergence speed when using the BP and QPSO algorithms and affecting the accuracy of the model [38]. Equation (9) was used [38] to normalize the data to obtain a range from 0 to 1:(9)xi∗=xi−xminxmax−xmin−1,wherexi is the original data, xmax is the maximum value in the original data sequence, xmin is the minimum value in the original data sequence, and xi∗ is the normalized data.Step 2. Determining the input and output sample sets of the neural network and the topology of the neural network.In the QPSO algorithm, the first-generation population is initialized within the range of0,1, and the size of the initial particle swarm population M, the scaling factor α, and the maximum number of iterations N are set.Step 3. Determining the coding scheme of the neural network weights and thresholds.We use a threshold or a weight to correspond to one dimension in the particle vector. For a specific neural network, if the input hasn dimensions, the number of hidden layers is l, the number of nodes in each corresponding hidden layer is ki, and the output dimension is m; then, the dimension of the particles can be obtained using the following equation [13, 15]:(10)D=nk1+∑i=1l−1kiki+1+kl−1m+∑i=1lki+m.Step 4. Recoding the population individuals of the quantum particle swarm into the weights and thresholds of the neural network.Step 5. Selecting the appropriate fitness function according to the neural network principle and input and output sample set.The objective is to minimize the difference between the expected output of the network and the actual output. The fitness function is expressed as[13](11)fx=∑i=1syi−yi^2,wheres is the number of samples, yi is the actual output of sample i, and yi^ is the predicted output of sample i.Step 6. Updating the particles according to equations (6) to (8).Step 7. Evaluation of the particles and calculation of the fitness value using equation (11).Step 8. If the end conditions are met, the optimization ends and the weights and thresholds at this time are used as the initial weights and thresholds of the neural network. Otherwise, return to Step (4).


## 3. Case Analysis

### 3.1. Obtaining and Processing Forecast Data

We selected four metro projects of the China State Construction Group Corporation for the case analysis. They include the Chengdu Metro Line 11 with 16 shield tunnels, Wuhan Metro Line 8 with 9 shield tunnels, Tianjin Metro Line 7 with 13 shield tunnels, and Wuhan Metro Line 21 with 12 tunnels. Table1 shows the 50 sets of sample data. The time factors and geographic factors have a significant impact on the cost of the building material, labor, and machinery; therefore, they also have a significant impact on the construction costs, and it is necessary to adjust the 50 sets of construction cost data based on geography and time. It is worth mentioning that, in Table 1, the actual construction costs represent the adjusted data.(1)
We use Chengdu City, Sichuan Province, China as a benchmark and use the regional adjustment coefficient for the adjustment. According to theStandard method of measurement for public utilities works (GB 50857-2013), it can be calculated by the following equation:(12)A=∑i=1nDri∗Pri∑i=1mD0i∗P0i,wheren is the total number of similar local projects, m is the total number of similar projects in other regions, Dri is the unit price of similar local projects, Pri is the weight of these projects. Pri=Gri/∑i=1nGri, Gri is the number of similar projects, and D0i and P0i are the unit price and weight of the same types of project in other regions.(2)
Based on the cost index of 2018, the time adjustment is carried out by using the cost index adjustment. According to theStandard method of measurement for public utilities works (GB 50857-2013), it can be calculated by the following equation:(13)B=∑i=1nDri∗Pri∑i=1mD0i∗P0i,wherem is the total number of similar projects in 2018, n is the total number of similar projects in other years, and Dri is the unit price of similar projects in the reporting period.


### 3.2. Forecast Results

Two sets of data were required for the BP neural network to establish the construction cost prediction model of metro shield engineering, i.e., the training set and test set. The training set data was used to train the neural network, and the test set data was used to validate the neural network. In this study, the training set consisted of 55 randomly selected datasets, and the rest of the data represented the test set. The ratio of training data to test data was 90% to 10%.The parameter settings of the BP neural network [15, 16] were as follows: the number of training times was 1000, the learning rate was 0.1, and the training target was 0.001. The calculation parameters of QPSO [20, 21] included the following: the number of iterations was 1000, the population size was 100, the scaling factor was 0.4, and the maximum error of the iteration termination was 0.00001.The convergence curve of the fitness function of the QPSO is shown in Figure3. The QPSO model converged after about 150 generations.Figure 3
Convergence curves of the QPSO, GA and PSO models.The results of the optimization calculation process of QPSO model are shown in Table3. The error between the 146th iteration and the 147th iteration is larger than the minimum acceptable precision (0.00001), and the error between the 147th iteration and the 148th iteration is smaller than the minimum acceptable precision. The error of the iterative calculation results after 148 generations is also less than 0.00001. The results in Figure 3 and Table 3 indicate that QPSO found the optimal neural network weights and thresholds in the 147th generation.Table 3
Optimization calculation process of the QPSO.
Iteration (n)Fitness (n − 1)Fitness (n)Fitness (n) − Fitness (n − 1)Result1472.3887992.3887990 < 0.00001Continue1482.3887992.3888070.000008 < 0.0001Continue1492.3888072.3888070 < 0.0001Continue10002.3888072.3888070 < 0.0001StopThe prediction results of different calculation models obtained from the Matlab software are shown in Figures4 and 5. In Figures 4 and 5, the predicted error and the relative error value were obtained according to the predicted and actual values. The formula is as follows:(14)E1=γp−γa,(15)E2=γp−γaγa,where E1 is the predicted error, E2 is the relative error, γp is the predicted value, and γa is the actual cost.Figure 4
The predicted errors of different algorithms.Figure 5
The relative errors of different algorithms.The calculation error of the QPSO-BP model is very small, the maximum predicted error is only 2.79, and the maximum relative error is only 1.39%. This proved that QPSO-BP had high calculation accuracy, and this paper detailedly compares the calculation accuracy of different calculation models inSection 4.1.The ratio of training set to test set might affect the prediction result. At present, the ratio of common training set and test set is 90% : 10%, 80% : 20%, or 70% : 30% [39, 40]. Therefore, this paper calculated the influence of these three proportions on the prediction results, and the calculation results are shown in Table 4.Table 4
Prediction results of different training and testing split ratio.
The training and testing split ratio90% : 10%80% : 20%70% : 30%Average absolute error5.2787.31312.763Average relative error (%)1.102.293.98According to the calculation results in Table4, it is not difficult to see that, with the decrease of the proportion of training sets, the calculation accuracy of the prediction results decreases continuously. However, when the proportion of training set is 70%, the calculation accuracy of prediction results is 3.98%. Compared with other methods (the multiple regression, BP, GA-BP, and PSO-BP), the QPSO-BP with 70% training set still has better calculation accuracy.


## 3.1. Obtaining and Processing Forecast Data

We selected four metro projects of the China State Construction Group Corporation for the case analysis. They include the Chengdu Metro Line 11 with 16 shield tunnels, Wuhan Metro Line 8 with 9 shield tunnels, Tianjin Metro Line 7 with 13 shield tunnels, and Wuhan Metro Line 21 with 12 tunnels. Table1 shows the 50 sets of sample data. The time factors and geographic factors have a significant impact on the cost of the building material, labor, and machinery; therefore, they also have a significant impact on the construction costs, and it is necessary to adjust the 50 sets of construction cost data based on geography and time. It is worth mentioning that, in Table 1, the actual construction costs represent the adjusted data.(1)
We use Chengdu City, Sichuan Province, China as a benchmark and use the regional adjustment coefficient for the adjustment. According to theStandard method of measurement for public utilities works (GB 50857-2013), it can be calculated by the following equation:(12)A=∑i=1nDri∗Pri∑i=1mD0i∗P0i,wheren is the total number of similar local projects, m is the total number of similar projects in other regions, Dri is the unit price of similar local projects, Pri is the weight of these projects. Pri=Gri/∑i=1nGri, Gri is the number of similar projects, and D0i and P0i are the unit price and weight of the same types of project in other regions.(2)
Based on the cost index of 2018, the time adjustment is carried out by using the cost index adjustment. According to theStandard method of measurement for public utilities works (GB 50857-2013), it can be calculated by the following equation:(13)B=∑i=1nDri∗Pri∑i=1mD0i∗P0i,wherem is the total number of similar projects in 2018, n is the total number of similar projects in other years, and Dri is the unit price of similar projects in the reporting period.


## 3.2. Forecast Results

Two sets of data were required for the BP neural network to establish the construction cost prediction model of metro shield engineering, i.e., the training set and test set. The training set data was used to train the neural network, and the test set data was used to validate the neural network. In this study, the training set consisted of 55 randomly selected datasets, and the rest of the data represented the test set. The ratio of training data to test data was 90% to 10%.The parameter settings of the BP neural network [15, 16] were as follows: the number of training times was 1000, the learning rate was 0.1, and the training target was 0.001. The calculation parameters of QPSO [20, 21] included the following: the number of iterations was 1000, the population size was 100, the scaling factor was 0.4, and the maximum error of the iteration termination was 0.00001.The convergence curve of the fitness function of the QPSO is shown in Figure3. The QPSO model converged after about 150 generations.Figure 3
Convergence curves of the QPSO, GA and PSO models.The results of the optimization calculation process of QPSO model are shown in Table3. The error between the 146th iteration and the 147th iteration is larger than the minimum acceptable precision (0.00001), and the error between the 147th iteration and the 148th iteration is smaller than the minimum acceptable precision. The error of the iterative calculation results after 148 generations is also less than 0.00001. The results in Figure 3 and Table 3 indicate that QPSO found the optimal neural network weights and thresholds in the 147th generation.Table 3
Optimization calculation process of the QPSO.
Iteration (n)Fitness (n − 1)Fitness (n)Fitness (n) − Fitness (n − 1)Result1472.3887992.3887990 < 0.00001Continue1482.3887992.3888070.000008 < 0.0001Continue1492.3888072.3888070 < 0.0001Continue10002.3888072.3888070 < 0.0001StopThe prediction results of different calculation models obtained from the Matlab software are shown in Figures4 and 5. In Figures 4 and 5, the predicted error and the relative error value were obtained according to the predicted and actual values. The formula is as follows:(14)E1=γp−γa,(15)E2=γp−γaγa,where E1 is the predicted error, E2 is the relative error, γp is the predicted value, and γa is the actual cost.Figure 4
The predicted errors of different algorithms.Figure 5
The relative errors of different algorithms.The calculation error of the QPSO-BP model is very small, the maximum predicted error is only 2.79, and the maximum relative error is only 1.39%. This proved that QPSO-BP had high calculation accuracy, and this paper detailedly compares the calculation accuracy of different calculation models inSection 4.1.The ratio of training set to test set might affect the prediction result. At present, the ratio of common training set and test set is 90% : 10%, 80% : 20%, or 70% : 30% [39, 40]. Therefore, this paper calculated the influence of these three proportions on the prediction results, and the calculation results are shown in Table 4.Table 4
Prediction results of different training and testing split ratio.
The training and testing split ratio90% : 10%80% : 20%70% : 30%Average absolute error5.2787.31312.763Average relative error (%)1.102.293.98According to the calculation results in Table4, it is not difficult to see that, with the decrease of the proportion of training sets, the calculation accuracy of the prediction results decreases continuously. However, when the proportion of training set is 70%, the calculation accuracy of prediction results is 3.98%. Compared with other methods (the multiple regression, BP, GA-BP, and PSO-BP), the QPSO-BP with 70% training set still has better calculation accuracy.


## 4. Discussion

In this paper, a rough set, the QPSO, and BP models were used to establish a prediction model of the construction cost of a subway shield project. This paper has the following limitations. (1) The primary indices of the index system in this study only refers to the engineering practices of the Chengdu Metro Line 11, Wuhan Metro Line 8, Tianjin Metro Line 7, and Wuhan Metro Line 21. If more projects are selected, the index system may change. (2) There is no complete analysis of the influence of input parameters on the prediction results in this paper. (3) Although the QPSO was successfully used to optimize the BP model, and its calculation accuracy was higher than that of common methods (the multiple regression, BP, PSO-BP, and GA-BP), many other optimization methods can also be used to optimize the BP model.


### 4.1. Error Analysis of the Different Computational Models

PSO and GA algorithms were used to optimize the BP model to determine the advantages of QPSO for optimizing the BP model. In the BP algorithm, the selected training function was traingda, the activation function of the hidden layer was logsig, and the activation function of the output layer was purelin [15–18]. We set the target error of training as 10−6 and the maximum number of iterations as 1000. The learning rule of the network was the gradient descent method. In the GA algorithm, the number of individuals in the population was 50, the maximum genetic algebra was 1000, the binary digits were 20, and the generation gap was 0.9 [15]. In the PSO algorithm, the swarm size was 200, the personal learning coefficient and global learning coefficient were both 2, the inertia weights decreased linearly from 0.9 to 0.4, the minimum acceptance accuracy was 0.00001, and the maximum number of iterations was 1000 [16]. In addition to population size and convergence conditions, the key parameter required for QPSO is only the scaling coefficient, whereas the GA algorithm and PSO require more parameters. Compared with the GA and PSO, the QPSO has the advantage of fewer parameters.In addition, in the subsequent error calculation and analysis, we also used the BP algorithm and multivariate linear programming, a commonly used method in engineering practice, to predict the construction cost of subway shield engineering in this case study.The convergence curves of the fitness function of the QPSO, GA, and PSO models are shown in Figure3. The QPSO model converged after about 150 generations, the PSO model converged after about 270 generations, and the GA model converged after about 350 generations.The results in Figure3 and Table 3 indicate that QPSO found the optimal neural network weights and thresholds in the 147th generation. However, the GA and PSO only found the optimal neural network weights and thresholds in the 353rd and 264th generations, respectively. Thus, QPSO has a faster convergence speed than the GA and PSO.According to Figures4 and 5 in Section 3.2, the calculation errors of the multiple regression, BP, GA-BP, and PSO-BP models are substantially larger than those of the QPSO-BP model. Therefore, the QPSO-BP model has the highest accuracy. In addition, the accuracy of the models follows the order QPSO-BP > PSO-BP > GA-BP > BP > multiple regression, which is similar to previous calculation results [22]. This phenomenon could prove the rationality of the case analysis results in this paper to a certain extent.In order to further analyse the calculation accuracy of different calculation methods, we selected four common error analysis tools for machine learning, including the determination coefficient (R2), the mean absolute percentage error (MAPE), the root mean square error (RMSE), and the mean absolute error (MAE) to evaluate the accuracy of several methods. The results of the four types of errors are shown in Figure 6.Figure 6
Error analysis results of the different calculation models.R2 indicates the degree of correlation between the actual value and the predicted value; the calculation method is shown in equation (16). The closer R2 is to 1, the higher the correlation is and vice versa [41]:(16)R2=1−∑i=1Nyiexp−yipre2∑i=1Nyiexp−yiexp¯2,where yiexp is the actual result, yipre is the predicted result, yiexp¯ is the average value of the actual results, and N is the number of test sets.As shown in Figure6, the R2 value of QPSO-BP is 0.9980, indicating that the predicted value obtained from the proposed model is very close to the actual value. Compared with the multiple regression analysis method, BP, GA-BP, and PSO-BP, the R2 value of QPSO-BP is higher than those of multiple regression analysis, BP, GA-BP, and PSO-BP, demonstrating that the proposed model has the best prediction results.TheMAPE represents the average relative error of the prediction results (equation (17)) [22, 41]. When the MAPE is 0%, the prediction result is perfect. The higher the value of MAPE, the worse the accuracy of the prediction result:(17)MAPE=1N∑i=1Nyiexp−yipreyiexp∗100%.TheMAPE of QPSO-BP is the smallest (1.29%), indicating that the prediction result of QPSO-BP is the best among all methods. The MAPE of multiple regression analysis, BP, GA-BP, and PSO-BP are 37.51%, 11.84%, 4.75%, and 2.83%, respectively.TheRMSE measures the deviation between the predicted value and the actual value (equation (18)) [22]. The smaller the RMSE value, the smaller the deviation between the predicted value and the actual value and the higher the prediction accuracy:(18)RMSE=1N∑i=1Nyiexp−yipre2.As shown in Figure6, the RMSE value of QPSO-BP is 6.4363, which is the lowest among all algorithms. The calculation results show that the proposed model provides the best prediction results.TheMAE is the average of the absolute values of the deviations between all predicted values and the arithmetic mean (equation (19)). The MAE avoids the problem of offset and describes the degree of data dispersion. The smaller the MAE value, the better the prediction accuracy:(19)MAE=1N∑i=1Nyiexp−yipre.TheMAE value of QPSO-BP is the smallest among all models, indicating good dispersion of the prediction results and high prediction accuracy of this model. All five test sets achieved good prediction results.We also conducted ten-fold cross-validation to test the accuracy of the algorithm [42] and verify whether the BP model was overfitted. The average accuracy of the ten-fold cross-validation is shown in Table 5.Table 5
Error results of the ten-fold cross-validation.
No.12345678910R20.99800.99810.99790.99750.99830.99800.99760.99760.99810.9982MAPE (%)1.291.291.301.341.271.291.331.331.281.27RMSE6.43636.67896.59736.79106.46726.47806.37896.47896.26716.6781MAE6.0786.38196.38916.47686.3056.26736.37316.27815.73786.2789The errors of ten calculation results are obviously smaller than the traditional BP and other three models, which shows the progressiveness of the calculation method. The accuracy is very high, and the average value of the maximum relative error calculated by the model proposed in this paper is 3.0009%. In addition, the values ofR2, MAPE, RMSE, and MAE of ten calculation results are stable, which proves the stability of the proposed algorithm. Based on the above analysis, we deduced that the BP model of case analysis was not overfitted.The number of input variables has a significant influence on the prediction results of the subway shield construction cost using an artificial intelligence method.InSection 2.1.3 of this paper, the rough set method was used to concentrate the input variables, and eight factors (X2, X3, X4, X5, X7, X8, X13, and X14) which had great correlation with other input variables were deleted. In order to verify the rationality and effectiveness of using rough sets, we investigated the impact of deleting this indicator on the calculation results. The average R2, MAPE, RMSE, and MAE of the ten prediction results changed slightly after these eight influencing factors were added, and the calculation time only increased by 54.7 s. According to the basic principle of the rough set, we thought that these eight factors had great correlation with other factors. It was reasonable to delete these eight factors in Section 2.1.3.


### 4.2. Stability Analysis of Different Computational Models

Stability determines the reliability of the calculation model in prediction. In this paper, the standard deviation was used to judge the stability of the calculation model. Among the 45 training samples in Section3.1, 30, 35, 40, and 45 samples were randomly selected as training sets, and the last 5 samples in Section 3.1 were used as test sets. The calculation was repeated 100 times, and the standard deviations of 5 test sets’ prediction errors are shown in Table 6.Table 6
Stability of calculation results under different sample numbers.
Model30354045Multiple regression33.89027.39020.1122.30BP1.8901.6311.0710.938GA-BP0.8170.6390.4510.381PSO-BP0.6910.5880.4320.340QPSO-BP0.3110.2730.1920.181It can be seen from Table6 that QPSO-BP has a low standard deviation, and the multiple regression had the largest standard deviation. With the increase of training sample size, the computational stability of QPSO-BP had been improved continuously. These results showed that QPSO-BP was more stable than other models.


### 4.3. Influence of the Topological Structure of the BP Model on the Calculation Results

To date, no unified and complete theoretical guidance has been developed for the selection of the BP neural network structure, and experience is commonly used for the selection. If the network structure is extensive, the training efficiency is not high, and there may be overfitting, resulting in low network performance and reduced fault tolerance. If the selection is too small, the network may not converge [43]. Therefore, we investigated the influence of the topological structure of different network models on the prediction results.In the case study, the topological structure of the BP model was determined using equation (10). In line with the research results in related fields [44], we investigated the topological structure of seven different network models based on the QPSO-BP model; the results are shown in Table 7.Table 7
The topological structures of the different network models.
ModelNumber of hidden layersNumber of hidden layer nodesR2MAPE (%)RMSEMAEBP1120.853137.89263.8791130.13150.886726.16153.381078.31180.915211.8454.981953.56212–40.899320.37139.311758.4812–60.920510.4848.528134.7612–80.92419.3846.318932.18GA-BP1120.95329.3731.381928.37150.96486.3927.904223.56180.97614.7522.304720.68212–40.97135.1727.313416.4912–60.98484.3120.844315.0112–80.99103.4817.361116.31PSO-BP1120.97988.3137.342826.44150.98314.4823.18915.31180.99202.8312.886412.76212–40.98475.3030.319717.5012–60.99454.3120.561012.8912–80.99212.309.38196.15QPSO-BP1120.98942.069.1508.367150.99311.477.6796.890180.99801.296.4366.078212–40.99421.377.0376.40112–60.99821.255.9365.56812–80.99851.235.7895.278The average error of the calculation results of the QPSO-BP neural network model was very small, regardless of the topological structure of the network model. The results demonstrate the excellent performance of the QPSO-BP neural network model for forecasting the cost of construction projects. As observed in [45, 46], since the QPSO-BP neural network model has high precision for prediction and cross-validation results, the probability of overfitting or underfitting is reduced. When using the same model, the average error of the prediction results decreases with an increase in the number of hidden layers and nodes, but the calculation time also increases.


## 4.1. Error Analysis of the Different Computational Models

PSO and GA algorithms were used to optimize the BP model to determine the advantages of QPSO for optimizing the BP model. In the BP algorithm, the selected training function was traingda, the activation function of the hidden layer was logsig, and the activation function of the output layer was purelin [15–18]. We set the target error of training as 10−6 and the maximum number of iterations as 1000. The learning rule of the network was the gradient descent method. In the GA algorithm, the number of individuals in the population was 50, the maximum genetic algebra was 1000, the binary digits were 20, and the generation gap was 0.9 [15]. In the PSO algorithm, the swarm size was 200, the personal learning coefficient and global learning coefficient were both 2, the inertia weights decreased linearly from 0.9 to 0.4, the minimum acceptance accuracy was 0.00001, and the maximum number of iterations was 1000 [16]. In addition to population size and convergence conditions, the key parameter required for QPSO is only the scaling coefficient, whereas the GA algorithm and PSO require more parameters. Compared with the GA and PSO, the QPSO has the advantage of fewer parameters.In addition, in the subsequent error calculation and analysis, we also used the BP algorithm and multivariate linear programming, a commonly used method in engineering practice, to predict the construction cost of subway shield engineering in this case study.The convergence curves of the fitness function of the QPSO, GA, and PSO models are shown in Figure3. The QPSO model converged after about 150 generations, the PSO model converged after about 270 generations, and the GA model converged after about 350 generations.The results in Figure3 and Table 3 indicate that QPSO found the optimal neural network weights and thresholds in the 147th generation. However, the GA and PSO only found the optimal neural network weights and thresholds in the 353rd and 264th generations, respectively. Thus, QPSO has a faster convergence speed than the GA and PSO.According to Figures4 and 5 in Section 3.2, the calculation errors of the multiple regression, BP, GA-BP, and PSO-BP models are substantially larger than those of the QPSO-BP model. Therefore, the QPSO-BP model has the highest accuracy. In addition, the accuracy of the models follows the order QPSO-BP > PSO-BP > GA-BP > BP > multiple regression, which is similar to previous calculation results [22]. This phenomenon could prove the rationality of the case analysis results in this paper to a certain extent.In order to further analyse the calculation accuracy of different calculation methods, we selected four common error analysis tools for machine learning, including the determination coefficient (R2), the mean absolute percentage error (MAPE), the root mean square error (RMSE), and the mean absolute error (MAE) to evaluate the accuracy of several methods. The results of the four types of errors are shown in Figure 6.Figure 6
Error analysis results of the different calculation models.R2 indicates the degree of correlation between the actual value and the predicted value; the calculation method is shown in equation (16). The closer R2 is to 1, the higher the correlation is and vice versa [41]:(16)R2=1−∑i=1Nyiexp−yipre2∑i=1Nyiexp−yiexp¯2,where yiexp is the actual result, yipre is the predicted result, yiexp¯ is the average value of the actual results, and N is the number of test sets.As shown in Figure6, the R2 value of QPSO-BP is 0.9980, indicating that the predicted value obtained from the proposed model is very close to the actual value. Compared with the multiple regression analysis method, BP, GA-BP, and PSO-BP, the R2 value of QPSO-BP is higher than those of multiple regression analysis, BP, GA-BP, and PSO-BP, demonstrating that the proposed model has the best prediction results.TheMAPE represents the average relative error of the prediction results (equation (17)) [22, 41]. When the MAPE is 0%, the prediction result is perfect. The higher the value of MAPE, the worse the accuracy of the prediction result:(17)MAPE=1N∑i=1Nyiexp−yipreyiexp∗100%.TheMAPE of QPSO-BP is the smallest (1.29%), indicating that the prediction result of QPSO-BP is the best among all methods. The MAPE of multiple regression analysis, BP, GA-BP, and PSO-BP are 37.51%, 11.84%, 4.75%, and 2.83%, respectively.TheRMSE measures the deviation between the predicted value and the actual value (equation (18)) [22]. The smaller the RMSE value, the smaller the deviation between the predicted value and the actual value and the higher the prediction accuracy:(18)RMSE=1N∑i=1Nyiexp−yipre2.As shown in Figure6, the RMSE value of QPSO-BP is 6.4363, which is the lowest among all algorithms. The calculation results show that the proposed model provides the best prediction results.TheMAE is the average of the absolute values of the deviations between all predicted values and the arithmetic mean (equation (19)). The MAE avoids the problem of offset and describes the degree of data dispersion. The smaller the MAE value, the better the prediction accuracy:(19)MAE=1N∑i=1Nyiexp−yipre.TheMAE value of QPSO-BP is the smallest among all models, indicating good dispersion of the prediction results and high prediction accuracy of this model. All five test sets achieved good prediction results.We also conducted ten-fold cross-validation to test the accuracy of the algorithm [42] and verify whether the BP model was overfitted. The average accuracy of the ten-fold cross-validation is shown in Table 5.Table 5
Error results of the ten-fold cross-validation.
No.12345678910R20.99800.99810.99790.99750.99830.99800.99760.99760.99810.9982MAPE (%)1.291.291.301.341.271.291.331.331.281.27RMSE6.43636.67896.59736.79106.46726.47806.37896.47896.26716.6781MAE6.0786.38196.38916.47686.3056.26736.37316.27815.73786.2789The errors of ten calculation results are obviously smaller than the traditional BP and other three models, which shows the progressiveness of the calculation method. The accuracy is very high, and the average value of the maximum relative error calculated by the model proposed in this paper is 3.0009%. In addition, the values ofR2, MAPE, RMSE, and MAE of ten calculation results are stable, which proves the stability of the proposed algorithm. Based on the above analysis, we deduced that the BP model of case analysis was not overfitted.The number of input variables has a significant influence on the prediction results of the subway shield construction cost using an artificial intelligence method.InSection 2.1.3 of this paper, the rough set method was used to concentrate the input variables, and eight factors (X2, X3, X4, X5, X7, X8, X13, and X14) which had great correlation with other input variables were deleted. In order to verify the rationality and effectiveness of using rough sets, we investigated the impact of deleting this indicator on the calculation results. The average R2, MAPE, RMSE, and MAE of the ten prediction results changed slightly after these eight influencing factors were added, and the calculation time only increased by 54.7 s. According to the basic principle of the rough set, we thought that these eight factors had great correlation with other factors. It was reasonable to delete these eight factors in Section 2.1.3.


## 4.2. Stability Analysis of Different Computational Models

Stability determines the reliability of the calculation model in prediction. In this paper, the standard deviation was used to judge the stability of the calculation model. Among the 45 training samples in Section3.1, 30, 35, 40, and 45 samples were randomly selected as training sets, and the last 5 samples in Section 3.1 were used as test sets. The calculation was repeated 100 times, and the standard deviations of 5 test sets’ prediction errors are shown in Table 6.Table 6
Stability of calculation results under different sample numbers.
Model30354045Multiple regression33.89027.39020.1122.30BP1.8901.6311.0710.938GA-BP0.8170.6390.4510.381PSO-BP0.6910.5880.4320.340QPSO-BP0.3110.2730.1920.181It can be seen from Table6 that QPSO-BP has a low standard deviation, and the multiple regression had the largest standard deviation. With the increase of training sample size, the computational stability of QPSO-BP had been improved continuously. These results showed that QPSO-BP was more stable than other models.


## 4.3. Influence of the Topological Structure of the BP Model on the Calculation Results

To date, no unified and complete theoretical guidance has been developed for the selection of the BP neural network structure, and experience is commonly used for the selection. If the network structure is extensive, the training efficiency is not high, and there may be overfitting, resulting in low network performance and reduced fault tolerance. If the selection is too small, the network may not converge [43]. Therefore, we investigated the influence of the topological structure of different network models on the prediction results.In the case study, the topological structure of the BP model was determined using equation (10). In line with the research results in related fields [44], we investigated the topological structure of seven different network models based on the QPSO-BP model; the results are shown in Table 7.Table 7
The topological structures of the different network models.
ModelNumber of hidden layersNumber of hidden layer nodesR2MAPE (%)RMSEMAEBP1120.853137.89263.8791130.13150.886726.16153.381078.31180.915211.8454.981953.56212–40.899320.37139.311758.4812–60.920510.4848.528134.7612–80.92419.3846.318932.18GA-BP1120.95329.3731.381928.37150.96486.3927.904223.56180.97614.7522.304720.68212–40.97135.1727.313416.4912–60.98484.3120.844315.0112–80.99103.4817.361116.31PSO-BP1120.97988.3137.342826.44150.98314.4823.18915.31180.99202.8312.886412.76212–40.98475.3030.319717.5012–60.99454.3120.561012.8912–80.99212.309.38196.15QPSO-BP1120.98942.069.1508.367150.99311.477.6796.890180.99801.296.4366.078212–40.99421.377.0376.40112–60.99821.255.9365.56812–80.99851.235.7895.278The average error of the calculation results of the QPSO-BP neural network model was very small, regardless of the topological structure of the network model. The results demonstrate the excellent performance of the QPSO-BP neural network model for forecasting the cost of construction projects. As observed in [45, 46], since the QPSO-BP neural network model has high precision for prediction and cross-validation results, the probability of overfitting or underfitting is reduced. When using the same model, the average error of the prediction results decreases with an increase in the number of hidden layers and nodes, but the calculation time also increases.


## 5. Conclusions

Prior to construction, it is of great significance to perform an accurate and science-based prediction of the construction cost for the smooth implementation of a subway shield project. Due to the complexity of subway shield construction, many factors affect the project construction cost, and the cost prediction is complex and involves many nonlinear factors. Therefore, we decomposed the project using a project analysis method, obtained the factors influencing the cost of the subway shield project according to the Chinese national standards and expert opinion, and used rough set theory to reduce the number of influencing factors, thus avoiding redundancy between the engineering features. The case analysis of 50 subway shield projects in China indicated that the key factors affecting the project construction cost which included the type of shield machine, geological features, number of times that the shield machine entered and left the site, segment thickness, construction period, environmental characteristics, tunnel structure, project management ability of the owner, project management ability of the construction company, and financing structure of the project. Thus, several factors affected the project construction cost of the subway shield. The BP neural network optimized by QPSO was used to predict the construction cost of the subway shield. This method integrated the global convergence of QPSO and the ability to perform nonlinear analysis of the BP algorithm, thus improving the accuracy and stability of BP network training. The proposed model was well suited for the complex nonlinear problem of cost prediction in subway shield construction. The results of the absolute error, relative error,R2, MAPE, RMSE, and MAE showed that the QPSO-BP model provided higher prediction accuracy than the multiple regression, BP, GA-BP, and PSO-BP methods.In the future, scholars should develop an information database of the construction cost of subway shield projects to investigate the influencing factors more comprehensively. In addition, the use of other artificial intelligence algorithms to optimize the BP neural network is also worthy of further study.



---

*Source: 6692130-2020-12-19.xml*