(57)【要約】  多相配電ネットワーク上の地絡故障に関係する相を決定するための装置（１００）であって、上記装置は、スパイク故障を判断することに基づき、相間故障電流の複数のＲＭＳ値のうちのどれが最小値であるかを判断するように構成され、上記最小ＲＭＳ値が、上記地絡故障に関係する上記多相配電ネットワークの特定の相を表している。【選択図】図１
2019507344.tif 000002
【特許請求の範囲】
【請求項１】
多相配電ネットワーク上の地絡故障に関係する相を決定するように構成された装置（１００）であり、前記装置は、スパイク故障を判断することに基づき、相間故障電流の複数のＲＭＳ値のうちのどれが最小値であるかを判断するように構成され、前記最小ＲＭＳ値が、前記地絡故障に関係する前記多相配電ネットワークの特定の相を表している、装置（１００）であって、
前記装置（１００）が、所定のサンプリングレートで求めた複数のサンプルセットおよび残留電流値、ｉＲ、に基づき前記スパイク故障の発生を判断するように構成され、前記サンプルセットが、特定のサンプリング時刻における、前記多相配電ネットワークの各相の相電流を含み、前記残留電流値、ｉＲ、が前記相電流の総和を含み、前記装置が次式のように前記残留電流のＲＭＳ値、ＩＲ（ｎ）、を決定するように構成され、
【数１】
2019507344.tif 000014
ここで、ｋが前記総和のインデックス値であり、ｎがスライディング窓の最新のサンプルを示し、Ｎがサンプルの数に関する前記スライディング窓サイズであり、前記スパイク故障の前記発生は、少なくともＩＲ（ｎ）が所定のスパイク時間しきい値（ＩＲｔｉｍｅ）よりも短い間所定の残留電流しきい値（ＩＲｓｅｔ）よりも大きいことに基づいて判断され、
前記装置（１００）は、各サンプルセットの前記サンプルのそれぞれの相電流の間の差違を含む相間電流ｉｘｙ（ｎ）および前記スライディング窓の始まりと終わりとにおける各相についての前記相間電流の間の差違を含む相間故障電流Δｉｘｙ（ｎ）に基づいて、前記相間故障電流Δｉｘｙ（ｎ）の複数のＲＭＳ値のうちのどれが最小値であるかを判断するように構成され、前記装置が次式のように前記相間故障電流の前記ＲＭＳ値ΔＩｘｙ（ｎ）を決定するように構成され、
【数２】
2019507344.tif 000015
ここで、ｘおよびｙが各相間電流を決定するために使用した２つの相を示す、
装置（１００）。
【請求項２】
前記装置（１００）が、相ａ、ｂおよびｃを含む３相配電ネットワークで使用されるように構成され、最小値である前記相間故障電流の前記ＲＭＳ値は、前記最小値ΔＩｍｉｎが、
ｉｖ．  ΔＩｂｃから成る場合には、前記地絡故障が相ａに関係付けられ、
ｖ．  ΔＩｃａから成る場合には、前記地絡故障が相ｂに関係付けられ、
ｖｉ．  ΔＩａｂから成る場合には、前記地絡故障が相ｃに関係付けられる
ことで地絡故障が存在する前記多相配電ネットワークの特定の相を表している、請求項１記載の装置。
【請求項３】
前記装置（１００）は、前記スパイク故障の前記発生を判断する際に、前記残留電流の前記ＲＭＳ値の変化率、および前記変化率が変化率のしきい値（２１４）よりも大きいか否かを決定するように構成され、
前記スパイク故障の前記発生は、少なくともＩＲ（ｎ）が前記所定のスパイク時間しきい値よりも短い間前記所定の残留電流しきい値（１１６）よりも大きいことおよびＩＲ（ｎ）が前記所定の残留電流しきい値（１１６）よりも大きい前記時間の少なくとも一部の間に前記変化率が前記変化率のしきい値（２１４）よりも大きいことに基づいて判断される、請求項１または２記載の装置。
【請求項４】
前記装置（１００）は、前記故障が、
複数の相の各々の組合せについて前記相間故障電流の前記ＲＭＳ値ΔＩｘｙ（ｎ）および前記残留故障電流のＲＭＳ値ΔＩＲ（ｎ）に基づいて順故障方向のものであるかどうかを判断するように構成され、前記残留故障電流が前記スライディング窓の前記始まりにおいて前記サンプルセットのうちの１つから導き出された前記残留電流と前記スライディング窓の前記終わりにおいて前記サンプルセットのうちの１つから導き出されたものとの間の差違を含み、前記残留故障電流の前記ＲＭＳ値ΔＩＲ（ｎ）が
【数３】
2019507344.tif 000016
から成り、
前記装置が、前記相間故障電流の前記ＲＭＳ値の最小値ΔＩｍｉｎおよび前記相間故障電流の前記ＲＭＳ値の最大値ΔＩｍａｘを特定するように構成され、前記装置は、
ΔＩｍｉｎが前記残留故障電流の前記ＲＭＳ値の第１の所定の割合未満であり、
ΔＩｍａｘが前記残留故障電流の前記ＲＭＳ値の第２の所定の割合よりも大きい
場合に順方向故障を判断するように構成される、請求項１乃至３のいずれか１項記載の装置。
【請求項５】
前記第１の所定の割合が０．１と０．２との間であり、および／または
前記第２の所定の割合が０．８と０．９との間である、
請求項４記載の装置。
【請求項６】
前記装置（１００）が、地絡故障に関係する相の前記決定に基づいて故障アラームを鳴らすように構成され、故障アラームを前記鳴らすことは、
前記スパイク故障の前記発生および前記スパイク故障中の前記地絡故障に関係する前記多相配電ネットワークの特定の相を前記特定すること、少なくとも所定の、アラーム、所定のアラーム期間内の回数（１１８）に依存する、請求項１乃至５のいずれか１項記載の装置。
【請求項７】
前記故障アラームを前記鳴らすことは、前記スパイク故障の各発生に対して請求項３記載の故障方向を前記決定することにさらに依存する、請求項６記載の装置。
【請求項８】
前記サンプルセットの前記サンプリングレートが２〜３ｋＨｚの範囲内である、請求項１乃至７のいずれか１項記載の装置。
【請求項９】
前記スライディング窓サイズが、前記多相配電ネットワークの前記基本周波数の単一の周期に対応する、請求項１乃至８のいずれか１項記載の装置。
【請求項１０】
前記所定の残留電流しきい値（１１６）が、前記配電ネットワークに印加される定格電流の１０％未満である、請求項１乃至９のいずれか１項記載の装置。
【請求項１１】
前記所定の残留電流しきい値（１１６）が、１００ｍＡ未満であり得る、請求項１乃至１０のいずれか１項記載の装置。
【請求項１２】
前記所定のスパイク時間しきい値が、前記多相配電ネットワークの前記基本周波数の１サイクルに基づく、請求項１乃至１１のいずれか１項記載の装置。
【請求項１３】
前記スライディング窓が、現在の時刻から現在の時刻に先立つ定格周波数の１サイクルの時間に及ぶ、請求項１乃至１２のいずれか１項記載の装置。
【請求項１４】
請求項１乃至１３のいずれか１項記載の地絡故障に関係する相を決定するための前記装置を含む、多相配電ネットワーク。
【請求項１５】
多相配電ネットワーク上の地絡故障に関係する相を決定するための方法であり、前記方法が、スパイク故障を判断することに基づき、相間故障電流の複数のＲＭＳ値のうちのどれが最小値であるかを判断するステップを含み、前記最小ＲＭＳ値が、前記地絡故障に関係する前記多相配電ネットワークの特定の相を表している、方法であって、
所定のサンプリングレートで求めた複数のサンプルセットおよび残留電流値、ｉＲ、に基づき前記スパイク故障の発生を判断するステップであり、前記サンプルセットが、特定のサンプリング時刻において、前記多相配電ネットワークの各相の相電流を含み、前記残留電流値、ｉＲ、が前記相電流の総和を含み、装置が次式のように前記残留電流のＲＭＳ値、ＩＲ（ｎ）、を決定するように構成され、
【数４】
2019507344.tif 000017
ここで、ｋが前記総和のインデックス値であり、ｎがスライディング窓の最新のサンプルを示し、Ｎがサンプルの数に関する前記スライディング窓サイズであり、前記スパイク故障の前記発生は、少なくともＩＲ（ｎ）が所定のスパイク時間しきい値（ＩＲｔｉｍｅ）よりも短い間所定の残留電流しきい値（ＩＲｓｅｔ）よりも大きいことに基づいて判断される、発生を判断するステップと、
各サンプルセットの前記サンプルのそれぞれの相電流の間の差違を含む相間電流ｉｘｙ（ｎ）および前記スライディング窓の始まりと終わりとにおける各相についての前記相間電流の間の差違を含む相間故障電流Δｉｘｙ（ｎ）に基づいて、前記相間故障電流Δｉｘｙ（ｎ）の複数のＲＭＳ値のうちのどれが最小値であるかを判断するステップであり、前記装置が次式のように前記相間故障電流の前記ＲＭＳ値ΔＩｘｙ（ｎ）を決定するように構成され、
【数５】
2019507344.tif 000018
ここで、ｘおよびｙが各相間電流を決定するために使用される２つの相を示す、どれが最小値であるかを判断するステップと
を含む方法。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、地絡故障を判断するための装置に関する。発明はまた、配電ネットワークおよび関連する方法にも関する。
【背景技術】
【０００２】
本開示は地絡故障を判断するための装置および方法に関する。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】欧州特許第２７３３４９８号
【発明の概要】
【課題を解決するための手段】
【０００４】
発明の第１の態様によれば、多相配電ネットワーク上の地絡故障に関係する相を決定するための装置が提供され、上記装置は、スパイク故障を判断することに基づき、相間故障電流の複数のＲＭＳ値のうちのどれが最小値であるかを判断するように構成され、上記最小ＲＭＳ値が、上記地絡故障に関係する上記多相配電ネットワークの特定の相を表しており、
上記装置が、所定のサンプリングレートで求めた複数のサンプルセットおよび残留電流値、ｉＲ、に基づき上記スパイク故障の発生を判断するように構成され、上記サンプルセットが、特定のサンプリング時刻における、上記多相配電ネットワークの各相の相電流を含み、上記残留電流値、ｉＲ、が上記相電流の総和を含み、上記装置が次式のように上記残留電流のＲＭＳ値、ＩＲ（ｎ）、を決定するように構成され、
【０００５】
【数１】
2019507344.tif 000003
ここで、ｋが上記総和のインデックス値であり、ｎがスライディング窓の最新のサンプルを示し、Ｎがサンプルの数に関する上記スライディング窓サイズであり、上記スパイク故障の上記発生は、少なくともＩＲ（ｎ）が所定のスパイク時間しきい値（ＩＲｔｉｍｅ）よりも短い間所定の残留電流しきい値（ＩＲｓｅｔ）よりも大きいことに基づいて判断され、
上記装置は、各サンプルセットの上記サンプルのそれぞれの相電流の間の差違を含む相間電流ｉｘｙ（ｎ）および上記スライディング窓の始まりと終わりとにおける各相についての上記相間電流の間の差違を含む相間故障電流Δｉｘｙ（ｎ）に基づいて、上記相間故障電流Δｉｘｙ（ｎ）の複数のＲＭＳ値のうちのどれが最小値であるかを判断するように構成され、上記装置が次式のように上記相間故障電流の上記ＲＭＳ値ΔＩｘｙ（ｎ）を決定するように構成され、
【０００６】
【数２】
2019507344.tif 000004
ここで、ｘおよびｙが各相間電流を決定するために使用した２つの相を示す。
【０００７】
１つまたは複数の例では、上記装置が、相ａ、ｂおよびｃを含む３相配電ネットワークで使用するために構成され、最小値である上記相間故障電流の上記ＲＭＳ値は、上記最小値ΔＩｍｉｎが、
ｉ．  ΔＩｂｃから成る場合には、上記地絡故障が相ａに関係付けられ、
ｉｉ．  ΔＩｃａから成る場合には、上記地絡故障が相ｂに関係付けられ、
ｉｉｉ．  ΔＩａｂから成る場合には、上記地絡故障が相ｃに関係付けられる
ことで地絡故障が存在する上記多相配電ネットワークの特定の相を表している。
【０００８】
１つまたは複数の例では、上記装置は、上記スパイク故障の上記発生を判断する際に、上記残留電流の上記ＲＭＳ値の変化率、および上記変化率が変化率のしきい値よりも大きいか否かを決定するように構成され、
上記スパイク故障の上記発生は、少なくともＩＲ（ｎ）が上記所定のスパイク時間しきい値よりも短い間上記所定の残留電流しきい値よりも大きいことおよびＩＲ（ｎ）が上記所定の残留電流しきい値よりも大きい上記時間の少なくとも一部の間に上記変化率が上記変化率のしきい値よりも大きいことに基づいて判断された。
【０００９】
１つまたは複数の例では、上記装置は、上記故障が、
複数の相の各々の組合せについて上記相間故障電流の上記ＲＭＳ値ΔＩｘｙ（ｎ）および上記残留故障電流のＲＭＳ値ΔＩＲ（ｎ）に基づいて順故障方向のものであるかどうかを判断するように構成され、上記残留故障電流が上記スライディング窓の上記始まりにおいて上記サンプルセットのうちの１つから導き出された上記残留電流と上記スライディング窓の上記終わりにおいて上記サンプルセットのうちの１つから導き出されたものとの間の差違を含み、上記残留故障電流の上記ＲＭＳ値ΔＩＲ（ｎ）が、
【００１０】
【数３】
2019507344.tif 000005
から成り、
上記装置が、上記相間故障電流の上記ＲＭＳ値の最小値ΔＩｍｉｎおよび上記相間故障電流の上記ＲＭＳ値の最大値ΔＩｍａｘを特定するように構成され、上記装置は、
ΔＩｍｉｎが上記残留故障電流の上記ＲＭＳ値の第１の所定の割合未満であり、
ΔＩｍａｘが上記残留故障電流の上記ＲＭＳ値の第２の所定の割合よりも大きい
場合に順方向故障を判断するように構成される。
【００１１】
１つまたは複数の例では、上記第１の所定の割合が０．２以下または０．１５以下または０．１以下である。１つまたは複数の例では、上記第１の所定の割合が０．１と０．２との間または０．０５と０．２５との間である。１つまたは複数の例では、上記第２の所定の割合が０．８以上または０．８５以上または０．９以上である。１つまたは複数の例では、上記第２の所定の割合が０．８と０．９との間または０．７と０．９５との間である。
【００１２】
１つまたは複数の例では、上記装置が、地絡故障に関係する相を上記決定することに基づいて故障アラームを鳴らすように構成され、故障アラームを上記鳴らすことは、
上記スパイク故障の上記発生および上記スパイク故障中の上記地絡故障に関係する上記多相配電ネットワークの特定の相を上記特定すること、少なくとも所定の、アラーム、所定のアラーム期間内の回数
に依存する。
【００１３】
１つまたは複数の例では、上記故障アラームを上記鳴らすことは、上記スパイク故障の各発生に対して請求項３記載の故障方向を決定することにさらに依存する。
【００１４】
１つまたは複数の例では、上記サンプルセットの上記サンプリングレートが、２４００Ｈｚなどの２〜３ｋＨｚの範囲内である、
１つまたは複数の例では、上記スライディング窓サイズが、上記多相配電ネットワークの上記基本周波数の単一の周期に対応する。このように、上記サンプリングレートが２４００Ｈｚであり、上記基本周波数が５０Ｈｚである場合には、Ｎ＝２４００／５０＝４８サンプルである。したがって、上記スライディング窓は一番最近に取り込んだサンプルからＮ個の連続するサンプルを包含するはずである。
【００１５】
１つまたは複数の例では、上記所定の残留電流しきい値が、定格電流の１０％未満または５％未満である。上記定格電流は、上記配電ネットワークに接続された変流器の二次側からなどの上記多相配電ネットワークに印加される上記電流であり得る。
【００１６】
１つまたは複数の例では、上記所定の残留電流しきい値が、１００ｍＡ未満、または５０ｍＡ未満であり得る。
【００１７】
１つまたは複数の例では、上記所定のスパイク時間しきい値が、上記多相配電ネットワークの上記基本周波数の１サイクルに基づいている。このように、例えば、５０Ｈｚネットワークに関して、上記スパイク時間しきい値は、１／５０＊１０００＝２０ｍｓであり得る、または６０Ｈｚネットワークに関して、上記スパイク時間しきい値は、１／６０＊１０００＝１６．６７ｍｓであり得る。
【００１８】
発明の第２の態様によれば、上記第１の態様に記載の地絡故障に関係する相を決定するための装置を含む配電ネットワークが提供される。
【００１９】
発明の第３の態様によれば、多相配電ネットワーク上の地絡故障に関係する相を決定するための方法が提供され、上記方法が、スパイク故障を判断することに基づき、相間故障電流の複数のＲＭＳ値のうちのどれが最小値であるかを判断するステップを含み、上記最小ＲＭＳ値が、上記地絡故障に関係する上記多相配電ネットワークの特定の相を表し、
所定のサンプリングレートで求めた複数のサンプルセットおよび残留電流値、ｉＲ、に基づき上記スパイク故障の発生を判断するステップであって、上記サンプルセットが、特定のサンプリング時刻における、上記多相配電ネットワークの各相の相電流を含み、上記残留電流値、ｉＲ、が上記相電流の総和を含み、装置が次式のように上記残留電流のＲＭＳ値、ＩＲ（ｎ）、を決定するように構成され、
【００２０】
【数４】
2019507344.tif 000006
ここで、ｋが上記総和のインデックス値であり、ｎがスライディング窓の最新のサンプルを示し、Ｎがサンプルの数に関する上記スライディング窓サイズであり、上記スパイク故障の上記発生は、少なくともＩＲ（ｎ）が所定のスパイク時間しきい値（ＩＲｔｉｍｅ）よりも短い間所定の残留電流しきい値（ＩＲｓｅｔ）よりも大きいことに基づいて判断される、発生を判断するステップと、
各サンプルセットの上記サンプルのそれぞれの相電流の間の差違を含む相間電流ｉｘｙ（ｎ）および上記スライディング窓の始まりと終わりとにおける各相についての上記相間電流の間の差違を含む相間故障電流Δｉｘｙ（ｎ）に基づいて、相間故障電流Δｉｘｙ（ｎ）の複数のＲＭＳ値のうちのどれが最小値であるかを判断するステップであり、上記装置が次式のように上記相間故障電流の上記ＲＭＳ値ΔＩｘｙ（ｎ）を決定するように構成され、
【００２１】
【数５】
2019507344.tif 000007
ここで、ｘおよびｙが各相間電流を決定するために使用した２つの相を示す、どれが最小値であるかを判断するステップと
を含む。
【００２２】
第４の態様によれば、プロセッサに関係するメモリを有するプロセッサで実行されると、上記第３の態様の方法を実行するように構成されるコンピュータプログラムコードを備えるコンピュータプログラムまたはコンピュータプログラム製品が提供される。
【００２３】
ここで、単に例として、下記の図を参照して発明の実施形態の詳細な説明を続ける。
【図面の簡単な説明】
【００２４】
【図１】多相配電ネットワーク上の地絡故障に関係する相を決定するための装置の一例の図である。
【図２】図１の装置のスパイク故障素子の一例の図である。
【図３】図１の装置の故障方向素子の一例の図である。
【図４】順方向の地絡故障の相の期間に装置により使用されるパラメータの記録の一例の図である。
【図５】逆方向の地絡故障の相の期間に装置により使用されるパラメータの記録の一例の図である。
【図６】図１の装置のアラーム素子の一例の図である。
【図７】装置に接続された多相配電ネットワーク上の地絡故障に関係する相を決定するための装置を有する配電ネットワークの図である。
【図８】方法の一例を図示するフローチャートである。
【発明を実施するための形態】
【００２５】
複数の相のうちの１つまたは複数が、ネットワークを形成するケーブル配線の損傷のためなどで地面に意図しないで接続すると、地絡故障が配電ネットワークに生じることがある。配電ネットワーク内の地絡故障の検出は、特に、地絡故障が断続的であるときには困難である。配電ネットワークの電気ケーブルの地絡故障が初期段階であるときには、電気的絶縁が完全には破壊していないことがある。これが、特別に高い故障インピーダンスをともなう断続的な地絡故障という結果になることがある。この種類の故障を検出することは非常に難しく、故障方向を判断することが特に難しい。
【００２６】
発明の装置および方法は、上に述べた初期段階にある超高インピーダンス地絡故障の検出を提供できることが有利である。発明の装置は、係る故障を判断することが電流の測定だけに基づくということでやはり有利であり得る。それはそうとして、電圧測定は必要ではない（とはいえ、いくつかの実装形態では使用されることがある）。装置は、特殊なハードウェアを何も追加せずに、故障方向検出および／または断続的な高インピーダンス故障検出を提供することができ、これを、大きな故障電流をともなう完全な地絡故障の今後の発生を予測するために使用することができる。このように、装置は、測定した電流だけを使用して故障方向検出を提供できることが有利であり、そして下記に説明するように、断続的な超高インピーダンス地絡故障検出を提供することができる。
【００２７】
図１は、多相配電ネットワーク上の地絡故障に関係する相を決定するための装置１００を示している。特に、装置１００は、故障スパイクと呼ばれるネットワーク内の電流の流れの断続的および／または過渡的な電気的乱れの間の地絡故障に関係する相を決定するように構成されることがある。装置には、多相配電ネットワークの相（Ａ、Ｂ、Ｃ）の各々の電流の入力サンプル１０１ａ、１０１ｂ、１０１ｃを供給する。この例では、多相配電ネットワークは、３相ネットワークであり、したがって装置は、３相で動作するように構成される。他の例では、異なる数の相を使用することができる。いくつかの例では、装置１００は、相の各々の電流を測定して電流サンプルのセットを生成するように構成されることがある。しかしながら、この例では、異なる装置（図示せず）が所定のサンプリングレートで相の各々の電流をサンプリングし、各サンプリング時刻における各相のサンプリングした電流ｉＡ、ｉＢ、ｉＣを含んでいる電流サンプルのストリームを装置に提供する。この例では、サンプリングレートは２４００Ｈｚであるが、２と３ｋＨｚとの間などの他のサンプリング周波数を使用することができる。
【００２８】
装置１００は、下記により詳細に説明するように、複数のサンプルセットに基づいてスパイク故障の発生を判断するためのスパイク故障素子１０２を含む。スパイク故障フラッグ信号１０３（ＤＤＢ＿ＳＴと呼ぶ）を、スパイク故障の検出でスパイク故障素子１０２により生成する。装置は、相間故障電流の複数のＲＭＳ値のうちのどれが最小値であるかを判断するように構成された故障相検出素子１０４をさらに備え、前記最小ＲＭＳ値は、地絡故障に関係する多相配電ネットワークのうちの特定の相を表す（その動作を下記により詳細に論じよう）。故障相検出素子１０４は、どの相に地絡故障があるかを判断することに応じて、相Ａ故障信号１０５（ＤＤＢ＿ＦＴＡと呼ぶ）、相Ｂ故障信号１０６（ＤＤＢ＿ＦＴＢと呼ぶ）または相Ｃ故障信号１０７（ＤＤＢ＿ＦＴＣと呼ぶ）のうちの１つを含む相インジケータ信号の生成を行う。信号１０３および信号１０５〜１０７のうちの１つは、装置１００の最終的な出力を与えることができる。
【００２９】
装置１００は、この例では、装置の精度および／または機能を高めることができるいくつかの任意選択の構成要素を含む。特に、装置１００は、故障方向を判断するように構成された方向素子１０８を含む。故障方向信号１１０（ＤＤＲ＿ＤＩＲと呼ぶ）は、順故障方向を指示することができる故障の判断した方向を指示する方向素子１０８により与えられる。装置１００は、さらなる装置へ故障アラーム指示信号（包括的に１１２、１１３、１１４）を提供するためのアラーム素子１１１をやはり含むことができ、さらなる装置は、修理／検査のために配電ネットワークの悪影響を受けた部分を停止させることができる。アラーム素子１１１は、アラームが関係する相の指示を相故障アラーム指示信号として与えることができ、そして相Ａアラーム１１２（ＤＤＢ＿Ａｌａｒｍ＿Ａと呼ぶ）、相Ｂアラーム１１３（ＤＤＢ＿Ａｌａｒｍ＿Ｂと呼ぶ）および相Ｃアラーム１１４（ＤＤＣＢ＿Ａｌａｒｍ＿Ｃと呼ぶ）のうちのいずれか１つをこのように含むことができる。
【００３０】
装置１００は、１１５において複数の所定のしきい値または設定値を受信するようにさらに示されている。設定値は、所定の残留電流しきい値１１６、所定のアラーム時間Ｔｓｅｔ１１７、所定のアラーム回数Ｎｓｅｔ１１８、およびＣＴ比１１９を含む。ＣＴ比は、配電ネットワークに接続された変流器の二次巻線に対する一次巻線の巻き数の比を含むことができる。
【００３１】
一般に、スパイク故障素子１０２は、各相の電流サンプルに基づいて、断続的／過渡的な「スパイク」または恒久的な地絡故障の間を区別するように構成される。スパイク故障素子１０２は、最新の時間窓をカバーしている複数の残留電流値、各残留電流値が特定の時刻における個々の相電流の総和を含む、の二乗平均平方根（ＲＭＳ）値を使用し、そしてスパイク故障の発生を判断するために所定のしきい値とのＲＭＳ値の比較を使用する。いくつかの例では、しきい値に対するＲＭＳ値の微分値の比較を加えて使用する。
【００３２】
一般に、方向素子１０８は、相間電流差から導き出されるＲＭＳ値の大きさをしきい値に対して比較して、故障方向を判断するように構成される。
【００３３】
一般に、故障相検出素子１０４は、相間故障電流から導き出されたＲＭＳ値の大きさを比較して、故障のある相を表す最小値を判断するように構成される。
【００３４】
一般に、アラーム素子１１１は、故障スパイクの時間において特定される故障のある相に応じて、所定のスライディング時間窓にわたって判断した故障スパイクの数をカウントし、カウント値が所定のしきい値よりも大きい場合にはアラームを鳴らすように構成される。
【００３５】
図２は、スパイク故障素子１０２のより詳細な図を示している。スパイク故障素子は、初期段階の地絡故障を表す場合がある超高インピーダンス故障に応答するように構成される。スパイク故障素子１０２は、２００において、各相１０１ａ、１０１ｂ、１０１ｃの電流サンプルのサンプルセットを受信する。スパイク故障素子１０２はまた、２０１において、残留電流しきい値の値ＩＲｓｅｔを受信するように示されており、残留電流しきい値の使用を下記に説明しよう。しきい値の値は、素子１０２の（ハードウェア、ソフトウェアまたはファームウェア内の）論理を用いて積分されてもよいことまたはユーザ設定もしくはユーザ調節可能であってもよいことが認識されるだろう。
【００３６】
スパイク故障素子１０２は、下記のステップを実行するように構成される。スパイク故障素子は、ブロック２０２において、相電流サンプルの各セットについて残留電流値を決定するように構成される。残留電流値ｉＲを、特定のサンプリング時刻における個々の相電流ｉＡ、ｉＢ、ｉＣの総和により決定する（式１）；
ｉＲ＝ｉＡ＋ｉＢ＋ｉＣ                                          （１）
ここで、ｉＲが残留電流サンプルであり、ｉＡが相Ａの電流サンプルであり、ｉＢが相Ｂの電流サンプルであり、そしてｉＣが相Ｃの電流サンプルである。
【００３７】
他の例では、残留電流値の計算を、例えば、電流サンプルのサンプリングレートよりも小さいレートで、受信した相電流サンプルのサブセットについて実行する。
【００３８】
残留電流のＲＭＳ値を、次いで、ブロック２０３において、現在の時刻から広がる窓サイズＮにより規定される最新の期間にわたって決定される複数の残留電流値を使用して計算する。式２は、ＲＭＳ値の計算を示している；
【００３９】
【数６】
2019507344.tif 000008
ここで、ｎはごく最近のインデックスでありしたがって最新の時刻における関係する値を表し、例えば、ＩＲ（ｎ）はｎにおけるごく最近のサンプルセットから導き出した残留電流値を表すことができ；ＮはＲＭＳ値を形成するために使用するサンプルの数に関する、すなわち相電流のサンプルセットの数に対応するスライディング窓サイズであり、そしてｋはサム関数のインデックスを表している。
【００４０】
窓サイズを、この例では、多相配電ネットワークの基本周波数の１サイクルに設定する。例えば、相電流のサンプリングレートが２４００Ｈｚであり、ネットワークの基本周波数が５０Ｈｚである場合には、Ｎ＝２４００／５０＝４８サンプルである。
【００４１】
このように、式２の各ＲＭＳ値は、使用する残留電流値に関して、次のＲＭＳ値とＮ−１サンプルの重なりを持つことができる。他の例では、ＲＭＳ値を、異なる重なりでまたは重なりなしに計算することができる。
【００４２】
故障スパイクの発生の判断は、多数のパラメータに基づくことがある。いくつかの例では、ＲＭＳ値の大きさをしきい値と比較することができ、そしてＲＭＳ値がしきい値を超える時間の長さを、故障スパイクの発生を判断するために使用することができる。この例では、下記に説明するように、さらなる任意選択のパラメータ、すなわちＲＭＳ値の微分値を使用する。
【００４３】
各ＲＭＳ値ＩＲ（ｎ）を、ブロック２０４において所定の残留電流しきい値（ＩＲｓｅｔ）と比較する。ＩＲ（ｎ）の値がＩＲｓｅｔよりも大きい場合には、２０５においてフラッグＤＤＢ＿ＩＲを立てる。ＩＲｓｅｔよりも大きいＩＲ（ｎ）は、電気的乱れを表すことができ、フラッグを立てることは、このパラメータが満足であることの指示を与える。故障が存在するかどうかの最終的な判断には、多数のパラメータが同時に満足されるべきであることを必要とすることがある。
【００４４】
所定の残留電流しきい値ＩＲｓｅｔを、電力配電ネットワークに接続された変流器の二次巻線の定格電流に基づいて決定することができる。典型的には、変流器の定格二次電流は１Ａである。例えば、所定の残留電流しきい値ＩＲｓｅｔを、前記定格「二次」電流の１０％未満または５％未満または約０．０５とすることができる。例えば、前記定格電流が１Ａである場合には、しきい値ＩＲｓｅｔを５０ｍＡ（すなわち５％）とすることができる。
【００４５】
ブロック２０６は、ＩＲ（ｎ）値がＩＲｓｅｔよりも大きい時間の長さの決定を行う。ブロック２０６はこのように、過渡なスパイクと恒久的な地絡故障との間の区別を行う。この例では、ブロック２０６は、フラッグＤＤＢ＿ＩＲで動作するが、他の例では、ブロック２０３の出力、しきい値２０１およびタイマを使用できる。ブロック２０６はこのように、ブロック２０７においてフラッグ信号ＤＤＢ＿ＩＲの時間的な幅を測定し、ブロック２０８において所定のスパイク時間しきい値（ＩＲｔｉｍｅ）２０９と比較する。所定のスパイク時間しきい値を、この例では、基本周波数の１サイクルに設定する。このように、ＤＤＢ＿ＩＲフラッグの時間的な幅がＩＲｔｉｍｅ＝２０ｍｓ（５０Ｈｚ方式に関して）よりも長い、または６０Ｈｚ方式に関してＩＲｔｉｍｅ＝１／６０＊１０００＝１６．６７ｍｓよりも長い場合には、故障がスパイク故障でないことを上記は表している。
【００４６】
ブロック２０８は、時間的な幅がＩＲｔｉｍｅしきい値を超えるときには「真」条件を生成する、このようにブロック２０６は、論理ＮＯＴブロック２１０を含み、その結果、フラッグ信号２１１はＩＲｔｉｍｅしきい値よりも短い間続くＤＤＢ＿ＩＲフラッグを表し、したがって「スパイク」を含む。
【００４７】
上に述べたように、フラッグ２０５およびフラッグ２１１に基づく論理ＡＮＤ関数は、スパイク故障フラッグ信号１０３の生成を行うことができる。しかしながら、この例は、ブロック２０３から微分ブロック２１２へ与えられるＲＭＳ値を示し、微分ブロック２１２は２つ以上のＲＭＳ値の時間に関する微分値、すなわち、残留電流のＲＭＳ値の変化率を決定できる。
【００４８】
ブロック２１３は、ＲＭＳ値の微分値を変化率のしきい値２１４と比較する。ＲＭＳ値の微分値が変化率のしきい値よりも大きい場合には、フラッグＤＤＢ＿ｄｅｌｔａＲ２１５を立てる。フラッグ２１５は、ブロック２１３の条件が満足されると、ブロック２１６により、ドウェル時間の間滞在するように与えられる。このように、フラッグ２１５は、ブロック２１３における条件がドウェル時間よりも短い間しか満足しないにもかかわらずドウェル時間の間与えられるであろう。この例では、ドウェル時間は１０ｍｓである。ドウェル時間は、基本周波数未満、または約０．５＊基本周波数以下などの基本周波数に関係することがある。
【００４９】
変化率のしきい値は、この例では、所定の残留電流しきい値ＩＲｓｅｔに基づくが、他の実施形態では、所定の残留電流しきい値とは無関係であってもよい。特に、変化率のしきい値は、この例では、０．５＊ＩＲｓｅｔを含む。ブロック２２０は、ＩＲｓｅｔ値のスケーリング調節を行って、ブロック２１３により使用される変化率のしきい値２１４を形成する。
【００５０】
ＤＤＢ＿ＩＲフラッグ２０５を、論理ＡＮＤブロック２１７によりＤＤＢ＿ｄｅｌｔａＲフラッグ２１５と組み合わせて、ＤＤＢ＿ＦＴフラッグ信号２１８を生成する。
【００５１】
さらなる論理ＡＮＤブロック２２１を、ＤＤＢ＿ＦＴフラッグ信号２１８と過渡的／恒久的地絡故障区別ブロック２０６からのフラッグ信号２１１とを組み合わせるために設けて、素子１０２の最終的な出力、スパイク故障フラッグ信号ＤＤＢ＿ＳＴ１０３を与える。遅延時間だけブロック２２１へのＤＤＢ＿ＦＴフラッグ信号２１８の提示を遅らせるために、遅延ブロック２２２を使用することができる。遅延時間は、この例では、２０ｍｓ、または基本周波数に関係するある他の値もしくは基本周期未満を含む。
【００５２】
スパイク故障の存在の推定を素子１０２により行うことができるそれぞれのフラッグを生成できる他のパラメータを使用することがでることが認識されるだろう。
【００５３】
図３は、方向素子１０８のより詳細な図を示している。方向素子１０８は任意選択の構成部品であり、故障の方向、および特に、故障が順方向のものであるかどうかを判断するように構成される。順方向故障は、保護デバイスが置かれている装置により保護された回線内の故障であることが認識されるだろう。
【００５４】
故障方向素子１０８は、３００において、各相１０１ａ、１０１ｂ、１０１ｃの電流サンプルのサンプルセットを受信する。ブロック３０１は、各サンプルセットのサンプルのそれぞれの相電流の間の差違を含んでいる相間電流ｉｘｙ（ｎ）、そして引き続いて、スライディング窓の始まりと終わりとにおける各相についての相間電流の間の差違を含んでいる相間故障電流Δｉｘｙ（ｎ）を決定するように構成される。
【００５５】
このように、ブロック３０１は、式３を使用して相間電流ｉｘｙ（ｎ）を決定するように構成され、下記の通りである；
ｉａｂ（ｎ）＝ｉａ（ｎ）−ｉｂ（ｎ）
ｉｂｃ（ｎ）＝ｉｂ（ｎ）−ｉｃ（ｎ）
ｉｃａ（ｎ）＝ｉｃ（ｎ）−ｉａ（ｎ）                              （３）
このように、各サンプルセット内のサンプリングした相電流の間の差違を式３により決定し、上記の差違を次いで、相間故障電流Δｉｘｙ（ｎ）を導き出すために使用する。相間故障電流は、第１の時刻における相間電流のうちの１つと第２の時刻における同じ相間電流との間の差違を含む。第１の時刻および第２の時刻は、この例では、上に説明したスライディング窓の始まりと終わりとを含み、ここではスライディング窓が、（この例では）定格周波数の１サイクルをカバーする。このように、相間故障電流は式４により決定され；
【００５６】
【数７】
2019507344.tif 000009
ここで、ｎが現在の時刻において関係付けられた期間を表し、Ｎがスライディング窓サイズである。スライディング窓サイズは、上に述べたように、基本周波数の１周期を含むことができ、そしてしたがってＮが基本周波数のサイクル当たりのサンプルの数である。相電流のサンプリングレートが２４００Ｈｚであり、基本周波数が５０Ｈｚである場合には、Ｎ＝２４００／５０＝４８である。
【００５７】
式４は、時刻ｎにおける残留電流と時刻ｎ−Ｎにおける残留電流との間の差違から残留故障電流ΔｉＲ（ｎ）を決定することを加えて示している。
【００５８】
装置は、ブロック３０２において、一般式５により下記のように、相間故障電流のＲＭＳ値ΔＩｘｙ（ｎ）を決定するように構成される；
【００５９】
【数８】
2019507344.tif 000010
ここで、ｘおよびｙは、各相間電流を決定するために使用される２つの相を示す。
【００６０】
このように、この３相の例では、残留故障電流のＲＭＳ値とともに３相間故障電流のＲＭＳ値を、より具体的に式６〜９により計算する：
【００６１】
【数９】
2019507344.tif 000011
ブロック３０３は、式６〜式９により決定した３相間故障電流のＲＭＳ値のうちのどれが最小値であり、どれが最大値であるかを判断するように構成される。したがって、ブロック３０３は、下記の式１０および式１１を使用する；
【００６２】
【数１０】
2019507344.tif 000012
ブロック３０３は、しきい値比較のための最小値および最大値を与える。最小ＲＭＳ値ΔＩｍｉｎを比較ブロック３０４に与える。最大ＲＭＳ値ΔＩｍａｘを比較ブロック３０５に与える。比較ブロック３０４、３０５は、残留故障電流のＲＭＳ値に基づいてしきい値を使用する。ブロック３０６は、残留故障電流のＲＭＳ値の第１の所定の割合を決定し、ブロック３０７は、残留故障電流のＲＭＳ値の第２の所定の割合を決定する。この例では、第１の所定の割合は、０．１から成り、第２の所定の割合は、０．９から成る。しかしながら、他のしきい値をネットワークの特性または決定した方向の所望の信頼度に合うように使用することができることが認識されるだろう。
【００６３】
論理ＡＮＤブロック３０８は、比較ブロック３０４および３０５の出力を統合する。このように、ブロック３０４〜３０７の組合せは、３相間故障成分電流のＲＭＳ値のうちの最小ＲＭＳ値を残留故障電流のＲＭＳ値の０．１倍と比較することにより、および最大値を残留故障電流のＲＭＳ値の０．９倍と比較することにより順方向故障があるか否かの判断を行う。すなわち、
【００６４】
【数１１】
2019507344.tif 000013
である場合には、順方向故障が想定され、そして故障方向フラッグ信号１１０を立てて、順方向故障を特定したことを示す。
【００６５】
このように、上に説明したように、相の各々の組合せについて相間故障電流Δｉｘｙ（ｎ）（ブロック３０１において計算される）のＲＭＳ値（ブロック３０２において計算される）および残留故障電流のＲＭＳ値Δｉｘｙ（ｎ）（ブロック３０２において計算される）に基づいて、故障が順故障方向のものであるかどうかを判断するように、ブロック１０８は構成され、残留故障電流は、スライディング窓の始まりにおけるサンプルセットのうちの１つの残留電流とスライディング窓の終わりにおけるサンプルセットのうちの１つとの間の差違（ブロック３０１において計算される）を含む。さらに、ブロック３０３において相間故障電流のＲＭＳ値のうちの最小値ΔＩｍｉｎおよび相間故障電流のＲＭＳ値のうちの最大値ΔＩｍａｘを特定し、そして
ΔＩｍｉｎが残留故障電流のＲＭＳ値の第１の所定の割合（ブロック３０４において計算される）よりも小さく、かつ
ΔＩｍａｘが残留故障電流のＲＭＳ値の第２の所定の割合（ブロック３０５において計算される）よりも大きい
場合には、ブロック１０８は、ブロック３０４、３０５において順方向故障を判断するように構成される。
【００６６】
図４は、ある相が順方向の地絡故障に至る間に式６〜式９により計算される３相間故障電流のＲＭＳ値の記録例４００を残留故障電流のＲＭＳ値とともに示している。特に、ΔＩａｂ（ｎ）の記録を４０１に示し、ΔＩｂｃ（ｎ）を４０２に示し、ΔＩｃａ（ｎ）を４０３に示し、ＩＲを４０４に示している。
【００６７】
図５は、ある相が逆方向の地絡故障に至る間に式６〜式９により計算される３相間故障電流のＲＭＳ値の記録例５００を残留故障電流のＲＭＳ値とともに示している。特に、ΔＩａｂ（ｎ）の記録を５０１に示し、ΔＩｂｃ（ｎ）を５０２に示し、ΔＩｃａ（ｎ）を５０３に示し、ＩＲを５０４に示している。
【００６８】
図１に戻って、故障相検出素子１０４は、１２０において、各相１０１ａ、１０１ｂ、１０１ｃの電流サンプルのサンプルセットを受信するように構成される。式３、式４および式６から式１０により決定したように、故障相検出素子１０４は、相間故障電流のＲＭＳ値のうちの最小値を決定するように構成される。したがって、故障相検出素子１０４は、最小ＲＭＳ値を受信するためにブロック３０３と通信することができる。あるいは、素子１０４は、ブロック３０１、３０２および３０３と同様に論理回路（ソフトウェア、ハードウェアまたはファームウェア）を組み込むことができる。式１０によって、最小値だけを決定するブロック３０３の修正版を与えることができることが認識されるだろう。
【００６９】
故障相検出素子１０４は、下記のように最小ＲＭＳ値から地絡故障に関係する相を決定する；
ｉ．  ΔＩｍｉｎ＝ΔＩｂｃである場合には、素子１０４は、故障相が相Ａであることを判断し、相Ａ故障信号１０５を与える。
【００７０】
ｉｉ．  ΔＩｍｉｎ＝ΔＩｃａである場合には、素子１０４は、故障相が相Ｂであることを判断し、相Ｂ故障信号１０６を与える。
【００７１】
ｉｉｉ．  ΔＩｍｉｎ＝ΔＩａｂである場合には、素子１０４は、故障相が相Ｃであることを判断し、相Ｃ故障信号１０７を与える。
【００７２】
いくつかの例では、装置１００の出力を、スパイク故障判断素子１０２の出力１０３および故障相検出素子１０４の出力１０５、１０６、１０７に基づいて地絡故障を有する相を表す信号とすることができる。このように、素子１０２および１０４からの出力は、同時に装置１００の出力を与えることができる。
【００７３】
この例では、装置は、故障相検出素子１０４、スパイク故障素子１０２および方向素子１０８から出力を受信するアラーム素子１１１を含み、アラームを鳴らす前に出力のさらなる処理を行う。アラームは、装置１００の出力を含むことができる。
【００７４】
図６は、アラーム素子１１１のより詳細な図を示している。素子１１１は、スパイク故障フラッグ信号１０３（ＤＤＢ＿ＳＴと呼ぶ）、故障方向信号１１０（ＤＤＢ＿ＤＩＲと呼ぶ）、およびどの相に地絡故障があるかの決定に応じて、相Ａ故障信号１０５（ＤＤＢ＿ＦＴＡと呼ぶ）、相Ｂ故障信号１０６（ＤＤＢ＿ＦＴＢと呼ぶ）または相Ｃ故障信号１０７（ＤＤＢ＿ＦＴＣと呼ぶ）のうちの１つを受信するように構成される。したがって、図６に示したように３つの配置がある場合があり、各相について１つそれゆえ相故障信号１０５、１０６、１０７のうちの１つを各々が受信する。素子１１１は、上に説明したフラッグ信号を受信するための論理ＡＮＤブロック６００を含む。ＤＤＢ＿ＳＴと、ＤＤＢ＿ＤＩＲと、ＤＤＢ＿ＦＴｘフラッグ信号のうちの１つとを同時に受信することが生じると、ブロック６００は故障フラッグ信号６０１を発生する。
【００７５】
特に、装置は、フラッグ信号ＤＤＢ＿ＳＴ、ＤＤＢ＿ＤＩＲおよびＤＤＢ＿ＦＴＡをブロック６００により受信すると、相Ａの順方向スパイク故障を判断するように構成される。装置は、フラッグ信号ＤＤＢ＿ＳＴ、ＤＤＢ＿ＤＩＲおよびＤＤＢ＿ＦＴＢをブロック６００により受信すると、相Ｂの順方向スパイク故障と判断するように構成される。最後に、装置は、フラッグ信号ＤＤＢ＿ＳＴ、ＤＤＢ＿ＤＩＲおよびＤＤＢ＿ＦＴＣをブロック６００により受信すると、相Ｃの順方向スパイク故障と判断するように構成される。各相についてのこれらの順方向スパイク故障信号は、装置１００の出力を形成することができる。
【００７６】
しかしながら、本例では、所定の最新の時間の長さ内に所定の数の故障フラッグ信号を受信するまで、アラーム信号は発生されない。カウンタブロック６０２は、６０３で示した所定のアラーム期間、Ｔｓｅｔにより規定される最新の期間内に受信される故障フラッグ信号６０１の数をカウントするように構成される。カウンタの出力６０４は、このように故障フラッグ信号がスライディング時間窓内に発生される回数の整数を含む。比較ブロック６０５は、出力６０４を６０６で示した所定の「アラーム」回数、Ｎｓｅｔと比較する。出力６０４がＮｓｅｔしきい値を超える場合には、アラーム信号６０７を発生する。図６の配置を相故障信号の各々に対して設けることが与えられると、アラーム信号６０７は、アラームを発生する相をやはり表している。
【００７７】
１１７および１１８において装置１００により受信される所定の値ＮｓｅｔおよびＴｓｅｔを、ユーザが規定したものとすることができる。
【００７８】
図１〜図３および図６に関連して上に論じた構成部品の機能を、ハードウェア、ソフトウェアまたはファームウェアに与えることができる。
【００７９】
図７は、相電流サンプリングデバイス（図示せず）を介してなどで多相配電ネットワークに接続された装置１００を有する多相配電ネットワーク７００を示している。装置１００は、ネットワーク７００の修理／検査を行うための故障モニタリング装置および／または故障アクション応答装置にさらに接続されても一部を形成してもよい。多相配電ネットワーク７００は、１つだけの電源を含むことができる。
【００８０】
図８は、スパイク故障の決定８０１に基づいて、相間故障電流の複数のＲＭＳ値のうちのどれが最小値であるかを判断するステップ８０２を示すフローチャートを示しており、前記最小ＲＭＳ値は、地絡故障に関係する多相配電ネットワークの特定の相を表している。
【００８１】
装置の利点は、
装置が配電ネットワークのケーブルを継続的にモニタすることを行うことができ、
装置が接地モードのすべてに利用可能であり、
（装置１００相により使用される配電ネットワークの唯一の測定値が電流サンプルであるので）電圧入力なしにフィーダリレーへと装置を容易に埋め込むことができ、
装置は単純であり、計算経費が小さく、
より優れた信頼性を提供でき、
インテリジェントなアルゴリズムを利用する方法よりもさらに実際的であり得る。
【符号の説明】
【００８２】
１００  装置
１０１ａ、１０１ｂ、１０１ｃ  入力サンプル
１０２  スパイク故障素子
１０３  スパイク故障フラッグ信号
１０４  故障相検出素子
１０５  相Ａ故障信号
１０６  相Ｂ故障信号
１０７  相Ｃ故障信号
１０８  方向素子
１１０  故障方向信号
１１１  アラーム素子
１１２  相Ａアラーム
１１３  相Ｂアラーム
１１４  相Ｃアラーム
１１５  所定のしきい値
１１６  所定の残留電流しきい値
１１７  所定のアラーム時間
１１８  所定のアラーム回数
１１９  ＣＴ比
２０２  ブロック
２０３  ブロック
２０４  ブロック
２０５  フラッグ
２０６  ブロック
２０７  ブロック
２０８  ブロック
２０９  所定のスパイク時間しきい値
２１０  論理ＮＯＴブロック
２１１  フラッグ信号
２１２  微分ブロック
２１３  ブロック
２１４  変化率のしきい値
２１５  ＤＤＢ＿ｄｅｌｔａＲフラッグ
２１６  ブロック
２１７  論理ＡＮＤブロック
２１８  ＤＤＢ＿ＦＴフラッグ信号２１８
２２０  ブロック
２２１  ブロック
２２２  遅延ブロック
３０１  ブロック
３０２  ブロック
３０３  ブロック
３０４  比較ブロック
３０５  比較ブロック
３０６  ブロック
３０７  ブロック
３０８  論理ＡＮＤブロック
４００  記録例
５００  記録例
６００  論理ＡＮＤブロック
６０１  故障フラッグ信号
６０２  カウンタブロック
６０４  出力
６０５  比較ブロック
６０７  アラーム信号
７００  多相配電ネットワーク
８０１  スパイク故障の決定
８０２  最小値の判断
【図１】
2019507344.tif 000019
【図２】
2019507344.tif 000020
【図３】
2019507344.tif 000021
【図４】
2019507344.tif 000022
【図５】
2019507344.tif 000023
【図６】
2019507344.tif 000024
【図７】
2019507344.tif 000025
【図８】
2019507344.tif 000026
【国際調査報告】
2019507344.tif 000027
2019507344.tif 000028