最大エントロピースペクトル法に基づく電力分配ケーブル欠陥局在技術

サマリー

このプロトコルは、周波数領域反射率測定を用いて長距離配電ケーブルにおける正確な欠陥局在化を最大エントロピースペクトル推定に基づく手法を提供します。高周波ケーブルモデリング、反射係数スペクトル解析、スペクトル外挿、シミュレーションや実験による検証を詳細に扱い、従来のFFT手法よりも優れた精度を実現しています。

要約

現代の都市配電網におけるケーブル敷設距離が拡大し続ける中、周波数領域反射率(FDR)による故障位置の精度は、長さが増すにつれて著しく低下しています。長ケーブルの精度向上のために、本研究は最大エントロピースペクトル推定に基づく手法を提案します。配電ケーブル用の高周波分布パラメータモデルを開発し、絶縁の劣化や機械的損傷が単位長の静電容量およびインダクタンスに与える影響を定量的に評価します。伝送線路理論を用いて、反射係数スペクトルと欠陥位置の関係を導き出し、スペクトル分解能を高めるためのステップ周波数最適化基準も導出します。高速フーリエ変換(FFT)ベースの周波数領域解析におけるスペクトル漏れや限られた分解能を軽減するため、高忠実度スペクトル推定には最大エントロピースペクトル法が用いられ、故障距離の位置特定精度を向上させています。シミュレーションの結果、従来のFFTと比較して、提案された手法は典型的な長ケーブルケースで位置誤差を2〜4.5倍削減できることが明らかになりました。実験結果により、最大エントロピーアプローチでは相対誤差が0.25%未満であることが確認され、従来の手法の0.55%未満の誤差を上回り、その有効性と工学的優位性が検証されました。

概要

現代の電力システムにおける電力供給の信頼性と安全性への需要が高まる中、ケーブルは都市配電網で広く採用されています。その中でも、著者の^{卓越した利点}により架橋ポリエチレン(XLPE)ケーブルが広く利用^{されています。}配電ケーブルは通常、過酷な運用環境の地下に設置されます。サービス時間が長くなり、負荷の変動や過負荷が重なると、ケーブルの劣化がより顕著になります。その結果、水木、電気木、部分的な放電、断熱湿気などの局所的な欠陥が徐々に発生することがあります。これらの欠陥が遅れて対処されなければ、深刻な故障に発展し、最終的にはグリッド^3,4の停電につながる可能性があります。したがって、配電ケーブルの欠陥の特定と故障位置特定は、電力網の運用と保守において不可欠な要素となっています。

配電ケーブルの故障位置特定方法は大きく分けてオンライン方式とオフライン方式の2種類に分類できます。オンライン方式は主に移動波位置特定法、部分放電(PD)監視法、共通モード漏れ電流法に代表されます。移動波位置測定は、断層で発生する過渡的な移動波を解析することで断層の位置を特定します。しかし、波面の正確な抽出の困難さ、非線形分布の断層特徴の重ね合い、設置の不便さなど、いくつかの課題に直面しています。さらに重要なのは、この方法は加齢に伴う欠陥の特定には適用できない^{ことです。}PDベースのオンライン監視はケーブル内で発生する欠陥信号を検出し、動作条件のリアルタイム評価を可能にします。それでもなお、老化の欠陥の認識や長期的なモニタリングにおいて技術的なボトルネックに直面しています。信号干渉、ノイズ抑制、複雑なデータ解析などの問題は、実^{用的な応用}においてしばしば誤判を招く^6,7,8。共通モード漏れ電流法は、ケーブル端子における共通モード漏れ電流の変動を測定することで、劣化欠陥の位置と重症度を決定します。この技術は高感度で非侵襲的な検出を提供します。しかし、その干渉防止能力は比較的弱く、可変周波数機器の動作時に発生する高周波差動モード電流は測定精度を大幅に低下させる可能性があります。さらに、通電ケーブルのオンライン監視は高電圧環境での信号絶縁にも対応しなければならず、これにより機器コストや保守の^複雑さが増加^します。

オフライン手法としては、時間領域反射測定(TDR)¹²、周波数領域反射測定(FDR)¹³、時間-周波数領域反射測定(TFDR)¹⁴など、いくつかのアプローチが提案されています。これらの中でも、TDRは注入されたパルス信号の高周波成分にエネルギーが低く、ケーブルの減衰やノイズ干渉に非常に敏感で、位置精度が限られ、初期の欠陥検出能力も低いですケーブル故障の初期段階では¹⁵。TFDRはガウスチープパルスを注入することで、TDR¹⁶における高周波成分の不足を緩和します。しかし、TFDRは複雑な信号変調回路と抽出アルゴリズムを必要とし、高周波変調信号は急速に減衰するため、長距離ケーブルでの欠陥^{検出が困難}です。一方、FDRはスイープ周波数信号に等しい出力分布で注入します。最も豊富な高周波成分を含み、より高い位置精度を実現し、TDRと比較して感度も高いため、急速な開発が進められています。取得信号の種類に応じて、FDRはさらに広帯域インピーダンススペクトラム(BIS)法^18,19,20と反射係数(RCS)法21,22,23,24に分類されます。

2003年、あるグループがケーブル²⁵の入力反射係数スペクトルに基づくケーブル局所欠陥局在法を初めて提案しました。しかし、このアプローチは分散パラメータの頻度依存効果を考慮しておらず、誤判断につながる可能性があります。後に逆高速フーリエ変換(IFFT)を用いてケーブル入力の反射係数スペクトルから情報を抽出し、熱老化、放射老化、押出変形などの局所欠陥の効果的な局在化が可能^{になりました。}また、スイープ信号の高周波数が高いほど、周波数領域反射率(FDR)の空間分解能も高くなることが観察されました。しかし、長いケーブルに適用すると、スイープ信号の上位周波数が大幅に低下し、定位精度が低下します。

FDR法の核心は、獲得した反射信号に対してFFTを実施することにあります。しかし、FFTに固有のスペクトル漏れ問題があるため、研究者たちはいくつかの改良案を提案しています。ある研究では、離散フーリエ変換(DFT)における距離ウィンドウとカイザーウィンドウ関数を組み合わせて、反射係数スペクトルの実部から欠陥成分を抽出し、ケーブル²⁷の緩んだ銅シールド欠陥の位置特定を可能にしました。カイザーウィンドウベースのDFT²⁸ は入力反射係数スペクトルの実数部分に適用され、オープンサーキット、ショートサーキット、高抵抗、低抵抗条件下の故障タイプ識別を実現しました。ある研究では、中間接合部の水分検出と局在化のために、ケーブル入力インピーダンスの虚部にブラックマンウィンドウベースのDFTを適用することが^{提案されました。}DFTのハニング自己畳み込みウィンドウを用いた補間アルゴリズムが導入され、反射係数スペクトル³⁰の実部における周期成分の周波数と位相を推定し、ケーブルインピーダンス異常の極性決定が可能となりました。その後、ナタルウィンドウ、4次3項ナタルウィンドウ、2次ナタル自己畳み込みウィンドウ、3点補間DFTアルゴリズム³¹の導入により改良が達成されました。全体として、既存の研究は主にウィンドウ関数設計や補間アルゴリズムの精緻化に焦点を当てており、これらは一定程度で欠陥局在精度の向上に寄与しています。

FFTの原理の観点から、周波数領域反射率(FDR)では、注入された信号の周波数が高いほど測定精度が向上します。理論的には、注入周波数が無限大に近づくと、ケーブルの欠陥を完璧な精度で特定できます。しかし、ケーブルを通じた超高周波信号の伝送は深刻な減衰を経験し、中長・長距離ケーブルには適していません^32,33。したがって、FDRにおける最大注入周波数と達成可能な欠陥位置精度との間には本質的なトレードオフが存在します。この問題に対処するため、本論文では最大エントロピースペクトル法に基づくケーブル欠陥局在技術を提案します。この手法は、測定された低周波データから最大エントロピー基準に従って反射係数スペクトルの高周波部分を予測し、その後反射係数の最大エントロピースペクトル(MES)を計算して欠陥位置を特定します。特に長距離ケーブル用途でスペクトルリーケージや分解能の制約に制約される従来のFFTベースのFDR法と比較して、提案された最大エントロピースペクトル法は、限られた測定データから高周波情報を効果的に外挿することで、より高いスペクトル分解能を提供します。ウィンドウ関数ベースの改良が主にサイドローブを抑制し、解像度を根本的に向上させるのとは異なり、MESアプローチはピークシャープネスと位置特定精度の両方を向上させます。さらに、補間法やパラメトリックスペクトル推定法と比較して、MES法は解像度と堅牢性のバランスが優れており、長い配電ケーブルの初期欠陥検出に特に適しています。このアプローチにより、局在精度が向上し、従来の手法では相対誤差が0.55%から0.25%に低減され、欠陥ピークの特定を容易にする滑らかな局在曲線が得られます。

プロトコル

注意:テストケーブルのすべての測定はオフライン状態(電源切れ)で行ってください。スイープ周波数信号注入および反射係数取得が可能なベクトルネットワークアナライザ(VNA)または同等の周波数領域反射計測機器を使用します。特に指定がない限り、ケーブルの端は必ずオープンサーキットにしてください。短く直接的な同軸接続と導電テープを用いて確実な接触をすることで、接続反射を最小限に抑えます。推奨機器やソフトウェアについては 材料表 を参照してください。

1. ケーブル試料と測定セットアップの準備

1. テスト用の電力分配ケーブル(例:YJLV22-8.7/10 kV-3×70 XLPEケーブルまたは同等品)を選択してください。
2. ケーブルの全長と既知の欠陥位置(人工的に導入された場合)を測定・記録してください。
3. ケーブル内の電磁波の伝播速度を測定してください。
   1. ケーブルを楽器の片端で接続し、反対側は開けたままにしてください。
   2. 適切な周波数帯域で反射係数スペクトルを取得します。
   3. ケーブル全長に対応する振動周期(往復)を特定します。
   4. 速度を (1)のように計算します。
4. ケーブルヘッドの接続を準備してください。
   1. ケーブルの端を剥がしてコアとシールドを露出させます。
   2. 両面導電性銅テープを使って、機器のリード線をコアとシールドにしっかりと固定します。
   3. 枝のリードの長さを最小限にしてヘッドエンドの反射を減らしましょう。

2. 周波数スイープパラメータの設定

1. 選択基準を用いて周波数ステップ(Δf)を決定します。
   1. ケーブル端の振動周期(2)を用います。
   2. サイクルあたり少なくとも200ポイントのサンプリングポイントを確保するように設定してください。
   3. 周波数ステップΔfは最大観測可能な距離範囲を決定する: (3)
   4. 全ケーブル長Lを曖昧さなく解決できるようにする:
      (4)
   5. 本研究では、解像度向上のためにより厳格な条件が採用されています:(5)
2. 周波数帯を設定します。
   1. 下限を100kHz(または楽器の最低限)に設定してください。
   2. 上限を過度な減衰なしに達成可能な最大値(ケーブルの場合は通常30〜50 MHz<200 m、長距離ケーブルの場合は減少)に設定してください。
   3. 上限周波数上限fmaxは空間分解能に直接影響します:(6)
   4. fmaxが高いほど解像度は細かくなりますが、ケーブルの減衰によって制約されます。
3. 必要なサンプリングポイントの数を計算します(7)
4. 計算された周波数帯域、ステップ、ポイント(例:100 kHz–50 MHz、高解像度の場合は~5,000ポイント)で機器を構成します。

3. 反射係数スペクトルを取得する

1. 計測器の校正を行います(VNAにはオープンショーロード校正が推奨されます)。
2. 準備済みのケーブルヘッドを計器ポートに接続します。
3. 掃走周波数の正弦波信号を注入し、各離散周波数点で入力反射係数Γ(ω)を測定します。
   1. 測定された反射係数は次のように表せます。(8)
   2. ここで: Γ_D: 欠陥での反射係数、 : 位相定数、 D: 欠陥距離。
   3. この関係から、反射係数の位相が周波数とともに線形に変化し、距離推定の基礎となります。
4. 複素反射係数スペクトル(大きさと位相、または実部と虚部)を記録します。
5. 周波数と実数・虚数の反射係数としてデータを保存します。

4. 従来のFFT法を用いたデータ処理(比較用)

1. 反射係数スペクトルの虚部を抽出します。
2. サイドローブ漏れを減らすためにウィンドウ関数(例:ブラックマンウィンドウ)を適用します。
3. ウィンドウ化された虚部に対して逆高速フーリエ変換(IFFT)を行います。
   1. 逆フーリエ変換は次のように導入されます: (9)
4. 欠陥距離は次のように計算されます:
   (10)
   1. ここで t は逆フーリエ変換後のピーク位置に対応する時間です。
