CN111817713A - 可快速同步对称故障下电压相位的高压直流锁相环及算法 - Google Patents
可快速同步对称故障下电压相位的高压直流锁相环及算法 Download PDFInfo
- Publication number
- CN111817713A CN111817713A CN202010516639.0A CN202010516639A CN111817713A CN 111817713 A CN111817713 A CN 111817713A CN 202010516639 A CN202010516639 A CN 202010516639A CN 111817713 A CN111817713 A CN 111817713A
- Authority
- CN
- China
- Prior art keywords
- phase
- voltage
- sampling
- fault
- current
- 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.)
- Granted
Links
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 33
- 238000006243 chemical reaction Methods 0.000 claims abstract description 22
- 238000001514 detection method Methods 0.000 claims abstract description 14
- 230000005540 biological transmission Effects 0.000 claims abstract description 11
- 230000008859 change Effects 0.000 claims abstract description 11
- 230000000737 periodic effect Effects 0.000 claims abstract description 3
- 238000005070 sampling Methods 0.000 claims description 90
- 238000000034 method Methods 0.000 claims description 30
- 230000010354 integration Effects 0.000 claims description 19
- 238000012546 transfer Methods 0.000 claims description 16
- 230000008569 process Effects 0.000 claims description 15
- 238000005259 measurement Methods 0.000 claims description 14
- 238000012545 processing Methods 0.000 claims description 14
- 230000009466 transformation Effects 0.000 claims description 12
- 238000001914 filtration Methods 0.000 claims description 11
- 230000001360 synchronised effect Effects 0.000 claims description 10
- 230000009471 action Effects 0.000 claims description 9
- 238000009499 grossing Methods 0.000 claims description 9
- 230000009191 jumping Effects 0.000 claims description 7
- 230000002194 synthesizing effect Effects 0.000 claims description 3
- 230000004044 response Effects 0.000 description 5
- 238000013178 mathematical model Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000013016 damping Methods 0.000 description 2
- 238000000113 differential scanning calorimetry Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- CNKHSLKYRMDDNQ-UHFFFAOYSA-N halofenozide Chemical compound C=1C=CC=CC=1C(=O)N(C(C)(C)C)NC(=O)C1=CC=C(Cl)C=C1 CNKHSLKYRMDDNQ-UHFFFAOYSA-N 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03L—AUTOMATIC CONTROL, STARTING, SYNCHRONISATION OR STABILISATION OF GENERATORS OF ELECTRONIC OSCILLATIONS OR PULSES
- H03L7/00—Automatic control of frequency or phase; Synchronisation
- H03L7/06—Automatic control of frequency or phase; Synchronisation using a reference signal applied to a frequency- or phase-locked loop
- H03L7/08—Details of the phase-locked loop
- H03L7/085—Details of the phase-locked loop concerning mainly the frequency- or phase-detection arrangement including the filtering or amplification of its output signal
- H03L7/093—Details of the phase-locked loop concerning mainly the frequency- or phase-detection arrangement including the filtering or amplification of its output signal using special filtering or amplification characteristics in the loop
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R29/00—Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
- G01R29/16—Measuring asymmetry of polyphase networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/24—Arrangements for preventing or reducing oscillations of power in networks
- H02J3/241—The oscillation concerning frequency
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/36—Arrangements for transfer of electric power between ac networks via a high-tension dc link
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03L—AUTOMATIC CONTROL, STARTING, SYNCHRONISATION OR STABILISATION OF GENERATORS OF ELECTRONIC OSCILLATIONS OR PULSES
- H03L7/00—Automatic control of frequency or phase; Synchronisation
- H03L7/06—Automatic control of frequency or phase; Synchronisation using a reference signal applied to a frequency- or phase-locked loop
- H03L7/08—Details of the phase-locked loop
- H03L7/085—Details of the phase-locked loop concerning mainly the frequency- or phase-detection arrangement including the filtering or amplification of its output signal
- H03L7/091—Details of the phase-locked loop concerning mainly the frequency- or phase-detection arrangement including the filtering or amplification of its output signal the phase or frequency detector using a sampling device
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
- Y02E60/60—Arrangements for transfer of electric power between AC networks or generators via a high voltage DC link [HVCD]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Inverter Devices (AREA)
- Stabilization Of Oscillater, Synchronisation, Frequency Synthesizers (AREA)
Abstract
本发明公开了一种可快速同步对称故障下电压相位的高压直流锁相环及算法,算法包括首先进行故障检测:在锁相环的测量系统中,根据高压直流输电系统的换流母线前三次的相位增量差、频率跟踪器的输出变化率和输入的换相电压幅值,以及换流母线当前a、b和c相电压的周期检测结果,检测是否发生对称故障;然后根据检测结果选定锁相方式;在发生对称故障的情况下计算相位跳变补偿值,并将对称故障下的相位跳变补偿值输入到锁相环的锁相控制系统中;最后按照选定的锁相方式对高压直流换相电压的相位进行同步。本发明解决了高压直流的换流母线发生对称故障时的相位同步问题。
Description
技术领域
本发明涉及高压直流输电技术领域,特别涉及一种可快速同步对称故障下电压相位的高压直流锁相环及算法。
背景技术
同步触发控制是整个直流输电控制系统的基础,通过改变换流阀触发脉冲的产生时刻或相位,实现对整个直流输电系统的调节,一般依赖于锁相环(PLL)对换相电压的相位进行同步。目前,我国高压直流输电工程中的控制保护系统都是基于高性能的控制平台,例如ABB公司的MACH2,SIEMENS公司的SIMADYN-D、SIMATIC-TDC。两家公司控制保护系统采用的锁相环有所不同。ABB-PLL本质上是一个静止坐标系锁相环,为了准确追踪交流测母线电压相位,通常将锁相环的带宽设置的很小。SIEMENS-PLL则是通过一个MXF128数字滤波器实现对基波电压的提取,利用反正切计算得到交流电压相角,最后利用PI控制器进行调节,使锁相环输出的角度值始终跟踪电网的实际相位,同时,附加的频率跟踪器动态调整采样周期。
实际上,当高压直流输电系统的换流母线发生故障时,换流母线的电压的幅值和相位都会发生跳变。ABB-PLL受限于其小带宽的锁相环调节器设置,在故障发生时刻及其后的0.1s内通常无法准确地跟踪换流母线电压的相位情况,使得触发控制系统无法及时地根据触发角指令产生触发脉冲,高压直流输电系统可能会因此而发生后续换相失败,且故障后的系统恢复性能受到影响。大多数情况下(特别是在故障工况下),实际电网中的频率都是恒定不变的,而SIEMENS-PLL所附加的频率跟踪器在故障工况下依然会检测到频率的偏移情况,这使得锁相环恢复零误差锁相的过程更为复杂。
因此,需要提供一种可以保持较大带宽和在故障下及时恒定电网频率使得频率与相位解耦的高压直流锁相环。
发明内容
本发明的第一目的在于克服现有技术的缺点与不足,提供一种可快速同步对称故障下电压相位的高压直流锁相环算法,该算法可以解决高压直流的换流母线发生对称故障时的相位同步问题。
本发明的第二目的在于提供一种可快速同步对称故障下电压相位的高压直流锁相环,该锁相环可以实现与对称故障下电压相位的快速同步。
本发明的第一目的通过下述技术方案实现:一种可快速同步对称故障下电压相位的高压直流锁相环算法,步骤如下:
S1、故障检测:在锁相环的测量系统中,根据高压直流输电系统的换流母线前三次的相位增量差、频率跟踪器的输出变化率和输入的换相电压幅值,以及换流母线当前a相电压va、b相电压vb和c相电压vc的周期检测结果,检测是否发生对称故障;
S2、选定锁相方式:若判定为对称故障,则将锁相环运行在第二锁相方式;若判定为额定工况或者非对称故障,则将锁相环运行在第一锁相方式;
S3、在发生对称故障的情况下计算相位跳变补偿值,并将对称故障下的相位跳变补偿值输入到锁相环的锁相控制系统中;
S4、按照步骤S2中选定的锁相方式对高压直流换相电压的相位进行同步。
优选的,步骤S1的故障检测过程具体如下:
S11、对于当前时刻及其前三个时刻,将这些时刻对应的相位增量输出按照时间序列依次进行比较,当相位增量差连续3次大于第一阈值,则生成第一控制信号;
将频率跟踪器的输出变化率与第二阈值进行比较,在输出变化率大于阈值的情况下生成第二控制信号;
获取换流母线当前a相电压va、b相电压vb和c相电压vc,并将这些电压分量在αβ坐标系下的电压幅值与第三阈值比较,当αβ坐标系下的电压幅值大于第三阈值时,生成第三控制信号;
S12、综合第一、二和第三控制信号,生成故障控制信号;
S13、将换流母线当前a相电压va、b相电压vb和c相电压vc分别与各自对应的前一周期的电压值进行比较,小于或等于则表示换流母线正常工作,处于额定工况;
大于则表示响应检测,其中,最先响应检测的相即为故障相,若三个相同时响应,则为对称故障,若只有最先响应的相响应,则为非对称故障。
更进一步的,第一阈值为0.3°,第二阈值为15000,第三阈值为0.9p.u.。
优选的,步骤S3中相位跳变补偿值的计算过程如下:
S31、将换流母线当前a相电压va、b相电压vb和c相电压vc转换至αβ坐标系下;
更进一步的,GDSC输入与输出的关系表示为:
式中,H(k)为单位阶跃函数;为相位跳变前的输入信号;为相位跳变后的信号;为相位跳变前的初相位;即为对称故障引起的相位跳变量;Ts为采样周期;ω为电网角频率;Vpre(k)为跳变前电压信号的幅值;为跳变前相角的指数表示;Vpost(k)为跳变后电压信号的幅值;为跳变后相角的指数表示;
最终,GDSC输入与输出的关系式为:
式中,γ(k)=[1+floor(32k/N)]/32,floor表示向下取整,那么△γ(k)=1/32;
步骤S32生成的比值向量为:
更进一步的,直流GDSC0输入与输出的关系表示为:
优选的,在步骤S4中,在额定工况或者非对称故障工况下,按照选定的第一锁相方式对高压直流换相电压的相位进行同步,过程如下:
S411、将上述的换流母线当前a相电压va、b相电压vb和c相电压vc,以采样周期Ts为采样基准,进行插补采样处理;
利用一个固定储存位置来储存采样周期Ts的累加时间进度值t0=kTs;
S412、将上述t0与运行时间t进行比较,运行时间t为运行步长tstep的累加时间,当t≥t0时,输出同步信号sync=1,作为使能信号,同时输出线性插值的百分数frac=(t-t0)/tstep;
S413、考虑到采样周期Ts不一定为运行步长的整数倍,需要对采样信号进行自定义插值处理:
Out′i=Outi-(Outi-Outi-1)×frac
=Outi×frac+Outi-1×(1-frac)
其中,Out′i表示最终的插值处理所得的信号结果;Outi表示最近一次的获取的电压信号;Outi-1表示上一次获取的电压信号;
插值采样后的获取的三相电压信号为:
S414、将上述插补采样后的三相电压信号转换到αβ坐标系下,然后进行Clarke变换:
Clarke变换后获得的电压信号为:
S415、利用FIR数字滤波器从上述电压信号中提取出αβ坐标系下的基频正序电压分量:
S417、将θ同时输入到频率跟踪器和锁相控制系统;
在频率跟踪器中:先将第k次的反正切运算结果与第(k-1)次的结果作差,并将得到的差值乘以N/4,再与π/2作差,经过频率跟踪器的调节器获取额定频率的变换量;
在锁相控制系统中:先将反正切运算结果θ与输出的相位进行比较,
优选的,在步骤S4中,在对称故障工况下,按照选定的第二锁相方式对高压直流换相电压的相位进行同步,过程如下:
S421、将换流母线当前a相电压va、b相电压vb和c相电压vc,以采样周期Ts为采样基准,进行插补采样处理;
利用一个固定储存位置来储存采样周期Ts的累加时间进度值t0=kTs;
S422、将上述t0与运行时间t进行比较,运行时间t为运行步长tstep的累加时间,当t≥t0时,输出同步信号sync=1,作为使能信号,同时输出线性插值的百分数frac=(t-t0)/tstep;
S423、考虑到采样周期Ts不一定为运行步长的整数倍,需要对采样信号进行自定义插值处理:
Out′i=Outi-(Outi-Outi-1)×frac
=Outi×frac+Outi-1×(1-frac)
其中,Out′i表示最终的插值处理所得的信号结果;Outi表示最近一次的获取的电压信号;Outi-1表示上一次获取的电压信号;
插值采样后的获取的三相电压信号为:
S424、将上述插补采样后的三相电压信号转换到αβ坐标系下,然后进行Clarke变换:
Clarke变换后获得的电压信号为:
S425、利用FIR数字滤波器从上述电压信号中提取出αβ坐标系下的基频正序电压分量:
其中,N对应于一个电网周期的采样次数,X(k-j)对应最新输入的j点数据,j为正整数,Y(k)为对应的输出,将所得的vα(k)和vβ(k)作为最新的输入数据X(k-j),所得的输出Y(k)即为在αβ坐标系下的基频正序电压分量和
S427、将θ同时输入到频率跟踪器和锁相控制系统;
在频率跟踪器中:故障控制信号动作前一个周期的估算频率是固定的,故保留故障控制信号动作前一个周期的输出再作取倒数和折算后可获得最新的采样周期Ts,Ts可继续应用于测量系统的三相电压的获取及在各模块的运行中进行采样;
优选的,在第一锁相方式的运行情况下,测量系统在离散域的开环传递函数为:
其中,z为开环传递函数的变量;ω为电网电压角频率;τ为平滑滤波时间常数;Ti为积分时间常数;
由该开环传递函数即可整理得到闭环传递函数的特征方程:
为保证测量系统稳定性,平滑滤波时间常数和积分时间常数的限制范围为:
在第一锁相方式的运行情况下,锁相控制系统的闭环传递函数为:
其中,ki为调节器的积分系数;kp为调节器的比例系数;T为第一积分环节的采样时间常数,T与运行步长相等;
为保证测量系统稳定性,比例系数和积分系数的可行域为:
本发明的第二目的通过下述技术方案实现:一种可快速同步对称故障下电压相位的高压直流锁相环,所述高压直流锁相环通过本发明第一目的所述的可快速同步对称故障下电压相位的高压直流锁相环算法,实现与对称故障下高压直流换相电压相位的快速同步。
本发明相对于现有技术具有如下的优点及效果:
(1)针对高压直流输电系统中换流母线发生对称性故障的情况,本发明高压直流锁相环算法通过计算相位跳变量来及时补偿响应,使得锁相环仍然能够快速准确地同步换流母线的电压相位,进而使得同步触发锁相控制系统尽可能地精准触发,有利于抑制逆变器发生后续换相失败。
(2)本发明高压直流锁相环算法能够基于检测结果选择对应的锁相方式,并进行对应的电压相位同步,在正常工作、非对称故障和对称故障工况下都能解决高压直流的换流母线的电压相位同步问题。
(3)本发明保留有SIEMENS锁相环的FIR数字滤波器功能模块,也即是保留了SIEMENS锁相环在严重谐波环境下的滤波功能,如此可以合理地整定锁相部分的调节器的参数,从而实现高带宽。
附图说明
图1是本发明可快速同步对称故障下电压相位的高压直流锁相环算法的流程图。
图2为应用图1算法的锁相环的整体结构框图。
图3(a)和图3(b)为图1算法中锁相环故障检测部分的原理图。
图4为频率跟踪器内部结构框图。
图5为第一锁相方式的数学模型。
图6(a)为第一锁相方式中频率跟踪器的参数整定范围。
图6(b)为第一锁相方式中调节器的参数整定范围。
图7为第二锁相方式的数学模型。
图8为图2锁相环与工程用SIEMENS锁相环在对称故障下的相位误差对比图。
具体实施方式
下面结合实施例及附图对本发明作进一步详细的描述,但本发明的实施方式不限于此。
实施例
本实施例公开了一种可快速同步对称故障下电压相位的高压直流锁相环算法,所述算法应用于高压直流锁相环,如图2所示,所述锁相环具有测量系统和锁相控制系统,锁相环通过该算法实现与对称故障下高压直流换相电压相位的快速同步。算法如图1所示,步骤如下:
S1、故障检测:在锁相环的测量系统中,根据高压直流输电系统的换流母线前三次的相位增量差、频率跟踪器的输出变化率和输入的换相电压幅值,以及换流母线当前a相电压va、b相电压vb和c相电压vc的周期检测结果,检测是否发生对称故障。
如图3(a)和图3(b)所示,故障检测过程具体如下:
S11、对于当前时刻及其前三个时刻,将这些时刻对应的相位增量输出按照时间序列依次进行比较,当相位增量差连续3次大于第一阈值,则生成第一控制信号;
将频率跟踪器的输出变化率与第二阈值进行比较,在输出变化率大于阈值的情况下生成第二控制信号;
获取换流母线当前a相电压va、b相电压vb和c相电压vc,并将这些电压分量在αβ坐标系下的电压幅值与第三阈值比较,当αβ坐标系下的电压幅值大于第三阈值时,生成第三控制信号。
S12、综合第一、二和第三控制信号,生成故障控制信号。
S13、将换流母线当前a相电压va、b相电压vb和c相电压vc分别与各自对应的前一周期的电压值进行比较,小于或等于则表示换流母线正常工作,处于额定工况;
大于则表示响应检测,其中,最先响应检测的相即为故障相,若三个相同时响应,则为对称故障,若只有最先响应的相响应,则为非对称故障。
在本实施例中,如图3(a)所示,第一阈值为0.3°,第二阈值为15000,第三阈值为0.9p.u.。
S2、选定锁相方式:若判定为对称故障,Ctrl信号输出高电平,则将锁相环运行在第二锁相方式;若判定为额定工况或者非对称故障,Ctrl信号输出低电平,则将锁相环运行在第一锁相方式。
S3、在发生对称故障的情况下计算相位跳变补偿值,并将对称故障下的相位跳变补偿值输入到锁相环的锁相控制系统中。
相位跳变补偿值的计算过程具体如下:
S31、将换流母线当前a相电压va、b相电压vb和c相电压vc转换至αβ坐标系下。
S32、基于GDSC(级联消去滤波法)对αβ坐标系下的电压分量进行滤波,获取GDSC的输出,基于该输出与其前一周期的输出,生成对应的比值向量前一周期的输出即是相位跳变前的信号,也就是指故障发生前的信号。
本实施例的级联消去滤波法具体是采用5个DSC(消去滤波法:将当前信号与1/n个周期前的信号相加并作平均,作平均后的信号作为最终的滤波输出;GDSC中,n=2,4,8,16,32),利用其输入与输出的关系进一步求取对称故障下的相位跳变量进行补偿,具有滤波效果好、原理清晰等优点。
式中,k表示第k次采样;为旋转因子,其中,e为自然常数,m为正整数,m的取值由GDSC的DSC滤波器决定,在本实施例中,由于5个级联的DSC滤波器中n最大为32,因此m分别选值为0,1,2,…,31;N为一个周期的采样次数,本实施例N取32;
式中,H(k)为单位阶跃函数;为相位跳变前的输入信号;为相位跳变后的信号;为相位跳变前的初相位;即为对称故障引起的相位跳变量;Ts为采样周期;ω为电网角频率;Vpre(k)为跳变前电压信号的幅值;为跳变前相角的指数表示;Vpost(k)为跳变后电压信号的幅值;为跳变后相角的指数表示;
最终,GDSC输入与输出的关系式为:
式中,γ(k)=[1+floor(32k/N)]/32,floor表示向下取整,那么△γ(k)=1/32;
步骤S32生成的比值向量为:
S4、按照步骤S2中选定的锁相方式对高压直流换相电压的相位进行同步。
具体来说,(1)在额定工况或者非对称故障工况下,按照选定的第一锁相方式对高压直流换相电压的相位进行同步的过程为:
S411、将上述的换流母线当前a相电压va、b相电压vb和c相电压vc,以采样周期Ts为采样基准,进行插补采样处理;这里,取Ts=0.625ms作为启动初始值;
这里是利用采样元件进行采样,其工作原理是设置一个固定储存位置来储存采样周期Ts的累加时间进度值t0=kTs;
S412、将上述t0与运行时间t进行比较,运行时间t为运行步长tstep的累加时间,当t≥t0时,输出同步信号sync=1,作为使能信号,同时输出线性插值的百分数frac=(t-t0)/tstep;使能信号用于对电压信号的采样,当有高电平时,就更新一次电压信号量,其后的值也跟着更新,如图4中,相位值的更新次数与使能脉冲次数一致;
S413、考虑到采样周期Ts不一定为运行步长的整数倍,需要对采样信号进行自定义插值处理:
Out′i=Outi-(Outi-Outi-1)×frac
=Outi×frac+Outi-1×(1-frac)
其中,Out′i表示最终的插值处理所得的信号结果;Outi表示最近一次的获取的电压信号;Outi-1表示上一次获取的电压信号;
插值采样后的获取的三相电压信号为:
S414、将上述插补采样后的三相电压信号转换到αβ坐标系下,然后进行Clarke变换:
Clarke变换后获得的电压信号为:
S415、利用FIR数字滤波器从上述电压信号中提取出αβ坐标系下的基频正序电压分量:
S417、将θ同时输入到频率跟踪器和锁相控制系统;
如图4所示,在频率跟踪器中:先将第k次的反正切运算结果与第(k-1)次的结果作差,并将得到的差值乘以N/4,再与π/2作差,经过频率跟踪器的调节器获取额定频率的变换量;
将上述变换量转为实际值后得到估算频率再作取倒数和折算(图2中的1/x和1/N)处理后获得最新的采样周期Ts,Ts可继续应用于测量系统的三相电压的获取及在各模块的运行中进行采样,采样是通过产生使能信号来实现。
在锁相控制系统中:先将反正切运算结果θ与输出的相位进行比较,
此步骤的频率跟踪器的调节器参数整定主要从稳定性、动态响应和稳定误差情况三方面进行考虑,由图5所示的第一锁相方式的数学模型可以得到,在第一锁相方式的运行情况下,测量系统在离散域的开环传递函数为:
其中,z为开环传递函数的变量;ω为电网电压角频率;τ为平滑滤波时间常数;Ti为积分时间常数;
由该开环传递函数即可整理得到闭环传递函数的特征方程:
为保证测量系统稳定性,平滑滤波时间常数和积分时间常数的限制范围为:
在本实施例中,采样周期稳定在0.625s附近,电网频率恒定为50Hz时,可行域为:
为使测量系统具有较快的响应速度,又避免产生严重的振荡,限制阻尼比介于0.4~0.8,此时闭环特征方程的根均较靠近坐标原点。当电网频率基本稳定在50Hz,采样周期介于0.5625~0.6875ms时,平滑滤波时间常数和积分时间常数的取值区域如图6(b)所示。结合误差传递函数和终值定理进行分析,优选出本实施例的参数:平滑滤波时间常数τ为0.02,积分时间常数Ti为0.05。
其中,ki为调节器的积分系数;kp为调节器的比例系数;T为第一积分环节的采样时间常数,T与运行步长相等,比Ts大得多;
根据劳斯判据,可得到稳定性可行域如图6(a)所示。在本实施例中,将该锁相控制系统的阻尼比选定为0.707,限制Ts在[0.5625ms,0.6875ms]区间取值,为保证测量系统稳定性,比例系数和积分系数的可行域为:
为实现即使在电网频率偏移变化的情况下,锁相控制系统的稳定误差较小,本实施例选定调节器的比例系数kp为1256,积分系数ki为5724,积分环节1的采样时间T为20us。
(2)在对称故障工况下,按照选定的第二锁相方式对高压直流换相电压的相位进行同步的过程为:
S421、将换流母线当前a相电压va、b相电压vb和c相电压vc,以采样周期Ts为采样基准,进行插补采样处理;
利用一个固定储存位置来储存采样周期Ts的累加时间进度值t0=kTs;
S422、将上述t0与运行时间t进行比较,运行时间t为运行步长tstep的累加时间,当t≥t0时,输出同步信号sync=1,作为使能信号,同时输出线性插值的百分数frac=(t-t0)/tstep;
S423、考虑到采样周期Ts不一定为运行步长的整数倍,需要对采样信号进行自定义插值处理:
Out′i=Outi-(Outi-Outi-1)×frac
=Outi×frac+Outi-1×(1-frac)
其中,Out′i表示最终的插值处理所得的信号结果;Outi表示最近一次的获取的电压信号;Outi-1表示上一次获取的电压信号;
插值采样后的获取的三相电压信号为:
S424、将上述插补采样后的三相电压信号转换到αβ坐标系下,然后进行Clarke变换:
Clarke变换后获得的电压信号为:
S425、利用FIR数字滤波器(finite impulse response digital filter)从上述电压信号中提取出αβ坐标系下的基频正序电压分量,FIR数字滤波器的计算公式为:
其中,N对应于一个电网周期的采样次数;X(k-j)对应最新输入的j点数据,j为正整数,Y(k)为对应的输出,将所得的vα(k)和vβ(k)作为最新的输入数据X(k-j),所得的输出Y(k)即为在αβ坐标系下的基频正序电压分量和
S427、将θ同时输入到频率跟踪器和锁相控制系统;
在频率跟踪器中:故障控制信号动作前一个周期的估算频率是固定的,故保留故障控制信号动作前一个周期的输出再作取倒数和折算后可获得最新的采样周期Ts,Ts可继续应用于测量系统的三相电压的获取及在各模块的运行中进行采样;
在本实施例中还对应用上述锁相环算法的锁相环进行仿真测试。首先,自定义三相电压源,并在2s时,使三相电压的相位同时跳变+30°,锁相环使用本实施例算法对电压进行同步,同时与工程用SIEMENS锁相环的相位同步结果进行对比,仿真结果如图8所示。工程用SIEMENS锁相环的参数设置:平滑滤波时间常数τ为0.02,积分时间常数Ti为0.1;比例系数kp为314,积分系数ki为3140。
从图8可以看到,本实施例的锁相环对于相位跳变的响应速度远远快于目前工程用的SIEMENS锁相环。
最后,为进一步验证本实施例算法对后续换相失败的抑制能力,将应用本实施例锁相环算法的锁相环应用于CIGRE高压直流标准测试模型中,同时对分别采用本实施例锁相环与上述工程用锁相环的在对称故障下抑制后续换相失败的效果进行统计,结果如表1所示。
表1
可以发现,本实施例在换流母线电压跌落至0.44p.u.时,仍然不会发生后续换相失败,而更严重的故障在电网实际中发生的概率极小。此结果得益于本实施例锁相环算法在对称故障下迅速补偿响应的相位跳变值,并将锁相环频率与输出相位解耦,使得本实施例锁相环仍然可以快速准确锁相。
上述实施例为本发明较佳的实施方式,但本发明的实施方式并不受上述实施例的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。
Claims (10)
1.一种可快速同步对称故障下电压相位的高压直流锁相环算法,其特征在于,步骤如下:
S1、故障检测:在锁相环的测量系统中,根据高压直流输电系统的换流母线前三次的相位增量差、频率跟踪器的输出变化率和输入的换相电压幅值,以及换流母线当前a相电压va、b相电压vb和c相电压vc的周期检测结果,检测是否发生对称故障;
S2、选定锁相方式:若判定为对称故障,则将锁相环运行在第二锁相方式;若判定为额定工况或者非对称故障,则将锁相环运行在第一锁相方式;
S3、在发生对称故障的情况下计算相位跳变补偿值,并将对称故障下的相位跳变补偿值输入到锁相环的锁相控制系统中;
S4、按照步骤S2中选定的锁相方式对高压直流换相电压的相位进行同步。
2.根据权利要求1所述的可快速同步对称故障下电压相位的高压直流锁相环算法,其特征在于,步骤S1的故障检测过程具体如下:
S11、对于当前时刻及其前三个时刻,将这些时刻对应的相位增量输出按照时间序列依次进行比较,当相位增量差连续3次大于第一阈值,则生成第一控制信号;
将频率跟踪器的输出变化率与第二阈值进行比较,在输出变化率大于阈值的情况下生成第二控制信号;
获取换流母线当前a相电压va、b相电压vb和c相电压vc,并将这些电压分量在αβ坐标系下的电压幅值与第三阈值比较,当αβ坐标系下的电压幅值大于第三阈值时,生成第三控制信号;
S12、综合第一、二和第三控制信号,生成故障控制信号;
S13、将换流母线当前a相电压va、b相电压vb和c相电压vc分别与各自对应的前一周期的电压值进行比较,小于或等于则表示换流母线正常工作,处于额定工况;
大于则表示响应检测,其中,最先响应检测的相即为故障相,若三个相同时响应,则为对称故障,若只有最先响应的相响应,则为非对称故障。
3.根据权利要求2所述的可快速同步对称故障下电压相位的高压直流锁相环算法,其特征在于,第一阈值为0.3°,第二阈值为15000,第三阈值为0.9p.u.。
5.根据权利要求4所述的可快速同步对称故障下电压相位的高压直流锁相环算法,其特征在于,GDSC输入与输出的关系表示为:
式中,H(k)为单位阶跃函数;为相位跳变前的输入信号;为相位跳变后的信号;为相位跳变前的初相位;即为对称故障引起的相位跳变量;Ts为采样周期;ω为电网角频率;Vpre(k)为跳变前电压信号的幅值;为跳变前相角的指数表示;Vpost(k)为跳变后电压信号的幅值;为跳变后相角的指数表示;
最终,GDSC输入与输出的关系式为:
式中,γ(k)=[1+floor(32k/N)]/32,floor表示向下取整,那么△γ(k)=1/32;
步骤S32生成的比值向量为:
7.根据权利要求1所述的可快速同步对称故障下电压相位的高压直流锁相环算法,其特征在于,在步骤S4中,在额定工况或者非对称故障工况下,按照选定的第一锁相方式对高压直流换相电压的相位进行同步,过程如下:
S411、将上述的换流母线当前a相电压va、b相电压vb和c相电压vc,以采样周期Ts为采样基准,进行插补采样处理;
利用一个固定储存位置来储存采样周期Ts的累加时间进度值t0=kTs;
S412、将上述t0与运行时间t进行比较,运行时间t为运行步长tstep的累加时间,当t≥t0时,输出同步信号sync=1,作为使能信号,同时输出线性插值的百分数frac=(t-t0)/tstep;
S413、考虑到采样周期Ts不一定为运行步长的整数倍,需要对采样信号进行自定义插值处理:
Out′i=Outi-(Outi-Outi-1)×frac
=Outi×frac+Outi-1×(1-frac)
其中,Out′i表示最终的插值处理所得的信号结果;Outi表示最近一次的获取的电压信号;Outi-1表示上一次获取的电压信号;
插值采样后的获取的三相电压信号为:
S414、将上述插补采样后的三相电压信号转换到αβ坐标系下,然后进行Clarke变换:
Clarke变换后获得的电压信号为:
S415、利用FIR数字滤波器从上述电压信号中提取出αβ坐标系下的基频正序电压分量:
S417、将θ同时输入到频率跟踪器和锁相控制系统;
在频率跟踪器中:先将第k次的反正切运算结果与第(k-1)次的结果作差,并将得到的差值乘以N/4,再与π/2作差,经过频率跟踪器的调节器获取额定频率的变换量;
在锁相控制系统中:先将反正切运算结果θ与输出的相位进行比较,
8.根据权利要求1所述的可快速同步对称故障下电压相位的高压直流锁相环算法,其特征在于,在步骤S4中,在对称故障工况下,按照选定的第二锁相方式对高压直流换相电压的相位进行同步,过程如下:
S421、将换流母线当前a相电压va、b相电压vb和c相电压vc,以采样周期Ts为采样基准,进行插补采样处理;
利用一个固定储存位置来储存采样周期Ts的累加时间进度值t0=kTs;
S422、将上述t0与运行时间t进行比较,运行时间t为运行步长tstep的累加时间,当t≥t0时,输出同步信号sync=1,作为使能信号,同时输出线性插值的百分数frac=(t-t0)/tstep;
S423、考虑到采样周期Ts不一定为运行步长的整数倍,需要对采样信号进行自定义插值处理:
Out′i=Outi-(Outi-Outi-1)×frac
=Outi×frac+Outi-1×(1-frac)
其中,Out′i表示最终的插值处理所得的信号结果;Outi表示最近一次的获取的电压信号;Outi-1表示上一次获取的电压信号;
插值采样后的获取的三相电压信号为:
S424、将上述插补采样后的三相电压信号转换到αβ坐标系下,然后进行Clarke变换:
Clarke变换后获得的电压信号为:
S425、利用FIR数字滤波器从上述电压信号中提取出αβ坐标系下的基频正序电压分量:
其中,N对应于一个电网周期的采样次数,X(k-j)对应最新输入的j点数据,j为正整数,Y(k)为对应的输出,将所得的vα(k)和vβ(k)作为最新的输入数据X(k-j),所得的输出Y(k)即为在αβ坐标系下的基频正序电压分量和
S427、将θ同时输入到频率跟踪器和锁相控制系统;
在频率跟踪器中:故障控制信号动作前一个周期的估算频率是固定的,故保留故障控制信号动作前一个周期的输出再作取倒数和折算后可获得最新的采样周期Ts,Ts可继续应用于测量系统的三相电压的获取及在各模块的运行中进行采样;
10.一种可快速同步对称故障下电压相位的高压直流锁相环,其特征在于,所述高压直流锁相环通过权利要求1~9中任一项所述的可快速同步对称故障下电压相位的高压直流锁相环算法,实现与对称故障下高压直流换相电压相位的快速同步。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010516639.0A CN111817713B (zh) | 2020-06-09 | 2020-06-09 | 可快速同步对称故障下电压相位的高压直流锁相环及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010516639.0A CN111817713B (zh) | 2020-06-09 | 2020-06-09 | 可快速同步对称故障下电压相位的高压直流锁相环及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111817713A true CN111817713A (zh) | 2020-10-23 |
CN111817713B CN111817713B (zh) | 2022-04-26 |
Family
ID=72845156
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010516639.0A Active CN111817713B (zh) | 2020-06-09 | 2020-06-09 | 可快速同步对称故障下电压相位的高压直流锁相环及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111817713B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113343456A (zh) * | 2021-05-31 | 2021-09-03 | 南方电网科学研究院有限责任公司 | Hvdc系统高精度闭环动态相量仿真计算、验证方法及装置 |
CN113960397A (zh) * | 2021-10-25 | 2022-01-21 | 矽力杰半导体技术(杭州)有限公司 | 一种故障检测方法和装置 |
CN114370813A (zh) * | 2021-12-28 | 2022-04-19 | 北京无线电计量测试研究所 | 一种u形微波腔电学对称性检测装置 |
CN117007863A (zh) * | 2023-08-03 | 2023-11-07 | 南方电网科学研究院有限责任公司 | 一种电力系统相位角的补偿方法及装置 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101604172A (zh) * | 2009-04-23 | 2009-12-16 | 上海交通大学 | 基于解耦多坐标系统锁相环的电压控制方法 |
CN101877587A (zh) * | 2010-06-30 | 2010-11-03 | 中国电力科学研究院 | 一种新型软锁相环 |
CN102545889A (zh) * | 2011-12-19 | 2012-07-04 | 西安赛博电气有限责任公司 | 基于瞬时无功功率理论的闭环锁相方法及锁相器 |
CN104600696A (zh) * | 2015-01-05 | 2015-05-06 | 上海电力学院 | 一种用于不平衡条件下的软件锁相环实现方法 |
CN104811188A (zh) * | 2015-05-06 | 2015-07-29 | 山东大学 | 基于滑动滤波器的锁相环动态性能改进方法 |
CN109067393A (zh) * | 2018-08-28 | 2018-12-21 | 南方电网科学研究院有限责任公司 | 一种电力系统的锁相方法、装置及设备 |
CN109412191A (zh) * | 2018-12-05 | 2019-03-01 | 华南理工大学 | 一种用于高压直流输电系统的锁相方法、装置和设备 |
CN110571839A (zh) * | 2019-08-05 | 2019-12-13 | 华南理工大学 | 一种基于换相电压相位检测的高压直流故障恢复方法 |
CN111082804A (zh) * | 2019-12-17 | 2020-04-28 | 中国石油大学(华东) | 一种频率补偿型数字锁相环实现方法 |
CN111082447A (zh) * | 2019-12-16 | 2020-04-28 | 华北电力大学 | 一种提升联接弱交流电网vsc-hvdc稳定性的附加相位变化速率负反馈控制方法 |
-
2020
- 2020-06-09 CN CN202010516639.0A patent/CN111817713B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101604172A (zh) * | 2009-04-23 | 2009-12-16 | 上海交通大学 | 基于解耦多坐标系统锁相环的电压控制方法 |
CN101877587A (zh) * | 2010-06-30 | 2010-11-03 | 中国电力科学研究院 | 一种新型软锁相环 |
CN102545889A (zh) * | 2011-12-19 | 2012-07-04 | 西安赛博电气有限责任公司 | 基于瞬时无功功率理论的闭环锁相方法及锁相器 |
CN104600696A (zh) * | 2015-01-05 | 2015-05-06 | 上海电力学院 | 一种用于不平衡条件下的软件锁相环实现方法 |
CN104811188A (zh) * | 2015-05-06 | 2015-07-29 | 山东大学 | 基于滑动滤波器的锁相环动态性能改进方法 |
CN109067393A (zh) * | 2018-08-28 | 2018-12-21 | 南方电网科学研究院有限责任公司 | 一种电力系统的锁相方法、装置及设备 |
CN109412191A (zh) * | 2018-12-05 | 2019-03-01 | 华南理工大学 | 一种用于高压直流输电系统的锁相方法、装置和设备 |
CN110571839A (zh) * | 2019-08-05 | 2019-12-13 | 华南理工大学 | 一种基于换相电压相位检测的高压直流故障恢复方法 |
CN111082447A (zh) * | 2019-12-16 | 2020-04-28 | 华北电力大学 | 一种提升联接弱交流电网vsc-hvdc稳定性的附加相位变化速率负反馈控制方法 |
CN111082804A (zh) * | 2019-12-17 | 2020-04-28 | 中国石油大学(华东) | 一种频率补偿型数字锁相环实现方法 |
Non-Patent Citations (2)
Title |
---|
吴秋媚 等: "一种高压直流锁相环的建模和分析", 《南方电网技术》 * |
郑睿娜 等: "高压直流输电系统交流侧故障检测方法", 《广东电力》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113343456A (zh) * | 2021-05-31 | 2021-09-03 | 南方电网科学研究院有限责任公司 | Hvdc系统高精度闭环动态相量仿真计算、验证方法及装置 |
CN113960397A (zh) * | 2021-10-25 | 2022-01-21 | 矽力杰半导体技术(杭州)有限公司 | 一种故障检测方法和装置 |
CN113960397B (zh) * | 2021-10-25 | 2024-04-16 | 矽力杰半导体技术(杭州)有限公司 | 一种故障检测方法和装置 |
CN114370813A (zh) * | 2021-12-28 | 2022-04-19 | 北京无线电计量测试研究所 | 一种u形微波腔电学对称性检测装置 |
CN114370813B (zh) * | 2021-12-28 | 2023-11-14 | 北京无线电计量测试研究所 | 一种u形微波腔电学对称性检测装置 |
CN117007863A (zh) * | 2023-08-03 | 2023-11-07 | 南方电网科学研究院有限责任公司 | 一种电力系统相位角的补偿方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN111817713B (zh) | 2022-04-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111817713B (zh) | 可快速同步对称故障下电压相位的高压直流锁相环及方法 | |
Rodriguez et al. | Multiple second order generalized integrators for harmonic synchronization of power converters | |
Freijedo et al. | Grid-synchronization methods for power converters | |
CN102305886B (zh) | 电网电压谐波畸变及不平衡时基波电压同步信号检测方法 | |
WO2018122391A1 (en) | Precise real-time advanced grid monitoring | |
Meral | Improved phase-locked loop for robust and fast tracking of three phases under unbalanced electric grid conditions | |
Patil et al. | Modified dual second-order generalised integrator FLL for synchronization of a distributed generator to a weak grid | |
Sinha et al. | A pre-filter based PLL for three-phase grid connected applications | |
Tang et al. | A novel repetitive controller assisted phase-locked loop with self-learning disturbance rejection capability for three-phase grids | |
CN104836255A (zh) | 一种基于隐式pi的数字锁相环与电网同步系统 | |
Pérez-Estévez et al. | Grid-tied inverter with AC voltage sensorless synchronization and soft start | |
JP2011055571A (ja) | 半導体電力変換装置 | |
CN107425743A (zh) | 一种基于预测偏差反馈矫正的并网逆变器mpc方法 | |
CN108763724A (zh) | 一种频率自适应延迟周期的锁相方法 | |
Nagliero et al. | Synchronization techniques for grid connected wind turbines | |
CN113014250B (zh) | 一种可消除直流偏移电压的锁相环及其锁相控制方法 | |
Angquist et al. | Auto-normalizing phase-locked loop for grid-connected converters | |
Agrawal et al. | A modernistic PLL based on feed forward frequency estimator with selective harmonic pre filter for grid imperfection | |
Silva et al. | A robust phase-locked loop against fundamental frequency deviations and harmonic distortions | |
Lee et al. | Grid synchronization PLL robust to frequency variation, unbalanced and distorted voltage | |
CN104184464B (zh) | 一种基于正负序快速识别的动态锁相同步方法 | |
Patil et al. | Modified dual second-order generalized integrator FLL for frequency estimation under various grid abnormalities | |
Bellini et al. | Performances of a PLL based digital filter for double-conversion UPS | |
Zhou et al. | Analysis and comparison of resonant-based current controllers implemented in stationary reference frame: A complex pole-zero placement perspective | |
CN112415266B (zh) | 一种有源电力滤波器负载谐波电流提取方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |