プロセスの状態予測方法
【課題】要求点ベクトルごとに最適な近傍データベクトルの数が自動的に決定され、予測精度の安定化が図られるプロセスの状態予測方法を提供する。
【解決手段】プロセスの操業状態を示す観測データから構成される入力ベクトル及び出力ベクトルが対となって蓄積されたデータベースを作成し、予測したい時点における出力ベクトルに対応する入力ベクトルからなる要求点ベクトルに類似する近傍データベクトルをデータベースから取得し、近傍データベクトルから局所モデルを構築して、予測したい時点における出力ベクトルを求めるプロセスの状態予測方法において、近傍データベクトルが格納された近傍データセットを近傍データ数を変えて複数作成し、複数の近傍データセットについて主成分分析を実施して各近傍データセットごとに前記要求点ベクトルに対するＱ統計量を算出して、Ｑ統計量が最小となる近傍データセットを選択して局所モデルを構築する。
【解決手段】プロセスの操業状態を示す観測データから構成される入力ベクトル及び出力ベクトルが対となって蓄積されたデータベースを作成し、予測したい時点における出力ベクトルに対応する入力ベクトルからなる要求点ベクトルに類似する近傍データベクトルをデータベースから取得し、近傍データベクトルから局所モデルを構築して、予測したい時点における出力ベクトルを求めるプロセスの状態予測方法において、近傍データベクトルが格納された近傍データセットを近傍データ数を変えて複数作成し、複数の近傍データセットについて主成分分析を実施して各近傍データセットごとに前記要求点ベクトルに対するＱ統計量を算出して、Ｑ統計量が最小となる近傍データセットを選択して局所モデルを構築する。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、プラント設備におけるプロセスの状態予測方法に関する。
【背景技術】
【０００２】
あるプロセスについて、その状態を把握する必要があるとき、通常の計測機器では分析に時間がかかり、リアルタイムにプロセス状態を把握できないことがある。また、計測環境や計測対象によっては、計測機器の設置自体が困難な場合もある。プロセスの状態を示す明確な物理モデルが得られる場合は、高精度な推定値を計算によって求めることができるが、プラント設備におけるプロセスは複雑な物理化学現象が複合した形で発現することが殆どであるため、物理モデルで表せない場合が多い。
【０００３】
そこで、近年、計算機ハードウェアやデータベースシステム技術の進歩に伴い、大量データの蓄積と高速検索が可能になったこと等を背景に、“Just-In-Time（ＪＩＴ）モデリング”と呼ばれる局所モデリング手法が注目されている。ＪＩＴモデリングでは、観測したデータをデータベースに蓄積しておき、システムの予測等の必要が生じるたびに、入力である“要求点ベクトル”と関連性の高いデータベクトルをデータベースから近傍データベクトルとして検索し、検索した近傍データベクトルの出力を補間する局所モデルを構成して、“要求点ベクトル”の出力を推定する。この手法では、観測データの更なる蓄積があるたびに既存の局所モデルを廃棄し、再び新たな局所モデルを構築する。
【０００４】
ＪＩＴモデリングでは、予測を行うたびに、データベースから要求点ベクトルと類似するデータベクトルを検索するため、データベースが大規模になると、計算負荷が大きくなりすぎるという問題がある。そのため、ＪＩＴモデリングにステップワイズ法を適用して変数の低次元化を行う大規模データベースオンラインモデリング（ＬＯＭ）という手法が開発されている。例えば、特許文献１、２では、熱反応炉の操業データからなる大規模データベースについて、ステップワイズ法を用いて炉頂ガス温度に対する寄与率が高い変数を選択して当該変数からなる新たなデータベースを作成し、新たなデータベースから取得した近傍データベクトルに基づいて構築した局所モデルを用いて炉頂ガス温度の予測を行っている。
【０００５】
なお、本明細書では、「要求点」と「近傍データ」がそれぞれベクトル量であることを明確にするため、「要求点」を「要求点ベクトル」、「近傍データ」を「近傍データベクトル」と記載する。また、データベクトルの集合である「データベクトル集合」を「データセット」と呼ぶことがある。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２００９−０７６０３６号公報
【特許文献２】特開２００９−０７６０３７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
ＪＩＴモデリングやＬＯＭは要求点の近傍データベクトルを基に局所モデルを構築するため、使用する近傍データ数が重要となる。従来、ＪＩＴモデリングやＬＯＭでは、対象によって試行錯誤的に最適と思われる近傍データベクトルの数を決定していた。しかし、同じ対象であっても要求点ベクトルごとに最適な近傍データベクトルの数は異なると考えられる。つまり、頻繁に発生するようなケースを要求点ベクトルとした場合、多くの近傍データベクトルを集めたほうが予測の安定度が増す。
一方、稀にしか発生しないようなケースを要求点ベクトルとした場合には、近傍データベクトルの数が多くなり過ぎると、関連性の低いデータベクトルまで使用して局所モデルを構築することになる。例えば、特許文献１、２に記載されたプロセスの状態予測方法では、ステップワイズ法で選択した変数について、予め設定した数（特許文献２では６個）だけ要求点の近傍データベクトルを取得しているため、予測精度が不安定であった。
【０００８】
本発明はかかる事情に鑑みてなされたもので、要求点ベクトルごとに最適な近傍データベクトルの数が自動的に決定され、予測精度の安定化が図られるプロセスの状態予測方法を提供することを目的とする。
【課題を解決するための手段】
【０００９】
上記目的を達成するため、本発明は、プラント設備におけるプロセスの操業状態を示す観測データから構成される入力ベクトル及び出力ベクトルが対となったデータベクトルが蓄積されたデータベースを作成し、予測したい時点における出力ベクトルに対応する入力ベクトルからなる要求点ベクトルに類似する前記データベクトルを近傍データベクトルとして前記データベースから少なくとも１つ以上取得し、前記近傍データベクトルから局所モデルを構築して、前記予測したい時点における出力ベクトルを推定するプロセスの状態予測方法において、
前記近傍データベクトルが格納された近傍データベクトル集合を該近傍データベクトルの数を変えて複数作成する工程と、前記複数の近傍データベクトル集合について主成分分析を実施して該各近傍データベクトル集合ごとに前記要求点ベクトルに対するＱ統計量を算出する工程と、前記Ｑ統計量が最小となる前記近傍データベクトル集合を選択して前記局所モデルを構築する工程とを備えることを特徴としている。
【００１０】
主成分分析では、変数間の相関関係を捉えるため、変数の線形結合によって主成分と呼ばれる新たな合成変数を作り出す。この主成分によって、対象とするデータベクトル集合の特徴を最も良く表現する部分空間を得ることができる。Ｑ統計量は、主成分によって張られる部分空間では表現できない部分を表している。つまり、Ｑ統計量は、対象とするデータベクトル集合と要求点ベクトルとの相関関係の非類似度を表し、Ｑ統計量が小さいほど、要求点ベクトルに類似するデータベクトル集合であると判断できる。
【００１１】
図１２は、要求点ベクトルと近傍データベクトルとの相関関係を表したものである。図１２（Ａ）は、ＪＩＴモデリングやＬＯＭの場合を示しており、ベクトル間距離に基づいて近傍データベクトルを選択するため、異なる相関関係を有する近傍データベクトルが選択されるおそれがある。一方、図１２（Ｂ）は、本発明に係るプロセスの状態予測方法の場合を示しており、Ｑ統計量を用いて、要求点ベクトルとデータベクトル集合（データセット）の相関関係を測るため、相関関係の高い○印のデータベクトル集合のみ選択される。
【００１２】
また、本発明に係るプロセスの状態予測方法では、前記近傍データベクトルの数が異なる前記複数の近傍データベクトル集合を作成する際、前記要求点ベクトルとのベクトル間距離が近い前記近傍データベクトルから順に格納して前記近傍データベクトル集合を作成することが好ましく、このようにすることで、より類似度の高い近傍データベクトル集合を選択することが可能となる。
【発明の効果】
【００１３】
本発明に係るプロセスの状態予測方法では、要求点ベクトルごとに作成される、近傍データベクトルの数が異なる複数の近傍データベクトル集合のなかから、要求点ベクトルに対するＱ統計量が最小となる近傍データベクトル集合を選択して局所モデルを構築するので、要求点ベクトルごとに最適な近傍データベクトルの数が自動的に決定され、予測精度の安定化を図ることができる。
【図面の簡単な説明】
【００１４】
【図１】本発明の一実施の形態に係るプロセスの状態予測方法を説明するためのフローチャートである。
【図２】データセットの構成を示すテーブルである。
【図３】要求点ベクトルの構成を示すテーブルである。
【図４】近傍データ数がＮＮＭＡＸである近傍データセットＡの構成を示すテーブルである。
【図５】近傍データ数がＮＮＭＩＮである近傍データセットＢ０の構成を示すテーブルである。
【図６】Ｑ統計量が格納されたＱ値テーブルである。
【図７】同実施の形態に係るプロセスの状態予測方法によって得られた予測値と実測値との相関度を示すグラフである。
【図８】近傍データ数を１２００個とした従来型ＬＯＭによって得られた予測値と実測値との相関度を示すグラフである。
【図９】同実施の形態に係るプロセスの状態予測方法によって１時間後の温度を連続予測した結果を示す時刻歴グラフである。
【図１０】近傍データ数を１２００個とした従来型ＬＯＭによって１時間後の温度を連続予測した結果を示す時刻歴グラフである。
【図１１】同実施の形態に係るプロセスの状態予測方法によって１時間後の温度を連続予測した際に使用した近傍データ数のグラフである。
【図１２】要求点ベクトルと近傍データベクトルとの相関関係を表す模式図であって、（Ａ）はＪＩＴモデリングやＬＯＭの場合、（Ｂ）は本発明に係るプロセスの状態予測方法の場合をそれぞれ示している。
【発明を実施するための形態】
【００１５】
続いて、添付した図面を参照しつつ、本発明を具体化した実施の形態に付き説明し、本発明の理解に供する。
【００１６】
［プロセスの状態予測方法の概略手順］
先ず、本発明の一実施の形態に係るプロセスの状態予測方法の概略手順を以下に示す。
（Ａ１）プラント設備におけるプロセスの操業状態を示す観測データから構成される入力ベクトル及び出力ベクトルが対となったデータベクトルが蓄積されたデータベースを作成する。
（Ａ２）予測したい時点における出力ベクトルに対応する入力ベクトルからなる要求点ベクトルに類似するデータベクトルを近傍データベクトルとしてデータベースから少なくとも１つ以上取得する。そして、近傍データベクトルが格納された近傍データセット（近傍データベクトル集合）を近傍データベクトルの数（以下では、単に「近傍データ数」と記載する。）を変えて複数作成する。
（Ａ３）複数の近傍データセットについて主成分分析を実施して各近傍データセットごとに要求点ベクトルに対するＱ統計量を算出し、Ｑ統計量が最小となる近傍データセットを選択して局所モデルを構築する。そして、当該局所モデルを用いて、予測したい時点における出力ベクトルの推定値を得る。
【００１７】
ここで、本実施の形態に係るプロセスの状態予測方法を構成する主要な手法について基本的な説明をしておく。
［ＪＩＴモデリング］
現在の挙動と近似した挙動が過去に観測されていたならば、現在の挙動が進展する様子は過去のものと近似したものになるであろうと考えることができる。この考え方を再現した予測手法の１つがJust-In-Time（ＪＩＴ）モデリングである。ＪＩＴモデリングは決まったモデルを持たない代わりに、過去のデータベクトルをそのままデータベースとして保持する。プロセスの予測が必要となったとき、過去データが蓄積されたデータベースから、要求点ベクトルと類似性の高いデータベクトルを検索し、局所モデルを構築して出力の推定を行う手法である。
【００１８】
対象とするプロセスが非線形かつ動的なプロセスであるとき、次式の回帰モデルでそのプロセスを表すことができる。
【００１９】
【数１】
【００２０】
ここで、プロセスの入力ベクトルｘｋと出力ベクトルｙｋを以下のように定義する。つまり、出力ベクトルｙｋは、ｋ時における入力ベクトルｘｋに対する（ｋ＋ｐ）時における出力、即ち予測値となる。
【００２１】
【数２】
【００２２】
時間の経過と共に、入力ベクトルｘｋと出力ベクトルｙｋのデータベクトルの組が（ｘ１，ｙ１），（ｘ２，ｙ２），…のように、対象とするプロセスから大量に得られ、データベクトル集合｛（ｘｋ，ｙｋ）｝（ｋ＝１，２，…）としてデータベースに蓄積される。ｋは離散化時間である。
【００２３】
予測したい時点における出力ベクトルｙｋｑに対応する入力ベクトルｘｋｑを要求点ベクトルとし、要求点ベクトルと類似性が高い近傍データベクトルを上記データベースから取得する。要求点ベクトルと類似性が高い近傍データベクトルを選択する際の指標としては、次式で示すようなベクトル間距離（ユークリッド距離）などを用いることができる。
【００２４】
【数３】
【００２５】
近傍データベクトル群｛（ｘｋｉ，ｙｋｉ）｝（ｉ＝１，２，…，ｍ）が取得されると、この近傍データベクトル群を用いて局所モデルの構築を行い、出力ベクトルｙｋｑの推定を行う。局所モデルとしては、重回帰モデルや、以下に示す相加平均法あるいは重み付き線形平均法などが用いられる。
【００２６】
【数４】
【００２７】
［ステップワイズ法］
ステップワイズ法は、目的変数に対する影響（寄与率）が小さい説明変数を除外するものであり、変数増加ステップと変数減少ステップとから構成される。以下、ステップワイズ法の手順について説明する。
【００２８】
（Ｂ１）目的変数に対する単寄与率が最大の説明変数を、先に選定した説明変数の中から決定する。具体的には、各説明変数ごとに単回帰モデルを作成して回帰係数を求めた後、（８）式によりＦ値を算出し、Ｆ値が最大となる説明変数を選択する。
【００２９】
【数５】
【００３０】
（Ｂ２）先の手順で決定された現モデルに説明変数を１つ追加することを考える。即ち、現モデルに含まれていない説明変数の中から、現モデルに一つ説明変数を追加した場合についてそれぞれ偏回帰係数を求め、（８）式によりＦ値を算出する。そして、Ｆ値が最大となる説明変数を探索する。
（Ｂ３）最大のＦ値が、前もって決めているＦｉｎ以上である場合は、その説明変数を現モデルに追加する。最大のＦ値がＦｉｎ未満の場合は、ステップワイズ法による選択手順を終了する。
【００３１】
（Ｂ４）新たに説明変数が追加された場合は、今度は逆に今までに取り込んだ説明変数が本当に有用な説明変数であるのかどうか調べる。即ち、現モデルを構成する説明変数の中で寄与率が一番低い説明変数を見つけるために、今までに取り込んだ説明変数を順番に一つずつ取り除いて当該説明変数が無い場合におけるＦ値を計算し、Ｆ値が最小となる説明変数を探索する。
（Ｂ５）最小のＦ値が、前もって決めているＦｏｕｔ（Ｆｉｎ≧Ｆｏｕｔ）未満である場合は、その説明変数を現モデルから削除する。説明変数を削除した場合は、（Ｂ４）、（Ｂ５）を繰り返し、さらに他の説明変数が削除できないか探索する。最小のＦ値がＦｏｕｔ以上の場合は、（Ｂ２）のステップに戻る。
【００３２】
［主成分分析］
主成分分析は、データの特徴抽出及び低次元化を目的とする多変量解析手法であり、変数間の相関関係を捉えるため、変数の線形結合によって得られる主成分と呼ばれる合成変数を使用する。主成分分析では、データを最も良く表現できる方向に第１主成分を設定し、第１主成分と直交する空間上で、第１主成分では表現できないデータの変動を最も良く表現できる方向に第２主成分を設定するという手順で、主成分を次々と設定していく。ここで、データを最も良く表現する方向というのは、主成分得点の分散が最大となる方向という意味である。また、主成分得点とは、主成分が張る部分空間へデータを射影した値である。
【００３３】
［Ｑ統計量］
Ｑ統計量は、データベクトルのうち、主成分によって張られる部分空間では表現できない部分を表す。Ｑ統計量は二乗予測誤差とも呼ばれ、以下のように定義されている。
Ｎ行×Ｍ列のデータ行列Ｘがあるものとする。ここで、Ｍは変数の数、Ｎはサンプル数であり、各変数は標準化されている。
データ行列Ｘを特異値分解すると次式のようになる。
【００３４】
【数６】
【００３５】
ＵとＶは直交行列であり、対角行列Ｓの対角要素には特異値ｓｒが降順に並んでいる。採用する主成分の数をＲとすると、第ｒ主成分は負荷量行列ＶＲの第ｒ列ｖｒで与えられる。
第ｒ主成分得点ｔｒは（１０）式で与えられ、第Ｒ主成分得点までをまとめて表現すると、（１１）式となる。
【００３６】
【数７】
【００３７】
ＴＲを元のＭ次元空間上の座標で表すと、再構築データ行列Ｘ＾は次のようになる。
【００３８】
【数８】
【００３９】
このとき、Ｑ統計量は次式で与えられる。
【００４０】
【数９】
【００４１】
［プロセスの状態予測方法の詳細手順］
続いて、図１のフローチャートに基づいて本実施の形態に係るプロセスの状態予測方法の手順について詳細に説明する。
（Ｃ１）プラント設備におけるプロセスの操業状態を示す観測データから構成される入力ベクトルｘｋ及び出力ベクトルｙｋのデータベクトルの組（ｘｋ，ｙｋ）（ｋ＝１，２，…）が蓄積された大規模データベース１０を作成する。
（Ｃ２）大規模データベース１０について、ステップワイズ法を用いて目的変数に対する寄与率が高い変数を選択して当該変数からなる新たなデータベース１１を作成する（ＳＴ１）。目的変数と説明変数の間に時間遅れが存在する可能性がある場合は、見込まれる最大の時間遅れ変数まで全て選択対象に加える。
作成されるデータベース１１（データセット）の構成を図２に示す。このデータベース１１では、入力変数の数がＭ個、出力変数の数がＬ個、各変数のサンプル数がＫ個とされている。各データは日時に応じたＩＤが付けられ、同じＩＤに属するデータは１つのデータベクトルとして扱われる。
【００４２】
（Ｃ３）予測したい時点における出力ベクトルＹｑに対応する入力ベクトルＸｑからなる要求点ベクトルを設定する（ＳＴ２）。図３に要求点ベクトルの構成を示す。
（Ｃ４）データベース１１に格納されている各データベクトルと要求点ベクトルとのベクトル間距離を（４）式や（５）式を用いて計算し、ベクトル間距離が小さいものから順にＮＮＭＡＸ個の近傍データベクトルを全て収集する。そして、収集した近傍データベクトルを、ベクトル間距離が近い順に近傍データセットＡとして保存する（ＳＴ３）。図４に近傍データセットＡの構成を示す。図４において「Ｎｏ．」が近傍データ数を表している。
（Ｃ５）近傍データセットＡの中から近傍データ数（Ｎｏ．）が１〜ＮＮＭＩＮまでの近傍データベクトルを選択して近傍データセットＢ０を作成する（ＳＴ４）。即ち、要求点ベクトルとのベクトル間距離が近いものから順にＮＮＭＩＮ個の近傍データベクトルを選択する。図５に近傍データセットＢ０の構成を示す。
【００４３】
（Ｃ６）近傍データセットＢ０に対して主成分分析を実施し、負荷量行列ＶＲを求める（ＳＴ５）。具体的には、近傍データセットＢ０をデータ行列Ｘとして特異値分解すればよい。
（Ｃ７）要求点ベクトルｘｑが（１４）式で表されるとすると、要求点ベクトルｘｑを再構築した再構築ベクトルｘ＾ｑは、負荷量行列ＶＲを用いて（１５）式により算出される。従って、近傍データセットＢ０に対するＱ統計量は、（１６）式より得ることができる（ＳＴ６）。算出されたＱ統計量は、図６に示すＱ値テーブルに保存される。
【００４４】
【数１０】
【００４５】
（Ｃ８）Ｑ統計量が算出された近傍データセットＢ０の近傍データ数がＮＮＭＡＸ以上であるかどうか判断される（ＳＴ７）。近傍データ数がＮＮＭＡＸ未満である場合は、近傍データセットＡの内、近傍データセットＢ０に含まれていない近傍データベクトルの中から、さらにＳ個の近傍データベクトルを、近傍データ数（Ｎｏ．）が小さいほうから（要求点ベクトルとのベクトル間距離が近いものから）選択し、近傍データセットＢ０に追加して新たな近傍データセットＢ１を作成する（ＳＴ９）。そして、ＳＴ５のステップに戻る。
（Ｃ９）一方、近傍データ数がＮＮＭＡＸ以上になった場合は、Ｑ値テーブルに基づいて、Ｑ統計量が最小となったデータセットＢｋをデータセットＡから選択する。そして、データセットＢｋに対応する出力ベクトルを、データセットＢｋのＩＤに基づいてデータベース１１から取得して、重回帰モデルや重み付き線形平均法などを用いて局所モデルを構築し、要求点ベクトルに対する出力の推定値を算出する（ＳＴ８）。
【００４６】
以上、本発明の一実施の形態について説明してきたが、本発明は何ら上記した実施の形態に記載の構成に限定されるものではなく、特許請求の範囲に記載されている事項の範囲内で考えられるその他の実施の形態や変形例も含むものである。例えば、上記実施の形態では、ＮＮＭＡＸ個の近傍データセット全てについてＱ統計量を求めているが、ある近傍データセットのＱ統計量が極小値であると判断された時点で、当該近傍データセットを選択してもよい。また、上記実施の形態では、変数を削減するためにステップワイズ法を使用しているが、変数が少ない場合はステップワイズ法を省略してもよい。
【実施例】
【００４７】
本実施の形態に係るプロセスの状態予測方法（以下では、「近傍データ数逐次可変型ＬＯＭ」と呼ぶ。）の効果について検証するため、ごみ処理プロセスにおけるガス化溶融炉の炉頂ガス温度の予測を行った。なお、従来の予測方法と比較するため、従来型ＬＯＭによるガス化溶融炉の炉頂ガス温度予測も併せて行った。
【００４８】
検証に使用したデータは、２年間に亘るごみ処理プロセスにおいて測定された観測データである。取り込んだデータは、ノイズ除去のため、１時間の移動平均フィルタを掛けて平滑化した。サンプリング時間は２０分、総データ数は３８８０９個である。
全ての変数に対してステップワイズ法を適用した場合、厖大な処理時間とコンピュータメモリが必要となると共に、不要な変数はモデリング精度低下の原因となる。そのため、予め炉頂ガス温度と関係があると考えられる説明変数を３７個に絞った後、各変数の遅れ時間を０から５０時間としてステップワイズ法による変数選択を実施した。その結果、２７個の説明変数を選択し、１時間後の炉頂ガス温度の予測を行った。
【００４９】
近傍データ数逐次可変型ＬＯＭにおける最大近傍データ数ＮＮＭＡＸは１２００個、最小近傍データ数ＮＮＭＩＮは１００個、近傍データ数の増加幅Ｓは１０個とした。また、主成分分析に使用した主成分数は１０個とした。
一方、従来型ＬＯＭにおける近傍データ数は、試行錯誤的に最も良い結果が得られる数値とした。
なお、近傍データ数逐次可変型ＬＯＭも従来型ＬＯＭも、局所モデルの構築には重回帰モデルを使用した。
【００５０】
要求点ベクトルを観測データからランダムに２００点選択して、各要求点ベクトルに対して１時間後の炉頂ガス温度の予測を行った。近傍データ数逐次可変型ＬＯＭによる予測値と実測値との相関度を図７に、近傍データ数を１２００個としたときの従来型ＬＯＭによる予測値と実測値との相関度を図８に示す。これらの図から、近傍データ数逐次可変型ＬＯＭは、従来型ＬＯＭに比べて予測値と実測値のバラツキが少なく、相関係数ｒが０．０８程度向上していることがわかる。
なお、図７以降において、目盛に表示されている「Ｔ」は基準温度を表している。即ち、基準温度Ｔからの変化量で表している。
【００５１】
次に、１時間ごとに１時間後の炉頂ガス温度を予測する連続予測を行った結果について説明する。近傍データ数逐次可変型ＬＯＭによる連続予測結果を示す時刻歴グラフを図９に、近傍データ数を１２００個としたときの従来型ＬＯＭによる連続予測結果を示す時刻歴グラフを図１０にそれぞれ示す。また、近傍データ数逐次可変型ＬＯＭによって連続予測した際に使用した近傍データ数のグラフを図１１に示す。
これらの図から、近傍データ数を１２００個とした場合には難しかった急激な温度上昇（時刻歴グラフの５時間目参照）が近傍データ数逐次可変型ＬＯＭでは捉えられていることがわかる。また、近傍データ数逐次可変型ＬＯＭでは、近傍データ数を１２００個とした場合に比べて全体的な誤差も小さくなっている。
【符号の説明】
【００５２】
１０：大規模データベース、１１：データベース
【技術分野】
【０００１】
本発明は、プラント設備におけるプロセスの状態予測方法に関する。
【背景技術】
【０００２】
あるプロセスについて、その状態を把握する必要があるとき、通常の計測機器では分析に時間がかかり、リアルタイムにプロセス状態を把握できないことがある。また、計測環境や計測対象によっては、計測機器の設置自体が困難な場合もある。プロセスの状態を示す明確な物理モデルが得られる場合は、高精度な推定値を計算によって求めることができるが、プラント設備におけるプロセスは複雑な物理化学現象が複合した形で発現することが殆どであるため、物理モデルで表せない場合が多い。
【０００３】
そこで、近年、計算機ハードウェアやデータベースシステム技術の進歩に伴い、大量データの蓄積と高速検索が可能になったこと等を背景に、“Just-In-Time（ＪＩＴ）モデリング”と呼ばれる局所モデリング手法が注目されている。ＪＩＴモデリングでは、観測したデータをデータベースに蓄積しておき、システムの予測等の必要が生じるたびに、入力である“要求点ベクトル”と関連性の高いデータベクトルをデータベースから近傍データベクトルとして検索し、検索した近傍データベクトルの出力を補間する局所モデルを構成して、“要求点ベクトル”の出力を推定する。この手法では、観測データの更なる蓄積があるたびに既存の局所モデルを廃棄し、再び新たな局所モデルを構築する。
【０００４】
ＪＩＴモデリングでは、予測を行うたびに、データベースから要求点ベクトルと類似するデータベクトルを検索するため、データベースが大規模になると、計算負荷が大きくなりすぎるという問題がある。そのため、ＪＩＴモデリングにステップワイズ法を適用して変数の低次元化を行う大規模データベースオンラインモデリング（ＬＯＭ）という手法が開発されている。例えば、特許文献１、２では、熱反応炉の操業データからなる大規模データベースについて、ステップワイズ法を用いて炉頂ガス温度に対する寄与率が高い変数を選択して当該変数からなる新たなデータベースを作成し、新たなデータベースから取得した近傍データベクトルに基づいて構築した局所モデルを用いて炉頂ガス温度の予測を行っている。
【０００５】
なお、本明細書では、「要求点」と「近傍データ」がそれぞれベクトル量であることを明確にするため、「要求点」を「要求点ベクトル」、「近傍データ」を「近傍データベクトル」と記載する。また、データベクトルの集合である「データベクトル集合」を「データセット」と呼ぶことがある。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２００９−０７６０３６号公報
【特許文献２】特開２００９−０７６０３７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
ＪＩＴモデリングやＬＯＭは要求点の近傍データベクトルを基に局所モデルを構築するため、使用する近傍データ数が重要となる。従来、ＪＩＴモデリングやＬＯＭでは、対象によって試行錯誤的に最適と思われる近傍データベクトルの数を決定していた。しかし、同じ対象であっても要求点ベクトルごとに最適な近傍データベクトルの数は異なると考えられる。つまり、頻繁に発生するようなケースを要求点ベクトルとした場合、多くの近傍データベクトルを集めたほうが予測の安定度が増す。
一方、稀にしか発生しないようなケースを要求点ベクトルとした場合には、近傍データベクトルの数が多くなり過ぎると、関連性の低いデータベクトルまで使用して局所モデルを構築することになる。例えば、特許文献１、２に記載されたプロセスの状態予測方法では、ステップワイズ法で選択した変数について、予め設定した数（特許文献２では６個）だけ要求点の近傍データベクトルを取得しているため、予測精度が不安定であった。
【０００８】
本発明はかかる事情に鑑みてなされたもので、要求点ベクトルごとに最適な近傍データベクトルの数が自動的に決定され、予測精度の安定化が図られるプロセスの状態予測方法を提供することを目的とする。
【課題を解決するための手段】
【０００９】
上記目的を達成するため、本発明は、プラント設備におけるプロセスの操業状態を示す観測データから構成される入力ベクトル及び出力ベクトルが対となったデータベクトルが蓄積されたデータベースを作成し、予測したい時点における出力ベクトルに対応する入力ベクトルからなる要求点ベクトルに類似する前記データベクトルを近傍データベクトルとして前記データベースから少なくとも１つ以上取得し、前記近傍データベクトルから局所モデルを構築して、前記予測したい時点における出力ベクトルを推定するプロセスの状態予測方法において、
前記近傍データベクトルが格納された近傍データベクトル集合を該近傍データベクトルの数を変えて複数作成する工程と、前記複数の近傍データベクトル集合について主成分分析を実施して該各近傍データベクトル集合ごとに前記要求点ベクトルに対するＱ統計量を算出する工程と、前記Ｑ統計量が最小となる前記近傍データベクトル集合を選択して前記局所モデルを構築する工程とを備えることを特徴としている。
【００１０】
主成分分析では、変数間の相関関係を捉えるため、変数の線形結合によって主成分と呼ばれる新たな合成変数を作り出す。この主成分によって、対象とするデータベクトル集合の特徴を最も良く表現する部分空間を得ることができる。Ｑ統計量は、主成分によって張られる部分空間では表現できない部分を表している。つまり、Ｑ統計量は、対象とするデータベクトル集合と要求点ベクトルとの相関関係の非類似度を表し、Ｑ統計量が小さいほど、要求点ベクトルに類似するデータベクトル集合であると判断できる。
【００１１】
図１２は、要求点ベクトルと近傍データベクトルとの相関関係を表したものである。図１２（Ａ）は、ＪＩＴモデリングやＬＯＭの場合を示しており、ベクトル間距離に基づいて近傍データベクトルを選択するため、異なる相関関係を有する近傍データベクトルが選択されるおそれがある。一方、図１２（Ｂ）は、本発明に係るプロセスの状態予測方法の場合を示しており、Ｑ統計量を用いて、要求点ベクトルとデータベクトル集合（データセット）の相関関係を測るため、相関関係の高い○印のデータベクトル集合のみ選択される。
【００１２】
また、本発明に係るプロセスの状態予測方法では、前記近傍データベクトルの数が異なる前記複数の近傍データベクトル集合を作成する際、前記要求点ベクトルとのベクトル間距離が近い前記近傍データベクトルから順に格納して前記近傍データベクトル集合を作成することが好ましく、このようにすることで、より類似度の高い近傍データベクトル集合を選択することが可能となる。
【発明の効果】
【００１３】
本発明に係るプロセスの状態予測方法では、要求点ベクトルごとに作成される、近傍データベクトルの数が異なる複数の近傍データベクトル集合のなかから、要求点ベクトルに対するＱ統計量が最小となる近傍データベクトル集合を選択して局所モデルを構築するので、要求点ベクトルごとに最適な近傍データベクトルの数が自動的に決定され、予測精度の安定化を図ることができる。
【図面の簡単な説明】
【００１４】
【図１】本発明の一実施の形態に係るプロセスの状態予測方法を説明するためのフローチャートである。
【図２】データセットの構成を示すテーブルである。
【図３】要求点ベクトルの構成を示すテーブルである。
【図４】近傍データ数がＮＮＭＡＸである近傍データセットＡの構成を示すテーブルである。
【図５】近傍データ数がＮＮＭＩＮである近傍データセットＢ０の構成を示すテーブルである。
【図６】Ｑ統計量が格納されたＱ値テーブルである。
【図７】同実施の形態に係るプロセスの状態予測方法によって得られた予測値と実測値との相関度を示すグラフである。
【図８】近傍データ数を１２００個とした従来型ＬＯＭによって得られた予測値と実測値との相関度を示すグラフである。
【図９】同実施の形態に係るプロセスの状態予測方法によって１時間後の温度を連続予測した結果を示す時刻歴グラフである。
【図１０】近傍データ数を１２００個とした従来型ＬＯＭによって１時間後の温度を連続予測した結果を示す時刻歴グラフである。
【図１１】同実施の形態に係るプロセスの状態予測方法によって１時間後の温度を連続予測した際に使用した近傍データ数のグラフである。
【図１２】要求点ベクトルと近傍データベクトルとの相関関係を表す模式図であって、（Ａ）はＪＩＴモデリングやＬＯＭの場合、（Ｂ）は本発明に係るプロセスの状態予測方法の場合をそれぞれ示している。
【発明を実施するための形態】
【００１５】
続いて、添付した図面を参照しつつ、本発明を具体化した実施の形態に付き説明し、本発明の理解に供する。
【００１６】
［プロセスの状態予測方法の概略手順］
先ず、本発明の一実施の形態に係るプロセスの状態予測方法の概略手順を以下に示す。
（Ａ１）プラント設備におけるプロセスの操業状態を示す観測データから構成される入力ベクトル及び出力ベクトルが対となったデータベクトルが蓄積されたデータベースを作成する。
（Ａ２）予測したい時点における出力ベクトルに対応する入力ベクトルからなる要求点ベクトルに類似するデータベクトルを近傍データベクトルとしてデータベースから少なくとも１つ以上取得する。そして、近傍データベクトルが格納された近傍データセット（近傍データベクトル集合）を近傍データベクトルの数（以下では、単に「近傍データ数」と記載する。）を変えて複数作成する。
（Ａ３）複数の近傍データセットについて主成分分析を実施して各近傍データセットごとに要求点ベクトルに対するＱ統計量を算出し、Ｑ統計量が最小となる近傍データセットを選択して局所モデルを構築する。そして、当該局所モデルを用いて、予測したい時点における出力ベクトルの推定値を得る。
【００１７】
ここで、本実施の形態に係るプロセスの状態予測方法を構成する主要な手法について基本的な説明をしておく。
［ＪＩＴモデリング］
現在の挙動と近似した挙動が過去に観測されていたならば、現在の挙動が進展する様子は過去のものと近似したものになるであろうと考えることができる。この考え方を再現した予測手法の１つがJust-In-Time（ＪＩＴ）モデリングである。ＪＩＴモデリングは決まったモデルを持たない代わりに、過去のデータベクトルをそのままデータベースとして保持する。プロセスの予測が必要となったとき、過去データが蓄積されたデータベースから、要求点ベクトルと類似性の高いデータベクトルを検索し、局所モデルを構築して出力の推定を行う手法である。
【００１８】
対象とするプロセスが非線形かつ動的なプロセスであるとき、次式の回帰モデルでそのプロセスを表すことができる。
【００１９】
【数１】
【００２０】
ここで、プロセスの入力ベクトルｘｋと出力ベクトルｙｋを以下のように定義する。つまり、出力ベクトルｙｋは、ｋ時における入力ベクトルｘｋに対する（ｋ＋ｐ）時における出力、即ち予測値となる。
【００２１】
【数２】
【００２２】
時間の経過と共に、入力ベクトルｘｋと出力ベクトルｙｋのデータベクトルの組が（ｘ１，ｙ１），（ｘ２，ｙ２），…のように、対象とするプロセスから大量に得られ、データベクトル集合｛（ｘｋ，ｙｋ）｝（ｋ＝１，２，…）としてデータベースに蓄積される。ｋは離散化時間である。
【００２３】
予測したい時点における出力ベクトルｙｋｑに対応する入力ベクトルｘｋｑを要求点ベクトルとし、要求点ベクトルと類似性が高い近傍データベクトルを上記データベースから取得する。要求点ベクトルと類似性が高い近傍データベクトルを選択する際の指標としては、次式で示すようなベクトル間距離（ユークリッド距離）などを用いることができる。
【００２４】
【数３】
【００２５】
近傍データベクトル群｛（ｘｋｉ，ｙｋｉ）｝（ｉ＝１，２，…，ｍ）が取得されると、この近傍データベクトル群を用いて局所モデルの構築を行い、出力ベクトルｙｋｑの推定を行う。局所モデルとしては、重回帰モデルや、以下に示す相加平均法あるいは重み付き線形平均法などが用いられる。
【００２６】
【数４】
【００２７】
［ステップワイズ法］
ステップワイズ法は、目的変数に対する影響（寄与率）が小さい説明変数を除外するものであり、変数増加ステップと変数減少ステップとから構成される。以下、ステップワイズ法の手順について説明する。
【００２８】
（Ｂ１）目的変数に対する単寄与率が最大の説明変数を、先に選定した説明変数の中から決定する。具体的には、各説明変数ごとに単回帰モデルを作成して回帰係数を求めた後、（８）式によりＦ値を算出し、Ｆ値が最大となる説明変数を選択する。
【００２９】
【数５】
【００３０】
（Ｂ２）先の手順で決定された現モデルに説明変数を１つ追加することを考える。即ち、現モデルに含まれていない説明変数の中から、現モデルに一つ説明変数を追加した場合についてそれぞれ偏回帰係数を求め、（８）式によりＦ値を算出する。そして、Ｆ値が最大となる説明変数を探索する。
（Ｂ３）最大のＦ値が、前もって決めているＦｉｎ以上である場合は、その説明変数を現モデルに追加する。最大のＦ値がＦｉｎ未満の場合は、ステップワイズ法による選択手順を終了する。
【００３１】
（Ｂ４）新たに説明変数が追加された場合は、今度は逆に今までに取り込んだ説明変数が本当に有用な説明変数であるのかどうか調べる。即ち、現モデルを構成する説明変数の中で寄与率が一番低い説明変数を見つけるために、今までに取り込んだ説明変数を順番に一つずつ取り除いて当該説明変数が無い場合におけるＦ値を計算し、Ｆ値が最小となる説明変数を探索する。
（Ｂ５）最小のＦ値が、前もって決めているＦｏｕｔ（Ｆｉｎ≧Ｆｏｕｔ）未満である場合は、その説明変数を現モデルから削除する。説明変数を削除した場合は、（Ｂ４）、（Ｂ５）を繰り返し、さらに他の説明変数が削除できないか探索する。最小のＦ値がＦｏｕｔ以上の場合は、（Ｂ２）のステップに戻る。
【００３２】
［主成分分析］
主成分分析は、データの特徴抽出及び低次元化を目的とする多変量解析手法であり、変数間の相関関係を捉えるため、変数の線形結合によって得られる主成分と呼ばれる合成変数を使用する。主成分分析では、データを最も良く表現できる方向に第１主成分を設定し、第１主成分と直交する空間上で、第１主成分では表現できないデータの変動を最も良く表現できる方向に第２主成分を設定するという手順で、主成分を次々と設定していく。ここで、データを最も良く表現する方向というのは、主成分得点の分散が最大となる方向という意味である。また、主成分得点とは、主成分が張る部分空間へデータを射影した値である。
【００３３】
［Ｑ統計量］
Ｑ統計量は、データベクトルのうち、主成分によって張られる部分空間では表現できない部分を表す。Ｑ統計量は二乗予測誤差とも呼ばれ、以下のように定義されている。
Ｎ行×Ｍ列のデータ行列Ｘがあるものとする。ここで、Ｍは変数の数、Ｎはサンプル数であり、各変数は標準化されている。
データ行列Ｘを特異値分解すると次式のようになる。
【００３４】
【数６】
【００３５】
ＵとＶは直交行列であり、対角行列Ｓの対角要素には特異値ｓｒが降順に並んでいる。採用する主成分の数をＲとすると、第ｒ主成分は負荷量行列ＶＲの第ｒ列ｖｒで与えられる。
第ｒ主成分得点ｔｒは（１０）式で与えられ、第Ｒ主成分得点までをまとめて表現すると、（１１）式となる。
【００３６】
【数７】
【００３７】
ＴＲを元のＭ次元空間上の座標で表すと、再構築データ行列Ｘ＾は次のようになる。
【００３８】
【数８】
【００３９】
このとき、Ｑ統計量は次式で与えられる。
【００４０】
【数９】
【００４１】
［プロセスの状態予測方法の詳細手順］
続いて、図１のフローチャートに基づいて本実施の形態に係るプロセスの状態予測方法の手順について詳細に説明する。
（Ｃ１）プラント設備におけるプロセスの操業状態を示す観測データから構成される入力ベクトルｘｋ及び出力ベクトルｙｋのデータベクトルの組（ｘｋ，ｙｋ）（ｋ＝１，２，…）が蓄積された大規模データベース１０を作成する。
（Ｃ２）大規模データベース１０について、ステップワイズ法を用いて目的変数に対する寄与率が高い変数を選択して当該変数からなる新たなデータベース１１を作成する（ＳＴ１）。目的変数と説明変数の間に時間遅れが存在する可能性がある場合は、見込まれる最大の時間遅れ変数まで全て選択対象に加える。
作成されるデータベース１１（データセット）の構成を図２に示す。このデータベース１１では、入力変数の数がＭ個、出力変数の数がＬ個、各変数のサンプル数がＫ個とされている。各データは日時に応じたＩＤが付けられ、同じＩＤに属するデータは１つのデータベクトルとして扱われる。
【００４２】
（Ｃ３）予測したい時点における出力ベクトルＹｑに対応する入力ベクトルＸｑからなる要求点ベクトルを設定する（ＳＴ２）。図３に要求点ベクトルの構成を示す。
（Ｃ４）データベース１１に格納されている各データベクトルと要求点ベクトルとのベクトル間距離を（４）式や（５）式を用いて計算し、ベクトル間距離が小さいものから順にＮＮＭＡＸ個の近傍データベクトルを全て収集する。そして、収集した近傍データベクトルを、ベクトル間距離が近い順に近傍データセットＡとして保存する（ＳＴ３）。図４に近傍データセットＡの構成を示す。図４において「Ｎｏ．」が近傍データ数を表している。
（Ｃ５）近傍データセットＡの中から近傍データ数（Ｎｏ．）が１〜ＮＮＭＩＮまでの近傍データベクトルを選択して近傍データセットＢ０を作成する（ＳＴ４）。即ち、要求点ベクトルとのベクトル間距離が近いものから順にＮＮＭＩＮ個の近傍データベクトルを選択する。図５に近傍データセットＢ０の構成を示す。
【００４３】
（Ｃ６）近傍データセットＢ０に対して主成分分析を実施し、負荷量行列ＶＲを求める（ＳＴ５）。具体的には、近傍データセットＢ０をデータ行列Ｘとして特異値分解すればよい。
（Ｃ７）要求点ベクトルｘｑが（１４）式で表されるとすると、要求点ベクトルｘｑを再構築した再構築ベクトルｘ＾ｑは、負荷量行列ＶＲを用いて（１５）式により算出される。従って、近傍データセットＢ０に対するＱ統計量は、（１６）式より得ることができる（ＳＴ６）。算出されたＱ統計量は、図６に示すＱ値テーブルに保存される。
【００４４】
【数１０】
【００４５】
（Ｃ８）Ｑ統計量が算出された近傍データセットＢ０の近傍データ数がＮＮＭＡＸ以上であるかどうか判断される（ＳＴ７）。近傍データ数がＮＮＭＡＸ未満である場合は、近傍データセットＡの内、近傍データセットＢ０に含まれていない近傍データベクトルの中から、さらにＳ個の近傍データベクトルを、近傍データ数（Ｎｏ．）が小さいほうから（要求点ベクトルとのベクトル間距離が近いものから）選択し、近傍データセットＢ０に追加して新たな近傍データセットＢ１を作成する（ＳＴ９）。そして、ＳＴ５のステップに戻る。
（Ｃ９）一方、近傍データ数がＮＮＭＡＸ以上になった場合は、Ｑ値テーブルに基づいて、Ｑ統計量が最小となったデータセットＢｋをデータセットＡから選択する。そして、データセットＢｋに対応する出力ベクトルを、データセットＢｋのＩＤに基づいてデータベース１１から取得して、重回帰モデルや重み付き線形平均法などを用いて局所モデルを構築し、要求点ベクトルに対する出力の推定値を算出する（ＳＴ８）。
【００４６】
以上、本発明の一実施の形態について説明してきたが、本発明は何ら上記した実施の形態に記載の構成に限定されるものではなく、特許請求の範囲に記載されている事項の範囲内で考えられるその他の実施の形態や変形例も含むものである。例えば、上記実施の形態では、ＮＮＭＡＸ個の近傍データセット全てについてＱ統計量を求めているが、ある近傍データセットのＱ統計量が極小値であると判断された時点で、当該近傍データセットを選択してもよい。また、上記実施の形態では、変数を削減するためにステップワイズ法を使用しているが、変数が少ない場合はステップワイズ法を省略してもよい。
【実施例】
【００４７】
本実施の形態に係るプロセスの状態予測方法（以下では、「近傍データ数逐次可変型ＬＯＭ」と呼ぶ。）の効果について検証するため、ごみ処理プロセスにおけるガス化溶融炉の炉頂ガス温度の予測を行った。なお、従来の予測方法と比較するため、従来型ＬＯＭによるガス化溶融炉の炉頂ガス温度予測も併せて行った。
【００４８】
検証に使用したデータは、２年間に亘るごみ処理プロセスにおいて測定された観測データである。取り込んだデータは、ノイズ除去のため、１時間の移動平均フィルタを掛けて平滑化した。サンプリング時間は２０分、総データ数は３８８０９個である。
全ての変数に対してステップワイズ法を適用した場合、厖大な処理時間とコンピュータメモリが必要となると共に、不要な変数はモデリング精度低下の原因となる。そのため、予め炉頂ガス温度と関係があると考えられる説明変数を３７個に絞った後、各変数の遅れ時間を０から５０時間としてステップワイズ法による変数選択を実施した。その結果、２７個の説明変数を選択し、１時間後の炉頂ガス温度の予測を行った。
【００４９】
近傍データ数逐次可変型ＬＯＭにおける最大近傍データ数ＮＮＭＡＸは１２００個、最小近傍データ数ＮＮＭＩＮは１００個、近傍データ数の増加幅Ｓは１０個とした。また、主成分分析に使用した主成分数は１０個とした。
一方、従来型ＬＯＭにおける近傍データ数は、試行錯誤的に最も良い結果が得られる数値とした。
なお、近傍データ数逐次可変型ＬＯＭも従来型ＬＯＭも、局所モデルの構築には重回帰モデルを使用した。
【００５０】
要求点ベクトルを観測データからランダムに２００点選択して、各要求点ベクトルに対して１時間後の炉頂ガス温度の予測を行った。近傍データ数逐次可変型ＬＯＭによる予測値と実測値との相関度を図７に、近傍データ数を１２００個としたときの従来型ＬＯＭによる予測値と実測値との相関度を図８に示す。これらの図から、近傍データ数逐次可変型ＬＯＭは、従来型ＬＯＭに比べて予測値と実測値のバラツキが少なく、相関係数ｒが０．０８程度向上していることがわかる。
なお、図７以降において、目盛に表示されている「Ｔ」は基準温度を表している。即ち、基準温度Ｔからの変化量で表している。
【００５１】
次に、１時間ごとに１時間後の炉頂ガス温度を予測する連続予測を行った結果について説明する。近傍データ数逐次可変型ＬＯＭによる連続予測結果を示す時刻歴グラフを図９に、近傍データ数を１２００個としたときの従来型ＬＯＭによる連続予測結果を示す時刻歴グラフを図１０にそれぞれ示す。また、近傍データ数逐次可変型ＬＯＭによって連続予測した際に使用した近傍データ数のグラフを図１１に示す。
これらの図から、近傍データ数を１２００個とした場合には難しかった急激な温度上昇（時刻歴グラフの５時間目参照）が近傍データ数逐次可変型ＬＯＭでは捉えられていることがわかる。また、近傍データ数逐次可変型ＬＯＭでは、近傍データ数を１２００個とした場合に比べて全体的な誤差も小さくなっている。
【符号の説明】
【００５２】
１０：大規模データベース、１１：データベース
【特許請求の範囲】
【請求項１】
プラント設備におけるプロセスの操業状態を示す観測データから構成される入力ベクトル及び出力ベクトルが対となったデータベクトルが蓄積されたデータベースを作成し、予測したい時点における出力ベクトルに対応する入力ベクトルからなる要求点ベクトルに類似する前記データベクトルを近傍データベクトルとして前記データベースから少なくとも１つ以上取得し、前記近傍データベクトルから局所モデルを構築して、前記予測したい時点における出力ベクトルを推定するプロセスの状態予測方法において、
前記近傍データベクトルが格納された近傍データベクトル集合を該近傍データベクトルの数を変えて複数作成する工程と、前記複数の近傍データベクトル集合について主成分分析を実施して該各近傍データベクトル集合ごとに前記要求点ベクトルに対するＱ統計量を算出する工程と、前記Ｑ統計量が最小となる前記近傍データベクトル集合を選択して前記局所モデルを構築する工程とを備えることを特徴とするプロセスの状態予測方法
【請求項２】
請求項１記載のプロセスの状態予測方法において、前記近傍データベクトルの数が異なる前記複数の近傍データベクトル集合を作成する際、前記要求点ベクトルとのベクトル間距離が近い前記近傍データベクトルから順に格納して前記近傍データベクトル集合を作成することを特徴とするプロセスの状態予測方法。
【請求項１】
プラント設備におけるプロセスの操業状態を示す観測データから構成される入力ベクトル及び出力ベクトルが対となったデータベクトルが蓄積されたデータベースを作成し、予測したい時点における出力ベクトルに対応する入力ベクトルからなる要求点ベクトルに類似する前記データベクトルを近傍データベクトルとして前記データベースから少なくとも１つ以上取得し、前記近傍データベクトルから局所モデルを構築して、前記予測したい時点における出力ベクトルを推定するプロセスの状態予測方法において、
前記近傍データベクトルが格納された近傍データベクトル集合を該近傍データベクトルの数を変えて複数作成する工程と、前記複数の近傍データベクトル集合について主成分分析を実施して該各近傍データベクトル集合ごとに前記要求点ベクトルに対するＱ統計量を算出する工程と、前記Ｑ統計量が最小となる前記近傍データベクトル集合を選択して前記局所モデルを構築する工程とを備えることを特徴とするプロセスの状態予測方法
【請求項２】
請求項１記載のプロセスの状態予測方法において、前記近傍データベクトルの数が異なる前記複数の近傍データベクトル集合を作成する際、前記要求点ベクトルとのベクトル間距離が近い前記近傍データベクトルから順に格納して前記近傍データベクトル集合を作成することを特徴とするプロセスの状態予測方法。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【公開番号】特開２０１２−２４７８５５（Ｐ２０１２−２４７８５５Ａ）
【公開日】平成２４年１２月１３日（２０１２．１２．１３）
【国際特許分類】
【出願番号】特願２０１１−１１７１９８（Ｐ２０１１−１１７１９８）
【出願日】平成２３年５月２５日（２０１１．５．２５）
【出願人】（３０６０２２５１３）新日鉄エンジニアリング株式会社 (897)
【出願人】（８９９００００６８）学校法人早稲田大学 (602)
【Ｆターム（参考）】
【公開日】平成２４年１２月１３日（２０１２．１２．１３）
【国際特許分類】
【出願日】平成２３年５月２５日（２０１１．５．２５）
【出願人】（３０６０２２５１３）新日鉄エンジニアリング株式会社 (897)
【出願人】（８９９００００６８）学校法人早稲田大学 (602)
【Ｆターム（参考）】
[ Back to top ]