谷志茹等

摘要：基于实际测量电力噪声特性参数的统计分布，推导非高斯分布的窄带电力噪声模型，主要包括背景噪声和脉冲噪声.前者应用升余弦滚降滤波法建模，其功率服从正态分布；后者利用随机分布的特性参数：脉冲的宽度、间隔和包络所构建的脉冲函数来建模.宽度、间隔和包络随机量均服从分隔马氏链.然后应用实际测量数据和马尔克夫链蒙特卡罗（Markov Chain Monte Carlo， MCMC）分析法比较评估了此模型，数值分析结果显示，模型能够很好地逼近实际的电力噪声.

关键词：电力线通信； 智能电网； 正交频分复用；分隔马氏链；马尔克夫链蒙特卡罗

中图分类号：TPN14 文献标识码：A

基于电力信道的通信技术是智能电网（Smart Grid）建设的主要通信方式.它不需要重新铺设专用的通信通路，通过在与传递电能同一路径的线路上确定通信链路来传递信息，并且电力线通信不会像无线通信那样容易被密集的建筑物所阻挡，从而减少安装和维护费用.电力线通信技术市场长期被单载波PLC技术所垄断.最近，被称为G3和PRIME标准[1]的OFDM窄带（10～95 kHz） PLC技术进入了市场.它提供高的通信速率，鲁棒的通信模式，这些是智能电网双向通信的关键.但是，作为信息传输媒介来进行数据或语音传输，低压配电网络具有负荷情况复杂、信号衰减大，信道容量小，并存在背景噪声、随机脉冲噪声、与工频同步的谐波噪声等.所以为满足智能电网双向通信的需求，以电力网络为信道的高速通信技术，必须有适应于电力信道的调制和编码方案，能够抵抗恶劣的信道属性等等，这都需要对实际电力线信道有深入的了解.

电力信道中影响高速通信的主要因素有因电缆损耗导致的信号衰减、多径传播和噪声等.其中多径传播在文献[2]中有详细论述，而应用于10～95 kHz的窄带电力线噪声的仿真模型还没有明确的推论.早期的研究显示电力环境中的噪声不是加性高斯白噪声（AWGN）[3].为了仿真噪声，许多文献提出了不同的噪声的统计模型，脉冲噪声的统计属性可以在文献[4]中发现；背景噪声的功率谱密度在文献[5]中建立；一种忽略脉冲噪声随机性的确定性窄带噪声数学模型在文献[6]中给出.然而，目前没有一种电力噪声模型，考虑了脉冲噪声在实际电力环境中包络和时间的分布.

由于电力噪声分布的多样性和随机性，用一种确定的数学模型很难体现其物理属性.本文通过构建真实电力线噪声的测量环境，记录并分析实测噪声，对其中的背景噪声和脉冲噪声进行了建模和仿真.前者的时域模型由高斯白噪声通过滚降系数为‘1的滚降滤波器而得到，其功率近似服从正态分布；后者的时域模型由服从确定统计规律的特性参数所描述，特性参数根据实际测量的数据统计得到，其时域的仿真波形由MCMC分析得到.

1噪声模型

根据背景噪声的平稳分布性和脉冲噪声的随机分布性，利用实际测量数据分别对两种噪声进行统计建模，再将两种噪声叠加[7]，得到实测的电力噪声模型.

1.1噪声测量

文中所述噪声的测量均按图1所示系统实现.电力信号首先通过容性耦合器，将50 Hz的工频电压谐波噪声滤除，然后通过10～500 kHz的带通滤波器，得到窄带电力噪声，最后由大容量数据存储的数字采集器记录保存.滤波器采用3阶巴特沃斯带通滤波器[8]，线性中心频率为255 kHz.数字采集器的型号为USB2085，实现A/D转换和数据采集的功能，设置其采样频率为400 kHz，采集时长为10 min/次.

现选湖南大学13舍实验楼作为测试点，对七楼配电箱A相进行分时段的连续测量记录，根据用电负荷和噪声干扰的不同，分成3个时段进行测量：轻量噪声时段，00：00-06：00；中量噪声时段，12：00-13：00和重量噪声时段，14：00-17：00.测量总时长为6 361 min，约106 h.

1.2背景噪声模型

电力背景噪声属于平稳分布的彩色噪声，其功率主要集中于低频段.图2为去除窄带及脉冲噪声后背景噪声的功率谱密度，3个时段的噪声在频率最低处，功率最大，集中于-10 dB附近，随着频率的增大，功率减小，所以具有低通特性.因此背景噪声的时域模型ηc（t），可以由加性高斯白噪声（AWGN）ξ（t），通过滚降系数为‘1的升余弦滚降滤波器而得到：

1.3脉冲噪声模型

脉冲噪声分为同步的周期性和异步的非周期性脉冲噪声.对于前者，因为具有确定性，所以建模相对简单；但是对于后者，因为时间和包络参数分布的随机性，所以必须用一个随机模型来描述，同时其特性参数服从确定的分布律.

1.3.1脉冲噪声模型

图3为实测脉冲噪声的时域波形，显然，每个脉冲可由3个特性参数描述：包络av，宽度tw和间隔td.其脉冲函数表述如下[5]：

时间/ms

1.3.2参数确定

这一节主要确定上节中所推得脉冲噪声的个数和特性参数的分布规律.

式中NA和Np分别是异步脉冲和周期脉冲的个数.结合式（6），表1可以通过脉冲的来源和数量来解释.在重量噪声时段，因为是工作时段，供电电源以同步于工频电压主频的速度持续动作，所以大量短时宽的周期脉冲占主导，即使脉冲的平均个数达到161个，但是所占时间仍然不足总时长的2%；在轻量噪声时段，实验楼内没有大量的用电设备，这时会有开关的瞬态过程，所以大部分是异步脉冲，周期脉冲可以忽略.因此在由“轻量噪声组”进入“重量噪声组”，周期性脉冲迅速增加.因为测量是在同一点位置进行的，所以3组时段异步脉冲的个数保持不变.基于此NA，NP和周期脉冲的tw_p可以推导如下：

2噪声仿真

电力噪声的仿真主要为背景噪声和脉冲噪声的仿真，再将两部分叠加，得：

其中背景噪声仿真主要是设计升余弦滚降滤波器，其截止频率由式（2）决定，功率谱密度由服从均值为-10，方差为2的正态分布的采样得到.

脉冲噪声的模型，因为参数是服从确定分布律的随机量，所以用MCMC能够得到快速收敛的仿真值.从图4可见，av和tw的测量统计具有快衰落的特点，在20个采样点后PDF（概率密度分布）已经低于10-4，所以用经典MCMC算法：MetropolisHastings （MH） 来仿真这两个量.画出仿真样点直方图，可以看出他们完全按照参数的概率分布模型抽取.

对于脉冲间隔td，概率密度分布于0～16 000 ms，变化缓慢，包括200 000个样点，用经典MCMC算法很难达到收敛，所以采MultipleTryMetropolis（MTM） 算法[11]，与经典MH算法比较，在不降低接受率的条件下，MTM具有更大步长的跳跃.这种算法首先产生几个相关样点，然后从中选择一个最优点.在这种算法中，目标分布和提议分布分别是概率密度函数Pdf_t（x）和正态分布函数N（x）.为得到新的采样，给定初始值为x0.一般，在MCMC算法的权重函数中给出更多的统计信息，能够改善仿真的性能.所以权重函数构造如下：

3模型评价

为了验证提议的模型，将实际电力噪声的测量值与提议模型的仿真值进行比较.在重量组选择5组具有不同（rimp，dr）的测量数据，仿真时长设置为20 s.根据每组（rimp，dr），仿真100组基于提议模型的随机噪声，统计仿真组的（rimp，dr）和均方功率[12]，并与测量值比较.

图5为相同（rimp，dr）条件下，测量值和仿真值的时域波形和功率谱密度，从图中可见，时频域波形比较接近.基于测量数据及仿真数据的干扰率和均

4结论

电力线环境作为通信信道时具有非高斯分布的噪声特性.因此应用高斯噪声作为信道干扰，进而设计和评价电力载波通信算法，其实质并不适应于电力线环境.为此本文根据实际测量电力噪声的统计分析，提出基于随机分布的特性参数所描述的随机模型.其中背景噪声模型由AWGN通过升余弦滚降滤波器得到；脉冲噪声模型由脉冲包络、宽度和间隔所限定的随机模型得到.通过MCMC算法仿真，得到模型的实现，最后通过与实测数据的比较证明了所提议的模型可以用来仿真实际的电力线噪声.

参考文献

[1]JAVIER M， SADOT A， RODRIGUEZ-MORCILLO C. Performance evaluation of two narrowband PLC systems： PRIME and G3 [J].Comput Stand Inter， 2013， 36（1）： 198-208.

[2]ZIMMERMANN M， DOSTERT K. A multipath model for the power line channel [J]. IEEE Trans Commun， 2002， 50（4）：555-559.

[3]GOTZ M， RAPP M，DOSTERT K. Power line channel characteristics and their effect on communication system design [J]. IEEE Communications Magazine， 2004， 42（4）：78-86.

[4]ZIMMERMANN M， DOSTERT K. Analysis and modeling of impulsive noise in broadband powerline communications [J]. IEEE Trans Electromagn Compat， 2002， 44（1）： 249-258.

[5]HOOIJEN O G. A channel model for the residential power circuit used as a digital communications medium [J]. IEEE Trans Power Del， 2004， 19（3）： 1057-1064.

[6]KATAYAMA M， YAMAZATO T， OKADA H. A mathematical model of noise in narrowband power line communication systems [J]. IEEE J Sel Area Comm， 2006， 24（7）： 1267-1276.

[7]ANDREADOU N， PAVLIDOU F N. Modeling the noise on the OFDM powerline communication system [J]. IEEE Trans Power Del， 2010， 25（1）： 150-157.

[8]YAO S N，COLLINS T，JAN P.Hybrid method for designing digital butter worth filters[J]. Comput Electr Eng， 2012， 38（4）： 811-818.

[9]MORRIS H DeGroot， MARK J Schervish. Probability and statistics[M]. 3rd ed. New York：Addison-Wesley， 2002： 219-237.

[10] THOMAS Kaijser. On Markov chains induced by partitioned transition probability matrices [J]. ActaMathematica Sinica， 2011， 27（3）： 221-227.

[11] MARTINO L， DEL Olmo V P， READ J. Amultipoint metropolis scheme with generic weight functions [J]. Statistics and Probability Letters， 2012， 82（7）： 121-134.

[12] 于善奇. 统计方法引论 [M]. 第2 版.北京：北京工业大学出版社，2014：251-267.

YU Shanqi. Introduction to statistical methods [M].2nd ed. Beijing： Press of Beijing University of Technology， 2014：251-267.（In Chinese）

[13] CHEN Shuchuan， LINDSAY Bruce. Improving mixture tree construction using better EM algorithms[J]. Computational Statistics & Data Analysis， 2014， 74（5）： 17-25.