JPH11294741A - Combustion control method in incineration plant with boiler - Google Patents

Combustion control method in incineration plant with boiler

Info

Publication number
JPH11294741A
JPH11294741A JP10670298A JP10670298A JPH11294741A JP H11294741 A JPH11294741 A JP H11294741A JP 10670298 A JP10670298 A JP 10670298A JP 10670298 A JP10670298 A JP 10670298A JP H11294741 A JPH11294741 A JP H11294741A
Authority
JP
Japan
Prior art keywords
pressure
boiler
set value
value
pressure set
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
JP10670298A
Other languages
Japanese (ja)
Inventor
Yoichi Kobayashi
洋一 小林
Kenichi Suematsu
健一 末松
Hisashi Okada
久 岡田
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.)
Ebara Corp
Original Assignee
Ebara Corp
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 Ebara Corp filed Critical Ebara Corp
Priority to JP10670298A priority Critical patent/JPH11294741A/en
Publication of JPH11294741A publication Critical patent/JPH11294741A/en
Pending legal-status Critical Current

Links

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E20/00Combustion technologies with mitigation potential
    • Y02E20/12Heat utilisation in combustion or incineration of waste

Landscapes

  • Incineration Of Waste (AREA)

Abstract

PROBLEM TO BE SOLVED: To achieve a stabilization and an efficiency in the generated output of electric power by stabilizing the boiler outlet pressure and the flow rate of steam. SOLUTION: A boiler pressure process value PPV (t) at a boiler outlet is measured cyclically and a pressure set value is corrected according to the deviation between the pressure process value and a pressure set value PSV (t) to obtain a corrected pressure set value PSV (t)'. The corrected pressure set value is supplied to a boiler pressure adjusting meter 2 in stead of the pressure set value and based on the value, the boiler pressure adjusting meter generates an output for adjusting manipulated variables K (t) such as number of revolutions of a dust feeder to adjust the amount of a fuel to be supplied to an incinerator. The correction of the set value is executed by performing a model forecasting control based on a model such as ARX model previously generated using the pressure set value of a boiler at the boiler outlet as input and the pressure process value as output. A more proper control is possible by correcting the coefficient of the model using the corrected pressure set value and the pressure process value obtained by the measurement.

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION

【0001】[0001]

【発明の技術分野】本発明は、一般廃棄物(都市ごみ)
や産業廃棄物等を焼却するためのボイラ付き焼却プラン
トにおける燃焼制御方法に関し、特に、定常運転時にボ
イラ圧力及びボイラ蒸気流量を安定化させ、それにより
発電量(発電電力)を安定化させるための燃焼制御方法
に関する。
TECHNICAL FIELD The present invention relates to general waste (municipal waste).
Control method in an incineration plant with a boiler for incineration of waste and industrial waste, etc., in particular, for stabilizing boiler pressure and boiler steam flow rate during steady operation, and thereby stabilizing power generation (power generation) The present invention relates to a combustion control method.

【0002】[0002]

【従来の技術】近年、都市ごみや産業廃棄物の焼却熱を
利用して蒸気タービンを駆動することにより、発電を行
うボイラ付き焼却プラントが開発され稼働している。こ
のような焼却プラントにおいては、燃料として用いる都
市ごみや産業廃棄物等の発熱量が変動するので、発電量
が所望の値となるように制御するためには、プラント自
体の応答遅延特性も相まって、経験則等に基づいた複雑
な制御を必要としている。このような制御を行うため
に、焼却プラントにおいては、定常時、分散制御システ
ム(DCS)等により、基本的に1入力1出力のPID
調節器を複数組み合わせたコントローラを構成して、ボ
イラ出口における蒸気圧および蒸気流量等の種々のプロ
セス値が設定範囲(設定値±許容誤差)内になるよう自
動燃焼制御(ACC)している。
2. Description of the Related Art In recent years, an incineration plant with a boiler for generating electricity by driving a steam turbine using incineration heat of municipal solid waste and industrial waste has been developed and operated. In such an incineration plant, the calorific value of municipal solid waste and industrial waste used as fuel fluctuates, and in order to control the amount of power generation to a desired value, the response delay characteristics of the plant itself are combined. And complicated control based on empirical rules and the like. In order to perform such control, in an incineration plant, a PID of one input and one output is basically used in a steady state by a distributed control system (DCS) or the like.
A controller is formed by combining a plurality of regulators, and automatic combustion control (ACC) is performed so that various process values such as a steam pressure and a steam flow rate at a boiler outlet fall within a set range (set value ± tolerance).

【0003】図5は、ボイラ付きごみ焼却プラントにお
ける従来例のPID自動燃焼制御を説明するためのブロ
ック図である。図5において、1はボイラ付きごみ燃焼
プラントであり、該プラント1は、ごみ焼却炉、廃熱ボ
イラ、給塵機等を含み、更に、ボイラ出口圧力センサ、
ボイラ出口蒸気流量センサ、焼却炉出口温度センサ等の
各種センサを含んでいる。2は、ボイラ出口圧力(蒸気
圧)調節計であり、オペレータ等により設定されるボイ
ラ出口圧力設定値Psv(t)と、プラント1中のボイラ
出口圧力センサからのボイラ出口圧力プロセス値Ppv
(t)とが入力され、それらの差をゼロとする方向に給
塵機の回転数等の操作量K(t)を変更してごみ供給量
を調節するよう構成されている。3は、ボイラ蒸気流量
調節計であり、オペレータ等により設定されるボイラ蒸
気流量設定値Qsv(t)と、ボイラ出口蒸気流量センサ
から蒸気流量プロセス値Qpv(t)とが入力され、それ
らの差をゼロとする方向に蒸気流量調節弁の開度R
(t)を調整するよう構成されている。4は炉頂温度調
節系等を含むその他の調節計であり、プラント1中の適
宜のセンサからのプロセス値とそれに対応する設定値と
が供給され、それらの差をゼロとする方向に適宜の操作
量を調整するように構成されている。
FIG. 5 is a block diagram for explaining a conventional PID automatic combustion control in a waste incineration plant with a boiler. In FIG. 5, reference numeral 1 denotes a boiler-equipped refuse combustion plant. The plant 1 includes a refuse incinerator, a waste heat boiler, a dust collector, and the like, and further includes a boiler outlet pressure sensor,
Various sensors such as a boiler outlet steam flow sensor and an incinerator outlet temperature sensor are included. Reference numeral 2 denotes a boiler outlet pressure (steam pressure) controller, and a boiler outlet pressure set value Psv (t) set by an operator or the like and a boiler outlet pressure process value Ppv from a boiler outlet pressure sensor in the plant 1.
(T) is input, and the amount of dust is adjusted by changing the operation amount K (t) such as the number of revolutions of the duster in a direction to make the difference between them zero. Reference numeral 3 denotes a boiler steam flow controller, which receives a boiler steam flow set value Qsv (t) set by an operator or the like and a steam flow process value Qpv (t) from a steam flow sensor at the boiler outlet. Of the steam flow control valve in the direction where
It is configured to adjust (t). Reference numeral 4 denotes another controller including a furnace top temperature control system, which is supplied with a process value from an appropriate sensor in the plant 1 and a set value corresponding to the process value. It is configured to adjust the operation amount.

【0004】これらの調節計2〜4は、所定の周期でP
ID燃焼制御を行うことにより、上記したように、プロ
セス値が目標とする設定範囲(設定値±許容誤差)とな
るよう操作量を自動調節し、さらに、プラント1を平衡
状態に保つように制御している。例えば、ボイラ出口圧
力プロセス値Ppv(t)がその設定値Psv(t)よりも
上昇したとき、ボイラ圧力調節計2の働きによって給塵
機の単位時間当たりの回転数を減少させるよう制御す
る。これに伴い、ボイラ入熱量が低下するのでボイラ出
口圧力プロセス値Ppv(t)は低下しボイラ出口蒸気流
量プロセス値Qpv(t)は減少しようとする。しかしな
がら、このとき、ボイラ蒸気流量調節計3の作用によっ
て、蒸気流量調節弁の開度が増大してボイラ出口圧力プ
ロセス値Ppv(t)の低下を助長すると共にボイラ出口
蒸気流量プロセス値Qpv(t)の減少を相殺して、プラ
ント1は平衡状態に達する。ボイラ出口圧力プロセス値
Ppv(t)が設定値より下降した場合は、上記と反対の
動作を行う。このようにして、各プロセス値が設定範囲
に入るようにするとともに、プラント1を平衡状態に保
つよう調整している。
[0004] These controllers 2 to 4 provide P
By performing the ID combustion control, as described above, the operation amount is automatically adjusted so that the process value falls within a target set range (set value ± permissible error), and further, control is performed so as to keep the plant 1 in an equilibrium state. doing. For example, when the boiler outlet pressure process value Ppv (t) is higher than the set value Psv (t), the boiler pressure controller 2 controls to reduce the number of rotations of the duster per unit time. As a result, the boiler heat input decreases, so that the boiler outlet pressure process value Ppv (t) decreases and the boiler outlet steam flow process value Qpv (t) tends to decrease. However, at this time, the opening of the steam flow control valve is increased by the action of the boiler steam flow controller 3 to promote a decrease in the boiler outlet pressure process value Ppv (t), and the boiler outlet steam flow process value Qpv (t). ), The plant 1 reaches an equilibrium state. When the boiler outlet pressure process value Ppv (t) falls below the set value, the operation opposite to the above is performed. In this way, adjustment is made so that each process value falls within the set range and the plant 1 is kept in an equilibrium state.

【0005】[0005]

【発明が解決しようとする課題】しかしながら、上記し
た従来例の燃焼制御方法においてもなお、ボイラ出口圧
力や蒸気流量が10分〜数10分程度の周期で変動し
て、その目標設定範囲から外れてしまうことがある。こ
のような変動原因の主要なものとして、焼却炉へ供給さ
れるごみ密度の変化等のごみ質変化を上げることができ
る。すなわち、ホッパ下部のごみ圧密度に着目すると、
ごみがクレーンによってホッパへ投入された時点ではご
み圧密度が最も大きく、その後徐々に小さくなり、次の
ごみ投入時点で再びごみ圧密度が大きくなる。このよう
なごみ圧密度の周期的変化に起因して、焼却炉へ供給さ
れるごみ密度が変化し、これにより、ボイラへの単位時
間当たりの入熱量が給塵装置の操作回数に比例しないこ
とになる。したがって、ごみ質を一定としてボイラ出口
圧力のPID制御を行っている従来例の燃焼制御方法に
おいては、ボイラ出口圧力及び蒸気流量が変動してしま
い、安定した発電量を得ることが困難な場合がある。な
お、ボイラ出口圧力の変動は、ごみ密度の変動以外の理
由、例えば、ごみ収集時の天候等に基づく水分含有量の
変化や、収集地域によるごみ質の違いによっても生じて
しまうものであり、更には、その他の外乱によっても生
じてしまうものである。外乱によるボイラ出口圧力の変
動は、ボイラ付きの産業廃棄物焼却プラントにおいても
生じてしまうものである。
However, even in the above-described conventional combustion control method, the boiler outlet pressure and the steam flow rate fluctuate at a period of about 10 minutes to several tens of minutes, and thus deviate from the target set range. Sometimes. As a main cause of such fluctuation, a change in waste quality such as a change in the density of the waste supplied to the incinerator can be raised. In other words, focusing on the dust pressure density at the bottom of the hopper,
The waste pressure density is highest when the refuse is put into the hopper by the crane, then gradually decreases, and then again increases at the next refuse input. Due to such a periodic change in the waste pressure density, the waste density supplied to the incinerator changes, and as a result, the heat input to the boiler per unit time is not proportional to the number of operations of the dust supply device. Become. Therefore, in the conventional combustion control method in which PID control of the boiler outlet pressure is performed while keeping the waste quality constant, the boiler outlet pressure and the steam flow rate may fluctuate, making it difficult to obtain a stable power generation amount. is there. The fluctuation of the boiler outlet pressure is also caused by a reason other than the fluctuation of the waste density, for example, a change in the water content based on the weather or the like at the time of collecting the waste, and a difference in the waste quality depending on the collection area. Furthermore, it is also caused by other disturbances. Fluctuations in boiler outlet pressure due to disturbances also occur in industrial waste incineration plants with boilers.

【0006】そして、ボイラ出口圧力が変動すると、ボ
イラ蒸気流量が変動してしまい、結局、発電量が変動し
てしまう。ボイラ出口圧力の変動による発電量の変動を
抑制するために、図5に示した従来例において、ボイラ
蒸気流量調節計3により蒸気流量調節弁の開度R(t)
を制御することも考えられるが、単に該調節弁の開度を
調節しても、ボイラ出口圧力は積分系であるため、入熱
量の変化が即座には圧力の変化として顕れない。このた
め、ボイラ出口圧力の変動がより増長されてしまい、結
局、発電量を安定化できない場合がある。また、燃料供
給量で蒸気流量を制御し、かつ調節弁でボイラ出口圧力
を調節する方式を採用した構成においても、上記と同様
な理由により、圧力を一定にすればするほど蒸気流量の
変化が大きくなってしまい、この場合も発電量を安定化
できない場合がある。本発明の目的は、ボイラ付き焼却
プラントにおいて、このような外乱によるボイラ出口圧
力の変動を低減させて焼却プラントの安全性及び信頼性
を向上させるとともに、それによって蒸気流量を安定化
して発電量の安定化及び効率化を図ることができる燃焼
制御方法を提供することである。
[0006] When the boiler outlet pressure fluctuates, the boiler steam flow rate fluctuates, and eventually, the power generation amount fluctuates. In order to suppress the fluctuation of the power generation amount due to the fluctuation of the boiler outlet pressure, in the conventional example shown in FIG. 5, the opening degree R (t) of the steam flow control valve is controlled by the boiler steam flow controller 3.
However, even if the opening of the control valve is simply adjusted, the change in heat input does not immediately appear as a change in pressure because the boiler outlet pressure is an integral system. For this reason, the fluctuation of the boiler outlet pressure is further increased, and eventually, the power generation amount may not be stabilized. Also, in a configuration in which the steam flow rate is controlled by the fuel supply amount and the boiler outlet pressure is adjusted by the control valve, for the same reason as above, the more the pressure is kept constant, the more the steam flow rate changes. In this case, the power generation amount may not be stabilized. An object of the present invention is to improve the safety and reliability of an incineration plant by reducing fluctuations in boiler outlet pressure due to such disturbances in an incineration plant with a boiler, thereby stabilizing the steam flow rate and thereby reducing the amount of power generation. An object of the present invention is to provide a combustion control method capable of stabilizing and improving efficiency.

【0007】[0007]

【課題を解決するための手段】上記した目的を達成する
ために、本発明のボイラ付き焼却プラントにおける燃焼
制御方法においては、ボイラ出口におけるボイラの圧力
プロセス値を周期的に測定するステップと、測定により
得られた圧力プロセス値と圧力設定値との偏差に応じて
圧力設定値を補正することにより、補正済圧力設定値を
得る設定値補正ステップと、得られた補正済圧力設定値
を圧力設定値の代わりにボイラ圧力調節計に供給するス
テップと、ボイラ圧力調節計の出力に基づいて、焼却炉
への燃料の供給量を調節する燃料調節ステップとを含む
ことを特徴としている。
In order to achieve the above object, a method for controlling combustion in a boiler-mounted incineration plant according to the present invention comprises the steps of periodically measuring a boiler pressure process value at a boiler outlet; A set value correction step of obtaining a corrected pressure set value by correcting the pressure set value according to the deviation between the pressure process value and the pressure set value obtained by the pressure setting, and setting the obtained corrected pressure set value to the pressure setting The method is characterized by including a step of supplying to the boiler pressure controller instead of the value, and a fuel adjusting step of adjusting a supply amount of fuel to the incinerator based on an output of the boiler pressure controller.

【0008】本発明の好適な実施例においては、設定値
補正ステップは、ボイラ出口におけるボイラの圧力設定
値を入力としかつ圧力プロセス値を出力として予め作成
された数式モデルに基づいて、モデル予測制御を行うこ
とにより、圧力設定値を補正して補正済圧力設定値を得
るよう構成されている。更に、本発明の好適な実施例に
おいては、設定値補正ステップはさらに、補正済圧力設
定値と測定により得られた圧力プロセス値とを用いて、
数式モデルの係数を補正するステップを含んでいる。ま
た、燃料調節ステップは、ボイラ圧力調節計の出力に基
づいて、焼却炉へ廃棄物を供給する給塵供給手段の駆動
回数を調節するよう構成されている。さらに、予め作成
された数式モデルはARXモデルである。
In a preferred embodiment of the present invention, the set value correcting step includes a model predictive control based on a mathematical model created in advance using a boiler pressure set value at the boiler outlet as an input and a pressure process value as an output. , The pressure set value is corrected to obtain a corrected pressure set value. Further, in a preferred embodiment of the present invention, the set value correcting step further comprises using the corrected pressure set value and the pressure process value obtained by the measurement,
Correcting the coefficients of the mathematical model. Further, the fuel adjusting step is configured to adjust the number of times of driving of the dust supply means for supplying waste to the incinerator based on the output of the boiler pressure controller. Further, the mathematical model created in advance is an ARX model.

【0009】[0009]

【発明の実施の態様】図1は、本発明の第1及び第2の
実施例の燃焼制御方法を実行するためのシステム機能ブ
ロック図を示しており、図1において、図5の構成要素
と同一の構成要素には、同一の参照番号が付されてい
る。本発明の燃焼制御方法においては、ボイラ圧力調節
計2に設定値補正量演算手段5及び加算手段6を付加
し、該設定値補正量演算手段5において、ボイラ出口に
設けられたボイラ出口圧力センサからのボイラ出口圧力
プロセス値(以下、単に「圧力プロセス値」と称する)
Ppv(t)及びボイラ出口圧力設定値(以下、単に「圧
力設定値」と称する)Psv(t)に基づいて、圧力設定
値Psv(t)に対する圧力補正値Paj(t)を演算し、
加算手段6において、これらの圧力設定値と圧力補正値
とを加算し、得られた値を補正済圧力設定値Psv′
(t)として、ボイラ圧力調節計2に供給することを特
徴としている。
FIG. 1 shows a functional block diagram of a system for executing a combustion control method according to first and second embodiments of the present invention. In FIG. 1, the components shown in FIG. The same components are denoted by the same reference numerals. In the combustion control method of the present invention, a set value correction amount calculating means 5 and an adding means 6 are added to the boiler pressure controller 2, and the set value correction amount calculating means 5 includes a boiler outlet pressure sensor provided at a boiler outlet. Boiler outlet pressure process value (hereinafter simply referred to as "pressure process value")
Based on Ppv (t) and a boiler outlet pressure set value (hereinafter simply referred to as “pressure set value”) Psv (t), a pressure correction value Paj (t) for the pressure set value Psv (t) is calculated,
The adding means 6 adds the pressure set value and the pressure correction value, and adds the obtained value to the corrected pressure set value Psv '.
As (t), it is supplied to the boiler pressure controller 2.

【0010】本発明の第1の実施例においては、設定値
補正量演算手段5において、圧力補正値Paj(t)及び
補正済圧力設定値Psv′(t)は、以下の式(1)又は
(2)に基づいて演算される。なお、「t」は、補正制
御が開始されてから所定の周期で実行されるPID制御
の実行時点(すなわち実行回)を表すものとする。
In the first embodiment of the present invention, in the set value correction amount calculating means 5, the pressure correction value Paj (t) and the corrected pressure set value Psv '(t) are calculated by the following equation (1) or It is calculated based on (2). Note that “t” represents the execution point (ie, the number of executions) of the PID control that is executed at a predetermined cycle after the correction control is started.

【数1】 Paj(t)=−{Ppv(t)−Psv(t)−Paj(t−1)} …(1) =−Σ(Ppv(t)−Psv(t)) …(2) ただし、Σは、t=1から現在実行時点についての加算 また、加算手段6においては、補正済圧力設定値Psv′
(t)
Paj (t) = − {Ppv (t) −Psv (t) −Paj (t−1)} (1) = −… (Ppv (t) −Psv (t)) (2) However, Σ is an addition from t = 1 to the current execution time.
(T)

【数2】 Psv′(t)=Psv(t)+Paj(t) …(3) が演算される。Psv ′ (t) = Psv (t) + Paj (t) (3)

【0011】式(1)から明らかなように、本発明の制
御方法においては、圧力プロセス値Ppv(t)の圧力設
定値Psv(t)からの偏差が今回の圧力補正値Paj
(t)と前回の圧力補正値Paj(t−1)との差となる
ように、今回の圧力補正値Paj(t)を設定しており、
Ppv(t)のPsv(t)との偏差に応じて圧力補正量を
増減している。また、式(2)及び(3)から、本発明
の制御方法においては、補正制御が開始されてから現時
点までの各実行時点における圧力プロセス値Ppv(t)
の圧力設定値Psv(t)からの偏差を求めてそれを積分
し、得られた積分値だけ該積分値と反対方向に、その時
点の圧力設定値Psv(t)を補正することにより、補正
済圧力設定値Psv′(t)を得ていると言うことができ
る。したがって、圧力プロセス値と圧力設定値との偏差
の履歴に基づいて圧力設定値を補正しているので、プロ
セス値の変動振幅を縮小することができる。
As is apparent from the equation (1), in the control method of the present invention, the deviation of the pressure process value Ppv (t) from the pressure set value Psv (t) is the current pressure correction value Paj.
The current pressure correction value Paj (t) is set so as to be the difference between (t) and the previous pressure correction value Paj (t-1).
The pressure correction amount is increased or decreased according to the deviation of Ppv (t) from Psv (t). Also, from the equations (2) and (3), in the control method of the present invention, the pressure process value Ppv (t) at each execution time from the start of the correction control to the present time is shown.
Of the pressure set value Psv (t) from the pressure set value Psv (t), and integrating the obtained value, and correcting the pressure set value Psv (t) at that point in the direction opposite to the integrated value by the obtained integrated value. It can be said that the completed pressure set value Psv '(t) has been obtained. Accordingly, since the pressure set value is corrected based on the history of the deviation between the pressure process value and the pressure set value, the fluctuation amplitude of the process value can be reduced.

【0012】上記した本発明の第1の実施例の燃焼制御
方法においては、過去の履歴に基づいて現在の圧力設定
値を補正することを意図しているが、本発明の第2の実
施例の燃焼制御方法においては、設定値補正量演算手段
5においてモデル予測制御を行うことにより、現在まで
の実測データに基づいて未来を予測し、現在の設定値の
補正量を最適となるようにするものである。したがっ
て、本発明の第2の実施例によれば、プラント1及び調
節計2〜4からなるシステムの変動周期と、圧力設定値
の変化に対するプロセス値の応答時間との比(変動周期
/応答時間)が比較的小さい場合にも、効果を発揮でき
る。以下に、本発明の第2の実施例におけるモデル予測
制御について説明する。
The above-described combustion control method according to the first embodiment of the present invention intends to correct the current pressure set value based on the past history, but the second embodiment of the present invention. In the combustion control method described above, the set value correction amount calculating means 5 performs model prediction control to predict the future based on actual measurement data up to the present time, and to optimize the correction amount of the current set value. Things. Therefore, according to the second embodiment of the present invention, the ratio of the fluctuation period of the system including the plant 1 and the controllers 2 to 4 to the response time of the process value to the change of the pressure set value (variation period / response time) The effect can be exerted even when () is relatively small. Hereinafter, the model prediction control according to the second embodiment of the present invention will be described.

【0013】モデル予測制御においては、まず事前に、
入力u(t)=圧力設定値Psv(t)とし、出力y
(t)=圧力プロセス値Ppv(t)とする、プラントと
従来制御系とからなるシステムを同定し、以下の演算式
(4)で表される予測制御用のARXモデルを決定す
る。
In the model predictive control, first,
Input u (t) = pressure set value Psv (t), output y
A system consisting of a plant and a conventional control system, where (t) = pressure process value Ppv (t) is identified, and an ARX model for predictive control represented by the following equation (4) is determined.

【数3】 A(q)y(t)=B(q)u(t)+e(t) …(4) ただし、A(q)=1+Σan-n …(5) [Σ:n=1〜naについての加算] B(q)=q-nkΣbn-n+1 …(6) [Σ:n=1〜nbについての加算] e(t):外乱 t:時間(t=……,−2,−1,0,1,2,……) na≧1,nb≧1,nk≧1 q-1:シフトオペレータ q-nx(t)=x(t−n) ここで、システムの応答時間がゼロでないことや、実時
間の位相遅れフィルタリングを前提にすれば、通常nk
≧1としてよいが、nk≧0として定式化することも可
能である。
Equation 3] A (q) y (t) = B (q) u (t) + e (t) ... (4) However, A (q) = 1 + Σa n q -n ... (5) [Σ: n = 1 to n added for a] B (q) = q -nk Σb n q -n + 1 ... (6) [Σ: n = 1~n addition of b] e (t): disturbance t: time ( t =..., -2, -1, 0, 1, 2,...) n a ≧ 1, n b ≧ 1, nk ≧ 1 q −1 : shift operator q −n x (t) = x ( t−n) Here, assuming that the response time of the system is not zero and that real-time phase lag filtering is assumed, nk is usually used.
Although it may be ≧ 1, it is also possible to formulate as nk ≧ 0.

【0014】次に、予測制御を行うが、現時点をt0
し、1回の演算に要する時間が無視できるほど小さいも
のとし、入力はt0を基点に操作周期TMで期間pTM
渡って変更可能であり、TS≧nk、TE≧nk、TP=p
M−TS+TEとしてt≧t0ではこの制御以外の手段に
よって入力が変更されないとの仮定をおく。これら時間
に関する関係は、図2に示されている通りである。そし
て、以下の式(7)で表されるARXモデルによる出力
の予測値列YPが、式(8)で表される予め定められた
参照軌道YRにできるだけ一致し、かつ、入力の変化が
できるだけ小さいような入力値列UFM[式(9)]を求
める。なお、参照軌道YRは、例えば、適当な時定数を
もつ一次遅れ曲線として定義しておく。
[0014] Next, performs predictive control, the current time and t 0, is assumed as the time required for one operation is negligible small, the input over a period pT M in operation period T M and t 0 to the base point T Snk , T Enk , T P = p
It is assumed that the input is not changed by means other than this control when t ≧ t 0 as T M −T S + T E. The relationship with respect to time is as shown in FIG. Then, the predicted value sequence Y P of the output by the ARX model represented by the following equation (7) matches the predetermined reference trajectory Y R represented by the equation (8) as much as possible, and the input change Is obtained as an input value sequence U FM [Equation (9)] that is as small as possible. The reference trajectory Y R is defined as a first-order lag curve having an appropriate time constant, for example.

【数4】 YP=[yP(t0+TS) yP(t0+TS+1) …… yP(t0+TS+TP)]T …(7) YR=[yR(t0+TS) yR(t0+TS+1) …… yR(t0+TS+TP)]T …(8) UFM=[uF(t0) uF(t0+TM) …… uF(t0+pTM)]T …(9) そして、入力値列UFM中の最初の値uF(t0)を実際に
プロセスに入力する。以後、操作周期TM毎に、この操
作を繰り返す。
Y P = [y P (t 0 + T S ) y P (t 0 + T S +1)... Y P (t 0 + T S + T P )] T (7) Y R = [y R ( t 0 + T S) y R (t 0 + T S +1) ...... y R (t 0 + T S + T P)] T ... (8) U FM = [u F (t 0) u F (t 0 + T M) ...... u F (t 0 + pT M)] T ... (9) and, actually input to process the first value u F in the input value column U FM (t 0). Thereafter, this operation is repeated for each operation cycle T M.

【0015】以下に、入力値UFMの求め方を具体的に示
す。まず、nk≧1の場合、ARXモデルによるt0+s
の出力予測値は以下の式(10)のように表される。
[0015] Hereinafter, concrete shows how to determine the input value U FM. First, if n k ≧ 1, t 0 + s based on the ARX model
Is represented as in the following equation (10).

【数5】 yP(t0+s) =−ΣasnO(t0−n+1)+ΣbsnΔuF(t0+s−n+1) +ΣGsnΔuO(t0−n)+dsu(t0−nb') …(10) ただし、△uF(t)=uF(t)−uF(t−1) △uO(t)=uO(t)−uO(t−1) Σ:asnに関しては、n=1〜naについての加算 bsnに関しては、n=1〜(s+1)についての加算 Gsnに関しては、n=1〜(nb'−1)についての加
算 ここで、ベクトルY0、△UF、△UO及びUHを、以下の
式(11)〜(14)のように定義すると、式(15)
〜(18)の形式で表される行列A〜Dに対して、式
(19)が成立する。ただし、UHは、要素数(TP
1)の列ベクトルとする。
Equation 5] y P (t 0 + s) = -Σa s, n y O (t 0 -n + 1) + Σb s, n Δu F (t 0 + s-n + 1) + ΣG s, n Δu O (t 0 -n) + d s u (t 0 -n b ') ... (10) However, △ u F (t) = u F (t) -u F (t-1) △ u O (t) = u O (t) - u O (t-1) Σ : a s, with respect to the n, n = 1 to n a sum b s for, for n, n = 1~ (s + 1) adding G s for, for n is, n = 1 to (n b '-1) Here, if the vectors Y 0 , △ U F , △ U O and U H are defined as in the following expressions (11) to (14), the expression ( 15)
Expression (19) holds for the matrices A to D expressed in the form of (18). Where U H is the number of elements (T P +
Let it be the column vector of 1).

【数6】 YO=[yO(t0) yO(t0−1) … yO(t0−na+1)]T …(11) △UF=[△uF(t0) △uF(t0+1) … △uF(t0+pTM)]T …(12) △UO =[△uO(t0−1) △uO(t0−2) … △uO(t0−nb'+1)]T …(13) UH =[uO(t0−nb') uO(t0−nb') … uO(t0−nb')]T …(14)[6] Y O = [y O (t 0) y O (t 0 -1) ... y O (t 0 -n a +1)] T ... (11) △ U F = [△ u F (t 0 ) △ u F (t 0 +1 ) ... △ u F (t 0 + pT M)] T ... (12) △ U O = [△ u O (t 0 -1) △ u O (t 0 -2) ... △ u O (t 0 -n b ' +1)] T ... (13) U H = [u O (t 0 -n b') u O (t 0 -n b ') ... u O (t 0 -n b ')] T ... (14)

【数7】 (Equation 7)

【0016】さらに、要素数p+1の列ベクトル△
FM、要素数TMの列ベクトルgM及び(pTM+1)x
(p+1)の行列GFを、以下の式(20)、(21)
及び(22)のようにおくと、式(23)が成り立つ。
Further, a column vector p of the number of elements p + 1
U FM, column vector g number of elements T M M and (pT M +1) x
The matrix G F of (p + 1) is calculated by the following equations (20) and (21).
And (22), equation (23) holds.

【数8】 したがって、B′=BGFと置くと、式(23)から式
(19)は、以下のように表すことができる。
(Equation 8) Thus, placing and B '= BG F, formula (19) from equation (23) can be expressed as follows.

【数9】 YP=−AYO+B′△UFM+C△UO+DUH …(24)Y P = −AY O + B ′ △ U FM + C △ U O + DU H (24)

【0017】今、モデル予測制御の評価関数Jを以下の
式(25)のように設定すると、該評価関数Jの偏微分
∂J/∂△UFMは、式(26)で表される。
Now, when the evaluation function J of the model predictive control is set as in the following equation (25), the partial differential ∂ △ J / ∂ △ U FM of the evaluation function J is expressed by the equation (26).

【数10】 J=‖YR−YP2+ξ‖△UFM2 =‖YR+AYO−B′△UFM−C△UO−DUH2+ξ‖△UFM2 …(25) ただし、‖X‖2=XTX ξ>0 ∂J/∂△UFM =[△∂J/∂△U1 ∂J/∂△U2 … ∂J/∂△Up+1T =−2(B′)T(YR+AYO−B′△UFM−C△UO−DUH)+2△ξUFM …(26) ただし、△Uiは、△UFMの第i成分Equation 10] J = ‖Y R -Y P2 + ξ‖ △ U FM 2 = ‖Y R + AY O -B '△ U FM -C △ U O -DU H ‖ 2 + ξ‖ △ U FM2 ... (25) However, ‖X‖ 2 = X T X ξ > 0 ∂J / ∂ △ U FM = [△ ∂J / ∂ △ U 1 ∂J / ∂ △ U 2 ... ∂J / ∂ △ U p + 1] T = -2 However (B ') T (Y R + AY O -B' △ U FM -C △ U O -DU H) + 2 △ ξU FM ... (26), △ U i is the △ U FM I-th component

【0018】ここで、∂J/∂△UFM=0とすると、Here, assuming that ∂J / ∂ △ U FM = 0,

【数11】 △UFM ={(B′)TB′+ξE}-1(B′)T(YR+AYO−C△UO−DUH) …(27) ただし、Eは、単位行列 が得られ、出力の予測値列YPが参照軌道YRに近く、か
つ入力値列の変化‖△UFM‖が小さくなるような△UFM
が求められ、したがって、入力値列UFMが求められる。
なお、得られた入力値列UFMは、式(4)〜(6)に関
連して説明したように、ボイラ圧力調節系2に直接供給
すべきボイラ圧力設定値であり、よって、得られた入力
値列UFMの個々の要素の値が、対応する時点での補正済
圧力設定値Psv′(t)となる。したがって、図1に示
された設定値演算手段5は、得られた入力値列の第1要
素の値から対応する時点の設定値Psv(t)を減算した
値を圧力補正値Paj(t)として、その時点で加算手段
6に供給する。また、加算手段6を介さずに、得られた
入力値列UFMの第1要素の値そのものをボイラ圧力制御
系2に直接供給してもよい。
Equation 11] △ U FM = {(B ' ) T B' + ξE} -1 (B ') T (Y R + AY O -C △ U O -DU H) ... (27) However, E is a unit matrix is obtained, the output of the predicted value sequence Y P reference trajectory Y R close to, and the change in the input value column ‖ △ U FM ‖ such decrease △ U FM
Is obtained, and thus the input value sequence UFM is obtained.
Note that the obtained input value sequence U FM is a boiler pressure set value to be directly supplied to the boiler pressure control system 2 as described in relation to the equations (4) to (6). the value of each element of the input value column U FM has becomes the corrected pressure setpoint Psv at the time the corresponding '(t). Therefore, the set value calculating means 5 shown in FIG. 1 calculates the pressure correction value Paj (t) by subtracting the set value Psv (t) at the corresponding time from the value of the first element of the obtained input value sequence. Is supplied to the adding means 6 at that time. Further, without using the addition means 6, the value itself of the first element of the obtained input value sequence U FM may be fed directly to the boiler pressure control system 2.

【0019】次に、本発明の第3の実施例における燃焼
制御方法を説明する。第3の実施例は、式(4)〜
(6)で表されるARXモデルに基づいて圧力補正量P
aj(t)を求める場合、ARXモデルの係数an及びbn
の調整を行うことができるようにしたものである。この
ため、図3に示すように、予測演算モデル係数補正手段
7を追加して、該係数補正手段7において、設定値補正
量演算手段5及び加算手段6を作動させた運転状態で得
られる補正済圧力設定値Psv′(t)及び圧力プロセス
値Ppv(t)に基づいて、ARXモデルの係数an及び
nを学習により補正するように構成している。該係数
の学習は、時点t0において計算されたモデル予測値を
p(t)と表し、t0≦t≦t0+Tの範囲で、圧力プ
ロセス値y(t)を教師信号として、モデル予測値yp
(t)が改善されるように、係数を補正する。いま、評
価関数を
Next, a description will be given of a combustion control method according to a third embodiment of the present invention. The third embodiment is based on Equations (4) to (4).
The pressure correction amount P based on the ARX model represented by (6)
When seeking aj (t), the coefficient of ARX model a n and b n
Can be adjusted. For this reason, as shown in FIG. 3, a prediction calculation model coefficient correction means 7 is added, and the correction value obtained in the operating state in which the set value correction amount calculation means 5 and the addition means 6 are operated in the coefficient correction means 7. based on the already pressure setpoint Psv '(t) and pressure process value Ppv (t), is configured to correct by learning the coefficients a n and b n of the ARX model. In the learning of the coefficient, the model prediction value calculated at the time point t 0 is represented as y p (t), and the pressure process value y (t) is used as a teacher signal within the range of t 0 ≦ t ≦ t 0 + T. Predicted value y p
Correct the coefficients so that (t) is improved. Now, the evaluation function

【数12】 J=Σ{y(t0+s)−yp(t0+s)}2/2 …(28) ただし、Σ:s=0〜Tについての加算 とすれば、係数an及びbnの補正代△an及び△bnは、
以下のように表される。
Equation 12] J = Σ {y (t 0 + s) -y p (t 0 + s)} 2/2 ... (28) However, sigma: if sum of s = 0 to T, the coefficient a n and b n of correction fee △ a n and △ b n is,
It is expressed as follows.

【数13】 △an=−ε∂J/∂an =εΣ{y(t0+s)−yp(t0+s)}pn(t0+s) …(29) △bn=−ε∂J/∂bn =εΣ{y(t0+s)−yp(t0+s)}qn(t0+s) …(30) ただし、ε:学習率 Σ:s=0〜Tについての加算 pn(t)=∂yp(t)/∂an (1≦n≦na) qn(t)=∂yp(t)/∂bn (1≦n≦nb) なお、実際には、複数の教師データの組から逐次1組を
選択し、予め設定した収束基準を満たすまで、係数の学
習を繰り返す。上記したモデル予測制御においてはAR
Xモデルについて説明したが、他の数式モデルを用いて
もよいことは言うまでもない。
Equation 13] △ a n = -ε∂J / ∂a n = εΣ {y (t 0 + s) -y p (t 0 + s)} p n (t 0 + s) ... (29) △ b n = - ε∂J / ∂b n = εΣ {y (t 0 + s) -y p (t 0 + s)} q n (t 0 + s) ... (30) However, epsilon: learning rate sigma: for s = 0 to t summing p n (t) = ∂y p (t) / ∂a n (1 ≦ n ≦ n a) q n (t) = ∂y p (t) / ∂b n (1 ≦ n ≦ n b) Actually, one set is sequentially selected from a plurality of sets of teacher data, and learning of coefficients is repeated until a preset convergence criterion is satisfied. In the above model predictive control, AR
Although the X model has been described, it goes without saying that another mathematical model may be used.

