CN105957703B - 电子式互感器的数字积分方法和数字积分器 - Google Patents
电子式互感器的数字积分方法和数字积分器 Download PDFInfo
- Publication number
- CN105957703B CN105957703B CN201610430451.8A CN201610430451A CN105957703B CN 105957703 B CN105957703 B CN 105957703B CN 201610430451 A CN201610430451 A CN 201610430451A CN 105957703 B CN105957703 B CN 105957703B
- Authority
- CN
- China
- Prior art keywords
- digital
- mutual inductor
- electronic mutual
- integrator
- present
- 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.)
- Active
Links
- 230000010354 integration Effects 0.000 title claims abstract description 21
- 238000000034 method Methods 0.000 title claims abstract description 16
- 238000012546 transfer Methods 0.000 claims abstract description 8
- 230000005540 biological transmission Effects 0.000 abstract description 13
- 238000005070 sampling Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000005350 ferromagnetic resonance Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000005291 magnetic effect Effects 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01F—MAGNETS; INDUCTANCES; TRANSFORMERS; SELECTION OF MATERIALS FOR THEIR MAGNETIC PROPERTIES
- H01F38/00—Adaptations of transformers or inductances for specific applications or functions
- H01F38/20—Instruments transformers
- H01F38/22—Instruments transformers for single phase ac
- H01F38/28—Current transformers
- H01F38/32—Circuit arrangements
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
本发明涉及一种电子式互感器数字积分方法和数字积分器,***传递函数为y(n)=A*[a*[x(n‑2)‑x(n)]+x(n‑1)+x(n)]‑B*y(n‑1)本发明提供的电子式互感器数字积分方法和积分器采用双线性算法,不易饱和;其中参数A和B的值由实际电子式互感器元器件的值决定,使本发明适用于各种电子式互感器;并且本发明的传递函数还设有用于相位补偿的a,使本发明不仅能进行积分运算,还能进行相位补偿,使得到的结果更加准确。
Description
技术领域
本发明属于电子式互感器的数据采集与处理领域,具体涉及一种用于电子式互感器的数字积分方法,即数字积分器的设计。
背景技术
随着电力***容量和电压等级的逐步提高,传统的电磁式电流互感器因体积大、重量重、动态范围窄。容易产生铁磁谐振,输出量为模拟量的缺点,无法满足当前智能电网的要求。电子式互感器绝缘性好,抗干扰能力强,不存在铁磁谐振和磁饱和、动态范围大等优点,适应智能电网的发展,在数字化变电站中得到了广泛应用。
罗氏线圈原理的电流互感器和阻容分压原理的电压互感器输出是对一次信号进行微分的信号,因此,要恢复与一次信号成正比例的信号就必须有相应的积分环节。
目前常用的方法包括模拟积分电路和数字积分器,模拟积分电路的方式,由于电子元器件不是理想器件,运放及其他模拟器件的时漂和温漂等因素都会对积分结果产生影响,造成积分误差,对电子互感器一次侧的信号还原不准确。
《电子测量技术》2014年9月第37卷第9期“Rogowski线圈数字积分器的分析与设计”披露了一种数字积分器。积分器的离散方程为:
y(n)=y(n-1)+1/2[x(n)+x(n-1)]T
但是这种数字积分器采用的是梯形算法,容易饱和,不能进行相位补偿,并且适用性差,只适用于特定的电子式互感器,当使用到不同的数字积分器时,会产生相应的误差。
发明内容
本发明提供了一种电子式互感器的数字积分方法和数字积分器,用于解决数字积分器适用性差的问题。
一种电子式互感器数字积分方法,包括通过以下***传递函数:y(n)=A*[a*[x(n-2)-x(n)]+x(n-1)+x(n)]-B*y(n-1)进行积分运算的步骤;
其中n为在二次侧采集信号的序号,x(n)为在二次侧采集到的信号,y(n)为积分后的数据;a在可调范围[0,Ts]上取值;A=α*Ts/R*C;B=β*Rf*C*Ts–1;
Ts、R、C是电子式互感器的参数,R为负载电阻,C为杂散电容,Rf是反馈电阻,Ts为一个采样周期,α和β为设定系数。
所述α=0.498,β=400。
一种电子式互感器数字积分器,包括通过以下***传递函数:y(n)=A*[a*[x(n-2)-x(n)]+x(n-1)+x(n)]-B*y(n-1)进行积分运算的模块;
其中n为在二次侧采集信号的序号,x(n)为在二次侧采集到的信号,y(n)为积分后的数据;a在可调范围[0,TS]上取值;A=α*Ts/R*C;B=β*Rf*C*Ts–1;
Ts、R、C是电子式互感器的参数,R为负载电阻,C为杂散电容,Rf是反馈电阻,Ts为一个采样周期,α和β为设定系数。
所述α=0.498,β=400。
本发明提供的方案采用的是双线性算法,得到的传递函数不易饱和;在传递函数中设有用于积分回路相位补偿的参数a,可以进行相位补偿,使得道的结果更加准确;并且传递函数中的参数A和B都是根据电子式互感器的实际参数进行变化的,当采用不同的电子式互感器时,只需根据电子式互感器中元器件的参数改变传递函数中A、B的值即可。
附图说明
图1为微分与积分的等效电路图;
图2为罗氏线圈的等效电路图;
图3为Simulink仿真模型;
图4为基波输入输出波形的对比;
图5为***幅频和相频特性曲线。
具体实施方式
本发明提供的电子式互感器数字积分方法,其电子式互感器的模拟微分和积分原型如图1所示,虚线前段是电子式互感器微分模型,虚线后端是模拟积分的模型,此处电子式互感器微分模型采用的是罗氏线圈。
图1中,电子式互感器的微分电路包括线圈的电感L,杂散电容C和负载电阻R,积分回路连续域的传递函数为:
将上述传递函数通过变换转换为离散域,得到积分回路在离散域的传递函数如下:
NUM=[A A]
DEN=[1 -B]
A=0.498*Ts/R*C;
B=400*Rf*C*Ts–1;
其中Ts为采样周期,R是连接放大器负载电阻,Rf是反馈电阻,Rf/R决定着***的增益系数,Rf*C决定***与电子式互感器的特性匹配。
在实际应用中,因为理想积分存在***饱和的问题,所以无论是模拟积分还是数字积分都不可能实现理想的积分特性,最终的输出结果存在相位差,因此需要对相位进行补偿。
采用二阶插值算法相位补偿,得到相位补偿传递函数
z(n)=a*x(n)+(1-a)*x(n-1)
其中0<a<1;
将相位补偿传递函数与积分传递函数进行卷积运算,得到最终的***传递函数:
y(n)=A*[a*[x(n-2)-x(n)]+x(n-1)+x(n)]-B*y(n-1)
其中0<a<1,0.001<A<0.1,0.9<B<1;
其中参数A,B的值由一次电子式互感器的参数决定,参数a用于积分回路相位补偿,可调范围为Ts,分辨率为1′。
以上为本发明的理论推导过程,本发明的方法为:通过以下***传递函数:y(n)=A*[a*[x(n-2)-x(n)]+x(n-1)+x(n)]-B*y(n-1)进行积分运算;
其中n为在二次侧采集信号的序号,x(n)为在二次侧采集到的信号,y(n)为积分后的数据;a在可调范围[0,Ts]上取值;A=α*Ts/R*C;B=β*Rf*C*Ts–1;
Ts、R、C是电子式互感器的参数,R为负载电阻,C为杂散电容,Rf是反馈电阻,Ts为一个采样周期,α和β为设定系数。
下面对本发明的方法进行验证。
如图2所示的电路图为220kV罗氏线圈的等效电路,其传递函数经过拉普拉斯变换后为:
其中,M为线圈的互感系数。
将参数代入,并经过双线性变化后,得到传递函数:
NUM=[-3e-3 0]
DEN=[1.875e-9 11.15e-6 1.006]
经过参数匹配得到积分算法的参数A=0.041654,B=0.99938,a=0.02;将电子式互感器和积分算法组成的模型通过图3的MATLAB内部simulink工具仿真。其基波的输入输出比较如图4所示。整个***的幅频和相频特性如图5所示。通过波形我们可以看出整个***能完整的还原波形,并且具有优秀的幅频和相频曲线。
在本方法实施例中,设定系数α=0.498,β=400,作为本发明的其他实施例,设定参数α可以采用其它与0.498相近的数值,β可以采用与400相近的数值。
依照上述方法,可以得到一种电子式互感器数字积分器,该积分器的***传递函数为:y(n)=A*[a*[x(n-2)-x(n)]+x(n-1)+x(n)]-B*y(n-1),参数解释等与上述方法相同。该积分器并不是硬件积分器,而是一种软件实现的数字积分器,依照上述方法进行编程即可实现,程序运行于采集设备中。
以上给出了本发明涉及的具体实施方式,但本发明不局限于所描述的实施方式。在本发明给出的思路下,采用对本领域技术人员而言容易想到的方式对上述实施例中的技术手段进行变换、替换、修改,并且起到的作用与本发明中的相应技术手段基本相同、实现的发明目的也基本相同,这样形成的技术方案是对上述实施例进行微调形成的,这种技术方案仍落入本发明的保护范围内。
Claims (4)
1.一种电子式互感器数字积分方法,其特征在于,包括通过以下***传递函数:y(n)=A*[a*[x(n-2)-x(n)]+x(n-1)+x(n)]-B*y(n-1)进行积分运算的步骤;
其中n为在二次侧采集信号的序号,x(n)为在二次侧采集到的信号,y(n)为积分后的数据;a在可调范围[0,Ts]上取值;A=α*Ts/R*C;B=β*Rf*C*Ts–1;
Ts、R、C是电子式互感器的参数,R为负载电阻,C为杂散电容,Rf是反馈电阻,Ts为一个采样周期,α和β为设定系数。
2.根据权利要求1所述的电子式互感器数字积分方法,其特征在于,所述α=0.498,β=400。
3.一种电子式互感器数字积分器,其特征在于,包括通过以下***传递函数:y(n)=A*[a*[x(n-2)-x(n)]+x(n-1)+x(n)]-B*y(n-1)进行积分运算的模块;
其中n为在二次侧采集信号的序号,x(n)为在二次侧采集到的信号,y(n)为积分后的数据;a在可调范围[0,Ts]上取值;A=α*Ts/R*C;B=β*Rf*C*Ts–1;
Ts、R、C是电子式互感器的参数,R为负载电阻,C为杂散电容,Rf是反馈电阻,Ts为一个采样周期,α和β为设定系数。
4.根据权利要求3所述的电子式互感器数字积分器,其特征在于,所述α=0.498,β=400。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610430451.8A CN105957703B (zh) | 2016-06-16 | 2016-06-16 | 电子式互感器的数字积分方法和数字积分器 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610430451.8A CN105957703B (zh) | 2016-06-16 | 2016-06-16 | 电子式互感器的数字积分方法和数字积分器 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105957703A CN105957703A (zh) | 2016-09-21 |
CN105957703B true CN105957703B (zh) | 2017-09-22 |
Family
ID=56906767
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610430451.8A Active CN105957703B (zh) | 2016-06-16 | 2016-06-16 | 电子式互感器的数字积分方法和数字积分器 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105957703B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112379323B (zh) * | 2020-11-05 | 2022-10-21 | 山东山大电力技术股份有限公司 | 一种基于罗氏线圈特性的等效模拟发生器及电子设备 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5327133A (en) * | 1993-02-16 | 1994-07-05 | Motorola, Inc. | Digital integrator with reduced circuit area and analog-to-digital converter using same |
CN1821794A (zh) * | 2006-03-31 | 2006-08-23 | 北京万工科技有限公司 | 用于洛高夫斯基(Rogowski)线圈的积分器及其实现方法 |
CN102128966A (zh) * | 2011-01-11 | 2011-07-20 | 沈阳互感器有限责任公司 | 基于罗氏线圈的电子式电流互感器用积分器电路 |
WO2013177002A1 (en) * | 2012-05-23 | 2013-11-28 | Finisar Corporation | Low power and compact area digital integrator for a digital phase detector |
CN104504272A (zh) * | 2014-12-26 | 2015-04-08 | 常州大学 | 基于理想重构滤波器冲激响应拟合的数字积分器实现方法 |
CN104503294A (zh) * | 2014-11-21 | 2015-04-08 | 三峡大学 | 一种基于直流负反馈原理的梯形-矩形加权数字积分器 |
CN104868900A (zh) * | 2015-05-25 | 2015-08-26 | 三峡大学 | 一种基于龙贝格算法的高精度数字积分器 |
-
2016
- 2016-06-16 CN CN201610430451.8A patent/CN105957703B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5327133A (en) * | 1993-02-16 | 1994-07-05 | Motorola, Inc. | Digital integrator with reduced circuit area and analog-to-digital converter using same |
CN1821794A (zh) * | 2006-03-31 | 2006-08-23 | 北京万工科技有限公司 | 用于洛高夫斯基(Rogowski)线圈的积分器及其实现方法 |
CN102128966A (zh) * | 2011-01-11 | 2011-07-20 | 沈阳互感器有限责任公司 | 基于罗氏线圈的电子式电流互感器用积分器电路 |
WO2013177002A1 (en) * | 2012-05-23 | 2013-11-28 | Finisar Corporation | Low power and compact area digital integrator for a digital phase detector |
CN104503294A (zh) * | 2014-11-21 | 2015-04-08 | 三峡大学 | 一种基于直流负反馈原理的梯形-矩形加权数字积分器 |
CN104504272A (zh) * | 2014-12-26 | 2015-04-08 | 常州大学 | 基于理想重构滤波器冲激响应拟合的数字积分器实现方法 |
CN104868900A (zh) * | 2015-05-25 | 2015-08-26 | 三峡大学 | 一种基于龙贝格算法的高精度数字积分器 |
Also Published As
Publication number | Publication date |
---|---|
CN105957703A (zh) | 2016-09-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Gobet et al. | Noise analysis of switched capacitor networks | |
CN102301602B (zh) | 用于使信号数字化的装置和方法 | |
DE112016006971T5 (de) | Verfahren und Vorrichtung zum Herausfinden des Phasenfehlers oder von Timing-Verzögerungen in einem Stromwandler und Leistungsmessvorrichtung einschließlich Stromwandlerfehlerkorrektur | |
CN100439926C (zh) | 用于洛高夫斯基(Rogowski)线圈的积分器及其实现方法 | |
CN106093547B (zh) | 基于辅助线圈修正的罗氏线圈电流互感器测量装置及方法 | |
CN103257271A (zh) | 一种基于stm32f107vct6的微电网谐波与间谐波检测装置及检测方法 | |
CN108710774B (zh) | 一种电流传感器频率特性建模方法 | |
CN101846699A (zh) | 电参数测量装置、***及方法 | |
CN105487034A (zh) | 一种0.05级电子式互感器校验方法及*** | |
CN103487648B (zh) | sigma-delta PLL频率测量电路及方法 | |
CN107478896A (zh) | 一种基于级联广义积分器的频率自适应谐波电流检测方法 | |
CN106546798A (zh) | 一种模拟积分器及电子电流互感器采集电路 | |
CN106291092A (zh) | 一种智能电表防外部恒磁场干扰方法 | |
CN114035142A (zh) | 基于分段拉格朗日插值的电能表误差补偿方法及*** | |
CN105957703B (zh) | 电子式互感器的数字积分方法和数字积分器 | |
Wang et al. | Nonlinear behavior immunity modeling of an LDO voltage regulator under conducted EMI | |
CN109546660A (zh) | 基于神经滑模控制策略的有源电力滤波电路和控制方法 | |
CN204439713U (zh) | 基于新型直流隔离积分器的空心线圈电流互感器 | |
CN206594219U (zh) | 一种基于罗氏线圈暂态特性的采集电路及模拟积分器 | |
CN106685376A (zh) | 一种分数阶电感电路 | |
WO2023124592A1 (zh) | 基于逆黑盒及逆电磁对偶模型的pt一次电压重构方法 | |
CN205725711U (zh) | 用于鼠笼式空芯电流互感器的高集成化数字积分器 | |
CN101872006B (zh) | 带升压器的电子式周期性非正弦波基准电压互感器 | |
CN113190949B (zh) | 一种cvt宽频仿真模型的时域构建方法 | |
CN104504272A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |