Abstract
Simulated moving bed (SMB) chromatographic separation technology is a new adsorption separation technology with strong separation ability. Based on the principle of the adaptive neural fuzzy inference system (ANFIS), a soft sensing modeling method was proposed for realizing the prediction of the purity of the extract and raffinate components in the SMB chromatographic separation process. The input data space of the established soft sensor model is divided, and the premise parameters are determined by utilizing the meshing partition method, subtractive clustering algorithm, and fuzzy -means (FCM) clustering algorithm. The gradient, Kalman, Kaczmarz, and PseudoInv algorithms were used to optimize the conclusion parameters of ANFIS soft sensor models so as to predict the purity of the extract and raffinate components in the SMB chromatographic separation process. The simulation results indicate that the proposed ANFIS soft sensor models can effectively predict the key economic and technical indicators of the SMB chromatographic separation process.
1. Introduction
SMB chromatographic separation technology is a new separation technology developed on the basis of traditional fixed bed adsorption operation and true moving bed (TMB) chromatographic separation technology [1]. As the main modern adsorption separation technology, SMB chromatography has been utilized more and more in the complex mixture separation process, such as petrochemical, fine chemical, biopharmaceutical, and food processing. SMB chromatography technology is the cutting-edge technology in separation science, which preserves the high separation rate, low energy consumption, and low material consumption of the chromatogram. It also introduces the continuous, countercurrent, rectification, reflux, and other mechanisms of moving bed technology. Compared with the existing chemical separation technology (distillation, extraction, single-column chromatography), this technology can achieve automatic continuous separation, which can increase separation capacity and improve product yield, yield, and efficiency. It is the key technology for the chemical separation process and can also reuse the stationary phase and mobile phase to reduce cost and energy consumption [2]. In the cyclical operation of continuous production, the SMB system fully reflects its nonlinear, nonideal, and nonequilibrium characteristics. This multidegree-of-freedom system finds it difficult to make the important performance indicators, such as product purity, yield, and mobile phase consumption, in periodicity optimal states [3, 4].
The adaptive neural fuzzy inference system is an artificial intelligence inference technology with the advantages of fuzzy logic and neural network [5], which is a good nonlinear mapping technique. ANFIS can construct a clear input-output mapping structure, which provides a powerful learning and decision-making process for target systems including nonlinear systems and has been successfully applied in the prediction of short-term electricity prices [6], clay expansion capacity [7], wind speed profile estimation [8], and the thermal error compensation on the CNC machine [9]. The ANFIS fuzzy neural network model was used to predict the cutting speed in WEDM [10]. The ANFIS model was used to predict the density of ionic liquids at various temperatures [11]. The particle swarm optimization (PSO) algorithm was used to obtain appropriate subtractive clustering (SC) parameter settings, and an adaptive ANFIS model was built to predict the business success rate [12]. ANFIS and artificial neural network (ANN) models were established for predicting the drying characteristics of potato, garlic, and cantaloupe in convection hot air dryers [13]. The polynomial equation of the bank section by utilizing the hybrid ANFIS was used to stabilize the design and implementation of the river section [14]. A data-driven model of the debutanizer based on ANFIS was established for realizing the real-time composition monitoring [15]. ANFIS was applied to develop an accurate temperature-dependent intelligent model for correlating the loading capacity of amino acid salt solutions for a variety of amino acids [16]. The combination of the ANFIS and CFD method provides the nondiscrete domain in various dimensions and makes a smart tool to locally predict multiphase flow [17]. The ANFIS prediction model was established for the response during the stainless steel 202 turning operation [18].
In this paper, an ANFIS-based soft modeling method for the SMB chromatography separation process was proposed, which lays a foundation for the SMB system to achieve the periodic optimum stable working condition. The structure of the paper is described as follows. Section 2 introduces the SMB chromatographic separation technique and the structure of the established soft sensor model. Three input data space partitioning methods and four conclusion parameter optimization algorithms of ANFIS are introduced in Section 3. Section 4 describes the experimental simulation and result analysis. Finally, the conclusion is listed in the last part.
2. SMB Chromatography Separation Technology and Soft Sensor Modeling
2.1. SMB Chromatography Separation Technology
SMB chromatographic separation technology simulates the reverse flow of a stationary phase adsorbent by continuously switching the positions of each feed and discharge port. The basic principle of typical SMB is described as follows. The several chromatographic columns are connected to form a loop. The relative countercurrent of the stationary phase and the mobile phase is simulated by orderly moving the eluent inlet, the extract outlet, the feed liquid inlet, and the raffinate outlet in the direction of the mobile phase so as to achieve the separation of two components [19]. Taking the separation of two components as an example, the working principle is briefly explained. The working principle of the typical SMB chromatographic separation process is shown in Figure 1 [20, 21]. It is assumed that two components to be separated are A and B, respectively. The adsorption capacity of component A is stronger than that of component B. D is the desorbent, E is the extract, F is the feed, and R is the raffinate. The entire bed can be divided into four zones (respectively, referred to as I, II, III, and IV or 1, 2, 3, and 4) depending on the location and function of the liquid inlet and outlet. Each zone performs a specific function.
Zone I. The strong adsorption component A was resolved from the stationary phase adsorbent by using desorbent D. The fresh desorbent enters from the bottom of zone I and contacts with the adsorbent, and component A is rinsed out and then discharged from the extract outlet at the top of the zone to obtain an extract.
Zone II. This is the analytical zone of weakly adsorbed component B. In the liquid entering zone II from zone I (including A and D), since component A has a stronger adsorption capacity than component B, component B in zone II is continuously displaced by component A and enters zone III along with the flowing liquid and the fresh feed (A and B).
Zone III. Strongly adsorbed component A is adsorbed by the adsorbent in this zone. B is discharged from the raffinate outlet through a liquid flowing downward to obtain the raffinate.
Zone IV. Its function is to achieve regeneration of eluent D and reduce the amount of eluent D. Since the concentration of D in zone IV is much higher than the concentration of D flowing into zone IV from the zone III region (because the region is zone I in the last switching cycle, eluent D is adsorbed by the adsorbent after being flushed with A, so its concentration is higher). Therefore, the adsorption equilibrium will be reestablished and eluent D will be resolved to achieve regeneration. The regenerated eluent D and the fresh eluent D flow into zone I together. In this way, the reverse flow of the stationary phase adsorbent is simulated by continuously switching the liquid inlet and outlet position in the switching direction shown in Figure 1, thereby achieving an effect equivalent to TMB.
The feeding and discharging ports divide the SMB flow system into different zones. For the SMB with three zones, the functions of each zone are described as follows.
Zone I. Elution zone: it desorbs the adsorbed substances on the stationary phase and flows out from the E port to complete the regeneration process of the stationary phase.
Zone II. Distillation zone: continuous elution with the mobile phase makes the fast and slow components further separated in the chromatographic column.
Zone III. Adsorption zone: the raw materials from the inlet flow into this zone and are adsorbed and separated, so that fast components flow out of the outlet (R port) of this zone and slow components are adsorbed in the stationary phase.
2.2. Soft Sensor Modeling for the SMB Chromatography Separation Process
In the SMB chromatography adsorption separation process, the purity of the product components in the extraction solution and raffinate is more important. However, the online measurement of the component purity as a quality index is difficult to achieve, and there are many factors affecting the product component purity in the SMB separation process. Therefore, the soft sensor modeling method for predicting the component purity in the SMB chromatography separation process has theoretical significance and application value [22]. The basic idea of soft sensing technology is to combine the automatic control theory with the knowledge of the production process. According to some optimal criteria, because some important variables are difficult to measure or temporarily unable to measure, the computer technology is used to select other variables that are easy to measure and to infer or estimate by constructing a mathematical relationship. The soft sensor model can be expressed as where is the dominant variable estimated value, is the measurable disturbance, is the control input variable, is the measurable output variable, and is the off-line sampled value or the analytically calculated value of the estimated variable. The optimal estimate of can be obtained by constructing the relationship between these variables. Equation (1) reflects the relationship between the dominant variable and the input variable , the auxiliary variable , and the measurable disturbance , while the off-line sample value is often used for self-correction of the model. The structure of the soft sensor model for the SMB chromatographic separation process established in this paper is shown in Figure 2.
The auxiliary variables in the soft sensor model are selected from the measurement variables that can be provided in the control system based on the technique mechanism of the production process. Based on the technique mechanism and prior knowledge of the SMB chromatography separation process, the following variables are selected as the auxiliary variables of the soft sensor model. (1)Flow rate of the injection pump (F pump) of the raw material liquid inlet, whose unit is ml/min(2)Flow rate of the flushing pump (D pump) of the flushing fluid inlet, whose unit is ml/min(3)Valve switching time, whose unit is min
Select the following variables as the output variable of the established soft sensor model. (1)The purity of the target in the effluent of the E port. If there is impurity at the E port, this purity is smaller than 1(2)The purity of the impurity in the outlet effluent of the R port. If there is target material at the R port, this purity is smaller than 1(3)Divide the mass of the target in the E port by the target injection quality so as to obtain the yield of the target at the E port(4)Divide the mass of the impurity in the R port by the impurity injection quality so as to obtain the yield of the impurity at the R port
The auxiliary variables are used as inputs, and the purity of the target in the E port and the purity of the impurities in the R port are output variables. ANFIS is used to fit the nonlinear relationship among them so as to establish a predictive model of the corresponding economic and technical indicators. Based on a SMB chromatographic separation process, a soft sensor modeling method on the purity and yield of E and R ports is proposed by utilizing ANFIS. The historical data of the SMB chromatographic separation process were collected, and 1000 sets of historical data with uniformity and representative shown in Table 1 are selected as the modeling data. Before carrying out the soft sensor modeling with these data, it is necessary to perform data normalization shown in Equation (2) in order to eliminate the influence of variable unit on modeling accuracy. where is the initial value of the data set, is the data after being normalized, is the maximum in the data set, and is the minimum in the data set.
3. ANFIS and Training Algorithms
3.1. Basic Principles of ANFIS Based on Mesh Generation
In the early 1990s, Jang proposed an ANFIS based on the T-S model [5], which is a new type of neural network structure combining fuzzy theory. It uses a neural network to realize three basic processes of fuzzification part, fuzzy inference part, and defuzzification in fuzzy control so that the established system has the advantages of both theories and realizes the adaptive adjustment on the research object [23]. The typical structure of the ANFIS is shown in Figure 3. The prototype ANFIS adopts the meshing partition method to linearly divide the input space so as to determine the number of fuzzy rules [24].
Assume that the fuzzy inference system under consideration has two inputs and and a single output . For the first-order Sugeno fuzzy model, the general rule set with two if-then rules is represented as follows.
Rule 1. If is and is , then
Rule 2. If is and is , then
Nodes in the same layer have the same function. Here, the output of the th node in the layer is .
Layer 1. This layer node obscures the input signal. where or is the input node and is the membership of the fuzzy set (, or ). Here, the membership function of takes a bell-shaped function, which can be represented as where {, , } is called the antecedent parameters.
Layer 2. Multiply all input signals passed in the previous layer as the output of this layer.
Layer 3. The output in this layer is the normalized excitation intensity. The principle is that the th node calculates the ratio of the excitation strength of the th rule to the sum of the excitation strengths of all the rules.
Layer 4. All nodes in this layer are used as adaptive nodes with the node functions. where is the normalized excitation intensity from the third layer and {, , } is the conclusion parameters or consequent parameters.
Layer 5. The last layer sums all the signals to calculate the system output.
It can be seen from the structure of ANFIS that the system has two adaptive layers (layer 1 and layer 4). The first layer has three adjustable antecedent parameters associated with the input membership functions. Layer 4 has three adjustable consequent parameters associated with the first-order polynomial. The linear combination of the consequent parameters can be used to obtain the output of ANFIS, which can be expressed as where the elements of the column vector constitute a set of conclusion parameters {}. If there are groups of input and output data pairs and given the antecedent parameters, the order of the matrix is , , and , respectively. In general, the number of sample data is much larger than the number of unknown parameters, that is to say, . The least squares method (LSM) is used to obtain the best estimate of the conclusion parameter vector in the sense of minimizing , namely,
Then, based on the identification result of the conclusion parameters , the estimated output of ANFIS is calculated by
The root mean square error (RMSE) under the current antecedent and conclusion parameters can be calculated by Equation (14).
3.2. Antecedent Parameter Determination Algorithm of ANFIS
In essence, the prototype structure of ANFIS uses the mesh generation method to divide the input space by default. If the input is complex nonlinear data, it will inevitably lead to the exponential growth of the number of fuzzy rules, which will inevitably bring about dimensionality disasters. In addition, the linear division of the grid partition method cannot accurately reflect the spatial distribution of input data. Therefore, the subtraction clustering and fuzzy -means clustering algorithms are utilized to realize the segmentation of data space and the determination of the antecedent parameters.
3.2.1. Subtractive Clustering Algorithm
In a certain group of data, each data point is regarded as a candidate point of the cluster center, and an independent and fast single-time clustering algorithm that can calculate the number of clusters and the class centers is named as the subtractive clustering method [25]. According to the subtractive clustering method, there is a simple linear relationship between the number of data points and the calculation amount, which is not necessarily related to the dimension of the research object. Assuming that the data points in the -dimensional space have been normalized to a hypercube space, the density index at the data point is set as where is a positive integer representing a neighborhood radius of such data points. It can be seen that a data point with multiple data points in the vicinity has a high density value, and a point outside a certain field has a minimal influence on its density index.
In the subtractive clustering algorithm, the class center may be any point. Therefore, the density index of the point and its surrounding points in the neighborhood must be calculated with each point as the class center. The possibility of this point as the class center is analyzed. After the above operations, the first class center is the point with the highest density index selected from it, which is written as and whose density value is . The corrected density index of any point can be obtained by Eq. (16). where .
Obviously, the closer the data point is to , the smaller the density value and then, the less likely it is a class center. In the neighborhood with radius , the density value shows a decreasing trend. In general, in order to prevent the situation of the closer class centers, usually , for example, . After the density values of all the points are modified, the calculations of other class centers are deduced by the same method, and the density value correction is performed for all the points of the nonclass center. In this paper, the subtractive clustering method is used for pattern decomposition. Suppose the degree of the input and output belonging to a certain kind of center is and the input set contains input and output data, wherein each point is expressed as where the input of the th input and output data pair is represented as and the corresponding output is .
The steps of pattern extraction based on the subtractive clustering algorithm are described as follows. Step 1.There are vectors , and the initial value is definedStep 2.Let be the reference vector and be the comparison vector. Then, calculate the relationship between the two vectorswhere , , and indicates the Euclidean distance between the two vectors. The variance of the Gaussian function is calculated by the mean square error of the performance index Step 3.Update the relationship degree between the reference vector and the comparison vector bywhere is a small positive number () Step 4.Calculate by Eq. (20).Step 5.When all and are the same, the algorithm continues; otherwise, let and go to Step 2Step 6.According to the final obtained in the above steps, the number of clusters is equal to the number of convergence vectors. The original data with the same convergence vector is divided into one class with the convergence vector as the center
3.2.2. FCM Clustering Algorithm
The FCM clustering algorithm [26] divides the feature points in the feature space into classes , and the cluster center of the th class is represented by , where any feature point belongs to the membership degree of the th class, and satisfies the following conditions.
The objective function of the FCM clustering algorithm is described as where is the index weight that affects the degree of fuzzification of the membership matrix. The clustering problem is to find the membership degree matrix and the category center to make Equation (3) minimized. When utilizing the FCM clustering algorithm, a parameter that needs to be determined in advance is the number of classifications of the data set. The clustering validity function used in this paper is described as follows [27]. where
This clustering criterion combines the tightness and separation of fuzzy partitioning. By combining the clustering ISODATA algorithm with the proposed criterion function to obtain the optimal fuzzy partition, the algorithm procedure is described as follows. Step 1.Specify the maximum number of cluster centers (), the maximum number of iterations , the index weight , and the stop threshold Step 2.Initialize the fuzzy clustering center for the number of cluster centers Step 3.For the number of iterations , calculate the membership matrix and cluster center according to Eq. (24) and Eq. (25).where and . If , go to the next step; otherwise, repeat Step 3 Step 4.Calculate by Equation (23). If , go to Step 2; otherwise, stop the iteration and output the optimal cluster number , which satisfies .
3.3. Consequent Parameter Optimization Algorithms of ANFIS
3.3.1. Gradient Algorithm
The gradient descent method, also known as the fastest descent method, is a multidimensional unconstrained optimization problem calculation method based on gradients [28]. Considering an unconstrained optimization problem, where has a first-order continuous partial derivative. Starting from any point , the Taylor formula shows that the decreases fastest along the direction of the negative gradient. Carry out the first-order Taylor spread on at along to obtain where is the searching step size and is a unit vector, that is to say, .
Carry out the transposition on Equation (28) to get
It can be seen from Equation (29) that if , there are and . The term of the square brackets in Equation (29) is negative, and thus, , where is the descending direction of at . Based on the principle of the Cauchy-Schwarz theory, obtain
Equation (30) is the only condition for Equation (31) to take the equal sign. For the iteration, get where . So, the iteration formula of the steepest descent method can be represented as where is the number of iterations and is a constant. The expression of the th component in Equation (33) can be described as
3.3.2. Kalman Algorithm
Information fusion is the core of the Kalman algorithm [29]. The probability distribution with position as the variable obeys the Gaussian distribution, which is described as
The probability of is proportional to , and the result of the product also obeys a new Gaussian distribution due to
The item in parentheses is 0, so there is
Substitute the mean into to get
The above description can convert two Gaussian distributions into one Gaussian distribution. When there are multiple Gaussian distributions, they can be merged one by one. So the and can be replaced with the optimal estimate, that is to say,
There are also and . Arrange the equations to obtain a new covariance formula:
Iteratively update the scale factor by Equation (41).
Feed into the mean to get
Bring into the covariance formula to obtain
3.3.3. Kaczmarz Algorithm
The basic principle of the Kaczmarz algorithm is described as follows [30]. Let the map from to be defined as
The mapping from to can be defined as
Optionally select and calculate the iterative sequence by utilizing the recursive formula. where is defined by Equations (44) and (45).
The orthogonal projection from to the th plane can be defined as
So the mapping shown in Equation (47) is then expressed as
Let , where , and matrix satisfies the following relationship:
According to the above equations, the expression of can be deduced. where , and matrix depends on matrix . Based on Equations (49) and (50), the Kaczmarz iterative may be described as
Let , where is the orthogonal projection from to . Then, define
Based on Equation (52), matrix has the following property: where is a projection along to and , which can be replaced by .
3.3.4. PseudoInv Algorithm
The PseudoInv algorithm is implemented using the Moore-Penrose generalized inverse [31]. It is often applied to obtain the least norm least squares solution (least squares method) on the nonuniform linear equations and makes the form of the solution simple. The Moore-Penrose generalized inverse of the matrix is unique in both the real and complex domains and can be obtained by the singular value decomposition. The pseudoinverse of matrix can be defined as
The actual algorithm for calculating the pseudoinverse adopts Eq. (55).
When the number of columns of matrix is more than the number of rows, solving the linear equation using pseudoinverse is one of the many possible solutions. Specifically, is with the smallest Euclidean norm in all feasible solutions for this equation. When the number of rows in matrix is more than the number of columns, there may be no solution. In this case, obtained by pseudoreverse is a solution that minimizes the Euclidean distance of and .
4. Simulation Experiment and Result Analysis
Based on the technique of the SMB chromatographic separation process, a soft sensor model on the target purity of E and R ports is proposed by utilizing ANFIS. The ANFIS network has three input variables and one output variable. The historical data of the SMB chromatographic separation process was collected, and 1000 sets of historical data with uniformity and representative were selected. Then, the processed data is divided into two parts. The first 800 sets of data are used as training data, and the last 200 sets of data are used to verify the performance of the soft sensor models. According to the reference data, it can be found that the target purity in the E effluent and the yield of the impurity at the R port are correlated, and the impurity purity in the R effluent and the yield of the target at the E port are also correlated. Therefore, in the experimental simulation stage, the target purity in the E effluent and the impurity purity in the R effluent are selected as the output. In order to measure the performance of the predictive models, several performance indicators are defined below, where is the estimated value and is the actual value [32]. The model performance indicators are defined as shown in Table 2.
Based on the meshing partition method, subtractive clustering algorithm, and FCM clustering algorithm, the input data space of the SMB chromatographic separation process is divided and the premise parameters are determined. The Sugeno-type ANFIS soft sensor models based on three algorithms are established. Then, these models are optimized by the gradient, Kalman, Kaczmarz, and PseudoInv algorithms to obtain the optimized conclusion parameters so as to achieve the prediction of the target purity of the E port and R port in the SMB chromatographic separation process.
4.1. ANFIS Soft Sensor Model Based on Mesh Partition
When the input data space of the soft sensor model is divided by the meshing partition method, the membership function selects the Gaussian function, and the number of membership functions is 5. The membership function curves of the input data of the Sugeno-type ANFIS model based on meshing partition are shown in Figure 4. The number of training iterations is 500, and the established soft sensor models are tested by the utilized 200 sets of data with average and representative. The predicted simulation results are shown in Figures 5–8. Figure 5 shows the output contrast of the target purity in the E port effluent in the SMB chromatography process under the gradient, Kalman, Kaczmarz, and PseudoInv algorithms to obtain the optimized conclusion parameters of ANFIS. Figure 6 shows the predictive error comparison curves. Figure 7 shows the output contrast of the impurity purity in the R port effluent in the SMB chromatography process under the gradient, Kalman, Kaczmarz, and PseudoInv algorithms to obtain the optimized conclusion parameters of ANFIS. Figure 8 shows the predictive error comparison curves. Table 3 compares the predicted performance indicators of the established soft sensor models.
According to the above simulation results, it can be seen that the input data spatial division and premise parameter determination are realized by the meshing partition method and ANFIS soft sensor models based on the gradient, Kalman, and PseudoInv algorithms for optimizing conclusion parameters have the better prediction results on the key economic and technical indicators of the SMB chromatographic separation process. In order to distinguish the optimization performance of the four algorithms, based on the adopted four performance indicators (RMSE, SSE, MAPE, and MPE), the ANFIS soft sensor model based on the meshing partition method and the gradient algorithm has higher prediction accuracy than the other three algorithms.
4.2. ANFIS Soft Sensor Model Based on the Subtractive Clustering Algorithm
When the input data space of the soft sensor model is divided by the subtractive clustering algorithm, the membership function selects the Gaussian function. The membership function curves of the input data of the Sugeno-type ANFIS model based on the subtractive clustering algorithm are shown in Figures 9–11. The number of training iterations is 500, and the established soft sensor models are tested by the utilized 200 sets of data with average and representative. The predicted simulation results are shown in Figures 12–15. Figure 12 shows the output contrast of the target purity in the E port effluent in the SMB chromatography process under the gradient, Kalman, Kaczmarz, and PseudoInv algorithms to obtain the optimized conclusion parameters of ANFIS. Figure 13 shows the predictive error comparison curves. Figure 14 shows the output contrast of the impurity purity in the R port effluent in the SMB chromatography process under the gradient, Kalman, Kaczmarz, and PseudoInv algorithms to obtain the optimized conclusion parameters of ANFIS. Figure 15 shows the predictive error comparison curves. Table 4 compares the predicted performance indicators of the established soft sensor models.
According to the above simulation results, it can be seen that the input data spatial division and premise parameter determination are realized by the subtractive clustering algorithm and ANFIS soft sensor models based on the gradient, Kalman, and PseudoInv algorithms for optimizing conclusion parameters have the better prediction results on the key economic and technical indicators of the SMB chromatographic separation process. In order to distinguish the optimization performance of the four algorithms, based on the adopted four performance indicators (RMSE, SSE, MAPE, and MPE), the ANFIS soft sensor model based on the meshing partition method and the PseudoInv algorithm has higher prediction accuracy than the other three algorithms.
4.3. ANFIS Soft Sensor Model Based on the FCM Clustering Algorithm
When the input data space of the soft sensor model is divided by the FCM clustering algorithm, the membership function selects the Gaussian function. The membership function curves of the input data of the Sugeno-type ANFIS model based on the FCM clustering algorithm are shown in Figures 16–18. The number of training iterations is 500, and the established soft sensor models are tested by the utilized 200 sets of data with average and representative. The predicted simulation results are shown in Figures 19–22. Figure 19 shows the output contrast of the target purity in the E port effluent in the SMB chromatography process under the gradient, Kalman, Kaczmarz, and PseudoInv algorithms to obtain the optimized conclusion parameters of ANFIS. Figure 20 shows the predictive error comparison curves. Figure 21 shows the output contrast of the impurity purity in the R port effluent in the SMB chromatography process under the gradient, Kalman, Kaczmarz, and PseudoInv algorithms to obtain the optimized conclusion parameters of ANFIS. Figure 22 shows the predictive error comparison curves. Table 5 compares the predicted performance indicators of the established soft sensor models.
According to the above simulation results, it can be seen that the input data spatial division and premise parameter determination are realized by utilizing the FCM clustering algorithm and ANFIS soft sensor models based on the Kalman, Kaczmarz, and PseudoInv algorithms for optimizing conclusion parameters have the better prediction results on the key economic and technical indicators of the SMB chromatographic separation process. In order to distinguish the optimization performance of the four algorithms, based on the adopted four performance indicators (RMSE, SSE, MAPE, and MPE), the ANFIS soft sensor model based on the FCM clustering algorithm and the PseudoInv algorithm has higher prediction accuracy than the other three algorithms.
5. Conclusions
In this paper, a soft sensing modeling method of the SMB chromatographic separation process based on ANFIS is proposed. Three input data space division and antecedent parameter determination methods are combined with four consequent parameter optimization algorithms to realize the ANFIS soft sensing models on the target purity of the E port and R port in the SMB chromatographic separation process. The simulation results verify that the proposed models can obtain better prediction results on the purity of components, which will supply the foundation of the quality closed-loop control.
Data Availability
There are no data available for this paper.
Conflicts of Interest
The authors declare no conflict of interest.
Authors’ Contributions
Dan Wang participated in the algorithm simulation and draft writing. Jie-Sheng Wang participated in the concept, design, and critical revision of this paper. Shao-Yan Wang participated in the interpretation and commented on the manuscript. Shou-Jiang Li participated in the SMB chromatographic separation process technique. Zhen Yan participated in the data collection and analysis. Wei-Zhen Sun participated in the simulation methods.
Acknowledgments
This work is partially supported by the project of National Natural Science Foundation of China (Grant No. 21576127), the Basic Scientific Research Project of Institution of Higher Learning of Liaoning Province (Grant No. 2017FWDF10), and the project of Liaoning Provincial Natural Science Foundation of China (Grant No. 20180550700).