一种基于一致性理论的孤岛微电网分层分布式控制策略的制作方法

本发明涉及微电网控制技术领域，特别涉及一种孤岛微电网的分层分布式控制方法。

背景技术：

随着用户对供电可靠性和电能质量的关注不断提高，以及太阳能、风能等各种形式的可再生能源大量利用，分布式电源(distributedgenerator，dg)因其充分利用各地丰富的清洁和可再生能源而得到大力发展。微电网作为由dg、储能装置、能量转换装置、可控负荷和监控、保护装置等单元汇集而成的小型发配电系统，能够实现自我控制、保护和管理，被公认为解决分布式电源接入配电网最有效的方案。特别地，在孤岛运行模式下，微电网可以给海岛和边防等偏远地区供电，也可以在外部电网故障时继续向重要负荷供电，具有很大的推广价值。

但由于孤岛微电网失去了大电网的支撑，网络结构的改变和负荷的扰动都易造成系统电压和频率的不稳定，并破坏负载分配机制的准确性，因此孤岛微电网的控制问题成为了亟待解决的难题。同时，由于现代电力系统中分布式电源特性各异，数量越来越多且接入具有分散特点。而现有的孤岛微电网分层控制方法多采用传统的集中式，二级控制和三级控制都需要依靠微电网中央控制器。该方法面临计算负担大、通信复杂以及灵活性不足等限制。因此，本领域技术人员需要提供一种新的孤岛微电网分层分布式控制方法，能够克服集中式控制的缺陷。

技术实现要素：

本发明的目的在于克服现有技术的上述不足，提供一种分层分布式协调控制策略，解决含分布式电源的孤岛微电网中电压和频率的控制问题，该方法不仅提高了孤岛微电网系统安全稳定和经济运行水平，而且不需要中央控制器，能够有效地解决集中式控制方式可靠性差及鲁棒性低等缺陷。

为实现上述目的，本发明采用了以下技术方案：一种基于一致性理论的孤岛微电网分层分布式控制策略，该方法的技术方案如下：

根据孤岛微电网系统的频率和电压调整以及经济优化的不同时间尺度，将控制目标划分为三个层次来实现，即：初级下垂控制、二级频率和电压控制、三级经济调度。

本发明的技术方案主要包括以下步骤：

(1)首先引入初级下垂控制策略，并设计基于下垂控制的分布式电源逆变器多环控制器，包括逆变器的硬件部分、外环功率控制器以及内环电压电流控制器，以实现分布式电源的快速响应，稳定系统电压和频率。

(2)然后基于有限时间一致性算法、平移下垂曲线法以及电压-无功下垂系数修正法，发明了分布式二级频率和电压控制策略，包含分布式二级频率控制和分布式二级电压-无功控制，以去中心化的方式修正初级下垂控制带来的系统频率和电压偏移，并实现无功功率的比例分配。

(3)最后以微电网系统总发电成本最优为目标建立经济调度模型，并采用拉格朗日对偶分解法对该模型进行求解，进一步基于有限时间一致性算法，发明孤岛微电网三级分布式经济调度方法，使各台分布式电源在本地求解最优输出有功功率，实现孤岛微电网的经济运行。

三层控制之间存在着通信联系，共同嵌入在每台dg中，其中初级下垂控制由各台dg自主执行，而分布式二级频率和电压控制以及分布式三级经济调度除了需要本地信息外，还需要通过通信网络与相邻的dg进行信息交换。

进一步的，步骤(1)所述的“初级下垂控制策略”的具体内容为：

下垂控制是利用与传统同步发电机相似的电压-无功和频率-有功的下垂特性将系统不平衡功率动态分配给并联运行的各台dg，从而获取稳定的电压和频率。

进一步的，步骤(2)所述的“分布式二级频率和电压控制策略”的具体内容为：

a.建立分布式二级频率控制

针对初级下垂控制带来的系统频率偏移，二级频率控制旨在实现频率的无差调节。当系统负荷出现正的波动，初级下垂控制动作，使dg输出有功功率增加，系统频率降低，为了使各台dg的输出频率恢复至额定值，采用下垂曲线平移法补偿频率。

传统的频率调节方法是通过实测输出角频率ωi与额定角频率ω^*比较得到频率补偿量，这种方法可能导致频率响应较慢，并且在干扰情况下会产生较大的频率偏差。而dg的实时输出有功功率pi可以很容易由功率计算单元计算得到，因此本发明将采集各台dg的实时输出有功功率pi和三级经济调度下发的最优有功功率参考值prefi进行计算得到频率补偿量。

此外，若各台dg仅利用本地信息各自计算出频率补偿量δωi进行分散式调节，由于功率计算单元可能产生计算误差，会造成并联的各台dg频率补偿量不相等，容易产生振荡和超调，导致系统频率失稳。因此，本发明将各台dg频率补偿量的平均值作为各台dg最终的频率补偿量。

