JP2001062246A - 脱硝制御方法およびその制御装置 - Google Patents

脱硝制御方法およびその制御装置

Info

Publication number
JP2001062246A
JP2001062246A JP24012399A JP24012399A JP2001062246A JP 2001062246 A JP2001062246 A JP 2001062246A JP 24012399 A JP24012399 A JP 24012399A JP 24012399 A JP24012399 A JP 24012399A JP 2001062246 A JP2001062246 A JP 2001062246A
Authority
JP
Japan
Prior art keywords
denitration
nox
data
amount
denitration process
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
Application number
JP24012399A
Other languages
English (en)
Other versions
JP4057750B2 (ja
Inventor
Yoshiko Shimizu
佳子 清水
Masashi Nakamoto
政志 中本
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.)
Toshiba Corp
Original Assignee
Toshiba 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 Toshiba Corp filed Critical Toshiba Corp
Priority to JP24012399A priority Critical patent/JP4057750B2/ja
Publication of JP2001062246A publication Critical patent/JP2001062246A/ja
Application granted granted Critical
Publication of JP4057750B2 publication Critical patent/JP4057750B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Exhaust Gas Treatment By Means Of Catalyst (AREA)
  • Feedback Control In General (AREA)
  • Treating Waste Gases (AREA)

Abstract

(57)【要約】 【課題】 最適化処理を用いて、グループ全体に課
される制約を満たしつつ、アンモニアの消費量を大幅に
低減化することにある。 【解決手段】 各脱硝プロセスの特性・制約条件その他
のデータを設定するパラメータ・条件設定ステップS1
と、この設定された特性データから得られるアンモニア
使用量の2乗を評価関数、前記グループ脱硝プロセス全
体のNOx排出量規制値データを制約条件とし、この制
約条件を満たしつつ前記評価関数を最小化する最適化演
算処理を実行し、前記各脱硝プロセスのNOx配分率を
決定する各脱硝プロセスNOx配分処理ステップS2
と、得られたNOx配分率に基づいて各脱硝プロセスの
制御用設定値を演算出力する設定値計算処理ステップS
3とを有する脱硝制御方法である。

Description

【発明の詳細な説明】
【0001】
【発明の属する技術分野】本発明は、一軸型・多軸型の
複合プラント、一台以上のガスタービン発電設備をもつ
シンプルサイクルプラントその他NOxを発生する各種
の発電プラントに設置される脱硝プロセスに利用される
脱硝制御方法およびその制御装置に関する。
【0002】
【従来の技術】この種の複合プラント、シンプルサイク
ルプラントその他の発電プラントは、燃料の種類に応じ
て、燃焼排ガス中に様々な有害物質が含まれており、大
気に放出する前にその物質を除去する必要がある。その
中でも、特にNOxの除去は非常に重要な技術的課題と
なっている。
【0003】NOxを取り除く最も一般的な方法は、触
媒に気化したアンモニアを吹きかけてアンモニアを活性
化し、NOxと反応させて無害な窒素ガスと水蒸気とに
分解させる方法である。このNOxを分解させる反応を
脱硝反応と言い、この脱硝反応を実現するプロセスを脱
硝プロセスと呼んでいる。
【0004】この脱硝プロセスには次のような特徴をも
っている。
【0005】(1) アンモニアの注入量を増加させる
と、残存するNOx量は減少する。
【0006】(2) アンモニアを多量に注入しても、
NOxを100%取り除くことはできないこと。
【0007】(3) アンモニアと触媒・NOxとの脱
硝反応は化学反応であり、排気ガスの温度・流量・圧力
やアンモニア・NOxの濃度等の変動により反応の特性
が変化する。つまり、非線形性が非常に強い反応系であ
ること。
【0008】(4) 脱硝触媒部の入・出口側NOx濃
度は、常時測定されているか、或いは推定値が計算され
ている。測定の場合には、検出部のハードウエア上の制
約から、非常に大きな無駄時間が発生する(数十秒程
度)。
【0009】一方、以上のような特徴とは別に、発電所
から排出されるNOxについて厳しい制約が存在する。
【0010】すなわち、発電所から排出されるNOx
は、地方自治体などにより様々な種類の制約値が課せら
れているが、大別すると、複数のNOx発生源ないしは
複数の脱硝プロセスに対して課される制約と、個別NO
x発生源ないしは脱硝プロセスに対して課される規制と
に分けられる。
【0011】以下、複数のNOx発生源ないしは複数の
脱硝プロセスをグループと総称し、単独のNOx発生源
ないしは脱硝プロセスを個別プロセスと総称する。
【0012】グループから排出されるNOxの制約は、
グループ全体の総流量あるいは総合的な濃度の制約値と
なるが、それは瞬間的な流量(瞬時値)であったり、あ
る一定時間の流量平均値(移動時間平均値)であったり
する。
【0013】また、個別NOx発生源に対する規制は、
各NOx発生源ごとの流量、濃度、実測値、計算値の何
れかであったり、あるいは瞬時値であったり、さらに移
動時間平均値である場合もある。
【0014】従って、脱硝プロセスは、以上のような特
徴をもつ反面、要求される種々の規制に対処しなければ
ならない。
【0015】そこで、従来、以上のような状況を踏まえ
つつ次のような脱硝制御方法ないし脱硝制御装置が提案
されている。
【0016】発電所から排出されるNOxの制約には、
グループ全体に対して課せられる制約と個別プロセスに
対して課せられる規制とがあり、それぞれに対処するた
めの制御がなされている。
【0017】グループ全体に対する制御は、系列脱硝制
御、グループ脱硝制御、ユニット脱硝制御、統括脱硝制
御などと呼ばれている。なお、発電所全体を1つのグル
ープと見なす場合もあり、発電所内に複数のグループが
存在する場合もある。また、個別NOx発生源に対処す
るための制御が別に行われるのが一般的であり、この場
合にはNOx制御、軸NOx制御、個別脱硝プロセス制
御などと呼ばれている。
【0018】ここで、取り上げる複数の脱硝プロセスに
関する制御は、グループ全体に対する制御を説明する
が、これは系列脱硝制御、統括脱硝制御およびユニット
脱硝制御を含む概念である。
【0019】従来、各ユニットごとの脱硝制御方法およ
びその制御装置に関しては、多数の提案がなされてい
る。例えば特開平10−116105号公報、特開平1
0−15354号公報、特開平9−870号公報、特開
平8−334214号公報、特開平8−281061号
公報、特開平8−257371号公報などである。しか
し、これらの方法及び装置は、系外へ排出されるNOx
量を計測し、この計測値と予め定められる制約値に相当
する設定値との偏差を算出し、その偏差を解消するため
に例えばPID制御器等のごとき補償器を用いて調節演
算を行って操作信号を取出すフィードバック制御系と、
発生したNOx量に応じて予測先行的にアンモニア注入
量を変化させるフィードフォワード制御系とを組合せる
ことにより、系外へ排出するNOx量を予め課せられて
いる規制値を越えないように制御するものである。つま
り、この種の先行技術は、一般的な制御手法を採用し、
排ガス中に含まれるNOx量に応じてアンモニア注入流
量を調整し、系外へ排出するNOx量の低減化を図るも
のである。
【0020】なお、このフィードバック制御系とフィー
ドフォワード制御系との組合せ方法にも種々の種類があ
り、またフィードバック制御部分、フィードフォワード
制御部分は、一般的な予測制御、ファジィ制御であった
り、シミュレータを用いたものなど、様々な制御の提案
が出されている。
【0021】次に、従来のグループ脱硝制御装置として
は、第1の先行技術(特公昭63−62251号公
報)、第2の先行技術(特開平6−257425号公
報)および第3の先行技術(特開平7−163835号
公報)等があり、これらは主としてガスタービンを用い
た発電プラントを対象とする。
【0022】(1) 第1の先行技術によるグループ全
体のNOx排出量制御について。
【0023】このグループ脱硝制御装置は、一台のガス
タービン発電機と一台ないし複数台の排熱回収ボイラ・
蒸気タービンとで構成されている複合発電システム、ま
たは一台以上のガスタービン発電機からなるシンプルサ
イクルなどを備えた発電プラントに適用され、これら各
システムを1ユニットと呼ぶことにする。
【0024】複合発電プラントやガスタービンシンプル
サイクルプラントは、数ユニットずつのグループ毎に取
り扱われる場合が多い。例えば発電出力要求値や有害ガ
スの排出量規制値などは1グループ当りの数値が問題と
なる。このグループを系列と呼ぶ。
【0025】ところで、ガスタービン起動中および停止
中は、排気ガス温度が脱硝反応プロセスに必要な温度よ
りも低いので、脱硝制御が不可能な状態にある。そこ
で、グループ内に起動・停止運転中のユニットが存在す
る場合、予めそれ以外のユニットのNOx排出量を減ら
しておけば、グループ全体のNOx排出量を規制値以下
に抑えることが可能である。
【0026】第1の先行技術は、起動・停止を考慮した
グループ全体の脱硝制御を行うグループ脱硝制御装置で
あり、その構成は図8に示す通りである。
【0027】この制御装置100は、ガスタービン起動
時および停止時の制御不能なNOx排出量パターンをそ
れぞれ記憶するパターン保持器101,102が設けら
れ、さらにグループ全体のNOx排出量設定値(規制
値)が入力され、設定値演算器103は、NOx排出量
検出値であるプロセス量(瞬時値,移動時間平均排出
量)とグループ全体に対して課せられるNOx排出量規
制値である設定値とを比較し、その偏差に基づいてプロ
セス量が規制値を越えないような操作信号を求めた後、
運転軸数信号nのもとに当該操作信号を除算し、各ユニ
ットのNOx排出量設定値を取得する構成である。
【0028】このとき、設定値演算器103は、起動・
停止スケジュール信号から起動・停止と判断されたと
き、パターン保持器101,102のパターンに基づ
き、グループ全体NOx排出量設定値から起動・停止中
のNO制御不可能なユニットから排出されるNOx量を
減算し、その減算結果のNOx量を、残りの制御可能な
ユニットに対して均等にNOx排出量設定値として送出
する。
【0029】なお、ガスタービンの起動・停止のスケジ
ュールは事前に計算され決定されていることが多い。そ
のため予めユニット毎のスケジュールを把握しておき、
現在運転中のユニットのNOx排出量を減らすのに必要
な時間Tだけ遡って、制御不能となるNOx量に基づい
て計算しNOx排出量設定値を得るようにすれば、スケ
ジュール通りの運転が可能となる。
【0030】(2) 第2の先行技術によるグループ全
体のNOx排出量制御について。
【0031】第1の先行技術では時間Tの算出方法が決
められていないが、第2の先行技術では、脱硝プロセス
の動特性モデルを用いることにより、将来の状態を予測
計算し、その計算結果に従って時間Tを決定する技術で
ある。
【0032】(3) 第3の先行技術によるグループ全
体のNOx排出量制御について。
【0033】このグループ脱硝制御装置は、図9に示す
ようにNOx移動時間平均排出量制御部111、NOx
瞬間排出量制御部112、アンモニア瞬間排出量制御部
113およびアンモニア移動時間平均排出量制御部11
4を設け、そのうちNOx移動時間平均排出量制御部1
11には第1軸および第2軸NOx瞬間排出量測定値11
5,116の加算出力117と系列NOx移動時間平均
排出量規制値118とが導入され、NOx瞬間排出量制
御部112には同じく加算出力117と系列NOx瞬間
排出量規制値119とが導入され、さらにアンモニア瞬
間排出量制御部113には第1軸および第2軸アンモニア
瞬間排出量測定値120,121の加算出力122と系
列アンモニア瞬間排出量規制値123とが導入され、ま
たアンモニア移動時間平均排出量制御部114には前記
加算出力122と系列アンモニア移動時間平均排出量規
制値124とが導入され、これら各制御部111〜11
4から移動時間平均および瞬間値によるNOx設定値1
25〜128を取出す。なお、各制御部111〜114
には濃度,流量が導入される場合もある。
【0034】そして、これら設定値125〜128はそ
れぞれ低位優先器130および高位優先器131,13
2で比較することにより、各設定値のうち一番厳しい値
を選択し、運転軸数信号nのもとに、各軸,ひいては各
ユニットのNOx排出量設定値を取出す構成である。
【0035】従って、このグループ脱硝制御装置は、昨
今の環境に対する関心の高まりとともにNOx関係の規
制が多様化されてきているが、これに対処すべく複数の
制約条件として、NOxの瞬時値、移動時間平均値、濃
度、流量などが用いられ、また未反応のまま系外へ放出
されるアンモニア分についても同様に瞬時値、移動時間
平均値、濃度、流量などをそれぞれ管理・規制する制御
技術である。
【0036】
【発明が解決しようとする課題】ところで、以上のよう
な先行技術のうち、例えば第1,第3の先行技術は、グ
ループ内の各脱硝プロセスおよびNOx発生源の特性が
等しいことを前提としているが、実際上,各脱硝プロセ
スやNOx発生源の特性はそれぞれ異なるのが一般的で
ある。
【0037】また、発電プラントにおける運用面の多様
化に伴い、個別ユニットを別な運転条件で運転する必要
性が増えているが、その場合には同じグループ内の各脱
硝プロセスを同じ条件で運転する必然性もない。
【0038】さらに、第2の先行技術は、時間的な応答
変化について考察されているものの、時間的な応答変化
が落ち着いている状態や定常状態等については検討され
ていない。
【0039】従って、これら脱硝制御装置は、個別脱硝
プロセスの特性が異なる場合や個別脱硝プロセスに入っ
てくるNOx流量ないし濃度が異なる場合、系外へ排出
されるNOxが最小になり、かつ、アンモニア消費量が
最小になる点は、従来技術により決定される等配分点で
なくなる。このことは、グループ全体から考えたとき、
過剰なアンモニアが消費されていることを意味し、経済
的に大きな損失となる。
【0040】しかも、アンモニアの注入に際し、余分と
なったアンモニア分がNOxに反応せずに大気中に放出
されるので、環境の面からも好ましくない。
【0041】本発明は上記事情に鑑みてなされたもの
で、脱硝プロセスの特性を有効に活用しつつ、グループ
全体でアンモニア消費量を最小化し、かつ、グループ全
体に課される制約を満たす制御を実現する脱硝制御方法
およびその制御装置を提供することにある。
【0042】また、本発明の他の目的は、脱硝プロセス
・ユニットを含むプラントの運転データを用いて、オフ
ラインまたはオンラインで脱硝プロセスの特性を求める
脱硝制御装置を提供することにある。
【0043】さらに、本発明の他の目的は、脱硝プロセ
スの特性を容易に確認可能とする脱硝制御装置を提供す
ることにある。
【0044】
【課題を解決するための手段】(1) 本発明に係る脱
硝制御方法は、上記課題を解決するために、グループ脱
硝プロセスを構成する複数の脱硝プロセスの制御用設定
値を演算する脱硝制御方法において、各脱硝プロセスの
排出NOx量の総和が予め定めるグループ脱硝プロセス
全体のNOx排出量規制値を越えないとする制約条件を
満たしつつ各脱硝プロセスのアンモニア使用量の2乗を
評価関数を最小化する最適化演算処理を行うことによ
り、各脱硝プロセスのNOx配分率を決定し、これらN
Ox配分率に応じた制御用設定値を得るものである。
【0045】(2) また、本発明に係る脱硝制御方法
は、少なくとも各脱硝プロセスの特性データ、前記グル
ープ脱硝プロセス全体のNOx排出量規制値データを設
定するパラメータ・条件設定ステップと、この設定ステ
ップによって設定される前記特性データから得られるア
ンモニア使用量の2乗を評価関数、前記グループ脱硝プ
ロセス全体のNOx排出量規制値データを制約条件とし
て用い、この制約条件を満たしつつ前記評価関数を最小
化する最適化演算処理を実行し、前記各脱硝プロセスの
NOx配分率を決定する各脱硝プロセスNOx配分処理
ステップと、このNOx配分処理ステップによって得ら
れるNOx配分率に基づいて前記各脱硝プロセスの制御
用設定値を演算出力する設定値計算処理ステップとを有
する脱硝制御方法である。
【0046】この発明は、以上のような手段を講じたこ
とにより、グループ脱硝プロセスから排出されるNOx
量を適切に保持しつつアンモニア消費量を最小化するこ
とが可能である。
【0047】(3) 前記各脱硝プロセスNOx配分処
理ステップとしては、NOx発生量に関する計算開始条
件値を設定した後、各脱硝プロセスごとにアンモニア使
用量を増加させたときのNOx減少変化分を求め、この
NOx減少変化分の大きさに応じてNOx配分率を修正
し、この各脱硝プロセスの修正NOx発生量の総和が前
記グループ脱硝プロセス全体のNOx排出量規制値デー
タより小さいとき、それらの修正NOx発生量に対応す
るNOx配分率を出力するようにすれば、各脱硝プロセ
スの能力を考慮しつつアンモニア消費量を最小化しつつ
NOxを抑制した状態で運転することが可能となる。
【0048】(4) さらに、本発明に係る脱硝制御装
置は、前記(1)項、(2)項および(3)項の制御方
法を実施する手段を備えた構成とすることもできる。
【0049】例えば前記(2)項の制御方法を実施する
脱硝制御装置としては、少なくとも各脱硝プロセスの特
性データ、前記グループ脱硝プロセス全体のNOx排出
量規制値データを設定するパラメータ・条件設定手段
と、この設定手段によって設定される前記特性データか
ら得られるアンモニア使用量の2乗を評価関数、前記グ
ループ脱硝プロセス全体のNOx排出量規制値データを
制約条件として用い、この制約条件を満たしつつ前記評
価関数を最小化する最適化演算処理を実行し、前記各脱
硝プロセスのNOx配分率を決定する各脱硝プロセスN
Ox配分処理手段と、このNOx配分処理手段によって
得られるNOx配分率に基づいて前記各脱硝プロセスの
制御用設定値を演算出力する設定値計算処理手段とを備
えた構成とする。
【0050】(5) さらに、本発明に係る脱硝制御装
置は、前記(4)項に記載する構成要素に新たに、各脱
硝プロセスを含むプラントの運転データ、運転開始前ま
たは運転中の運転データから求められる前記各脱硝プロ
セスの特性データ等を記憶するデータ記憶手段と、この
データ記憶手段に記憶される特性データの中から同一の
脱硝プロセスの過去特性データ及び前記運転データから
得られる今回特性データとを読み出して比較表示し、或
いは複数の脱硝プロセスの特性データを読み出して比較
表示し、或いは種々の異なる運転データであるパラメー
タを考慮して得られた脱硝プロセスの特性データを読み
出して二次元または三次元にて比較表示する表示制御手
段とを設けた構成とすることにより、脱硝プロセスの特
性状態を一目瞭然に把握でき、監視,制御、さらに経済
性の観点から有効なものとなる。
【0051】さらに、予め運転データの中の少なくとも
アンモニア注入量とNOx発生量の第1のモル比とNO
x排出量と前記NOx発生量との第2のモル比との関係
の近似式データを記憶し、入力される第1または第2の
モル比に基づき、前記近似式により前記第2または第1
のモル比を順次計算しテーブル化し、そのテーブル内容
を表示してもよい。
【0052】なお、前記各脱硝プロセスの特性データ
は、アンモニア注入量、NOx発生量およびNOx排出
量以外に種々のプロセス量を用いて求めることができ
る。
【0053】
【発明の実施の形態】以下、本発明に係る脱硝制御方法
およびその制御装置の実施の形態について図面を参照し
て説明する。
【0054】(第1の実施の形態)図1は本発明に係る
脱硝制御方法の一実施の形態を説明する図である。
【0055】この脱硝制御方法は、例えば複数のユニッ
トにそれぞれ設置される各脱硝プロセスからなるグルー
プ脱硝プロセスに適用するグループ脱硝制御であって、
ここで言うグループ脱硝制御とは、前述したように系列
脱硝制御、ユニット脱硝制御、統括脱硝制御、さらには
単独のNOx発生源であるにも拘らず、複数に分岐され
る煙道にそれぞれ脱硝プロセスが設置されるプラントの
ユニット脱硝制御などを含む概念である。
【0056】このグループ脱硝制御方法は、各個別脱硝
プロセスに対応する軸制御装置に与える最適な設定値を
演算し、得られた設定値を各軸制御装置に設定する方法
であり、大きく分けると図1に示すように、個別脱硝プ
ロセスに必要な特性及び条件などを設定するパラメータ
・条件データ設定ステップS1と、これらパラメータお
よび条件設定データのもとに最適化処理によってグルー
プから排出されるNOx量を適切な状態に抑制し、か
つ、アンモニア使用量を最小化するために、各脱硝プロ
セスに対するNOx配分率を求める各脱硝プロセスNO
x配分処理ステップS2と、各脱硝プロセスに対するN
Ox配分率に基づいて、各個別脱硝プロセスに対する設
定値を計算処理する設定値計算処理ステップS3とから
なり、以下、各ステップについて詳細に説明する。
【0057】A. パラメータ・条件データ設定ステッ
プS1について このデータ設定ステップS1は個別脱硝プロセスの特
性、条件などを設定する処理ステップであって、個別脱
硝プロセスの特性は図1に示す如く例えばNOx排出量
NOxoutとNOx発生量NOxinとのモル比x=NO
xout/NOxinを横軸とし、アンモニア注入量NH3in
とNOx発生量NOxinとのモル比y=NH3in/NO
xinを縦軸とする,いわゆる各軸脱硝反応特性が設定さ
れる。モル比に代えて濃度比または体積流量比であって
もよい。
【0058】また、同図においてf1(x)〜fn(x)は各
脱硝プロセスの特性を表し、nは個別脱硝プロセスの数
であって、1≦i≦nの関係にある。その他、グループ
全体のNOx排出量設定値r、ユニット・脱硝プロセス
を含むプラントの運転データであるプロセス量(NOx
発生量,排ガス流量その他多くの物理量を含む)および
追加制約条件などが入力設定される。なお、グループ全
体のNOx排出量設定値rは規制値に相当するものであ
って、規制の種類は特に限定されるものでなく、例えば
流量,濃度,実測値,計算値,或いは瞬時値,移動時間
平均値の何れであっても、有効性を損なわない限りにお
いて全て包括する概念である。また、追加制約条件とし
ては例えば各ユニットから排出されるNOx量の最小値
などが挙げられる。
【0059】従って、パラメータ・条件データ設定ステ
ップS1は、以上述べたように最適化処理に必要な各種
のデータを設定し、各脱硝プロセスNOx配分処理ステ
ップS2に入力する。
【0060】B. 各脱硝プロセスNOx配分処理ステ
ップS2について この各脱硝プロセスNOx配分処理ステップS2は、パ
ラメータ・条件データ設定ステップS1から与えられる
パラメータ,条件データなどを用いて、n個のユニッ
ト、ひいてはn個の脱硝プロセスに対するNOxの配分
を決定する処理であって、次のような最適化処理を行
う。
【0061】以下、機能を分かりやすく説明するため
に、例えば各ユニットから発生するNOx発生量が同一
であり、各脱硝プロセスの特性が異なる場合の例につい
て説明する。
【0062】b10.仮定 説明を簡単にするために、以下の仮定を満たすグループ
脱硝プロセスの最適化処理の例を対象とする。
【0063】b11:グループ脱硝プロセス全体のNO
x排出量に対する瞬時値の流量は制約条件として課せら
れている。
【0064】b12:脱硝プロセスの特性は、図1に示
すようにNOxのモル比とアンモニアのモル比のみによ
って決まること。
【0065】b13:追加制約条件は、各ユニットから
排出されるNOxの最小値oのみとする。つまり、各ユ
ニットから排出されるNOx量oiは最小値oよりも大
きくなくてはならないこと。このことは、o=0の場合
にはどんなに多量のアンモニアを注入しても、NOx排
出量は負にならないことを意味する。ここで、oiはi
番目の脱硝プロセスから排出されるNOxの体積流量を
意味する。
【0066】b14:各ユニット・各脱硝プロセスを含
むプラントのプロセス量は、例えばNOx発生量および
排ガス流量であるとする。
【0067】b15:各脱硝プロセスの配分計算は、非
線形計画法を用いる。
【0068】b16:配分処理結果、最適なアンモニア
とNOxとのモル比が得られるものとする。
【0069】b17:各ユニットの脱硝制御では、排出
されるNOx流量設定値が必要であるとする。
【0070】b18:脱硝プロセスn=3とする。
【0071】b20.物理的な関係式 今、i番目の脱硝プロセスで使用するアンモニアの体積
流量hi、i番目の脱硝プロセスの触媒入口NOx体積
流量gi、i番目の脱硝プロセスから排出されるNOx
の体積流量oi、i番目の脱硝プロセスの特性xi,y
iとすると、これらの間には次のような関係が成立する
ものとする。
【0072】 yi=oi/gi ……(1) ∴oi=yi*gi ……(1a) xi=hi/gi ……(2) ∴hi=xi*gi ……(2a) これとは別に、fi(x)=aioi2+bioi+cと前記(2)式とから、 ∂fi/∂hi=2(a/gi2)*oi+bi/gi ……(3) なる式が得られる。
【0073】b30.最適化問題の定式化今、最適化の
指標である評価関数Jは、 J=h12+h22+h32 ……(4) なる式で表わす。これは3個の脱硝プロセスに対するア
ンモニア使用量の2乗の和をもって表わされる。
【0074】一方、制約条件は、グループ脱硝プロセス
全体のNOx排出量設定値(規制値)がr、各脱硝プロ
セス1〜3から排出されるNOxの体積流量がo1,o
2,o3であるので、 o1+o2+o3≦r ……(5) の関係になければならない。但し、 o1≧0、o2≧0、o3≧0 ……(5a) である必要がある。すなわち、グループ脱硝制御の目的
は、制約条件を満たしつつ評価関数を最小化する排出N
Oxの体積流量o1,o2,o3を決定することにあ
る。
【0075】b40.最適な配分を得るアルゴリズム b41:このアルゴリズムの実行に際し、前記(5)式
を満足させるための最初の値(計算開始条件値)o1
(0),o2(0),o3(0)を適当に定めて設定す
る必要がある。ここでは、制約条件に基づき、o1
(0)=o2(0)=o3(0)=(1/3)rを計算
開始条件値として設定する。
【0076】b42:Jの∂J/∂o(変数に対する傾
きないし変化分)を計算する。つまり、ここでは、以上
のような計算開始条件値のもとに、各脱硝プロセスに対
してアンモニアの注入量を少しずつ増加させたときの各
軸ないし各脱硝プロセスのNOx発生量がどの程度減少
したかを表わす変化分を演算するものであって、下記す
る(6)式によって計算する。
【0077】
【数1】
【0078】b43:各脱硝プロセスのNOx流量oi
の更新 ここでは、計算開始条件値のもとにアンモニアを少しず
つ増加注入させたときの各脱硝プロセスのNOx減少傾
きないし減少変化分を求めた後、このNOx減少変化分
の大きい脱硝プロセスには大きなNOx配分率、NOx
減少変化分の小さい脱硝プロセスには小さいNOx配分
率を求めて、下式に従って次回の修正NOx流量∧(ハ
ット)oi(1)(但し、i=1,2,3)を求める。
【0079】
【数2】
【0080】すなわち、この(7a)式、(7b)式、
(7c)式は、制約条件(5)式及び(5a)式を満た
すように、最小値との比較を行って∧(ハット)oi
(1)を求めている。ここで、αは予め決めた値ないし
直線探索アルゴリズムにより決定される係数である。
【0081】b44:制約条件の確認と条件不可時の処
理 以上のような(7a)式〜(7c)式によって修正され
た∧(ハット)oi(1)が制約条件(5)式を満たし
ているかどうかを確認する。ここで、制約条件を満たし
ている場合にはその修正値であるNOx流量oi(1)
であるNOx配分率を採用し、満たしていない場合には
下式により∧(ハット)oi(1)を修正しNOx流量
oi(1)を求める。
【0082】
【数3】
【0083】この(8a)式〜(8c)式から得られる
値は各脱硝プロセスに対して排出してもよいNOx量o
1(1)〜o3(1)であって、NOx配分率w1,w
2、w3に相当する値となる。
【0084】b45:収束判定 さらに、計算終了条件を満足したか否かを判定し、満足
していない場合には前記b42〜b44の処理を繰り返
し行う。計算終了条件としては、予め定めた計算回数を
越えた場合、oi(k+1)とoi(k)との大きさの
差が予め定めた値よりも小さくなった場合、oi(k+
1)とoi(k)とを使ってJの値を計算し、その差が
予め定めた値よりも小さくなった場合などが挙げられ
る。なお、何れの条件にも該当しない場合は、oi
(1)を初期値としてd43に戻り、d45まで繰り返
し、各脱硝プロセスのNOx配分率を求める。kは処理
回数を意味する。なお、図1においてxaiはi番目の
脱硝プロセスに課される制約値である。
【0085】C. 設定値計算処理ステップS3につい
て 以上のようにして各脱硝プロセスのNOx配分率w1〜
w3=o1〜o3(o1〜o3:体積流量)を求めた後、各
脱硝プロセスの制御装置15の設定値を求めるが、これ
らの設定値sv1,sv2,sv3は各脱硝プロセスの
制御装置15で何の設定値を必要とするかによって異な
る。以下、異なる2つの方法について述べる。
【0086】イ.各脱硝プロセスの制御装置15で使用
しているNOx制御設定値が体積流量に関するものであ
れば、o1〜o3をそのままsv1〜sv3として使用す
る。
【0087】ロ.各脱硝プロセスの制御装置15で使用
しているNOx制御設定値がモル比=(アンモニア注入
量/NOx発生量)→sv1〜sv3である時、各脱硝
プロセスの特性関数を用いてo1〜o3→sv1〜sv3
を計算する。
【0088】その他、各脱硝プロセスに対するNOx配
分率w1,w2,w3と対応するユニット・脱硝プロセ
スを含むプラントのプロセス量とから、個別脱硝プロセ
スを制御する軸制御装置の設定値sv1〜sv3を計算
処理し、該当軸制御装置に設定する。
【0089】従って、以上のような実施の形態によれ
ば、アンモニア注入量を少しずつ増やしながら、制約条
件を満たしつつ各脱硝プロセスのNOx減少変化分を計
算し、減少変化分が大きい脱硝プロセスほどNOx配分
率を大きくし、このNOx配分率に従って脱硝プロセス
を制御する制御装置のアンモニア注入設定値を決定する
ので、図2(a)に示すようにNOx量oから決まる各
脱硝プロセスのNOx最小発生量x1,x2,x3とし
つつ、最小のアンモニア注入量ないし使用量y1,y
2,y3を用いて、各脱硝プロセスを制御できる。
【0090】これに対し、図8に示す脱硝制御装置で
は、図2(b)に示すような脱硝制御がなされている。
すなわち、従来のNOx排出流量設定値の等配分による
制御について述べる。従来の制御方法では、個別脱硝プ
ロセスの特性を考慮せずに各軸に同じ設定値o10=o
20=o30=(1/3)*o0=(1/3)*rを出
力する。このとき、o0=rであり、グループ全体のN
Ox排出量は課せられた値を満たしている。この時のx
10,x20,x30は、g1=g2=g3=gと仮定
すると、 x10=o10/g=o0/3g ……(9a) x20=o20/g=o0/3g ……(9b) x30=o30/g=o0/3g ……(9c) から、 ∴x10=x20=x30 ……(9d) が得られる。
【0091】このときの全体のアンモニア使用量y0
は、 y0=g*(y10+y20+y30) ……(10) である。
【0092】本発明では、NOx総排出量oは、 o=g*(x1+x2+x3)=r ……(11) であって、例えば従来と同じであるにも拘らず、全体のアンモニア使用量yは、 y=g*(y1+y2+y3) ……(12) であり、この例では明らかにy〈y0の関係となる。す
なわち、本発明による方法は、より少ないアンモニア消
費量を用いてNOxの低減化を達成できることを意味す
る。
【0093】次に、この脱硝制御方法の他の実施形態例
について説明する。
【0094】この他の実施の形態は、上記実施の形態と
は逆に個別発生源に応じてNOx発生量が異なるが、脱
硝プロセスの特性が同一となる場合の例である。なお、
ここでは、特に各ユニットNOx配分処理ステップS2
のみについて説明する。
【0095】b10:仮定 説明を簡単に行うためのb11〜b18は前記B項で述
べたと同等の仮定を満たす脱硝プロセスについての最適
化の例である。
【0096】b20:物理的な関係式 前記B項で述べた前記(1)式〜(3)式の関係式が成
立する。
【0097】b30:最適化問題の定式化 前記B項で述べた通りであり、評価関数は(4)式、制
約条件は前記(5)式及び(5a)式とする。
【0098】b40:最適な配分率を得るアルゴリズム b41:計算開始の条件値は前記B項で述べたのと同様
である。
【0099】b42:変化分の計算,つまりJの∂J/
∂o(変数に対する傾き)の計算は前記B項で述べたの
と同じ処理である。この実施の形態では、個別脱硝プロ
セスの特性が等しいと仮定しているので、関数全体で
は、f1=f2=f3となる。
【0100】しかし、一般的には、x1(k),x2
(k),x3(k)は等しくないので、f1{x1
(k)}≠f2{x2(k)}≠f3{x3(k)}で
あるので、前記B項のb42で述べた計算式をそのまま
使用する。
【0101】b43〜b45についても前記B項で述べ
たのと同じ処理となる。
【0102】従って、この実施の形態によれば、前述と
同様に最適化処理により、アンモニアを少しずつ増やし
たときに制約条件を満たしつつNOx変化分の大きな脱
硝プロセスに大きなNOx配分率を設定することによ
り、図3(a)に示すような制御を行うことが出来るの
に対し、従来方法では、各脱硝プロセスに対してNOx
排出流量設定値を等配分としたものであって、同図
(b)に示すような脱硝制御となる。
【0103】以下、効果を明らかにするために、具体的
に数値をもって説明する。
【0104】今、g1=100、g2=70、g3=5
0、r=20とする。本発明方法では、最適化処理によ
り異なるxnの値を決める必要がある。例えば本例で
は、最適値は等モル比配分となるので、x1=x2=x
3=0.1が得られ、その結果、 o=g1x1+g2x2+g3x3 =(g1+g2+g3)x1 =220x1=20 となり、グループ全体のNOx排出量はr=20とな
る。このとき、f1=f2=f3により、全体のアンモ
ニア使用量yは、 y=g1y1+g2y2+g3y3 =(g1+g2+g3)y1 =220*1.0=220である。
【0105】これに対し、従来方法においては、個別脱
硝プロセスの特性を考慮に入れずに各軸に同じ設定値と
してr/3を出力する。このときのNOx総排出量o0
は、g1=100、g2=70、g3=50、r=20
とすると、 100x10=70x20=50x3=20/3 ∴x10=0.0666 x20=0.0952 x30=0.1333 である。一方、全体のアンモニア使用量yoは、 y1=1.5、y2=1.05、y3=0.7から、 yo=g1y10+g2y20+g3y30 =100*1.5+70*1.05+50*0.7 =150+73.5+35 =258.5 となる。
【0106】以上の例からも明らかなように、y〈yo
の関係となり、本発明方法では、より少ないアンモニア
消費量を用いて、NOxの低減化を達成できる。
【0107】(第2の実施の形態)図4は本発明に係る
脱硝制御装置の一実施の形態を説明する構成図である。
【0108】この実施の形態は、グループ脱硝プロセス
10と、このグループ脱硝プロセス10を制御するため
の設定値を与えるためのグループ脱硝制御装置20とか
らなる。ここで、グループ脱硝制御とは、第1の実施の
形態で述べたように、系列脱硝制御、ユニット脱硝制
御、統括脱硝制御、さらには単独のNOx発生源に対し
て複数に分岐される煙道にそれぞれ脱硝プロセスが設置
されているプラントのユニット脱硝制御などを含む概念
である。
【0109】このグループ脱硝プロセス10は、複合プ
ラント、シンプルサイクルプラントその他NOxを発生
する発電プラント等の複数のユニットと複数の個別脱硝
プロセスとからなり、各ユニットにそれぞれ個別脱硝プ
ロセスが設置されている。
【0110】各ユニットの一例は、例えば燃焼器11、
ガスタービン12、燃料の燃焼によって発生する排気ガ
スを導く排気ダクト13および排気ガスを大気に放出す
る煙突14その他必要な機器(図示せず)等によって構
成されている。
【0111】各個別脱硝プロセスは、グループ脱硝制御
装置20から設定される設定値svi(i=1,2,
…,n)のもとに所要の調節演算を実行し弁操作信号を
取出す軸制御装置15、排気ダクト13内に設置される
アンモニア噴霧器16、この噴霧器16下流側に設置さ
れる脱硝触媒17、軸制御装置15からの弁操作信号を
受けて弁開度が調整されて所要量のアンモニアを注入す
るアンモニア注入弁18a、この注入弁18aを介して
注入されるアンモニアと空気ブロア18bから送り込ま
れてくる空気とを混合しアンモニア噴霧器16に供給す
る混合器19等によって構成され、このアンモニア噴霧
器16から排気ダクト13内の脱硝触媒17に吹きかけ
ることにより、排気ガス中のNOx発生量を低減する構
成となっている。
【0112】前記グループ脱硝制御装置20は、本発明
の要部となる部分であって、個別脱硝プロセスに必要な
特性及び条件などのデータを設定するパラメータ・条件
データ設定手段21と、これらパラメータおよび条件設
定データのもとに最適化処理によってグループ脱硝プロ
セス全体から排出されるNOx量を適切な状態に抑制
し、かつ、アンモニア使用量を最小化するために、各個
別脱硝プロセスに対するNOx配分率を求める各脱硝プ
ロセスNOx配分処理手段22と、このNOx配分率に
基づき、各個別脱硝プロセスの軸制御装置16に与える
設定値を計算処理する設定値計算処理手段23とによっ
て構成されている。
【0113】次に、以上のような装置の動作について説
明する。
【0114】先ず、パラメータ・条件データ設定手段2
1は、第1の実施の形態におけるパラメータ・条件デー
タ設定ステップS1で説明したように、各脱硝プロセス
の特性f1(x),f2(x),…,fn(x)、グル
ープ全体のNOx排出量設定値(規制値)rその他の条
件データ等を設定し、各脱硝プロセスNOx配分処理手
段22に入力する。
【0115】この各脱硝プロセスNOx配分処理手段2
2は、CPUで構成され、予め定める図5に示すプログ
ラムに従って所定の処理を実行する。先ず、最適化処理
のための仮定や最適化問題の定式化は第1の実施の形態
と同様であるとすれば、前記(5)式、(5a)式の制
約条件を満たしつつ前記(4)式の評価関数Jを最小化
するNOx排出量o1,o2,o3を決定する必要があ
るが、そのためには計算開始の点,つまり計算開始の条
件値設定する(S11)。制約条件は、グループ全体の
NOx排出量設定値(規制値)がr、各脱硝プロセス1
〜3に対する排出NOxの体積流量がo1,o2,o3
とすれば、o1+o2+o3≦rであるので、計算開始
条件値は下式の値を設定することになる。
【0116】 o1(0)=o2(0)=o3(0)=(1/3)r この状態においてアンモニア注入量を少しずつ増やして
いき、各脱硝プロセスのNOx量がどの程度減少するか
のNOx変化分を計算する(S12)。具体的には前記
(6)式に基づいて演算を実行する。なお、ステップS
12において変数に付けたkはk回目に得られた値(変
数に対する傾き)であることを表わす。
【0117】以上のようにしてNOx変化分に対し、制
約条件である前記(5)式、(5a)式を満たすように
最小値と比較を行いつつ、前記(7a)式〜(7c)式
に基づき、各脱硝プロセスの修正NOx量∧(ハット)
oi(1)(i=1,2,3)を計算する(S13)。
この修正NOx量∧(ハット)oi(1)は修正された
NOx配分率に相当するものである。
【0118】しかる後、各脱硝プロセスの修正NOx量
∧(ハット)oi(1)を加算し、その総和がグループ
全体のNOx排出量設定値rを越えているか否かを判断
し(S14)、越えていない場合にはその修正NOx量
をそのまま採用する(S15)。ステップS14におい
て修正NOx量の総和がNOx排出量設定値rを越えて
いる場合には越えないような処理を実行する(S16)
(前記(8a)式〜(8c)式参照)。
【0119】以上のようにして修正NOx量の総和が制
約条件を満たすような処理を行った後、計算終了条件を
満足したか否かを判断し(S17)。ここで、計算終了
条件を満足していない場合、kを+1インクリメントし
た後(S18)、ステップS12に戻り、同様の処理を
繰り返し実行し、各脱硝プロセスに対するアンモニア使
用量が少なく、かつ、NOx排出量を低減化できる設定
値sv1,sv2,sv3を見つけ出し、設定値計算処
理手段23に送出する。
【0120】なお、設定値計算処理手段23は第1の実
施の形態で説明した設定値計算処理ステップS3と同様
であるので、ここではその説明を省略する。
【0121】従って、以上のような実施の形態によれ
ば、アンモニア注入量を少しずつ増やしながら、制約条
件を満たしつつ各脱硝プロセスのNOx減少変化分を計
算し、減少変化分が大きい脱硝プロセスほどNOx配分
率を大きくし、このNOx配分率に従って脱硝プロセス
を制御する制御装置のアンモニア注入設定値を決定する
ので、図2(a)及び図3(a)に示すような脱硝制御
を実現できる。すなわち、本発明装置は、個別ユニット
のNOx発生量及び各脱硝プロセスの特性を用いて最適
化演算を行うことにより、各脱硝プロセスのNOx排出
量設定値を算出でき、従来の制御装置に比べて、例えば
同じ性能を維持しつつアンモニア消費量を格段に少なく
でき、環境的および経済的にも優れた脱硝制御装置を実
現できる。
【0122】(第3の実施の形態)第1及び第2の実施
の形態では、予め各脱硝プロセスの特性f1(x)〜f
n(x)を設定することを前提としている。その理由
は、例えば発電プラントの排煙脱硝プロセスの特性は、
通常その運転開始前に確認することが要求されているの
で、その値を用いることが現実的とされているためであ
る。
【0123】しかし、脱硝プロセスの特性は、時間が経
過するに伴って変化することが経験的に知られている。
その結果、運転開始から時間が経つと、)既に得られて
いる特性と実際の特性との間に差異が生じる。このこと
は、実際の特性とは異なる過去の特性を用いて最適化処
理を行ったことになり、ひいては最適な運転とは言えな
くなり、本制御方法及び制御装置のメリットがなくな
る。
【0124】そこで、本実施の形態は、経時的な変化に
対応させるために、ユニット・脱硝プロセスを含むプラ
ントの運転データから脱硝プロセスの特性を取得する機
能を設けた構成例である。
【0125】以下、特性の取得手順を説明するが、当該
取得手順の説明はi軸の個別脱硝プロセスの適用例であ
り、グループ全体の場合には、以下の取得手順が第1軸
〜第n軸に対して個別に実施することになる。
【0126】この実施の形態における各脱硝プロセスの
特性は、第1の実施の形態と同様に、アンモニア注入量
とNOx発生量とのモル比をy、NOx排出量とNOx
発生量とのモル比をxとすると、yとxの関係は次式の
ような多項式で近似できる。
【0127】 y=a0+a1x+a2x2+…+amxm …(21) 先ず、運転データからy,xの組を収集する。今、全部
でn組の収集データ(x,y)を得たものとすると、任
意のk番目の収集データに関し前記(21)式の関係が
成立すると仮定したとき、n組の収集データの誤差を最
小にするパラメータ(a0,…,am)は以下の式から
計算できる。すなわち、
【0128】
【数4】
【0129】とおくと、前記(21)式が常に成立する
ならば、 Y=Фθ ……(23) となるので、 θ=(ФTФ)-1ФTY ……(24) となり、この(24)式から得られるθは、二乗誤差を
最小にする値となる。よって、(24)式の計算式から
第i軸の関数パラメータを取得できる。各軸の関数パラ
メータの取得に対しても同様の計算を用いて求めること
ができる。
【0130】従って、この実施の形態によれば、プラン
トの運転データであるアンモニア注入量とNOx発生量
とを収集し所定の演算を実行することにより、経時的な
変化に影響されない各脱硝プロセスの特性を容易に取得
でき、当該特性を用いることにより最適運転に必要な設
定値を取出すことができる。
【0131】(第4の実施の形態)以上述べた第3の実
施の形態はオフラインによる計算によって脱硝プロセス
の特性を求める例であるが、本実施の形態は、プラント
運転中にプランから運転データを取り込んでオンライン
による計算により各脱硝プロセスの特性を求める例であ
る。つまり、プラントの運転中に特性曲線のパラメータ
を求める制御方法及び制御装置である。
【0132】本実施の形態は、第1および第2の実施の形
態に関わる運転データからオンラインにて脱硝プロセス
の特性を取得する機能を設けた構成例である。
【0133】すなわち、この実施の形態は、第3の実施
の形態と同様に、アンモニア注入量とNOx発生量との
モル比をy、NOx排出量とNOx発生量とのモル比を
xとすると、yとxとの関係は前述したように(21)
式で近似できる。
【0134】今、パラメータa0,…,anに関し、前
回の計算結果a0(t−1),…,an(t−1)が与
えられているとする。計算開始時に適宜なa0(0),
…,an(0)を与えることにより、常に前回の計算結
果が利用できる。
【0135】
【数5】
【0136】とすると、前回の制御周期の計算結果θ
(t−1)+P(t)と今回の制御周期に得られたデー
タy(t),x(t)とを用いて、下記式のごとく新し
いパラメータθ(t)と補助行列P(t)が算出でき
る。
【0137】 P(t)=P(t−1)−P(t−1)ψ(t){I+ψ(t)T P(t−1)ψ(t)}-1ψT(t)P(t−1) ……(26a) θ(t)=θ(t−1)+P(t)ψ(t) {y(t)−ψ(t)Tθ(t−1)}……(26b) 従って、以上のような実施の形態によれば、オンライン
にて脱硝プロセスの特性曲線のパラメータを取得でき、
運転を停止することなくパラメータを取得しつつ最適化
処理に利用できる。
【0138】(第5の実施の形態)この実施の形態は、
第1ないし第4の実施の形態で得られる過去の脱硝プロセ
スの特性と第3、第4の実施の形態で取得されるプラン
トの運転データおよびこれら運転データから求める脱硝
プロセスの特性を表示部に表示し、その特性の違いを一
目で分かるようにした例である。
【0139】図6はかかる実施の形態に関する脱硝制御
装置の一例を示す構成図ある。
【0140】このグループ脱硝制御装置20は、前述し
たパラメータ・条件データ設定手段21、各ユニットN
Ox配分処理手段22、設定値計算処理手段23の他、
各脱硝プロセスの特性関数を求めるために各ユニット・
各脱硝プロセスを含むプラントから必要なデータ例えば
アンモニア注入量,NOx発生量,NOx排出量その他
特性を求めるために必要なプロセス量などの運転データ
を取り込む入力ポートまたはインターフェース31と、
データメモリ32と、インターフェース31を介して取
り込んだ運転データやパラメータ・条件データ設定手段
21で既に設定されている各脱硝プロセスの特性データ
等を取り込んでデータメモリ33に記憶するデータ記憶
手段33と、運転データから各脱硝プロセスの特性を演
算するCPUで構成された特性演算処理手段34と、表
示すべき画像データを記憶する画像メモリ35と、出力
制御手段36と、表示部37と、プリンタ38等によっ
て構成されている。
【0141】前記データメモリ32には、運転データ、
パラメータ・条件データ設定手段21で既に設定されて
いる各脱硝プロセスの特性データの他、各脱硝プロセス
の特性を求めるための近似式データ、演算処理上必要な
固定データ、演算処理途中及び演算処理後の特性データ
その他必要なデータを記憶するものである。なお、少な
くとも運転データは各脱硝プロセスごとに区分して記憶
するものとする。
【0142】前記特性演算処理手段34は、プログラム
データに基づいて所定周期ごとに運転データを取り込ん
でデータメモリ32に記憶するとともに、このデータメ
モリ32に記憶される運転データを用いて第3,第4の
実施の形態で説明したように各脱硝プロセスごとの特性
を演算し、外部からの表示指示に基づいて所要とする形
式で表示するためにデータメモリ32から各脱硝プロセ
スごとの特性等を取出して画像メモリ35に書込み、出
力制御手段36を介して表示部37に表示し、必要に応
じてプリンタ38に出力する等の一連の処理を実行す
る。
【0143】なお、外部からの表示指示は如何なる形式
(図7(a)〜(f))で表示するかによって異なる。
【0144】図7(a)はある同一の脱硝プロセスにお
ける第1〜第4の実施の形態に係る特性f1(x)と、当
該脱硝プロセスに係るプラントから取り込んだ運転デー
タ(イ)とこれら運転データから求めた今回の特性f1
(x)newとの表示比較例を示す図である。
【0145】なお、運転データから今回の特性f1
(x)newを求める例は、特性演算手段34においてデ
ータメモリ32に記憶される近似式或いは第3,第4の
実施の形態で求めた演算によって計算できる。
【0146】図7(b)はデータメモリ32から複数の
脱硝プロセスの特性f1(x)とf2(x)とを取出し
て比較表示した例を示す図である。
【0147】図7(c)はデータメモリ32に記憶され
る特性演算の近似式を読み出して表示し、入力指示のも
とにオペレータから入力されるNOx排出量とNOx発
生量とのモル比xまたはアンモニア注入量とNOx発生
量とのモル比yに基づき、近似式を用いて自動的にモル
比yまたはモル比xを演算し表示する例である。そし
て、オペレータから入力される複数のNOx排出量とN
Ox発生量とのモル比xまたはアンモニア注入量とNO
x発生量とのモル比yに基づいて求めたモル比yまたは
モル比xはデータメモリ32にテーブル化される。
【0148】図7(d)は同図(c)によって入力され
た複数のxまたはyに基づいて計算されたyまたはxの
一覧表を表示する例である。この表示の場合、表示部3
7の一画面に図7(c)と同図(d)とを上下または左
右に並べて表示してもよい。
【0149】図7(e)はある脱硝プロセスに対し、種
々のパラメータを考慮したときの特性関数の違いを同時
表示した表示例を示す図である。
【0150】図7(f)は同図(e)を発展させた表示
例であって、パラメータの他、例えば触媒温度T等を考
慮した演算された特性関数を三次元表示した表例を示す
図である。
【0151】(その他の実施の形態) (1) 以上述べた実施の形態における個別脱硝プロセ
スの特性及び第5の実施の形態における表示例は、専ら
排出されるNOxのモル数と発生したNOxのモル数と
の比、アンモニアのモル数と発生したNOxのモル数と
の比によって表現したものに限定したが、プラントから
得られるNOx量,アンモニア注入量以外の運転データ
であるプロセス量を用いて各脱硝プロセスの特性を求め
ることが出来る。ここで、言うプロセス量とは、例えば
ガスタービンの回転数・発電出力・燃料弁開度・起動時
間・排ガス温度・排ガス流量・排ガス圧力、排ガスを構
成する成分としての窒素ガス(N2)・酸素ガス(O2
・窒素酸化物ガス(NOx)・炭素酸化物(COx)・
硫黄酸化物(SOx)・水分(H2O)のそれぞれの濃
度及び流量、触媒温度、蒸気発生器が設置されている場
合には周囲の熱交換器の蒸気温度・蒸気圧力、蒸気ター
ビンの回転数・発電出力・加減弁開度・起動時間・給水
流量等を意味し、これらプロセス量のうち少なくとも1
つ以上を用いて、各脱硝プロセスの特性を算出でき、ま
た表示をすることも可能である。つまり、各脱硝プロセ
スの特性を求めるられるNOx量,アンモニア注入量以
外の他のプロセス量を用いて、各脱硝プロセスの特性を
演算し、最適化処理に利用することが可能である。これ
は、容易に類推可能であることから、特に上記実施の形
態では記述されていない。
【0152】(2) 上記実施の形態では、最適化のア
ルゴリズムとして、非線形計画法を使用したが、従来か
ら知られている例えばニューラルネットワーク、遺伝的
アルゴリズムなどを利用して最適化処理を行ってもよ
い。
【0153】(3) また、上記実施の形態では、最適
化処理に関し、説明を簡単にするために複数の仮定(b
11〜b18)を述べたが、これらの仮定を設けない場
合でも最適化処理が可能なアルゴリズムに容易に拡張で
きることは言うまでもない。また、各種の制約条件が存
在する場合でも、最適化処理を行うことが出来るように
アルゴリズムを拡張することは容易である。従って、こ
れらの拡張は容易に類推可能であることから、特に上記
実施の形態では記述されていない。さらに、グループ内
において特性の大きく異なるNOx発生源ないしは脱硝
プロセスが混在する場合であっても、最適化処理に有効
であることは言うまでもない。
【0154】(4) また、第1及び第2の実施の形態
では、複数の脱硝プロセスからなるプラントについて有
効な制御アルゴリズムであるのに対し、第3〜第5の実
施の形態では、グループ及び単独のNOx発生源および
単独の脱硝プロセスからなるプラントについて適用可能
な技術である。
【0155】(5) 脱硝プロセスの特性として、第1
〜第5の実施の形態では、関数の形で記述できるものと
したが、例えば図7(d)に示すようにテーブル形式で
特性が与えられている場合でも、同様に適用可能であ
る。なお、テーブル形式で特性データが与えられた場合
には、例えば第1,第2の実施の形態における最適化処
理の中の傾き計算を行う部分について、何らかの関数で
データ間を補間して傾きを計算するとか、微分可能な関
数で近似して傾き計算を行う手法を使用できる。
【0156】その他、本発明は、上記実施の形態に限定
されず、その要旨を逸脱しない範囲で種々変形して実施
可能であることは言うまでもない。
【0157】
【発明の効果】以上説明したように本発明によれば、各
脱硝プロセスの特性を有効に活用しつつ、最適化処理に
よって決定される各脱硝プロセスのNOx配分から定ま
る設定値を用いて制御することにより、グループ全体に
課される制約を満たしつつ、グループ全体でアンモニア
の消費量を大幅に低減化できる脱硝制御方法およびその
制御装置を提供できる。
【0158】また、プラントの運転データを用いて、オ
フラインまたはオンラインで容易に各脱硝プロセスの特
性を求めることができる脱硝制御装置を提供できる。
【0159】さらに、脱硝プロセスの特性その他の状態
を容易に確認でき、脱硝プロセスの健全性や特性変化を
把握でき、監視,制御及び経済性等の判断に広く利用で
き、非常に有益なものとなる。
【図面の簡単な説明】
【図1】 本発明に係る脱硝制御方法の一実施の形態を
説明する処理の流れ図。
【図2】 本発明方法と従来方法と脱硝制御の状態を表
わす図。
【図3】 同じく本発明方法と従来方法と脱硝制御の状
態を表わす図。
【図4】 本発明に係る脱硝制御装置の一実施の形態を
示す構成図。
【図5】 本発明装置の動作を説明するフローチャー
ト。
【図6】 本発明に係る脱硝制御装置の他の実施形態を
示す構成図。
【図7】 図6に示す装置を用いて脱硝プロセスの特性
を表示する一表示例図。
【図8】 従来の脱硝制御装置を説明する構成図。
【図9】 従来の別の脱硝制御装置を説明する構成図。
【符号の説明】
S1…パラメータ・条件データ設定ステップ S2…各脱硝プロセスNOx配分処理ステップ S3…設定値計算処理ステップ 10…グループ脱硝プロセス 15…軸制御装置 16…アンモニア噴霧器 17…脱硝触媒 18a…アンモニア注入弁 21…パラメータ・条件データ設定手段 22…各脱硝プロセスNOx配分処理手段 23…設定値計算処理手段 36…表示部
───────────────────────────────────────────────────── フロントページの続き Fターム(参考) 4D002 AA12 AC01 AC10 BA03 BA06 BA14 CA01 CA11 CA20 DA07 EA02 GA02 GA03 GB01 GB02 GB06 GB08 HA10 4D048 AA06 AB02 AC04 CC38 DA01 DA10 5H004 GB04 HA02 HA04 HB02 HB04 KC09 KC10 KC14

Claims (10)

    【特許請求の範囲】
  1. 【請求項1】 グループ脱硝プロセスを構成する複数の
    脱硝プロセスの制御用設定値を演算する脱硝制御方法に
    おいて、 前記各脱硝プロセスのアンモニア使用量の2乗を評価関
    数とし、かつ、前記各脱硝プロセスの排出NOx量の総
    和が予め定める前記グループ脱硝プロセス全体のNOx
    排出量規制値を越えないことを制約条件とし、 この制約条件を満たしつつ前記評価関数を最小化する最
    適化演算処理を行い、前記制御用設定値に相当する前記
    各脱硝プロセスのNOx配分率を決定することを特徴と
    する脱硝制御方法。
  2. 【請求項2】 グループ脱硝プロセスを構成する複数の
    脱硝プロセスの制御用設定値を演算する脱硝制御方法に
    おいて、 少なくとも各脱硝プロセスの特性データ、前記グループ
    脱硝プロセス全体のNOx排出量規制値データを設定す
    るパラメータ・条件設定ステップと、 この設定ステップによって設定される前記特性データか
    ら得られるアンモニア使用量の2乗を評価関数、前記グ
    ループ脱硝プロセス全体のNOx排出量規制値データを
    制約条件として用い、この制約条件を満たしつつ前記評
    価関数を最小化する最適化演算処理を実行し、前記各脱
    硝プロセスのNOx配分率を決定する各脱硝プロセスN
    Ox配分処理ステップと、 このNOx配分処理ステップによって得られるNOx配
    分率に基づいて前記各脱硝プロセスの制御用設定値を演
    算出力する設定値計算処理ステップとを有することを特
    徴とする脱硝制御方法。
  3. 【請求項3】 前記各脱硝プロセスNOx配分処理ステ
    ップは、前記特性データから得られる各脱硝プロセスの
    NOx発生量の総和が前記グループ脱硝プロセス全体の
    NOx排出量規制値データと等しくなるようなNOx発
    生量に関する計算開始条件値を設定する第1のステップ
    と、各脱硝プロセスごとにアンモニア使用量を増加させ
    たときのNOx減少変化分を求める第2のステップと、
    このNOx減少変化分の大きさに応じてNOx配分率を
    修正する第3のステップと、この第3のステップによる
    配分率に基づく各脱硝プロセスの修正NOx発生量の総
    和が前記グループ脱硝プロセス全体のNOx排出量規制
    値データより小さいとき、それらの修正NOx発生量に
    対応するNOx配分率を出力する第4のステップとを実
    行することを特徴とする請求項2に記載の脱硝制御方
    法。
  4. 【請求項4】 グループ脱硝プロセスを構成する複数の
    脱硝プロセスの制御用設定値を演算する脱硝制御装置に
    おいて、 前記各脱硝プロセスのアンモニア使用量の2乗を評価関
    数とし、また前記各脱硝プロセスの排出NOx量の総和
    が予め定める前記グループ脱硝プロセス全体のNOx排
    出量規制値を越えないことを制約条件とする手段と、 この制約条件を満たしつつ前記評価関数を最小化する最
    適化演算処理を行い、前記各脱硝プロセスのNOx配分
    率を決定する最適化問題解法手段とを備え、 この各脱硝プロセスのNOx配分率に対応する前記制御
    用設定値を得るようにしたことを特徴とする脱硝制御装
    置。
  5. 【請求項5】 グループ脱硝プロセスを構成する複数の
    脱硝プロセスの制御用設定値を演算する脱硝制御装置に
    おいて、 少なくとも各脱硝プロセスの特性データ、前記グループ
    脱硝プロセス全体のNOx排出量規制値データを設定す
    るパラメータ・条件設定手段と、 この設定手段によって設定される前記特性データから得
    られるアンモニア使用量の2乗を評価関数、前記グルー
    プ脱硝プロセス全体のNOx排出量規制値データを制約
    条件として用い、この制約条件を満たしつつ前記評価関
    数を最小化する最適化演算処理を実行し、前記各脱硝プ
    ロセスのNOx配分率を決定する各脱硝プロセスNOx
    配分処理手段と、 このNOx配分処理手段によって得られるNOx配分率
    に基づいて前記各脱硝プロセスの制御用設定値を演算出
    力する設定値計算処理手段とを備えたことを特徴とする
    脱硝制御装置。
  6. 【請求項6】 前記各脱硝プロセスNOx配分処理手段
    は、前記特性データから得られる各脱硝プロセスのNO
    x発生量の総和が前記グループ脱硝プロセス全体のNO
    x排出量規制値データと等しくなるようなNOx発生量
    に関する計算開始条件値を設定する設定手段と、各脱硝
    プロセスごとにアンモニア使用量を増加させたときのN
    Ox減少変化分を求めるNOx減少変化分演算手段と、
    このNOx減少変化分の大きさに応じてNOx配分率を
    修正する配分率修正手段と、この配分率に基づく各脱硝
    プロセスの修正NOx発生量の総和が前記グループ脱硝
    プロセス全体のNOx排出量規制値データより小さいか
    否かを判断し、小さいときそれらの修正NOx発生量に
    対応するNOx配分率を出力する判断手段とを設けたこ
    とを特徴とする請求項5に記載の脱硝制御装置。
  7. 【請求項7】 請求項5に記載の脱硝制御装置におい
    て、 前記各脱硝プロセスを含むプラントの運転データ、運転
    開始前または運転中の運転データから求められる前記各
    脱硝プロセスの特性データ等を記憶するデータ記憶手段
    と、このデータ記憶手段に記憶される特性データの中か
    ら同一の脱硝プロセスの過去特性データ及び前記運転デ
    ータから得られる今回特性データとを読み出して比較表
    示し、或いは複数の脱硝プロセスの特性データを読み出
    して比較表示し、或いは種々の異なる運転データである
    パラメータを考慮して得られた脱硝プロセスの特性デー
    タを読み出して二次元または三次元にて比較表示する表
    示制御手段とを設けたことを特徴とする脱硝制御装置。
  8. 【請求項8】 請求項5に記載の脱硝制御装置におい
    て、 前記各脱硝プロセスを含むプラントの運転データ、運転
    開始前または運転中の運転データから求められる前記各
    脱硝プロセスの特性データ、前記運転データの少なくと
    もアンモニア注入量とNOx発生量の第1のモル比とN
    Ox排出量と前記NOx発生量との第2のモル比との関
    係の近似式データ等を記憶するデータ記憶手段と、入力
    指示のもとに入力される前記第1または第2のモル比に
    基づき、前記データ記憶手段に記憶される前記近似式に
    より前記第2または第1のモル比を計算する手段と、こ
    の計算結果を順次テーブル化し表示する表示制御手段と
    を備えたことを特徴とする脱硝制御装置。
  9. 【請求項9】 前記各脱硝プロセスの特性データとして
    は、プラント運転前に予め確認する各脱硝プロセスの特
    性データまたは前記各脱硝プロセスを含むプラントの運
    転データであるオフライン時またはオンライン時のアン
    モニア注入量とNOx発生量のモル比とNOx排出量と
    前記NOx発生量とのモル比との多項式の関係から求め
    る各脱硝プロセスの特性データを用いることを特徴とす
    る請求項5、請求項7及び請求項8の何れか1つに記載
    の脱硝制御装置。
  10. 【請求項10】 前記各脱硝プロセスの特性データは、
    プラントから得られる運転データのうち、アンモニア注
    入量,NOx発生量およびNOx排出量以外のプロセス
    量から求める各脱硝プロセスの特性データを用いること
    を特徴とする請求項5、請求項7及び請求項8の何れか
    1つに記載の脱硝制御装置。
JP24012399A 1999-08-26 1999-08-26 脱硝制御方法およびその制御装置 Expired - Fee Related JP4057750B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP24012399A JP4057750B2 (ja) 1999-08-26 1999-08-26 脱硝制御方法およびその制御装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP24012399A JP4057750B2 (ja) 1999-08-26 1999-08-26 脱硝制御方法およびその制御装置

Publications (2)

Publication Number Publication Date
JP2001062246A true JP2001062246A (ja) 2001-03-13
JP4057750B2 JP4057750B2 (ja) 2008-03-05

Family

ID=17054848

Family Applications (1)

Application Number Title Priority Date Filing Date
JP24012399A Expired - Fee Related JP4057750B2 (ja) 1999-08-26 1999-08-26 脱硝制御方法およびその制御装置

Country Status (1)

Country Link
JP (1) JP4057750B2 (ja)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014211102A (ja) * 2013-04-18 2014-11-13 株式会社日立製作所 コンバインドサイクル発電プラント
KR101571994B1 (ko) 2012-04-02 2015-11-25 미쓰비시덴키 가부시키가이샤 파라미터 설정 장치
CN109583077A (zh) * 2018-11-27 2019-04-05 南京杰科丰环保技术装备研究院有限公司 一种集成多场状态信息的环保装置三维协同建模方法
CN115245735A (zh) * 2022-06-16 2022-10-28 浙江省生态环境科学设计研究院 基于预测控制模型的水泥窑烟气iSNCR控制方法
WO2023123577A1 (zh) * 2021-12-30 2023-07-06 深圳华明环保科技有限公司 一种干法脱硝工艺
CN116832614A (zh) * 2023-05-24 2023-10-03 华能国际电力股份有限公司上海石洞口第二电厂 一种scr脱硝系统喷氨量控制方法及系统

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101571994B1 (ko) 2012-04-02 2015-11-25 미쓰비시덴키 가부시키가이샤 파라미터 설정 장치
JP2014211102A (ja) * 2013-04-18 2014-11-13 株式会社日立製作所 コンバインドサイクル発電プラント
CN109583077A (zh) * 2018-11-27 2019-04-05 南京杰科丰环保技术装备研究院有限公司 一种集成多场状态信息的环保装置三维协同建模方法
CN109583077B (zh) * 2018-11-27 2022-12-27 南京杰科丰环保技术装备研究院有限公司 一种集成多场状态信息的环保装置三维协同建模方法
WO2023123577A1 (zh) * 2021-12-30 2023-07-06 深圳华明环保科技有限公司 一种干法脱硝工艺
CN115245735A (zh) * 2022-06-16 2022-10-28 浙江省生态环境科学设计研究院 基于预测控制模型的水泥窑烟气iSNCR控制方法
CN116832614A (zh) * 2023-05-24 2023-10-03 华能国际电力股份有限公司上海石洞口第二电厂 一种scr脱硝系统喷氨量控制方法及系统
CN116832614B (zh) * 2023-05-24 2024-05-24 华能国际电力股份有限公司上海石洞口第二电厂 一种scr脱硝系统喷氨量控制方法及系统

Also Published As

Publication number Publication date
JP4057750B2 (ja) 2008-03-05

Similar Documents

Publication Publication Date Title
US7219040B2 (en) Method and system for model based control of heavy duty gas turbine
US8185216B2 (en) Plant controlling device and method, thermal power plant, and its control method
JPH06205935A (ja) 脱硝制御装置
KR20130089283A (ko) 산소 연료 연소 파워 플랜트를 위한 최적화된 통합 제어
JP6173840B2 (ja) 脱硝制御装置
CN111841276B (zh) 循环流化床机组sncr脱硝控制方法及装置、存储介质
CN111365733B (zh) 一种燃煤锅炉经济氧量动态寻优方法、系统及设备
Rafiei et al. Integration of design and control for industrial-scale applications under uncertainty: A trust region approach
JP2001062246A (ja) 脱硝制御方法およびその制御装置
US20060121616A1 (en) Method and system for increasing efficiency of FGD operation in fossil fuel boilers
Henao et al. Surrogate-based process synthesis
JP3299531B2 (ja) 発電プラント
JP4792696B2 (ja) 脱硝制御方法、脱硝制御装置及びそのプログラム
JP5117232B2 (ja) ボイラを備えたプラントの制御装置、及びボイラを備えたプラントの制御方法
JP7075194B2 (ja) 廃棄物処理計画演算方法および廃棄物処理計画演算装置
CN113578008A (zh) 一种垃圾焚烧锅炉脱硝控制方法及系统
JP2000163102A (ja) 非線形最適状態フィードバック制御方法及び装置
CN114089636A (zh) Scr脱硝外挂式智能喷氨闭环控制方法与设备
CN112926765A (zh) 一种脱硫系统运行优化方法以及信息物理融合系统
CN206848713U (zh) 一种模糊规则在线可调的自适应模糊控制器
JPH08168639A (ja) 脱硝触媒を内蔵した脱硝装置へのアンモニア注入量制御方法および制御装置
Li et al. A data‐based optimal setting method for the coking flue gas denitration process
JPH06257425A (ja) 系列脱硝制御装置
CN109395549A (zh) 一种喷氨控制方法、装置及计算机可读存储介质
JP3492744B2 (ja) 系列脱硝制御装置

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20040908

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20060329

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20070109

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20070307

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20070828

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20070928

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20071211

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20071214

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20101221

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20111221

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20121221

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20121221

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20131221

Year of fee payment: 6

LAPS Cancellation because of no payment of annual fees