張志鑫 楊緒龍 徐賽梅 朱云武 任楠楠

?

Allan方差在互感器精度評測中的應用研究

張志鑫1楊緒龍2徐賽梅2朱云武2任楠楠1

（1. 北京航天時代光電科技有限公司，北京 100094； 2. 云南電網(wǎng)公司楚雄供電局，云南楚雄 675000）

按照與一次輸入信號的相關性，可將噪聲分為兩大類，即與一次輸入相關的乘性噪聲和與一次輸入無關的加性噪聲。本文在互感器精度評測中引入Allan方差分析法，充分發(fā)揮Allan方差在隨機誤差辨識方面的性能，實現(xiàn)了互感器乘性噪聲和加性噪聲的分離，數(shù)據(jù)仿真結果證明了方法的有效性。

互感器；乘性噪聲；加性噪聲；比差；Allan方差

互感器在傳變一次信號、信號傳輸、數(shù)據(jù)處理的過程中，不可避免的會引入噪聲干擾。其中，加性噪聲包含可建模補償?shù)牧阄灰约皼]有固定模型的隨機噪聲，隨機噪聲表現(xiàn)為測量值圍繞真實值的上下波動，其圍繞真實值的波動的范圍與一次輸入信號無關，可用隨機噪聲的方差描述；而乘性噪聲則表現(xiàn)為測量值偏離真實值的趨勢項，是影響互感器性能指標的關鍵，對互感器在電能計量和繼電保護中的應用帶來潛在的風險和隱患。

目前，互感器精度評測方法無法準確區(qū)分乘性噪聲和加性噪聲。本文中采用Allan方差分析法，應用在互感器精度評測中。

1 精度評測中的噪聲項

研究兩類噪聲信號，當一次側輸入信號有效值為in時，互感器的輸出有效值可表示為

式中，0為乘性噪聲信號；1為加性噪聲信號，互感器的零位已經(jīng)過標定扣除。

由標準GB/T 20840.8—2007中互感器比差的計算方法可得[1]

式中，A為比差測量值，（-1）為比差真實值；0為乘性噪聲引起的比差測量誤差，表現(xiàn)為比差測量值的趨勢項；1/in為加性噪聲引起的比差測量誤差，表現(xiàn)為比差測量值圍繞真實值的上下波動。

2 Allan方差法分離噪聲項

Allan方差法的主要特點是能夠比較容易地對各種誤差源及其對整個噪聲統(tǒng)計特性進行細致地表征和辨識，而且便于計算、易于分離。

Allan方差的定義如式（3）所示[2-4]。

對于數(shù)據(jù)的間隔時間為0的比差測量值數(shù)據(jù)A()（=1,2,…,）的Allan方差計算步驟如下[5-6]：

1）對數(shù)據(jù)進行二次采樣，采樣時間為=0，分別取等于0, 20, …,0（＜/2，求出每一時間長度為的數(shù)組平均值，即

對數(shù)組平均值求差，即

（5）

2）計算相關時間時的Allan方差，即

3）取不同的相關時間，重復步驟1至步驟3，可得()與的雙對數(shù)坐標曲線，該曲線稱為Allan方差曲線，如圖1所示。

圖1 Allan方差曲線

圖1中，量化噪聲、比差隨機游走、馬爾可夫噪聲、正弦噪聲與一次輸入信號無關，共同構成了精度評測中的加性噪聲部分；比差趨勢項隨機游走、比差趨勢項噪聲與一次輸入信號有關，共同構成了精度評測中的乘性噪聲部分。圖1中斜率為0的部分為噪聲項分離之后的比差不確定度。

在互感器精度評測中，采用Allan方差法，可方便地獲取互感器的加性噪聲方差、乘性噪聲方差以及比差真實值的方差（不確定度），為全面評測互感器的性能指標提供了數(shù)據(jù)基礎。

3 仿真分析

無噪聲信號=6sin(2p50)的波形如圖2所示。

圖2 無噪聲50Hz信號

在此基礎上構造了50Hz正弦信號與加性噪聲疊加信號、50Hz正弦信號與乘性噪聲疊加信號、50Hz正弦信號與乘性噪聲和加性噪聲疊加信號等仿真數(shù)據(jù)，并采用Allan方差法進行不同種類噪聲的分離。

3.1 50Hz正弦信號與加性噪聲疊加

50Hz正弦信號與均值為0正態(tài)分布的加性噪聲疊加如式（7）所示。時域波形如圖3所示。

圖3 時域波形

式中，加性噪聲1()的概率密度函數(shù)滿足

（8）

每10個周波（200ms）計算一次比差測量值，結果如圖4所示。

比差測量值的Allan方差如圖5所示。

圖5表明，比差測量值中只含有比差隨機游走等加性噪聲，不含趨勢項噪聲和趨勢項隨機游走等乘性噪聲。圖5中Allan方差模型可表示為

圖4 比差測量值

圖5 比差測量值的Allan方差

式中，當=-1時，由最小二乘法計算可得比差隨機游走，即

圖5中，加性噪聲Allan方差不超過4%（1），代表了圖4中比差測量值縱軸的波動范圍。隨著相關時間的延長，加性噪聲對比差測量值的影響逐漸減小，因此，可通過對比差測量值取平均的方法剔除加性噪聲的影響，圖5中，取50s比差測量值的平均值可使加性噪聲的影響不超過0.1%（1）。