为了使各台dg在本地估算出平均频率补偿量实现分布式二级频率控制，本发明采用有限时间一致性算法，通过通信网络获取相邻dg的频率补偿量进行迭代计算以获得全局所有dg频率补偿量的平均值。

各台dg将最终估算出的频率补偿量反馈至p-ω初级下垂控制环节，即可实现分布式的二级频率控制，使系统频率恢复至额定值。并且由于频率是全局变量，因此采用上述分布式频率控制方法不会破坏各台dg之间有功功率的比例分配。

b.建立分布式二级电压—无功控制

传统的分散式电压调节通过各台dg各自改变励磁调压，将输出电压均恢复至额定值，这种方法将恶化各台dg之间本就不成比例的无功功率分配。因此，为了同时能够实现合理的无功分配，可允许小范围的电压偏差，将dg的输出电压幅值调节至额定值附近即可。

各台dg与相邻dg交换输出相电压幅值信息，采用有限时间一致性算法，各自在本地进行迭代计算，得到所有dg输出相电压幅值平均值。

各台dg将最终估算得到的输出相电压幅值平均值与输出相电压额定值un相减，得到电压偏差量，再将电压偏差量经过pi控制器gi(s)得到电压幅值补偿量δei，最后将δei反馈至q-e初级下垂控制环节，即可实现分布式的二级电压控制，使各台dg的输出电压调节至额定电压值附近的允许范围之内，并且给无功功率的比例分配提供了可能性。

进而采用q-e下垂系数修正法实现无功功率的比例分配，修正过程一直持续到下垂系数偏差量等于零时，q-e下垂系数便稳定在恒定值，此时各台dg的输出无功功率均调节至与初始下垂系数相匹配，实现了dg之间输出无功功率的比例分配。

进一步的，步骤(3)所述的“三级分布式经济调度方法”的具体内容为：

在初级控制和二级控制保证孤岛微电网安全稳定运行的基础上，三级经济调度旨在实现孤岛微电网的经济运行。

假设微电网系统中含有n台dg，建立第i台dg的发电成本ci与其有功功率pi的二次函数：ci(pi)＝aipi²+bipi+ci，式中：ai、bi和ci为相关的燃料消耗系数；以其所有dg发电成本最小为经济调度的目标函数同时建立满足供需平衡的等式约束和dg的发电容量上下限的不等式约束pi,min≤pi≤pi,max，其中：pl是系统负荷需求功率，pi,max和pi,min分别为第i台dg的发电功率上下限。

引入拉格朗日对偶分解法对上述经济调度模型进行求解，可以得到交替求解方程：式中：i＝1,2,…,n，k为迭代次数；ρ>0为迭代步长；λ为拉格朗日乘子；prefi^k+1代表第k+1次迭代下dgi计算出的有功功率参考值。系统总不平衡功率和拉格朗日乘子λ均为全局变量，即需要中央控制器收集和处理系统中的负荷信息以及每台dg的有功功率参考值信息。为了消除中央控制器，得到完全分布式的经济调度方法，本发明将利用有限时间一致性算法实现每个dg本地独立迭代求解。

与现有技术相比，本发明具有以下优点：

本发明根据孤岛微电网系统控制优化的不同时间尺度，将控制目标划分初级下垂控制、二级频率和电压控制、三级经济调度三个层次，并区别于传统的集中式控制方式，将采取基于有限时间一致性算法的分布式控制方式以实现上述分层控制结构。

附图说明

图1是本发明所提的分层分布式控制策略整体控制框图；

图2是二级频率控制原理图；

图3是三级分布式经济调度方法的流程图。

具体实施方式

下面结合附图对本发明做进一步说明：

如图1所示，一种基于一致性理论的孤岛微电网分层分布式控制策略，该控制策略包括以下步骤：

(1)通过下垂控制策略向dg逆变器内环控制器提供频率和电压幅值参考值，使dg快速响应系统中的负荷波动，稳定系统频率和电压，实现微电网实时控制，建立初级下垂控制；

(2)调整由初级下垂控制带来的系统频率和电压偏差，同时在dg之间保持精确的有功和无功功率比例分配，建立分布式二级频率控制和分布式二级电压—无功控制；

(3)确定各台dg的最优发电有功功率参考值，并下发给初级下垂控制层，实现孤岛微电网系统经济运行，建立分布式三级经济调度。

具体的讲，步骤(1)中初级下垂控制的具体方法为：

下垂控制是利用与传统同步发电机相似的电压-无功和频率-有功的下垂特性将系统不平衡功率动态分配给并联运行的各台dg，从而获取稳定的电压和频率。初级下垂控制表达式为：式中：ωi和ei分别为dgi输出电压角频率和幅值参考值，用于产生逆变器内环控制器的控制信号；ωi^*和ei^*分别为系统的额定角频率和额定电压幅值；mi和ni分别为p-ω和q-e下垂系数；pi和qi分别为dgi实时输出的有功和无功功率；prefi为三级经济调度下达的dgi最优有功功率参考值。

另外，由于频率为全局变量，即并联运行的各台dg输出角频率ωi相等，且各台dg的额定角频率ωi^*均为100π(rad/s)。为了能够实现二级频率控制，在一级下垂控制中应使各台dg的空载频率ωi0^*相等(其中ωi0^*＝ωi^*+miprefi)，需选取p-ω下垂系数mi，使其满足：m1pref1＝m2pref2＝…＝mnprefn＝c，式中：c为常数，即各台dg的p-ω下垂系数需根据三级经济调度下达的最优有功功率参考值确定。为了实现dg间无功功率的比例分配，各台dg的q-e下垂系数将由二级电压-无功控制确定。

具体的讲，步骤(2)中有限时间一致性算法的具体方法为：

一致性算法的本质是通过本地节点与相邻节点的信息交互，更新本地节点的状态量，使得网络中所有节点的状态量收敛于稳定的共同值。离散时间一致性算法可描述如下：式中：xi代表第i个节点vi的状态量，如系统的电压、频率和功率等；k为迭代次数；wij为节点vi和vj之间的权重因子。上式也可以写成矩阵形式：x(k+1)＝wx(k)，式中：w为权重矩阵。

根据一致性算法逐步迭代到当k→∞时，所有节点的状态量都会趋于一致，且最终的收敛结果为节点状态量初始值的平均值，即：因此一致性算法只需各节点和相邻节点进行通信，便可进行系统状态量平均值的估算，无需中央控制器进行全局状态量的收集和计算，大大减少了计算负担和通信成本，提高了计算速度。

将改写为其中：in为n阶单位矩阵，jn＝(1/n)inin^t。由公式x(k+1)＝wx(k)和可以得出因此，要想实现在有限的步数之内达到精确收敛，只需找到一组与网络拓扑一致的权重矩阵{wm}m＝1,2,…,k，满足：其中k为通信网络的拉普拉斯矩阵的不同非零特征值个数。

因此由有限时间一致性算法可知，通过设置一组权重矩阵wm可以使得系统中所有节点的状态量在k步迭代之内达到平均一致性收敛，即：

具体的讲，步骤(2)中分布式二级频率控制的具体方法为：

通过两台并联dg的初级下垂控制特性曲线来说明二级频率调节过程，如附图2所示。假设初始时刻dg1和dg2均分别运行于功率参考点即a点和b点，此时两台dg的输出频率均稳定在额定频率ω^*。当系统负荷出现正的波动，初级下垂控制动作，使dg输出有功功率增加，dg1和dg2运行点分别偏移至c点和d点，系统频率降低至ω。为了使各台dg的输出频率恢复至额定值，采用下垂曲线平移法，使dg1和dg2的运行点分别向上平移至e点和f点，即系统的运行频率被修正至额定值。dg1和dg2的p-ω下垂特性曲线斜率分别为m1和m2，则下垂曲线的平移量即第i台dg的频率补偿量δωi可根据下式计算：

传统的频率调节方法是通过实测输出角频率ωi与额定角频率ω^*比较得到频率补偿量，这种方法可能导致频率响应较慢，并且在干扰情况下会产生较大的频率偏差。而dg的实时输出有功功率pi可以很容易由功率计算单元计算得到，因此本文将采集各台dg的实时输出有功功率pi和三级经济调度下发的最优有功功率参考值prefi进行计算得到频率补偿量。

此外，若各台dg仅利用本地信息各自计算出频率补偿量δωi进行分散式调节，由于功率计算单元可能产生计算误差，会造成并联的各台dg频率补偿量不相等，容易产生振荡和超调，导致系统频率失稳。因此，本文将各台dg频率补偿量的平均值作为各台dg最终的频率补偿量，即：

根据有限时间一致性算法的原理，获得权重矩阵，则上式可以在k步迭代之内达到平均一致性收敛，即：

各台dg将最终估算出的频率补偿量反馈至p-ω初级下垂控制环节，即可实现分布式的二级频率控制，使系统频率恢复至额定值。并且由于频率是全局变量，因此采用上述分布式频率控制方法不会破坏各台dg之间有功功率的比例分配。

具体的讲，步骤(2)中分布式二级电压—无功控制的具体方法为：

