CN106877309A - 主从控制模式下直流微电网的稳定性判定方法 - Google Patents

主从控制模式下直流微电网的稳定性判定方法 Download PDF

Info

Publication number
CN106877309A
CN106877309A CN201710005641.XA CN201710005641A CN106877309A CN 106877309 A CN106877309 A CN 106877309A CN 201710005641 A CN201710005641 A CN 201710005641A CN 106877309 A CN106877309 A CN 106877309A
Authority
CN
China
Prior art keywords
node
admittance
matrix
stability
power supply
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
CN201710005641.XA
Other languages
English (en)
Other versions
CN106877309B (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.)
State Grid Zhejiang Xiaoshan District Power Supply Co ltd
Zhejiang Zhongxin Electric Power Engineering Construction Co Ltd
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201710005641.XA priority Critical patent/CN106877309B/zh
Publication of CN106877309A publication Critical patent/CN106877309A/zh
Application granted granted Critical
Publication of CN106877309B publication Critical patent/CN106877309B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J1/00Circuit arrangements for dc mains or dc distribution networks
    • H02J1/10Parallel operation of dc sources
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P80/00Climate change mitigation technologies for sector-wide applications
    • Y02P80/10Efficient use of energy, e.g. using compressed air or pressurized fluid as energy carrier
    • Y02P80/14District level solutions, i.e. local energy networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Power Engineering (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明涉及直流微电网技术中的稳定性分析技术,旨在提供一种主从控制模式下直流微电网的稳定性判定方法。包括:利用单元连接法对各电源、恒功率负载的功率模块进行戴维南/诺顿等效,得到相应的阻抗/导纳,并获得主电源的阻抗矩阵、从电源和恒功率负载的导纳矩阵;根据已知的各节点电阻性负载大小得到阻性负载导纳矩阵,根据节点间线缆参数得到系统节点导纳矩阵;将各矩阵带入稳定性判据并计算其极点;如果极点全部在左半平面,则系统稳定;如果有右半平面极点,则系统不稳定。本发明能兼容考虑系统节点间阻抗,对系统结构无特定要求,适用范围广。分析系统级稳定性时简化计算,可有效分析复杂直流微电网的稳定性。

Description

主从控制模式下直流微电网的稳定性判定方法
技术领域
本发明涉及主从控制模式下直流微电网的稳定性判定方法,属于直流微电网技术中的稳定性分析技术。
背景技术
随着可再生新能源的大量接入,直流微电网作为一种高效利用的系统被广泛应用。直流微电网一般由多个电源、恒功率负载等功率模块组成,多个电源之间采用主从控制,即一个电源是电压控制型作为主电源,其余电源是电流控制型作为从电源。通常单个功率模块在控制设计时能确保独立运行时稳定,但组成系统时将会互相影响,降低系统稳定性。另一方面,恒功率负载的负阻抗外特性会进一步加剧系统不稳定。如何分析直流微电网的稳定性是很重要、很有意义的问题。
针对直流微电网稳定性分析问题,以往大量研究是针对单源单负载、多源并联多负荷并联情形的稳定性,文献《Impedance-based local stability criterion for DCdistributed power systems》、《Impedance-based stability criterion for grid-connected inverters》、《Impedance specifications for stable DC distributedpower systems》给出了很多有益结论,但忽略了系统节点间的阻抗。文献《直流微电网节点阻抗特性与系统稳定性分析》考虑了系统节点间阻抗,利用系统节点导纳求取系统任意节点阻抗,但只适用于放射型结构,对环网结构系统无法有效判定系统稳定性。对于环网结构系统,状态空间法可有效分析。文献《Reduced-order model and stability analysis oflow-voltage DC microgrid》通过建立环网直流微电网的完整状态空间模型,求取系统矩阵的特征根分析系统稳定性,但是需要系统各部分的详细模型,得到的系统矩阵也较复杂,每个矩阵元素对应的物理意义不清晰。文献《Modeling and analysis of harmonicstability in an AC power-electronics-based power system》将单元连接法应用于主从控制的交流微电网中分析系统稳定性,用具有明确物理意义、可通过阻抗分析仪测量得到的阻抗和导纳对各功率模块进行建模,但并未将系统节点导纳矩阵与其他部分阻抗分离开来,很难直观得到系统的一般化稳定性判据,不利于扩展。
发明内容
本发明要解决的技术问题是,克服现有技术中的不足,提供一种主从控制模式下直流微电网的稳定性判定方法,用于解决直流微电网的稳定性分析问题。
为解决技术问题,本发明的解决方案是:
提供一种主从控制模式下直流微电网的稳定性判定方法,包括下述步骤:
(1)利用单元连接法对各电源、恒功率负载的功率模块进行戴维南/诺顿等效,得到相应的阻抗/导纳,并获得主电源的阻抗矩阵Zv、从电源和恒功率负载的导纳矩阵Yc
(2)根据已知的各节点电阻性负载大小得到阻性负载导纳矩阵YL,根据节点间线缆参数得到系统节点导纳矩阵Ynet
(3)将主电源的阻抗矩阵Zv、从电源和恒功率负载的导纳矩阵Yc、阻性负载导纳矩阵YL和系统节点导纳矩阵Ynet带入稳定性判据T,如下式所示:
T=(Yc+Zv(YL+Ynet))-1 (1)
计算稳定性判据T的极点;如果稳定性判据T的极点全部在左半平面,则系统稳定;如果稳定性判据T有右半平面极点,则系统不稳定。
本发明中,所述主电源的阻抗矩阵Zv是通过下述方式计算获得:
其中,Zv1是主电源的输出阻抗,通过计算或者阻抗分析仪进行测量获得;n是系统节点数。
本发明中,所述从电源和恒功率负载的导纳矩阵Yc是通过下述方式计算获得:
其中,m是从电源个数,不失一般性,假设m个从电源是接在节点2至节点m+1;其余节点(m+2)至节点n所接为恒功率负载;Ycj(j=2~m+1)是第j节点从电源的输出导纳;Yck(k=m+2~n)是第k节点恒功率负载的输入导纳;通过计算或者阻抗分析仪进行测量获得。
本发明中,所述阻性负载导纳矩阵YL是通过下述方式计算获得:
其中,YLi(i=1~n)是第i节点上的电阻性负载导纳,通过计算或者阻抗分析仪进行测量获得。
本发明中,所述系统节点导纳矩阵Ynet是通过下述方式计算获得:
其中,系统节点导纳矩阵Ynet的各个元素yij(i,j=1~n)通过下述方式获得:
其中,rcablejk(j,k=1~n,j≠k)是节点j和节点k之间线缆采用集中等效模型中的电阻,Lcablejk(j,k=1~n,j≠k)是节点j和节点k之间线缆采用集中等效模型中的电感,Ccablejk(j,k=1~n,j≠k)是节点j和节点k之间线缆采用集中等效模型中的电容,平均分布在二边。
与现有技术相比,本发明的优点有:
(1)本发明提出的主从控制模式下直流微电网稳定性判定方法能兼容考虑系统节点间阻抗;稳定性判据对系统结构无特定要求,可适用于放射型系统,也适用于环形系统,适用范围广。
(2)该判定方法在各功率模块自身稳定的前提下,只与各功率模块的阻抗/导纳相关,相应的阻抗/导纳是有明确物理意义的量,分析系统级稳定性时不用再顾及单个模块的稳定性,简化计算;且与系统节点导纳矩阵与其他部分解开。
(3)该判定方法中的主电源阻抗矩阵Zv、从电源和恒功率负载的导纳矩阵Yc都是对角矩阵,易于扩展。
(4)使用本发明提出的主从控制模式下直流微电网稳定性判定方法,可有效分析含电压控制型主电源、电流控制型从电源、恒功率负载、系统节点间阻抗等复杂直流微电网的稳定性。
附图说明
图1为4节点直流微电网。
图2为直流微电网的等效电路。
具体实施方式
下面结合附图对本发明作进一步详细描述:
考虑节点间阻抗的4节点直流微电网如图1所示。系统中有1个主电源,2个从电源,1个恒功率负载;主电源所接的节点为节点1,2个从电源分别接在节点2和节点3,恒功率负载接在节点4;各节点间通过线缆连接。主电源采用恒电压控制,维持系统电压在额定点;从电源采用电流控制,控制从电源输入系统的电流大小;恒功率负载通过控制从系统中吸收恒定功率。其中v1~v4是4个节点的电压;i1~i4是4个功率模块向节点注入的电流,流出功率模块为正方向;icablejk(j,k=1~4,j≠k)是节点j与节点k之间的线缆电感电流,方向以附图1中所示为正方向,节点之间不相连时无电流,如节点1和节点4、节点2和节点3。
从阻抗分析的角度来看,即使各个电源及恒功率负载单独设计时控制性能很好,能保证各自的稳定性。但接入系统后将会产生复杂的相互作用,系统稳定性将会降低。本发明采用单元连接法对各电源、恒功率负载等功率模块进行戴维南/诺顿等效来分析系统稳定性。其中主电源采用戴维南等效,从电流源采用诺顿等效,恒功率负载也采用诺顿等效。将单元连接法应用于附图1所示的系统,得到如附图2所示的等效电路。其中,Zv1是主电源输出阻抗,Yc2、Yc3是从电源的输出导纳,Yc4是恒功率负载的输入导纳。YLi(i=1~4)是各个节点上的阻性负载导纳。线缆全部采用π型集中等效模型,rcable12、Lcable12、Ccable12分别是节点1和节点2之间线缆的电阻参数、电感参数、分布电容参数;rcable13、Lcable13、Ccable13分别是节点1和节点3之间线缆的电阻参数、电感参数、分布电容参数;rcable24、Lcable24、Ccable24分别是节点2和节点4之间线缆的电阻参数、电感参数、分布电容参数;rcable34、Lcable34、Ccable34分别是节点3和节点4之间线缆的电阻参数、电感参数、分布电容参数。v1ref是主电源的参考电压,i2ref、i3ref是从电源的参考电流,v4ref是恒功率负载的参考电压。Gi(i=1~4)是各个功率模块参考到输出的闭环传递函数。
通过计算或者阻抗分析仪对节点1处的主电源进行测量获得主电源的输出阻抗Zv1,构建主电源的阻抗矩阵Zv
通过计算或者阻抗分析仪获得节点2、节点3的从电源输出导纳Yc2、Yc3,通过计算或者阻抗分析仪获得节点4的恒功率负载的输入导纳Yc4,构建从电源和恒功率负载的导纳矩阵Yc
通过计算或者阻抗分析仪进行测量获得各节点上的电阻性负载导纳YLi(i=1~4),构建阻性负载导纳矩阵YL
构建系统节点导纳矩阵Ynet,如下式所示:
其中,系统节点导纳矩阵Ynet的各个元素yij(i,j=1~4)通过下述方式获得:
将上述得到的主电源的阻抗矩阵Zv、从电源和恒功率负载的导纳矩阵Yc、阻性负载导纳矩阵YL、系统节点导纳矩阵Ynet带入稳定性判据T=(Yc+Zv(YL+Ynet))-1。计算稳定性判据T的极点,如果稳定性判据T的极点全部在左半平面,则系统稳定;如果稳定性判据T有右半平面极点,则系统不稳定。

Claims (5)

1.主从控制模式下直流微电网的稳定性判定方法,其特征在于,包括下述步骤:
(1)利用单元连接法对各电源、恒功率负载的功率模块进行戴维南/诺顿等效,得到相应的阻抗/导纳,并获得主电源的阻抗矩阵Zv、从电源和恒功率负载的导纳矩阵Yc
(2)根据已知的各节点电阻性负载大小得到阻性负载导纳矩阵YL,根据节点间线缆参数得到系统节点导纳矩阵Ynet
(3)将主电源的阻抗矩阵Zv、从电源和恒功率负载的导纳矩阵Yc、阻性负载导纳矩阵YL和系统节点导纳矩阵Ynet带入稳定性判据T,如下式所示:
T=(Yc+Zv(YL+Ynet))-1 (1)
计算稳定性判据T的极点;如果稳定性判据T的极点全部在左半平面,则系统稳定;如果稳定性判据T有右半平面极点,则系统不稳定。
2.根据权利要求1所述的方法,其特征在于,所述主电源的阻抗矩阵Zv是通过下述方式计算获得:
其中,Zv1是主电源的输出阻抗,通过计算或者阻抗分析仪进行测量获得;n是系统节点数。
3.根据权利要求1所述的方法,其特征在于,所述从电源和恒功率负载的导纳矩阵Yc是通过下述方式计算获得:
其中,m是从电源个数,不失一般性,假设m个从电源是接在节点2至节点m+1;其余节点(m+2)至节点n所接为恒功率负载;Ycj(j=2~m+1)是第j节点从电源的输出导纳;Yck(k=m+2~n)是第k节点恒功率负载的输入导纳;通过计算或者阻抗分析仪进行测量获得。
4.根据权利要求1所述的方法,其特征在于,所述阻性负载导纳矩阵YL是通过下述方式计算获得:
其中,YLi(i=1~n)是第i节点上的电阻性负载导纳,通过计算或者阻抗分析仪进行测量获得。
5.根据权利要求1所述的方法,其特征在于,所述系统节点导纳矩阵Ynet是通过下述方式计算获得:
其中,系统节点导纳矩阵Ynet的各个元素yij(i,j=1~n)通过下述方式获得:
y i i = Σ j ≠ i ( C c a b l e i j s / 2 + 1 / ( r c a b l e i j + L c a b l e i j s ) ) ; i = 1 ~ n y j k = - 1 / ( r c a b l e j k + L c a b l e j k s ) ; j , k = 1 ~ n , j ≠ k r c a b l e j k = r c a b l e k j ; L c a b l e j k = L c a b l e k j ; C c a b l e j k = C c a b l e k j ; j , k = 1 ~ n , j ≠ k - - - ( 6 )
其中,rcablejk(j,k=1~n,j≠k)是节点j和节点k之间线缆采用集中等效模型中的电阻,Lcablejk(j,k=1~n,j≠k)是节点j和节点k之间线缆采用集中等效模型中的电感,Ccablejk(j,k=1~n,j≠k)是节点j和节点k之间线缆采用集中等效模型中的电容,平均分布在二边。
CN201710005641.XA 2017-01-04 2017-01-04 主从控制模式下直流微电网的稳定性判定方法 Active CN106877309B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710005641.XA CN106877309B (zh) 2017-01-04 2017-01-04 主从控制模式下直流微电网的稳定性判定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710005641.XA CN106877309B (zh) 2017-01-04 2017-01-04 主从控制模式下直流微电网的稳定性判定方法

Publications (2)

Publication Number Publication Date
CN106877309A true CN106877309A (zh) 2017-06-20
CN106877309B CN106877309B (zh) 2019-05-17

Family

ID=59165532

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710005641.XA Active CN106877309B (zh) 2017-01-04 2017-01-04 主从控制模式下直流微电网的稳定性判定方法

Country Status (1)

Country Link
CN (1) CN106877309B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107544249A (zh) * 2017-09-14 2018-01-05 浙江大学 对等控制模式下直流分布式系统的稳定性判定方法
CN107947237A (zh) * 2017-12-14 2018-04-20 合肥工业大学 一种多类型逆变器孤岛微网稳定性分析方法
CN112034251A (zh) * 2020-08-28 2020-12-04 武汉大学 评价诺顿型间谐波源接入点激励并联谐振能力的方法
CN112217191A (zh) * 2019-07-11 2021-01-12 华北电力大学(保定) 一种基于节点阻抗矩阵的直流配电网稳定性分析方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105305402A (zh) * 2015-10-20 2016-02-03 天津大学 一种直流微网母线电压鲁棒自治控制方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105305402A (zh) * 2015-10-20 2016-02-03 天津大学 一种直流微网母线电压鲁棒自治控制方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
刘晓东等: "直流微电网节点阻抗特性与系统稳定性分析", 《电网技术》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107544249A (zh) * 2017-09-14 2018-01-05 浙江大学 对等控制模式下直流分布式系统的稳定性判定方法
CN107947237A (zh) * 2017-12-14 2018-04-20 合肥工业大学 一种多类型逆变器孤岛微网稳定性分析方法
CN112217191A (zh) * 2019-07-11 2021-01-12 华北电力大学(保定) 一种基于节点阻抗矩阵的直流配电网稳定性分析方法
CN112217191B (zh) * 2019-07-11 2023-07-28 华北电力大学(保定) 一种基于节点阻抗矩阵的直流配电网稳定性分析方法
CN112034251A (zh) * 2020-08-28 2020-12-04 武汉大学 评价诺顿型间谐波源接入点激励并联谐振能力的方法
CN112034251B (zh) * 2020-08-28 2021-11-23 武汉大学 评价诺顿型间谐波源接入点激励并联谐振能力的方法

Also Published As

Publication number Publication date
CN106877309B (zh) 2019-05-17

Similar Documents

Publication Publication Date Title
Xu et al. Robust dispatch of high wind power-penetrated power systems against transient instability
Beerten et al. VSC MTDC systems with a distributed DC voltage control-A power flow approach
CN103296677B (zh) 一种在线大电网恢复辅助决策系统
CN103973203B (zh) 一种适用于安全稳定分析的大型光伏电站在线等值建模方法
CN103825267B (zh) 一种mmc-mtdc直流侧短路电流的计算方法
CN103049617B (zh) 保留无源性的大规模配电网络电磁暂态仿真模型化简方法
CN101719182A (zh) 一种交直流电力系统分割并行电磁暂态数字仿真方法
CN103632313A (zh) 一种基于pmu数据的电网动态可观方法
CN103077480B (zh) 一种电力系统的安全校核方法
CN106877309A (zh) 主从控制模式下直流微电网的稳定性判定方法
CN105140907A (zh) 直流微网多智能体自适应下垂一致性协调控制方法及装置
CN103595040A (zh) 含直流分布式电源的配电网及其综合负荷特性的分析方法
Sanchez et al. Stability evaluation of a DC micro-grid and future interconnection to an AC system
CN103345162A (zh) 一种功率级数模混合仿真系统
CN106786529B (zh) 一种分布式静态安全分析方法
CN103795057A (zh) 配电网潮流计算中基于搜索的配电网拓扑编号生成方法
CN110797874A (zh) 含电力电子变压器的交直流混合配电网状态估计方法
CN109494719A (zh) 一种中低压混合配电网层次化阻抗分析方法
CN105281325A (zh) 一种提升受端电网外受电能力的网架优化方法
Feng et al. Graph computing based distributed parallel power flow for AC/DC systems with improved initial estimate
CN104283208B (zh) 大型电网概率性可用输电能力的分解协调计算方法
CN103390890B (zh) 基于电流分布系数的配电网潮流分析方法
Tao et al. Power connection technology for full-digital real-time simulation and analog simulator
Chen et al. An improved load-shedding model based on power flow tracing
CN113315100A (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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20220518

Address after: 311201 No. 678, chenduan Road, Xintang street, Xiaoshan District, Hangzhou City, Zhejiang Province

Patentee after: ZHEJIANG ZHONGXIN ELECTRIC POWER ENGINEERING CONSTRUCTION Co.,Ltd. AUTOMATION BRANCH

Patentee after: ZHEJIANG ZHONGXIN ELECTRIC POWER ENGINEERING CONSTRUCTION Co.,Ltd.

Patentee after: STATE GRID ZHEJIANG XIAOSHAN DISTRICT POWER SUPPLY Co.,Ltd.

Address before: 310058 Yuhang Tang Road, Xihu District, Hangzhou, Zhejiang 866

Patentee before: ZHEJIANG University