特許 6148229 遷移信号の動的クラスタリング

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

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

(11)【特許番号】6148229

(24)【登録日】2017年5月26日

(45)【発行日】2017年6月14日

(54)【発明の名称】遷移信号の動的クラスタリング

(51)【国際特許分類】

   G06F 17/30 20060101AFI20170607BHJP

【ＦＩ】

   G06F17/30 210D

【請求項の数】7

【全頁数】12

(21)【出願番号】特願2014-513017(P2014-513017)

(86)(22)【出願日】2012年5月25日

(65)【公表番号】特表2014-517398(P2014-517398A)

(43)【公表日】2014年7月17日

(86)【国際出願番号】CA2012050343

(87)【国際公開番号】WO2012162825

(87)【国際公開日】20121206

【審査請求日】2015年4月3日

(31)【優先権主張番号】2,741,202

(32)【優先日】2011年5月27日

(33)【優先権主張国】CA

【前置審査】

(73)【特許権者】

【識別番号】507381363

【氏名又は名称】ハイドロ−ケベック

(74)【代理人】

【識別番号】100158920

【弁理士】

【氏名又は名称】上野 英樹

(72)【発明者】

【氏名】レオナール， フランソワ

【審査官】 齊藤 貴孝

(56)【参考文献】

【文献】 特開平０７−００６１４９（ＪＰ，Ａ）

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

【文献】 特表２００８−５３９０３８（ＪＰ，Ａ）

【文献】 米国特許第０５４７９５７０（ＵＳ，Ａ）

【文献】 米国特許出願公開第２０１１／０１９１０７６（ＵＳ，Ａ１）

【文献】 山本 貴大、外２名，クラスタリングを用いた直交変調信号の同定，電子情報通信学会技術研究報告，日本，社団法人電子情報通信学会，２００８年 ５月２２日，第１０８巻，第６２号，ｐ．６５−６９

【文献】 殿村 正延、外１名，多層パーセプトロンにおける内部情報最適化アルゴリズムと汎化能力の解析，電子情報通信学会論文誌，日本，社団法人電子情報通信学会，２００１年 ５月 １日，第Ｊ８４−Ｄ−ＩＩ巻，第５号，ｐ．８３０−８４２

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

Ｇ０６Ｆ １７／３０

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

【請求項1】

遷移信号をクラスタリングするための方法であって、
獲得手段によって前記遷移信号を獲得するステップと、
獲得された新たな遷移信号がそれぞれ、獲得された時点で前に獲得された同様の遷移信号と共にメモリに記憶されるいずれかのクラスタに入るように、前記メモリに記憶された比較およびクラスタリング規則に基づいて、処理装置により、それぞれの前記獲得された新たな前記遷移信号を前に獲得された遷移信号によって定義されるクラスタシグネチャと比較することによって、同様の遷移信号からなるクラスタを超空間内で動的に生成するステップと、
前記メモリに記憶された前記クラスタを前記処理装置によって解析して、前記メモリに記憶された前記クラスタに集められた前記遷移信号によって定義されるそれぞれのシグネチャを決定するステップと、
前記処理装置によって前記シグネチャを処理して、前記遷移信号の固有の属性に関連付けることができる現象を検出するステップと
を含み、前記比較およびクラスタリング規則は計量のモデリングを用い、各クラスタは前記クラスタに集められた前記遷移信号における雑音の半径特性を有する超球に近い、方法。

【請求項2】

前記処理装置による前記クラスタの生成が、前記獲得された遷移信号を再整合するステップと、遷移と遷移、遷移とシグネチャ、およびシグネチャとシグネチャの距離の比較の少なくとも１つを行うステップとを含む請求項１に記載の方法。

【請求項3】

前記クラスタの生成が、獲得されたいくつかの初期遷移信号から生成されるクラスタから始まる請求項１に記載の方法。

【請求項4】

前記処理装置による前記クラスタの生成が、同様の遷移信号および同様のクラスタをマージするステップと、獲得された新たな遷移信号のクラスタリングに適用可能な前記比較およびクラスタ規則の新たなパラメータを計算するステップとを含む請求項１に記載の方法。

