Two-Dimensional, Two-Layer Quality Regression Model Based Batch Process Monitoring

In this paper, a two-dimensional, two-layer quality regression model is established to monitor multi-phase, multi-mode batch processes. Firstly, aiming at the multi-phase problem and the multi-mode problem simultaneously, the relations among modes and phases are captured through the analysis between process variables and quality variables by establishing a two-dimensional, two-layer regression partial least squares (PLS) model. The two-dimensional regression traces the intra-batch and inter-batch characteristics, while the two-layer structure establishes the relationship between the target process and historical modes and phases. Consequently, online monitoring is carried out for multi-phase, multi-mode batch processes based on quality prediction. In addition, the online quality prediction and monitoring results based on the proposed method and those based on the traditional phase mean PLS method are compared to prove the effectiveness of the proposed method.


Introduction
Batch process is a way of production closely related to people's life in the modern process industry. It is now widely used in the production and preparation of small-batch and high value-added products to meet the rapidly changing market demand, such as in the production of fine chemical industry, food, polymer reaction, metal processing, biopharmaceutical, etc. Batch operation process characteristics are more complex and have richer data statistical characteristics than the continuous industrial process. Firstly, a batch cycle usually consists of several fixed phases, and different phases in each batch may show different potential behaviors, which is called multi-phase characteristic. Multi-phase characteristic is one kind of intra-batch characteristic evolution; it is one characteristic evolution along the time direction within each batch. Secondly, in the process of industrial production, due to the change of operating conditions, there are many different stable working points. Certain number of batches within the same working points belong to the same mode, which is called multi-mode characteristic. Multi-mode characteristic is one kind of inter-batch characteristic evolution; it is one characteristic evolution along the batch direction through the process. Multi-phase characteristic and multi-mode characteristic exist in batch processes simultaneously, which makes the batch processes complex and interesting for researchers.
At the end of the 1980s, principal component analysis (PCA) [1] and partial least squares (PLS) [2], as the focus of multivariate statistical modeling methods, began to be applied to statistical modeling, online monitoring, quality control, and fault diagnosis based on process data [3][4][5]. In order to maximize their effectiveness, it is necessary to extend these methods to batch processes. In the mid-1990s, Nomikos and MacGregor proposed multilinear principal component analysis (MPCA) and multiway partial least squares regression (MPLS) [6,7]. After, many international research groups invested a lot of manpower and material resources to carry out the research work for batch processes. In 1998, Wold et al. described an approach to multivariate batch process modeling and comprehensive explanation of the process characteristics under the two irrelevant methods. In addition, when the methods in both directions need to be jointed, the method structure is relatively complex and difficult to comprehend and implement. Therefore, it is necessary to further improve the existing methods by proposing a novel united strategy framework according to the two-dimensional evolution characteristics of the batch processes. That is, it is necessary to study the method of solving the process evolution in both directions from the overall perspective of the two-dimensional evolution of batch processes. Based on the above analysis, the research focus in this paper is to establish the overall framework of tracking the two-dimensional evolution based on a two-dimensional regression model. Moreover, the basic one-dimensional regression model used in the above framework should be selected appropriately to have the ability to establish the relationship between the target process and historical process (modes or phases). The between-mode quality analysis [20], which has a two-layer mode analysis structure, is adopted, and the object is changed to processes (modes or phases), rather than modes. By the two-layer regression analysis strategy, the relationship between the target process and historical processes will be established.
Therefore, in this work, by expanding two-layer regression analysis strategy from dealing with one direction process evolution into two direction process evolutions under the two-dimensional framework, both the intra-batch and inter-batch characteristics will be traced for multi-phase and multi-mode processes, based on which a process monitoring strategy is developed. The two-dimensional regression traces the intra-batch and interbatch characteristics, while the two-layer structure establishes the relationship between the target process and historical modes and phases. Firstly, aiming at the multi-mode problem and the multi-phase problem simultaneously, the relations among modes and phases are captured through the analysis of the characteristics of each mode and each phase in the batch process by establishing a two-dimensional regression model. This model is expanded from the two-layer between-mode analyzing model which can extract useful quality-related information from historical modes. In the first layer of relation analysis, for all historical processes, different quality regression models are developed by PLS and according regression parameters are obtained. In addition, by applying these regression parameters to the process variables of the target process, a series of assumed quality predictions would be obtained. Then, in the second layer, PLS is conducted between the assumed quality predictions and real qualities of the target process to judge the relationship between the target process and the historical processes. Based on the two-dimensional, two-layer regression model, for the target process, quality prediction can be conducted using all available quality information underling both the intra-batch and inter-batch evolutions. Consequently, online monitoring is carried out for multi-phase, multi-mode batch processes based on quality prediction. In addition, the online quality prediction and monitoring results based on the proposed method and those based on the traditional method are compared to prove the effectiveness of the proposed method. In the traditional strategy [16,23], historical modes are considered together and modeled by PLS method to include more mode information and avoid the overfitting problem. However, for the multi-phase problem, several phase mean models are built for different phases, respectively. To be fair, all concerned modes are available for modeling in both the proposed method and the traditional method.
The remaining work of this paper includes the following aspects: firstly, the Section 2 introduces the establishment of two-dimensional, two-layer regression model and the principle of online monitoring. The Section 3 briefly introduces the characteristics of the injection molding process, and the application of the proposed method is illustrated by the online monitoring of the injection molding process. The Section 4 is conclusions.

