Next Article in Journal
Low-Frequency Seismic Noise Properties in the Japanese Islands
Next Article in Special Issue
Sensor and Actuator Fault Diagnosis for Robot Joint Based on Deep CNN
Previous Article in Journal
Empirical Mode Decomposition-Derived Entropy Features Are Beneficial to Distinguish Elderly People with a Falling History on a Force Plate Signal
Previous Article in Special Issue
Constrained Parameter Estimation for a Mechanistic Kinetic Model of Cobalt–Hydrogen Electrochemical Competition during a Cobalt Removal Process
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimation of Feeding Composition of Industrial Process Based on Data Reconciliation

School of Automation, Central South University, Changsha 410083, China
*
Author to whom correspondence should be addressed.
Entropy 2021, 23(4), 473; https://doi.org/10.3390/e23040473
Submission received: 1 March 2021 / Revised: 11 April 2021 / Accepted: 13 April 2021 / Published: 16 April 2021
(This article belongs to the Special Issue Complex Dynamic System Modelling, Identification and Control)

Abstract

:
For an industrial process, the estimation of feeding composition is important for analyzing production status and making control decisions. However, random errors or even gross ones inevitably contaminate the actual measurements. Feeding composition is conventionally obtained via discrete and low-rate artificial testing. To address these problems, a feeding composition estimation approach based on data reconciliation procedure is developed. To improve the variable accuracy, a novel robust M-estimator is first proposed. Then, an iterative robust hierarchical data reconciliation and estimation strategy is applied to estimate the feeding composition. The feasibility and effectiveness of the estimation approach are verified on a fluidized bed roaster. The proposed M-estimator showed better overall performance.

1. Introduction

Complete, accurate and reliable data measurements are important for process modeling, model identification, real-time online optimization and process control. However, the measured values of variables in the actual measurement process are inevitably contaminated by random errors or even gross errors. Therefore, the process measurements will deviate from the real values. In addition, the feeding composition, which is important in making control decisions, cannot be measured online due to the high economic cost and technical limitations. Therefore, it is necessary to process the measurements to guarantee the reliability and observability of the system. Based on this, the estimation of feeding composition can be better performed. Data reconciliation is an effective approach for this aim [1]. The principle of data reconciliation is to minimize the sum of square deviations between the coordinated value of variables and their measurements under the condition of satisfying material balances, heat balances and boundary constraints of process variables. From the mathematical point of view, data reconciliation is a constrained optimization problem [2], which has been widely used in many important procedures such as in-line process monitoring [3,4], enhanced control strategy performance [5,6], parameter estimation [7], soft sensor applications [8], quality control [9] and industrial optimization [10,11] among others [3,4,5,6,7,8,9,10,11].
Traditional data reconciliation typically assumes that measurements contain only random errors with an average value of zero and normal distributions. Nonetheless, the measurements will be contaminated by gross errors. If gross errors are not processed, the accuracy of data reconciliation results will be reduced. The gross error detection is therefore particularly important for data reconciliation. There are currently three main methods for handling measurements with gross errors. The first strategy is sequential gross error detection and data reconciliation. A variety of statistical tests, including the global test [12], the measurement test [13], the nodal test [14], generalized likelihood ratio test [15] and principal component analysis test [16], are used for gross error detection. Following the detection and elimination of gross errors, the traditional data reconciliation is carried out. The second strategy is simultaneous gross error detection and data reconciliation. Yu et al. [17] proposed a support vector regression method for recursive simultaneous data reconciliation and gross error detection. Zhang et al. [18] used a novel particle filter (PF) algorithm based on the measurement test (MT) to solve the dynamic simultaneous gross error detection and data reconciliation. Yuan et al. [19] established a new hierarchical Bayesian framework, which can simultaneously estimate the real value of process variables and obtain the magnitudes of gross errors. The last strategy is robust data reconciliation. The estimations close to the true value can be obtained by reducing the influence of gross errors on data reconciliation results. Tjoa and Biegler [20] constructed an M-estimator based on mixed distribution, which takes both random errors and gross errors into account and minimized the M-estimator as an objective function. Johnson and Kramer [21] presented a new robust method for data reconciliation based on a probability bootstrapping theory. Arora and Biegler [22] evaluated the Fair and Hampel estimators in stationary and dynamical constrained problems. Hampel’s three-part estimator presented the best performance. Wang and Romagnoli [23] designed a partially adaptive estimator based on a generalized T distribution and a fully adaptive estimator based on density estimation by studying the robustness of a set of existing estimators. Özyurt and Pike [24] analyzed the performance of a number of robust estimators such as Hampel, Cauchy and Fair. Ragot et al. [25] adopted a cost function that has been improved based on the contaminated probability density distribution model for the design of robust estimator. Prata et al. [26] applied the proposed strategy based on the particle swarm optimization (PSO) and the robust Welsch estimator to real-time online detection of Bulk polymerization of propylene. Zhang et al. [27] constructed a robust estimator by using the quasi-weighted least squares. Through the comparative analysis of various methods like Welsch, quasi-weighted least squares and comentropy M-estimators, the feasibility and effectiveness of robust estimators in simultaneous gross error detection and data reconciliation were demonstrated [28]. Alighardashi et al. [29] proposed a maximum likelihood framework for simultaneous data reconciliation and gross error detection for steady-state data. Xie et al. [30] utilized a novel robust estimator to improve the robustness of data reconciliation. Forty-eight different robust M-estimators can be founded in a recent review study, including Fair, Cauchy, Biweight, Jin, Welsch estimator and Xie [31].
The key to estimating feeding composition based on data reconciliation is the robust estimator and the design of the estimation strategy. An excellent robust estimator is embodied in two aspects: On the one hand, the objective function is bounded, and the bound value is small; on the other hand, the influence function can converge quickly and converge to zero. The objective functions of Cauchy [32] and Fair [33] estimators increase with the increase in standard measurement residuals, indicating that the error will greatly affect the reconciled results. Besides, the convergence speed of their influence functions is very slow, and the bound values are large, which shows that the two estimators cannot suppress the influence of gross errors on reconciled results. The objective function of Welsch [26] estimator is limited, but the bound value is also large, and the influence function does not converge rapidly. Although the influence function of Xie [30] estimator may converge to zero and the convergence speed is improved in comparison to the estimators above, it still cannot achieve the fast standard. Based on the robust estimation theory, a novel robust estimator is proposed to enhance the performance of the robust estimators. The objective function is bounded, the bound value is smaller, and the influence function converges toward zero faster. At the same time, aiming at the problem that the change of feeding composition caused by the long test period of industrial processes is unknown, an iterative robust hierarchical data reconciliation and estimation strategy based on the heat balance is proposed. In order to further verify the robustness and effectiveness of the proposed robust estimator and estimation of feeding composition strategy, a series of comparative experiments are carried out through two numerical examples and the real data from a fluidized bed roaster for zinc smelting.
After initial remarks, motivation aspects and applications of the data reconciliation issue, the paper is organized as follows. Section 1 presents preliminary approaches to traditional data reconciliation and robust estimator. The novel robust estimator is proposed, and two numerical examples are used to demonstrate the effectiveness of the proposed estimator in Section 2. Section 3 provides an iterative robust hierarchical data reconciliation and estimation strategy to estimate feeding composition. The feeding composition of fluidized bed roaster is then estimated in Section 4. Finally, Section 5 provides the conclusions.

2. Preliminaries

2.1. Data Reconciliation

The purpose of data reconciliation is to obtain the reconciled data by processing the raw measurements. The reconciled data not only satisfy the constraints of the process model but are also closer to the true value. There are three main assumptions for the system using data reconciliation: The system is in a steady state, the measurement errors follow a normal distribution with zero mean, and each measured variable is independent of each other. Taking into account the presence of measurement errors, the measurement model may be expressed as follows:
x = x * + ε
where x denotes the vector of raw process measurements, x * denotes the vector of true values of the process variables and ε is the random measurement errors.
Based on the principle of data reconciliation, the general steady-state data reconciliation problem can be stated as a form of solving the weighted least squares solution satisfying the process model and boundary constraints:
min   ( x x ˜ ) T 1 ( x x ˜ ) s . t .       F ( x ˜ , u ) = 0                   G ( x ˜ , u ) 0
where x ˜ is the vector of reconciled data, is the diagonal covariance matrix of measurement errors, u is the vector of unmeasured variables; F represents the process model, which is used as equality constraints in the optimization problem, and G denotes the inequality constraints indicating variable boundaries.

2.2. M-Estimator

In the classical weighted least square of data reconciliation, Equation (2), it is usually assumed that the process measurements contain only random errors. However, in the actual measurement process, the measured values may contain gross errors. The presence of gross errors has a serious effect on the conventional data reconciliation, resulting in the propagation of gross errors. This will cause the reconciled data to not satisfy the process model and deviate from the true value of variables seriously. As a result, the robust estimator is used to account for gross errors in measurements.
Currently, there are numerous approaches for robust data reconciliation, the majority of which are based on the theory of M-estimator. The function of the measurement residuals is constructed as the objective function of the robust estimator. By reducing the weight of the measurements with gross errors, it can prevent gross errors smearing other measurements. Hence, the robust estimator can well reconcile the data with gross errors. The problem of robust data reconciliation can be expressed as follows:
min   i = 1 n ρ ( ξ i ) = min   i = 1 n ρ ( x i x ˜ i σ i ) s . t .       F ( x ˜ , u ) = 0                   x ˜ i min x ˜ i x ˜ i max ,   i = 1 , 2 , , n                   u l min u l u l max ,   l = 1 , 2 , , N n
where ρ is the robust estimator, ξ i = ( x i x ˜ i ) / σ i is the standardized residual for the i th measured variable, x i ,   x ˜ i and u l are, respectively, the measured data, the reconciled data for the i th measured variable and the estimated value for the j th unmeasured variable, x ˜ i min ,   x ˜ i max ,   u l min and u l max are, respectively, the lower limit, the upper limit for the i th reconciled variable, the lower limit, the upper limit for the l th unmeasured variable. N ,   n ,   N n denote the total number of variables, the total number of measured variables and the total number of variables not measured, respectively. σ i represents the standard deviation for variable measurement errors.
Most M-estimators are not based on a clearly defined probability distribution known in advance. Most of them are based on a simple mathematical structure. The mathematical structure of M-estimator, i.e., ρ , has some general characteristics as follows:
ρ is is continuous;
ρ ( ξ ) = ρ ( ξ ) ;
ρ ( ξ ) 0 and ρ is integrable
ρ ( ξ 1 ) ρ ( ξ 2 ) , for | ξ i | < | ξ j | ;
ρ ( 0 ) = 0 .
As a crucial index in assessing the robustness of the M-estimator, the influence function (IF) [34] is directly proportional to the derivative of the objective function of the estimator. The influence function value refers to the effect of different deviation on the estimator. In general, the function may be defined as the first derivative of the objective function to standardized residuals, which can be expressed as follows:
IF ( ξ ) = d ρ ( ξ ) d ξ
Some general characteristics of IF [35] are:
IF is limited;
IF is continuous or piecewise continuous;
IF ( ξ ) = IF ( ξ ) ;
IF is nearly linear near the origin ( IF ( ξ ) k ξ , k 0 , for small ξ ), but this characteristic is not necessary;
the rejection point of IF (the point where IF is zero) is finite to suppress large deviations.
On the basis of a robust M-estimation theory, the influence function needs to be continuous and bounded. When ξ is small, the influence function is proportional to ξ ; when ξ is infinite, the influence function can converge to a constant, which indicates that the robust estimator can suppress the effect of gross errors in the reconciled results. The influence function of the weighted least squares estimator is the standardized residual ξ , which indicates that the influence function will increase with the measurement errors. It can be seen that the weighted least squares estimator is not robust because its influence function is boundless, and the reconciled results can be easily affected by gross errors.

3. Data Reconciliation Based on a Novel Robust Estimator

3.1. A Novel Robust Estimator

The most robust estimators are constructed according to the following requirements: The objective function of the estimator is bounded, and the bound value is small; the influence function can rapidly converge to a constant. From the observation of the mathematical structure of 48 different robust estimators in the literature [31], it has been shown that the objective functions of these M-estimators comprise limitless function, limited function, non-piecewise function, piecewise function, multi parameters and a single parameter. The performance of the M-estimator whose objective function is limited, non-piecewise and has single parameter is analyzed. It is found that the M-estimator with a relatively rapid convergence rate basically contains the mathematical expression with e as the base number. The derivative of exponential function based on e is equal to itself. Moreover, e ξ decreases with increasing ξ and tends towards 0. When ξ tends towards 0, e ξ tends towards 1. Therefore, the objective function can be limited if the structure of 1 e ξ is included in. For instance, the Welsch estimator and the objective function is presented in Equation (5). In order to improve the robustness of the M-estimator, Xie divided 1 e ξ with 1 + e ξ . By the decline of 1 + e ξ , the overall rise speed of the objective function is improved. Furthermore, the converging speed of the influence function is accelerated. The objective function of Xie estimator is illustrated under Equation (6).
Welsch:
ρ ( ξ i ) = c w 2 2 ( 1 exp ( ( ξ i c w ) 2 ) )
Xie:
ρ ( ξ i ) = 1 exp ( ( ξ i c x ) 2 ) 1 + exp ( ( ξ i c x ) 2 )
where c w and c x are, respectively, tuning parameters of Welsch and Xie estimators.
As discussed above, in order to limit the effect of gross errors on the reconciled results, a novel robust estimator is created based on the general characteristics of the M-estimator and influence function. The objective function of the novel robust estimator is defined as follows:
ρ ( ξ i ) = c p 2 [ 1 exp ( ( ξ i / c p ) 2 ) ] 4 [ 1 + exp ( ( ξ i / c p ) 4 ) ]
where c p is the tuning parameter; the definition of ξ i is the same as above. By increasing the index of ξ , the downward velocity of 1 + e ξ is increased. Then, the objective function increases faster, and the influence function converts to 0 faster. At the same time, in order to satisfy the characteristic that the influence function is almost linear near the origin, the objective function is multiplied by c p 2 / 4 . The influence function of the novel robust estimator function expressed in Equation (7) is as follows:
IF ( ξ i ) = d ρ ( ξ i ) d ξ i = ξ i exp ( ( ξ i c p ) 2 ) [ 1 + exp ( ( ξ i c p ) 4 ) + 2 ξ i 2 c p 2 ( 1 exp ( ( ξ i c p ) 2 ) ) ] 2 [ 1 + exp ( ( ξ i c p ) 4 ) ] 2
It can be seen from Equations (7) and (8) that when the standardized residual is small, the objective function of the novel robust estimator tends to be gradually zero, and the influence function is proportional to the standardized residual. It shows that the weighted least squares method can be used to reconcile the data where measurements contain only random errors. When the standardized residual is significant, the objective function of the novel robust estimator tends to be constant little by little, and the influence function quickly converges to zero. It is demonstrated that the novel robust estimator can adjust the weight of the estimator based on the standardized residual size when the measurements contain gross errors. The larger the gross errors, the lower the weights, and the greater the effects of restraining the gross errors. Therefore, the proposed estimator is more robust.

3.2. The Tuning Parameter of the Novel Robust Estimator

The use of robust estimators in data regression typically involves the inverse relationship between relative efficiency and robustness [36]. Relative efficiency refers to the fit quality of the estimated value reconciled by M-estimator, when errors follow another distribution relative to the reference distribution (usually assumed to be Normal distribution). The so-called robustness refers to the performance of the M-estimator under various nonnormal error distributions. The more robust an M-estimator is, the less its relative efficiency is. The relationship between robustness and relative efficiency is addressed by the tuning parameters. When comparing the performance of different M-estimators, the level of relative efficiency and the reference distribution should be consistent. In general, the tuning parameters for each M-estimator are obtained at the relative efficiency level of 95%. Then the robustness of distinct M-estimators is judged.
The mathematical definition of relative efficiency is illustrated in Equation (9):
E f f [ IF ( ξ ) , f ( ξ ) ] = V f [ IF f ( ξ ) , f ( ξ ) ] V [ IF ( ξ ) , f ( ξ ) ]
where V f is the asymptotic variance of the reference estimator, V is the asymptotic variance of the specified estimator, f ( ξ ) is the error reference probability distribution, IF f ( ξ ) is the influence function of the reference estimator and IF ( ξ ) is the influence function of the specified estimator. The asymptotic variance V is defined as follows [37]:
V [ IF ( ξ ) , f ( ξ ) ] = + IF 2 ( ξ ) f ( ξ ) d ξ [ + IF ( ξ ) f ( ξ ) d ξ ] 2
where IF ( ξ ) is the derivative of the influence function IF ( ξ ) . Since IF ( ξ ) may be discontinuous, Equation (10) can be further expressed by Equation (11):
V [ IF ( ξ ) , f ( ξ ) ] = + IF 2 ( ξ ) f ( ξ ) d ξ [ + IF ( ξ ) f ( ξ ) d ξ ] 2 = + IF 2 ( ξ ) f ( ξ ) d ξ [ + IF ( ξ ) f ( ξ ) d ξ ] 2
where f ( ξ ) is the derivative of f ( ξ ) .
The tuning parameter of the proposed estimator at the relative efficiency level of 95% in respect to the Normal distribution is calculated as follows:
(1) The Normal distribution is selected as the reference distribution, and its probabilistic distribution and the first derivative are respectively expressed as follows:
f ( ξ ) = 1 σ 2 π exp ( ξ 2 2 )
f ( ξ ) = ξ 1 σ 2 π exp ( ξ 2 2 ) = ξ f ( ξ )
(2) The least square estimator is selected as the reference estimator, and its influence function is expressed as follows:
IF f ( ξ ) = d ρ ( ξ ) d ξ = d d ξ ( ξ 2 2 ) = ξ
(3) The influence function of the proposed estimator is shown in Equation (8). The expressions of V f and V can be obtained by using Equation (11):
V f = 2 0 + IF f 2 ( ξ ) f ( ξ ) d ξ ( 2 0 + IF f ( ξ ) ξ f ( ξ ) d ξ ) 2
V = 2 0 + IF 2 ( ξ ) f ( ξ ) d ξ ( 2 0 + IF ( ξ ) ξ f ( ξ ) d ξ ) 2
(4) The expression of relative efficiency can be obtained from Equation (9). Making the expression equal to 0.95, a univariate equation on the tuning parameter c p can be constructed, as shown in Equation (17). By numerical calculation, the tuning parameter at the relative efficiency level of 95% in respect to the Normal distribution is 1.5424.
V f V = 2 0 + ξ 2 1 σ 2 π exp ( ( ξ 2 / 2 ) ) d ξ [ 2 0 + ξ 2 1 σ 2 π exp ( ( ξ 2 / 2 ) ) d ξ ] 2 2 0 + ξ 2 exp ( 2 ( ξ / c p ) 2 ) [ 1 + exp ( ( ξ / c p ) 4 ) + 2 ξ 2 ( 1 exp ( ( ξ / c p ) 2 ) ) / c p 2 ] 2 4 [ 1 + exp ( ( ξ / c p ) 4 ] 4 1 σ 2 π exp ( ( ξ 2 / 2 ) ) d ξ [ 2 0 + ξ 2 exp ( ( ξ / c p ) 2 ) [ 1 + exp ( ( ξ / c p ) 4 ) + 2 ξ 2 ( 1 exp ( ( ξ / c p ) 2 ) ) / c p 2 ] 2 [ 1 + exp ( ( ξ / c p ) 4 ] 2 1 σ 2 π exp ( ( ξ 2 / 2 ) ) d ξ ] 2 = 0.95
To more accurately reflect the competitive relationship between the relative efficiency and robustness of the M-estimator, the comparative images of the objective function and influence function of the proposed, Xie and Welsch estimators are shown in Figure 1. Relative efficiency levels are 90%, 95%, 98% and 99%, respectively. The tuning parameters of the three estimators at different relative efficiency levels are presented in Table 1. Table 1 and Figure 1 show that the higher the relative efficiency level, the higher the values of the tuning parameters for the three estimators. The higher the tuning parameter, the slower the convergence speed of the influence function to 0 and the weaker the robustness. Hence, when comparing the robustness of distinct estimators, the comparative analysis has to be performed at the same level of relative efficiency.
In order to further analyze the effectiveness of the proposed estimator, the objective function and influence function of the proposed estimator are compared with the other four robust estimators. The four types of robust estimators are Fair, Cauchy, Welsch and Xie. The objective functions of Fair and Cauchy are shown in Equations (18) and (19):
Fair:
ρ ( ξ i ) = c f 2 ( | ξ i | c f ln ( 1 + | ξ i | c f ) )
Cauchy:
ρ ( ξ i ) = c c 2 ln ( 1 + ξ i 2 c c 2 )
where c f and c c are, respectively, tuning parameters of Fair and Cauchy estimators. The tuning parameters of the five M-estimators under the level of 95% relative efficiency are c f = 1.3998 , c w = 2.9846 , c c = 2.3849 , c x = 1.9597 and c p = 1.5424 . The results of the comparison are illustrated in Figure 2 and Figure 3.
As can be seen from Figure 2 and Figure 3, the objective function and influence function of the proposed M-estimator satisfy their general characteristics. When the standardized residual is low, the objective and influence functions of several robust estimators are all close, which indicates that low measurement errors have little effect on these estimators. When the standardized residual is large, the objective function of the Fair and Cauchy estimator increases at a higher rate. Their influence functions do not converge rapidly and eventually converge to a nonzero value. It shows that both Fair and Cauchy estimators are sensitive to gross errors. Although the objective function of Welsch does not diverge, it tends to be a greater constant than Xie and the novel robust estimator. The objective functions of the novel robust estimator and Xie estimator are less than those of the other three estimators in the case of gross errors in measurements. Furthermore, the proposed M-estimator tends to be a constant lower than Xie. It can be seen further from Figure 2 that when there are gross errors in measurements, the influence function of the proposed M-estimator decreases much faster than that of Xie. When the standardized residual exceeds 4, the influence function of the proposed M-estimator converges to 0. Nevertheless, the influence function of Xie converges to 0 only when the standardized residual is above 5. Hence one can see that the convergence speed of influence function of the novel robust estimator is faster than that of Xie, and the effect of suppressing gross errors is excellent. The proposed method is therefore more robust and efficient than other estimators.
Next, two numerical examples are used to further verify the effectiveness of the proposed robust estimator, and the proposed robust estimator is compared with several estimators. The data reconciliation procedure can be interpreted mathematically as an optimization problem. In this paper, the state transition algorithm (STA) is applied as the optimization method for data reconciliation problems [38]. The local search operator, the global search operator and the heuristic search operator are adopted as the state transformation operators. This can prevent the solution process from falling into the local optimum and reduce the search time for the global optimum.

3.3. Linear Case

In this part, the measurement network [26] as shown in Figure 4 is adopted. The network is a linear structure consisting of four nodes and seven streams. As can be seen from Figure 4, the seven streams meet the following material balance
x 1 x 2 + x 4 + x 6 = 0 x 2 x 3 = 0 x 3 x 4 x 5 = 0 x 5 x 6 x 7 = 0
where x = [ x 1 ,   x 2 ,     ,   x 7 ] T is set of stream variables. In this linear example, all variables are assumed to be measured, and the true value of each variable is X = [5, 15, 15, 5, 10, 5, 5]. The true value of each variable is added with some random errors to get the corresponding measurements. Let the standard deviation of random errors be 2.5% of the true value of each stream. The diagonal matrix of the variance is as follows:
= 0.01 × d i a g ( 1.562 , 14.062 , 14.062 , 1.562 , 6.250 , 1.562 , 1.562 )
For a precise analysis of the performance of each robust estimator, SSE (sum of squares due to error), TER (total error reduction) and RER (relative error reduction) [27] are used as indicators to assess the precision of the reconciled results. The definitions included are:
SSE = i = 1 n ( x ˜ i x i * ) 2 + j = 1 m ( u ˜ j u j * ) 2
TER = i = 1 n ( x i x i * σ i ) 2 i = 1 n ( x ˜ i x i * σ i ) 2 i = 1 n ( x i x i * σ i ) 2
RER = i = 1 n ( | x i * x i | x i * | x i * x ˜ i | x i * ) i = 1 n | x i * x i | x i *
where x ˜ i represents the reconciled data for the i th measurable variable; x i * represents the true value for the i th measurable variable; x i represents the measured data for the i th measurable variable; u ˜ j is the estimated value of the j th unmeasurable variable; u j * is the true value of the j th unmeasurable variable; σ i denotes the standard deviation for the i th measurable variable; n is the total number of measurable variables; m is the total number of unmeasurable variables. The smaller the SSE, the closer the reconciled results are to the true value. The larger the TER, the more accurate the reconciled results. The larger RER shows that the reconciled results are less affected by gross errors.

3.3.1. There Are Two Gross Errors in the Measurement Variables

Streams 2 and 5 are selected to add gross errors in the magnitude of 2 and 1, respectively, and other streams contain only random errors. The results of the data reconciliation and indicators for several methods are presented in Table 2. For convenience, "Proposed" is used to represent the proposed estimator in this table and subsequent tables.
Table 2 shows that the SSE of the proposed robust estimator (Proposed) is much lower than that of Fair with 1.2117 and Cauchy with 0.1287 and is less than that of Xie and Welsch estimator. The TER of the novel robust estimator is much larger than that of Fair with 0.2953 and Cauchy with 0.7606, and higher than that of Xie and Welsch estimator. In addition, the RER of Proposed is also higher than that of the other four robust estimators. The results show that in the case of two gross errors, the reconciled results of the novel robust estimator are closer to the true value, and the proposed estimator performs better than other methods.

3.3.2. There Are Three Gross Errors in the Measurement Variables

Streams 2, 5 and 7 are selected to add gross errors with the magnitudes of 2, 1 and −1.5, respectively, and other streams contain only random errors. The results of the data reconciliation and indicators for several methods are presented in Table 3.
Table 3 shows that the SSE of the proposed robust estimator (Proposed) is much lower than that of Fair with 1.2866 and Cauchy with 0.1071 and is less than that of Xie and Welsch estimator. The TER of the novel robust estimator is much larger than that of Fair with 0.3474, and higher than that of Xie, Welsch and Cauchy estimator. In addition, the RER of Proposed is also higher than that of the other four robust estimators. The results show that in the case of three gross errors, the reconciled results of the novel robust estimator are closer to the true value, and the proposed estimator performs better than other methods.
In order to better illustrate the effectiveness of the proposed estimator, the linear numerical example is divided into three groups for comparative experiments. Each group comprises 100 data samples. In the first group, there is gross error in stream 2 ( x 2 ). In group two, there are gross errors in stream 2 ( x 2 ) and stream 5 ( x 5 ). The third group contains gross errors in stream 2 ( x 2 ), stream 5 ( x 5 ) and stream 7 ( x 7 ). In each group of experiments, the selected variables are added with gross errors ranging from 3 times to 9 times the standard deviation. The measurement test (MT) is used to detect gross errors. Two criteria are selected to assess the performance of different methods, including the observed power (OP) and the Average Number of Type I errors (ATVI) [2]. OP and ATVI are defined as follows:
O P = The   number   of   gross   errors   detected   correctly The   total   number   of   gross   errors
A T V I = The   number   of   gross   errors   detected   mistakenly The   total   number   of   simulations
OP indicates the number of gross errors correctly detected. ATVI indicates the number of gross errors mistakenly detected. The larger value of the OP is, the stronger ability of the estimator to identify gross errors. The lower value of the ATVI is, the less frequently gross errors are detected incorrectly. The statistical graph of the process in which gross errors are detected by distinct methods in each group is shown in Figure 5. The abscissa of each graph corresponds to seven stream variables. The y-axis is the number of times each variable is detected to contain gross errors. The OP and ATVI of the proposed estimator, Xie, Welsch, Cauchy and Fair are obtained by statistical analysis of the process diagram, as illustrated in Figure 6.
As depicted in Figure 6, when one variable, two variables and three variables contain gross errors, the OP obtained by the proposed estimator is higher than that obtained by other estimators. The result shows that the proposed estimator has a higher probability of identifying gross errors correctly. In addition, the AVTI obtained by the proposed estimator is lower than that obtained by other methods. It can also be seen from Figure 5 that the proposed estimator has a lower number of gross errors mistakenly detected than other estimators. It is shown that gross errors have less influence on the data reconciliation based on the Proposed estimator.

3.4. Nonlinear Case

A nonlinear numerical example [22] is used in this part. Within the nonlinear system, there are five measurable variables, three unmeasurable variables and six nonlinear constraint equations, which are described as Equation (27):
{ 0.5 x 1 2 0.7 x 2 + x 3 u 1 + x 2 2 u 1 u 2 + 2 x 3 u 3 2 255.8 = 0 x 1 2 x 2 + 3 x 1 x 3 2 x 2 u 1 x 2 u 2 u 3 + 111.2 = 0 x 3 u 1 x 1 + 3 x 2 + x 1 u 2 x 3 u 3 33.57 = 0 x 4 x 1 x 3 2 + u 2 + 3 u 3 = 0 x 5 2 x 3 u 2 u 3 = 0 2 x 1 + x 2 x 3 u 1 + u 2 u 3 126.6 = 0
where x i ( i = 1 , 2 , , 5 ) is the measured variable and u j ( j = 1 , 2 , 3 ) is the unmeasured variable. Measurements of the measured variables are obtained from the true value added, some random errors and gross errors. The true value, the measured value, the standard deviation of the measured variables and the true value of the unmeasured variables are given in Table 4. Amongst these,   x 2   ,   x 3 and x 5 contain gross errors. The results of the data reconciliation and indicators for several methods are presented in Table 4.
Table 4 shows that the SSE calculated by the proposed robust estimator is significantly less than that calculated by Fair, Cauchy and Welsch estimator, and it is less than that calculated by Xie estimator. Furthermore, the TER and RER of the proposed estimator are larger than those of the other four estimators. The results show tha for the nonlinear case, the reconciled results of the new robust estimator are also near the true value. Consequently, the proposed robust estimator can effectively suppress the effect of gross errors on the reconciled results and has greater robustness and excellent effectiveness.

4. Feeding Composition Estimation Based on Iterative Data Reconciliation

In order to address the problem caused by the unknown change in feeding composition, an iterative robust hierarchical data reconciliation and estimation strategy is proposed. The method is applied to a fluidized bed roaster for zinc smelting. Firstly, the industrial process of fluidized bed roaster is introduced. At the same time, difficulties in reconciling the fluidized bed roaster are discussed in Section 3.1. Then, on the basis of the establishment of the mechanism balance model of the fluidized bed roaster, specific solutions are presented in accordance with the existing difficulties.

4.1. Industrial Process Description

The fluidized bed roaster is a type of thermal equipment that applies fluidization technology to make materials desulphurized by oxidation roasting in the metallurgical industry. The roasting process of the fluidized bed roaster is a gas-solid reaction procedure. By sending a great deal of air into the oven from the bottom, a very intense exothermic reaction occurs in the material layer of the sulfur ore under the agitation of the air. Oxygen combines with sulfide to form sulfur dioxide and precious metals are converted into oxides or sulphates. The structure diagram for the fluidized bed roaster is shown in Figure 7.
In the fluidized roasting process, accurate and reliable measurements are the basis of process control, functional analysis and production management. However, in the real process, due to inaccurate sensors, equipment leaks and measurement bias, real-time measurements are inevitably affected by random errors and gross errors. That has a serious effect on the accuracy of measurements. In addition, some important variables, such as smoke volume and calcine quality, cannot be measured due to limitations in measurement technology and the environment.
The ZnS content of the feed during the roasting of the zinc concentrate will affect the indices of operation, such as the resistance to roasting, the temperature of the boiling layer and the content in the gaseous effluents. Then, the technical-economic index of the level of zinc soluble in calcine is affected, which field staff can judge from according to roasting. The compositions of mixed zinc concentrate and calcine are sampled only twice a day, once in the morning and once at night. Then the percentage of Zn, S and other elements in the mixed zinc concentrate and that of soluble zinc and soluble sulfur in the calcine are detected by the laboratory. Nevertheless, the sampling frequency of the variables measured in real time, like feed volume and blast volume, is 5 seconds, which is much greater than that of the laboratory data. Hence, due to the long test cycle of the mixed zinc concentrate composition in the fluidized bed roaster, the laboratory data cannot reflect the change in the real-time feeding composition.
Concerning the problems of measurements contaminated by gross errors, irredundant process data and a long test cycle of the feeding composition, an iterative robust hierarchical data reconciliation and estimation strategy for feeding composition is proposed.

4.2. Iterative Robust Hierarchical Data Reconciliation and Composition Estimation Framework

Figure 8 illustrates the framework of iterative robust hierarchical data reconciliation and estimation strategy for feeding composition. The major steps are as follows:
Step 1: Establish the steady-state mechanism model of the fluidized bed roaster The establishment of a reasonable process model is the basis for data reconciliation. However, the actual industrial process is complex and needs to be simplified accordingly. The following assumptions are put forward to establish a credible and comprehensible mechanism model:
(1) The only exothermal component of the mixed zinc concentrate is ZnS.
(2) Since the specific heat of calcine and soot is close and the mass ratio is close to 1:1, soot can be classified as calcine in the calculation.
Therefore, the material and heat balance at steady-state are expressed by Equations (28) and (29):
Material balance:
16 97 α ( M 1 + M 2 ) 1000 + V 0 ( ρ 0 ρ 1 ) = 0
Heat balance:
1000 ( M 1 + M 2 ) T 1 c p 1 + V 0 T 0 c p 2 + 443508 97 α ( M 1 + M 2 ) 1000 = 1000 ( 1 16 97 α ) ( M 1 + M 2 ) T 3 c p 3 + V 0 T 2 c p 0 + Q
where α denotes the percentage of ZnS in zinc concentrate; M 1 and M 2 represent the first and second feeding quantity, respectively; V 0 is the blast velocity; T 0 , T 1 , T 2 and T 3 denote the temperature variables for blast, zinc concentrate, gas and calcine, respectively; c p 0 , c p 1 , c p 2 and c p 3 are specific heat for gas, zinc concentrate, blast and calcine, respectively; Q describe the heat loss; ρ 0 and ρ 1 are, respectively, the density of air and gas. Among these variables, M 1 , M 2 , V 0 , T 0 , T 2 and T 3 are measured variables. ρ 1 , Q and α are unmeasured variables. The other parameters are all fixed values.
Step 2: Preliminary estimation of α based on measurements and heat balance. Estimation of α can be considered as an optimization problem. Combined with the heat balance and measurements, the ZnS estimation model before data reconciliation could be expressed as follows:
min Y 1 = ( 1000 ( M 1 + M 2 ) T 1 c p 1 + V 0 T 0 c p 2 + 443 , 508 97 α ( M 1 + M 2 ) 1000 1000 ( 1 16 97 α ) ( M 1 + M 2 ) T 3 c p 3 V 0 T 2 c p 0 Q ) 2
s . t .           α min α α max
where Y 1 is the objective of the optimization problem; α min and α max are the top and bottom bounds of α , respectively. In this case, Q can be roughly considered as a fixed value, which can be obtained by expert knowledge. The preliminary estimation of α will be used to construct the hierarchical model in Step 3.
Step 3: Robust hierarchical data reconciliation. First, the robust reconciliation model of material balance layer is established to obtain the reconciled values of the first and second feeding quantity and blast velocity. Then, the reconciled values of the material balance layer are used as precise values to construct the robust reconciliation model of heat balance layer. Finally, the reconciled results of the temperature for the blast, gas and calcine can be obtained. Therefore, the hierarchical data reconciliation model based on the novel robust estimation is expressed as follows:
Material balance layer:
min f 1 = i = 1 n ρ ( x i , x ^ i ) = i = 1 n ( 1 + ( ( x i x ^ i ) 2 / σ i 2 c p 14 1 ) exp ( ( ( x i x ^ i ) / σ i c p ) 2 ) ) = 1 + ( ( M 1 M ^ 1 ) 2 / σ M 1 2 c p 14 1 ) exp ( ( ( M 1 M ^ 1 ) / σ M 1 c p ) 2 ) + 1 + ( ( M 2 M ^ 2 ) 2 / σ M 2 2 c p 14 1 ) exp ( ( ( M 2 M ^ 2 ) / σ M 2 c p ) 2 ) + 1 + ( ( V 0 V ^ 0 ) 2 / σ V 0 2 c p 14 1 ) exp ( ( ( V 0 V ^ 0 ) / σ V 0 c p ) 2 ) + 1 + ( ρ 1 ρ ˜ 1 2 / c p 14 1 ) exp ( ( ρ 1 ρ ˜ 1 / c p ) 2 )
s . t .   { 16 97 α ( M ^ 1 + M ^ 2 ) + V ^ 0 ( ρ 0 ρ 1 ) = 0 M ^ 1 min M ^ 1 M ^ 1 max M ^ 2 min M ^ 2 M ^ 2 max V ^ 0 min V ^ 0 V ^ 0 max ρ 1 min ρ 1 ρ 1 max
where f 1 is the objective function of the robust estimator; x i and x i ^ represent the measurement and reconciled data for the i th measured variable respectively; M ^ 1 min , M ^ 2 min , V ^ 0 min and ρ 1 min denote the bottom bounds of the reconciled results for the first and second feeding quantity, blast velocity and the density of the air, respectively. M ^ 1 max , M ^ 2 max , V ^ 0 max and ρ 1 max denote the top bounds of the reconciled results for the first and second feeding quantity, blast velocity and the density of the air, respectively.
Heat balance layer:
min f 2 = i = 1 n ρ ( x i , x ^ i ) = i = 1 n ( 1 + ( ( x i x ^ i ) 2 / σ i 2 c p 14 1 ) exp ( ( ( x i x ^ i ) / σ i c p ) 2 ) ) = 1 + ( ( T 0 T ^ 0 ) 2 / σ T 0 2 c p 14 1 ) exp ( ( ( T 0 T ^ 0 ) / σ T 0 c p ) 2 ) + 1 + ( ( T 2 T ^ 2 ) 2 / σ T 2 2 c p 14 1 ) exp ( ( ( T 2 T ^ 2 ) / σ T 2 c p ) 2 ) + 1 + ( ( T 3 T ^ 3 ) 2 / σ T 3 2 c p 14 1 ) exp ( ( ( T 3 T ^ 3 ) / σ T 3 c p ) 2 ) + 1 + ( Q Q ˜ 2 / c p 14 1 ) exp ( ( Q Q ˜ / c p ) 2 )
s . t . { 1000 ( M ^ 1 + M ^ 2 ) T 1 c p 1 + V ^ 0 T ^ 0 c p 2 + 443508 97 α ( M ^ 1 + M ^ 2 ) 1000 = 1000 ( 1 16 97 α ) ( M ^ 1 + M ^ 2 ) T ^ 3 c p 3 + V ^ 0 T ^ 2 c p 0 + Q T ^ 0 min T ^ 0 T ^ 0 max T ^ 2 min T ^ 2 T ^ 2 max T ^ 3 min T ^ 3 T ^ 3 max Q min Q Q max
where f 2 is the objective function of the robust estimator; T ^ 0 min , T ^ 2 min , T ^ 3 min and Q min denote the bottom bounds of the reconciled results for the temperature of the blast, gas, calcine and heat loss, respectively; T ^ 0 max , T ^ 2 max , T ^ 3 max and Q max denote the top bounds of the reconciled results for the temperature of the blast, gas, calcine and heat loss, respectively.
Step 4: Reconciled estimation of α ^ based on reconciled results and heat balance. The reconciled results obtained in Step 3 are substituted for the objective function of the optimization problem based on heat balance. Then the estimated value of the percentage for ZnS after reconciliation could be obtained. The ZnS estimation model after data reconciliation can be represented as follows:
min Y 2 = ( 1000 ( M ^ 1 + M ^ 2 ) T 1 c p 1 + V ^ 0 T ^ 0 c p 2 + 443 , 508 97 α ^ ( M ^ 1 + M ^ 2 ) 1000 1000 ( 1 16 97 α ^ ) ( M ^ 1 + M ^ 2 ) T ^ 3 c p 3 V ^ 0 T ^ 2 c p 0 Q ) 2
s . t . α ^ min α ^ α ^ max
where Y 2 is the objective of the optimization problem; α ^ represents the estimated value of percentage for ZnS after reconciliation; α ^ min and α ^ max denote the bottom and top bounds of α ^ respectively. Through the iterative process of the above steps, greater accuracy and a smaller range of ZnS fluctuation in the feed can be obtained. Therefore, staff can better understand changes of the feeding composition and make adjustments as required.

5. Results from Real Industrial Data

To further demonstrate the effectiveness of the proposed method, 100 steady-state samples are collected from the real industry of a fluidized bed roaster. In addition, the data reconciliation results of the proposed estimator are compared with those of Xie estimator. The results of the comparison between the two methods are presented in Figure 9, Figure 10, Figure 11, Figure 12, Figure 13 and Figure 14, which respectively present the reconciled results for the first feeding quantity, the second feeding quantity, blast velocity, blast temperature, gas temperature and calcine temperature. From Figure 9, Figure 10, Figure 11, Figure 12, Figure 13 and Figure 14, it is known that the reconciled results by the novel robust estimator have a smaller range of fluctuation than that obtained by Xie estimator. It shows that whatever the uncertainty of the fluctuation of the measurements and the number of gross errors they contain, the data reconciliation approach based on the novel robust estimator is highly robust.
In addition, the standard deviations of the reconciled results and measurements for six variables obtained from the two methods are calculated, and the results of the comparison are presented in Figure 15. Figure 15 shows that, for the first feeding quantity, the seconding feed quantity, the blast velocity and the temperature of blast, the standard deviation of the reconciled results based on the novel robust estimator is less than Xie. For the temperature of gas and calcine, the standard deviations of reconciled results for the two methods are closer, because the measurement fluctuation ranges for both two variables are all very small. The results demonstrate that the data reconciliation based on the novel robust estimator provides better reconciliation results for variables with a large fluctuation range.
For the change in feeding composition, the above 100 steady-state samples are also collected to compare the percentage of ZnS before and after data reconciliation. Experimental results are shown in Figure 16 and Figure 17. It is evident that the ZnS fluctuation range obtained by the proposed method is smaller than that before data reconciliation. It illustrates that the percentage of ZnS estimated by the reconciled results is more accurate and could better reflect the change of feeding composition over a period of time.

6. Conclusions

The feeding composition in an industrial process can better reflect the state of production. In this paper, estimation of feeding composition based on robust data reconciliation is proposed to account for the unknown change in feeding composition. According to the robust estimation theory, a novel robust estimator is presented to handle the measurements with random errors and gross errors. In comparison with the objective function and influence function of other robust estimators, the proposed estimator is proved to be more robust. Besides, the linear and nonlinear numerical examples are used to further verify the effectiveness of the estimator. Then, an iterative robust hierarchical data reconciliation and estimation of feeding composition strategy is proposed and used for a fluidized bed roaster. The reconciled and estimated results show that the proposed strategy is beneficial for the estimation of feeding composition in the actual industrial process.

Author Contributions

Conceptualization, Y.L. and B.S.; methodology, Y.L.; software, Y.L.; validation, Y.L.; formal analysis, Y.L.; investigation, Y.L., Z.F. and M.J.; resources, Y.L.; data curation, Y.L. and Z.F.; writing—original draft preparation, Y.L.; writing—review and editing, Y.L. and B.S.; visualization, Y.L.; supervision, Y.L.; project administration, Y.L.; funding acquisition, B.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (61973321) and the Natural Science Foundation of Hunan Province (Grant No. 2019JJ50823).

Data Availability Statement

Publicly available datasets were analyzed in this study. It can be found here: https://open.library.ubc.ca/cIRcle/collections/ubctheses/831/items/1.0078522 (accessed on 1 February 2021).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kuehn, D.R.; Davidson, H. Computer control II. Mathematics of control. Chem. Eng. Prog. 1961, 57, 44–47. [Google Scholar]
  2. Narasimhan, S.; Jordache, C. Data Reconciliation and Gross Error Detection: An Intelligent Use of Process Data; Elsevier: Amsterdam, The Netherlands, 1999. [Google Scholar]
  3. Prata, D.M.; Lima, E.L.; Pinto, J.C. In-Line Monitoring of Bulk Polypropylene Reactors Based on Data Reconciliation Procedures. Macromol. Symp. 2008, 271, 26–37. [Google Scholar] [CrossRef]
  4. Prata, D.M.; Lima, E.L.; Pinto, J.C. Nonlinear Dynamic Data Reconciliation in Real Time in Actual Processes. Comput. Aided Chem. Eng. 2009, 27, 47–54. [Google Scholar]
  5. Abu-El-Zeet, Z.H.; Roberts, P.D.; Becerra, V.M. Enhancing model predictive control using dynamic data reconciliation. Alche J. 2002, 48, 324–333. [Google Scholar] [CrossRef]
  6. Sun, B.; Yang, C.; Wang, Y.; Gui, W.; Craig, I.; Olivier, L. A comprehensive hybrid first principles/machine learning modeling framework for complex industrial processes. J. Process Control 2020, 86, 30–43. [Google Scholar] [CrossRef]
  7. Faber, R.; Arellano-Garcia, H.; Li, P.; Wozny, G. An optimization framework for parameter estimation of large-scale systems. Chem. Eng. Process. 2007, 46, 1085–1095. [Google Scholar] [CrossRef]
  8. Schladt, M.; Hu, B. Soft Sensor Based on Nonlinear Steady-State Data Reconciliation in the Process Industry. Chem. Eng. Process. 2007, 46, 1107–1115. [Google Scholar] [CrossRef]
  9. Su, Q.; Bommireddy, Y.; Shah, Y.; Ganesh, S.; Moreno, M.; Liu, J.; Gonzalez, M.; Yasdanpanah, N.; O’Connor, T.; Reklaitis, G.V.; et al. Data reconciliation in the Quality-byDesign (QbD) implementation of pharmaceutical continuous tablet manufacturing. Int. J. Pharm. 2019, 563, 259–272. [Google Scholar] [CrossRef]
  10. Chiari, M.; Bussani, G.; Grotolli, M.G.; Pierucci, S. On-line Data Reconciliation and Optimisation: Refinery Applications. Comput. Chem. Eng. 1997, 21, s1185–s1190. [Google Scholar] [CrossRef]
  11. Lee, M.H.; Lee, S.J.; Han, C.; Chang, K.S.; Kim, S.H.; Sang, H.Y. Hierarchical on-line data reconciliation and optimization for an industrial utility plant. Comput. Chem. Eng. 1998, 22, s247–s254. [Google Scholar] [CrossRef]
  12. Almásy, G.A.; Sztanó, T. Checking and correction of measurements on the basis of linear system model. Probl. Control Inf. Theory 1975, 4, 57–69. [Google Scholar]
  13. Mah, R.S.H.; Tamhane, A.C. Detection of gross errors in process data. Aiche J. 2010, 28, 828–830. [Google Scholar] [CrossRef]
  14. Mah, R.S.; Stanley, G.M.; Downing, D.M. Reconciliation and rectification of process flow and inventory data. Ind. Eng. Chem. Process Des. Dev. 1976, 15, 175–183. [Google Scholar] [CrossRef]
  15. Narasimhan, S.; Mah, R.S.H. Generalized likelihood ratio method for gross error identification. Aiche J. 1987, 33, 1514–1521. [Google Scholar] [CrossRef]
  16. Tong, H.; Crowe, C.M. Detection of gross erros in data reconciliation by principal component analysis. Aiche J. 2010, 41, 1712–1722. [Google Scholar] [CrossRef]
  17. Yu, M.; Hong-Ye, S.U.; Jian, C.H.U. A support vector regression approach for recursive simultaneous data reconciliation and gross error detection in nonlinear dynamical systems. Acta Autom. Sin. 2009, 35, 707–716. [Google Scholar]
  18. Zhang, Z.; Chen, J. Simultaneous data reconciliation and gross error detection for dynamic systems using particle filter and measurement test. Comput. Chem. Eng. 2014, 69, 66–74. [Google Scholar] [CrossRef]
  19. Yuan, Y.; Khatibisepehr, S.; Huang, B.; Li, Z. Bayesian method for simultaneous gross error detection and data reconciliation. Aiche J. 2015, 61, 3232–3248. [Google Scholar] [CrossRef]
  20. Tjoa, I.B.; Biegler, L.T. Simultaneous strategies for data reconciliation and gross error detection of nonlinear systems. Comput. Chem. Eng. 1991, 15, 679–690. [Google Scholar] [CrossRef]
  21. Johnson, L.P.M.; Kramer, M.A. Maximum likelihood data rectification: Steady-state systems. Alche J. 1995, 41, 2415–2426. [Google Scholar] [CrossRef]
  22. Arora, N.; Biegler, L.T. Redescending estimators for data reconciliation and parameter estimation. Comput. Chem. Eng. 2001, 25, 1585–1599. [Google Scholar] [CrossRef]
  23. Wang, D.; Romagnoli, J.A. A Framework for Robust Data Reconciliation Based on a Generalized Objective Function. Ind. Eng. Chem. Res. 2003, 42, 3075–3084. [Google Scholar] [CrossRef]
  24. Özyurt, D.B.; Pike, R.W. Theory and practice of simultaneous data reconciliation and gross error detection for chemical processes. Comput. Chem. Eng. 2004, 28, 381–402. [Google Scholar] [CrossRef]
  25. Ragot, J.; Chadli, M.; Maquin, D. Mass balance equilibration: A robust approach using contaminated distribution. Aiche J. 2005, 51, 1569–1575. [Google Scholar] [CrossRef]
  26. Prata, D.M.; Schwaab, M.; Lima, E.L.; Pinto, J.C. Simultaneous robust data reconciliation and gross error detection through particle swarm optimization for an industrial polypropylene reactor. Chem. Eng. Sci. 2010, 65, 4943–4954. [Google Scholar] [CrossRef]
  27. Zhang, Z.; Shao, Z.; Chen, X.; Wang, K.; Qian, J. Quasi-weighted least squares estimator for data reconciliation. Comput. Chem. Eng. 2010, 34, 154–162. [Google Scholar] [CrossRef]
  28. Llanos, C.E.; Sanchéz, M.C.; Maronna, R.A. Robust estimators for data reconciliation. Ind. Eng. Chem. Res. 2015, 54, 5096–5105. [Google Scholar] [CrossRef]
  29. Alighardashi, H.; Magbool Jan, N.; Huang, B. Expectation Maximization Approach for Simultaneous Gross Error Detection and Data Reconciliation Using Gaussian Mixture Distribution. Ind. Eng. Chem. Res. 2017, 56, 14530–14544. [Google Scholar] [CrossRef]
  30. Xie, S.; Yang, C.; Yuan, X.; Wang, X.; Xie, Y. A novel robust data reconciliation method for industrial processes. Control Eng. Pract. 2019, 83, 203–212. [Google Scholar] [CrossRef]
  31. de Menezes, D.Q.F.; Prata, D.M.; Secchi, A.R.; Pinto, J.C. A review on robust M-estimators for regression analysis. Comput. Chem. Eng. 2021, 147, 107254. [Google Scholar] [CrossRef]
  32. Rey, W.J.J. Introduction to Robust and Quasi-Robust Statistical Methods; Springer Science & Business Media: Berlin, Germany, 2012. [Google Scholar]
  33. Fair, R.C. On the robust estimation of econometric models. Ann. Econ. Soc. Meas. 1974, 3, 667. [Google Scholar]
  34. Jin, S.Y.; Li, X.W.; Huang, Z.J.; Meng, L. A new target function for robust data reconciliation. Ind. Eng. Chem. Res. 2012, 51, 10220–10224. [Google Scholar] [CrossRef]
  35. Hoaglin, D.C.; Mosteller, F.; Tukey, J.W. Understanding Robust and Exploratory Data Analysis; Wiley: New York, NY, USA, 1983; Volume 3. [Google Scholar]
  36. Albuquerque, J.S.; Biegler, L.T. Data reconciliation and gross-error detection for dynamic systems. Aiche J. 1996, 42, 2841–2856. [Google Scholar] [CrossRef]
  37. Shevlyakov, G.; Morgenthaler, S.; Shurygin, A. Redescending m-estimators. J. Stat. Plann. Inference 2008, 138, 2906–2917. [Google Scholar] [CrossRef] [Green Version]
  38. Zhou, X.J.; Yang, C.H.; Gui, W.H. State transition algorithm. J. Ind. Manag. Optim. 2012, 8, 1039–1056. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Characteristicfunctions of three estimators at distinct relative efficiency levels.
Figure 1. Characteristicfunctions of three estimators at distinct relative efficiency levels.
Entropy 23 00473 g001aEntropy 23 00473 g001b
Figure 2. Objective functions for Xie, Welsch, Cauchy, Fair and Proposed estimator.
Figure 2. Objective functions for Xie, Welsch, Cauchy, Fair and Proposed estimator.
Entropy 23 00473 g002
Figure 3. Influence functions for Xie, Welsch, Cauchy, Fair and Proposed estimator.
Figure 3. Influence functions for Xie, Welsch, Cauchy, Fair and Proposed estimator.
Entropy 23 00473 g003
Figure 4. Diagram of the measurement network.
Figure 4. Diagram of the measurement network.
Entropy 23 00473 g004
Figure 5. Statistical diagram of gross errors detection process for distinct methods.
Figure 5. Statistical diagram of gross errors detection process for distinct methods.
Entropy 23 00473 g005aEntropy 23 00473 g005b
Figure 6. Results of OP and AVTI with different methods.
Figure 6. Results of OP and AVTI with different methods.
Entropy 23 00473 g006
Figure 7. Structure diagram of fluidized bed roaster.
Figure 7. Structure diagram of fluidized bed roaster.
Entropy 23 00473 g007
Figure 8. The framework of iterative robust hierarchical data reconciliation and composition estimation.
Figure 8. The framework of iterative robust hierarchical data reconciliation and composition estimation.
Entropy 23 00473 g008
Figure 9. Reconciled results for the first feed material quantity.
Figure 9. Reconciled results for the first feed material quantity.
Entropy 23 00473 g009
Figure 10. Reconciled results for the second feed material quantity.
Figure 10. Reconciled results for the second feed material quantity.
Entropy 23 00473 g010
Figure 11. Reconciled results for the blast velocity.
Figure 11. Reconciled results for the blast velocity.
Entropy 23 00473 g011
Figure 12. Reconciled results for the temperature of blast.
Figure 12. Reconciled results for the temperature of blast.
Entropy 23 00473 g012
Figure 13. Reconciled results for the temperature of gas.
Figure 13. Reconciled results for the temperature of gas.
Entropy 23 00473 g013
Figure 14. Reconciled results for the temperature of calcine.
Figure 14. Reconciled results for the temperature of calcine.
Entropy 23 00473 g014
Figure 15. The standard deviation of reconciled results and measurements.
Figure 15. The standard deviation of reconciled results and measurements.
Entropy 23 00473 g015
Figure 16. The percentage of ZnS before data reconciliation.
Figure 16. The percentage of ZnS before data reconciliation.
Entropy 23 00473 g016
Figure 17. The percentage of ZnS after data reconciliation.
Figure 17. The percentage of ZnS after data reconciliation.
Entropy 23 00473 g017
Table 1. Tuning parameters of three estimators at distinct relative efficiency levels.
Table 1. Tuning parameters of three estimators at distinct relative efficiency levels.
M-EstimatorTuning Parameter
1Welsch { c w = 2.3828           E f f = 90 % c w = 2.9846           E f f = 95 % c w = 3.9077           E f f = 98 % c w = 4.7343           E f f = 99 %
2Xie { c x = 1.6705           E f f = 90 % c x = 1.9597           E f f = 95 % c x = 2.3409           E f f = 98 % c x = 2.6359           E f f = 99 %
3Proposed { c p = 1.3082           E f f = 90 % c p = 1.5424           E f f = 95 % c p = 2.0942           E f f = 98 % c p = 2.6060           E f f = 99 %
Table 2. Reconciled results of different methods with two gross errors.
Table 2. Reconciled results of different methods with two gross errors.
StreamTrueOriginal Meas. Meas. with Gross Error ProposedXieWelschCauchyFair
x 1 54.9954.9955.01115.01525.03465.05195.1777
x 2 1514.9116.9115.042015.054015.120215.197015.6413
x 3 1515.0115.0115.042015.054015.120215.197015.6413
x 4 55.0025.0024.99874.99855.00535.02305.1795
x 5 109.9810.9810.043310.055510.114910.174010.4618
x 6 55.0195.0195.03225.04035.08035.12215.2841
x 7 55.0145.0145.01115.01525.03465.05195.1777
SSE------0.00670.01100.05100.12871.2117
TER------0.94240.92630.84570.76060.2953
RER------0.91010.88380.75010.6008−0.2626
Table 3. Reconciled results of different methods with three gross errors.
Table 3. Reconciled results of different methods with three gross errors.
StreamTrueOriginal Meas. Meas. with Gross Error ProposedXieWelschCauchyFair
x 1 54.9954.9955.00805.01645.05745.03234.3533
x 2 1514.9116.9115.039015.055215.142515.180215.1833
x 3 1515.0115.0115.039015.055215.142515.180215.1833
x 4 55.0025.0024.99914.99845.00375.02445.3066
x 5 109.9810.9810.039910.056810.138810.15589.8767
x 6 55.0195.0195.03205.04045.08145.12355.5233
x 7 55.0143.5145.00805.01645.05745.03234.3533
SSE------0.00580.01150.07310.10711.2866
TER------0.97430.96420.91110.89540.3474
RER------0.96410.94700.86210.84460.1268
Table 4. Reconciled results of different methods for nonlinear constraints.
Table 4. Reconciled results of different methods for nonlinear constraints.
StreamStandard DeviationTrueMeas.ProposedXieWelschCauchyFair
x 1 0.54.51244.53604.53784.52804.53794.57964.4727
x 2 0.65.58195.90705.57545.56555.53315.53605.6514
x 3 0.21.92601.80741.92211.92231.92001.91531.9321
x 4 0.21.45601.46531.46531.48421.49241.50961.4655
x 5 0.54.85454.54914.81564.80104.80834.78824.7870
u 1 --11.070--11.098811.117811.193911.207910.9076
u 2 --0.6147--0.61430.61600.61870.61680.6104
u 3 --2.0504--2.04692.03492.03172.03452.0372
SSE------0.00310.00670.02220.03330.0377
TER------0.89500.81920.77510.66310.7994
RER------0.88050.80080.73240.59300.7693
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Luan, Y.; Jiang, M.; Feng, Z.; Sun, B. Estimation of Feeding Composition of Industrial Process Based on Data Reconciliation. Entropy 2021, 23, 473. https://doi.org/10.3390/e23040473

AMA Style

Luan Y, Jiang M, Feng Z, Sun B. Estimation of Feeding Composition of Industrial Process Based on Data Reconciliation. Entropy. 2021; 23(4):473. https://doi.org/10.3390/e23040473

Chicago/Turabian Style

Luan, Yusi, Mengxuan Jiang, Zhenxiang Feng, and Bei Sun. 2021. "Estimation of Feeding Composition of Industrial Process Based on Data Reconciliation" Entropy 23, no. 4: 473. https://doi.org/10.3390/e23040473

APA Style

Luan, Y., Jiang, M., Feng, Z., & Sun, B. (2021). Estimation of Feeding Composition of Industrial Process Based on Data Reconciliation. Entropy, 23(4), 473. https://doi.org/10.3390/e23040473

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop