(57)【要約】  ネットワーク化されたマイクログリッドのシステム、方法、技術および装置が開示される。例示的な一実施形態は、ネットワーク化されたマイクログリッドシステムから偽データを除去するための方法であって、マイクログリッド制御システムで第１のローカル状態推定値を計算するステップと、マイクログリッド制御システムおよび複数のマイクログリッドで複数のグローバル状態推定値を計算するステップと、複数のグローバル状態推定値が収束しなかったと判断するステップと、マイクログリッド制御システムおよび複数のマイクログリッドを用いて、マイクログリッド制御システムで計算された複数のグローバル状態推定値のうちの第１のグローバル状態推定値における偽データを検出するステップと、複数のグローバル状態推定値から第１グローバル状態推定値を除去するステップと、複数のグローバル状態推定値のうち残りのグローバル状態推定値が収束するまで、マイクログリッド制御システムおよび複数のマイクログリッドで、複数のグローバル状態推定値のうち残りのグローバル状態推定値を反復的に更新するステップとを含む。
2022501994.tif 000002
【特許請求の範囲】
【請求項１】
第１のマイクログリッドおよび複数のマイクログリッドを含むネットワーク化されたマイクログリッドシステムから偽データを除去するための方法であって、
第１の複数のローカル測定値を用いてマイクログリッド制御システムで第１のローカル状態推定値を計算するステップと、
前記マイクログリッド制御システムで、前記複数のマイクログリッドから複数のローカル状態推定値を受信するステップと、
前記第１のローカル状態推定値および前記複数のローカル状態推定値を用いて、前記マイクログリッド制御システムおよび前記複数のマイクログリッドで複数のグローバル状態推定値を計算するステップと、
前記複数のグローバル状態推定値が収束しなかったと判断するステップと、
前記マイクログリッド制御システムおよび前記複数のマイクログリッドを用いて、前記マイクログリッド制御システムで計算された前記複数のグローバル状態推定値のうちの第１のグローバル状態推定値における偽データを検出するステップと、
前記複数のグローバル状態推定値から前記第１のグローバル状態推定値を除去するステップと、
前記複数のグローバル状態推定値のうち残りのグローバル状態推定値が収束するまで、前記マイクログリッド制御システムおよび前記複数のマイクログリッドで、前記複数のグローバル状態推定値のうち前記残りのグローバル状態推定値を反復的に更新するステップとを含む、方法。
【請求項２】
前記複数のマイクログリッドの第１の部分に関連する直接通信の際に、前記第１のマイクログリッド制御システムで、前記第１の部分から情報を受信し、その後、前記第１のマイクログリッド制御システムとの間接通信の際に、前記第１の部分から前記複数のマイクログリッドの第２の部分に関連する情報を受信することによって、通信ネットワークトポロジを決定するステップを含む、請求項１に記載の方法。
【請求項３】
複数のグローバル状態推定値は、各グローバル状態推定値が閾値未満の標準偏差を有するときに収束する、請求項１に記載の方法。
【請求項４】
前記複数のグローバル状態推定値を計算するステップは、複数回の反復を用いて行われ、前記複数のグローバル状態推定値が収束しなかったと判断するステップは、前記複数回の反復の後に行なわれる、請求項１に記載の方法。
【請求項５】
前記第１のグローバル状態推定値における偽データを検出するステップは、前記複数のマイクログリッドのうち１つのマイクログリッドから第１のマイクログリッド制御システムによって受信されるグローバル状態推定値に適用される重みに対応する第１の信頼因数と、前記第１のマイクログリッド制御システムによって計算される前記グローバル状態推定値に適用される重みに対応する第２の信頼因数とを含む複数の信頼因数を含む動的信頼行列を更新するステップを含む、請求項１に記載の方法。
【請求項６】
前記第１のグローバル推定値における偽データを検出するステップは、以下の式セットを用いて前記複数のグローバル状態推定値を更新するステップを含み、
【数１】
2022501994.tif 000030
ここで、ｋはエージェントマイクログリッド制御システムであり、ｌは近隣マイクログリッド制御システムであり、
【数２】
2022501994.tif 000031
は、反復ｉでのエージェントマイクログリッド制御システムｋのグローバル状態推定値を示し、ωｋ，ｉは、反復ｉでのエージェントマイクログリッド制御システムｋに関する中間変数であり、μｋは、エージェントマイクログリッド制御システムｋの非負更新パラメータであり、∇ωＪｌ（ωｋ，ｉ−１）は、反復ｉ−１での中間状態ωのエージェントマイクログリッド制御システムｋに関する確率勾配であり、νｋは忘却因数であり、
【数３】
2022501994.tif 000032
は瞬間誤差メトリックである、請求項１に記載の方法。
【請求項７】
前記複数のグローバル状態推定値の各グローバル状態推定値は、前記第１のマイクログリッドのバスの電圧大きさに対応する、請求項１に記載の方法。
【請求項８】
マイクログリッドシステムにおける偽データを検出および軽減するための方法であって、
第１のマイクログリッドの第１のマイクログリッド制御システム、第２のマイクログリッドの第２のマイクログリッド制御システム、および第３のマイクログリッドの第３のマイクログリッド制御システムを動作させるステップと、
第１の複数の測定値を用いて前記第１のマイクログリッド制御システムで第１のローカル状態推定値を計算し、第２の複数の測定値を用いて前記第２のマイクログリッド制御システムで第２のローカル状態推定値を計算し、第３の複数の測定値を用いて前記第３のマイクログリッド制御システムで第３のローカル状態推定値を計算するステップと、
前記第２のマイクログリッド制御システムおよび第前記３のマイクログリッド制御システムで前記第１のローカル状態推定値を受信し、前記第１のマイクログリッド制御システムおよび前記第３のマイクログリッド制御システムで前記第２のローカル状態推定値を受信し、前記第１のマイクログリッド制御システムおよび前記第２のマイクログリッド制御システムで前記第３のローカル状態推定値を受信するステップと、
前記第１のマイクログリッド制御システムで前記第１のマイクログリッドの電気的特性の第１のグローバル状態推定値を計算し、前記第２のマイクログリッド制御システムで電気的特性の第２のグローバル状態推定値を計算し、前記第３のマイクログリッド制御システムで電気的特性の第３のグローバル状態推定値を計算するステップと、
前記第１のマイクログリッド制御システムで、前記第１のグローバル状態推定値、前記第２のグローバル状態推定値、および前記第３のグローバル状態推定値が収束しないと判断するステップと、
動的信頼行列を用いて、前記第１のマイクログリッド制御システムで前記第１のグローバル状態推定値を反復的に更新し、前記第２のマイクログリッド制御システムで前記第２のグローバル状態推定値を反復的に更新し、前記第３のマイクログリッド制御システムで前記第３のグローバル状態推定値を反復的に更新するステップと、
前記第１のグローバル状態推定値、前記第２のグローバル状態推定値、および前記第３のグローバル状態推定値を反復的に更新するステップを用いて、前記第１のグローバル状態推定値内の偽データを検出するステップと、
前記第２のグローバル状態推定値および前記第３のグローバル状態推定値が収束するまで、前記第２のマイクログリッド制御システムでの前記第２のグローバル状態推定値の更新および前記第３のマイクログリッド制御システムでの前記第３のグローバル状態推定値の更新を続けるステップとを含む、方法。
【請求項９】
第１の反復中に、前記第１のマイクログリッド制御システムで、前記第２のマイクログリッドに関する情報を前記第２のマイクログリッド制御システムから受信し、前記第１の反復中に、前記第２のマイクログリッド制御システムで、前記第３のマイクログリッドに関する情報を前記第３のマイクログリッド制御システムから受信し、第２の反復中に、前記第１のマイクログリッド制御システムで、前記第３のマイクログリッドに関する情報を前記第２のマイクログリッド制御システムから受信することによって、通信ネットワークトポロジを決定するステップを含む、請求項８に記載の方法。
【請求項１０】
前記第１のローカル状態推定値、前記第１のグローバル状態推定値、前記第２のグローバル状態推定値、および前記第３のグローバル状態推定値は、前記第１のマイクログリッドのバスの推定電圧である、請求項８に記載の方法。
【請求項１１】
前記第１のマイクログリッド制御システムは、第２の複数のローカル測定値または第３の複数のローカル測定値を受信しない、請求項８に記載の方法。
【請求項１２】
前記第１のグローバル状態推定値、前記第２のグローバル状態推定値、および前記第３のグローバル状態推定値を計算するステップは、拡散アルゴリズムを用いるステップを含み、前記第２のグローバル状態推定値および前記第３のグローバル状態推定値の更新を続けるステップは、前記拡散アルゴリズムを用いるステップを含む、請求項８に記載の方法。
【請求項１３】
前記動的信頼行列は、第１のマイクログリッド制御システムによって前記第１のグローバル状態推定値に適用される重みに対応する第１の信頼因数と、前記第１のマイクログリッド制御システムによって前記第２のグローバル状態推定値に適用される重みに対応する第２の信頼因数とを含む複数の信頼因数を含む、請求項８に記載の方法。
【請求項１４】
前記動的信頼行列を用いて前記第１のグローバル状態推定値、前記第２のグローバル状態推定値、および前記第３のグローバル状態推定値を反復的に更新するステップは、前記第１のマイクログリッド制御システムによって前記第２のグローバル状態推定値内のノイズを検出することに応答して前記第２の信頼因数をゼロに低減するステップを含む、請求項１３に記載の方法。
【請求項１５】
前記第１のグローバル状態推定値、前記第２のグローバル状態推定値、および前記第３のグローバル状態推定値を反復的に更新するステップは、以下の式セットを用いて前記動的信頼行列を更新するステップを含み、
【数４】
2022501994.tif 000033
ここで、ｋはエージェントマイクログリッド制御システムであり、ｌは近隣マイクログリッド制御システムであり、
【数５】
2022501994.tif 000034
は、反復ｉでのエージェントマイクログリッド制御システムｋのグローバル状態推定値を示し、ωｋ，ｉは、反復ｉでのエージェントマイクログリッド制御システムｋに関する中間変数であり、μｋは、エージェントマイクログリッド制御システムｋの非負更新パラメータであり、νｋは忘却因数であり、
【数６】
2022501994.tif 000035
は瞬間誤差メトリックである、請求項８に記載の方法。
【請求項１６】
非一時的コンピュータ可読媒体であって、
第１のマイクログリッドの第１のマイクログリッド制御システムの処理デバイスによって実行されるように構成された命令のセットを含み、前記命令のセットは、
第１の複数の測定値を用いて第１のローカル状態推定値を計算し、
前記第１のマイクログリッドに結合された複数のネットワーク化されたマイクログリッドから複数のローカル状態推定値を受信し、
前記第１のローカル状態推定値および前記複数のローカル状態推定値を用いて第１のグローバル状態推定値を計算し、
前記複数のネットワーク化されたマイクログリッドから複数のグローバル状態推定値を受信し、
前記第１のグローバル状態推定値および前記複数のグローバル状態推定値が収束しないと判断し、
前記第１のグローバル状態推定値および前記複数のグローバル状態推定値における検出されたノイズに応じて更新されるように構成された動的信頼行列を用いて、前記第１のグローバル状態推定値および前記複数のグローバル状態推定値を反復的に更新することによって、前記複数のグローバル状態推定値のうち少なくとも１つのグローバル状態推定値が偽データを含むと判断し、
前記複数のグローバル状態推定値のうち偽データを含むと判断されなかったグローバル状態推定値で前記第１のグローバル状態推定値が収束するまで前記第１のグローバル状態推定値を反復的に更新する、
のに有効である、非一時的コンピュータ可読媒体。
【請求項１７】
前記複数のフィールドデバイスは、前記第１のマイクログリッドの複数のバス電圧を測定するように構成された複数の電圧センサを含む、請求項１６に記載の非一時的メモリデバイス。
【請求項１８】
処理デバイスによって実行されるように構成された前記命令のセットは、前記複数のネットワーク化されたマイクログリッドの第１の部分に関連する直接通信の際に前記第１の部分から情報を受信し、その後、前記第１のマイクログリッド制御システムとの間接通信の際に、前記第１の部分から、前記複数のネットワーク化されたマイクログリッドの第２の部分に関連する情報を受信することによって、通信ネットワークトポロジを決定するのに有効である、請求項１６に記載の非一時的メモリデバイス。
【請求項１９】
前記動的信頼行列は、複数のマイクログリッドのうちの１つのマイクログリッドから前記第１のマイクログリッド制御システムによって受信されるグローバル状態推定値に適用される重みに対応する第１の信頼因数と、前記第１のグローバル状態推定値に適用される重みに対応する第２の信頼因数とを含む複数の信頼因数を含む、請求項１６に記載の非一時的メモリデバイス。
【請求項２０】
前記複数のグローバル状態推定値のうち少なくとも１つのグローバル状態推定値が偽データを含むと判断することは、拡散アルゴリズムを用いることと、以下の式セットを用いて前記動的信頼行列を更新することとを含み、
【数７】
2022501994.tif 000036
ここで、ｋはエージェントマイクログリッド制御システムであり、ｌは近隣マイクログリッド制御システムであり、
【数８】
2022501994.tif 000037
は、反復ｉでのエージェントマイクログリッド制御システムｋに関するグローバル状態推定値を示し、ωｋ，ｉは、反復ｉでのエージェントマイクログリッド制御システムｋに関する中間変数であり、μｋは、エージェントマイクログリッド制御システムｋの非負更新パラメータであり、νｋは忘却因数であり、
【数９】
2022501994.tif 000038
は瞬間誤差メトリックである、請求項１６に記載の非一時的メモリデバイス。
【発明の詳細な説明】
【技術分野】
【０００１】
背景
本開示は、概して、少なくとも１つの電源と１つの負荷とを含む電力ネットワークの相互接続された部分であるネットワーク化されたマイクログリッドに関するものであって、当該１つの負荷は、電力ネットワークの残りの部分から切断された場合に独立して動作することができる。ネットワーク化されたマイクログリッドシステムにおける状態推定は、オンライン監視およびネットワーク認識を可能にする。いくつかのシステムでは、状態推定は一元化された場所で実行される。すべてのデータを１つの中央コントローラに送信することで、状態推定が実行可能となる。電力ネットワーク内のマイクログリッドおよび通信デバイスの数が増加するのに応じて、リアルタイムでの集中型状態推定のための計算性能制約が指数関数的に増加する。分散型状態推定は、集中型状態推定の一代替例であって、１つの集中型制御システムではなく複数のローカル制御システムによって状態推定が実行される。
【背景技術】
【０００２】
連邦政府による資金提供を受けた研究開発
本発明は、契約番号ＤＥ−ＯＥ００００８３１に関して米国エネルギ省からの政府支援を受けて行われた。米国政府は本発明について一定の権利を有する。
【発明の概要】
【課題を解決するための手段】
【０００３】
開示の概要
本開示の例示的な実施形態は、ネットワーク化されたマイクログリッドのための固有のシステム、方法、技術および装置を含む。本開示のさらなる実施形態、形態、目的、特徴、利点、局面および利益が以下の記載および添付の図面から明らかになるだろう。
【図面の簡単な説明】
【０００４】
【図１】例示的なネットワーク化されたマイクログリッドシステムを示す。
【図２】偽データ検出および軽減のための例示的なプロセスを示すフローチャートである。
【図３】図１のネットワーク化されたマイクログリッドシステムでのグローバル状態推定を示すグラフのセットである。
【図４】図１のネットワーク化されたマイクログリッドシステムでの偽データ軽減を示すグラフのセットである。
【図５】図１のネットワーク化されたマイクログリッドシステムについての例示的な初期化プロセスを示す表である。
【図６Ａ】例示的な分散型状態推定プロセスに関する収束を示すグラフである。
【図６Ｂ】代替的な分散型状態推定プロセスに関する収束を示すグラフである。
【図７】例示的な分散型状態推定プロセスおよび従来の集中型状態推定プロセスに関する計算性能を示す表である。
【図８】例示的なマイクログリッド制御システムを示すブロック図である。
【発明を実施するための形態】
【０００５】
具体的な実施形態の詳細な説明
図１を参照すると、ネットワーク化されたマイクログリッドシステム１００が示されている。システム１００が、いくつか例を挙げると、ユーティリティグリッド配電システム、産業プラント配電システム、および車両配電システムを含むさまざまな用途で実現され得ることを認識されたい。システム１００のトポロジは、説明のために例示されるものであって本開示を限定するものとして意図されるものではないことを認識されたい。たとえば、システム１００が含み得るマイクログリッドは、より多くてもよくまたはより少なくてもよく、システム１００のマイクログリッドは、各マイクログリッドが少なくとも１つの他のマイクログリッドに結合されている限り、如何なる構成で配置されてもよい。さらに、各マイクログリッドのトポロジは、説明のために例示されるものであって本開示を限定するものとして意図されていないことを認識されたい。たとえば、マイクログリッド１１０が含み得るバス、電源、エネルギ貯蔵デバイス、負荷およびフィールドデバイスはそれぞれ、より多くてもよくまたはより少なくてもよい。システム１００は単線図で示されているが、当該システム１００は、単相交流（alternating current：ＡＣ）電力または多相ＡＣ電力、または直流（direct current：ＤＣ）電力を伝送するように構成されてもよい。
【０００６】
システム１００は、マイクログリッド１１０、１２０、１３０、１４０および１５０を含む。各マイクログリッドは、少なくとも１つの負荷と、複数のバスのうちの１つに結合される少なくとも１つの電源とを含む。各マイクログリッドは、マイクログリッドの電気的特性または物理的特性を測定するように構成された複数のフィールドデバイスを含む。たとえば、マイクログリッド１５０は、バス１５３を含む複数のバスと、負荷１５５を含む複数の負荷と、ＡＣ電源１５７と、複数のフィールドデバイス１５２とを含む。
【０００７】
ＡＣ電源１５７は、ソーラーパネルアレイ、風力タービン、天然ガス発生器、または、電力を発生するように構成された他の任意のデバイスもしくはシステムを含み得る。負荷１５５は、電力を消費するように構成された任意の種類のデバイスまたはシステムであり得る。複数のフィールドデバイス１５２は、いくつか例を挙げると、電圧センサ、インテリジェント電子デバイス（intelligent electronic device：ＩＥＤ）、遠隔端末ユニット（remote terminal unit：ＲＴＵ）、リレー、再閉器、電流センサ、変圧器、および変流器を含み得る。
【０００８】
各マイクログリッドはまた、複数の通信チャネルを介して複数のフィールドデバイスに結合されたマイクログリッド制御システムを含む。たとえば、マイクログリッド１５０は、複数の通信チャネルを介して複数のフィールドデバイス１５２に結合されるマイクログリッド制御システム１５１を含む。各マイクログリッド制御システムは、マイクログリッドの電気的特性および物理的特性を監視し、マイクログリッド内の制御可能なデバイスを動作させ、マイクログリッドに結合されたデバイスを故障および他の状況から保護するように構成される。
【０００９】
マイクログリッド制御システム１１１、１２１、１３１、１４１および１５１は、１つ以上の近隣のマイクログリッド制御システムと結合された通信チャネルを含む通信ネットワーク１６０を介して通信するように構成される。図示した実施形態では、マイクログリッド制御システム１１１はマイクログリッド制御システム１２１および１３１と通信し、マイクログリッド制御システム１２１はマイクログリッド制御システム１１１のみと通信し、マイクログリッド制御システム１３１はマイクログリッド制御システム１１１、１４１および１５１と通信し、マイクログリッド制御システム１４１はマイクログリッド制御システム１３１のみと通信し、マイクログリッド制御システム１５１はマイクログリッド制御システム１３１のみと通信する。
【００１０】
フィールドデバイスとマイクログリッド制御システムとの間、ならびにマイクログリッド間の通信チャネルは有線でもあってもよく、または無線であってもよい。各通信チャネルは、いくつか例を挙げると、ＩＥＣ６１８５０、ＯＰＣ  ＵＡ、Ｐｕｂ−Ｓｕｂ、ＸＭＰＰまたはＤＤＳなどの標準化された通信プロトコルを用いてもよい。いくつかの実施形態では、標準化された通信プロトコルが分散型状態推定のためにマイクログリッドの挙動を定義するのに不十分である場合、やり取りされるメッセージのセマンティクスが共通の知識表現を用いて定義され得る。
【００１１】
たとえば、インテリジェント物理エージェント財団（the Foundation for Intelligent Physical Agents：ＦＩＰＡ）の規格は、やり取りされたメッセージのセマンティクスを記述するためのエージェント通信言語（Agent Communication Language：ＡＣＬ）を組込んでいる。ＡＣＬメッセージの内容は、当該メッセージを通じて渡されるパラメータおよびアクションのオブジェクトを含む。送信者パラメータおよび受信者パラメータは、送信者および意図された受信者エージェントの名前を指定する。ＡＣＬメッセージは、いくつか例を挙げると、通知、要求、拒否、加入および伝搬などの通信動作を表わす。ＩＥＣ６１８５０通信プロトコルは、ＡＣＬ通信動作に基づく抽象通信サービスインターフェイス（Abstract Communication Service Interface：ＡＣＳＩ）として非標準ＩＥＣ６１８５０論理ノードを用いて拡張され得るとともに、ＧＯＯＳＥメッセージにマッピングすることができる。
【００１２】
既存の分散型状態推定システムには多くの欠点および不利点がある。サイバーセキュリティの増強の要求や、プライバシー問題の増加に伴う要求が依然として満たされていないままである。たとえば、マイクログリッドは、分散型状態推定を実行するために、他のマイクログリッドのセキュリティだけでなく通信システムのセキュリティにも依拠しなければならない。カスケード式電力ネットワークの障害を引起こすには、１つのマイクログリッドのセキュリティの違反だけで十分であるだろう。さらに、独立して動作するマイクログリッドは、プライバシー問題のために、未処理の測定値などの他のマイクログリッドと共有されるデータの種類を制限する可能性がある。当技術分野におけるこれらおよび他の欠点に鑑みると、本明細書に開示される固有の装置、方法、システムおよび技術が大いに必要とされている。以下でより詳細に説明するように、システム１００のマイクログリッド制御システムは各々、他のマイクログリッド制御システムから受信した状態推定値およびエージェントマイクログリッド制御システムに関連付けられたマイクログリッドのローカル測定値のみを用いてグローバル状態推定を実行するように構成されている。
【００１３】
図２を参照すると、完全分散型状態推定を用いて、ネットワーク化されたマイクログリッドシステムにおいて偽データを検出および軽減するためのプロセス２００が示されている。図１に示されるネットワーク化されたマイクログリッドシステム１００を参照して、以下においてプロセス２００を説明する。具体的には、プロセス２００は、１つのマイクログリッド制御システムを参照して説明されるが、ネットワーク化されたマイクログリッドシステム１００の各マイクログリッド制御システムが同時にプロセス２００を実行することが理解されるはずである。プロセス２００は、他の形態のネットワーク化されたマイクログリッドシステムと組合わせて用いられてもよい。たとえば、プロセス２００についての１つ以上の局面の省略、さらなる条件および動作の追加、ならびに／または、動作および条件の別個のプロセスへの再編成または分離を含む、プロセス２００に対するいくつかの変形およ変更が企図されることがさらに認識されるだろう。
【００１４】
プロセス２００は開始動作２０１を含み、動作２０３へと進む。動作２０３で、マイクログリッド制御システムは、他のマイクログリッド制御システムを発見するとともにマイクログリッド間の位相角に関連する情報を受信することによって、マイクログリッド通信ネットワークのトポロジを決定する。マイクログリッド制御システムは、状態推定のための電力ネットワークモデルを必要としない。代わりに、ネットワーク化されたマイクログリッドシステムは、システム１００内のマイクログリッド制御システムの数およびマイクログリッド間の位相角差を識別するために完全分散型の反復プロセスを用いる。特定の実施形態では、システム１００の各マイクログリッドは、フェーザ測定ユニットを含んでおり、マイクログリッド間の位相角差を決定する必要はない。
【００１５】
マイクログリッド制御システムがマイクログリッド通信ネットワークトポロジを決定すると、プロセス２００は動作２０５に進む。動作２０５で、マイクログリッド制御システムは、マイクログリッド制御システムと同じマイクログリッドの複数のフィールドデバイスから受信したデータを用いてローカル状態推定を実行する。ローカル状態推定は、推定すべきマイクログリッドの電気的特性または物理的特性に関する入力データを用いてマイクログリッドの電気的特定または物理的特性が推定されるプロセスである。複数のフィールドデバイスからの入力データは、いくつか例を挙げると、電圧、電流、有効電力流、無効電力流の測定値、または他の種類の情報、たとえば回路遮断器状態などを含み得る。各々の測定値または推定値は、大きさ、位相角、またはこれら両方を含み得る。
【００１６】
特定の実施形態では、ローカル状態推定は加重最小二乗法を用いて実行される。各マイクログリッドが知っているのは、それ自体のローカル測定ベクトルｚ、ローカル共分散行列Ｒ、およびローカルヤコビ行列Ｈだけである。ローカル測定ベクトルｚは、バス電圧、電流注入および電流の流れなどの、マイクログリッドにおけるフィールドデバイスからのデータを含む。行列Ｈは、ｈ（ｘ）のヤコビ行列であって、各マイクログリッドについての非線形測定関数
【００１７】
【数１】
2022501994.tif 000003
【００１８】
である。行列Ｒは、各測定値の標準偏差を含む測定共分散行列
【００１９】
【数２】
2022501994.tif 000004
【００２０】
である。反復プロセスにおいて以下の式のセットを用いて、各マイクログリッドがローカル状態推定を実行する。ここで、Ｇ（ｘ）は利得行列であり、Ｗは共分散行列Ｒの逆数であり、ｘはローカル状態推定ベクトルであり、ｉは反復数である。
【００２１】
【数３】
2022501994.tif 000005
【００２２】
プロセス２００が動作２０７に進む。動作２０７で、マイクログリッド制御システムは、発見されたマイクログリッドトポロジを用いるグローバル状態推定値、それ自体の計算されたローカル状態推定値、およびネットワーク化されたマイクログリッドの他のマイクログリッド制御システムによって計算されるローカル状態推定値を初期化する。グローバル状態推定は、ネットワーク化されたマイクログリッドシステムの全マイクログリッドの電気的特定または物理的特性が、推定されるべきネットワーク化されたマイクログリッドシステムの電気的特定または物理的特性に関連する入力データを用いて推定されるプロセスである。たとえば、グローバル状態推定は、ネストされたマイクログリッドシステム内のすべてのバスのバス電圧を推定することと、ネストされたマイクログリッドシステム内のすべてのマイクログリッドを流れる電流を推定することとを含み得る。プロセス２００は、各マイクログリッド制御システムが知っているのがそれ自体の測定値だけであるので、ローカル情報ポリシーおよびプライバシーを考慮しながらネットワーク化されたマイクログリッドシステムのグローバル状態推定を達成するために分散型の２層レベルアルゴリズムを実現する。グローバル状態推定中にマイクログリッド間でやり取りされるデータは、以下に説明されるように、反復的に計算される初期グローバル状態推定値およびグローバル状態推定値を決定するためのローカル状態推定値のみを含む。各マイクログリッド制御システムは、グローバル状態推定を実行するために、別のマイクログリッドまたはネットワークモデルからの未処理のフィールドデバイスデータを必要としない。
【００２３】
マイクログリッド制御システムは、ネットワーク化されたマイクログリッドシステムにおける他のマイクログリッド制御システムの各々からベクトルの形でローカル状態推定値を受信する。各ローカル状態推定ベクトルを用いて、マイクログリッド制御システムは、ローカル状態推定ベクトルを１つのベクトルに組合わせることによってグローバル推定ベクトルを生成する。
【００２４】
動作２０７中にグローバル状態推定ベクトルを生成することに加えて、マイクログリッド制御システムは、マイクログリッド間の関係の各々に重み付けされた値を割当てるために信頼行列を生成する。ネットワーク化されたマイクログリッドシステムは、複数のマイクログリッド制御システムとして動的な無向グラフで表現され得る。信頼行列または近隣行列である
【００２５】
【数４】
2022501994.tif 000006
【００２６】
は、エージェントマイクログリッド制御システムｋがネットワーク化されたマイクログリッドシステムにおける近隣マイクログリッド制御システムｌから受信した通信に伝達する非負重みを表わす。
【００２７】
係数ａｋｌは、各マイクログリッド制御システムが互いに対して課す信頼因数である。係数は対称的である必要はない（すなわち、ａｋｌ≠ａｌｋ）が、各係数は０以上でなければならず、１つのマイクログリッド制御システムによって近隣に割当てられる信頼因数の合計は１に等しくなければならず、近隣ではないマイクログリッド制御システムに割当てられる信頼因数は０でなければならない。
【００２８】
【数５】
2022501994.tif 000007
【００２９】
信頼因数は、エージェントマイクログリッド制御システムとその近隣との間の信頼因数がエージェントマイクログリッド制御システムまたは近隣マイクログリッド制御システムのいずれかに対する接続されるマイクログリッド制御システムの最大数の逆数に等しくなるように、かつ、マイクログリッド制御システムとマイクログリッド制御システム自体との間の信頼因数がマイクログリッド制御システムとその近隣との間の信頼因数の和を引いたものに等しくなるように、割当てられてもよい。式セット（３）は、上述した割当てを示している。ここで、ｎｋは、マイクログリッド制御システムｋ自体を含むエージェントマイクログリッド制御システムｋのうち接続されたマイクログリッド制御システムの数であり、ｎｌは、マイクログリッド制御システムｌ自体を含む近隣マイクログリッド制御システムｌのうち接続されたマイクログリッド制御システムの数であり、Ｎｌは近隣ｌの近傍である。
【００３０】
【数６】
2022501994.tif 000008
【００３１】
言い換えれば、ｋからｌまでの信頼因数は、いずれかのマイクログリッド制御システムが有する近隣の最大数に依存しており、自己重みは、Ａ内の各行（および列）の合計を１にする。信頼行列は、システムトポロジの変化に応じて再計算されてもよいが、それ以外の場合には静的なままである。
【００３２】
プロセス２００が動作２０８に進む。動作２０８で、動作２０７中に生成された信頼行列およびグローバル推定ベクトルを用いて、反復グローバル状態推定手順が実行される。グローバル状態推定ベクトルおよび信頼行列を用いることで、ネットワーク化されたマイクログリッドシステムのマイクログリッド制御システムは、拡散戦略およびピアツーピア反復通信を用いてグローバル状態推定を実行する。拡散アルゴリズムは、同じパラメータベクトルであるω０、Ｎ×１のサイズを推定するＮ個のマイクログリッド制御システムの集合を考慮に入れる。拡散アルゴリズムは、Combine-then-Adapt（ＣＴＡ）およびAdapt-then-Combine（ＡＴＣ）という２種類の戦略を用いる。ＡＴＣは以下のように要約することができる。すなわち、すべての反復瞬間ｉにおいて、ＡＴＣ戦略は２つのステップを実行する。第１のステップは、マイクログリッド制御システムｋがその近隣からそれらのモーメント｛Ｒｕ，ｌ，ｒｄｕ，ｌ｝を受信する情報交換ステップである。マイクログリッド制御システムｋは、この情報を組合わせ、これを用いて、その既存の推定値ωｋ，ｉ−１を中間値
【００３３】
【数７】
2022501994.tif 000009
【００３４】
に更新する。ネットワーク化されたマイクログリッドシステムにおける他のすべてのマイクログリッド制御システムは、同様のステップを実行するとともに、それらの近隣からの情報を用いることによってそれらの推定値｛ωｋ，ｉ−１｝を中間推定値
【００３５】
【数８】
2022501994.tif 000010
【００３６】
に更新している。第２のステップは、マイクログリッド制御システムｋがその近隣の中間推定値同士を組合せてその更新推定値ωｋ，ｉを得る集約ステップである。
【００３７】
グローバル状態推定は、エージェントマイクログリッド制御システムｋによって計算されるグローバル状態推定値が、式セット（４）を用いて反復ごとに更新される反復プロセスにおいて実行される。この場合、
【００３８】
【数９】
2022501994.tif 000011
【００３９】
は、反復ｉでのエージェントマイクログリッド制御システムｋの状態を示しており、ωｋ，ｉは、反復ｉでのエージェントｋに関する中間変数であり、μｋは、エージェントマイクログリッド制御システムｋの非負更新パラメータであり、∇ωＪｌ（ωｋ，ｉ−１）は、反復ｉでの中間状態ωのエージェントマイクログリッド制御システムｋに関する確率勾配である。
【００４０】
【数１０】
2022501994.tif 000012
【００４１】
特定の実施形態では、μｋは０．２であり、νｋは．０１である。他の実施形態では、μｋは０〜１の範囲内であり、νｋは０〜０．１の範囲内である。確率勾配は、ある反復から他の反復までのωの差である。各マイクログリッド制御システムは、その勾配を追跡し、それに応じてそのすべての近隣の組合わせを更新する。
【００４２】
各々の反復中、マイクログリッド制御システムは、それ自体の更新済みグローバル状態推定値を他のマイクログリッド制御システムに送信するとともに、他のマイクログリッド制御システムの各々から更新済みグローバル状態推定値を受信する。マイクログリッド制御システムは、それ自体の更新済みグローバル状態推定値を他のマイクログリッド制御システムによって計算される更新済みグローバル状態推定値と比較して、すべてのグローバル状態推定値が収束したかどうか、すなわち、グローバル状態ベクトルの標準偏差が閾値を下回っているかどうかを判断する。たとえば、閾値は、いくつか例を挙げると、σ＜０．００１またはσ＜０．１であり得る。動作２０８は、すべてのグローバル状態推定値が収束してしまうまで、またはすべてのグローバル状態推定値が収束していないが規定された計算期間が終了するまで、続けられる。規定された計算期間は、いくつか例を挙げると、一定回数の反復または一定量の反復であり得る。
【００４３】
プロセス２００が条件２０９に進む。条件２０９で、マイクログリッド制御システムは、ネットワーク化されたマイクログリッドシステムの全マイクログリッド制御システムによって計算された全グローバル状態推定値が収束したかどうかを判断する。
【００４４】
プロセス２００は、全グローバル状態推定値が収束した場合、条件２０９から動作２１５に進む。動作２１５で、他のアプリケーションによる使用のためにグローバル状態推定値が利用可能にされる。たとえば、グローバル状態推定値は、いくつか例を挙げると、市場シミュレーション、安定性分析、および需要応答のために用いられてもよい。
【００４５】
グローバル状態推定値が収束していない場合、プロセス２００は条件２０９から動作２１１に進む。動作２１１で、マイクログリッド制御システムは、通常の推定モードから偽データ軽減モードに移行する。通常の条件下では、グローバル状態推定は、最良の収束精度結果を達成するための最適なポリシーとして静的信頼行列を適用する。グローバル状態推定値が収束していないとき、マイクログリッド制御システムは、動的信頼行列を含む修正済みグローバル状態推定式セットを用いて、損なわれたマイクログリッド制御システムを検出する。動的信頼行列を用いることで、マイクログリッド制御システムは、複数の近隣の信頼因数を低減させるとともに、信頼因数が閾値未満に低減した場合にリンクを完全に切断し得る。近隣のグローバル状態推定値が、当該近隣のグローバル状態推定値のノイズのせいで他のグローバル状態推定値とクラスタ化できない場合、近隣の信頼因数が低下する。
【００４６】
偽データ軽減モードでは、損なわれたマイクログリッド制御システムを示すノイズが、受信したグローバル状態推定値において検出されるかどうかに応じて、信頼因数に対して低減された重みが割当てられる。損なわれたマイクログリッド制御システムを検出するために、式セット（３）および（４）の代わりに以下の式セット（５）が反復的に計算される。ここで、ｋはエージェントマイクログリッド制御システムであり、ｌは近隣マイクログリッド制御システムであり、
【００４７】
【数１１】
2022501994.tif 000013
【００４８】
は反復ｉでのエージェントｋの状態を示し、ωｋ，ｉは反復ｉでのエージェントｋに関する中間変数であり、μｋはエージェントｋの非負更新パラメータであり、∇ωＪｌ（ωｋ，ｉ−１）は、反復ｉ−１での中間状態ωのエージェントｋに関する確率勾配であり、νｋは忘却因数であり
【００４９】
【数１２】
2022501994.tif 000014
【００５０】
は瞬間誤差メトリックである。
【００５１】
【数１３】
2022501994.tif 000015
【００５２】
特定の実施形態では、μｋは０．２であり、νｋは．０１である。他の実施形態では、μｋは０〜１の範囲内であり、νｋは０〜０．１の範囲内である。確率勾配は、ある反復から他の反復までのωの差である。各マイクログリッド制御システムは、その勾配を追跡し、それに応じてそのすべての近隣の組合わせを更新する。
【００５３】
プロセス２００は、動作２１１から動作２１３に進む。ここで、プロセス２００の残りの部分に関して、各マイクログリッド制御システムによるさらなる考慮要素から偽データのソースが除去される。たとえば、各マイクログリッド制御システムは、偽データに関連付けられたグローバル状態推定値を除去することによって、そのグローバル状態推定ベクトルを更新してもよい。いくつかの実施形態では、１つ以上のマイクログリッド制御システムは、将来のグローバル状態推定値から偽データソースを除去するために、再び動作２０３、２０５および２０７を実行する。いくつかの実施形態では、マイクログリッド制御システムは、アラームをシステムオペレータに送信してもよく、または、状態推定のためにデータがデータソースから収集されるのを阻止してもよい。
【００５４】
プロセス２００が動作２１３から動作２０８に進む。ここで、マイクログリッド制御システムは、動作２１１において除去されたグローバル状態推定値なしで、式セット（４）を用いるグローバル状態推定および式セット（３）を用いて計算される静的信頼行列の反復手順を続ける。動作２０８は、残りのグローバル状態推定値が収束するまで続く。次いで、プロセス２００は条件２０９、動作２１５に進み、最後に動作２１７で終了する。
【００５５】
図３を参照すると、プロセス２００に従ってシステム１００によって実行される、グローバル状態推定値の収束を例示する複数のグラフ３００が示されている。システム１００の各マイクログリッド制御システムは、それ自体のマイクログリッドのフィールドデバイスからのローカル測定値を用いて、ローカル状態推定値を計算する。ベクトルセット（６）は、すべてのマイクログリッドに関するローカル状態推定ベクトルを含み、各ベクトルは、マイクログリッド内の各バスに関する初期電圧および位相角
【００５６】
【数１４】
2022501994.tif 000016
【００５７】
の推定値を含む。０回目の反復で、各マイクログリッドは、初期推定ベクトルｘＭＧｎ（０）を設定する。
【００５８】
【数１５】
2022501994.tif 000017
【００５９】
各マイクログリッド制御システムが他のマイクログリッド制御システムの各々からローカル状態推定ベクトルを受信すると、各マイクログリッド制御システムはベクトルセット（６）を用いてグローバル推定ベクトルｘＭＧｇｌｏｂａｌを生成する。
【００６０】
【数１６】
2022501994.tif 000018
【００６１】
式セット（３）を用いて、各マイクログリッド制御システムは信頼行列を計算する。以下の行列は、式セット（３）を用いて信頼行列に変換されているシステム１００のトポロジ行列を示す。
【００６２】
【数１７】
2022501994.tif 000019
【００６３】
各マイクログリッド制御システムは、信頼行列（８）および式セット（４）を用いて、図２の動作に関して上述したようにグローバル状態推定値を反復的に計算する。グラフ３１０は、各反復ごとにマイクログリッド１１０〜１５０によって計算されるバス１５３の電圧大きさに関するそれぞれのグローバル状態推定値３１１〜３１５を含む。グラフ３２０は、各反復ごとにマイクログリッド１１０〜１５０によって計算されるバス１５３の電圧位相角に関するそれぞれのグローバル状態推定値３２１〜３２５を含む。４０回目の反復において、１×１０−５未満の標準偏差を有する、各マイクログリッド制御システムのグローバル状態推定値が収束する。
【００６４】
図４を参照すると、図２のプロセス２００に従ってシステム１００によって実行されるグローバル状態推定中の偽データ検出および軽減を例示する複数のグラフ４００が示されている。グラフ４１０は、各反復ごとにマイクログリッド１１０〜１５０によって実行されるバス１５３の電圧大きさに関するそれぞれのグローバル状態推定４１１〜４１５の値を含む。グラフ４２０は、各反復ごとにマイクログリッド１１０〜１５０によって実行されるバス１５３の電圧大きさに関するそれぞれのグローバル状態推定４１１〜４１５の値を含む。マイクログリッド制御システム１４１によって計算されて推定値４１４および４２４として示されるグローバル状態推定値は、グローバル状態推定値に１０％バイアス因数を導入する攻撃者によって損なわれたものであって、これにより、近隣エージェントに不正確な推定値が送られることとなった。
【００６５】
いくつかの実施形態では、偽データは、１つのマイクログリッドからのグローバル状態推定値の一部または全てを損なうように、グローバル状態推定ベクトル内の１つ以上のグローバル状態推定値に適用され得る。いくつかの実施形態では、偽データは、２つ以上のマイクログリッド制御システムからグローバル状態推定値に適用され得る。
【００６６】
０回目の反復で、マイクログリッド制御システムの動作は通常推定モードである。１０回目の反復から始まって、１０％バイアス因数がマイクログリッド制御システム１４１のグローバル推定値４１４、４２４に追加される。サイバー攻撃が開始されると、マイクログリッド制御システムのグローバル状態推定値が互いから逸脱し始める。マイクログリッド制御システムは、３０回目の反復まで、通常の推定モードで動作し続ける。
【００６７】
３０回目の反復で、マイクログリッド制御システムは、グローバル状態推定値が収束していないと判断し、それに応じて偽データ軽減モードに進む。３０回目の反復の後、偽データなしの推定値は、最初にそれらの信頼重みに基づいてクラスタを形成する傾向がある。グラフ４１０に関して、マイクログリッド制御システム１４１の推定値４１４は、マイクログリッド制御システム１１１および１２１の推定値４１１および４１２ならびにマイクログリッド制御システム１３１および１５１の推定値４１３および４１５によって形成される２つのクラスタから大きく逸脱している。反復が進むにつれて、推定値４１１、４１２、４１３および４１５は互いに向かって収束するが、損なわれた推定値４１４は逸れている。
【００６８】
反復を繰り返した後の信頼行列の適応が行列（９）〜（１２）によって示される。行列（９）は、反復３１での動的信頼行列であり、行列（１０）は、反復３２での動的信頼行列であり、行列（１１）は、反復４０での動的信頼行列であり、行列（１２）は、反復５０での動的信頼行列である。３０回目の反復から反復が進むにつれて、マイクログリッド制御システム１４１から受信したグローバル状態推定値に対するマイクログリッド制御システム１３１の信頼因数である行列Ａ（４，３）が減少する。いくつかの実施形態では、信頼因数が所定の閾値よりも減少すると分離が発生するだろう。たとえば、予め定義された閾値は、いくつか例を挙げると、ａｋｌ＜０．００００１、ａｋｌ＜．００５、または、ａｋｌ＜．０５となり得る。
【００６９】
【数１８】
2022501994.tif 000020
【００７０】
図５を参照すると、図１のシステム１００のための例示的なマイクログリッド制御システムおよび位相角発見プロセスを例示する表が示されている。マイクログリッド制御システム間の通信のたびに、１つのマイクログリッド制御システムから別のマイクログリッド制御システムに送られるすべてのメッセージが限られた時間内に破損なしで到着するが、必ずしも同じ順序で公開されるわけではない。さらに、各マイクログリッド制御システムは、その近隣と、受信したメッセージの発信元である近隣のアイデンティティとを認識している。
【００７１】
一例では、マイクログリッド１１０の共通結合の格子点は、位相角計算のための基準バスとなるように決定される。マイクログリッド制御システムは、それ自体と近隣マイクログリッドとの間で位相角差を決定するための情報をやり取りすることで、各マイクログリッドが基準バスに対するそれ自体の位相角変位を計算できるようにする。
【００７２】
反復（０）で、各マイクログリッド制御システムはそれ自体の固有のアイデンティティを認識している。０回目の反復後の反復ごとに、各マイクログリッド制御システムは、その近隣からデータを要求することによって情報を取得する。反復（１）で、近隣マイクログリッド制御システムが、エージェントとしても知られるエージェントマイクログリッド制御システムと情報をやり取りし、新しいマイクログリッド制御システムおよび位相角差が発見される。位相角は、直接通信されてもよく、または、やり取りされた電圧および電力流などの、近隣間でやり取りされた情報を用いて計算されてもよい。たとえば、以下の式セット（１）を用いて、２つの隣接するマイクログリッド間の位相角を決定してもよい。ここで、ｉは１つのマイクログリッドの相互接続バスであり、ｊは隣接するマイクログリッドの相互接続バスであり、θｉ−ｊはバスｉとバスｊとの間の位相角差であり、Ｐｉｊはバスｉとバスｊとの間のライン上を流れる有効電力であり、Ｖｉはバスｉにおける電圧であり、Ｖｊはバスｊにおける電圧であり、Ｒｉｊはバスｉとバスｊとの間のライン抵抗であり、Ｘｉｊはバスｉとバスｊとの間のラインリアクタンスである。
【００７３】
【数１９】
2022501994.tif 000021
【００７４】
反復（２）で、エージェントマイクログリッド制御システムは、反復（１）中にその近隣によって事前に取得されていた情報、すなわち近隣の近隣からの情報、を取得する。当該プロセスは、マイクログリッド制御システムがいずれも新しいマイクログリッド制御システムを発見しなくなるまで続く。図５に示されるように、発見プロセスは、ネットワークの直径Ψである、４回目の反復の終了時に完了する。ネットワーク内のマイクログリッド制御システムの数は、Ｎ＝Ψ＋１として決定される。分散型マイクログリッド制御システム発見は、各推定プロセスの後に適用されて、次回の状態推定に何個のマイクログリッドが関与するかを決定する。マイクログリッド制御システム発見の終了時に、すべてのマイクログリッド制御システムは、θ１−２、θ１−３、θ３−４、θ３−５についての情報を認識するようになる。決定された位相角差を用いて、各マイクログリッド制御システムは、マイクログリッド１１０の基準バスに対するそれ自体の位相角差を計算する。たとえば、マイクログリッド制御システム１５１は、θ１−３−θ３−５としてθ１−５を計算する。マイクログリッド１１０以外のシステム１００内の各マイクログリッドは、ローカル状態推定の位相角結果をシフトさせるために、この位相角差を用いるだろう。
【００７５】
図６Ａおよび図６Ｂを参照すると、例示的なマイクログリッド状態推定プロセスと、コンセンサスアルゴリズムを用いる代替的なマイクログリッド状態推定プロセスとの計算性能を例示するグラフが示されている。２つのプロセスは、２．７０ＧＨｚでクロックするＣｏｒｅ  ｉ７と８ＧＢ  ＲＡＭとを含むコンピュータにおいて実装されている。各プロセスの計算性能は、３０マイクログリッドの３００バスネットワークにおいて分析される。図６Ａのグラフ６１０は図２のプロセス２００の通常動作モードの実行を示す。図６Ｂのグラフ６２０は、コンセンサスアルゴリズムを用いた別のマイクログリッド分散型状態推定プロセスの実行を示す。グラフ６２０の複数の推定値６１２は、１００回反復しても、４．９×１０−３の標準偏差を有する正確な推定値に収束しない。しかしながら、グラフ６１０の複数の推定値６１１は、１．０４６１×１０−５の標準偏差を有する正確な推定値に収束する。
【００７６】
さらに、図７の表に示されるように、例示的なグローバル状態推定プロセスのための収束時間は、同じグローバル状態推定値を計算するための従来の集中型プロセスと比べて著しく短い。２つのプロセスは、２．７０ＧＨｚでクロックするＣｏｒｅ  ｉ７と８ＧＢのＲＡＭとを含むコンピュータにおいて実装されていた。各プロセスの計算性能は、３０マイクログリッドの３００バスネットワークにおいて分析される。集中型状態推定に必要な反復はより少ないものの、集中型状態推定において用いられる行列がより大きいので各反復の経過時間は著しくより長くなる。
【００７７】
図８を参照すると、図１のシステム１００のマイクログリッド制御システムのうちの１つなどの例示的なマイクログリッド制御システム８００の概略ブロック図が例示されている。マイクログリッド制御システム８００は、処理デバイス８０２と、入出力デバイス８０４と、メモリデバイス８０６と、動作論理８０８とを含む。さらに、コンピューティングデバイス８００は、ネストされたマイクログリッドシステムまたはローカルフィールドデバイスにおける他のマイクログリッド制御システムを含む１つ以上の外部デバイス８１０と通信する。マイクログリッド制御システム８００は、独立型デバイス、埋込み型システム、またはシステム８００に関して説明した機能を実行するように構成された複数のデバイスであってもよい。たとえば、システム１５１は、マイクログリッドＳＣＡＤＡゲートウェイに組込まれてもよい。
【００７８】
入出力デバイス８０４は、マイクログリッド制御システム８００がローカルフィールドデバイスまたは他のマイクログリッドと通信することを可能にする。入出力デバイス８０４は、いくつか例を挙げると、ネットワークアダプタ、ネットワーククレデンシャル、インターフェイス、またはポート（たとえば、ＵＳＢポート、シリアルポート、パラレルポート、アナログポート、デジタルポート、ＶＧＡ、ＤＶＩ、ＨＤＭＩ（登録商標）、ＦｉｒｅＷｉｒｅ、ＣＡＴ５、イーサネット（登録商標）、ファイバ、または任意の他の種類のポートもしくはインターフェイス）を含み得る。入出力デバイス８０４は、これらのアダプタ、クレデンシャル、またはポート、たとえば、データを受信するための第１のポートおよびデータを送信するための第２のポートなどのうち２つ以上を含み得る。
【００７９】
処理デバイス８０２は、いくつか例を挙げると、１つまたは複数のプロセッサ、演算論理装置（arithmetic-logic unit：ＡＬＵ）、中央処理装置（central processing unit：ＣＰＵ）、デジタル信号プロセッサ（digital signal processor：ＤＳＰ）、またはフィールドプログラマブルゲートアレイ（field-programmable gate array：ＦＰＧＡ）を含み得る。複数の処理ユニットを有する処理デバイスの形態に関して、分散処理、パイプライン処理、または並列処理が用いられ得る。処理デバイス８０２は、本明細書で説明される動作のみを実行するためだけのものであってもよく、または１つ以上の追加のアプリケーションにおいて用いられてもよい。処理デバイス８０２は、メモリ８０６に格納されたプログラミング命令（ソフトウェアまたはファームウェアなど）によって定義されるように、動作論理８０８に従ってアルゴリズムを実行するとともにデータを処理するプログラム可能な種類のものであってもよい。代替的にまたは付加的には、処理デバイス１２０２のための動作論理８０８は、ハードワイヤード論理または他のハードウェアによって少なくとも部分的に定義されている。処理デバイス８０２は、入出力デバイス８０４または他の箇所から受信した信号を処理するとともに所望の出力信号を提供するのに適した任意の種類の１つ以上の構成要素を備え得る。このような構成要素は、デジタル回路、アナログ回路、またはこれら両方の組合わせを含み得る。
【００８０】
メモリデバイス８０６は、コンピュータ可読媒体としても公知であり、いくつか例を挙げると、ソリッドステート型、電磁型、光学型、またはこれらの形態の組合わせなどの１つ以上の種類のメモリであり得る。さらに、メモリデバイス８０６は、揮発性、不揮発性、一時的、非一時的、またはこれらの種類の組合わせであってもよく、メモリデバイス８０６の一部または全体は、いくつか例を挙げると、ディスク、テープ、メモリスティック、またはカートリッジなどの携帯型であってもよい。加えて、メモリデバイス８０６は、一例を挙げると、動作論理８０８を定義するプログラミング命令を格納することに加えて、またはその代わりに、入出力デバイス８０４から受信した信号および／または入出力デバイス８０４に送信される信号を表わすデータなどの、処理デバイス８０２の動作論理８０８によって操作されるデータを格納していてもよい。メモリデバイス８０６は、処理デバイス８０２に含まれてもよく、および／または、処理デバイス８０２に結合されてもよい。
【００８１】
いくつかの例示的な実施形態をここでさらに説明する。一実施形態は、第１のマイクログリッドおよび複数のマイクログリッドを含むネットワーク化されたマイクログリッドシステムから偽データを除去するための方法であって、第１の複数のローカル測定値を用いてマイクログリッド制御システムで第１のローカル状態推定値を計算するステップと、マイクログリッド制御システムで、複数のマイクログリッドから複数のローカル状態推定値を受信するステップと、第１のローカル状態推定値および複数のローカル状態推定値を用いて、マイクログリッド制御システムおよび複数のマイクログリッドで複数のグローバル状態推定値を計算するステップと、複数のグローバル状態推定値が収束しなかったと判断するステップと、マイクログリッド制御システムおよび複数のマイクログリッドを用いて、マイクログリッド制御システムで計算された複数のグローバル状態推定値のうちの第１のグローバル状態推定値における偽データを検出するステップと、複数のグローバル状態推定値から第１のグローバル状態推定値を除去するステップと、複数のグローバル状態推定値のうち残りのグローバル状態推定値が収束するまで、マイクログリッド制御システムおよび複数のマイクログリッドで、複数のグローバル状態推定値のうち残りのグローバル状態推定値を反復的に更新するステップとを含む。
【００８２】
上述の方法のいくつかの例では、当該方法は、第１のマイクログリッド制御システムで、複数のマイクログリッドの第１の部分に関連する直接通信の際に当該第１の部分から情報を受信し、その後、第１のマイクログリッド制御システムとの間接通信の際に、当該第１の部分から複数のマイクログリッドの第２の部分に関連する情報を受信することによって、通信ネットワークトポロジを決定するステップを含む。いくつかの例では、グローバル状態推定値は、各グローバル状態推定値が閾値未満の標準偏差を有するときに収束する。いくつかの例では、複数のグローバル状態推定値を計算するステップは、複数回の反復を用いて行われ、複数のグローバル状態推定値が収束しなかったと判断するステップは、複数回の反復の後に行なわれる。いくつかの例では、第１のグローバル状態推定値における偽データを検出するステップは、複数のマイクログリッドのうち１つのマイクログリッドから第１のマイクログリッド制御システムによって受信されるグローバル状態推定値に適用される重みに対応する第１の信頼因数と、第１のマイクログリッド制御システムによって計算されるグローバル状態推定値に適用される重みに対応する第２の信頼因数とを含む複数の信頼因数を含む動的信頼行列を更新するステップを含む。いくつかの例では、第１のグローバル推定値における偽データを検出するステップは、式セット（５）を用いて複数のグローバル状態推定値を更新するステップを含み、ここで、ｋはエージェントマイクログリッド制御システムであり、ｌは近隣マイクログリッド制御システムであり、
【００８３】
【数２０】
2022501994.tif 000022
【００８４】
は、反復ｉでのエージェントマイクログリッド制御システムｋのグローバル状態推定値を示し、ωｋ，ｉは、反復ｉでのエージェントマイクログリッド制御システムｋに関する中間変数であり、μｋは、エージェントマイクログリッド制御システムｋの非負更新パラメータであり、∇ωＪｌ（ωｋ，ｉ−１）は、反復ｉ−１での中間状態ωのエージェントマイクログリッド制御システムｋに関する確率勾配であり、νｋは忘却因数であり、
【００８５】
【数２１】
2022501994.tif 000023
【００８６】
は瞬間誤差メトリックである。いくつかの例では、複数のグローバル状態推定値の各グローバル状態推定値は、第１のマイクログリッドのバスの電圧大きさに対応する。
【００８７】
例示的な実施形態は、マイクログリッドシステムにおける偽データを検出および軽減するための方法であって、第１のマイクログリッドの第１のマイクログリッド制御システム、第２のマイクログリッドの第２のマイクログリッド制御システム、および第３のマイクログリッドの第３のマイクログリッド制御システムを動作させるステップと、第１の複数の測定値を用いて第１のマイクログリッド制御システムで第１のローカル状態推定値を計算し、第２の複数の測定値を用いて第２のマイクログリッド制御システムで第２のローカル状態推定値を計算し、第３の複数の測定値を用いて第３のマイクログリッド制御システムで第３のローカル状態推定値を計算するステップと、第２のマイクログリッド制御システムおよび第３のマイクログリッド制御システムで第１のローカル状態推定値を受信し、第１のマイクログリッド制御システムおよび第３のマイクログリッド制御システムで第２のローカル状態推定値を受信し、第１のマイクログリッド制御システムおよび第２のマイクログリッド制御システムで第３のローカル状態推定値を受信するステップと、第１のマイクログリッド制御システムで第１のマイクログリッドの電気的特性の第１のグローバル状態推定値を計算し、第２のマイクログリッド制御システムで電気的特性の第２のグローバル状態推定値を計算し、第３のマイクログリッド制御システムで電気的特性の第３のグローバル状態推定値を計算するステップと、第１のマイクログリッド制御システムで、当該第１のグローバル状態推定値、当該第２のグローバル状態推定値、および当該第３のグローバル状態推定値が収束しないと判断するステップと、動的信頼行列を用いて、当該第１のマイクログリッド制御システムで第１のグローバル状態推定値を反復的に更新し、第２のマイクログリッド制御システムで第２のグローバル状態推定値を反復的に更新し、第３のマイクログリッド制御システムで第３のグローバル状態推定値を反復的に更新するステップと、当該第１のグローバル状態推定値、当該第２のグローバル状態推定値、および当該第３のグローバル状態推定値を反復的に更新するステップを用いて、当該第１のグローバル状態推定値内の偽データを検出するステップと、当該第２のグローバル状態推定値および当該第３のグローバル状態推定値が収束するまで、当該第２のマイクログリッド制御システムでの当該第２のグローバル状態推定値の更新および当該第３のマイクログリッド制御システムでの当該第３のグローバル状態推定値の更新を続けるステップとを含む。
【００８８】
上記の方法のいくつかの例では、当該方法は、第１の反復中に、第１のマイクログリッド制御システムで、第２のマイクログリッドに関する情報を第２のマイクログリッド制御システムから受信し、第１の反復中に、第２のマイクログリッド制御システムで、第３のマイクログリッドに関する情報を第３のマイクログリッド制御システムから受信し、第２の反復中に、第１のマイクログリッド制御システムで、第３のマイクログリッドに関する情報を第２のマイクログリッド制御システムから受信することによって、通信ネットワークトポロジを決定するステップを含む。いくつかの例では、第１のローカル状態推定値、第１のグローバル状態推定値、第２のグローバル状態推定値、および第３のグローバル状態推定値は、第１のマイクログリッドのバスの推定電圧である。いくつかの例では、第１のマイクログリッド制御システムは、第２の複数のローカル測定値または第３の複数のローカル測定値を受信しない。いくつかの例では、第１のグローバル状態推定値、第２のグローバル状態推定値、および第３のグローバル状態推定値を計算するステップは、拡散アルゴリズムを用いるステップを含み、第２のグローバル状態推定値および第３のグローバル状態推定値の更新を続けるステップは、拡散アルゴリズムを用いるステップを含む。いくつかの例では、動的信頼行列は、第１のマイクログリッド制御システムによって第１のグローバル状態推定値に適用される重みに対応する第１の信頼因数と、第１のマイクログリッド制御システムによって第２のグローバル状態推定値に適用される重みに対応する第２の信頼因数とを含む複数の信頼因数を含む。いくつかの例では、動的信頼行列を用いて第１のグローバル状態推定値、第２のグローバル状態推定値、および第３のグローバル状態推定値を反復的に更新するステップは、第１のマイクログリッド制御システムによって第２のグローバル状態推定値内のノイズを検出することに応答して第２の信頼因数をゼロに低減するステップを含む。いくつかの例では、第１のグローバル状態推定値、第２のグローバル状態推定値、および第３のグローバル状態推定値を反復的に更新するステップは、式セット（５）のうち以下の部分を用いて動的信頼行列を更新するステップを含む。
【００８９】
【数２２】
2022501994.tif 000024
【００９０】
ここで、ｋは、エージェントマイクログリッド制御システムであり、ｌは近隣マイクログリッド制御システムであり、
【００９１】
【数２３】
2022501994.tif 000025
【００９２】
は、反復ｉでのエージェントマイクログリッド制御システムｋのグローバル状態推定値であり、ωｋ，ｉは、反復ｉでのエージェントマイクログリッド制御システムｋに関する中間変数であり、μｋは、エージェントマイクログリッド制御システムｋの非負更新パラメータであり、νｋは忘却因数であり、
【００９３】
【数２４】
2022501994.tif 000026
【００９４】
は瞬間誤差メトリックである。
例示的な実施形態は非一時的メモリデバイスである。当該非一時的メモリデバイスは、第１のマイクログリッドの第１のマイクログリッド制御システムの処理デバイスによって実行されるように構成された命令のセットを含み、当該命令のセットは、第１の複数の測定値を用いて第１のローカル状態推定値を計算し、第１のマイクログリッドに結合された複数のネットワーク化されたマイクログリッドから複数のローカル状態推定値を受信し、第１のローカル状態推定値および複数のローカル状態推定値を用いて第１のグローバル状態推定値を計算し、複数のネットワーク化されたマイクログリッドから複数のグローバル状態推定値を受信し、第１のグローバル状態推定値および複数のグローバル状態推定値が収束しないと判断し、第１のグローバル状態推定値および複数のグローバル状態推定値における検出されたノイズの大きさに応じて更新されるように構成された動的信頼行列を用いて、第１のグローバル状態推定値および複数のグローバル状態推定値を反復的に更新することによって、複数のグローバル状態推定値のうち少なくとも１つのグローバル状態推定値が偽データを含むと判断し、複数のグローバル状態推定値のうち偽データを含むと判断されなかったグローバル状態推定値で第１のグローバル状態推定値が収束するまで、第１のグローバル状態推定値を反復的に更新する、のに有効である。
【００９５】
前述の非一時的メモリデバイスのいくつかの例では、複数のフィールドデバイスは、第１のマイクログリッドの複数のバス電圧を測定するように構成された複数の電圧センサを含む。いくつかの例では、処理デバイスによって実行されるように構成された命令のセットは、複数のネットワーク化されたマイクログリッドの第１の部分に関連する直接通信の際に当該第１の部分から情報を受信し、その後、第１のマイクログリッド制御システムとの間接通信の際に、第１の部分から、複数のネットワーク化されたマイクログリッドの第２の部分に関連する情報を受信することによって、通信ネットワークトポロジを決定するのに有効である。いくつかの例では、動的信頼行列は、複数のマイクログリッドのうちの１つのマイクログリッドから第１のマイクログリッド制御システムによって受信されるグローバル状態推定値に適用される重みに対応する第１の信頼因数と、第１のグローバル状態推定値に適用される重みに対応する第２の信頼因数とを含む複数の信頼因数を含む。いくつかの例では、複数のグローバル状態推定値のうち少なくとも１つのグローバル状態推定値が偽データを含むと判断することは、拡散アルゴリズムを用いることと、式セット（５）のうち以下の部分を用いて動的信頼行列を更新することとを含む。
【００９６】
【数２５】
2022501994.tif 000027
【００９７】
ここで、ｋはエージェントマイクログリッド制御システムであり、ｌは近隣マイクログリッド制御システムであり、
【００９８】
【数２６】
2022501994.tif 000028
【００９９】
は、反復ｉでのエージェントマイクログリッド制御システムｋに関するグローバル状態推定値を示し、ωｋ，ｉは、反復ｉでのエージェントマイクログリッド制御システムｋに関する中間変数であり、μｋは、エージェントマイクログリッド制御システムｋの非負更新パラメータであり、νｋは忘却因数であり、
【０１００】
【数２７】
2022501994.tif 000029
【０１０１】
は瞬間誤差メトリックである。
さまざまな実施形態からのさまざまな局面、特徴、プロセスおよび動作は、そうではないと明示的に述べられない限り、他の実施形態のいずれにも使用され得ることが企図されている。例示される特定の動作は、コンピュータが非一時的なコンピュータ可読記憶媒体上でコンピュータプログラム製品を実行することによって実装されてもよく、この場合、コンピュータプログラム製品は、コンピュータに動作のうちの１つ以上を実行させるか、または、当該コンピュータに、１つ以上の動作を実行するためのコマンドを他のデバイスに発行させる命令を含む。
【０１０２】
本開示の非限定的な例示的実施形態と、その作製および使用方法およびプロセスとを明確、簡潔かつ正確に説明するために、さらに、その実施、作製および使用を可能にするために、図に示されるものを含む特定の例示的実施形態が参照されており、これを説明するのに特定の言語が用いられている。それにも関わらず、本開示の範囲の限定がこれによってなされるものではなく、本開示が、本開示の利益を伴って当業者に想起され得るように例示的な実施形態のこのような変更例、変形例およびさらなる適用例を包含および保護することが理解されるはずである。
【０１０３】
本開示が添付の図面および上述の記載において詳細に図示および説明されてきたが、本開示は、特徴を限定するのではなく例示するものと見なされるべきであり、特定の例示的実施形態のみが図示および説明されていること、ならびに、本開示の精神の範囲内に収まるすべての変更例および変形例が保護されることが望ましいことを理解されたい。上述の記載において利用される「好まれる」、「好ましくは」、「好ましい」、または「より好ましい」などの語の使用は、そのように記載される特徴がより望ましい可能性があることを示しているが、必ずしも必要でなくてもよく、それが欠如している実施形態も、添付の特許請求の範囲によって定義されている本開示の範囲内となるよう企図されていると理解されるべきである。請求項を読む際に、「a」、「an」といった不定冠詞、「少なくとも１つ（at least one）」、または「少なくとも１つの部分（at least one portion）」などの単語が用いられる場合、請求項において特段の記載がない限り、請求項を１つのアイテムだけに限定する意図がないことが意図されている。前置詞「の（of）」は、それが用いられる文脈によって示されるように、別のアイテムとの関連付けまたは別のアイテムへのつながり、さらに、他のアイテムへの帰属または他のアイテムとのつながりを意味し得る。語句「に結合される（coupled to）」、「と結合される（coupled with）」などは、間接的な接続および結合を含んでおり、さらに、そうでないと明記されない限り、直接的な結合または接続を、必ずしも必要とはしないが含んでいる。「少なくとも１つの部分（at least a portion）」および／または「１つの部分（a portion）」という語が用いられる場合、そうでないと明記されない限り、アイテムは１つの部分および／またはアイテム全体を含み得る。
【図１】
2022501994.tif 000039
【図２】
2022501994.tif 000040
【図３】
2022501994.tif 000041
【図４】
2022501994.tif 000042
【図５】
2022501994.tif 000043
【図６Ａ】
2022501994.tif 000044
【図６Ｂ】
2022501994.tif 000045
【図７】
2022501994.tif 000046
【図８】
2022501994.tif 000047
【国際調査報告】
2022501994.tif 000048