一种低压微电网小扰动稳定性分析方法与流程

本申请所描述的实施例一般涉及微电网的静态稳定领域。尤其是，涉及一种低压微电网小扰动稳定性分析方法。

背景技术

近年来随着分布式电源接入配电网的数量和容量的增加，其带来的扰动对配电网的稳定运行的影响日益严重。设想一个优良的、稳固的控制装置，肯定可以提供更为良好的相对高效和可靠的产品服务。由于逆变型的分布式电源的本身惯性较小，虽然通过电力电子变换器接入电网，但受到扰动时，仍易发生振荡失稳，分析微电网系统小扰动稳定性，对电力工业有重要的意义。

小扰动法是电力系统稳定判据的依据，将非线性的动态行为转变成线性动态行为，已有较为成熟的理论。但现代的电力系统并不是简单的只有一个发电机元件需要列写状态方程，而是通过电力电子变换器接入电网，响应时间较长，网络及负荷将影响系统的稳定性。分布式电源的开关频率较高忽略分布式电源侧直流母线的动态过程，看做理想电源，以逆变器的模型替代各分布式电源。

传统下垂控制p-f/q-u下微电网的小扰动稳定性已有价位成熟的研究。但对于低压微电网，即为用户侧的小型微网，线路阻抗主要是电阻，传统的下垂已不再适用。因此基于阻性下垂控制的低压微电网小扰动稳定性的研究是极为必要的。

技术实现要素：

相应地，本发明的实施例提供一种低压微电网小扰动稳定性分析方法，解决了以传统下垂控制应用在低压微电网的电压降落较大的问题。

本发明的一种低压微电网小扰动稳定性分析方法，将低压微电网分为网络负荷模块和逆变器模块，分别构建网络负荷模块小扰动模型和逆变器模块小扰动模型，根据网络负荷模块小扰动模型和逆变器模块小扰动模型得到微电网小扰动模型，采用小扰动特征分析法，把描述电力系统动态行为的非线性微分方程组在平衡点做线性化处理。求解模型的状态空间的系数矩阵，通过计算一些提供的初始参数和系统参数状态下的特征值分析低压微电网系统的稳定性；

所述逆变器模块采用阻性下垂曲线p-u/q-f作为控制方式，p-u下垂曲线为有功功率-电压下垂曲线，q-f下垂曲线为无功功率-频率下垂曲线。

进一步的，上述p-u下垂曲线方程为u＝un-mp(p-pn)，un为逆变器的额定的输出电压，u为逆变器输出的电压，p为逆变器输出的有功功率，pn为逆变器输出的额定有功功率，mp为有功功率下垂系数；

q-f下垂曲线方程为f＝fn+nq(q-qn)，f为逆变器的输出频率，fn为逆变器的额定的输出频率，q为逆变器输出的无功功率，qn为逆变器输出的额定无功功率。nq为无功功率下垂系数。

进一步的，所述逆变器模块的小扰动模型的构建方法包括构建控制单元小扰动模型、电压电流环小扰动模型和接口电路小扰动模型，将控制单元小扰动模型、电压电流环小扰动模型和接口电路小扰动模型联立得到逆变器模块的小扰动模型。

进一步的，所述下垂曲线方程线性化处理后得到下垂控制线性化的表达式为：

式中，分别表示电压和频率偏差的求导。分别表示逆变器有功功率和无功功率偏差的求导值。

进一步的，建立控制单元的小扰动模型，表达式为：

式中，表示控制器在旋转坐标系下的输出到电压电流环的电压参考值偏差的导数，δudqref表示表示控制器输出到电压电流环的电压参考值的偏差δp，δq分别表示逆变器有功功率和无功功率的偏差，δil1dq表示滤波电路的电感电流偏差，δucdq、δil2dq表示输入功率计算模块的电压偏差和电流的偏差,a1，a2，a3，a4表示系数矩阵，其中

进一步的，所述电压电流环小扰动模型的表达式为：

上式中，δuidqref表示经过电压电流环输出到逆变单元的电压偏差，δidq，δudq分别表示正交派克变换后所得的d/q轴电流偏差和电压偏差，δf表示系统频率偏差，a5，a6，a7，a8，a9表示系数矩阵。

进一步的，接口电路的小扰动模型为：

上式中，和分别表示滤波电容的电压偏差的导数、电流偏差的导数和耦合电感的电流偏差的导数，δuidq为逆变器输出的电压偏差，δufdq为馈线上的电压偏差，a10，a11，a12，a13表示系数矩阵，令δuidqref＝δuidq。

进一步的，所述逆变器模块的小扰动模型的表达式如下：

式中，δx为状态变量，为状态变量δx的导数，

其中：其中k表示逆变器的数量，δpk、δqk分别代表第k个逆变器有功功率和无功功率的偏差，δfk表示第k个逆变器的频率偏差，δukdqref表示第k个控制器在旋转坐标系下的输出到电压电流环的电压参考值偏差,δikdq和δukdq分别表示第k个控制器正交派克变换后所得的d/q轴电压和电流偏差值，δuckdq和δil2kdq分别表示第k个逆变器输入功率计算模块的电压偏差和电流的偏差，δil1kdq表示第k个逆变器的滤波电路的电感电流偏差。

δufkdq表示第k个逆变器馈线上的电压偏差，下标1,2…k表示第k个逆变器与网络馈线的连接节点，a、b表示系数矩阵。

进一步的，所述网络负荷模块的小扰动模型包括：

负荷小扰动模型：

式中，状态变量δilkdq为第k条负荷支路流过负荷的电流偏差，δufdq为馈线电压的偏差，an，bn1，bn2表示系数矩阵；

网络小扰动模型：

δufdq＝cn1×δildq+cn2×δilinedq；

式中，δildq为负荷支路流过负荷的电流偏差，δilinedq为分布式电源之间的电流偏差，cn1，cn2表示系数矩阵。

进一步的，将负荷小扰动模型、网络小扰动模型和逆变器模块小扰动模型综合，得到低压微电网小扰动模型，表达式为：

式中，am为系统状态矩阵；

计算状态矩阵am的特征值，根据特征值的分布对微电网稳定性进行分析，具体分析方法包括：

如果特征值的实部均为负，表明微网在小扰动下呈现衰减振荡，能保持小扰动下的稳定；

若有特征值的实部为正，则系统不稳定，对于不稳定的特征根，对其进行灵敏度分析；通过改变相应的参数，再次计算特征根分析系统的稳定性。

本发明所述逆变模块包括p-u/q-f下垂控制器、电压电流环控制器、逆变单元、lc滤波器和耦合电抗。

本发明针对线路阻感较大的低压微电网，逆变模块的控制系统采用阻性下垂控制p-u/q-f的控制方法，即有功功率调节逆变器输出电压的幅值，无功功率调节逆变器输出电压的频率，解决了以传统下垂控制应用在低压微电网的电压降落较大的问题，通过有功-电压/无功-频率下垂曲线来调节电压和频率。通过电力电子变换器接入电网，响应时间较长，网络及负荷将影响系统的稳定性。通过建立p-u/q-f控制下的微电网的小扰动数学模型，通过正交派克变换和线性化处理得系统的状态矩阵，求解特征值，观察特征值的分布来分析系统的稳定性，对于不稳定的特征根，对其进行灵敏度分析，通过改变相应的参数，再次计算特征根分析系统的稳定性，为低压微电网的控制提供依据。

附图说明

图1示出了根据本发明的实施例的逆变器模块控制系统结构示意图；

图2示出了根据本发明的实施例的逆变型分布式电源功率传输等效模型。

具体实施方式

分布式电源的开关频率较高忽略分布式电源侧直流母线的动态过程，看作理想电源，以逆变器的模型替代各分布式电源，因此该发明中逆变器的小扰动模型即为分布式电源的小扰动模型，本实施例的基于阻性下垂控制的低压微电网小扰动稳定性分析方法，其基本原理是将低压微电网分为网络负荷模块和逆变器模块，分别构建网络负荷模块小扰动模型和逆变器模块小扰动模型，根据网络负荷模块小扰动模型和逆变器模块小扰动模型得到微电网小扰动模型，通过微电网小扰动模型得到状态特征值，根据特征值分布确定系统的稳定性；

表一不同电压等级线路阻抗参数

以单个逆变器与微电网交流母线相连为例，u为逆变电源输出的电压幅值,e为微电网公共交流母线电压幅值,z为线路阻抗的幅值,θ为线路阻抗的相角,δ是以微电网公共交流母线电压矢量作为参考时逆变电源输出电压矢量的相角。

由图2可得下列公式：

(10)代入(11)可得逆变器输出的有功功率和无功功率：

传统的下垂控制，以高压线路为主，忽略了r，则z＝x,θ＝90°，并联运行时，δ很小，sinδ＝δ，cosδ＝1。

则有对(14)、(15)进行微分运算得又ω＝2πf，可得传统的下垂控制数学表达式为：f＝fn-kp(p-pn)和u＝un-kq(q-qn)；

同理，低压线路下，由表一可以发现阻抗比为7.7，因此新型的下垂控制忽略x,z＝r，θ＝0°,则有

对(16)、(17)进行微分运算得

又ω＝2πf，可得传统的下垂控制数学表达式为

u＝un-mp(p-pn)和f＝fn+nq(q-qn)；

本实施例具体包括以下步骤：

步骤1、确定微电网的线路参数以及逆变器的元件参数，包括：线路参数，负荷参数，各个逆变器的控制器和接口电路的参数，滤波电容和耦合电感的值l1、c、l2、有功功率下垂系数mp、无功功率下垂系数nq、电流和电压环的比例和积分系数，初始稳态运行的参数，电网的额定电压，额定频率，逆变器本地坐标与参考坐标的夹角、各个馈线的电压、逆变器的输出电压和输出电流、耦合电感电流。

步骤2、对负荷建立小扰动模型，即：

式中，状态变量δilkdq为第k条负荷支路流过负荷的电流偏差，δufdq为馈线电压的偏差，an，bn1，bn2表示系数矩阵。

步骤3、对网络建立小扰动模型，即：

δufdq＝cn1×δildq+cn2×δilinedq(2)

式中，δildq为负荷支路流过负荷的电流偏差，δilinedq为分布式电源之间的电流偏差，cn1，cn2表示系数矩阵。

步骤4、逆变器模块的控制系统采用基于阻性的下垂控制，功角或频率通过无功功率来进行调节，逆变器的电压通过调节有功功率来进行调节；

p-u下垂曲线即有功功率—电压下垂曲线，方程如下：

u＝un-mp(p-pn)；

其中，u为逆变器输出的电压，un为逆变器的额定的输出电压，p为逆变器输出的有功功率，pn为逆变器输出的额定有功功率，mp为有功功率下垂系数。

q-f下垂曲线即无功功率—频率下垂曲线，方程如下：

f＝fn+nq(q-qn)；

其中，f为逆变器的输出频率，fn为逆变器的额定的输出频率，q为逆变器输出的无功功率，qn为逆变器输出的额定无功功率,nq为无功功率下垂系数，当系统阻抗呈现为阻性时，无功功率和有功功率可以实现解耦。传统的下垂控制，得到有功功率和无功功率的表达式，取其微分形式，改变频率对有功产生较大影响，改变电压对无功产生较大影响；本实施例的新型下垂控制，可以通过逆变器向母线注入的有功和无功的表达式可以发现，取其微分形式，可以发现改变频率对无功功率产生影响，改变电压对有功功率产生影响，这种改变一个量只对一个量产生影响，实现解耦过程。

本实施例的逆变器模块的控制系统包括p-u/q-f下垂控制器、电压电流环控制环节及接口电路的lc滤波器和耦合电抗；

步骤5、通过park变换将三相静止坐标系转化为dq两相同步旋转坐标系，将各模块在各自的坐标系上建立小扰动模型，选择参考坐标系，再通过坐标变换，得到逆变器完整的小扰动模型，包括：

