特許 7133154 運転異常検出システム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】

(24)【登録日】2022-08-31

(45)【発行日】2022-09-08

(54)【発明の名称】運転異常検出システム

(51)【国際特許分類】

   G05B 19/418 20060101AFI20220901BHJP

   G05B 23/02 20060101ALI20220901BHJP

【ＦＩ】

G05B19/418 Z

G05B23/02 302Z

【請求項の数】 6

(21)【出願番号】P 2019006970

(22)【出願日】2019-01-18

(65)【公開番号】P2020119009

(43)【公開日】2020-08-06

【審査請求日】2021-08-23

(73)【特許権者】

【識別番号】000183266

【氏名又は名称】住友大阪セメント株式会社

(74)【代理人】

【識別番号】110002734

【氏名又は名称】特許業務法人藤本パートナーズ

(72)【発明者】

【氏名】稲津 和喜

【審査官】石川 健一

(56)【参考文献】

【文献】特開２００２－３１８６１７（ＪＰ，Ａ）

【文献】特開２００８－１１７３８０（ＪＰ，Ａ）

【文献】特開２０１５－２１５７０９（ＪＰ，Ａ）

【文献】特開２０１４－０３５２８２（ＪＰ，Ａ）

【文献】特開２０１８－００１９９７（ＪＰ，Ａ）

【文献】特開２０１７－０９１０５６（ＪＰ，Ａ）

【文献】特開平０６－１５０１７８（ＪＰ，Ａ）

【文献】特開２０１８－１１６５４５（ＪＰ，Ａ）

(58)【調査した分野】(Int.Cl.，ＤＢ名)

Ｇ０５Ｂ １９／４１８

Ｇ０５Ｂ ２３／０２

(57)【特許請求の範囲】

【請求項1】

製造ラインにおける運転異常を検出する運転異常検出システムであって、
製造ラインにおける運転状態に対応する運転データを取得する運転データ取得手段と、
該運転データ取得手段により取得した運転データを所定の期間で分割し、分割した所定の期間の運転データから運転異常の判断に用いる指標データを該所定の期間毎に抽出する指標データ抽出手段と、
該指標データ抽出手段により抽出された指標データに対応した運転により製造された製造物の品質の分析を運転異常検出システム外で行った結果である既知の品質分析データを取得する品質分析データ取得手段と、
該品質分析データ取得手段で取得した品質分析データに基づいて前記指標データ抽出手段により所定の期間毎に抽出した指標データを正常時と異常時に仕分けする仕分手段と、
該仕分手段により仕分けられた指標データのうちの正常時及び異常時の指標データから正常時及び異常時のマハラノビス距離をそれぞれ算出する算出手段と、
正常時と異常時の指標データから前記算出手段が算出したマハラノビス距離の度数分布を求めて、該マハラノビス距離の度数分布から指標データの正常異常を識別するための識別度を求め、該識別度により異常検出に用いる最適な指標データを選択する指標データ選択手段と、
該指標データ選択手段で選択された最適な指標データのうちの正常時の指標データを用いて、正常運転時に対応する単位空間を構築し、構築された単位空間により算出されたマハラノビス距離に基づいて前記指標データ抽出手段により抽出される指標データに対応した運転が運転異常であるか否かを判断する運転異常判断手段と、
を備えていることを特徴とする運転異常検出システム。

【請求項2】

前記運転異常判断手段により運転異常であると判断された場合に、前記運転データに対応した項目であって、運転異常を引き起こす関連性の高い少なくとも一つの運転異常項目を、該運転異常項目の有無とこれに対応したマハラノビス距離の変化に基づいて抽出する運転異常項目抽出手段と、
前記運転異常項目が複数であった場合に、前記指標データを抽出する前記所定の期間の長さを変更し、変更した所定の期間の運転データから該指標データを前記指標データ抽出手段により再度抽出し、抽出した指標データから前記算出手段によりマハラノビス距離を再度算出し、構築した前記単位空間により算出されたマハラノビス距離に基づいて前記再度抽出された指標データに対応した運転が運転異常であるか否かを再度判断することにより前記運転異常項目抽出手段により抽出した運転異常項目の数を減らすように絞り込みを行う運転異常項目絞り込み手段と、
を備えていることを特徴とする請求項１に記載の運転異常検出システム。

【請求項3】

前記運転異常であると判断するための第１閾値並びに前記運転異常項目を抽出するための第２閾値を設定する閾値設定手段を更に備え、前記運転異常判断手段は、前記閾値設定手段により設定された第１閾値により前記指標データ抽出手段により抽出される指標データに対応した運転が運転異常であるか否かを判断し、前記運転異常判断手段により運転異常であると判断された場合に、運転異常を引き起こす関連性の高い少なくとも一つの運転異常項目を前記閾値設定手段により設定された第２閾値により抽出するように構成されていることを特徴とする請求項２に記載の運転異常検出システム。

