王碩禾， 許繼勇， 許惠敏， 蔡清亮

(1．石家莊鐵道大學(xué) 電氣與電子工程學(xué)院，河北 石家莊 050043;

2．上海理工大學(xué) 光電信息與計(jì)算機(jī)工程學(xué)院，上海 200093;3．北京鐵路局 石家莊供電段，河北 石家莊 050000)

0 引言

我國(guó)交流電氣化鐵路的牽引負(fù)荷，全部采用單相工頻25 kV 交流制，電氣化鐵路客運(yùn)專線牽引供電制式為單相工頻25 kV( 直供) 或2 ×25 kV 交流制( AT) 。隨著電力電子技術(shù)的迅速發(fā)展，電力機(jī)車已經(jīng)實(shí)現(xiàn)了由直流傳動(dòng)向現(xiàn)代交流傳動(dòng)的轉(zhuǎn)變。目前國(guó)際上已有數(shù)千臺(tái)的交流傳動(dòng)電力機(jī)車和電動(dòng)車組投入了實(shí)際運(yùn)營(yíng);CRH 型電動(dòng)車組的投入運(yùn)行標(biāo)志著我國(guó)的交流傳動(dòng)技術(shù)已經(jīng)走向?qū)嵱没２捎媒涣鱾鲃?dòng)模式，大幅度提高了牽引網(wǎng)的功率因數(shù)，并顯著改善了牽引取流波形，降低了電氣化鐵路牽引取流對(duì)電力系統(tǒng)電能質(zhì)量的干擾。即便如此，電氣化鐵路是仍然是電力系統(tǒng)的主要負(fù)序和諧波干擾源之一，世界各國(guó)對(duì)于電氣化鐵路對(duì)公用電網(wǎng)的干擾問題一直十分關(guān)注。

目前，在電能質(zhì)量檢測(cè)產(chǎn)品領(lǐng)域，高端產(chǎn)品基本上被美國(guó)Fluke、日本日置和瑞士萊姆等國(guó)際知名公司占據(jù)，這些公司的產(chǎn)品性能高但價(jià)格昂貴［1］。而國(guó)產(chǎn)電能質(zhì)量檢測(cè)產(chǎn)品還存在一定問題，如檢測(cè)指標(biāo)不全面、質(zhì)量不高、使用不便等。鐵路牽引供電系統(tǒng)一般直接供電方式( 如圖1( a) ) 具有4 條饋電線;AT供電方式( 如圖1( b) ) 每條供電臂含有接觸網(wǎng)( T) 和正饋線( F) ，共計(jì)8 路電流饋線。在牽引供電系統(tǒng)中，需要同時(shí)檢測(cè)高壓側(cè)和牽引側(cè)電壓、電流參數(shù)，以便分析各電量參數(shù)之間的相關(guān)性特點(diǎn)?？v觀整個(gè)電能質(zhì)量檢測(cè)領(lǐng)域，大多數(shù)是針對(duì)三相四線制的電力系統(tǒng)，不能滿足牽引供電系統(tǒng)中直供、AT 等不同供電方式下多通道實(shí)時(shí)同步采集的要求，因此研制專用于牽引供電系統(tǒng)的電能質(zhì)量檢測(cè)與分析系統(tǒng)具有現(xiàn)實(shí)意義和經(jīng)濟(jì)意義。利用LabVIEW 在測(cè)控領(lǐng)域的先進(jìn)技術(shù)和強(qiáng)大的數(shù)字信號(hào)處理能力的特點(diǎn)，將其作為系統(tǒng)的軟件平臺(tái)，提高了系統(tǒng)的通用性和檢測(cè)的精確度。

1 系統(tǒng)硬件設(shè)計(jì)

圖1 供電方式

檢測(cè)電網(wǎng)電能質(zhì)量即檢測(cè)電網(wǎng)的電壓及電流需要達(dá)到一定標(biāo)準(zhǔn)。為了實(shí)現(xiàn)實(shí)時(shí)分析電網(wǎng)電能質(zhì)量的要求，系統(tǒng)采樣的速率及精度必須很高。作為牽引供電系統(tǒng)的研究，非常關(guān)心不同供電臂上同一時(shí)刻電能質(zhì)量參數(shù)的變化，以此來分析不同載運(yùn)條件下供電系統(tǒng)的特征;而且系統(tǒng)需測(cè)量負(fù)序電流分量，這就要求采樣時(shí)必須多路同步。需要分析的電能質(zhì)量指標(biāo)多，計(jì)算復(fù)雜，運(yùn)算速度要求較快，現(xiàn)有的DSP、ARM 很難滿足。而且，系統(tǒng)要求具有普遍適用性，即在采用不同供電方式的牽引變電所都能使用該系統(tǒng)。基于以上做出的分析，并從性價(jià)比及實(shí)用性方面考慮，系統(tǒng)選用如圖2 所示的總體設(shè)計(jì)方案。

牽引變電所電壓互感器高壓側(cè)的電壓等級(jí)為220 kV 或110 kV，牽引側(cè)為55 kV 或27．5 kV，經(jīng)PT 取樣為100 V 電壓，再經(jīng)DVDI 系列互感器變?yōu)? ～1 V，實(shí)現(xiàn)牽引變電所電壓信號(hào)采集，并實(shí)現(xiàn)高壓與采集分析電路的隔離;牽引變電所電流互感器將高壓側(cè)和牽引側(cè)的電流變?yōu)? A，采用電流鉗將電流信號(hào)轉(zhuǎn)為0 ～2 V 電壓信號(hào)以便采集。電流鉗選用FLUKE i200s 交流電流鉗，可將電流信號(hào)轉(zhuǎn)為電壓信號(hào)以便采集。信號(hào)測(cè)量頻率范圍為40 Hz ～10 kHz，安全級(jí)別為CAT III 600 V。

圖2 總體設(shè)計(jì)方案

將由調(diào)理、采集電路采集的牽引供電系統(tǒng)的高壓側(cè)三相電壓和電流、牽引側(cè)兩牽引臂電壓和最多8 路饋線電流共16 路信號(hào)，送入高速多路同步采集卡進(jìn)行采集。采集卡選用阿爾泰公司生產(chǎn)的PCI2008 數(shù)據(jù)采集卡，該卡是12 位16 路同步采樣的AD 卡，采用PCI2 ．2 總線標(biāo)準(zhǔn)，支持多卡同步采集，通過率為25K/通道，單卡總通過率可達(dá)400 K，多卡采集總通過率可達(dá)1 M以上，并實(shí)現(xiàn)連續(xù)存盤，滿足系統(tǒng)設(shè)計(jì)要求。

2 系統(tǒng)軟件設(shè)計(jì)

系統(tǒng)的軟件設(shè)計(jì)由LabVIEW 開發(fā)平臺(tái)編寫。LabVIEW( Laboratory Virtual Instrument Engineering Workbench) 是一種廣泛被工業(yè)界、學(xué)術(shù)界和研究實(shí)驗(yàn)室所接受的圖形化的編程語言的開發(fā)環(huán)境，且被視為一個(gè)標(biāo)準(zhǔn)的數(shù)據(jù)采集和儀器控制軟件，也是電能質(zhì)量檢測(cè)裝置中常用的軟件開發(fā)平臺(tái)［2-3］?；谔摂M儀器的電能質(zhì)量檢測(cè)系統(tǒng)分為檢測(cè)和分析兩大部分。

檢測(cè)部分的任務(wù)是將原始電壓電流信號(hào)采集并變換成數(shù)字信號(hào)，然后顯示出波形和頻譜，存儲(chǔ)并輸出數(shù)據(jù)文件，其包含的功能模塊有數(shù)據(jù)傳輸模塊、數(shù)據(jù)分析模塊、圖形顯示模塊、數(shù)據(jù)文件輸出模塊，軟件結(jié)構(gòu)流程圖如圖3 所示。

分析部分是對(duì)監(jiān)測(cè)部分輸出的數(shù)據(jù)文件進(jìn)行系統(tǒng)的深入分析，得出電能質(zhì)量各項(xiàng)指標(biāo)參數(shù)的統(tǒng)計(jì)結(jié)果、變化趨勢(shì)等，并根據(jù)需要打印出牽引供電系統(tǒng)的電能質(zhì)量分析報(bào)表，軟件結(jié)構(gòu)流程圖如圖4 所示。

2．1 軟件設(shè)計(jì)中頻率測(cè)量算法的實(shí)現(xiàn)

常用的頻率測(cè)量算法有周期法、DFT ( FFT) 類算法［4］、線性調(diào)頻變換( CZT)［5］和小波分析［6］等。周期法通過測(cè)量信號(hào)波形相繼過零點(diǎn)間的時(shí)間寬度來計(jì)算頻率，該方法物理概念清晰、易于實(shí)現(xiàn)，但精度低，受諧波、噪聲和非周期分量的影響，實(shí)時(shí)性不好。在理想的模型下，選擇合適的參數(shù)，DFT ( FFT) 類和線性調(diào)頻變換( CZT) 、小波分析算法能準(zhǔn)確的計(jì)算出頻率，但存在需分析數(shù)據(jù)長(zhǎng)度較大、增大計(jì)算量和時(shí)滯等缺點(diǎn)。三點(diǎn)法是一種建立在三角函數(shù)變換基礎(chǔ)上的數(shù)據(jù)擬合方法。在被測(cè)函數(shù)是正弦函數(shù)，且采樣信號(hào)為等間隔采樣的前提下( 含非整周期采樣) ，利用相鄰連續(xù)的3 個(gè)采樣數(shù)據(jù)點(diǎn)，構(gòu)造求解信號(hào)頻率的線性方程，從而擬合出方程的系數(shù)，最終求出頻率［7］。

圖3 電能質(zhì)量在線監(jiān)測(cè)部分軟件結(jié)構(gòu)圖

設(shè)電壓或電流信號(hào)為u( t) = Umsin( ωt+φ) ，系統(tǒng)的采樣頻率為fs，則有兩相鄰采樣點(diǎn)的相角差為

圖4 電能質(zhì)量監(jiān)測(cè)分析部分軟件結(jié)構(gòu)圖

令α = ωt +φ，則任意3 個(gè)相鄰的采樣點(diǎn)均可表示為

因?yàn)閡( k) + u( k +2) = Um［sin αk+ sin( αk+2θ) ］= 2Umsin( αk+ θ) cos θ = 2u( k +1) cos θ。令y( n) =u( k) + u( k +2) ，x( n) = 2u( k +1) ，可有

通過構(gòu)造線性方程(3) ，只要對(duì)采樣數(shù)據(jù)進(jìn)行上述處理，通過最小二乘法擬合就可以得到一個(gè)較為準(zhǔn)確的斜率，即cosθ，根據(jù)式(1) 可以求出正弦波u( t) 的頻率f。LabVIEW 程序如圖5 所示。

圖5 三點(diǎn)法測(cè)頻率程序

實(shí)際應(yīng)用中，特別是用該方法進(jìn)行測(cè)量電力系統(tǒng)電壓( 電流) 基波頻率時(shí)，由于電力系統(tǒng)中電網(wǎng)電壓( 電流) 含有大量諧波成分，不可能為一標(biāo)準(zhǔn)正弦波，應(yīng)用該算法時(shí)，需對(duì)采樣信號(hào)進(jìn)行濾波。

2．2 軟件設(shè)計(jì)中負(fù)序測(cè)量算法的實(shí)現(xiàn)

在三相電路中，對(duì)于任意一組不對(duì)稱的三相電壓( 或電流) ，可以分解為三組對(duì)稱的相量，若以A 相作為基準(zhǔn)相時(shí)，三相相量與其對(duì)稱分量之間的關(guān)系為( 以電流為例)

式中，a = ej120°= -

三相電流不平衡度

牽引供電系統(tǒng)中，牽引變壓器主要采用的四種接線方式:單相結(jié)線、三相V/v 結(jié)線、三相YN/d11 結(jié)線和斯科特結(jié)線，牽引變壓器負(fù)荷都不會(huì)在電力系統(tǒng)中產(chǎn)生零序電流［8］，僅存在正序和負(fù)序電流，三個(gè)電流向量之和為零，可以構(gòu)成一個(gè)封閉的三角形［9］，如圖6 所示。

假定˙IA、˙IB、˙IB初相角分別為0°，θ1，θ2，根據(jù)余弦定理有

圖6 封閉三角形

由式(4) 和式(6) 就可以算出不接地系統(tǒng)中負(fù)序電流I2、相角θ2和正序電流I1、相角θ1，而在中性點(diǎn)接地的系統(tǒng)中，可以根據(jù)式(4) 求出各相中的負(fù)序電流后，減去零序電流，再運(yùn)用上述公式，同樣也可以求出系統(tǒng)中負(fù)序電流I2、相角θ2和正序電流I1、相角θ1。利用LabVIEW 的Matlab Script 節(jié)點(diǎn)，調(diào)用Matlab 程序。

2．3 諧波分析算法的實(shí)現(xiàn)

目前，國(guó)際上諧波檢測(cè)的理論主要有離散傅里葉變換、小波分析理論［10-11］和神經(jīng)網(wǎng)絡(luò)理論［12-13］。文獻(xiàn)［10］利用小波變換來評(píng)價(jià)經(jīng)支持向量機(jī)( SVM) 訓(xùn)練后的數(shù)據(jù)，計(jì)算出各頻段相應(yīng)的分頻系數(shù)，然后利用支持向量機(jī)分類的方法對(duì)信號(hào)諧波進(jìn)行分類;文獻(xiàn)［11］采用不同小波濾波器對(duì)電力系統(tǒng)中不同頻段中的特定的擾動(dòng)進(jìn)行多重濾波，以檢測(cè)電力系統(tǒng)中的擾動(dòng);文獻(xiàn)［12］采用前饋型遞歸神經(jīng)網(wǎng)絡(luò)訓(xùn)練樣本，將各諧波含量不需經(jīng)過傅里葉變換就可以分解出來，以實(shí)現(xiàn)檢測(cè)諧波;文獻(xiàn)［13］將經(jīng)過模塊神經(jīng)網(wǎng)絡(luò)訓(xùn)練后的樣本直接做S 變換，從而將電能質(zhì)量的11 種擾動(dòng)檢測(cè)出來。綜合比較上述兩種理論，小波變換理論在諧波檢測(cè)中常被用作濾波器使用，相比數(shù)字濾波器，具有良好的自適應(yīng)性，而神經(jīng)網(wǎng)絡(luò)理論雖然同樣可以實(shí)現(xiàn)較高精度的電能質(zhì)量分析，但是需要對(duì)樣本進(jìn)行多次訓(xùn)練，且訓(xùn)練精度與選擇樣本有關(guān)，計(jì)算量較大，故在一般電能質(zhì)量檢測(cè)中很少應(yīng)用。為便于討論，僅以離散傅里葉變換作為諧波檢測(cè)的方案為例，改進(jìn)算法將在另文中詳細(xì)討論。離散傅里葉變換是當(dāng)今諧波與間諧波檢測(cè)中應(yīng)用最多的一種方法。依據(jù)IEC 標(biāo)準(zhǔn)( IEC 61000-4-30) ，諧波分析DFT 譜線的計(jì)算，要求頻譜分析間隔為5 Hz，采用S 級(jí)計(jì)算諧波共11 條譜線，間諧波9 條譜線［14］，LabVIEW 程序如圖7 所示。

圖7 S 級(jí)DFT 諧波分析LabVIEW 程序

2．4 故障錄波功能

牽引接觸網(wǎng)是一個(gè)置于大自然中龐大的供電裝置，不僅工作條件惡劣，而且結(jié)構(gòu)復(fù)雜。牽引供電系統(tǒng)絕大部分故障發(fā)生在接觸網(wǎng)系統(tǒng)中，說明接觸網(wǎng)系統(tǒng)是牽引供電系統(tǒng)的薄弱環(huán)節(jié)。以2003 年和2004年接觸網(wǎng)故障統(tǒng)計(jì)數(shù)據(jù)為例，2003 年接觸網(wǎng)系統(tǒng)發(fā)生故障340 件，累計(jì)停電31 148 min。故障平均停電時(shí)間為91．61 min，占牽引供電系統(tǒng)故障的91%:2004 年接觸網(wǎng)系統(tǒng)發(fā)生故障274 件，累計(jì)停電25 160 min，故障平均停電時(shí)間為91．82 min，占牽引供電系統(tǒng)故障的95%［15］。綜合分析這些故障的共同點(diǎn)，可以發(fā)現(xiàn)，在故障發(fā)生前，電網(wǎng)會(huì)出現(xiàn)電壓跌落或電流驟升的現(xiàn)象，依據(jù)這一共同特點(diǎn)，本系統(tǒng)對(duì)牽引網(wǎng)中出現(xiàn)的電壓跌落或電流上升過程不加區(qū)別的采集樣本信息，實(shí)現(xiàn)以下程序功能:設(shè)定一定閾值，當(dāng)大于或小于閾值后，程序采集從故障前10 個(gè)周波始設(shè)定錄波時(shí)間內(nèi)的故障數(shù)據(jù)，用以分析故障信號(hào)的頻率分布，找出故障發(fā)生前的特征信號(hào)。

2．5 參數(shù)的設(shè)置與數(shù)據(jù)統(tǒng)計(jì)功能

牽引變電所參數(shù)選擇、增加、修改、刪除功能。由于不同牽引變電所采用的供電方式，結(jié)線方式的不同，各變電所都有其不同之處，需要設(shè)置不同的參數(shù)，在滿足不同牽引變電所的需要的同時(shí)，也要保證電能質(zhì)量分析系統(tǒng)的正確運(yùn)行。牽引供電系統(tǒng)電能質(zhì)量統(tǒng)計(jì)與分析。電能質(zhì)量分析系統(tǒng)不僅要求能實(shí)時(shí)顯示各種電能質(zhì)量指標(biāo)、電力參數(shù)，還需要具有數(shù)據(jù)統(tǒng)計(jì)分析功能，如計(jì)算指定時(shí)間段內(nèi)的報(bào)表數(shù)據(jù)。系統(tǒng)設(shè)計(jì)中利用LabVIEW 與數(shù)據(jù)庫(kù)技術(shù)，進(jìn)行數(shù)據(jù)分析并形成報(bào)表［16］。

3 結(jié)束語

本牽引供電系統(tǒng)電能質(zhì)量檢測(cè)分析系統(tǒng)，充分利用了虛擬儀器技術(shù)在測(cè)控領(lǐng)域的優(yōu)勢(shì)，硬件上采用Fluke i200s 電流鉗、阿爾泰科技PCI2008 12 位16 路數(shù)據(jù)采集卡和便攜式工控機(jī)實(shí)現(xiàn)數(shù)據(jù)的采集，軟件上利用NI 公司的虛擬儀器軟件LabVIEW 編制了適合牽引供電系統(tǒng)特點(diǎn)的電能質(zhì)量檢測(cè)與分析系統(tǒng)程序的界面，具備信號(hào)采集、分析計(jì)算、存儲(chǔ)、統(tǒng)計(jì)分析、報(bào)表等功能模塊，實(shí)現(xiàn)了穩(wěn)態(tài)5 項(xiàng)電能質(zhì)量指標(biāo)的監(jiān)測(cè)與分析功能，并且初步具有暫態(tài)電能質(zhì)量事件捕獲、分析的功能。