3.2 50Hz正弦信號與乘性噪聲疊加

50Hz正弦信號與幅值線性變化的乘性噪聲疊加如式（10）所示。時域波形如圖6所示。

每10個周波（200ms）計算一次比差測量值，結果如圖7所示。

比差測量值的Allan方差如圖8所示。

圖8表明，比差測量值中不含量化噪聲、比差隨機游走等加性噪聲，只含趨勢項噪聲和趨勢項隨機游走等乘性噪聲。由于趨勢項的存在，導致比差測量值與真實值的偏離隨時間逐漸增加，如圖7所示。

圖6 時域波形

圖7 比差測量值

圖8 比差測量值的Allan方差

3.3 50Hz正弦信號與乘性噪聲和加性噪聲疊加

50Hz正弦信號與幅值線性變化的乘性噪聲和均值為0正態(tài)分布的加性噪聲疊加如式（11）所示。時域波形如圖9所示。

圖9 時域波形

每10個周波（200ms）計算一次比差測量值，結果如圖10所示。

圖10 比差測量值

比差測量值的Allan方差如圖11所示。

圖11 比差測量值的Allan方差

圖11中，比差測量值既含有加性噪聲又含有乘性噪聲。其中加性噪聲的Allan方差不超過0.4%，代表了圖10中比差測量值的波動范圍（1）；乘性噪聲反映了圖10中比差測量值隨時間偏離比差真實值的情況。

50Hz正弦信號與幅值曲線變化的乘性噪聲和均值為0正態(tài)分布的加性噪聲疊加如式（12）所示。時域波形如圖12所示。

圖12 時域波形

每10個周波（200ms）計算一次比差測量值，結果如圖13所示。

圖13 比差測量值

比差測量值的Allan方差如圖14所示。

圖14 比差測量值的Allan方差

圖14中，比差測量值既含有加性噪聲又含有乘性噪聲。其中加性噪聲的Allan方差不超過0.4%，代表了圖13中比差測量值的波動范圍（1）；乘性噪聲反映了圖13中比差測量值的曲線變化趨勢。

4 結論

互感器精度評測中比差測量值的總誤差由乘性噪聲引起的誤差和加性噪聲引起的誤差構成。一文在互感器精度評測中引入Allan方差分析法，試驗中充分發(fā)揮Allan方差在隨機誤差辨識方面的性能，數(shù)據(jù)仿真結果實現(xiàn)了互感器乘性噪聲和加性噪聲的分離。

[1] GB/T 20840.8—2007. 互感器(第8部分): 電子式電流互感器[S].

[2] 王巍. 干涉型光纖陀螺儀技術[M]. 北京: 中國宇航出版社, 2010.

[3] 王立輝, 伍雪峰, 孫健, 等. 光纖電流互感器噪聲特性及建模方法研究[J]. 電力系統(tǒng)保護與控制, 2011, 39(1): 62-66.

[4] 張桂才. 光纖陀螺原理與技術[M]. 北京: 國防工業(yè)出版社, 2008.

[5] Allan DW. Statistics of atomic frequency standards[J]. Proceedings of the IEEE, 1966, 54(2): 221-230.

[6] IEEE Std. 1139—2008. IEEE Standard Definitions of Physical Quantities for Fundamental Frequency and Time Metrology[S].

[7] 關巍. 光纖傳感技術在智能電網(wǎng)中的應用[J]. 電氣技術, 2014, 15(3): 115-119, 127.

[8] 馮昊. 諧波對電流互感器誤差的影響[J]. 電氣技術, 2015, 16(8): 45-49.

[9] 產(chǎn)焰萍, 繆希仁. 低壓配電系統(tǒng)短路故障早期檢測的信號濾波方法研究[J]. 電氣技術, 2016, 17(8): 9-13.

[10] 李深旺, 郭志忠, 張國慶, 等. 集磁環(huán)式光學電流互感器的自校正測量方法[J]. 電工技術學報, 2016, 31(9): 197-202.

[11] 李振華, 李紅斌, 張秋雁, 等. 一種高壓電子式電流互感器在線校驗系統(tǒng)[J]. 電工技術學報, 2014, 29(7): 229-236.

[12] 馮利民, 王曉波, 吳聯(lián)梓, 等. 500kV GIS變電站VFTO對于電子式互感器的電磁騷擾研究[J]. 電工技術學報, 2016, 31(1): 85-90.

Research on the Application of Allan Variance in the Accuracy Evaluation of Transformer

Zhang Zhixin1Yang Xulong2Xu Saimei2Zhu Yunwu2Ren Nannan1

(1. Beijing Aerospace Times Optical-electronic Technology Co., Ltd, Beijing 100094; 2. Chuxiong Power Supply Bureau of Yunnan Power Grid Corporation, Chuxiong, Yunnan 675000)

According to the correlation with primary signal, the noise can be divided into two categories: multiplicative noise that is correlated with primary signal and additive noise that is independent of primary signal. A method based on Allan variance is presented in this article. Multiplicative noise and additive noise of transformer are separated by the method. The effectiveness of the method is verified by the results of computer simulation.

transformer; multiplicative noise; additive noise; ratio error; allan variance

張志鑫（1964-），男，工學碩士，研究員，主要研究方向為光學傳感技術的理論研究和工程應用。