CN103472866B - The pesticide waste liquid incinerator furnace temperature optimization system of intelligent fuzzy system and method - Google Patents

Abstract

The invention discloses a kind of pesticide waste liquid incinerator furnace temperature optimization system of intelligent fuzzy system.The method, based on fuzzy system, is carried out Fuzzy processing to the training sample Output rusults mapping process through error backward propagation method, is obtained last furnace temperature predicted value and the performance variable value making furnace temperature the best.In the present invention, from DCS database, the input matrix of data as training sample of described variable when producing normal is gathered, through the process of standardization module, as the input of fuzzy system module; The furnace temperature predicted value that fuzzy system obtains and make the performance variable value of furnace temperature the best pass to result display module, simultaneously model modification module, by the sampling time interval of setting, collection site intelligent instrument signal.The present invention finally achieves accurately and in real time controlling of furnace temperature, has carried out effective suppression, avoid and occur that furnace temperature is too low or too high the noise in system.

Description

The pesticide waste liquid incinerator furnace temperature optimization system of intelligent fuzzy system and method

Technical field

The present invention relates to pesticide producing liquid waste incineration field, especially, relate to pesticide waste liquid incinerator furnace temperature optimization system and the method for intelligent fuzzy system.

Background technology

Along with developing rapidly of pesticide industry, the problem of environmental pollution of emission has caused the great attention of national governments and corresponding environmental administration.The qualified discharge of research and solution agricultural chemicals organic liquid waste controls and harmless minimization, not only becomes difficult point and the focus of various countries' scientific research, is also the science proposition of the national active demand being related to social sustainable development simultaneously.

Burning method be process at present agricultural chemicals raffinate and waste residue the most effectively, thoroughly, the most general method of application.In burning process, incinerator furnace temperature must remain on a suitable temperature, and too low furnace temperature is unfavorable for the decomposition of poisonous and harmful element in discarded object; Too high furnace temperature not only increases fuel consumption, increases equipment operating cost, and easily damages inboard wall of burner hearth, shortens equipment life.In addition, excessive temperature may increase the volatile quantity of metal and the generation of nitrogen oxide in discarded object.Special in chloride waste water, suitable furnace temperature more can reduce the corrosion of inwall.But the factor affecting furnace temperature in actual burning process is complicated and changeable, easily there is the phenomenon that furnace temperature is too low or too high.

Artificial neural network in recent years, especially error backward propagation method, achieve good effect in system optimization.Neural network has very strong self-adaptation, self-organization, the ability of self study and the ability of large-scale parallel computing.But in actual applications, neural network also exposes some self intrinsic defects: the initialization of weights is random, is easily absorbed in local minimum; In learning process, the interstitial content of hidden layer and the selection of other parameters can only rule of thumb be selected with experiment; Convergence time is long, poor robustness etc.Secondly, the DCS data that industry spot collects also because noise, manual operation error etc. are with certain uncertain error, so use the general Generalization Ability of model of the artificial neural network that determinacy is strong or not.

Nineteen sixty-five U.S. mathematician L.Zadeh first proposed the concept of fuzzy set.Subsequently fuzzy logic with its problem closer to daily people and the meaning of one's words statement mode, start the classical logic replacing adhering to that all things can represent with binary item.Fuzzy logic so far successful Application among multiple fields of industry, the such as field such as household electrical appliances, Industry Control.2003, Demirci proposed the concept of fuzzy system, by using fuzzy membership matrix and the input matrix new with its distortion structure one, then in local equation, showed that analytic value is as last output using the gravity model appoach in Anti-fuzzy method.For pesticide waste liquid incinerator furnace temperature optimization system and method, consider the noise effect in industrial processes and operate miss, the fuzzy performance of fuzzy logic can be used to reduce error to the impact of precision.

Particle cluster algorithm, i.e. Particle Swarm Optimization, being a kind of a kind of biological intelligence optimizing algorithm seeking global optimum by imitating birds flight data dispose put forward by Kennedy and Eberhart professor, being called for short PSO.This algorithm is influenced each other by interparticle in colony, decreases the risk that searching algorithm is absorbed in locally optimal solution, has good global search performance.Particle cluster algorithm is used to the best parameter group searching for error back propagation fuzzy system, to reach the object of Optimized model.

Summary of the invention

Being difficult to control, easily occur the deficiency that furnace temperature is too low or too high to overcome existing incinerator furnace temperature, the invention provides and a kind ofly realize that furnace temperature accurately controls, the pesticide waste liquid incinerator furnace temperature optimization system avoiding occurring that furnace temperature is too low or too high and method.

The technical solution adopted for the present invention to solve the technical problems is:

The pesticide waste liquid incinerator furnace temperature optimization system of intelligent fuzzy system, comprise incinerator, intelligent instrument, DCS system, data-interface and host computer, described DCS system comprises control station and database; Described field intelligent instrument is connected with DCS system, and described DCS system is connected with host computer, and described host computer comprises:

Standardization module, carries out pre-service for the model training sample will inputted from DCS database, to training sample centralization, namely deducts the mean value of sample, then carries out standardization to it:

Computation of mean values: $\stackrel{\‾}{\mathrm{TX}}=\frac{1}{N}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{\mathrm{TX}}_i---\left(1\right)$

Calculate variance: ${\mathrm{\σ}}_x^2=\frac{1}{N-1}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}({\mathrm{TX}}_i-\stackrel{\‾}{\mathrm{TX}})---\left(2\right)$

Standardization: $X=\frac{\mathrm{TX}-\stackrel{\‾}{\mathrm{TX}}}{{\mathrm{\σ}}_x}---\left(3\right)$

Wherein, TX _ibe i-th training sample, be the production that gathers from DCS database normal time key variables, furnace temperature and make the data of the optimized performance variable of furnace temperature, N is number of training, for the average of training sample, X is the training sample after standardization.σ _xrepresent the standard deviation of training sample, σ ² _xrepresent the variance of training sample.

Fuzzy system module, to the training sample X passed from data preprocessing module after the standardization of coming, carries out obfuscation.If have c in fuzzy system ^*individual fuzzy group, the center of fuzzy group k, j is respectively v _k, v _j, then the training sample X after i-th standardization _ifor the degree of membership μ of fuzzy group k _ikfor:

${\mathrm{\μ}}_{\mathrm{ik}}={\left(\underset{j=1}{\overset{c^{*}}{\mathrm{\Σ}}}{\left(\frac{\left|\right|X_i-v_k\left|\right|}{\left|\right|X_i-v_j\left|\right|}\right)}^{\frac{2}{n-1}}\right)}^{-1}---\left(4\right)$

In formula, n is the partitioned matrix index needed in fuzzy classification process, is usually taken as 2, || || be norm expression formula.

Use and be subordinate to angle value or its distortion to obtain new input matrix above, for fuzzy group k, its input matrix is deformed into:

${\mathrm{\Φ}}_{\mathrm{ik}}(X_i,{\mathrm{\μ}}_{\mathrm{ik}})=\left[\begin{array}{ccc}1& \mathrm{func}\left({\mathrm{\μ}}_{\mathrm{ik}}\right)& X_i\end{array}\right]---\left(5\right)$

Wherein func (μ _ik) for being subordinate to angle value μ _ikwarping function, generally get exp (μ _ik) etc., Φ _ik(X _i, μ _ik) represent i-th input variable X _iand the degree of membership μ of fuzzy group k _ikcorresponding new input matrix.

Error backward propagation method, as the local equation of fuzzy system, if the prediction of a kth error back propagation fuzzy system output layer exports is be input as net, in the hidden layer that layer is adjacent therewith, the output of arbitrary neuron l is s _l, then have:

$\mathrm{net}=\underset{l}{\mathrm{\Σ}}{w}_{\mathrm{lk}}\×s_l---\left(6\right)$

${\hat{y}}_{\mathrm{ik}}=f\left(\mathrm{net}\right)---\left(7\right)$

In formula, w _lkbe the connection weight between hidden neuron l and output layer neuron, f () is neuronic output function, is usually taken as Sigmoid function, is expressed as:

$f\left(\mathrm{net}\right)=1/(1+e^{-(\mathrm{net}+h)/{\mathrm{\θ}}_0})---\left(8\right)$

In formula, h is the neuronic threshold value of output layer, θ ₀for steepness parameter, in order to regulate the steepness of Sigmoid function; Finally, the output of last fuzzy system is obtained by the gravity model appoach in Anti-fuzzy method:

${\hat{y}}_i=\frac{{\mathrm{\Σ}}_{k=1}^{c^{*}}{\mathrm{\μ}}_{\mathrm{ik}}{\hat{y}}_{\mathrm{ik}}}{{\mathrm{\Σ}}_{k=1}^{c^{*}}{\mathrm{\μ}}_{\mathrm{ik}}}---\left(9\right)$

the prediction being fuzzy system exports.

Intelligent optimization module, for adopting particle cluster algorithm to the w of fuzzy system medial error reverse transmittance nerve network local equation _lkbe optimized, specific implementation step is as follows:

1. determine that the Optimal Parameters of population is the w of error backward propagation method local equation _lk, population individual amount popsize, largest loop optimizing number of times iter _max, a p particle initial position r _p, initial velocity v _p, local optimum Lbest _pand the global optimum Gbest of whole population.

2. set optimization object function, be converted into fitness, each On Local Fuzzy equation is evaluated; Calculate fitness function by corresponding error function, and think that the large particle fitness of error is little, the fitness function of particle p is expressed as:

f _p=1/(E _p+1) （10）

In formula, E _pbe the error function of fuzzy system, be expressed as:

$E_p=\sqrt{\frac{1}{N}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{({\hat{y}}_i-O_i)}^2}---\left(11\right)$

In formula, the prediction output of fuzzy system, O _ifor the target of fuzzy system exports, N is number of training;

3. according to following formula, circulation upgrades speed and the position of each particle,

v _p(iter+1)=ω×v _p(iter)+m ₁a ₁(Lbest _p-r _p(iter))+m ₂a ₂(Gbest-r _p(iter))

（12）

r _p(iter+1)=r _p(iter)+v _p(iter+1) （13）

In formula, v _prepresent the speed of more new particle p, r _prepresent the position of more new particle p, Lbest _prepresent the individual optimal value of more new particle p, Gbest is the training sample X after corresponding to i-th standardization _ifurnace temperature predicted value and make the performance variable value of furnace temperature the best, iter represents cycle index, and ω is the inertia weight in particle cluster algorithm, m ₁, m ₂corresponding accelerator coefficient, a ₁, a ₂it is the random number between [0,1];

4. for particle p, if new fitness is greater than original individual optimal value, the more individual optimal value of new particle:

Lbest _p=f _p（14）

If the 5. individual optimal value Lbest of particle p _pbe greater than original population global optimum Gbest, upgrade original population global optimum Gbest:

Gbest=Lbest _p（15）

6. judge whether to meet performance requirement, if so, terminate optimizing, obtain the local equation parameter of the fuzzy system that a group is optimized; Otherwise return step 3., continue iteration optimizing, until reach maximum iteration time iter _max.

Gbest is the training sample X after corresponding to i-th standardization _ifurnace temperature predicted value and make the performance variable value of furnace temperature the best.

As preferred a kind of scheme: described host computer also comprises: model modification module, for the sampling time interval by setting, collection site intelligent instrument signal, the actual measurement furnace temperature function predicted value obtained is compared, if relative error be greater than 10% or furnace temperature exceed produce normal bound scope, then the new data of furnace temperature the best that makes when producing normal in DCS database is added training sample data, upgrade soft-sensing model.

Further, described host computer also comprises: result display module, for by furnace temperature predicted value and make the performance variable value of furnace temperature the best pass to DCS system, shows at the control station of DCS, and is delivered to operator station by DCS system and fieldbus and shows; Meanwhile, obtained makes the performance variable value of furnace temperature the best as new performance variable setting value by DCS system, automatically performs the operation of furnace temperature optimization.

Signal acquisition module, for the time interval of each sampling according to setting, image data from database.

Further again, described key variables comprise enter incinerator waste liquid flow, enter the air mass flow of incinerator and enter the fuel flow rate of incinerator; Described performance variable comprises the air mass flow entering incinerator and the fuel flow rate entering incinerator.

The furnace temperature optimization method of the pesticide waste liquid incinerator furnace temperature optimization system realization of intelligent fuzzy system, described furnace temperature optimization method specific implementation step is as follows:

1), determine key variables used, from DCS database, gather the input matrix of data as training sample TX of described variable when producing normal, gather corresponding furnace temperature and make the optimized performance variable data of furnace temperature as output matrix O;

2), by the model training sample inputted from DCS database carry out pre-service, to training sample centralization, namely deduct the mean value of sample, then carry out standardization to it, make its average be 0, variance is 1.This process adopts following formula process:

2.1) computation of mean values: $\stackrel{\‾}{\mathrm{TX}}=\frac{1}{N}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{\mathrm{TX}}_i---\left(1\right)$

2.2) variance is calculated: ${\mathrm{\σ}}_x^2=\frac{1}{N-1}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}({\mathrm{TX}}_i-\stackrel{\‾}{\mathrm{TX}})---\left(2\right)$

2.3) standardization: $X=\frac{\mathrm{TX}-\stackrel{\‾}{\mathrm{TX}}}{{\mathrm{\σ}}_x}---\left(3\right)$

Wherein, TX _ibe i-th training sample, be the production that gathers from DCS database normal time key variables, furnace temperature and make the data of the optimized performance variable of furnace temperature, N is number of training, for the average of training sample, X is the training sample after standardization.σ _xrepresent the standard deviation of training sample, σ ² _xrepresent the variance of training sample.

3), to passing the training sample of coming from data preprocessing module, obfuscation is carried out.If have c in fuzzy system ^*individual fuzzy group, the center of fuzzy group k, j is respectively v _k, v _j, then the training sample X after i-th standardization _ifor the degree of membership μ of fuzzy group k _ikfor:

${\mathrm{\μ}}_{\mathrm{ik}}={\left(\underset{j=1}{\overset{c^{*}}{\mathrm{\Σ}}}{\left(\frac{\left|\right|X_i-v_k\left|\right|}{\left|\right|X_i-v_j\left|\right|}\right)}^{\frac{2}{n-1}}\right)}^{-1}---\left(4\right)$

In formula, n is the partitioned matrix index needed in fuzzy classification process, is usually taken as 2, || || be norm expression formula.

Use and be subordinate to angle value or its distortion to obtain new input matrix above, for fuzzy group k, its input matrix is deformed into:

${\mathrm{\Φ}}_{\mathrm{ik}}(X_i,{\mathrm{\μ}}_{\mathrm{ik}})=\left[\begin{array}{ccc}1& \mathrm{func}\left({\mathrm{\μ}}_{\mathrm{ik}}\right)& X_i\end{array}\right]---\left(5\right)$

Wherein func (μ _ik) for being subordinate to angle value μ _ikwarping function, generally get exp (μ _ik) etc., Φ _ik(X _i, μ _ik) represent i-th input variable X _iand the degree of membership μ of fuzzy group k _ikcorresponding new input matrix.

Error backward propagation method, as the local equation of fuzzy system, if the prediction of a kth error back propagation fuzzy system output layer exports is be input as net, in the hidden layer that layer is adjacent therewith, the output of arbitrary neuron l is s _l, then have:

$\mathrm{net}=\underset{l}{\mathrm{\Σ}}{w}_{\mathrm{lk}}\×s_l---\left(6\right)$

${\hat{y}}_{\mathrm{ik}}=f\left(\mathrm{net}\right)---\left(7\right)$

In formula, w _lkbe the connection weight between hidden neuron l and output layer neuron, f () is neuronic output function, is usually taken as Sigmoid function, is expressed as:

$f\left(\mathrm{net}\right)=1/(1+e^{-(\mathrm{net}+h)/{\mathrm{\θ}}_0})---\left(8\right)$

In formula, h is the neuronic threshold value of output layer, θ ₀for steepness parameter, in order to regulate the steepness of Sigmoid function; Finally, the output of last fuzzy system is obtained by the gravity model appoach in Anti-fuzzy method:

${\hat{y}}_i=\frac{{\mathrm{\Σ}}_{k=1}^{c^{*}}{\mathrm{\μ}}_{\mathrm{ik}}{\hat{y}}_{\mathrm{ik}}}{{\mathrm{\Σ}}_{k=1}^{c^{*}}{\mathrm{\μ}}_{\mathrm{ik}}}---\left(9\right)$

the prediction being fuzzy system exports.

4), adopt particle cluster algorithm to the w of fuzzy system medial error reverse transmittance nerve network local equation _lkbe optimized, specific implementation step is as follows:

1. determine that the Optimal Parameters of population is the w of error backward propagation method local equation _lk, population individual amount popsize, largest loop optimizing number of times iter _max, a p particle initial position r _p, initial velocity v _p, local optimum Lbest _pand the global optimum Gbest of whole population.

2. set optimization object function, be converted into fitness, each On Local Fuzzy equation is evaluated; Calculate fitness function by corresponding error function, and think that the large particle fitness of error is little, the fitness function of particle p is expressed as:

f _p=1/(E _p+1) （10）

In formula, E _pbe the error function of fuzzy system, be expressed as:

$E_p=\sqrt{\frac{1}{N}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{({\hat{y}}_i-O_i)}^2}---\left(11\right)$

In formula, the prediction output of fuzzy system, O _ifor the target of fuzzy system exports, N is number of training;

3. according to following formula, circulation upgrades speed and the position of each particle,

v _p(iter+1)=ω×v _p(iter)+m ₁a ₁(Lbest _p-r _p(iter))+m ₂a ₂(Gbest-r _p(iter))

（12）

r _p(iter+1)=r _p(iter)+v _p(iter+1) （13）

In formula, v _prepresent the speed of more new particle p, r _prepresent the position of more new particle p, Lbest _prepresent the individual optimal value of more new particle p, Gbest is the training sample X after corresponding to i-th standardization _ifurnace temperature predicted value and make the performance variable value of furnace temperature the best, iter represents cycle index, and ω is the inertia weight in particle cluster algorithm, m ₁, m ₂corresponding accelerator coefficient, a ₁, a ₂it is the random number between [0,1];

4. for particle p, if new fitness is greater than original individual optimal value, the more individual optimal value of new particle:

Lbest _p=f _p（14）

If the 5. individual optimal value Lbest of particle p _pbe greater than original population global optimum Gbest, upgrade original population global optimum Gbest:

Gbest=Lbest _p（15）

6. judge whether to meet performance requirement, if so, terminate optimizing, obtain the local equation parameter of the fuzzy system that a group is optimized; Otherwise return step 3., continue iteration optimizing, until reach maximum iteration time iter _max.

Gbest is the training sample X after corresponding to i-th standardization _ifurnace temperature predicted value and make the performance variable value of furnace temperature the best.

As preferred a kind of scheme: described method also comprises: 5), by the sampling time interval set, collection site intelligent instrument signal, the actual measurement furnace temperature function predicted value obtained is compared, if relative error be greater than 10% or furnace temperature exceed produce normal bound scope, then the new data of furnace temperature the best that makes when producing normal in DCS database is added training sample data, upgrade soft-sensing model.

Further, the furnace temperature predicted value obtained in described step 4) and make the performance variable value of furnace temperature the best, passes to DCS system by result, shows at the control station of DCS, and is delivered to operator station by DCS system and fieldbus and shows; Meanwhile, obtained makes the performance variable value of furnace temperature the best as new performance variable setting value by DCS system, automatically performs the operation of furnace temperature optimization.

Further again, described key variables comprise enter incinerator waste liquid flow, enter the air mass flow of incinerator and enter the fuel flow rate of incinerator; Described performance variable comprises the air mass flow entering incinerator and the fuel flow rate entering incinerator.

Technical conceive of the present invention is: pesticide waste liquid incinerator furnace temperature optimization system and the method for having invented intelligent fuzzy system, search out the performance variable value making furnace temperature the best.

Beneficial effect of the present invention is mainly manifested in: the online soft sensor model 1, establishing quantitative relationship between system core variable and furnace temperature; 2, the operating conditions making furnace temperature the best is found rapidly.

Accompanying drawing explanation

Fig. 1 is the hardware structure diagram of system proposed by the invention;

Fig. 2 is the functional structure chart of host computer proposed by the invention.

Embodiment

Below in conjunction with accompanying drawing, the invention will be further described.The embodiment of the present invention is used for explaining and the present invention is described, instead of limits the invention, and in the protection domain of spirit of the present invention and claim, any amendment make the present invention and change, all fall into protection scope of the present invention.

Embodiment 1

With reference to Fig. 1, Fig. 2, the pesticide waste liquid incinerator furnace temperature optimization system of intelligent fuzzy system, comprise field intelligent instrument 2, DCS system and the host computer 6 be connected with incinerator object 1, described DCS system comprises data-interface 3, control station 4 and database 5, described field intelligent instrument 2 is connected with data-interface 3, described data-interface is connected with control station 4, database 5 and host computer 6, and described host computer 6 comprises:

Standardization module 7, carries out pre-service for the model training sample will inputted from DCS database, to training sample centralization, namely deducts the mean value of sample, then carries out standardization to it:

Computation of mean values: $\stackrel{\‾}{\mathrm{TX}}=\frac{1}{N}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{\mathrm{TX}}_i---\left(1\right)$

Calculate variance: ${\mathrm{\σ}}_x^2=\frac{1}{N-1}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}({\mathrm{TX}}_i-\stackrel{\‾}{\mathrm{TX}})---\left(2\right)$

Standardization: $X=\frac{\mathrm{TX}-\stackrel{\‾}{\mathrm{TX}}}{{\mathrm{\σ}}_x}---\left(3\right)$

Wherein, TX _ibe i-th training sample, be the production that gathers from DCS database normal time key variables, furnace temperature and make the data of the optimized performance variable of furnace temperature, N is number of training, for the average of training sample, X is the training sample after standardization.σ _xrepresent the standard deviation of training sample, σ ² _xrepresent the variance of training sample.

Fuzzy system module 8, to the training sample X passed from data preprocessing module after the standardization of coming, carries out obfuscation.If have c in fuzzy system ^*individual fuzzy group, the center of fuzzy group k, j is respectively v _k, v _j, then the training sample X after i-th standardization _ifor the degree of membership μ of fuzzy group k _ikfor:

${\mathrm{\μ}}_{\mathrm{ik}}={\left(\underset{j=1}{\overset{c^{*}}{\mathrm{\Σ}}}{\left(\frac{\left|\right|X_i-v_k\left|\right|}{\left|\right|X_i-v_j\left|\right|}\right)}^{\frac{2}{n-1}}\right)}^{-1}---\left(4\right)$

In formula, n is the partitioned matrix index needed in fuzzy classification process, is usually taken as 2, || || be norm expression formula.

Use and be subordinate to angle value or its distortion to obtain new input matrix above, for fuzzy group k, its input matrix is deformed into:

${\mathrm{\Φ}}_{\mathrm{ik}}(X_i,{\mathrm{\μ}}_{\mathrm{ik}})=\left[\begin{array}{ccc}1& \mathrm{func}\left({\mathrm{\μ}}_{\mathrm{ik}}\right)& X_i\end{array}\right]---\left(5\right)$

Wherein func (μ _ik) for being subordinate to angle value μ _ikwarping function, generally get exp (μ _ik) etc., Φ _ik(X _i, μ _ik) represent i-th input variable X _iand the degree of membership μ of fuzzy group k _ikcorresponding new input matrix.

Error backward propagation method, as the local equation of fuzzy system, if the prediction of a kth error back propagation fuzzy system output layer exports is be input as net, in the hidden layer that layer is adjacent therewith, the output of arbitrary neuron l is s _l, then have:

$\mathrm{net}=\underset{l}{\mathrm{\Σ}}{w}_{\mathrm{lk}}\×s_l---\left(6\right)$

${\hat{y}}_{\mathrm{ik}}=f\left(\mathrm{net}\right)---\left(7\right)$

In formula, w _lkbe the connection weight between hidden neuron l and output layer neuron, f () is neuronic output function, is usually taken as Sigmoid function, is expressed as:

$f\left(\mathrm{net}\right)=1/(1+e^{-(\mathrm{net}+h)/{\mathrm{\θ}}_0})---\left(8\right)$

In formula, h is the neuronic threshold value of output layer, θ ₀for steepness parameter, in order to regulate the steepness of Sigmoid function; Finally, the output of last fuzzy system is obtained by the gravity model appoach in Anti-fuzzy method:

${\hat{y}}_i=\frac{{\mathrm{\Σ}}_{k=1}^{c^{*}}{\mathrm{\μ}}_{\mathrm{ik}}{\hat{y}}_{\mathrm{ik}}}{{\mathrm{\Σ}}_{k=1}^{c^{*}}{\mathrm{\μ}}_{\mathrm{ik}}}---\left(9\right)$

the prediction being fuzzy system exports.

Intelligent optimization module 9, for adopting particle cluster algorithm to the w of fuzzy system medial error reverse transmittance nerve network local equation _lkbe optimized, specific implementation step is as follows:

1. determine that the Optimal Parameters of population is the w of error backward propagation method local equation _lk, population individual amount popsize, largest loop optimizing number of times iter _max, a p particle initial position r _p, initial velocity v _p, local optimum Lbest _pand the global optimum Gbest of whole population.

2. set optimization object function, be converted into fitness, each On Local Fuzzy equation is evaluated; Calculate fitness function by corresponding error function, and think that the large particle fitness of error is little, the fitness function of particle p is expressed as:

f _p=1/(E _p+1) （10）

In formula, E _pbe the error function of fuzzy system, be expressed as:

$E_p=\sqrt{\frac{1}{N}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{({\hat{y}}_i-O_i)}^2}---\left(11\right)$

In formula, the prediction output of fuzzy system, O _ifor the target of fuzzy system exports, N is number of training;

3. according to following formula, circulation upgrades speed and the position of each particle,

v _p(iter+1)=ω×v _p(iter)+m ₁a ₁(Lbest _p-r _p(iter))+m ₂a ₂(Gbest-r _p(iter))

（12）

r _p(iter+1)=r _p(iter)+v _p(iter+1) （13）

In formula, v _prepresent the speed of more new particle p, r _prepresent the position of more new particle p, Lbest _prepresent the individual optimal value of more new particle p, Gbest is the training sample X after corresponding to i-th standardization _ifurnace temperature predicted value and make the performance variable value of furnace temperature the best, iter represents cycle index, and ω is the inertia weight in particle cluster algorithm, m ₁, m ₂corresponding accelerator coefficient, a ₁, a ₂it is the random number between [0,1];

4. for particle p, if new fitness is greater than original individual optimal value, the more individual optimal value of new particle:

Lbest _p=f _p（14）

If the 5. individual optimal value Lbest of particle p _pbe greater than original population global optimum Gbest, upgrade original population global optimum Gbest:

Gbest=Lbest _p（15）

6. judge whether to meet performance requirement, if so, terminate optimizing, obtain the local equation parameter of the fuzzy system that a group is optimized; Otherwise return step 3., continue iteration optimizing, until reach maximum iteration time iter _max.

Gbest is the training sample X after corresponding to i-th standardization _ifurnace temperature predicted value and make the performance variable value of furnace temperature the best.

Described host computer 6 also comprises: signal acquisition module 11, for the time interval of each sampling according to setting, and image data from database.

Described host computer 6 also comprises: model modification module 12, by the sampling time interval of setting, collection site intelligent instrument signal, the actual measurement furnace temperature function predicted value obtained is compared, if relative error be greater than 10% or furnace temperature exceed produce normal bound scope, then the new data of furnace temperature the best that makes when producing normal in DCS database is added training sample data, upgrade soft-sensing model.

Described key variables comprise enter incinerator waste liquid flow, enter the air mass flow of incinerator and enter the fuel flow rate of incinerator; Described performance variable comprises the air mass flow entering incinerator and the fuel flow rate entering incinerator.

Described system also comprises DCS system, and described DCS system is made up of data-interface 3, control station 4, database 5; Intelligent instrument 2, DCS system, host computer 6 are connected successively by fieldbus; Host computer 6 also comprises result display module 10, for forecast result is passed to DCS system, and in the control station procedure for displaying state of DCS, by DCS system and fieldbus, process state information is delivered to operator station and shows; Meanwhile, obtained makes the performance variable value of furnace temperature the best as new performance variable setting value by DCS system, automatically performs the operation of furnace temperature optimization.

When liquid waste incineration process has been furnished with DCS system, the detection of sample real-time dynamic data, the real-time of memory DCS system and historical data base, obtained furnace temperature predicted value and the function of the performance variable value of furnace temperature the best mainly completed on host computer.

When liquid waste incineration process is not equipped with DCS system, data-carrier store is adopted to carry out the data storage function of the real-time of alternative DCS system and historical data base, and do not rely on the independently complete SOC (system on a chip) of of DCS system by what obtain furnace temperature predicted value and make the function system of the performance variable value of furnace temperature the best manufacture to comprise I/O element, data-carrier store, program storage, arithmetical unit, several large component of display module, when whether being equipped with DCS regardless of burning process, independently can both use, more be of value to and promoting the use of.

Embodiment 2

With reference to Fig. 1, Fig. 2, the pesticide waste liquid incinerator furnace temperature optimization method of intelligent fuzzy system, described method specific implementation step is as follows:

1), determine key variables used, from DCS database, gather the input matrix of data as training sample TX of described variable when producing normal, gather corresponding furnace temperature and make the optimized performance variable data of furnace temperature as output matrix O;

2), by the model training sample inputted from DCS database carry out pre-service, to training sample centralization, namely deduct the mean value of sample, then carry out standardization to it, make its average be 0, variance is 1.This process adopts following formula process:

2.1) computation of mean values: $\stackrel{\‾}{\mathrm{TX}}=\frac{1}{N}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{\mathrm{TX}}_i---\left(1\right)$

2.2) variance is calculated: ${\mathrm{\σ}}_x^2=\frac{1}{N-1}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}({\mathrm{TX}}_i-\stackrel{\‾}{\mathrm{TX}})---\left(2\right)$

2.3) standardization: $X=\frac{\mathrm{TX}-\stackrel{\‾}{\mathrm{TX}}}{{\mathrm{\σ}}_x}---\left(3\right)$

Wherein, TX _ibe i-th training sample, be the production that gathers from DCS database normal time key variables, furnace temperature and make the data of the optimized performance variable of furnace temperature, N is number of training, for the average of training sample, X is the training sample after standardization.σ _xrepresent the standard deviation of training sample, σ ² _xrepresent the variance of training sample.

3), to passing the training sample of coming from data preprocessing module, obfuscation is carried out.If have c in fuzzy system ^*individual fuzzy group, the center of fuzzy group k, j is respectively v _k, v _j, then the training sample X after i-th standardization _ifor the degree of membership μ of fuzzy group k _ikfor:

${\mathrm{\μ}}_{\mathrm{ik}}={\left(\underset{j=1}{\overset{c^{*}}{\mathrm{\Σ}}}{\left(\frac{\left|\right|X_i-v_k\left|\right|}{\left|\right|X_i-v_j\left|\right|}\right)}^{\frac{2}{n-1}}\right)}^{-1}---\left(4\right)$

In formula, n is the partitioned matrix index needed in fuzzy classification process, is usually taken as 2, || || be norm expression formula.

Use and be subordinate to angle value or its distortion to obtain new input matrix above, for fuzzy group k, its input matrix is deformed into:

${\mathrm{\Φ}}_{\mathrm{ik}}(X_i,{\mathrm{\μ}}_{\mathrm{ik}})=\left[\begin{array}{ccc}1& \mathrm{func}\left({\mathrm{\μ}}_{\mathrm{ik}}\right)& X_i\end{array}\right]---\left(5\right)$

Wherein func (μ _ik) for being subordinate to angle value μ _ikwarping function, generally get exp (μ _ik) etc., Φ _ik(X _i, μ _ik) represent i-th input variable X _iand the degree of membership μ of fuzzy group k _ikcorresponding new input matrix.

Error backward propagation method, as the local equation of fuzzy system, if the prediction of a kth error back propagation fuzzy system output layer exports is be input as net, in the hidden layer that layer is adjacent therewith, the output of arbitrary neuron l is s _l, then have:

$\mathrm{net}=\underset{l}{\mathrm{\Σ}}{w}_{\mathrm{lk}}\×s_l---\left(6\right)$

${\hat{y}}_{\mathrm{ik}}=f\left(\mathrm{net}\right)---\left(7\right)$

In formula, w _lkbe the connection weight between hidden neuron l and output layer neuron, f () is neuronic output function, is usually taken as Sigmoid function, is expressed as:

$f\left(\mathrm{net}\right)=1/(1+e^{-(\mathrm{net}+h)/{\mathrm{\θ}}_0})---\left(8\right)$

In formula, h is the neuronic threshold value of output layer, θ ₀for steepness parameter, in order to regulate the steepness of Sigmoid function; Finally, the output of last fuzzy system is obtained by the gravity model appoach in Anti-fuzzy method:

${\hat{y}}_i=\frac{{\mathrm{\Σ}}_{k=1}^{c^{*}}{\mathrm{\μ}}_{\mathrm{ik}}{\hat{y}}_{\mathrm{ik}}}{{\mathrm{\Σ}}_{k=1}^{c^{*}}{\mathrm{\μ}}_{\mathrm{ik}}}---\left(9\right)$

the prediction being fuzzy system exports.

4), for adopting particle cluster algorithm to the w of fuzzy system medial error reverse transmittance nerve network local equation _lkbe optimized, specific implementation step is as follows:

1. determine that the Optimal Parameters of population is the w of error backward propagation method local equation _lk, population individual amount popsize, largest loop optimizing number of times iter _max, a p particle initial position r _p, initial velocity v _p, local optimum Lbest _pand the global optimum Gbest of whole population.

2. set optimization object function, be converted into fitness, each On Local Fuzzy equation is evaluated; Calculate fitness function by corresponding error function, and think that the large particle fitness of error is little, the fitness function of particle p is expressed as:

f _p=1/(E _p+1) （10）

In formula, E _pbe the error function of fuzzy system, be expressed as:

$E_p=\sqrt{\frac{1}{N}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{({\hat{y}}_i-O_i)}^2}---\left(11\right)$

In formula, the prediction output of fuzzy system, O _ifor the target of fuzzy system exports, N is number of training;

3. according to following formula, circulation upgrades speed and the position of each particle,

v _p(iter+1)=ω×v _p(iter)+m ₁a ₁(Lbest _p-r _p(iter))+m ₂a ₂(Gbest-r _p(iter))

（12）

r _p(iter+1)=r _p(iter)+v _p(iter+1) （13）

In formula, v _prepresent the speed of more new particle p, r _prepresent the position of more new particle p, Lbest _prepresent the individual optimal value of more new particle p, Gbest is the training sample X after corresponding to i-th standardization _ifurnace temperature predicted value and make the performance variable value of furnace temperature the best, iter represents cycle index, and ω is the inertia weight in particle cluster algorithm, m ₁, m ₂corresponding accelerator coefficient, a ₁, a ₂it is the random number between [0,1];

4. for particle p, if new fitness is greater than original individual optimal value, the more individual optimal value of new particle:

Lbest _p=f _p（14）

If the 5. individual optimal value Lbest of particle p _pbe greater than original population global optimum Gbest, upgrade original population global optimum Gbest:

Gbest=Lbest _p（15）

6. judge whether to meet performance requirement, if so, terminate optimizing, obtain the local equation parameter of the fuzzy system that a group is optimized; Otherwise return step 3., continue iteration optimizing, until reach maximum iteration time iter _max.

Gbest is the training sample X after corresponding to i-th standardization _ifurnace temperature predicted value and make the performance variable value of furnace temperature the best.

Described method also comprises: 5), by the sampling time interval set, collection site intelligent instrument signal, the actual measurement furnace temperature obtained is compared with predicts value, if relative error be greater than 10% or furnace temperature exceed produce normal bound scope, then the new data of furnace temperature the best that makes when producing normal in DCS database is added training sample data, upgrade soft-sensing model.

Optimum Operation variate-value is calculated in described step 4), by the furnace temperature predicted value obtained and make the performance variable value of furnace temperature the best pass to DCS system, show at the control station of DCS, and be delivered to operator station by DCS system and fieldbus and show; Meanwhile, obtained makes the performance variable value of furnace temperature the best as new performance variable setting value by DCS system, automatically performs the operation of furnace temperature optimization.

Described key variables comprise enter incinerator waste liquid flow, enter the air mass flow of incinerator and enter the fuel flow rate of incinerator; Described performance variable comprises the air mass flow entering incinerator and the fuel flow rate entering incinerator.

Claims (2)

1. a pesticide waste liquid incinerator furnace temperature optimization system for intelligent fuzzy system, comprise incinerator, field intelligent instrument, DCS system, data-interface and host computer, described DCS system comprises control station and database; Described field intelligent instrument is connected with DCS system, and described DCS system is connected with host computer, it is characterized in that: described host computer comprises:

Standardization module, carries out pre-service for the model training sample will inputted from DCS database, to training sample centralization, namely deducts the mean value of sample, then carries out standardization to it:

Computation of mean values:

Calculate variance:

Standardization:

Wherein, TX _ibe i-th training sample, be the production that gathers from DCS database normal time key variables, furnace temperature and make the data of the optimized performance variable of furnace temperature, N is number of training, for the average of training sample, X is the training sample after standardization; σ _xrepresent the standard deviation of training sample, σ ² _xrepresent the variance of training sample;

Fuzzy system module, to the training sample X passed from standardization module after the standardization of coming, carries out obfuscation; If have c in fuzzy system ^*individual fuzzy group, the center of fuzzy group k, j is respectively v _k, v _j, then the training sample X after i-th standardization _ifor the degree of membership μ of fuzzy group k _ikfor:

In formula, n is the partitioned matrix index needed in fuzzy classification process, is usually taken as 2, || || be norm expression formula;

Use and be subordinate to angle value or its distortion to obtain new input matrix above, for fuzzy group k, its input matrix is deformed into:

Φ _ik(X _i,μ _ik)＝[1 func(μ _ik) X _i] (5)

Wherein func (μ _ik) for being subordinate to angle value μ _ikwarping function, generally get exp (μ _ik), Φ _ik(X _i, μ _ik) represent i-th input variable X _iand the degree of membership μ of fuzzy group k _ikcorresponding new input matrix;

Error backward propagation method, as the local equation of fuzzy system, if the prediction of a kth error back propagation fuzzy system output layer exports is be input as net, in the hidden layer that layer is adjacent therewith, the output of arbitrary neuron l is s _l, then have:

In formula, w _lkbe the connection weight between hidden neuron l and output layer neuron, f () is neuronic output function, is usually taken as Sigmoid function, is expressed as:

In formula, h is the neuronic threshold value of output layer, θ ₀for steepness parameter, in order to regulate the steepness of Sigmoid function; Finally, the output of last fuzzy system is obtained by the gravity model appoach in Anti-fuzzy method:

the prediction being fuzzy system exports;

Intelligent optimization module, for adopting particle cluster algorithm to the w of fuzzy system medial error reverse transmittance nerve network local equation _lkbe optimized, specific implementation step is as follows:

1. determine that the Optimal Parameters of population is the w of error backward propagation method local equation _lk, population individual amount popsize, largest loop optimizing number of times iter _max, a p particle initial position r _p, initial velocity v _p, local optimum Lbest _pand the global optimum Gbest of whole population;

2. set optimization object function, be converted into fitness, each On Local Fuzzy system is evaluated; Calculate fitness function by corresponding error function, and think that the large particle fitness of error is little, the fitness function of particle p is expressed as:

f _p＝1/(E _p+1) (10)

In formula, E _pbe the error function of fuzzy system, be expressed as:

In formula, the prediction output of fuzzy system, O _ifor the target of fuzzy system exports, N is number of training;

3. according to following formula, circulation upgrades speed and the position of each particle,

v _p(iter+1)＝ω×v _p(iter)+m ₁a ₁(Lbest _p-r _p(iter))+m ₂a ₂(Gbest-r _p(iter))

(12)

r _p(iter+1)＝r _p(iter)+v _p(iter+1) (13)

In formula, v _prepresent the speed of more new particle p, r _prepresent the position of more new particle p, Lbest _prepresent the individual optimal value of more new particle p, Gbest is the training sample X after corresponding to i-th standardization _ifurnace temperature predicted value and make the performance variable value of furnace temperature the best, iter represents cycle index, and ω is the inertia weight in particle cluster algorithm, m ₁, m ₂corresponding accelerator coefficient, a ₁, a ₂it is the random number between [0,1];

4. for particle p, if new fitness is greater than original individual optimal value, the more individual optimal value of new particle:

Lbest _p＝f _p(14)

If the 5. individual optimal value Lbest of particle p _pbe greater than original population global optimum Gbest, upgrade original population global optimum Gbest:

Gbest＝Lbest _p(15)

6. judge whether to meet performance requirement, if so, terminate optimizing, obtain the local equation parameter of the fuzzy system that a group is optimized; Otherwise return step 3., continue iteration optimizing, until reach largest loop optimizing number of times iter _max;

Gbest is the training sample X after corresponding to i-th standardization _ifurnace temperature predicted value and make the performance variable value of furnace temperature the best;

Described host computer also comprises:

Model modification module, for the sampling time interval by setting, collection site intelligent instrument signal, the actual measurement furnace temperature obtained is compared with predicts value, if relative error be greater than 10% or furnace temperature exceed produce normal bound scope, then the new data of furnace temperature the best that makes when producing normal in DCS database is added training sample data, upgrade soft-sensing model;

Result display module, for by the furnace temperature obtained predicted value and make the performance variable value of furnace temperature the best pass to DCS system, shows at the control station of DCS, and is delivered to operator station by DCS system and fieldbus and shows; Meanwhile, obtained makes the performance variable value of furnace temperature the best as new performance variable setting value by DCS system, automatically performs the operation of furnace temperature optimization;

