頻域反射法中阻抗變化類(lèi)型判斷技術(shù)

饒顯杰 周 凱 黃永祿 李明志 李 蓉

頻域反射法中阻抗變化類(lèi)型判斷技術(shù)

饒顯杰 周 凱 黃永祿 李明志 李 蓉

（四川大學(xué)電氣工程學(xué)院 成都 610065）

利用行波的反射可以有效地探測(cè)電纜的局部缺陷和故障等阻抗異變點(diǎn)，反射波的傳播距離和極性可分別用于定位阻抗異變點(diǎn)和判斷阻抗異變點(diǎn)的類(lèi)型。頻域反射（FDR）技術(shù)是識(shí)別電纜中反射波的有效手段之一，但是現(xiàn)有FDR技術(shù)僅能計(jì)算反射波的傳播距離，而不能判斷波形極性，因此，該文提出一種FDR的波形極性判斷技術(shù)。首先闡述FDR中反射波的識(shí)別原理，將Hanning自卷積窗的快速傅里葉變換插值算法用于FDR中反射系數(shù)譜分析，對(duì)FDR測(cè)試下限頻率和頻變波速造成的相位偏差進(jìn)行修正，得到頻率值和修正后的相位值用于計(jì)算反射波的傳播距離和判斷波形極性；然后對(duì)單反射波和多反射波的電纜模型進(jìn)行仿真，檢驗(yàn)該技術(shù)的可行性；最后在實(shí)驗(yàn)室105m長(zhǎng)的10kV交聯(lián)聚乙烯電力電纜上制作阻抗異變點(diǎn)對(duì)該方法的準(zhǔn)確性進(jìn)行驗(yàn)證。仿真和實(shí)驗(yàn)結(jié)果表明，該文所提方法可以有效獲取FDR中反射波的傳播距離和極性。

頻域反射 反射波 極性判斷 Hanning自卷積窗 相位修正

0 引言

近年來(lái)，伴隨著我國(guó)城市的高速擴(kuò)張和發(fā)展，城市用電量逐年攀升，對(duì)城市配電網(wǎng)的供電可靠性提出了更高的要求[1-2]。為了滿足城市發(fā)展的需求，交聯(lián)聚乙烯電力電纜開(kāi)始逐漸在城市配電網(wǎng)的輸電系統(tǒng)中占據(jù)主導(dǎo)地位[3-4]。惡劣復(fù)雜的服役環(huán)境和錯(cuò)誤的施工操作可能破壞電力電纜，迫使其出現(xiàn)局部缺陷，隨著服役年限的增加，局部缺陷可能會(huì)演變?yōu)楣收希瑢?dǎo)致城市電網(wǎng)意外性停電事故，造成經(jīng)濟(jì)損失[5-6]。因此電纜中缺陷或故障的主動(dòng)定位探測(cè)技術(shù)極具研究意義。

電纜中局部缺陷或故障會(huì)迫使電纜局部特性阻抗發(fā)生異常變化，產(chǎn)生阻抗異變點(diǎn)，同時(shí)阻抗異變點(diǎn)會(huì)引起行波在該位置的反射現(xiàn)象，因此大量學(xué)者通過(guò)研究反射波的識(shí)別技術(shù)來(lái)探測(cè)電纜的阻抗異變點(diǎn)。其中，反射波的傳播距離可用于定位阻抗異變點(diǎn)；反射波的極性可用于判斷阻抗異變點(diǎn)的類(lèi)型[7]。

時(shí)域反射（Time Domain Reflectometry, TDR）技術(shù)[8-9]在電纜首端位置注入脈沖信號(hào)，然后直接在時(shí)域中識(shí)別反射波形，由于該方法可以直觀地得到反射波的傳播距離和極性，因此被廣泛用于電纜中反射波的識(shí)別。但是該方法受環(huán)境噪聲影響大，同時(shí)該方法中脈沖信號(hào)的高頻能量少而且高頻能量在電纜中衰減明顯，因此，TDR難以識(shí)別微弱反射波。文獻(xiàn)[10]在TDR基礎(chǔ)上提出了擴(kuò)展頻譜時(shí)域反射法，該方法利用單頻點(diǎn)的反射信息結(jié)合互相關(guān)算法完成反射波的識(shí)別，提高了傳統(tǒng)TDR方法的抗噪能力，但是該方法的頻域信息過(guò)少，因此也難以識(shí)別微弱反射波。