不同于二级频率控制，由于线路阻抗不等及负荷分布不均，将导致各台dg的输出电压水平不相等，dg间较小的电压差就会造成较大的无功偏差，进而dg输出无功功率并不能按q-e下垂系数比例分配，容易造成dg过载。传统的分散式电压调节通过各台dg各自改变励磁调压，将输出电压均恢复至额定值，这种方法将恶化各台dg之间本就不成比例的无功功率分配。因此，为了同时能够实现合理的无功分配，可允许小范围的电压偏差，将dg的输出电压幅值调节至额定值附近即可。分布式电压控制过程如下：

各台dg与相邻dg交换输出相电压幅值信息uoj，采用有限时间一致性算法，各自在本地进行迭代计算，即：与分布式频率控制采用相同的通信网络，则上式同样可在k步迭代之内达到平均一致性收敛，即：

各台dg将最终估算得到的输出相电压幅值平均值与输出相电压额定值un相减，得到电压偏差量，再将电压偏差量经过pi控制器gi(s)得到电压幅值补偿量δei，最后将δei反馈至q-e初级下垂控制环节，即可实现分布式的二级电压控制，使各台dg的输出电压调节至额定电压值附近的允许范围之内，并且给无功功率的比例分配提供了可能性。

进而采用q-e下垂系数修正法实现无功功率的比例分配，分布式无功控制过程如下：

各台dg与相邻dg交换q-e下垂系数与无功功率的乘积信息，获取下垂系数偏差量δni，即：式中：b为电压幅值与q-e下垂系数之间的耦合增益，lij为图的拉普拉斯矩阵中的元素。将求得的偏差量δni经过pi控制器hi(s)得到下垂系数修正项δni，并将δni反馈至q-e初级下垂控制环节，更新q-e下垂系数：ni(t)＝ni0+δni(t)。

上述修正过程一直持续到下垂系数偏差量δni等于零时，q-e下垂系数便稳定在恒定值，此时各台dg的输出无功功率均调节至与初始下垂系数相匹配，即满足n10q1＝n20q2＝…＝nn0qn，实现了dg之间输出无功功率的比例分配。

具体的讲，步骤(3)中分布式三级经济调度的具体方法为：

在初级控制和二级控制保证孤岛微电网安全稳定运行的基础上，三级经济调度旨在实现孤岛微电网的经济运行。假设微电网系统中含有n台dg，dgi的发电成本一般是其输出有功功率的二次函数：ci(pi)＝aipi²+bipi+ci，式中：ai、bi和ci为相关的燃料消耗系数。那么微电网系统经济调度目标函数为：

同时，微电网的经济运行需要满足供需平衡约束以及dg的发电容量不等式约束如下(忽略网损)：式中pl表示系统总有功负荷，pi,max和pi,min分别表示dgi的出力上下限。

为了得到分布式的经济调度方法，首先引入拉格朗日对偶分解法对上述经济调度模型进行求解，可以得到交替求解方程如下：式中：i＝1,2,…,n，k为迭代次数；ρ>0为迭代步长；λ为拉格朗日乘子；prefi^k+1代表第k+1次迭代下dgi计算出的有功功率参考值。令：代表第k+1次迭代下系统的总不平衡功率，进一步可以得到如下迭代求解方程：考虑到dg输出功率上下限，实际的输出有功功率参考值应满足：重复以上迭代步骤，当满足r^k+1＝λ^k+1-λ^k＝ρpe^k+1＜ε时，dg有功功率参考值将获得最优解，其中ε是收敛精度。

然而上述过程中使用到的系统总不平衡功率pe和拉格朗日乘子λ均为全局变量，即需要中央控制器收集和处理系统中的负荷信息以及每台dg的有功功率参考值信息。为了消除中央控制，得到完全分布式的经济调度方法，本发明利用有限时间一致性算法使每台dg在本地独立估算出pe。估算方法如下：

a.定义第k次迭代下dgi的不平衡功率为在初级下垂控制的作用下，各台dg的实时输出有功功率能够满足：则可以推导出即第k次迭代下所有dg的不平衡功率总和就等于系统总不平衡功率。

b.各台dg与相邻dg交换不平衡功率按照有限时间一致性算法各自在本地进行迭代计算，则dgi的不平衡功率可以在k步迭代之内达到平均一致性收敛，即：

c.则每一台dg都可以估算出第k次迭代下的系统总不平衡功率：

按照如上估算方法，每一台dg只需要本地信息和相邻dg的信息便可估算出系统总不平衡功率，pe不再是一个全局变量，从而拉格朗日乘子的更新可以由各dg自主完成，则公式可以本地进行独立迭代计算，分布式求解出每一台dg的最优有功功率参考值。

综上所述，附图3给出了所发明的分布式经济调度方法具体流程。

以上所述的实施例仅仅是对本发明的优选实施方式进行描述，并非对本发明的范围进行限定，在不脱离本发明设计精神的前提下，本领域普通技术人员对本发明的技术方案作出的各种变形和改进，均应落入本发明权利要求书确定的保护范围内。