JP2003216205A - プラント運用制御装置 - Google Patents

プラント運用制御装置

Info

Publication number
JP2003216205A
JP2003216205A JP2002010377A JP2002010377A JP2003216205A JP 2003216205 A JP2003216205 A JP 2003216205A JP 2002010377 A JP2002010377 A JP 2002010377A JP 2002010377 A JP2002010377 A JP 2002010377A JP 2003216205 A JP2003216205 A JP 2003216205A
Authority
JP
Japan
Prior art keywords
plant
unit
operation plan
plan
control device
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
JP2002010377A
Other languages
English (en)
Other versions
JP4099334B2 (ja
Inventor
Yoshiyuki Sakamoto
本 義 行 坂
Futoshi Kurokawa
川 太 黒
Kenichi Yamazaki
崎 謙 一 山
Tatsuo Ashiki
木 達 雄 芦
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 JP2002010377A priority Critical patent/JP4099334B2/ja
Publication of JP2003216205A publication Critical patent/JP2003216205A/ja
Application granted granted Critical
Publication of JP4099334B2 publication Critical patent/JP4099334B2/ja
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Landscapes

  • Complex Calculations (AREA)
  • Feedback Control In General (AREA)

Abstract

(57)【要約】 【課題】 制御対象プラントが複雑もしくは大規模なも
のであったとしても、最適な運用計画についての演算を
高速且つ精度良く行うことができ、効率的で安定したプ
ラント運用を行えるようにすること。 【解決手段】 プロセスデータ記憶部4は、プロセスコ
ントローラ2を介して得られるプロセスデータの計測値
や種々のパラメータ設定値等のデータを保存する。需要
予測部5は、プロセスデータ記憶部4から得られる過去
の実績需要地、あるいは気象情報その他のデータに基づ
いて、運転該当日の単位時間あたりの需要量を1日分予
測する。運用計画部6は、需要予測部5が予測した単位
時間あたりの1日分需要量とプロセスの計測値に基づい
て該当日の単位時間あたりの機器起動停止計画を整数値
遺伝子でコード化した遺伝的アルゴリズムで演算する。

Description

【発明の詳細な説明】
【0001】
【発明の属する技術分野】本発明は、上下水道プラント
をはじめとし、鉄鋼プラントなどの産業プラント、地域
冷暖房プラントや各ビル毎の熱源プラント、あるいは発
電所プラントなどの各種プラントにおいて用いられるプ
ラント運用制御装置に関するものである。
【0002】
【従来の技術】プラントは、一般に制御対象となる機器
の数が多く、また、種々の制御条件を考慮しなければな
らないために制御パラメータ数が非常に多くなる。した
がって、プラントの運用計画を理論的、数学的手法に基
づき立案するのは困難である。更に、上下水道プラント
などいくつかのプラントでは、天候や季節等によっても
少なからず需要量が左右されることがあり、理論的、数
学的手法に基づく立案は一層困難なものとなっている。
そのため、熟練したオペレータの経験側を考慮しつつ、
プラントの基本的運用を行っている例が多くなってい
る。
【0003】
【発明が解決しようとする課題】しかし、近年はコンピ
ュータ技術の発展に伴って、所謂「動的計画法」や「分
岐限定法」といった数学的最適化手法が著しく発展して
きており、ある程度数学的に立案することが可能になっ
てきている。ところが、このような数学的最適化手法で
は、プラント規模や機器数の増大あるいは制御内容の複
雑化に伴い、演算時間が指数関数的に急激に増大し、実
用性が失われることになる。例えば、ある日の夜に、明
日の朝以降のプラント運用計画を立案しようとしている
場合に、演算に要する時間が十数時間以上では、得られ
る結果が最適なものであったとしても何ら意味をなさな
くなる。
【0004】本発明はこのような事情に鑑みてなされた
ものであり、制御対象プラントが複雑もしくは大規模な
ものであったとしても、最適な運用計画についての演算
を高速且つ精度良く行うことができ、効率的で安定した
プラント運用を実現することができるプラント運用制御
装置を提供することを目的としている。
【0005】
【課題を解決するための手段】上記課題を解決するため
の手段として、請求項1記載の発明は、プラント内のプ
ラント機器に対する制御を予め作成してある運用計画に
基づき行うプラント運用制御装置において、前記プラン
トからプロセスコントローラを介して入力したプロセス
データを記憶するプロセスデータ記憶部と、前記プロセ
スデータ記憶部に記憶された過去の実績値、及び所定の
情報に基づき、需要を予測する需要予測部と、前記需要
予測部で予測された需要に基づき、遺伝的アルゴリズム
を用いて最適な運用計画についての演算を行う運用計画
演算手段を有する運用計画部と、を備えたことを特徴と
する。
【0006】請求項2記載の発明は、請求項1記載の発
明において、前記遺伝的アルゴリズムは、コード化され
た整数値遺伝子を用いるものである、ことを特徴とす
る。
【0007】請求項3記載の発明は、請求項1又は2記
載の発明において、前記運用計画部は、前記遺伝的アル
ゴリズム以外の最適化手法により得られた暫定解を、前
記遺伝的アルゴリズムの初期値として演算し、これを前
記運用計画演算部に入力させる初期値演算手段を有す
る、ことを特徴とする。
【0008】請求項4記載の発明は、請求項1又は2記
載の発明において、前記運用計画部は、過去のプラント
運用における制御で実際に用いた値を前記プロセスデー
タ記憶部から取り出し、これを前記遺伝的アルゴリズム
の初期値として、前記運用計画演算手段に入力させるプ
ラント実績初期値入力手段を有する、ことを特徴とす
る。
【0009】請求項5記載の発明は、請求項1乃至4の
いずれかに記載の発明において、前記運用計画部は、前
記運用計画演算手段が前記遺伝的アルゴリズムを用いて
演算した運用計画について、更に、遺伝的アルゴリズム
以外の最適化手法に基づき最適化演算を行うハイブリッ
ド最適化手段を有する、ことを特徴とする。
【0010】請求項6記載の発明は、請求項1乃至5の
いずれかに記載の発明において、前記運用計画部は、前
記運用計画演算手段が演算を行う際の制約条件を緩和す
る制約条件緩和手段を有する、ことを特徴とする。
【0011】請求項7記載の発明は、請求項1乃至6の
いずれかに記載の発明において、前記運用計画部は、前
記運用計画演算手段又は前記ハイブリッド最適化手段が
演算した運用計画に係る経時的な特性曲線中に凹部又は
凸部が存在する場合、最適化のための評価値を低下させ
ることなく、これらの凹部又は凸部を平滑化する凹凸平
滑化手段を有する、ことを特徴とする。
【0012】請求項8記載の発明は、請求項1乃至7の
いずれかに記載の発明において、前記運用計画部は、前
記運用計画演算手段又は前記ハイブリッド最適化手段が
演算した運用計画に係る経時的な特性曲線中に上り階段
状部分又は下り階段状部分が存在する場合、最適化のた
めの評価値を低下させることなく、これらの上り階段状
部分又は下り階段状部分を平滑化する上下階段平滑化手
段を有する、ことを特徴とする。
【0013】請求項9記載の発明は、請求項1乃至8の
いずれかに記載の発明において、プラント内にローカル
制御コントローラのみにより制御されるローカル機器が
含まれている場合に、そのローカル制御の模擬を行うロ
ーカル制御模擬計画部を備えており、前記運用計画部の
運用計画演算手段は、このローカル制御模擬計画部の模
擬結果を考慮して前記運用計画についての演算を行うも
のである、ことを特徴とする。
【0014】請求項10記載の発明は、請求項1乃至9
のいずれかに記載の発明において、前記需要予測部は、
プラント機器同士間の離間距離が大きく、あるプラント
機器の運転開始後その運転結果が他のプラント機器に反
映されるまでに一定レベル以上の時間遅れを生じる場合
に、その時間遅れを考慮した需要予測を行う時間遅れ修
正手段を有するものである、ことを特徴とする。
【0015】
【発明の実施の形態】以下、本発明の実施形態を制御対
象プラントが上水道プラントである場合を例に取り説明
する。図2は、プラント1の概略構成図である。この図
において、複数台の取水ポンプ51は河川等から原水を
取水し、これを着水井52に送り出す。着水井52に貯
溜された原水は、流量計を介して混和池53に流入す
る。混和池53に流入した原水は、薬品の注入と共に攪
拌機による攪拌が行われた後フロック形成池54に流入
する。そして、ここで徐々に沈殿するフロックが形成さ
れた後、沈殿池55において汚物が沈殿し、上澄み水が
ろ過池56に送られる。ろ過池56では、このろ過水が
ろ過された後、塩素の注入により殺菌処理が施され、浄
水池57に流入する。浄水池57に貯溜された処理水
は、複数台の送水ポンプ58により配水池59に送水さ
れ、更に、配水池59に貯溜された処理水は複数台の配
水ポンプ60により配水ポンプ61に対する配水が行わ
れる。
【0016】図1は、上記のプラント1に対する運用を
行うプラント運用制御装置の構成を示すブロック図であ
る。プラント運用制御装置3は、プロセスコントローラ
2を介してプラント1との間でデータの授受を行うよう
になっており、プロセスデータ記憶部4、需要予測部
5、運用計画部6、及びローカル制御模擬計画手段7を
有している。そして、需要予測部5は時間遅れ修正手段
8を有しており、また、運用計画部6は、初期値入力手
段9、初期値演算手段10、制約条件緩和手段11、運
用計画演算手段12、ハイブリッド最適化手段13、凹
凸平滑化手段14、及び上下平滑化手段15を有してい
る。
【0017】上水道プラントの自動制御において、計算
機が作成した運用計画に基づいてプラントの運用及び制
御を行う場合、浄水を断水させることなく、需要家に対
して安定した供給を行うことが絶対に必要とされる条件
となる。それと同時に、プラント運用上、運用コストや
設備保守費等を極力低減するために、効率的なプラント
の運用及び制御を行うことが要求される。したがって、
このような観点からの要望を満たすように、最適又は近
似的に最適に近いプラント運用計画の作成を可能にする
ことが重要となるが、近年のプラントの複雑化及び大規
模化に伴い、上記のような最適又は近似的に最適に近い
プラント運用計画を迅速に演算することが益々困難な状
況となってきている。
【0018】図1に示したプラント運用制御装置3は、
このような課題を解決する機能を有するものである。プ
ロセスデータ記憶部4は、プロセスコントローラ2を介
して得られるプロセスデータの計測値や種々のパラメー
タ設定値等のデータを保存するものである。需要予測部
5は、プロセスデータ記憶部4から得られる過去の実績
需要地、あるいは気象情報その他のデータに基づいて、
運転該当日の単位時間あたりの需要量を1日分予測する
ものである。運用計画部6は、需要予測部5が予測した
単位時間あたりの1日分需要量とプロセスの計測値に基
づいて該当日の単位時間あたりの機器起動停止計画を整
数値遺伝子でコード化した遺伝的アルゴリズムで演算す
るものである。
【0019】需要予測部5が行う需要予測の方法として
は、例えば、統計的な手法や最小自乗法、GMDH(Gr
ouping Method of Data Handling)などの各種同定手
法、ニューラルネットワークを用いた手法などが考えら
れるが、特に手法が限定されるわけではなく、いずれの
手法を用いてもよい。
【0020】需要予測と運用計画は1日1回以上、定刻
に起動する。まず、定刻前までに、需要予測に必要なデ
ータを手動又は自動で入力する。例えば、需要予測した
い日の天気予報や最高又は最低気温予報といった天候情
報や、それまでに得られた天候情報の実績値や需要量実
績値などである。需要予測部5による需要予測の結果は
ある単位時間毎に少なくとも1日の区切りの時刻分まで
を出力する。
【0021】需要予測部5による需要側結果とプラント
の現在の計測値(浄水の送水計画であれば、配水池59
の水位や送水ポンプ58の流量、起動台数、配水量等)
と送水ポンプ流量特性、配水池59の容量(運用水位上
下限値)などのパラメータに基づいて、配水池59の運
用水位上下限の逸脱や排水量予測値に不足なく且つ急激
な送水量の変化がない送水計画を、運用計画演算手段1
2が演算する。この際、運用コストをできるだけ低減化
するような考慮が要求されることもある。この場合に
は、必要であれば、電力消費のピーク時間帯にポンプの
起動可能な台数を制約することも考えられる。こうし
て、整数値遺伝子でコード化した遺伝的アルゴリズムを
持つ運用計画演算手段12は、プラント機器の最適な運
用計画を単位時間毎に、少なくとも1日の区切りの時刻
分までを出力する。浄水の送水計画であれば、送水ポン
プ58の起動停止計画と、単位時間毎の総水量、配水池
59の単位時間毎の水位となる。
【0022】いま、図2に示した上水道プラントにおい
て、浄水池57から送水ポンプ58によって、配水池5
9を軽油して配水区61へ浄水を配水するプロセスを想
定して説明することとする。もちろん、自然流下で浄水
池57から配水池59へ送水したり、配水池59から配
水区61へ配水することもあるが、ここでは、自然流下
ではなく、送水ポンプ58で送水し、配水ポンプ60で
配水することとしても一般性は失われない。
【0023】或る時刻kで浄水池57から配水池59へ
の総水量Qp(k)は、起動する配水ポンプ60や配管に
取り付けられたバルブの開度によって制御される。その
目標となる流量を離散的に何段階かに決めることとす
る。例えば、バルブであれ、ポンプであれ、とにかく送
水量の物理的最大値が100〔m3/h〕だとするなら、
仮にこれを6段階に分けるとすると、0,20,40,
60,80,100〔m3/h〕が送水計画として取り得
る送水量となる。送水ポンプ58は、回転数固定の固定
速ポンプn台のみだとすると、固定速ポンプの台数と送
水計画としてとり得る送水量の離散値は1対1で対応す
るが、説明を簡単にするため以後この場合を想定するこ
とにする。
【0024】図3は、ある1日における1時間毎のポン
プ運転台数を示した図表であるが、遺伝的アルゴリズム
ではこのポンプ台数の欄に記載された整数値を遺伝子と
呼び、この図表にに相当するものを整数値遺伝子(個
体)と呼ぶ。このように、整数値遺伝子でコード化した
遺伝的アルゴリズムでは、遺伝子が整数値で表現され
る。この図3の内容は、時刻0時台(0:00〜1:00)では
ポンプを1台起動する計画を表しており、1時台(1:00
〜2:00)では2台、2時台(2:00〜3:00)では3台を起
動する計画を表している。以降の時間帯も同様である。
【0025】上記のように、ポンプの台数(自然流下の
場合はバルブの開度によって制御される総水量を表し、
また、可変速ポンプで送水の場合は送水ポンプの送水流
量を離散的に決めた何段階かのステップ数を表す)を表
す整数値を遺伝的アルゴリズムでは遺伝子と呼んでいる
が、固定速ポンプのみ5台からなる送水プロセスのとき
の遺伝子は、0,1,2,…,5のいずれかを取ること
になる。
【0026】いま、考えている送水計画の最適化問題を
次のように定式化することとする。定式化の方法は、ど
のような送水計画を最適化したいかによって変わるので
一意ではないが、どのように定式化したとしても、以下
に示すような組み合わせ最適化問題と呼ばれる定式化で
あれば、整数値遺伝子でコード化された遺伝的アルゴリ
ズムによって近似的に最適化することができる。以下に
示す最適化問題は、目的関数fを最小化し得るx,yを
求めようとするものである。目的関数fは(1)式によ
り表すことができ、(1)式におけるf1は(2)式で
表され、また、f2は(3)式又は(4)式で表され
る。また、このときの制約条件は(5),(6)式によ
り表される。
【数1】
【0027】上記の説明では、例えば、送水プロセスの
送水ポンプの数を5台と想定しているので、k=5とな
る。このような問題は、一般に組み合わせ最適化問題と
呼ばれ、いわゆる分岐限定法(BBM:Branch and Bounded
Method)によって最適なx(i)を求めることができ
る。図4は、この分岐限定法についての説明図である。
分岐限定法は、上述した最適送水計画問題の変数x
(i)を決定するために、以下のような動作をする。ま
ず、計画開始時(時刻0)の実際の離散的に決められた
送水量ステップ(固定速ポンプのみならポンプ運転台
数)が初期値となり、次の時刻(時刻1,2,3,…)
に進む。そのときのポンプ運転台数で制約条件を満足で
きるか否かをチェックしながら時刻毎のポンプ運転台数
を次々に分岐していく。計画最終時刻まで到達できたと
きに、目的関数値を評価する。目的関数値がこれまでの
解より良好であれば、その解を暫定解として保存し、さ
らに分岐探索を行うというものである。
【0028】このような分岐限定法は列挙法の一種であ
り、整数計画問題のような組み合わせ問題を解く場合、
最悪では全組み合わせ通りの探索を行ってしまうが、演
算時間に制限がなく、可解であれば必ず最適解が得られ
る。しかし、最適送水計画問題のような組み合わせ問題
は所謂NP-困難な問題の1つとして知られており、変
数の数に対して指数関数的に増大する演算回数を要す
る。そのため、送水ポンプを多数含むような大規模プラ
ントを対象とする場合、オンライン制御に有効な時間内
に最適解が得られない虞があり得る。
【0029】そこで、本実施形態で示すような、整数値
遺伝子でコード化した遺伝的アルゴリズムによる高速解
法が有効である。つまり、図3で示したように、整数値
遺伝子でコード化したら、図5のフローチャートに示す
ような遺伝的アルゴリズムを適用することによって、最
適又は近似的に最適な運用計画を得ることができる。
【0030】すなわち、まず、ステップ1では、ランダ
ムに遺伝子を割り当てて生成した固体をそれぞれ予め定
義した個数n個だけ発生させる(初期個体群の生成)。
ここで、個体とは、図3に示したような1つの遺伝子列
のことである。
【0031】ステップ2では、制約条件を満足していな
い場合には満足するように遺伝子を変更し、制約条件を
満足する固体がn個生成できたら、各固体の適応度f′
及びその世代での平均適応度を計算する(各個体の評
価)。なお、一般に、遺伝的アルゴリズムは最大問題に
適用されるものであるが、本実施形態では最小の目的関
数fを求めることになる。したがって、本実施形態にお
ける適応度f′とは、(1)式における目的関数fの逆
数に相当するものとなる。
【0032】ステップ3では、今回の世代の平均適応度
と前回以前数世代の平均適応度とを比較し、両者の差が
予め定めてある設定値ε以下であるか否かを判別する。
そして、判別結果が「YES」であれば処理を終了し、
「NO」であればステップ4に進む。
【0033】ステップ4では、制約条件を満足しない固
体や予め定義した固体数以上の固体が存在する場合に、
適応度の悪い(適応度の小さい)固体や、余分な個体を
その定義した個数になるまで淘汰(削除)する(淘汰処
理)。
【0034】ここで、ステップ2〜4につき具体的に説
明すると、ステップ2で述べた個体が(5)式の制約条
件を満足していない場合、満足していない時間の遺伝子
を次のようにして変更する。まず、(5)式の制約条件
における上限値を満足していない場合、配水池59にお
いては配水池59から出ていく量よりも入ってくる量
(送水量)が多く、水位が上昇するので上限値を下回る
ようになるまで送水量を減らせばよい。例えば、当初、
或る時刻の送水量のステップ数(ポンプ台数)が5であ
ったとすると、5→4→3… と、上限値を超えなくな
るまで送水量を減らしていく。そして、もし送水量のス
テップ数を0にしても上限値を超えている場合は、時間
を1時間分だけ遡って同様に送水量を減らしていく。結
局、計画のスタート時刻に相当する時間まで遡っても上
限値を超えてしまう場合には、その個体を淘汰対象と
し、次の個体の制約条件をチェックすることとする。
(5)式の制約条件における下限値を満足していない場
合も、上限値の場合と同様にして処理することができ
る。但し、下限値を上回るように、送水量を増やして行
く点が上限値の場合と異なる。
【0035】ステップ5では、ステップ4の淘汰処理の
結果、個体数が予め定義した数より少ない場合に、残っ
た個体の中から適応度が最良の個体を増殖(コピー)し
て、個体数が定義した数になるようにする(増殖処
理)。
【0036】ステップ6では、ランダムにペアリングを
行う(交叉処理)。ペアリングは全個体数に対する割合
(交叉率)分だけ行い、ペア毎にランダムに遺伝子座
(遺伝子の場所)を選び、一点交叉(選んだ遺伝子の場
所から交互に遺伝子のセットを交換)させる。例えば、
図6に示すように、2つの個体(親)が選ばれたとき、
3時と4時との間で遺伝子を交換させることとすると、
同図右に示すように、新たに2つの遺伝子列(個体)が
子として生成される。これを全体の個体数に対して、決
められた割合(交叉率)にしたがって操作する。このよ
うな処理が交叉処理である。
【0037】ステップ7では、全個体数に割合(突然変
異率)分だけランダムに個体を選び、各個体の任意の遺
伝子座(ランダムに決定する)を変更させる(突然変異
処理)。例えば、図7に示すように、ある個体(親)が
選ばれたとき、23時の遺伝子が同図左(遺伝子0)か
ら同図右(遺伝子1)のように、突然変更されることが
あるが、これを突然変異と呼ぶ。この突然変異は、全体
の遺伝子のうち、決められた割合(突然変位率)にした
がって操作される。
【0038】ステップ8では、世代数が予め定めてある
規定値以上に達したか否かを判別し、規定値以上に達し
ていれば処理を終了し、達していなければステップ2に
戻って同様の処理を繰り返す。
【0039】運用計画部6の運用計画演算手段12は、
需要予測部5の需要予測に基づいて、上記したステップ
1〜8の処理を行う。ここで、運用計画部6の最も簡単
な構成は、その構成要素が運用計画演算手段12のみの
場合であるが、本実施形態では手段9〜11、及び手段
13〜15が付加されている。次に、これら各手段の機
能につき順次説明する。
【0040】初期値演算手段9は、運用計画演算手段1
2が上記のステップ1で生成する初期個体すなわち初期
値を演算するものである。遺伝的アルゴリズムでは、初
期値として乱数などの任意の値を用いても最終的には最
適解もしくは最適解に近いものを求めることができる
が、ステップ2の制約条件のチェックに多くの処理時間
を要することとなり、結果的には演算実行時間のリミッ
トを超えてしまい、何らの運用計画を得られない場合も
あり得る。そこで、初期値演算手段9により、遺伝的ア
ルゴリズム以外の最適化手法(例えば、前述した分岐限
定法)を用い、最適運用計画ではないけれども制約条件
を一応全て満足する運用計画の候補を求め、この候補を
初期値として運用計画演算手段12に入力させるように
している。
【0041】初期値入力手段10は、初期個体をランダ
ム生成するだけでなく、過去の実績値から遺伝子列を予
め定義した個数だけ生成する。つまり、ランダムに生成
した個体やその他の方法(例えば分岐限定法)で生成し
た個体と合わせてn個となるように過去の実績値から遺
伝子列(個体)を生成するものである。この初期値入力
手段10を用いた構成では、過去に実際に用いた運用計
画を個体に含んでいるので、初期値演算手段9のみで演
算した運用計画を用いる場合よりも、より実用的なもの
であると考えられる。なお、本実施形態では、初期値演
算手段9又は初期値入力手段10のいずれか一が選択的
に用いられるようになっている。
【0042】制約条件緩和手段11は、運用計画演算手
段12の演算が所定時間内に終了しない場合や、所定プ
ラント状態の上下限範囲を満たさない場合に、この所定
の上下限範囲を制約条件から除外し、代わりに物理的な
上下限範囲を新たな制約条件に設定することで制約条件
を緩和し、所定時間内に最適運用計画を得られるように
するためのものである。但し、この場合、運用計画の最
適化演算における目的関数に、プラント状態のもとの上
下限範囲からの逸脱量に応じたペナルティを課すように
して、運用計画を演算する。
【0043】この制約条件緩和の具体例につき説明する
と、最適化問題を解く演算が予め定めた時間、例えば1
分などで完了しなかった場合、又は予め定めたプラント
の状態の上下限範囲(例えば、配水池の運用水位上下限
値)を満たさない場合、予め定めたプラントの状態の上
下限範囲を緩和し、運用計画の最適化演算における目的
関数にプラント状態の上下限範囲からの逸脱量にペナル
ティを課すよう、以下のように最適運用計画を再演算す
る。
【数2】
【0044】ハイブリッド最適化手段13は、運用計画
演算手段12で得られた近似的最適運用計画を入力し、
これを初期値にして遺伝的アルゴリズム以外の最適化手
法(例えば、分岐限定法や動的計画法)により最適運用
計画を得ようとするものである。遺伝的アルゴリズムの
みでは得られた運用計画が真に最適なものであるか否か
が分からず、また、分岐限定法などの最適化手法では演
算に多くの時間を要することになるが、このように遺伝
的アルゴリズムと他の手法とを組み合わせることによ
り、短時間で最適な運用計画を得ることが可能になる。
【0045】凹凸平滑化手段14は、運用計画演算手段
12又はハイブリッド最適化手段13により得られた運
用計画が、図8(a),(b)の各左側に示すように、
凹状の時間帯を含む場合(凹状計画)あるいは凸状の時
間帯を含む場合(凸状計画)に、そのときの最適化のた
めの評価値(目的関数や適応度)が悪化せず、且つ制約
条件が遵守されるという制限の下に、図8(a),
(b)の各右側に示すように平滑化を行うものである。
このような処理は、整数値遺伝子でコード化された遺伝
的アルゴリズムを、より効果的に作用させるために導入
したヒューリスティックな(経験に基づく発見的な)も
のであり、主に(2)式のy(i)を含む項の改善に効
果を発揮する。特に、最適化のための評価値が悪化しな
いときのみ平滑化し、そうでない場合は平滑化しないと
いう点が重要である。
【0046】上下平滑化手段15も凹凸平滑化手段14
と同様の機能を有するものであり、運用計画演算手段1
2又はハイブリッド最適化手段13により得られた運用
計画が、図9(a),(b)の各左側に示すように、上
り階段状の時間帯を含む場合(上り階段状計画)あるい
は下り階段状の時間帯を含む場合(下り階段状計画)
に、そのときの最適化のための評価値(目的関数や適応
度)が悪化せず、且つ制約条件が遵守されるという制限
の下に、図9(a),(b)の各右側に示すように平滑
化を行うものである。
【0047】次に、ローカル制御模擬計画手段7の機能
につき説明する。これは、プラントにローカル制御コン
トローラがあり、プラント運用制御装置3からは制御で
きない機器が含まれている場合に、ローカル制御を模擬
し、その結果を運用計画演算手段12に入力して運用計
画演算手段12がプラントの運用計画を演算できるよう
にするためのものである。
【0048】例えば、図10は図2とは異なり、2つの
送水ポンプ58A,B、及び2つの配水池59A,Bを
有するものであるが、浄水池57から配水池59Aへの
送水計画を最適化する場合、時刻iでの需要量に相当す
る配水区61での配水需要予測値qd(i)に対し、いま
配水池59Aから配水池59Bへの送水を行う送水ポン
プ58Bの起動・停止といった制御は、送水計画とは無
関係にローカル制御がかかっているものとする。
【0049】この場合、配水池59Bにとっての需要予
測値は需要予測を行えば得られるが、配水池59Aにと
っての需要予測に相当する量(配水池59Aから配水池
59Bへの送水量)が分からない。そこで、送水ポンプ
58Bの制御シーケンスはいろいろ考えられるが、いま
は、ごく単純に配水池59Bの水位h2(i)が、適切に
定義した上下限値Hu,Hlに対し、Hu<h2(i)なら
送水ポンプ58Bが全台停止し、Hl>h2(i)なら送
水ポンプ58Bが全台起動するものとする。iは1分周
期であったり、5分周期であったり、適切な計算周期を
設定するものとする。なお、配水池59Bのh2(i)の
計算には(6)式を用いればよい。
【0050】すると、どのタイミングで送水ポンプ58
Bが起動又は停止するかが分かるので、運用計画演算手
段12で演算する最適運用計画の計画刻み(1時間単位
で24時間分、あるいは30分単位で24時間分など、
その計画における単位時間のこと)に模擬して得られた
結果を整えて、運用計画演算手段12へ入力すれば、送
水ポンプ58Aによる最適又は近似的な送水計画を演算
することができる。
【0051】例えば、1:34に100〔m3/h〕の送水
ポンプ58Bが3:00まで起動するとき、運用計画演
算手段12への入力として、1時台は、100×(60
−34)/60=43.3 となり、2時台は100と
なる。このような演算を必要な時間分行えば、配水池5
9Aから見た需要予測値にあたる配水池59Bへの送水
量がわかり、最適又は近似的に最適な送水計画を演算す
ることができる。
【0052】このように、プラント1にローカル制御コ
ントローラがあり、プラント運用制御装置3からは制御
できない機器が含まれている場合でも、ローカル制御模
擬計画手段7によりローカル制御を模擬し、その結果を
運用計画演算手段12に入力することにより最適な運用
計画を得ることができる。
【0053】次に、需要予測部5が有する時間遅れ修正
手段8の機能につき説明する。この時間遅れ修正手段8
は、プラント1の設置場所がプロセスコントローラ2に
よる制御応答が一定レベル以上悪くなるような遠距離で
ある場合に、需要予測部5がその制御応答の時間遅れを
考慮した需要予測を行えるようにするためのものであ
る。
【0054】すなわち、図2において、送水ポンプ58
の設置されている場所と配水池59とが数kmにもわた
って離れており、送水ポンプ58の運転開始の結果が配
水池59に反映されるまでに時間が大きく遅れる場合に
は、この時間遅れを考慮しないと制御上問題を生じるこ
とになる。例えば、送水ポンプ58の設置場所から浄水
が配水池59に到達するのに10分かかる場合、送水量
を11:00に100〔m3/h〕から150〔m3/h〕に増
やさないと配水池59の水位が下限値を下回ってしまう
ような場合、少なくとも10分前にその量を想定して、
送水量を増やしてやる必要がある。
【0055】図11は、時間遅れ修正手段8が行う修正
内容を説明するためのグラフ図であり、(a)は修正前
の需要予測値、(b)は修正後の需要予測値である。ま
た、図12は、図11のグラフ図に対応する図表であ
る。図11及び図12に示すように、需要予測を10分
単位で行っているとすると、送水ポンプ58の設置場所
から浄水が配水池59へ到達するのに10分かかるな
ら、単純に10分、すなわち1ステップ先送りした量を
需要予測値とする。もし、送水ポンプ58の設置場所か
ら浄水が配水池59へ到達するのに20分かかるなら、
2ステップ先送りすればよい。そして、運用計画演算手
段12で演算する最適運用計画の計画刻み(1時間単位
で24時間分、あるいは30分単位で24時間分など、
その計画における単位時間のこと)に合わせて、需要予
測値を加工し使用する。
【0056】例えば、10:00〜11:00までの1
時間単位の需要予測値qd(11)につき時間遅れを考慮
して求める場合、時間遅れを考慮しない場合の配水量
は、 110+130+126+128+110+124=7
28 となる。そして、時間遅れが10分なので1ステップ先
送りすると、その場合の配水量は、 130+126+128+110+124+126=7
44 となる。需要予測部5は、このようにして時間遅れ修正
手段8により時間遅れ分が修正された需要予測を運用計
画演算手段12に出力し、運用計画演算手段12はこの
入力に基づき運用計画を演算する。
【0057】
【発明の効果】以上のように、本発明によれば、プロセ
スデータ記憶部がプラントからプロセスコントローラを
介して入力したプロセスデータを記憶し、この記憶され
た過去の実績値、及び所定の情報に基づき、需要予測部
が需要を予測し、この予測された需要に基づき、運用計
画演算手段が遺伝的アルゴリズムを用いて最適な運用計
画についての演算を行う構成としているので、制御対象
プラントが複雑もしくは大規模なものであったとして
も、最適な運用計画についての演算を高速且つ精度良く
行うことができ、効率的で安定したプラント運用を実現
することができる。
【図面の簡単な説明】
【図1】本発明の実施形態に係るプラント運用制御装置
の構成を示すブロック図。
【図2】本発明の実施形態に係るプラント運用制御装置
の制御対象となるプラントの概略構成図。
【図3】本発明の実施形態における整数値遺伝子の例を
示す図表。
【図4】運用計画演算の一手法である分岐限定法につい
ての説明図。
【図5】本発明の実施形態で用いる遺伝的アルゴリズム
を説明するためのフローチャート。
【図6】図5のステップ6における交叉処理を説明する
ための図表。
【図7】図5のステップ7における突然変位処理を説明
するための図表。
【図8】図1における凹凸平滑化手段14の機能につい
ての説明図であり、(a)は凸状計画に対する平滑化、
(b)は凹状計画に対する平滑化を示している。
【図9】図1における上下平滑化手段15の機能につい
ての説明図であり、(a)は上り階段状計画に対する平
滑化、(b)は下り階段状計画に対する平滑化を示して
いる。
【図10】ローカル制御が行われるプラント機器を有す
るプラントの例を示すぶろっく図である。
【図11】図1の時間遅れ修正手段8が行う修正内容を
説明するためのグラフ図であり、(a)は修正前の需要
予測値、(b)は修正後の需要予測値を示している。
【図12】図11のグラフ図に対応する図表。
【符号の説明】 1 プラント 2 プロセスコントローラ 3 プラント運用制御装置 4 プロセスデータ記憶部 5 需要予測部 6 運用計画部 7 ローカル制御模擬計画手段 8 時間遅れ修正手段 9 初期値入力手段 10 初期値演算手段 11 制約条件緩和手段 12 運用計画演算手段 13 ハイブリッド最適化手段 14 凹凸平滑化手段 15 上下平滑化手段 51 取水ポンプ 52 着水井 53 混和池 54 フロック形成池 55 沈殿池 56 ろ過池 57 浄水池 58 送水ポンプ 59 配水池 60 配水ポンプ 61 配水区
───────────────────────────────────────────────────── フロントページの続き (72)発明者 山 崎 謙 一 東京都府中市東芝町1番地 株式会社東芝 府中事業所内 (72)発明者 芦 木 達 雄 東京都港区芝浦一丁目1番1号 株式会社 東芝本社事務所内 Fターム(参考) 5B056 BB91 HH00 5H004 GA34 GA36 GB01 GB08 HA05 HB05 KA16 KC07 KC10 KC23 KD67

Claims (10)

    【特許請求の範囲】
  1. 【請求項1】プラント内のプラント機器に対する制御を
    予め作成してある運用計画に基づき行うプラント運用制
    御装置において、 前記プラントからプロセスコントローラを介して入力し
    たプロセスデータを記憶するプロセスデータ記憶部と、 前記プロセスデータ記憶部に記憶された過去の実績値、
    及び所定の情報に基づき、需要を予測する需要予測部
    と、 前記需要予測部で予測された需要に基づき、遺伝的アル
    ゴリズムを用いて最適な運用計画についての演算を行う
    運用計画演算手段を有する運用計画部と、 を備えたことを特徴とするプラント運用制御装置。
  2. 【請求項2】前記遺伝的アルゴリズムは、コード化され
    た整数値遺伝子を用いるものである、 ことを特徴とする請求項1記載のプラント運用制御装
    置。
  3. 【請求項3】前記運用計画部は、前記遺伝的アルゴリズ
    ム以外の最適化手法により得られた暫定解を、前記遺伝
    的アルゴリズムの初期値として演算し、これを前記運用
    計画演算部に入力させる初期値演算手段を有する、 ことを特徴とする請求項1又は2記載のプラント運用制
    御装置。
  4. 【請求項4】前記運用計画部は、過去のプラント運用に
    おける制御で実際に用いた値を前記プロセスデータ記憶
    部から取り出し、これを前記遺伝的アルゴリズムの初期
    値として、前記運用計画演算手段に入力させるプラント
    実績初期値入力手段を有する、 ことを特徴とする請求項1又は2記載のプラント運用制
    御装置。
  5. 【請求項5】前記運用計画部は、前記運用計画演算手段
    が前記遺伝的アルゴリズムを用いて演算した運用計画に
    ついて、更に、遺伝的アルゴリズム以外の最適化手法に
    基づき最適化演算を行うハイブリッド最適化手段を有す
    る、 ことを特徴とする請求項1乃至4のいずれかに記載のプ
    ラント運用制御装置。
  6. 【請求項6】前記運用計画部は、前記運用計画演算手段
    が演算を行う際の制約条件を緩和する制約条件緩和手段
    を有する、 ことを特徴とする請求項1乃至5のいずれかに記載のプ
    ラント運用制御装置。
  7. 【請求項7】前記運用計画部は、前記運用計画演算手段
    又は前記ハイブリッド最適化手段が演算した運用計画に
    係る経時的な特性曲線中に凹部又は凸部が存在する場
    合、最適化のための評価値を低下させることなく、これ
    らの凹部又は凸部を平滑化する凹凸平滑化手段を有す
    る、 ことを特徴とする請求項1乃至6のいずれかに記載のプ
    ラント運用制御装置。
  8. 【請求項8】前記運用計画部は、前記運用計画演算手段
    又は前記ハイブリッド最適化手段が演算した運用計画に
    係る経時的な特性曲線中に上り階段状部分又は下り階段
    状部分が存在する場合、最適化のための評価値を低下さ
    せることなく、これらの上り階段状部分又は下り階段状
    部分を平滑化する上下階段平滑化手段を有する、 ことを特徴とする請求項1乃至7のいずれかに記載のプ
    ラント運用制御装置。
  9. 【請求項9】プラント内にローカル制御コントローラの
    みにより制御されるローカル機器が含まれている場合
    に、そのローカル制御の模擬を行うローカル制御模擬計
    画部を備えており、 前記運用計画部の運用計画演算手段は、このローカル制
    御模擬計画部の模擬結果を考慮して前記運用計画につい
    ての演算を行うものである、 ことを特徴とする請求項1乃至8のいずれかに記載のプ
    ラント運用制御装置。
  10. 【請求項10】前記需要予測部は、プラント機器同士間
    の離間距離が大きく、あるプラント機器の運転開始後そ
    の運転結果が他のプラント機器に反映されるまでに一定
    レベル以上の時間遅れを生じる場合に、その時間遅れを
    考慮した需要予測を行う時間遅れ修正手段を有するもの
    である、 ことを特徴とする請求項1乃至9のいずれかに記載のプ
    ラント運用制御装置。
JP2002010377A 2002-01-18 2002-01-18 プラント運用制御装置 Expired - Lifetime JP4099334B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2002010377A JP4099334B2 (ja) 2002-01-18 2002-01-18 プラント運用制御装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2002010377A JP4099334B2 (ja) 2002-01-18 2002-01-18 プラント運用制御装置

Publications (2)

Publication Number Publication Date
JP2003216205A true JP2003216205A (ja) 2003-07-31
JP4099334B2 JP4099334B2 (ja) 2008-06-11

Family

ID=27648133

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2002010377A Expired - Lifetime JP4099334B2 (ja) 2002-01-18 2002-01-18 プラント運用制御装置

Country Status (1)

Country Link
JP (1) JP4099334B2 (ja)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009048524A (ja) * 2007-08-22 2009-03-05 Toshiba Corp プラント最適運転制御システム
JP2011220668A (ja) * 2010-03-24 2011-11-04 Metawater Co Ltd 循環式流動焼却炉の温度制御装置及びその温度制御方法
JP2013250992A (ja) * 2007-05-29 2013-12-12 Palo Alto Research Center Inc 機器環境においてジョブの計画を作成する方法及び機器のジョブのためのオフライン及びオンラインジョブの計画作成を実行するためのシステム

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH01277902A (ja) * 1988-04-28 1989-11-08 Fujitsu Ltd スケジューリング問題解決支援システム
JPH04131600A (ja) * 1990-09-19 1992-05-06 Hitachi Ltd 都市エネルギーシステム
JPH0755469A (ja) * 1993-08-11 1995-03-03 Nippon Telegr & Teleph Corp <Ntt> 多視点ステレオ画像計測法
JPH07134603A (ja) * 1992-09-14 1995-05-23 Neuralware Inc プロセス制御装置及び方法
JPH08317559A (ja) * 1995-05-12 1996-11-29 Toshiba Corp 電力系統運用訓練シミュレータ
JPH1161893A (ja) * 1997-08-20 1999-03-05 Toshiba Corp 上水道プラントの制御装置
JPH11126102A (ja) * 1997-10-23 1999-05-11 Nkk Corp ガス需要量の予測方法
JP2000132204A (ja) * 1998-10-21 2000-05-12 Toshiba Corp プラント制御装置
JP2001022437A (ja) * 1999-07-13 2001-01-26 Toshiba Corp プラント制御装置およびプラント制御プログラムを記録したコンピュータ読み取り可能な記録媒体
JP2001055763A (ja) * 1999-08-18 2001-02-27 Toshiba Corp 上水道プラントの広域最適水運用装置
JP2001083107A (ja) * 1999-09-14 2001-03-30 Asahi Kasei Corp 膜の物質特性解析装置
JP2001084004A (ja) * 1999-09-14 2001-03-30 Mitsubishi Electric Corp 機器最適運転支援方法および装置
JP2001214470A (ja) * 2000-01-31 2001-08-07 Toshiba Corp 最適水運用制御装置

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH01277902A (ja) * 1988-04-28 1989-11-08 Fujitsu Ltd スケジューリング問題解決支援システム
JPH04131600A (ja) * 1990-09-19 1992-05-06 Hitachi Ltd 都市エネルギーシステム
JPH07134603A (ja) * 1992-09-14 1995-05-23 Neuralware Inc プロセス制御装置及び方法
JPH0755469A (ja) * 1993-08-11 1995-03-03 Nippon Telegr & Teleph Corp <Ntt> 多視点ステレオ画像計測法
JPH08317559A (ja) * 1995-05-12 1996-11-29 Toshiba Corp 電力系統運用訓練シミュレータ
JPH1161893A (ja) * 1997-08-20 1999-03-05 Toshiba Corp 上水道プラントの制御装置
JPH11126102A (ja) * 1997-10-23 1999-05-11 Nkk Corp ガス需要量の予測方法
JP2000132204A (ja) * 1998-10-21 2000-05-12 Toshiba Corp プラント制御装置
JP2001022437A (ja) * 1999-07-13 2001-01-26 Toshiba Corp プラント制御装置およびプラント制御プログラムを記録したコンピュータ読み取り可能な記録媒体
JP2001055763A (ja) * 1999-08-18 2001-02-27 Toshiba Corp 上水道プラントの広域最適水運用装置
JP2001083107A (ja) * 1999-09-14 2001-03-30 Asahi Kasei Corp 膜の物質特性解析装置
JP2001084004A (ja) * 1999-09-14 2001-03-30 Mitsubishi Electric Corp 機器最適運転支援方法および装置
JP2001214470A (ja) * 2000-01-31 2001-08-07 Toshiba Corp 最適水運用制御装置

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013250992A (ja) * 2007-05-29 2013-12-12 Palo Alto Research Center Inc 機器環境においてジョブの計画を作成する方法及び機器のジョブのためのオフライン及びオンラインジョブの計画作成を実行するためのシステム
JP2009048524A (ja) * 2007-08-22 2009-03-05 Toshiba Corp プラント最適運転制御システム
JP2011220668A (ja) * 2010-03-24 2011-11-04 Metawater Co Ltd 循環式流動焼却炉の温度制御装置及びその温度制御方法

Also Published As

Publication number Publication date
JP4099334B2 (ja) 2008-06-11

Similar Documents

Publication Publication Date Title
Makaremi et al. Optimization of pump scheduling program in water supply systems using a self-adaptive NSGA-II; a review of theory to real application
JP4365598B2 (ja) 広域プラントの最適運用制御装置
JP2013059743A (ja) 海水淡水化プラントシステム
CN108805329B (zh) 一种梯级水库实现实时调度的方法和系统
Wardlaw et al. Comparison of genetic algorithm and linear programming approaches for lateral canal scheduling
JP5703965B2 (ja) 水運用システム及び水運用方法
CN103761385A (zh) 一种多热源环状管网的优化设计方法
CN101789053A (zh) 一种供水型水库生态调度技术
JP2007070829A (ja) 送水運用制御装置
JP2003216205A (ja) プラント運用制御装置
JP2001214470A (ja) 最適水運用制御装置
Fang et al. Optimal operation of multi-storage tank multi-source system based on storage policy
JP2003113739A (ja) エネルギ供給設備の運転計画システム
CN115983489A (zh) 一种随机模糊环境下的热泵与集中供热互补优化运行方法
CN116025026A (zh) 一种城市水箱联动控制方法
CN116415824A (zh) 一种考虑不同需求的多层级梯级水库模拟调度系统及方法
CN116307094A (zh) 一种基于多目标粒子群算法的城市供水优化调度方法
KR101993449B1 (ko) 수처리시설 계장계측제어를 통한 에너지최적화시스템
Kazantzis et al. A new methodology for optimizing the daily operations of a pumping plant
JP3940571B2 (ja) プラント自動制御装置
Fang et al. Optimising pump system with constant and variable speed pumps: case study
JP2006018363A (ja) プラント運用最適化装置
JP4300092B2 (ja) プラント運用計画導出装置
JP6603156B2 (ja) 下水道設備運転支援装置、そのシステムおよび方法
JP7292559B1 (ja) 水処理場運転支援システムおよび水処理場運転支援方法

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20040310

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20060203

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20060207

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20060410

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20060609

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20060720

A911 Transfer to examiner for re-examination before appeal (zenchi)

Free format text: JAPANESE INTERMEDIATE CODE: A911

Effective date: 20060925

A912 Re-examination (zenchi) completed and case transferred to appeal board

Free format text: JAPANESE INTERMEDIATE CODE: A912

Effective date: 20061020

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20080116

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20080317

R151 Written notification of patent or utility model registration

Ref document number: 4099334

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151

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

Free format text: PAYMENT UNTIL: 20110321

Year of fee payment: 3

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

Free format text: PAYMENT UNTIL: 20120321

Year of fee payment: 4

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

Free format text: PAYMENT UNTIL: 20130321

Year of fee payment: 5

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

Free format text: PAYMENT UNTIL: 20130321

Year of fee payment: 5

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

Free format text: PAYMENT UNTIL: 20140321

Year of fee payment: 6

EXPY Cancellation because of completion of term