CN110768299A - 负荷虚拟同步机的序阻抗建模与稳定性分析方法 - Google Patents

负荷虚拟同步机的序阻抗建模与稳定性分析方法 Download PDF

Info

Publication number
CN110768299A
CN110768299A CN201911213771.8A CN201911213771A CN110768299A CN 110768299 A CN110768299 A CN 110768299A CN 201911213771 A CN201911213771 A CN 201911213771A CN 110768299 A CN110768299 A CN 110768299A
Authority
CN
China
Prior art keywords
synchronous machine
virtual synchronous
load virtual
disturbance
sequence
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
Application number
CN201911213771.8A
Other languages
English (en)
Other versions
CN110768299B (zh
Inventor
周小平
刘一锋
于浩祺
洪乐荣
夏海涛
朱仁龙
舒钰成
尹寒航
邓凌峰
陈斌超
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hunan University
Original Assignee
Hunan University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hunan University filed Critical Hunan University
Priority to CN201911213771.8A priority Critical patent/CN110768299B/zh
Publication of CN110768299A publication Critical patent/CN110768299A/zh
Application granted granted Critical
Publication of CN110768299B publication Critical patent/CN110768299B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种负荷虚拟同步机的序阻抗建模与稳定性分析方法,综合考虑了负荷虚拟同步机的直流电压控制器、无功控制器、有功控制器、交流电流控制器以及正、负序扰动频率耦合的影响,并采用谐波线性化的方法,建立了负荷虚拟同步机精确的序阻抗模型,解决了多环控制以及频率耦合影响下的负荷虚拟同步机小信号阻抗建模的难题。利用所建的负荷虚拟同步机的序阻抗模型、电网阻抗模型以及奈奎斯特稳定性判据,分析了不同电网阻抗对负荷虚拟同步机并网系统稳定性的影响。分析结果表明:在电网阻抗较大即弱电网的情况下,负荷虚拟同步机依然能稳定运行,不易与电网交互产生振荡。本发明为负荷虚拟同步机接入电网的小信号稳定性分析提供了模型和方法,有利于负荷虚拟同步机的推广与应用。

Description

负荷虚拟同步机的序阻抗建模与稳定性分析方法
技术领域
本发明涉及电力电子变换器与电网的交互稳定性领域,特别是一种负荷虚拟同步机的序阻抗建模与稳定性分析方法。
背景技术
随着以电力电子变换器为接口的可再生能源和负荷大规模、高渗透率地接入电网,由于其缺乏传统同步电机的阻尼与惯性,给电力系统的安全稳定运行带来了不可忽视的影响。为了增强高比例电力电子配电网的惯性,不仅需要新能源并网逆变器具有传统电机的惯性和阻尼,同样需要负荷侧的电力电子变换器能为电网提供惯性支撑并自主参与电网调节。负荷虚拟同步机可以使负荷并网接口具有同步电机的外特性,为负荷友好并网提供了新的技术手段,成为了当下研究的热点。
为了使负荷虚拟同步机能够稳定并网运行,需要合理设计系统模型参数。目前已有文献将负荷虚拟同步机技术应用到电动汽车充电接口、柔性直流输电等领域,并分析了其功率内环以及直流电压外环的小信号模型,给出了系统控制参数优化设计方法,但尚未有文献分析负荷虚拟同步机与电网的交互稳定性问题。
大部分负荷连接在配电网的末端,远离新能源发电,使得电网呈现高感抗弱电网的特性,有必要建立负荷虚拟同步机精确的阻抗模型,进而研究与电网的交互稳定性,为负荷虚拟同步机技术的推广和应用提供理论技术支撑。
发明内容
本发明所要解决的技术问题是,针对现有技术不足,提供一种负荷虚拟同步机的序阻抗建模与稳定性分析方法,解决考虑直流电压控制器、无功控制器、有功控制器、交流电流控制器以及正、负序扰动频率耦合影响下的负荷虚拟同步机小信号阻抗建模的难题。
为解决上述技术问题,本发明所采用的技术方案是:一种负荷虚拟同步机的序阻抗建模与稳定性分析方法,该方法主要实现过程如下:在三相静止坐标系下建立负荷虚拟同步机精确的小信号序阻抗模型,利用所建的负荷虚拟同步机小信号序阻抗模型、电网阻抗模型以及奈奎斯特稳定性判据,分析不同电网阻抗对负荷虚拟同步机并网系统稳定性的影响。
在三相静止坐标系下建立负荷虚拟同步机精确的小信号序阻抗模型的具体实现过程包括:
1)在时域中,向负荷虚拟同步机交流侧施加正序小信号电压扰动,则正序小信号电压扰动通过负荷虚拟同步机的控制作用后,产生正序调制波扰动、耦合的负序调制波扰动、正序交流电流扰动、耦合的负序交流电流扰动和零序的直流电压扰动,将上述扰动转化到频域,以A相和频域中的正频率分量为例,得到负荷虚拟同步机的交流电压va正序扰动分量Vp、调制波ma正序扰动分量Mp、耦合的负序调制波扰动分量Mp2、交流电流ia正序扰动分量Ip、耦合的负序交流电流扰动分量Ip2以及直流电压vdc零序分量Vpdc;根据负荷虚拟同步机的主电路拓扑结构,利用频域卷积定理和谐波线性化方法,得到负荷虚拟同步机的主电路小信号模型在频域中的表达式如下:
sLfIp=Vp-KmVdcMp-KmM1Vpdc
(s-j2ω1)LfIp2=-KmVdcMp2-KmM-1Vpdc
[(s-jω1)Cdc+1/RL]Vpdc=3Km(MpI-1+Mp2I1+M-1Ip+M1Ip2)
式中:Lf表示交流滤波电感;Km表示调制度;RL为直流侧负载等效电阻;Cdc为直流侧电容;ω1表示基波角频率;Vdc表示直流侧电压参考值;M1和M-1分别对应调制波基波分量在频域中的正频率分量和负频率分量;Ι1和I-1分别对应交流电流基波分量在频域中的正频率分量和负频率分量;
2)根据瞬时功率理论可得负荷虚拟同步机输入的有功功率Pe和无功功率Qe,再根据Vp、Ip、Ip2和频域卷积定理,可得在交流电压正序小信号扰动下Pe和Qe扰动的频域表达式如下:
Pe[fp-f1]=3(V-1Ip+VpI-1+V1Ip2)
Qe[fp-f1]=j3(V-1Ip-VpI-1-V1Ip2)
式中:fp、f1分别表示基波频率和正序扰动频率;V1和V-1分别对应交流电压基波分量在频域中的正频率分量和负频率分量;
3)根据负荷虚拟同步机的无功控制器可得其输出内电势幅值em,在交流电压正序小信号扰动下,结合无功功率扰动的频域表达式,可得em扰动在频域中的表达式如下:
Εm[fp-f1]=j3kq(V-1Ip-VpI-1-V1Ip2)
式中:kq表示无功调节系数;
4)根据负荷虚拟同步机的直流电压控制器和有功控制器可得负荷虚拟同步机调制波的相角θ,在交流电压正序小信号扰动下,结合直流电压扰动和有功功率扰动的频域表达式,可得相角扰动Δθ在频域中的表达式如下:
Δθ[fp-f1]=T(s){Pe[fp-f1]/ω1-Hv(s)Vpdc}
式中:
Figure BDA0002298911330000031
J表示虚拟惯性系数;D表示阻尼系数;Pr表示有功功率参考值;Hv(s)表示直流侧电压PI控制器的传递函数;
5)负荷虚拟同步机控制生成的正弦信号cosθ≈cosθ1-Δθsinθ1,其中θ1为负荷虚拟同步机调制波相角的基波分量,结合相角扰动Δθ的频域表达式,利用频域卷积定理并考虑频率耦合影响,可得在交流电压正序小信号扰动下cosθ的频域表达式如下:
Figure BDA0002298911330000032
式中:
Figure BDA0002298911330000033
表示负荷虚拟同步机内电势与并网点电压之间的相位差;Em为内电势幅值em的参考值。
6)根据负荷虚拟同步机电气部分控制可得调制波表达式,结合内电势幅值em、正弦信号cosθ、交流电压、交流电流和直流电压扰动的频域表达式,利用频域卷积定理并考虑电压电流采样延时、PWM延时和频率耦合的影响,可得在交流电压正序小信号扰动下调制波正序扰动分量Mp、耦合的负序扰动分量Mp2在频域中的表达式如下:
Figure BDA0002298911330000034
Figure BDA0002298911330000035
式中:Gv(s)、Gi(s)分别表示电压和电流采样延时传递函数;Hr(s)表示交流电流PR控制传递函数;Gm(s)表示PWM延时传递函数;
7)将负荷虚拟同步机调制波扰动表达式代入主电路小信号模型频域表达式中,可以求得负荷虚拟同步机在交流电压正序小信号扰动Vp作用下产生交流电流正序扰动Ip对应的导纳Ypp(s)以及耦合出的交流电流负序扰动Ip2对应的导纳Ypn(s);
8)根据正序阻抗与负序阻抗的共轭关系,可以得到负荷虚拟同步机在交流电压负序小信号扰动Vp2作用下产生交流电流负序扰动Ip2对应的导纳Ynn(s)以及耦合产生交流电流正序扰动Ip对应的导纳Ynp(s);
9)考虑电网阻抗影响,则交流电压正序扰动Vp通过负荷虚拟同步机耦合产生的交流电流负序扰动Ip2会通过电网阻抗产生负序电压扰动Vp2,Vp2又会施加在负荷虚拟同步机上产生Ip2和Ip,根据这一扰动传递过程和正、负序阻抗的共轭关系即可求得负荷虚拟同步机的正、负序阻抗模型Zp(s)和Zn(s)计算公式为:
Figure BDA0002298911330000041
Figure BDA0002298911330000042
式中:符号*表示对复数的共轭;ω1表示基波角频率;
Figure BDA0002298911330000044
分别表示对Ypn(s)和Ypp(s)取共轭并用s-j2ω1代替复变量s;
Yg(s-j2ω1)表示电网导纳;
Figure BDA0002298911330000045
表示对Zp(s)取共轭;
利用所建的负荷虚拟同步机小信号序阻抗模型、电网阻抗模型以及奈奎斯特稳定性判据,分析不同电网阻抗对负荷虚拟同步机并网系统稳定性的影响的具体实现过程包括:
1)根据电网的阻抗特性,建立电网的正、负序阻抗模型分别为Zgp(s)和Zgn(s);
2)分别计算负荷虚拟同步机并网系统的正序阻抗比Zgp(s)/Zp(s)和负序阻抗比Zgn(s)/Zn(s);Zp(s)和Zn(s)分别为负荷虚拟同步机的正序阻抗模型、负序阻抗模型;
3)绘制正序阻抗比与负序阻抗比的奈奎斯特曲线,利用奈奎斯特稳定性判据进行稳定性分析,在电网阻抗变化时,只有当正序阻抗比和负序阻抗比都满足奈奎斯特稳定性判据时,系统才稳定。
与现有技术相比,本发明所具有的有益效果为:本发明全面考虑了负荷虚拟同步机的控制环,在三相静止坐标系下建立了负荷虚拟同步机小信号序阻抗模型,所建的阻抗模型精度高,具有明确的物理意义,稳定性分析方法简便有效等优点;本发明为负荷虚拟同步机接入电网的小信号稳定性分析提供了模型和方法,有利于负荷虚拟同步机的推广与应用。
附图说明
图1为本发明一实施例负荷虚拟同步机的主电路拓扑结构;
图2为本发明一实施例负荷虚拟同步机的控制框图;
图3为本发明一实施例负荷虚拟同步机并网系统的扰动传递信号流图;
图4为本发明一实施例负荷虚拟同步机的正、负序阻抗特性曲线及其仿真测量结果;
图5为本发明一实施例不同电网阻抗下系统的判稳正、负序阻抗比的奈奎斯特图。
具体实施方式
图1为负荷虚拟同步机的主电路拓扑结构。其中:Q1-Q6构成了负荷虚拟同步机三相整流桥;vdc、idc分别为负荷虚拟同步机的直流侧电压和电流;ea、eb和ec为负荷虚拟同步机的内电势;va、vb和vc为负荷虚拟同步机的并网点电压;ia、ib、和ic为负荷虚拟同步机的交流输入电流;vga、vgb和vgc为电网电压,iga、igb和igc为电网电流;RL为直流侧负载等效电阻;Cdc为直流侧电容;Lf、Cf和Rf分别为交流滤波电感、滤波电容和阻尼电阻;Lg和Rg分别为电网的等效线路电感和电阻。
图2为负荷虚拟同步机的控制框图。其通过控制模拟传统同步机电机的机电暂态方程,使得变换器具有类似的惯性、有功调频和无功调压的特性。同时,为了保证并网电流的质量,根据变换器输出电压指令和电气关系计算得到输出电流指令,再对输出电流进行比例谐振控制。负荷虚拟同步机的有功控制器、无功控制器和电气部分控制的数学方程为:
D(ω-ω0)+Pe/ω-(vdc-vdcr)Hv(s)=Jsω
(1)
ω=sθ
(2)
em=kq(Qe-Qr)+E0
(3)
ma=[(va-emcosθ)/(sLf)-ia]Hr(s)
(4)
式中:ω和ω0分别表示负荷虚拟同步机输出角频率和电网额定角频率;J表示虚拟惯性系数;D表示阻尼系数;θ表示负荷虚拟同步机调制波的相角;vdcr表示负荷虚拟同步机直流侧电压指令值;em表示负荷虚拟同步机输出内电势幅值;Qr表示负荷虚拟同步机无功功率指令值;E0表示负荷虚拟同步机空载电动势;kq表示无功调节系数;ma表示负荷虚拟同步机a相调制波;Hv(s)表示直流侧电压PI控制器的传递函数;Hr(s)表示交流电流PR控制传递函数;
一种负荷虚拟同步机的序阻抗建模与稳定性分析方法,其特征在于,在三相静止坐标系下建立负荷虚拟同步机精确的小信号序阻抗模型,利用所建的负荷虚拟同步机小信号序阻抗模型、电网阻抗模型以及奈奎斯特稳定性判据,分析了不同电网阻抗对负荷虚拟同步机并网系统稳定性的影响。
所述在三相静止坐标系下建立负荷虚拟同步机精确的小信号序阻抗模型的具体实现过程包括:
1)在时域中,向负荷虚拟同步机交流侧施加正序小信号电压扰动,则正序小信号电压扰动通过负荷虚拟同步机的控制作用后,产生正序调制波扰动、耦合的负序调制波扰动、正序交流电流扰动、耦合的负序交流电流扰动和零序的直流电压扰动,将上述扰动转化到频域,以A相和频域中的正频率分量为例,得到负荷虚拟同步机的交流电压va正序扰动分量Vp、调制波ma正序扰动分量Mp、耦合的负序调制波扰动分量Mp2、交流电流ia正序扰动分量Ip、耦合的负序交流电流扰动分量Ip2以及直流电压vdc零序分量Vpdc;根据负荷虚拟同步机的主电路拓扑结构,利用频域卷积定理和谐波线性化方法,得到负荷虚拟同步机的主电路小信号模型在频域中的表达式如下:
sLfIp=Vp-KmVdcMp-KmM1Vpdc
(5)
(s-j2ω1)LfIp2=-KmVdcMp2-KmM-1Vpdc
(6)
[(s-jω1)Cdc+1/RL]Vpdc=3Km(MpI-1+Mp2I1+M-1Ip+M1Ip2)
(7)
式中:Lf表示交流滤波电感;Km表示调制度;RL为直流侧负载等效电阻;Cdc为直流侧电容;ω1表示基波角频率;Vdc表示直流侧电压参考值;M1和M-1分别对应调制波基波分量在频域中的正频率分量和负频率分量;Ι1和I-1分别对应交流电流基波分量在频域中的正频率分量和负频率分量;
2)根据瞬时功率理论可得负荷虚拟同步机输入的有功功率Pe和无功功率Qe,再根据Vp、Ip、Ip2和频域卷积定理,可得在交流电压正序小信号扰动下Pe和Qe扰动的频域表达式如下:
Pe[fp-f1]=3(V-1Ip+VpI-1+V1Ip2) (8)
Qe[fp-f1]=j3(V-1Ip-VpI-1-V1Ip2) (9)
式中:fp、f1分别表示基波频率和正序扰动频率;V1和V-1分别对应交流电压基波分量在频域中的正频率分量和负频率分量;
3)根据负荷虚拟同步机的无功控制器可得其输出内电势幅值em,在交流电压正序小信号扰动下,结合无功功率扰动的频域表达式,可得em扰动在频域中的表达式如下:
Εm[fp-f1]=j3kq(V-1Ip-VpI-1-V1Ip2) (10)
式中:kq表示无功调节系数;
4)根据负荷虚拟同步机的直流电压控制器和有功控制器可得负荷虚拟同步机调制波的相角θ,在交流电压正序小信号扰动下,结合直流电压扰动和有功功率扰动的频域表达式,可得相角扰动Δθ在频域中的表达式如下:
Δθ[fp-f1]=T(s){Pe[fp-f1]/ω1-Hv(s)Vpdc} (11)
式中:
Figure BDA0002298911330000081
J表示虚拟惯性系数;D表示阻尼系数;Pr表示有功功率参考值;Hv(s)表示直流侧电压PI控制器的传递函数;
5)负荷虚拟同步机控制生成的正弦信号cosθ≈cosθ1-Δθsinθ1,其中θ1为负荷虚拟同步机调制波相角的基波分量,结合相角扰动Δθ的频域表达式,利用频域卷积定理并考虑频率耦合影响,可得在交流电压正序小信号扰动下cosθ的频域表达式如下:
Figure BDA0002298911330000082
式中:表示负荷虚拟同步机内电势与并网点电压之间的相位差;Em为内电势幅值em的参考值。
6)根据负荷虚拟同步机电气部分控制可得调制波表达式,结合内电势幅值em、正弦信号cosθ、交流电压、交流电流和直流电压扰动的频域表达式,利用频域卷积定理并考虑电压电流采样延时、PWM延时和频率耦合的影响,可得在交流电压正序小信号扰动下调制波正序扰动分量Mp、耦合的负序扰动分量Mp2在频域中的表达式如下:
Figure BDA0002298911330000084
Figure BDA0002298911330000085
式中:Gv(s)、Gi(s)分别表示电压和电流采样延时传递函数;Hr(s)表示交流电流PR控制传递函数;Gm(s)表示PWM延时传递函数。
7)将负荷虚拟同步机调制波扰动表达式代入主电路小信号模型频域表达式中,可以求得负荷虚拟同步机在交流电压正序小信号扰动Vp作用下产生交流电流正序扰动Ip对应的导纳Ypp(s)以及耦合出的交流电流负序扰动Ip2对应的导纳Ypn(s);
8)根据正序阻抗与负序阻抗的共轭关系,可以得到负荷虚拟同步机在交流电压负序小信号扰动Vp2作用下产生交流电流负序扰动Ip2对应的导纳Ynn(s)以及耦合产生交流电流正序扰动Ip对应的导纳Ynp(s);
9)考虑电网阻抗影响,则交流电压正序扰动Vp通过负荷虚拟同步机耦合产生的交流电流负序扰动Ip2会通过电网阻抗产生负序电压扰动Vp2,Vp2又会施加在负荷虚拟同步机上产生Ip2和Ip,根据这一扰动传递过程和正、负序阻抗的共轭关系即可求得负荷虚拟同步机的正、负序阻抗模型Zp(s)和Zn(s)计算公式为:
Figure BDA0002298911330000091
Figure BDA0002298911330000092
式中:符号*表示对复数的共轭;ω1表示基波角频率;分别表示对Ypn(s)和Ypp(s)取共轭并用s-j2ω1代替复变量s;Yg(s-j2ω1)表示电网导纳;
Figure BDA0002298911330000095
表示对Zp(s)取共轭;
所述基于所建的负荷虚拟同步机小信号序阻抗模型、电网阻抗模型以及奈奎斯特稳定性判据,分析不同电网阻抗对负荷虚拟同步机并网系统稳定性的影响的具体实现过程包括:
1)根据电网的阻抗特性,建立电网的正、负序阻抗模型分别为Zgp(s)和Zgn(s);
2)分别计算负荷虚拟同步机并网系统的正序阻抗比Zgp(s)/Zp(s)和负序阻抗比Zgn(s)/Zn(s);Zp(s)和Zn(s)分别为负荷虚拟同步机的正序阻抗模型、负序阻抗模型;
3)绘制正序阻抗比与负序阻抗比的奈奎斯特曲线,利用奈奎斯特稳定性判据进行稳定性分析,在电网阻抗变化时,只有当正序阻抗比和负序阻抗比都满足奈奎斯特稳定性判据时,系统才稳定。
图3为负荷虚拟同步机并网系统的扰动传递信号流图。如图3所示,考虑电网阻抗的影响,交流电压正序扰动Vp通过负荷虚拟同步机耦合产生的交流电流负序扰动Ip2会通过电网阻抗产生负序电压扰动Vp2,Vp2又会施加在负荷虚拟同步机上产生Ip2和Ip
图4为负荷虚拟同步机的正、负序阻抗特性曲线及其仿真测量结果。图中Zg(s)表示电网阻抗实际值,其正、负序阻抗是一样的,Zp(s)和Zn(s)分别表示负荷虚拟同步机的正、负序阻抗(实线为理论值,圆圈代表正序阻抗的测量值,星号代表负序阻抗的测量值)。从图中可以看出:阻抗测量的结果和所建的阻抗模型吻合,验证了所推导的负荷虚拟同步机的序阻抗模型的正确性。
图5为不同电网阻抗下系统的判稳正、负序阻抗比的奈奎斯特图。正序阻抗比用实线表示,负序阻抗比用虚线表示,Lg为电网等效电感。负荷虚拟同步机的正序阻抗无右半平面零点,而负序阻抗有一对共轭的右半平面零点。因此,当其正序阻抗比奈奎斯特曲线不包围(-1,j0)点以及负序阻抗比奈奎斯特曲线从(-1,j0)左侧正穿越一次时,系统稳定。从图5中可以看出,随着电网阻抗的增大,负荷虚拟同步机始终能保持稳定。因此,负荷虚拟同步机具有较强的适应弱电网的能力。

