基于改進ITD-LS的介質(zhì)損耗角檢測方法

徐　娟,孫大偉,亓　亮,張亞鵬

(1.國網(wǎng)寧夏供電公司,寧夏 靈武 750411; 2.國網(wǎng)寧夏電力科學研究院,銀川 750001; 3.國網(wǎng)寧夏經(jīng)濟技術研究院,銀川 750001)

?

基于改進ITD-LS的介質(zhì)損耗角檢測方法

徐娟1,孫大偉2,亓亮2,張亞鵬3

(1.國網(wǎng)寧夏供電公司,寧夏 靈武 750411; 2.國網(wǎng)寧夏電力科學研究院,銀川 750001; 3.國網(wǎng)寧夏經(jīng)濟技術研究院,銀川 750001)

為準確檢測介質(zhì)損耗角,提出基于改進ITD算法與最小二乘法的介質(zhì)損耗角檢測方法。首先構造協(xié)方差矩陣進行數(shù)據(jù)預處理,利用改進ITD法提取基波電壓和電流分量的頻率特性,然后采用最小二乘法獲取基波電壓、電流的相位信息,實現(xiàn)高精度檢測介質(zhì)損耗角。仿真實驗結果表明了在介損角真實值、基波電壓頻率波動、3次諧波比例、直流分量比例、采樣點數(shù)變化時該方法的可行性與有效性。

介質(zhì)損耗角;改進ITD法;最小二乘法

電網(wǎng)中電容型電氣設備占據(jù)很大的比重,而且它們的絕緣狀態(tài)直接關系電網(wǎng)的安全穩(wěn)定運行。因此,準確有效地測試出電氣設備的絕緣狀況對電網(wǎng)的安全運行有著重要的意義[1-2]。目前通常采用介質(zhì)損耗角的檢測方法,該方法包括硬件法和軟件法。硬件方法有西林電橋法和過零點比較法,二者都存在處理環(huán)節(jié)多、累計誤差大、補償困難和抗干擾能力差等缺點。文獻[3]提出了一種用小波變換和正弦波參數(shù)法去除電磁干擾的數(shù)學原理,但若小波基參數(shù)選取不當會導致頻帶混疊;文獻[4-5]采用的FFT法在模態(tài)分解時無法避免模態(tài)混疊;文獻[6-7]中原理復雜,運算量大,無法進行快速檢測。對此,為了準確檢測介質(zhì)損耗角,本文提出一種改進ITD-LS方法。首先應用改進ITD法提取基波電壓和電流分量的頻率,然后利用最小二乘法對基波電壓、電流的相位進行檢測,完成對介質(zhì)損耗角的高精度檢測,并通過仿真實驗驗證了該方法的可行性與有效性。

1　改進ITD法

1.1ITD的基本理論

ITD算法的基本思路是:通過采集實時量測數(shù)據(jù),將量測數(shù)據(jù)等間隔采樣后,構造衰減響應矩陣,利用響應矩陣和其時間延遲矩陣的關系,求取量測信號的廣義特征值[8-9]。該算法以量測數(shù)據(jù)表示成復指數(shù)關系的數(shù)學模型,依據(jù)系統(tǒng)的特征值與振動系統(tǒng)特征值的關系,利用最小二乘法求解其狀態(tài)矩陣,求解出系統(tǒng)的模態(tài)參數(shù)[10]。

自由振動響應微分方程為

{x(t)}N×t=[Φ]N×2N{est}2N×t

{x(t)}=[x1(t),x2(t),…，xN(t)]T

[Φ]=[{Φ1},{Φ2},…，{Φ2N}]

{est}={es1t,es2t，…，es2Nt}

式中:{x(t)}為系統(tǒng)自由振動響應;[Φ]為系統(tǒng)的振型矩陣;Sn為系統(tǒng)第n階特征值;N為系統(tǒng)的模態(tài)階數(shù)。

1.2改進ITD算法

隨機子空間中協(xié)方差矩陣(Toeplitz)具有保留全部原始數(shù)據(jù)信息、去除噪聲干擾、減小運算量等優(yōu)點,故選擇建立協(xié)方差矩陣進行數(shù)據(jù)預處理。首先根據(jù)原數(shù)據(jù)構造Hankel矩陣,然后變換Hankel矩陣構造Toeplitz 矩陣,即

(1)

本文采用單邊最小二乘法替換式(1)中[A]的求法,以進一步降低噪聲對算法的干擾:

對[A]進行特征值分解:

[A][Φ]=[Φ][α]

可以得到矩陣[A]的第r階的特征值為esrΔt,設所求的特征值Vr為

Rr=lnVr=srΔt

(2)

從而求出系統(tǒng)的模態(tài)頻率為

2　最小二乘法

定義電壓、電流的數(shù)學模型為[5]

式中:A0、Ak、fk、θk分別是信號中的直流分量、幅值、頻率、相位,n(t)為噪聲信號。則

mk=Akcosθk

nk=Aksinθk

由公式可以推出基波信號的幅值和相位為

X=[1,cos 2πf1t,…,cos 2πfnt,sin 2πf1,…,sin 2πfnt]

(3)

W=[A0e-λt,M1,…,Mn,N1,…,Nn]T

Y=XW+e

式中,e為誤差向量。

W=(XTX)-1XTY

由參數(shù)W可以求得mk、nk,然后由式(2)求得基波電壓、電流信號的幅值和相位。由式(3)可以求出絕緣電氣設備的介質(zhì)損耗角δ[5],其中φi、φu為基波電流電壓的初相位。

δ=π/2-(φi-φu)

(4)

3改進ITD法和最小二乘法測量介損角步驟

介質(zhì)損耗角檢測流程圖如圖1所示。具體步驟如下:

1) 對采集到的電壓、電流信號進行預處理,構造Hankle矩陣。

2) 由Hankle矩陣來構造Toeplitz矩陣,將Toeplitz矩陣作為改進ITD法的輸入數(shù)據(jù)進行信號的基波頻率參數(shù)提取。

3) 利用最小二乘法得到基波電壓、電流信號的相位φi、φu。

4) 利用式(4)求出介質(zhì)損耗角。

圖1　介質(zhì)損耗角檢測流程圖

4　仿真算例

4.1數(shù)值仿真信號

電氣信號模型為

為了更好地驗證算法的有效性,仿真參數(shù)取值如下:采樣頻率為1 kHz,采樣點數(shù)為1000,基波頻率取為50 Hz,取0.003 rad為介質(zhì)損耗角的實際值。在求解介質(zhì)損耗角時,信號的初始相角選取方式為:在0～π/2之間隨機采取100個數(shù)據(jù)點,取平均值做為初始相角,3次諧波和5次諧波設定為基頻信號的10%、1%。

分別在介損角實際值、基頻頻率波動、3次諧波比例、非周期分量、采樣頻率、采樣個數(shù)以及白噪聲等情況下,對改進ITD-LS算法進行驗證。

4.2介質(zhì)損耗角實際值變化對仿真的影響

當介質(zhì)損耗角實際值從0.004到0.02范圍內(nèi)變動時,求取的結果誤差如表1所示。

表1　真實值對檢測結果的影響

由表1仿真結果可知:當介損角的實際值隨機變動時,誤差絕對值都在10-5rad數(shù)量級,對改進ITD-LS算法的計算精度影響微弱。

4.3信號頻率波動對仿真的影響

設定基波信號頻率取值為49.5、49.7、50、50.3、50.5 Hz,應用改進ITD-LS算法求取介損角,仿真結果如表2所示。

表2　基頻變化對檢測結果的影響

由表2仿真結果可知:當檢測信號中含有諧波成分,信號頻率的波動在電能質(zhì)量允許范圍內(nèi)變化時,取得的介質(zhì)損耗角最大誤差絕對值均低于1×10-5rad,誤差的標準差均低于0.3×10-5rad,檢測信號精度高,且頻率的波動對改進ITD-LS算法影響微弱。

4.43次諧波對仿真的影響

在電氣測量信號中3次諧波占較大比重,檢驗3次諧波在檢測信號的比重變化對算法的影響,仿真結果如表3所示。

表3　3次諧波對檢測結果的影響

由表3仿真結果可知:當調(diào)整3次諧波在檢測信號的比重時,檢測結果誤差絕對值都低于1×10-5rad,誤差的標準差低于0.5×10-5rad,遠遠滿足介質(zhì)損耗角的檢測精度要求,對改進ITD-LS算法檢測結果影響微弱。

4.5非周期分量對仿真的影響

電氣信號的測量中經(jīng)常發(fā)生零漂現(xiàn)象,將會在檢測信號中產(chǎn)生非周期分量。當非周期分量在檢測信號比重變化時,仿真得到的介質(zhì)損耗角誤差情況如表4所示。

表4　非周期分量對檢測結果的影響

由表4仿真結果可知:當檢測信號中非周期分量占比重從0逐漸增加為60%時,檢測結果誤差的絕對值都低于1×10-5rad,誤差標準差都低于0.5×10-5rad,誤差相當微小,對改進ITD-LS算法檢測結果影響微弱。

4.6采樣點數(shù)的影響

設定采樣頻率為1000 Hz,采樣點數(shù)從200點增加到1000點,采樣間隔為200點,仿真結果如表5所示。

表5　采樣點數(shù)對檢測結果的影響

由表5仿真結果可知:采樣點的變化范圍為400點至1000點,計算最大誤差絕對值低于0.2×10-5rad, 標準差都低于10-4rad,滿足介質(zhì)損耗角的誤差要求。為進一步節(jié)省運算時間,采樣400點可獲精確的介損角。

4.7與SSI算法的對比

目前,針對信號檢測方法中,SSI具有較高精度,構造含有噪聲的信號為

y(t)=cos(2πf1t)+0.15cos(2πf2t)+

0.25sin(2πf3t)+n(t)

算法對比結果如表6所示。

表6　算法對比結果

由表6仿真結果可知:在沒有降低檢測精度的前提下,改進ITD算法運算速度約為隨機子空間的50%,表明本方法能夠滿足電氣信號檢測精度要求,且運算速度快。

4.8仿真結果分析

電力行業(yè)標準DLT 1154-2012《高壓電氣設備額定電壓下介質(zhì)損耗因數(shù)測試導則》中對高壓介質(zhì)損耗角的測量精度要求為

ΔD=±arctan[1%×tanD+0.0005]

(5)

式中:ΔD為介損角的測量誤差上限。

將上述各種情況的檢測真實值代入式(5)中求得誤差精度要求,然后與最大檢測誤差對比,其結果如表7所示。

表7　檢測結果與精度要求

由表7可知,基于本文方法的介損角檢測方法誤差比精度要求小兩個數(shù)量級,完全滿足電力行業(yè)標準,由此進一步驗證了本文算法的有效性與準確性。

5　結　論

通過對ITD算法的研究,提出了一種采用Toeplitz矩陣代替Hankel矩陣作為端部數(shù)據(jù)量來提升ITD算法性能的一種高精度的模態(tài)參數(shù)提取方法,并結合LS算法將其運用到介損角的檢測中,成功檢測出基波電壓、電流信號的初相位。仿真結果表明,ITD-LS算法具有較強的免疫力和可靠性,同時該檢測方法的快速性和準確性具有一定的工程實用價值,為電力系統(tǒng)介損角的檢測提供了一種新的思路。

[1] 談克雄,李福祺,張會平.提高電容型設備介損監(jiān)測裝置性能的意見[J]. 高電壓技術, 2002,28(11): 21-23.

TAN Kexiong, LI Fuqi, ZHANG Huiping. Advice for improving performance of on-line tanδ monitoring devices for capacitive type equipment[J]. High Voltage Engineering, 2002,28(11): 21-23.

[2] 趙秀山,談克雄,朱德恒,等.介質(zhì)損耗角的數(shù)字化測量[J].清華大學學報(自然科學版), 1996,36(9):51-56.

ZHAO Xiushan, TAN Kexiong, ZHU Deheng, et al. Study on digital measurements for dielectric loss angle[J].Journal of Tsinghua University (Sci&Tech), 1996,36(9):51-56.

[3] 任佩余, 周健平.計算局部放電電容增量和損耗因數(shù)增量的方法[J].電工技術學報, 1986, 2(2):21-25.

REN Peiyu, ZHOU Jianping. A method for calculation of the capacitance increment and dissipation factor tip-up from parallelogram figures[J]. Transactions of China Electrotechnical Society, 1986, 2(2):21-25.

[4] 蔡國雄,甄為紅,楊曉洪,等. 測量介質(zhì)損耗的數(shù)字化過零點電壓比較法[J]. 電網(wǎng)技術, 2002, 26(7):15-18.CAI Guoxiong, ZHEN Weihong, YANG Xiaohong, et al. A digitized zero cross point voltage method of comparison for measurement of dielectric losses[J]. Power System Technology, 2002, 26(7):15-18.

[5] 柴旭崢,關根志,文習山,等. tanδ高準確度測量的加權插值FFT算法[J]. 高電壓技術, 2003, 29(2):32-33.

CHAI Xuzheng, GUAN Genzhi, WEN Xishan, et al. The hanning windowed interpolated algorithm for dielectric loss measurement with high accuracy[J]. High Voltage Engineering, 2003, 29(2):32-33.

[6] 左自強,徐陽,曹曉瓏,等.計算電容型設備介質(zhì)損耗因數(shù)的相關函數(shù)法的改進[J].電網(wǎng)技術, 2004,28(18):53-57.ZUO Ziqiang, XU Yang, CAO Xiaolong, et al. Improvement to correlation algorithm for dielectric dissipation factor calculation of capacitance equipments[J]. Power System Technology, 2004,28(18):53-57.

[7] 李天云,王靜,郭躍霞,等. 基于經(jīng)驗模態(tài)分解和正弦波參數(shù)法的介損角測量算法[J]. 電網(wǎng)技術,2007,31(4):76-80.

LI Tianyun, WANG Jing, GUO Yuexia, et al. Dielectric loss angle measuring algorithm based on EMD and sin-wave parameter method[J]. Power System Technology, 2007,31(4):76-80.

[8] 楊佑發(fā),李帥,李海龍. 環(huán)境激勵下結構模態(tài)參數(shù)識別的改進ITD法[J]. 振動與沖擊, 2014, 33(1):194-199.YANG Youfa, LI Shuai, LI Hailong. Improved ITD method for structural modal parameter identification under ambient excitation[J]. Journal of Vibration and Shock, 2014, 33(1):194-199.

[9] 王建元,亓亮,張亞鵬,等. 改進ITD-LS結合奇異值差分譜的間諧波快速檢測方法[J]. 電測與儀表,2014,51(20):65-69.

WANG Jianyuan, QI Liang, ZHANG Yapeng, et al. A quicky harmonic detection combined improved ITD-LS with singular value difference spectrum[J]. Electrical Measurement & Instrumentation, 2014,51(20):65-69.

[10] 周幫友,胡紹全,杜強.特征系統(tǒng)實現(xiàn)算法中的模型定階方法研究[J]. 科學技術與工程, 2009,9(10): 2715-2722.ZHOU Bangyou, HU Shaoquan, DU Qiang. Study about calculating the order of model in eigen system realization algorithm[J]. Science Technology and Engineering, 2009,9(10): 2715-2722.

(責任編輯郭金光)

Dielectric loss angle measurement based on improved ITD-LS

XU Juan1, SUN Dawei2, QI Liang2, ZHANG Yapeng3

(1.State Grid Ningxia Power Supply Company, Lingwu,750411,China; 2.State Grid Ningxia Electeic Power Research Institute,Yinchuan,750001,China; 3.State Grid Ningxia Electeic Economic Institute of Technology, Yinchuan,750001,China)

In order to accurately detect dielectric loss, this paper proposed the dielectric loss angle measurement based on the improved ITD algorithm and least square method. It is a method that firstly conforms covariance matrix as data preprocessing, then extracts the frequency characteristics of fundamental wave voltage and current component based on the improved ITD algorithm, and finally acquires their phase information by the least square method so as to realize the dielectric loss angle measurement with high accuracy. The simulation result shows that the method is feasible and effective when the true value of dielectric loss angle, frequency fluctuations of fundamental voltage, the three harmonic ratio, the DC component proportion, and sampling points change.

dielectric loss angle; improved ITD; least square

2016-01-15；

2016-02-25。

徐娟(1982—),女,工程師,從事電網(wǎng)運行工作。

TM854

A

2095-6843(2016)03-0209-05