CN110703599A - Organic Rankine cycle system control performance optimization method based on dynamic data correction - Google Patents
Organic Rankine cycle system control performance optimization method based on dynamic data correction Download PDFInfo
- Publication number
- CN110703599A CN110703599A CN201910877359.XA CN201910877359A CN110703599A CN 110703599 A CN110703599 A CN 110703599A CN 201910877359 A CN201910877359 A CN 201910877359A CN 110703599 A CN110703599 A CN 110703599A
- Authority
- CN
- China
- Prior art keywords
- data
- output
- dynamic data
- entropy
- control performance
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
Abstract
The invention discloses an organic Rankine cycle system control performance optimization method based on dynamic data correction, which comprises the following steps of: the first stage, dynamic data correction: describing a system model under non-Gaussian disturbance; solving a likelihood function of the correction data according to a maximum posterior distribution principle; supplementing complete data by utilizing a maximum expectation algorithm idea to obtain a likelihood function of the complete data, and taking a logarithm of the likelihood function; finding a tight lower bound of the logarithm of the likelihood function to obtain an optimal estimation of the unknown parameter; and a second stage, optimizing the performance based on the minimum entropy control of dynamic data correction: selecting performance indexes based on statistical information; solving for control input; and constructing a performance evaluation standard. The invention considers the defect that the Renyi entropy does not meet the consistency, provides the minimum rational entropy control performance optimization method based on dynamic data correction, and can obviously improve the control performance of the organic Rankine cycle system containing non-Gaussian noise.
Description
Technical Field
The invention discloses an organic Rankine cycle system control performance optimization method based on dynamic data correction, and belongs to the technical field of organic Rankine cycle systems.
Background
An Organic Rankine Cycle (ORC) waste heat utilization process is a thermodynamic Cycle process that converts waste heat into electrical energy. The waste heat utilization has important significance for improving the energy utilization efficiency, reducing the industrial operation cost and reducing the pollution emission. The superheated steam temperature is one of the most important parameters in an organic rankine cycle process throughout the thermodynamic description of the organic rankine cycle. In order to improve the energy conversion efficiency of the organic rankine cycle process, the temperature of the organic working medium at the outlet of the evaporator cannot be too high from the aspects of safety and economy, and the temperature cannot be too low in order to avoid some critical conditions, namely, the temperature at the outlet of the evaporator should be controlled within a reasonable range. The temperature sensor is used for accurately measuring the outlet temperature of the evaporator and effectively controlling the outlet temperature, so that the organic Rankine cycle energy conversion efficiency is improved. If a bad control loop is not discovered in time, the energy conversion efficiency and the service life of the equipment are affected. Therefore, Control Performance Assessment (CPA) is very important in the industry.
Because factors such as aging and unstable voltage inevitably exist in the use process of the sensor, the measured data of the sensor are inaccurate, and the control effect is poor. The influence brought by measurement noise is not considered in the conventional control strategy and control performance evaluation method applied to the organic Rankine cycle waste heat utilization process. Dynamic data correction (DDR) can improve the quality of measurement data in a transient state, so that the measurement data including measurement noise is dynamically corrected to operate in an optimal state, which is significant for improving the economy and stability of the organic rankine cycle system.
At present, the dynamic data correction method mainly comprises a Kalman filtering method, an extended Kalman filtering method and an unscented Kalman filtering method, the methods are derived under the assumption of Gauss, and when the measurement noise distribution does not meet the assumption of Gauss, the data correction method cannot obtain a satisfactory result. The patent of non-gaussian system dynamic data correction and system control performance optimization (patent application number: 201811199510.0) is published, and combines a dynamic data correction method with a minimum Renyi entropy control strategy to provide a control performance optimization method, but does not consider the importance of control system performance evaluation, and the Renyi entropy in the patent does not meet the requirement of consistency, so that a new non-gaussian system control performance optimization and evaluation method is needed.
Disclosure of Invention
In order to overcome the defects of the prior art, the invention provides the organic Rankine cycle system control performance optimization method based on dynamic data correction, and the method can obviously improve the control performance of the organic Rankine cycle system containing non-Gaussian noise.
The invention is realized by the following technical scheme:
the organic Rankine cycle system control performance optimization method based on dynamic data correction comprises the following steps:
the first stage, dynamic data correction:
step 1.1, describing a system model under non-Gaussian disturbance;
step 1.2, solving correction data y according to the maximum posterior distribution principler(t) a likelihood function;
step 1.3, supplementing complete data by utilizing the thought of a maximum expectation algorithm to obtain a likelihood function of the complete data, and taking a logarithm of the likelihood function;
step 1.4. find the logarithm of the likelihood functionObtaining the optimal estimation of the unknown parameter theta by the tight lower bound of the parameter theta;
and a second stage, optimizing the performance based on the minimum entropy control of dynamic data correction:
2.1, selecting performance indexes based on statistical information;
step 2.1, solving a control input u (t);
and 2.3, constructing a performance evaluation standard.
Preferably, step 1.1 is obtained by:
in a non-linear system, the output at time t is decomposed into two parts: prediction outputAnd an unpredictable output δ (t), the process output y (t) being expressed as:
wherein u (t) is a control input, d (t) is random process interference, and f (-) and g (-) are arbitrary nonlinear relations which do not affect each other; true measurement output y due to the presence of measurement noise ε (t)m(t) is:
ym(t)=y(t)+ε(t) (2)
probability density functions of non-gaussian features have arbitrary shapes that are difficult to directly represent and process. As is well known, any probability density function can be approximated by a mixture of a finite number of gaussian probability density functions, and therefore, using a gaussian mixture model to derive a general dynamic data correction method, the probability density functions for the measured noise, epsilon (t) and delta (t), are as follows:
wherein K and N are the number of Gaussian mixture models, lambdakAnd munRespectively the mean value of each gaussian model,andvariance, a, of each Gaussian mixture modelkAnd bnAre weighted and satisfy φ1kProbability density function, phi, representing the kth Gaussian fraction model of epsilon (t)2nA probability density function representing the nth gaussian component model of delta (t),θ1and theta2All unknown parameters in the probability density function representing ε (t) and δ (t), respectively, i.e. θ1=(a1,a2,…,aK;θ11,θ12,…,θ1K),θ2=(b1,b2,…,bN;θ21,θ22,…,θ2N)。
Preferably, step 1.2 is obtained by:
measurement output ym(t) likelihood function L [ y [)m(t)|yr(t),θ2]And predicted outputLikelihood function ofIs represented as follows:
wherein, yr(t) denotes the correction output, and all unknown parameters are denoted as θ ═ yr,θ1,θ2) Output y from the measurement according to Bayesian criterionm(t) and prediction outputGiven of yrThe posterior probability of (t) is proportional to the product of the likelihood functions of the measured output and the predicted output, i.e.
According to the maximum a posteriori distribution principle, yr(t) is obtained by the following formula:
preferably, step 1.3 is obtained by:
to obtain the unknown parameter θ, a set of dataObviously, the method of Monte Carlo is adopted to generate J group dataThe estimation of the optimal parameters is obtained by:
let observation data ym(ym1,ym2,Λ,ymJ) Is generated as follows:
first according to the probability akSelecting the kth Gaussian distribution partial model g (y)m|θk) (ii) a Then according to the probability distribution g (y) of the k-th partial modelm|θk) Generating data ym(ii) a At this time, observation data ymjJ ═ 1,2, Λ, J, known; observation data of reaction ymjData from the K-th partial model is unknown, K1, 2jkRepresenting observation data ymjData from the kth partial model, defined as follows:
data are obtained by the same methodCorresponding hidden variablen=1,2,...,N,Representing observed data of reactionData from the nth partial model is defined as follows:
to pairTaking the logarithm to obtain:
the process of solving for θ directly by equation (14) is still complicated if foundAnd then continually optimize this tight lower bound, an optimal estimate of the unknown parameter theta will be obtained.
Said step 1.4 findingThe optimal estimate of the unknown parameter θ is obtained by the tight lower bound of:
according to formula (15), obtaining
Order toThen
namely find outThe equal sign of the formula (20) is ensured to be established, and the lower bound is further ensured to be tight, then
Wherein, C1And C2Is a constant; due to Q (gamma)jk) Andare respectively gammajkAndis a distribution function of
Q(γjk) Andhas obtained, Q (gamma)jk) Represents gammajkThe probability of 1 is given as the probability of,representsThe probability of (d); define Q (γ) respectivelyjk) Andis estimated asAnd
So far, a function which is easy to solve for the unknown parameter theta is found, and the maximization is realizedThe iteration is carried out until the required precision is reached,
θ(i)representing the value of theta at the i-th iteration, theta(i)The unknown parameters contained in the data are passed through the pairsAnd respectively calculating partial derivatives, and respectively setting the partial derivatives to be 0 to obtain:
repeating the iterative equations (27) - (33) until convergence to obtain the corrected output yr。
Preferably, in the second stage, the minimum entropy control performance optimization method based on dynamic data correction is as follows:
step 2.1. Performance index selection based on statistical information
Defining the tracking error e (t) as: e (t) ═ r (t) -yr(t) since the process is affected by non-gaussian random disturbances, more general statistical information besides mean, variance is usedA measure of row randomness, i.e., entropy, replaces the variance in a gaussian system; considering that the Renyi entropy does not meet the requirement of consistency, the rational entropy has the advantages of the Renyi entropy and meets the requirement of consistency, and therefore the entropy measure selected by the method is the rational entropy.
wherein γ (e) is the probability density function of e (t);
in order to characterize tracking performance and interference rejection performance, the average value and entropy of the tracking error are minimized; furthermore, the control energy is also minimized; thus, the controller is obtained by minimizing the following performance indicators:
wherein R is1Weights, R, representing rational entropy2Weight representing mean value of error, R3A weight representing the input square;
step 2.1. solving for control input u (t):
obtaining an optimal control input u (t) by adopting a random gradient method:
wherein λ is the learning rate, and the iteration is stopped when the required precision is reached;
the method for constructing the performance evaluation criterion in the step 2.3 comprises the following steps:
the control performance evaluation criterion is defined as the ratio of the optimal control performance index to the current control performance index:
wherein the content of the first and second substances,is the entropy value of the system output in the minimum rational entropy control strategy,is the actual entropy value of the current controller, and the control performance evaluation reference etaMRECA scalar value, eta, is given for evaluating the performance of the currently used controllerMRECHas a value between 0 and 1, when ηMRECWhen the value of (A) is close to 0, the current control state of the system is poor; when etaMRECWhen the value of (b) is close to 1, it indicates that the current process has a good control effect.
Compared with the prior art, the invention has the following beneficial effects:
because the actual industrial process is inevitably influenced by random noise and is generally non-Gaussian random noise, the influence caused by the random noise and the non-Gaussian random noise is considered when the distribution of the measurement noise and the process noise is determined, and compared with the existing method in which the influence of the noise is ignored or the noise is assumed to obey Gaussian distribution, the method has more generality and practical significance. In addition, the invention considers the defect that the Renyi entropy does not meet the consistency, provides the minimum rational entropy control performance optimization method based on dynamic data correction, is suitable for all non-Gaussian nonlinear processes, has a larger application value in practice, and has a larger significance for improving the safety and the stability of a control system.
Drawings
FIG. 1 is a structural diagram of temperature control of an organic Rankine cycle system according to the present invention.
FIG. 2 is a schematic diagram of a conventional organic Rankine cycle system.
FIG. 3 is a graph of the process noise and measurement noise distribution of the present invention.
FIG. 4 is a graph comparing control results according to the present invention.
Detailed Description
The present invention will be described in further detail with reference to specific examples, but the scope of the present invention is not limited to these examples, and all changes or equivalent substitutions that do not depart from the spirit of the present invention are intended to be included within the scope of the present invention.
Taking an input and output model of a superheated steam temperature system in an organic Rankine cycle process as an example, the ORC system mainly comprises the following components: working medium pump, evaporator, expander, condenser. As shown in the attached drawing 2, in the thermodynamic cycle process, a low-temperature low-pressure liquid organic working medium R245fa is pressurized by a working medium pump and enters an evaporator, in the evaporator, the working medium exchanges heat with a low-temperature waste heat source for evaporation, the working medium is heated by the heat source into superheated steam with higher temperature and high pressure, and then the steam enters an expander and expands in the expander to work to drive a generator to generate power. The low-temperature low-pressure organic working medium (namely exhaust gas) after acting enters a condenser, the organic working medium is condensed into a liquid state by releasing heat to a low-temperature heat source, and finally, the working medium pump conveys the liquid working medium back to the evaporator again to complete the whole cycle.
The rotating speed of a working medium pump in the model is used as control input, namely u (t), the temperature of the superheated steam in the ORC process is used as control output, namely y (t), and d (t) represents non-Gaussian random disturbance in the control process, namely non-Gaussian random disturbance of inlet temperature and smoke mass flow.
The organic Rankine cycle system temperature control method based on dynamic data correction comprises two stages of performance optimization of dynamic data correction and minimum entropy control performance evaluation based on the dynamic data correction, wherein the dynamic data correction performance optimization stage comprises the following steps:
the method comprises the following steps: and (4) describing a system model under non-Gaussian disturbance.
In an organic rankine cycle process superheated steam temperature system, the output at time t can be broken down into two parts: predictable portionAnd an unpredictable portion δ (t), the true measurement output, i.e. the measurement y of the superheated steam temperaturem(t) is:
ym(t)=y(t)+ε(t) (37)
ε (t) is the temperature sensor measurement noise. The Probability Density Function (PDF) of the measurement noise ∈ (t) and δ (t) is as follows:
the distribution diagram is shown in fig. 3, wherein K-N-3,θ1and theta2All unknown parameters in the probability density function representing ε (t) and δ (t), respectively, i.e. θ1=(a1,a2,…,aK;θ11,θ12,…,θ1K),θ2=(b1,b2,…,bN;θ21,θ22,…,θ2N)。
Step two: obtaining solving correction data y according to maximum posterior distribution principlerLikelihood function of (t)
all unknown parameters are noted as θ ═ yr,θ1,θ2). According to Bayes' rule, by ym(t) andgiven of yrThe posterior probability of (t) is proportional to the product of the likelihood functions of the measured output and the predicted output, i.e.
According to the maximum a posteriori distribution principle, yr(t) can be obtained by the following formula:
obviously, the calculation amount of the above formula is too large for directly solving the unknown parameter θ ═ yr,θ1,θ2)。
Step three: full data is supplemented using the most desirable algorithm concept.
Generating J-group data by adopting Monte Carlo methodIn this example J is 200. The estimate of the optimal parameters can be obtained by:
imagine observation data ym(ym1,ym2,Λ,ymJ) Is generated as follows: first according to the probability akSelecting the kth Gaussian distribution partial model g (y)m|θk) (ii) a Then according to the probability distribution g (y) of the k-th partial modelm|θk) Generating data ym(ii) a At this time, observation data ymjJ ═ 1,2, Λ, J, known; observation data of reaction ymjData from the K-th partial model is unknown, K1, 2jkWhich is defined as follows:
data are obtained by the same methodRelative to each otherImplicit variable of responseN ═ 1,2,., N is defined as follows:
then it is determined that,
the process of solving for theta directly by the above equation is still complex if foundAnd then continually optimize this tight lower bound, an optimal estimate of the unknown parameter theta will be obtained.
So far, a function which is easy to solve for the unknown parameter theta is found, and the maximization is realizedThe iteration is carried out until the required precision is reached,
θ(i)representing the value of theta at the ith iteration. Theta(i)The unknown parameters contained in the data can be processed by the pairThe partial derivatives were calculated and made to be 0, respectively. The results were:
repeating the iterative equations (50) - (56) until convergence to obtain the corrected output yr。
The minimum entropy control performance optimization stage based on dynamic data correction comprises the following steps:
the method comprises the following steps: selecting performance indexes based on statistical information, which is specifically as follows:
defining a tracking error as e (t) ═ r (t) -yr(t), the entropy measure selected by the invention is rational entropyWherein γ (e) is the probability density function of e (t); in order to characterize tracking performance and interference rejection performance, the average value and entropy of the tracking error are minimized; furthermore, the control energy is also minimized; thus, the controller is obtained by minimizing the following performance indicators:
wherein R is1,R2And R3Are weighted and are respectively set as R1=0.06,R20.02 and R3=0.001。
Step two: solving for the control input u (t).
If a target set point for the superheated steam temperature is given, then minimizing the objective function allows control and regulation of the pump speed u (t). Obtaining an optimal control input u (t) by adopting a random gradient method:
where λ is the learning rate, and empirically, λ is 0.001 and the accuracy is set to 0.001.
Step three: and constructing a performance evaluation standard.
The control performance evaluation criterion is defined as the ratio of the optimal control performance index to the current control performance index, and is a convenient method for measuring relative performance, and is as follows:
whereinIs the entropy value of the system output in the minimum rational entropy control strategy,is the actual entropy value, η, of the current controllerMRECScalar values for evaluating the performance of the currently used controller are given. Obviously, etaMRECIs between 0 and 1. When approaching 0, it is indicated that the current control state of the system is very poor. If the value is close to 1, the control effect of the current process is good.
Table 1 and FIG. 4 show the performance optimization results of the minimum entropy control based on the dynamic data correction in this embodiment, and table 1 shows the entropy values of the tracking errors
Output of | Noiseless output | Measuring output | Corrected output of the invention |
Entropy of the entropy | 0.004 | 0.018 | 0.006 |
From table 1, it can be found that the present invention can significantly improve the control performance of the organic rankine cycle system containing non-gaussian noise, and the stabilized correction output is close to the case without measurement noise, with a good effect. Based on the proposed performance evaluation criterion etaMRECThe control performance of the PID controller in the organic rankine cycle system was evaluated. When the system control performance is good, the parameters of the PID controller are respectively kp=0.035、ki0.15 and k d0. When the system control performance is poor, the parameters of the PID controller are respectively kp=0.11、ki0.09 and kd=0。
TABLE 2 control Performance evaluation criteria based on dynamic data correction
Table 2 shows the minimum entropy control performance evaluation result based on dynamic data correction in this embodiment, and it can be found that the present invention can accurately evaluate the actual control effect of the organic rankine cycle system.
The present invention is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.
Claims (6)
1. The method for optimizing the control performance of the organic Rankine cycle system based on dynamic data correction is characterized by comprising the following steps of:
the first stage, dynamic data correction:
step 1.1, describing a system model under non-Gaussian disturbance;
step 1.2, solving correction data y according to the maximum posterior distribution principler(t) a likelihood function;
step 1.3, supplementing complete data by utilizing the thought of a maximum expectation algorithm to obtain a likelihood function of the complete data, and taking a logarithm of the likelihood function;
step 1.4. find the logarithm of the likelihood functionObtaining the optimal estimation of the unknown parameter theta by the tight lower bound of the parameter theta;
and a second stage, optimizing the performance based on the minimum entropy control of dynamic data correction:
2.1, selecting performance indexes based on statistical information;
step 2.1, solving a control input u (t);
and 2.3, constructing a performance evaluation standard.
2. The dynamic data correction based organic rankine cycle system control performance optimization method according to claim 1, characterized in that the step 1.1 is obtained by the following method:
in a non-linear system, the output at time t is decomposed into two parts: prediction outputAnd an unpredictable output δ (t), the process output y (t) being expressed as:
wherein u (t) is a control input, d (t) is random process interference, and f (-) and g (-) are arbitrary nonlinear relations which do not affect each other; true measurement output y due to the presence of measurement noise ε (t)m(t) is:
ym(t)=y(t)+ε(t) (2)
the probability density functions of the measured noise epsilon (t) and delta (t) are as follows:
wherein K and N are the number of Gaussian mixture models, lambdakAnd munRespectively the mean value of each gaussian model,andvariance, a, of each Gaussian mixture modelkAnd bnAre weighted and satisfy φ1kProbability density function, phi, representing the kth Gaussian fraction model of epsilon (t)2nA probability density function representing the nth gaussian component model of delta (t),θ1and theta2All unknown parameters in the probability density function representing ε (t) and δ (t), respectively, i.e. θ1=(a1,a2,…,aK;θ11,θ12,…,θ1K),θ2=(b1,b2,…,bN;θ21,θ22,…,θ2N)。
3. The dynamic data correction based organic rankine cycle system control performance optimization method according to claim 1, characterized in that the step 1.2 is obtained by the following method:
measurement output ym(t) likelihood function L [ y [)m(t)|yr(t),θ2]And predicted outputLikelihood function ofIs represented as follows:
wherein, yr(t) denotes the correction output, and all unknown parameters are denoted as θ ═ yr,θ1,θ2) Output y from the measurement according to Bayesian criterionm(t) and prediction outputGiven of yrThe posterior probability of (t) is proportional to the product of the likelihood functions of the measured output and the predicted output, i.e.
According to the maximum posterior distribution principle, y is obtainedr(t) is:
4. the dynamic data correction based organic rankine cycle system control performance optimization method according to claim 1, characterized in that the step 1.3 is obtained by the following method:
generating J-group data by adopting Monte Carlo methodThe estimation of the optimal parameters is obtained by:
let observation data ym(ym1,ym2,Λ,ymJ) Is generated as follows:
first according to the probability akSelecting the kth Gaussian distribution partial model g (y)m|θk) (ii) a Then according to the probability distribution g (y) of the k-th partial modelm|θk) Generating data ym(ii) a At this time, observation data ymjJ ═ 1,2, Λ, J, known; observation data of reaction ymjData from the K-th partial model is unknown, K1, 2jkRepresenting observation data ymjData from the kth partial model, defined as follows:
data are obtained by the same methodCorresponding hidden variable Representing observed data of reactionData from the nth partial model is defined as follows:
to pairTaking the logarithm to obtain:
5. the dynamic data correction based organic rankine cycle system control performance optimization method according to claim 1, characterized in that the step 1.4 is obtained by the following method:
an equal sign holds if and only if x is a constant;
order to
according to formula (15), obtaining
namely find outThe equal sign of the formula (20) is ensured to be established, and the lower bound is further ensured to be tight, then
Wherein, C1And C2Is a constant; due to Q (gamma)jk) Andare respectively gammajkAndis a distribution function of
Q(γjk) Andhas been obtained to represent respectively gammajk1 andthe probability of (d); define Q (γ) respectivelyjk) Andis estimated asAnd
So far, a function which is easy to solve for the unknown parameter theta is found, and the maximization is realizedThe iteration is carried out until the required precision is reached,
θ(i)representing the value of theta at the i-th iteration, theta(i)The unknown parameters contained in the data are passed through the pairsAnd respectively calculating partial derivatives, and respectively setting the partial derivatives to be 0 to obtain:
repeating the iterative equations (27) - (33) until convergence to obtain the corrected output yr。
6. The dynamic data correction-based organic Rankine cycle system control performance optimization method according to claim 1, wherein the second stage, minimum entropy control performance optimization based on dynamic data correction is obtained by:
step 2.1. Performance index selection based on statistical information
Defining the tracking error e (t) as: e (t) ═ r (t) -yr(t),
wherein γ (e) is the probability density function of e (t);
the controller is obtained by minimizing the following performance indicators:
wherein R is1Weights, R, representing rational entropy2Weight representing mean value of error, R3A weight representing the input square;
step 2.1. solving for control input u (t):
obtaining an optimal control input u (t) by adopting a random gradient method:
wherein λ is the learning rate, and the iteration is stopped when the required precision is reached;
the method for constructing the performance evaluation criterion in the step 2.3 comprises the following steps:
the control performance evaluation criterion is defined as the ratio of the optimal control performance index to the current control performance index:
wherein the content of the first and second substances,is the entropy value of the system output in the minimum rational entropy control strategy,is the actual entropy value of the current controller, and the control performance evaluation reference etaMRECA scalar value, eta, is given for evaluating the performance of the currently used controllerMRECHas a value between 0 and 1, when ηMRECWhen the value of (A) is close to 0, the current control state of the system is poor; when etaMRECWhen the value of (b) is close to 1, it indicates that the current process has a good control effect.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910877359.XA CN110703599B (en) | 2019-09-17 | 2019-09-17 | Organic Rankine cycle system control performance optimization method based on dynamic data correction |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910877359.XA CN110703599B (en) | 2019-09-17 | 2019-09-17 | Organic Rankine cycle system control performance optimization method based on dynamic data correction |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110703599A true CN110703599A (en) | 2020-01-17 |
CN110703599B CN110703599B (en) | 2022-06-07 |
Family
ID=69194983
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910877359.XA Active CN110703599B (en) | 2019-09-17 | 2019-09-17 | Organic Rankine cycle system control performance optimization method based on dynamic data correction |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110703599B (en) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103927412A (en) * | 2014-04-01 | 2014-07-16 | 浙江大学 | Real-time learning debutanizer soft measurement modeling method on basis of Gaussian mixture models |
CN104392136A (en) * | 2014-11-28 | 2015-03-04 | 东南大学 | High-precision data fusion method oriented to high-dynamic non-Gaussian-model robustness measurement |
CN107608939A (en) * | 2017-08-16 | 2018-01-19 | 北京师范大学 | TRMM precipitation data NO emissions reduction methods based on high resolution satellite remote sensing image |
CN107942682A (en) * | 2017-12-22 | 2018-04-20 | 太原理工大学 | The dynamic economy performance optimization of non-Gaussian filtering and control integrated design method |
CN108446529A (en) * | 2018-06-22 | 2018-08-24 | 太原理工大学 | Organic rankine cycle system fault detection method based on broad sense cross-entropy-DPCA algorithms |
CN108520107A (en) * | 2018-03-19 | 2018-09-11 | 山西大学 | System state estimation method based on maximum-likelihood criterion Robust Kalman Filter |
CN108920868A (en) * | 2018-07-23 | 2018-11-30 | 国网江西省电力有限公司电力科学研究院 | Load based on intelligent electric meter data recognizes model parameter correction method |
CN109002835A (en) * | 2018-06-19 | 2018-12-14 | 西安电子科技大学 | A kind of particle filter data correlation method based on maximum entropy fuzzy clustering |
CN109240085A (en) * | 2018-10-15 | 2019-01-18 | 太原理工大学 | Non-Gaussian filtering dynamic data rectification and system control performance optimization method |
CN109447272A (en) * | 2018-09-29 | 2019-03-08 | 西安交通大学 | A kind of extreme learning machine method based on center of maximum cross-correlation entropy criterion |
-
2019
- 2019-09-17 CN CN201910877359.XA patent/CN110703599B/en active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103927412A (en) * | 2014-04-01 | 2014-07-16 | 浙江大学 | Real-time learning debutanizer soft measurement modeling method on basis of Gaussian mixture models |
CN104392136A (en) * | 2014-11-28 | 2015-03-04 | 东南大学 | High-precision data fusion method oriented to high-dynamic non-Gaussian-model robustness measurement |
CN107608939A (en) * | 2017-08-16 | 2018-01-19 | 北京师范大学 | TRMM precipitation data NO emissions reduction methods based on high resolution satellite remote sensing image |
CN107942682A (en) * | 2017-12-22 | 2018-04-20 | 太原理工大学 | The dynamic economy performance optimization of non-Gaussian filtering and control integrated design method |
CN108520107A (en) * | 2018-03-19 | 2018-09-11 | 山西大学 | System state estimation method based on maximum-likelihood criterion Robust Kalman Filter |
CN109002835A (en) * | 2018-06-19 | 2018-12-14 | 西安电子科技大学 | A kind of particle filter data correlation method based on maximum entropy fuzzy clustering |
CN108446529A (en) * | 2018-06-22 | 2018-08-24 | 太原理工大学 | Organic rankine cycle system fault detection method based on broad sense cross-entropy-DPCA algorithms |
CN108920868A (en) * | 2018-07-23 | 2018-11-30 | 国网江西省电力有限公司电力科学研究院 | Load based on intelligent electric meter data recognizes model parameter correction method |
CN109447272A (en) * | 2018-09-29 | 2019-03-08 | 西安交通大学 | A kind of extreme learning machine method based on center of maximum cross-correlation entropy criterion |
CN109240085A (en) * | 2018-10-15 | 2019-01-18 | 太原理工大学 | Non-Gaussian filtering dynamic data rectification and system control performance optimization method |
Non-Patent Citations (3)
Title |
---|
ZHENGJIANG ZHANG等: "Dynamic data reconciliation for enhancing performance of minimum variance control in univariate and multivariate systems", 《INDUSTRIAL&ENGINEERING CHEMISTRY RESEARCH》 * |
成婷等: "基于中心误差熵准则的非高斯***滤波器设计", 《太原理工大学学报》 * |
赵雅兰等: "数据驱动框架下的非高斯批次过程最小熵性能评估算法", 《太原理工大学学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN110703599B (en) | 2022-06-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112085277B (en) | SCR denitration system prediction model optimization method based on machine learning | |
TWI714757B (en) | Eco-efficiency characterization tool | |
Prasad et al. | Plant-wide predictive control for a thermal power plant based on a physical plant model | |
CN112818595A (en) | Method and system for correcting digital twin model data of evaporation zone of thermal power plant | |
Shi et al. | Dual-mode fast DMC algorithm for the control of ORC based waste heat recovery system | |
CN115510904B (en) | Boiler heating area ash monitoring method based on time sequence prediction | |
Yu et al. | Hybrid modelling and simulation of thermal systems of in-service power plants for digital twin development | |
CN103235512B (en) | A kind of method that genset runs | |
CN108182553B (en) | Coal-fired boiler combustion efficiency online measurement method | |
CN115077640A (en) | Real-time prediction method for main steam flow in urban solid waste incineration process | |
Shajahan et al. | Improvement in efficiency of thermal power plant using optimization and robust controller | |
CN110703599B (en) | Organic Rankine cycle system control performance optimization method based on dynamic data correction | |
CN107247994B (en) | Fuzzy modeling method for desulfurization efficiency of tray tower desulfurization device | |
CN116865343B (en) | Model-free self-adaptive control method, device and medium for distributed photovoltaic power distribution network | |
CN115081327A (en) | Partition wall type heat exchanger dynamic simulation model modeling method based on reinforcement learning | |
Li et al. | Operation Data based Modelling and Optimization of Thermal Power Units under Full Working Conditions | |
Chen et al. | Optimization of turbine cold-end system based on BP neural network and genetic algorithm | |
CN112670997A (en) | Electric heating energy source system time sequence probability load flow calculation method considering photovoltaic uncertainty | |
Ou et al. | Performance analysis of a liquid desiccant air conditioning system based on a data-driven model | |
Zheng et al. | Improved prediction of nitrogen oxides using GRNN with k-means clustering and EDA | |
Guo et al. | Neural Network Modeling and Predictive Control of Low Oxygen Combustion System in Petrochemical Heating Furnace | |
Chen et al. | Dynamic modeling of biotechnical process based on online support vector machine | |
CN114320828B (en) | Heat accumulating type compressed air energy storage system and control method | |
Jitchaiyapoom et al. | Production Capacity Prediction and Optimization in the Glycerin Purification Process: A Simulation-Assisted Few-Shot Learning Approach | |
Guo et al. | Nox prediction of gas turbine based on dual attention and lstm |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |