基于多元線性回歸的阻性和容性電流分解

韓永森， 李忠華， 鄭歡， 郭文敏

(哈爾濱理工大學(xué) 教育部工程電介質(zhì)及其應(yīng)用技術(shù)重點實驗室，黑龍江 哈爾濱 150080)

?

基于多元線性回歸的阻性和容性電流分解

韓永森， 李忠華， 鄭歡， 郭文敏

(哈爾濱理工大學(xué) 教育部工程電介質(zhì)及其應(yīng)用技術(shù)重點實驗室，黑龍江 哈爾濱 150080)

為了研究交流電壓作用下非線性半導(dǎo)體器件和非線性絕緣電介質(zhì)的絕緣狀態(tài)和介電性能，提出一種阻性和容性電流分解算法。以非線性電阻和非線性電容構(gòu)成的并聯(lián)等效電路為研究對象，推導(dǎo)響應(yīng)電流關(guān)于激勵電壓的非線性方程。通過坐標變換，將其轉(zhuǎn)化成多元線性方程。利用多元線性回歸方法，獲得等效電路參數(shù)且實現(xiàn)了阻性和容性電流的分解。定性分析該算法的抗干擾能力和對非標準正弦波電壓的適應(yīng)能力。仿真結(jié)果表明：該算法可以準確地實現(xiàn)阻性和容性電流的分解；當響應(yīng)電流含有55 dB的噪聲時，電路參數(shù)的求解誤差較小；激勵電壓諧波分量對電路參數(shù)的求解幾乎沒有影響。實驗結(jié)果表明：該算法可以實現(xiàn)MOA閥片在交流電壓作用下全泄露電流的分解。

電流分解；非線性等效電路；電路參數(shù)估計；多元線性回歸

0 引 言

阻性和容性電流的分解是非線性半導(dǎo)體器件和絕緣電介質(zhì)在交流電壓作用下的狀態(tài)監(jiān)測和電氣性能測試的關(guān)鍵技術(shù)環(huán)節(jié)。

金屬氧化物避雷器(MOA)以其優(yōu)異的非線性特性和良好的浪涌吸收能力已逐漸廣泛應(yīng)用于電力系統(tǒng)的過電壓保護中[1-2]。但是，MOA在運行期間會受到各種過電壓以及外界環(huán)境因素(污濁、潮濕等)的影響而使其電氣性能劣化，主要表現(xiàn)在泄漏電流中的阻性電流分量明顯增大，同時功率損耗的增加會進一步加劇MOA的老化，甚至?xí)贡芾灼鲹p壞[3]。因此，MOA阻性電流分量是MOA狀態(tài)監(jiān)測與診斷的重要依據(jù)。

MOA阻性電流提取的方法主要有：基波法[4]、容性電流補償法[5]、三次諧波法[6]、相位比較法[7-8]等。這些方法均是頻域分析法，是以標準正弦波激勵電壓為基本假設(shè)條件的。

非線性絕緣電介質(zhì)系指其介電參數(shù)隨電場變化而變化的絕緣材料[9]，具有在不均勻電場中自行調(diào)控電場分布的特點[10]。應(yīng)用于電纜終端的電場應(yīng)力控制片(或管)[11]和高壓電機定子線棒端部防電暈結(jié)構(gòu)中的防電暈帶[12]就屬于非線性絕緣電介質(zhì)，其在高壓直流電纜和套管中也具有潛在的應(yīng)用前景[13-14]。

采用研究線性絕緣電介質(zhì)的測試設(shè)備(比如，交流高壓電橋、阻抗分析儀和介電譜儀等)測量非線性絕緣電介質(zhì)，測試結(jié)果不能真實反映其介電特性[15]。國內(nèi)學(xué)者通過仿真和實驗手段研究了不同復(fù)合絕緣電介質(zhì)的直流非線性電導(dǎo)特性[16-18]。國外Obrzut等對非線性絕緣電介質(zhì)的交流激勵電壓和響應(yīng)電流波形進行了測試，并用響應(yīng)電流的諧波分量表征了非線性絕緣電介質(zhì)的非線性特性[19]。文獻[20]利用非線性最小二乘算法實現(xiàn)了非線性絕緣電介質(zhì)的阻性電流和容性電流的分解。

交聯(lián)聚乙烯(XLPE)電力電纜已經(jīng)成為城市供電網(wǎng)絡(luò)的主體，為確保城市電網(wǎng)供電安全，電纜水樹老化狀態(tài)的診斷技術(shù)備受關(guān)注。XLPE電纜水樹老化后，由于水分和低分子產(chǎn)物的存在使電纜絕緣表現(xiàn)出場致增強型電導(dǎo)特性，標準正弦波電壓作用下其阻性電流將含有諧波分量，因此阻性電流及其諧波的測量[21]已成為XLPE電纜絕緣狀態(tài)診斷的有效途徑之一。

綜上所述，準確地實現(xiàn)在交流電壓作用下的非線性半導(dǎo)體器件和絕緣電介質(zhì)的阻性電流和容性電流的分解是研究器件和絕緣電介質(zhì)的電氣性能的關(guān)鍵。本文建立了非線性絕緣電介質(zhì)的并聯(lián)等效電路模型，得到了響應(yīng)電流關(guān)于激勵電壓的時域表達式，利用多元線性回歸方法得到了等效電路參數(shù)，實現(xiàn)了阻性和容性電流的分解。仿真和實驗驗證了該方法的可行性和實用性。

1 理論基礎(chǔ)

1.1 非線性電路模型

為了實現(xiàn)普遍意義上的阻性電流和容性電流的分解，以非線性絕緣電介質(zhì)為例進行分析。參照線性絕緣電介質(zhì)的表示方式，用非線性電導(dǎo)與非線性電容的并聯(lián)等效電路來描述非線性絕緣介質(zhì)，如圖1所示。其中，u(t)是施加在此介質(zhì)上的交流激勵電壓；i(t)是流過此介質(zhì)的響應(yīng)電流；iG(t)是流過此介質(zhì)的阻性電流；iC(t)是流過此介質(zhì)的容性電流；Gx(u)是此介質(zhì)的非線性等效電導(dǎo)；Cx(u)是此介質(zhì)的非線性等效電容，t為時間。

圖1 非線性絕緣電介質(zhì)的并聯(lián)等效電路Fig.1 Parallel equivalent circuit for nonlinear insulating dielectrics

用電壓u(t)的多項式形式來表示非線性電導(dǎo)Gx(u)和非線性電容Cx(u)?？紤]到非線性絕緣電介質(zhì)中不存在整流和濾波過程，電流中只會含有奇次諧波分量，因此多項式中只采用了u(t)的偶次項。

非線性電導(dǎo)Gx(u)與電壓u(t)的關(guān)系為

(1)

式中，n=0，1，2，…；Gx,2k為Gx(u)對應(yīng)的第2k個電路參數(shù)。

Gx(u)在u(t)作用下產(chǎn)生的阻性電流iG(t)為

(2)

非線性電容Cx(u)與電壓u(t)的關(guān)系為

(3)

式中，n=0，1，2，…；Cx,2k為Cx(u)對應(yīng)的第2k個電路參數(shù)。

Cx(u)在u(t)作用下產(chǎn)生的容性電流iC(t)為

(4)

響應(yīng)電流i(t)為

i(t)= iG(t)+iC(t)=

(5)

由此可知，響應(yīng)電流i(t)是關(guān)于交流激勵電壓u(t)的非線性方程，共包含2n+2個電路參數(shù)。

1.2 阻性、容性電流的分解

由式(2)和式(4)可知，上述2n+2個電路參數(shù)是實現(xiàn)阻性電流和容性電流分解的關(guān)鍵。這些電路參數(shù)通常都是未知的，需要通過求解非線性方程(5)來得到這些參數(shù)。然而，非線性方程的求解算法復(fù)雜，且收斂結(jié)果與初始值的選擇有關(guān)。為了方便地獲得準確的電路參數(shù)，本文首先將非線性方程轉(zhuǎn)化成多元線性方程，然后用多元線性回歸方法求解上述電路參數(shù)。

令

y=i(t)，

(6)

xk=u(t)2k+1，

(7)

bk=Gx,2k。

(8)

式中，k=0，1，…，n。

(9)

bk=Cx,2(k-n-1)。

(10)

式中，k=n+1，n+2，…，2n+1。

根據(jù)式(6)～式(10)便可以將非線性方程(5)轉(zhuǎn)化為多元線性方程

(11)

通過實驗或仿真可以得到m組觀測數(shù)據(jù)，為

xj,0，xj,1，…，xj,2n+1，yj。

其中：xj,0為自變量x0的第j個觀測數(shù)據(jù)；xj,1為自變量x1的第j個觀測數(shù)據(jù)；…；xj, 2n+1為自變量x2n+1的第j個觀測數(shù)據(jù)；yj為因變量y的第j個觀測數(shù)據(jù)；j=1，2，…，m。

以此m組觀測數(shù)據(jù)和多元線性方程(11)為基礎(chǔ)，用多元線性回歸方法[22]得到正則方程組，即

(XTX)B=XTY。

(12)

式中，Y=(y1，y2，…，ym)T；B=(b0，b1，…，b2n+1)T；

考慮到方程組求解的穩(wěn)定性和運算速度，用QR分解法[23]求解此正則方程組，得到回歸系數(shù)向量B，即電路參數(shù)向量。

為了便于分析，在公式(1)和式(3)中采用相同數(shù)量的電路參數(shù)來描述非線性電導(dǎo)Gx(u)和非線性電容Cx(u)。利用F檢驗可以對各電路參數(shù)bk的顯著程度進行定量的判斷。對于作用不顯著的電路參數(shù)bk應(yīng)舍去，此時F檢驗滿足：

(13)

對于給定的顯著水平α，若Fk>F1-α，則認為電路參數(shù)bk作用顯著。將作用顯著的各電路參數(shù)代入到公式(2)和式(4)便可以得到阻性電流和容性電流。

1.3 分解算法的抗干擾能力和對非標準正弦波電壓的適應(yīng)能力分析

通過電壓互感器、電容分壓器或電阻分壓器測得的電壓u(t)的相對誤差通常較小，此測試誤差幾乎不影響阻性、容性電流的準確分解。但是，由于響應(yīng)電流i(t)的數(shù)量級較小(mA級以下)，在測試過程中響應(yīng)電流i(t)容易受到外界電磁場的干擾而混有一定信噪比SNR的白噪聲。信噪比SNR可按下式計算：

(14)

式中，PS為信號能量；PN為噪聲能量。

相對于響應(yīng)電流i(t)而言，白噪聲具有較寬的頻率范圍(包含響應(yīng)電流的頻率)。當白噪聲中與響應(yīng)電流不相關(guān)的頻率的振幅遠低于響應(yīng)電流中各頻率的振幅時，求解到的電路參數(shù)的相對誤差較小，對阻性、容性電流的分解準確度的影響也?。划敯自肼曋信c響應(yīng)電流不相關(guān)的頻率的振幅接近或高于響應(yīng)電流中某頻率的振幅時，部分電路參數(shù)的求解結(jié)果的相對誤差較大，甚至遠遠偏離真實值，以至于不能準確地實現(xiàn)阻性、容性電流的分解。

在響應(yīng)電流i(t)不含噪聲的前提下，當激勵電壓u(t)中含有諧波分量(主要是奇次諧波分量)時，即

(15)

式中：U2k-1為第2k-1次諧波分量的有效值；φ2k-1為第2k-1次諧波分量的初相位；ω為角頻率；n1=1，2，…。

對應(yīng)的多元線性方程為

y+Δy=BT×(x+Δx)。

(16)

式中：y為激勵電壓的基波單獨作用時對應(yīng)的響應(yīng)；Δy為激勵電壓的各次諧波協(xié)同作用時對應(yīng)的響應(yīng)；x=(x0，x1，…，x2n+1)T為激勵電壓的基波單獨作用時對應(yīng)的激勵；Δx=(Δx0，Δx1，…，Δx2n+1)T為激勵電壓的各次諧波協(xié)同作用時對應(yīng)的激勵。

由公式(16)可知，當u(t)含有諧波分量時，響應(yīng)y上會疊加Δy而影響到響應(yīng)電流i(t)的波形。但是，對于給定的非線性絕緣電介質(zhì)，其對應(yīng)的電路參數(shù)向量B是確定的，即為常系數(shù)向量。因此，通過多元線性方程(16)和(11)對應(yīng)的正則方程組求得的電路參數(shù)相同，即u(t)中有無諧波分量不影響阻性、容性電流的準確分解。

2 分解算法的實現(xiàn)

利用Matlab軟件編寫的程序來實現(xiàn)電路參數(shù)的求解和阻性電流與容性電流的分解，圖2給出了分解算法實現(xiàn)的程序流程圖。

圖2 分解算法流程圖Fig.2 Flow chart of decomposition algorithm

首先，輸入相關(guān)的數(shù)據(jù)和參數(shù)，根據(jù)相關(guān)公式得到矩陣X和向量Y。然后，用QR分解法求解正規(guī)方程組(12)，得到電路參數(shù)向量B，并用F檢驗對電路參數(shù)bk的顯著性進行檢驗，將作用不顯著的電路參數(shù)舍去。最后，將作用顯著的各電路參數(shù)代入式(2)和式(4)就可以得到阻性電流和容性電流。

3 仿真驗證

仿真所用的仿真參數(shù)為：G0x=4×10-8S，G2x=6×10-15S/V2，G4x=3×10-19S/V4，C0x=1×10-9F，C2x=2×10-17F/V2，C4x=5×10-22F/V4。

3.1 原理性驗證

采用的激勵電壓

u(t)=Umsin(2πft+φ)。

式中：Um=500 V；f=50 Hz；φ=0°。

根據(jù)仿真生成的激勵電壓u(t)和響應(yīng)電流i(t)，分別在n=2，α=0.05和n=3，α=0.05時用本文算法對響應(yīng)電流i(t)(如圖3所示)進行分解。表1給出了電路參數(shù)的計算結(jié)果。

圖3 響應(yīng)電流波形Fig.3 Simulation waveform of responding current

電路參數(shù)n=2n=3bk顯著性bk顯著性G0x4.0×10-8*4.0×10-8*G2x6.0×10-15*6.0×10-15*G4x3.0×10-19*3.0×10-19*G6x5.2×10-37C0x1.0×10-9*1.0×10-9*C2x2.0×10-17*2.0×10-17*C4x5.0×10-22*5.0×10-22*C6x4.3×10-38

由表1可知，n=2和n=3時分解得到的電路參數(shù)G0x～G4x和C0x～C4x與仿真生成響應(yīng)電流i(t)時所用到的相應(yīng)電路參數(shù)相同，而且通過顯著性檢驗可以將n=3時分解得到的作用不顯著的電路參數(shù)G6x和C6x舍去。仿真結(jié)果說明：在激勵電壓u(t)不含諧波分量和響應(yīng)電流i(t)不含噪聲的情況下，本文算法能夠?qū)Ψ蔷€性并聯(lián)等效電路的電路參數(shù)進行準確地求解，這為阻性、容性電流的準確分解提供了基礎(chǔ)。

將n=2時分解得到的作用顯著的各電路參數(shù)代入到式(2)和式(4)，得到阻性電流和容性電流波形，如圖4所示。

3.2 響應(yīng)電流含有噪聲

在仿真生成的響應(yīng)電流i(t)上分別疊加了信噪比SNR為55 dB和35 dB的白噪聲，且在n=2，α=0.05時用本文算法分別對這兩種情況下的i(t)進行分解，所得到的電路參數(shù)如表2所示。

圖4 分解得到的阻性電流和容性電流波形Fig.4 Resistive current and capacitive current waveforms gotten by decomposing

電路參數(shù)55dB35dBbk相對誤差/%bk相對誤差/%G0x3.996×10-8-0.093.964×10-8-0.90G2x6.292×10-154.873.427×10-15-42.89G4x2.997×10-19-0.093.231×10-197.69C0x1.000×10-90.011.001×10-90.05C2x1.841×10-17-7.931.704×10-17-14.83C4x4.686×10-22-6.285.757×10-2215.14

由表2可以得出，當噪聲的信噪比為55 dB時，所得到的電路參數(shù)的相對誤差較?。欢斣肼暤男旁氡葹?5 dB時，電路參數(shù)的相對誤差較大，特別地，G2x的相對誤差已經(jīng)達到-42.89%。結(jié)果表明：當噪聲有較大的信噪比時，比如55 dB，該算法能較準確地實現(xiàn)阻性、容性電流的分解。

3.3 激勵電壓含有諧波分量

考慮到實際電網(wǎng)中諧波分量主要為奇次，且7次以上的高次諧波含量占總諧波比例較小，這里僅對激勵電壓u(t)含有不同比例的諧波分量(7次及以下奇次諧波成分)的情況進行討論。

此時，

其中，U2k-1為第2k-1次諧波分量的有效值；ω=100π(rad/s)；k=1，2，3，4。

在n=2，α=0.05時，用本文算法對激勵電壓含有諧波分量時產(chǎn)生的響應(yīng)電流進行分解。表3列出了相關(guān)參數(shù)的計算結(jié)果。

表3 激勵電壓分別含不同比例的奇次諧波時相關(guān)參數(shù)的計算結(jié)果

通過表3可以看出，電路參數(shù)的相對誤差幾乎均為0%，這說明含有諧波分量的激勵電壓u(t)對電路參數(shù)的準確求解幾乎沒有影響。因此，在非標準正弦交流激勵電壓下，該算法也能夠準確地實現(xiàn)阻性和容性電流的分解。

4 實驗驗證

為了清晰說明本文算法的合理性與實用價值，將市售的MOA閥片切片成厚度為5 mm，直徑為42 mm的試樣。首先，對該試樣的直流伏安特性進行了測試，得到其電流密度J與電場強度E的關(guān)系曲線，如圖5所示。

圖5 電流密度與電場強度關(guān)系曲線Fig.5 Curve of current density against electric field strength

由圖5可知，該試樣的J-E特性曲線具有典型的非線性特性，因此采用該試樣作為本文實驗驗證的測試對象。

用文獻[15]所建立的非線性絕緣電介質(zhì)交流測試系統(tǒng)，對該試樣的交流激勵電壓u(t)和響應(yīng)電流i(t)的時域波形進行了測試，測試結(jié)果見圖6。測試條件為：室溫；標準正弦波激勵電壓，幅值為1 300 V，頻率為50 Hz；平板三電極系統(tǒng)。

圖6 激勵電壓和響應(yīng)電流的實測波形Fig.6 Test waveforms of excitation voltage and responding current

用本文算法在n≥3的情況下對此實測的響應(yīng)電流i(t)進行分解時，由于矩陣XTX近似于奇異矩陣，所得到的分解結(jié)果會不準確。為了得到準確的分解結(jié)果，用本文算法在n=2，α=0.05時對實測的響應(yīng)電流i(t)進行了分解，得到的電路參數(shù)見表4。圖7給出了響應(yīng)電流的擬合曲線和測試曲線對比結(jié)果。

表4 電路參數(shù)的計算結(jié)果

圖7 響應(yīng)電流回歸曲線和測試曲線Fig.7 Regression and measurement waveforms of responding current

由表4可以發(fā)現(xiàn)，這些電路參數(shù)的作用均顯著。通過圖7中響應(yīng)電流的擬合曲線和測試曲線的對比，可以得知該算法能夠?qū)崿F(xiàn)MOA閥片在交流電壓作用下的響應(yīng)電流的分解。利用該算法能夠同時得到MOA閥片的阻性和容性電流波形，如圖8所示。結(jié)合試樣的幾何尺寸，可以得到圖9所示的阻性電流密度與電場強度以及容性電流密度與電場強度導(dǎo)數(shù)的關(guān)系曲線。

圖8 阻性和容性電流波形Fig.8 Resistive and capacitive current waveforms

由圖5和圖9(a)的對比可知，在相同電場強度下，交流電壓下的阻性電流密度要比直流電壓下的大兩個數(shù)量級，這是由于MOA閥片中晶界處的界面極化[24]使得交流電導(dǎo)率較大的緣故。通過圖9(a)還可以發(fā)現(xiàn)，當電場強度足夠高時阻性電流密度jG會表現(xiàn)出非線性。

由圖9(b)可知，在交流電壓作用下，當電場強度導(dǎo)數(shù)足夠大時容性電流密度jC也會表現(xiàn)出非線性。因此，可以用本文算法得到的阻性、容性電流密度來研究MOA閥片在交流電壓作用下的非線性特性。

圖9 阻性電流密度與電場強度以及容性電流密度與電場強度導(dǎo)數(shù)關(guān)系曲線Fig.9 Curves of resistive current density vs electrical field strength and capacitive current density vs electrical field strength derivative

5 結(jié) 論

通過對非線性絕緣電介質(zhì)并聯(lián)等效電路的分析，以并聯(lián)等效電路參數(shù)為基礎(chǔ)推導(dǎo)了響應(yīng)電流關(guān)于激勵電壓的關(guān)系式，提出了一種基于多元線性回歸的阻性電流和容性電流分解方法。理論分析表明：該方法的原理簡單，易于實現(xiàn)，得到的阻性電流和容性電流中包含豐富的頻域信息，具有時域和頻域相結(jié)合的特點。仿真結(jié)果表明：該方法不僅可以準確的實現(xiàn)阻性、容性電流的分解，還可以得到非線性電路參數(shù)，其受較大信噪比噪聲的影響較小但幾乎不受激勵電壓諧波分量的影響。實驗結(jié)果表明：該方法可以實現(xiàn)MOA閥片在交流電壓作用下的響應(yīng)電流分解，同時發(fā)現(xiàn)其阻性、容性電流密度會表現(xiàn)出非線性特性，該方法具有潛在的實用價值。

[1] 譚向宇，許學(xué)勤，孫福，等．新型自適應(yīng)廣義形態(tài)濾波在 MOA 在線監(jiān)測數(shù)據(jù)處理中的應(yīng)用[J]．中國電機工程學(xué)報，2008，28(19):25-29．

TAN Xiangyu, XU Xueqin, SUN Fu, et al. Application of novel self-adaptive generalized morphological filter to MOA on-line monitoring data processing[J]. Proceedings of the CSEE, 2008, 28(19):25-29.

[2] 段大鵬，江秀臣，孫才新，等．基于正交分解的MOA 漏電有功分量提取方法[J]．電工技術(shù)學(xué)報，2008，23(7):56-61．

DUAN Dapeng, JIANG Xiuchen, SUN Caixin, et al. Extracting algorithm of active power component of MOA leakage current based on orthogonal decomposition[J]. Transactions of China ElectroTechnical Society, 2008, 23(7):56-61.

[3] ZHAO Tong, LI Qingmin, QIAN Jiali. Investigation on digital algorithm for online monitoring and diagnostics of metal oxide surge arrester based on an accurate Mode[J]. IEEE Transactions on Power Delivery, 2005, 20(2):51-756.

[4] LUNDQUIST O, STENSTROM L, SCHEI A, et al. New method for measurement of the resistive current of metal oxide surge arresters in service[J]. IEEE Transactions on Power Delivery, 1990, 5(4):1811-1822.

