CN104659810B - 一种用于不确定通信拓扑的微电网协同控制方法 - Google Patents
一种用于不确定通信拓扑的微电网协同控制方法 Download PDFInfo
- Publication number
- CN104659810B CN104659810B CN201510044890.0A CN201510044890A CN104659810B CN 104659810 B CN104659810 B CN 104659810B CN 201510044890 A CN201510044890 A CN 201510044890A CN 104659810 B CN104659810 B CN 104659810B
- Authority
- CN
- China
- Prior art keywords
- agency
- represent
- micro
- control
- capacitance sensor
- 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.)
- Expired - Fee Related
Links
- 238000004891 communication Methods 0.000 title claims abstract description 69
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000004146 energy storage Methods 0.000 claims abstract description 37
- 230000008859 change Effects 0.000 claims abstract description 22
- 230000008878 coupling Effects 0.000 claims abstract description 17
- 238000010168 coupling process Methods 0.000 claims abstract description 17
- 238000005859 coupling reaction Methods 0.000 claims abstract description 17
- 239000011159 matrix material Substances 0.000 claims description 63
- 230000005611 electricity Effects 0.000 claims description 5
- 230000008901 benefit Effects 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 230000001360 synchronised effect Effects 0.000 description 3
- 101100499229 Mus musculus Dhrsx gene Proteins 0.000 description 2
- 238000011960 computer-aided design Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- 238000005303 weighing 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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Electrical Variables (AREA)
Abstract
本发明公开了一种用于不确定通信拓扑的微电网协同控制方法，包括：步骤10）微电网储能系统运行在恒压/恒频控制模式下，维持微电网的功率平衡，发生事故时，储能系统代理自动进行一次控制；步骤20）在不确定通信拓扑下，测算牵制控制的预定义一致值；步骤30）多代理系统中的除牵制代理外的各代理通过通信耦合与牵制代理寻求同步；步骤40）各分布式电源代理根据预定义一致值增发功率，完成系统二次控制，实现总功率缺额在分布式电源之间共享。该控制方法以牵制控制为基础，消除了中央控制器的需求和复杂的通信拓扑，能够适应微电网通信拓扑的变化，满足微电网中通信线路开断和分布式电源即插即用操作的需求。
Description
技术领域
本发明属于微电网运行控制领域，具体来说，涉及一种用于不确定通信拓扑的微电网协同控制方法。
背景技术
未来的电力系统中，将会有越来越多的可能包含风力发电机、光伏、微型燃气轮机、小水电、储能系统、电动汽车等可控单元的微电网。因此，自治微电网的稳定控制，近来得到了特别的关注。由于运行模式的转变和分布式发电单元的间歇性，微电网中经常发生频率和电压的扰动或故障。
基于多智能体系统的分布式控制策略被认识到能够在维持自治微电网稳定方面发挥重要作用。分布式控制的优点包括抗不确定干扰和分布式信息更新能力，从而使信息有效共享，最终使得决策制定和实施更加迅速。牵制控制是一种多智能体系统有效的分布式控制方式，引入牵制控制后，需控制的代理个数将明显减少，可以应用于大型复杂系统。同时，牵制控制可以在变通信拓扑结构下进行分布式信息交互，能够适应微电网通信拓扑的变化，满足微电网中通信线路开断和分布式电源即插即用操作的需求。
发明内容
技术问题：本发明所要解决的技术问题是：提供一种用于不确定通信拓扑的微电网协同控制方法，该控制方法以牵制控制为基础，是一种分布式控制方法，消除了中央控制器的需求和复杂的通信拓扑，适用于不确定通信拓扑的微电网，能够适应微电网通信拓扑的变化，满足微电网中通信线路开断和分布式电源即插即用操作的需求。
技术方案：为解决上述技术问题，本发明采用的技术方案如下：
一种用于不确定通信拓扑的微电网协同控制方法，微电网采用多代理系统，包括分布式电源代理、储能系统代理、负荷代理，共有n个代理，其中一部分代理是被控制的牵制代理，其余代理通过与牵制代理的通信耦合，以分布的方式跟踪同步；该控制方法包括以下步骤：
步骤10)进行一次控制，获取功率缺额：在孤岛模式下，微电网储能系统运行在恒压/恒频控制模式下，维持微电网的功率平衡，发生事故时，储能系统代理自动进行一次控制，在一次控制中，储能系统输出有功功率的参考值为微电网总有功功率的缺额，储能系统输出无功功率的参考值为微电网总无功功率的缺额；
步骤20)测算牵制控制的预定义一致值：在不确定通信拓扑下，测算牵制控制的预定义一致值，该预定义一致值包括有功功率预定义一致值ΔPc和无功功率预定义一致值ΔQc；
步骤30)代理通过通信耦合与牵制代理寻求同步：多代理系统中的除牵制代理外的各代理通过通信耦合与牵制代理寻求同步，使得参与牵制控制的分布式电源代理的有功功率增发量和无功功率增发量达到预定义一致值；
步骤40)进行二次控制，实现总功率缺额的共享：各分布式电源代理的有功功率增发量和无功功率增发量达到预定义一致值时，各分布式电源代理根据预定义一致值增发功率，完成系统二次控制，实现总功率缺额在分布式电源之间共享。
进一步，所述的步骤20)中，依照式(1)测算牵制控制的预定义一致值：
其中，mΔ(t)表示参与牵制控制的代理总数，xΔ(t)表示牵制代理的数量。
进一步，所述的步骤20)中，基于受到分布式电源的发电容量限制，依照式(2)测算牵制控制的预定义一致值：
其中，YΔ(t)表示达到最大容量限制的分布式电源代理集合，ΔPMh表示第h个分布式电源的最大有功功率增加量，mΔ(t)表示参与牵制控制的代理总数，xΔ(t)表示牵制代理的数量，yΔ(t)表示达到最大容量限制的分布式电源代理数量，ΔQMh表示第h个分布式电源的最大无功功率增加量。
进一步，所述的步骤30)包括以下过程：
定义一个时变矩阵ΘΦ和一个切换序列S表示不确定通信拓扑，切换序列 k为整数，且k≥0；表示切换序列S中第1个切换区间，表示切换序列S中第2个切换区间，表示切换序列S中第k+1个切换区间；表示在时间t0时的通信拓扑，表示在时间t1时的通信拓扑，表示在时间tk时的通信拓扑；Φ表示通信拓扑集合，时变矩阵ΘΦ用于表征微电网的在不同时刻的拓扑结构；ΘΦ基于切换序列S，随着切换区间的改变而改变：
第i个代理在不确定通信拓扑下的牵制控制，依照式(3)确定：
其中，Ni(t)表示与第i个代理相邻的代理的集合，Ni(t)为一个时变的集合；Δ(t)是介于时间tk和tk+1之间的时间区间，Δ(t)＝[tk,tk+1)；表示切换序列S中切换区间上的拓扑矩阵；表示第j个代理增发有功功率的参考值，表示第i个代理增发有功功率的参考值，di,Δ(t)表示第i个代理的牵制控制增益，di,Δ(t)≥0；di,Δ(t)＝0表示第i个代理没有被牵制控制，表示第i个代理增发无功功率的参考值，表示第j个代理增发无功功率的参考值，表示的微分，表示的微分；WΔ(t)表示时变的通信耦合矩阵，wij,Δ(t)是WΔ(t)矩阵的第(i，j)个元素，若代理i和代理j通过通信线路连接，wij,Δ(t)≠0；否则，wij,Δ(t)＝0；
如式(4)所示，定义控制误差ei：
式(3)根据ei，改写为式(5)：
其中，ei表示第i个代理的控制误差，ej表示第j个代理的控制误差，ep,i表示第i个代理的有功功率控制误差，eq,i是第i个代理的无功功率控制误差；δei表示ei的微分；表示切换序列S中切换区间上的拓扑矩阵；
整个微电网的牵制控制，使用如式(6)所示的矩阵形式表示：
其中，E表示微电网控制误差矩阵，δE是E的微分，WΔ(t)表示时变的通信耦合矩阵，EP表示微电网有功功率控制误差矩阵，EP＝[ep,1、ep,2、…、ep,n]，ep,1表示第1个代理的有功功率控制误差，ep,2表示第2个代理的有功功率控制误差，ep,n表示第n个代理的有功功率控制误差；δEP表示EP的微分；DΔ(t)表示微电网牵制矩阵，DΔ(t)＝[d1,Δ(t)、d2,Δ(t)、…、dn,Δ(t)]，d1,Δ(t)表示第1个代理的牵制控制增益，d2,Δ(t)表示第2个代理的牵制控制增益，dn,Δ(t)表示第n个代理的牵制控制增益；Im表示单位矩阵；EQ表示微电网无功功率控制误差矩阵，EQ＝[eq,1、eq,2、…、eq,n]，eq,1表示第1个代理的无功功率控制误差，eq,2表示第2个代理的无功功率控制误差，eq,n表示第n个代理的无功功率控制误差，δEQ表示EQ的微分；
当微电网的拓扑结构发生变化，式(6)改写成式(7)：
ΘΦδE＝WP,Δ(t)ΘΦE 式(7)
式中，WP,Δ(t)表示微电网牵制控制的更新矩阵。
进一步，所述的步骤30)中，wij,Δ(t)的数值按照式(8)计算，适应通信拓扑的改变：
其中，θ为常数，0<θ<1；ni表示与第i个代理相邻的代理的数量，nj表示与第j个代理相邻的代理的数量。
进一步，所述的步骤30)中，WP,Δ(t)如式(9)所示：
其中，N1(t)表示与第1个代理相邻的代理集合，Ni(t)表示与第i个代理相邻的代理集合，Nn(t)表示与第n个代理相邻的代理集合，N1(t)、Ni(t)、Nn(t)均为时变的集合；w1j,Δ(t)表示WΔ(t)矩阵的第(1，j)个元素，w1i,Δ(t)表示WΔ(t)矩阵的第(1，i)个元素，w1n,Δ(t)表示WΔ(t)矩阵的第(1，n)个元素，wi1,Δ(t)表示WΔ(t)矩阵的第(i，1)个元素，wij,Δ(t)表示WΔ(t)矩阵的第(i，j)个元素，win,Δ(t)表示WΔ(t)矩阵的第(i，n)个元素，wn1,Δ(t)表示WΔ(t)矩阵的第(n，1)个元素，wni,Δ(t)表示WΔ(t)矩阵的第(n，i)个元素，wnj,Δ(t)表示WΔ(t)矩阵的第(n，j)个元素，d1,Δ(t)表示第1个代理的牵制控制增益，di,Δ(t)表示第i个代理的牵制控制增益，dn,Δ(t)表示第n个代理的牵制控制增益。
有益效果：与现有技术相比，本发明具有以下有益效果：本发明的控制方法是一种适用于不确定通信拓扑微电网协同控制的方法。本发明以牵制控制为基础，是一种分布式控制方法，消除了中央控制器的需求和复杂的通信拓扑，能够适应微电网通信拓扑的变化，满足微电网中通信线路开断和分布式电源即插即用操作的需求。应用牵制控制，计算牵制控制的预定义一致值，实施过程中考虑考虑分布式电源的容量限制和发电容量约束；多代理系统中的分布式电源代理通过通信耦合与牵制代理寻求同步，当微电网拓扑结构发生变化时，该方法可以对通信耦合矩阵和牵制控制矩阵进行更新，能够适应微电网通信拓扑的变化；当牵制一致达到预定义值时，完成总功率缺额在分布式电源之间共享，从而实现微电网的分布式协同控制。
附图说明
图1是本发明的流程框图。
图2是本发明实施例中微电网仿真系统的结构示意图。
图3是本发明实施例中微电网多代理系统的通信拓扑图。
图4是本发明实施例中光伏退出运行，储能系统和各分布式电源输出的有功功率曲线图。
图5是本发明实施例中光伏退出运行，储能系统和各分布式电源输出的无功功率曲线图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施案例对本发明进行深入地详细说明。应当理解，此处所描述的具体实施案例仅仅用以解释本发明，并不用于限定发明。
如图1所示，本发明中，微电网采用多代理系统，包括分布式电源代理、储能系统代理、负荷代理，共有n个代理，其中一部分代理是被控制的牵制代理，其余代理通过与牵制代理的通信耦合，以分布的方式跟踪同步。本发明的用于不确定通信拓扑的微电网协同控制方法，包括以下步骤：
步骤10)进行一次控制：在孤岛模式下，微电网储能系统运行在恒压/恒频控制(V/f控制)模式下，维持微电网的功率平衡，发生事故时，储能系统代理自动进行一次控制，在一次控制中，储能系统输出有功功率的参考值为微电网总有功功率的缺额，储能系统输出无功功率的参考值为微电网总无功功率的缺额。
步骤20)在不确定通信拓扑下，测算牵制控制的预定义一致值，该预定义一致值包括有功功率预定义一致值ΔPc和无功功率预定义一致值ΔQc。
作为一种优选方案，步骤20)依照式(1)测算牵制控制的预定义一致值：
其中，mΔ(t)表示参与牵制控制的代理总数，xΔ(t)表示牵制代理的数量。mΔ(t)和xΔ(t)根据通信拓扑的改变而改变，例如，在通信拓扑插入操作中，mΔ(t)变为(mΔ(t)+1)，当通信线路断开时xΔ(t)也会改变。
作为另一种优选方案，考虑受到分布式电源的发电容量限制，所述的步骤20)依照式(2)测算牵制控制的预定义一致值：
其中，YΔ(t)表示达到最大容量限制的分布式电源代理集合，ΔPMh表示第h个分布式电源的最大有功功率增加量，mΔ(t)表示参与牵制控制的代理总数，xΔ(t)表示牵制代理的数量，yΔ(t)表示达到最大容量限制的分布式电源代理数量，ΔQMh表示第h个分布式电源的最大无功功率增加量。
步骤30)多代理系统中的除牵制代理外的各代理通过通信耦合与牵制代理寻求同步，使得参与牵制控制的分布式电源代理的有功功率增发量和无功功率增发量达到预定义一致值。步骤30)具体包括以下过程：
定义一个时变矩阵ΘΦ和一个切换序列S表示不确定通信拓扑，切换序列 k为整数，且k≥0；表示切换序列S中第1个切换区间，表示切换序列S中第2个切换区间，表示切换序列S中第k+1个切换区间；表示在时间t0时的通信拓扑，表示在时间t1时的通信拓扑，表示在时间tk时的通信拓扑；Φ表示通信拓扑集合，时变矩阵ΘΦ用于表征微电网的在不同时刻的拓扑结构；ΘΦ基于切换序列S，随着切换区间的改变而改变。
第i个代理在不确定通信拓扑下的牵制控制，依照式(3)确定：
其中，Ni(t)表示与第i个代理相邻的代理的集合，Ni(t)为一个时变的集合；Δ(t)是介于时间tk和tk+1之间的时间区间，Δ(t)＝[tk,tk+1)；表示切换序列S中切换区间上的拓扑矩阵；表示第j个代理增发有功功率的参考值，表示第i个代理增发有功功率的参考值，di,Δ(t)表示第i个代理的牵制控制增益，di,Δ(t)≥0；di,Δ(t)＝0表示第i个代理没有被牵制控制，表示第i个代理增发无功功率的参考值，表示第j个代理增发无功功率的参考值，表示的微分，表示的微分；WΔ(t)表示时变的通信耦合矩阵，wij,Δ(t)是WΔ(t)矩阵的第(i，j)个元素，若代理i和代理j通过通信线路连接，wij,Δ(t)≠0；否则，wij,Δ(t)＝0。wij,Δ(t)的数值按照式(8)计算，适应通信拓扑的改变：
其中，θ为常数，0<θ<1；ni表示与第i个代理相邻的代理的数量，nj表示与第j个代理相邻的代理的数量。
如式(4)所示，定义控制误差ei：
式(3)根据ei，改写为式(5)：
其中，ei表示第i个代理的控制误差，ej表示第j个代理的控制误差，ep,i表示第i个代理的有功功率控制误差，eq,i是第i个代理的无功功率控制误差；δei表示ei的微分；表示切换序列S中切换区间上的拓扑矩阵。
整个微电网的牵制控制，使用如式(6)所示的矩阵形式表示：
其中，E表示微电网控制误差矩阵，δE是E的微分，WΔ(t)表示时变的通信耦合矩阵，EP表示微电网有功功率控制误差矩阵，EP＝[ep,1、ep,2、…、ep,n]，ep,1表示第1个代理的有功功率控制误差，ep,2表示第2个代理的有功功率控制误差，ep,n表示第n个代理的有功功率控制误差；δEP表示EP的微分；DΔ(t)表示微电网牵制矩阵，DΔ(t)＝[d1,Δ(t)、d2,Δ(t)、…、dn,Δ(t)]，d1,Δ(t)表示第1个代理的牵制控制增益，d2,Δ(t)表示第2个代理的牵制控制增益，dn,Δ(t)表示第n个代理的牵制控制增益；Im表示单位矩阵；EQ表示微电网无功功率控制误差矩阵，EQ＝[eq,1、eq,2、…、eq,n]，eq,1表示第1个代理的无功功率控制误差，eq,2表示第2个代理的无功功率控制误差，eq,n表示第n个代理的无功功率控制误差，δEQ表示EQ的微分。
当微电网的拓扑结构发生变化，式(6)改写成式(7)：
ΘΦδE＝WP,Δ(t)ΘΦE 式(7)
式中，WP,Δ(t)表示微电网牵制控制的更新矩阵。WP,Δ(t)如式(9)所示：
其中，N1(t)表示与第1个代理相邻的代理集合，Ni(t)表示与第i个代理相邻的代理集合，Nn(t)表示与第n个代理相邻的代理集合，N1(t)、Ni(t)、Nn(t)均为时变的集合；w1j,Δ(t)表示WΔ(t)矩阵的第(1，j)个元素，w1i,Δ(t)表示WΔ(t)矩阵的第(1，i)个元素，w1n,Δ(t)表示WΔ(t)矩阵的第(1，n)个元素，wi1,Δ(t)表示WΔ(t)矩阵的第(i，1)个元素，wij,Δ(t)表示WΔ(t)矩阵的第(i，j)个元素，win,Δ(t)表示WΔ(t)矩阵的第(i，n)个元素，wn1,Δ(t)表示WΔ(t)矩阵的第(n，1)个元素，wni,Δ(t)表示WΔ(t)矩阵的第(n，i)个元素，wnj,Δ(t)表示WΔ(t)矩阵的第(n，j)个元素，d1,Δ(t)表示第1个代理的牵制控制增益，di,Δ(t)表示第i个代理的牵制控制增益，dn,Δ(t)表示第n个代理的牵制控制增益。
当微电网的拓扑结构发生变化时，利用式(8)对wij,Δ(t)进行更新，进而利用式(9)对WP,Δ(t)进行更新。
步骤40)各分布式电源代理的有功功率增发量和无功功率增发量达到预定义一致值时，各分布式电源代理根据预定义一致值增发功率，完成系统二次控制，实现总功率缺额在分布式电源之间共享。其他分布式电源代理不具备功率增发能力，只是得到功率增发参考信息，并不参与系统二次控制，不进行功率增发。
本发明的控制方法可以适用于不确定通信拓扑的微电网，以牵制控制为基础，在实施过程中考虑考虑分布式电源的容量限制和发电容量约束，是一种分布式控制方法，消除了中央控制器的需求和复杂的通信拓扑，能够适应微电网通信拓扑的变化，满足微电网中通信线路开断和分布式电源即插即用操作的需求。
下面例举一个实施例。
图2为本实施例采用的仿真微电网结构图，图3为本实施例采用的仿真微电网通信拓扑图。该仿真微电网由微型燃气轮机(简称MT)、风力发电机(简称WT)、光伏系统(简称PV)、储能系统(ESS)及5个负荷单元(Load1、Load2、Load3、Load4、Load5)组成，各分布式单元由电力电子元件接入0.38kV低压配电网。如图2和图3所示，系统有8个分布式代理(Agent)，分别用A1，A2，A3，A4，A5，A6，A7，A8表示，其中A3为牵制代理。一个代理只能与其在通信拓扑上直接相邻的代理进行通信。系统有8个本地控制器(LC)，分别用LC1，LC2，LC3，LC4，LC5，LC6，LC7，LC8表示。基于电力系统计算机辅助设计/含直流电磁暂态仿真(英文简称：PSCAD/EMTDC)平台搭建仿真微电网模型，在矩阵实验室(英文简称MATLAB)中模拟多代理系统，建立多智能体系统基于牵制的微电网分布式协同控制算法程序，利用PSACD中的用户自定义接口(英文简称UDI)模型将MATLAB中的算法与电力系统计算机辅助设计(英文简称PSCAD)模型联合运行，从而利用联合仿真技术实现本发明的控制方法的仿真验证。
针对孤岛运行模式下的微电网进行分布式电源拔出操作的情况进行了仿真，验证本发明方法的控制效果，即适用于不确定通信拓扑的微电网协同控制方法，能够适应微电网通信拓扑的变化，满足微电网中通信线路开断和分布式电源即插即用操作的需求。微电网运行在孤岛模式，储能系统运行在V/f控制(对应中文为：恒压/恒频控制)模式下，其余分布式电源运行于PQ(对应中文为：恒功率控制)控制模式下。
6秒时，微电网进行分布式电源拔出操作，光伏由于故障退出运行，微电网功率平衡被打破，微电网通信拓扑发生改变。首先，储储能系统进行一次控制，增发功率维持孤岛微电网的功率平衡。由于光伏代理A1切除，参与牵制控制的代理总数nΔ(t)变为(nΔ(t)-1)，预定义一致值做出相应的改变，通信耦合矩阵和牵制控制更新矩阵进行相应的更新。风力发电机组和微型燃气轮机进行协同增发有功功率和无功功率，寻求与牵制代理同步。风力发电机组增发后达到容量限制0.06MW，这是预定义一致值按照式(2)进行调整，风力发电机增发至最大限制功率，微型燃气轮机增发与预定义一致值跟踪同步。当分布式电源输出达到预定义值时，储能系统功率输出自动变为0，即完成总功率缺额在分布式电源之间共享，由分布式电源取代储能系统弥补系统功率缺额，实现微电网分布式协同控制。仿真结果如图4至图5所示。
图4中，PESS表示储能系统的有功功率输出曲线。储能系统运行于V/f控制模式，6秒时，光伏退出运行，微电网出现功率缺额，储能系统增发有功功率维持系统功率平衡。在风力发电机和微型燃气轮机协同增发后，储能系统的有功功率输出逐渐变为0。PMT表示微型燃气轮机有功功率输出曲线，PWT表示风力发电机有功功率输出曲线。微型燃气轮机和风力发电机运行于PQ控制模式下，6秒时微型燃气轮机和风力发电机协同增发有功功率，代替储能系统平衡系统功率缺额，同时使储能系统有功功率输出逐渐变为0。其中，由于风力发电机发电容量限制，其增发后最大输出有功功率只能达到0.06MW。PPV表示光伏有功功率输出曲线。光伏运行于PQ控制模式，维持功率稳定输出，6秒时，光伏由于故障退出运行，输出有功功率变为0。
图5中，QESS表示储能系统的无功功率输出曲线。储能系统运行于V/f控制模式，6秒时，光伏退出运行，储能系统增发无功功率维持系统功率平衡。在风力发电机和微型燃气轮机协同增发后，储能系统的无功功率输出逐渐变为0。QMT表示微型燃气轮机无功功率输出曲线，QWT表示风力发电机无功功率输出曲线。微型燃气轮机和风力发电机运行于PQ控制模式下，6秒时微型燃气轮机和风力发电机协同增发无功功率，代替储能系统平衡系统功率缺额，同时使储能系统无功功率输出逐渐变为0。QPV表示光伏无功功率输出曲线。光伏运行于PQ控制模式，维持功率稳定输出，6秒时，光伏由于故障退出运行，输出无功功率变为0。
从本实施例可以看出，采用本发明的控制方法后，微电网协同控制可以有效的适应通信拓扑的变化，满足微电网中分布式电源即插即用操作的需求，各分布式电源协同控制维持系统功率平衡，同时在实施工程中考虑了分布式电源发电容量的限制。本发明提出的方法有很好地控制效果。
Claims (6)
1.一种用于不确定通信拓扑的微电网协同控制方法，其特征在于，微电网采用多代理系统，包括分布式电源代理、储能系统代理、负荷代理，共有n个代理，其中一部分代理是被控制的牵制代理，其余代理通过与牵制代理的通信耦合，以分布的方式跟踪同步；该控制方法包括以下步骤：
步骤10)进行一次控制，获取功率缺额：在孤岛模式下，微电网储能系统运行在恒压/恒频控制模式下，维持微电网的功率平衡，发生事故时，储能系统代理自动进行一次控制，在一次控制中，储能系统输出有功功率的参考值为微电网总有功功率的缺额，储能系统输出无功功率的参考值为微电网总无功功率的缺额；
步骤20)测算牵制控制的预定义一致值：在不确定通信拓扑下，测算牵制控制的预定义一致值，该预定义一致值包括有功功率预定义一致值ΔPc和无功功率预定义一致值ΔQc；
步骤30)代理通过通信耦合与牵制代理寻求同步：多代理系统中的除牵制代理外的各代理通过通信耦合与牵制代理寻求同步，使得参与牵制控制的分布式电源代理的有功功率增发量和无功功率增发量达到预定义一致值；
步骤40)进行二次控制，实现总功率缺额的共享：各分布式电源代理的有功功率增发量和无功功率增发量达到预定义一致值时，各分布式电源代理根据预定义一致值增发功率，完成系统二次控制，实现总功率缺额在分布式电源之间共享。
2.按照权利要求1所述的用于不确定通信拓扑的微电网协同控制方法，其特征在于，所述的步骤20)中，依照式(1)测算牵制控制的预定义一致值：
其中，mΔ(t)表示参与牵制控制的代理总数，xΔ(t)表示牵制代理的数量。
3.按照权利要求1所述的用于不确定通信拓扑的微电网协同控制方法，其特征在于，所述的步骤20)中，基于受到分布式电源的发电容量限制，依照式(2)测算牵制控制的预定义一致值：
其中，YΔ(t)表示达到最大容量限制的分布式电源代理集合，ΔPMh表示第h个分布式电源的最大有功功率增加量，mΔ(t)表示参与牵制控制的代理总数，xΔ(t)表示牵制代理的数量，yΔ(t)表示达到最大容量限制的分布式电源代理数量，ΔQMh表示第h个分布式电源的最大无功功率增加量。
4.按照权利要求1所述的用于不确定通信拓扑的微电网协同控制方法，其特征在于，所述的步骤30)包括以下过程：
定义一个时变矩阵ΘΦ和一个切换序列S表示不确定通信拓扑，切换序列 k为整数，且k≥0；表示切换序列S中第1个切换区间，表示切换序列S中第2个切换区间，表示切换序列S中第k+1个切换区间；表示在时间t0时的通信拓扑，表示在时间t1时的通信拓扑，表示在时间tk时的通信拓扑；Φ表示通信拓扑集合，时变矩阵ΘΦ用于表征微电网的在不同时刻的拓扑结构；ΘΦ基于切换序列S，随着切换区间的改变而改变：
第i个代理在不确定通信拓扑下的牵制控制，依照式(3)确定：
其中，Ni(t)表示与第i个代理相邻的代理的集合，Ni(t)为一个时变的集合；Δ(t)是介于时间tk和tk+1之间的时间区间，Δ(t)＝[tk,tk+1)；表示切换序列S中切换区间上的拓扑矩阵；表示第j个代理增发有功功率的参考值，表示第i个代理增发有功功率的参考值，di,Δ(t)表示第i个代理的牵制控制增益，di,Δ(t)≥0；di,Δ(t)＝0表示第i个代理没有被牵制控制，表示第i个代理增发无功功率的参考值，表示第j个代理增发无功功率的参考值，表示的微分，表示的微分；WΔ(t)表示时变的通信耦合矩阵，wij,Δ(t)是WΔ(t)矩阵的第(i，j)个元素，若代理i和代理j通过通信线路连接，wij,Δ(t)≠0；否则，wij,Δ(t)＝0；
如式(4)所示，定义控制误差ei：
式(3)根据ei，改写为式(5)：
其中，ei表示第i个代理的控制误差，ej表示第j个代理的控制误差，ep,i表示第i个代理的有功功率控制误差，eq,i是第i个代理的无功功率控制误差；δei表示ei的微分；
整个微电网的牵制控制，使用如式(6)所示的矩阵形式表示：
其中，E表示微电网控制误差矩阵，δE是E的微分，EP表示微电网有功功率控制误差矩阵，EP＝[ep,1、ep,2、…、ep,n]，ep,1表示第1个代理的有功功率控制误差，ep,2表示第2个代理的有功功率控制误差，ep,n表示第n个代理的有功功率控制误差；δEP表示EP的微分；DΔ(t)表示微电网牵制矩阵，DΔ(t)＝[d1,Δ(t)、d2,Δ(t)、…、dn,Δ(t)]，d1,Δ(t)表示第1个代理的牵制控制增益，d2,Δ(t)表示第2个代理的牵制控制增益，dn,Δ(t)表示第n个代理的牵制控制增益；Im表示单位矩阵；EQ表示微电网无功功率控制误差矩阵，EQ＝[eq,1、eq,2、…、eq,n]，eq,1表示第1个代理的无功功率控制误差，eq,2表示第2个代理的无功功率控制误差，eq,n表示第n个代理的无功功率控制误差，δEQ表示EQ的微分；
当微电网的拓扑结构发生变化，式(6)改写成式(7)：
ΘΦδE＝WP,Δ(t)ΘΦE 式(7)
式中，WP,Δ(t)表示微电网牵制控制的更新矩阵。
5.按照权利要求4所述的用于不确定通信拓扑的微电网协同控制方法，其特征在于，所述的步骤30)中，wij,Δ(t)的数值按照式(8)计算，适应通信拓扑的改变：
其中，θ为常数，0<θ<1；ni表示与第i个代理相邻的代理的数量，nj表示与第j个代理相邻的代理的数量。
6.按照权利要求4所述的用于不确定通信拓扑的微电网协同控制方法，其特征在于，所述的步骤30)中，WP,Δ(t)如式(9)所示：
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510044890.0A CN104659810B (zh) | 2015-01-28 | 2015-01-28 | 一种用于不确定通信拓扑的微电网协同控制方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510044890.0A CN104659810B (zh) | 2015-01-28 | 2015-01-28 | 一种用于不确定通信拓扑的微电网协同控制方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104659810A CN104659810A (zh) | 2015-05-27 |
| CN104659810B true CN104659810B (zh) | 2017-04-05 |
Family
ID=53250634
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510044890.0A Expired - Fee Related CN104659810B (zh) | 2015-01-28 | 2015-01-28 | 一种用于不确定通信拓扑的微电网协同控制方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104659810B (zh) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104901313B (zh) * | 2015-06-16 | 2017-07-25 | 北京亿利智慧能源科技有限公司 | 基于Gossip通信算法的微网电压不平衡补偿方法 |
| CN105138006B (zh) * | 2015-07-09 | 2017-08-04 | 哈尔滨工程大学 | 一种时滞非线性多智能体系统的协同追踪控制方法 |
| CN106129999B (zh) * | 2016-07-01 | 2018-06-15 | 东南大学 | 基于有限时间一致性的直流微电网分布式协同控制方法 |
| CN106169756B (zh) * | 2016-07-15 | 2018-06-05 | 河海大学 | 一种计及经济特性的孤岛微网频率协同控制方法 |
| CN106410808B (zh) * | 2016-09-27 | 2018-04-24 | 东南大学 | 包含恒功率和下垂控制的通用型微电网群分布式控制方法 |
| CN106972512B (zh) * | 2017-05-16 | 2021-03-05 | 中国石油大学(华东) | 一种电力系统同步稳定运行的分散控制器设计方法 |
| CN110611333B (zh) * | 2019-11-11 | 2022-12-13 | 燕山大学 | 一种考虑通信问题的孤岛微电网分布式协调控制方法 |
| CN111082989B (zh) * | 2020-03-05 | 2021-11-12 | 东南大学 | 基于网络镜像和全局传播率的微电网分布式通信拓扑设计方法 |
| CN114629168B (zh) * | 2022-03-30 | 2024-03-01 | 合肥工业大学 | 一种考虑通信延时及拓扑变换的微电网次级控制方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE102008037575A1 (de) * | 2008-11-21 | 2010-07-29 | EnBW Energie Baden-Württemberg AG | Computergestütztes Verfahren zur Optimierung der Energienutzung in einem lokalen System |
| CN104218681B (zh) * | 2014-09-28 | 2016-01-13 | 东南大学 | 一种用于降低孤岛微电网切负荷成本的控制方法 |
-
2015
- 2015-01-28 CN CN201510044890.0A patent/CN104659810B/zh not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN104659810A (zh) | 2015-05-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104659810B (zh) | 一种用于不确定通信拓扑的微电网协同控制方法 | |
| CN105634020B (zh) | 基于有限时间一致性的孤岛微电网分布式协调控制方法 | |
| CN108363306B (zh) | 基于线性二次型优化的微电网分布式控制器参数确定方法 | |
| CN104659811B (zh) | 一种基于牵制的微电网分布式协同控制方法 | |
| WO2018058804A1 (zh) | 通用型包含恒功率和下垂控制的微电网群分布式控制方法 | |
| CN106549384B (zh) | 一种含upfc电力系统的通用潮流计算方法 | |
| CN106532715B (zh) | 一种基于非线性状态观测器的微电网分散式电压控制方法 | |
| CN107664954A (zh) | 一种电力系统电磁暂态可视化仿真模型自动生成方法 | |
| CN104600749B (zh) | 一种适用于微电网的分层同步控制方法 | |
| CN103532140B (zh) | 一种含分布式电源的配电网故障后供电恢复的方法及系统 | |
| CN107294085A (zh) | 基于临界特征根跟踪的微电网延时裕度计算方法 | |
| CN105119283B (zh) | 面向对等模式下微电网的基于牵制的分布式协同控制方法 | |
| CN106849112A (zh) | 基于非支配邻域免疫算法的配电网多目标无功优化方法 | |
| CN105186500B (zh) | 一种基于重加权加速Lagrangian的配电网能量分散协调优化方法 | |
| CN106129999A (zh) | 基于有限时间一致性的直流微电网分布式协同控制方法 | |
| CN106374498A (zh) | 一种考虑二次电压频率控制的微电网潮流计算方法 | |
| Zhang et al. | A distributed multi-control-center dynamic power flow algorithm based on asynchronous iteration scheme | |
| CN106208099A (zh) | 一种基于二层规划的电力系统无功优化方法及其应用 | |
| Wang et al. | A test bed for self-regulating distribution systems: Modeling integrated renewable energy and demand response in the GridLAB-D/MATLAB environment | |
| CN108879667A (zh) | 一种电网闭环控制潮流仿真方法 | |
| CN104218681B (zh) | 一种用于降低孤岛微电网切负荷成本的控制方法 | |
| CN114676569B (zh) | 电网仿真分析算例及其生成方法、生成系统、设备、介质 | |
| CN103996325B (zh) | 一种培训系统中多级调度一体化模式下的教案制作方法 | |
| CN110120673B (zh) | 基于戴维南等值参数辨识的分布式输配协同无功优化方法及系统 | |
| CN108599974A (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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170405 |