【0020】図4は、本発明の第2の実施例の燃焼制御
方法を実機テスト(圧力設定値Psv(t)=一定)によ
り実行した結果を、圧力設定値を補正しない従来例の燃
焼制御方法のテスト結果と対比して表したグラフであ
る。図4において、(A)はボイラ出口におけるボイラ
圧力を、(B)はボイラ蒸気流量を示しており、時間T
0〜T1及びT2〜T3は従来例により得られた圧力プ
ロセス値を、時間T1〜T2は本発明により得られた圧
力プロセス値を示している。これらのグラフから明らか
なように、本発明によれば、従来例に比べて、圧力プロ
セス値の設定値からの変動を極めて小さくでき、しかも
ボイラ蒸気流量の変動も抑圧できるので、安定した発電
量を提供できることが実証された。
FIG. 4 shows the result of executing the combustion control method according to the second embodiment of the present invention by an actual machine test (set value of pressure Psv (t) = constant). Fig. 3 is a graph showing the test results of the method in comparison with the test results. 4A shows the boiler pressure at the boiler outlet, FIG. 4B shows the boiler steam flow rate, and the time T
0 to T1 and T2 to T3 indicate pressure process values obtained by the conventional example, and times T1 to T2 indicate pressure process values obtained by the present invention. As is clear from these graphs, according to the present invention, the fluctuation from the set value of the pressure process value can be extremely reduced and the fluctuation of the steam flow rate of the boiler can be suppressed as compared with the conventional example. It has been demonstrated that

【0021】本発明は以上のように構成されているの
で、ボイラ出口圧力プロセス値と圧力設定値との偏差に
応じて設定値を補正することにより、ボイラ圧力の変動
を十分に低減することができるので、プラントの安全性
及び信頼性を向上させることができる。さらに、これに
よりボイラ蒸気流量を安定化させることができるので、
発電量の安定化及び効率化をはかることができる。ま
た、ボイラ圧力設定値を補正するための補正量の決定に
モデル予測制御を採用すれば、システムの変動周期とボ
イラ圧力設定値の変化に対するプロセス値の応答時間と
の比が比較的小さい場合にも、対処することが可能とな
り、ボイラ圧力及び蒸気流量を安定化させる等の制御性
能が向上する。さらにまた、モデル予測制御におけるモ
デル係数の補正機能を備えれば、焼却炉の経時変化等に
拘わらず、制御性能を悪化させることなく発電量を安定
化することができる。
Since the present invention is configured as described above, it is possible to sufficiently reduce the fluctuation of the boiler pressure by correcting the set value according to the deviation between the boiler outlet pressure process value and the pressure set value. Therefore, the safety and reliability of the plant can be improved. In addition, this can stabilize the boiler steam flow,
The amount of power generation can be stabilized and efficiency can be improved. In addition, if model predictive control is adopted to determine the correction amount for correcting the boiler pressure set value, the ratio between the system fluctuation cycle and the response time of the process value to the change in the boiler pressure set value is relatively small. Can be dealt with, and the control performance such as stabilizing the boiler pressure and the steam flow rate is improved. Furthermore, if a function of correcting the model coefficient in the model predictive control is provided, the power generation amount can be stabilized without deteriorating the control performance irrespective of the aging of the incinerator.

【図面の簡単な説明】[Brief description of the drawings]

【図1】本発明の第1及び第2の燃焼制御方法を説明す
るためのブロック図である。
FIG. 1 is a block diagram for explaining first and second combustion control methods of the present invention.

【図2】本発明の第2の実施例の燃焼制御方法における
時間ファクタの関係を説明するための図である。
FIG. 2 is a diagram for explaining a time factor relationship in a combustion control method according to a second embodiment of the present invention.

【図3】本発明の第3の実施例の燃焼制御方法を説明す
るためのブロック図である。
FIG. 3 is a block diagram for explaining a combustion control method according to a third embodiment of the present invention.

【図4】本発明の第2の実施例の燃焼制御方法を用いて
実機テストを行った結果を、従来例の結果と対比して表
したグラフである。
FIG. 4 is a graph showing a result of an actual machine test performed using the combustion control method according to the second embodiment of the present invention in comparison with a result of a conventional example.

【図5】従来例の燃焼制御方法を説明するためのブロッ
ク図である。
FIG. 5 is a block diagram for explaining a conventional combustion control method.

Claims (5)

【特許請求の範囲】[Claims] 【請求項1】 ボイラ付き焼却プラントにおける燃焼制
御方法において、 ボイラ出口におけるボイラの圧力プロセス値を周期的に
測定するステップと、 測定により得られた圧力プロセス値と圧力設定値との偏
差に応じて圧力設定値を補正することにより、補正済圧
力設定値を得る設定値補正ステップと、 得られた補正済圧力設定値を圧力設定値の代わりにボイ
ラ圧力調節計に供給するステップと、 ボイラ圧力調節計の出力に基づいて、焼却炉への燃料の
供給量を調節する燃料調節ステップとを含むことを特徴
とする燃焼制御方法。
1. A method for controlling combustion in a boiler-equipped incineration plant, comprising the steps of: periodically measuring a boiler pressure process value at a boiler outlet; A set value correction step of obtaining a corrected pressure set value by correcting the pressure set value; a step of supplying the obtained corrected pressure set value to the boiler pressure controller instead of the pressure set value; A fuel adjusting step of adjusting the amount of fuel supplied to the incinerator based on the output of the meter.
【請求項2】 請求項1記載の燃焼制御方法において、
設定値補正ステップは、 ボイラ出口におけるボイラの圧力設定値を入力としかつ
圧力プロセス値を出力として予め作成された数式モデル
に基づいて、モデル予測制御を行うことにより、圧力設
定値を補正して補正済圧力設定値を得るステップからな
ることを特徴とする燃焼制御方法。
2. The combustion control method according to claim 1, wherein
The set value correction step corrects and corrects the pressure set value by performing model predictive control based on a mathematical model created in advance with the boiler pressure set value at the boiler outlet as input and the pressure process value as output. A method for controlling the combustion pressure, the method comprising: obtaining a set pressure value.
【請求項3】 請求項2記載の燃焼制御方法において、
設定値補正ステップはさらに、 補正済圧力設定値と測定により得られた圧力プロセス値
とを用いて、数式モデルの係数を補正するステップを含
んでいることを特徴とする燃焼制御方法。
3. The combustion control method according to claim 2, wherein
The combustion control method, further comprising the step of correcting the coefficient of the mathematical model using the corrected pressure set value and the pressure process value obtained by the measurement.
【請求項4】 請求項2又は3記載の燃焼制御方法にお
いて、予め作成された数式モデルは、ARXモデルであ
ることを特徴とする燃焼制御方法。
4. The combustion control method according to claim 2, wherein the mathematical model created in advance is an ARX model.
【請求項5】 請求項1〜5いずれかに記載の燃焼制御
方法において、燃料調節ステップは、 ボイラ圧力調節計の出力に基づいて、焼却炉へ廃棄物を
供給する給塵供給手段の駆動回数を調節するステップで
あることを特徴とする燃焼制御方法。
5. The combustion control method according to claim 1, wherein the fuel adjusting step includes the step of driving the dust supply means for supplying waste to the incinerator based on the output of the boiler pressure controller. Controlling the combustion.
JP10670298A 1998-04-16 1998-04-16 Combustion control method in incineration plant with boiler Pending JPH11294741A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP10670298A JPH11294741A (en) 1998-04-16 1998-04-16 Combustion control method in incineration plant with boiler

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP10670298A JPH11294741A (en) 1998-04-16 1998-04-16 Combustion control method in incineration plant with boiler

Publications (1)

Publication Number Publication Date
JPH11294741A true JPH11294741A (en) 1999-10-29

Family

ID=14440344

Family Applications (1)

Application Number Title Priority Date Filing Date
JP10670298A Pending JPH11294741A (en) 1998-04-16 1998-04-16 Combustion control method in incineration plant with boiler

Country Status (1)

Country Link
JP (1) JPH11294741A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014013093A (en) * 2012-07-03 2014-01-23 Takuma Co Ltd Method for controlling air-cooled wall surface temperature for combustion furnace, and combustion furnace using the same

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014013093A (en) * 2012-07-03 2014-01-23 Takuma Co Ltd Method for controlling air-cooled wall surface temperature for combustion furnace, and combustion furnace using the same

Similar Documents

Publication Publication Date Title
US6496749B1 (en) Semiconductor producing apparatus and temperature control method therefor
JP3775753B2 (en) Model-based temperature controller for semiconductor thermal processing machines
US6207936B1 (en) Model-based predictive control of thermal processing
US6278899B1 (en) Method for on-line optimization of a plant
US4768143A (en) Apparatus and method using adaptive gain scheduling algorithm
JP3949164B2 (en) A feedback method for controlling nonlinear processes.
US20030014131A1 (en) Method for optimizing a plant with multiple inputs
EP0788626A1 (en) A variable horizon predictor for controlling dead time dominant processes, multivariable interactive processes, and processes with time variant dynamics
WO1997028669A9 (en) Model-based predictive control of thermal processing
US6711531B1 (en) Temperature control simulation method and apparatus
JP4112561B2 (en) Predictive control device, predictive control method, predictive control program, and storage medium
CN110618706A (en) Multistage intelligent denitration online optimization control system based on data driving
JP3061450B2 (en) Model predictive controller
JPH11294741A (en) Combustion control method in incineration plant with boiler
JPH0751202B2 (en) Computer-controlled spray drying method
JP2021152702A (en) Apparatus for assisting plant-operation optimization, and apparatus and method for controlling plant-operation optimization
JP4381628B2 (en) Concentration control device and concentration control method for concentration equipment
JP3653599B2 (en) Apparatus and method for controlling ammonia injection amount of flue gas denitration equipment
JP2010218007A (en) Disturbance estimation device, control object model estimation device, feedforward amount estimation device, and controller
KR100221231B1 (en) Model-based predictive control system and method
JPH07332602A (en) Steam temperature prediction control device
JPH07167402A (en) Automatic controller for thermal power plant
JPH06222808A (en) Automatic controller for thermoelectric power plant
JP3124872B2 (en) Control equipment for thermal power plant
JP3448823B2 (en) Power plant control equipment