Signal acquisition module, for the time interval of each sampling according to setting, image data from database;

Described key variables comprise enter incinerator waste liquid flow, enter the air mass flow of incinerator and enter the fuel flow rate of incinerator; Described performance variable comprises the air mass flow entering incinerator and the fuel flow rate entering incinerator.

2. a pesticide waste liquid incinerator furnace temperature optimization method for intelligent fuzzy system, is characterized in that: described furnace temperature optimization method comprises the following steps:

1), determine key variables used, from DCS database, gather the input matrix of data as training sample TX of described variable when producing normal, gather corresponding furnace temperature and make the optimized performance variable data of furnace temperature as output matrix O;

2), by the model training sample inputted from DCS database carry out pre-service, to training sample centralization, namely deduct the mean value of sample, then carry out standardization to it, make its average be 0, variance is 1; This process adopts following formula process:

2.1) computation of mean values:

2.2) variance is calculated:

2.3) standardization:

Wherein, TX _ibe i-th training sample, be the production that gathers from DCS database normal time key variables, furnace temperature and make the data of the optimized performance variable of furnace temperature, N is number of training, for the average of training sample, X is the training sample after standardization; σ _xrepresent the standard deviation of training sample, σ ² _xrepresent the variance of training sample;

3), to passing the training sample of coming from standardization module, obfuscation is carried out; If have c in fuzzy system ^*individual fuzzy group, the center of fuzzy group k, j is respectively v _k, v _j, then the training sample X after i-th standardization _ifor the degree of membership μ of fuzzy group k _ikfor:

In formula, n is the partitioned matrix index needed in fuzzy classification process, is usually taken as 2, || || be norm expression formula;

Use and be subordinate to angle value or its distortion to obtain new input matrix above, for fuzzy group k, its input matrix is deformed into:

Φ _ik(X _i,μ _ik)＝[1 func(μ _ik) X _i] (5)

Wherein func (μ _ik) for being subordinate to angle value μ _ikwarping function, generally get exp (μ _ik), Φ _ik(X _i, μ _ik) represent i-th input variable X _iand the degree of membership μ of fuzzy group k _ikcorresponding new input matrix;

Error backward propagation method, as the local equation of fuzzy system, if the prediction of a kth error back propagation fuzzy system output layer exports is be input as net, in the hidden layer that layer is adjacent therewith, the output of arbitrary neuron l is s _l, then have:

In formula, w _lkbe the connection weight between hidden neuron l and output layer neuron, f () is neuronic output function, is usually taken as Sigmoid function, is expressed as:

In formula, h is the neuronic threshold value of output layer, θ ₀for steepness parameter, in order to regulate the steepness of Sigmoid function; Finally, the output of last fuzzy system is obtained by the gravity model appoach in Anti-fuzzy method:

the prediction being fuzzy system exports;

4), adopt particle cluster algorithm to the w of fuzzy system medial error reverse transmittance nerve network local equation _lkbe optimized, specific implementation step is as follows:

1. determine that the Optimal Parameters of population is the w of error backward propagation method local equation _lk, population individual amount popsize, largest loop optimizing number of times iter _max, a p particle initial position r _p, initial velocity v _p, local optimum Lbest _pand the global optimum Gbest of whole population;

2. set optimization object function, be converted into fitness, each On Local Fuzzy system is evaluated; Calculate fitness function by corresponding error function, and think that the large particle fitness of error is little, the fitness function of particle p is expressed as:

f _p＝1/(E _p+1) (10)

In formula, E _pbe the error function of fuzzy system, be expressed as:

In formula, the prediction output of fuzzy system, O _ifor the target of fuzzy system exports, N is number of training;

3. according to following formula, circulation upgrades speed and the position of each particle,

v _p(iter+1)＝ω×v _p(iter)+m ₁a ₁(Lbest _p-r _p(iter))+m ₂a ₂(Gbest-r _p(iter))

(12)

r _p(iter+1)＝r _p(iter)+v _p(iter+1) (13)

In formula, v _prepresent the speed of more new particle p, r _prepresent the position of more new particle p, Lbest _prepresent the individual optimal value of more new particle p, Gbest is the training sample X after corresponding to i-th standardization _ifurnace temperature predicted value and make the performance variable value of furnace temperature the best, iter represents cycle index, and ω is the inertia weight in particle cluster algorithm, m ₁, m ₂corresponding accelerator coefficient, a ₁, a ₂it is the random number between [0,1];

4. for particle p, if new fitness is greater than original individual optimal value, the more individual optimal value of new particle:

Lbest _p＝f _p(14)

If the 5. individual optimal value Lbest of particle p _pbe greater than original population global optimum Gbest, upgrade original population global optimum Gbest:

Gbest＝Lbest _p(15)

6. judge whether to meet performance requirement, if so, terminate optimizing, obtain the local equation parameter of the fuzzy system that a group is optimized; Otherwise return step 3., continue iteration optimizing, until reach largest loop optimizing number of times iter _max;

Gbest is the training sample X after corresponding to i-th standardization _ifurnace temperature predicted value and make the performance variable value of furnace temperature the best;

Described method also comprises:

5), by the sampling time interval of setting, collection site intelligent instrument signal, the actual measurement furnace temperature function calculated value obtained is compared, if relative error be greater than 10% or furnace temperature exceed produce normal bound scope, then the new data of furnace temperature the best that makes when producing normal in DCS database is added training sample data, upgrade soft-sensing model;

6), in described step 4) in calculate best performance variable value, by the furnace temperature predicted value obtained and make the performance variable value of furnace temperature the best pass to DCS system, show at the control station of DCS, and be delivered to operator station by DCS system and fieldbus and show; Meanwhile, obtained makes the performance variable value of furnace temperature the best as new performance variable setting value by DCS system, automatically performs the operation of furnace temperature optimization;

Described key variables comprise enter incinerator waste liquid flow, enter the air mass flow of incinerator and enter the fuel flow rate of incinerator; Described performance variable comprises the air mass flow entering incinerator and the fuel flow rate entering incinerator.