目前，整個(gè)系統(tǒng)所具備的功能在傳統(tǒng)電能質(zhì)量分析功能的基礎(chǔ)上，研制出適合于牽引供電系統(tǒng)特點(diǎn)的電能質(zhì)量檢測(cè)裝置。利用LabVIEW 的信號(hào)處理工具包，實(shí)現(xiàn)高級(jí)算法。在現(xiàn)有電能質(zhì)量分析平臺(tái)上，可以方便的實(shí)現(xiàn)對(duì)電力機(jī)車在車網(wǎng)耦合時(shí)，產(chǎn)生的過電壓而引起牽引變電所故障的原因進(jìn)行分析;實(shí)現(xiàn)對(duì)故障信號(hào)的分析，找出故障特征信號(hào)，從而為今后調(diào)整鐵路運(yùn)行方案提供依據(jù)，規(guī)避故障;實(shí)現(xiàn)不同車型負(fù)載對(duì)牽引網(wǎng)電能產(chǎn)生的影響的分析;電力機(jī)車型的識(shí)別;實(shí)現(xiàn)其他電能質(zhì)量檢測(cè)與分析，如電壓波動(dòng)與閃變等功能。系統(tǒng)經(jīng)監(jiān)督檢測(cè)后，先后應(yīng)用于石太客專、京廣線、石徳線等多個(gè)牽引變電所的測(cè)試中，使用結(jié)果顯示達(dá)到了設(shè)計(jì)的預(yù)期目的和要求。

［1］粟時(shí)平，劉桂英．現(xiàn)代電能質(zhì)量檢測(cè)技術(shù)［M］．北京:中國(guó)電力出版社，2008．

［2］雷振山，趙晨光．LabVIEW 8．2 基礎(chǔ)教程［M］．北京:中國(guó)鐵道出版社，2008．

［3］Julio Barros，Matilde de Apraiz，Ramón I． Diego． A virtual measurement instrument for electrical power quality analysis using wavelets［J］．Measurement，2009，42:298-307．

［4］李紹銘，紀(jì)萍，彭玉龍．基于加窗插值FFT 和動(dòng)態(tài)頻率的諧波檢測(cè)算法［J］．自動(dòng)化與儀器儀表，2009(3) :93-95．

［5］和昆英，郭虹，劉洛琨，等．一種FFT 和CZT 聯(lián)合的快速高精度頻率估計(jì)算法［J］．?dāng)?shù)字電視與數(shù)字視頻，2006( 8) :18-20．

［6］趙成勇，胥國(guó)毅，何明鋒，等．基于改進(jìn)遞歸小波的電力系統(tǒng)頻率測(cè)量［J］．電工技術(shù)學(xué)報(bào)，2006，20(6) :62-65．

［7］陳昊，曹福凱．基于LabVIEW 的三點(diǎn)法頻率測(cè)量技術(shù)［J］．試驗(yàn)技術(shù)與試驗(yàn)機(jī)，2007，11(4) :75-77．

［8］譚秀炳．交流電氣化鐵道牽引供電系統(tǒng)［M］．2 版．成都:西南交通大學(xué)出版社，2007．

［9］張文生．基于MATLAB 語言的負(fù)序電流通用計(jì)算法［J］．東北電力技術(shù)，2003(11) :37-38．

［10］S Tuntisak，S Premrudeepreechacharn．Harmonic Detection in Distribution Systems Using Wavelet Transform and Support Vector Machine［C］//Power Electronics and Motion Control Conference，2006 5th International．Lausanne:Switzerland，2006:1540-1545．

［11］C Sharmeela，M R Mohan，G Uma ，et al．A Novel Detection and Classification Algorithm for Power Quality Disturbances using Wavelets［J］．American Journal of Applied Sciences，2006，3(10) :2049-2053．

［12］Fevzullah Temurtas，Rustu Gunturkun，Nejat Yumusak，et al． An Application of Elman＇s Recurrent Neural Networks to Harmonic Detection［J］．Neural Networks，2004，3029:1043-1052．

［13］C N Bhende，S Mishra ，B K Panigrahi．Detection and classification of power quality disturbances using S-transform and modular neural network［J］． Electric Power Systems Research，2008，78:122-128．

［14］美國(guó)福祿克公司．現(xiàn)代電能質(zhì)量測(cè)量技術(shù)［M］．華盛頓:美國(guó)福祿克公司，2008．

［15］林飛，羅君，李志鋒．2003、2004 年全國(guó)電力牽引供電系統(tǒng)故障統(tǒng)計(jì)分析［J］．工業(yè)科技，2006，35(3) :27-28．

［16］楊樂平，李海濤，趙勇，等．LabVIEW 高級(jí)程序設(shè)計(jì)［M］．北京:清華大學(xué)出版社，2003．