一種基于實際工況的數字化電能表校驗方法及其誤差分析

寇英剛, 范 潔， 楊世海，陳 剛， 胡 琛， 穆小星， 徐敏銳

(1. 國網江蘇省電力公司, 江蘇 南京 210024；2. 國家電網公司電能計量重點實驗室(國網江蘇省電力公司電力科學研究院)，江蘇 南京 211103；3. 華中科技大學電氣與電子工程學院，湖北 武漢 430074)

一種基于實際工況的數字化電能表校驗方法及其誤差分析

寇英剛1, 范 潔1， 楊世海2，陳 剛2， 胡 琛3， 穆小星2， 徐敏銳2

(1. 國網江蘇省電力公司, 江蘇 南京 210024；2. 國家電網公司電能計量重點實驗室(國網江蘇省電力公司電力科學研究院)，江蘇 南京 211103；3. 華中科技大學電氣與電子工程學院，湖北 武漢 430074)

數字化電能表在實際工況下受到頻率波動、諧波、輸入噪聲等因素影響，常出現誤差超差的現象。為了研究現場實際工況下數字化電能表的計量性能，提出了一種基于實際工況的校驗方法，研制了工況復現裝置，推出了Blackman離散傅里葉變換(DFT)+自適應線性(Adaline)神經網絡算法，實現了標準電能的計算，并將該數字化電能表校驗方法和瓦秒法進行了比較分析。誤差分析結果表明基于實際工況的數字化電能表校驗方法和瓦秒法均能用于校驗數字化電能表，但是前者的測試結果波動更小，更加穩定,且能夠為現場復雜工況下電能表性能評估提供參考。

數字化電能表；實際工況；校驗；傅里葉變換；自適應線性神經網絡

0 引言

智能變電站利用電磁式互感及模擬量輸入合并單元或電子式互感器作為信號傳感單元， IEC 61850標準作為通信協議，數字化電能表作為計量單元，實現了基于IEC 61850協議的數字計量系統[1-3]。數字化電能表作為數字電能計量系統的重要組成部分，可應用于數字化變電站的貿易結算，因此研究其計量性能具有重要意義。

理論上，數字化電能表的準確性和可靠性較高。但實際上，通過實驗室檢定的數字化電能表在現場運行時多次出現誤差超差現象，且數字化電能表的溯源方案暫無定論，因此其還不具備貿易結算方面的應用條件[4-6]。由于數字化電能表是工作在頻率波動、諧波、輸入噪聲等復雜工況下的高精度儀器，實驗室的模擬環境不能完全還原復雜的現場工況，因此實驗室測試得到的不超差數字化電能表，在現場工況下的準確度還有待考核，有必要對數字化電能表在現場工況下的誤差狀態進行評估。

目前廣泛采用的數字化電能表校驗方法均為虛負荷校驗法，包括標準數字表法、標準數字源法、標準模擬表法和瓦秒法[7-13]。標準表法將標準數字表和被校數字表的電能誤差進行比較；標準源法利用輸出協議數據和高頻標準電能脈沖對數字化電能表進行校驗；標準模擬表法利用信號調理單元、協議轉換單元和標準模擬表，計算數字化表誤差；瓦秒法采用計算被校數字化電能表輸出電能脈沖數并測量時間的原理，計算被校表誤差。然而，上述數字化電能表校驗方案尚未建立現場實際工況和被校數字化電能表誤差之間的聯系，數字化電能表在現場實際工況下的誤差特性研究還存在空白。

文中提出了一種基于實際工況的數字化電能表校驗方法，研制了工況復現裝置，并提出了Blackman離散傅里葉變換(DFT)+自適應線性(Adaline)神經網絡算法，實現了標準電能的計算。

1 校驗方法的原理

校驗方法的基本原理如圖1所示。

圖1 基于實際工況的校驗方案原理Fig.1 Calibration scheme for digital watt-hour meter using real load

工況讀取模塊接收合并單元輸出的數據，存儲為自定義格式數據包；工況復現模塊將自定義格式數據包轉換為標準格式數據包，向被校表發包，同時記錄當前已發送數據包的數量；被校表接收標準格式數據包后輸出電能脈沖，集成于校驗儀的標準表，同時計算標準電能數據，兩路數據進行比較，得到被校表的誤差。校驗儀一方面接收實際工況讀取模塊的數據文件，另一方面利用反饋的數據包值截取數據文件進行電能計算。

電能脈沖計算原理如圖2所示。在T0時刻收到第一個脈沖，記錄當前發包數量N0；在T1時刻收到第二個脈沖，記錄當前發包裝置N1；在Ti時刻收到第i+1個脈沖，記錄當前發包裝置Ni。利用數據文件中N0～Ni的原始數據進行計算，實現對數據文件的精確截取，標準通道和被校表用同一組數據進行電能計算，得出二者計算電能的誤差，測試時間設為5 min，數字化電能表的相對誤差為：

(1)

式中：m0為標準數字化電能表算定的脈沖數；m為被校數字化電能表脈沖數。

圖2 實際工況下電能脈沖計算原理Fig.2 Calculation principle of energy pulse under real load

2 校驗方法的實現

2.1 工況讀取模塊

工況讀取模塊如圖3所示，包括光纖收發器和計算機(PC機)。首先由光纖收發器(TP-LINK TR-932D)將合并單元光口發出的數據轉化為電信號；其次與PC機的以太網口進行通信；最后PC機進行IEC 61850-9-2 LE協議數據幀的抓包、解析和提取工作，以自定義格式文件形式存儲。為保證抓包不丟包和數據解析的完整性，工況讀取模塊軟件基于WinPcap和C語言開發。模塊工作流程為：初始化網卡并進行配置，建立合并單元和工況讀取模塊之間的通信；基于pcap_next_ex函數接收報文，讀取報文中的Ethertype字段，判斷報文類型；當報文類型為88BA時，讀取報文中應用協議數據單元(APDU)內容，并對其中的計量用通道進行識別；最后將計量用通道的數據保存為自定義格式數據。

圖3 工況讀取模塊Fig.3 Real load reading module

2.2 工況復現模塊

工況復現模塊如圖4所示。自定義格式數據在該模塊中添加目的地址、源地址、采樣值標識ID(SVID)等IEC 61850-9-2 LE協議相關參數，轉換為標準格式后等間隔發包(4 kHz)。

圖4 工況復現模塊Fig.4 Real load reconstruction module

該模塊設計基于現場可編程門陣列(FPGA)板和PC機。FPGA板為Altera Cyclone Ⅱ EP4CE15F17C8N系列。FPGA板外部擴展一個電可編程只讀存儲器(EPROM)，FPGA板和PC機之間通過一個串口和一個100 M網口連接，實現雙向通信。在PC機組幀完畢后，百兆網口發送符合IEC 61850-9-2 LE協議的數據包，串口用于FPGA向PC機反饋存儲狀態和發包數量。當EPROM儲存數據包不足時，FPGA向PC機反饋儲存不足，PC機停止發送IEC 61850-9-1/2協議數據包進行儲存，否則繼續數據包的發送。

2.3 校驗儀

校驗儀如圖5所示。由PC機完成被校數字電能表電能和標準數字電能之間的計算與比較工作。工況讀取模塊、復現模塊、校驗儀共用一臺PC機。 PC機的CPU為Intel core i5-4460，存儲容量為500 GB。PC機除了工況讀取模塊和工況復現模塊的接口外，還擴展脈沖采集板，獲取被校數字化電能表的電能脈沖，并計算被校通道電能。當合并單元數據速率為8 Mbit/s時，數據存儲時間可達138 d。

圖5 校驗儀Fig. 5 Calibration instrument

3 標準數字電能計量算法

為實現標準電能的計算，文中提出Blackman DFT+ Adaline神經網絡算法。 Blackman DFT算法應用較廣[14]，因此文中主要對Adaline神經網絡進行分析。

在實際工況下電網中含有各次諧波，電網波形可以表示為周期信號：

(2)

式中:Ak為信號幅值;k為諧波次數;m為最高諧波次數;f0為電網基波頻率;φk為初相位。將式(2)離散化后得到表達式(3)：

(3)

式中：n在0～(N-1)間取值，N為采樣點數；d0為直流分量；ak=Aksinφk；bk=Akcosφk；ω0=2πf0。只要得到ak和bk的值，就可以求出各次諧波的幅值和相角。

