CN103488206B - The optimizing temperature of pesticide production waste liquid incinerator system and method for intelligence radial basis - Google Patents

The optimizing temperature of pesticide production waste liquid incinerator system and method for intelligence radial basis Download PDF

Info

Publication number
CN103488206B
CN103488206B CN201310436865.8A CN201310436865A CN103488206B CN 103488206 B CN103488206 B CN 103488206B CN 201310436865 A CN201310436865 A CN 201310436865A CN 103488206 B CN103488206 B CN 103488206B
Authority
CN
China
Prior art keywords
mrow
msub
fuzzy
mfrac
math
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201310436865.8A
Other languages
Chinese (zh)
Other versions
CN103488206A (en
Inventor
刘兴高
李见会
张明明
孙优贤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201310436865.8A priority Critical patent/CN103488206B/en
Publication of CN103488206A publication Critical patent/CN103488206A/en
Application granted granted Critical
Publication of CN103488206B publication Critical patent/CN103488206B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Feedback Control In General (AREA)

Abstract

The invention discloses a kind of optimizing temperature of pesticide production waste liquid incinerator system and method for intelligent radial basis.The method uses radial base neural net as the local equation of fuzzy system, carrying out Fuzzy processing, reducing the impact of noise on Control for Kiln Temperature precision by exporting the mapping of training sample.In the present invention, standardization module, for training sample centralization; Fuzzy system module, for soft sensor modeling; Intelligent optimization module, for optimization system parameter; Model modification module presses the sampling time interval of setting, collection site intelligent instrument signal; 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; Signal acquisition module, for the time interval according to setting, image data from database.The real-time control, the control accuracy that present invention achieves furnace temperature are high, by avoiding to the suppression of noise the overshoot occurring furnace temperature.

Description

Intelligent radial-based pesticide production waste liquid incinerator temperature optimization system and method
Technical Field
The invention relates to the field of pesticide production waste liquid incineration, in particular to an intelligent radial-based pesticide production waste liquid incinerator temperature optimization system and method.
Background
With the rapid development of the pesticide industry, the environmental pollution problem of the emissions has attracted high attention from governments and corresponding environmental protection departments of various countries. The research and the solution of the standard-reaching discharge control and the harmless minimization treatment of the pesticide organic waste liquid not only become the difficulty and the hot point of scientific research of various countries, but also are the scientific proposition of the national urgent need related to the sustainable development of society.
The incineration method is the most effective and thorough method for treating pesticide residue and waste residue at present and is the most common method for application. The temperature of the incinerator must be kept at a proper temperature in the incineration process, and the excessively low incinerator temperature is not beneficial to the decomposition of toxic and harmful components in the waste; the overhigh furnace temperature not only increases the fuel consumption and the equipment operation cost, but also easily damages the inner wall of the hearth and shortens the service life of the equipment. In addition, excessive temperatures may increase the amount of metal volatilization and nitrogen oxide formation in the waste. Particularly for chlorine-containing wastewater, the corrosion of the inner wall can be reduced by proper furnace temperature. However, factors influencing the furnace temperature in the actual incineration process are complex and changeable, and the phenomenon that the furnace temperature is too low or too high is easy to occur.
In recent years, artificial neural networks, especially radial basis function neural networks, have achieved good results in terms of system optimization. The neural network has strong self-adaption, self-organization and self-learning capabilities and large-scale parallel operation capabilities. In practical applications, however, neural networks also expose some inherent drawbacks: the initialization of the weight is random and is easy to fall into local minimum; the selection of the number of nodes and other parameters of the hidden layer in the learning process can be selected only according to experience and experiments; too long convergence time, poor robustness, etc. Secondly, the DCS data collected in the industrial field also has certain uncertain errors due to noise, manual operation errors, and the like, so that a model using an artificial neural network with strong certainty is generally not strong in popularization capability.
Zadeh first proposed the concept of fuzzy aggregation in 1965 by american mathematician l. Fuzzy logic then begins to replace the classical logic that persists in that everything can be represented in terms of binary terms in a way that it more closely resembles the question and semantic statement of everyday people. Fuzzy logic has been successfully applied in various fields of industry, such as home appliances, industrial control, etc. In 2003, Demirci proposed the concept of fuzzy system, which constructs a new input matrix by using fuzzy membership matrix and its variants, and then derives the analytic value as the final output in the local equation by the centroid method in the inverse fuzzy method. For the system and the method for optimizing the temperature of the pesticide production waste liquid incinerator, the influence of noise and operation errors in the industrial production process are considered, and the influence of fuzzy performance of fuzzy logic on the accuracy can be reduced.
Particle Swarm Optimization, namely Particle Swarm Optimization, is a biological intelligent Optimization algorithm proposed by professors Kennedy and Eberhart to seek global optimality by simulating bird flight behavior, called PSO for short. The algorithm reduces the risk that the search algorithm is trapped in the local optimal solution through the mutual influence among particles in the group, and has good global search performance. The particle swarm algorithm is used for searching the optimal parameter combination of the radial basis fuzzy system so as to achieve the aim of optimizing the model.
Disclosure of Invention
In order to overcome the defects that the furnace temperature of the conventional incinerator is difficult to control and is easy to be too low or too high, the invention provides the system and the method for optimizing the furnace temperature of the pesticide production waste liquid incinerator, which realize accurate control of the furnace temperature and avoid the occurrence of too low or too high furnace temperature.
The technical scheme adopted by the invention for solving the technical problems is as follows:
the intelligent radial-base system for optimizing the furnace temperature of the pesticide production waste liquid incinerator comprises the incinerator, an intelligent instrument, a DCS (distributed control system), a data interface and an upper computer, wherein the DCS comprises a control station and a database; on-spot intelligent instrument and DCS headtotail, the DCS system is connected with the host computer, the host computer include: the standardization processing module is used for preprocessing the model training samples input from the DCS database, centralizing the training samples, namely subtracting the average value of the samples, and then standardizing the training samples:
calculating an average value: <math> <mrow> <mover> <mi>TX</mi> <mo>&OverBar;</mo> </mover> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>TX</mi> <mi>i</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </math>
calculating the variance: <math> <mrow> <msubsup> <mi>&sigma;</mi> <mi>x</mi> <mn>2</mn> </msubsup> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </mfrac> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <mrow> <mo>(</mo> <msub> <mi>TX</mi> <mi>i</mi> </msub> <mo>-</mo> <mover> <mi>TX</mi> <mo>&OverBar;</mo> </mover> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow> </math>
and (3) standardization: <math> <mrow> <mi>X</mi> <mo>=</mo> <mfrac> <mrow> <mi>TX</mi> <mo>-</mo> <mover> <mi>TX</mi> <mo>&OverBar;</mo> </mover> </mrow> <msub> <mi>&sigma;</mi> <mi>x</mi> </msub> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein, TXiThe ith training sample is data of key variables, furnace temperature and operating variables for optimizing the furnace temperature during normal production collected from the DCS database, N is the number of training samples,is the mean of the training samples, and X is the normalized training sample. σ x denotes the standard deviation of the training samples, σ2x represents the variance of the training samples.
And the fuzzy system module is used for fuzzifying the standardized training sample X transmitted from the data preprocessing module. Let c in fuzzy system*A fuzzy group, the centers of the fuzzy groups k and j are vk、vjThe ith normalized training sample XiMembership mu for fuzzy group kikComprises the following steps:
<math> <mrow> <msub> <mi>&mu;</mi> <mi>ik</mi> </msub> <mo>=</mo> <msup> <mrow> <mo>(</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msup> <mi>c</mi> <mo>*</mo> </msup> </munderover> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>v</mi> <mi>k</mi> </msub> <mo>|</mo> <mo>|</mo> </mrow> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>v</mi> <mi>j</mi> </msub> <mo>|</mo> <mo>|</mo> </mrow> </mfrac> <mo>)</mo> </mrow> <mfrac> <mn>2</mn> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </mfrac> </msup> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow> </math>
in the formula, n is a blocking matrix index required in the fuzzy classification process, and is usually taken as 2, | | · |, as a norm expression.
Using the above membership values or its variants to obtain a new input matrix, for the fuzzy group k, its input matrix variant is:
Φik(Xiik)=[1 func(μik) Xi] (5)
wherein func (. mu.)ik) Is a membership value muikDeformation function of, in general, takeexp(μik) Equal, phiik(Xiik) Denotes the ith input variable XiAnd membership mu of its fuzzy group kikThe corresponding new input matrix.
The radial basis function neural network is used as a local equation of the fuzzy system, and optimal fitting is carried out on each fuzzy group. The output of the kth particle swarm optimization radial basis fuzzy system is:
<math> <mrow> <msub> <mover> <mi>y</mi> <mo>^</mo> </mover> <mi>ik</mi> </msub> <mo>=</mo> <munder> <mi>&Sigma;</mi> <mi>l</mi> </munder> <msub> <mi>w</mi> <mi>lk</mi> </msub> <msub> <mi>&Phi;</mi> <mi>lk</mi> </msub> <mrow> <mo>(</mo> <mo>|</mo> <mo>|</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>C</mi> <mi>lk</mi> </msub> <mo>|</mo> <mo>|</mo> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein C islkAnd wlkOptimizing the center and output weight, phi, of the ith node of the radial basis fuzzy system by the kth particle swarm algorithmlk(||Xi-Clk| |) is the radial basis function of the ith node of the k particle swarm optimization radial basis fuzzy system, and is determined by the following formula:
<math> <mrow> <msub> <mi>&Phi;</mi> <mi>lk</mi> </msub> <mrow> <mo>(</mo> <mo>|</mo> <mo>|</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>C</mi> <mi>lk</mi> </msub> <mo>|</mo> <mo>|</mo> <mo>)</mo> </mrow> <mo>=</mo> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mfrac> <msup> <mrow> <mo>(</mo> <mo>|</mo> <mo>|</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>C</mi> <mi>lk</mi> </msub> <mo>|</mo> <mo>|</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mrow> <mn>2</mn> <mo>&times;</mo> <msub> <mi>&sigma;</mi> <mi>lk</mi> </msub> </mrow> </mfrac> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein sigmalkIs the gaussian width of the corresponding radial basis function, | | | · | |, is a norm expression. Finally, the center of gravity method in the anti-fuzzy method is used to obtain the final prediction output of the fuzzy equation system
<math> <mrow> <msub> <mover> <mi>y</mi> <mo>^</mo> </mover> <mi>i</mi> </msub> <mo>=</mo> <mfrac> <mrow> <msubsup> <mi>&Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <msup> <mi>c</mi> <mo>*</mo> </msup> </msubsup> <msub> <mi>&mu;</mi> <mi>ik</mi> </msub> <msub> <mover> <mi>y</mi> <mo>^</mo> </mover> <mi>ik</mi> </msub> </mrow> <mrow> <msubsup> <mi>&Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <msup> <mi>c</mi> <mo>*</mo> </msup> </msubsup> <msub> <mi>&mu;</mi> <mi>ik</mi> </msub> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow> </math>
In the formula,is the predicted output of the fuzzy equation system.
An intelligent optimization module for adopting particle swarm optimization to perform C of the radial basis function neural network local equation in the fuzzy systempk、σpk、wpkThe optimization is carried out, and the specific steps are as follows:
firstly, determining optimized parameters of particle swarm as C of radial basis function neural network local equationlk、σlk、wlkNumber of individual particle groups popsize, maximum number of cyclic optimization itermaxInitial position r of the p-th particlepInitial velocity vpLocal optimum value LbestpAnd a global optimum Gbest for the entire population of particles.
Secondly, setting an optimization objective function, converting the optimization objective function into fitness, and evaluating each local fuzzy equation; calculating a fitness function through a corresponding error function, considering that the fitness of the particle with large error is small, and expressing the fitness function of the particle p as follows:
fp=1/(Ep+1) (9)
in the formula, EpIs the error function of the fuzzy system, expressed as:
<math> <mrow> <msub> <mi>E</mi> <mi>p</mi> </msub> <mo>=</mo> <msqrt> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msup> <mrow> <mo>(</mo> <msub> <mover> <mi>y</mi> <mo>^</mo> </mover> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>O</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </msqrt> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow> </math>
in the formula,is the predicted output of a system of fuzzy equations, OiThe target output of the fuzzy equation system is shown, and N is the number of training samples;
thirdly, according to the following formula, the speed and the position of each particle are circularly updated,
vp(iter+1)=ω×vp(iter)+m1a1(Lbestp-rp(iter))+m2a2(Gbest-rp(iter))
(11)
rp(iter+1)=rp(iter)+vp(iter+1) (12)
in the formula, vpIndicates the velocity, r, of the update particle ppIndicating the position of the update particle p, LbestpRepresenting the individual optimum value of the update particle p, Gbest being that corresponding to the ith normalizedTraining sample XiIs a value of an operating variable for optimizing the furnace temperature, iter represents the number of cycles, ω is an inertial weight in a particle swarm algorithm, m1、m2Is the corresponding acceleration factor, a1、a2Is [0,1 ]]A random number in between;
and fourthly, for the particle p, if the new fitness is larger than the original individual optimal value, updating the individual optimal value of the particle:
Lbestp=fp (13)
if the individual optimum value Lbest of the particle ppAnd (3) updating the original particle swarm global optimum value Gbest if the particle swarm global optimum value Gbest is larger than the original particle swarm global optimum value Gbest:
Gbest=Lbestp (14)
sixthly, judging whether the performance requirement is met, if so, finishing the optimization to obtain a group of optimized local equation parameters of the fuzzy system; otherwise, returning to the step III, continuing the iteration optimization until reaching the maximum iteration number itermax
Gbest is the training sample X corresponding to the ith standardiAnd an operating variable value for optimizing the furnace temperature.
As a preferred solution: the host computer still include: and the model updating module is used for acquiring field intelligent instrument signals according to a set sampling time interval, comparing the obtained actually-measured furnace temperature with a system forecast value, and adding new data which enables the furnace temperature to be optimal when the furnace temperature is normally produced in the DCS database into the training sample data to update the soft measurement model if the relative error is more than 10% or the furnace temperature exceeds the upper and lower normal production limit ranges.
Further, the host computer still include: the result display module is used for transmitting the obtained furnace temperature forecast value and the operation variable value which enables the furnace temperature to be optimal to the DCS, displaying the values at a control station of the DCS and transmitting the values to a field operation station for displaying through the DCS and a field bus; at the same time, the DCS system automatically executes the furnace temperature optimization operation by using the obtained value of the operation variable that optimizes the furnace temperature as a new operation variable set value.
And the signal acquisition module is used for acquiring data from the database according to the set time interval of each sampling.
Still further, the key variables include the flow of waste liquid into the incinerator, the flow of air into the incinerator, and the flow of fuel into the incinerator; the manipulated variables include air flow into the incinerator and fuel flow into the incinerator.
The method for optimizing the furnace temperature of the pesticide production waste liquid incinerator by using the intelligent radial basis system comprises the following specific implementation steps of:
1) determining used key variables, collecting data of the variables in normal production from a DCS (distributed control system) database as an input matrix of a training sample TX, and collecting corresponding furnace temperature and operation variable data for optimizing the furnace temperature as an output matrix O;
2) preprocessing a model training sample input from a DCS database, centralizing the training sample, namely subtracting the average value of the sample, and then normalizing the training sample so that the average value is 0 and the variance is 1. The processing is accomplished using the following mathematical process:
2.1) calculating the mean value: <math> <mrow> <mover> <mi>TX</mi> <mo>&OverBar;</mo> </mover> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>TX</mi> <mi>i</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </math>
2.2) calculating the variance: <math> <mrow> <msubsup> <mi>&sigma;</mi> <mi>x</mi> <mn>2</mn> </msubsup> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </mfrac> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <mrow> <mo>(</mo> <msub> <mi>TX</mi> <mi>i</mi> </msub> <mo>-</mo> <mover> <mi>TX</mi> <mo>&OverBar;</mo> </mover> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow> </math>
2.3) standardization: <math> <mrow> <mi>X</mi> <mo>=</mo> <mfrac> <mrow> <mi>TX</mi> <mo>-</mo> <mover> <mi>TX</mi> <mo>&OverBar;</mo> </mover> </mrow> <msub> <mi>&sigma;</mi> <mi>x</mi> </msub> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow> </math> wherein, TXiThe ith training sample is data of key variables, furnace temperature and operating variables for optimizing the furnace temperature during normal production collected from the DCS database, N is the number of training samples,is the mean of the training samples, and X is the normalized training sample. σ x denotes the standard deviation of the training samples, and σ 2x denotes the variance of the training samples.
3) And fuzzifying the training sample transmitted from the data preprocessing module. Let c in fuzzy system*One fuzzy group, among fuzzy groups k, jHeart is v respectivelyk、vjThe ith normalized training sample XiMembership mu for fuzzy group kikComprises the following steps:
<math> <mrow> <msub> <mi>&mu;</mi> <mi>ik</mi> </msub> <mo>=</mo> <msup> <mrow> <mo>(</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msup> <mi>c</mi> <mo>*</mo> </msup> </munderover> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>v</mi> <mi>k</mi> </msub> <mo>|</mo> <mo>|</mo> </mrow> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>v</mi> <mi>j</mi> </msub> <mo>|</mo> <mo>|</mo> </mrow> </mfrac> <mo>)</mo> </mrow> <mfrac> <mn>2</mn> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </mfrac> </msup> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow> </math>
in the formula, n is a blocking matrix index required in the fuzzy classification process, and is usually taken as 2, | | · |, as a norm expression.
Using the above membership values or its variants to obtain a new input matrix, for the fuzzy group k, its input matrix variant is:
Φik(Xiik)=[1 func(μik) Xi] (5)
wherein func (. mu.)ik) Is a membership value muikDeformation function of, in general, takeexp(μik) Equal, phiik(Xiik) Denotes the ith input variable XiAnd membership mu of its fuzzy group kikThe corresponding new input matrix.
The radial basis function neural network is used as a local equation of the fuzzy system, and optimal fitting is carried out on each fuzzy group. The output of the kth particle swarm optimization radial basis fuzzy system is:
<math> <mrow> <msub> <mover> <mi>y</mi> <mo>^</mo> </mover> <mi>ik</mi> </msub> <mo>=</mo> <munder> <mi>&Sigma;</mi> <mi>l</mi> </munder> <msub> <mi>w</mi> <mi>lk</mi> </msub> <msub> <mi>&Phi;</mi> <mi>lk</mi> </msub> <mrow> <mo>(</mo> <mo>|</mo> <mo>|</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>C</mi> <mi>lk</mi> </msub> <mo>|</mo> <mo>|</mo> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein C islkAnd wlkOptimizing the center and output weight, phi, of the ith node of the radial basis fuzzy system by the kth particle swarm algorithmlk(||Xi-Clk| |) is the radial basis function of the ith node of the k particle swarm optimization radial basis fuzzy system, and is determined by the following formula:
<math> <mrow> <msub> <mi>&Phi;</mi> <mi>lk</mi> </msub> <mo>=</mo> <mrow> <mo>(</mo> <mo>|</mo> <mo>|</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>C</mi> <mi>lk</mi> </msub> <mo>|</mo> <mo>|</mo> <mo>)</mo> </mrow> <mo>=</mo> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mfrac> <msup> <mrow> <mo>(</mo> <mo>|</mo> <mo>|</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>C</mi> <mi>lk</mi> </msub> <mo>|</mo> <mo>|</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mrow> <mn>2</mn> <mo>&times;</mo> <msub> <mi>&sigma;</mi> <mi>lk</mi> </msub> </mrow> </mfrac> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein sigmalkIs the gaussian width of the corresponding radial basis function, | | | · | |, is a norm expression. And finally, obtaining the final prediction output of the fuzzy equation system by a gravity center method in an anti-fuzzy method:
<math> <mrow> <msub> <mover> <mi>y</mi> <mo>^</mo> </mover> <mi>i</mi> </msub> <mo>=</mo> <mfrac> <mrow> <msubsup> <mi>&Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <msup> <mi>c</mi> <mo>*</mo> </msup> </msubsup> <msub> <mi>&mu;</mi> <mi>ik</mi> </msub> <msub> <mover> <mi>y</mi> <mo>^</mo> </mover> <mi>ik</mi> </msub> </mrow> <mrow> <msubsup> <mi>&Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <msup> <mi>c</mi> <mo>*</mo> </msup> </msubsup> <msub> <mi>&mu;</mi> <mi>ik</mi> </msub> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow> </math>
in the formula,is the predicted output of the fuzzy equation system.
4) C for radial basis function neural network local equation in fuzzy system by adopting particle swarm optimizationpk、σpk、wpkThe optimization is carried out, and the specific implementation steps are as follows:
firstly, determining optimized parameters of particle swarm as C of radial basis function neural network local equationlk、σlk、wlkNumber of individual particle groups popsize, maximum number of cyclic optimization itermaxInitial position r of the p-th particlepInitial velocity vpLocal optimum value LbestpAnd a global optimum Gbest for the entire population of particles.
Secondly, setting an optimization objective function, converting the optimization objective function into fitness, and evaluating each local fuzzy equation; calculating a fitness function through a corresponding error function, considering that the fitness of the particle with large error is small, and expressing the fitness function of the particle p as follows:
fp=1/(Ep+1) (9)
in the formula, EpIs the error function of the fuzzy system, expressed as:
<math> <mrow> <msub> <mi>E</mi> <mi>p</mi> </msub> <mo>=</mo> <msqrt> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msup> <mrow> <mo>(</mo> <msub> <mover> <mi>y</mi> <mo>^</mo> </mover> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>O</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </msqrt> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow> </math>
in the formula,is the prediction output of the fuzzy system, OiThe target output of the fuzzy equation system is shown, and N is the number of training samples;
thirdly, according to the following formula, the speed and the position of each particle are circularly updated,
vp(iter+1)=ω×vp(iter)+m1a1(Lbestp-rp(iter))+m2a2(Gbest-rp(iter))
(11)
rp(iter+1)=rp(iter)+vp(iter+1) (12)
in the formula, vpIndicates the velocity, r, of the update particle ppIndicating the position of the update particle p, LbestpRepresenting the individual optimum value of the update particle p, Gbest being the training sample X corresponding to the ith normalizationiIs a value of an operating variable for optimizing the furnace temperature, iter represents the number of cycles, ω is an inertial weight in a particle swarm algorithm, m1、m2Is the corresponding acceleration factor, a1、a2Is [0,1 ]]A random number in between;
and fourthly, for the particle p, if the new fitness is larger than the original individual optimal value, updating the individual optimal value of the particle:
Lbestp=fp (13)
if the individual optimum value Lbest of the particle ppAnd (3) updating the original particle swarm global optimum value Gbest if the particle swarm global optimum value Gbest is larger than the original particle swarm global optimum value Gbest:
Gbest=Lbestp (14)
sixthly, judging whether the performance requirement is met, if so, finishing the optimization to obtain a group of optimized local equation parameters of the fuzzy system; otherwise, returning to the step III, continuing the iteration optimization until reaching the maximum iteration number itermax
Gbest is the training sample X corresponding to the ith standardiAnd an operating variable value for optimizing the furnace temperature.
As a preferred solution: the method further comprises the following steps: 5) acquiring on-site intelligent instrument signals according to a set sampling time interval, comparing the obtained actually-measured furnace temperature with a system forecast value, and if the relative error is more than 10% or the furnace temperature exceeds the upper and lower production normal limit ranges, adding new data which enables the furnace temperature to be optimal when the furnace temperature is normally produced in the DCS database into training sample data, and updating the soft measurement model.
6) Further, calculating to obtain an optimal operating variable value in the step 4), transmitting the obtained furnace temperature forecast value and the operating variable value which enables the furnace temperature to be optimal to a DCS (distributed control system), displaying the values on a control station of the DCS, and transmitting the values to a field operating station through the DCS and a field bus for displaying; at the same time, the DCS system automatically executes the furnace temperature optimization operation by using the obtained value of the operation variable that optimizes the furnace temperature as a new operation variable set value.
Still further, the key variables include the flow of waste liquid into the incinerator, the flow of air into the incinerator, and the flow of fuel into the incinerator; the manipulated variables include air flow into the incinerator and fuel flow into the incinerator.
The technical conception of the invention is as follows: the invention discloses an intelligent radial basis system and method for optimizing the furnace temperature of a pesticide production waste liquid incinerator, and aims to find a furnace temperature forecast value and an operation variable value for optimizing the furnace temperature.
The invention has the following beneficial effects: 1. establishing an online soft measurement model of a quantitative relation between a system key variable and furnace temperature; 2. the operating conditions that optimize the furnace temperature are quickly found.
Drawings
FIG. 1 is a hardware block diagram of the system proposed by the present invention;
fig. 2 is a functional structure diagram of the upper computer according to the present invention.
Detailed Description
The invention is further described below with reference to the accompanying drawings. The examples are intended to illustrate the invention, but not to limit the invention, and any modifications and variations of the invention within the spirit and scope of the claims are intended to fall within the scope of the invention.
Example 1
Referring to fig. 1 and 2, the system for optimizing the furnace temperature of the intelligent radial-based pesticide production waste liquid incinerator comprises a field intelligent instrument 2 connected with an incinerator object 1, a DCS system and an upper computer 6, wherein the DCS system comprises a data interface 3, a control station 4 and a database 5, the field intelligent instrument 2 is connected with the data interface 3, the data interface is connected with the control station 4, the database 5 and the upper computer 6, and the upper computer 6 comprises:
a normalization processing module 7, configured to pre-process the model training samples input from the DCS database, centralize the training samples, that is, subtract the average value of the samples, and then normalize them:
calculating an average value: <math> <mrow> <mover> <mi>TX</mi> <mo>&OverBar;</mo> </mover> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>TX</mi> <mi>i</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </math>
calculating the variance: <math> <mrow> <msubsup> <mi>&sigma;</mi> <mi>x</mi> <mn>2</mn> </msubsup> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </mfrac> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <mrow> <mo>(</mo> <msub> <mi>TX</mi> <mi>i</mi> </msub> <mo>-</mo> <mover> <mi>TX</mi> <mo>&OverBar;</mo> </mover> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow> </math>
and (3) standardization: <math> <mrow> <mi>X</mi> <mo>=</mo> <mfrac> <mrow> <mi>TX</mi> <mo>-</mo> <mover> <mi>TX</mi> <mo>&OverBar;</mo> </mover> </mrow> <msub> <mi>&sigma;</mi> <mi>x</mi> </msub> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein, TXiThe ith training sample is data of key variables, furnace temperature and operating variables for optimizing the furnace temperature during normal production collected from the DCS database, N is the number of training samples,is the mean of the training samples, and X is the normalized training sample. σ x denotes the standard deviation of the training samples, and σ 2x denotes the variance of the training samples.
And the fuzzy system module 8 is used for fuzzifying the standardized training sample X transmitted from the data preprocessing module. Let c in fuzzy system*A fuzzy group, the centers of the fuzzy groups k and j are vk、vjThe ith normalized training sample XiMembership mu for fuzzy group kikComprises the following steps:
<math> <mrow> <msub> <mi>&mu;</mi> <mi>ik</mi> </msub> <mo>=</mo> <msup> <mrow> <mo>(</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msup> <mi>c</mi> <mo>*</mo> </msup> </munderover> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>v</mi> <mi>k</mi> </msub> <mo>|</mo> <mo>|</mo> </mrow> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>v</mi> <mi>j</mi> </msub> <mo>|</mo> <mo>|</mo> </mrow> </mfrac> <mo>)</mo> </mrow> <mfrac> <mn>2</mn> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </mfrac> </msup> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow> </math>
in the formula, n is a blocking matrix index required in the fuzzy classification process, and is usually taken as 2, | | · |, as a norm expression.
Using the above membership values or its variants to obtain a new input matrix, for the fuzzy group k, its input matrix variant is:
Φik(Xiik)=[1 func(μik) Xi] (5)
wherein func (. mu.)ik) Is a membership value muikDeformation function of, in general, takeexp(μik) Equal, phiik(Xiik) Denotes the ith input variable XiAnd membership mu of its fuzzy group kikThe corresponding new input matrix.
The radial basis function neural network is used as a local equation of the fuzzy system, and optimal fitting is carried out on each fuzzy group. The output of the kth particle swarm optimization radial basis fuzzy system is:
<math> <mrow> <msub> <mover> <mi>y</mi> <mo>^</mo> </mover> <mi>ik</mi> </msub> <mo>=</mo> <munder> <mi>&Sigma;</mi> <mi>l</mi> </munder> <msub> <mi>w</mi> <mi>lk</mi> </msub> <msub> <mi>&Phi;</mi> <mi>lk</mi> </msub> <mrow> <mo>(</mo> <mo>|</mo> <mo>|</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>C</mi> <mi>lk</mi> </msub> <mo>|</mo> <mo>|</mo> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein C islkAnd wlkOptimizing the center and output weight, phi, of the ith node of the radial basis fuzzy system by the kth particle swarm algorithmlk(||Xi-Clk| |) is the radial basis function of the ith node of the k particle swarm optimization radial basis fuzzy system, and is determined by the following formula:
<math> <mrow> <msub> <mi>&Phi;</mi> <mi>lk</mi> </msub> <mrow> <mo>(</mo> <mo>|</mo> <mo>|</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>C</mi> <mi>lk</mi> </msub> <mo>|</mo> <mo>|</mo> <mo>)</mo> </mrow> <mo>=</mo> <mi>exp</mi> <mrow> <mo>(</mo> <mfrac> <msup> <mrow> <mo>(</mo> <mo>|</mo> <mo>|</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>C</mi> <mi>lk</mi> </msub> <mo>|</mo> <mo>|</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mrow> <mn>2</mn> <mo>&times;</mo> <msub> <mi>&sigma;</mi> <mi>lk</mi> </msub> </mrow> </mfrac> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein sigmalkIs the gaussian width of the corresponding radial basis function, | | | · | |, is a norm expression. And finally, obtaining the final prediction output of the fuzzy equation system by a gravity center method in an anti-fuzzy method:
<math> <mrow> <msub> <mover> <mi>y</mi> <mo>^</mo> </mover> <mi>i</mi> </msub> <mo>=</mo> <mfrac> <mrow> <msubsup> <mi>&Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <msup> <mi>c</mi> <mo>*</mo> </msup> </msubsup> <msub> <mi>&mu;</mi> <mi>ik</mi> </msub> <msub> <mover> <mi>y</mi> <mo>^</mo> </mover> <mi>ik</mi> </msub> </mrow> <mrow> <msubsup> <mi>&Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <msup> <mi>c</mi> <mo>*</mo> </msup> </msubsup> <msub> <mi>&mu;</mi> <mi>ik</mi> </msub> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow> </math>
in the formula,is the predicted output of the fuzzy equation system.
An intelligent optimization module 9 for adopting the particle swarm algorithm to perform C on the radial basis function neural network local equation in the fuzzy systempk、σpk、wpkThe optimization is carried out, and the specific implementation steps are as follows:
firstly, determining optimized parameters of particle swarm as C of radial basis function neural network local equationlk、σlk、wlkNumber of individual particle groups popsize, maximum number of cyclic optimization itermaxInitial position r of the p-th particlepInitial velocity vpLocal optimum value LbestpAnd a global optimum Gbest for the entire population of particles.
Secondly, setting an optimization objective function, converting the optimization objective function into fitness, and evaluating each local fuzzy equation; calculating a fitness function through a corresponding error function, considering that the fitness of the particle with large error is small, and expressing the fitness function of the particle p as follows:
fp=1/(Ep+1) (9)
in the formula, EpIs the error function of the fuzzy system, expressed as:
<math> <mrow> <msub> <mi>E</mi> <mi>p</mi> </msub> <mo>=</mo> <msqrt> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msup> <mrow> <mo>(</mo> <msub> <mover> <mi>y</mi> <mo>^</mo> </mover> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>O</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </msqrt> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow> </math>
in the formula,is the prediction output of the fuzzy system, OiThe target output of the fuzzy equation system is shown, and N is the number of training samples;
thirdly, according to the following formula, the speed and the position of each particle are circularly updated,
vp(iter+1)=ω×vp(iter)+m1a1(Lbestp-rp(iter))+m2a2(Gbest-rp(iter))
(11)
rp(iter+1)=rp(iter)+vp(iter+1) (12)
in the formula, vpIndicates the velocity, r, of the update particle ppIndicating the position of the update particle p, LbestpTo representUpdating the individual optimal value of the particle p, wherein Gbest is the training sample X corresponding to the ith normalizationiIs a value of an operating variable for optimizing the furnace temperature, iter represents the number of cycles, ω is an inertial weight in a particle swarm algorithm, m1、m2Is the corresponding acceleration factor, a1、a2Is [0,1 ]]A random number in between;
and fourthly, for the particle p, if the new fitness is larger than the original individual optimal value, updating the individual optimal value of the particle:
Lbestp=fp (13)
if the individual optimum value Lbest of the particle ppAnd (3) updating the original particle swarm global optimum value Gbest if the particle swarm global optimum value Gbest is larger than the original particle swarm global optimum value Gbest:
Gbest=Lbestp (14)
sixthly, judging whether the performance requirement is met, if so, finishing the optimization to obtain a group of optimized local equation parameters of the fuzzy system; otherwise, returning to the step III, continuing the iteration optimization until reaching the maximum iteration number itermax
Gbest is the training sample X corresponding to the ith standardiAnd an operating variable value for optimizing the furnace temperature.
The upper computer 6 further comprises: and the signal acquisition module 11 is used for acquiring data from the database according to a set time interval of each sampling.
The upper computer 6 further comprises: and the model updating module 12 is used for acquiring field intelligent instrument signals according to a set sampling time interval, comparing the obtained actually-measured furnace temperature with a system forecast value, and adding new data which enables the furnace temperature to be optimal when the furnace temperature is normally produced in the DCS database into the training sample data to update the soft measurement model if the relative error is more than 10% or the furnace temperature exceeds the upper and lower normal production limit ranges.
The key variables include the flow of waste liquid into the incinerator, the flow of air into the incinerator and the flow of fuel into the incinerator; the manipulated variables include air flow into the incinerator and fuel flow into the incinerator.
The system also comprises a DCS (distributed control system), wherein the DCS is composed of a data interface 3, a control station 4 and a database 5; the intelligent instrument 2, the DCS system and the upper computer 6 are sequentially connected through a field bus; the upper computer 6 further comprises a result display module 10, which is used for transmitting the calculation optimal result to the DCS system, displaying the process state at the control station of the DCS, and transmitting the process state information to the field operation station for displaying through the DCS system and the field bus.
When the waste liquid incineration process is provided with the DCS system, the functions of obtaining the furnace temperature forecast value and the operation variable value for optimizing the furnace temperature are mainly completed on the upper computer by utilizing the real-time and historical databases of the DCS system to detect and store the real-time dynamic data of the sample.
When the waste liquid incineration process is not provided with the DCS system, the data memory is adopted to replace the data storage function of a real-time and historical database of the DCS system, and the functional system for obtaining the furnace temperature forecast value and the operation variable value for optimizing the furnace temperature is manufactured into an independent complete system-on-chip which comprises an I/O element, a data memory, a program memory, an arithmetic unit and a display module and does not depend on the DCS system, so that the system-on-chip can be independently used regardless of whether the incineration process is provided with the DCS or not, and is more beneficial to popularization and use.
Example 2
Referring to fig. 1 and 2, the method for optimizing the temperature of the pesticide production waste liquid incinerator with the intelligent radial basis specifically comprises the following implementation steps:
1) determining used key variables, collecting data of the variables in normal production from a DCS (distributed control system) database as an input matrix of a training sample TX, and collecting corresponding furnace temperature and operation variable data for optimizing the furnace temperature as an output matrix O;
2) preprocessing a model training sample input from a DCS database, centralizing the training sample, namely subtracting the average value of the sample, and then normalizing the training sample so that the average value is 0 and the variance is 1. The processing is accomplished using the following mathematical process:
2.1) calculating the mean value: <math> <mrow> <mover> <mi>TX</mi> <mo>&OverBar;</mo> </mover> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>TX</mi> <mi>i</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </math>
2.2) calculating the variance: <math> <mrow> <msubsup> <mi>&sigma;</mi> <mi>x</mi> <mn>2</mn> </msubsup> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </mfrac> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <mrow> <mo>(</mo> <msub> <mi>TX</mi> <mi>i</mi> </msub> <mo>-</mo> <mover> <mi>TX</mi> <mo>&OverBar;</mo> </mover> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow> </math>
2.3) standardization: <math> <mrow> <mi>X</mi> <mo>=</mo> <mfrac> <mrow> <mi>TX</mi> <mo>-</mo> <mover> <mi>TX</mi> <mo>&OverBar;</mo> </mover> </mrow> <msub> <mi>&sigma;</mi> <mi>x</mi> </msub> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein, TXiThe ith training sample is data of key variables, furnace temperature and operating variables for optimizing the furnace temperature during normal production collected from the DCS database, N is the number of training samples,is the mean of the training samples, and X is the normalized training sample. σ x denotes the standard deviation of the training samples, and σ 2x denotes the variance of the training samples.
3) And fuzzifying the training sample transmitted from the data preprocessing module. Let c in fuzzy system*A fuzzy group, the centers of the fuzzy groups k and j are vk、vjThe ith normalized training sample XiMembership mu for fuzzy group kikComprises the following steps:
<math> <mrow> <msub> <mi>&mu;</mi> <mi>ik</mi> </msub> <mo>=</mo> <msup> <mrow> <mo>(</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <msup> <mi>c</mi> <mo>*</mo> </msup> </munderover> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>v</mi> <mi>k</mi> </msub> <mo>|</mo> <mo>|</mo> </mrow> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>v</mi> <mi>j</mi> </msub> <mo>|</mo> <mo>|</mo> </mrow> </mfrac> <mo>)</mo> </mrow> <mfrac> <mn>2</mn> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </mfrac> </msup> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow> </math>
in the formula, n is a blocking matrix index required in the fuzzy classification process, and is usually taken as 2, | | · |, as a norm expression.
Using the above membership values or its variants to obtain a new input matrix, for the fuzzy group k, its input matrix variant is:
Φik(Xiik)=[1 func(μik) Xi] (5)
wherein func (. mu.)ik) Is a membership value muikDeformation function of, in general, takeexp(μik) Equal, phiik(Xiik) Denotes the ith input variable XiAnd membership mu of its fuzzy group kikThe corresponding new input matrix.
The radial basis function neural network is used as a local equation of the fuzzy system, and optimal fitting is carried out on each fuzzy group. The output of the kth particle swarm optimization radial basis fuzzy system is:
<math> <mrow> <msub> <mover> <mi>y</mi> <mo>^</mo> </mover> <mi>ik</mi> </msub> <mo>=</mo> <munder> <mi>&Sigma;</mi> <mi>l</mi> </munder> <msub> <mi>w</mi> <mi>lk</mi> </msub> <msub> <mi>&Phi;</mi> <mi>lk</mi> </msub> <mrow> <mo>(</mo> <mo>|</mo> <mo>|</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>C</mi> <mi>lk</mi> </msub> <mo>|</mo> <mo>|</mo> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein C islkAnd wlkOptimizing the center and output weight, phi, of the ith node of the radial basis fuzzy system by the kth particle swarm algorithmlk(||Xi-Clk| |) is the radial basis function of the ith node of the k particle swarm optimization radial basis fuzzy system, and is determined by the following formula:
<math> <mrow> <msub> <mi>&Phi;</mi> <mi>lk</mi> </msub> <mo>=</mo> <mrow> <mo>(</mo> <mo>|</mo> <mo>|</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>C</mi> <mi>lk</mi> </msub> <mo>|</mo> <mo>|</mo> <mo>)</mo> </mrow> <mo>=</mo> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mfrac> <msup> <mrow> <mo>(</mo> <mo>|</mo> <mo>|</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>C</mi> <mi>lk</mi> </msub> <mo>|</mo> <mo>|</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mrow> <mn>2</mn> <mo>&times;</mo> <msub> <mi>&sigma;</mi> <mi>lk</mi> </msub> </mrow> </mfrac> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein sigmalkIs the gaussian width of the corresponding radial basis function, | | | · | |, is a norm expression. And finally, obtaining the final prediction output of the fuzzy equation system by a gravity center method in an anti-fuzzy method:
<math> <mrow> <msub> <mover> <mi>y</mi> <mo>^</mo> </mover> <mi>i</mi> </msub> <mo>=</mo> <mfrac> <mrow> <msubsup> <mi>&Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <msup> <mi>c</mi> <mo>*</mo> </msup> </msubsup> <msub> <mi>&mu;</mi> <mi>ik</mi> </msub> <msub> <mover> <mi>y</mi> <mo>^</mo> </mover> <mi>ik</mi> </msub> </mrow> <mrow> <msubsup> <mi>&Sigma;</mi> <mrow> <mi>k</mi> <mo>-</mo> <mn>1</mn> </mrow> <msup> <mi>c</mi> <mo>*</mo> </msup> </msubsup> <msub> <mi>&mu;</mi> <mi>ik</mi> </msub> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow> </math>
in the formula,is the predicted output of the fuzzy equation system.
4) C for radial basis function neural network local equation in fuzzy system by adopting particle swarm optimizationpk、σpk、wpkThe optimization is carried out, and the specific implementation steps are as follows:
firstly, determining optimized parameters of particle swarm as C of radial basis function neural network local equationlk、σlk、wlkNumber of individual particle groups popsize, maximum number of cyclic optimization itermaxInitial position r of the p-th particlepInitial velocity vpLocal optimum value LbestpAnd a global optimum Gbest for the entire population of particles.
Secondly, setting an optimization objective function, converting the optimization objective function into fitness, and evaluating each local fuzzy equation; calculating a fitness function through a corresponding error function, considering that the fitness of the particle with large error is small, and expressing the fitness function of the particle p as follows:
fp=1/(Ep+1) (9)
in the formula, EpIs the error function of the fuzzy system, expressed as:
<math> <mrow> <msub> <mi>E</mi> <mi>p</mi> </msub> <mo>=</mo> <msqrt> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msup> <mrow> <mo>(</mo> <msub> <mover> <mi>y</mi> <mo>^</mo> </mover> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>O</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </msqrt> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow> </math>
in the formula,is the prediction output of the fuzzy system, OiThe target output of the fuzzy equation system is shown, and N is the number of training samples;
thirdly, according to the following formula, the speed and the position of each particle are circularly updated,
vp(iter+1)=ω×vp(iter)+m1a1(Lbestp-rp(iter))+m2a2(Gbest-rp(iter))
(11)
rp(iter+1)=rp(iter)+vp(iter+1)(12)
in the formula, vpIndicates the velocity, r, of the update particle ppIndicating the position of the update particle p, LbestpRepresenting the individual optimum value of the update particle p, Gbest being the training sample X corresponding to the ith normalizationiIs a value of an operating variable for optimizing the furnace temperature, iter represents the number of cycles, ω is an inertial weight in a particle swarm algorithm, m1、m2Is the corresponding acceleration factor, a1、a2Is [0,1 ]]A random number in between;
and fourthly, for the particle p, if the new fitness is larger than the original individual optimal value, updating the individual optimal value of the particle:
Lbestp=fp (13)
if the individual optimum value Lbest of the particle ppAnd (3) updating the original particle swarm global optimum value Gbest if the particle swarm global optimum value Gbest is larger than the original particle swarm global optimum value Gbest:
Gbest=Lbestp (14)
sixthly, judging whether the performance requirement is met, if so, finishing the optimization to obtain a group of optimized local equation parameters of the fuzzy system; otherwise, returning to the step III, continuing the iteration optimization until reaching the maximum iteration number itermax
Gbest is the training sample X corresponding to the ith standardiAnd an operating variable value for optimizing the furnace temperature.
The method further comprises the following steps: 5) acquiring on-site intelligent instrument signals according to a set sampling time interval, comparing the obtained actually-measured furnace temperature with a system forecast value, and if the relative error is more than 10% or the furnace temperature exceeds the upper and lower production normal limit ranges, adding new data which enables the furnace temperature to be optimal when the furnace temperature is normally produced in the DCS database into training sample data, and updating the soft measurement model.
Calculating to obtain an optimal operation variable value in the step 4), transmitting the obtained furnace temperature forecast value and the operation variable value which enables the furnace temperature to be optimal to the DCS, displaying the values at a control station of the DCS, and transmitting the values to a field operation station for displaying through the DCS and a field bus; at the same time, the DCS system automatically executes the furnace temperature optimization operation by using the obtained value of the operation variable that optimizes the furnace temperature as a new operation variable set value.
The key variables include the flow of waste liquid into the incinerator, the flow of air into the incinerator and the flow of fuel into the incinerator; the manipulated variables include air flow into the incinerator and fuel flow into the incinerator.

Claims (2)

1. An intelligent radial-based pesticide production waste liquid incinerator temperature optimization system comprises an incinerator, a field intelligent instrument, a DCS (distributed control system), a data interface and an upper computer, wherein the DCS comprises a control station and a database; on-spot intelligent instrument and DCS headtotail, the DCS system is connected with the host computer, its characterized in that: the host computer include:
the standardization processing module is used for preprocessing the model training samples input from the DCS database, centralizing the training samples, namely subtracting the average value of the samples, and then standardizing the training samples:
calculating an average value: <math> <mrow> <mover> <mrow> <mi>T</mi> <mi>X</mi> </mrow> <mo>&OverBar;</mo> </mover> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>TX</mi> <mi>i</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </math>
calculating the variance: <math> <mrow> <msubsup> <mi>&sigma;</mi> <mi>x</mi> <mn>2</mn> </msubsup> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </mfrac> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <mrow> <mo>(</mo> <msub> <mi>TX</mi> <mi>i</mi> </msub> <mo>-</mo> <mover> <mrow> <mi>T</mi> <mi>X</mi> </mrow> <mo>&OverBar;</mo> </mover> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow> </math>
and (3) standardization: <math> <mrow> <mi>X</mi> <mo>=</mo> <mfrac> <mrow> <mi>T</mi> <mi>X</mi> <mo>-</mo> <mover> <mrow> <mi>T</mi> <mi>X</mi> </mrow> <mo>&OverBar;</mo> </mover> </mrow> <msub> <mi>&sigma;</mi> <mi>x</mi> </msub> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein, TXiThe ith training sample is data of key variables, furnace temperature and operating variables for optimizing the furnace temperature during normal production collected from the DCS database, N is the number of training samples,is the mean value of the training sample, and X is the training sample after standardization; sigmaxRepresenting the standard deviation, σ, of the training samples2 xRepresenting the variance of the training samples;
the fuzzy system module is used for fuzzifying the standardized training sample X transmitted from the standardization processing module; let c in fuzzy system*A fuzzy group, the centers of the fuzzy groups k and j are vk、vjThe ith normalized training sample XiMembership mu for fuzzy group kikComprises the following steps:
<math> <mrow> <msub> <mi>&mu;</mi> <mrow> <mi>i</mi> <mi>k</mi> </mrow> </msub> <mo>=</mo> <msup> <mrow> <mo>(</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mrow> <mi>c</mi> <mo>*</mo> </mrow> </munderover> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>v</mi> <mi>k</mi> </msub> <mo>|</mo> <mo>|</mo> </mrow> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>v</mi> <mi>j</mi> </msub> <mo>|</mo> <mo>|</mo> </mrow> </mfrac> <mo>)</mo> </mrow> <mfrac> <mn>2</mn> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </mfrac> </msup> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow> </math>
in the formula, n is a blocking matrix index required in the fuzzy classification process, and is taken as 2, | | · | | is a norm expression;
using the above membership values or its variants to obtain a new input matrix, for the fuzzy group k, its input matrix variant is:
Φik(Xiik)=[1 func(μik) Xi] (5)
wherein func (. mu.)ik) Is a membership value muikIs taken as a deformation function ofOr exp (μ)ik),Φik(Xiik) Denotes the ith input variable XiAnd membership mu of its fuzzy group kikThe corresponding new input matrix;
the radial basis function neural network is used as a local equation of the fuzzy equation system, and optimization fitting is carried out on each fuzzy group; let the output of the kth radial basis blurring system be:
<math> <mrow> <msub> <mover> <mi>y</mi> <mo>^</mo> </mover> <mrow> <mi>i</mi> <mi>k</mi> </mrow> </msub> <mo>=</mo> <munder> <mo>&Sigma;</mo> <mi>l</mi> </munder> <msub> <mi>w</mi> <mrow> <mi>l</mi> <mi>k</mi> </mrow> </msub> <msub> <mi>&Phi;</mi> <mrow> <mi>l</mi> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <mo>|</mo> <mo>|</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>C</mi> <mrow> <mi>l</mi> <mi>k</mi> </mrow> </msub> <mo>|</mo> <mo>|</mo> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein,is the output of the kth radial basis blur System, ClkAnd wlkIs the center and output weight, phi, of the ith node of the kth radial basis ambiguity systemlk(||Xi-Clk| |) is the radial basis function of the kth node of the kth radial basis fuzzy system, determined by the following equation:
<math> <mrow> <msub> <mi>&Phi;</mi> <mrow> <mi>l</mi> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <mo>|</mo> <mo>|</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>C</mi> <mrow> <mi>l</mi> <mi>k</mi> </mrow> </msub> <mo>|</mo> <mo>|</mo> <mo>)</mo> </mrow> <mo>=</mo> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mfrac> <msup> <mrow> <mo>(</mo> <mo>|</mo> <mo>|</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>C</mi> <mrow> <mi>l</mi> <mi>k</mi> </mrow> </msub> <mo>|</mo> <mo>|</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mrow> <mn>2</mn> <mo>&times;</mo> <msub> <mi>&sigma;</mi> <mrow> <mi>l</mi> <mi>k</mi> </mrow> </msub> </mrow> </mfrac> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein sigmalkIs the gaussian width of the corresponding radial basis function, | | | · | |, is a norm expression; and finally, obtaining the final prediction output of the fuzzy equation system by a gravity center method in an anti-fuzzy method:
<math> <mrow> <msub> <mover> <mi>y</mi> <mo>^</mo> </mover> <mi>i</mi> </msub> <mo>=</mo> <mfrac> <mrow> <msubsup> <mo>&Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <msup> <mi>c</mi> <mo>*</mo> </msup> </msubsup> <mrow> <msub> <mi>&mu;</mi> <mrow> <mi>i</mi> <mi>k</mi> </mrow> </msub> <msub> <mover> <mi>y</mi> <mo>^</mo> </mover> <mrow> <mi>i</mi> <mi>k</mi> </mrow> </msub> </mrow> </mrow> <mrow> <msubsup> <mo>&Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <msup> <mi>c</mi> <mo>*</mo> </msup> </msubsup> <msub> <mi>&mu;</mi> <mrow> <mi>i</mi> <mi>k</mi> </mrow> </msub> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow> </math>
in the formula,is the predicted output of the fuzzy equation system;
an intelligent optimization module for adopting particle swarm optimization to perform C on the radial basis function neural network local equation in the fuzzy equationpk、σpk、wpkThe optimization is carried out, and the specific steps are as follows:
firstly, determining optimized parameters of particle swarm as C of radial basis function neural network local equationlk、σlk、wlkNumber of individual particle groups popsize, maximum number of cyclic optimization itermaxInitial position r of the p-th particlepInitial velocity vpLocal optimum value LbestpAnd the global optimum value Gbest of the whole particle swarm;
setting an optimization objective function, converting the optimization objective function into fitness, and evaluating each local fuzzy equation; calculating a fitness function through a corresponding error function, considering that the fitness of the particle with large error is small, and expressing the fitness function of the particle p as follows:
fp=1/(Ep+1) (9)
in the formula, EpIs the error function of the fuzzy equation system, expressed as:
<math> <mrow> <msub> <mi>E</mi> <mi>p</mi> </msub> <mo>=</mo> <msqrt> <mrow> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msup> <mrow> <mo>(</mo> <msub> <mover> <mi>y</mi> <mo>^</mo> </mover> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>O</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow> </math>
in the formula,is the predicted output of a system of fuzzy equations, OiThe target output of the fuzzy equation system is shown, and N is the number of training samples;
thirdly, according to the following formula, the speed and the position of each particle are circularly updated,
vp(iter+1)=ω×vp(iter)+m1a1(Lbestp-rp(iter))+m2a2(Gbest-rp(iter))
(11)
rp(iter+1)=rp(iter)+vp(iter+1) (12)
in the formula, vpIndicates the velocity, r, of the update particle ppIndicating the position of the update particle p, LbestpRepresenting the individual optimum of the update particle p, iter representing the number of cycles, ω being the inertial weight in the particle swarm optimization, m1、m2Is the corresponding acceleration factor, a1、a2Is [0,1 ]]A random number in between;
and fourthly, for the particle p, if the new fitness is larger than the original individual optimal value, updating the individual optimal value of the particle:
Lbestp=fp (13)
if the individual optimum value Lbest of the particle ppAnd (3) updating the original particle swarm global optimum value Gbest if the particle swarm global optimum value Gbest is larger than the original particle swarm global optimum value Gbest:
Gbest=Lbestp (14)
sixthly, judging whether the performance requirement is met, if so, finishing the optimization to obtain a set of optimized local equation parameters of the fuzzy equation; otherwise, returning to the step III, continuing iterative optimization until reaching the maximum cyclic optimization times itermax
The host computer still include:
the model updating module is used for acquiring field intelligent instrument signals according to a set sampling time interval, comparing the obtained actually-measured furnace temperature with a system forecast value, and if the relative error is more than 10% or the furnace temperature exceeds the upper and lower production normal limit ranges, adding new data which enables the furnace temperature to be optimal when the furnace temperature is normally produced in the DCS database into training sample data, and updating the soft measurement model;
the result display module is used for transmitting the obtained furnace temperature forecast value and the operation variable value which enables the furnace temperature to be optimal to the DCS, displaying the values at a control station of the DCS and transmitting the values to a field operation station for displaying through the DCS and a field bus;
the signal acquisition module is used for acquiring data from the database according to the set time interval of each sampling;
the key variables include the flow of waste liquid into the incinerator, the flow of air into the incinerator and the flow of fuel into the incinerator; the manipulated variables include air flow into the incinerator and fuel flow into the incinerator.
2. An intelligent radial-based pesticide production waste liquid incinerator temperature optimization method is characterized by comprising the following steps: the furnace temperature optimization method comprises the following specific implementation steps:
1) determining used key variables, collecting data of the variables in normal production from a DCS (distributed control system) database as an input matrix of a training sample TX, and collecting corresponding furnace temperature and operation variable data for optimizing the furnace temperature as an output matrix O;
2) preprocessing a model training sample input from a DCS database, centralizing the training sample, namely subtracting the average value of the sample, and then standardizing the training sample to ensure that the average value is 0 and the variance is 1; the processing is accomplished using the following mathematical process:
2.1) calculating the mean value: <math> <mrow> <mover> <mrow> <mi>T</mi> <mi>X</mi> </mrow> <mo>&OverBar;</mo> </mover> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>TX</mi> <mi>i</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </math>
2.2) calculating the variance: <math> <mrow> <msubsup> <mi>&sigma;</mi> <mi>x</mi> <mn>2</mn> </msubsup> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </mfrac> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <mrow> <mo>(</mo> <msub> <mi>TX</mi> <mi>i</mi> </msub> <mo>-</mo> <mover> <mrow> <mi>T</mi> <mi>X</mi> </mrow> <mo>&OverBar;</mo> </mover> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow> </math>
2.3) standardization: <math> <mrow> <mi>X</mi> <mo>=</mo> <mfrac> <mrow> <mi>T</mi> <mi>X</mi> <mo>-</mo> <mover> <mrow> <mi>T</mi> <mi>X</mi> </mrow> <mo>&OverBar;</mo> </mover> </mrow> <msub> <mi>&sigma;</mi> <mi>x</mi> </msub> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein, TXiThe ith training sample is data of key variables, furnace temperature and operating variables for optimizing the furnace temperature during normal production collected from the DCS database, N is the number of training samples,is the mean value of the training sample, and X is the training sample after standardization; sigmaxRepresenting the standard deviation, σ, of the training samples2 xRepresenting the variance of the training samples;
3) fuzzification is carried out on the training samples transmitted from the standardization processing module; let c in fuzzy system*A fuzzy group, the centers of the fuzzy groups k and j are vk、vjThe ith normalized training sample XiMembership mu for fuzzy group kikComprises the following steps:
<math> <mrow> <msub> <mi>&mu;</mi> <mrow> <mi>i</mi> <mi>k</mi> </mrow> </msub> <mo>=</mo> <msup> <mrow> <mo>(</mo> <munderover> <mo>&Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mrow> <mi>c</mi> <mo>*</mo> </mrow> </munderover> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>v</mi> <mi>k</mi> </msub> <mo>|</mo> <mo>|</mo> </mrow> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>&nu;</mi> <mi>j</mi> </msub> <mo>|</mo> <mo>|</mo> </mrow> </mfrac> <mo>)</mo> </mrow> <mfrac> <mn>2</mn> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </mfrac> </msup> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow> </math>
in the formula, n is a blocking matrix index required in the fuzzy classification process, and is taken as 2, | | · | | is a norm expression;
using the above membership values or its variants to obtain a new input matrix, for the fuzzy group k, its input matrix variant is:
Φik(Xiik)=[1 func(μik) Xi] (5)
wherein func (. mu.)ik) Is a membership value muikIs taken as a deformation function ofOr exp (μ)ik),Φik(Xiik) Denotes the ith input variable XiAnd membership mu of its fuzzy group kikThe corresponding new input matrix;
the radial basis function neural network is used as a local equation of the fuzzy equation system, and optimization fitting is carried out on each fuzzy group; let the output of the kth radial basis blurring system be:
<math> <mrow> <msub> <mover> <mi>y</mi> <mo>^</mo> </mover> <mrow> <mi>i</mi> <mi>k</mi> </mrow> </msub> <mo>=</mo> <munder> <mo>&Sigma;</mo> <mi>l</mi> </munder> <msub> <mi>w</mi> <mrow> <mi>l</mi> <mi>k</mi> </mrow> </msub> <msub> <mi>&Phi;</mi> <mrow> <mi>l</mi> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <mo>|</mo> <mo>|</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>C</mi> <mrow> <mi>l</mi> <mi>k</mi> </mrow> </msub> <mo>|</mo> <mo>|</mo> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein,is the output of the kth radial basis blur System, ClkAnd wlkIs the center and output weight, phi, of the ith node of the kth radial basis ambiguity systemlk(||Xi-Clk| |) is the radial basis function of the kth node of the kth radial basis fuzzy system, determined by the following equation:
<math> <mrow> <msub> <mi>&Phi;</mi> <mrow> <mi>l</mi> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <mo>|</mo> <mo>|</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>C</mi> <mrow> <mi>l</mi> <mi>k</mi> </mrow> </msub> <mo>|</mo> <mo>|</mo> <mo>)</mo> </mrow> <mo>=</mo> <mi>exp</mi> <mrow> <mo>(</mo> <mo>-</mo> <mfrac> <msup> <mrow> <mo>(</mo> <mo>|</mo> <mo>|</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>C</mi> <mrow> <mi>l</mi> <mi>k</mi> </mrow> </msub> <mo>|</mo> <mo>|</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mrow> <mn>2</mn> <mo>&times;</mo> <msub> <mi>&sigma;</mi> <mrow> <mi>l</mi> <mi>k</mi> </mrow> </msub> </mrow> </mfrac> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow> </math>
wherein sigmalkIs the gaussian width of the corresponding radial basis function, | | | · | |, is a norm expression; and finally, obtaining the final prediction output of the fuzzy equation system by a gravity center method in an anti-fuzzy method:
<math> <mrow> <msub> <mover> <mi>y</mi> <mo>^</mo> </mover> <mi>i</mi> </msub> <mo>=</mo> <mfrac> <mrow> <msubsup> <mi>&Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <msup> <mi>c</mi> <mo>*</mo> </msup> </msubsup> <msub> <mi>&mu;</mi> <mrow> <mi>i</mi> <mi>k</mi> </mrow> </msub> <msub> <mover> <mi>y</mi> <mo>^</mo> </mover> <mrow> <mi>i</mi> <mi>k</mi> </mrow> </msub> </mrow> <mrow> <msubsup> <mi>&Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <msup> <mi>c</mi> <mo>*</mo> </msup> </msubsup> <msub> <mi>&mu;</mi> <mrow> <mi>i</mi> <mi>k</mi> </mrow> </msub> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow> </math>
in the formula,is the predicted output of the fuzzy equation system;
4) c for radial basis function neural network local equation in fuzzy equation by adopting particle swarm optimizationpk、σpk、wpkThe optimization is carried out, and the specific steps are as follows:
firstly, determining optimized parameters of particle swarm as C of radial basis function neural network local equationlk、σlk、wlkNumber of individual particle groups popsize, maximum number of cyclic optimization itermaxInitial position r of the p-th particlepInitial velocity vpLocal optimum value LbestpAnd the global optimum value Gbest of the whole particle swarm;
setting an optimization objective function, converting the optimization objective function into fitness, and evaluating each local fuzzy equation; calculating a fitness function through a corresponding error function, considering that the fitness of the particle with large error is small, and expressing the fitness function of the particle p as follows:
fp=1/(Ep+1) (9)
in the formula, EpIs a fuzzy squareThe error function of the equation system, expressed as:
<math> <mrow> <msub> <mi>E</mi> <mi>p</mi> </msub> <mo>=</mo> <msqrt> <mrow> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <munderover> <mo>&Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msup> <mrow> <mo>(</mo> <msub> <mover> <mi>y</mi> <mo>^</mo> </mover> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>O</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow> </math>
in the formula,is the predicted output of a system of fuzzy equations, OiThe target output of the fuzzy equation system is shown, and N is the number of training samples;
thirdly, according to the following formula, the speed and the position of each particle are circularly updated,
vp(iter+1)=ω×vp(iter)+m1a1(Lbestp-rp(iter))+m2a2(Gbest-rp(iter))
(11)
rp(iter+1)=rp(iter)+vp(iter+1) (12)
in the formula, vpIndicates the velocity, r, of the update particle ppIndicating the position of the update particle p, LbestpRepresenting the individual optimum of the update particle p, iter representing the number of cycles, ω being the inertial weight in the particle swarm optimization, m1、m2Is the corresponding acceleration factor, a1、a2Is [0,1 ]]A random number in between;
and fourthly, for the particle p, if the new fitness is larger than the original individual optimal value, updating the individual optimal value of the particle:
Lbestp=fp (13)
if the individual optimum value Lbest of the particle ppAnd (3) updating the original particle swarm global optimum value Gbest if the particle swarm global optimum value Gbest is larger than the original particle swarm global optimum value Gbest:
Gbest=Lbestp (14)
sixthly, judging whether the performance requirement is met, if so, finishing the optimization to obtain a set of optimized local equation parameters of the fuzzy equation; otherwise, returning to the step III, continuing iterative optimization until reaching the maximum cyclic optimization times itermax
The method further comprises the following steps:
5) acquiring field intelligent instrument signals according to a set sampling time interval, comparing the obtained actually-measured furnace temperature with a system forecast value, and if the relative error is more than 10% or the furnace temperature exceeds the upper and lower production normal limit ranges, adding new data which enables the furnace temperature to be optimal when the furnace temperature is normally produced in a DCS database into training sample data, and updating a soft measurement model;
6) calculating to obtain an optimal operating variable value in the step 4), transmitting the obtained furnace temperature forecast value and the operating variable value which enables the furnace temperature to be optimal to the DCS, displaying the values at a control station of the DCS, and transmitting the values to a field operating station through the DCS and a field bus for displaying; meanwhile, the DCS system takes the obtained operating variable value which enables the furnace temperature to be optimal as a new operating variable set value, and automatically executes furnace temperature optimization operation;
the key variables include the flow of waste liquid into the incinerator, the flow of air into the incinerator and the flow of fuel into the incinerator; the manipulated variables include air flow into the incinerator and fuel flow into the incinerator.
CN201310436865.8A 2013-09-22 2013-09-22 The optimizing temperature of pesticide production waste liquid incinerator system and method for intelligence radial basis Expired - Fee Related CN103488206B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310436865.8A CN103488206B (en) 2013-09-22 2013-09-22 The optimizing temperature of pesticide production waste liquid incinerator system and method for intelligence radial basis

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310436865.8A CN103488206B (en) 2013-09-22 2013-09-22 The optimizing temperature of pesticide production waste liquid incinerator system and method for intelligence radial basis

