金增杰，衛(wèi)永琴，張琦，王瀟龍

（1.山東科技大學(xué)電氣與自動(dòng)化工程學(xué)院，山東 青島 266590；2.國家電網(wǎng)山東省電力公司檢修公司，山東 濟(jì)南 272000）

目前，全球用電量增長迅速，我國作為能源大國對電能的利用遠(yuǎn)超其他國家[1]，這就要求我國的輸配電系統(tǒng)具有較高的運(yùn)行可靠性及持續(xù)穩(wěn)定性[2-3]。而輸電線路上的電線電纜有如下優(yōu)點(diǎn)[4]：1）具有外護(hù)層，使得電纜芯線使用壽命長；2）電纜一般埋于地下，不占地上面積，自然也不會(huì)有影響環(huán)境美觀的問題。但同時(shí)也帶來了一些實(shí)際問題，最明顯的是其在地下，走線復(fù)雜曲折，一旦發(fā)生故障，很難定位故障點(diǎn)，并且長時(shí)間的停電檢修，浪費(fèi)人力物力，對社會(huì)生產(chǎn)和經(jīng)濟(jì)都會(huì)造成較大的損失[5-6]。那么如何迅速準(zhǔn)確地找到故障點(diǎn)位置已經(jīng)成為了研究的熱點(diǎn)問題。

故障測距主要是為了快速發(fā)現(xiàn)故障點(diǎn)，這樣工作人員就可以及時(shí)排除故障和修復(fù)線路，所以我們研究的首要問題是如何在最短時(shí)間內(nèi)把故障范圍縮小到最小。現(xiàn)有的電纜故障測距方法主要分為離線和在線測距兩種[7]。離線測距是我國目前主要的電纜故障測距方法[8-9]。而行波法是研究最多的一種離線測距法，已有的故障測距的計(jì)算公式均與故障點(diǎn)相位有關(guān)，不易實(shí)現(xiàn)。本文推導(dǎo)出了易操作的測距計(jì)算公式，采用多頻測相法在Matlab平臺(tái)上仿真并驗(yàn)證。最后設(shè)計(jì)了電纜故障測距硬件系統(tǒng)，完成調(diào)試并通過實(shí)驗(yàn)對系統(tǒng)可行性進(jìn)行了驗(yàn)證。

1 電纜故障行波等效模型

下面研究行波在電纜故障中傳播的等效模型，在工程計(jì)算中，假設(shè)行波傳播速度的經(jīng)驗(yàn)值等于行波的真實(shí)傳播速度[10]。

圖1為行波在故障電纜中傳播的示意圖，電纜兩端點(diǎn)分別用h和i表示，f為電纜中出現(xiàn)的故障點(diǎn)。信號Vh從h端進(jìn)入電纜，當(dāng)其到達(dá)f點(diǎn)時(shí)，便會(huì)產(chǎn)生折射和反射。Vfh為反射信號，從f點(diǎn)往回傳播到端點(diǎn)h；Vfi為折射信號，向前傳播至端點(diǎn)i。

圖1 故障電纜行波傳播示意圖Fig.1 Traveling wave propagation schematic of fault cables

下面建立電纜故障行波傳播的數(shù)學(xué)模型，假設(shè)X1為原始信號，其表達(dá)式為

式中：H為信號幅值；ω為角頻率；t為傳播時(shí)間；θ為初始相位。

測試信號在電纜中傳播時(shí)，到達(dá)故障點(diǎn)產(chǎn)生的反射信號Xfh與電纜首端的輸入信號Xfi疊加，用X2表示如下：

X2(t)=H0sin(ωt+θ)+H1sin(ωt+θ+φ)（2）

式中：H0為注入信號幅值；H1為故障點(diǎn)反射信號的幅值；φ為反射信號Vfh和折射信號Vfi的相角差。

2 推導(dǎo)測距公式及相位差計(jì)算

相位法測距是基于行波測距原理提出的一種間接測距方法，測距基本公式為

式中：λ為正弦信號的波長；N為信號在入射點(diǎn)和故障點(diǎn)之間完整傳播的次數(shù)；Δφ為不足1個(gè)周期的相位差；2π為1個(gè)完整的傳播周期。

多頻測相法是一種依據(jù)電纜總長來逐步增加信號頻率的測距方法，這種方法可以使測距更精確。根據(jù)已有文獻(xiàn)的研究可知，電壓的折射和反射系數(shù)都會(huì)受到頻率大小的影響[11]。下面針對斷路和短路兩種情況，推導(dǎo)出不同的測距公式。因斷路時(shí)故障信號具有正的反射系數(shù)[12]，φ的取值范圍為0～2π，則測距公式可以表示為

短路故障時(shí)，信號的反射系數(shù)為負(fù)[12]，φ的取值范圍為-π～π，測距公式為

定義入射信號為X1，經(jīng)故障反射后的疊加信號為X2，對比信號為X3。令X1與X3頻率相同，均為ω。令X1的初相位為0，則X2，X3兩信號經(jīng)過乘法器運(yùn)算之后，再通過低通濾波器將高頻信號過濾掉，得到輸出信號：

為了計(jì)算簡便，設(shè)定信號X3的初始相位θ分別為0和π/2，得到φ的表達(dá)式為

式中：Xπ/2為θ=π/2時(shí)X（t）的值；X0為θ=0時(shí)X（t）的值。

通過tanφ的值即可求出所需相位差。

3 Matlab建模仿真分析

通過Matlab/Simulink搭建電纜故障模型，首先依據(jù)實(shí)際情況對電纜參數(shù)進(jìn)行設(shè)定，然后通過多頻測相原理和相位法對電纜故障進(jìn)行仿真分析，模型如圖2所示（只進(jìn)行了A，B相電纜檢測）。

圖2 電纜故障模型Fig.2 Cable fault model

圖2中，受控電壓源接收正弦波信號，分別作為檢測的入射信號X1和對比信號X3，X1通過電阻模塊注入電纜，電纜分布參數(shù)設(shè)置如表1所示。

表1 電纜分布參數(shù)Tab.1 The parameters of the cable distribution

3.1 短路故障

首先，設(shè)置電纜總長度S為23 km，A段長度即故障距離為17.65 km，故障設(shè)為接地電阻為10 Ω的單相接地故障。采樣頻率fs=100 MHz，采樣時(shí)間t=0.005 s。信號1和信號2的幅值均為5 V，X1的初相位為0，由此可得入射信號X1，疊加信號X2以及對比信號X3。

然后，1）根據(jù)電纜長度選擇4個(gè)測距信號的頻率f1～f4，要求滿足λ＞2S且留有一定裕度，取波長λ1=50 km，假設(shè)選擇的第一頻率f1=3.993 8 kHz，確定故障距離的十位數(shù)字；2）確定故障距離的個(gè)位數(shù)字，此時(shí)λ2/2=10 km，f2=9.967 kHz；3）確定故障距離小數(shù)點(diǎn)后第1位數(shù)字，此時(shí)λ3/2=1 km，f3=99.67 kHz；4）確定故障距離小數(shù)點(diǎn)后第2位數(shù)字，此時(shí)λ4/2=0.1 km，f4=996.7 kHz。

f1=3.993 8 kHz時(shí)，可以得到θ=0和θ=π/2兩種情況下示波器V4的波形，如圖3所示。

圖3 短路故障f1=3.993 8 kHz時(shí)V4波形Fig.3 V4waveforms at f1=3.993 8 kHz when short circuit fault

圖3中高頻分量被濾掉，只剩直流分量，X0=1.991 V，Xπ/2=-4.528 V，可得tanφ=7.114，相位差φ為-0.544π或0.456π，由于故障反射信號Vf總滯后于由電壓測量元件1測得的發(fā)射信號VM1，因此φ的值為正，取φ=0.456 π，故障類型為短路故障。根據(jù)式（5）可得故障距離S1=18.2 km，因此故障距離十位數(shù)字為1。此時(shí)誤差為0.55 km，相對誤差為2.391%。運(yùn)用同樣的方法可分別計(jì)算出其他3個(gè)頻率測距時(shí)故障距離的個(gè)位、小數(shù)點(diǎn)后第1位和小數(shù)點(diǎn)后第2位數(shù)字，最終得到故障距離S4為17.613 2 km，誤差為-0.036 8 km，相對誤差從第1個(gè)頻率測距時(shí)的2.391%精確到了0.16%。

為了進(jìn)一步驗(yàn)證該方法測距的準(zhǔn)確性，分別設(shè)置故障距離為2.59 km，9.68 km，13.24 km，17.32 km和22.05 km。仿真結(jié)果如表2所示，其中，S為測量的故障距離，e為相對誤差。

表2 短路故障仿真結(jié)果Tab.2 Simulation calculation results when short circuit fault

3.2 斷路故障

將故障類型設(shè)置為A相單相斷路，采樣時(shí)間t=0.005 s，采樣頻率fs=100 MHz。由于故障距離沒有改變，因此四組測量頻率與短路故障時(shí)相同。f1=3.993 8 kHz時(shí)，可以得到θ=0和θ=π/2兩種情況下示波器V4的波形如圖4所示。

圖4 斷路故障f1=3.993 8 kHz時(shí)V4波形Fig.4 V4waveforms at f1=3.993 8 kHz when open circuit fault

