Research of an EPB shield pressure and depth prediction model based on deep neural network and its control device

<!-- metadata commented in wiki content


<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
'''Research of an EPB shield pressure and depth prediction model based on deep neural network and its control device'''</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Jiacheng Shao<sup>12</sup>,Jingxiu Ling<sup>123a</sup>,Rongchang zhang<sup>12</sup>,Xiaoyuan Cheng<sup>12</sup>,Hao Zhang<sup>3</sup></div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<sup>1</sup>The Key Laboratory of Intelligent Machining Technology and Equipment, Fujian University of Technology, China</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<sup>2</sup> School of Mechanical and Automotive Engineering, Fujian University of Technology, China</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<sup>3</sup>CSCEC Strait Construction and Development Co, Ltd, Fuzhou, China</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
<sup>a</sup>Corresponding Autor，E-mail：[mailto:ljxyxj@fjut.edu.cn ljxyxj@fjut.edu.cn]</div>
-->
==Abstract==

Based on the construction data of Fuzhou Metro Line 4 in Fujian Province, China, this paper proposes a soil pressure prediction model that combines Long Short-Term Memory (LSTM) and Particle Swarm Optimization (PSO). The values of Mean Absolute Error, Mean Squared Error, and Coefficient of Determination are 0.007MPa, 0.007%, and 0.93, respectively, indicating an improvement in accuracy.Wang-Mendel algorithm is used to establish fuzzy rules. The Mean Absolute Error and Mean Squared Error of the rotating speed of the screw machine are 0.065rpm and 1.528%, and the Coefficient of Determination is 0.82. The calculation accuracy of this algorithm is high.A set of knob intelligent control device is developed.The Mean Absolute Error and Mean Squared Error of 0.015rpm and 0.392%, respectively, and the Coefficient of Determination of 0.95, indicating a small execution error of the device. This paper provides a new and effective method for the control of EPB shield pressure.

'''Keywords''': EPB, soil pressure prediction, LSTM, PSO

==1. Introduction==

EPB uses dynamic balance control of water and soil pressure in front of the excavation face to avoid surface deformation. Soil pressure control involves multi-field coupling of multiple subsystems, such as the cutterhead system, propulsion system, and screw conveyor system, as well as the comprehensive effects of multiple factors such as excavation parameters and muck properties. Therefore, establishing a soil pressure control model is one of the key challenges in shield tunneling technology.Currently, the control of soil pressure in shield tunneling is generally achieved through two approaches: the establishment of a mechanism model based on mathematical and physical relationships or the use of intelligent methods such as neural networks based on on-site data. With an accurate model, it is possible to make correct predictions of future states based on past and current conditions and then optimize and control controllable variables to achieve effective control of soil pressure.

Previous studies have conducted extensive research on earth pressure balance.Hu et al. used electro-hydraulic proportional control technology to control the propulsion speed or screw machine speed in real-time on a simulation experimental platform, demonstrating that controlling the screw machine speed can achieve control of soil pressure. Based on this, they further used self-tuning PID to regulate the screw conveyor speed to achieve control of soil pressure [1]. Wang and Fu derived three basic equations for EPB shield based on model test results, and then obtained two total balance equations for EPB shield, establishing a mathematical model for EPB shield excavation [2]. Wang et al. proposed a soil pressure balance control system based on the feedback of advance speed and soil pressure, building upon previous mechanism models [3]. In recent years, neural networks have been widely used in soil pressure prediction and shield tunneling. Yeh pioneered the use of neural networks in the control of soil pressure by using a backpropagation neural network to establish a control model with the current soil pressure, advancing speed, screw conveyor speed, next moment advancing speed, next moment screw conveyor speed as inputs and next moment soil pressure as output [4]. Wei et al. combined Monte Carlo method and BP neural network model to predict the TBM's advancing speed [5]. Gao et al. used three types of recurrent neural networks (RNNs) to process real-time prediction data based on on-site TBM operating parameters, and showed that the three proposed RNN-based predictors outperformed classical regression models [6]. Liu et al. developed a long short-term memory (LSTM) network based on a global attention mechanism to predict the lithology of the excavation face [7]. Cachim and Bezuijen used a multi-layer feedforward artificial neural network to predict the cutterhead torque of an EPB-TBM, and used a time-series neural network to predict the torque as a function of the measured torque and injected foam volume at the previous time step [8]. Song et al. proposed a time series segmentation method based on fuzzy c-means clustering for TBM time series data [9]. Shi et al. proposed a new clustering algorithm based on support vector regression and fuzzy c-means algorithm (SVR-FCM) [10]. Sun et al. developed a prediction model based on integrated heterogeneous in-situ data using data-driven techniques. This method was applied in the field of non-homogeneous tunnel boring machine, which not only accurately predicted the dynamic characteristics of the load, but also accurately estimated the statistical characteristics of the load [11].

Some scholars have made research on control and particle swarm optimization. Wen et al. proposed a soil discharge prediction method based on an adaptive adjustment strategy and parameters to achieve optimal regression. The study analyzed the factors affecting soil flow data models using actual engineering data and found that the screw conveyor speed and total thrust were the main factors affecting the soil flow rate [12]. Liu et al. proposed an optimized control method for earth pressure based on ADHDP, which realized the online autonomous learning and adaptive control of shield tunneling chamber pressure [13]. Zhang used artificial experience to train and learn the membership functions and fuzzy rule parameters of a fuzzy system using neural networks. They established an ANFIS model for predicting the soil pressure in the shield tunneling machine and used genetic algorithms for optimizing control [14]. Zhu  used a fuzzy PID control algorithm to control the factors that affect the earth pressure, and verified the correctness of the algorithm through experimental simulation [15]. Li and Gong proposed a model predictive control (MPC) system for predicting and controlling mud pressure balance during construction. The system consists of a diagonal recurrent neural network (DRNN), an optimizer that generates optimal air pressure and slurry level based on multi-step prediction, and an evolved particle swarm optimization (EPSO) algorithm. It can effectively adjust the mud circulation and air pressure maintenance system according to geological conditions [16]. Li and Li used a BP neural network for prediction and adopted the PSO algorithm to calculate the balanced control of the screw conveyor speed required for the next moment[17]. Zhang et al. proposed a method for controlling the shield tunneling machine's earth pressure using a fuzzy PID control approach based on a mechanistic model. The method was compared with the traditional PID control method, and the results showed that the proposed approach was more effective [18]. Yang  developed a PSO-BP model for predicting the soil pressure in the tunnel and obtained the relationship between the predicted soil pressure and the screw conveyor speed using polynomial fitting to achieve control of soil pressure [19].

After reviewing the research status above, it was found that the input parameters considered in the mechanism control model are limited, and it is difficult to obtain some parameters such as the screw conveyor excavation efficiency and the tangential modulus of deformation of the spoil. In intelligent control models such as neural networks, after the pressure prediction model is obtained, the corresponding screw conveyor speed needs to be calculated using algorithms such as PSO and genetic algorithms, which is time-consuming. In order to facilitate engineering application, this paper will establish a PSO-LSTM deep neural network prediction model, use the WM algorithm to establish fuzzy rules for easy calculation, and develop an intelligent control device. The feasibility of the device will be verified through actual engineering data.

==2. Theoretical background==

===2.1 LSTM neural network===