頻域反射（Frequency Domain Reflectometry, FDR）技術(shù)[11-12]中信號(hào)高頻能量更多，頻域信息豐富，能夠識(shí)別更加微弱的反射波，因此逐漸取代了TDR技術(shù)。該方法采用掃頻信號(hào)對(duì)電纜首端進(jìn)行測(cè)試，然后通過(guò)特定算法進(jìn)行時(shí)頻等效變換，最后完成反射波的識(shí)別。文獻(xiàn)[13]提出將FDR中寬頻阻抗譜借助積分變換來(lái)識(shí)別電纜中反射波，并取得較好的識(shí)別效果，但是該方法要求被測(cè)電纜正常時(shí)的多組原始數(shù)據(jù)，因此工程上的實(shí)用性較差。文獻(xiàn)[14-15]提出利用快速傅里葉反變換處理阻抗譜數(shù)據(jù)完成反射波識(shí)別，但是該方法要求掃頻信號(hào)的上限頻率更高，同時(shí)采集點(diǎn)數(shù)也要求更多，因此測(cè)試成本較高；另一方面，測(cè)試頻率過(guò)高時(shí)測(cè)量夾具和電力電纜的連接效果將嚴(yán)重影響測(cè)試結(jié)果。文獻(xiàn)[16]提出將FDR中反射系數(shù)譜在加Kaiser窗后進(jìn)行離散傅里葉變換實(shí)現(xiàn)微弱反射波的識(shí)別，該方法對(duì)測(cè)試頻帶和采集點(diǎn)數(shù)的要求不高，同時(shí)不需要被測(cè)電纜完好時(shí)的測(cè)試數(shù)據(jù)，因此具有更好的應(yīng)用前景，但是該方法存在參數(shù)難以確定和計(jì)算效率低的問(wèn)題。另一方面，目前基于FDR的研究?jī)H能計(jì)算反射波的傳播距離，即定位阻抗異變點(diǎn)，而無(wú)法確定反射波的極性，即無(wú)法確定阻抗異變點(diǎn)的類(lèi)型。

綜上所述，本文在前人研究的基礎(chǔ)上提出一種頻域反射的波形極性判斷技術(shù)，該方法可以利用FDR中反射系數(shù)譜確定反射波的傳播距離和極性，進(jìn)而定位阻抗異變點(diǎn)和判斷阻抗異變點(diǎn)的類(lèi)型。本文首先介紹了FDR的反射波識(shí)別原理，將FDR中反射波的傳播距離計(jì)算問(wèn)題和極性判斷問(wèn)題轉(zhuǎn)化為周期分量的頻率和相位估計(jì)問(wèn)題；接著將Hanning自卷積窗（Hanning Self-Convolution Windows, HSCW）的快速傅里葉變換（Fast Fourier Transform, FFT）插值算法用于估計(jì)周期分量的頻率和相位，以解決原有算法中參數(shù)難以確定和計(jì)算效率低的問(wèn)題，估計(jì)的頻率值可直接用于計(jì)算反射波的傳播距離；最后考慮FDR測(cè)試下限頻率和頻變波速的影響，對(duì)估計(jì)的相位值進(jìn)行修正，利用修正后相位值確定反射波的極性。仿真和實(shí)驗(yàn)驗(yàn)證了該方法的可行性與有效性。

1 FDR的反射波識(shí)別原理

1.1 傳輸線基本理論

式中，為角頻率，f為信號(hào)的頻率；rc為纜心的半徑；rs為金屬屏蔽層的半徑；為纜心的電阻率；為金屬屏蔽層的電阻率；為真空磁導(dǎo)率；為交聯(lián)聚乙烯的電導(dǎo)率；e 為交聯(lián)聚乙烯的介電常數(shù)。

其中

式中，Z為距電纜首端處的等效阻抗；0為特性阻抗。

當(dāng)傳輸信號(hào)的頻率較高時(shí)，0為

1.2 反射波的傳播距離計(jì)算原理

圖2 反射系數(shù)與負(fù)載阻抗關(guān)系

圖3 反射波傳播示意圖

將式（11）代入式（9）得到

借助歐拉公式分解式（12）得到

觀察式（14）可知，當(dāng)注入行波為高頻信號(hào)時(shí)，將看作自變量，反射系數(shù)實(shí)部中產(chǎn)生了頻率為2/的等效周期分量，且該頻率的數(shù)值正好為行波傳播2所需時(shí)間的數(shù)值。該周期分量是由電纜末端的反射現(xiàn)象產(chǎn)生的，因此該周期分量對(duì)應(yīng)電纜末端的反射波，即圖3中的反射波P。

觀察式（15）可知，對(duì)應(yīng)的等效周期分量的頻率為2/，該頻率的數(shù)值正好為行波傳播2所需時(shí)間的數(shù)值。以此類(lèi)推，當(dāng)FDR測(cè)試結(jié)果中存在多個(gè)反射波時(shí)，可以通過(guò)對(duì)反射系數(shù)實(shí)部中各等效周期分量進(jìn)行頻率估計(jì)來(lái)確定FDR中各反射波的傳播距離。

1.3 反射波的極性判斷原理

目前已有的FDR方法僅能計(jì)算反射波的傳播距離，借此確定阻抗異變點(diǎn)位置，但是無(wú)法判斷FDR中反射波極性，即無(wú)法確定阻抗異變點(diǎn)類(lèi)型，由此本文側(cè)重于研究FDR中反射波極性判斷技術(shù)。

利用圖2改寫(xiě)式（15），得到

需要說(shuō)明的是，當(dāng)=0時(shí)，此時(shí)電纜距首端處沒(méi)有阻抗異變點(diǎn)，該位置不存在反射信號(hào)，因此本文不對(duì)該情況進(jìn)行討論。對(duì)式（16）利用三角函數(shù)變換公式進(jìn)行進(jìn)一步改寫(xiě)，得到

式中，為整數(shù)。

綜上所述，借助反射系數(shù)譜的實(shí)部，F(xiàn)DR中反射波的傳播距離計(jì)算問(wèn)題和極性的判斷問(wèn)題轉(zhuǎn)化為周期分量的頻率和相位估計(jì)問(wèn)題。

在實(shí)際工程中進(jìn)行FDR測(cè)試時(shí)，綜合考慮測(cè)試的成本和效率后，會(huì)設(shè)定一個(gè)測(cè)試下限頻率min和測(cè)試上限頻率max，同時(shí)為了準(zhǔn)確地識(shí)別FDR中各反射波，F(xiàn)DR的測(cè)試頻率間隔s應(yīng)滿足奈奎斯特采樣準(zhǔn)則，如式（18）所示，保證其能有效測(cè)試到電纜末端的反射信號(hào)。

2 Hanning自卷積窗的FFT插值算法

本文采用傅里葉變換算法分析反射系數(shù)實(shí)部的周期分量，由于實(shí)際采集的信號(hào)存在非同步采樣和數(shù)據(jù)截?cái)嗟默F(xiàn)象，用傅里葉變換算法進(jìn)行周期分量識(shí)別時(shí)會(huì)產(chǎn)生頻譜泄漏和柵欄效應(yīng)[21]，影響FDR中反射波的辨識(shí)能力。因此需要進(jìn)行離散頻譜的校正，本文采用FFT插值[22]的方法來(lái)校正離散頻譜，該方法可以在頻譜上增加數(shù)據(jù)的頻率分辨率，進(jìn)而提高周期分量的頻率和相位計(jì)算精度。對(duì)FDR中反射波的識(shí)別而言，該方法可以在有效的測(cè)試頻帶中，更好地確定反射波形的特征。同時(shí)借助FFT插值算法可以有效提升計(jì)算效率。

由于本文需要精確估計(jì)反射系數(shù)實(shí)部中各周期分量的頻率和相位，且要提升對(duì)微弱反射波的識(shí)別能力，因此需降低各周期分量的譜間干擾。當(dāng)增加卷積階數(shù)時(shí)，Hanning自卷積窗[23]的旁瓣衰減速率快速增加，此時(shí)各周期分量的相互干擾能夠得到有效抑制，頻譜泄漏現(xiàn)象得到改善。因此本文利用該技術(shù)處理反射系數(shù)譜以獲得更好的反射波識(shí)別效果。

Hanning窗離散形式下可表示[23]為

式中，為Hanning窗的窗長(zhǎng)。

把個(gè)Hanning窗自卷積得到階HSCW為

當(dāng)構(gòu)造窗長(zhǎng)為的階HSCW時(shí)，各Hanning窗的為

式中，floor為向下取整。

利用式（20）和式（21）構(gòu)造HSCW后，在序列末尾補(bǔ)零便得到階HSCW。圖4展示了=126的1～3階HSCW幅頻響應(yīng)曲線，當(dāng)卷積階數(shù)增加時(shí)，HSCW的旁瓣電平在減小，旁瓣衰減速率在增加，這說(shuō)明HSCW可以有效抑制頻譜泄漏，減少譜間干擾，因此將其用于反射系數(shù)譜處理有助于精確估計(jì)各周期分量的頻率和相位。當(dāng)選取太大時(shí)，不僅計(jì)算復(fù)雜，同時(shí)主瓣寬度會(huì)過(guò)大，不適合頻譜密集情況下的識(shí)別，因此本文采用2階HSCW，以解決原有Kaiser窗[16]中需要確定參數(shù)的難題。

圖4 HSCW的旁瓣性能

3 FDR的波形極性判斷技術(shù)

3.1 FDR測(cè)試下限頻率的相位修正

圖5 q的相位函數(shù)波形

3.2 頻變波速的相位修正

當(dāng)傳輸信號(hào)的頻率變化時(shí)，電纜單位長(zhǎng)度的電容、電感參數(shù)值會(huì)產(chǎn)生不同程度的變化，導(dǎo)致不同頻率的波速略有不同，即電纜中的色散現(xiàn)象[18]。圖6為某型號(hào)10kV電纜的實(shí)測(cè)波速，可以看出，電纜波速在較低頻段區(qū)間變化較大，這種變化會(huì)導(dǎo)致反射系數(shù)實(shí)部中各周期分量發(fā)生相位偏移，為了直觀說(shuō)明這種現(xiàn)象，構(gòu)建周期分量1為

式中，為圖6所示實(shí)測(cè)電纜相速度。f=0.01～80MHz，f的間隔為0.01MHz。

圖7 y1和y2的對(duì)比結(jié)果

3.3 初始相位判斷波形極性

結(jié)合式（17）可以得到，當(dāng)=0°時(shí)，該反射波形為正極性；當(dāng)=180°時(shí)，該反射波形為負(fù)極性。在實(shí)際測(cè)試中，綜合考慮測(cè)量誤差和計(jì)算誤差的影響，提出FDR的波形極性判斷條件如圖8所示。當(dāng)處于一、四象限，實(shí)軸正半軸或虛軸正半軸時(shí)，判斷該反射波為正極性；當(dāng)處于二、三象限，實(shí)軸負(fù)半軸或虛軸負(fù)半軸時(shí)，判斷該反射波為負(fù)極性。

4 仿真算例

為了驗(yàn)證本文所提方法判斷FDR中反射波極性的準(zhǔn)確性，本文分別對(duì)單反射波和多反射波兩類(lèi)情況進(jìn)行驗(yàn)證。需要說(shuō)明的是，本文仿真僅討論傳播距離小于電纜末端反射波的反射波。本文采用Matlab軟件仿真10kV的交聯(lián)聚乙烯電纜，具體的電纜參數(shù)[18-19]見(jiàn)表1。

表1 電纜模型參數(shù)

Tab.1 The parameters of cable model

4.1 單反射波

本文采用常見(jiàn)的過(guò)渡電阻接地故障模擬電纜中單反射波的情況，其模型如圖9所示。圖9中，g為過(guò)渡電阻，本文分別將其設(shè)置為20W、200W和2kW。FDR中頻段選擇為0.15～80MHz，s= 0.005MHz。利用本文算法識(shí)別FDR中反射波如圖10所示，得到極性判斷結(jié)果見(jiàn)表2。從圖10可以看出，F(xiàn)DR的反射波識(shí)別可以確定電纜中反射波的傳播距離，但是無(wú)法直接判斷其極性。表2的數(shù)據(jù)證明了本文方法判斷FDR中反射波極性的準(zhǔn)確性。

對(duì)圖9中電纜模型進(jìn)行TDR測(cè)試分析，得到其反射波識(shí)別結(jié)果如圖11所示。從圖11可以看出，當(dāng)g=20W時(shí)，接地故障位置阻抗變化程度較大，1反射波能量較強(qiáng)，TDR能夠成功識(shí)別該反射波并判斷其極性。當(dāng)g增加時(shí)，接地故障位置阻抗變化程度變小，1反射波能量變?nèi)酰琓DR中1的識(shí)別效果變差。當(dāng)g=2kW時(shí)，TDR已經(jīng)無(wú)法識(shí)別1，而FDR仍能夠準(zhǔn)確識(shí)別1。因此對(duì)于接地電阻過(guò)高的高阻接地而言，F(xiàn)DR可以準(zhǔn)確地識(shí)別該反射波，證明了FDR相比于TDR能更好地識(shí)別電纜中微弱的反射波。

圖9 單反射波的電纜模型

圖10 單反射波的FDR識(shí)別結(jié)果

表2 單反射波的極性判斷結(jié)果

Tab.2 The polarity judgment results of single reflected wave

圖11 單反射波的TDR識(shí)別結(jié)果

4.2 多反射波

本文通過(guò)改變電纜中一段區(qū)域的特性阻抗來(lái)制作電纜中多反射波的情況，其模型如圖12所示，圖中，阻抗變化段的特性阻抗設(shè)置為原來(lái)的1.005倍。FDR中頻段選擇為0.15～80MHz，s=0.005MHz。利用本文算法識(shí)別該模型的FDR反射波如圖13所示，得到反射波極性判斷結(jié)果，見(jiàn)表3。圖13同樣說(shuō)明了FDR的波形識(shí)別圖雖然可以識(shí)別反射波1、2和3并判斷其傳播距離，但是無(wú)法判斷其極性。表3表明即使存在多個(gè)反射波時(shí)，本文方法也可以有效判斷FDR中反射波極性。

圖12 多反射波的電纜模型

圖13 多反射波的FDR識(shí)別結(jié)果

表3 多反射波的極性判斷結(jié)果

Tab.3 The polarity judgment results of multiple reflected waves

對(duì)圖12中的電纜模型進(jìn)行TDR測(cè)試，得到TDR測(cè)試結(jié)果如圖14所示。對(duì)比圖13和圖14可以看出，對(duì)于反射能量較強(qiáng)的反射波1和2，TDR和FDR均能有效識(shí)別，并且都能準(zhǔn)確地判斷其波形極性；但是對(duì)于反射能量較弱的反射波3而言，F(xiàn)DR仍可以有效識(shí)別并準(zhǔn)確地判斷其波形極性，但是TDR無(wú)法進(jìn)行識(shí)別，該仿真同樣證明了FDR能夠有效識(shí)別電纜中微弱的反射波。

圖14 多反射波的TDR識(shí)別結(jié)果

5 實(shí)驗(yàn)驗(yàn)證

為了實(shí)際驗(yàn)證本文所提方法判斷FDR中反射波極性的有效性，在實(shí)驗(yàn)室中將105m 10kV交聯(lián)聚乙烯電力電纜在距首端43m處進(jìn)行鹽水腐蝕破壞。在進(jìn)行一段時(shí)間的腐蝕破壞后，該位置的銅屏蔽層出現(xiàn)了銅綠，因此判斷該電纜中出現(xiàn)了阻抗異變點(diǎn)，對(duì)該電纜進(jìn)行FDR測(cè)試，如圖15所示。

圖15 FDR的測(cè)試系統(tǒng)示意圖

圖15中，在控制PC的控制下，調(diào)頻信號(hào)發(fā)生器產(chǎn)生調(diào)頻信號(hào)，該信號(hào)借助功分器實(shí)現(xiàn)信號(hào)能量的平分，平分后的一半信號(hào)成為入射信號(hào)經(jīng)過(guò)耦合器之后傳播到電纜中；平分后的另一半信號(hào)成為參考信號(hào)直接被數(shù)據(jù)采集模塊捕獲。當(dāng)入射信號(hào)傳播到電纜中后，會(huì)引發(fā)電纜的行波反射現(xiàn)象，進(jìn)而產(chǎn)生反射信號(hào)，反射信號(hào)會(huì)由電纜向耦合器方向傳播，在耦合器的作用下反射信號(hào)被捕獲到數(shù)據(jù)采集裝置中。此時(shí)數(shù)據(jù)采集裝置同時(shí)采集得到參考信號(hào)和反射信號(hào)，計(jì)算得到反射系數(shù)譜并將其返回給控制PC，進(jìn)而完成FDR中反射系數(shù)譜的測(cè)試。

圖16 實(shí)測(cè)的FDR識(shí)別結(jié)果

圖17 實(shí)測(cè)的TDR識(shí)別結(jié)果

6 結(jié)論

本文針對(duì)傳統(tǒng)FDR方法無(wú)法判斷反射波極性，即無(wú)法判斷阻抗異變點(diǎn)類(lèi)型的問(wèn)題，提出了一種新的FDR反射波極性判斷方法，該方法首先對(duì)FDR中反射系數(shù)實(shí)部進(jìn)行周期分量分析，利用其頻率值和修正后的相位值分別計(jì)算反射波的傳播距離和判斷其極性。通過(guò)對(duì)仿真和實(shí)測(cè)分析，得到以下結(jié)論：

1）本文通過(guò)分析傳輸線理論，證明了利用FDR中反射系數(shù)實(shí)部進(jìn)行反射波識(shí)別的可行性。

2）HSCW的FFT插值算法改善了原有FDR方法中參數(shù)難以確定和計(jì)算效率低的問(wèn)題，同時(shí)其周期分量的譜間干擾抑制效果較好，利于FDR的反射波識(shí)別。

3）在利用相位判斷反射波極性之前，需要對(duì)FDR測(cè)試下限頻率和頻變波速造成的相位偏差進(jìn)行修正，以提高極性判斷的準(zhǔn)確率。

4）仿真和實(shí)際測(cè)試結(jié)果表明，本文所提方法可以有效計(jì)算FDR中反射波的傳播距離和判斷其極性，可以進(jìn)一步定位電纜局部缺陷和故障等阻抗異變點(diǎn)并判斷其類(lèi)型。

[1] 周凱, 謝敏, 趙世林, 等. 基于改進(jìn)FastICA的局部放電在線監(jiān)測(cè)窄帶干擾高保真性抑制方法[J]. 電工技術(shù)學(xué)報(bào), 2018, 33(11): 2604-2612.

Zhou Kai, Xie Min, Zhao Shilin, et al. Periodic narrowband noise rejection with high fidelity of partial discharge on line monitoring based on improved FastICA algorithm[J]. Transactions of China Elec- trotechnical Society, 2018, 33(11): 2604-2612.

[2] 饒顯杰, 周凱, 謝敏, 等. 穩(wěn)定圖法在極化等效電路參數(shù)辨識(shí)中的應(yīng)用[J]. 電工技術(shù)學(xué)報(bào), 2020, 35(10): 2248-2256.

Rao Xianjie, Zhou Kai, Xie Min, et al. Application of stabilization diagram method for solving polarization equivalent circuit parameters[J]. Transactions of China Electrotechnical Society, 2020, 35(10): 2248- 2256.

[3] 李長(zhǎng)明, 伍國(guó)方, 李春陽(yáng), 等. XLPE絕緣高壓直流電纜終端內(nèi)缺陷對(duì)電場(chǎng)分布的影響[J]. 電機(jī)與控制學(xué)報(bào), 2018, 22(12): 62-67.

Li Changming, Wu Guofang, Li Chunyang, et al. Effect of the defects inside XLPE insulated HVDC cable termination on the electric field distribution[J].

Electric Machines and Control, 2018, 22(12): 62-67.

[4] 樂(lè)彥杰, 汪洋, 鄭新龍, 等. 交聯(lián)聚乙烯與浸漬紙絕緣直流電纜接頭電場(chǎng)分布[J]. 電機(jī)與控制學(xué)報(bào), 2019, 23(2): 75-86.

Le Yanjie, Wang Yang, Zheng Xinlong, et al. Electric field distributions in the joint of XLPE and mass impregnated HVDC cables[J]. Electric Machines and Control, 2019, 23(2): 75-86.

[5] 李露露, 雍靜, 曾禮強(qiáng), 等. 基于系統(tǒng)電力擾動(dòng)的交叉互聯(lián)電纜絕緣整體老化在線監(jiān)測(cè)[J]. 電工技術(shù)學(xué)報(bào), 2018, 33(14): 3396-3405.

Li Lulu, Yong Jing, Zeng Liqiang, et al. On-line monitoring of insulation overall aging for cross- bonded cables based on system power disturbances[J]. Transactions of China Electrotechnical Society, 2018, 33(14): 3396-3405.

[6] 秦雪, 錢(qián)勇, 許永鵬, 等. 基于2D-LPEWT的特征提取方法在電纜局部放電分析中的應(yīng)用[J]. 電工技術(shù)學(xué)報(bào), 2019, 34(1): 170-178.

Qin Xue, Qian Yong, Xu Yongpeng, et al. Appli- cation of feature extraction method based on 2D- LPEWT in cable partial discharge analysis[J]. Transactions of China Electrotechnical Society, 2019, 34(1): 170- 178.

[7] 李露. 船用電力電纜故障診斷的仿真研究[D]. 大連: 大連理工大學(xué), 2017.

[8] 李保生. 基于時(shí)域脈沖反射原理的電線電纜精確測(cè)長(zhǎng)技術(shù)研究[D]. 西安: 西安電子科技大學(xué), 2010.

[9] 張俊民, 魏娟, 謝華博, 等. 基于時(shí)域反射法的航空導(dǎo)線絕緣故障檢測(cè)與分析[J]. 航空學(xué)報(bào), 2009, 30(4): 706-712.

Zhang Junmin, Wei Juan, Xie Huabo, et al. Detection and analysis of aerospace wire insulation faults based on TDR[J]. Acta Aeronouticaet Astronautica Sinica, 2009, 30(4): 706-712.

