Abstract
In order to improve the control on sulphur content at the endpoint of Kanbara Reactor (KR) desulphurization process, the case-based reasoning (CBR) method based on mechanistic model correction is used to predict the endpoint sulphur content of molten iron. First, the KR desulphurization process is analysed to determine its kinetic mechanistic model, and partial derivatives are obtained for different attributes. Then, according to the analysis ofthe attributes, the corrected model is determined by the method of selecting the attributes, determining the reasonable weights and fitting the calculation results. Finally, the CBR method, the Back Propagation Neural Network (BPNN) and the corrected model are used to predict the endpoint sulphur content at KR desulphurization. The experiment results indicate that the prediction accuracy of corrected model is significantly higher than that of BPNN and mechanistic model.
Introduction
As the sulphur is generally harmful to most steel grades in the steel materials, it is necessary to strictly control the sulphur content in steel so as to ensure the mechanical and process properties of steel. At present, the sulphur is mainly removed by molten iron pretreatment for most steel grades, and Ladle Furnace (LF) refining is required by only a few of steel grades for deep desulphurization. In the desulphurization process of molten iron pretreatment, Kanbara Reactor (KR) stirring method and injecting magnesium method are the most widely used. With the increased power of the stirring equipment and the improvement of the process, the KR stirring method has become the most mainstream method for the desulphurization of molten iron in China and Japan.
The purpose of the KR stirring method is to remove the sulphur element in the molten iron and to supply molten iron that has qualified sulphur content for the converter steelmaking, so it is very important to effectively improve the control on the endpoint sulphur content in molten iron in KR desulphurization. However, in the actual production, the endpoint sulphur content in molten iron in the KR desulphurization process is obtained by sampling and testing, which not only prolongs the process time, but also increases the sampling and testing cost. Therefore, more and more researches attempt to solve the control problem of the endpoint sulphur content by means of prediction.
The case-based reasoning (CBR) method based on mechanistic model correction is proposed to control the endpoint sulphur content in KR desulphurization. First, the mechanistic model of KR desulphurization process is established. Then, the influence degree of each attribute on the objective function is analysed by obtaining the partial derivative for the attribute in the mechanistic model, and the corrected model is determined by selecting the attribute, determining the reasonable weight and fitting the calculation results. Finally, the prediction accuracy of the corrected model is verified using actual production data and shows good result. This model can provide a reference for the actual production to improve the control accuracy of the endpoint sulphur content.
Section ‘literature review’ reviews the related literatures. CBR method based on Mechanistic Model correction is proposed in Section ‘CBR method based on mechanistic model correction’. Experiment and results are obtained in Section ‘Experiment and results’, and conclusion is made in Section ‘Conclusions’.
Literature review
In the study of steel production, the models for sulphur content prediction at the endpoint of the smelting process include mechanistic model and data model.
As for mechanistic model, its research on desulphurization reaction mainly focuses on process optimization. Shougang research institute of technology established a kinetic model of the 60% Mg and 40% CaO injection desulphurization. By analysing the model, the desulphurization process of high sulphur molten iron had three stages, incubation stage, rapid desulphurization stage and slow desulphurization stage [1]. Wuhan University of Science and Technology established a mathematical model about the variation of sulphur content in liquid steel during a 300 t RH refining process based on the consideration of thermodynamics, kinetics of desulphurization and conditions of practical production. The effects of initial sulphur content in the liquid steel, powder injection rate, particle size of powder, flow rate of RH driving gas and initial temperature of the bath on the temperature drop were taken into account. Sulphur content in the steel bath during desulphurization process and the final sulphur content in the liquid steel were analysed. The model can be used to predict desulphurization process [2]. The accurate parameter setting and the data collection of all influencing factors are required in the application of the mechanistic model, which greatly limits the application of this model in the actual production.
As for the data model, there are many studies on the endpoint compositions and temperature prediction of the monomer process. To improve the control precision of endpoint because of the constraints of economy and technology, Yue proposed a method to solve the problem that the small and middle converters unable to introduce the sublance detection technology. The method combined the pedigree cluster and neural network. Simulation results show that the multi-neural network model has better prediction results [3]
Extensive attention has been given to the application of CBR methods to solve related problems in the research field of data model. The CBR is a machine learning method in which the similarity among the influencing factors is used to solve the problem. A better solution can always be obtained for cases that have too complicated/weak correlation between influencing factors, or cases whose influencing factors are seriously insufficient. An integrated CBR model composed of rough search and delicate search, is proposed to predict the endpoint temperature of molten steel in AOD [11]. Combining with the ultra-fast cooling process, a developed CBR model was proposed, which mainly improves the case representation, similarity relation and retrieval module. Retrieval process was simplified and retrieval efficiency is improved apparently by the windmill retrieval algorithm. The proposed CBR model was used for predicting the case of cooling strategy and its capability is superior to the traditional process model [12]. A combined method of CBR and Bayesian belief network (BBN) has been proposed. The evaluation of the reliability of cases is conducted by applying BBN for the assessment which the CBR method lacks in [13]. The establishment based on CBR and radial basis function neural network of coke oven flue temperature intelligent prediction model, realize the real-time prediction of the temperature, and help to realize the coke oven production process of intelligent optimization control [14]. A two-step CBR method using genetic algorithm to find the optimal attributes subset based on the evaluation method of Correlation-based Feature Selection was proposed for predicting the endpoint phosphorus content in BOF efficiently [15]. To be simulate dynamics of cold tandem rolling processes, a novel hybrid intelligent dynamic modelling approach is proposed based on the combination of a linearized state space model derived from various mechanism equations, a CBR algorithm for multi-state space models selection, a genetic algorithm for optimization of case attributes, an adaptive fractal filtering algorithm for the identification of state space model parameters, a neural network-based simulation error compensation model for the strip exit velocity [16]. CBR based on two-step retrieval approach and the correlation-based feature weighting method was proposed for predicting end temperature of molten steel in LF [17].
At present, the research of CBR mainly focuses on the case classification processing of case store and the similarity algorithm improvement, but the necessary mechanistic analysis for the researched problems is absent. Better prediction accuracy can be obtained using the mechanistic model to correct the similarity of influencing factors in the CBR method. The mechanistic model is established to characterize the correlation between the influencing factors and the target factors. If the mechanistic model characterizing the correlation between the influencing factors and the target factors is established to correct the solution provided in the CBR method for the similarity among the influencing factors, a better prediction accuracy will be obtained. The research in this paper is based on this idea.
CBR method based on the mechanistic model correction
Arithmetic process of CBR
The CBR is a method to obtain solutions using the similarity between current problems and existing cases in the case store. The standard CBR model often consists of four processes: case representation, case retrieval, case reuse and case retain. Among them, case retrieval is the core of CBR.
Before the case retrieval, it is generally necessary to pre-process the case store and the current problem. The 0–1 normalization is used for data preprocess in this paper.
A typical retrieval method is to perform retrieval with the similarity between attributes. The grey distance similarity calculation method is used in this paper. Grey distance is an analysis method based on grey relational space theory, which integrates the characteristics of distance space and topological space. Its basic idea is to judge the degree of correlation between factors according to the similarity of the magnitude changes between the factor curves. Because it does not need the typical distribution law of data, and can combine qualitative analysis with quantitative analysis, it has been widely used in different fields.
Given that the problem case is
Then the formula for calculating the similarity of the grey distance between the problem case
In formula (1),
The case results obtained by the similarity calculation are reused to solve the problem case. The cases can be sorted by similarity which is calculated using the similarity calculation method. If there are k cases with the maximum similarities, the result of the problem case is computed as follows:
In which Gi is the similarity between the new problem and the i-th case, and Ti is the result in the i-th stored case. The value of m will be 4 in the paper.
The CBR method is based on the similarity between the case attributes to obtain the result of the problem case, but the consideration of the relationship between the attribute and the result is absent. As far as the most practical problems are concerned, this correlation cannot be ignored. Therefore, a method is proposed in this paper to correct CBR calculation process using the mechanistic model.
The correction on the CBR method using the mechanistic model
The CBR is to realize case reuse to solve new problems by analysing and calculating the attribute similarity between existing case and new problem. Only the influence degree of different attributes on the target in the attribute weight of the similarity algorithm is considered in this method, but the inspection on the degree of necessity of attributes is absent. At the same time, it is impossible to compare and analyse the coupling relationship between attributes. Therefore, it is necessary to study the degree of necessity of related attributes and the coupling relationship between different attributes by means of the mechanistic model of the studied problem, and then correct the CBR method to obtain more accurate results.
The objective function
The correction by mechanistic model on the CBR method includes the following three aspects: Attribute selection: By comparing the size of each attribute Determining the weight: the weight determination method based on the mechanistic model is to determine the corresponding weight according to the size of each attribute Fitting of calculation results: for the new problem, the CBR method and the mechanistic model are used for calculation respectively, and then the result of the two methods are fitted. The related fitting parameters are determined by the method of case store verification, then more accurate results are obtained. That is: the correction model is constructed as
After the correction by the mechanistic model, the specific process of the CBR method is shown in Figure 1: Process of CBR model corrected by mechanistic model.
Experiment and results
Experiment data
Production performance data collected from a steel plant in the process of KR desulphurization production were used as the data source for CBR calculation and verification. The KR mechanical stirring method is a method in which the refractory stirrer is immersed into the molten iron pool by a certain depth, and then the desulphurizing agent is uniformly dispersed into the molten iron by the vortex arising from stirring, so that the solid and liquid are fully contacted, and high-efficiency desulphurization is achieved. The lime-based desulphurizer is used in KR desulphurization in this plant (CaO content is 80%, SiO2 content is 3%, fluorite content is 8%, other ingredients content is 9%), and particles with diameter of 0.3–1.2 mm account for 80%. The operation flow of the KR desulphurization process is: entry → first sampling and temperature measurement → first slag removal → desulphurizer is added and stirred for desulphurization → stop stirring → second slag removal → second sampling and temperature measurement → exit. In this paper, 5701 sets of production data for 3 months of this plant are selected to conduct experiments. Four thousand sets of data are randomly selected as case store data, and the remaining 1701 sets of data are used as test data.
Data distribution of influencing factors on endpoint sulphur content at KR desulphurization.
Note: The [S] content in molten iron is measured as mass percentage.
Based on this data, the sulphur content at KR endpoint is predicted.
The predictive results are obtained based on back-propagation (BP) neural network. This neural network is implemented in Matlab. It is a three-layer network with input, hidden, and output layer. The input layer consists of 5 nodes representing the 5 attributes. The hidden layer also has 5 nodes with transig as transfer function. The output layer is composed of just one node representing the predicted the sulphur content with purelin as transfer function. The max epoch is 1500. The hit rate of predicting is illustrated in Figure 2(a). Hit rate of endpoint sulphur content prediction in KR process ((a) BPNN; (b) the ordinary CBR).
According to the CBR method described in Section ‘Arithmetic process of CBR’. The final predictive results were obtained based on ordinary CBR that used GRD with average weight and on a case base that contained 4000 sets of data. The distribution of its hit rate is illustrated in Figure 2(b).
Hit rates of endpoint sulphur content prediction in KR processwith BPNN and ordinary CBR.
It can be seen from the comparison between BPNN and ordinary CBR that, BPNN, as a method of identifying data correlation between attributes and targets, has high prediction accuracy when the prediction error range is large, such as the prediction results within the error range of [−4, 4], [ −5,5] and [−6,6]. However, CBR, as a method of identifying attribute similarity between cases, has a high prediction accuracy when the prediction error range is small, such as the prediction results within the error range of [−1, 1] and [−2, 2]. This indicates that, for the prediction of endpoint sulphur content in KR desulphurization, only when both the relationship between the attribute and the target and the similarity of the attributes between the cases have been considered, overall better prediction result can be obtained.
Kinetic mechanistic model of KR desulphurization
In order to obtain a better prediction result, it is necessary to analyse the relationship between the endpoint sulphur content and its influencing factors by means of mechanistic analysis.
The KR stirring desulphurization reaction equation and the corresponding standard free energy are as follows:
Relationship between molten iron temperature and endpoint sulphur content at the point of desulphurization reaction equilibrium.
In the actual production, the average endpoint sulphur content in KR desulphurization is 6.44×10−5, which indicates that the reaction is far from equilibrium, so the endpoint sulphur content depends on the kinetic conditions of the reaction process.
Complicated mass transfer process and interfacial reaction are involved in the KR desulphurization process. According to the slag ion structure theory, the control element of the desulphurization reaction is the diffusion process of sulphur from molten iron to slag. The desulphurization rate is:
Considering the lime particles as spherical solid particles, the mass transfer coefficient is
Arrange the above formula, bring in the corresponding parameters:
The slagging agent composition and the slagging process of the KR desulphurization process in this plant are relatively stable, so
By fitting the actual production data, we obtain
Prediction result of CBR based on the mechanistic model correction
Based on the KR desulphurization kinetic model, as shown in formula (6), the partial derivatives of
According to the method in Section ‘The correction on CBR method using mechanistic model’ of this paper, the CBR method is corrected based on the above results: Attribute selection:
According to the KR desulphurization model, if the threshold is set
The influence of the molten iron temperature T on the endpoint sulphur content is negligible, because the temperature change of the molten iron has little influence on the kinetic conditions of the desulphurization reaction during the KR desulphurization process. The desulphurizer consumption
Through the attribute selection, the five influencing factors can be reduced to three influencing factors, namely the entry sulphur content Determine the weight
Relevant factors weight for five and three factors scenarios.
The CBR method is used to predict the endpoint sulphur content in KR desulphurization for the mechanistic model weight of five factors, the average weight of three factors and the mechanistic model weight of three factors, and the prediction results are as shown in Figure 3 (a–c). The mechanistic model is used to predict the endpoint sulphur content in KR desulphurization for the five factors and three factors, respectively. When the endpoint sulphur content is predicted by the three-factor mechanistic model, the other parameters in formula (6) are taken at the average value, and the prediction results are as shown in Figure 3(d,e). Hit rate of predicting endpoint sulphur content in KR process ((a) CBR.-5 factor-mechanistic model weight; (b) CBR.-3 factor-average weight; (c) CBR.-3 factor-mechanistic model weight; (d) mechanistic model-5 factor; (e) mechanistic model-3 factor).
Hit rates of predicting endpoint sulphur content in KR process under five scenarios.
It can be seen from Table 5 that, the prediction results of the CBR.−3 factor-mechanistic model weight are overall better than that of other CBR. At the same time, the prediction accuracy of the mechanistic model-3 factors is also higher than that of the mechanistic model-5 factors comprehensively. This indicates that the unnecessary influencing factors can be removed from the attribute selection and attribute weight determination based on the mechanistic model, thus avoiding the influence of fluctuation of this factor on prediction accuracy. At the same time, the relative importance of each attribute can be determined, thus effectively improving the accuracy of CBR method. Fitting of the calculation results
In order to further improve the prediction accuracy, a corrected model of the prediction results is constructed
The prediction results of the corrected model (7) on the endpoint sulphur content are as shown in Figure 4: Hit rate of predicting endpoint sulphur content in KR process using the corrected model.
Comparison of prediction hit rates of endpoint sulphur content in KR process using the corrected model and other methods.
It can be seen from the above test that, compared with other CBR methods, although the prediction hit rate of corrected model is lowered in smaller error range, its prediction hit rates in other error ranges are greatly improved. For example, in the error range from −3×10−5 to 3×10−5, the prediction hit rate of the corrected model has improved 5.18% than CBR.-3 factor-mechanistic model weight, and 7.41% higher than the ordinary CBR method. In the error range from −4×10−5 to 4×10−5, the corrected model has improved 11.11% than CBR.-3 factor-mechanistic model weight, and 13.05% higher than the ordinary case reasoning method. In the error range from −5×10−5 to 5×10−5, the corrected model has improved 10.23% than CBR.-3 factor-mechanistic model weight, and 12.93% higher than the ordinary case reasoning method. Considering that the actual measured accuracy of sulphur content of molten iron in KR desulphurization process is 1×10−5, and the reduction of prediction accuracy in error range from −1×10−5 to 1×10−5 and from −2×10−5 to 2×10−5 is acceptable.
As compared with the mechanistic model, although basically there is no change in the prediction hit rate within small error range of the corrected model, its hit rate in other error ranges is significantly improved. For example, in the error range from −3×10−5 to 3×10−5, the corrected model has improved 3.53% than mechanistic model-3 factors. In the error range from −4×10−5 to 4×10−5, the corrected model has improved 9.34% than mechanistic model-3 factors. In the error range from −5×10−5 to 5×10−5, the corrected model has improved 8.17% than mechanistic model -3 factors
As compared with BPNN, the prediction hit rate of corrected model in all error ranges has improved. For example, in the error range from −1×10−5 to 1×10−5, the corrected model has improved 4.23% than BPNN. In the error range from −3×10−5 to 3×10−5, the corrected model has improved 6.41% than BPNN. In the error range from −5×10−5 to 5×10−5, the corrected model has improved 3.82% than BPNN.
Conclusions
In this paper, a CBR method based on mechanistic model correction is proposed to predict the endpoint sulphur content of molten iron in the KR desulphurization process.
Based on the analysis of the kinetic mechanism of KR desulphurization process, and using the actual data to fit the relevant parameters, the kinetic mechanistic model of KR desulphurization process is established, and the partial derivative of each influencing factor is obtained to characterize its influence degree on the endpoint sulphur content.
Based on the kinetic mechanistic model and the analysis of the necessity of influencing factors, the CBR is corrected by the method of attribute selection, determining the reasonable attribute weight and the fitting of the calculation result, and finally, the corrected model of KR desulphurization is determined.
The CBR method (for the three factors, five factors, average weighting and mechanistic model weight, respectively), the mechanistic model (for the three factors and five factors respectively), BPNN and the corrected model are separately used to predict the endpoint sulphur content in KR desulphurization. The results show that the prediction accuracy of the corrected model is significantly higher than that of the BPNN and the mechanistic model. At the same time, the prediction accuracy of the corrected model is significantly improved as compared with the prediction results of all CBR methods in respect of the error ranges of the actual production process. The result is can be reference in instructing the production process.
Footnotes
Disclosure statement
No potential conflict of interest was reported by the authors.
