Flicker Parameters Detection Technology Based on Improved Energy Operator of Optimized Sampling Interval and MSLD-SCW Function Spectrum Correction in and Its Application in Wind Farms
-
摘要: 风力发电并网和电力电子技术的广泛发展与应用,造成的电压波动与闪变已成为当前智能电网不容忽视的问题。针对人眼更为敏感的矩形波调制下的电压闪变模型,该文提出改进能量算子和加截断窗谱线插值快速傅里叶变换(FFT) 的闪变包络参数识别法。通过优化能量算子采样间隔实现电压波动分量准确提取,以改进的最大旁瓣衰减速率(MSLD)6项组合余弦窗为母窗构建频域性能优良的MSLD自卷积窗(MSLD-SCW)函数,推导基于新型2阶 MSLD-SCW函数谱线插值校正公式,据此实现矩形方波调制的闪变参数识别与分析。仿真结果表明,优化改进算法在单频矩形方波调制、多频方波调幅波调制、含有谐波与次/超同步间谐波干扰、基频变动及含有噪声等干扰源背景下,相较于传统检测算法均保持较高准确性。最后,将优化算法应用于新疆某地区电网电压闪变识别,验证其有效性。Abstract: With the wide development and application of grid-connected wind power generation and power electronics technologies, the voltage fluctuation and flicker problem can not be ignored in smart grid. Considering the voltage flicker model under rectangular wave modulation which is more sensitive to human eyes, a method for detecting voltage flicker envelope parameters based on improved energy operator and windowed interpolation Fast Fourier Transform (FFT) is proposed. By optimizing the sampling interval of energy operator, the voltage fluctuation components can be accurately extracted. A Maximum Side-Lobe Decay(MSLD) Self-Convolution Window(MSLD-SCW) function with excellent frequency domain performance can be obtained by improved six-term combined cosine window with MSLD, the rectification formula of spectral line interpolation based on new MSLD-SCW is derived, and the detection and analysis of flicker parameters are realized accordingly. The simulation results indicate that the optimization and improvement algorithm can effectively maintain higher detection accuracy than the traditional methods in the case of single frequency rectangular wave modulation, multi-frequency modulation, containing harmonics and sub/supersynchronous inter-harmonics interference, power grid fundamental frequency deviation and noise interference. Finally, the optimization algorithm is applied to voltage flicker envelope parameters detection in a certain area of XinJiang, and the effectiveness of the algorithm is verified.
-
1. 引言
随着大数据和5G技术的快速发展,数据安全变得越来越重要,而保障数据安全的最好方式之一是运用密码学对数据进行加密处理。在密码学应用中,无论是密码算法中密钥的生成或是密码协议中特定变量的随机初始化,都需要用到真随机数源。真随机数在统计学上具有随机性,包括时间上的独立性和空间上的均匀性,还具有不可重复性和不可预测性[1-3]。随机数主要应用于密码算法协处理器中的密钥、身份认证和数字签名等,这些都需要真随机数发生器输出的随机系列具有熵源高,功耗低,面积小等特点[4,5]。
目前用于产生随机数的硬件实现方法:(1)离散时间混沌的方法[6,7];(2)采用相位抖动的方法[8];(3)采用SRAM的方法[9];(4)热噪声的方法[10]。
方法(1)所产生的随机数本质上是伪随机数,是通过数学算法实现的,不是真随机数;方法(2)利用相位抖动或振荡器的漂移作为随机源,但靠这一点不能达到足够的随机性能;方法(3)产生的随机数性能与工艺相关性很大,随机性能不能保证;方法(4)利用热噪声是由大量自由电子运动产生的,其统计特性服从高斯分布,产生的随机数性能好。
本文通过对一种低功耗高噪声源真随机数发生器的研究,设计了一种新型的低频时钟电路,可以把电阻热噪声放大100倍以上,从而减少低频时钟电路的带宽和电阻值,使电路的面积和功耗减少,并且使低频时钟的jitter到达58.2 ns。随机数输出满足AIS31真随机数熵源测试要求,且通过了国密2安全测试。
2. 真随机数电路结构
真随机数电路结构如图1所示,电路由低频时钟、高频时钟、触发器、后处理等电路构成[11-13]。在真随机数电路设计过程中,通常要求低频时钟的抖动标准差(jitter)在高振荡器周期的10~20倍之间,以提高真随机数发生器的抗干扰能力以及输出序列的随机性能,低频时钟的jitter值与随机数的随机性能正相关。要使低频时钟的jitter值变大,一般情况下,会增大低频时钟的噪声电阻和噪声带宽,增大噪声电阻会增加电路面积,噪声带宽会增大电路的功耗,故低频时钟的jitter值与电路的面积和功耗相互制约;高频时钟频率输出频率也会影响随机数的随机性能,比如低频时钟的抖动标准差为20 ns,根据通常要求慢振荡器抖动标准差在高振荡器周期的10~20倍之间,以15倍计算,则要求高速时钟输出频率为0.75 GHz,这样会大大增加高频时钟的功耗,为了减少功耗,增大低频时钟的jitter值是必然要求,此时又要增大电路的面积和功耗[14]。设计随机数会在面积和功耗与输出序列的随机性之间进行折中。
3. 新型低频时钟设计
3.1 低频时钟原理和jitter理论分析
低频时钟带有抖动的慢振荡器(clkslow),电路如图2所示。低频时钟模块包括跨阻放大器、跨导放大器、迟滞比较器、电荷泵、基准电压VREF等电路。
电路使能开启后,偏置电路开始提供偏置电压和电流。当SCLK_OUT为高电平时,电荷泵放电使跨导放大器正端电压下降,跨导放大器输出电流降低,输出的电流通过跨阻放大器转换为电压,当跨阻放大器输出电压到达低阈值-VTL, SCLK_OUT变为低电平;当SCLK_OUT为低电平时,电荷泵冲电使跨导放大器正端电压上升,跨导放大器输出电流增大,输出的电流通过跨阻放大器转换为电压,当跨阻放大器输出电压到达低阈值+VTL,SCLK_OUT变为高电平,这样就可以得到低频时钟输出信号SCLK_OUT。跨阻放大器输出信号是一个在迟滞比较器高低阈值间来回摆动的三角波;电阻R1和R2上的热噪声经过运放放大后叠加在三角波上,得到SCLK_OUT 的时钟沿抖动与热噪声一样,满足正态分布。
跨阻放大器输出的三角波信号如图3所示,其中S是运放输出三角波的斜率,Tcl是慢时钟信号的周期。可以得到
Tcl=t1+t2 (1) 其中t1和t2是随噪声电压而独立随机变化的。同时可得
V(t)=−VTL+S⋅t+Vn(t) (2) S·t是三角波的电压,V(t)是输出电压随时间变化的函数,Vn(t)是放大后的电阻热噪声,可以推出
t1=(VTH+VTL−Vn(t))/S (3) 因为Vn(t)的平均值为0,那么
E{Tcl}=2(VTH+VTL)/S (4) σ{Tcl}=√2σ{Vn}/S (5) 其中E{Tcl}是clkslow的平均周期,σ{Tcl}和σ{Vn}分别是clkslow 的jitter和放大后噪声电压的标准差。由上式可知,σ{Tcl}与σ{Vn}成正比,与S成反比。σ{Vn}和S这两个值都是与电路的某些参数相关的,可以得到
S=±(Ich⋅A1)/C (6) A1=GmI1IR3 (7) Ich是电荷泵充放电电流,A1是电荷泵的输出信号到跨阻放大器的输出信号之间增益,即信号增益。Gm是跨导放大器的跨导,I是跨导放大器的输出电流,I1是流过电阻R3的电流,R3是跨阻放大器的增益,C是电荷泵到地的电容。
σ{Vn}=√(4kT⋅Bw⋅2Rno⋅A22) (8) Rno是噪声电阻,A2是噪声电压增益,Bw是噪声带宽(保守计算时可用主极点频率代替)。
由式(4)、式(6)可得低频时钟周期公式
E{Tcl}=2(VTH+VTL)⋅CIchA1 (9) 从式(9)可知,低频时钟频率与比较器的迟滞范围VTH+VTL、电容C、信号增益A1、电荷泵充放电电流Ich相关,如果低频时钟频率确定,其相应的参数也会确定。
由式(5)、式(6)、式(8)可得低频时钟jitter的σ{Tcl}的具体公式
σ{Tcl}=√2×√4kTBW⋅2Rno⋅A2⋅CIchA1 (10) 3.2 跨阻放大器的设计
从图2可知,噪声电阻在跨阻放大器的两端,噪声是通过跨阻放大器放大的,跨阻放大器的带宽即噪声带宽。从式(10)可知,σ{Tcl}值和噪声电阻Rno、噪声带宽Bw、噪声电压增益A2、迟滞范围VTH+VTL、电容C、信号增益A1、电荷泵充放电电流Ich相关,当低频时钟确定之后,其参数迟滞范围VTH+VTL、电容C、信号增益A1、电荷泵充放电电流Ich会跟着确定,要改变σ{Tcl}值,只有噪声电阻Rno、噪声带宽Bw、噪声电压增益A2可以改变,这3个参数都和跨阻放大器相关,故跨阻放大器的设计是本文的关键。
在图2有跨阻放大器的原理图,电流信号通过电阻R3转换成电压信号,跨阻放大器的闭环增益为R3,噪声电阻为R1, R2,电流信号Iin分成两路,一路通过R3转换成电压信号,一路通过电阻R4。
根据跨阻放大器的原理图求噪声电压增益,噪声电阻R1的噪声电压
vn2=4kT⋅Bw⋅R1 (11) 噪声电流
in2=4kT⋅Bw⋅R1R42 (12) 噪声电流加载在Iin信号上,通过R3进行放大,即输出噪声电压为
vnout2=4kT⋅Bw⋅R1R42R32 (13) 跨阻放大器放大后噪声电压的标准差
σ{Vnout}=√4kT⋅Bw⋅2Rno⋅R3R4 (14) 即噪声电压增益
A2=R3R4 (15) 在这里噪声电压增益可以取100~200倍之间,这样可以在噪声电阻和噪声带宽在比较小的情况下,得到比较大的噪声电压标准差,即可以设计出比较大的低频时钟jitter,增大随机数的性能,并且可以减少功耗和面积。
图4是普通的电压放大器放大电阻噪声运原理图,这种结构信号增益和噪声增益相同,根据式(10)可知,要增大低频时钟jitter,必须增大噪声电阻和噪声带宽,这样会使电路面积和功耗大大增加。
表1为这两种结构下低频时钟jitter为60 ns,时钟频率为2 MHz下,需要的噪声电阻值和功耗。
表 1 两种结构下噪声电阻值和功耗指标 噪声电阻值
(Ω)信号增益A1(倍) 信号增益A2(倍) 噪声带宽(MHz) 功耗(mW) 跨阻放大器结构 64 k 1 100 1 0.081 电压放大器结构 2 M 5 5 80 0.220 通过对两种结构的比较可知,新结构用较小的功耗和面积实现了较大的低频时钟jitter。
3.3 低频时钟仿真
本文在设计真随机数产生器时,使输出吞吐率范围为1.38~3.34 Mbit/s,典型值为2.13 Mbit/s。
表2为低频时钟输出频率对应的σ{Tcl}仿真结果,从表1中可知,σ{Tcl}的典型值为58.2 ns,功耗为74 uA。
表 2 低频时钟频率仿真结果指标 仿真结果 MIN TYP MAX 输出频率(MHz) 1.38 2.13 3.34 Jitter(σ{Tcl})(ns) 77.89 58.2 40 功耗(mW) 0.055 0.081 0.110 图5 为低频时钟频率在2.13 MHz时的jitter的分布情况;图6是对jitter值的分布情况做正态分布处理,从图6中可知低频时钟的jitter值服从正态分布,标准差σ{Tcl}值为58.2 ns。
4. 低功耗高频时钟设计
在上面的章节中可知时钟jitter值比较大,其值为58.2 ns,这样就可以减少高频时钟频率,从而减少功耗。
为了提高真随机数发生器的抗干扰能力以及输出序列的随机性能,快时钟振荡器的周期是慢振荡器jitter的1/20-1/10。取高频时钟的周期是σ{Tcl}的1/15,快时钟振荡器的周期为3.88 ns,时钟频率为257 MHz,这样可以减少高频时钟的功耗。
由于高频时钟的输出中心频率要达到257 MHz,本文中高频时钟的电路采取环路振荡器的结构,晶体管M1控制环路的总电流,用来减少功耗,为保证在PVT变化时输出高频时钟的占空比为50%,环振的输出信号需经过一个高速二分频电路。二分频电路采取高速二分频电路(TSPC)的结构。
高频时钟的仿真结果表3,其时钟频率的中心频率为250 MHz,占空比为50.28%。
表 3 高频时钟频率仿真结果仿真 HOSC频率 MIN TYP MAX 频率(GHz) 0.186 0.25 0.313 功耗(mW) 0.017 0.024 0.035 占空比(%) 50.09 50.28 50.43 5. 整体版图设计和测试
本文设计的真随机数发生器电路采用SMIC 40 nm CMOS工艺,核心电路版图面积小于0.00789 mm2,图7为整体电路的版图设计,分为低频时钟模块、高频时钟模块、基准电压模块。
基于噪声放大的真随机数发生器整体电路测试结果为输出时钟频率为2.4 MHz左右,功耗为0.11 mW,实现了低成本低功耗真随机数的设计,输出的随机系列通过了AIS31测试,并且过了国密2安全检测。
表4是流片测试结果和文献调研中参考的国外相关真随机数发生器的性能参数的对比和具体参考文献,本文的结果在最后一行。表的最后一列是归一化参数S的结果比较,参数S的定义为
S=面积×功耗速度×1000 (16) 如表4所示,利用相位抖动原理产生的随机数虽然面积和功耗都较低,但是其随机数的速度较慢,因此其应用场景会受限。在很多实际应用场景中,对随机数的速度都有要求,本文所设计的随机数主要运用在高速安全芯片和主控芯片中,因此随机数的速度至上要达到1 MHz以上,速度太慢,会影响芯片的性能。
6. 结束语
本文设计基于电阻热噪声振荡器的真随机数产生器,电路包括了低频时钟电路、高频时钟电路、D触发器等模块。设计了一种新型的低频时钟电路,可以把电阻热噪声放大100倍以上,从而减少低频时钟电路的带宽和电阻值,使电路的面积和功耗减少。该电路结构可以保证获得较大的周期抖动从而减少电路的面积和功耗。在本设计中使低频时钟的jitter到达58.2 ns,从而减少高频时钟频率,进一步减少功耗,最终使整体电路功耗为0.11 mW,面积为0.00789 mm2。低频时钟输出噪声符合白噪声分布,随机性较好,随机数输出结果满足AIS31随机性测试,并且过了国密2安全检测。
本文设计的随机数在面积和功耗方面有了很大的提升,可以应用在信息安全、计算随机模拟、数字系统内置的检测性能和电子商务系统等领域。
-
表 1 多频率包络调制仿真参数设置
闪变幅值(%) 闪变频率(Hz) 0.10 22 0.08 12 0.04 5 -
[1] 李加升, 杨金辉, 吴顺秋, 等. 风能发电并网时引起的电压波动与闪变检测的仿真研究[J]. 电测与仪表, 2010, 47(10): 32–35. doi: 10.3969/j.issn.1001-1390.2010.10.008LI Jiasheng, YANG Jinhui, WU Shunqiu, et al. Simulation for voltage fluctuation and flicker when wind power grid detection[J]. Electrical Measurement &Instrumentation, 2010, 47(10): 32–35. doi: 10.3969/j.issn.1001-1390.2010.10.008 [2] 张超, 王维庆, 邱衍江, 等. 大规模风电并网地区次同步谐波检测方法[J]. 高电压技术, 2019, 45(7): 2194–2202. doi: 10.13336/j.1003-6520.hve.20181207008ZHANG Chao, WANG Weiqing, QIU Yanjiang, et al. Detection method of subsynchronous harmonic in regions with large scale wind power paralleled in grid[J]. High Voltage Engineering, 2019, 45(7): 2194–2202. doi: 10.13336/j.1003-6520.hve.20181207008 [3] 李芊芊, 张江滨, 李亚琳. 适用于双馈风机的并网电压波动与闪变检测系统的研究[J]. 电网与清洁能源, 2015, 31(12): 126–131. doi: 10.3969/j.issn.1674-3814.2015.12.023LI Qianqian, ZHANG Jiangbin, and LI Yalin. Research on grid-connected voltage fluctuation and flicker detection system suitable for doubly-fed fan[J]. Power System and C1ean Energy, 2015, 31(12): 126–131. doi: 10.3969/j.issn.1674-3814.2015.12.023 [4] 张超, 王维庆, 王海云, 等. 风火打捆外送系统220kV电网次同步振荡监控策略研究[J]. 电力系统保护与控制, 2018, 46(11): 138–144. doi: 10.7667/PSPC170755ZHANG Chao, WANG Weiqing, WANG Haiyun, et al. Study on monitoring and control strategy of subsynchronous oscillation for 220 kV power grid in wind-thermal-bundled transmission system[J]. Power System Protection and Control, 2018, 46(11): 138–144. doi: 10.7667/PSPC170755 [5] 倪良华, 肖李俊, 历馨, 等. 基于HHT与闪变功率流向法的闪变源定位识别[J]. 电网与清洁能源, 2018, 34(1): 26–31. doi: 10.3969/j.issn.1674-3814.2018.01.005NI Lianghua, XIAO Lijun, LI Xin, et al. Flicker source tracing based on HHT and flicker power flow direction method[J]. Power System and Clean Energy, 2018, 34(1): 26–31. doi: 10.3969/j.issn.1674-3814.2018.01.005 [6] 倪良华, 肖李俊, 吕干云, 等. 基于改进HHT的电压闪变检测[J]. 电测与仪表, 2017, 54(18): 72–78. doi: 10.3969/j.issn.1001-1390.2017.18.012NI Lianghua, XIAO Lijun, LÜ Ganyun, et al. Voltage flicker parameters detection based on modified HHT[J]. Electrical Measurement &Instrumentation, 2017, 54(18): 72–78. doi: 10.3969/j.issn.1001-1390.2017.18.012 [7] 石黄霞, 付涛. 基于小波变换的单频时变电压闪变检测仿真研究[J]. 电气自动化, 2020, 42(3): 52–54. doi: 10.3969//j.issn.1000-3886.2020.03.016SHI Huangxia and FU Tao. Simulation study on single-frequency time-varying voltage flicker detection based on wavelet transform[J]. Electrical Automation, 2020, 42(3): 52–54. doi: 10.3969//j.issn.1000-3886.2020.03.016 [8] 贾峰, 孙曼, 杨勰颖. 基于改进小波变换的电力谐波检测方法研究[J]. 计算机仿真, 2017, 34(1): 99–103. doi: 10.3969/j.issn.1006-9348.2017.01.021JIA Feng, SUN Man, and YANG Xieying. Power harmonic detection method research basedon improved wavelet transform[J]. Computer Simulation, 2017, 34(1): 99–103. doi: 10.3969/j.issn.1006-9348.2017.01.021 [9] 吕干云, 武阳, 时翔, 等. 一种基于S变换的闪变检测和源定位方法[J]. 电力电容器与无功补偿, 2018, 39(5): 87–93. doi: 10.14044/j.1674-1757.pcrpc.2018.05.016LÜ Ganyun, WU Yang, SHI Xiang, et al. Flicker detection and source location method based on S-transform[J]. Power Capacitor &Reactive Power Compensation, 2018, 39(5): 87–93. doi: 10.14044/j.1674-1757.pcrpc.2018.05.016 [10] 王粉花, 谢斌, 王华涛. 基于FPGA的快速差频测量系统设计[J]. 电子与信息学报, 2019, 41(1): 187–194. doi: 10.11999/JEIT180243WANG Fenhua, XIE Bin, and WANG Huatao. Design of fast differential frequency measurement system based on FPGA[J]. Journal of Electronics &Information Technology, 2019, 41(1): 187–194. doi: 10.11999/JEIT180243 [11] 古庭赟, 高云鹏, 吴聪, 等. 基于改进能量算子和六项余弦窗频谱校正的电压闪变包络参数检测[J]. 电力系统保护与控制, 2019, 47(23): 44–51. doi: 10.19783/j.cnki.pspc.181600GU Tingyun, GAO Yunpeng, WU Cong, et al. Voltage flicker envelope parameter detection based on improved energy operator and six-term cosine window spectrum correction[J]. Power System Protection and Control, 2019, 47(23): 44–51. doi: 10.19783/j.cnki.pspc.181600 [12] 罗瑜, 张珍珍. 一种方向插值预测变长编码的帧存有损压缩算法[J]. 电子与信息学报, 2019, 41(10): 2495–2500. doi: 10.11999/JEIT181195LUO Yu and ZHANG Zhenzhen. A lossy frame memory compression algorithm using directional interpolation prediction variable length coding[J]. Journal of Electronics &Information Technology, 2019, 41(10): 2495–2500. doi: 10.11999/JEIT181195 [13] 陈至豪, 王立德, 王冲, 等. 基于组合余弦优化窗四谱线插值FFT的电力谐波分析方法[J]. 电网技术, 2020, 44(3): 1105–1113. doi: 10.13335/j.1000-3673.pst.2019.0075CHEN Zhihao, WANG Lide, WANG Chong, et al. An approach for electrical harmonic analysis based on optimized composite cosine window four-spectrum-line interpolation FFT[J]. Power System Technology, 2020, 44(3): 1105–1113. doi: 10.13335/j.1000-3673.pst.2019.0075 [14] 曾博, 唐求, 卿柏元, 等. 基于Nuttall自卷积窗的改进FFT谱分析方法[J]. 电工技术学报, 2014, 29(7): 59–65. doi: 10.3969/j.issn.1000-6753.2014.07.009ZENG Bo, TANG Qiu, QING Baiyuan, et al. Spectral analysis method based on improved FFT by nuttall self-convolution window[J]. Transactions of China Electrotechnical Society, 2014, 29(7): 59–65. doi: 10.3969/j.issn.1000-6753.2014.07.009 [15] 孙仲民, 何正友, 臧天磊. 一种混合卷积窗及其在谐波分析中的应用[J]. 电工技术学报, 2016, 31(16): 207–214. doi: 10.3969/j.issn.1000-6753.2016.16.024SUN Zhongmin, HE Zhengyou, and ZANG Tianlei. A kind of hybrid convolution window and its application in harmonic analysis[J]. Transactions of China Electrotechnical Society, 2016, 31(16): 207–214. doi: 10.3969/j.issn.1000-6753.2016.16.024 [16] 崔嘉, 杨俊友, 李连富, 等. 含分布式风电场的配电网电压波动分析与抑制[J]. 电网技术, 2015, 39(12): 3414–3420. doi: 10.13335/j.1000-3673.pst.2015.12.013CUI Jia, YANG Junyou, LI Lianfu, et al. Voltage fluctuation analysis and mitigation of distribution network containing distributed wind farm[J]. Power System Technology, 2015, 39(12): 3414–3420. doi: 10.13335/j.1000-3673.pst.2015.12.013 期刊类型引用(4)
1. 王鹏,邹彬,刘金枝,周丹阳. 基于Xilinx型FPGA系统单粒子效应评估方法研究. 电子学报. 2022(11): 2716-2721 . 百度学术
2. 蔡阳阳,陶伟,郭刚,魏敬和. CAN控制器单粒子效应测试系统的研制. 中国电子科学研究院学报. 2019(04): 387-392 . 百度学术
3. 王鹏,芦浩,刘金枝,薛茜男,金志威. 基于协同进化的航空高度单粒子翻转故障生成方法研究. 现代电子技术. 2019(16): 112-116+121 . 百度学术
4. 陈晓娟,陈东阳,吴洁. CMOS反相器低频噪声模型及可靠性表征研究. 电子学报. 2016(11): 2646-2652 . 百度学术
其他类型引用(5)
-