CN107355767A - Control method and apparatus, the storage medium of separator temperature - Google Patents

Control method and apparatus, the storage medium of separator temperature Download PDF

Info

Publication number
CN107355767A
CN107355767A CN201710446637.7A CN201710446637A CN107355767A CN 107355767 A CN107355767 A CN 107355767A CN 201710446637 A CN201710446637 A CN 201710446637A CN 107355767 A CN107355767 A CN 107355767A
Authority
CN
China
Prior art keywords
mrow
msub
msup
fuel
variable
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.)
Pending
Application number
CN201710446637.7A
Other languages
Chinese (zh)
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.)
China Shenhua Energy Co Ltd
Beijing Guohua Electric Power Co Ltd
Zhejiang Guohua Zheneng Power Generation Co Ltd
Original Assignee
China Shenhua Energy Co Ltd
Beijing Guohua Electric Power Co Ltd
Zhejiang Guohua Zheneng Power Generation Co Ltd
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 China Shenhua Energy Co Ltd, Beijing Guohua Electric Power Co Ltd, Zhejiang Guohua Zheneng Power Generation Co Ltd filed Critical China Shenhua Energy Co Ltd
Priority to CN201710446637.7A priority Critical patent/CN107355767A/en
Publication of CN107355767A publication Critical patent/CN107355767A/en
Pending legal-status Critical Current

Links