【請求項5】

前記パラメータが、距離計算計量を備える請求項４に記載の方法。

【請求項6】

請求項１〜５のいずれか一項に記載の方法のステップを行うための、処理装置と、前記処理装置によって実行される実行可能命令を記憶するメモリとを有するコンピュータシステム。

【請求項7】

請求項１〜５のいずれか一項に記載の方法のステップを行うためのコンピュータシステムによって実行される実行可能命令を記憶する非一時的な有形のコンピュータ可読記憶媒体。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、信号解析方法に関し、より詳細には、遷移信号の動的クラスタリングのためのシステムおよび方法に関する。

【背景技術】

【0002】

遷移信号は、レーダー、部分放電、アークノイズ（プラズマでの不定の放電）、株価変動、流体キャビテーション、音響放出、地球波（ｔｅｌｌｕｒｉｃ ｗａｖｅｓ）、および画像（ｉｍａｇｅｒｙ）など様々な分野で見られる。

【0003】

これらの分野の多くにおける継続的な問題は、各遷移の別個の処理が、多大な計算労力を必要とし、雑音の多い信号を対象とすることである。

【0004】

米国特許第６，０８８，６５８号明細書（Ｙａｚｉｃｉ他）、米国特許第６，８６８，３６５号明細書（Ｂａｌａｎ他）、米国特許第７，５７９，８４３号明細書（Ｙｏｕｎｓｉ他）、および米国特許出願第２００８／００８８３１４号明細書（Ｙｏｕｎｓｉ他）が、信号を解析するための従来技術のシステムおよび方法の例を提供しているが、これらは、時間がかかり、大きなリソースを消費する計算およびコンピュータ関連タスクを含む。

【発明の概要】

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

【0005】

本発明の目的は、遷移信号をクラスタリングするシステムおよび方法であって、同様の遷移を１つの特性シグネチャに再グループ化して、各グループに関してより少数のシグネチャ、例えば１グループ当たり１つのシグネチャを処理するシステムおよび方法を提供することである。シグネチャがより少ないので、計算時間が短縮される。個々の成分よりも雑音の少ないシグネチャにより、処理結果がはるかに正確になり、また結果が予め分類されている。

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

【0006】

本発明の一態様によれば、遷移信号をクラスタリングするための方法であって、
遷移信号が来たときに遷移信号を獲得するステップと、
獲得された新たな遷移信号がそれぞれ、前に獲得された同様の遷移信号と共に１つのクラスタに入るように、比較およびクラスタリング規則に基づいて、同様の遷移信号からなるクラスタを超空間内で動的に生成するステップと、
クラスタを解析して、クラスタに集められた遷移信号によって定義されるそれぞれのシグネチャを決定するステップと、
シグネチャを処理して、遷移信号の固有の属性に関連付けることができる現象を検出するステップと
を含む方法が提供される。

【0007】

本発明の別の態様によれば、処理装置と、上記の方法のステップを行うために処理装置によって実行すべき実行可能命令を記憶するメモリとを有するコンピュータシステムも提供される。

【0008】

本発明の別の態様によれば、上記の方法のステップを行うためにコンピュータシステムが実行すべき実行可能命令を記憶する非一時的な有形のコンピュータ可読記憶媒体が提供される。

【0009】

以下、本明細書において、以下の図面を参照しながら好ましい実施形態の詳細な説明を述べる。

【図面の簡単な説明】

【0010】

【図1】本発明による方法を示す流れ図である。

【図2】２つのクラスタを示す概略図である。

【図3】測定値のクラスタ部分放電シグネチャ（灰色のライン）と、１つの測定値スナップショット（細い黒色のライン）とを示すグラフである。

【図4】解像力を示すヒストグラムである。

【図5】クラスタコヒーレンシーを示すヒストグラムである。

【図6】シーケンス当たりのクラスタの数を示すヒストグラムである。

【図7】ｃ_distパラメータに伴う処理時間の減少を示すグラフである。

【図8】重畳されたクラスタの位相分解された部分放電図である。

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

