李磊，關(guān)海鷗，杜松懷，祁廣云，李春蘭，蘇娟

（1.黑龍江八一農(nóng)墾大學(xué)信息技術(shù)學(xué)院，大慶163319；2.中國(guó)農(nóng)業(yè)大學(xué)信息與電氣工程學(xué)院；3.新疆農(nóng)業(yè)大學(xué)機(jī)械交通學(xué)院）

農(nóng)村低壓電網(wǎng)觸電生物電流頻譜特性分析

李磊1，關(guān)海鷗1，杜松懷2，祁廣云1，李春蘭3，蘇娟2

（1.黑龍江八一農(nóng)墾大學(xué)信息技術(shù)學(xué)院，大慶163319；2.中國(guó)農(nóng)業(yè)大學(xué)信息與電氣工程學(xué)院；3.新疆農(nóng)業(yè)大學(xué)機(jī)械交通學(xué)院）

針對(duì)農(nóng)村低壓電網(wǎng)發(fā)生觸電故障時(shí)觸電生物電流及其成分分量的復(fù)雜性和多樣性，采用快速傅里葉變換對(duì)觸電生物實(shí)驗(yàn)獲取的電源電壓、剩余電流以及觸電生物電流進(jìn)行頻譜分析，明確了剩余電流與觸電生物電流頻譜特性及各成分分量變換過(guò)程。仿真實(shí)驗(yàn)表明：在觸電故障發(fā)生時(shí)剩余電流與觸電生物電流相關(guān)系數(shù)為0.999，其中觸電生物電流中直流分量的最大畸變率最大可達(dá)0.044 5，2～5次諧波的最大畸變率分別為0.013 4，0.035 4，0.006 9，0.014 2，能夠?yàn)樘崛∩镉|電故障過(guò)程中剩余電流的變化規(guī)律提供依據(jù)。

剩余電流保護(hù)；快速傅里葉變換；觸電生物電流；頻譜分析

目前我國(guó)農(nóng)村用電水平發(fā)生了根本性的變化，并廣泛采用剩余電流保護(hù)裝置，農(nóng)村低壓電網(wǎng)的設(shè)施狀況得到了明顯的改善［1-2］。但農(nóng)村觸電傷亡事故仍然時(shí)有發(fā)生，在線運(yùn)行的剩余電流保護(hù)器［3］，普遍將剩余電流有效值的大小作為發(fā)生觸電故障的動(dòng)作依據(jù)［4-5］，其動(dòng)作電流的整定值與觸電生物電流無(wú)關(guān)。由于我國(guó)農(nóng)村地區(qū)實(shí)行的是大雜居小聚居，居民居住分散，運(yùn)行環(huán)境復(fù)雜，以及開(kāi)關(guān)切換和諧波干擾因素影響［6］，使得剩余電流保護(hù)器仍然存在誤脫扣現(xiàn)象，為了保證人身安全和低壓電網(wǎng)的正常運(yùn)行，需要開(kāi)發(fā)和研制基于人體觸電電流而動(dòng)作的剩余電流保護(hù)裝置。目前電力系統(tǒng)暫態(tài)故障的噪聲信號(hào)數(shù)字特征分析主要研究了衡量信號(hào)相似的度量包含重要信息［7-8］，文獻(xiàn)［9］利用高階統(tǒng)計(jì)量分析檢測(cè)信號(hào)的數(shù)字特征，提出了抽取總漏電電流變化過(guò)程表征指標(biāo)的定量描述方法。但由于生物體觸電故障發(fā)生過(guò)程的隨機(jī)性及復(fù)雜性，難以用一種時(shí)域方法完成準(zhǔn)確判斷。因此，研究針對(duì)發(fā)生觸電故障時(shí)觸電生物電流及其頻譜復(fù)多樣特點(diǎn)，采用快速傅里葉變換對(duì)觸電生物實(shí)驗(yàn)獲取的電源電壓、剩余電流以及觸電生物電流的時(shí)域信號(hào)進(jìn)行頻譜分析，明確剩余電流及觸電生物電流各成分分量的變化特性及相關(guān)性，為研發(fā)出新一代基于人體觸電電流的剩余電流保護(hù)器提供進(jìn)一步依據(jù)。

1 觸電生物電流原始信號(hào)

選取實(shí)驗(yàn)數(shù)據(jù)是2012年5月以來(lái)利用觸電物理實(shí)驗(yàn)平臺(tái)，以四肢哺乳動(dòng)物為實(shí)驗(yàn)對(duì)象，通過(guò)故障錄波器獲取其觸電故障過(guò)程的電氣信號(hào)及數(shù)據(jù)波形，其中包括：電源電壓、觸電生物電流、剩余電流。研究中截取其發(fā)生觸電故障時(shí)刻前5個(gè)工頻周期和后10個(gè)工頻周期，共計(jì)15個(gè)采樣周期（3000個(gè)采樣點(diǎn)）數(shù)據(jù)，如圖1所示。

圖1 原始信號(hào)Fig.1 Original signal

圖1原始信號(hào)中（a）為電源電壓的原始信號(hào)波形圖，橫坐標(biāo)為采樣點(diǎn)數(shù)，縱坐標(biāo)為其幅值（V）；（b）為觸電生物電流的原始波形圖，橫坐標(biāo)為采樣點(diǎn)數(shù)，縱坐標(biāo)為其幅值（mA）；（c）為剩余電流的原始信號(hào)波形圖，橫坐標(biāo)為采樣點(diǎn)數(shù)，縱坐標(biāo)為其幅值（mA）。

2 基于FFT的觸電生物電流

2.1 FFT的基本原理

研究選取工頻信號(hào)采樣周期長(zhǎng)度為輸入序列點(diǎn)數(shù)N=200，滑動(dòng)窗長(zhǎng)度為1，提取其直流分量幅值和9次以內(nèi)諧波分量幅值，并利用各個(gè)分量幅值與原始信號(hào)做相關(guān)性分析。依據(jù)文獻(xiàn)［10］結(jié)合實(shí)際觸電生物電流信號(hào)，采用一種快速離散傅里葉變換算法進(jìn)行分析。

設(shè)N點(diǎn)序列的x（n）點(diǎn)DFT表示成：

依據(jù)文獻(xiàn)［10］設(shè)離散信號(hào)具有各次諧波的周期信號(hào)可表示為：

式中，c0為直流成分，cm表示m次諧波的幅值，mω為角頻率，φm為m次諧波角相位，N為最高諧波次數(shù)。

2.2 觸電生物電流信號(hào)分析

在觸電故障發(fā)生時(shí)刻前5個(gè)工頻周期，觸電生物電流的幅值為零，剩余電流幅值為電網(wǎng)穩(wěn)定時(shí)對(duì)地泄漏電流的幅值，因此進(jìn)行觸電生物電流與剩余電流在發(fā)生觸電故障后10個(gè)工頻周期內(nèi)原始信號(hào)相關(guān)性分析。

設(shè)兩個(gè)信號(hào)x（t）和y（t）的相關(guān)程度可以采用相關(guān)系數(shù)ρ來(lái)衡量，其定義［11］為：

其中D（·）表示信號(hào)方差，E（·）表示數(shù)學(xué)期望。當(dāng)ρ=0時(shí)，兩信號(hào)互不相關(guān)；當(dāng)0＜|ρ|＜1時(shí)，兩信號(hào)相關(guān)；當(dāng)|ρ|=1時(shí)，兩信號(hào)相干，即完全相關(guān)。

電源電壓是決定剩余電流和觸電生物電流的根本條件。剩余電流保護(hù)技術(shù)中普遍以剩余電流為判斷依據(jù)，與觸電支路產(chǎn)生的觸電生物電流無(wú)關(guān)，如何明確剩余電流與觸電生物電流的關(guān)系成為主要問(wèn)題。所以研究主要針對(duì)剩余電流與觸電生物電流的相關(guān)性進(jìn)行分析，利用公式（3）觸電生物電流與剩余電流在觸電故障發(fā)生后的原始信號(hào)進(jìn)行相關(guān)性分析其相關(guān)性系數(shù)范圍為0.997 4～0.999，觸電生物電流和剩余電流在觸電后的原始信號(hào)具有高度相關(guān)性。

2.3 觸電生物電流頻譜特性分析

由于觸電前5個(gè)工頻周期沒(méi)有發(fā)生觸電故障，所以選取觸電故障發(fā)生后10個(gè)工頻周期數(shù)據(jù)進(jìn)行分析。剩余電流各成分分量與基波的比值，即剩余電流直流分量、2～9次諧波的畸變率［12］如圖2所示。

圖2 剩余電流各成分分量畸變率Fig.2 Aberration rate of each component of the residual current

圖2中橫坐標(biāo)為采樣點(diǎn)數(shù)，縱坐標(biāo)為各次諧波的幅值與基波幅值的比值即畸變率。直流分量的畸變率最小值為0.005 4，最大值為0.034 1；2～9次諧波的畸變率幅值范圍分別為：0.002 5～0.010 2，0.0136～0.029 5，0.000 020 064～0.005 3，0.005 7～0.012 0，0.000 016 421～0.003 2，0.002 0～0.006 7，0.000 004 458 3～0.002 5，0.004 4～0.008 3。剩余電流各成分分量的畸變率隨諧波頻率的增加而逐漸減小，各次諧波的畸變率具有隨采樣點(diǎn)數(shù)的增加，其幅值比例逐漸減小，各次諧波的畸變率存在具有先增加到極值而后逐漸衰減至平穩(wěn)狀態(tài)。

同時(shí)選取在觸電故障發(fā)生后10個(gè)工頻周期的數(shù)據(jù)計(jì)算觸電生物電流各成分分量畸變率，如圖3所示。

圖3中橫坐標(biāo)為采樣點(diǎn)數(shù)，縱坐標(biāo)為畸變率。剩余電流各成分分量的畸變率隨諧波頻率的增加而逐漸減小，各次諧波的畸變率大都具有隨著時(shí)間的增加，其幅值比例逐漸減小。其中直流分量的畸變率最小值為0.007 3，最大值為0.044 5；2～9次諧波的畸變率的值域分別為：0.003 2～0.013 4，0.014 7～0.035 4，0.000006 000 4～0.006 9，0.006～0.014 2，0.000 005 691 4～0.004 4，0.002 5～0.008 3，0.000 003 185 5～0.003 3，0.003 0～0.008 2，并且各次諧波的畸變率普遍具有先增加到極值而后逐漸衰減至平穩(wěn)。

由于現(xiàn)有剩余保護(hù)裝置普遍以剩余電流的大小作為動(dòng)作依據(jù)，無(wú)法分離剩余電流中觸電生物電流信號(hào)，供電回路中只能檢測(cè)剩余電流。所以，研究剩余電流與觸電生物電流的各次諧波畸變率之間存在的規(guī)律特性，對(duì)于研發(fā)基于觸電生物電流的剩余電流保護(hù)裝置具有重大意義。將發(fā)生觸電故障時(shí)刻后10個(gè)工頻周期的剩余電流與觸電生物電流各成分分量的畸變率變化過(guò)程，如圖4所示。

圖3 觸電生物電流各成分分量畸變率Fig.3 Aberration rate of each component of the electric shock bioelectric current

圖4 剩余電流與觸電生物電流各成分分量畸變率比較圖Fig.4 Comparison of each component’s aberration rate of residual current with that of electric shock bioelectric current

圖4橫坐標(biāo)為采樣點(diǎn)數(shù)，縱坐標(biāo)為各次諧波畸變率幅值，其中紅色和藍(lán)色分別代表觸電生物電流與剩余電流的畸變率曲線。觸電生物電流在上說(shuō)明其畸變率的幅值變化較為明顯，剩余電流較觸電生物電流變化頻率緩慢。二者在4、6、8次諧波重合率異常顯著。其中直流分量，2～9次諧波畸變率的極大值點(diǎn)分別在：89，40，1，17，30，7，14，6，1附近。2～9次諧波分別在：1 016，1 127，826，830，159，1 095，186，582點(diǎn)附近整體趨勢(shì)趨于平穩(wěn)狀態(tài)。從剩余電流與觸電生物電流各成分分量畸變率比較，觸電生物電流與剩余電流各成分分量畸變率的變化規(guī)律趨于一致。

3 結(jié)果分析

通過(guò)對(duì)觸電生物電流與剩余電流在觸電故障發(fā)生后10個(gè)工頻周期的原始信號(hào)及快速傅里葉變換提取直流分量、2～9次諧波并計(jì)算諧波畸變率及相關(guān)性結(jié)果分析如下：

（1）觸電生物電流原始信號(hào)分析結(jié)果

在生物發(fā)生觸電故障時(shí)，電源電壓對(duì)觸電支路在觸電發(fā)生時(shí)刻起到了決定性的作用，剩余電流為低壓電網(wǎng)發(fā)生觸電故障時(shí)流過(guò)剩余電流保護(hù)裝置主回路電流瞬時(shí)值的矢量和。剩余電流與觸電生物電流在觸電發(fā)生時(shí)刻二者變化趨勢(shì)相同，達(dá)到穩(wěn)態(tài)的點(diǎn)基本一致，并結(jié)合觸電生物電流信號(hào)的分析可得，觸電發(fā)生時(shí)刻剩余電流與觸電生物電流二者具有高度相關(guān)性。所以剩余電流在一定程度可以反映工程上不可測(cè)量的觸電支路中觸電生物電流的特性及變化規(guī)律。

（2）觸電生物電流頻譜分析結(jié)果

觸電生物電流與剩余電流各成分分量變化規(guī)律趨勢(shì)，可以計(jì)算n次諧波的相關(guān)性系數(shù)跨度為r（n）為：

依據(jù)公式（4）可計(jì)算剩余電流與觸電生物電流在觸電故障發(fā)生時(shí)刻后10個(gè)工頻周期的各成分分量的相關(guān)性系數(shù)跨度幅值，如圖5所示。

圖5 觸電生物電流與剩余電流觸電后10個(gè)工頻周期各成分分量相關(guān)性跨度Fig.5 Correlation span of each component of electric shock bioelectric current and residual current after the electric shock 10 power frequency cycle

圖5中橫坐標(biāo)為諧波次數(shù)，縱坐標(biāo)為相關(guān)性系數(shù)跨度幅值。觸電生物電流與剩余電流在觸電故障發(fā)生時(shí)刻后10個(gè)工頻周期內(nèi)各成分分量的相關(guān)性系數(shù)跨度普遍趨近于零點(diǎn)，直流分量相關(guān)性跨度為極大值0.125 6，5次諧波相關(guān)性跨度為極大值0.009；2次諧波相關(guān)性跨度為極小值0.001 2。奇次諧波相關(guān)性跨度比偶次諧波相關(guān)性跨度大，說(shuō)明奇次諧波相關(guān)性有一定波動(dòng)；偶次諧波相對(duì)較小表明其偶次諧波比奇次諧波的高相關(guān)性變化平穩(wěn)。

通過(guò)以上分析可知，剩余電流與觸電生物電流在觸電故障發(fā)生時(shí)刻后10個(gè)工頻周期內(nèi)各成分分量具有高度相關(guān)性，其中觸電生物電流與剩余電流在觸電故障發(fā)生時(shí)刻后10個(gè)工頻周期各成分分量的相關(guān)性系數(shù)平均值為0.989 5，并且觸電生物電流與剩余電流各成分的畸變率變化規(guī)律趨于一致；剩余電流與觸電生物電流各成分的畸變率的極值點(diǎn)、趨于平穩(wěn)狀態(tài)的點(diǎn)相對(duì)統(tǒng)一，說(shuō)明二者的整體變化趨勢(shì)具有高度相關(guān)性。所以，剩余電流在一定程度可以反映觸電生物電流的頻譜特性，為新一代剩余電流保護(hù)技術(shù)的研究提供參考基礎(chǔ)，以解決農(nóng)村低壓電網(wǎng)運(yùn)行環(huán)境復(fù)雜、干擾較多、自然條件多變等因素影響下，剩余電流保護(hù)裝置存在著誤動(dòng)和拒動(dòng)現(xiàn)象。為全面有效地支持我國(guó)新一輪農(nóng)村電網(wǎng)改造升級(jí)工程實(shí)施，充分發(fā)揮科技生產(chǎn)力，并提供有效的科學(xué)技術(shù)支持作用。

［1］張學(xué)玲.淺析農(nóng)村低壓電網(wǎng)觸電事故防范對(duì)策及措施［J］.現(xiàn)代農(nóng)業(yè)，2009（10）：104-105.

［2］董赫，李維凱.農(nóng)用高壓電網(wǎng)圓柱形殷鋼電容分壓器邊緣效應(yīng)的研究［J］.黑龍江八一農(nóng)墾大學(xué)學(xué)報(bào)，2015，27（2）：92-96.