Classifications

    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F22—STEAM GENERATION
    • F22B—METHODS OF STEAM GENERATION; STEAM BOILERS
    • F22B35/00—Control systems for steam boilers
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F22—STEAM GENERATION
    • F22B—METHODS OF STEAM GENERATION; STEAM BOILERS
    • F22B37/00—Component parts or details of steam boilers
    • F22B37/02—Component parts or details of steam boilers applicable to more than one kind or type of steam boiler
    • F22B37/26—Steam-separating arrangements

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Thermal Sciences (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Fuel Cell (AREA)

Abstract

The embodiment of the present invention provides a kind of method for controlling separator temperature, belongs to generator set control field.The method of the control separator temperature includes：Merge separator temperature variable and Fuel- Water Rate variable, to generate optimized variable；The optimized variable is brought into the optimal objective function of the separator temperature variable and the Fuel- Water Rate variable, to obtain quadratic programming problem；Derive the KKT conditions of the quadratic programming problem；The KKT conditions are solved, to obtain the optimal Fuel- Water Rate of the Fuel- Water Rate；And separator temperature is controlled according to the optimal Fuel- Water Rate.

Description

Control method and apparatus, the storage medium of separator temperature
Technical field
The present invention relates to generator set control field, method and apparatus, storage more particularly to control separator temperature Medium.
Background technology
In generating set AGC control systems, the response of boiler and turbine control loop determines different generating set productions The response difference of process control methods.The response that separator temperature inputs change to boiler oil more lags, and this is by fuel What a series of links such as the hysteresis of supply, heat interaction, heat transmission, the generation of steam, transmission were determined.In general, divide It is in Higher-order inertia link for fuel input response from device temperature, time constant is in or so a few minutes.Change of the steam turbine for pitch rings Should be very fast, only in low pressure stage, there is certain hysteresis.And existing Optimal Control System is still adopted in overall control structure With the control model of feedforward+feedback, and utilize traditional PID control.Therefore, steam turbine is easily controlled, and boiler implosion has necessarily Difficulty.
For the gen-set control system in thermoelectricity electricity generation system, tradition is coordinated to control most weak link to be, has There is the boiler combustion system of Higher-order inertia link, but to control main vapour pressure.Boiler heat load instruction by pulverized coal preparation system (several minutes Inertia time) heat energy is converted to, heating surface heat exchange, working medium transmission, last boiler boosting are passed through among this.There is high-order at this In inertia system, opened loop control can cause fluctuation, it is necessary to can just make system stable by the regulation of controller.But existing In AGC control systems, there is the problem of excessively hysteresis in separator temperature, thus be unfavorable for the overall stability of AGC control systems.
Fundamentally to solve the above problems, should be by advanced control technology such as：PREDICTIVE CONTROL, ANN Control, from The technologies such as suitable solution, fuzzy control are applied in the optimal control of generating set.Advanced AGC real-time optimal controls system is melted Closed a variety of control technologies state-of-the-art in the world, be exclusively for the difficulties solved in above-mentioned generating set AGC control and The advanced control platform of research and development.
The content of the invention
The purpose of the embodiment of the present invention is to provide a kind of method and apparatus, storage medium for controlling separator temperature, the party Method and device can improve the stability of separator temperature control and do disturbance ability.
To achieve these goals, embodiments of the invention provide a kind of method for controlling separator temperature, this method bag Include：Merge separator temperature variable and Fuel- Water Rate variable, to generate optimized variable；Bring the optimized variable into the separator In the optimal objective function of temperature variable and the Fuel- Water Rate variable, to obtain quadratic programming problem；Derive the quadratic programming The KKT conditions of problem；The KKT conditions are solved, to obtain the optimal Fuel- Water Rate of the Fuel- Water Rate；And according to the optimal combustion Water is than control separator temperature.
Wherein, it is described to solve the KKT conditions, it can be included with obtaining the optimal Fuel- Water Rate of the Fuel- Water Rate：By described in KKT conditions are converted into the KKT conditions with Min-Algebra type functions；With Min-Algebra type functions described in solving KKT conditions, to obtain optimal sequence；First value is taken from the optimal sequence as optimal solution；And according to described optimal Solution determines the optimal Fuel- Water Rate.
Wherein, the quadratic programming problem is：
Wherein,
X=[x (k+1) ' ... x (k+N) '] ', U=[u (k) ' ... x (k+N-1) '] ', Z=[U'X'] '
C=[0 ... 0, (- 2Qr) ' ... (- 2Qr) '] '
B=[umax′,umax′…umax′]′
Wherein, A, B are the coefficient of the Space admittance of the separator temperature variable and the Fuel- Water Rate variable, X tables Show the separator temperature matrix of variables in prediction time domain N, U represents the Fuel- Water Rate variable square in prediction time domain N Battle array, Z are optimized variable matrix, and b is bound variable, and Q is X diagonal positive definite matrix, and R is U diagonal positive definite matrix, and G is constraint The coefficient matrix of condition, I are unit matrix, and r is the setting value of separator temperature, and x (k) is the separator temperature variable, u (k) the Fuel- Water Rate variable, umaxFor the higher limit of controlled quentity controlled variable Fuel- Water Rate.
Wherein, the KKT conditions can be：
MZ+c+G′y+Aeq' γ=0
Wherein, γ, yiFor Lagrange coefficient, i is vector order.
Wherein, the Min-Algebra type functions are：
Φmin(y)=0
Wherein,
According to another aspect of the present invention, a kind of device for controlling separator temperature is also provided, the device includes：Secondary rule Problem generation module is drawn, for merging separator temperature variable and Fuel- Water Rate variable, to generate optimized variable, and by the optimization Variable is brought into the optimal objective function of the separator temperature variable and the Fuel- Water Rate variable, is asked with obtaining quadratic programming Topic；Optimal Fuel- Water Rate acquisition module, for deriving the KKT conditions of the quadratic programming problem, and the KKT conditions are solved, with Obtain the optimal Fuel- Water Rate of the Fuel- Water Rate；And control module, for controlling separator temperature according to the optimal Fuel- Water Rate.
Wherein, it is described to ask optimal Fuel- Water Rate solution module to include：KKT condition modular converters, for the KKT conditions to be turned Turn to the KKT conditions with Min-Algebra type functions；Module is solved, described Min-Algebra type functions are carried for solving KKT conditions, to obtain optimal sequence；Determining module, for determining the optimal Fuel- Water Rate according to the optimal solution.
Wherein, the quadratic programming problem is：
Wherein,
X=[x (k+1) ' ... x (k+N) '] ', U=[u (k) ' ... x (k+N-1) '] ', Z=[U'X'] '
C=[0 ... 0, (- 2Qr) ' ... (- 2Qr) '] '
B=[umax′,umax′…umax′]′
Wherein, A, B are the coefficient of the Space admittance of the separator temperature variable and the Fuel- Water Rate variable, X tables Show the separator temperature matrix of variables in prediction time domain N, U represents the Fuel- Water Rate variable square in prediction time domain N Battle array, Z are optimized variable matrix, and b is bound variable, and Q is X diagonal positive definite matrix, and R is U diagonal positive definite matrix, and G is constraint The coefficient matrix of condition, I are unit matrix, and r is the setting value of separator temperature, and x (k) is the separator temperature variable, u (k) it is the Fuel- Water Rate variable, umaxFor the higher limit of controlled quentity controlled variable Fuel- Water Rate.
Wherein, the KKT conditions are：
MZ+c+G′y+Aeq' γ=0
Wherein, γ, yiFor Lagrange coefficient, i is vector order.
Wherein, the Min-Algebra type functions are：
Φmin(y)=0
Wherein,
On the other hand, the present invention provides a kind of machinable medium, and finger is stored with the machinable medium Order, the instruction are used for the method for make it that machine performs control separator temperature described above.
Pass through above-mentioned technical proposal, being capable of look-ahead regulated variable using Prediction and Control Technology --- separator temperature Future trends, then it is controlled according to the change in future amount of regulated variable, effectively adjusts process in advance, so as to greatly improve The closed loop stability and Ability of Resisting Disturbance of unit AGC control systems.At the same time, for realistic model dimension feature, use Retain equality condition without removing by way of iteration erased condition variable to obtain the solution to decision variable so that matrix meter The complexity of calculation reduces, and the quadratic programming problem of equality constraint and inequality constraints is contained while obtaining by KKT conditions i.e. Equality condition can be removed and be integrated into Min-Algebra type object functions, realize the increase for no matter calculating dimension whether (during prediction Length of field) it can use relatively low iterations realization quickly to restrain and obtain optimal output, effective simplify control program, improve Control efficiency.
The further feature and advantage of the embodiment of the present invention will be described in detail in subsequent specific embodiment part.
Brief description of the drawings
Accompanying drawing is that the embodiment of the present invention is further understood for providing, and a part for constitution instruction, with The embodiment in face is used to explain the embodiment of the present invention together, but does not form the limitation to the embodiment of the present invention.Attached In figure：
Fig. 1 is the flow chart of the method for control separator temperature according to an embodiment of the invention；
Fig. 2 is the flow chart of the method for control separator temperature according to another embodiment of the present invention；
Fig. 3 is the structured flowchart of the device of control separator temperature according to another embodiment of the present invention；
Fig. 4 be control separator temperature according to another embodiment of the present invention device in optimal Fuel- Water Rate solve module Structured flowchart；And
Fig. 5 is that the method for the control separator temperature of the present invention is applied to a kind of feelings in generating set AGC control systems The control signal of shape is thought.
Description of reference numerals
310：Quadratic programming problem generation module 320：Optimal Fuel- Water Rate acquisition module
330：Control module 410：KKT condition modular converters
420：Solve module 430：Determining module
Embodiment
The embodiment of the embodiment of the present invention is described in detail below in conjunction with accompanying drawing.It should be appreciated that this The embodiment of place description is merely to illustrate and explain the present invention embodiment, is not intended to limit the invention embodiment.
Fig. 1 is the flow chart of the method for control separator temperature according to an embodiment of the invention.As shown in figure 1, the party Method comprises the following steps：
In step s 110, separator temperature variable and Fuel- Water Rate variable are merged, to generate optimized variable.
In the step s 120, the optimized variable is brought into the separator temperature variable and the Fuel- Water Rate variable most In excellent object function, to obtain quadratic programming problem.
In step s 130, the KKT conditions of the quadratic programming problem are derived.
In step S140, the KKT conditions are solved, to obtain the optimal Fuel- Water Rate of the Fuel- Water Rate.
Separator temperature is controlled in step S150, and according to the optimal Fuel- Water Rate.
Fig. 2 is the flow chart of the method for control separator temperature according to another embodiment of the present invention.As shown in Fig. 2 should Method may include following steps：
Wherein, step S210-S230 is with above-mentioned steps S110-S130, and here is omitted.
In step S240, the KKT conditions are converted into the KKT conditions with Min-Algebra type functions.
In step s 250, the KKT conditions for carrying Min-Algebra type functions are solved, to obtain optimal sequence.
In step S260, for determining the optimal Fuel- Water Rate according to the optimal solution.
In step S270, the optimal Fuel- Water Rate is determined according to the optimal solution.
The method of the optimal Fuel- Water Rate of solution in the preferred embodiments of the present invention described further below.In following content In, symbol " ' " represent that transposed matrix, such as U' represent U transposed matrix.
The quadratic programming problem can be represented by the mathematical expression one：
Mathematical expression 1：
Wherein, in above-mentioned mathematical expression 1,
X=[x (k+1) ' ... x (k+N) '] ', U=[u (k) ' ... x (k+N-1) '] ', Z=[U'X'] '
C=[0 ... 0, (- 2Qr) ' ... (- 2Qr) '] '
B=[umax′,umax′…umax′]′
Wherein, A, B are the coefficient of the Space admittance of the separator temperature variable and the Fuel- Water Rate variable, X tables Show the separator temperature matrix of variables in prediction time domain N, U represents the Fuel- Water Rate variable square in prediction time domain N Battle array, Z are optimized variable matrix, and b is bound variable, and Q is X diagonal positive definite matrix, and R is U diagonal positive definite matrix, and G is constraint The coefficient matrix of condition, I are unit matrix, and r is the setting value of separator temperature, and x (k) is the separator temperature variable, u (k) it is the Fuel- Water Rate variable, umaxFor the higher limit of controlled quentity controlled variable Fuel- Water Rate, i.e. constraints.
The Space admittance is represented by as shown in mathematical expression 2：
Mathematical expression 2：
X (k+1)=Ax (k)+Bu (k), u (k)≤umax
Wherein, A, B in mathematical expression 2 are the Space admittance of the separator temperature variable and the Fuel- Water Rate variable Coefficient.
The object function can be expressed as shown in mathematical expression 3：
Mathematical expression 3：
Wherein, xrefFor the desired value of separator temperature.The Space admittance and the object function are this areas Known routine, therefore no longer this is described in detail for this specification.
Above-mentioned mathematical expression 1 can be converted into KKT conditions, and the KKT conditions can be expressed as shown in mathematical expression 4：
Mathematical expression 4：
MZ+c+G′y+Aeq' γ=0
Wherein, γ, yiFor Lagrange coefficient, i is vector order.
In the KKT conditions described in mathematical expression 4, M and AeqAll it is sparse diagonal matrix, therefore calculates simple.Can be further Equation simultaneous in equality constraint and above-mentioned mathematical expression 4 is eliminated into Lagrange coefficient γ, so as to can obtain Z optimum control sequence Row, can be expressed as shown in mathematical expression 5：
Mathematical expression 5：
Z=- (DG'y+d)
Wherein,
Meanwhile the inequality constraints condition in mathematical expression 4 can be converted to as shown in mathematical expression 6.
Mathematical expression 6：
GZ-min (b, GZ+y)=0
Further mathematical expression 5 can be substituted into mathematical expression 6, it can thus be concluded that to the Min-Algebra type functions on y, should Min-Algebra type functions are represented by as shown in mathematical expression 7.
Mathematical expression 7：
Φmin(y)=0
Thus, the KKT conditions with Min-Algebra type functions, the KKT bars of band Min-Algebra type functions can be drawn Part is represented by as shown in mathematical expression 8：
Mathematical expression 8：
MZ+c+G′y+Aeq' γ=0
Φmin(y)=0
The KKT conditions of the band Min-Algebra type functions as shown in mathematical expression 8 after conversion calculate simple, this turn of solution KKT after change, which obtains part, can draw Lagrange multiplier y optimal sequence, because obtained optimal sequence is a pre- sequencing Row, so the value in the optimal sequence more rearward is due to simply predicting solution so deviation is bigger, therefore generally chooses optimal sequence First value be used as optimal solution.
After the optimal solution is substituted into mathematical expression 5, optimized variable Z optimal solution can must be obtained, because Z is on separator temperature The optimized variable of variable and Fuel- Water Rate variable, therefore, optimal Fuel- Water Rate can be from which further followed that.
And then separator temperature can be controlled according to the optimal Fuel- Water Rate in the DCS control systems of AGC control systems, Specific control process can be illustrated as shown in Figure 5.
Fig. 5 is that the method for the control separator temperature of the present invention is applied to a kind of feelings in generating set AGC control systems The control schematic diagram of shape.In Figure 5, Tsp0 represents separator temperature setting value, and Tsp represents the value of feedback of separator temperature, FWR Represent Fuel- Water Rate.Boiler sends boiler master instruction, and master control instruction generation confluent instruction is to control total confluent, wherein always Confluent can be known, and generate coal-supplying amount instruction, and control system is always given according to coal-supplying amount instruction and Fuel- Water Rate FWR controls Coal amount.Fuel- Water Rate is by separator temperature controller according to the optimal Fuel- Water Rate control that obtains in the method for the control separator temperature Fuel- Water Rate FWR processed output valve.
Fig. 3 is the structured flowchart of the device of control separator temperature according to another embodiment of the present invention.As shown in figure 3, The device includes：Quadratic programming problem generation module 310, for merging separator temperature variable and Fuel- Water Rate variable, with generation Optimized variable, and the optimized variable is brought into the optimal objective function of the separator temperature variable and the Fuel- Water Rate variable In, to obtain quadratic programming problem；Optimal Fuel- Water Rate acquisition module 320, for deriving the KKT bars of the quadratic programming problem Part, and the KKT conditions are solved, to obtain the optimal Fuel- Water Rate of the Fuel- Water Rate；And control module 330, for according to institute State optimal Fuel- Water Rate control separator temperature.
Fig. 4 be control separator temperature according to another embodiment of the present invention device in optimal Fuel- Water Rate solve module Structured flowchart.As shown in figure 4, described ask optimal Fuel- Water Rate solution module to include：KKT conditions modular converter 410, for by described in KKT conditions are converted into the KKT conditions with Min-Algebra type functions；Module 420 is solved, described Min- is carried for solving The KKT conditions of Algebra type functions, to obtain optimal sequence；Determining module 430, described in being determined according to the optimal solution Optimal Fuel- Water Rate.
Wherein, the device of the control separator temperature obtains the process of optimal Fuel- Water Rate and above-mentioned control separator temperature Method described in content it is identical, as space is limited, here is omitted.
Another embodiment of the present invention also provides machinable medium, and finger is stored with the machinable medium Order, the instruction are used for the method for make it that machine performs control separator temperature described above.
The optional embodiment of example of the present invention, still, the embodiment of the present invention and unlimited are described in detail above in association with accompanying drawing Detail in above-mentioned embodiment, can be to the embodiment of the present invention in the range of the technology design of the embodiment of the present invention Technical scheme carry out a variety of simple variants, these simple variants belong to the protection domain of the embodiment of the present invention.
It is further to note that each particular technique feature described in above-mentioned embodiment, in not lance In the case of shield, it can be combined by any suitable means.In order to avoid unnecessary repetition, the embodiment of the present invention pair Various combinations of possible ways no longer separately illustrate.
Prediction and Control Technology is incorporated into the mode of traditional feedforward plus feedback by the present invention, and generating set AGC is controlled The temperature control of separator is optimized in system, and AGC control systems are in overall control structure still using feedforward plus feedback Control model, but with due to applying Predictive Control System, solving the problems, such as control hysteresis.Using the solution of the present invention energy It is enough to be controlled according to the change in future amount of Fuel- Water Rate, process is effectively adjusted in advance, so as to which unit AGC controls greatly improved The closed loop stability and Ability of Resisting Disturbance of separator temperature in system.
It will be appreciated by those skilled in the art that realize that all or part of step in above-described embodiment method is to pass through Program instructs the hardware of correlation to complete, and the program storage is in the storage medium, including some instructions are causing one Individual (can be single-chip microcomputer, chip etc.) or processor (processor) perform the whole of each embodiment methods described of the application Or part steps.And foregoing storage medium includes：USB flash disk, mobile hard disk, read-only storage (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), magnetic disc or CD etc. are various can store journey The medium of sequence code.
In addition, it can also be combined between a variety of embodiments of the embodiment of the present invention, as long as it is not The thought of the embodiment of the present invention is run counter to, it should equally be considered as disclosure of that of the embodiment of the present invention.

Claims (11)

  1. A kind of 1. method for controlling separator temperature, it is characterised in that this method includes：
    Merge separator temperature variable and Fuel- Water Rate variable, to generate optimized variable；
    The optimized variable is brought into the optimal objective function of the separator temperature variable and the Fuel- Water Rate variable, with To quadratic programming problem；
    Derive the KKT conditions of the quadratic programming problem；
    The KKT conditions are solved, to obtain the optimal Fuel- Water Rate of the Fuel- Water Rate；And
    The separator temperature is controlled according to the optimal Fuel- Water Rate.
  2. 2. the method for control separator temperature according to claim 1, it is characterised in that it is described to solve the KKT conditions, Included with obtaining the optimal Fuel- Water Rate of the Fuel- Water Rate：
    The KKT conditions are converted into the KKT conditions with Min-Algebra type functions；
    The KKT conditions for carrying Min-Algebra type functions are solved, to obtain optimal sequence；
    First value is taken from the optimal sequence as optimal solution；And
    The optimal Fuel- Water Rate is determined according to the optimal solution.
  3. 3. the method for control separator temperature according to claim 1 or 2, it is characterised in that the quadratic programming problem For：
    <mrow> <mi>min</mi> <mo>{</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msup> <mi>Z</mi> <mo>&amp;prime;</mo> </msup> <mi>M</mi> <mi>Z</mi> <mo>+</mo> <msup> <mi>c</mi> <mo>&amp;prime;</mo> </msup> <mi>Z</mi> <mo>|</mo> <mi>G</mi> <mi>Z</mi> <mo>&amp;le;</mo> <mi>b</mi> <mo>,</mo> <msub> <mi>A</mi> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mi>Z</mi> <mo>=</mo> <msub> <mi>b</mi> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mo>}</mo> </mrow>
    Wherein,
    X=[x (k+1) ' ... x (k+N) '] ', U=[u (k) ' ... x (k+N-1) '] ', Z=[U'X'] '
    <mrow> <mi>M</mi> <mo>=</mo> <mn>2</mn> <mi>d</mi> <mi>i</mi> <mi>a</mi> <mi>g</mi> <mrow> <mo>(</mo> <mover> <mi>Q</mi> <mo>&amp;OverBar;</mo> </mover> <mover> <mi>R</mi> <mo>&amp;OverBar;</mo> </mover> <mo>)</mo> </mrow> <mo>,</mo> <mover> <mi>Q</mi> <mo>&amp;OverBar;</mo> </mover> <mo>=</mo> <mi>d</mi> <mi>i</mi> <mi>a</mi> <mi>g</mi> <mrow> <mo>(</mo> <mi>Q</mi> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mi>Q</mi> <mo>,</mo> <mi>Q</mi> <mo>)</mo> </mrow> <mo>,</mo> <mover> <mi>R</mi> <mo>&amp;OverBar;</mo> </mover> <mo>=</mo> <mi>d</mi> <mi>i</mi> <mi>a</mi> <mi>g</mi> <mrow> <mo>(</mo> <mi>R</mi> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mi>R</mi> <mo>)</mo> </mrow> </mrow>
    C=[0 ... 0, (- 2Qr) ' ... (- 2Qr) '] '
    B=[umax′,umax′…umax′]′
    <mrow> <msub> <mi>b</mi> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mi>A</mi> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mi>x</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> </mrow>
    Wherein, A, B are the coefficient of the Space admittance of the separator temperature variable and the Fuel- Water Rate variable, and X is represented The separator temperature matrix of variables in time domain N is predicted, U represents the Fuel- Water Rate matrix of variables in prediction time domain N, Z For optimized variable matrix, b is bound variable, and Q is X diagonal positive definite matrix, and R is U diagonal positive definite matrix, and G is constraints Coefficient matrix, I is unit matrix, and r is the setting value of separator temperature, and x (k) is the separator temperature variable, and u (k) is The Fuel- Water Rate variable, umaxFor the higher limit of controlled quentity controlled variable Fuel- Water Rate.
  4. 4. the method for control separator temperature according to claim 3, it is characterised in that the KKT conditions are：
    MZ+c+G′y+Aeq' γ=0
    <mrow> <msub> <mi>G</mi> <mi>i</mi> </msub> <mo>.</mo> <mi>U</mi> <mo>&lt;</mo> <msub> <mi>b</mi> <mi>i</mi> </msub> <mo>&amp;DoubleRightArrow;</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>=</mo> <mn>0</mn> </mrow>
    <mrow> <msub> <mi>b</mi> <mi>i</mi> </msub> <mo>=</mo> <msub> <mi>G</mi> <mi>i</mi> </msub> <mo>.</mo> <mi>U</mi> <mo>&amp;DoubleRightArrow;</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>&amp;GreaterEqual;</mo> <mn>0</mn> </mrow>
    Wherein, γ, yiFor Lagrange coefficient, i is vector order.
  5. 5. the method for control separator temperature according to claim 3, it is characterised in that the Min-Algebra types letter Number is：
    Φmin(y)=0
    Wherein,
    <mrow> <mi>D</mi> <mover> <mo>=</mo> <mi>&amp;Delta;</mi> </mover> <mi>G</mi> <mrow> <mo>(</mo> <mo>-</mo> <msup> <mi>M</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>+</mo> <msup> <mi>M</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msup> <msub> <mi>A</mi> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mo>&amp;prime;</mo> </msup> <msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>A</mi> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> <msup> <mi>M</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msup> <msub> <mi>A</mi> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mo>&amp;prime;</mo> </msup> </mrow> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msub> <mi>A</mi> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> <msup> <mi>M</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>)</mo> </mrow> <msup> <mi>G</mi> <mo>&amp;prime;</mo> </msup> <mo>,</mo> </mrow>
    <mrow> <mi>d</mi> <mover> <mo>=</mo> <mi>&amp;Delta;</mi> </mover> <msup> <mi>GM</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mrow> <mo>(</mo> <mo>-</mo> <mi>c</mi> <mo>+</mo> <msup> <msub> <mi>A</mi> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mo>&amp;prime;</mo> </msup> <msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>A</mi> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> <msup> <mi>M</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msup> <msub> <mi>A</mi> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mo>&amp;prime;</mo> </msup> </mrow> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msub> <mi>A</mi> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> <msup> <mi>M</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mi>c</mi> <mo>+</mo> <msup> <msub> <mi>A</mi> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mo>&amp;prime;</mo> </msup> <msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>A</mi> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> <msup> <mi>M</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msup> <msub> <mi>A</mi> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mo>&amp;prime;</mo> </msup> </mrow> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msub> <mi>b</mi> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>.</mo> </mrow>
  6. 6. a kind of device for controlling separator temperature, it is characterised in that the device includes：
    Quadratic programming problem generation module, for merging separator temperature variable and Fuel- Water Rate variable, to generate optimized variable, and The optimized variable is brought into the optimal objective function of the separator temperature variable and the Fuel- Water Rate variable, to obtain two Secondary planning problem；
    Optimal Fuel- Water Rate acquisition module, for deriving the KKT conditions of the quadratic programming problem, and the KKT conditions are solved, with Obtain the optimal Fuel- Water Rate of the Fuel- Water Rate；And
    Control module, for controlling separator temperature according to the optimal Fuel- Water Rate.
  7. 7. the device of control separator temperature according to claim 6, it is characterised in that described to ask optimal Fuel- Water Rate to solve Module includes：
    KKT condition modular converters, for the KKT conditions to be converted into the KKT conditions with Min-Algebra type functions；
    Module is solved, for solving the KKT conditions for carrying Min-Algebra type functions, to obtain optimal sequence；
    Determining module, for determining the optimal Fuel- Water Rate according to the optimal solution.
  8. 8. the device of the control separator temperature according to claim 6 or 7, it is characterised in that the quadratic programming problem For：
    <mrow> <mi>min</mi> <mo>{</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msup> <mi>Z</mi> <mo>&amp;prime;</mo> </msup> <mi>M</mi> <mi>Z</mi> <mo>+</mo> <msup> <mi>c</mi> <mo>&amp;prime;</mo> </msup> <mi>Z</mi> <mo>|</mo> <mi>G</mi> <mi>Z</mi> <mo>&amp;le;</mo> <mi>b</mi> <mo>,</mo> <msub> <mi>A</mi> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mi>Z</mi> <mo>=</mo> <msub> <mi>b</mi> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mo>}</mo> </mrow>
    Wherein,
    X=[x (k+1) ' ... x (k+N) '] ', U=[u (k) ' ... x (k+N-1) '] ', Z=[U'X'] '
    <mrow> <mi>M</mi> <mo>=</mo> <mn>2</mn> <mi>d</mi> <mi>i</mi> <mi>a</mi> <mi>g</mi> <mrow> <mo>(</mo> <mover> <mi>Q</mi> <mo>&amp;OverBar;</mo> </mover> <mover> <mi>R</mi> <mo>&amp;OverBar;</mo> </mover> <mo>)</mo> </mrow> <mo>,</mo> <mover> <mi>Q</mi> <mo>&amp;OverBar;</mo> </mover> <mo>=</mo> <mi>d</mi> <mi>i</mi> <mi>a</mi> <mi>g</mi> <mrow> <mo>(</mo> <mi>Q</mi> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mi>Q</mi> <mo>,</mo> <mi>Q</mi> <mo>)</mo> </mrow> <mo>,</mo> <mover> <mi>R</mi> <mo>&amp;OverBar;</mo> </mover> <mo>=</mo> <mi>d</mi> <mi>i</mi> <mi>a</mi> <mi>g</mi> <mrow> <mo>(</mo> <mi>R</mi> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mi>R</mi> <mo>)</mo> </mrow> </mrow>
    C=[0 ... 0, (- 2Qr) ' ... (- 2Qr) '] '
    B=[umax′,umax′…umax′]′
    <mrow> <msub> <mi>b</mi> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mi>A</mi> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mi>x</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> </mrow>
    Wherein, A, B are the coefficient of the Space admittance of the separator temperature variable and the Fuel- Water Rate variable, and X is represented The separator temperature matrix of variables in time domain N is predicted, U represents the Fuel- Water Rate matrix of variables in prediction time domain N, Z For optimized variable matrix, b is bound variable, and Q is X diagonal positive definite matrix, and R is U diagonal positive definite matrix, and G is constraints Coefficient matrix, I is unit matrix, and r is the setting value of separator temperature, and x (k) is the separator temperature variable, and u (k) is The Fuel- Water Rate variable, umaxFor the higher limit of controlled quentity controlled variable Fuel- Water Rate.
  9. 9. the device of control separator temperature according to claim 8, it is characterised in that the KKT conditions are：
    MZ+c+G′y+Aeq' γ=0
    <mrow> <msub> <mi>G</mi> <mi>i</mi> </msub> <mo>.</mo> <mi>U</mi> <mo>&lt;</mo> <msub> <mi>b</mi> <mi>i</mi> </msub> <mo>&amp;DoubleRightArrow;</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>=</mo> <mn>0</mn> </mrow>
    <mrow> <msub> <mi>b</mi> <mi>i</mi> </msub> <mo>=</mo> <msub> <mi>G</mi> <mi>i</mi> </msub> <mo>.</mo> <mi>U</mi> <mo>&amp;DoubleRightArrow;</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>&amp;GreaterEqual;</mo> <mn>0</mn> </mrow>
    Wherein, γ, yiFor Lagrange coefficient, i is vector order.
  10. 10. the device of control separator temperature according to claim 8, it is characterised in that the Min-Algebra types letter Number is：
    Φmin(y)=0
    Wherein,
    <mrow> <mi>D</mi> <mover> <mo>=</mo> <mi>&amp;Delta;</mi> </mover> <mi>G</mi> <mrow> <mo>(</mo> <mo>-</mo> <msup> <mi>M</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>+</mo> <msup> <mi>M</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msup> <msub> <mi>A</mi> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mo>&amp;prime;</mo> </msup> <msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>A</mi> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> <msup> <mi>M</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msup> <msub> <mi>A</mi> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mo>&amp;prime;</mo> </msup> </mrow> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msub> <mi>A</mi> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> <msup> <mi>M</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>)</mo> </mrow> <msup> <mi>G</mi> <mo>&amp;prime;</mo> </msup> <mo>,</mo> </mrow>
    <mrow> <mi>d</mi> <mover> <mo>=</mo> <mi>&amp;Delta;</mi> </mover> <msup> <mi>GM</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mrow> <mo>(</mo> <mo>-</mo> <mi>c</mi> <mo>+</mo> <msup> <msub> <mi>A</mi> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mo>&amp;prime;</mo> </msup> <msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>A</mi> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> <msup> <mi>M</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msup> <msub> <mi>A</mi> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mo>&amp;prime;</mo> </msup> </mrow> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msub> <mi>A</mi> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> <msup> <mi>M</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mi>c</mi> <mo>+</mo> <msup> <msub> <mi>A</mi> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mo>&amp;prime;</mo> </msup> <msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>A</mi> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> <msup> <mi>M</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msup> <msub> <mi>A</mi> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mo>&amp;prime;</mo> </msup> </mrow> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <msub> <mi>b</mi> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>.</mo> </mrow>
  11. 11. a kind of machinable medium, instruction is stored with the machinable medium, the instruction is used to cause machine Perform claim requires the method that separator temperature is controlled any one of 1-5.
CN201710446637.7A 2017-06-14 2017-06-14 Control method and apparatus, the storage medium of separator temperature Pending CN107355767A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710446637.7A CN107355767A (en) 2017-06-14 2017-06-14 Control method and apparatus, the storage medium of separator temperature

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710446637.7A CN107355767A (en) 2017-06-14 2017-06-14 Control method and apparatus, the storage medium of separator temperature

Publications (1)

Publication Number Publication Date
CN107355767A true CN107355767A (en) 2017-11-17

Family

ID=60272936

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710446637.7A Pending CN107355767A (en) 2017-06-14 2017-06-14 Control method and apparatus, the storage medium of separator temperature

Country Status (1)

Country Link
CN (1) CN107355767A (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102109172A (en) * 2009-12-24 2011-06-29 江苏方天电力技术有限公司 Fuel-water ratio control method for supercritical and ultra supercritical unit
CN105955023A (en) * 2016-05-16 2016-09-21 浙江大学 Model prediction control method used for rapid motion control system

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102109172A (en) * 2009-12-24 2011-06-29 江苏方天电力技术有限公司 Fuel-water ratio control method for supercritical and ultra supercritical unit
CN105955023A (en) * 2016-05-16 2016-09-21 浙江大学 Model prediction control method used for rapid motion control system

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
郑卫东 等: "先进AGC及汽温控制系统在1000MW超超临界机组的应用", 《浙江电力》 *
郭仁生: "《机械工程设计分析和MATLAB应用》", 30 November 2014, 机械工业出版社 *
麦好: "《机器学习实践指南：案例应用解析》", 31 July 2016, 机械工业出版社 *

Similar Documents

Publication Publication Date Title
Dong et al. Model-free adaptive control law for nuclear superheated-steam supply systems
CN103544551B (en) A kind of operation optimization method and apparatus of steam power system
Wu et al. Reinforced coordinated control of coal-fired power plant retrofitted with solvent based CO2 capture using model predictive controls
Wu et al. Flexible operation of supercritical coal-fired power plant integrated with solvent-based CO2 capture through collaborative predictive control
Li et al. Flexibility control and simulation with multi-model and LQG/LTR design for PWR core load following operation
CN109507910A (en) A kind of fired power generating unit modeling and control method based on data-driven
CN104036334A (en) Cascade hydropower station multi-object optimal scheduling mixed searching method coupling peak shaving and navigation demands
Johansen et al. Gain-scheduled control of a solar power plant
Jiang et al. Practical dynamic matrix control of MHTGR-based nuclear steam supply systems
Dong et al. Coordination control of SMR-based NSSS modules integrated by feedwater distribution
AboRas et al. Voltage and frequency regulation in smart grids via a unique Fuzzy PIDD2 controller optimized by Gradient-Based Optimization algorithm
Tripathy et al. Optimal design of PI/PD dual mode controller based on quasi opposition based learning for power system frequency control
Li et al. Study on the dynamic characteristics, control strategies and load variation rates of the concentrated solar power plant
Wang et al. Simulation study of frequency control characteristics of a generation III+ nuclear power plant
Rosado-Tamariz et al. Optimization of a drum boiler startup using dynamic simulation and a micro-genetic algorithm
CN107355767A (en) Control method and apparatus, the storage medium of separator temperature
Zhao et al. A novel coordinated control for Integrated Pressurized Water Reactor
Kumar et al. Genetic fuzzy PID controller based on adaptive gain scheduling for load frequency control
Dong et al. Power-level control of nuclear reactors based on feedback dissipation and backstepping
Varshney et al. Fast frequency response using model predictive control for a hybrid power system
CN114547983A (en) Improved multi-population genetic algorithm-based reactor operation optimization method
Nguyen Study on a novel speed control strategy of hydraulic turbines applying artificial neural network and metaheuristic optimization technique
Yang et al. Module coordination control of multi-modular lead-based reactor systems
Malik et al. Fractional order modelling and robust multi-model intelligent controllers' synthesis for ACP1000 nuclear power plant
Oluseyi et al. Optimal load frequency control of two area power system

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20171117

RJ01 Rejection of invention patent application after publication