CN103293363B - 一种互感器采样值延时补偿方法 - Google Patents
一种互感器采样值延时补偿方法 Download PDFInfo
- Publication number
- CN103293363B CN103293363B CN201310276475.9A CN201310276475A CN103293363B CN 103293363 B CN103293363 B CN 103293363B CN 201310276475 A CN201310276475 A CN 201310276475A CN 103293363 B CN103293363 B CN 103293363B
- Authority
- CN
- China
- Prior art keywords
- value
- delay
- mutual inductor
- sampled
- electrical angle
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Complex Calculations (AREA)
Abstract
本发明提出了一种互感器采样值延时补偿方法。所述方法首先通过采样值映射关系测量出延时,再利用三点乘积法得出采样点幅值的平方值,然后利用CORDIC算法双曲向量模式进行开方得出采样值的虚部,最后通过CORDIC算法圆周旋转模式根据测算出的数据延时电角度对原采样值进行旋转补偿,得到补偿电角度后的采样值。本发明方法能够在合并单元数据处理时间的不确定的情况下,对互感器数据采样延时进行实时动态补偿,得到补偿电角度后的采样点参数,方法相对简单使用,在工程应用中取得了较好的效果。
Description
技术领域
本发明属于电工技术领域,尤其是涉及一种互感器采样值延时补偿方法。
背景技术
在数字化变电站与智能变电站的建设中,电子式互感器作为一次测量设备起着关键作用,采样数据的准确有效是整个变电站稳定运行的基础。电子式互感器的输出精度主要包括相位和幅值两方面,幅值精度取决于变比系数,设计时很容易考虑到,而相位精度较为复杂,与***的电路设计、软件处理都有关系。因此有效改善并提高相位精度,对电网的稳定与发展具有十分重要的意义。
传统的短数据窗移项算法是根据额定数据延时Tdr,对采样电流/电压数据补偿固定的额定延时位移θ。
设t时刻的电流或电压采样值可表示为:
x(t)=Xcos(ωt+δ0) (1)
设两个采样点之间的时间间隔为TS,则上一个采样点t-TS时刻电流采样值为
x(t-TS)=Xcos(ωt+δ0-ωTS) (2)
补偿额定延时Tdr,即向后旋转电角度补偿后的采样值为
xtdr(t)=Xcos(ωt+δ0+θ) (3)
设
其中a、b为系数,三个时刻的采样值对应的向量分别为 其中,∠表示角度符号,∠后面是角度。
在向量三角形中,由正弦定理可推算出
由于电流/电压采样值在补偿前与补偿后只是经过θ电角度的旋转,幅值均为X,因此系数a、b可直接由式(5)算出
而ωTS由每工频周期采样点数N及谐波次数n确定,ωTS=2πn/N。针对基波进行补偿时,n=1,即ωTS=2π/N。因此系数a、b的值可预先计算出来,在上位机对处理芯片进行配置时,动态配置到EEPROM中,以便芯片每次启动时读取系数a、b的值进行延时补偿计算。
由式(4)得:xtdr(t)=a·x(t)-b·x(t-TS) (7)
因此处理芯片在进行计算时,只需用到上一个点的采样值x(t-TS)与本次采样值x(t),经过式(7)的两次乘法计算,一次减法计算就可计算出补偿后的电流采样值xtdr(t)。
现在的补偿算法都是根据预先估计的额定延时位移进行定制补偿,但当线路改造时传感器与合并单元的光纤通信线发生改动、传感器的晶振由于室外温度变化而产生误差、合并单元数据处理时间的不确定等因素造成数据延时产生扰动时,补偿算法便会造成较大的误差。而且每次重新安装时都需要对合并单元的额定数据延时重新测量配置,故该算法有很大的局限性。
因此,发明一种更为有效地提高互感器的相位精度的新方法成为亟需解决的课题。
发明内容
本发明所要解决的技术问题在于克服现有技术的不足,针对静态的额定延时补偿方法带来的补偿算法造成较大误差的问题,本发明提出了一种互感器采样值延时补偿方法,该方法基于坐标旋转数字计算CORDIC算法,通过实时测量数据延时,并对延时位移进行动态补偿的方法,有效提高了互感器的相位精度。
为解决上述技术问题,本发明采用的技术方案如下:
一种互感器采样值延时补偿方法,所述方法针对互感器采样延时的不确定性,首先通过采样值映射关系测量出延时,再利用三点乘积法得出采样点幅值的平方值,然后利用CORDIC算法双曲向量模式进行开方得出采样值的虚部,最后通过CORDIC算法圆周旋转模式根据测算出的数据延时电角度对原采样值进行旋转补偿,可以得到补偿电角度后的采样值;具体步骤如下:
步骤A,利用高压端传感器与互感器合并单元之间的采样值映射关系,实时测量采集数据的延时Td,进而得出数据延时电角度θ;
步骤B,计算采样值幅值的平方值X2;
其中,x为当前时刻的采样值,x1,x2分别为前两个时刻的采样值,Ts为采样周期,sin2ωTS为定值;
步骤C,计算采样值虚部值y的平方值y2:
y2=X2-x2
步骤D,计算采样值虚部y:
将(y2+1/4,y2-1/4)代入CORDIC算法双曲向量模式中进行迭代,得到采样值虚部y;
步骤E,将(x,y)代入CORDIC算法圆周旋转模式中进行迭代,得经延时补偿θ电角度后的坐标(x',y'),从而得到补偿电角度后的采样值。
所述步骤A中,所述数据延时电角度θ,其计算公式为:
其中,L'为在Td时间内互感器合并单元晶振振动次数,L为360°电气角度对应的互感器合并单元晶振振动次数。
所述步骤D中,所述CORDIC算法双曲向量模式,经n次迭代后坐标关系为
式中,Sn=-sgn(Yn),CORDIC算法双曲向量模式比例因子在迭代次数n确定后K′n为常数;预先对迭代初始值(x0,y0)进行比例修正,乘以比例因子K′n得到修正后的初始值(x0′,y0′)再进行迭代计算,ti满足tanh(ti)=2-i,i为自然数;
其中,sgn(Yn)为符号运算,Yn表示n次迭代后采样值的虚部,Yn为正实数时,sgn(Yn)值为1;Yn为负实数时,sgn(Yn)值为-1;Yn为0时,sgn(Yn)值为0;
逐次迭代后y趋于0,最终收敛于x轴上的向量(r,0),有
所述步骤E中,所述CORDIC算法圆周旋转模式,经n次旋转后坐标关系为
其中,CORDIC算法圆周旋转模式比例因子Kn在迭代次数n确定后为常数,且n→∞时,Kn→0.858784;Sn=sgn(Zn),
Z0=θ;θi满足tanθi=2-i,i为自然数。
本发明的有益效果是:本发明提出了一种互感器采样值延时补偿方法,所述方法采用的CORDIC算法只用到加法器和移位器,电路实现极为简单,而且CORDIC算法迭代计算采用流水线结构,整个处理器被精简成一个内部互连的加法器阵列,提高了***工作的速度。基于CORDIC算法采样值延时补偿方法能够在合并单元数据处理时间的不确定的情况下,对互感器数据采样延时进行实时动态补偿,得到补偿电角度后的采样点参数,方法相对简单使用,在工程应用中取得了较好的效果。
附图说明
图1是CORDIC算法圆周旋转示意图。
图2是CORDIC算法双曲向量示意图。
图3是CORDIC算法迭代结构图。
图4是CORDIC算法流水线结构图。
图5是延时补偿算法流程图。
图6是编程计算流程图。
具体实施方式
下面结合附图,对本发明提出的一种互感器采样值延时补偿方法进行详细说明:
1、CORDIC算法
为了更好的理解本发明,首先对CORDIC算法中所涉及到的部分原理进行介绍。
1)CORDIC算法圆周旋转模式
在二维向量空间中,给定向量(x0,y0),则根据公式(1)可以得到按照顺时针旋转θ电角度后到达坐标(x1,y1),如图1所示。
式(1)可改写为
选择一系列满足tanθi=2-i的特殊角θi(i=1,2,3…)进行不断的旋转,则乘以tanθi就转化为乘以2-i,即可以通过简单的移位运算逐步旋转±θi逼近θ角。由于是对延时造成的相位移动进行补偿,因此补偿的θ角较小,可忽略θ0=45°时的旋转,虽然缩小了θ角的补偿范围(θ<54.827°),仍可满足补偿要求,在硬件资源和处理时间上省略了一次迭代。
经n次旋转后转角误差为
式中Sn=sgn(Zn),转角误差Zn小于arctan2-n-1,代表向量的旋转方向,1代表顺时针旋转,-1代表逆时针旋转。由式(3)可以看出,迭代的次数越多,转角误差越小,转角越趋近θ。
经n次旋转后坐标关系为
上述计算可分为两部分,第一部分为比例因子第二部分为迭代计算,可记为
比例因子Kn在迭代次数n确定后为常数,且n→∞时,Kn→0.858784,即旋转后造成了半径的增大,需要对坐标进行比例因子Kn的缩小补偿。由于合并单元的数字信号处理模块中有定标模块,通过数据定标将采集到的数字化电气量重新标定以便与一次侧的电气参量相符,即使得保护用电流互感器二次侧额定输出为0x01CF,测量用电子式电流互感器和电子式电压互感器的二次侧额定输出均为0x2D41。因此该比例因子的补偿计算可以整合到定标模块中实现。
2)CORDIC算法双曲向量模式
设向量(x0,y0)=(r·cosht,r·sinht)沿双曲线x2-y2=r2移动到另一向量(x1,y1),如图2所示。双曲旋转定义为
式(6)可改写为
选择一系列满足tanh(ti)=2-i的ti进行迭代计算(i=1,2,3,…),则乘以tanh(ti)就转化为乘以2-i,可通过移位运算实现。
经n次迭代后坐标关系为
式中Sn=-sgn(Yn),Yn表示n次迭代后采样值的虚部,比例因子
在迭代次数n确定后K′n为常数。因此可预先对初始值(x0,y0)进行比例修正,乘以比例因子K′n得到修正后的初始值(x′0,y′0)再进行迭代计算。
令(xn,yn)=(r,0),则逐次迭代后y→0,最终收敛于x轴上的向量(r,0),有
若令x0=W+14,y0=W-14,则可得
2、三点乘积法
取任意三个连续的采样时刻t-2TS、t-TS、t,其采样值分别对应为x1、x2、x。
有x1=Xcos(ωt+δ0-2ωTS)=Xcos(δ1-2ωTS) (10)
x2=Xcos(ωt+δ0-ωTS)=Xcos(δ1-ωTS) (11)
x=Xcos(ωt+δ0)=Xcosδ1 (12)
由于
=X2(cos2δ1cos2ωTS+sin2δ1sin2ωTS+sinδ1cosδ1sin2ωTS (13)
xx1=X2cosδ1cos(δ1-2ωTS) (14)
=X2(cos2δ1cos2ωTS-cos2δ1sin2ωTS+sinδ1cosδ1sin2ωTS
x2-xx1=X2(sin2δ1sin2ωTS+cos2δ1sin2ωTS)=X2sin2ωTS (15)
所以有
其中sin2ωTS为定值,可预先求出。
3、方法介绍
本发明的整个延时补偿算法的流程图如图5所示,首先接收采样数据,测出数据延时Td,并计算出数据延时电角度θ。然后用三点乘积法求出信号的幅值平方X2,接着计算出y2=X2-x2,随后将(y2+1/4,y2-1/4)代入CORDIC算法双曲向量模式中迭代,逐步令纵坐标趋于0,迭代n次后可得收敛于x轴上的向量(y,0)。接着将(x,y)代入CORDIC算法圆周旋转模式中进行迭代,取Z0=θ,逐步令Z趋于0,迭代n次后可得经延时补偿θ电角度后的坐标(x',y')。最后将补偿后的采样值写入IEC61850-9-1/2报文中进行报文发送。
CORDIC算法FPGA中的迭代结构图如图3所示,根据所需的工作模式,决定旋转方向控制因子Sn。旋转模式中Sn=sgn(Zn),向量模式中Sn=-sgn(Yn)。实际工作时,x、y、z寄存器中的初始值分别为x0、y0、θ。在第n个时钟周期时,x、y寄存器内的值右移n位后进入加/减法器,计算出的结果再分别存入寄存器,进行迭代计算。而ROM存储器中存储了第n次迭代时补偿的固定角度值,将第n个地址中的角度值与z寄存器内的值一起调入减法器,计算出剩余待补偿的角度再存入z寄存器,进行迭代计算。但该迭代结构中,移位寄存器的移位参数每次迭代都要修改,ROM的查找地址每次迭代也要递增,因此需要较大的扇入能力,输出时间由迭代次数决定,造成延时较大。
将上述迭代计算处理器进行展开,如图4所示,n个处理单元中每个处理单元都可以同时并行处理一个相同的迭代算法。该流水线方式结构中各个移位寄存器上的移位次数固定,角度值作为各个角度累加器中的一个常量输入,常量可以用硬件连线来代替存储空间,因此整个处理器被精简成一个内部互连的加/减法器阵列,输出延时为组合逻辑电路的硬件延时,提高了***工作的速度。
实施例
本发明采用Xilix的Spartan3系列芯片进行编程测试,选择FPGA作为采样主控芯片,FPGA丰富的I/O口和快速的并行处理速度,行程高速实时数据流。整个编程计算运行流程图如图6所示。
将本发明所提出的算法经编程下载到电子式电流互感器合并单元中,采用江苏凌创NT700电子式互感器稳态校验***进行了运行测试,额定测量电流为5A,实验部分数据结果如表1和表2所示。其中绝对延时Tad为额定延时Tdr与实际延时Td的差值。
Tad=Td-Tdr (17)
表1传统延时补偿试验数据
表2动态延时补偿试验数据
表1试验数据采用的是传统的延时补偿算法,对额定延时进行定值补偿。考虑到***相位偏移的补偿误差与延时补偿的计算误差,其数据结果基本满足式推导所得的1us延时对应1.08′相位差的结论,因此补偿后的相位差主要与绝对延时有关。
表2为同等条件下,采用了CORDIC动态延时补偿算法的试验数据。可以看出,由于采用了动态补偿,绝对延时也计算并补偿在内,因此相位差得到了明显的改善,这时相位差主要与采样精度和计算误差有关,其精度可以达到0.1级要求。
Claims (3)
1.一种互感器采样值延时补偿方法,其特征在于,所述方法首先通过采样值映射关系测量出延时,再利用三点乘积法得出采样点幅值的平方值,然后利用CORDIC算法双曲向量模式进行开方得出采样值的虚部,最后通过CORDIC算法圆周旋转模式根据测算出的数据延时电角度对原采样值进行旋转补偿,得到补偿电角度后的采样值;具体步骤如下:
步骤A,利用高压端传感器与互感器合并单元之间的采样值映射关系,实时测量采集数据的延时Td,进而得出数据延时电角度θ;
步骤B,计算采样值幅值的平方值X2;
其中,x为当前时刻的采样值,x1,x2分别为前两个时刻的采样值,Ts为采样周期,ω为角频率,ω=2πf,f为工频等于50Hz,sin2ωTS为定值;
步骤C,计算采样值虚部值y的平方值y2:
y2=X2-x2
步骤D,计算采样值虚部y:
将(y2+1/4,y2-1/4)代入CORDIC算法双曲向量模式中进行迭代,得到采样值虚部y;
步骤E,将(x,y)代入CORDIC算法圆周旋转模式中进行迭代,得经延时补偿θ电角度后的坐标(x',y'),从而得到补偿电角度后的采样值。
2.根据权利要求1所述的一种互感器采样值延时补偿方法,其特征在于,在步骤A中,所述数据延时电角度θ,其计算公式为:
其中,L'为在Td时间内互感器合并单元晶振振动次数,L为360°电气角度对应的互感器合并单元晶振振动次数。
3.根据权利要求1所述的一种互感器采样值延时补偿方法,其特征在于,在步骤D中,所述CORDIC算法双曲向量模式,经n次迭代后坐标关系为:
式中,Sn=-sgn(Yn),CORDIC算法双曲向量模式比例因子在迭代次数n确定后K′n为常数;预先对迭代初始值(x0,y0)进行比例修正,乘以比例因子K′n得到修正后的初始值(x′0,y′0)再进行迭代计算,ti满足tanh(ti)=2-i,i为自然数;
其中,sgn(Yn)为符号运算,Yn是n次迭代后采样值的虚部,Yn为正实数时,sgn(Yn)值为1;Yn为负实数时,sgn(Yn)值为-1;Yn为0时,sgn(Yn)值为0;
逐次迭代后y趋于0,最终收敛于x轴上的向量(r,0),有
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310276475.9A CN103293363B (zh) | 2013-07-02 | 2013-07-02 | 一种互感器采样值延时补偿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310276475.9A CN103293363B (zh) | 2013-07-02 | 2013-07-02 | 一种互感器采样值延时补偿方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103293363A CN103293363A (zh) | 2013-09-11 |
CN103293363B true CN103293363B (zh) | 2015-10-07 |
Family
ID=49094580
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310276475.9A Expired - Fee Related CN103293363B (zh) | 2013-07-02 | 2013-07-02 | 一种互感器采样值延时补偿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103293363B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103983828B (zh) * | 2014-05-26 | 2016-09-07 | 国家电网公司 | 一种电子式互感器数字化相位补偿方法 |
CN107942280A (zh) * | 2018-01-03 | 2018-04-20 | 中国电力科学研究院有限公司 | 一种用于对绝对延迟时间进行校准的方法及*** |
CN110672974A (zh) * | 2019-09-27 | 2020-01-10 | 长沙理工大学 | 一种新型配电网电压跌落检测方法 |
CN111261374B (zh) * | 2020-03-26 | 2021-08-27 | 河北中开明泰电气设备有限公司 | 一种环保型智能互感器及互感器补偿方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6483501B1 (en) * | 1999-08-18 | 2002-11-19 | Tai-Her Yang | Load current compensated type voltage sampling circuit |
CN101556318A (zh) * | 2009-04-20 | 2009-10-14 | 江苏省电力试验研究院有限公司 | 自适应合并器延时同步的电子式互感器稳态校验装置 |
CN101738495A (zh) * | 2009-12-18 | 2010-06-16 | 浙江大学 | 基于cordic算法的电容式微加速度计信号检测装置 |
CN102914695A (zh) * | 2012-10-11 | 2013-02-06 | 南京亚派科技实业有限公司 | 具有实时延时补偿功能的谐波检测方法 |
-
2013
- 2013-07-02 CN CN201310276475.9A patent/CN103293363B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6483501B1 (en) * | 1999-08-18 | 2002-11-19 | Tai-Her Yang | Load current compensated type voltage sampling circuit |
CN101556318A (zh) * | 2009-04-20 | 2009-10-14 | 江苏省电力试验研究院有限公司 | 自适应合并器延时同步的电子式互感器稳态校验装置 |
CN101738495A (zh) * | 2009-12-18 | 2010-06-16 | 浙江大学 | 基于cordic算法的电容式微加速度计信号检测装置 |
CN102914695A (zh) * | 2012-10-11 | 2013-02-06 | 南京亚派科技实业有限公司 | 具有实时延时补偿功能的谐波检测方法 |
Non-Patent Citations (1)
Title |
---|
电子式互感器采样***固有延时测试研究;张晓华等;《电测与仪表》;20110731;第28卷(第547期);42-45 * |
Also Published As
Publication number | Publication date |
---|---|
CN103293363A (zh) | 2013-09-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105738696B (zh) | 全相位时移相位差频率估计方法及装置 | |
CN108535679B (zh) | 基于插值法的电容式电压互感器谐波电压精准测量方法 | |
US8198904B2 (en) | Synchrophasor measuring device and inter-bus-line phase angle difference measurement unit using the same | |
CN103293363B (zh) | 一种互感器采样值延时补偿方法 | |
CN105652151B (zh) | 基于线路参数检测与数据非同步校验的双端测距方法 | |
CN104038134B (zh) | 一种基于线性霍尔的永磁同步电机转子位置误差校正方法 | |
CN109633262A (zh) | 基于组合窗多谱线fft的三相谐波电能计量方法、装置 | |
CN107632236B (zh) | 一种基于对端母线反射波识别的单出线输电线路单端故障测距方法 | |
CN103399203A (zh) | 一种基于复合迭代算法的谐波参数高精度估计方法 | |
CN102435844A (zh) | 一种频率无关的正弦信号相量计算方法 | |
CN106777836A (zh) | 一种变压器直流偏磁仿真模拟方法及装置 | |
CN116359605B (zh) | 一种基于二次加权的谐波信号分析方法 | |
CN108562820A (zh) | 小电流接地***单相接地故障选线方法及装置 | |
CN110967658B (zh) | 一种基于数字微差法的模拟量输入合并单元校验仪溯源的方法 | |
CN103217578A (zh) | 基于锁相环技术测量信号间相位差的数字相位计及其方法 | |
CN105510706A (zh) | 一种高精度欠采样测频方法 | |
CN111308198B (zh) | 一种基于Hanning窗的加窗插值DFT的谐波测量方法 | |
CN104360156A (zh) | 一种电力***频率偏移时信号相角的测量方法 | |
CN114280526B (zh) | 一种电子式互感器校验仪数字微差溯源***及方法 | |
CN107037324B (zh) | 基于单端电气量的不受过渡电阻影响的故障测距方法 | |
CN103383443A (zh) | 智能电网用模块化三相仪表检定装置 | |
CN105067882A (zh) | 一种电量幅值的测定方法 | |
CN104330079A (zh) | 一种多陀螺仪的角速度测量方法及*** | |
CN108982954B (zh) | 适用于馈线终端的计算相电压幅值与相位的方法及*** | |
CN104135284B (zh) | 一种鉴相方法、装置,以及锁相方法和锁相环 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20151007 Termination date: 20180702 |
|
CF01 | Termination of patent right due to non-payment of annual fee |