Adaline神經元結構原理如圖6所示。其中x1n～xkn為Adaline神經元在n時刻的輸入信號，取其向量形式為Xin=[x0n,x1n,…,xkn]T，對應同一時刻每一個輸入信號都有一個相應的權值，取其向量形式為Win=[w0n,w1n,…,wkn]T。Adaline神經元的輸出為：

(4)

圖6 Adaline神經元結構原理Fig.6 Adaline neuron structure

Wi(n+1)=Win+ηenXin

(5)

式中：η稱為步長，也稱學習速率，0<η<1。

在進行諧波分析時取基波頻率f0為已知，即基波角頻率ω0=2πf0已知。初始值取：

(6)

校驗儀采集數字電能表輸出的脈沖時間間隔為t，采用時間-數字轉換(TDC)技術，將時間t與數字功率源輸出數據頻率f作乘積，取整得到時間間隔t內數字功率源輸出的點數N，將抓包程序解出的數據截取N個數據點，對N個數據點利用文中算法解析出數據波形中基波、諧波、頻率數據及N點的功率。

將文中算法與常規DFT算法和插值4階卷積窗DFT算法[15-17]進行比較，計算數字電能計量的誤差計算。電壓/電流有效值均取1，功率因素為1.0，基波頻率為50 Hz，對輸入噪聲、頻率偏差、諧波等因素進行誤差仿真，仿真次數為10，表1為仿真結果的平均值。

表1 算法誤差仿真結果Table 1 Simulation results of algorithm error

可知：(1) 輸入噪聲誤差仿真時，在基波上疊加了高斯白噪聲，信噪比為30 dB，誤差為隨機性誤差，仿真結果中最大誤差為10-4數量級，誤差幾乎為零，可以忽略；(2) 頻率偏差誤差仿真時，頻率值取50.5 Hz，頻率偏差為0.5 Hz，誤差為10-7數量級，誤差較小，可以忽略；(3) 諧波誤差仿真時，取2～11次諧波共同疊加至基波中，按照電能質量公用電網諧波國家標準中的規定對諧波允許值進行取值，取偶次諧波為基波的0.8%，奇次諧波為基波的1.6%，電壓總諧波畸變率為：

(7)

式中：h表示基波以上的各次諧波次數，電流以相同方法進行取值，仿真誤差為10-7數量級。

由結果可見：常規DFT算法在頻率偏差、諧波情況下誤差較大；插值DFT算法和文中算法的計算結果比較接近，但是信號疊加噪聲工況下的誤差較大；文中提出的數字電能計量算法在各種復雜實際工況下，均能準確的進行電能計量。

4 校驗方法的驗證

為驗證文中校驗方法的可行性，采用瓦秒法和文中方法對同一塊0.2 s級的數字電能表進行校驗。電能計算均采用點積和算法，取連續測試的10個脈沖作為測試結果，均值計算結果如表2所示。

表2 測試結果Table 2 The experimental results %

采用瓦秒法測試數字化電能表的誤差和采用基于實際工況的校驗方法測試誤差的均值基本相同。但是瓦秒法不確定度較大，而基于實際工況的校驗方法測試精度高，誤差波動小。試驗表明了該校驗方法檢定數字化電能表的有效性和可行性。

5 結語

文中分析了數字化電能表校驗方案，研究了基于實際工況的數字化電能表校驗方法，研制了工況復現裝置，提出了Blackman DFT+Adaline神經網絡算法計算標準電能，并和瓦秒法進行了比較。研究和試驗結果表明：(1) 基于實際工況的校驗方法可用于重建現場工況，評估復雜環境下數字化電能表的計量性能；(2) 與瓦秒法測試結果相比，該校驗方案測試結果的波動較小，穩定性較高。

[1] 艾 兵，江 波. 數字電能計量及其電能表檢測技術[J]. 四川電力技術，2011，34(2)：10-13，17.

AI Bing, JIANG Bo. Digital electric energy metering and its detection technique[J]. Sichuan Electric Power Technology，2011，34(2)：10-13，17.