Claims (3)

1.一种负荷虚拟同步机的序阻抗建模与稳定性分析方法,其特征在于,该方法主要实现过程如下:在三相静止坐标系下建立负荷虚拟同步机的小信号序阻抗模型,利用所建的负荷虚拟同步机的小信号序阻抗模型、电网阻抗模型以及奈奎斯特稳定性判据,分析不同电网阻抗对负荷虚拟同步机并网系统稳定性的影响。
2.根据权利要求1所述的负荷虚拟同步机的序阻抗建模与稳定性分析方法,其特征在于,在三相静止坐标系下建立负荷虚拟同步机的小信号序阻抗模型的具体实现过程包括:
1)在负荷虚拟同步机交流侧施加正序小信号电压扰动,获取负荷虚拟同步机的交流电压va正序扰动分量Vp、调制波ma正序扰动分量Mp、耦合的负序调制波扰动分量Mp2、交流电流ia正序扰动分量Ip、耦合的负序交流电流扰动分量Ip2以及直流电压vdc零序分量Vpdc,得到负荷虚拟同步机的主电路小信号模型在频域中的表达式如下:
sLfIp=Vp-KmVdcMp-KmM1Vpdc
(s-j2ω1)LfIp2=-KmVdcMp2-KmM-1Vpdc
[(s-jω1)Cdc+1/RL]Vpdc=3Km(MpI-1+Mp2I1+M-1Ip+M1Ip2);
式中:Lf表示交流滤波电感;Km表示调制度;RL为直流侧负载等效电阻;Cdc为直流侧电容;ω1表示基波角频率;Vdc表示直流侧电压参考值;M1和M-1分别对应调制波基波分量在频域中的正频率分量和负频率分量;Ι1和I-1分别对应交流电流基波分量在频域中的正频率分量和负频率分量;
2)根据瞬时功率理论得到负荷虚拟同步机输入的有功功率Pe和无功功率Qe,再根据Vp、Ip、Ip2和频域卷积定理,得到在交流电压正序小信号扰动下Pe和Qe扰动的频域表达式如下:
Pe[fp-f1]=3(V-1Ip+VpI-1+V1Ip2);
Qe[fp-f1]=j3(V-1Ip-VpI-1-V1Ip2);
式中:fp、f1分别表示基波频率和正序扰动频率;V1和V-1分别对应交流电压基波分量在频域中的正频率分量和负频率分量;
3)根据负荷虚拟同步机的无功控制器得到其输出内电势幅值em,在交流电压正序小信号扰动下,结合无功功率扰动的频域表达式,得到em在频域中的表达式如下:
Εm[fp-f1]=j3kq(V-1Ip-VpI-1-V1Ip2);
式中:kq表示无功调节系数;Em为小信号扰动下,em在频域中的表达式;
4)根据负荷虚拟同步机的直流电压控制器和有功控制器得到负荷虚拟同步机调制波的相角θ,在交流电压正序小信号扰动下,结合直流电压扰动和有功功率扰动的频域表达式,得到相角扰动Δθ在频域中的表达式如下:
Δθ[fp-f1]=T(s){Pe[fp-f1]/ω1-Hv(s)Vpdc};
式中:
Figure FDA0002298911320000021
J表示虚拟惯性系数;D表示阻尼系数;Pr表示有功功率参考值;Hv(s)表示直流侧电压PI控制器的传递函数;
5)负荷虚拟同步机控制生成的正弦信号cosθ≈cosθ1-Δθsinθ1,其中θ1为负荷虚拟同步机调制波相角的基波分量,结合相角扰动Δθ的频域表达式,利用频域卷积定理并考虑频率耦合影响,得到在交流电压正序小信号扰动下cosθ的频域表达式如下:
Figure FDA0002298911320000022
式中:
Figure FDA0002298911320000023
表示负荷虚拟同步机内电势与并网点电压之间的相位差;Em为内电势幅值em的参考值;
6)根据负荷虚拟同步机电气部分控制得到调制波表达式,结合内电势幅值em、正弦信号cosθ、交流电压、交流电流和直流电压扰动的频域表达式,利用频域卷积定理并考虑电压电流采样延时、PWM延时和频率耦合的影响,得到在交流电压正序小信号扰动下调制波正序扰动分量Mp、耦合的负序扰动分量Mp2在频域中的表达式如下:
Figure FDA0002298911320000024
Figure FDA0002298911320000025
式中:Gv(s)、Gi(s)分别表示电压和电流采样延时传递函数;Hr(s)表示交流电流PR控制传递函数;Gm(s)表示PWM延时传递函数;
7)将负荷虚拟同步机调制波扰动表达式代入主电路小信号模型频域表达式中,求得负荷虚拟同步机在交流电压正序小信号扰动Vp作用下产生交流电流正序扰动Ip对应的导纳Ypp(s)以及耦合出的交流电流负序扰动Ip2对应的导纳Ypn(s);
8)根据正序阻抗与负序阻抗的共轭关系,得到负荷虚拟同步机在交流电压负序小信号扰动Vp2作用下产生交流电流负序扰动Ip2对应的导纳Ynn(s)以及耦合产生交流电流正序扰动Ip对应的导纳Ynp(s);
9)求得负荷虚拟同步机的正、负序阻抗模型Zp(s)和Zn(s)计算公式为:
Figure FDA0002298911320000031
式中:符号*表示对复数的共轭;ω1表示基波角频率;
Figure FDA0002298911320000032
Figure FDA0002298911320000033
分别表示对Ypn(s)和Ypp(s)取共轭并用s-j2ω1代替复变量s;Yg(s-j2ω1)表示电网导纳;表示对Zp(s)取共轭。
3.根据权利要求1或2所述的负荷虚拟同步机的序阻抗建模与稳定性分析方法,其特征在于,利用所建的负荷虚拟同步机的小信号序阻抗模型、电网阻抗模型以及奈奎斯特稳定性判据,分析不同电网阻抗对负荷虚拟同步机并网系统稳定性的影响的具体实现过程包括:
1)根据电网的阻抗特性,建立电网的正、负序阻抗模型分别为Zgp(s)和Zgn(s);
2)分别计算负荷虚拟同步机并网系统的正序阻抗比Zgp(s)/Zp(s)和负序阻抗比Zgn(s)/Zn(s);Zp(s)和Zn(s)分别为负荷虚拟同步机的正序阻抗模型、负序阻抗模型;
3)绘制正序阻抗比与负序阻抗比的奈奎斯特曲线,利用奈奎斯特稳定性判据进行稳定性分析,在电网阻抗变化时,只有当正序阻抗比和负序阻抗比都满足奈奎斯特稳定性判据时,系统才稳定。
CN201911213771.8A 2019-12-02 2019-12-02 负荷虚拟同步机的序阻抗建模与稳定性分析方法 Active CN110768299B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911213771.8A CN110768299B (zh) 2019-12-02 2019-12-02 负荷虚拟同步机的序阻抗建模与稳定性分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911213771.8A CN110768299B (zh) 2019-12-02 2019-12-02 负荷虚拟同步机的序阻抗建模与稳定性分析方法