【0011】

本開示において、時間または空間は、任意の他の性質を有する任意の他の次元に置き換えることができる。

【0012】

本発明は、母集団の一部に関して反復的な特性を有する遷移信号に関する。用語「反復」は、その時間または空間内で同じ遷移を複数回観察することができることを意味するものと理解されたい。これらの遷移は、変化することがある振幅を有し、また、雑音、測定誤差、遷移キャリアの時間的もしくは空間的な歪み、または任意の他のモデル化可能な現象によって（デジタルで、解析的に、または統計的に）説明することができる低い非類似特性も有する。

【0013】

本発明は、遷移が測定システムによって捕捉されるか、または適切なデバイスによって獲得されるときに、遷移を動的に再グループ化することを提案する。用語「再グループ化」によって、同様の遷移が同じグループに集められることを意味するものと理解されたい。したがって、所与のグループが少なくとも１つの遷移を含み、その結果、複数の遷移が関連付けられる少なくとも１つのグループが生じる。用語「非動的なグループ化」は、得られるすべての遷移との比較が行われることを意味する。これは、より最適に近いが、この手法の計算時間は膨大である。

【0014】

空間

【数1】

（ここで、Ｎは、遷移を特徴付ける時間点または空間点の数）内で表すと、グループ化は、この超空間内に点の雲のように現れる。グループ化の遷移の平均に対応するグループ化の質量中心をシグネチャと呼ぶ。

【0015】

遷移と遷移、遷移とシグネチャ、またはシグネチャとシグネチャの比較は、比較する両要素間の相関を最大にするため、または距離を最小にするために、時間または空間シフトを必要とする。距離ベースの比較基準の場合、このシフトは、

【数2】

が遷移ＸとシグネチャＳとの間の距離を表すように、ゼロ次で、１ブロックで実現される。遷移（またはシグネチャ）を補間して伸縮させることによって、１次シフトを実現することが可能である。２次についても同じである。さらに、動的時間伸縮型の方法を、距離計算のためのシフティング手段とみなすことができる。基本的には、比較する両要素間の適切なシフトは、比較時に行われる。

【0016】

図１を参照すると、本発明による方法は以下のように進む。ブロック２によって示されるように、遷移が捕捉され、ブロック４によって示されるように、遷移は、クラスタシグネチャを有する動的クラスタ構成に関してクラスタリングを受ける。ブロック６によって示されるように、シグネチャが解析され、ブロック８によって示される精査中の遷移の特性に関係付けられるさらなる処理に有用なシグネチャ特性を決定する。この方法は、処理装置と、上述したステップを行うために処理装置によって実行すべき実行可能命令を記憶するメモリとを有するコンピュータシステムで実装することができる。また、この方法は、この方法のステップを行うためにコンピュータシステムが実行すべき実行可能命令を記憶する非一時的な有形のコンピュータ可読記憶媒体の形態を取ることもできる。

【0017】

以下に、部分放電検出、位置特定、および解析に関連した本発明の一実施形態の一例を与える。本発明は、そのような実施形態および適用例に限定されず、本発明から逸脱することなく変更および修正を施すことができることを理解されたい。

【0018】

本発明による方法は、動的時間クラスタリングによる電圧遷移分類のために使用することができる。地下配電網において、数秒で数百の信号遷移を捕捉することができ、それらの多くが部分放電（ＰＤ）である。各遷移ごとの信号処理およびパターン認識には時間がかかる。Ｍ個の遷移をＩ個のクラスタにグループ化することにより、処理時間が劇的に短縮され、対応するＩ個のシグネチャの信号対雑音比が大幅に高められる。クラスタリングは、数百次元Ｎにわたって行うことができ、各次元が信号時間サンプルに対応する。遷移の時間位置は時間ジッタによって損なわれるので、Ｔ回の異なる時間整合に関して距離関数が計算される。ｋ−ｍｅａｎｓアルゴリズムと同様のヒューリスティックが、「スフェアハードニング（ｓｐｈｅｒｅ ｈａｒｄｅｎｉｎｇ）」現象に基づいて説明され、Ｉ個のクラスタに関してＯ（Ｔ×Ｎ×Ｍ×Ｉ）の複雑さを有する。クラスタリングプロセスの精度を査定するため、およびこの方法のいくつかのパラメータを最適化するために、様々なツールが提案される。