Publications (2)

Publication Number Publication Date
CN103488206A CN103488206A (en) 2014-01-01
CN103488206B true CN103488206B (en) 2015-09-09

Family

ID=49828512

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310436865.8A Expired - Fee Related CN103488206B (en) 2013-09-22 2013-09-22 The optimizing temperature of pesticide production waste liquid incinerator system and method for intelligence radial basis

Country Status (1)

Country Link
CN (1) CN103488206B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106931453B (en) * 2017-02-27 2018-02-16 浙江大学 The forecasting system and method for circulating fluid bed domestic garbage burning emission of NOx of boiler

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6507774B1 (en) * 1999-08-24 2003-01-14 The University Of Chicago Intelligent emissions controller for substance injection in the post-primary combustion zone of fossil-fired boilers
CN101457264A (en) * 2008-12-29 2009-06-17 杭州电子科技大学 Blast furnace temperature optimization control method
CN101763085A (en) * 2009-12-29 2010-06-30 浙江大学 System and method for optimizing temperature of pesticide production waste liquid incinerator
CN101763084A (en) * 2009-12-29 2010-06-30 浙江大学 System and method for minimizing chemical oxygen demand (COD) discharge of pesticide production waste liquid incinerator
CN102176221A (en) * 2011-03-16 2011-09-07 中南大学 Coke furnace temperature predicting method based on dynamic working conditions in coke furnace heating and burning process

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040143149A1 (en) * 2002-08-02 2004-07-22 Decourcy Michael Stanley Method for reducing nitrogen oxide emissions in industrial process

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6507774B1 (en) * 1999-08-24 2003-01-14 The University Of Chicago Intelligent emissions controller for substance injection in the post-primary combustion zone of fossil-fired boilers
CN101457264A (en) * 2008-12-29 2009-06-17 杭州电子科技大学 Blast furnace temperature optimization control method
CN101763085A (en) * 2009-12-29 2010-06-30 浙江大学 System and method for optimizing temperature of pesticide production waste liquid incinerator
CN101763084A (en) * 2009-12-29 2010-06-30 浙江大学 System and method for minimizing chemical oxygen demand (COD) discharge of pesticide production waste liquid incinerator
CN102176221A (en) * 2011-03-16 2011-09-07 中南大学 Coke furnace temperature predicting method based on dynamic working conditions in coke furnace heating and burning process

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
农药废液焚烧建模与优化研究;刘兴高等;《仪器仪表学报》;20090630;第30卷(第6期);第266-269页 *
有机废液焚烧炉控制系统设计;潘华丰等;《控制工程》;20090731;第16卷;第53-56页 *