[2] 王忠東，李紅斌，程含渺，等. 模擬量輸入合并單元計量性能測試研究[J]. 電網技術，2014，38(12): 3522-3527.

WANG Zhongdong, LI Hongbin, CHENG Hanmiao, et al. Metering performance test of analog-input merging unit[J]. Power System Technology，2014，38(12)：3522-3527.

[3] 程含渺，徐 晴，紀 峰，等. 一種用于電子式互感器相位補償的數字移相方法[J]. 電力工程技術，2017，36(2)：82-87.

CHENG Hanmiao, XU Qing, JI Feng, et al. A digital phase-shift method for phase compensation of electronic instrument transformer[J]. Electric Power Engineering Technology，2017，36(2)：82-87.

[4] 謝 東, 吳 濤, 吳伯華，等. 數字量輸入電能表誤差等級問題的探討[J]. 電測與儀表, 2011(8):87-91.

XIE Dong, WU Tao, WU Bohua, et al. Discussion on accuracy class issue for digital interface electricity meter[J]. Electrical Measurement & Instrumentation, 2011(8):87-91.

[5] 張秋雁，魏 偉，李紅斌. 基于IEC 61850標準的數字化電能表丟幀誤差分析[J]. 電測與儀表，2016，53(24)：56-60.

ZHANG Qiuyan, WEI Wei, LI Hongbin. Frame loss error analysis of digital Watt-hour meter based on IEC 61850[J]. Electrical Measurement & Instrumentation，2016，53(24)：56-60.

[6] 段梅梅，吳偉將，紀 峰，等. 頻率偏移對數字化電能表計量誤差影響分析[J]. 電測與儀表，2016，53(S1)：63-66.

DUAN Meimei, WU Weijiang, JI Feng, et al． Research on influence of frequency deviation on digital input electricity meters[J]. Electrical Measurement & Instrumentation. 2016，53( S1)：63-66.

[7] 李 前，章述漢，陸以彪，等. 數字電能計量系統現場檢定技術研究[J]. 電測與儀表，2010，47(10)：25-28.

LI Qian，ZHANG Shuhan，LU Yibiao，et al． The Study of on-site meteringand calibration of fully digital power metering system[J]. Electrical Measurement & Instrumentation. 2010, 47( 10) : 25-28.

[8] 楊華云，艾 兵，江 波，等. 數字電能表現場校驗儀的研制[J]. 電測與儀表，2013，50(6)：53-57.

YANG Huayun，AI Bing，JIANG Bo，et al． Development of the field tester for digital electricity meter[J]. Electrical Measurement & Instrumentation，2013，50(6)：53-57.

[9] 張 健，穆小星，段梅梅. 一種小模擬量輸入電能表檢定裝置的研制[J]. 江蘇電機工程，2015，34(2)：55-57.

ZHANG Jian, MU Xiaoxing, DUAN Meimei. Development of a small analog input testing equipment for electricity meter [J]. Jiangsu Electrical Engineering，2015，34(2)：55-57.

[10] 唐 毅，江 波，李紅斌，等. 數字電能計量系統檢定方法綜述[J]. 電工技術學報, 2013，28(S2)：372-377.

TANG Yi, JIANG Bo, LI Hongbin, et al. Survey on the verification methods for digital electrical energy metering[J]. Transactions of China Electrotechnical Society，2013，28(S2)：372-377.

[11] 林偉斌，陳 壘，肖 勇，等. 電能表現場校驗技術研發應用現狀及其發展趨勢(上)[J]. 電測與儀表, 2016, 53(12): 1-6,41.

LIN Weibin, CHEN Lei, XIAO Yong, et al. Current research status and trend of on-site calibration of Watt-hour meters[J]. Electrical Measurement & Instrumentation. 2016，53(12)：1-6,41.

[12] 林國營，周尚禮，孫衛明，等. 數字化變電站電能計量裝置檢驗技術[J]. 電力系統及其自動化學報，2011，23(3)：145-148.

LIN Guoying, ZHOU Shangli, SUN Weiming, et al. Calibration technology of non-traditional electric energy measuring equipment in digitalized substation[J]. Proceedings of the CSU-EPSA. 2011，23(3)：145-148.

[13] GB/T 17215.303—2013 交流電測量設備 特殊要求 第3部分：數字化電能表[S]. 2013.

GB/T 17215.303—2013 Electric metering equipment (AC) particular requirements dart 3: digital input electricity meters[S]. 2013.

[14] 許 珉，張文強，王兆強，等. 基于加窗遞推DFT算法的快速相位差校正法研究[J]. 電力系統保護與控制，2010，38(14)：1-4.

XU Min, ZHANG Wenqiang, WANG Zhaoqiang,et al. Fast phase difference correction method by using recursive DFT with blackman-harris window[J]. Power System Protection and Control，2010，38( 14) : 1-4.

[15] ZHANG F, GENG Z, YUAN W. The algorithm of interpolating windowed FFT for harmonic analysis of electric power system[J]. IEEE Transactions on Power Delivery, 2001, 16(2)：

160-164.

[16] QIAN H, ZHAO R, CHEN T. Interharmonics analysis based on interpolating windowed FFT algorithm[J]. IEEE Transactions on Power Delivery，2007，22(2)：1064-1069.

[17] 王 平，丁 力. 加窗插值FFT的電網諧波分析算法研究[J]. 江蘇電機工程，2008，27(2)：46-50.

WANG Ping, DING Li． Power system harmonic measurement analysis based on windows and interpolated FFT [J]． Jiangsu Electrical Engineering，2008，27(2)：46-50.

寇英剛

寇英剛(1974—)，男，天津人，碩士，高級工程師，從事電能計量管理工作(E-mail：Yinggang_kou@qq.com)；

范 潔(1977—)，女，江蘇南通人，碩士，教授級高級工程師，從事電能計量檢測與管理工作(E-mail：fanjie125@163.com)；

楊世海(1976—)，男，江蘇徐州人，碩士，教授級高級工程師，從事電子式互感器、數字化電能計量研究工作(E-mail：ysh.young@163.com)；

陳 剛(1984—)，男，湖北黃梅人，碩士，高級工程師，從事電子式互感器、數字化電能計量研究工作(E-mail：cg09_jseprc@sina.com)；

胡 琛(1986—)，男，江西南昌人，博士研究生，研究方向為電子式互感器狀態評估、數字化電能計量系統仿真平臺研發(E-mail: huchena@126.com)；

穆小星(1963—)，男，江蘇南京人，碩士，教授級高級工程師，從事電氣自動化測試與研發工作(E-mail：m139129@163.com)；

徐敏銳(1976—)，男，江蘇淮安人，碩士，教授級高級工程師，從事互感器現場檢測、數字化電能計量研究工作(E-mail：xuminrui18@163.com)。

(編輯方 晶)

ACalibrationMethodandErrorAnalysisofDigitalEnergyMeterBasedonActualWorkingCondition

KOU Yinggang1, FAN Jie1, YANG Shihai2,CHEN Gang2,,HU Chen3,MU Xiaoxin2，XU Mingrui2

(1. State Grid Jiangsu Electric Power Company, Nanjing 210024,China; 2. State Grid Key Laboratory ofEnergy Metering(State Power Research Jiangsu Electric Power Company Institute), Nanjing 211103,China;3. SEEE of Huazhong University of Science and Technology,Wuhan 430074, China)

Digital energy meters often appear error phenomenon in the actual working conditions such as frequency fluctuations, harmonics, input noise and other factors. In order to study the measurement performance of the digital energy meter under the actual working conditions, a calibration method based on the actual working condition is presented. A condition recovery device is developed. A Blackman discrete Fourier transform (DFT) + adaptive linear (Adaline) neural network algorithm is proposed. The calculation of standard electrical energy is achieved. The digital energy meter calibration method and the wattage-second method are compared and analyzed. The results of error analysis show that the digital energy meter calibration method based on the actual working conditions and the wattage-second method can be used to verify the digital energy meter. But the test results of the former are smaller and more stable, and can provide reference for the performance evaluation of the energy meter under complex conditions.

digital watt-hour meter; actual working condition; calibration; DFT; Adaline neural network

TM933

A

2096-3203(2017)06-0053-05

2017-06-27；

2017-07-26

國家重點研發計劃資助項目(2016YFB0901104);國家電網公司科技項目(5210EF17001M)