Batch Process Monitoring Based on Phase Mean PLS
Generally, process data of batch processes are stored in data matrix X(I × J x × K), where I refers to the number of batches; J x refers to the number of process variables; and K where K c is the data length of phase c, and X c,k is the time-slice process data matrix. The PLS algorithm is used to decompose the average process data matrix X c and quality data matrix Y linearly, and the model is as follows: where T c is the score matrix of X c , U c is the score matrix of Y, P c is the load matrix of X c , Q c is the load matrix of Y, E c is the fitting error matrix of X c , and F c is the fitting error matrix of Y.
The previous model can be expressed by the regression model as: where B c is the regression parameter matrix, andŶ c is the predicted quality. When considering a single quality variable y(I × 1), the regression model can be simply expressed as: where β c is the regression parameter.
In this paper, PLS model is used for online monitoring. The commonly used statistics are Hotelling-T 2 statistics and square prediction error SPE statistics. Hotelling-T 2 statistics reflects the deviation degree of latent variables from the established model in amplitude and process data development trend. SPE describes the deviation degree of the measured value of the input variable from the latent variable space in the batch process [13].
In the online monitoring of batch process, the current operating condition can be determined by observing the Hotelling-T 2 and the control limit of the square prediction error SPE. The definitions of T 2 statistics and SPE statistics are as follows: where x c is the residual vector, R c = W c (P T c W c ) −1 , and W c is the weight matrix. The detailed properties and calculations can be found in reference [8].
The corresponding control limits are: where After obtaining the x c,k vector at the k-th moment, the online T 2 statistics and online SPE statistics are calculated, and the calculation formula is as follows: where x c,k is the residual vector at the k-th moment.