[10] 吳春華, 胡雅, 李智華, 等. 基于SSTDR的光伏系統(tǒng)直流母線電弧故障在線檢測(cè)與定位[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2020, 40(8): 2725-2735.

Wu Chunhua, Hu Ya, Li Zhihua, et al. On-line detection and location of DC bus arc faults in PV systems based on SSTDR[J]. Proceedings of the CSEE, 2020, 40(8): 2725-2735.

[11] Shi Qinghai, Kanoun O. Wire fault diagnosis in the frequency domain by impedance spectroscopy[J]. IEEE Transactions on Instrumentation and Measure- ment, 2015, 64(8): 2179-2187.

[12] Hirai N, Yamada T, Ohki Y. Comparison of broad- band impedance spectroscopy and time domain refle- ctometry for locating cable degradation[C]//IEEE International Conference Condition Monitoring and Diagnosis, Bali, Indonesia, 2012: 229-232.

[13] Zhou Zhiqiang, Zhang Dandan, He Junjia, et al. Local degradation diagnosis for cable insulation based on broadband impedance spectroscopy[J]. IEEE Transa- ctions on Dielectrics and Electrical Insulation, 2015, 22(4): 2097-2107.

[14] Ohki Y, Yamada T, Hirai N. Precise location of the excessive temperature points in polymer insulated cables[J]. IEEE Transactions on Dielectrics and Electrical Insulation, 2013, 20(6): 2099-2106.

[15] Ohki Y, Yamada T, Hirai N. Diagnosis of cable aging by broadband impedance spectroscopy[C]//IEEE Conference on Electrical Insulation and Dielectric Phenomena, Cancun, Mexico, 2011: 24-27.

[16] 謝敏, 周凱, 趙世林, 等. 新型基于反射系數(shù)譜的電力電纜局部缺陷定位方法[J]. 電網(wǎng)技術(shù), 2017, 41(9): 3083-3089.

Xie Min, Zhou Kai, Zhao Shilin, et al. A new location method of local defects in power cables based on reflection coefficient spectrum[J]. Power System Technology, 2017, 41(9): 3083-3089.

[17] 周志強(qiáng). 基于寬頻阻抗譜的電纜局部缺陷診斷方法研究[D]. 武漢: 華中科技大學(xué), 2015.

[18] 謝敏, 周凱, 趙世林, 等. 考慮相速度頻變效應(yīng)的改進(jìn)互相關(guān)算法局部放電定位[J]. 電網(wǎng)技術(shù), 2018, 42(5): 1661-1667.

Xie Min, Zhou Kai, Zhao Shilin, et al. Research on partial discharge location using modified cross- correlation method considering frequency characteri- stic of phase velocity[J]. Power System Technology, 2018, 42(5): 1661-1667.

[19] 謝敏, 周凱, 何珉, 等. 基于時(shí)間反演技術(shù)的電力電纜局部放電定位方法[J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2018, 38(11): 3402-3409.

Xie Min, Zhou Kai, He Min, et al. Partial discharge location for power cable based on the time reversal technique[J]. Proceedings of the CSEE, 2018, 38(11): 3402-3409.

[20] 洪博, 王莉, 毛健美, 等. 一種航天器一次母線故障在線檢測(cè)定位方法[J]. 電工技術(shù)學(xué)報(bào), 2016, 31(5): 118-125.

Hong Bo, Wang Li, Mao Jianmei, et al. An on-line detection and location method for primary bus fault of spacecraft[J]. Transactions of China Electro- technical Society, 2016, 31(5): 118-125.

[21] 王永, 滕召勝, 李建閩, 等. 基于采樣序列重構(gòu)的高精度介質(zhì)損耗角測(cè)量方法[J]. 電工技術(shù)學(xué)報(bào), 2018, 33(23): 5607-5615.

Wang Yong, Teng Zhaosheng, Li Jianmin, et al. High accuracy dielectric loss angle measurement method based on sampling sequence reconstruction[J]. Transactions of China Electrotechnical Society, 2018, 33(23): 5607-5615.

[22] 段虎明, 秦樹(shù)人, 李寧. 離散頻譜的校正方法綜述[J]. 振動(dòng)與沖擊, 2007(11): 138-145, 189.

Duan Huming, Qin Shuren, Li Ning. Review of correction methods for discrete spectrum[J]. Journal of Vibration and Shock, 2007(11): 138-145, 189.

[23] 溫和, 滕召勝, 卿柏元. Hanning自卷積窗及其在諧波分析中的應(yīng)用[J]. 電工技術(shù)學(xué)報(bào), 2009, 24(2): 164-169.

Wen He, Teng Zhaosheng, Qing Baiyuan. Hanning self-convolution windows and its application to harmonic analysis[J]. Transactions of China Electro- technical Society, 2009, 24(2): 164-169.

Type Judgement Technology of Impedance Variation in Frequency Domain Reflection Method

（School of Electrical Engineering Sichuan University Chengdu 610065 China）

The reflection of the traveling wave can be used to detect the impedance variation of cable including local defects and faults. The propagation distance and polarity of the reflected wave are respectively used to locate the impedance variation position and judge the type of the impedance variation. Frequency domain reflection (FDR) is one of the effective means to identify the reflected wave in the cable. However, the traditional FDR method can only obtain the propagation distance of the reflected wave, but cannot judge the polarity of the reflected wave. Therefore, this paper presents a polarity judgement technology for the reflected wave in FDR. Firstly, this paper describes the principle of reflected wave identification in FDR, and analyzes the reflection coefficient in FDR by the fast Fourier transform interpolation algorithm of Hanning self-convolution windows. The phase deviations caused by lower limit frequency of FDR and frequency-dependent wave velocity are corrected, and the frequency value and the modified phase value are respectively used to obtain the propagation distance and the polarity of the reflected wave. Then, the cable model of single reflected wave and multiple reflected wave is simulated. Finally, the proposed method is verified by experiments on a 105m 10kV cross-linked polyethylene cable. The results show that the method proposed can effectively obtain the propagation distance and polarity of the reflected wave in FDR.

Frequency domain reflectometry (FDR), reflected wave, polarity judgement, Hanning self-convolution windows, phase correction

10.19595/j.cnki.1000-6753.tces.200463

TM41

饒顯杰 男，1996年生，碩士，研究方向?yàn)殡娏υO(shè)備狀態(tài)監(jiān)測(cè)。E-mail: 1953520244@qq.com

周 凱 男，1975年生，教授，博士生導(dǎo)師，研究方向?yàn)殡娎|絕緣狀態(tài)檢測(cè)與修復(fù)等。E-mail: zhoukai_scu@163.com（通信作者）

國(guó)家自然科學(xué)基金資助項(xiàng)目（51477106）。

2020-05-06

2020-06-16

（編輯 崔文靜）