プラントの異常を診断する異常診断装置であって、所定のタイミングで、グループデータ記憶部から異なる複数の変数で構成される複数のグループに関するグループデータを読み出し、各グループについて、読み出したグループデータに含まれる変数の値と、プラントから新たに入力した当該プラントの状態を表す複数の変数の値とを利用してマハラノビス距離を算出する算出部と、各グループについて、前記算出部で算出されたマハラノビス距離と、グループ毎に予め設定される異常の判定に利用する判定値とを比較して異常を診断する診断部と、前記診断部における全てのグループについての診断結果を統合し、異常の原因を特定する特定部と、を備えることを特徴とする異常診断装置。所定のタイミングで、異常と特定されたグループの判定値を変更する変更部をさらに備えることを特徴とする請求項1記載の異常診断装置。上述したように、従来の異常診断では、高精度で異常診断することは困難であった。上記課題に鑑み、高精度で異常診断する異常診断装置を提供することを目的としている。上記目的を達成するために、請求項1記載の発明は、所定のタイミングで、グループデータ記憶部から異なる複数の変数で構成される複数のグループに関するグループデータを読み出し、各グループについて、読み出したグループデータに含まれる変数の値と、プラントから新たに入力した当該プラントの状態を表す複数の変数の値とを利用してマハラノビス距離を算出する算出部と、各グループについて、算出部で算出されたマハラノビス距離と、グループ毎に予め設定される異常の判定に利用する判定値とを比較して異常を診断する診断部と、診断部における全てのグループについての診断結果を統合し、異常の原因を特定する特定部とを備える。また、請求項2の発明は、所定のタイミングで、異常と特定されたグループの判定値を変更する変更部をさらに備える。本発明によれば、プラントの異常診断の精度を向上することができる。実施形態に係る異常診断装置は、プラントの運転の異常を診断する異常診断装置である。例えば、異常診断装置が診断するプラントは、発電プラントである。発電プラントは、複数の機器ポンプ、バルブ等を備えており、これらの機器を制御する値が目標値として設定されている。この目標値は、例えば、ポンプの圧力、バルブの開閉等である。また、発電プラントでは、複数のセンサを備えており、各センサで吸気量、排気量、発電量等が計測されている。さらに、発電プラントでは、様々な運転状態が設定可能であって、運転状態によって各機器に設定する目標値、各センサで計測される計測値が異なる。例えば、発電プラントでは、常に一定量の電力を発電する必要はなく、電力需要に合わせて発電する。したがって、時間帯、曜日、季節等に合わせて運転状態を変えることができ、運転状態で各機器の目標値が異なる。また、各機器の目標値が異なることにより、運転状態で各センサの計測値も異なる。図1に示すように、実施形態に係る異常診断装置1は、プラントから入力した新たな変数で蓄積データ21を更新する更新部11と、異常の診断に利用するグループを生成し、各グループに診断で利用する判定値を設定してグループデータ22として記憶させる生成部12と、プラントから新たに入力する変数とグループデータ22の変数を利用して複数のグループについてマハラノビス距離を算出する算出部13と、算出部13が算出したマハラノビス距離を利用してグループ毎に異常を診断する診断部14と、診断部14で得られた各グループに関する診断結果を利用して特定データ23から異常の原因を特定する特定部15と、グループデータ22の判定値を変更する変更部16とを備えている。異常診断装置1は、例えば、中央処理装置CPU10や記憶装置20を備える情報処理装置であって、記憶装置20に記憶される異常診断プログラムPが読み出されて実行されることで、図1に示すように、CPU10に更新部11、生成部12、算出部13、診断部14、特定部15及び変更部16が実装される。記憶装置20では、異常診断プログラムPの他、蓄積データ21、グループデータ22及び特定データ23を記憶している。また、異常診断装置1は、操作を入力するキーボード、マウス、操作ボタン、タッチパネル等の入力装置2と接続されており、異常診断の処理過程や結果を出力するディスプレイ、スピーカ等の出力装置3と接続されている。更新部11は、プラントから新たに各変数の値を入力すると、入力した各変数の値を含む新たなレコードを生成し、生成したレコードを追加して記憶装置20に記憶される蓄積データ21を更新する。蓄積データ21は、過去のプラントの状態を表す変数のレコードを蓄積したデータである。蓄積データ21の変数は、プラントの運転の条件値及び当該条件値の場合にプラントで計測された計測値である。図2に示す一例では、蓄積データ21は、プラントの運転の条件値であるプラントの各機器に設定する目標値変数1、2と、このモード及び目標値の場合に計測された計測値変数4、5等を関連づけたレコードを有している。また、図2に示す蓄積データ21の一例では、各レコードに時刻等のサンプル番号を付している。生成部12は、グループを生成する所定のタイミングで、記憶装置20から蓄積データ21を読み出し、変数が所定の抽出条件を満たすレコードを抽出して複数のグループを生成する。この抽出条件は、例えば、入力装置2を介して入力される。例えば、変数1や変数2目標値の値、変数4や変数5計測値の値のいずれか一つを抽出条件としても良いし、複数の変数の値の組み合わせを抽出条件としても良い。生成部12は、生成した複数のグループにプラントの異常診断に利用する判定値を設定し、グループの各変数の値と判定値とを含むグループデータ22を生成して記憶装置20に記憶する。すなわち、グループデータ22は、図3aや図3bに示すように、蓄積データ21の選択されたレコードのみを有するデータである。なお、図1の一例では、記憶装置20には1つのグループデータ22を記憶するものとして示しているが、記憶装置20では、複数のグループデータ22、すなわち、各グループについてのグループデータ22を記憶するものとする。ここで、生成部12がグループを生成する所定のタイミングは、例えば、定期的なタイミングや、入力装置2を介して入力されるタイミングである。また、生成部12で判定値とする値は、このグループについて求められるマハラノビス空間を表す値である。この判定値は、例えば、入力装置2を介して入力される。算出部13は、異常を診断する所定のタイミングで、記憶装置20から各グループのグループデータ22を読み出し、読み出したグループデータ22が有する各変数の値を抽出する。また、算出部13は、グループデータ22から抽出した各変数の値とプラントから入力した新たな変数とを利用して、各グループのマハラノビス距離を求め、求めた値を診断部14に出力する。ここで、異常を診断する所定のタイミングとは、例えば、定期的なタイミングである。具体的には、1まず、算出部13は、各変数に対して平均値と標準偏差を求める。2その後、算出部13は、各変数の値と、各変数に対して求めた平均値及び標準偏差を利用してデータを基準化し、各変数に対する基準化値を求める。3続いて、算出部13は、各変数に対して求めた基準化値を利用して、各変数についての相関行列を求めるとともに、相関行列の逆行列を求める。4最後に、算出部13は、求めた逆行列を利用してマハラノビス距離を求める。このようにして求められたマハラノビス距離により、対象のグループの値について現在の状態といつもの状態との違いを特定することができ、プラントにおける現在の状態といつもの状態との違いを求めることができる。診断部14は、算出部13から各グループのマハラノビス距離を入力すると、各グループのグループデータ22を読み出し、読み出したグループデータ22が有する判定値を抽出する。また、診断部14は、グループ毎に、算出部13から入力したマハラノビス距離を、グループデータ22から抽出した判定値と比較し、プラントの運転状態が正常であるか異常であるかを判定し、各グループのマハラノビス距離とともに診断結果として特定部15及び出力装置3に出力する。例えば、診断部14から診断結果を入力する出力装置3は、診断部14による診断結果を出力して異常の発生をユーザ等に通知することができる。例えば、図5に示すように、MT法を利用して異常を判定する場合、算出部13で求めたマハラノビス距離が設定された判定値より小さいとき、すなわちマハラノビス空間内にあるとき、診断部14は、現在のプラントの状態はいつもの状態と同じであるとし、プラントは正常に運転していると診断する。一方、算出部13で求めたマハラノビス距離が設定された判定値より大きいとき、すなわち、マハラノビス空間から外れているとき、診断部14は、現在のプラントの状態がいつもの状態とは異なる状態であるとし、プラントで異常が発生していると診断する。特定部15は、診断部14から診断結果を入力すると、特定データ23を読み出し、診断部14の各グループに関する診断結果を利用して、異常の原因を特定し、特定した異常の原因を変更部16及び出力装置3に出力する。例えば、特定部15から異常の原因を入力する出力装置3は、特定部15によって特定された異常の原因を出力して異常の原因を通知する。特定データ23は、各グループの診断結果と、異常原因のパターンとが関連づけられるデータである。ここで、例えば、特定データ23では、各グループのマハラノビス距離の変化パターンと異常原因とを関連付けていてもよい。すなわち、あるグループに関するマハラノビス距離が増加しているときには、プラントの所定箇所が異常原因であると特定することもできる。このようにマハラノビス距離の変化パターンから異常の原因を特定する場合、異常診断装置1では、記憶装置20等にマハラノビス距離の変化の履歴を記憶しておき、特定部15は、新たに入力する診断結果に加え、記憶される所定期間のマハラノビス距離の履歴を利用して異常の原因を特定する。変更部16は、異常診断に利用する判定値を更新するタイミングで、新たな判定値を求めて、グループデータ22の判定値を変更する。異常診断に利用する判定値を更新するタイミングとは、例えば、特定部15で異常の原因が特定されたタイミングである。または、変更部16は、定期的なタイミングで判定値を更新してもよい。例えば、変更部16は、過去の各状態でのマハラノビス距離の平均値を解析した統計値から対象のグループについての新たな判定値を求める。具体的には、正常時のマハラノビス距離の平均値と、異常時のマハラノビス距離の平均値とを比較することで、異常の判定に利用する判定値を求めることができる。なお、異常診断装置1は、複数の情報処理装置から構成されていてもよく、例えば、更新部11のみ他の処理部12~16とは異なる情報処理装置に含まれていてもよい。また、記憶装置20に記憶されるデータの一部のみ外部の記憶装置に記憶されていてもよい。続いて、図4に示すフローチャートを用いて、異常診断装置1において最適なグループで異常を診断する処理について説明する。生成部12が、グループデータ22を生成するタイミングで、記憶装置20から蓄積データ21を読み出して抽出条件を満たすレコードの値を抽出し、グループデータ22を生成して記憶装置20に記憶させるS1。その後、算出部13は、各グループについてそれぞれ、記憶装置20からグループデータ22を読み出し、プラントから入力した変数の値とグループデータ22に含まれる変数の値を利用して、マハラノビス距離を算出するS2。算出部13から複数のグループのマハラノビス距離を入力した診断部14は、各グループについて、入力したマハラノビス距離をグループデータ22において設定されている判定値と比較し、プラントの運転状態が異常であるか診断し、診断結果を出力装置3及び特定部15に出力するS3。特定部15は、各グループの診断結果を統合して、異常があるとき、特定データを利用して異常箇所を特定するS4。異常診断装置1では、判定値の変更のタイミングでない場合S5でNO、ステップS2乃至S4の処理を繰り返す。また、異常診断装置1では、判定値の変更のタイミングになるとS5でYES、変更部16で判定値を変更した後S6、ステップS2乃至S4の処理を繰り返す。なお、図4に示すフローチャートでは、グループデータ22の生成を繰り返すことについては示していないが、グループデータ22も、例えば定期的等の所定のタイミングで繰り返して生成しても良い。上述したように、本発明に係る異常診断装置では、グループ毎に設定した診断基準判定値でそれぞれ異常を診断したうえで、複数のグループの診断結果を統合してプラントの状態を診断するため、最適な判断基準により、高精度で異常診断することができる。また、異常診断装置では、異常診断をしながら変更した判断基準判定値を利用してプラントの異常を診断するため、高精度で異常診断することができる。以上、実施形態を用いて本発明を詳細に説明したが、本発明は本明細書中に説明した実施形態に限定されるものではない。本発明の範囲は、特許請求の範囲の記載及び特許請求の範囲の記載と均等の範囲により決定されるものである。