Framework of Two-Dimensional, Two-Layer Regression Modeling Strategy
To perform process monitoring based on the three-dimensional data matrix of the batch process, in this paper, a novel two-dimensional, two-layer regression modeling strategy is established, where the process evolution that exists in both directions within and between batches can be tracked, and important quality information in historical processes can be extracted.
In the process of batch process modeling and online application, the most direct embodiment of the process characteristics is the data, and the main characteristic of the batch process is the three-dimensional data matrix. Extraction and analysis of process characteristics from the data is the basis of subsequent process modeling and process analysis. On the basis of understanding the characteristics of the actual three-dimensional data matrix of the batch process, appropriate scientific methods should be adopted to extract the process evolution information intelligently and accurately, so as to deeply understand the characteristics of the process and further establish the process model. The analysis of process evolution characteristics reflects the respective characteristics in both intra-batch direction and inter-batch direction but also the consistency of the two directions from the three-dimensional data matrix.
The intra-batch evolution problem and inter-batch evolution problem are shown in Figure 1. In early work, the researchers proposed different methods to track the intrabatch evolution and inter-batch evolution, respectively. In general, these methods are individual and can be selected simultaneously in evolution analysis and used for modeling in one of the two directions, as shown in Figure 1a. However, it is relatively complex to choose different methods when processing jointly in both directions. In addition, the most important is that the evolution nature in the both directions is the same. Although the evolutions in the two directions may be raised because of different physical significances, they are both process characteristic evolutions. It is a straightforward idea to adopt one method to handle both problems in the two directions. Therefore, a thinking framework is proposed, that is, expanding the method to inter-batch evolution based on intra-batch evolution, and maintaining the unity of the method under a reasonable premise, so as to achieve the method simplification. It should be noted that the model building method should follow the physical significances of both the inter-batch evolution and the intrabatch evolution during application. Following this framework, the proposed methods in both directions are consistent in core ideas and will be more convenient when uniting.
Secondly, from the overall perspective of batch process three-dimensional data matrix, the other kind of method can solve the process evolution in two directions simultaneously. As shown in Figure 1b, this kind of method can be split into two directions, and performs under the thought of handling the whole first, and dealing with the parts later. This kind of strategy adopts a unified core algorithm to extract and analyze the evolution in two directions. The core method here should be as suitable for the tracking of process evolution characteristics as much as possible. Attention is paid to the use of intelligent analysis method to avoid computational redundancy and reflect the superiority of the algorithm starting from the 3-dimensional data matrix.
neously, as shown in Figure 1a, while, for the strategy shown in Figure 1b, the research will be talked about in later works. Next, the modeling method based on two-layer multi mode regression analysis will be expanded to deal with the multi-mode and multi-phase problem at the same time, and a two-dimensional, two-layer regression model will be proposed.
(a) Unity of methods for intra-batch evolution and inter-batch evolution  In the multi-mode batch process, when a new mode occurs, the single mode mode will no longer match the target process, resulting in serious overrun of T 2 and SPE statis tics when monitoring. Therefore, the statistical modeling and online monitoring of the batch process should not only consider a single mode but also carry out the quality anal ysis of multiple modes.
The main strategy to deal with the multi-mode problem is to extract the relationship between the historical modes and the target process. By doing this, more between-mode information would be captured by the model, and the multi-mode process can be better predicted and monitored. In this two-layer regression strategy, first, the PLS regression models are built for historical modes. Then, using those regression models and the process variables of the target process, a series of assumed quality predictions can be obtained in the first regression layer. In the second layer, based on these assumed quality predictions and real quality values of the target process, the PLS regression model is established fo quality prediction. Thus, this work is to explore more effective data modeling methods on the basis of preliminary work and trace the intra-batch evolution and inter-batch evolution simultaneously, as shown in Figure 1a, while, for the strategy shown in Figure 1b, the research will be talked about in later works. Next, the modeling method based on two-layer multi-mode regression analysis will be expanded to deal with the multi-mode and multi-phase problem at the same time, and a two-dimensional, two-layer regression model will be proposed.
In the multi-mode batch process, when a new mode occurs, the single mode model will no longer match the target process, resulting in serious overrun of T 2 and SPE statistics when monitoring. Therefore, the statistical modeling and online monitoring of the batch process should not only consider a single mode but also carry out the quality analysis of multiple modes.
The main strategy to deal with the multi-mode problem is to extract the relationship between the historical modes and the target process. By doing this, more between-mode information would be captured by the model, and the multi-mode process can be better predicted and monitored. In this two-layer regression strategy, first, the PLS regression models are built for historical modes. Then, using those regression models and the process variables of the target process, a series of assumed quality predictions can be obtained in the first regression layer. In the second layer, based on these assumed quality predictions and real quality values of the target process, the PLS regression model is established for quality prediction.
To adopt the same strategy to handle the multi-mode problem and the multi-phase problem at the same time, the strategy above will be applied to the intra-batch batch direction and handle the multi-phase problem by analyzing the intra-batch evolution. The framework of the two-dimensional, two-layer regression modeling strategy is shown in Figure 2. In the first regression layer, first, the PLS regression models are built between the process variables and quality variables of the historical processes with phase and mode indexes to get regression parameters. Then, using those regression parameters and the process variables of the target process, a series of assumed quality predictions can be obtained. Here, it is assumed that the target process belongs to each modeled process. This step connects all historical process information with the target process. Whether the information captured is related to the final quality of the target process will be judged in next regression layer. Consequently, in the second regression layer, based on these assumed quality predictions and real quality values of the target process, the PLS regression model is established, which captures the valuable information existing in the historical processes related to the final quality of the target process. Finally, by applying the two-dimensional, two-layer regression model, the final prediction quality can be obtained. problem at the same time, the strategy above will be applied to the intra-batch batch direction and handle the multi-phase problem by analyzing the intra-batch evolution. The framework of the two-dimensional, two-layer regression modeling strategy is shown in Figure 2. In the first regression layer, first, the PLS regression models are built between the process variables and quality variables of the historical processes with phase and mode indexes to get regression parameters. Then, using those regression parameters and the process variables of the target process, a series of assumed quality predictions can be obtained. Here, it is assumed that the target process belongs to each modeled process. This step connects all historical process information with the target process. Whether the information captured is related to the final quality of the target process will be judged in next regression layer. Consequently, in the second regression layer, based on these assumed quality predictions and real quality values of the target process, the PLS regression model is established, which captures the valuable information existing in the historical processes related to the final quality of the target process. Finally, by applying the two-dimensional, two-layer regression model, the final prediction quality can be obtained. Figure 2. Illustration of two-dimensional, two-layer regression modeling.
In the proposed method, the two-dimensional regression, which is mainly reflected by the regression parameters of the historical processes with both phase and mode indexes, traces the intra-batch and inter-batch characteristics, while the two-layer structure establishes the relationship between the target process and historical processes. Compared with the between-mode model, the two-dimensional, two-layer regression model handles the inter-batch evolution and intra-batch evolution simultaneously, and, at the same time, maintaining the unity of the method. However, the between-mode method deals with only the inter-batch multi-mode problem, which means other methods dealing with the intra-batch multi-phase problem or phase information are necessary. Compared with the traditional strategy [16,23], where historical modes are considered together and modeled by phase mean PLS method, and several phase mean models are built for different phases, respectively, the proposed method mixes all process mode and phase information together in the first regression layer, then using the second regression layer to ex- In the proposed method, the two-dimensional regression, which is mainly reflected by the regression parameters of the historical processes with both phase and mode indexes, traces the intra-batch and inter-batch characteristics, while the two-layer structure establishes the relationship between the target process and historical processes. Compared with the between-mode model, the two-dimensional, two-layer regression model handles the inter-batch evolution and intra-batch evolution simultaneously, and, at the same time, maintaining the unity of the method. However, the between-mode method deals with only the inter-batch multi-mode problem, which means other methods dealing with the intra-batch multi-phase problem or phase information are necessary. Compared with the traditional strategy [16,23], where historical modes are considered together and modeled by phase mean PLS method, and several phase mean models are built for different phases, respectively, the proposed method mixes all process mode and phase information together in the first regression layer, then using the second regression layer to extract valuable information for the target process, which is more reasonable to use regression method to judge the relationship between the processes, rather than excluding the information by phase division.
The details of the proposed two-dimensional, two-layer quality regression modeling strategy are introduced in the next section.

Two-Dimensional, Two-Layer Regression Model
In this section, after establishing the models representing the quality related information of each historical processes, the regression relationship between the process variables and the quality of the target process is analyzed through the regression model of the historical phases of historical modes, and, on this basis, online monitoring is carried out.
Firstly, the models representing the quality related information of each historical processes are established. Within phase c of mode m, X c,m,k (I c,m × J x ) and Y c,m are the normalized time-slice process variables and quality variables. The PLS algorithm is used to capture the relationship of the process data matrix X c,m,k and quality data matrix Y c,m , and the model is as follows: where T c,m,k and U c,m,k are the score matrices, P c,m,k and Q c,m,k are the loading matrices, and E c,m,k and F c,m,k are the residual matrices. The prediction model can be expressed as below: where B c,m,k is the regression parameter matrix. When only one final quality variable is considered, the prediction model can be written as:ŷ where β c,m,k is the regression parameter vector. So, the regression coefficient of historical phases β c,m is obtained: where K c stands for the number of time-slices within phase c of mode m, m stands for the number of the historical modes, m = 1, 2, . . . , M, and c stands for the historical phases in each batch, c = 1, 2, . . . , C.
For the target process, X t,k (I t × J x ) and y t are the normalized time-slice process variables and quality variable.
The regression model and the process variables of the target process are used to obtain the assumed quality prediction of the target process: where t stands for the target process, m stands for the number of the historical modes, m = 1, 2, . . . , M, and c stands for the number of historical phases, c = 1, 2, . . . , C. It should be noted thatŷ t,c,m,k is called the assumed quality prediction, which means it is obtained by assuming that the process variables belong to the certain phase of the certain mode with the regression parameter β c,m . By obtaining the assumed quality predictions, the quality information of historical phases and modes is shared by the target process. Further, the quality information of historical phases and modes will be judged and extracted by the next regression. Then, the relationship between the assumed quality predictions and the quality data of the target process will be established. All these assumed predictions of the historical modes can comprise a new matrix Z t,k (I t × MC), Z t,k = ŷ t,1,1,k , . . . ,ŷ t,m,1,k , . . . ,ŷ t,M,1,k , . . . ,  y t,1,c,k , . . . ,ŷ t,m,c,k , . . . ,ŷ t,M,c,k , . . . ,ŷ t,1,C,k , . . . ,ŷ t,m,C,k , . . . ,ŷ t,M,C,k . Then, the kth time-slice PLS regression model is built between Z t,k and y t as follows: where T t,k and u t,k are the score matrix and vector of the target process, P t,k and q t,k are the loading matrix and vector of the target process, and E t,k and f t,k are the residual matrix and vector of the target process. Then, novel predictions are obtained: where α t,k is the regression parameter of the kth time-slice model. The mean regression parameters of the target process can be obtained from the regression parameters of the time-slice models, where K t is the number of the time intervals within phase c. Then, the predictions based on the regression parameter of the whole phase, α t , are obtained: Corresponding coefficients can be obtained based on the assumed prediction quality obtained from the above analysis: After obtaining the z t,k vector at the k-th moment, the online T 2 statistics and online SPE statistics are calculated, and the calculation formula is as follows: where T 2 t,k is the T 2 statistic of the current k-th moment, and SPE t,k is the SPE statistic of the current k-th moment; z t,k is the residual vector at the k-th moment.
The corresponding control limits are calculated similarly as Equations (8) and (9).

