(57)【要約】      （修正有）【課題】パラメータ調整をできる限り簡単に行え、また与えられたホワイトボックスモデルのパラメータの可同定性を判断しながらパラメータを調整するようにする。【解決手段】プロセスは、ダイナミクスが複数の物質間の相互の変化で表現される様な限定された生物学的・化学プロセスであり、プロセスにおけるプロセス外部入力と複数の物質の計測データを収集してデータを保持するプロセス外部入力収集手段２１と、プロセス出力収集手段２２とから構成されるプロセスデータ収集手段２と、物質収支と反応速度によって構成されるプロセスモデルに対して、プロセス外部入力データを入力して、シミュレーションを行うプロセスシミュレーション手段３と、プロセス出力表示手段１０と、プロセス状態遷移図表示手段１１と、パラメータ値を入力するパラメータ値入力手段と、を具えたプロセスモデルのパラメータ調整支援装置とその方法。【選択図】図４
2005346714.tif 000002
【特許請求の範囲】
【請求項１】
外乱および操作量を含むプロセス外部入力と、プロセス状態量と、前記プロセス外部入力と前記プロセス状態量との関数であるプロセス観測出力と、前記プロセス外部入力を計測する外部入力センサーと、前記プロセス観測出力を計測するプロセス出力センサーと、を有する任意のプロセスにおいて、
前記プロセスにおける前記プロセス外部入力と前記プロセス状態が複数の物質から構成され、前記プロセスは、ダイナミクスが前記複数の物質間の相互の変化で表現される様な限定された生物学的・化学プロセスであり、
前記外部入力センサーの計測データを収集して所定の期間に亘ってデータを保持するプロセス外部入力収集手段と、前記プロセス出力センサーの計測データを収集して所定の期間に亘ってデータを保持するプロセス出力収集手段とから構成されるプロセスデータ収集手段と、
複数のパラメータからなるパラメータセット（パラメータベクトル）θεＲｐ（ｐ：パラメータ数）を持ち、物質収支と反応速度によって構成されるプロセスモデルに対して、前記プロセスデータ外部入力収集手段を通して供給されるプロセス外部入力データを入力することにより、プロセスの時間的および空間的な挙動のシミュレーションを行うプロセスシミュレーション手段と、
前記プロセスシミュレーション手段によるプロセス出力シミュレーション値と前記プロセス出力値とを表示するプロセス出力表示手段と、
前記プロセスモデルの状態の相互作用と状態の変化を表すプロセス状態遷移図表示手段と、
前記プロセス状態遷移図に付属する前記パラメータセットのパラメータ値を入力するパラメータ値入力手段と、
を具えたことを特徴とするプロセスモデルのパラメータ調整支援装置。
【請求項２】
前記パラメータセット入力手段は、前記パラメータセットの要素である各パラメータの物理化学的意味を説明するヘルプ機能を有する
ことを特徴とする請求項１に記載のプロセスモデルのパラメータ調整支援装置。
【請求項３】
前記プロセス状態遷移表示手段から数式へリンクを張り、前記プロセス状態遷移表示手段と前記プロセスモデルを表す数式とを対応づける機能を有する
ことを特徴とする請求項１に記載のプロセスモデルのパラメータ調整支援装置。
【請求項４】
前記プロセス状態遷移図表示手段は、表示されるプロセス状態変数を示す図やプロセス状態間を結合する線の大きさや色が、前記パラメータセット入力手段によって入力されたパラメータ値の大きさに依存して可変となるように修正されている
ことを特徴とする請求項１記載のプロセスモデルのパラメータ調整支援装置。
【請求項５】
前記プロセス状態遷移表示手段と前記パラメータセット入力手段を前記プロセスシミュレーション手段とは異なる装置に実装し、ＷＷＷ（World  Wide  Web  ）ブラウザ上から遠隔でパラメータセットを入力することができるように修正されている
ことを特徴とする請求項１に記載のプロセスモデルのパラメータ調整支援装置。
【請求項６】
外乱および操作量を含むプロセス外部入力と、プロセス状態量と、前記プロセス外部入力と前記プロセス状態量との関数であるプロセス観測出力と、前記プロセス外部入力を計測する外部入力センサーと、前記プロセス観測出力を計測するプロセス出力センサーと、を有する任意のプロセスにおいて、
前記プロセスにおける前記プロセス外部入力と前記プロセス状態が複数の物質から構成され、前記プロセスは、ダイナミクスが前記複数の物質間の相互の変化で表現される様な限定された生物学的・化学プロセスであり、
前記外部入力センサーの計測データを収集して所定の期間に亘ってデータを保持するプロセス外部入力収集工程と、前記プロセス出力センサーの計測データを収集して所定の期間に亘ってデータを保持するプロセス出力収集工程とから構成されるプロセスデータ収集工程と、
複数のパラメータからなるパラメータセット（パラメータベクトル）θεＲｐ（ｐ：パラメータ数）を持ち、物質収支と反応速度によって構成されるプロセスモデルに対して、前記プロセスデータ外部入力収集工程を通して供給されるプロセス外部入力データを入力することにより、プロセスの時間的および空間的な挙動のシミュレーションを行うプロセスシミュレーション工程と、
前記プロセスシミュレーション手段によるプロセス出力シミュレーション値と前記プロセス出力値とを表示するプロセス出力表示工程と、
前記プロセスモデルの状態の相互作用と状態の変化を表すプロセス状態遷移図表示工程と、
前記プロセス状態遷移図に付属する前記パラメータセットのパラメータ値を入力するパラメータ値入力工程と、
を具えたことを特徴とするプロセスモデルのパラメータ調整支援方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、下水処理プロセスや食品加工プロセス等の生物学的プロセスおよび石油精製等の化学プロセスなどのプロセス系を中心とする一般的なプロセスの挙動を模擬するプロセスモデルのパラメータ調整支援装置、及び方法に関する。
【背景技術】
【０００２】
上下水道の水質プロセスなどの水処理プロセスや食品加工プロセス等の微生物の代謝を利用した生物学的プロセスや石油精製プロセスなどの化学プロセスでは、近年プロセスの挙動を模擬するためのプロセスモデル／プロセスシミュレータを利用することが多くなってきている。（なお、ここで、プロセスシュミレータという言葉は、微分方程式や偏微分方程式などの数式で記述されるプロセスモデルをコンピュータなどで計算できるようにコンピュータ上に実装したもののことを示すことが多い。しかし、ここではプロセスモデルとプロセスシュミレータを同一視し、文脈に応じてより適切な言葉を用いる。）
これは、現実のプロセスを詳細に模擬することのできるプロセスシミュレータがあれば、例えば、以下の様なことを行うことができるからである。
【０００３】
（１）プロセスを設計する場合の事前検討にプロセスシミュレータを用いて、最適なプロセスを設計する（設計支援）。
【０００４】
（２）プロセスの運転管理や維持管理をする場合に、プロセスシミュレータによる挙動解析により、より適切な運転管理／維持管理を行うための検討を行う（運転管理支援）。
【０００５】
（３）様々な状況を想定したシナリオシミュレーションを行うことによって、予め起こり得る危険な現象を予測する（シナリオシミュレーション）。
【０００６】
（４）プロセスシミュレータを利用して、オンラインで計測できる変数から、オンラインで計測できない変数の値を推定する（ソフトウエアセンサー、オブザーバ）。
【０００７】
（５）プロセスシミュレータを組み込んだ自動制御コントローラを設計する（モデルベース制御）。
【０００８】
（６）プロセスシミュレータを用いて、経済コストと品質（化学プロセス等の場合）、経済コストと環境コスト（水処理プロセスなどの環境プロセスの場合）の最適化を図る（短期のプロセス最適化）。
【０００９】
（７）プロセスシミュレータを用いて、プロセスの設計段階からプロセスを減価償却し廃棄する段階までの間のライフサイクルコスト（ＬＣＣ）の計算や、ライフサイクルアセスメント（ＬＣＡ）を行う（長期のプロセス最適化）。
【００１０】
（８）プロセスに熟知した専門家以外には理解しづらいプロセスの内部現象をシミュレータによって模擬して可視化することにより、プロセス技術者の教育や訓練に用いる（教育／運転訓練）。
【００１１】
このように、プロセスモデル／プロセスシミュレータは様々な利用用途を潜在的に持っているが、プロセスモデルが実際に現実の様々な問題に適用できるか否かは、プロセスモデルの持つパラメータを現実のプロセスに適合する様にうまく調整できるか否かにかかっている。すなわち、実世界（リアルワールド）で現実に起こっている現象の何らかの形の縮約であるプロセスモデルは、リアルワールドを精度よく近似することができて初めてその存在意義がある。
【００１２】
ところで、プロセスモデルの作り方は、ブラックボックス的なアプローチとホワイトボード的なアプローチの２つに大別される。
【００１３】
ブラックボックス的アプローチは、実際のプロセスから入手したデータを解析することによって数学的にモデルを構築する方法であり、このモデルをブラックボックスモデルと呼ぶことが多い。ブラックボックス的アプローチの代表的なものには、制御理論の一分野であるシステム同定法や、人間の脳の働きを数学的にモデル化したニューラルネットワーク、統計学の分野で用いられる多変量解析などがある。ブラックボックス的なアプローチは、化学プロセスや生物学的プロセスなどのように、機械システムや電磁気システムと比較して複雑な内部構造を持つプロセスをモデル化する場合には有望なアプローチである。このアプローチでは、実在のデータから実際の現象に適合する様にモデルを構築するため、リアルワールドを精度よく近似することが期待できるが、一方、以下の様な問題点がある。
【００１４】
（Ａ）実プロセスの時系列データが存在しない場合はモデル化ができない、そのため、プロセス設計などのプロセス構築前の事前検討に用いることができない。（Ｂ）基本的には実データに適合する様にモデル化を行うため、過去に生じていない現象を模擬することが困難な場合が多い、そのため、シナリオシミュレーションや異常診断に用いる場合には限界がある（汎化能力や予測能力の精度に限界がある）。
【００１５】
（Ｃ）ブラックボックスモデルが持つパラメータは、（数学的な意味やシステム論的な意味は持つが）、通常物理化学的な意味を持たないため、物理量や化学量との対応付けが難しい。そのため、専門家以外の人へのモデルの説明が困難であり、専門家でない一般の人々には受け入れられがたい。そのため、普及させることが難しい。
【００１６】
一方、ホワイトボックス的なアプローチは、物理化学的な洞察に基づいてプロセスモデルを構築する方法であり、この方法によって構築されたモデルはホワイトボックスモデル、物理モデル、物理化学モデルなどと呼ばれる。ホワイトボックスモデルは、ブラックボックスモデルの持つ欠点である上記（Ａ）〜（Ｃ）を解決できるという意味では魅力的な方法である。しかし、ホワイトボックスモデルは物理化学法則に従って構築されているにもかかわらず、これを現実のプロセスに適合させることは、（特にモデルが複雑な場合には）しばしば極めて困難なものとなる。この理由は、ホワイトボックスモデルの持つ複数のパラメータを現実のプロセスに適合する様に調整する統一的な手順が明確でないことが多いためである。そのため、しばしば、プロセス技術者が経験的知識を用いて試行錯誤で調整を行うことが多い、しかし、このような方法は、以下の様な問題点を持つ。
【００１７】
（Ａ’）プロセスの内部構造に関する詳細な知識がない人にはパラメータ調整の指針を立てることができない。また、指針を示されてもその意味を理解することが難しい。従って、結果的にホワイトボックスモデルは「絵に書いた餅」となる可能性が高い。
【００１８】
（Ｂ’）具体的にパラメータの調整を行う以前に、そもそも原理的にホワイトボックス的アプローチで構築されたプロセスモデルのパラメータの値を一通りに決めることができるか否かが明らかでない（構造的可同定性の問題）。（なお、ここで、原理的にパラメータの値を一意に決めることができるか否かという概念を可同定性という。特に、現実のプロセスを完全（厳密）に表現できるモデルとノイズなどを含まない完全なデータが与えられた場合にパラメータの値を一意に決めることができるモデルを構造的に可同定なモデルという。構造的可同定性の問題は、一般には、モデルのもつパラメータ空間Θ⊆Ｒｐ 、（ｐ：パラメータ数）からモデルの出力（センサーによる観測量）空間Ｙ⊆Ｒｐ 、（ｐ：出力数）への写像の逆写像が一意に存在するか否かという問題に帰着される。しかし、個々のプロセスモデルに対して、構造的可同定性を調べることは容易ではない。）
そのため、構造的に可同定でないモデルでは、複数のパラメータの異なる値の組み合わせに対しても同じ応答（出力）が得られることがあり、パラメータ同定（調整）の指針を示すことが難しくなることがある。
【００１９】
（Ｃ’）実際に対象とするプロセスから得られる応答データがノイズに乱されていたり、（プロセスの入力が変化しないために）応答データがあまり変化しなかったりする場合には、パラメータ調整を行うための十分な情報が得られない。そのため、調整できるパラメータの数が限られてしまう（実用的可同定性の問題） 。（なお、ここで、構造的に可同定なモデルであっても、得られるデータが限られていたり、データにノイズが含まれている場合には、そのモデルのパラメータを同定できなくなる場合がある。このように、モデル構造のみに依存する可同定性ではなく（モデル構造と）データの質に依存する可同定性を実用的可同定性という。実際のプロセスでは、実用的可同定であることが望まれる。）
ホワイトボックスモデルをリアルワールドに適合させるためには、（Ａ’）〜（Ｃ’）の様な問題があるにも関わらず、しばしば、ホワイトボックスモデルを扱う人（特にホワイトボックスモデルの作成者）は、物質収支、エネルギー収支、反応過程や反応速度等の物理化学的な諸性質にのみに興味を集中させ、パラメータ同定／調整の問題を単なる最終的なチューニングの問題として軽視する場合があった。そのため、ホワイトボックスモデルを利用する人はこのパラメータ調整のための試行錯誤を繰り返し、これに多大の時間を費やしてしまうということがあった。そして、場合によってはパラメータの調整に失敗し、結果的にホワイトボックスモデルが使われなくなることがしばしば生じていた。
【発明の開示】
【発明が解決しようとする課題】
【００２０】
本発明は、上記の様な問題点を鑑みて、ホワイトボックスモデルをリアルワールドに適合させることのできる実際的な調整方法を提供することを目的とする。
【００２１】
特に、上記（Ａ’）に関連して、プロセスの内部構造を熟知していないプロセス技術者がパラメータ調整をできる限り簡単に行える様な調整の支援装置を提供すること、および、上記（Ｂ’）と（Ｃ’）に関連して、与えられたホワイトボックスモデルのパラメータの可同定性を判断しながら、パラメータを調整するパラメータ調整装置を提供することである。
【課題を解決するための手段】
【００２２】
本発明に係るプロセスモデルのパラメータ調整支援装置は、外乱および操作量を含むプロセス外部入力と、プロセス状態量と、前記プロセス外部入力と前記プロセス状態量との関数であるプロセス観測出力と、前記プロセス外部入力を計測する外部入力センサーと、前記プロセス観測出力を計測するプロセス出力センサーと、を有する任意のプロセスにおいて、前記プロセスにおける前記プロセス外部入力と前記プロセス状態が複数の物質から構成され、前記プロセスは、ダイナミクスが前記複数の物質間の相互の変化で表現される様な限定された生物学的・化学プロセスであり、
前記外部入力センサーの計測データを収集して所定の期間に亘ってデータを保持するプロセス外部入力収集手段と、前記プロセス出力センサーの計測データを収集して所定の期間に亘ってデータを保持するプロセス出力収集手段とから構成されるプロセスデータ収集手段と、
複数のパラメータからなるパラメータセット（パラメータベクトル）θεＲｐ（ｐ：パラメータ数）を持ち、物質収支と反応速度によって構成されるプロセスモデルに対して、前記プロセスデータ外部入力収集手段を通して供給されるプロセス外部入力データを入力することにより、プロセスの時間的および空間的な挙動のシミュレーションを行うプロセスシミュレーション手段と、
前記プロセスシミュレーション手段によるプロセス出力シミュレーション値と前記プロセス出力値とを表示するプロセス出力表示手段と、
前記プロセスモデルの状態の相互作用と状態の変化を表すプロセス状態遷移図表示手段と、
前記プロセス状態遷移図に付属する前記パラメータセットのパラメータ値を入力するパラメータ値入力手段と、
を具えたことを特徴とする。
【００２３】
本発明によれば、パラメータ調整者が、前記プロセスモデルの持つパラメータの物理化学的な意味をプロセス状態遷移図に基づいて考えながらパラメータの調整を行うことができ、パラメータ調整を容易なものとすることができる。
【００２４】
また、本発明は、前記パラメータセット入力手段は、前記パラメータセットの要素である各パラメータの物理化学的意味を説明するヘルプ機能を有することを特徴とする。
【００２５】
本発明によれば、パラメータ調整者に対して、前記プロセスモデルのパラメータの物理化学的意味を詳細かつ分かり易く提示することにより、パラメータの調整の指針を支援することができる。そして、これにより、プロセスモデルをより現実のプロセスにより近いものとすることができる。
【００２６】
また、本発明は、前記プロセス状態遷移表示手段から数式へリンクを張り、前記プロセス状態遷移表示手段と前記プロセスモデルを表す数式とを対応づける機能を有することを特徴とする。
【００２７】
本発明によれば、パラメータ調整者が前記プロセスモデルの数学的な記述と前記プロセス状態遷移図との対応関係を容易に理解することができ、プロセスモデルの持つ意味に対する理解を深めることができる。そのため、プロセスモデルに関する数学的素養のあるパラメータ調整者に対して、パラメータの持つ意味をより深く理解させることができ、ひいてはパラメータ調整者がプロセスモデルを現実のプロセスにより精密に適合させることを可能にすることができる。
【００２８】
また、本発明は、前記プロセス状態遷移図表示手段は、表示されるプロセス状態変数を示す図やプロセス状態間を結合する線の大きさや色が、前記パラメータセット入力手段によって入力されたパラメータ値の大きさに依存して可変となるように修正されていることを特徴とする。
【００２９】
本発明によれば、前記パラメータセットの各パラメータの持つ意味を視覚的に理解することが可能になる。そのため、例えばパラメータのデフォルト値（初期値）を与えた場合のプロセス状態遷移図とパラメータ調整後のプロセス状態遷移図と同時に表示することなどにより、パラメータ調整者がどのパラメータをどの程度調整したのかを視覚的に容易に理解することが可能になる。そのため、特に数多くのパラメータを含むようなプロセスシミュレータのパラメータ調整をより簡単なものとすることができる。
【００３０】
また、本発明は、前記プロセス状態遷移表示手段と前記パラメータセット入力手段を前記プロセスシミュレーション手段とは異なる装置に実装し、ＷＷＷ（World  Wide  Web  ）ブラウザ上から遠隔でパラメータセットを入力することができるように修正されていることを特徴とする。
【００３１】
本発明によれば、プロセスシミュレータを組み込んだプロセスコントローラ（例えばＩＭＣ：Internal  Model  Control  ）を設計する場合のように、プロセスコントローラを実際のプロセスが存在する現場に設置しプロセスモデルのパラメータを監視室などから設定するような状況において、インターネット技術を利用して遠隔でパラメータ調整を容易に行うことができる。これにより、遠隔操作によるプロセスコントローラチューニングが可能となる。
【００３２】
また、ロセスモデルのパラメータ調整支援方法は、プロセス外部入力（外乱および操作量）と、プロセス状態量と、前記プロセス外部入力と前記プロセス状態量との関数であるプロセス観測出力と、前記プロセス外部入力を計測する外部入力センサーと、前記プロセス観測出力を計測するプロセス出力センサーと、を有する任意のプロセスにおいて、
前記プロセスにおける前記プロセス外部入力と前記プロセス状態が複数の物質から構成され、前記プロセスは、ダイナミクスが前記複数の物質間の相互の変化で表現される様な限定された生物学的・化学プロセスであり、
前記外部入力センサーの計測データを収集して所定の期間に亘ってデータを保持するプロセス外部入力収集工程と、前記プロセス出力センサーの計測データを収集して所定の期間に亘ってデータを保持するプロセス出力収集工程とから構成されるプロセスデータ収集工程と、
複数のパラメータからなるパラメータセット（パラメータベクトル）θεＲｐ（ｐ：パラメータ数）を持ち、物質収支と反応速度によって構成されるプロセスモデルに対して、前記プロセスデータ外部入力収集工程を通して供給されるプロセス外部入力データを入力することにより、プロセスの時間的および空間的な挙動のシミュレーションを行うプロセスシミュレーション工程と、
前記プロセスシミュレーション手段によるプロセス出力シミュレーション値と前記プロセス出力値とを表示するプロセス出力表示工程と、
前記プロセスモデルの状態の相互作用と状態の変化を表すプロセス状態遷移図表示工程と、
前記プロセス状態遷移図に付属する前記パラメータセットのパラメータ値を入力するパラメータ値入力工程と、
を具えたことを特徴とする。
【発明の効果】
【００３３】
以上のように、本発明の構成により得られる効果を要約すると以下の少なくとも２点である。
【００３４】
（１）常微分方程式や偏微分方程式で表される複数のパラメータを持つプロセスシミュレータにおいて、入手可能なデータの質に応じて自動的に同定可能なパラメータセットを抽出し、その同定を行うことができる。
【００３５】
（２）下水処理プロセス、上水処理プロセス、石油化学プロセスなどの生物学的プロセスや化学プロセスを表すプロセスモデルの持つパラメータ調整を、視覚情報を提供することによって容易にすることができる。
【００３６】
これら２点の効果により、最終的には、対象とするプロセスの物理化学モデルを現実の対象プロセスに適合することができ、［従来の技術］で述べた様に、プロセスモデルを用いた様々な検討が可能になる。
【発明を実施するための最良の形態】
【００３７】
図１は、本発明の実施形態の基本的な構成を、対象とするプロセスを下水処理プロセスである場合を想定して、下水処理プロセスと共に示したものである。
【００３８】
図１に示す下水処理プロセスパラメータ調整装置は、制御対象としての下水処理プロセス１と、下水処理プロセス１への流入下水量や流入下水水質等の外乱データと曝気風量や返送汚泥量などの操作量データから成るプロセス外部入力データと、各反応槽に設置されたセンサーによって計測している溶存酸素（ＤＯ）や浮遊固形物量（ＭＬＳＳ）やｐＨ、あるいはＣＯＤ（化学的酸素要求量）などの観測出力データを、予め決められた周期で収集し、プロセスデータサーバへ保持するプロセスデータ収集手段２と、下水処理プロセス１の物理化学的な知見に基づいて構成されるプロセスモデルを用いて下水処理プロセス１のＣＯＤ（化学的酸素要求量）等の指標で計測される有機物量、ＤＯ、ＭＬＳＳ、アンモニア態窒素（ＮＨ４−Ｎ）、硝酸態窒素（ＮＯ３−Ｎ）やリン酸態リン（ＰＯ４−Ｐ）などの各種水質の挙動をシミュレーションすることのできるプロセスシミュレーション手段３と、プロセスシミュレーション手段３で用いているプロセスモデルの持つ、反応速度に関連する最大比増殖速度、半飽和定数、加水分解速度定数などの動力学パラメータや物質の収支に関連する収率などの化学量論パラメータなどのパラメータの初期値を設定するパラメータ初期値設定手段４と、プロセスシミュレーション手段３で用いているプロセスモデルから導出することのできるパラメータ感度関数モデルにパラメータ初期値設定手段４によって前記動力学パラメータや化学量論パラメータなどの初期値を設定して、初期設定値近傍でのパラメータの感度を解析するパラメータ感度解析手段５と、パラメータ感度解析手段５によって解析されたパラメータの感度の時間変化を表示するパラメータ感度表示手段６と、パラメータ感度解析手段５によって解析された全てのパラメータの中から高感度なパラメータセットを抽出する高感度パラメータセット抽出手段７と、高感度パラメータセット抽出手段７によって抽出された高感度なパラメータセットの各パラメータを、プロセスデータ収集手段２によるプロセス出力データとプロセスシミュレーション手段３によるプロセス出力シミュレーションデータの値の誤差が小さくなる様に調整／同定するパラメータ同定手段８と、プロセスデータ収集手段２のプロセス外部入力データとプロセス出力データと、プロセスシミュレーション手段３によって計算されるプロセス出力シミュレーションデータの誤差の計算と保存を行いパラメータ同定手段８の評価に用いる誤差計算および保存手段９と、プロセス出力データ収集手段２２から供給されるプロセス出力データとプロセスシミュレーション手段３によって計算されるプロセス出力シミュレーションデータを同時に表示するプロセス値／プロセスシミュレーション値表示手段１０と、から構成される。
【００３９】
また、下水処理プロセス１は、最初沈殿池１１と、生物反応槽１２と、最終沈殿池１３と、下水処理プロセス１の各施設への操作量であるブロワ１４１と返送ポンプ１４２と余剰汚泥引き抜きポンプ１４３の３つのアクチュエータと、プロセス外乱と考えられる流入下水の量および各種流入下水水質を計測するプロセス外乱センサー１５１と曝気風量を計測する風量センサー１５２と返送汚泥量を計測するセンサー１５３と余剰汚泥引抜量を計測するセンサー１５４と、下水処理プロセス１の各設備のＤＯ、ＮＨ４−Ｎ、ＮＯ３−Ｎ、ＰＯ４−Ｐなどの各種水質値を計測するプロセス出力センサー１６１〜１６３と、から構成される。
【００４０】
また、プロセスデータ収集手段２は、下水処理プロセスの外乱と操作量データからなる外部入力データを収集するプロセス外部入力データ収集手段２１と、下水処理プロセスの各種水質データを収集するプロセス出力データ収集手段２２と、から構成される。
【００４１】
図１を用いて、本実施例の作用を説明する。  
本実施例では、説明を簡単にするため、下水処理プロセスの挙動を模擬する単純なプロセスモデルを用いて説明を行うが、本発明の適用対象がこのような単純なプロセスモデルに限定されるわけではない。
【００４２】
まず、予め、プロセスシミュレーション手段３で用いるプロセスモデルから、パラメータ感度解析手段５に用いるパラメータ感度解析モデルを導出する。
【００４３】
最初沈殿池１１と最終沈殿池１３の存在を無視して生物反応槽１２に直接下水が流入するような攪拌式連続反応槽（ＣＳＴＲ）である場合のプロセスモデルは、例えば以下のように与えられる。
【００４４】
【数１】
2005346714.tif 000003
ここで、Ｙhは収率、μmaxは最大比増殖速度、Ｋは酸素濃度に関する半飽和定数、Ｋhは加水分解定数、Ｋxは微生物−難分解性有機物の濃度比に関する半飽和定数、ＫＬは総括酸素移動容量に関する係数（ＫＬ・Ｑbが総括酸素移動容量係数）、ｂは死滅速度、を表すパラメータである。Ｖは反応槽の水量であり一定値であるとする。また、Ｓo2(t)は溶存酸素濃度（ＤＯ）、Ｓc(t)は溶解性有機物濃度、Ｘc(t)は難分解性有機物濃度、Ｘn(t)は微生物濃度、を表す状態変数であり、Ｑ(t)、ＱB(t)は、各々流入量と曝気量を表す入力変数である。各状態変数に添え字のＩＮをつけたものは、各状態変数に対応する流入水質である。また、Ｙ(t)は、溶存酸素濃度（ＤＯ）と溶解性のＣＯＤに対応する出力変数であるが、ＤＯとＭＬＳＳを出力変数とした場合は、以下のように変更される。
【００４５】
【数２】
2005346714.tif 000004
なお、参考のため、このプロセスモデルの反応部分を参考文献［１］［LAWQ Task Group on Mathematical Modeling for Design and Operation of Biological Wastewater Treatment"Activated Sludge Model No.2", IAWQ Scientific Technical Report No.3,(1995)］の記法に従って記述すると、表１のようになる。
【００４６】
【表１】
2005346714.tif 000005
また、表記の簡単化のため、次式の様にプロセスモデルをベクトル形式で表現しておく。
【００４７】
【数３】
2005346714.tif 000006
ここで、調整の対象となるパラメータの候補は、θ＝［Ｙh，μmax，Ｋ，Ｋh，Ｋｘ，ＫＬ，ｂ］の７つである。
【００４８】
さて、状態変数に対するパラメータ感度関数Ｚθ(t)を、次式の様に定義する。
【００４９】
【数４】
2005346714.tif 000007
ここで、Ｚθ(t)のｉ行ｊ列の要素（Ｚθ(t)）ｉｊは、次式で定義される。
【００５０】
【数５】
2005346714.tif 000008
また、センサーで計測している計測値に対してのみパラメータの感度を考える場合は、出力変数に対するパラメータ感度関数Ｙθ(t)を次式で定義する。
【００５１】
【数６】
2005346714.tif 000009
また、パラメータ調整によってプロセスシミュレーションの出力値と実際にセンサーで計測している観測値とを適合させたい出力変数が複数ある場合には、出力変数に対するパラメータ感度関数（１９）の汎関数（Ｒｑ−Ｒ１への写像、ｑは出力変数の数）に対するパラメータ感度関数を定義する。
【００５２】
例えば、（１）〜（４）式の例において、Ｓo2(t)とＳc(t)のプロセスシミュレーション値を実際の値に適合させたい場合は、前記汎関数として、以下の様な評価関数Ｊを設定し、この評価関数に対する感度関数を定義することができる。
【００５３】
【数７】
2005346714.tif 000010
ここで、（・）ｒｅａｌは実際にセンサーなどで計測した計測値を表しｗ１やｗ２は重みを表す。評価関数Ｊに関する感度は、次式で定義できる。
【００５４】
【数８】
2005346714.tif 000011
他の汎関数の選び方として、例えば、
【数９】
2005346714.tif 000012
とすることもできる。この場合の感度関数は、
【数１０】
2005346714.tif 000013
となる。
【００５５】
（２４）式や（２８）式は、（１５）式と重みと観測出力とシミュレーション出力から計算されることに注意する。従って、（１５）式が計算できれば、既知の情報を用いて（２４）式や（２８）式を計算することができる。
【００５６】
なお、パラメータベクトルθの要素の中のいくつかを予め調整の対象としない場合は、θの中のいくつかの要素、例えばθ：＝［μmax  Ｋh  ｂ  ＫＬ］などを新たにパラメータベクトルと見なして、パラメータ感度関数を定義してもよい。
【００５７】
さて、以上から、（１５）式を計算できれば、（２４）式や（２８）式は計算できるので、（１５）式の計算方法を考えればよい。Ｚθを時間ｔで微分し、（７）式を用いると、パラメータ感度関数Ｚθ(t)に関する以下の微分方程式を導くことができる。
【００５８】
【数１１】
2005346714.tif 000014
（３０）式は、プロセスモデルとパラメータ（のデフォルト）値が与えられれば計算可能である。従って、（２４）式や（２８）式も計算可能である。
【００５９】
以下で説明する３パラメータ感度解析手段で用いられるパラメータ感度解析モデルは、（３０）式を（１９），（２４），（２８）式に代入したものや（３０）式自身である。
【００６０】
（３０）式の具体例として、例えば、パラメータベクトルθ：＝［μmax  Ｋh  ｂ  ＫＬ］とした場合の（３０）式の感度関数は、次式のようになる。
【００６１】
【数１２】
2005346714.tif 000015
このような計算を行うことによって、パラメータ感度解析モデルを導出することができる。
【００６２】
なお、このパラメータ感度解析モデルに関して特に注意すべき重要なことは、以下の通りである。
【００６３】
（Ｉ）（３０）式で用いられる∂ｆ／∂Ｚと∂ｆ／∂θは、３プロセスシミュレーション手段で用いるプロセスモデルが存在すれば、そのモデルから解析的に直接計算できる。そのため、パラメータの感度計算において、各パラメータの値を変化させて繰り返しシミュレーションを行うというような面倒な手続きを行う必要がない。
【００６４】
（II）パラメータ感度関数を（２４）式や（２８）式で定義する場合は、（２４）式や（２８）式の感度関数はパラメータの数と同じ要素数を持つベクトルとなっている。そのため、このように定義された感度関数を用いて高感度パラメータセットを抽出することを考える場合には、必ずパラメータ数以下の高感度パラメータセットを抽出することができる。
【００６５】
（III）（３０）で用いられる∂ｆ／∂Ｚと∂ｆ／∂θは、プロセスの状態変数や入力に依存する関数であるため、プロセス入力値が変化すれば、それに伴いこれらの値も変化する。その結果、パラメータ感度関数はプロセス入力値に依存した関数となり、プロセス入力の種類や入力の波形に応じて変化する。以下で説明する高感度パラメータセット抽出手段では、このような入力値の変化に応じて高感度パラメータセットとして抽出するパラメータ数を変化させることができる。そのため、プロセスモデルの構造的および実用的可同定性を陽に意識することなく、可同定なパラメータを抽出できる。
【００６６】
以上の準備のもとで、実施例の作用と効果を説明する。
【００６７】
まず、下水処理プロセス１における外乱センサー１５１で計測される流入量Ｑ(t)と複数の流入水質成分、例えば、溶存酸素濃度ＤＯ、溶解性有機物Ｓｃ(t)（溶解性のＣＯＤとして計測される）、浮遊性有機物Ｘｃ(t)（全ＣＯＤ−溶解性ＣＯＤの成分として計測されたものにα＜１を乗じたもの）、従属栄養性微生物Ｘｈ(t)（全ＣＯＤ−溶解性ＣＯＤの成分として計測されたものに１−αを乗じたもの）をプロセス外乱入力データ収集手段２１によってプロセスデータサーバの所定の領域に所定の計測周期で格納する。同時に下水処理プロセス１におけるの操作量のセンサー１５２〜１５４で収集した操作量の計測値をプロセス外部入力データ収集手段２１によってプロセスデータサーバの所定の領域に所定の計測周期で格納する。また、センサー１６１〜１６３で計測したＤＯ、ＭＬＳＳ、等の水質値をプロセス出力データ収集手段２２によってプロセスデータサーバの所定の領域に格納する。
【００６８】
同時あるいは予め、パラメータ初期値設定手段４によって、プロセスシミュレーション手段３とパラメータ感度解析手段５で用いられる、プロセスモデルおよびパラメータ感度解析モデルが持っているパラメータ初期値を設定する。ここでは、（パラメータの感度はパラメータ初期値にも依存するため）、予め物理化学的な考察によって想定される範囲のパラメータ値を設定しておくことが好ましい。例えば、下水処理プロセスモデルの生物反応モデルとして、ＩＷＡ（International  Water  Association ）は、活性汚泥モデルNo.１〜No.３（ＡＳＭ１〜ＡＳＭ３）を公表しているが、これらのモデルでは、パラメータのデフォルト値が与えられており、このような場合には、デフォルト値を設定しておけば良い。また、このようなデフォルト値が与えられていない場合には、例えば、物質Ａが物質Ｂ内に含まれる含有率というパラメータは０〜１内の値であるというように、理論的に可能な範囲の値を設定しておく。
【００６９】
次に、プロセス外部入力データ収集手段２１によって収集された外部入力データをプロセスシミュレーション手段３のプロセスモデルに入力することによって、プロセスの状態値を計算する。例えば、（１）〜（５）式のプロセスモデルを考える場合は、Ｓo2IN(t)、ＳcIN(t)、ＸcIN(t)、ＸhIN(t)やＱ(t)、などを入力し、Ｓo2(t)、Ｓc(t)、Ｘc(t)、Ｘh(t)などの値を計算する。
【００７０】
次に、プロセスシミュレーション手段３によって計算されたプロセスの状態値とプロセス外部入力データ収集手段２１によって収集された外部入力データをパラメータ感度解析手段５に入力する。この際、プロセスの状態値の中でオンライン計測可能なデータを、プロセスシミュレーション手段３によって計算されたプロセス状態値の代わりに代入してもよい。パラメータ感度解析手段５では、このようなデータが入力されると、（１５）式を計算することができるので、これを計算し、さらに（２４）式や（２８）式を計算する。例えば、（１）〜（５）式のプロセスモデルを考えた場合は、まず、（３１）〜（３３）式の様な感度解析モデルが得られる。ここで、（３１）式は、（３２）〜（３３）式に、μmax、Ｋｈ、ｂ、ＫＬのパラメータ初期値と、Ｓo2(t)、Ｓc(t)、Ｘc(t)、Ｘh(t)のプロセス状態の計算値あるいは計測値、およびＱB(t)などのプロセス外部入力値を与えることにより、計算することができる。そして、（３１）式が計算できれば、（２４）式や（２８）式を計算することもできる。このようにして、５パラメータ感度解析手段が実行される。
【００７１】
次に、パラメータ感度表示手段６では、パラメータ感度解析手段５によって、解析されたパラメータ感度の時系列データを表示する。（１５）式の感度解析結果を全て表示する場合は、（感度解析を行う）パラメータの個数×状態変数の個数分の感度解析結果が表示される。また、（２４）式や（２８）式の感度解析結果を表示する場合は、パラメータの個数分の感度解析結果が表示される。この際、外部入力データや計算あるいは計測された状態変数データを同時に表示してもよい。例えば、（１）〜（５）式のプロセスモデルを考えた場合、
【数１３】
2005346714.tif 000016
に対する感度解析結果は、図２に示すようなグラフとして表示される。
【００７２】
次に、高感度パラメータセット抽出手段７では、パラメータ感度解析手段５とパラメータ感度表示手段６の結果に基づき、高感度なパラメータセットを抽出する。
【００７３】
具体例として、例えば親出願出願時の請求項２に記載の発明を用いる場合には、まず、図２のパラメータ感度関数を予め決めた周期Ｔでサンプリングし、時系列データに変換する。この概念図を図３に示す。図３に示す様に（２４）式や（２８）式を用いるとパラメータの個数分の感度関数ベクトルを生成することができる。例えば、（１）〜（５）式のプロセスモデルを考えた場合は、μmax、Ｋh、ｂ、ＫＬに対応する４つの感度関数ベクトルを生成する。これを、それぞれｅμmax、ｅKn、ｅb、ｅKLとする。次にこれらの生成されたベクトルを並べることによって、次式の様な感度関数行列Ｅを作成する。
【００７４】
【数１４】
2005346714.tif 000017
次にこの感度関数行列Ｅに対して、主成分分析、あるいはほぼ同じことであるが、ＥγＥの特異値分解を行う。（なお、ここで、ＥγＥに対する特異値分解は主成分分析の計算途中で現れ、この特異値分解が計算できれば、主成分分析に必要な情報は全て得られる。）
すると次式の様な形に変換することができる。
【００７５】
【数１５】
2005346714.tif 000018
このような変換を行った場合、もしパラメータ感度関数同士のお互いの相関が強い場合、例えば図２において、μmaxとｂとＫＬのように相関が強い場合は、λｉ２のいくつかは、０に近い非常に小さな値になる。図２に示したような感度関数グラフのような場合は、λ３２とλ４２は非常に小さな値になる。従って、この場合λ１２とλ２２の２つに対応する２個の高感度パラメータを抽出することができる。この場合主成分分析により、パラメータベクトル自身も以下の様に変換される。
【００７６】
【数１６】
2005346714.tif 000019
そのため、この場合は、高感度パラメータセットはθmoaの第１成分と第２成分の２つのパラメータからなるベクトルとして選ばれる。つまり、
【数１７】
2005346714.tif 000020
が高感度パラメータとして選択される。自動的にこの操作を行う場合は、λｉ２に対する閾値λ２を設定しておき、この閾値以上のものを高感度パラメータセットとして抽出することにすればよい。これが、高感度パラメータセット抽出手段７の一例である。
【００７７】
このように高感度パラメータセットを抽出すると以下の様な効果が得られる。（１）数値的に計算できる定量的な基準によって高感度パラメータセットを抽出できる。
【００７８】
（２）プロセスモデルの構造的可同定性や実用的可同定性を陽に意識することなく高感度パラメータセットを抽出できる。もし、プロセスモデルが構造的に可同定でない場合には、上記の特異値分解（主成分分析）の特異値のいくつかは必ず０になり、必然的にパラメータ数未満の高感度パラメータセットしか得られない。また、実用的可同定性は得られるデータの質に依存するが、上記の手続から分かるように、データの質が向上すると、閾値λ２以上の値を持つ特異値の数が多くなり、結果として調整できるパラメータの数が増加する。逆にデータの質が悪い場合には、それに応じて少ない数のパラメータのみが調整可能になる。このような操作を定量的かつ自動的に行うことができる。
【００７９】
（３）高感度パラメータセットを抽出する場合に選ばれる変換パラメータは、変換前の全てのパラメータの線形結合で表され、選ばれる変換パラメータの数は変換前のパラメータの数以下である。選ばれる変換パラメータの数が変換前のパラメータの数より少ない場合には、変換前のパラメータの値は一意に決定されない。そのため、高感度パラメータセット抽出後のパラメータ同定において、物理化学的な考察を合わせて行うことにより、変換前のパラメータの値が一意に決定されないという自由度を利用して、あるパラメータ値だけが物理化学的にありえないような極端な値にならないようにパラメータ調整を行うことができる。
【００８０】
その他の高感度パラメータセット抽出手段７の具体例として、例えば請求項３に記載の発明を用いる場合は、Ｅの要素であるｅμmax、ｅKh、ｅb、ｅKLを抽出するところまでは、請求項２に記載の発明の場合と同様のことを行う。以降のアルゴリズムは、例えば次の様にして実行できる。
【００８１】
Step１
Ｅの要素、ｅμmax  、ｅＫh、ｅb、ｅKLを感度の高いもの順に並べる。例えば、‖ｅ*‖i（＊＝μmax、Ｋh、ｂ、ＫＬ、‖・‖iはｉノルム、ｉ＝１，…，∞）などのノルムを用いて感度の高さを評価し、このノルムの大きなものから順に並べる様にすればよい。そして、例えば、感度の高いものから順にｅ1，ｅ2，ｅ3，ｅ4とする。
【００８２】
Step２
Ｅ2：＝［ｅ1，ｅ2］として特異値分解を行う。この特異値の値が２つとも閾値λ２以上であれば、これを高感度パラメータとして採用する。特異値のうち１つが閾値λ２以下であれば、ｅ2は高感度パラメータセットとして採用せず、Ｅ2＝ｅ1と再定義する。
【００８３】
Step３
Ｅi+1：＝［Ｅi，ｅi+1］としてステップ２へ戻る。ｉがパラメータ数と等しくなったら終了する。
【００８４】
但し、親出願出願時の請求項３に記載の発明の場合のアルゴリズムは、上記のアルゴリズムに限られるわけではなく、高感度パラメータの定義が正確になされており、その中で互いの独立性が強いものを選ぶアルゴリズムであれば任意のものでよい。
【００８５】
このように高感度パラメータセットを抽出すると、親出願出願時の請求項２に記載の発明の場合による高感度パラメータセット抽出手段の効果の（１）と（２）に加えて、物理的に意味を持つパラメータを直接調整対象のパラメータとすることができる。これにより、上記の（３）の効果は失われるが、相互に独立性の高いパラメータでかつ感度の高いものが高感度パラメータセットとして抽出されることになり、従属的なパラメータや感度の低いパラメータは選択されなくなり、パラメータ調整を効率良く行える高感度パラメータセットを抽出できる。例えば、微生物の増殖速度と死滅速度の様な逆相関に近い関係を持つパラメータが共に採用されることは少なくなり、通常は感度の高い増殖速度を表すパラメータが高感度パラメータセットの要素として選択されることになる。
【００８６】
また、高感度パラメータセット抽出手段７の他の例として、親出願出願時の請求項４に記載の発明を用いる場合には、（２４）式や（２８）式のようなパラメータと同じ個数の感度関数を６パラメータ感度表示手段によってディスプレイ上に表示し、その中から感度の高いものを抽出する。あるいは、（１５）式の（パラメータ数×状態変数の数）分の感度関数を６パラメータ感度表示手段によってディスプレイ上に表示し、その中から、パラメータ個数を越えない範囲で、適切なパラメータを抽出する。例えば、状態変数の数がパラメータ数以下である場合は、各状態変数に対して最も高感度なパラメータを１つずつ選択するというようなことが考えられる。具体例として、（１）〜（４）式のモデルを考える場合、例えば、Ｓo2(t)に対してKL，Ｓo(t)に対してｂ，Ｘh(t)に対してμmax，Ｘc(t)に対してＫhなどとして、高感度パラメータを選択することができる。そして、例えばＸh(t)の値が計測されておらず分からない場合などには、ＫＬとｂとＫhを調整するという風にすればよい。
【００８７】
このような高感度パラメータセットの抽出方法は、人間系を介している為効率は悪いが、人間の持つ知識を高感度パラメータセットの抽出に反映できるという利点を持つ。例えば、「ＫＬというパラメータは総括酸素移動容量係数を表し、直接溶存酸素濃度Ｓo2(t)に影響を与え、その他の状態変数には直接的な影響を与えない」という知識を人間は持っている。そのため、Ｓo2(t)のみの調整を行う場合にＫＬを高感度パラメータセットの要素として選択するというような合理的な高感度パラメータセット抽出を行うことが可能になる。
【００８８】
次に、パラメータ同定手段８では、高感度パラメータセット抽出手段７によって抽出した複数の高感度パラメータセットに対して、その値の調整（同定）を行う。例えば、親出願出願時の請求項５に記載の発明を用いる場合には、例えば評価関数（あるいは適合度）として（２４）式や（２８）式の誤差の２乗平方和をとり、その値が最小となるようにヒューリスティックな探索を行う。（なお、ここで、適合度の場合は誤差の２乗平方和の逆数などとして定義する場合が多いため、その場合は最大。）
具体例として、μmaxとＫhのパラメータ同定に免疫的アルゴリズムを用いる場合は、以下の様な手順で行う。
【００８９】
Step１
μmaxとＫhを離散データとみなし、複数の離散データからなる初期抗体群を発生させる。具体的には、μmaxとＫhに対して物理化学的に意味がある範囲内の実数値を適当に離散化したデータ集合を記憶細胞とみなし、この記憶細胞から初期抗体群を生成する。
【００９０】
Step２
抗体同士の親和度と抗体と抗原の親和度を計算する。抗体同士の親和度の計算において、各抗体（パラメータ値の候補）に対するシミュレーション出力値と実際のセンサーで計測した値との二乗誤差などの評価値を用いる。これは９誤差計算および保存手段から供給された値に基づいて行う。
【００９１】
Step３
抗原との親和度の高い抗体を記憶細胞に追加する。
【００９２】
Step４
次世代に残る抗体の期待値を計算し、期待値の低い抗体を消滅させる。
【００９３】
Step５
消滅させた抗体に代わる新しい抗体を算出し、適当な確率で突然変異させる。
【００９４】
Step６
基準を満たせば終了。満たさなければ、Step２へ戻る。
【００９５】
このような操作により、パラメータ同定手段８を実行することができる。
【００９６】
このようなメタヒューリスティクスを用いるパラメータ同定手段８は、比較的高速に実行でき、数学的に非凸な探索問題にも適用できるため、パラメータに関して複雑な関数であるプロセスモデルのパラメータをオンラインで調整する場合には、特に有効な方法である。
【００９７】
パラメータ同定手段８の他の実行方法として、例えば親出願出願時の請求項６に記載の発明を適用する場合には、以下の様に行うことができる。まず、ｈ（Ｚ，θ）：＝Ｙ(t)と定義する。これは、（５）式や（６）式で定義されるＹ(t)は、θを含む（１）〜（４）式の出力であるためＹ(t)自身θの関数であり、形式的にｈ（Ｚ，θ）と書くことが可能である。この時、評価関数を（２０）式で定義すると、例えば非線形パラメータ同定の繰り返しアルゴリズムとして、以下の様なものを採用することができる。
【００９８】
【数１８】
2005346714.tif 000021
ここで、ｉは繰り返し回数を示し、ｉ＝１，２，…である。初期値であるθ０は、物理化学的に意味のある値を設定しておく。また、Ｋiは繰り返し回数がｉ回目である場合の勾配ゲインを表す行列であり、正定行列である。Ｋiの選び方には、多くの自由度があるが、例えば、Gauss-Newton法を利用する場合は、
【数１９】
2005346714.tif 000022
の様なアルゴリズムとなる。但し、Ｑは正定行列である。この中で利用する式は（１）〜（５）式を用いれば、全て計算可能であるため、これは実行可能なアルゴリズムである。
【００９９】
このようなパラメータ同定手段は、非線形最小二乗法や最大推定法などを根拠にしているため、二乗誤差最小の意味や光度の意味で最適なパラメータ値を求めることができる。また、例えばＡＳＭ２などの様にパラメータのデフォルト値が予め与えられている場合には、その近傍を探索していくため、物理化学的に有り得ないような異常な値になる可能性が少なくなるという効果がある。
【０１００】
他の方法として、パラメータ同定手段８の他の実行に親出願出願時の請求項７に記載の発明を採用する場合には、高感度パラメータセット抽出手段７で抽出したパラメータの最も高感度なものから順にプロセス値／プロセスシミュレーション値表示手段１０の表示結果を参考にして調整を行う。この場合、パラメータの変化量に対する状態変数の変化量の割合を示す（１５）式などを用いれば、パラメータの値を変えた場合に出力がどのように変化するかを容易に計算することが可能であるため、パラメータ値を変化させた場合の１０プロセス値／プロセスシミュレーション値表示手段の表示を容易に変更できる。そのため、パラメータの値の変更に伴い、（１５）式を用いて状態変数の値を修正し、その結果を表示しながらパラメータの調整（同定）を行うことができる。
【０１０１】
このようなパラメータ同定手段は、物理化学的な知見と照らし合わせながら、高感度なパラメータのみをなるべく効率良く調整することができるという効果がある。
【０１０２】
［その他の実施例］
図４に示す下水処理プロセスパラメータ調整装置は、制御対象としての下水処理プロセス１と、下水処理プロセスへの流入下水量や流入下水水質等の外乱データと曝気風量や返送汚泥量などの操作量データから成るプロセス外部入力データと、各反応槽に設置されたセンサーによって計測している溶存酸素（ＤＯ）や浮遊固形物量（ＭＬＳＳ）やｐＨ、あるいはＣＯＤなどの観測出力データを、予め決められた周期で収集し、プロセスデータサーバへ保持するプロセスデータ収集手段２と、下水処理プロセス１の物理化学的な知見に基づいて構成されるプロセスモデルを用いて下水処理プロセス１のＣＯＤ（化学的酸素要求量）等の指標で計測される有機物量、ＤＯ、ＭＬＳＳ、アンモニア態窒素（ＮＨ４−Ｎ）、硝酸態窒素（ＮＯ３−Ｎ）やリン酸態リン（ＰＯ４−Ｐ）などの各種水質の挙動をシミュレーションすることのできるプロセスシミュレーション手段３と、プロセスシミュレーション手段３で計算され出力されるプロセスシミュレーション値とプロセスデータ収集手段２によって収集された観測出力データを表示するプロセス値／プロセスシミュレーション値表示手段１０と、プロセスシミュレーション手段３で用いているプロセスモデルの反応に従って、プロセスの状態が遷移していく様子を示すプロセス状態遷移図表示手段１１と、プロセス状態遷移図表示手段１１に付属するパラメータ値設定手段１２と、から構成される。
【０１０３】
また、下水処理プロセス１は、図１と同様に構成されている。プロセスデータ収集手段２も図１と同様に構成されている。
【０１０４】
次に、本実施例の作用と効果について説明する。  
図４を用いて、本実施例の作用を説明する。本実施例では、前述の参考文献［１］に記述されている活性汚泥モデルＡＳＭ２をプロセスモデルの例として説明を行う。もちろん、その他の生物学的プロセスや化学プロセスにおいても同様な考え方を用いることができる。
【０１０５】
まず、図１を用いた場合の実施例の作用と同様に、下水処理プロセス１における外乱センサー１５１で計測される流入量Ｑ(t)と複数の流入水質成分、例えば、溶存酸素濃度ＤＯ、溶解性有機物Ｓｃ(t)（溶解性のＣＯＤとして計測される）。浮遊性有機物Ｘc(t)（全ＣＯＤ−溶解性ＣＯＤの成分として計測されたものにα＜１を乗じたもの）、従属栄養性微生物Ｘh(t)（全ＣＯＤ−溶解性ＣＯＤの成分として計測されたものに１−αを乗じたもの）、リン酸態リン（ＰＯ４−Ｐセンサーで計測される）、アンモニア態窒素（ＮＨ４−Ｎセンサーで計測される）、硝酸態窒素（ＮＯ３−Ｎセンサーで計測される）などをプロセス外部入力データ収集手段２１によってプロセスデータサーバの所定の領域に所定の計測周期で格納する。同時に下水処理プロセス１における操作量のセンサー１５２〜１５４で収集した操作量の計測値をプロセス外部入力データ収集手段２１によってプロセスデータサーバの所定の領域に所定の計測周期で格納する。また、センサー１６１〜１６３で計測したＤＯ、ＭＬＳＳ、ＰＯ４−Ｐ、ＮＨ４−Ｎ、ＮＯ３−Ｎ等の水質値をプロセス出力データ収集手段２２によってプロセスデータサーバの所定の領域に格納する。
【０１０６】
次に、プロセスシミュレーション手段３にプロセス外部入力データ収集手段２１からプロセス外部入力データを供給することによって、プロセス状態の時間的および空間的変化のシミュレーションを行う。
【０１０７】
次に、プロセス出力データ収集手段２２から供給されるプロセス出力値とプロセスシミュレーション手段３によって計算したプロセス出力シミュレーション値を同時にプロセス値／プロセスシミュレーション値表示手段１０に表示する。
【０１０８】
一方、図５〜図７に示すようなプロセス状態遷移図をプロセス状態遷移図表示手段１１に表示しておく 。（なお、ここで、図５〜図７のプロセス状態遷移図は、ＡＳＭ２の消化、脱窒、脱リン、有機物除去に関連する状態遷移を示したものである。以降の説明には図５のみを用いるが、ＡＳＭ２で考慮しているプロセスの反応の複雑さを強調するため、参考にために図６と図７も載せておく。）
同時に、図５に示す様にプロセス状態遷移図上にパラメータ値を設定することのできる入力部を設けておく。これが、パラメータ値設定手段１２の設定を行うための入力部である。プロセス状態遷移図表示手段１１をパソコンのWindows上に構築しておくような場合には、このパラメータ値設定手段１２の化学量論パラメータに関する入力部は、例えば、テキストボックスに入力できるようにしておいてもよいし、メニュー形式で値を選択できるようにしておいてもよい。また、図が複雑に成りすぎるため全ての表示を行っていないが、図５の右に示す様に、全ての状態遷移を表す矢印部分には、動力学パラメータを設定できるようにしておく。この動力学パラメータ設定部もテキストボックスやメニュー形式で設定できるようにしておけばよい。パソコンなどのWindows上の表示面積が小さい場合には、状態遷移の矢印部分をクリックすると、動力学パラメータの設定を行える画面が例えばポップアップウィンドウとして表示されるようにしておいてもよい（図５の右部分参照）。
【０１０９】
さて、以上の準備のもとで、まず、プロセス値／プロセスシミュレーション値表示手段１０に表示されたプロセス出力値とプロセス出力シミュレーション値の差（誤差）を確認する。例えば、ＡＳＭ２を用いて下水処理プロセスのシミュレーションを行った場合は、リン酸態リンＰＯ４−Ｐやアンモニア態窒素ＮＨ４−Ｎや硝酸態窒素ＮＯ３−Ｎなどの状態変数の濃度の出力値とシミュレーション値の誤差を確認する。そして、最も誤差の大きい状態変数を選ぶ。ここでは、例えば、仮に硝酸態窒素濃度の誤差が最も大きかったとしておく。
【０１１０】
次に最も誤差の大きい状態変数、この場合硝酸態窒素濃度に影響を与えるパラメータをプロセス状態遷移図表示手段１１によって探し出す。例えば、化学量論パラメータであるＹautやＹhが硝酸態窒素濃度に影響を与えることがわかる。また、図５には表示されていないが、ＡＳＭ２では、図５の硝化の部分（Ａの部分）および脱窒の部分（Ｂ１とＢ２の部分）の反応速度は、次式で表されている。
【０１１１】
【数２０】
2005346714.tif 000023
ここで、動力学パラメータは、μaut，μh，ηno3，Ｋo2，Ｋnh4，Ｋno3，Ｋpo4，Ｋalk，Ｋf，Ｋaであるため、これらのパラメータは硝酸態窒素濃度に影響を与えることが分かる。
【０１１２】
そこで、１１プロセス状態遷移図表示手段を用いると、硝酸態窒素濃度の値を調整する場合は、Ｙaut，Ｙh，μaut，μh，ηno3，Ｋo2，Ｋnh4，Ｋno3，Ｋpo4，Ｋalk，Ｋf，Ｋaの値のいくつかを調整すればよいことが分かる。ＡＳＭ２などの化学・生物反応モデルの物理化学的な意味と数学的な意味を理解しているパラメータ調整者であれば、この中から最も影響を与えそうなものから順に試行錯誤で調整を行えばよい。
【０１１３】
プロセスモデルに対する物理化学的な知識がないパラメータ調整者のためには、本願に記載の発明によって、パラメータの意味を表示するヘルプ機能を用意する。ヘルプ機能のヘルプ情報の表示例を数例以下に示す。Ｙautについては、「硝化菌に対する収率を表す化学量論パラメータ  この値を大きくすると、硝酸態窒素濃度は上昇し、アンモニア態窒素濃度と酸素濃度は低下する。動力学パラメータを調整しても調整が困難な場合にこのパラメータ値を変更することが望ましい。パラメータ値は０〜１の範囲で設定すること、」などの表示をしておくとよい。Ｋo2については、「溶存酸素濃度に関する半飽和定数を表す動力学パラメータ、この値を小さくすると溶存酸素濃度が比較的低い場合でも好気条件になり易くなり好気条件の反応が進行しやすくなる。逆に嫌気条件での反応は進行しにくくなる。好気条件／嫌気条件に関わる全ての状態変数に影響を与えるため、調整には注意を要する。」などとしておくとよい。また、ηｎｏ３については、「無酸素条件での脱窒反応の最大速度を補正する動力学パラメータ、無酸素条件での脱窒反応のみに関係するため、脱窒に関連する硝酸態窒素濃度などの調整に利用するとよい。」などとしておく。このようなヘルプメッセージがあれば、例えば、脱窒反応に関わる硝酸態窒素濃度ＮＯ３−ＮとＳＦの濃度とＳＡの濃度のプロセス出力シミュレーション値とプロセス出力値が大きくずれており、他の状態変数の値はほぼプロセス出力シミュレーション値とプロセス出力値が一致している場合には、ηｎｏ３を用いて調整すればよいことが分かる。
【０１１４】
このようなヘルプ機能は、調整の対象となるパラメータセットの中から最も適したパラメータを選択するために有効である。
【０１１５】
また、モデルの数学的な知識や物理化学的な知識をある程度持っているパラメータ調整者のためには、本願に記載の発明により、図５〜図７のプロセス状態遷移図と図８の様なプロセスモデルを表す数式とを対応づけておくことが役に立つ。例えば、図５に示す様に、「モデル表示」をクリックすると図８の様なプロセスモデルを表す数式が表示されるようにしておくとよい。
【０１１６】
このようプロセス状態遷移図と数式を対応づける機能を持つことにより、化学反応や生物反応に関する物理化学的および数学的知識を持つパラメータ調整者が調整を行う場合には、数式とプロセス状態遷移図の両方を参考にしながら、調整を行うためのパラメータを選択していくことができる。
【０１１７】
また、プロセスモデルで表現しているプロセス反応の仕組みと現実のプロセスの対応関係をより視覚的に表示したい場合には、本願に記載の発明により、設定したパラメータの値とプロセス状態遷移図の図が連動するようにしておくとよい。例えば、化学量論パラメータに関しては、図９に示す様に、プロセス状態遷移図において状態変数の表示の大きさと化学量論パラメータの値とを連動させればよい。図９では、状態変数である従属栄養性微生物について示しているが、従属栄養性微生物の収率Ｙhを大きくする場合には、それに伴って従属栄養性微生物を表す状態変数の大きさを大きくすればよい。具体的には、収率Ｙhの大きさと従属栄養性微生物を表す円の半径の間に比例関係を持たせておくとすればよい。動力学パラメータに関しては、通常複数の状態変数の値に依存して反応速度が変わるため、例えば各状態変数に対して色を決めた矢印を複数用意し、その矢印の太さや矢印の「矢」の部分の大きさを動力学パラメータと対応させればよい。例えば図１０に示す様に、「酸素濃度に関しては黒、硝酸態窒素濃度に対しては赤、…」などと色を決めておき、これらの濃度が関連する状態遷移については、図５〜図７に示すプロセス状態遷移図の矢印を１つの反応に対して複数並列に表示するようにしておく。そして、例えば、動力学パラメータの一つである半飽和定数は反応がどれぐらい速やかに進行するのかを示すパラメータであるため、半飽和定数が大きい場合には、図１１に示す様に矢印の「矢」の部分を大きくし、半飽和定数が小さい場合には「矢」の部分を小さくするなどの対応づけを行っておけばよい。例えば、酸素濃度に関係する半飽和定数Ｋo2の値は、黒の矢印の「矢」の部分と対応づけられる。また、動力学パラメータの一つである最大比増殖速度は、反応速度の最大値を与えるパラメータであるため、最大比増殖速度が大きい場合には、図１２に示す様に矢印を太く表示し、逆に最大比増殖速度が小さい場合には、矢印を細く表示するなどの対応づけを行っておけばよい。
【０１１８】
このようにパラメータ値とプロセス状態遷移図の視覚的な表示との対応づけを行うことにより、各パラメータの持つ意味を視覚的に容易に理解できる。さらに、パラメータのデフォルト値を与えた場合のプロセス状態遷移図と調整を行ったプロセス状態遷移図を同時に表示するなどの工夫を行えば、特にパラメータの数が多いプロセスモデルの場合において、パラメータ調整者自身が、数多くのパラメータの中のどのパラメータをデフォルト値から既に変更したのかを視覚的に容易に理解できる。
【０１１９】
以上の様なパラメータ調整の支援情報をもとにして、調整すべきパラメータを選択して、実際にパラメータ値を決定した後、これを１１プロセス状態遷移図表示手段に付属する１２パラメータ値設定手段を通して入力し、この値を３プロセスシミュレーション手段のプロセスモデルのパラメータに反映させる。そして、新たに設定されたパラメータ値を用いて３プロセスシミュレーション手段によって更新されたプロセス出力シミュレーション値を計算する。
【０１２０】
以上の操作を繰り返して、最も誤差の大きい状態変数に対するパラメータ調整を行った後、プロセス出力値とプロセス出力シミュレーション値の誤差の大きい順に、同様の操作を繰り返すことによって、（主に化学・生物反応に関係する）プロセスのパラメータ調整（支援）を行うことができる。
【０１２１】
また、プロセスシミュレーションの目的が、例えばInternal  Model  Control等の様に最終的にプロセスを適切に制御する操作量を演算するような場合には、３プロセスシミュレーション手段をコントローラの内部に持っておくことが好ましい。このような場合、本願に記載の発明により、図５〜図７の様なプロセス５状態遷移図表示手段を、例えば対象とするプロセスを監視する監視室あるいは対象プロセスから遠く離れた場所に設置されたオフィスなどにあるパソコンやワークステーション上のＷＷＷブラウザ上に構築し、その値をインターネットや専用回線経由で対象プロセスが存在する場所のコントローラ内部のプロセスシミュレータのパラメータ値に反映できるようにしておくとよい。
【０１２２】
これにより、遠隔操作でプロセスシミュレータのパラメータを調整することができ、プロセスシミュレータがコントローラに組み込まれているため、コントローラのパラメータを遠隔操作で調整していることにもなる。
【０１２３】
なお、以上の実施例は、多くの化学プロセスや生物学的プロセスに対して適用できる。なぜなら、化学プロセスや生物学的プロセスの多くのものは、図５に示す様な複数の物質（の量や濃度）に対応する複数の状態変数間の相互作用を表すCompartmentalモデルと呼ばれるモデルで記述することができ、これらのモデルのパラメータは、基本的に物質の量的な関係を表す化学量論パラメータと、状態変数間の相互作用（化学・生物反応）の反応速度に関連する動力学パラメータから構成されるためである。従って、上記で説明した方法は、化学プロセスや生物学的プロセスにおいては、かなり汎用的に適用することができる。
【図面の簡単な説明】
【０１２４】
【図１】パラメータ調整装置を実施する場合の構成図と対象プロセスの例である下水処理プロセスを併記した図。
【図２】パラメータ感度関数の時間変化の例を示す図。
【図３】パラメータ感度関数ベクトルを構成する方法を示す概念図。
【図４】パラメータ調整支援装置を実施する場合の構成図と対象プロセスの例である下水処理プロセスを併記した図。
【図５】プロセス状態遷移図表示手段とパラメータ値設定手段の例を示す図。
【図６】プロセス状態遷移図表示手段の例２を示す図。
【図７】プロセス状態遷移図表示手段の例３を示す図。
【図８】プロセス状態遷移図表示に関連付けられたプロセスモデル（数学モデル）の例を示す図。
【図９】化学量論パラメータとプロセス状態遷移図表示手段の表示方法の対応付けの例を示す図。
【図１０】反応速度に関連する状態変数とプロセス状態遷移図表示手段の表示方法の対応付けの例を示す図。
【図１１】動力学パラメータとプロセス状態遷移図表示手段の表示方法の対応付けの例１（半飽和定数の場合）を示す図。
【図１２】動力学パラメータとプロセス状態遷移図表示手段の表示方法の対応付けの例２（最大比増殖速度の場合）を示す図。
【符号の説明】
【０１２５】
１  下水道プロセス
１１  最初沈殿池
１２  生物反応槽
１３  最終沈殿池
１４１  ブロワ（アクチュエータ）
１４２  返送ポンプ（アクチュエータ）
１４３  余剰汚泥引き抜きポンプ（アクチュエータ）
１５１  プロセス外乱センサー（プロセス入力センサー）
１５２  曝気風量センサー（プロセス入力センサー）
１５３  返送汚泥量センサー（プロセス入力センサー）
１５４  余剰汚泥引抜量センサー（プロセス入力センサー）
１６１  最初沈殿池状態センサー（プロセス出力センサー）
１６２  生物反応槽状態センサー（プロセス出力センサー）
１６３  最終沈殿池状態センサー（プロセス出力センサー）
２  プロセスデータ収集手段
２１  プロセス外部入力データ収集手段
２２  プロセス出力データ収集手段
３  プロセスシミュレーション手段
４  パラメータ初期値設定手段
５  パラメータ感度解析手段
６  パラメータ感度表示手段
７  高感度パラメータセット抽出手段
８  パラメータ同定手段
９  誤差計算および保存手段
１０  プロセス値／プロセスシミュレーション値表示手段
１１  プロセス状態遷移図表示手段
１２  パラメータ値設定手段
【図１】
2005346714.tif 000024
【図２】
2005346714.tif 000025
【図３】
2005346714.tif 000026
【図４】
2005346714.tif 000027
【図５】
2005346714.tif 000028
【図６】
2005346714.tif 000029
【図７】
2005346714.tif 000030
【図８】
2005346714.tif 000031
【図９】
2005346714.tif 000032
【図１０】
2005346714.tif 000033
【図１１】
2005346714.tif 000034
【図１２】
2005346714.tif 000035