Abstract

Annular furnaces have multivariate, nonlinear, large time lag, and cross coupling characteristics. The prediction and control of the exit temperature of a tube billet are important but difficult. We establish a prediction model for the final temperature of a tube billet through OS-ELM-DRPLS method. We address the complex production characteristics, integrate the advantages of PLS and ELM algorithms in establishing linear and nonlinear models, and consider model update and data lag. Based on the proposed model, we design a prediction control algorithm for tube billet temperature. The algorithm is validated using the practical production data of Baosteel Co., Ltd. Results show that the model achieves the precision required in industrial applications. The temperature of the tube billet can be controlled within the required temperature range through compensation control method.

1. Introduction

A seamless tube should be heated to a given temperature in an annular furnace prior to piercing. The heating quality of the tube billet directly influences the quality of the seamless tube. In seamless tube production, the main evaluation criterion for tube billet heating quality is the final exit temperature of the tube billet. The heating reaction mechanism of an annular furnace is complicated. It has nonlinear, large time lag, and uncertainty characteristics. The temperature distribution on the surface of a tube billet within the furnace cannot be directly measured and controlled [1]. Thus, controlling the outlet temperature of the billet is difficult.

Measuring the temperature of the billet in the furnace is also difficult. A common method is to establish a temperature prediction model for the billet. Several scholars have proposed mechanism models for billet temperature in furnaces. Chen et al. [2] established a temperature drop model of a tube billet during the transfer of the tube billet from the heating furnace to the rolling mill based on the partial differential equation of heating transfer. Jaklič et al. [3] established a mathematical model of deformation and heat flow during rough rolling. Considering that mechanism modeling is complex and modeling consumes much time, a model is usually established with a single furnace type and under many restricting conditions. Zhang et al. [4, 5] estimated the tube billet exit temperature through dynamic modeling of furnace temperature; however, the prediction error was too large to meet the requirements of the heating process. Wick [6] applied Kalman filter technique to estimate the temperature distribution of a tube billet inside a heating furnace. The disadvantage of this method is that the surface temperature of the tube billet inside the furnace must be measured, which is difficult to perform in practical production. Xiao et al. [7] employed production data and applied PCR method to establish a soft measurement model of tube billet temperature inside a furnace. However, this model has poor prediction precision because of the strong nonlinearity of production data. Chen and Chai [8] designed a preprocessing system for production process data. This system can predict several variables that are difficult to measure through the use of a self-adapting fuzzy-neural network. Cui and Ding [9] established a soft measurement model of tube billet temperature based on RBF neural network; however, model update was not considered. Iwamoto et al. [10] designed an automatic control system for a tube billet reheating rotary hearth furnace. The system consists of a component that calculates the tube billet temperature and a component that calculates the optimal furnace temperature set point. In large Chinese steel companies, such as Baosteel Co., Panzhihua Iron and Steel Co., Ltd., Anshan Iron and Steel Co., Ltd., and Capital Iron and Steel Co., Ltd., temperature prediction models are utilized in several heating furnaces. However, their prediction models are almost entirely engineering models imported from abroad. Therefore, these models are difficult to maintain and transplant, and the costs of doing so are high.

With the development of configuration software and database technique, increasing amounts of production data are being collected and stored. Therefore, data-driven modeling and control methods are eliciting more and more attention. He et al. [11] and Lv et al. [12] established data models for the Ladle furnace through data-driven methods. In the present work, the production data of an annular furnace were obtained from the seamless tube subcompany of Baosteel. Industrial process data contain noise, which reduces the modeling accuracy of the extreme learning machine (ELM) algorithm. By contrast, the capability of the partial least square (PLS) algorithm to process linear relevant data is suitable. Moreover, some cyclical changes occur in production. Thus, the model should be updated in real time. Online sequential (OS) ELM and recursive PLS algorithm can realize model update. Online sequential extreme learning machine dynamic recursive partial least square (OS-ELM-DRPLS) algorithm was proposed in this study. A tube billet temperature prediction model based on this algorithm was established, and a strategy for the optimization and control of tube billet temperature was proposed based on this model. OS-ELM-DRPLS not only has the advantages of the OS-ELM algorithm (e.g., rapid nonlinear modeling and update) but also has the capability of the RPLS algorithm to process linear relevant data. The large time lag and reduced model precision were solved through dynamic processing. The algorithm is easy to implement with advanced computer language. Configuration software, such as WINCC6.1, can be utilized to compile the modeling and control algorithms into special modules for use in industrial sites. Simulation and actual experiments prove that tube billet temperature can be predicted and controlled within the scope of the production process requirements with the established temperature prediction model and the proposed strategy of optimization and control based on OS-ELM-DRPLS algorithm.

2. Annular Furnace

An annular furnace is a type of rotary hearth furnace utilized for tube billet heating. It consists of the furnace shaft and auxiliary equipment for charging and discharging. The furnace shaft consists of a fixed furnace roof, a ring-type tunnel surrounded by a fixed furnace roof wall, and a circular ring-like rotary hearth, as shown in Figure 1.

External and inner ring seams are placed between the fixed furnace wall and rotary hearth in an annular furnace. Internal and external water seal tanks are arranged beneath the external and inner ring seams to maintain normal temperature and pressure in the furnace cavity and prevent external cold air from entering the furnace cavity. Fuel gas and combustion-supporting air are blown into the furnace through burning nozzles mounted on the external and internal walls or furnace roof to make the gas burn within the furnace and heat the tube billet. The fumes produced by gas burning within the furnace move conversely through the rotary hearth to the tail end of the soaking zone, enter the flue and chimney outside the furnace, and exit to the atmosphere. The external wall of the furnace has charging and discharging furnace doors, in which a charger and discharger are placed, respectively. Charging and discharging proceed simultaneously. When a tube billet is placed in the furnace, the bottom rotates at a certain angle. Tube billets follow a radial layout inside the furnace and are arranged either in a single row or in multirows.

The furnace cavity is divided into preheating, heating, and soaking zones according to the heating process of the tube billet in the annular furnace. Burning nozzles are not mounted in the preheating zone. A flue opening is arranged on the side wall near the charging furnace door in this zone. High-temperature exhaust gas flows toward the opposite direction of hearth rotation and exits into the atmosphere through the flue opening in the heating and soaking zones. During the flow process of high-temperature exhaust gas, the tube billets in the preheating zone are mainly convection heated. The length of the preheating zone accounts for approximately one-fourth of the peripheral length of the annular furnace. Temperature differences between the surface and center and between both ends exist in the tube billet rapidly heated in the heating zone. To reduce the temperature difference of the tube billets and eliminate their male-female faces, the tube billets must be heated in the soaking zone. The length of the soaking zone is approximately three-twentieths of the peripheral length of the annular furnace. In addition, no tube billet and burning nozzle are present between the charging and discharging furnace doors. A partition wall is placed in the middle. The distance between the charging and discharging furnace doors is approximately one-tenth of the peripheral length of the annular furnace.

3. Dynamic Nonlinear PLS Method

Given that a linear PLS model cannot correctly describe the nonlinear relationship between independent variable and dependent variable ( is the variable matrix that affects the heating temperature of the tube billet and is the variable matrix of the heating temperature), nonlinear PLS (NLPLS) method is required. Wold et al. extended the PLS method to the nonlinear field [13, 14]. Two feasible methods exist in nonlinear PLS methods. One method is to perform array extension for the input matrix, introduce several nonlinear terms of the original variable (e.g., the square term), and then regress the extended input and output matrix through PLS method. If prior knowledge on the relationship of original input variable does not exist, this method cannot be utilized as a reference in the selection of the combined mode and may lead to an oversized dimension of the input matrix and difficulties in processing. The other method is to reserve the linear external model of the PLS method. The internal model is nonlinear. The effect of various input variables on the final tube billet temperature has a different time lag because of the large time lag characteristics of tube billet heating. Accurately predicting the tube billet temperature through traditional PLS method is difficult. In this study, dynamic PLS method was utilized to calculate the lag time of various input channels and significantly improve the model’s precision. The algorithm is as follows:where are the ratio of lag time to sampling period for sampling variables .

(1) The external relation model iswhere is the number of reserved eigenvector; and are the score vectors of and , respectively; and are the load vectors of and , respectively; and are the score matrixes of and , respectively; and are the load matrixes of and , respectively; and and are the fit residual matrixes of and , respectively.

(2) The internal relation model iswhere is the nonlinear function and is the residual.

Given that a neural network is capable of nonlinearity fitting during the modeling of the batch process, nonlinear MPLS method, where the internal model adopts a neural network, has gained extensive applications. Considering that a traditional feed-forward neural network adopts a gradient learning algorithm during training, the parameters in the network need iteration and updating. Training not only consumes much time but also easily results in issues of local minimum and excessive training [15].

4. OS-ELM-DRPLS Algorithm

4.1. ELM Algorithm

In supervised batch learning, the learning algorithms employ a finite number of input-output samples for training [1622]. For arbitrary distinct samples , where is a input vector and is a target vector, if a single hidden layer feed-forward neural network (SLFN) [2326] with hidden nodes can approximate these samples with zero error, then , , and exist such thatIn the expression above, , , and are the learning parameters of the hidden nodes (the weight vector connecting the input node to the hidden node and the threshold of the hidden node) randomly selected according to the proof provided by Huang et al. is the weight connecting the th hidden node to the output node. Error term is added to avoid overfitting the noise in the data. is the output of the th hidden node with respect to input , and is the number of hidden nodes that can be determined by trial and error or prior experience. Equation (4) can then be written compactly aswhereIn the expressions above, is the hidden layer output matrix of the network; the th column of is the th hidden node’s output vector with respect to inputs , and the th row of is the output vector of the hidden layer with respect to input . The hidden node parameters and need not be tuned during training and may simply be assigned with random values. Equation (5) then becomes a linear system, and the output weights are estimated aswhere is the Moore-Penrose generalized inverse of hidden layer output matrix [27, 28].

4.2. OS-ELM Algorithm

In actual applications, training data may arrive chunk by chunk or one by one. Hence, the batch ELM algorithm has to be modified and made online sequential for this case [29, 30].

Output weight matrix provided in (7) is a least-squares solution of (5). We consider the case, where is the number of hidden nodes. Under this condition, of (7) is provided by

If is singular, one can make it nonsingular by selecting a small network size or increasing data number in the initialization phase of OS-ELM. Substituting (8) to (7), becomes

Equation (9) is the least-squares solution to . Sequential implementation of (9) results in OS-ELM [31].

Given a chunk of initial training set and , if the batch ELM algorithm is employed, the solution of minimizing , which is given by , where , must be considered.

We consider another chunk of data , where is the number of samples in this chunk. The problem involves minimizing

Considering both and , output weight becomes

For sequential learning, should be expressed as a function of and not as a function of dataset . can be written as

Combining (11) and (13), is obtained withwhere .

When the th chunk of datasetis received, where and denotes the number of samples in the (k+1)th chunk, we have

rather than is utilized to compute from in (16). The update formula for is derived with the Woodbury formula

We let . The equation for updating can be written as

Equation (18) provides the recursive formula for .

4.3. OS-ELM-DRPLS Modeling Steps

The difference of nonlinear DRPLS modeling method based on OS-ELM from linear PLS method is that the former employs ELM to establish the internal nonlinear model and updates the internal and external models. This method reserves the linear external model, extracts the attributive information of the process through PLS, eliminates the colinearity of data, reduces the dimension of the input variable, and then adopts ELM to establish a nonlinear internal model between the input score vector matrix and the output score vector; the nonlinear processing capability of the internal model is enhanced. Thus, OS-ELM-DRPLS method has the advantages of PLS and ELM (i.e., the robustness and feature extraction capability of PLS method and quick nonlinear processing capability of ELM as well as precision accuracy through real-time model update).

The modeling and testing steps of nonlinear DRPLS method based on OS-ELM are as follows.

(1) Two standardized data matrices, and , are assigned. The dynamic nonlinear PLS regression model can be expressed as follows: where are the ratio of lag time to the sampling period for sampling variables .

(2) The batch data of the batch process are deployed, cross-validation method is implemented to determine the number of latent variables, and linear PLS method is applied to calculate score vector matrices and and load vector matrices and for modeling samples and :

(3) A node number is assigned to the ELM hidden layer and activation function (e.g., sigmoid function), ELM is employed to establish a nonlinear model between internal models and , and is obtained, where is the nonlinear function indicated by ELM. The hidden nodes in SLFN transform the feature space into another feature space. The original ELM regards the number of nodes as a parameter to be defined. We increase the number of hidden nodes until stop criteria (e.g., residual error reduction) are reached. Meanwhile, the number of hidden nodes is less than .

(4) When one new batch of data , is obtained, PLS decomposition is performed, and score and load vectors are obtained:

According to (18), the OS-ELM algorithm is adopted to update the output layer weight value and the internal model. Weighted mean is conducted on the load matrix of the external model, and external RPLS update, where is the weight value factor, is achieved. The above steps are repeated, and model update is conducted for every batch:

(5) Testing data are utilized to verify the model’s precision. PLS decomposition is conducted on testing data , and score vector is obtained:

is introduced into the OS-ELM model. is obtained, and the model prediction value is determined through .

(6) A system error is obtained by comparing with the practical output. can vary within . After each variation, are substituted back to (19) for calculation and to obtain an estimation error. Finally, one is obtained by exhausting a group of optimal values to minimize the model estimation error:

Model parameters and of the OS-ELM-DRPLS model are then obtained through the aforementioned calculation.

5. Prediction and Control of Tube Billet Heating Quality Based on OS-ELM-DRPLS Model

5.1. Introduction of the Site and Selection of Measuring Points

In the seamless tube subcompany of Baosteel, the designed output of an annular furnace was 160 t/h. Its intermediate diameter was 35 m, and the effective width of hearth was 4.5 m. The hearth was divided into six burning control sections. The diameter of the heated tube billet was 178 mm. The temperature upon entering the furnace was 20°C, and the maximum temperature upon leaving the furnace was 1280°C. In the annular furnace, mixed gas that consists of 52% blast furnace gas, 13% converter gas, and 34.8% coke oven gas was utilized. The composition of the blast furnace gas was 23.5% CO, 2% H2, 19.5% CO2, and 35% N; the composition of coke oven gas was 53% H2, 29.2 CH4, 2.8% weight carbon hydride, 7.5% CO, 2.0% CO2, 0.6% O2, and 4.4% N2. The composition of converter gas was 56% CO, 24% N2, and 19.7% CO2. The specific technical parameters are shown in Table 1.

The final exit temperature of the tube billet was predicted through OS-ELM-DRPLS method. First, the variation in the tube billet temperature was reflected, and the measured variables were easily obtained. On one hand, gas could not be obtained through the peep holes because the peep holes in the furnace were closed. On the other hand, opening of the furnace door to obtain gas affects the testing precision because of the absorption of cold air. Therefore, the measuring points in the site were set at the lighting holes of burning nozzles in the external surrounding furnace walls. Six flow rate detecting points were set for the burning nozzles. Nine thermocouples were mounted in the six working sections to measure the temperature inside the furnace cavity. The specific positions of flow rate and furnace temperature are shown in Figure 2. Fifteen measuring variables were selected to predict the final tube billet exit temperature. were measuring points for numbers 1–6 burning nozzle flow rate, and were measuring points for numbers 1–9 furnace cavity temperature. The variable table is shown in Table 2. After selecting the measuring variables and gathering site production data, OS-ELM-DRPLS method was applied to the prediction model of tube billet temperature.

5.2. Establishment and Checking of the Tube Billet Final Temperature Prediction Model

The production data for 70 pieces of tube billets produced by Baosteel in March, 2013, were utilized. The first forty samples were utilized as training data to establish the prediction model of tube billet final temperature. The last thirty samples were used for model update. Lump update was employed. Every group of five was considered a lump. The model was updated. The last thirty samples acted as the testing samples to check the precision of model prediction. Prior to modeling, data were expanded, they were standardized-processed, and they underwent cross checking. The number of PLS potential variables was determined to be 4. The number of ELM hidden layer nodes was 10. The excitation function was a sigmoid function. The ratio of lag time of in (19) was calculated in formulas equation (25).

The same data were tested with RPLS, RBF-PLS, and OS-ELM-DRPLS methods. The predicted mean square error and modeling time are shown in Table 3. Although the three methods meet the requirements of industrial application, OS-ELM-RPLS method exhibits better expansion capability, prediction precision, and nonlinear fitting capability for industrial application than RPLS method. Compared with nonlinear RBF-PLS method, the training time in OS-ELM-RPLS method is shorter. OS-ELM-RPLS method can achieve rapid modeling and model update and is significant to the intermittent production processes, such as tube billet heating:

The unit of [] is the sample time. Figure 3 shows a comparison between regression data and practical modeling data using RPLS and OS-ELM-DRPLS models. The maximum error was 6.9°C and the mean error was 2.3°C, which meet the requirements of the production site. To further verify the accuracy of the model, new data were introduced into the model and substituted into the following equation to obtain estimation value of the new data. The comparison with is shown in Figure 4. The maximum error was 9.8°C and the mean error was 3.1°C, which meet the requirements of the production site:

5.3. Predicted Control of Tube Billet Final Temperature

The aforementioned data indicate that tube billet exit temperature often fluctuates in the temperature range of 1200°C to 1300°C and often deviates from the ideal piercing temperature (1270°C). Such condition degrades the quality of the tube. The tube billet exit temperature should be controlled within the temperature range of 1255°C to 1295°C. The gas flow rate can be adjusted according to the prediction, practical measuring, and target temperatures. Its control period was 1 s. An ELM model predicted controller (EPC) was designed for the annular furnace system with the OS-ELM-DRPLS model predictor (EMP), as shown in Figure 5.

The basic operating principle of predictive control is to generate a sequence of control signals at each sample interval that optimize the control effort to follow the reference trajectory exactly [32, 33]. The ELM model predictive control law was obtained by minimizing the following predictive performance criterion:where

is the predictive output horizon, is the input reference signal at discrete time , and is the step-ahead prediction of . In general, is selected to include all responses that are significantly affected by the present control. In this study, is .

The control, , was obtained from the optimization of the cost function (29) based on gradient descent method; that is,whereTo reduce the computational load of EPC, we let . The EPC controller is expressed in the formwhere

A schematic of the proposed PLC-based temperature control system is shown in Figure 6. The actual temperature control system of the annular furnace is depicted in Figure 7. SIMATIC S7-400 was selected as the PLC of the control system. The entire system is mainly composed of a PLC master station, a remote I/O station, an operator station, a programmer and communication bus, and other components. The main modules of PLC include nine slot bases (UR2), a 4 A power supply module (PS407), a central processor (CPU416-2DP), 1 M memory card, and a network communication module (CP443-1). The main modules of I/O expansion include a power supply module (PS307), an interface module (IM153-1), a digital input module (SM321 DC24V × DI16), a digital output module (SM322 DC24V × DO16), a counter function module (8CH FM350-2), an eight-thermocouple input module (SM331), an eight-RTD input module (SM331), and a four-output module (SM332). The main modules of the workstation include a CPU, (Intel Core i7-930, 2.8 GHz × 4), hard disk (WD 2TB), memory (Kingston 8 GB), color LED (24′′, 1280 × 1024 resolution), and a net card (Siemens 10/100 MB). The main module of communication includes Ethernet SINEC H1 and field bus PROFIBUS-DP. The main Software programs are Windows 2003 Prof, STEP7 V5.4, and WINCC6.1.

The tube billet exit temperature should be controlled as best as possible within the temperature range of 1255°C to 1295°C. Thirty tube billets were controlled by ELM model predicted control. A thermocouple was “buried” in a tube billet. The temperature course of the tube billet with the buried thermocouple is shown in Figure 8. In position 30, the predicted temperature of the tube billet is 1211. OS-ELM-DRPLS-based predictive control algorithm was employed to make the tube billet reach the lowest required temperature (1255°C). By adjusting the input, the temperature of the tube billet reached 1265°C. The variation in tube billet temperature after introducing temperature compensation control is shown in Figure 9. The variation in tube billet temperature after introducing PID temperature control is shown in Figure 10. The effect of predicted control is better than that of the PID method.

Figure 9 shows that the tube billet exit temperature basically fluctuates in the range of [1255°C, 1295°C]; the tube billet heating quality is better than that before prediction control and meets the requirements of piercing production for tubes.

6. Conclusion

Measuring and controlling tube billet heating temperature are difficult because of the complex reaction mechanism during the heating process in an annular furnace. A tube billet final temperature prediction model was established in this study through OS-ELM-DRPLS modeling method. An OS-ELM-DRPLS-based predictive controller for the control of tube billet temperature was also systematically developed. The tube billet heating quality increased to a certain extent. This finding lays the foundation for the improvement of seamless tube quality. After the developed model was compiled into a universal module through the advanced computer language of the configuration software, the modules not only assisted in production by guiding front line workers to operate manually but also formed a perfect close loop control circuit together with the heating furnace model and controller. Hence, tube billet heating quality was improved effectively. Experimentation proved that this method is feasible. This modeling method is also versatile and can be extended to other processes with a large time lag.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This research is supported by the National Natural Science Foundation of China (Grant nos. 61203214, 41371437, and 61304121) and Provincial Science and Technology Department of Education Projects, the General Project (L2013101).