Introduction of Injection Molding Process
Injection molding a typical batch process, which is one of the important technologies of plastic processing. To obtain plastic products, an injection molding process mainly consists of mold closing, injection, packing-holding, plasticizing, cooling, mold opening, part ejection, and other processes. Among those phases, injection, packing-holding, plasticizing, and cooling are the most important four operation phases to determine the quality of products: in the injection phase, the molten plastic is injected into the mold; then, in the packing-holding phase, a certain pressure is maintained to fulfill the hollows which may occur due to the temperature decrease; consequently, in the plasticizing phase, the raw material, i.e., plastic particles, is transported forward, plasticized, and melted; finally, in the cooling phase, the plastic is cooled in the mold until the product becomes sufficiently rigid for ejection. As the development of computers, all important process variables can be obtained online by high-precision sensors.
The material used in this experiment is high density polyethylene (HDPE). The quality analyzed in this experiment is the weight of injection molded parts. The operating conditions are shown in Table 1. The selected process variables are shown in Table 2, which are used to establish the model. The variable data shown in the table can be collected by sensors. The mode information of this experiment includes the packing pressure (PP) and the barrel temperature (BT). Five different modes are obtained. The experimental conditions are shown in Table 3. The data used in the modeling process are all real experiment data.

Normal Batch Monitoring
The two-dimensional, two-layer regression model proposed in this paper can track the process evolution that exists in both directions within and between batches and extract important quality information from historical processes. To illustrate the superiority of this method, real injection molding process data with both multi-mode and multi-phase characteristics are used for analyzing.
In this simulation, 18 batches in mode 1, mode 2, and mode 3, 22 batches in mode 4, and 17 batches in mode 5 are selected, respectively, as normal batches of historical modes. Mode 3, with 18 batches, is used as the target mode for modeling. Five batches of mode 3 are selected as the normal test batches. During regression modeling, the four-fold cross-validation method is used to determine the number of reserved latent variables of the traditional method and the proposed method, which is 2 after analyzation. The confidence level of α is set to 0.99.
To illustrate the advantages of the proposed method, it is compared with the traditional multi-mode and multi-phase strategy. In this kind of traditional strategy, historical modes are considered together and modeled by phase mean PLS method, and several models are built for different phases, respectively. To be fair, all concerned modes are available for modeling in both the proposed method and the traditional method.
The mean RMSE of the quality predictions of the test batches from mode 3 by the traditional method is 0.0077, while the mean RMSE obtained by the proposed method is 0.0071, so the proposed method provides a more accurate prediction. The RMSE values of online quality predictions of the test batches using the two methods are shown in Figure 3. The dot dashed line is the RMSE value of the traditional method, the full line is the RMSE value of the proposed method, and vertical dotted lines are the phase boundaries. It can be seen from Figure 3 that the RMSE values of the traditional method have obvious phase variation, which indicates that the precision of the traditional method may be affected by different phase models. Therefore, the proposed method provides better quality prediction results than the proposed method. The RMSE values of the traditional method reflect obvious different phase characteristics because, in the traditional method, four models have been established for the four phases, while, for the proposed method, the RMSE values are not only lower than those of the traditional method within each phase but also keep a continuous low value throughout the four phases. This is because the proposed method mixes all process mode and phase information together in the first regression layer and then uses the second regression layer to extract valuable information for the target process. The results prove that using the two-layer regression method to judge the relationship between the target process and the historical processes is better than the traditional method based on phase division.  Further, one of the five test batches from mode 3 is selected to show the online qualit prediction and monitoring results. The online quality prediction results of one test batc are shown in Figure 4. Similar conclusions can be obtained from Figure 4 with those from Figure 3, which proves the correctness of the conclusion from Figure 3.  Further, one of the five test batches from mode 3 is selected to show the online quality prediction and monitoring results. The online quality prediction results of one test batch are shown in Figure 4. Similar conclusions can be obtained from Figure 4 with those from Figure 3, which proves the correctness of the conclusion from Figure 3. Further, one of the five test batches from mode 3 is selected to show the online qualit prediction and monitoring results. The online quality prediction results of one test batc are shown in Figure 4. Similar conclusions can be obtained from Figure 4 with those from Figure 3, which proves the correctness of the conclusion from Figure 3. The monitoring results of the four phases of one normal test batch from mode 3 ar shown in Figures 5, 6, 7 and 8, respectively. This batch, which has some fluctuation, ma be identified as anormal batch by the traditional method, and false alarm will occur du to the process fluctuation. While, by the propose method, this normal batch can be nor mally monitored. The details are as follows: In Figure 5 and 6, it can be seen that both T and SPE of the proposed method and the traditional method do not exceed the contro limits in the injection phase and the packing-holding phase. In Figure 7, it can be seen tha in the plasticizing phase, T 2 and SPE of the proposed method do not exceed their respec tive control limits, while, for the traditional method, although T 2 does not exceed its con trol limit, there is an obvious period at the beginning of the phase during which the SPE value exceeds its control limit. In fact, the fluctuation of this batch is caused by the transi tion between the packing-holding phase and the plasticizing phase. During the transition  The monitoring results of the four phases of one normal test batch from mode 3 are shown in Figures 5-8, respectively. This batch, which has some fluctuation, may be identified as anormal batch by the traditional method, and false alarm will occur due to the process fluctuation. While, by the propose method, this normal batch can be normally monitored. The details are as follows: In Figures 5 and 6, it can be seen that both T 2 and SPE of the proposed method and the traditional method do not exceed the control limits in the injection phase and the packing-holding phase. In Figure 7, it can be seen that, in the plasticizing phase, T 2 and SPE of the proposed method do not exceed their respective control limits, while, for the traditional method, although T 2 does not exceed its control limit, there is an obvious period at the beginning of the phase during which the SPE value exceeds its control limit. In fact, the fluctuation of this batch is caused by the transition between the packing-holding phase and the plasticizing phase. During the transition, process characteristic change from the packing-holding phase characteristic to the plasticizing phase characteristic. This characteristic fluctuation is natural progress and does not affect the production, so this batch is normal. Because the traditional method which uses the phase mean model to represent the phase characteristic is more easily affected by the characteristic fluctuation away from the phase mean characteristic, this fluctuation is identified as a fault by the traditional method. However, this is not consistent with the conclusion from the real production. Because the proposed method comprehensively deals with the multi-phase and multi-mode problems by involving all quality-related phase and mode information in historical processes in the regression model, it is not easily affected by phase fluctuation, and it can offer the right monitoring result, in which this batch is normal during the plasticizing phase. In Figure 8, there is one point of SPE value of the proposed method exceeds the control limit. Generally, this point will not be considered as a fault. Therefore, it can be concluded from the monitoring results of the normal test batch that the proposed method is better than the traditional method because, by the two-dimensional, two-layer quality regression, all quality-related phase and mode information in historical processes have been extracted and utilized for quality prediction and monitoring.
proposed method exceeds the control limit. Generally, this point will not be considered as a fault. Therefore, it can be concluded from the monitoring results of the normal test batch that the proposed method is better than the traditional method because, by the two-dimensional, two-layer quality regression, all quality-related phase and mode information in historical processes have been extracted and utilized for quality prediction and monitoring.
(a) T 2 statistics ( b) SPE statistics a fault. Therefore, it can be concluded from the monitoring results of the normal test batch that the proposed method is better than the traditional method because, by the two-dimensional, two-layer quality regression, all quality-related phase and mode information in historical processes have been extracted and utilized for quality prediction and monitoring.