5. 反射係数スペクトル内の欠陥ピーク(ヘッドエンドおよび終端反射を除く)を特定します。

5. 欠陥局在化に最大エントロピースペクトル法を適用する

1. 取得した反射係数スペクトルデータを読み込みます。
2. デトレンドと想像上の部分の配列を中心に据えてください。
3. 自己回帰モデルの順序を選択します
   (11)
   1. ここで N はデータポイントの数です。
4. Burg再帰アルゴリズムを実装します。
   1. 前方および後方予測誤差を初期化します。
   2. 再帰的に反射係数と自己回帰パラメータを計算しつつ、順方向および後方誤差のパワーを最小化します。
   3. 最大エントロピー基準に基づいて自己相関列を外挿します。
5. 外挿された配列からパワースペクトルを計算します。
6. 最大エントロピーパワースペクトルのピークを特定します。
7. ピーク位置を欠陥距離に変換するには(12)が使われます。
8. ピークの鮮明さと局在誤差を従来のFFT結果と比較してください。

6. 結果の解釈

1. 正のピーク(欠陥検出に成功)と、偽のピークや最適でないトレースを区別します。
2. 相対局在誤差を(13)のように計算します。
3. 相対誤差<0.25%と、明確で滑らかなピークを大きなサイドローブなしに確認してください。

結果

反射係数スペクトルおよび欠陥局在の結果は、最大エントロピースペクトル法を用いたシミュレーションおよび実験の両方から得られ、従来のFFTベースの手法と比較されました。成功した欠陥局在は、正確な欠陥位置に鋭く目立つピークを持つ滑らかなパワースペクトル、サイドローブや漏れによる擬ピーク(偽ピーク)が最小限、そして低い相対誤差(通常<0.3%)を特徴とします。これに対し、従来のFFT手法による最適でない結果は、多数の偽ピーク、解像度の制限やスペクトル漏れによるより広い欠陥ピーク、そして最大0.55%に達する高い相対誤差を示し、ピーク識別を複雑にし精度を低下させることがあります。

シミュレーションでは、PSCADで構造パラメータを表 1に示した配線ケーブルモデルが確立されました。長さと欠陥位置の異なる3つのケーブルサンプルが構成されました(表2)。より短いケーブル(サンプル1)では、従来のウィンドウなしFFTは複雑なスペクトルを生成し、リークが大きく、相対誤...

ディスカッション

本研究では、最大エントロピースペクトル法(MES)を用いて周波数領域反射率測定(FDR)を用いた電力配電ケーブルの欠陥を位置特定するための詳細なプロトコルを提示します。これまでの研究は、スペクトル漏れ、長いケーブルでの限界、欠陥同定を複雑にする偽ピークなど、従来の高速フーリエ変換(FFT)に基づく解析に頼ることが多く、問題が解決しにくい。さらに、高度な信号処理の専門知識を持たない研究者は、周波数スイープパラメータの最適化、高周波データの外挿、複雑なスペクトル^35、36、37の解釈に課題を抱えることがあります。これらの課題に対処するため、著者のプロトコルは高周波分散パラメータケーブルモデル、反射係数スペクトルの導出、最適化されたステップ周波数選択、BurgアルゴリズムによるMES推定を統合しています。従来のFFT(ウィンドウの有無にかかわら...

材料

この記事で使用された材料の一覧

名前	会社	カタログ番号	コメント
導電性銅テープ、両面	詳細は不明です	該当なし	ヘッドエンドの反射干渉を最小限に抑えるために、ケーブルコアとシールドにテストリードを接続するために使用されます https://www.shcenxin.com/metaltape/93.html
YJLV22-8.7/10 kV-3×70中圧ケーブル	詳細は不明です	YJLV22-8.7/10 kV-3×70	架橋ポリエチレン(XLPE)アルミニウムコア鋼テープ装甲PVC被覆電源ケーブル;全長171メートル;外部から誘発される欠陥 http://www.xmdl518.com/jspn01-Products-6234877/ を試験対象として使用
10 kVケーブル欠陥定位装置	詳細は不明です	詳細は不明です	周波数領域反射計測装置;周波数帯100 kHz–50 MHz;反射係数の取得 https://www.163.com/dy/article/K6OHFJOO0552X57P.html のために5,001のサンプリングポイントを配置しています