In many studies that use neural networks, it is necessary to deal with sequential data and utilize its historical information. Therefore, it is necessary to enable neural networks to have memory function. The long short-term memory (LSTM) network is selected, and its structure is shown in [[#img-1|Figure 1]].

<div id='img-1'></div>
{| class="wikitable" style="margin: 0em auto 0.1em auto;border-collapse: collapse;width:auto;" 
|-style="background:white;"
|style="text-align: center;padding:10px;"| [[Image:Draft_Shao_683360558-image1.png|342px]]
|-
| style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 1'''. LSTM network structure
|}


LSTM network introduces an internal state <math>c_t</math> specifically for the linear recurrent information transfer, as shown in Eq.(1). At the same time, it transfers the non-linear output information to the external state <math>h_t</math> of the hidden layer, as shown in Eq.(2):

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>c_t=f_t\odot c_{t-1}+i_t\odot g_t</math> 
|}
| style="width: 5px;text-align: right;white-space: nowrap;" |(1)
|}

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>h_t=o_t\odot</math> 
|}
| style="width: 5px;text-align: right;white-space: nowrap;" |(2)
|}

In Eqs.(1) and (2), <math>c_t</math> represents, at each time step <math> t </math>, the internal state <math>c_t</math> stores the historical information from the beginning up to the current time step; ⊙denotes element-wise multiplication of vectors; <math>c_{t-1}</math> is the memory cell from the previous time step; <math>h_t</math> is both the current time step <math>t</math>'s output of the hidden layer and the input of the hidden layer for the next time step <math>t+1</math>; <math>\tanh(x)</math> is the activation function, defined by Eq.(3), with a range of <math>(-1,1)</math>；

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>\tanh=\frac{\exp\left(x\right)-\exp\left(-x\right)}{\exp\left(x\right)+ \exp\left(-x\right)}</math> 
|}
| style="width: 5px;text-align: right;white-space: nowrap;" |(3)
|}

where <math display="inline">g_t  </math> is the candidate state obtained through a nonlinear function, as shown in Eq.(4)

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>g_t=\tanh\left(W_{ig}x_t+b_{ig}+W_{hg}h_{t-1}+b_{hg}\right)</math> 
|}
| style="width: 5px;text-align: right;white-space: nowrap;" |(4)
|}

where <math display="inline">f_t \in [0,1]</math> is the forget gate, which controls how much information from the previous internal state <math>c_{t-1}</math> should be forgotten, as shown in Eq.(5)

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>i_t=\sigma \left(W_{ii}x_t+b_{ii}+W_{hi}h_{t-1}+b_{hi}\right)</math> 
|}
| style="width: 5px;text-align: right;white-space: nowrap;" |(5)
|}

where <math display="inline">i_t \in [0,1]</math>  represents the input gate, which controls how much information the candidate state <math display="inline">g_t  </math> should store at the current time step, as shown in Eq.(6)

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>f_t=\sigma (W_{if}x_t+b_{if}+W_{hf}h_{t-1}+b_{hf})</math> 
|}
| style="width: 5px;text-align: right;white-space: nowrap;" |(6)
|}

where <math display="inline">o_t \in [0,1]</math> represents the output gate, which controls how much information from the current internal state <math display="inline">c_t  </math> is output to the external state <math display="inline">h_t  </math>, as shown in Eq.(7)

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>o_t=\sigma \left(W_{io}x_t+b_{io}+W_{ho}h_{t-1}+b_{ho}\right)</math> 
|}
| style="width: 5px;text-align: right;white-space: nowrap;" |(7)
|}

where <math display="inline">x_t</math>  is the input at the current time step, and <math>h_{t-1}</math>  is the external state at the previous time step; <math display="inline">w_{ii}</math>, <math display="inline">b_{ii}</math>, <math display="inline">w_{if}</math>, <math display="inline">b_{if}</math>, <math display="inline">w_{io}</math>, <math display="inline">b_{io}</math>, <math display="inline">w_{ig}</math>, <math display="inline">b_{ig}</math> are parameter matrices for the linear transformation of input <math display="inline">x_t</math>; <math display="inline">w_{hi}</math>, <math display="inline">b_{hi}</math>, <math display="inline">w_{hf}</math>, <math display="inline">b_{hf}</math>, <math display="inline">w_{ho}</math>, <math display="inline">b_{ho}</math>, <math display="inline">w_{hg}</math>, <math display="inline">b_{hg}</math> are the parameter matrices for the linear transformation of <math>h_{t-1}</math>; <math display="inline">\sigma (x)</math> is the activation function, whose range is <math>(0,1)</math>, and is defined by Eq.(8)

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>\sigma \left(x\right)=\frac{1}{1+exp\left(-x\right)}</math> 
|}
| style="width: 5px;text-align: right;white-space: nowrap;" |(8)
|}


The computation process of LSTM neural network is as follows:

(1) Calculate three gates and a candidate state <math display="inline">g_t  </math> using the current input <math display="inline">x_t</math>  and the previous external state <math>h_{t-1}</math>.

(2) Update the memory cell <math display="inline">c_t  </math> through the forget gate <math display="inline">f_t  </math> and the input gate <math display="inline">o_t  </math>.

(3) Use the output gate <math display="inline">o_t  </math> to pass the information from the internal state to the external state <math display="inline">h_t  </math>.

===2.2 PSO (Particle Swarm Optimization)===

Particle Swarm Optimization (PSO) is mainly inspired by the collective behavior of bird flocking and the cognitive mechanism of humans.Assuming the position of particle i randomly placed in the search space is ''x''<sub>i</sub>''=''（''x''<sub>1</sub>'',x''<sub>2</sub>'',···,x''<sub>id</sub>'',···,x''<sub>in</sub>）and its velocity is ''v<sub>i</sub>=''(''v<sub>i</sub>''<sub>1</sub>'',v<sub>i</sub>''<sub>2</sub>'',···,v<sub>i</sub>''<sub>d</sub>'',···,v<sub>i</sub>''<sub>n</sub>), then the individual best value of the particle based on the fitness function is ''P<sub>i</sub>=''(''P<sub>i</sub>''<sub>1</sub>'',P<sub>i</sub>''<sub>2</sub>'',···,P<sub>i</sub>''<sub>d</sub>'',···,P<sub>i</sub>''<sub>n</sub>), and the global best value of the swarm is ''P<sub>gi</sub>=''(''P<sub>g</sub>''<sub>1</sub>'',P<sub>g</sub>''<sub>2</sub>'',···,P<sub>g</sub>''<sub>d</sub>'',···,P<sub>g</sub>''<sub>n</sub>). The evolution equations for the particle swarm algorithm are given by Eqs. (9) and (10)

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>  V_{id}^{k+1}=\omega V_{id}^k+c_1r_1(P_{id}^k-X_{id}^k)+c_2r_2(P_{gd}^k-X_{gd}^k) </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" |(9)
|}

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math> X_{id}^{k+1}=X_{id}^k+V_{id}^{k+1}  </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" |(10)
|}

In the equation, <math>\omega </math> is the weight coefficient: <math>d = 1 , 2, 3... N </math>; <math>k</math> is the current iteration number; <math>c_1</math>, <math>c_2</math>, <math>r_1</math>, <math>r_2</math> are non-negative constants randomly distributed in <math>[0,1]</math>; to avoid large-scale blind search, the range of particle positions and velocities are limited to <math>[-X_{\max}, X_{\max}]</math>  and <math>[-V_{\max}, V_{\max}]</math> , respectively.

In this paper, PSO algorithm is mainly used to optimize the hyperparameters in the soil pressure prediction model. Therefore, the error between the predicted values obtained by the prediction model and the actual values is used as a reference value for the fitness function, and the optimization goal is to find the best combination of hyperparameters for the deep neural network. The relationship between the model to be optimized and the optimization goal can be expressed as Eq.(11):

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math> loss_{\min}=\frac{\displaystyle\sum {\left(\mbox{y}_{\mbox{test}}-\mbox{y}_{\mbox{true}}\right)}^2}{\mbox{N}_{\mbox{test}}}  </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" |(11)
|}

where <math>{y}_{\mbox{test}}</math> is the predicted result of the model; <math>{y}_{\mbox{true}}</math> is the actual value; <math>N_{\mbox{test}}</math> is the number of test samples selected.

The process of optimizing the prediction model using PSO can be summarized as follows:

(1) Initialize PSO parameters and set initial values for hyperparameters of the deep neural network prediction model. Calculate the initial position for each particle, denoted as PBest, and set the best value in PBest as the global optimal value gBest.

(2) Calculate the fitness value of the prediction model by using the error between the predicted values obtained by each particle and the actual values, and record the current particle position and corresponding fitness function value.

(3) Compare the fitness value of the particle with the best fitness value it has reached, update the particle's velocity and position, and optimize.

(4) If the number of iterations is not completed, return to step 2 to continue the operation, and finally obtain the globally optimal hyperparameter combination.

===2.3 Wang-Mendel algorithm===

The Wang-Mendel (WM) method proposed by Wang and Mendel can effectively generate fuzzy rules from sample data [20]. This method partitions the input space into fuzzy regions and then uses a lookup table to extract rules for each fuzzy subspace. The WM method does not require repeated learning and can use a large amount of input-output data to obtain a complete fuzzy rule base that covers the entire input domain [21]. The WM method has become a classic method for fuzzy rule generation and has been widely used in various applications [22].

Therefore, in this study, the WM (Wang-Mendel) algorithm was used to establish a fuzzy rule base, which mainly includes three steps: rule extraction, rule merging, and rule inference [21].

Rule extraction refers to extracting fuzzy rules from input factor data. Firstly, the input factors <math>x_j</math> should be divided into <math>m_j</math> fuzzy subsets. In this paper, the input data is in the form of n inputs and one output, which can be expressed as (<math>x^{(I)},y)</math>,  <math>I=1,2,\ldots ,n</math>. Where <math>x^{(I)}\in R^n</math> represents the input factors and <math>y\in R</math> represents the output factor. The input factor data is fuzzified and divided into L fuzzy subspaces using the triangular membership function Eq.(14), as shown in Eq.(13). Then, one rule can be extracted from each of the <math>  k</math> data points, as shown in Eq.(12). In this way, <math>  k</math> rules are obtained

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>  if{\mbox{ x}}_1\mbox{ is }A_1^{l_1}\mbox{ and}\cdots {\mbox{and x}}_n\mbox{ is }A_n^{l_n},\mbox{then y is }B^l </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" |(12)
|}

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>  L=\prod_{j=1}^nm_j </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" |(13)
|}

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math> u(y,y_c,\sigma )=\frac{\sigma -\vert y-y_c\vert }{\sigma } </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" |(14)
|}

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math> y_c{}^l=\frac{\displaystyle\sum_{k=1}^{N_k}y_k^{\left(i\right)}{\omega }_k^{\left(i\right)}}{\displaystyle\sum_{k=1}^{N_k}{\omega }_k^{\left(i\right)}} </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" |(15)
|}

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math> \sigma =\frac{\displaystyle\sum_{k=1}^{Nk}\vert y_k^{\left(i\right)}-y_c^l\vert {\omega }_k^{\left(i\right)}}{\displaystyle\sum_{k=1}^{Nk}{\omega }_k^{\left(i\right)}} </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" |(16)
|}

where <math>l=1,\ldots,L</math>; there are <math>m_n</math> fuzzy sets <math>\{A^{lb}_n|l_b=1,\ldots ,m_n\}</math> on the input domain of input factor <math>x_n</math>. <math>A^{lb}_n</math> represents the fuzzy set on input <math>x_b</math> with index <math>l_d</math>; <math>B^l</math> represents the fuzzy set on the output domain <math>U_B^l(y)=u(y,y^l_c,\, \sigma^l)</math>, where <math>\sigma</math> is the width of the fuzzy set, <math>y_c</math> is the center point of the fuzzy set, defined by Eqs.(15) and (16), where <math>y\in [y_c -\sigma,y_c+\sigma]</math>. <math>l</math> is the rule index corresponding to it in the rule base. <math>L</math> fuzzy sets correspond to <math>L</math> fuzzy rules. Then, the rules with the same premise conditions are merged, and the confidence degrees of the duplicate and conflicting rules are calculated. By comparing the confidence degrees, the rules with low confidence are removed. Suppose there are <math display="inline"> H </math> groups of duplicate rules, each containing <math display="inline"> N_h </math> rules. The weight of each rule is calculated using Eq.17. The larger the weight, the closer the rule is to the center of the fuzzy space. Then, <math display="inline"> N_h </math> rules are merged into a single rule in the form of Eq.(12), and the reliability of the rule is calculated using Eq.(18). Finally, the reliable rules are selected and added to the rule base.

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;" 
|-
| style="text-align: center;" |<math>{\omega}_k^{\left(i\right)}=\prod_{j=1}^nU_{A_J}^{l_j}(x_j^{\left(i\right)})</math>                                                                                                                                               
|}
| style="width: 5px;text-align: right;white-space: nowrap;" |(17)
|}


{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;" 
|-
| style="text-align: center;" |<math>doc^l=1-\frac{{\sigma }^l}{{\max}_{k,t=1}^{N_h}\vert y_k^{\left(i\right)}-y_t^{\left(i\right)}\vert }</math> 
|}
| style="width: 5px;text-align: right;white-space: nowrap;" |(18)
|}

==3. Process and analysis=

===3.1 Data Acquisition and Cleaning===

The research data in this paper is based on the construction site of Subway Line 4 in Fuzhou City. Data is collected and saved to the database in real-time at intervals of one minute.

In the research process, the correct selection of factor data determines the correct match between the analysis object and the results. Due to the large number of excavation factor data generated by the shield machine during the excavation process, and the existence of missing values and outliers to some extent in each factor, it is necessary to clean the raw data to improve the quality of factor data.

The effectiveness of data cleaning directly affects the accuracy of subsequent data analysis and model predictions. Through investigation and monitoring of the construction site of Fuzhou Metro Line 4, based on traditional data preprocessing methods and fully considering factors such as the shield tunneling construction process, construction conditions, and the operator's habits, this article formulated the following five data cleaning principles:

(1) Based on detailed geological survey data, select shield tunneling data under the same geological conditions as the research object, excluding the influence of different geological conditions on the soil pressure caused by the shield machine.

(2) Due to the large amount of data, delete the real-time factor data corresponding to missing or abnormal values in each factor to ensure the correspondence between the data of each factor.

(3) Data corresponding to abnormal tunneling states such as stand-by, stoppage, cutterhead replacement, and mechanical failure are deleted to ensure that the data only includes those from normal tunneling processes of the shield machine.

(4) Data for monitoring factors such as horizontal deviation, front and rear mileage, mechanical and electrical parameters, and vertical deviation with values of 0 are removed.

(5) If the data of a certain factor is 0, the corresponding real-time data of other factors will be deleted to ensure the correspondence between each factor data.

After applying the above data cleaning principles, a total of 27 rings, 3525 sets of data, and 32 excavation factors were selected in the excavation factors data ranging from ring 407 to 434.

===3.2 Spearman correlation analysis===

In this paper, to investigate whether there is a relationship between the shield tunneling factors and the earth pressure in the cutterhead chamber, as well as the degree of correlation, it is necessary to perform a correlation analysis.Two commonly used methods for correlation analysis are the Pearson correlation coefficient and the Spearman correlation coefficient. The Pearson correlation coefficient belongs to parametric statistical methods, while the Spearman correlation coefficient belongs to non-parametric statistical methods.  Therefore, it is necessary to perform normality test on the 32 data factors obtained from the previous data cleaning. This article lists the histograms and Q-Q plots of the following three excavation factors, as shown in [[#img-2|Figures 2]] to [[#img-4|4]]. As shown in the figures, some factors follow normal distribution while others do not, indicating that the Spearman correlation coefficient should be chosen for analysis.

<div id='img-2'></div>
{| class="wikitable" style="margin: 0em auto 0.1em auto;border-collapse: collapse;width:auto;" 
|-style="background:white;"
|style="text-align: center;padding:10px;"| [[Image:Draft_Shao_683360558-image21-c.png|600px]]
|-
| style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 2'''. Boring pressure of Group B. (a) Histogram. (b) Q-Q plot
|}


<div id='img-3'></div>
{| class="wikitable" style="margin: 0em auto 0.1em auto;border-collapse: collapse;width:auto;" 
|-style="background:white;"
|style="text-align: center;padding:10px;"| [[Image:Draft_Shao_683360558-image22-c.png|600px]]
|-
| style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 3'''. Grouting pressure of grouting hole A3. (a) Histogram. (b) Q-Q plot
|}


<div id='img-4'></div>
{| class="wikitable" style="margin: 0em auto 0.1em auto;border-collapse: collapse;width:auto;" 
|-style="background:white;"
|style="text-align: center;padding:10px;"| [[Image:Draft_Shao_683360558-image23-c.png|600px]]
|-
| style="background:#efefef;text-align:left;padding:10px;font-size: 85%;"| '''Figure 4'''. Average oil pressure of the tail seal before the shield closure of the front chamber. (a) Histogram. (b) Q-Q plot
|}


The Spearman correlation coefficient, also known as the rank correlation coefficient, is calculated using the following Eq.(19):

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math> r=\frac{\sum_{i=1}^N(R_i-\overline{R})(S_i-\overline{S})}{{\left[\sum_{i=1}^N(R_i-\overline{R})^2\sum_{i=1}^N(S_i-\overline{S})^2\right]}^{\frac{1}{2}}}  </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" |(19)
|}

where <math>R_i</math> and <math>S_i</math> represent the rank order of the observed value <math>i</math>; <math>R</math> and <math>S</math> are the average rank orders of <math>x</math> and <math>y</math>, respectively; <math>N</math> is the total number of observed values.

Apart from the factor of silo pressure, this study conducted correlation analysis between 31 other factors and silo pressure using SPSS. The results are shown in [[#tab-1|Table 1]] and [[#img-5|Figure 5]]. The value of the Spearman correlation coefficient <math>r</math> ranges from -1 to 1, where values between 0 and 1 indicate a positive correlation and values between -1 and 0 indicate a negative correlation. The absolute value of the correlation coefficient <math>\vert r\vert</math>  represents the degree of correlation between the two factors, and the correlation strength can generally be judged and described using [[#tab-2|Table 2]] [23].

<div class="center" style="font-size: 75%;">'''Table 1'''.  Spearman correlation coefficients between tunneling factors and earth pressure for each excavation factor ID</div>

<div id='tab-1'></div>
{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;" 
|-style="text-align:center"
! Factor number !! Factor name !! Spearman correlation coefficient <math>r</math> !! Degree of correlation
|-
|  style="text-align: center;"|1
|  style="text-align: center;"|Advance speed（mm/min）
|  style="text-align: center;"|0.297**
|  style="text-align: center;"|Weak
|-
|  style="text-align: center;"|2
|  style="text-align: center;"|Cutter speed（r/min）
|  style="text-align: center;"|-0.239**
|  style="text-align: center;"|Weak
|-
|  style="text-align: center;"|3
|  style="text-align: center;"|Cutter Torque（kn.m）
|  style="text-align: center;"|0.378**
|  style="text-align: center;"|Weak
|-
|  style="text-align: center;"|4
|  style="text-align: center;"|Total thrust（kn）
|  style="text-align: center;"|0.683**
|  style="text-align: center;"|Strong
|-
|  style="text-align: center;"|5
|  style="text-align: center;"|Group A propulsion pressure（MPa）
|  style="text-align: center;"|0.07**
|  style="text-align: center;"|Unrelated
|-
|  style="text-align: center;"|6
|  style="text-align: center;"|Group B propulsion pressure（MPa）
|  style="text-align: center;"|0.532**
|  style="text-align: center;"|Moderate
|-
|  style="text-align: center;"|7
|  style="text-align: center;"|Group C propulsion pressure（MPa）
|  style="text-align: center;vertical-align: top;"|0.062**
|  style="text-align: center;vertical-align: top;"|Unrelated
|-
|  style="text-align: center;"|8
|  style="text-align: center;"|Group D propulsion pressure（MPa）
|  style="text-align: center;vertical-align: top;"|0.151**
|  style="text-align: center;vertical-align: top;"|Unrelated
|-
|  style="text-align: center;"|9
|  style="text-align: center;"|Group A advance displacement（mm/m）
|  style="text-align: center;vertical-align: top;"|-0.125**
|  style="text-align: center;vertical-align: top;"|Unrelated
|-
|  style="text-align: center;"|10
|  style="text-align: center;"|Group B advance displacement（mm/m）
|  style="text-align: center;vertical-align: top;"|-0.122**
|  style="text-align: center;vertical-align: top;"|Unrelated
|-
|  style="text-align: center;"|11
|  style="text-align: center;"|Group C advance displacement（mm/m）
|  style="text-align: center;vertical-align: top;"|-0.122**
|  style="text-align: center;vertical-align: top;"|Unrelated
|-
|  style="text-align: center;"|12
|  style="text-align: center;"|Group D advance displacement（mm/m）
|  style="text-align: center;vertical-align: top;"|-0.129**
|  style="text-align: center;vertical-align: top;"|Unrelated
|-
|  style="text-align: center;"|13
|  style="text-align: center;"|Grouting hole A1 grouting pressure（kn/m<sup>2</sup>）
|  style="text-align: center;vertical-align: top;"|0.33**
|  style="text-align: center;vertical-align: top;"|Weak
|-
|  style="text-align: center;"|14
|  style="text-align: center;"|Grouting hole A2 grouting pressure（kn/m<sup>2</sup>）
|  style="text-align: center;vertical-align: top;"|0.049**
|  style="text-align: center;vertical-align: top;"|Unrelated
|-
|  style="text-align: center;"|15
|  style="text-align: center;"|Grouting hole A3 grouting pressure（kn/m<sup>2</sup>）
|  style="text-align: center;vertical-align: top;"|0.33**
|  style="text-align: center;vertical-align: top;"|Weak
|-
|  style="text-align: center;"|16
|  style="text-align: center;"|Grouting hole A4 grouting pressure（kn/m<sup>2</sup>）
|  style="text-align: center;vertical-align: top;"|0.025**
|  style="text-align: center;vertical-align: top;"|Unrelated
|-
|  style="text-align: center;"|17
|  style="text-align: center;"|screw machine pressure（MPa）
|  style="text-align: center;vertical-align: top;"|0.232**
|  style="text-align: center;vertical-align: top;"|Weak
|-
|  style="text-align: center;"|18
|  style="text-align: center;"|Screw machine torque（KN<sup>2</sup>m）
|  style="text-align: center;vertical-align: top;"|0.232**
|  style="text-align: center;vertical-align: top;"|Weak
|-
|  style="text-align: center;"|19
|  style="text-align: center;"|Screw speed（rpm）
|  style="text-align: center;vertical-align: top;"|0.326**
|  style="text-align: center;vertical-align: top;"|Weak
|-
|  style="text-align: center;"|20
|  style="text-align: center;"|Spiral machine oil temperature（℃）
|  style="text-align: center;vertical-align: top;"|-0.188**
|  style="text-align: center;vertical-align: top;"|Unrelated
|-
|  style="text-align: center;"|21
|  style="text-align: center;"|The average oil pressure of the shield front and tail seals（MPa）
|  style="text-align: center;vertical-align: top;"|0.278**
|  style="text-align: center;vertical-align: top;"|Weak
|-
|  style="text-align: center;"|22
|  style="text-align: center;"|The average oil pressure of the shield tail seal of the sealed warehouse（MPa）
|  style="text-align: center;vertical-align: top;"|0.274**
|  style="text-align: center;vertical-align: top;"|Weak
|-
|  style="text-align: center;"|23
|  style="text-align: center;"|Average pressure of foam system（MPa）
|  style="text-align: center;vertical-align: top;"|0.464**
|  style="text-align: center;vertical-align: top;"|Moderate
|-
|  style="text-align: center;"|24
|  style="text-align: center;"|Foam Circuit Average Air Flow（L/min）
|  style="text-align: center;vertical-align: top;"|0.16**
|  style="text-align: center;vertical-align: top;"|Unrelated
|-
|  style="text-align: center;"|25
|  style="text-align: center;"|Foam Circuit Average Liquid Flow（L/min）
|  style="text-align: center;vertical-align: top;"|-0.021
|  style="text-align: center;vertical-align: top;"|Unrelated
|-
|  style="text-align: center;"|26
|  style="text-align: center;"|Articulation cylinder pressure（MPa）
|  style="text-align: center;vertical-align: top;"|-0.119**
|  style="text-align: center;vertical-align: top;"|Unrelated
|-
|  style="text-align: center;"|27
|  style="text-align: center;"|Group A articulated cylinder stroke（mm）
|  style="text-align: center;vertical-align: top;"|-0.153**
|  style="text-align: center;vertical-align: top;"|Unrelated
|-
|  style="text-align: center;"|28
|  style="text-align: center;"|Group B articulated cylinder stroke（mm）
|  style="text-align: center;vertical-align: top;"|-0.211**
|  style="text-align: center;vertical-align: top;"|Weak
|-
|  style="text-align: center;"|29
|  style="text-align: center;"|Group C articulated cylinder stroke（mm）
|  style="text-align: center;vertical-align: top;"|0.165**
|  style="text-align: center;vertical-align: top;"|Unrelated
|-
|  style="text-align: center;"|30
|  style="text-align: center;"|Group D articulated cylinder stroke（mm）
|  style="text-align: center;vertical-align: top;"|0.332**
|  style="text-align: center;vertical-align: top;"|Weak
|-
|  style="text-align: center;"|31
|  style="text-align: center;"|Foam water real-time cumulative average（m<sup>3</sup>/h）
|  style="text-align: center;vertical-align: top;"|-0.104**
|  style="text-align: center;vertical-align: top;"|Unrelated
|}

**<nowiki>**.exhibits significant correlation at the 0.01 level (two-tailed)</nowiki>


<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Table 2 |r| Values and Corresponding Degrees of Correlation</div>

{| style="width: 100%;border-collapse: collapse;" 
|-
|<nowiki>  |r|</nowiki>||r|correlation coefficient
|-
|  style="border-top: 1pt solid black;"|(0,0.2]
|  style="border-top: 1pt solid black;"|Extremely Weak Correlation

/Unrelated
|-
| (0.2,0.4]
| Weak Correlation
|-
| (0.4,0.6]
| Moderate Correlation
|-
| (0.6,0.8]
| Strong Correlation
|-
|  style="border-bottom: 2pt solid black;"|(0.8,1]
|  style="border-bottom: 2pt solid black;"|Extremely Strong Correlation
|}


<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Shao_683360558-image25-c.png|270px]] </div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Fig.5 Results of Spearman correlation analysis</div>

The correlation between total thrust and soil pressure is strongly correlated. The average pressure of the foam system and the push pressure of group B have a moderate correlation. The correlation between cutterhead torque, advance speed, cutterhead speed, grouting hole A1 grouting pressure, grouting hole A3 grouting pressure, screw conveyor pressure, screw conveyor torque, screw conveyor speed, average oil pressure of the front shield tail seal, average oil pressure of the rear shield tail seal, the stroke of the hinge oil cylinder of group B and group D have a weak correlation with soil pressure. Other factors have a very weak correlation and can be considered to have no correlation with soil pressure.

Therefore, this study finally determined 15 excavation factors including total thrust, B-group advance pressure, foam system average pressure, cutterhead torque, advance speed, cutterhead speed, grouting hole A1 grouting pressure, grouting hole A3 grouting pressure, screw conveyor pressure, screw conveyor torque, screw conveyor speed, average oil pressure of shield tail sealing at front chamber, average oil pressure of shield tail sealing at rear chamber, B-group articulated oil cylinder stroke, and D-group articulated oil cylinder stroke as inputs for the neural network prediction model.

==3.3 Evaluation Indicators==

The evaluation indicators selected in this paper are mean absolute error (''M''<sub>AE</sub>), mean squared error (''M''<sub>SE</sub>), and coefficient of determination (R<sup>2</sup>).The mean absolute error (''M''<sub>AE</sub>) is used in this paper to reflect the error between the actual and predicted values, and its Eq. is shown in Eq.20. The mean squared error (''M''<sub>SE</sub>) is used to evaluate the degree of data variation, and the smaller the value, the better the predictive model's ability to accurately describe the data, as shown in Eq.21 [24].The coefficient of determination R<sup>2</sup> is used to evaluate the goodness-of-fit of a regression model. The closer its value is to 1, the better the performance of the model, as shown in Eq.22.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 <math display="inline">M_{\mbox{AE}}=\frac{1}{n}\sum_{i=1}^n\vert {\overset{\mbox{ˆ}}{y}}_l-y_i\vert</math> （20）</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 <math display="inline">M_{\mbox{SE}}=\frac{1}{n}\sum_{i=1}^n(y_i-{\overset{\mbox{ˆ}}{y}}_l)^2</math> （21）</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 <math display="inline">R^2=1-\frac{\sum_i{\left({\overline{y}}_i-{\overset{\mbox{ˆ}}{y}}_i\right)}^2}{\sum_i{\left(y_i-{\overline{y}}_i\right)}^2}</math> （22）</div>

In the Eq., ''n ''represents the sample size; ''ŷ<sub>i</sub>'' represents the predicted value; ''y<sub>i</sub>'' represents the actual value; and ''ȳ<sub>i</sub>'' represents the mean value.

==3.4 Model training==

This paper builds an LSTM deep neural network prediction model using PyTorch deep learning framework, based on mini-batch gradient descent method and further optimized by Adam algorithm. In order to improve the accuracy of the LSTM prediction model, PSO is used to optimize LSTM parameters, and a PSO-BP prediction model is constructed using the classical BP neural network for comparison.The loss function of the prediction model is the root mean square error (RMSE) loss function. The training process is iteratively trained using the training set, and the prediction results are compared and analyzed using the test set. The hyperparameters of the four comparative models are tuned by parameter tuning, as shown in Table 3.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Table 3 Parameter Tuning of Four Comparative Models</div>

{| style="width: 76%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|Model parameters
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|learning_rate
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|lstm_layers
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|hidden_size
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|batch_size
|-
|  style="border-top: 1pt solid black;text-align: center;"|LSTM
|  style="border-top: 1pt solid black;text-align: center;"|0.001
|  style="border-top: 1pt solid black;text-align: center;"|3
|  style="border-top: 1pt solid black;text-align: center;"|52
|  style="border-top: 1pt solid black;text-align: center;"|612
|-
|  style="text-align: center;"|BP
|  style="text-align: center;"|0.001
|  style="text-align: center;"|2
|  style="text-align: center;"|62
|  style="text-align: center;"|640
|-
|  style="text-align: center;"|PSO-LSTM
|  style="text-align: center;"|0.005
|  style="text-align: center;"|3
|  style="text-align: center;"|46
|  style="text-align: center;"|1024
|-
|  style="border-bottom: 2pt solid black;text-align: center;"|PSO-BP
|  style="border-bottom: 2pt solid black;text-align: center;"|0.001
|  style="border-bottom: 2pt solid black;text-align: center;"|3
|  style="border-bottom: 2pt solid black;text-align: center;"|48
|  style="border-bottom: 2pt solid black;text-align: center;"|820
|}


<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Shao_683360558-image29-c.png|270px]] </div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Fig.6 Graphs of the Four Comparative Models</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Table 4. Comparison results</div>

{| style="width: 79%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|test group
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|''M''<sub>AE</sub>(MPa)
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|''M''<sub>SE</sub>(%)
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|R<sup>2</sup>
|-
|  style="border-top: 1pt solid black;border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|PSO-BP
|  style="border-top: 1pt solid black;border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|0.010
|  style="border-top: 1pt solid black;border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|0.018
|  style="border-top: 1pt solid black;border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|0.84
|-
|  style="border-top: 2pt solid black;border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|PSO-LSTM
|  style="border-top: 2pt solid black;border-bottom: 2pt solid black;text-align: center;"|0.007
|  style="border-top: 2pt solid black;border-bottom: 2pt solid black;text-align: center;"|0.007
|  style="border-top: 2pt solid black;border-bottom: 2pt solid black;text-align: center;"|0.93
|-
|  style="border-top: 2pt solid black;border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|LSTM
|  style="border-top: 2pt solid black;border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|0.011
|  style="border-top: 2pt solid black;border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|0.015
|  style="border-top: 2pt solid black;border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|0.87
|-
|  style="border-top: 2pt solid black;border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|BP
|  style="border-top: 2pt solid black;border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|0.015
|  style="border-top: 2pt solid black;border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|0.028
|  style="border-top: 2pt solid black;border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|0.75
|}


According to Figure 6 and Table 4, it can be seen that compared to the BP prediction model, the LSTM prediction model has reduced ''M''<sub>AE</sub> from 0.015 MPa to 0.011 MPa by 26.7%, reduced ''M''<sub>SE</sub> from 0.028% to 0.015% by 46.4%, and increased the absolute coefficient R<sup>2</sup> from 0.75 to 0.87 by 16%. It can be seen that the LSTM prediction model is better than the BP prediction model in predicting the pressure of the soil bin.Compared with PSO-BP, PSO-LSTM reduced ''M''<sub>AE</sub> from 0.010 MPa to 0.007 MPa, a decrease of 30%, reduced ''M''<sub>SE</sub> from 0.018% to 0.007%, a decrease of 61.1%, and increased R<sup>2</sup> from 0.84 to 0.93, an improvement of 10.7%. It can be seen that PSO-optimized LSTM prediction model performs better than PSO-optimized BP prediction model in predicting the pressure in the soil bin.Compared with the LSTM prediction model, the PSO-optimized model has reduced ''M''<sub>AE</sub> from 0.011 MPa to 0.007 MPa, a decrease of 36.4%, ''M''<sub>SE</sub> from 0.015% to 0.007%, a decrease of 53.3%, and R<sup>2</sup> from 0.87 to 0.93, an increase of 6.9%, indicating that using PSO to optimize the LSTM prediction model can improve its prediction accuracy.

==3.5 Fuzzy control==

In order to achieve the application of intelligent control of shield machines in practical engineering, this paper needs to obtain the corresponding relationship between the soil pressure and the screw machine speed. As the six factors of advancing speed, cutterhead speed, total thrust, cutterhead torque, screw machine speed, and soil pressure can be obtained from the monitoring room of a shield machine, historical data of these six factors from ring 407 to 433 were selected as the training set, and the data from ring 434 was used as the test set to construct a 5-input single-output fuzzy control rule.

The maximum and minimum values of each factor are taken as the domain range, and the triangular membership function is used to automatically partition the fuzzy set based on the specified number of partition intervals for each factor. The data is then fuzzified, and the membership degree of each data point on the fuzzy set is calculated. The highest membership degree value is selected as the basis for a fuzzy rule. Due to the large number of repeated rules in these rules, their confidence levels are calculated separately and filtered to establish a fuzzy rule base.

The predicted value of the 434th ring shield pressure obtained by the LSTM prediction model, together with the real-time read push speed, cutterhead speed, total thrust, and cutterhead torque that can be obtained from the shield machine monitoring room, are used as inputs to the fuzzy rules. The corresponding target value of the screw conveyor speed is obtained, as shown in Figure 7.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Shao_683360558-image30-c.png|288px]] </div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Fig.7 Graph of the variation in screw machine speed</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Table 5 Evaluation Indicators of WM Calculated Values and Motor Operation Values</div>

{| style="width: 79%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|Evaluation index
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|''M''<sub>AE</sub>(rpm)
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|''M''<sub>SE</sub>(%)
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|R<sup>2</sup>
|-
|  style="border-top: 1pt solid black;border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|Fuzzy rule calculation value
|  style="border-top: 1pt solid black;border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|0.065
|  style="border-top: 1pt solid black;border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|1.528
|  style="border-top: 1pt solid black;border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|0.82
|}


According to Table 5, the errors between the extracted fuzzy rules for screw machine speed and the actual values, obtained through the WM algorithm from historical data, are 0.065rpm for ''M''<sub>AE</sub>, 1.528% for ''M''<sub>SE</sub>, and 0.82 for absolute coefficient R<sup>2</sup>, indicating that the accuracy of the fuzzy rule base established by the algorithm is relatively high and basically meets the requirements.

==3.6 intelligent control device==

To achieve precise control of the screw conveyor, it is essentially precise control of the drive motor. Based on the actual situation of the Fuzhou Metro Line 4 shield machine control room, this paper uses an Arduino main control board, Simple FOC to drive a BLDC brushless DC motor, and an AS5600 magnetic encoder as the motor's position sensor to create an intelligent control device for the shield screw machine control knob, as shown in Figure 8.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Shao_683360558-image31-c.png|438px]] </div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Fig.8 Intelligent Control Device for the Screw Machine Speed Control Knob.</div>

The working principle of the intelligent control device for the shield spiral conveyor is as follows: When the PC sends a position target command to the Arduino main control board through USB serial communication, the Simple FOC drives the BLDC brushless DC motor to start. The motor gear fixed on the top of the BLDC brushless DC motor rotates, and through the gear meshing with the knob gear, the motor torque is transmitted to the knob gear, thereby driving the speed knob of the spiral conveyor fixed by screws under the knob gear. The AS5600 magnetic sensor obtains the position data of the BLDC brushless DC motor rotation by fixing a radial magnet on the top of the cylindrical center of the motor gear, feeds it back to the Arduino main control board through I2C and transmits it to the PC through USB serial communication. The brushless DC motor stops rotating when it reaches the target position.

The intelligent control device established in this article mainly includes two parts: hardware and software. The key to establishing the overall control program is to realize the communication between Python and Arduino microcontroller. The structure of the overall control program established in this article is shown in Figure 9, and the control flow of the overall program is shown in Figure 10.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Shao_683360558-image32-c.png|402px]] </div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Fig.9 Structure of the Central Control Program</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Shao_683360558-image33-c.png|600px]] </div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Fig.10 Overall program control flow</div>

The performance of the intelligent control device was validated using actual data from the 434-ring earth pressure balance shield machine of Fuzhou Metro Line 4. The predicted silo pressure values obtained from the PSO-LSTM model were used as inputs to the fuzzy rules, and the corresponding spiral machine speed target control values were obtained. The results are shown in Figure 11.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
 [[Image:Draft_Shao_683360558-image34-c.png|372px]] </div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Fig.11 Comparison between control target instructions and actual spiral machine speed and feedback values of the control device</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Table 6 Evaluation indicators of the comparison between the control target instruction, device operation feedback value, and actual engineering data</div>

{| style="width: 79%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|contrast
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|''M''<sub>AE</sub>(rpm)
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|''M''<sub>SE</sub>(%)
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|R<sup>2</sup>
|-
|  style="border-top: 1pt solid black;border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|Control target command value
|  style="border-top: 1pt solid black;border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|0.065
|  style="border-top: 1pt solid black;border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|1.528
|  style="border-top: 1pt solid black;border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|0.82
|-
|  style="border-top: 2pt solid black;border-bottom: 2pt solid black;text-align: center;vertical-align: top;"|Device operation feedback value
|  style="border-top: 2pt solid black;border-bottom: 2pt solid black;text-align: center;"|0.015
|  style="border-top: 2pt solid black;border-bottom: 2pt solid black;text-align: center;"|0.392
|  style="border-top: 2pt solid black;border-bottom: 2pt solid black;text-align: center;"|0.95
|}


According to Table 6, the ''M''<sub>AE</sub> and ''M''<sub>SE</sub> errors of the fuzzy rules calculated from the WM algorithm, based on the PSO-LSTM model's predicted value of silo pressure, and the corresponding spiral speed value compared with the actual spiral speed value are 0.065 rpm and 1.528%, respectively. The absolute coefficient R<sup>2</sup> can reach 0.82, indicating that the calculation effect is good and basically meets the requirements. Comparing the real-time data of the motor's rotation position feedback by the intelligent control device's magnetic encoder and the control target instruction value, there are ''M''<sub>AE</sub> and ''M''<sub>SE</sub> errors of 0.015 rpm and 0.392%, respectively, between the effect executed by the intelligent control device and the target input instruction. The absolute coefficient R<sup>2</sup> is 0.95, indicating that the execution error of the device is small. Therefore, the feasibility and reliability of the device have been verified.

=4. Conclusion=

This paper relies on the construction site data of Fuzhou Metro Line 4 in Fujian Province, China, and proposes a soil bin pressure prediction model combining LSTM and PSO. Its mean absolute error, mean squared error, and coefficient of determination values are 0.007MPa, 0.007%, and 0.93, respectively, which show improved accuracy compared to other prediction models.Meanwhile, the Wang-Mendel algorithm was used to extract fuzzy rules from historical data. The mean absolute error and mean square error of the rotating speed of the screw machine were 0.065rpm and 1.528%, and the coefficient of determination was 0.82, indicating that the fuzzy rule base established by the algorithm had high calculation accuracy.The development of an intelligent control device aimed at controlling the screw machine speed knob showed the mean absolute error and mean square error of 0.015rpm and 0.392%, respectively, with the absolute coefficient of 0.95, indicating that the device had low execution error. The research content of this paper is reliable and provides a new and effective method for controlling the EPB soil pressure.

==Acknowledgments==

This work was supported by the Natural Science Foundation of Fujian Province (Grant No.2022J01390), the China Postdoctoral Science Foundation (Grant No.2020M671956) ,Fujian Province Foreign Cooperation Industrialization Project(Grant No.2021I1006) and Fujian Province technology innovation key research and industrialization projects(Grant No.2023XQ002）

==References==

<span id='_Ref130840916'></span>[[#cite-_Ref130840916|[1]]] HU Guo-liang,GONG Guo-fang,YANG Hua-yong.Realization of earth pressure balance for shield tunnelling machine[J].JOURNAL OF ZHEJIANG UNIVERSITY (ENGINEERING SCIENCE),2006(05):874-877.(in Chinese)

<span id='_Ref130840930'></span>[[#cite-_Ref130840930|[2]]] Wang Hongxin,Fu Deming.A mathematical model and the related parameters for EPB shield tunneling[J].CHINA CIVIL ENGINEERING JOURNAL,2006(09):86-90.(in Chinese)

<span id='_Ref130840987'></span>[[#cite-_Ref130840987|[3]]] WANG Lintao,GONG Guofang,YANG Huayong, HOU Dianqing.Earth pressure balance control based on feedforward-feedback compound control[J].Journal of Central South University (Science and Technology),2013,44(07):2726-2735.(in Chinese)

<span id='_Ref130841013'></span>[[#cite-_Ref130841013|[4]]] Yeh I C . Application of neural networks to automatic soil pressure balance control for shield tunneling[J]. Automation in Construction, 1997, 5(5):421-426.

[5] Wei, M., Wang, Z., Wang, X. et al. Prediction of TBM penetration rate based on Monte Carlo-BP neural network. Neural Comput & Applic 33, 603–611 (2021).

[6] Xianjie Gao, Maolin Shi, Xueguan Song, Chao Zhang, Hongwei Zhang,Recurrent neural networks for real-time prediction of TBM operating parameters,Automation in Construction,Volume 98,2019,Pages 225-235,ISSN [0926-5805，https://doi.org/10.1016/j.autcon.2018.11.013. 0926-5805，https://doi.org/10.1016/j.autcon.2018.11.013.]

[7] Zaobao Liu, Long Li, Xingli Fang, Wenbiao Qi, Jimei Shen, Hongyuan Zhou, Yulong Zhang,Hard-rock tunnel lithology prediction with TBM construction big data using a global-attention-mechanism-based LSTM network,Automation in Construction,Volume 125,2021,103647,ISSN 0926-5805,[https://doi.org/10.1016/j.autcon.2021.103647. https://doi.org/10.1016/j.autcon.2021.103647.]

[8] Cachim, P., Bezuijen, A. Modelling the Torque with Artificial Neural Networks on a Tunnel Boring Machine. KSCE J Civ Eng 23, 4529–4537 (2019). [https://doi.org/10.1007/s12205-019-0302-0 https://doi.org/10.1007/s12205-019-0302-0].

[9] Xueguan Song, Maolin Shi, Jianguo Wu, Wei Sun,A new fuzzy c-means clustering-based time series segmentation approach and its application on tunnel boring machine analysis,Mechanical Systems and Signal Processing,Volume 133,2019,106279,ISSN 0888-3270,[https://doi.org/10.1016/j.ymssp.2019.106279. https://doi.org/10.1016/j.ymssp.2019.106279.]

[10] Maolin Shi, Tianci Zhang, Liyong Zhang, Wei Sun, Xueguan Song,A fuzzy c-means algorithm based on the relationship among attributes of data and its application in tunnel boring machine,Knowledge-Based Systems,Volume 191,2020,105229,ISSN 0950-7051.

[11] Wei Sun, Maolin Shi, Chao Zhang, Junhong Zhao, Xueguan Song,Dynamic load prediction of tunnel boring machine (TBM) based on heterogeneous in-situ data,Automation in Construction,Volume 92,2018,Pages 23-34,ISSN 0926-5805,[https://doi.org/10.1016/j.autcon.2018.03.030. https://doi.org/10.1016/j.autcon.2018.03.030.]

[12] Wen, Z., Wang, Z., Rong, X. et al. Prediction of the Amount of Soil Discharged by an Earth Pressure Balanced Shield Machine Based on Feature Engineering. KSCE J Civ Eng 25, 4868–4886 (2021).

<span id='_Ref130840963'></span>[[#cite-_Ref130840963|[13]]]<nowiki> Liu X , Xu S , Huang Y . Optimal control for earth pressure balance of shield machine based on action-dependent heuristic dynamic programming, ISA Transactions ,Volume 94,2019,Pages 28-35,ISSN 0019-0578,https://doi.org/10.1016/j.isatra.2019.04.007.</nowiki>

<span id='_Ref130841043'></span>[[#cite-_Ref130841043|[14]]]<nowiki> Zhang Yan. Multi-point Earth Pressure Prediction and Optimal Control Method for EPB Shield Capsule[D].Dalian University of Technology,2011.(in Chinese)，https://kns.cnki.net/kcms2/article/abstract?v=A67obWUfw0lxTK8BeOwngvPL3MO_OuZ_T4CcFElsylpkirzW2IvOdmfAD7ItlhsZVvGzKIL80xofvgBmjQKlbTXqRHxg6-CbC1-oyJ5P5S0LbPccE9PpoWB3ePDxgo42&uniplatform=NZKPT&language=CHS</nowiki>

<span id='_Ref130841054'></span>[[#cite-_Ref130841054|[15]]] ZHU Hong-tang.Research of Earth Pressure Balance Control of Shield Machine Based on the Fuzzy PID Algorithm[J].CHINESE HYDRAULICS & PNEUMATICS,2013,No.262(06):103-107.(in Chinese)

<span id='_Ref130841087'></span>[[#cite-_Ref130841087|[16]]] Xiaofei Li, Guofang Gong,Predictive control of slurry pressure balance in shield tunneling using diagonal recurrent neural network and evolved particle swarm optimization,Automation in Construction,Volume 107,2019,102928,ISSN 0926-5805.

[17] Li Qiang,Li Jingxuan.Chamber earth pressure balance control for shield machine based on PSO-BP prediction model[J].INFORMATION TECHNOLOGY AND NETWORK SECURITY,2018,37(04):94-99.DOI:10.19358/j.issn.2096-5133.2018.04.022.(in Chinese)

<span id='_Ref130841098'></span>[[#cite-_Ref130841098|[18]]] Zhang Chao,Guo Jingbo,Zhang Haidong,Zhou Qingxiang.Research of Earth Pressure Balance Control Based on the Fuzzy PID[J].JOURNALOF' SHIJIAZHUANG TIEDAO UNIVERSITY(NaturalScienceEdition),2019,32(03):82-88.DOI:10.13319/j.cnki.sjztddxxbzrb.20180024. (in Chinese)

<span id='_Ref130841149'></span>[[#cite-_Ref130841149|[19]]] Yang Xiaoda. Intelligent Prediction of In-chamber Earth Pressure for EPB Shield and Optimization of Tunneling Parameters Based on Machine Learning[D].Shandong University,2022.DOI:10.27272/d.cnki.gshdu.2022.002796.(in Chinese)

<span id='_Ref130841744'></span>[[#cite-_Ref130841744|[20]]] L.X. Wang, J.M. Mendel, Generating fuzzy rules by learning from examples,IEEE Trans. Syst. Man Cybern. 22 (6) (1992) 1414–1427.

<span id='_Ref130841756'></span>[[#cite-_Ref130841756|[21]]] SHEN Jun，LI Shao-yuan.Extraction and application of the fuzzy control rules based on Wang-Mendel fuzzy model[J].Journal of Systems Englneering and Electronics,2007,No.336(09):1524-1528.(in Chinese)

<span id='_Ref130841766'></span>[[#cite-_Ref130841766|[22]]] Y. Wang, D. Wang, T. Chai, Extract of fuzzy rules with completeness and robustness, Acta Autom. Sin. 36 (9) (2010) 1337–1342 (in Chinese).

<span id='_Ref130840836'></span>[[#cite-_Ref130840836|[23]]] Ll Kun,TIAN Guan-feng,MA Hong-wei,YANG Xiao-lin,PU Bao-jun,JIANG Yun-chen.Correlation Analysis and Prediction Model of Boring Parameters of Earth Pressure Balance Shield[J].Science Technology and Engineering,2021,21(09):3814-3821.(in Chinese)

<span id='_Ref130841576'></span>[[#cite-_Ref130841576|[24]]] YANG Kaihong,ZHANG Qian,ZHOU Siyang.Advance Rate Prediction by Recurrent Neural Networks Analysis on In-situ Data of Tunnel Boring Machines[J].MECHANICAL SCIENCE AND TECHNOLOGY FOR AEROSPACE ENGINEERING，2021，40（6）：835-839.(in Chinese)