CN106709197A - Molten iron silicon content predicting method based on slide window T-S fuzzy neural network model - Google Patents

Abstract

The invention discloses a molten iron silicon content predicting method based on a slide window T-S fuzzy neural network model and belongs to the field of industrial process monitoring, modeling and simulation. The method includes: selecting a T-S fuzzy neural network to serve as the basic model of the predicting; adding a slide window model on the basis of the neural network, wherein a training sample set can be updated constantly to well track the variation trend of silicon content; selecting 11 parameters, which affect the molten iron silicon content most, according to practical experience and mutual information calculation to serve as the input of the model, and using the molten iron silicon content as the output of the model; normalizing training samples to train the model, and using the trained model to predict the molten iron silicon content during production. The molten iron silicon content fluctuates fiercely and is unpredictable due to the features of time varying, dynamicity, non-linearity, high inertia and multi-dimensionality of a blast furnace during iron making. Compared with the prior art, the molten iron silicon content predicting method is high in precision, low in mean-square error and applicable to online real-time predicting.

Description

Molten iron silicon content Forecasting Methodology based on sliding window T-S fuzzy neural network models

Technical field

The invention belongs to industrial process monitoring, modeling and simulation field, more particularly to a kind of modified EMD-Elman nerves The method of neural network forecast molten iron silicon content.

Background technology

Complicated mass-and heat-transfer, heterogeneous reaction and seal in blast furnace so that blast furnace has complicated time-varying, dynamic, non- Linearly, strong inertia and multiple dimensioned characteristic, the ironmaking processes allowed in blast furnace turn into one of most complicated industrial processes.Blast furnace Interior HTHP, deep-etching, strongly disturbing environment so that we are difficult directly to measure the hot situation in stove.But, in molten iron Silicone content is linearly related to furnace temperature, can reflect the quality of molten iron, and people are generally represented with the size of silicone content in molten iron The height of furnace temperature.Silicone content is too high to represent that furnace temperature is too high, can consume extra fuel, and can reduce the yield of iron；And silicon contains The low expression in-furnace temperature of amount is relatively low, may trigger and freeze the accidents such as cylinder.Therefore, for the stable smooth operation of blast furnace, people are needed stove In certain zone of reasonableness, the prediction of silicone content is just particularly important temperature control system.

The change of silicon is main in blast furnace is made up of three below reaction：

1/2CO+O₂=1/2CO₂

SiO₂+ CO=SiO (g)+CO₂

SiO (g)+C=Si+CO

By Arrhenius equation, temperature and concentration have a great impact to chemical reaction rate, in the phase of silicon Close in reacting, it can be seen that the influence of temperature, oxygen concentration and carbonomonoxide concentration to molten iron silicon content is maximum.Therefore, have Scholar establishes the mechanism model of molten iron silicon content prediction by dynamics and thermodynamics, and they pay close attention to the heat in course of reaction Amount and the conservation of mass.But, because mass-and heat-transfer complicated in blast furnace, phase change and chemical reaction so that modelling by mechanism is little The content of silicon can accurately be predicted.

Nowadays, the development of detection means allows that we measure substantial amounts of data, and the fast development of computer technology makes Obtaining us can in a short time carry out substantial amounts of computing, and these technological progresses cause that the modeling based on data-driven becomes more Easily, the also modeling method as main flow.The model based on data-driven for having existed has neutral net, linear regression, mixes Ignorant and supporting vector machine model etc., they have oneself respective strong point in some aspects.For example, the chaos grain that Jiang is proposed Subgroup optimized algorithm can well predict the temperature in successive reaction kettle in pharmaceuticals industry.But, these models are all to set up It is determined that data set on, be not suitable for industrial on-line prediction.

T-S fuzzy neural networks possess very strong adaptive ability, can automatically update model structure parameter, and can correct The membership function of fuzzy subset, can well be used for the molten iron silicon content in PREDICTIVE CONTROL ironmaking processes.Slided by combining Dynamic window model, model can more new training sample set at any time, and then update the parameter and coefficient of T-S fuzzy neural networks.It is sliding The characteristics of dynamic window T-S fuzzy neural networks can be well adapted for dynamic, the non-linear and strong inertia of ironmaking processes, in molten iron Good performance is shown in the prediction of silicone content.

The content of the invention

For the weak point of existing Silicon Content Prediction in Process of Iron model, it is proposed that one kind is based on sliding window T-S fuzznets The molten iron silicon content Forecasting Methodology of network.The method is modeled from sliding window model and T-S fuzzy neural network models, and is chosen 11 major parameters as model input, using silicone content as model output.The method have hit rate higher and compared with Small mean square error, can provide accurately prediction for the operating personnel of blast furnace, help their advance operation blast furnaces, make blast furnace steady Determine direct motion.The method is comprised the steps of：

Step one：T-S fuzzy neural network models are chosen, and combines sliding window model, for the prediction of silicone content.

Step 2：The input that 11 parameters are chosen as model, silicone content conduct are calculated by practical experience and mutual information Output.

Step 3：After by model initialization, with normalized training sample training pattern, the model that will be trained is used for silicon The prediction of content.

The structure of the T-S fuzzy neural networks described in step one is as follows：

T-S fuzzy neural networks are constituted by four layers, are respectively input layer, obfuscation layer, fuzzy rule computation layer and output Layer.Wherein input is fuzzy, and exports what is be to determine, and this represents that output is the linear combination of input.T-S fuzzy neural networks It is defined as follows：

WhereinIt is fuzzy subset, y_iIt is the calculating output of fuzzy rule.

(1) obfuscation layer is that it is defined as follows based on probability density function μ：

X in formula_jIt is input variable,WithIt is center and the width of probability density function, k is the dimension of |input paramete, n It is the quantity of fuzzy subset.

(2) fuzzy rule computation layer is made up of following formula：

(3) output layer is calculated by following formula：

The learning algorithm of the T-S fuzzy neural networks described in step one is as follows：

(1) error calculation：

Wherein y_dIt is actual value, y_cIt is predicted value, e is both differences.

(2) coefficient amendment：

In formulaIt is the coefficient of T-S fuzzy neural networks, and α is its learning rate.

(3) parameters revision：

Sliding window model principle described in step one is as follows：

Sliding window model is built upon in a kind of hypothesis, i.e., current output depends on current input, and is input into defeated Mapping ruler between going out can be obtained by historical data.According to this it is assumed that we preset a certain amount of training set Sample, is then continuously updated sample data and gives up earliest data point.With the slip of window, T-S fuzzy neural networks Its structural parameters can be constantly updated and newest predicted value is given.

The selection process of the input variable described in step 2 is as follows：

Mutual information is a kind of important method of test variable correlation, and Kraskov proposes a kind of k-NN methods can be very It is convenient to be used for calculating mutual information, comprise the following steps that described：

K is the number of neighbour given at the beginning in formula, and ψ is that Digamma functions can be expressed as：

ψ (x)=Γ (x)^-1dΓ(x)/dx

It obeys following iterative relation：

ψ (x+1)=ψ (x)+1/x

Ψ (1)=- C, C=0.5772156...

In order to obtain n_xAnd n_y, it is necessary to calculate sample z_iAnd z_jThe distance between d_i,j：

d_i,j=| | z_i-z_j||:d_i,j1≤d_i,j2≤d_i,j3...

||z_i-z_j| |=max | | x_i-x_j||,||y_i-y_j||}

As ε (i)=max { ε_x(i),ε_y(i) }, ε (i)/2 are taken as z_iWith the distance of k rank neighbours.Obviously, n_xI () is to x_i Distance is less than the number of the point of ε (i)/2, n_yI () is to y_iNumber of the distance less than the point of ε (i)/2.

Advise that we have chosen input of 11 variables as model by the practical experience of execute-in-place engineer, it Be respectively top pressure, top temperature, gas permeability, coal powder injection, oxygen enrichment percentage, full tower pressure difference, hot-blast pressure, hot blast temperature, hot air flow, Air humidity and previous stove silicone content.

Method for normalizing described in step 3 is as follows：

The present invention has advantages below：

1st, for the time-varying of blast furnace, dynamic, non-linear, strong inertia and multiple dimensioned characteristic in ironmaking processes, tool has been selected There are the T-S fuzzy neural networks of very strongly-adaptive, it has very strong learning ability, can find out latent between input and output In contact.Additionally, by adding sliding window, model can well track the variation tendency of molten iron silicon content, improve pre- The precision of survey.

2nd, calculated by operating experience and mutual information, selected top pressure, top temperature, gas permeability, coal powder injection, oxygen enrichment percentage, complete Tower pressure difference, hot-blast pressure, hot blast temperature, hot air flow, air humidity and previous stove silicone content etc. influence most on current silicone content 11 big parameters can make full use of the respective advantage of modelling by mechanism and data-driven modeling as the input variable of model.

Brief description of the drawings

Fig. 1 is the structural representation of T-S fuzzy neural networks,

Fig. 2 is the schematic diagram of sliding window,

Fig. 3 is the molten iron silicon content of 1000 stoves,

Fig. 4 is that this method predicts the outcome to molten iron silicon content.

Specific embodiment

The present invention proposes a kind of molten iron silicon content Forecasting Methodology based on sliding window T-S fuzzy neural networks, the party Method is comprised the steps of：

Step one：T-S fuzzy neural network models (such as Fig. 1) is chosen, and combines sliding window model (such as Fig. 2), be used for The prediction of silicone content.

Step 2：The input that 11 parameters are chosen as model, silicone content conduct are calculated by practical experience and mutual information Output.

Step 3：After by model initialization, with normalized training sample training pattern, the model that will be trained is used for silicon The prediction (such as Fig. 4) of content.

The structure of the T-S fuzzy neural networks described in step one is as follows：

T-S fuzzy neural networks are constituted by four layers, are respectively input layer, obfuscation layer, fuzzy rule computation layer and output Layer.Wherein input is fuzzy, and exports what is be to determine, and this represents that output is the linear combination of input.T-S fuzzy neural networks It is defined as follows：

WhereinIt is fuzzy subset, y_iIt is the calculating output of fuzzy rule.

(1) obfuscation layer is that it is defined as follows based on probability density function μ：

X in formula_jIt is input variable,WithIt is center and the width of probability density function, k is the dimension of |input paramete, n It is the quantity of fuzzy subset.

(2) fuzzy rule computation layer is made up of following formula：

(3) output layer is calculated by following formula：

The learning algorithm of the T-S fuzzy neural networks described in step one is as follows：

(1) error calculation：

Wherein y_dIt is actual value, y_cIt is predicted value, e is both differences.

(2) coefficient amendment：

In formulaIt is the coefficient of T-S fuzzy neural networks, and α is its learning rate.

(3) parameters revision：

Sliding window model principle described in step one is as follows：

Sliding window model is built upon in a kind of hypothesis, i.e., current output depends on current input, and is input into defeated Mapping ruler between going out can be obtained by historical data.According to this it is assumed that we preset a certain amount of training set Sample, is then continuously updated sample data and gives up earliest data point.With the slip of window, T-S fuzzy neural networks Its structural parameters can be constantly updated and newest predicted value is given.

The selection process of the input variable described in step 2 is as follows：

Mutual information is a kind of important method of test variable correlation, and Kraskov proposes a kind of k-NN methods can be very It is convenient to be used for calculating mutual information, comprise the following steps that described：

K is the number of neighbour given at the beginning in formula, and ψ is that Digamma functions can be expressed as：

ψ (x)=Γ (x)^-1dΓ(x)/dx

It obeys following iterative relation：

ψ (x+1)=ψ (x)+1/x

Ψ (1)=- C, C=0.5772156...

In order to obtain n_xAnd n_y, it is necessary to calculate sample z_iAnd z_jThe distance between d_i,j：

d_i,j=| | z_i-z_j||:d_i,j1≤d_i,j2≤d_i,j3...

||z_i-z_j| |=max | | x_i-x_j||,||y_i-y_j||}

As ε (i)=max { ε_x(i),ε_y(i) }, ε (i)/2 are taken as z_iWith the distance of k rank neighbours.Obviously, n_xI () is to x_i Distance is less than the number of the point of ε (i)/2, n_yI () is to y_iNumber of the distance less than the point of ε (i)/2.

Method for normalizing described in step 3 is as follows：

Embodiment

In steel manufacture process, blast furnace ironmaking is all good link of crucial importance, and the energy consumption of its consumption accounts for whole flow process 70%, thus blast furnace stable smooth operation be whole production process safe and highly efficient operation guarantee.Because the environment in blast furnace is extremely Badly, HTHP deep-etching so that conventional measurement means are difficult to carry out, and operating personnel are difficult to know the actual heat in blast furnace Situation, when tapping a blast furnace, molten iron loss amount of heat can not react actual furnace temperature.People are generally using molten iron silicon content come reacting furnace Interior actual state, thus the prediction of molten iron silicon content just seems of crucial importance, accurately prediction not only assists in operating personnel Rational regulation operating parameter, moreover it is possible to instruct blast furnace stable smooth operation.

We verify the accuracy of the model for proposing by studying No. 2 the 1000 of blast furnace groups of data (shown in Fig. 3) of Liu Gang. Below, we are explained in detail with reference to detailed process to implementation steps：

Step one：T-S fuzzy neural network models are chosen, and combines sliding window model, for the prediction of silicone content.

Step 2：The input that 11 parameters are chosen as model, silicone content conduct are calculated by practical experience and mutual information Output.

Step 3：After by model initialization, with normalized training sample training pattern, the model that will be trained is used for silicon The prediction of content.

Advise that we have chosen input of 11 variables as model by the practical experience of execute-in-place engineer, it Be respectively top pressure, top temperature, gas permeability, coal powder injection, oxygen enrichment percentage, full tower pressure difference, hot-blast pressure, hot blast temperature, hot air flow, Air humidity and previous stove silicone content.Their association relationships with current silicone content are calculated, following result is obtained：

Numbering	Variable	Unit	Mutual information
				1	Top pressure	kPa	0.12
2	Top temperature	℃	0.22
				3	Gas permeability	m3/min·kPa	0.14
4	Coal powder injection	t/h	0.29
				5	Oxygen enrichment percentage	Vol%	0.21
6	Full tower pressure difference	kPa	0.10
				7	Hot-blast pressure	kPa	0.15
8	Hot blast temperature	℃	0.32
				9	Hot air flow	m3/min	0.13
10	Air humidity	Vol%	0.08
				11	Previous stove silicone content	Wt%	0.45

Association relationship is between 0 to 1, and two correlation of variables of bigger expression are stronger.Previous stove silicon contains as can be seen from the table Amount contacts maximum with current silicone content, but the influence of remaining variable can not be ignored.

Method for normalizing described in step 3 is as follows：

The sample size that we set training data is 400, and test data set is 50.With prediction hit rate J and mean square error MSE two indices verify the precision of model prediction：

In actual production process, predicated error can meet requirement less than 0.1.It is proposed that the hit rate of model reach To 90%, mean square error is 0.0023.With precision very high, it is entirely capable of meeting the demand of actual production.

Above-described embodiment is used for illustrating the present invention, rather than limiting the invention, in spirit of the invention and In scope of the claims, any modifications and changes made to the present invention belong to protection scope of the present invention.

Claims (6)

1. a kind of molten iron silicon content Forecasting Methodology based on sliding window T-S fuzzy neural networks, it is characterised in that step is such as Under：

Step one：T-S fuzzy neural network models are chosen, and combines sliding window model, for the prediction of silicone content；

Step 2：The input for choosing 11 parameters as model is calculated by practical experience and mutual information, silicone content as output, Described parameter be respectively top pressure, top temperature, gas permeability, coal powder injection, oxygen enrichment percentage, full tower pressure difference, hot-blast pressure, hot blast temperature, Hot air flow, air humidity and previous stove silicone content；

Step 3：After by model initialization, with normalized training sample training pattern, the model that will be trained is used for silicone content Prediction；

The structure of the T-S fuzzy neural networks described in step one is as follows：

T-S fuzzy neural networks are constituted by four layers, are respectively input layer, obfuscation layer, fuzzy rule computation layer and output layer, its Middle input is fuzzy, and exports and be to determine, this represents that output is the linear combination of input, and T-S fuzzy neural networks are determined Justice is as follows：

$R^i:If{x}_{1}is{A}_1^i,...,x_{k}is{A}_k^i,then{y}_i=p_0^i+p_1^i{x}_1+...+p_k^i{x}_k$

WhereinIt is fuzzy subset, y_iIt is the calculating output of fuzzy rule；

(1) obfuscation layer is that it is defined as follows based on probability density function μ：

${\mathrm{\μ}}_{A_j^i}=e^{-{(x_j-c_j^i)}^2/b_j^i},(i=1,...,n;j=1,...,k)$

X in formula_jIt is input variable,WithIt is center and the width of probability density function, k is the dimension of |input paramete, and n is mould Paste the quantity of subset；

(2) fuzzy rule computation layer is made up of following formula：

${\mathrm{\ω}}^i={\mathrm{\μ}}_{A_j^1}\left(x_1\right)\×{\mathrm{\μ}}_{A_j^2}\left(x_2\right)\×...\×{\mathrm{\μ}}_{A_j^k}\left(x_k\right),(i=1,...,n)$

(3) output layer is calculated by following formula：

$y_i=\frac{\underset{i=1}{\overset{n}{\Σ}}{\mathrm{\ω}}^i{y}_i}{\underset{i=1}{\overset{n}{\Σ}}{\mathrm{\ω}}^i}.$

2. method according to claim 1, it is characterised in that the study of the T-S fuzzy neural networks described in step one is calculated Method is as follows：

(1) error calculation：

$e=\frac{1}{2}{(y_d-y_c)}^2$

Wherein y_dIt is actual value, y_cIt is predicted value, e is both differences；

(2) coefficient amendment：

$p_j^i=p_j^i(k-1)-\mathrm{\α}\frac{\∂e}{\∂p_j^i}$

$\frac{\∂e}{\∂p_j^i}=\frac{(y_d-y_c){\mathrm{\ω}}^i}{\underset{i=1}{\overset{m}{\Σ}}{\mathrm{\ω}}^i\·x_j}$

In formulaIt is the coefficient of T-S fuzzy neural networks, and α is its learning rate；

(3) parameters revision：

$\begin{array}{c}{c}_j^i\left(k\right)=c_j^i(k-1)-\mathrm{\β}\frac{\∂e}{\∂c_j^i}\\ b_j^i\left(k\right)=b_j^i(k-1)-\mathrm{\β}\frac{\∂e}{\∂b_j^i}\end{array}.$

3. method according to claim 1, it is characterised in that the sliding window model principle described in step one is as follows：

Sliding window model be built upon it is a kind of assume on, i.e., current output depends on current input, and input and output it Between mapping ruler can be obtained by historical data, according to this it is assumed that presetting a certain amount of training set sample, then It is continuously updated sample data and gives up earliest data point, with the slip of window, T-S fuzzy neural networks can be constantly updated Its structural parameters simultaneously provides newest predicted value.

4. method according to claim 1, it is characterised in that the selection process of the input variable described in step 2 is as follows：

Mutual information is a kind of method of test variable correlation, is comprised the following steps that described：

$\begin{array}{c}I(X;Y)\\ =\underset{x\&Element;X}{\&Sigma;}\underset{y\&Element;Y}{\&Sigma;}p(x,y){\log }_b\frac{p(x,y)}{p\left(x\right)p\left(y\right)}\\ =\mathrm{\&psi;}\left(k\right)-<\mathrm{\&psi;}(n_x+1)+\mathrm{\&psi;}(n_y+1)>+\mathrm{\&psi;}\left(N\right)\end{array}$

K is the number of neighbour given at the beginning in formula, and ψ is that Digamma functions can be expressed as：

ψ (x)=Γ (x)^-1dΓ(x)/dx

It obeys following iterative relation：

ψ (x+1)=ψ (x)+1/x

Ψ (1)=- C, C=0.5772156...

$<...>=N^{-1}\underset{i=1}{\overset{N}{\&Sigma;}}E\&lsqb;...\left(i\right)\&rsqb;$

In order to obtain n_xAnd n_y, it is necessary to calculate sample z_iAnd z_jThe distance between d_i,j：

d_i,j=| | z_i-z_j||:d_i,j1≤d_i,j2≤d_i,j3…

||z_i-z_j| |=max | | x_i-x_j||,||y_i-y_j||}

As ε (i)=max { ε_x(i),ε_y(i) }, ε (i)/2 are taken as z_iWith the distance of k rank neighbours, n_xI () is to x_iDistance is less than ε The number of the point of (i)/2, n_yI () is to y_iNumber of the distance less than the point of ε (i)/2；

Advised by the practical experience of execute-in-place engineer, have chosen input of 11 variables as model, calculate they with The association relationship of current silicone content, association relationship is between 0 to 1, and two correlation of variables of bigger expression are stronger.

5. method according to claim 1, it is characterised in that the method for normalizing described in step 3 is as follows：

$y=\frac{(y_{\max }-y_{\min })(x-x_{\min })}{x_{\max }-x_{\min }}+y_{\min },(y_{min}=-1,y_{max}=1).$

6. method according to claim 1, it is characterised in that model suitable for ironmaking processes the time-varying of blast furnace, dynamic, Non-linear, strong inertia and multiple dimensioned characteristic.