【0019】

本発明により処理すべき信号は、例えば、１Ｇｓ／ｓでサンプリングし、フィルタし、補間し、切り捨てることができる。数百個の時間サンプルＮが遷移パターンを表す。ＰＤ放出位置「ｉ」に関して、様々な振幅にわたる一定の正規化時間シグネチャｓ_i（ｔ）を仮定すると、
ｘ_m（ｔ）＝ａ_m・ｓ_i（ｔ−ｔ_m）＋ｎ_m（ｔ） （１）
が得られる。これは、遷移シグネチャの測定「ｍ」の実現であり、ここで、ａ_mは実現振幅であり、ｔ_mは実現遅延であり、ｎ_mは加法性雑音である。対応するモデリングは、打切りにより、
Ｘ_mn＝ａ_m・Ｓ_n-dm＋ｎ_mn （２）
となる。１つのアナログ入力スケールにわたって取られる連続する測定値を、シーケンスと呼ぶ。ａ_mダイナミックレンジは、１シーケンスに関して１０ｄＢ未満でよく、すなわち、トリガ設定レベルに対するクリッピングレベルの比である。第１のクラスタリングステップは、固定スケールに関して行うことができる。次いで、第２のステップで、様々なＡ／Ｄスケールから得られるクラスタをマージすることによって、フルダイナミックレンジを得ることができる。いくつかの場合には、ａ_mダイナミックレンジは、３０ｄＢを超えることがある。

【0020】

ガウス雑音を仮定し、ａ_mダイナミクスを無視すると、

【数3】

投射は、「ｉ」シグネチャ
Ｓ_i＝｛Ｓ_i,1，Ｓ_i,2，．．．，Ｓ_i,N｝ （３）
を中心とする超球を示す。ここで、測定値
Ｘ_m＝｛Ｘ_m,1，Ｘ_m,2，．．．，Ｘ_m,N｝ （４）
は、超球境界に近い。境界の厚さは、計量、測定信号対雑音比（ＳＮＲ）、および時間サンプルＮの数の関数である。ユークリッド計量では、距離

【数4】

は、ｘ_m∈クラスタ「ｉ」に関して、期待平均値

【数5】

と、標準偏差

【数6】

とを有する。

【0021】

図２を参照すると、半径および超球境界厚さは、

【数7】

において、それぞれｒ_iおよび２・σ_iで示される。超球半径に対する境界厚さの比は、Ｎ→∞のときには０となる傾向がある。この現象は、スフェアハードニングと呼ばれる。多数の雑音サンプルを使用して計算すると、距離Ｘ_m−Ｓ_iは、ほとんど一定でない。超球の境界付近以外では、超球内に測定値はない。

【0022】

【数8】

では、クラスタ確率密度は、同様の半径および厚さを有する分散されたシェルのように見える。大きなａ_mダイナミクスが存在する場合、単一点のシグネチャが、軸の原点に向いたロッドに置き換えられる。対応するシェルは、ロッドの軸に沿って拡げられる。シェル厚さは、ロッドの方向で増加する。シグネチャは、異なる振幅ダイナミクス伸長関数を有する分散されたシェルのように見える。適切な計量の使用により、このシェル歪みを一部克服することができる。

【0023】