【請求項4】

前記識別度は、前記正常時の指標データから求めたマハラノビス距離の度数分布における前記第１閾値未満の割合を第１割合（Ａ％）とし、前記異常時の指標データから求めたマハラノビス距離の度数分布における前記第１閾値以上の割合を第２割合（Ｂ％）としたとき、正常時の指標データが誤って異常であると判定される第１判定割合（１００－Ａ％）と異常時の指標データが誤って正常であると判定される第２判定割合（１００－Ｂ％）との和の１／２で算出される誤判定確率であり、前記第１閾値は、運転データ項目数で決定されるχ２（カイ二乗）分布における１％棄却域値又は予め設定される２．０の高い方の値が選択されることを特徴とする請求項３に記載の運転異常検出システム。

【請求項5】

前記指標データ選択手段は、前記指標データ抽出手段により抽出された指標データを、期間平均された第１指標データと、移動平均の傾きで表された第２指標データと、移動平均からの変動で表された第３指標データと、期間平均からの変動で表された第４指標データと、変動周波数成分に変換された第５指標データと、で構成し、第１指標データから第５指標データのそれぞれの識別度を前記誤判定確率で算出し、算出された５個の誤判定確率のうちの低い誤判定確率の指標データを単数又は複数選択することで前記最適な指標データを選択する手段であることを特徴とする請求項４に記載の運転異常検出システム。

【請求項6】

請求項１～５に記載の運転異常検出システムを備えたセメント製造装置。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、製造ラインにおける運転異常を検出する運転異常検出システムに関する。

【背景技術】

【0002】

従来から、オペレータが、製造ラインにおいて検出される代表的な運転データ（温度、圧力等の検出値）の挙動を監視操作画面を通して見ることで監視している。しかしながら、オペレータによる監視では、オペレータ毎の経験やスキルの差等により、異常が発生しているか否かの判断が異なってしまうという不都合がある。

【0003】

そこで、コンピュータを用いて、算出されたマハラノビス距離に基づいて製造工程の異常を数値化することにより判断するものが既に提案されている。具体的には、製造ラインの正常状態下における複数の製造管理パラメータのデータをサンプリングし、サンプリングした第１のサンプリングデータに基づいて複数の製造管理パラメータ毎のマハラノビス空間を作成し、製造ラインの稼働時における複数の製造管理パラメータのデータを所定の周期でサンプリングし、サンプリングした第２のサンプリングデータ及び前記マハラノビス空間からマハラノビス距離を算出し、算出したマハラノビス距離と所定値（予め設定されている閾値）とを比較することによって、製造ラインが異常であるか否かを判断するようにしている（例えば、特許文献１）。

【先行技術文献】

【特許文献】

【0004】

【文献】特開２００２－３１８６１７号公報

【発明の概要】

【発明が解決しようとする課題】

【0005】

特許文献１の構成では、マハラノビス距離を用いて、製造ラインの異常を判断する有効な手法の一つであるものの、製造ラインの正常状態下における複数の製造管理パラメータのデータをサンプリングした第１のサンプリングデータに基づいて、複数の製造管理パラメータ毎のマハラノビス空間を単に作成しているだけである。そのため、作成されたマハラノビス空間が、各製造管理パラメータに対応した有効なマハラノビス空間になっているとは限らない。そのため、製造ラインの異常を精度よく判断することができないことがあり、改善の余地があった。

【0006】

そこで本発明は、製造ラインの異常を精度よく検出することができる運転異常検出システムを提供することを目的とする。

【課題を解決するための手段】

【0007】

本発明の運転異常検出システムは、製造ラインにおける運転異常を検出する運転異常検出システムであって、製造ラインにおける運転状態に対応する運転データを取得する運転データ取得手段と、該運転データ取得手段により取得した運転データを所定の期間で分割し、分割した所定の期間の運転データから運転異常の判断に用いる指標データを該所定の期間毎に抽出する指標データ抽出手段と、該指標データ抽出手段により抽出された指標データに対応した運転により製造された製造物の品質の分析を運転異常検出システム外で行った結果である既知の品質分析データを取得する品質分析データ取得手段と、該品質分析データ取得手段で取得した品質分析データに基づいて前記指標データ抽出手段により所定の期間毎に抽出した指標データを正常時と異常時に仕分けする仕分手段と、該仕分手段により仕分けられた指標データのうちの正常時及び異常時の指標データから正常時及び異常時のマハラノビス距離をそれぞれ算出する算出手段と、正常時と異常時の指標データから前記算出手段が算出したマハラノビス距離の度数分布を求めて、該マハラノビス距離の度数分布から指標データの正常異常を識別するための識別度を求め、該識別度により異常検出に用いる最適な指標データを選択する指標データ選択手段と、該指標データ選択手段で選択された最適な指標データのうちの正常時の指標データを用いて、正常運転時に対応する単位空間を構築し、構築された単位空間により算出されたマハラノビス距離に基づいて前記指標データ抽出手段により抽出される指標データに対応した運転が運転異常であるか否かを判断する運転異常判断手段と、を備えていることを特徴としている。