Also Published As

Publication number Publication date
CN103488206A (en) 2014-01-01

Similar Documents

Publication Publication Date Title
CN103472865B (en) The pesticide waste liquid incinerator furnace temperature optimization system of intelligence least square and method
Cheng et al. A convolutional neural network based degradation indicator construction and health prognosis using bidirectional long short-term memory network for rolling bearings
CN115169479A (en) Remote monitoring method, system and storage medium for sewage treatment process
CN101763085B (en) System and method for optimizing temperature of pesticide production waste liquid incinerator
CN103226741A (en) Urban water supply network tube explosion prediction method
CN106199174A (en) Extruder energy consumption predicting abnormality method based on transfer learning
CN111814956A (en) Multi-task learning air quality prediction method based on multi-dimensional secondary feature extraction
CN103472866B (en) The pesticide waste liquid incinerator furnace temperature optimization system of intelligent fuzzy system and method
He et al. Fault prediction method for wireless sensor network based on evidential reasoning and belief-rule-base
Zhang et al. A behavior prediction method for complex system based on belief rule base with structural adaptive
CN105183659A (en) Software system behavior anomaly detection method based on multi-level mode predication
CN118228606A (en) Online quick pollution prediction method for polluted land block based on machine learning, electronic equipment and storage medium
CN103472867B (en) The optimizing temperature of pesticide production waste liquid incinerator system and method for support vector machine
CN103488206B (en) The optimizing temperature of pesticide production waste liquid incinerator system and method for intelligence radial basis
Yuan et al. Prediction of PM2. 5 time series by seasonal trend decomposition-based dendritic neuron model
CN116861256A (en) Furnace temperature prediction method, system, equipment and medium for solid waste incineration process
CN103472721B (en) The pesticide waste liquid incinerator furnace temperature optimization system of self-adaptation machine learning and method
CN103488209B (en) The pesticide waste liquid incinerator furnace temperature optimization system of Intelligent Support vector machine and method
Xu et al. Prediction of the Wastewater's pH Based on Deep Learning Incorporating Sliding Windows.
CN103675012B (en) The industrial melt index soft measurement instrument of BP particle group optimizing and method
CN103472727B (en) The pesticide incinerator hazardous emission controls up to par system and method for gunz weighting
CN103488208B (en) The optimizing temperature of pesticide production waste liquid incinerator system and method for least square
CN103472729A (en) Crowd-sourcing system and method for controlling harmful emissions of pesticide waste liquid incinerator to reach standard
CN103472728B (en) The pesticide incinerator hazardous emission controls up to par system and method for machine learning
CN103488145B (en) The incinerator hazardous emission controls up to par system and method for gunz FUZZY NETWORK

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20150909

Termination date: 20180922