Abnormal Batch Monitoring
In the injection molding process, the sensor fault is one of the common faults. Because of the sensor fault, some important process variable cannot be detected, the control system cannot know the real progress, and wrong control commands may be delivered, leading to serious dangerousness of production.
To simulate the sensor fault, a test batch with the pressure variable removed during the packing-holding phase is monitored. The T 2 and SPE monitoring effects of the traditional method and the proposed method are shown in Figure 9, respectively. Compared with the traditional method, the amplitudes of the statistics of the proposed method are relatively larger. It can be concluded that the sensitivity of the proposed method is not affected, although it involves all quality-related phase and mode information in historical processes in the regression model, and it can identify the fault batch as fast as the traditional method, with an even larger alarm signal. So, the proposed method is better than the traditional method.

Abnormal Batch Monitoring
In the injection molding process, the sensor fault is one of the common faults. Because of the sensor fault, some important process variable cannot be detected, the control system cannot know the real progress, and wrong control commands may be delivered, leading to serious dangerousness of production.
To simulate the sensor fault, a test batch with the pressure variable removed during the packing-holding phase is monitored. The T 2 and SPE monitoring effects of the traditional method and the proposed method are shown in Figure 9, respectively. Compared with the traditional method, the amplitudes of the statistics of the proposed method are relatively larger. It can be concluded that the sensitivity of the proposed method is not affected, although it involves all quality-related phase and mode information in historical processes in the regression model, and it can identify the fault batch as fast as the traditional method, with an even larger alarm signal. So, the proposed method is better than the traditional method.