【0008】

まず、運転データ取得手段により取得した運転データを所定の期間で複数に分割する。これら分割した複数の所定の期間の運転データから運転異常の判断に用いる指標データを所定の期間毎に抽出する。次に、品質分析データ取得手段で取得した品質分析データに基づいて前記所定の期間毎に抽出された指標データを正常時と異常時とに仕分ける。仕分けられた指標データのうちの正常時及び異常時の指標データから正常時及び異常時のマハラノビス距離をそれぞれ算出手段により算出する。指標データ選択手段により、正常時と異常時の指標データから算出したマハラノビス距離の度数分布を求める。求めたマハラノビス距離の度数分布から指標データの正常異常を識別するための識別度を求め、該識別度により異常検出に用いる最適な指標データを選択する。選択された最適な指標データのうちの正常時の指標データを用いて、正常運転時に対応する単位空間を構築し、構築された単位空間により算出されたマハラノビス距離に基づいて前記指標データ抽出手段により抽出される指標データに対応した運転が運転異常であるか否かを運転異常判断手段により判断する。このように、識別度により選択された最適な指標データのうちの正常時の指標データを用いて、正常運転時に対応する単位空間を構築することによって、その構築された単位空間により算出されたマハラノビス距離は、信頼性の高いものとなる。その信頼性の高いマハラノビス距離に基づいて抽出される指標データに対応した運転が運転異常であるか否かを精度よく判断することができる。

【0009】

また、本発明の運転異常検出システムは、前記運転異常判断手段により運転異常であると判断された場合に、前記運転データに対応した項目であって、運転異常を引き起こす関連性の高い少なくとも一つの運転異常項目を、該運転異常項目の有無とこれに対応したマハラノビス距離の変化に基づいて抽出する運転異常項目抽出手段と、前記運転異常項目が複数であった場合に、前記指標データを抽出する前記所定の期間の長さを変更し、変更した所定の期間の運転データから該指標データを前記指標データ抽出手段により再度抽出し、抽出した指標データから前記算出手段によりマハラノビス距離を再度算出し、構築した前記単位空間により算出されたマハラノビス距離に基づいて前記再度抽出された指標データに対応した運転が運転異常であるか否かを再度判断することにより前記運転異常項目抽出手段により抽出した運転異常項目の数を減らすように絞り込みを行う運転異常項目絞り込み手段と、を備えていてもよい。

【0010】

上記のように、運転異常判断手段により運転異常であると判断された場合に、運転異常を引き起こす関連性の高い少なくとも一つの運転異常項目を、該運転異常項目の有無とこれに対応したマハラノビス距離の変化に基づいて抽出し、運転異常項目が抽出された場合には、その運転異常項目による運転異常となる原因を解消することができる。また、前記運転異常項目が複数であった場合に、指標データ抽出手段による所定の期間の長さを変更して、マハラノビス距離を再度算出し、算出したマハラノビス距離に基づいて前記再度抽出された指標データに対応した運転が運転異常であるか否かを再度判断することにより、運転異常項目の数を減らすように運転異常項目数の絞り込みを行うことで、運転異常項目の運転異常となる原因を早期に見つけ出して運転異常となる原因を解消することができる。

【0011】

また、本発明の運転異常検出システムは、前記運転異常であると判断するための第１閾値並びに前記運転異常項目を抽出するための第２閾値を設定する閾値設定手段を更に備え、前記運転異常判断手段は、前記閾値設定手段により設定された第１閾値により前記指標データ抽出手段により抽出される指標データに対応した運転が運転異常であるか否かを判断し、前記運転異常判断手段により運転異常であると判断された場合に、運転異常を引き起こす関連性の高い少なくとも一つの運転異常項目を前記閾値設定手段により設定された第２閾値により抽出するように構成されていてもよい。

【0012】

上記のように、運転異常判断手段は、第１閾値により運転異常であるか否かを的確に判断し、かつ、第２閾値により運転異常項目を確実に抽出することができる。そして、運転異常判断手段により運転異常であると判断された場合に、運転異常を引き起こす関連性の高い少なくとも一つの運転異常項目を抽出するので、運転異常を早期に解消することができる。

【0013】

また、本発明の運転異常検出システムは、前記識別度は、前記正常時の指標データから求めたマハラノビス距離の度数分布における前記第１閾値未満の割合を第１割合（Ａ％）とし、前記異常時の指標データから求めたマハラノビス距離の度数分布における前記第１閾値以上の割合を第２割合（Ｂ％）としたとき、正常時の指標データが誤って異常であると判定される第１判定割合（１００－Ａ％）と異常時の指標データが誤って正常であると判定される第２判定割合（１００－Ｂ％）との和の１／２で算出される誤判定確率であり、前記第１閾値は、運転データ項目数で決定されるχ２（カイ二乗）分布における１％棄却域値又は予め設定される２．０の高い方の値が選択されてもよい。

【0014】

上記のように、指標データの正常異常を識別するための識別度を誤判定確率から算出することによって、運転が正常か異常かを判定する精度を更に高めることができる。

【0015】

また、本発明の運転異常検出システムは、前記指標データ選択手段は、前記指標データ抽出手段により抽出された指標データを、期間平均された第１指標データと、移動平均の傾きで表された第２指標データと、移動平均からの変動で表された第３指標データと、期間平均からの変動で表された第４指標データと、変動周波数成分に変換された第５指標データと、で構成し、第１指標データから第５指標データのそれぞれの識別度を前記誤判定確率で算出し、算出された５個の誤判定確率のうちの低い誤判定確率の指標データを単数又は複数選択することで前記最適な指標データを選択する手段であってもよい。

【0016】

上記のように、指標データ選択手段により算出された５個の誤判定確率のうちの誤判定する確率の低い誤判定確率の指標データを単数又は複数選択して、指標データを誤判定確率の低いものに絞り込むことによって、運転異常をより一層早期に解消することができる。

【0017】

また、本発明は、前記運転異常検出システムを備えたセメント製造装置であってもよい。

【発明の効果】

【0018】

本発明は、選択された最適な指標データのうちの正常時の指標データを用いて、正常運転時に対応する単位空間を構築し、構築された単位空間により算出されたマハラノビス距離に基づいて抽出される指標データに対応した運転が運転異常であるか否かを運転異常判断手段により判断することによって、製造ラインの異常を精度よく検出することができる運転異常検出システムを提供することができる。

【図面の簡単な説明】

【0019】

【図1】運転異常検出システムを構成するブロック図である。

【図2】同システムの構築を行うためのフローチャートである。

【図3】構築されたシステムに基づいて異常検出・分析を行うためのフローチャートである。

【図4】指標データの生データ、移動平均、期間平均を示すグラフである。

【図5】期間平均からの変動、移動平均からの変動を示すグラフである。

【発明を実施するための形態】

【0020】

以下、本発明の一実施形態に係る運転異常検出システムについて説明する。

【0021】

運転異常検出システムは、この実施形態では、石灰石などの粉粒状の原料を高温ガスで加熱することによって、製造物としての焼成物（セメントクリンカ）を製造するセメント製造ライン（以下において単に製造ラインという）における運転異常を検出するためのシステムとして使用する。つまり、この実施形態では、運転異常検出システムを備えたセメント製造装置について説明する。より具体的には、図１に示すように、運転異常検出システムに備えている制御装置Ｓに、運転データ取得手段１と、指標データ抽出手段２と、品質分析データ取得手段３と、仕分手段４と、算出手段５と、指標データ選択手段６と、運転異常判断手段７と、運転異常項目抽出手段８と、閾値設定手段９と、運転異常項目絞り込み手段１０と、を備えている。

【0022】

運転データ取得手段１は、製造ラインにおける運転状態に対応する運転データを取得する手段であり、取得した運転データを、例えば制御装置Ｓに備えているハードディスク（記憶手段）に記憶しておく。運転データとしては、例えば表４に示すように、キルン石炭ビン重量、キルン石炭焚量、キルン焼点温度、キルン電力、キルン排ガスＣＯ濃度や、表５に示すように、仮焼炉入口ガス温度、バーナ吹込ファン電流、最下段ＣｙＣＯ濃度、最下段Ｃｙ入口ガス温度、表４にもあるキルン電力等が挙げられる。ここでいう運転データとは、運転状態により変動する測定可能な数値からなるデータである。

【0023】

指標データ抽出手段２は、運転データ取得手段１により取得した運転データを所定の期間で分割し、分割した所定の期間の運転データから運転異常の判断に用いる指標データを所定の期間毎に抽出する手段であり、この場合も、抽出した指標データを、例えば制御装置Ｓに備えているハードディスク（記憶手段）に記憶しておく。指標データとしては、後述する期間平均、移動平均からの傾き、移動平均からの変動、期間平均からの変動、変動周波数成分等のデータである。指標データとは、運転データの数値変動に係るデータである。

【0024】

品質分析データ取得手段３は、指標データ抽出手段２により抽出された指標データに対応した運転により製造された製造物の品質の分析を運転異常検出システム外で行った結果である既知の品質分析データを前記システム外から取得する手段であり、この場合も、取得した既知の品質分析データを、例えば制御装置Ｓに備えているハードディスク（記憶手段）に記憶しておく。

【0025】

仕分手段４は、品質分析データ取得手段３で取得した品質分析データに基づいて指標データ抽出手段２により所定の期間毎に抽出された指標データを正常時と異常時に仕分けする手段である。

【0026】