[5] 唐炬，張新強，孫才新，等．金屬氧化物避雷器阻性電流檢測的多元補償法研究[J]．高電壓技術(shù)，1999，25(1):36-41．

TANG Ju, ZHANG Xinqiang, SUN Caixin, et al. Study of multi-coefficient compensation method on resistive current passing through MOA[J]. High Voltage Engineering, 1999, 25(1):36-41.

[6] 臧殿紅．基于DSP 的氧化鋅避雷器在線監(jiān)測裝置的研究[D]．青島:中國石油大學(xué)，2009:12-13．

[7] 祝銘，鄒彬，萬善良，等．采用測量阻性電流提高MOA狀態(tài)監(jiān)測的可行性研究[J]．電瓷避雷器，2007，(6):31-34．

ZHU Ming, ZOU Bin, WAN Shanliang, et al. Study on feasibility of monitoring MOA with measurement of resistive current[J]. Insulator and Surge Arresters, 2007, (6):31-34.

[8] 朱寶森，宋璠，魏新勞．基于電壓過零點電流幅值分析的MOA阻性電流提取方法[J]．中國電機工程學(xué)報， 2012，32(16):165-171．

ZHU Baosen, SONG Fan, WEI Xinlao. An extracting method based on analysis of current magnitude at zero-point of voltage for resistive current of MOA[J]. Proceedings of the CSEE, 2013, 32(16):165-171.

[9] QI X, ZHENG Z, BOGGS S. Engineering with nonlinear dielectrics[J]. Electrical Insulation Magazine, 2004, 20(6):27-34.

[10] AUCKLAND D W,SU W,VARLOW B R.Nonlinear fillers in electrical insulation[J]. IEE Pro-Sci Meas Technol, 1997, 144(3):127-133.

[11] POKRYVAILO A,CARP C,SCAPELLATI C.Comparative testing of simple terminations of high-voltage cables[J]. Electrical Insulation Magazine, 2010, 26(1):7-14.

[12] DONZEL L,GREUTER F,CHRISTEN T.Nonlinear resistive electric field grading Part 2: materials and applications[J]. Electrical Insulation Magazine, 2011,27(2):18-29.

[13] BOGGS S A,DAMMON D H.Effect of insulation properties on the field grading of solid dielectric DC cable[J]. IEEE Transactions on Power Delivery, 2011, 16(4):456-461.

[14] MONGA S,GORUR R S,HANSEN P,et al.Design optimization of high voltage bushing using electric field computations[J]. IEEE Transactions on Dielectric and Electrical Insulation, 2006, 13(6):1217-1224.

[15] 鄭歡，李忠華，黃志鵬，等．非線性絕緣電介質(zhì)AC介電特性參數(shù)及其測量[J]．中國電機工程學(xué)報，2013，33(4):201-208．

ZHENG Huan, LI Zhonghua, HUANG Zhipeng, et al. AC dielectric characterization parameters of nonlinear insulating dielectrics and their measurement[J]. Proceedings of the CSEE, 2013, 33(4):201-208.

[16] 韓寶忠，高思明，李春陽，等．納米石墨/聚乙烯復(fù)合材料電導(dǎo)特性模擬仿真[J]．哈爾濱理工大學(xué)學(xué)報，2014，19(1):90-94．

HAN Baozhong, GAO Siming, LI Chunyang, et al. Simulation on electrical conductivity of nano graphite/PE composite[J]. Journal of Harbin University of Science and Technology, 2014, 19(1):90-94.

[17] 王飛風(fēng)，張沛紅，高銘澤．納米碳化硅/硅橡膠復(fù)合物介電性研究[J]．哈爾濱理工大學(xué)學(xué)報，2015，20(3):82-85．

WANG Feifeng, ZHANG Peihong, GAO Mingze. Research on the dielectric characteristics of nano-sic/silicone rubber composites[J]. Journal of Harbin University of Science and Technology, 2015, 20(3):82-85.