Conclusions
In this work, a two-dimensional, two-layer quality regression model is established to monitor multi-phase, multi-mode batch processes. First, traditional batch process monitoring method based on phase mean PLS is introduced. Then, aiming at the multi-mode problem and the multi-phase problem simultaneously, the framework of two-dimensional, two-layer quality regression modeling strategy is proposed by expanding the multi-mode analysis strategy into the multi-phase problem and dealing with the intrabatch evolution and inter-batch evolution at the same time. Accordingly, the two-dimen-

Conclusions
In this work, a two-dimensional, two-layer quality regression model is established to monitor multi-phase, multi-mode batch processes. First, traditional batch process monitoring method based on phase mean PLS is introduced. Then, aiming at the multi-mode problem and the multi-phase problem simultaneously, the framework of two-dimensional, two-layer quality regression modeling strategy is proposed by expanding the multi-mode analysis strategy into the multi-phase problem and dealing with the intra-batch evolution and inter-batch evolution at the same time. Accordingly, the two-dimensional, two-layer regression model is established, based on which the quality predictions are obtained, and the process is monitored. Through the application to an injection molding process, it is proved that the proposed methodology can provide better quality prediction and process monitoring results for multi-phase, multi-mode batch processes than the traditional method.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are openly available in https://share. weiyun.com/Jbc9pbWF.