［3］杜松懷，張莜慧，王淑平.電力系統(tǒng)接地技術(shù)［M］.北京：中國(guó)電力出版社，2011.

［4］Liu Jianfeng，Sun Yongguo.Research on residual current protective device based on microcontroller of MSP430［J］. Advanced Materials Research，2013，662：838-841.

［5］李由，付衛(wèi)強(qiáng)，包涵，等.農(nóng)機(jī)車(chē)載電子設(shè)備用戶界面框架設(shè)計(jì)［J］.黑龍江八一農(nóng)墾大學(xué)學(xué)報(bào)，2014，26（6）：84-87.

［6］李春蘭，杜松懷，廖從研，等.剩余電流保護(hù)裝置及其性能綜述［J］.電網(wǎng)技術(shù)，2009，33（11）：180-183.

［7］王濤，山拜，達(dá)拉拜，等.基于高階統(tǒng)計(jì)量的信號(hào)特征描述［J］.通信技術(shù)，2011，44（2）：151-154.

［8］田金超，張興周，韓志學(xué).基于高階統(tǒng)計(jì)量信號(hào)檢測(cè)的研究［J］.應(yīng)用科技，2006，33（1）：24-26.

［9］關(guān)海鷗，郭元建，杜松懷，等.農(nóng)村低壓電網(wǎng)觸電電流數(shù)字特征描述與分析［J］.中國(guó)農(nóng)業(yè)大學(xué)學(xué)報(bào)，2014，19（3）：202-207.

［10］易立強(qiáng)，鄺繼順.一種基于FFT的實(shí)時(shí)諧波分析算法［J］.電力系統(tǒng)及自動(dòng)化學(xué)報(bào)，2007，19（2）：98-102.

［11］劉春靜，張錦中.多徑效應(yīng)下的信號(hào)相關(guān)性分析［J］.火控雷達(dá)技術(shù)，2011，40（3）：16-20.

［12］潘紹明，羅功坤，勞有蘭.基于快速傅里葉變換的電力諧波分析與程序?qū)崿F(xiàn)［J］.廣西工學(xué)院學(xué)報(bào)，2012，23（1）：10-13.

Analysis of Electric Shock Bioelectric Current Spectrum Characteristic in Rural Low-voltage Power Grid

Li Lei1，Guan Haiou1，Du Songhuai2，Qi Guangyun1，Li Chunlan3，Su Juan2
（1.College of Information and Technology，Heilongjiang Bayi Agricultural University，Daqing 163319；2.College of Information and Electrical Engineering，China Agricultural University；3.College of Mechanical and Traffic，Xinjiang Agricultural University）

In terms of the complexity and variety of electric shock bioelectric current and its composition component when electric shock fault occurring in rural low-voltage power grid，residual current and electric shock bioelectric current spectrum characteristics and transformation process of each component were clarified by using fast Fourier transform to analyze the spectrum of the power supply voltage，residual current and electric shock bioelectric current obtained from electric shock bio-experiment.Simulation results showed that when the electric shock fault occurred，the residual current and electric shock bioelectric current correlation coefficient was 0.999，the greatest component in the electric shock bioelectric current aberration rate was up to 0.044 5，2-5 times of the highest harmonic aberration rate was 0.013 4，0.035 4，0.006 9 and 0.014 2 respectively，which could provide the basis for the extraction in the process of biological electric fault caused by variation of residual current.

residual current protection；fast fourier transform；electric shock bioelectric current；spectrum analysis

TM77

A

1002-2090（2016）03-0129-05

10.3969/j.issn.1002-2090.2016.03.025

2015-04-28

新疆維吾爾自治區(qū)高等學(xué)校科研計(jì)劃項(xiàng)目（XJEDU2012I15）；黑龍江八一農(nóng)墾大學(xué)校博士啟動(dòng)基金項(xiàng)目（XDB2015-10）；國(guó)家自然科學(xué)基金項(xiàng)目（51177165）。

李磊（1989-），男，黑龍江八一農(nóng)墾大學(xué)信息技術(shù)學(xué)院2014級(jí)碩士研究生。

杜松懷，教授，博士研究生導(dǎo)師，E-mail：songhuaidu@cau.edu.cn。