[18] 范勇，謝艷紅，趙偉，等．PAA短鏈處理納米粒子改性PI薄膜的介電性能[J]．哈爾濱理工大學(xué)學(xué)報，2015，20(5):56-60．

FAN Yong, XIE Yanhong, ZHAO Wei, et al. Dielectric properties of polyimide films modified by nano-particles treated with short-chain PAA molecules[J]. Journal of Harbin University of Science and Technology, 2015, 20(5):56-60.

[19] OBRZUT J, KANO K. Impedance and nonlinear dielectric testing at high AC voltages using waveforms[J]. IEEE Transaction on Instrumentation and Measurement, 2005, 54(4):1570-1574.

[20] 許恒迎，夏榮，李忠華．基于非線性最小二乘擬合的容性、阻性電流信號分解算法[J]．哈爾濱理工大學(xué)學(xué)報，2006，11(1):120-123．

XU Hengying, XIA Rong, LI Zhonghua. The capacitive current and resistive current decomposition algorithm using non linear least squares method[J]. Journal of Harbin University of Science and Technology, 2006, 11(1):120-123.

[21] TANAKA A., YAGI Y., TANAKA H. On-site diagnostic method for water treed XLPE cable by harmonics in AC loss current[J]. Electrical Insulation News in Asia, 2003(10):29-30.

[22] 梁馮珍，宋占杰，張玉環(huán)．應(yīng)用概率統(tǒng)計[M]．天津： 天津大學(xué)出版社，2004:363-369．

[23] Gerald R．數(shù)值方法和MATLAB實現(xiàn)與應(yīng)用[M]．北京:機械工業(yè)出版社，2004:337-344．

[24] 陳季丹，劉子玉．電介質(zhì)物理學(xué)[M]．北京:機械工業(yè)出版社，1982:178-183．

(編輯：劉素菊)

Decomposition method of resistive current and capacitive current based on multiple linear regression

HAN Yong-sen， LI Zhong-hua， ZHENG Huan， GUO Wen-min

(Key Laboratory of Engineering Dielectrics and Its Application, Ministry of Education, Harbin University of Science and Technology, Harbin 150080, China)

To study the insulation status and dielectric performance of nonlinear semiconductor equipment and nonlinear insulating dielectrics under AC voltage, a new decomposition method of resistive and capacitive current was proposed. Nonlinear parallel equivalent circuit which consists of nonlinear resistance and nonlinear capacitance was chosen as the research object. The nonlinear expression of responding current on excitation voltage was deduced and then was transformed to multiple linear expression with coordinate transformation. The equivalent circuit parameters were got and the decomposition of the resistive and capacitive currents was achieved using multiple linear regression. The anti-interference performance and the adaptability to the nonstandard sinusoidal voltage were analyzed qualitatively. The simulation results show that the proposed method can decompose the resistive and capacitive current accurately. The relative errors of circuit parameters are smaller when the signal to noise ratio of noise is 55 dB.The harmonics of excitation voltage have little effect on the solution of circuit parameters. The experimental results indicate that this method can decompose the total leakage current for MOA varistor under AC voltage.

current decomposition；nonlinear equivalent circuit；circuit parameters estimation；multiple linear regression

2013-08-07

國家重點基礎(chǔ)研究發(fā)展計劃項目(973計劃)(2014CB239504)

韓永森(1985—)，男，博士研究生，研究方向為非線性電介質(zhì)介電特性測試；

李忠華(1962—)，男，博士，教授，博士生導(dǎo)師，研究方向為非線性絕緣電介質(zhì)基礎(chǔ)理論和測試技術(shù)；

李忠華

10.15938/j.emc.2016.11.008

TM 85

A

1007-449X(2016)11-0053-08

鄭 歡(1981—)，男，博士，副教授，研究方向電氣絕緣監(jiān)測與診斷技術(shù)；

郭文敏(1977—)，女，博士，副教授，研究方向為非線性絕緣電介質(zhì)基礎(chǔ)理論。