控制单元的小扰动模型为：

其中，下垂控制线性化的表达式为：

其中,分别表示电压和频率偏差的求导，分别表示逆变器有功功率和无功功率偏差的求导值。

步骤6、建立下垂控制器的小扰动模型，表达式为：

式中，表示控制器在旋转坐标系下的输出到电压电流环的电压参考值偏差的导数，δudqref表示表示控制器输出到电压电流环的电压参考值的偏差，δp，δq分别表示逆变器有功功率的偏差和无功功率的偏差，δil1dq表示滤波电路的电感电流偏差，δucdq、δil2dq表示输入功率计算模块的电压偏差和电流的偏差,a1，a2，a3，a4表示系数矩阵，其中下文下标为dq的量均表示类似列向量。

步骤7、构建电压电流环的小扰动模型为：

式中，δuidqref表示经过电压电流环输出到逆变单元的电压偏差，δidq，δudq分别表示正交派克变换后所得的d/q轴电压偏差和电流偏差，δf表示系统频率偏差。a5，a6，a7，a8，a9表示系数矩阵。

步骤8、构建接口电路的小扰动模型为：

式中，和分别表示滤波电容的电压偏差的导数、电流偏差的导数和耦合电感的电流偏差的导数，δuidq为逆变器输出的电压偏差。δufdq为馈线上的电压偏差，a10，a11，a12，a13表示系数矩阵。令δuidqref＝δuidq，模型(6)可化简。

步骤9、联立模型(3)、(4)、(5)、(6)可得逆变器的整体小扰动模型，即为：

式中，δx为状态变量，代表状态变量δx的导数

其中

其中k表示逆变器的数量，δpk、δqk分别代表第k个逆变器有功功率和无功功率的偏差，δfk表示第k个逆变器的频率偏差，δukdqref表示第k个控制器在旋转坐标系下的输出到电压电流环的电压参考值偏差,δikdq和δukdq分别表示第k个控制器正交派克变换后所得的d/q轴电压和电流偏差值，δuckdq和δil2kdq分别表示第k个逆变器输入功率计算模块的电压和电流的偏差，δil1kdq表示第k个逆变器的滤波电路的电感电流偏差。δufkdq表示第k个逆变器馈线上的电压偏差，下标1,2…k表示第k个逆变器与网络馈线的连接节点，a、b表示系数矩阵。

步骤10、综合网络及负荷的小扰动模型(1)、(2)和低压下垂控制的逆变器的小扰动数学模型(7)，可得微电网的整体小扰动数学模型，为：

式中，am为系统状态矩阵。

状态矩阵am特征值对矩阵元素的灵敏度可表示为：

对于特征值的灵敏度抽取典型分析对系统的影响。

步骤11、根据所提出的微电网的小扰动数学模型，计算一些提供的初始参数和系统参数状态下状态矩阵am的特征值，根据系统的状态特征值分布；

如果特征值的实部均为负，表明微网在小扰动下呈现衰减振荡，能保持小扰动下的稳定；

若有特征值的实部为正，则系统不稳定，对于不稳定的特征根，对其进行灵敏度分析。

特征值的灵敏度可以反映参数变化引起相应特征值变化的大小，从而反映参数变化对系统扰动相应特性的影响程度，通过计算特征根灵敏度，得出电压电流环的比例系数对系统小扰动的稳定性的影响较大，积分系数影响较小，其结果与特征值分析一致。通过特征值灵敏度的分析结果可为控制器参数的优化提供依据，通过改变相对应的参数大小，再次计算特征根分析系统的稳定性。

本发明的实施例的上述描述是为了示例和说明的目的而给出的。它们并不是穷举性，也不意于将本发明限制于这些精确描述的内容，在上述教导的指引下，还可以有许多改动和变化。这些实施例被选中和描述仅是为了最好解释本发明的原理以及它们的实际应用，从而使得本领域技术人员能够更好地在各种实施例中并且使用适合于预期的特定使用的各种改动来应用本发明。因此，应当理解的是，本发明意欲覆盖在下面权利要求范围内的所有改动和等同。