Publications (2)

Publication Number Publication Date
CN110768299A true CN110768299A (zh) 2020-02-07
CN110768299B CN110768299B (zh) 2021-06-25

Family

ID=69340779

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911213771.8A Active CN110768299B (zh) 2019-12-02 2019-12-02 负荷虚拟同步机的序阻抗建模与稳定性分析方法

Country Status (1)

Country Link
CN (1) CN110768299B (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111525555A (zh) * 2020-04-29 2020-08-11 北京交通大学 一种虚拟同步机黑箱模型的建模方法
CN111555339A (zh) * 2020-02-25 2020-08-18 暨南大学 用于稳定性分析的变流器并网通用序阻抗模型及建模方法
CN111931343A (zh) * 2020-05-21 2020-11-13 北京理工大学 用于永磁同步电机驱动系统的稳定性分析方法及系统
CN113013926A (zh) * 2021-03-22 2021-06-22 国网山东综合能源服务有限公司 一种分布式并网发电系统序阻抗聚合方法及系统
CN113285447A (zh) * 2021-05-21 2021-08-20 上海电力大学 一种风储黑启动长线路参数变化谐振分析方法
CN113419588A (zh) * 2021-07-01 2021-09-21 湘潭大学 一种基于阻抗特性的两级级联变换器系统稳定性分析方法
CN113890096A (zh) * 2021-10-20 2022-01-04 东南大学 一种整流逆变模式下的变流器序阻抗建模方法
CN115632401A (zh) * 2022-12-21 2023-01-20 浙江大学 一种计及负载和电网阻抗影响的sapf参数设计方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016015046A1 (en) * 2014-07-25 2016-01-28 Axsun Technologies Llc Real time fpga resampling for swept source optical coherence tomography
CN107994606A (zh) * 2017-11-27 2018-05-04 浙江大学 一种多频率耦合因素共存情况下并网逆变器系统稳定性分析方法
CN108281986A (zh) * 2018-03-26 2018-07-13 湖南大学 电压控制型虚拟同步发电机的阻抗建模与稳定性分析方法
CN108418253A (zh) * 2018-03-26 2018-08-17 湖南大学 电流控制型虚拟同步发电机的阻抗建模与稳定性分析方法
CN109921413A (zh) * 2019-02-28 2019-06-21 湖南大学 三相并网逆变器稳定性分析方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016015046A1 (en) * 2014-07-25 2016-01-28 Axsun Technologies Llc Real time fpga resampling for swept source optical coherence tomography
CN107994606A (zh) * 2017-11-27 2018-05-04 浙江大学 一种多频率耦合因素共存情况下并网逆变器系统稳定性分析方法
CN108281986A (zh) * 2018-03-26 2018-07-13 湖南大学 电压控制型虚拟同步发电机的阻抗建模与稳定性分析方法
CN108418253A (zh) * 2018-03-26 2018-08-17 湖南大学 电流控制型虚拟同步发电机的阻抗建模与稳定性分析方法
CN109921413A (zh) * 2019-02-28 2019-06-21 湖南大学 三相并网逆变器稳定性分析方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
ATLE RYGG, MARTA MOLINAS: "A Modified Sequence-Domain Impedance Definition and Its Equivalence to the dq-Domain Impedance Definition for the Stability Analysis of AC Power Electronic Systems", 《IEEE JOURNAL OF EMERGING AND SELECTED TOPICS IN POWER ELECTRONICS》 *
MAURICIO CESPEDES: "Impedance Modeling and Analysis of Grid-Connected Voltage-Source Converters", 《IEEE TRANSACTIONS ON POWER ELECTRONICS》 *
伍文华: "序阻抗视角下虚拟同步发电机与传统并网逆变器的稳定性对比分析", 《中国电机工程学报》 *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111555339A (zh) * 2020-02-25 2020-08-18 暨南大学 用于稳定性分析的变流器并网通用序阻抗模型及建模方法
CN111555339B (zh) * 2020-02-25 2022-01-18 暨南大学 用于稳定性分析的变流器并网通用序阻抗模型及建模方法
CN111525555A (zh) * 2020-04-29 2020-08-11 北京交通大学 一种虚拟同步机黑箱模型的建模方法
CN111931343B (zh) * 2020-05-21 2023-03-24 北京理工大学 用于永磁同步电机驱动系统的稳定性分析方法及系统
CN111931343A (zh) * 2020-05-21 2020-11-13 北京理工大学 用于永磁同步电机驱动系统的稳定性分析方法及系统
CN113013926A (zh) * 2021-03-22 2021-06-22 国网山东综合能源服务有限公司 一种分布式并网发电系统序阻抗聚合方法及系统
CN113013926B (zh) * 2021-03-22 2024-02-13 国网山东综合能源服务有限公司 一种分布式并网发电系统序阻抗聚合方法及系统
CN113285447A (zh) * 2021-05-21 2021-08-20 上海电力大学 一种风储黑启动长线路参数变化谐振分析方法
CN113419588A (zh) * 2021-07-01 2021-09-21 湘潭大学 一种基于阻抗特性的两级级联变换器系统稳定性分析方法
CN113890096A (zh) * 2021-10-20 2022-01-04 东南大学 一种整流逆变模式下的变流器序阻抗建模方法
CN113890096B (zh) * 2021-10-20 2023-12-22 东南大学 一种整流逆变模式下的变流器序阻抗建模方法
CN115632401B (zh) * 2022-12-21 2023-03-07 浙江大学 一种计及负载和电网阻抗影响的sapf参数设计方法
CN115632401A (zh) * 2022-12-21 2023-01-20 浙江大学 一种计及负载和电网阻抗影响的sapf参数设计方法

Also Published As

Publication number Publication date
CN110768299B (zh) 2021-06-25

Similar Documents

Publication Publication Date Title
CN110768299B (zh) 负荷虚拟同步机的序阻抗建模与稳定性分析方法
CN108418253B (zh) 电流控制型虚拟同步发电机的阻抗建模与稳定性分析方法
JP6265826B2 (ja) 単相系統に接続される電力変換装置
CN108154315B (zh) 一种考虑锁相环影响的并网变流器次同步振荡风险分析方法
CN108923463B (zh) 考虑锁相环的单相lcl型并网逆变器的频率耦合建模方法
CN108281986B (zh) 电压控制型虚拟同步发电机的阻抗建模与稳定性分析方法
US9729085B2 (en) Observer-based control system for grid-connected DC/AC converters with LCL-filter
CN108280271B (zh) 基于开关周期平均原理的统一潮流控制器等效建模方法
CN110676874B (zh) 计及频率耦合效应的直驱式风机次同步振荡电气量分析方法
CN109217371B (zh) 考虑锁相环影响的电压源型变流器并网系统稳定性分析方法、装置及系统
CN103066876B (zh) 一种具有同步发电机特性的逆变控制器
CN106532749B (zh) 一种微电网不平衡功率和谐波电压补偿系统及其应用
CN106786647A (zh) 一种三相四线制并联apf双闭环非线性复合控制方法
CN110739678A (zh) 一种并网换流器串联虚拟阻抗的控制方法
CN107611971A (zh) 针对电网电压谐波畸变工况的网侧逆变器谐振全阶滑模控制方法
CN107276091A (zh) Npc型三电平三相四线制sapf非线性无源控制方法
CN105406741B (zh) 一种三相电网电压不平衡时pwm整流器模糊滑模变结构控制方法
CN103678827B (zh) 一种用于逆变器的电磁暂态建模方法
CN106712115A (zh) 一种不用电流反馈的虚拟同步发电机控制器
CN107947237A (zh) 一种多类型逆变器孤岛微网稳定性分析方法
CN109802385B (zh) 电压源型逆变器的阻抗建模方法
CN113346513A (zh) 一种辨识直驱风机强迫次同步振荡的方法
CN103117562A (zh) 一种高压级联能量回馈变频器功率模块的控制方法
CN104362675A (zh) 光伏并网逆变器控制方法、装置及逆变器系统
CN103973151A (zh) 电感不平衡条件下三相pwm并网逆变器的解耦控制方法

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