算出手段５は、仕分手段４により仕分けられた指標データのうちの正常時及び異常時の指標データから正常時及び異常時のマハラノビス距離を算出する手段である。

【0027】

指標データ選択手段６は、算出手段５により、正常時と異常時の指標データから算出したマハラノビス距離の（に基づいて）度数分布を求めて、マハラノビス距離の度数分布から指標データの正常異常を識別するための識別度を求め、識別度により異常検出に用いる最適な指標データを選択する手段である。

【0028】

また、指標データ選択手段６は、指標データ抽出手段２により抽出された指標データを、期間平均された第１指標データと、移動平均の傾きにされた第２指標データと、移動平均からの変動で表された第３指標データと、期間平均からの変動で表された第４指標データと、変動周波数成分に変換された第５指標データと、に変換し、第１指標データから第５指標データのそれぞれの識別度を後述する誤判定確率で算出し、算出された５個の誤判定確率のうちの低い誤判定確率の指標データを単数又は複数（ここでは５個のうちの３個）選択することで前記最適な指標データを選択する手段である。

【0029】

期間平均とは、期間毎の生データ値を足し算して期間で割った値であり、図４のグラフに破線で示している。また、移動平均とは、時系列においてある一定間毎の平均値を区間をずらしながら求めた値であり、グラフ上に変化傾向を表すことができ、図４の細い実線で示している。また、移動平均からの変動とは、移動平均値と生データ値との差（乖離している量）を示したものであり、図５に示すグラフに太い実線で示している。また、期間平均からの変動とは、期間平均値と生データ値との差（乖離している量）を示したものであり、図５に示すグラフに細い実線で示している。

【0030】

前記識別度は、前記正常時の指標データから求めたマハラノビス距離の度数分布における前記第１閾値未満の第１割合（Ａ％）と、前記異常時の指標データから求めたマハラノビス距離の度数分布における前記第１閾値以上の第２割合（Ｂ％）としたとき、正常時の指標データが誤って異常であると判定される第１判定割合（１００－Ａ％）と異常時の指標データが誤って正常であると判定される第２判定割合（１００－Ｂ％）との和の１／２で算出される誤判定確率であり、前記第１閾値は、運転データ項目数で決定されるχ２（カイ二乗）分布における１％棄却域値又は予め設定される２．０の高い方の値が選択される。ここでは、運転データ項目数で決定されるχ２（カイ二乗）分布における１％棄却域値が２．０よりも低いため、第１閾値を２．０にしている。

【0031】

前記１％棄却域値について説明する。まず、データ項目数をｎ’とし、時系列データ数をnとし、マハラノビス距離をMDとすると、MD×ｎ’は、自由度（＝データ項目数）ｎ’のχ2分布に従う。次に、（１）～（５）の順に沿ってMDの算出を行う。
（１）データをD_ij(i=1～n、j=1～ｎ’)、データ項目jの平均をU_j、標準偏差をσ_jとする。;D_ij：行列
（２）D_ijをU_j,σ_jで規格化すると、D'_ij=(D_ij-U_j)/σ_jとなる。
（３）D'_ijの相関行列を求めると、A_kl=Σ_i(D'_ik×D'_il)/n ; i:1～n k,l:1～ｎ’となる。
Σ_i全てのiについて和をとる
（４）A_klの逆行列を求めると、B_kl=(A_kl)^-1となる。
（５）時系列i番目のマハラノビス距離は、MD_i=Σ_k,l (B_kl×D'_ik×D'_il)/ｎ’ ; i:1～n k,l:1～ｎ’
Σ_k,l全てのk,lの組み合わせにおいて和をとる。
そして、例えば、MD１%棄却域を算出すると、
データ項目数ｎ’=10の場合を例にすると、
χ2分布表（表７）より1%棄却域の値は、自由度n’=10、確率p=0.010を読み取り、23.2093であり、
MD×ｎ’=23.2093より、MD=23.2093/10=2.32093がMDの1%棄却域となる。
また、例えばデータ項目数ｎ’=40の場合では、
χ2分布表（表７参照）より1%棄却域の値は、自由度n’=40、確率p=0.010を読み取り、63.6907であり、
MD×ｎ’=63.6907より、MD=63.6907/40=1.59227がMDの1%棄却域となる。
この実施形態では、データ項目数ｎ’=30としており、上記より求まる1%棄却域は、50.8922/30=1.69641となり、1.69641＜2.0のため、MD閾値は2.0に設定している。なお、表７は、必要な部分のみの数値を取り出した表としている。つまり、確率ｐの０．９９０～０．０５０の間の確率についての自由度ｎ’の値を省き、また、自由度２～９、１１～２９、３１～３９、４１以降を省いた表にしている。

【0032】

【表1】

【0033】

表３に、表１及び表２の指標データの度数分布表の値から演算された識別度を示している。ここで、期間平均の正常時の第１割合（Ａ％）は、１１０／（１１０＋２３）×１００＝８３％である。なお、１１０は、第１閾値の２．０未満の正常時の指標データの合計数であり、２３は、第１閾値の２．０以上の正常時の指標データの合計数である。また、期間平均の異常時の第２割合（Ｂ％）は、５７／（５７＋１９）×１００＝７５％である。なお、５７は、第１閾値の２．０以上の異常時の指標データの合計数であり、１９は、第１閾値の２．０未満の異常時の指標データの合計数である。他の移動平均の傾き、移動平均からの変動、期間平均からの変動、変動周波数成分も同様に算出できる。これらの算出結果に基づいて識別度が前述した式により算出される。算出された結果が、表３に示すように、期間平均された第１指標データにおける識別度が２１％、移動平均の傾きにされた第２指標データにおける識別度が１４％、移動平均からの変動で表された第３指標データにおける識別度が１７％、期間平均からの変動で表された第４指標データにおける識別度が３７％、変動周波数成分に変換された第５指標データにおける識別度が４３％となっている。第４指標データにおける識別度及び第５指標データにおける識別度が悪いため、除外し、第１指標データと第２指標データと第３指標データとを用いて運転データの異常検出及び分析を行った。

【0034】

【表2】

【0035】

【表3】

【0036】

【表4】

【0037】

運転異常判断手段７は、前記指標データ選択手段６により選択された最適な指標データのうちの正常時の指標データを用いて、正常運転時に対応する単位空間を構築する。単位空間の構築とは、前記指標データ選択手段６により選択された最適な指標データの個数をα個とした場合、段落番号［００３０］の手順におけるデータ項目数ｎ’をｎ’×αに読み替えた上で同様な手順に従い、マハラノビス距離の算出式（［００３０］（５）の式）を導出することを意味する。このようにして構築された単位空間により算出されるマハラノビス距離に基づいて前記指標データ抽出手段２により抽出される指標データに対応した運転（この運転は異常検出・分析を開始して実際に運転しているリアルタイムの運転である）が運転異常であるか否かを判断する手段である。

【0038】

運転異常項目抽出手段８は、後述する閾値設定手段により設定された第１閾値により前記指標データ抽出手段により抽出される指標データに対応した運転が運転異常であるか否かを判断し、前記運転異常判断手段７により運転異常であると判断された場合に、運転異常を引き起こす関連性の高い少なくとも一つの運転異常項目を、該運転異常項目の有無とこれに対応したマハラノビス距離の変化に基づいて後述する閾値設定手段９により設定された第２閾値により抽出する手段である。運転異常項目とは、前記運転異常判断手段７によって運転異常と判断された指標データの属する運転データに対応した項目である。

【0039】

閾値設定手段９は、運転異常であると判断するための第１閾値並びに前記運転異常項目を抽出するための第２閾値を設定する手段である。この第２閾値は、０を越える値でマハラノビス距離の平均値である１．０を越えない値、つまり、０．１～１．０の間の任意の値に設定することが好ましく、より好ましい範囲が０．５～１．０であり、この実施形態では、０．５に設定している。

【0040】

運転異常項目絞り込み手段１０は、前記運転異常項目が複数であった場合に、前記指標データを抽出する前記所定の期間の長さを変更（図３では短縮であるが、長くしてもよい）し、変更した所定の期間の運転データから該指標データを指標データ抽出手段２により再度抽出し、抽出した指標データから算出手段５によりマハラノビス距離を再度算出し、構築した前記単位空間により算出されたマハラノビス距離に基づいて前記再度抽出された指標データに対応した運転が運転異常であるか否かを再度判断することにより前記運転異常項目抽出手段により抽出した運転異常項目の数を減らすように絞り込みを行う手段である。

【0041】

次に、運転データの異常検出・分析を行うためのマハラノビス距離の算出式を構築することを、図２のシステム構築開始のフローチャートに基づいて説明する。

【0042】

まず、既存の運転データと品質分析データとを取得し、前述したように記憶させる（ステップＳ１）。取得した運転データを所定の期間で複数分割する（ステップＳ２）。分割した複数の所定の期間のうちの１つの所定の期間の運転データから指標データｉを所定の期間毎に抽出する（ステップＳ３）。これらステップＳ２とステップＳ３とが、前記制御装置Ｓに備えている指標データ抽出手段２により実行される。抽出した指標データｉを品質分析データ取得手段３で取得した品質分析データに基づいて正常時と異常時に仕分けする（ステップＳ４）。このステップＳ４が、前記制御装置Ｓに備えている仕分手段４により実行される。前記算出手段５により仕分けられた指標データのうちの正常時及び異常時の指標データから正常時及び異常時のマハラノビス距離をそれぞれ算出し、算出したマハラノビスの度数分布を求め、その度数分布から指標データの正常異常を識別するための識別度を求める（ステップＳ５）。指標データｉの全候補（１～ｎ個）を網羅したかどうかを判定する（ステップＳ６）。網羅していなければ、次の指標データに移ってステップＳ３に戻り、ステップＳ４で正常時と異常時とに仕分けて、ステップＳ５で正常異常を識別するための２番目の識別度ｉを求める。指標データｉの全候補を網羅すると、ｉ＝ｎとなり、ステップＳ７に移行する。ステップＳ７では、求めたｎ個の識別度（ｎ＞１）のうちの上位α個の指標データを選択する。つまり、前述したように、前記指標データ選択手段６により最適な識別度を選択する。選択した指標データのうちの正常時の指標データを用いて、単位空間とマハラノビス距離の算出式を構築して（ステップＳ８）、システム構築を終了する。

【0043】

システム構築が終了すると、製造ラインの異常検出・分析を行うことになり、図３に示す異常検出・分析開始のフローチャートに基づいて説明する。

【0044】

まず、焼成物（セメントクリンカ）を製造する製造ラインの運転を開始し、新規の所定の期間の運転データを前記運転データ取得手段１により取得する（ステップＳ１１）。取得した所定の期間の運転データから前記指標データ抽出手段２により指標データ（ここでは、期間平均、移動平均からの傾き、移動平均からの変動）を抽出する（ステップＳ１２）。抽出した指標データより前記構築したマハラノビス距離の算出式を用いてマハラノビス距離（ＭＤ）を算出する（ステップＳ１３）。算出したマハラノビス距離（ＭＤ）が前記第１閾値（図３ではＭ１）よりも大きいか小さいかを判定する（ステップＳ１４）。マハラノビス距離（ＭＤ）が第１閾値（Ｍ１）よりも大きいと判定されると、運転異常であると判断する。表４では、Ｎｏ３－６の４つのマハラノビス距離（ＭＤ）が第１閾値（Ｍ１＝２．０）よりも大きいため異常であると判断する。これら異常における運転異常項目を抽出するためにステップＳ１５に移行する。ステップＳ１５では、所定の期間の運転データの指標データの運転項目について、運転項目を含めた場合と含めない場合のマハラノビス距離（ＭＤ）の差（ΔＭＤ）を算出する。この算出した差（ΔＭＤ）が前記第２閾値（図３ではＭ２であり、前記のように０．５と定めた）よりも大きいか小さいかを判定する（ステップＳ１６）。差（ΔＭＤ）が前記第２閾値（Ｍ２）よりも大きいと判定されると、要監視項目（運転異常項目）として記録される（ステップＳ１７）。差（ΔＭＤ）が前記第２閾値（Ｍ２）よりも小さいと判定されると、要監視項目（運転異常項目）が存在しないと判断し、ステップ１８に移行する。なお、ステップＳ１４でＭＤがＭ１よりも小さいと判定されると、ステップＳ１９へ移行する。また、表４に、算出した前記差（ΔＭＤ）の値を示している。この表４では、キルン石炭ビン重量のＮｏ３～Ｎｏ６の４か所において差（ΔＭＤ）が第２閾値（Ｍ２＝０．５）よりも大きい値になっており、キルン石炭ビン重量の変動が異常発生の発端となり、この異常によってキルン焼点温度（Ｎｏ５、Ｎｏ６において差（ΔＭＤ）が第２閾値（Ｍ２＝０．５）よりも大きい値になっている）及びキルン排ガスＣＯ濃度（Ｎｏ６において差（ΔＭＤ）が第２閾値（Ｍ２＝０．５）よりも大きい値になっている）に波及しているものと考えられ、これら３個がステップＳ１７で要監視項目（運転異常項目）として記録される。

【0045】

【表5】

【0046】

前記３個の要監視項目（運転異常項目）を記録した後は、全項目の分析が終了したか否かを確認し（ステップＳ１８）、終了していない場合には、ステップＳ１５へ戻り、終了している場合には、全所定期間の分析が終了しているか否かを確認し（ステップＳ１９）、分析が終了しない場合には、ステップＳ１１に戻り、分析が終了している場合には、同時間帯で要監視項目（運転異常項目）が複数あるか否かを確認し（ステップＳ２０）、要監視項目（運転異常項目）が１個の場合には、異常検出・分析を終了する。また、同時間帯で要監視項目（運転異常項目）が複数ある場合には、運転データの取得期間を短縮し（ステップＳ２１）、ステップＳ１１へ戻り、再度、運転ラインの異常検出・分析を行う。

【0047】

同時間帯で要監視項目（運転異常項目）が複数ある場合を、例えば表５（表４とは異なる項目も含んだものである）に５個の項目におけるマハラノビス距離の差を示している。この場合、指標データの抽出期間が１時間になっている。表５では、Ｎｏ２－６の５つのマハラノビス距離（ＭＤ）が第１閾値（Ｍ１＝２．０）よりも大きいため異常であると判断する。Ｎｏ２で仮焼炉入口ガス温度（移動平均からの変動）、バーナ吹込ファン電流（移動平均からの変動）、最下段CyCO濃度（移動平均傾き）の３項目が同期間において異常項目と判定されているため、異常原因の項目を絞ることができない。そこで、指標データの抽出期間を短縮した１５分にしてマハラノビス距離の差を演算したものを表６に示している。表６では、Ｎｏ５－７の３つのマハラノビス距離（ＭＤ）が第１閾値（Ｍ１＝２．０）よりも大きいため異常であると判断する。そして、Ｎｏ５においてバーナ吹込ファン電流の移動平均からの変動のみが差（ΔＭＤ）が第２閾値（Ｍ２＝０．５）異常となっている。よって、バーナ吹込ファン電流の変動が異常発生の発端となり、この異常によって仮焼炉入口ガス温度の移動平均からの変動において、Ｎｏ７において差（ΔＭＤ）が第２閾値（Ｍ２＝０．５）よりも大きい値になっており、仮焼炉入口ガス温度の変動に波及しているものと考えられる。このように異常項目がバーナ吹込ファン電流であると特定でき、図３に示すステップＳ２０で同時間帯で要監視項目（運転異常項目）が１つと判断され、異常検出・分析を終了する。

【0048】

【表6】

【0049】

【表7】

【0050】

表４のNo３のキルン石炭ビン重量の３つのΔMDを計算する計算方法について（１）～（８）の手順に沿って説明する。
（１）規格化した正常時データ（表８－１～表８－３のD'_ij）から相関行列の逆行列（表９のB_ij）を求める。
（２）No３の規格化したデータ（表１０のD'3_i）を求める。
（３）No３のマハラノビス距離を求める（この例では項目数15のため自由度15である）。
MD₃＝(Σ_ij D'3_i×B_ij×D'3_j)/15＝2.14となる。
（４）規格化した正常時データより１列目（キルン石炭ビン重量期間平均）のみ削除したもの（表１１－１～表１１－３）、その相関行列の逆行列（表１２）、No３の規格化したデータより１列目のみ削除したもの（表１３）を作成する。
（５）（３）と同手順で1列目のみ削除した場合のNo３のマハラノビス距離を求める。
（１列削除したので自由度は１４となる）
MD_3-1 = 2.16
（６）（４），（５）と同手順で2列目（キルン石炭ビン重量移動平均傾き）のみ削除した場合のNo.3のマハラノビス距離を求める。
MD_3-2 = 1.90
（７）（４），（５）と同手順で3列目（キルン石炭ビン重量移動平均からの変動）のみ削除した場合のNo.3のマハラノビス距離を求める。
MD_3-3 = 1.43
（８）元のMDと一列削除した場合のMDとの差を求める。このＭＤの差がΔMDとなる。
MD₃-MD_3-1 ＝2.14 - 2.16 ＝-0.02→小数1位に丸めて 0.0
MD₃-MD_3-2 ＝2.14 - 1.90 ＝0.24→小数1位に丸めて 0.2
MD₃-MD_3-3 ＝2.14 - 1.43 ＝0.71→小数1位に丸めて 0.7
上記に示す０．０が、Ｎｏ３のキルン石炭ビン重量期間平均のΔＭＤであり、０．２が、Ｎｏ３のキルン石炭ビン重量移動平均傾きのΔＭＤ、０．７が、Ｎｏ３のキルン石炭ビン重量移動平均からの変動のΔＭＤである。

【0051】

【表8-1】

【0052】

【表8-2】

【0053】

【表8-3】

【0054】

【表9】

【0055】

【表10】

【0056】

【表11-1】

【0057】

【表11-2】

【0058】

【表11-3】

【0059】

【表12】

【0060】

【表13】

【0061】

以上、本発明につき一実施形態を取り上げて説明したが、本発明は前記実施形態に限定されない。

【0062】

前記実施形態では、石灰石などの粉粒状の原料を高温ガスで加熱することによって、焼成物（セメントクリンカ）を製造する製造ラインにおける運転異常を検出するためのシステムを示したが、抵抗やコンデンサ等の回路部品及び半導体等の電子部品並びに各種の機械を製造する製造ラインにおける運転異常を検出するためのシステムにも本発明は適用できる。

【0063】

また、前記実施形態では、運転異常項目絞り込み手段により一回のみ運転異常項目の絞り込みを行ったが、一回で運転異常項目を良好に絞り込めなかった場合には、良好に絞り込めるように複数回絞り込むようにしてもよい。

【符号の説明】

【0064】

１…運転データ取得手段、２…指標データ抽出手段、３…品質分析データ取得手段、４…仕分手段、５…算出手段、６…指標データ選択手段、７…運転異常判断手段、８…運転異常項目抽出手段、９…閾値設定手段、１０…運転異常項目絞り込み手段、Ｓ…制御装置

【図1】

【図2】

【図3】

【図4】

【図5】