CN107294085B - 基于临界特征根跟踪的微电网延时裕度计算方法 - Google Patents
基于临界特征根跟踪的微电网延时裕度计算方法 Download PDFInfo
- Publication number
- CN107294085B CN107294085B CN201710456420.4A CN201710456420A CN107294085B CN 107294085 B CN107294085 B CN 107294085B CN 201710456420 A CN201710456420 A CN 201710456420A CN 107294085 B CN107294085 B CN 107294085B
- Authority
- CN
- China
- Prior art keywords
- distributed power
- power supply
- small signal
- voltage
- delay
- 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
- 238000004364 calculation method Methods 0.000 title claims abstract description 15
- 230000006854 communication Effects 0.000 claims abstract description 41
- 238000004891 communication Methods 0.000 claims abstract description 40
- 238000000034 method Methods 0.000 claims abstract description 22
- 230000003068 static effect Effects 0.000 claims abstract description 10
- 239000011159 matrix material Substances 0.000 claims description 47
- 230000008859 change Effects 0.000 claims description 37
- 230000001934 delay Effects 0.000 claims description 10
- 230000005611 electricity Effects 0.000 claims 1
- 238000013461 design Methods 0.000 abstract description 6
- 230000000087 stabilizing effect Effects 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 15
- 230000009471 action Effects 0.000 description 12
- 230000000694 effects Effects 0.000 description 9
- 150000001875 compounds Chemical class 0.000 description 8
- 230000010355 oscillation Effects 0.000 description 8
- 239000011541 reaction mixture Substances 0.000 description 6
- 208000000419 Chronic Hepatitis B Diseases 0.000 description 5
- 208000002672 hepatitis B Diseases 0.000 description 5
- 238000011084 recovery Methods 0.000 description 5
- 238000004088 simulation Methods 0.000 description 5
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000009827 uniform distribution Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- 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/12—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
- H02J3/16—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load by adjustment of reactive power
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- 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/001—Methods to deal with contingencies, e.g. abnormalities, faults or failures
-
- 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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/381—Dispersed generators
-
- 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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
- H02J3/50—Controlling the sharing of the out-of-phase component
-
- 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
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- 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
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/30—Reactive power compensation
-
- 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
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P80/00—Climate change mitigation technologies for sector-wide applications
- Y02P80/10—Efficient use of energy, e.g. using compressed air or pressurized fluid as energy carrier
- Y02P80/14—District level solutions, i.e. local energy networks
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Business, Economics & Management (AREA)
- Power Engineering (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- Health & Medical Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Human Resources & Organizations (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了基于临界特征根跟踪的微电网延时裕度计算方法,属于微电网运行控制的技术领域。本发明基于静态输出反馈建立包含通讯延时电压反馈控制量的微电网闭环小信号模型,从而获取含有超越项的特征方程,对***特征方程的超越项进行临界特征根轨迹跟踪,搜寻可能的纯虚特征根进而计算使微电网稳定的最大延时时间,对控制器参数与延时裕度间的关系进行研究,从而指导控制参数的设计,有效提高微网稳定性及动态性能。
Description
技术领域
本发明公开了基于临界特征根跟踪的微电网延时裕度计算方法,尤其涉及一种微电网二次电压控制延时裕度的计算方法,属于微电网运行控制的技术领域。
背景技术
随着地球资源的日渐衰竭以及人们对环境问题的关注,可再生能源的接入越来越受到世界各国的重视。微电网是一种在能量供应***中增加可再生能源和分布式能源渗透率的新兴能量传输模式,其组成部分包括微型燃气轮机、风力发电机、光伏、燃料电池、储能设备等不同种类的分布式能源(Distributed Energy Resources,DER)、各种电负荷和/或热负荷的用户终端以及相关的监控保护装置。
微电网内部的电源主要由电力电子器件转换能量并提供必须的控制。微电网相对于外部大电网表现为单一的受控单元,可同时满足用户对电能质量和供电安全等要求。微电网与大电网之间通过公共连接点进行能量交换,双方互为备用,从而提高了供电的可靠性。微电网是规模较小的分散***且负荷的距离较近,在增加本地供电可靠性的同时降低了网损,这大大增加了能源利用效率,因此微电网是一种符合未来智能电网发展要求的新型供电模式。
下垂控制因可以实现无通讯的功率均分而受到关注,但各分布式电源输出电压会出现稳态偏差,同时,由于各分布式电源输出阻抗不同,无功功率均分很难达到满意效果,因此,需要采用微电网二次电压控制以提高无功均分效果及电压性能。目前,设计的协同电压控制为集中式控制结构,微电网集中式电压控制器产生控制信号并下发至各分布式电源本地控制器,该集中式控制结构依赖于通讯技术,但是通讯过程通常受到信息延时、数据丢包的影响,信息延时、数据丢包等影响导致微电网动态性能不佳甚至危及***稳定性。基于以上原因,有必要研究一套微电网二次电压控制延时裕度计算方法,分析使微电网稳定的最大通讯延时时间,有必要对微网集中控制器参数与延时裕度的关系进行分析,从而指导控制参数的设计,有效提高微网稳定性及动态性能。
发明内容
本发明的发明目的是针对在微电网无功功率均分和电压恢复控制中通常忽略通讯延时对动态性能影响的现象,充分考虑了电力电子接口型微电网惯性小从而导致通讯延时对***稳定性不可忽视的实际情况,提供了基于临界特征根跟踪的微电网延时裕度计算方法,通过求取微电网特征方程的所有可能纯虚特征根进而计算使微电网稳定的最大延时时间,通过对控制器参数与稳定裕度间的关系进行研究为控制参数的设计提供指导意见,解决了现有微电网***的稳定性受通讯技术影响的技术问题。
本发明为实现上述发明目的采用如下技术方案:
基于临界特征根跟踪的微电网延时裕度计算方法,根据静态反馈输出建立包含通讯延时电压反馈控制量的逆变器闭环小信号模型及分布式电源闭环小信号模型,结合连接网络、负载阻抗的动态方程及分布式电源闭环小信号模型建立微电网小信号模型,从微电网小信号模型获取含有超越项的特征方程,对超越项进行临界特征根轨迹跟踪进而确定满足***稳定性要求的延时裕度。
进一步地,基于临界特征根跟踪的微电网延时裕度计算方法中,根据静态反馈输出建立的包含通讯延时电压反馈控制量的逆变器闭环小信号模型为:Δxinv、分别为逆变器的闭环小信号状态变量及其变化率,Δxinv1、Δxinv2、Δxinvi、Δxinvn分别为第1个、第2个、第i个、第n个分布式电源的小信号状态变量,分别为第1个、第2个、第i个、第n个分布式电源的无功功率辅助小信号状态变量,第i个分布式电源的无功功率辅助小信号状态变量由表达式:确定,为第i个分布式电源无功功率辅助小信号状态变量的变化率,Qi为第i个分布式电源实际输出的无功功率,nQi为第i个分布式电源的电压下垂特性系数,n为分布式电源的数目,Δγ为分布式电源的电压辅助小信号状态变量,分布式电源的电压辅助小信号状态变量Δγ由表达式:确定,为分布式电源的电压辅助小信号状态变量的变化率,为第i个分布式电源平均电压的期望值,Vodi为在第i个分布式电源输出电压在其自身参考坐标系dq下的d轴分量,Ainv为分布式电源的状态矩阵,ΔVbDQ为母线电压在公共参考坐标系DQ中的小信号状态变量,ΔVbDQ=[ΔVbDQ1,ΔVbDQ2,…,ΔVbDQl,…,ΔVbDQm]T,ΔVbDQ1、ΔVbDQ2、ΔVbDQl、ΔVbDQm分别为第1根、第2根、第l根、第m根母线的电压在公共参考坐标系DQ中的小信号状态变量,m为母线的数目,Binv为分布式电源对母线电压的输入矩阵,Δu为分布式电源的二次电压小信号控制量,Δu=[Δu1,Δu2,…,Δui,…,Δun]T,Δu1、Δu2、Δui、Δun分别为第1个、第2个、第i个、第n个分布式电源的二次电压小信号控制量,Bu为分布式电源对二次电压小信号控制量的输入矩阵,Δui=KQiΔyinvQi(t-τi)+KViΔyinvV(t-τi),t为当前时刻,τi为第i个分布式电源本地控制器与微网二次电压集中控制器间的通讯时延,KQi、KVi分别为第i个分布式电源的无功功率控制系数、电压控制系数,ΔyinvQi为第i个分布式电源的无功功率输出小信号状态变量,ΔyinvQ、ΔyinvV分别为分布式电源的无功功率输出小信号状态变量、电压输出小信号状态变量,CinvQ、CinvV分别为分布式电源的无功功率输出矩阵、电压输出矩阵。
再进一步地,基于临界特征根跟踪的微电网延时裕度计算方法中,根据静态反馈输出建立的包含通讯延时电压反馈控制量的分布式电源闭环小信号模型为: 为第i个分布式电源的延时状态矩阵,Bui为第i个分布式电源对二次电压小信号控制量的输入矩阵,CinvQi为第i个分布式电源的无功功率输出矩阵,ΔioDQ为公共参考坐标系DQ中分布式电源输出电流的小信号状态变量,Cinvc为分布式电源的电流输出矩阵。
更进一步地,基于临界特征根跟踪的微电网延时裕度计算方法中,微电网小信号模型为x、分别为微电网小信号状态变量及其变化率,x=[ΔxinvΔilineDQΔiloadDQ]T,ΔilineDQ为公共参考坐标系DQ中分布式电源所连接母线间的连接线路的电流的小信号状态变量,公共参考坐标系DQ中第i个分布式电源所连接母线和第j个分布式电源所连接母线之间的连接线路ij的电流的小信号状态变量为:ΔilineDij、分别为连接线路ij的电流在公共参考坐标系DQ下的D轴小信号分量及其变化率,ΔilineQij、分别为连接线路ij的电流在公共参考坐标系DQ下的Q轴小信号分量及其变化率,rlineij、Llineij分别为连接线路ij的线路电阻和线路电感,ω0为微电网额定角频率,ΔVbusDi、ΔVbusQi分别为第i个分布式电源所连接母线的电压在公共参考坐标系DQ下的D轴分量、Q轴分量,ΔVbusDj、ΔVbusQj分别为第j个分布式电源所连接母线的电压在公共参考坐标系DQ下的D轴分量、Q轴分量,ΔiloadDQ为公共参考坐标系DQ中母线所连接负载的电流的小信号状态变量,公共参考坐标系DQ中第l根母所连接负载的电流的小信号状态变量为:ΔiloadDl、分别为第l根母线所连接负载的电流在公共参考坐标系DQ下的D轴分量及其变化率,ΔiloadQl、分别为第l根母线所连接负载的电流在公共参考坐标系DQ下的Q轴分量及其变化率,Rloadl、Lloadl分别为第l根母线所连接负载的负载电阻、负载电感,ΔVbusDl、ΔVbusQl分别为第l根母线的电压在公共参考坐标系DQ下的D轴分量、Q轴分量,Adi、τi分别为第i个分布式电源的延时状态矩阵和延时。
作为基于临界特征根跟踪的微电网延时裕度计算方法的进一步优化方案,从微电网小信号模型获取含有超越项的特征方程的方法为:在分布式电源的延时一致时得到微电网小信号模型的的特征方程:CEτ(s,τ)=det(sI-A-Ade-τs),s为时域复平面参数,τ为各分布式电源的一致时延时间,CEτ(·)表示各分布式电源一致时延τ时得到的微电网小信号模型的特征方程,det(·)为矩阵行列式,I为单位矩阵,Ad为分布式电源的延时状态矩阵,e-τs为超越项。
作为基于临界特征根跟踪的微电网延时裕度计算方法的再进一步优化方案,对超越项进行临界特征根轨迹跟踪进而确定满足***稳定性要求的延时裕度,具体方法为:以延时时间辅助变量作为特征方程的变量,求解特征方程在延时时间辅助变量变化周期内的所有纯虚特征根,从所有纯虚特征根对应的临界延时时间中选取最小值作为满足***稳定性要求的延时裕度,所述延时时间辅助变量为分布式电源延时和虚特征根幅值的乘积。
本发明采用上述技术方案,具有以下有益效果:
(1)本发明提出了一种微电网二次电压控制延时裕度的计算方法,该方法基于静态输出反馈建立包含通讯延时电压反馈控制量的微电网闭环小信号模型,从而获取含有超越项的特征方程,对***特征方程的超越项进行临界特征根轨迹跟踪,搜寻可能的纯虚特征根进而计算使微电网稳定的最大延时时间,该方法能够有效降低通讯延时对微网动态性能的影响,有效提高微网稳定性及动态性能;
(2)通过对不同控制器参数下的***稳定裕度进行求取,对控制器参数与延时裕度间的关系进行研究,从而指导控制参数的设计,有效提高微网稳定性及动态性能。
附图说明
图1是本发明实施例的流程图;
图2是本发明实施例中微电网一次、二次控制框图;
图3是本发明实施例中采用的微电网仿真***图;
图4是在某一组控制参数kIQ=0.02,kIV=20下,临界特征根轨迹跟踪示意图;
图5是本发明实施例中控制器参数与***延时裕度的关系;
图6(a)是本发明实例在某一组控制参数kIQ=0.02,kIV=20下,3种不同通讯延时对平均电压动态性能的影响;
图6(b)是本发明实例在某一组控制参数kIQ=0.02,kIV=20下,3种不同通讯延时对分布式电源1无功功率动态性能的影响;
图6(c)是本发明实例在某一组控制参数kIQ=0.02,kIV=20下,3种不同通讯延时对分布式电源2无功功率动态性能的影响;
图7(a)是本发明实例在某一组控制参数kIQ=0.04,kIV=40下,3种不同通讯延时对平均电压动态性能的影响;
图7(b)是本发明实例在某一组控制参数kIQ=0.04,kIV=40下,3种不同通讯延时对分布式电源1无功功率动态性能的影响;
图7(c)是本发明实例在某一组控制参数kIQ=0.04,kIV=40下,3种不同通讯延时对分布式电源2无功功率动态性能的影响。
具体实施方式
下面结合附图对发明的技术方案进行详细说明。
如图1所示,本发明公开的基于临界特征根跟踪的微电网延时裕度计算方法,包括下述步骤:
步骤10)基于静态输出反馈建立包含通讯时延电压反馈控制量的逆变器闭环小信号模型
各分布式电源通过本地控制器中的下垂控制环设置逆变器输出电压及频率参考指令,如式(1)所示:
式(1)中,ωi表示第i个分布式电源的本地角频率;ωn表示分布式电源本地角频率的参考值,单位:弧度/秒;mPi表示第i个分布式电源的频率下垂特性系数,单位:弧度/秒·瓦;P表示第i个分布式电源实际输出的有功功率,单位:瓦;kVi表示第i个分布式电源的下垂控制增益;表示第i个分布式电源输出电压的变化率,单位:伏/秒;Vn表示分布式电源输出电压的参考值,单位:伏;Vo,magi表示第i个分布式电源实际输出的电压,单位:伏;nQi表示第i个分布式电源的电压下垂特性系数,单位:伏/乏;Qi表示第i个分布式电源实际输出的无功功率,单位:乏。
第i个分布式电源实际输出的有功功率Pi、无功功率Qi通过低通滤波器获得,如式(2)所示:
式(2)中,表示第i个分布式电源实际输出有功功率的变化率,单位:瓦/秒;ωci表示第i个分布式电源所连接低通滤波器的剪切频率,单位:弧度/秒;Vodi表示在第i个分布式电源的dq参考坐标系中,第i个分布式电源输出电压的d轴分量,单位:伏;Voqi表示在第i个分布式电源的dq参考坐标系中,第i个分布式电源输出电压的q轴分量,单位:伏;iodi表示在第i个分布式电源的dq参考坐标系中,第i个分布式电源输出电流的d轴分量,单位:安;ioqi表示在第i个分布式电源的dq参考坐标系中,第i个分布式电源输出电压的q轴分量,单位:安;表示第i个分布式电源实际输出无功功率的变化率,单位:乏/秒。
微电网一次、二次控制框图如图2所示,各分布式电源一次控制通过锁相环控制使输出电压q轴分量为0,基于分布式电源电压的二次控制得到式(3):
式(3)中,表示在第i个分布式电源的dq参考坐标系下,第i个分布式电源输出电压的d轴分量的变化率,单位:伏/秒;Vni表示第i个分布式电源输出电压的参考值,ui表示二次电压控制量,单位:伏。
分布式电源输出电流的动态方程如式(4)所示:
式(4)中,表示在第i个分布式电源的dq参考坐标系中,第i个分布式电源输出电流的d轴分量的变化率,单位:安/秒;Rci表示第i个分布式电源至其所连接母线的连接电阻,单位:欧姆;Lci表示第i个分布式电源至其所连接母线的连接电感,单位:亨利;Vbusdi表示在第i个分布式电源的dq参考坐标系中,第i个分布式电源所连接母线的电压d轴分量;表示在第i个分布式电源的dq参考坐标系中,第i个分布式电源输出电流的q轴分量的变化率,单位:安/秒;Vbusqi表示在第i个分布式电源的dq参考坐标系中,第i个分布式电源所连接母线的电压q轴分量,单位:伏。
各分布式电源基于本地的dq参考坐标系建立模型,为建立含多个分布式电源的微电网整体模型,设定其中一个分布式电源的dq参考坐标系为公共参考坐标系DQ,则其它分布式电源dq参考坐标系下的输出电流需要转换到公共参考坐标系DQ下,转换方程如式(5)所示:
式(5)中,ioDi表示在公共参考坐标系DQ中,第i个分布式电源输出电流在D轴的分量,ioQi表示在公共参考坐标系DQ中,第i个分布式电源输出电流在Q轴的分量,单位:安;Ti表示第i个分布式电源输出电流从第i个分布式电源dq参考坐标系到公共参考坐标系DQ的转换矩阵,δi表示第i个分布式电源dq参考坐标系旋转角度与公共参考坐标系DQ旋转角度之间的静态差值,单位:度,δi可以由式(6)求得:
式(6)中,ωcom表示公共参考坐标系DQ的角频率;表示δi的变化率。
线性化式(1)~式(6)得到如式(7)所示的第i个分布式电源的开环小信号模型:
式(7)中,表示第i个分布式电源的小信号状态变量的变化率, Δxinvi表示第i个分布式电源的小信号状态变量,Δxinvi=[Δδi,ΔPi,ΔQi,ΔVodi,Δiodi,Δioqi]T;ΔVbDQi表示在公共参考坐标系DQ中第i个分布式电源所连接母线的电压的小信号状态变量;ΔVsDQi=[ΔVbDi,ΔVbQi]T,ΔVbDi表示在公共参考坐标系DQ中第i个分布式电源所连接母线的电压在D轴的小信号分量,ΔVbQi表示在公共参考坐标系DQ中第i个分布式电源所连接母线的电压在Q轴的小信号分量,单位:伏;Δωcom表示公共参考坐标系DQ角频率的小信号状态变量,单位:弧度/秒;Δui表示第i个分布式电源二次电压的小信号控制量,单位:伏;Ainvi表示第i个分布式电源的状态矩阵;Binvi表示第i个分布式电源对其所连接母线电压的输入矩阵;Biwcom表示第i个分布式电源对公共参考坐标系角频率的输入矩阵;Bui表示第i个分布式电源对其二次电压小信号控制量的输入矩阵;ΔioDQi表示在公共参考坐标系DQ中,第i个分布式电源输出电流的小信号状态变量,ΔioDQi=[ΔioDi,ΔioQi]T,单位:安;Cinvci表示第i个分布式电源的电流输出矩阵。
根据式(7),ΔVbusDQi和Δωcom作为第i个分布式电源的扰动变量,其中一般选取第1个分布式电源的参考坐标系作为公共参考坐标系DQ,则,
Δωcom=[0 -mP1 0 0 0 0]Δxinv1 式(8),
式(8)中,mP1表示第1个分布式电源的频率下垂特性系数,单位:弧度/秒·瓦;Δxinv1表示第1个分布式电源的小信号状态变量,Δxinv1=[Δδ1,ΔP1,ΔQ1,ΔVod1,Δiod1,Δioq1]T。
根据式(7)和式(8),可以得到n个分布式电源所组成***的小信号模型:
式(9)中,Δxinv1表示第1个分布式电源的小信号状态变量,Δxinv2表示第2个分布式电源的小信号状态变量,Δxinvn表示第n个分布式电源的小信号状态变量;ΔVbDQ=[ΔVbDQ1ΔVbDQ2...ΔVbusDQm]T,ΔVbDQ1=[ΔVbD1ΔVbQ1]T,ΔVbD1表示在公共参考坐标系DQ中母线1的电压在D轴的小信号分量,ΔVbQ1表示在公共参考坐标系DQ中母线1的电压在Q轴的小信号分量,ΔVbDQ2=[ΔVbD2ΔVbQ2]T,ΔVbD2表示在公共参考坐标系DQ中母线2的电压在D轴的小信号分量,ΔVbQ2表示在公共参考坐标系DQ中母线2的电压在Q轴的小信号分量,ΔVbDQm=[ΔVbDmΔVbQm]T,ΔVbDm表示在公共参考坐标系DQ中母线m的电压在D轴的小信号分量,ΔVbQm表示在公共参考坐标系DQ中母线m的电压在Q轴的小信号分量;Δu=[Δu1Δu2....Δun]T,Δu1表示分布式电源1的二次电压小信号控制量,Δu2表示分布式电源2的二次电压小信号控制量,Δun表示分布式电源n的二次电压小信号控制量;ΔioDQ=[ΔioDQ1ΔioDQ2...ΔioDQn]T,ΔioDQ1=[ΔioD1,ΔioQ1]T,ΔioD1表示在公共参考坐标系DQ中第1个分布式电源输出电流在D轴的小信号分量,ΔioQ1表示在公共参考坐标系DQ中第i个分布式电源输出电流在Q轴的小信号分量,ΔioDQ2=[ΔioD2,ΔioQ2]T,ΔioD2表示在公共参考坐标系DQ中第2个分布式电源输出电流在D轴的小信号分量,ΔioQ2表示在公共参考坐标系DQ中第2个分布式电源输出电流在Q轴的小信号分量;ΔioDQn=[ΔioDn,ΔioQn]T,ΔioDn表示在公共参考坐标系DQ中第n个分布式电源输出电流在D轴的小信号分量,ΔioQn表示在公共参考坐标系DQ中第n个分布式电源输出电流在Q轴的小信号分量,为n个分布式电源的状态矩阵;为n个分布式电源对母线电压的输入矩阵;为n个分布式电源对二次电压小信号控制量的输入矩阵;为n个分布式电源的电流输出矩阵。
本发明基于无功功率均分和电压恢复的控制要求实现微电网电压控制。无功功率均分是指各分布式电源输出无功功率按功率容量进行分配,电压恢复指各分布式电源输出电压平均值恢复至额定值,首先定义如下动态方程:
式(10)中,为第i个分布式电源的无功功率辅助小信号状态变量的变化率,单位:乏;为第i个分布式电源期望输出的无功功率,单位:乏;nQi表示第i个分布式电源的电压下垂特性系数,单位:伏/乏;为分布式电源的电压辅助小信号状态变量的变化率,单位:伏;为各分布式电源的平均输出电压,为第i个分布式电源平均电压的期望值,单位:伏。
因此,基于输出反馈的逆变器闭环小信号模型为:
式(11)中,Δxinv表示n个逆变器的闭环小信号状态变量, 为第1个分布式电源的无功功率辅助小信号状态变量,为第2个分布式电源的无功功率辅助小信号状态变量,为第i个分布式电源的无功功率辅助小信号状态变量,为第n个分布式电源的无功功率辅助小信号状态变量,Δγ为各分布式电源的电压辅助小信号状态变量;ΔyinvQ为无功功率输出小信号状态变量 为第1个分布式电源的无功功率辅助小信号状态变量的变化率,为第2个分布式电源的无功功率辅助小信号状态变量的变化率,为第n个分布式电源的无功功率辅助小信号状态变量的变化率;ΔyinvV为分布式电源的电压输出小信号状态变量, 为各分布式电源的电压辅助小信号状态变量的变化率;CinvQ表示各分布式电源的无功功率输出矩阵;CinvV表示各分布式电源的电压输出矩阵。
定义分布式电源控制量为:
式(12)中,δQi表示第i个分布式电源的无功功率控制信号;kPQ表示无功功率比例积分控制器中的比例项系数;kIQ表示无功功率比例积分控制器中的积分项系数;δVi表示第i个分布式电源的平均电压恢复控制信号;kPV表示平均电压比例积分控制器中的比例项系数;kIV表示平均电压比例积分控制器中的积分项系数。
当微网电压集中控制器与各分布式电源间存在通讯延时时,电压控制量为:
Δui=ΔδQi(t-τi)+ΔδVi(t-τi)=KQiΔyinvQi(t-τi)+KViΔyinvV(t-τi) 式(13),
式(13)中,τi为第i个分布式电源本地控制器与微网二次电压集中控制器间的通讯时延,单位:秒;KQi表示第i个分布式电源的无功功率控制器,KQi=[kPQikIQi];KVi表示第i个分布式电源的电压控制器,KVi=[kPVi kIVi]。
结合式(11)~式(13),得到n个分布式电源的闭环小信号模型为:
式(14)中,为第i个分布式电源的延时状态矩阵,Bui为第i个分布式电源对二次电压小信号控制量的输入矩阵,CinvQi为第i个分布式电源的无功功率输出矩阵,Cinvc为分布式电源的电流输出矩阵。
步骤20)结合连接网络、负载型阻抗的动态方程,建立微电网小信号模型
公共参考坐标系DQ中第i个分布式电源所连接母线和第j个分布式电源所连接母线之间的连接线路ij的电流小信号动态方程如式(15)所示:
式(15)中,表示在公共参考坐标系DQ中,第ij条连接线路电流D轴小信号分量的变化率,单位:安/秒;rlineij表示第ij条连接线路的线路电阻,单位:欧姆;Llineij表示第ij条连接线路的线路电感,单位:亨利;ΔilineDij表示在公共参考坐标系DQ中,第ij条连接线路电流的D轴小信号分量,ΔilineQij表示在公共参考坐标系DQ中,第ij条连接线路的电流的Q轴小信号分量,单位:安;ω0表示微网额定角频率,单位:弧度/秒;ΔVbusDi表示在公共参考坐标系DQ中,第i个分布式电源所连接母线的电压在D轴的小信号分量;ΔVbusDj表示在公共参考坐标系DQ中,第j个分布式电源所连接母线的电压在D轴的小信号分量;表示在公共参考坐标系DQ中,第ij条连接线路电流的Q轴小信号分量的变化率,单位:安/秒;ΔVbusQi表示在公共参考坐标系DQ中,第i个分布式电源所连接母线的电压在Q轴的小信号分量,ΔVbusQj表示在公共参考坐标系DQ中,第j个分布式电源所连接母线的电压在Q轴的小信号分量,单位:伏。
公共参考坐标系DQ中第l根母所连接负载的电流动态方程,如式(16)所示:
式(16)中,表示在公共参考坐标系DQ中,第l根母线所连接负载的电流在D轴的小信号分量变化率,单位:安/秒;Rloadl表示第l根母线所连接负载的负载电阻,单位:欧姆;Lloadl表示第l根母线所连接负载的负载电感,单位:亨利;ΔiloadDl为在公共参考坐标系DQ中,第l根母线所连接负载的电流在D轴的小信号分量,ΔiloadQl为在公共参考坐标系DQ中,第l根母线所连接负载的电流在Q轴的小信号分量,单位:安;表示在公共参考坐标系DQ中,第l根母线所连接负载的电流在Q轴的小信号分量变化率,单位:安/秒。
设定连接于第i个分布式电源所连接母线和第j个分布式电源所连接母线之间的连接线路的小信号方程如式(17)所示:
式(17),
式(17)中,Rloadj、Lloadj分别为第j个分布式电源所连接母线上负载的阻值和电感值;ΔioDj、ΔioQj分别为第j个分布式电源输出电流在公共参考坐标系DQ中的D轴小信号分量和Q轴小信号分量。
将式(17)代入式(14)~式(16),可得包含n个分布式电源、s条支路、p个负载的微电网小信号模型为:
式(18)中,x为微电网小信号状态变量,x=[ΔxinvΔilineDQΔiloadDQ]T,ΔilineDQ为公共参考坐标系DQ中分布式电源所连接母线间的连接线路的电流的小信号状态变量,ΔiloadDQ为公共参考坐标系DQ中母线所连接负载的电流的小信号状态变量;为微电网小信号状态变量的变化率;A为微电网状态矩阵;Adi为第i个分布式电源的延时状态矩阵;τi为第i个分布式电源的延时。
步骤30)获取微电网闭环小信号模型含有超越项的特征方程
在各分布式电源的延时一致时,式(18)的特征方程为式(19):
CEτ(s,τ)=det(sI-A-Ade-τs) 式(19),
式(19)中,s为时域复平面参数;τ为各分布式电源的一致时延时间,τ1=τ2=...=τn,单位:秒;det(·)表示矩阵行列式;I表示单位矩阵;Ad表示分布式电源的延时状态矩阵,e-τs为超越项。
步骤40)对***特征方法的超越项进行临界特征根轨迹跟踪以计算***稳定裕度
对式(19),当***特征根都在复平面左半平面时,***稳定;当存在特征根在复平面右半平面时,***不稳定;当特征根在复平面左半平面或者虚轴上时,***临界稳定。由于***特征根随着时延时间τ连续变化,因此要确定***稳定裕度τd,即,τ<τd时***稳定,τ>τd时***不稳定,需要确定***可能存在的纯虚特征根和对应的延时裕度。
定义ξ=τω,代入式(19),则,
CEξ(s,ξ)=det(sI-A-Ade-iξ) 式(20),
其中,ξ为时延时间辅助变量,ω为虚特征根幅值;这里i为虚数单位,i2=-1。
ξ在[0,2π]的周期内进行变化,获取式(20)的相应特征根。如果对应于某个ξ存在纯虚特征根,则临界延时时间为:
τc=ξc/abs(ωc) 式(21),
式中,ξc为使***存在纯虚特征根的延时时间辅助变量,abs(ωc)表示对应的纯虚特征根的幅值,τc为临界延时时间。
当ξ在[0,2π]周期内变化时,***可能存在多个临界延时时间,即τc1,τc2...τcL,延时裕度取最小值τd:
τd=min(τc1 τc2 … τcL) 式(22),
在上述实施例中,所述的公共参考坐标系DQ是指第1个分布式电源的dq参考坐标系,其余分布式电源、支路电流、负载电流的状态变量通过坐标变换转换到公共参考坐标系DQ中。在步骤10)中无功功率比例积分控制器和电压比例积分控制器中,由于比例项系数比较小,实际中可以分别简化为无功功率积分控制器和电压积分控制器。在步骤20)中,负载为阻抗型负载。
本实施例通过引入信号通讯延时时间的微电网闭环小信号模型,建立含有超越项的***特征方程,从而实现基于临界特征根跟踪的微电网延时裕度计算方法。针对常规的忽略通讯时延对***动态性能影响的微网二次控制方法,本实施例充分考虑了电力电子接口型微电网惯性小从而导致通讯延时对***稳定性不可忽视的实际情况,计算出***维持稳定的最大延时时间。本实施例的延时裕度计算方法,通过对不同控制器参数与延时裕度间关系的分析,指导控制器设计,从而提升了***稳定性和动态性能。
本发明实施例中的微电网控制***框图如2所示,该控制框图主要包括两层:第一层为各分布式电源的本地控制器,由功率计算、下垂控制和电压电流双环组成;第二层为二次电压控制层,实现无功功率均分和平均电压恢复。二次电压集中控制器采集各分布式电源输出电压、输出无功功率,计算出各二次电压控制量后,将控制指令下发至各分布式电源的本地控制器中。在控制指令下发过程中,通讯时延存在于二次电压集中控制器与各分布式电源本地控制器间,该时延对***动态性能产生影响。
下面例举一个实施例。
仿真***如图3所示,微电网由2个分布式电源,2条连接线路和3个负载组成,负载1连接于母线1,负载2连接于母线2,负载3连接于母线3。***中负载采用阻抗型负载。假设分布式电源1,分布式电源2容量比为1:1,则设计相应的频率下垂系数、电压下垂系数使各分布式电源期望输出有功功率、无功功率比值为1:1。研究在不同控制器参数下的微电网理论延时裕度,并基于MATLAB/Simulink平台搭建微电网仿真模型对理论延时裕度进行仿真验证。
图4为在控制器参数kIQ=0.02,kIV=20下,与***稳定性相关的临界特征根轨迹跟踪示意图。通讯延时辅助变量ξ在[0,2π]变化,2对共轭特征根与***稳定性密切相关,记录下4个经过复平面虚轴的临界特征根A(jωc1),A'(-jωc1),B(jωc2)and B'(-jωc2)及相应的ξ,根据式(21)和式(22)计算出延时裕度τd=0.0588s。
图5是本发明实施例中,在控制器参数0.005≤kIQ≤0.06,5≤kIV≤60下,基于临界特征根跟踪计算的微电网延时裕度与控制器参数的关系。由图可知,随着无功功率控制器积分系数kIQ或电压控制器积分系数kIV的增加,***延时裕度减少,也就是***鲁棒稳定性降低。因此当不同组合控制器参数达到相似的动态性能时,延时裕度将作为附加的鲁棒稳定性指标,指导控制器参数设计,提供***稳定性及动态性能。
图6为微电网采用本发明实施例在某一组控制器参数kIQ=0.02,kIV=20下,3种不同通讯延时对***动态性能的影响中的分散式控制方法的仿真结果。开始运行时,各分布式电源运行于下垂控制模式,0.5s时二次电压控制投入。仿真结果如图6所示,图6(a)为微电网中分布式电源平均电压曲线图,横坐标表示时间,单位:秒,纵坐标表示平均电压,单位:伏。瓦。如图6(a)所示,最初在下垂控制作用下,分布式电源平均电压存在稳态偏差,0.5s后在二次控制作用下,电压幅值提升。由图6(a)可知:***不存在通讯延时时,平均电压较平滑得到达额定值,当延时时间为53ms时,电压曲线经过衰减振荡恢复,当延时时间为61ms时,曲线增幅振荡,***不稳定。图6(b)为分布式电源1无功功率输出曲线图,单位:秒,纵坐标表示无功功率,单位:乏。由图6(b)可知,最初在下垂作用下无功功率均分效果并不理想(少于分布式电源1期望无功功率输出值),0.5s后在二次控制作用下,无功功率输出增加。由图6(b)可知,***不存在通讯延时时,无功功率较平滑得到达期望值,当延时时间为53ms时,功率曲线经过衰减振荡达到控制目标,当延时时间为61ms时,曲线增幅振荡,***不稳定。在二次控制作用下,微电网无功功率均分的效果得到显著改善。图6(c)为分布式电源2无功功率输出曲线图,单位:秒,纵坐标表示无功功率,单位:乏。由图6(c)可知,最初在下垂作用下无功功率均分效果并不理想(高于分布式电源2期望无功功率输出值),0.5s后在二次控制作用下,无功功率输出减少。由图6(c)可知,***不存在通讯延时时,无功功率较平滑得到达期望值,当延时时间为53ms时,功率曲线经过衰减振荡达到控制目标,当延时时间为61ms时,曲线增幅振荡,***不稳定。由图6可知,在此控制器参数下的***延时裕度介于53ms和61ms间,与理论计算值一致。
图7为微电网采用本发明实施例在某一组控制器参数kIQ=0.04,kIV=40下,3种不同通讯延时对***动态性能的影响中的分散式控制方法的仿真结果。开始运行时,各分布式电源运行于下垂控制模式,0.5s时二次电压控制投入。仿真结果如图7所示,图7(a)为微电网中分布式电源平均电压曲线图,横坐标表示时间,单位:秒,纵坐标表示平均电压,单位:伏。瓦。如图7(a)所示,最初在下垂控制作用下,分布式电源平均电压存在稳态偏差,0.5s后在二次控制作用下,电压幅值提升。由图7(a)可知:***不存在通讯延时时,平均电压较平滑得到达额定值,当延时时间为25ms时,电压曲线经过衰减振荡恢复,当延时时间为33ms时,曲线增幅振荡,***不稳定。图7(b)为分布式电源1无功功率输出曲线图,单位:秒,纵坐标表示无功功率,单位:乏。由图7(b)可知,最初在下垂作用下无功功率均分效果并不理想(少于分布式电源1期望无功功率输出值),0.5s后在二次控制作用下,无功功率输出增加。由图6(b)可知,***不存在通讯延时时,无功功率较平滑得到达期望值,当延时时间为25ms时,功率曲线经过衰减振荡达到控制目标,当延时时间为33ms时,曲线增幅振荡,***不稳定。在二次控制作用下,微电网无功功率均分的效果得到显著改善。图7(c)为分布式电源2无功功率输出曲线图,单位:秒,纵坐标表示无功功率,单位:乏。由图7(c)可知,最初在下垂作用下无功功率均分效果并不理想(高于分布式电源2期望无功功率输出值),0.5s后在二次控制作用下,无功功率输出减少。由图7(c)可知,***不存在通讯延时时,无功功率较平滑得到达期望值,当延时时间为25ms时,功率曲线经过衰减振荡达到控制目标,当延时时间为33ms时,曲线增幅振荡,***不稳定。由图6可知,在此控制器参数下的***延时裕度介于25ms和33ms间,与理论计算值一致。
本发明实施例的方法是基于临界特征根跟踪的微电网延时裕度计算方法,基于输出反馈建立含有通讯时延的微电网闭环小信号模型,分析使***稳定的最大延时时间,即延时裕度。针对常规的忽略通讯时延对***动态性能影响的微网二次控制方法,本实施例充分考虑了通讯延时对***稳定性的影响,此外通过研究不同控制器参数与延时裕度间关系,指导控制器设计,从而提升了微电网的鲁棒稳定性和动态性能。
Claims (5)
1.基于临界特征根跟踪的微电网延时裕度计算方法,其特征在于,根据静态反馈输出建立包含通讯延时电压反馈控制量的逆变器闭环小信号模型及分布式电源闭环小信号模型,结合连接网络、负载阻抗的动态方程及分布式电源闭环小信号模型建立微电网小信号模型,从微电网小信号模型获取含有超越项的特征方程,以延时时间辅助变量作为特征方程的变量,求解特征方程在延时时间辅助变量变化周期内的所有纯虚特征根,从所有纯虚特征根对应的临界延时时间中选取最小值作为满足***稳定性要求的延时裕度,所述延时时间辅助变量为分布式电源延时和虚特征根幅值的乘积。
2.根据权利要求1所述基于临界特征根跟踪的微电网延时裕度计算方法,其特征在于,根据静态反馈输出建立的包含通讯延时电压反馈控制量的逆变器闭环小信号模型为:Δxinv、分别为逆变器的闭环小信号状态变量及其变化率,Δxinv1、Δxinv2、Δxinvi、Δxinvn分别为第1个、第2个、第i个、第n个分布式电源的小信号状态变量,分别为第1个、第2个、第i个、第n个分布式电源的无功功率辅助小信号状态变量,第i个分布式电源的无功功率辅助小信号状态变量由表达式:确定,为第i个分布式电源无功功率辅助小信号状态变量的变化率,Qi为第i个分布式电源实际输出的无功功率,nQi为第i个分布式电源的电压下垂特性系数,n为分布式电源的数目,Δγ为分布式电源的电压辅助小信号状态变量,分布式电源的电压辅助小信号状态变量Δγ由表达式:确定,为分布式电源的电压辅助小信号状态变量的变化率,Vi *为第i个分布式电源平均电压的期望值,Vodi为在第i个分布式电源输出电压在其自身参考坐标系dq下的d轴分量,Ainv为分布式电源的状态矩阵,ΔVbDQ为母线电压在公共参考坐标系DQ中的小信号状态变量,ΔVbDQ=[ΔVbDQ1,ΔVbDQ2,…,ΔVbDQl,…,ΔVbDQm]T,ΔVbDQ1、ΔVbDQ2、ΔVbDQl、ΔVbDQm分别为第1根、第2根、第l根、第m根母线的电压在公共参考坐标系DQ中的小信号状态变量,m为母线的数目,Binv为分布式电源对母线电压的输入矩阵,Δu为分布式电源的二次电压小信号控制量,Δu=[Δu1,Δu2,…,Δui,…,Δun]T,Δu1、Δu2、Δui、Δun分别为第1个、第2个、第i个、第n个分布式电源的二次电压小信号控制量,Bu为分布式电源对二次电压小信号控制量的输入矩阵,Δui=KQiΔyinvQi(t-τi)+KViΔyinvV(t-τi),t为当前时刻,τi为第i个分布式电源本地控制器与微网二次电压集中控制器间的通讯时延,KQi、KVi分别为第i个分布式电源的无功功率控制系数、电压控制系数,ΔyinvQi为第i个分布式电源的无功功率输出小信号状态变量,ΔyinvQ、ΔyinvV分别为分布式电源的无功功率输出小信号状态变量、电压输出小信号状态变量,CinvQ、CinvV分别为分布式电源的无功功率输出矩阵、电压输出矩阵。
3.根据权利要求2所述基于临界特征根跟踪的微电网延时裕度计算方法,其特征在于,根据静态反馈输出建立的包含通讯延时电压反馈控制量的分布式电源闭环小信号模型为: 为第i个分布式电源的延时状态矩阵,Bui为第i个分布式电源对二次电压小信号控制量的输入矩阵,CinvQi为第i个分布式电源的无功功率输出矩阵,ΔioDQ为公共参考坐标系DQ中分布式电源输出电流的小信号状态变量,Cinvc为分布式电源的电流输出矩阵。
4.根据权利要求3所述基于临界特征根跟踪的微电网延时裕度计算方法,其特征在于,所述微电网小信号模型为x、分别为微电网小信号状态变量及其变化率,x=[ΔxinvΔilineDQΔiloadDQ]T,ΔilineDQ为公共参考坐标系DQ中分布式电源所连接母线间的连接线路的电流的小信号状态变量,公共参考坐标系DQ中第i个分布式电源所连接母线和第j个分布式电源所连接母线之间的连接线路ij的电流的小信号状态变量为:ΔilineDij、分别为连接线路ij的电流在公共参考坐标系DQ下的D轴小信号分量及其变化率,ΔilineQij、分别为连接线路ij的电流在公共参考坐标系DQ下的Q轴小信号分量及其变化率,rlineij、Llineij分别为连接线路ij的线路电阻和线路电感,ω0为微电网额定角频率,ΔVbusDi、ΔVbusQi分别为第i个分布式电源所连接母线的电压在公共参考坐标系DQ下的D轴分量、Q轴分量,ΔVbusDj、ΔVbusQj分别为第j个分布式电源所连接母线的电压在公共参考坐标系DQ下的D轴分量、Q轴分量,ΔiloadDQ为公共参考坐标系DQ中母线所连接负载的电流的小信号状态变量,公共参考坐标系DQ中第l根母所连接负载的电流的小信号状态变量为:ΔiloadDl、分别为第l根母线所连接负载的电流在公共参考坐标系DQ下的D轴分量及其变化率,ΔiloadQl、分别为第l根母线所连接负载的电流在公共参考坐标系DQ下的Q轴分量及其变化率,Rloadl、Lloadl分别为第l根母线所连接负载的负载电阻、负载电感,ΔVbusDl、ΔVbusQl分别为第l根母线的电压在公共参考坐标系DQ下的D轴分量、Q轴分量,Adi为第i个分布式电源的延时状态矩阵。
5.根据权利要求4所述基于临界特征根跟踪的微电网延时裕度计算方法,其特征在于,从微电网小信号模型获取含有超越项的特征方程的方法为:在分布式电源的延时一致时得到微电网小信号模型的特征方程:
CEτ(s,τ)=det(sI-A-Ade-τs),s为时域复平面参数,τ为各分布式电源的一致时延时间,CEτ(·)表示各分布式电源一致时延τ时得到的微电网小信号模型的特征方程,det(·)为矩阵行列式,I为单位矩阵,Ad为分布式电源的延时状态矩阵,e-τs为超越项。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710456420.4A CN107294085B (zh) | 2017-06-16 | 2017-06-16 | 基于临界特征根跟踪的微电网延时裕度计算方法 |
PCT/CN2018/084937 WO2018228068A1 (zh) | 2017-06-16 | 2018-04-27 | 基于临界特征根跟踪的微电网延时裕度计算方法 |
US16/618,378 US20200293703A1 (en) | 2017-06-16 | 2018-04-27 | Microgrid delay margin calculation method based on critical characteristic root tracking |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710456420.4A CN107294085B (zh) | 2017-06-16 | 2017-06-16 | 基于临界特征根跟踪的微电网延时裕度计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107294085A CN107294085A (zh) | 2017-10-24 |
CN107294085B true CN107294085B (zh) | 2019-12-17 |
Family
ID=60096713
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710456420.4A Active CN107294085B (zh) | 2017-06-16 | 2017-06-16 | 基于临界特征根跟踪的微电网延时裕度计算方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US20200293703A1 (zh) |
CN (1) | CN107294085B (zh) |
WO (1) | WO2018228068A1 (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107294085B (zh) * | 2017-06-16 | 2019-12-17 | 东南大学 | 基于临界特征根跟踪的微电网延时裕度计算方法 |
CN107994564A (zh) * | 2017-10-27 | 2018-05-04 | 东南大学 | 基于特征根聚类的多重时滞微电网电压稳定性分析方法 |
CN108363306B (zh) * | 2018-03-20 | 2020-04-24 | 东南大学 | 基于线性二次型优化的微电网分布式控制器参数确定方法 |
WO2020141009A1 (en) * | 2019-01-04 | 2020-07-09 | Vestas Wind Systems A/S | A hybrid renewable power plant |
CN109787234B (zh) * | 2019-01-25 | 2021-04-23 | 国网上海市电力公司 | 含vsc接口的分布式电源超高次谐波稳定模式获取方法 |
CN109946963B (zh) * | 2019-04-23 | 2021-10-15 | 北京航天飞腾装备技术有限责任公司 | 一种判断多回路控制***裕度的方法 |
CN110443302B (zh) * | 2019-08-02 | 2023-06-09 | 天津相和电气科技有限公司 | 基于特征融合与深度学习的负荷辨识方法及其应用 |
CN110649642B (zh) * | 2019-09-29 | 2021-09-17 | 山东理工大学 | 交直流配电***电压协调控制方法及交直流配电*** |
CN114069718B (zh) * | 2020-08-03 | 2024-03-22 | 北京机械设备研究所 | 一种并联变换器的同步控制装置和方法 |
CN112260251B (zh) * | 2020-10-12 | 2022-07-05 | 国网河北省电力有限公司经济技术研究院 | 一种微电网控制周期稳定性分析方法、*** |
CN112670992B (zh) * | 2021-01-22 | 2023-11-07 | 上海交通大学 | 含能量路由器的配电网稳定性分析和失稳校正方法及*** |
CN112865094B (zh) * | 2021-03-11 | 2022-12-06 | 南方电网科学研究院有限责任公司 | 多端直流输电***低压线路重启的协调控制方法及装置 |
CN113078645B (zh) * | 2021-05-20 | 2022-09-27 | 合肥工业大学 | 一种考虑延时与拓扑切换的微电网参数自适应控制方法 |
CN117353396B (zh) * | 2023-12-06 | 2024-03-08 | 国网浙江省电力有限公司信息通信分公司 | 一种基于启停曲线的火电机组调度优化方法和装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102623992A (zh) * | 2012-04-12 | 2012-08-01 | 山东大学 | 基于旋转坐标虚拟阻抗的孤岛微电网控制及优化方法 |
CN103472731A (zh) * | 2013-09-24 | 2013-12-25 | 南方电网科学研究院有限责任公司 | 一种微电网小信号稳定性分析并参数协调整定的方法 |
CN104578097A (zh) * | 2014-12-28 | 2015-04-29 | 国网山东省电力公司日照供电公司 | 一种链式svg控制器的电压增益调节装置及控制方法 |
CN105162134A (zh) * | 2015-08-26 | 2015-12-16 | 电子科技大学 | 新型微电网***及其功率均衡控制策略和小信号建模方法 |
CN106532715A (zh) * | 2016-12-30 | 2017-03-22 | 东南大学 | 一种基于非线性状态观测器的微电网分散式电压控制方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7508224B2 (en) * | 2006-04-07 | 2009-03-24 | Michael Lamar Williams | Method for measuring stability margin at a node of a polyphase power grid |
CN101408908A (zh) * | 2008-11-26 | 2009-04-15 | 天津大学 | 基于优化的电力***实用时滞稳定裕度计算方法 |
CN107294085B (zh) * | 2017-06-16 | 2019-12-17 | 东南大学 | 基于临界特征根跟踪的微电网延时裕度计算方法 |
-
2017
- 2017-06-16 CN CN201710456420.4A patent/CN107294085B/zh active Active
-
2018
- 2018-04-27 WO PCT/CN2018/084937 patent/WO2018228068A1/zh active Application Filing
- 2018-04-27 US US16/618,378 patent/US20200293703A1/en not_active Abandoned
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102623992A (zh) * | 2012-04-12 | 2012-08-01 | 山东大学 | 基于旋转坐标虚拟阻抗的孤岛微电网控制及优化方法 |
CN103472731A (zh) * | 2013-09-24 | 2013-12-25 | 南方电网科学研究院有限责任公司 | 一种微电网小信号稳定性分析并参数协调整定的方法 |
CN104578097A (zh) * | 2014-12-28 | 2015-04-29 | 国网山东省电力公司日照供电公司 | 一种链式svg控制器的电压增益调节装置及控制方法 |
CN105162134A (zh) * | 2015-08-26 | 2015-12-16 | 电子科技大学 | 新型微电网***及其功率均衡控制策略和小信号建模方法 |
CN106532715A (zh) * | 2016-12-30 | 2017-03-22 | 东南大学 | 一种基于非线性状态观测器的微电网分散式电压控制方法 |
Non-Patent Citations (2)
Title |
---|
具有同步发电机特性的并联逆变器小信号稳定性分析;颜湘武 等;《电网技术》;20160331;第40卷(第3期);第910-917页 * |
逆变型分布式电源微网小信号稳定性动态建模分析;张建华;《电力***自动化》;20101125;第34卷(第22期);第97-102页 * |
Also Published As
Publication number | Publication date |
---|---|
US20200293703A1 (en) | 2020-09-17 |
WO2018228068A1 (zh) | 2018-12-20 |
CN107294085A (zh) | 2017-10-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107294085B (zh) | 基于临界特征根跟踪的微电网延时裕度计算方法 | |
CN108363306B (zh) | 基于线性二次型优化的微电网分布式控制器参数确定方法 | |
Hu et al. | A coordinated control of hybrid ac/dc microgrids with PV-wind-battery under variable generation and load conditions | |
WO2018121732A1 (zh) | 一种基于非线性状态观测器的微电网分散式电压控制方法 | |
Zhao et al. | Distributed frequency control for stability and economic dispatch in power networks | |
Ravinder et al. | Investigations on shunt active power filter in a PV-wind-FC based hybrid renewable energy system to improve power quality using hardware-in-the-loop testing platform | |
CN109256803A (zh) | 虚拟同步机孤岛运行小信号建模及参数灵敏度计算方法 | |
Behera et al. | A novel resilient control of grid-integrated solar PV-hybrid energy storage microgrid for power smoothing and pulse power load accommodation | |
CN113131521A (zh) | 虚拟同步机多机并联稳定控制及其惯量匹配方法 | |
Kou et al. | Stable and optimal load sharing of multiple PMSGs in an islanded DC microgrid | |
Dong et al. | Output control method of microgrid VSI control network based on dynamic matrix control algorithm | |
Singh et al. | Control of single-phase distributed PV-battery microgrid for smooth mode transition with improved power quality | |
Fahad et al. | A coordination control strategy for power quality enhancement of an active distribution network | |
Beus et al. | A model predictive control approach to operation optimization of an ultracapacitor bank for frequency control | |
Wang et al. | Power quality enhancement in islanded microgrids via closed-loop adaptive virtual impedance control | |
Singh et al. | Control strategy for multiple residential solar PV system in distribution network with improved power quality | |
CN108879797A (zh) | 一种主动配电网端口pq控制方法 | |
CN109066765A (zh) | 基于一致性策略的孤岛微电网自适应同步频率控制方法 | |
Wong et al. | An optimal secondary multi-bus voltage and reactive power sharing control based on non-iterative decoupled linearized power flow for islanded microgrids | |
Arteaga et al. | Control of energy storage and photovoltaic systems using model predictive control | |
Pradhan et al. | Adaptive position observer for multimode wind-BES based microgrid interfaced to distribution network | |
Huang et al. | Coupling characteristic analysis and synchronization stability control for Multi-Paralleled VSCs system under symmetric faults | |
Yang et al. | A coordinated optimization method considering time-delay effect of islanded photovoltaic microgrid based on modified backtracking search algorithm | |
CN107994564A (zh) | 基于特征根聚类的多重时滞微电网电压稳定性分析方法 | |
Hossain et al. | Distributed control scheme to regulate power flow and minimize interactions in multiple microgrids |
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 |