技术领域
[0001] 本
发明涉及
智能电网同步控制领域的电压参数鲁棒估计方法，具体地，涉及一种基于现代控制理论、滑模观测器和线性矩阵不等式的电压参数鲁棒估计方法。
背景技术
[0002] 根据国家
能源战略，大规模分布式
可再生能源得到大
力发展，并网速度加快。分布式能源并网的核心问题是同步控制，而实时、准确地获取电网电压参数(幅值、
频率和
相位)则是实现同步控制的关键。此外，实时电压参数还是多逆变器并联运行、电压
质量控制和柔性输电所必需的信息。
[0003] 电网电压参数估计，是指利用电压
信号的高频
采样值和适当的辨识
算法，实时获取正弦电压信号的3个参数，即幅值、频率和相位。由于实际电网运行中存在多种干扰，如电压跌落和低频振荡，电压参数也会随之发生变化。滑模观测器具有响应速度快、鲁棒性强、计算量小等优点，已成功应用于
电机转速和锂
电池电量估计等领域。
[0004] 对
现有技术的检索发现，中国
专利申请号CN201310602875.4，公开日2014-2-19，提出了一种基于参考频率的电网
同步信号实时检测方法，该方法对经典
锁相环方法进行了改进，在低通滤波环节引入
负反馈和比例控制，解决了同步信号由于频率突变导致的相
角偏移问题。然而，该方法本质上还是采用的
锁相环原理，该原理的优点是设计简单、易于实现，但其本征缺点是响应速度慢(100毫秒级)，不能准确估计低频振荡情况下所需的、响应速度为毫秒级的电压参数。在电力系统实际应用中，大多采用经典锁相环技术(PLL)，但其响应速度慢，估计
精度易受电压幅值和频率
波动的影响。因此，现有专利技术和经典锁相环方法均不能满足电压参数估计在快速性、准确性和鲁棒性等方面的高要求。
发明内容
[0005] 针对现有技术
缺陷，本发明提供一种基于滑模观测器的单相电网电压参数鲁棒估计方法，可为提高单相电网电压参数的估计精度和鲁棒性等提供技术
支撑，亦可用于电网电压跌落和低频振荡故障的在线检测。
[0006] 为实现以上目的，本发明提供一种基于滑模观测器的单相电网电压参数的鲁棒估计方法，该方法包括以下步骤：
[0007] 步骤一、根据单相电网正弦电压信号的动态特性，构建其二阶
状态空间模型；
[0008] 步骤二、在上述状态空间模型的
基础上，采用坐标变换，设计滑模观测器，其中，单相电网正弦电压信号的角频率包括时变的情况；
[0009] 步骤三、针对步骤二中角频率时变的情况，对频率系数的观测误差进行重构，获取频率系数的动态观测方程；
[0010] 步骤四、确定滑模观测器参数，以保证滑模观测器的收敛性和良好动态特性；
[0011] 步骤五、依据单相电网电压信号的实时采样数据和步骤四得到的滑模观测器，进行电压参数的在线估计。
[0012] 优选地，步骤一中，将正弦电压信号等效为一个时变的线性二阶系统，并考虑电网频率的变化，构建电压信号的状态空间模型。
[0013] 优选地，所述步骤一，具体如下：
[0014] 不考虑谐波干扰时，单相电网的正弦电压信号表示为：
[0015]
[0016] 其中，t表示时间，Vm为电压幅值，ω为角频率，φ0和 分别表示电压信号的初相位和实时相位，其范围为[0,2π]；
[0017] 根据正弦电压信号的动态特性，构建其状态空间模型如下：
[0018]
[0019] 其中，x和y分别表示状态空间模型的状态变量和可实时测量的输出变量，表示x的一阶倒数，A和C分别表示状态空间模型的系统矩阵和输出矩阵，它们的具体表达式如下：
[0020]
[0021]
[0022] C＝[1 0]
[0023] 式中，ω0表示额定角频率，θ表示频率偏离额定值时的变化系数即频率系数。
[0024] 优选地，步骤二中，采用坐标变换和滑模观测器方法，同时引入观测误差的线性和非线性反馈项，以提高电压参数估计的响应速度和鲁棒性。
[0025] 优选地，所述步骤二，具体如下：
[0026] 采用如下坐标变换：
[0027] z＝Tx (4)
[0028] 其中，T为可逆的坐标变换矩阵，
[0029]
[0030] 式中，x表示状态空间模型的状态变量；k1为待设计的常数参数，ω0表示额定角频率，θ表示频率偏离额定值时的变化系数即频率系数；
[0031] 采用式(4)的坐标变换后，正弦电压信号的状态空间模型变换如下：
[0032]
[0033] 其中，y表示状态空间模型的可实时测量的输出变量，z表示新
坐标系下的状态变量，新坐标系下的系统矩阵和输出矩阵表示如下：
[0034]
[0035]
[0036] 定义观测误差ey为：
[0037]
[0038] 设计滑模观测器为：
[0039]
[0040] v＝-ρ1sign(ey)
[0041]
[0042] 其中， 和 分别表示z、y和θ的观测值，表示 的一阶导数，Gl和Gn分别表示滑模观测器的线性增益和非线性增益矩阵；v为非线性反馈项，ρ1为待设计的反馈系数，sign表示符号函数。
[0043] 优选地，所述步骤三，具体如下：
[0044] 定义式(5)所描述状态空间模型的观测误差如下：
[0045]
[0046] 结合式(5)和式(7)，可得滑模观测器的观测误差方程为：
[0047]
[0048] 式(9)中，
[0049]
[0050]
[0051]
[0053]
[0054] 上述公式中：Ac11、Ac21、Ac22分别表示矩阵Ac的三个分块矩阵， 分别表示观测误差向量z的两个分量，表示 的第一个分量；
[0055] 选择增益矩阵Gl和Gn保证误差观测方程(10)的收敛性，等状态变量进入滑模态时，有
[0056]
[0057] 则通过重构得到频率系数的观测误差：
[0058]
[0059] 其中，Gn2表示非线性增益Gn2的第二个分量，veq表示非线性项的等效控制输入；设计频率系数的观测方程为：
[0060]
[0061] 其中，表示 的一阶导数；ρ2和Г为待设计的滑模观测器参数，用于调节观测器的动态性能，由步骤四确定；
[0062] 式(7)和式(13)一同构成了电压参数的滑模观测器。
[0063] 优选地，步骤四中，采用Lyapunov理论和线性矩阵不等式的方法，获得可保证滑模观测器收敛性和动态性能的滑模观测器参数。
[0064] 优选地，所述步骤四，具体如下：
[0065] 滑模观测器的线性增益矩阵取为：
[0066]
[0067] 其中，k2为待定参数，由步骤四确定，且有k2>k1>0，ω0表示额定角频率；
[0068] 滑模观测器的非线性增益为：
[0069]
[0070] 其中，P-1表示可逆对称阵P的逆矩阵，P由下式决定：
[0071]
[0072] 其中，Ac表示滑模观测器的观测误差矩阵， 为Ac的转置，B为输入矩阵，Q为待求解的对称正定矩阵，利用线性矩阵不等式方法数值求解式(16)，从而确定Gn。
[0073] 优选地，步骤五中，根据滑模观测器的观测结果来计算电网电压的参数，其中：
[0074] 通过坐标逆变换得到状态变量的估计值：
[0075]
[0076] 其中，T为可逆的坐标变换矩阵，T-1为T的逆矩阵：
[0077]
[0078] 电网电压的幅值、频率和相角按下式计算：
[0079]
[0080] 其中：k1为待设计的常数参数，ω0表示额定角频率，θ表示频率偏离额定值时的变化系数即频率系数；Vm为电压幅值，ω为电压角频率，表示电压信号的实时相位，其范围为[0,2π]；f表示电压频率，x1、x2分别表示式(2)中状态向量x的两个分量。
[0081] 与现有技术相比，本发明具有如下的有益效果：
[0082] 本发明基于滑模观测器进行单相电网电压参数的鲁棒估计，采用了滑模观测器、Lyapunov和线性矩阵不等式等方法，鲁棒估计具有响应速度快、鲁棒性、计算量小等优点，可快速准确地实时估计电压参数，还可检测出电网电压跌落和低频振荡等故障。
[0083] 本发明可以为提高单相电网电压参数的实时估计精度和鲁棒性等提供理论方法，对智能电网的实时控制和保护具有重要意义。
附图说明
[0084] 通过阅读参照以下附图对非限制性
实施例所作的详细描述，本发明的其它特征、目的和优点将会变得更明显：
[0085] 图1为本发明一实施例基于滑模观测器的单相电网电压参数估计方法原理图；
[0086] 图2为本发明一实施例在电压幅值阶跃条件下，基于滑模观测器(SMO)和锁相环法(PLL)的电压参数估计结果对比。
[0087] 图3为本发明一实施例在电压相位阶跃条件下，基于滑模观测器(SMO)和锁相环法(PLL)的电压参数估计结果对比。
[0088] 图4为本发明一实施例在低频振荡条件下，基于滑模观测器(SMO)和锁相环法(PLL)的电压参数估计结果对比。
具体实施方式
[0089] 下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干
变形和改进。这些都属于本发明的保护范围。
[0090] 如图1所示，本发明所述的基于滑模观测器的单相电网电压参数鲁棒估计方法原理图，其中所述方法步骤包括：
[0091] 步骤1、根据单相电网电压信号的动态特性，构建其状态空间模型；
[0092] 步骤2、采用适当的坐标变换，进行滑模观测器的设计,其中,单相电网正弦电压信号的角频率可以是时变的；
[0093] 步骤3、针对步骤2中角频率时变的情况，对频率系数的观测误差进行重构，获取频率系数的动态观测方程；
[0094] 步骤4、确定滑模观测器参数以保证滑模观测器的收敛性和良好动态特性；
[0095] 步骤5、依据单相电网电压信号的高频采样数据和所设计的滑模观测器，进行电压参数(幅值、频率和相角)的在线估计。
[0096] 上述步骤1中，与信号滤波和频域分析等常规方法不同，本发明将正弦电压信号等效为一个时变的线性二阶系统，并考虑了电网频率的变化，构建电压信号的状态空间模型。
[0097] 上述步骤2中，采用坐标变换和滑模观测器方法，同时引入观测误差的线性和非线性反馈项，以提高电压参数估计的响应速度和鲁棒性。
[0098] 上述步骤3中，采用故障重构获得频率系数的观测误差动态方程，以解决电网低频振荡的实时检测问题、提高频率的估计精度。
[0099] 上述步骤4中，采用Lyapunov理论和线性矩阵不等式的方法，获得可保证滑模观测器收敛性和动态性能的滑模观测器参数。
[0100] 上述步骤5中，根据滑模观测器的观测结果来计算电网电压的3个参数(电网电压的幅值、频率和相角)，与常规的频域计算方法不同，无需附加任何滤波环节。
[0101] 本发明方法在状态空间模型的基础上进行坐标变换和滑模观测器设计，采用Lyapunov理论和线性矩阵不等式方法，获得滑模观测器参数，采用故障重构获得频率系数观测误差的动态方程，从而得到高精度的电压参数估计。
[0102] 以下提供实施例，对上述方法的具体实施进行详细的说明，以便进一步理解本发明的技术方案。
[0103] 本实施例涉及以民用220V单相配电网为例，额定频率为50Hz，并设初相角为零。
[0104] 本实施例提供基于滑模观测器的单相电网电压参数鲁棒估计方法，该方法具体包括以下步骤：
[0105] 步骤一、根据单相电网正弦电压信号的动态特性，构建其二阶状态空间模型。
[0106] 不考虑谐波干扰时，单相电网的正弦电压信号可表示为：
[0107]
[0108] 其中，t表示时间，Vm为电压幅值，ω为角频率,φ0和 分别表示电压信号的初相位和实时相位，其范围为[0,2π]；
[0109] 根据正弦电压信号的动态特性，构建其状态空间模型如下：
[0110]
[0111] 其中，x和y分别表示状态空间模型的状态变量和可实时测量的输出变量，表示x的一阶倒数，A和C分别表示状态空间模型的系统矩阵和输出矩阵，它们的具体表达式如下：
[0112]
[0113]
[0114] C＝[1 0]
[0115] 式中，ω0表示额定角频率，θ表示频率偏离额定值时的变化系数即频率系数。
[0116] 步骤二、在上述状态空间模型的基础上，采用坐标变换，进行滑模观测器设计。
[0117] 首先，为了便于分析，采用如下坐标变换：
[0118] z＝Tx (4)
[0119] 其中，T为可逆的坐标变换矩阵：
[0120]
[0121] 式中，k1为待设计的常数参数。
[0122] 采用式(4)的坐标变换后，式(1)所描述正弦电压信号的状态空间模型变换如下：
[0123]
[0124] 其中，z表示新坐标系下的状态变量，新坐标系下的系统矩阵和输出矩阵表示如下：
[0125]
[0126]
[0127] 定义滑模观测器的观测误差ey为：
[0128]
[0129] 设计滑模观测器为：
[0130]
[0131]
[0132] 其中， 和 分别表示z、y和θ的观测值，表示 的一阶导数，Gl和Gn分别表示滑模观测器的线性增益和非线性增益矩阵；v为非线性反馈项，ρ1为待设计的反馈系数，sign表示符号函数。
[0133] 步骤三、针对步骤二中频率时变的情况，对频率系数的观测误差进行重构，获取频率系数 动态观测方程。
[0134] 定义滑模观测器的观测误差如下：
[0135]
[0136] 结合式(5)和式(7)，可得滑模观测器的观测误差方程为
[0137]
[0138] 式(9)中，
[0139]
[0140]
[0141]
[0142] 根据矩阵分块理论，式(9)可改写为：
[0143]
[0144] Ac11、Ac21、Ac22分别表示矩阵Ac的三个分块矩阵， 分别表示观测误差向量z的两个分量，表示 的第一个分量。
[0145] 选择合适增益矩阵Gl和Gn可以保证误差观测方程(10)的收敛性，等状态变量进入滑模态时，有
[0146]
[0147] 则可通过重构得到频率系数的观测误差：
[0148]
[0149] 设计频率系数的观测方程为：
[0150]
[0151] 其中，表示 的一阶导数，ρ2和Г为待设计的滑模观测器参数，用于调节观测器的动态性能，由步骤四确定。
[0152] 综上可知，式(7)和式(13)一同构成了电压参数的滑模观测器。
[0153] 步骤四、确定滑模观测器参数以保证滑模观测器的收敛性和良好动态特性。
[0154] 滑模观测器的线性增益矩阵取为：
[0155]
[0156] 其中，k2为待定参数，且有k2>k1>0.
[0157] 滑模观测器的非线性增益为：
[0158]
[0159] 其中，P-1表示可逆对称阵P的逆矩阵，P由下式决定：
[0160]
[0161] 其中，Ac表示滑模观测器的观测误差矩阵， 为Ac的转置，B为输入矩阵，Q为待求解的对称正定矩阵。利用线性矩阵不等式方法(LMI)数值求解式(16)，从而确定Gn＝[0,1]T。
[0162] 实施例中其它参数：ρ1、ρ2和Г均为正数，以可保证滑模观测器的收敛性。
[0163] 步骤五、依据单相电网电压信号的实时采样数据和所设计的滑模观测器，进行电压参数的在线辨识。
[0164] 通过坐标逆变换得到如式(2)所描述的系统状态变量的估计值：
[0165]
[0166] 其中，T-1为T的逆矩阵：
[0167]
[0168] 电网电压的3个参数可按下式计算：
[0169]
[0170] 设计滑模观测器的参数为k1＝5,k2＝8,Г＝2/ω0,ρ1＝1.7×105,ρ2＝3000；作为对比，锁相环的参数设为kP＝20，kI＝5，kD＝0。
[0171] 在Matlab/Simulink
软件环境下，构建单相电网的仿真模型，采用ode45算法进行仿真。图2、图3和图4分别是在不同条件下，基于锁相环和滑模观测器的仿真实验结果对比。
[0172] 本发明基于滑模观测器进行单相电网电压参数的鲁棒估计，采用了滑模观测器、Lyapunov和线性矩阵不等式等方法，鲁棒估计具有响应速度快、鲁棒性、计算量小等优点，可快速准确地实时估计电压参数，还可检测出电网电压跌落和低频振荡等故障；本发明可以为提高单相电网电压参数的实时估计精度和鲁棒性等提供理论方法，对智能电网的实时控制和保护具有重要意义。
[0173] 以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在
权利要求的范围内做出各种变形或
修改，这并不影响本发明的实质内容。