由圖4可知，X0=1.996 V，Xπ/2=-4.529 V，同短路故障相同，利用式（8）可以求出tanφ=7.171 8，所以相位差φ=0.456π或 1.456π，因?yàn)閄π/2＜0，而φ∈(π～2π），可得φ的值為1.456π，利用式（4）可得故障距離S1=18.2 km，故障距離十位數(shù)字為1，計(jì)算可得誤差為0.55 km，相對誤差為2.4%。同樣的方法可得到其他3個(gè)頻率測距時(shí)測量故障距離的個(gè)位、小數(shù)點(diǎn)后一位和小數(shù)點(diǎn)后第二位數(shù)字。最終計(jì)算出故障距離為17.726 6 km，誤差為0.076 6 km，相對誤差從第一個(gè)頻率測距時(shí)的2.4%精確到了0.33%。

同樣，分別設(shè)置故障距離2.59 km，9.68 km，13.24 km，17.32 km和22.05 km，仿真計(jì)算結(jié)果如表3所示。

表3 斷路故障仿真結(jié)果Tab.3 Simulation calculation results when open circuit fault

依據(jù)以上結(jié)果可知，當(dāng)電纜故障發(fā)生在不同的距離時(shí)，此方法不受故障類型的影響，能檢測出準(zhǔn)確的故障距離，誤差不超過0.6%。

4 系統(tǒng)設(shè)計(jì)及實(shí)驗(yàn)

4.1 系統(tǒng)總結(jié)構(gòu)

根據(jù)以上的公式推導(dǎo)和仿真結(jié)果理論分析，設(shè)計(jì)了一種多頻測相法電纜故障測距系統(tǒng)，如圖5所示。

圖5 多頻測相法電纜故障測距系統(tǒng)Fig.5 Ranging system with multi-frequency phase method when cable fault

本文提出的多頻測相法電纜故障測距系統(tǒng)工作原理如下：首先將故障數(shù)據(jù)輸入系統(tǒng)，中央處理器檢測信號頻率有無超出臨界值，若低于臨界值便通過信號發(fā)生器發(fā)出兩個(gè)頻率同步的信號，分別輸入電纜首端和乘法器中，進(jìn)入電纜的信號到達(dá)末端再次沿電纜返回后同樣注入乘法器。之后乘法器輸出信號，通過低通濾波器處理該信號變?yōu)橹绷?，最后通過A/D轉(zhuǎn)換器將信號輸入中央處理器。第一頻率完成運(yùn)行計(jì)算之后，中央處理器自主選取下一頻率，繼續(xù)重復(fù)上述過程，最終由上位機(jī)顯示距離數(shù)值。

4.2 系統(tǒng)組裝及實(shí)驗(yàn)測試

系統(tǒng)包括STM32系統(tǒng)、信號發(fā)生器AD9959及2個(gè)運(yùn)算放大器、2個(gè)功率放大器、1個(gè)乘法器和1個(gè)4階巴特沃斯低通濾波器等，組裝焊接并調(diào)試。

下面利用這套測距系統(tǒng)做實(shí)驗(yàn)，實(shí)驗(yàn)中的電纜絕緣材料為聚氯乙烯，選其波速為130 m/μs，將電纜在1 km處斷開，也就是斷路故障。將入射信號的幅值設(shè)置為5 V，并注入電纜；將對比信號的幅值設(shè)置為2.5 V，并注入乘法器。運(yùn)行系統(tǒng)，根據(jù)系統(tǒng)提示向下位機(jī)發(fā)送電纜長度和故障類型（0表示短路故障、1表示斷路故障），隨后系統(tǒng)會(huì)根據(jù)電纜故障信息自動(dòng)選擇頻率進(jìn)行測距，如圖6所示。

圖6 測距結(jié)果Fig.6 Ranging result

第一頻率為66 kHz，分別設(shè)置對比信號的相位θ為0和π/2，用示波器觀測相應(yīng)的濾波器輸出直流信號X0和Xπ/2，如圖7所示.

圖7 第一頻率為66 kHz時(shí)的直流分量Fig.7 The DC components at the first frequency of 66 kHz

圖7中直流分量X0的值為1.93 V，直流分量Xπ/2的值為-0.8 V，可以計(jì)算出X0的初始值為6.655 V，Xπ/2的初始值為-2.759 V，相位差φ為1.843 6π或0.843 6π，最后算得故障距離為0.921 8 km，與系統(tǒng)運(yùn)行結(jié)果一致。

第2次選頻后的直流信號中直流分量X0=1.053 V，直流分量Xπ/2=-0.4 V。經(jīng)計(jì)算可知故障距離為0.990 8 km。

將電纜末端設(shè)置為短路，通過本系統(tǒng)測距得到的故障距離為0.997 6 km。通過以上結(jié)果可知，該測距系統(tǒng)能夠根據(jù)設(shè)計(jì)方案測量故障距離，測得的斷路故障和短路故障相對誤差分別為0.92%和0.24%，符合預(yù)期目標(biāo)。

5 結(jié)論

本文通過多頻相位法對電纜故障測距進(jìn)行研究，推導(dǎo)出短路故障和斷路故障下不同的測距公式，然后用Matlab建立電纜故障模型并利用多頻測相的相位法測量故障距離。

通過仿真結(jié)果分析可知，本文提出的多頻相位對電纜故障測距方法，可測得故障距離，且誤差均小于0.6%，充分論證了此方法的可行性，同時(shí)此測距方法具有不受故障點(diǎn)位置和故障類型的影響的優(yōu)越性。

但是因?yàn)闂l件不足，系統(tǒng)雖能完成測距，但設(shè)計(jì)上依然存在有待改進(jìn)之處：

1）存在部分元件帶寬不足的問題，限制了頻率的選取。

2）系統(tǒng)只根據(jù)線路絕緣材料選擇了一種行波速度，如果通過頻率的差異設(shè)置不同的波速會(huì)進(jìn)一步提高測距的精度。