動的クラスタリングでは、クラスタの数およびクラスタ重心位置を動的に調節することができる。作業次元は、Ｉ_max（最大許容シグネチャ）およびＩ_min（最小許容シグネチャ）によって制限されることがある。シーケンスのプロセス中に、シグネチャに対する平均距離測定値の移動推定（ｒｕｎｎｉｎｇ ｅｓｔｉｍａｔｉｏｎ）を行うことができる。平均距離は、雑音と振幅ダイナミクスの寄与を含む。最大許容距離は、平均距離ｒに距離係数ｃ_distを乗算した値として定義される。この係数は、約１．５に設定することができる。プロセスは、比較およびクラスタリング規則に基づいて以下のように進むことがある。
− 新たなシグネチャとして、第１のＩ_min測定値を記憶する；
− 次の測定値に関して、シグネチャに対する最も近い距離測定値を見出す；
− 最も近い距離が最大許容距離未満である場合、
− 測定値を最も近いシグネチャにマージして、平均距離ｒを再推定する、
− そうでない場合、シグネチャ対シグネチャ距離を計算して、最も近いシグネチャ対シグネチャ距離を見出す；
− 距離が最大許容距離未満である場合、
− 最も近い２つのシグネチャをマージして、平均距離ｒを再推定する、
− そうでない場合には、測定値を新たなシグネチャとするか、あるいは、Ｉ_maxを超えている場合には、測定値を最も近いシグネチャとマージさせる；
− 最後の測定後、２Ｔ→Ｔに関して、時間整合のスパンが倍増され：
１．シグネチャ対シグネチャ距離を計算し；
２．最も近いシグネチャ対シグネチャ距離を見出し；
３．最も近い距離が最大許容距離未満である場合：
４．最も近い２つのシグネチャをマージし、新たなシグネチャに対するシグネチャ距離を計算し、結果を、シグネチャ対シグネチャ距離三角行列に挿入し、ステップ２に戻る：
５．そうでない場合には、停止する。

【0024】

スフェアハードニングを使用して結果を精密化することができる。すなわち、球面境界から離れた測定値を再処理して、別のクラスタに割り振りし直すことができる。しきい値は、

【数9】

に比例的に固定することができ、ここで、ｖａｒ（σ_i）は、σ_i推定値に対する分散の不確かさである。

【0025】

ｒ_i≒ｒ_jと仮定すると、平均距離

【数10】

は、クラスタ生成中に計算された最小距離Ｄ_m,iから再帰的に推定される。この平均距離

【数11】

は、測定値の雑音と、シグネチャ位置誤差とを含む。右側の因子、すなわちクラスタ母集団Ｐ_iの関数は、シグネチャ位置の分散を考慮に入れる。係数ｂは、重み付け係数＞１である。

【0026】

新たなシグネチャ計算に関して、マージプロセス

【数12】

は、

【数13】

であり、ここで、Ｐ_iおよびＰ_jは、クラスタ母集団である。マージ前、シグネチャは、最小距離に対して時間整合される。シグネチャ時間整合も、クラスタ母集団によって重み付けされる。

【数14】

ここで、ｄは、サンプルの数で表現されるシグネチャ間の距離である（式５）。測定値の到着の順序は、最終結果へわずかに影響を及ぼすが、最後には、

【数15】

となる。

【0027】

様々な計量のうち、試験されたＴ回の時間整合にわたって計算された二乗距離

【数16】

の最小値が、最良の結果を生み出すことができる。右側の項を最小にすることが、相関を最大にすることに相当することに留意されたい。多くの小さな雑音パターンがいくつかのＰＤシグネチャに相関するので、相関の最大化のみでは効果がないことがある。この計量では、ＰＤ振幅のばらつきの影響は、第２の項によって減少される。

【0028】

様々な測定スケールによって生成されるクラスタをマージするために行われる第２のクラスタリングで、相関（すなわち、ブラベ−ピアソン（Ｂｒａｖａｉｓ−Ｐｅａｒｓｏｎ）係数）を使用することができる。この第２のステップでは、クラスタシグネチャＳＮＲは高く、ＰＤクラスタシグネチャと雑音シグネチャを間違うことはあり得ない。

【0029】

計算の複雑さは、シグネチャに対する距離測定に関しては、Ｏ（Ｎ×Ｔ×Ｍ×Ｉ）であり、シグネチャ対シグネチャ距離三角行列に関しては、Ｏ（Ｎ×Ｔ×Ｍ×Ｉ（Ｉ−１）／２）である。係数ｃ_dist、Ｉ_max、およびＩ_minは、計算時間と、誤ったマージが生じる確率との折衷点を見出す。多数の測定値に関して、クラスタの数が安定化されるとき、シグネチャ対シグネチャ距離の計算はもはや必要ない。すなわち全体的な計算の複雑さは、Ｏ（Ｎ×Ｔ×Ｍ×Ｉ）になる傾向がある。

【0030】

精度測定値は、クラスタ分散およびクラスタ重畳に関するいくつかの情報を含むことがある。後者の情報はデータセットにわたって一定であるので、クラスタリングの誤りを表す１つまたは複数の精度推定値を使用してプロセスを最適化することがその狙いである。提案される推定値は、２つの相反する方向に基づく。一方で、クラスタ半径に対するクラスタ間の距離の比として定義される解像力が、クラスタ同士の重畳に関係する。他方で、コヒーレンシー、すなわち合計のクラスタエネルギーに対するコヒーレントエネルギーの比が、クラスタ自体に関係する。

【0031】

解像力

【数17】

は、クラスタ解像能力の尺度である。解像力

【数18】

は、信号対雑音比に対応し、ここで、ＲＭＳクラスタ間距離が信号であり、ＲＭＳクラスタ半径が雑音である。等式

【数19】

が、第１の項の分母を説明し、ΣＰ_i＝Ｍである。

【0032】

クラスタ「ｉ」のコヒーレンシー

【数20】

は、このクラスタのＰ_i要素から計算される。平均コヒーレンシー

【数21】

は、寄与するすべてのクラスタを考慮するシーケンスに関して定義される。最良の解に関して、解像力と平均コヒーレンシーはそれらの最大レベルであり、クラスタリングプロセスでの誤りの累積と共に減少するという仮説を使用することができる。

【0033】

図３は、典型的なＰＤシグネチャを示す。これらの結果は、地下配電網で取られた１７３０個のシーケンス（１１８９３２個の測定値）に関して得られている。最初のローブは、最も鋭く、ここではＰＤ極性を負の値に設定する。解像力の値の範囲１〜１００００は、対数（ｄＢ）スケーリングを必要とする。距離係数ｃ_distは、解像力未満に設定しなければならない。図４の左で、ｃ_dist＝１０に関する０ｄＢでの２０カウントは、ｃ_dist＞解像力であるので、クラスタリングの誤りとして説明される。

【0034】

クラスタコヒーレンシーは、主に、小さな解像力を有するシーケンスに関して変えられる。図５で、シーケンスのほとんどは、大きな解像力を有する。すなわち、ｃ_distの値を１０にシフトしても、コヒーレンシーに対してほとんど影響がない。しかし、同じシフトが、いくつかのシーケンスのクラスタ母集団を大幅に減少させる（図６）。図６は、最大許容シグネチャＩ_maxを３０に設定することができることを示す。図７は、ｃ_distと共に減少する処理時間を示す。処理時間とクラスタリングの誤りの確率との折衷点は、データ自体と、診断目的での望ましい信頼度とによって設定される。後者および他の結果に基づいて、ｃ_dist＝１．５の設定が、データ母集団の９９％超を対象とし、測定時間（約８秒）に近い処理時間を生み出す。

【0035】

時間領域での遷移のクラスタリングの寄与は、ＰＲＰＤ図でクラスタ区別によって認識される。図８は、合計１１８２個の測定値に関する９個の他のクラスタのうち、６５個の測定値からなる１つのクラスタを強調する。

【0036】

クラスタシグネチャに対して適用されるさらなる信号処理による、最初のステップの解析としての時間領域での遷移信号のクラスタリングは、少なくとも以下の点において有利である。すなわち、（１）情報が、多数の測定値ではなくいくつかのシグネチャに減少され、（２）ＳＮＲシグネチャが、クラスタ母集団と共に増加し、また（３）後処理の時間が短縮される。さらに、重畳されたクラスタが、ＰＲＰＤ図において区別される。説明した次善のヒューリスティックも、高速で正確である。複数のフィールドデータを使用する試験を、パラメータを調整するため、および本発明による方法での計量の選択を行うために実現することができる。

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】