本发明公开了一种NPC三相三电平并网逆变器的模型预测控制方法及装置。方法为：对三相入网电流和三相电网电压采样值进行Clark变换；利用DSOGI‑SPLL模块和电网电压特征值提取模块得到计算入网参考电流时所必须的电网电压特征值；采用两步预测法对三相入网电流进行跟踪控制和直流侧中点电压抑制，以补偿因采样和计算造成的延迟。该装置包括主功率电路、控制电路和检测电路，主功率电路包括NPC三相三电平逆变器和L型低通滤波器，检测电路包括入网电流检测电路、电网电压检测电路和直流侧中点电压检测电路。本发明能够在非理想电网条件下实现对NPC三相三电平逆变器良好的并网控制，大大提高入网电能品质，同时有效抑制直流侧中点电压的波动。
1.一种NPC三相三电平并网逆变器的模型预测控制方法，其特征在于，采用两步预测法进行模型预测计算，过程包括以下步骤：
步骤S1，输出当前即第k个采样周期最优开关状态Sopt(k)控制逆变器，Sopt(k)由上一个即第k-1个采样周期计算得到；
步骤S2，采样及变换：对入网电流进行检测得到a、b、c三相入网电流ia(k)、ib(k)、ic(k)，并对ia(k)、ib(k)、ic(k)进行Clark变换得到iα(k)、iβ(k)；对电网电压进行检测得到a、b、c三相电网电压ea(k)、eb(k)、ec(k)，并对ea(k)、eb(k)、ec(k)进行Clark变换得到eα(k)、eβ(k)；根据检测所得三相电网电压ea(k)、eb(k)、ec(k)，提取出电网电压特征值并进行入网参考电流计算，得到αβ坐标系下的入网参考电流i*α(k)、i*β(k)；对电容电压进行检测得到直流侧中点电压Δvc(k)；
步骤S3，相角补偿：对eα(k)、eβ(k)做一个采样周期Ts相角补偿得到eα(k+1)、eβ(k+1)，对i*α(k)、i*β(k)做两个采样周期2Ts相角补偿得到步骤S4，计算入网电流反馈值iα(k+1)、iβ(k+1)：结合Sopt(k)对应的逆变器交流输出侧电压矢量值uα(k)、uβ(k)与步骤S2的采样及变换结果，根据逆变器的离散数学模型计算出第k+1个采样周期入网电流反馈值iα(k+1)、iβ(k+1)；
步骤S5，计算直流侧中点电压反馈值Δvc(k+1)：结合Sopt(k)对应的相开关函数状态Sa、Sb、Sc与步骤S2的采样及变换结果，根据逆变器的离散数学模型计算出第k+1个采样周期直流侧中点电压反馈值Δvc(k+1)；
步骤S6，对第k+1个采样周期入网电流反馈值iα(k+1)、iβ(k+1)进行反Clark变换，得到ia(k+1)、ib(k+1)、ic(k+1)；
步骤S7，遍历计算：结合iα(k+1)、iβ(k+1)、eα(k+1)、eβ(k+1)、Δvc(k+1)、ia(k+1)、ib(k+1)、ic(k+1)、逆变器交流输出侧电压矢量值以及该电压矢量值对应的相开关函数状态，根据逆变器的离散数学模型预测计算出第k+2个采样周期的入网电流iα(k+2)(i)、iβ(k+2)(i)和直流侧中点电压Δvc(k+2)(i)；
步骤S8，建立目标函数g，作为选出27组相开关函数状态中的最优开关状态Sopt(k+1)的依据；
步骤S9，存储目标函数g最小值对应的开关状态，该状态即最优开关状态Sopt(k+1)，将Sopt(k+1)做为第k+1个采样周期的最优输出；
步骤S10，等待本次采样周期时间结束，返回步骤S1，进入下一个循环。
2.根据权利要求1所述的NPC三相三电平并网逆变器的模型预测控制方法，其特征在于，步骤S2所述根据检测所得三相电网电压ea(k)、eb(k)、ec(k)，提取出电网电压特征值并进行参考电流计算，得到αβ坐标系下的入网参考电流i*α(k)、i*β(k)，具体如下：+ - 5-
(2.1)非理想电网电压中含有基频正序分量e、基频负序分量e、5倍频负序分量e 、7倍频正序分量e7+，因此假设：
其中，下标a、b、c表示a相、b相、c相，上标+、-、5-、7+分别表示基频正序分量、基频负序分量、5倍频负序分量、7倍频正序分量；
(2.2)定义：基频正序旋转坐标系dq+以相对于αβ坐标系中α轴的角速度ω逆时针旋转；
基频负序旋转坐标系dq-以相对于αβ坐标系中α轴的角速度ω顺时针旋转；5倍频负序旋转坐标系dq5-以相对于αβ坐标系中α轴的角速度5ω顺时针旋转；7倍频正序旋转坐标系dq7+以相对于αβ坐标系中α轴的角速度7ω逆时针旋转；其中ω为电网电压基波角频率；
(2.3)对步骤(2.1)中公式两端进行Clark变换，得到下式：将α轴量视为向量实部，将β轴量视为向量虚部，上式转化为：根据Clark变换的特性， 以相对于αβ坐标系中α轴的角速度ω逆时针旋转， 以相对于αβ坐标系中α轴的角速度ω顺时针旋转， 以相对于αβ坐标系中α轴的角速度5ω顺时针旋转， 以相对于αβ坐标系中α轴的角速度7ω逆时针旋转；
结合步骤(2.2)和步骤(2.3)， 与基频正序旋转坐标系dq+相对静止， 与基频负序旋转坐标系dq-相对静止， 与5倍频负序旋转坐标系dq5-相对静止， 与7倍频正序旋转坐标系dq7+相对静止；
(2.4)定义：电网电压特征值 分别为 在基频正序旋转坐标系dq+中d+轴、q+轴上的投影，电网电压特征值 分别为 在基频负序旋转坐标系dq-中d-轴、q-轴上的投影，电网电压特征值 分别为 在5倍频负序旋转坐标系dq5-中d5-轴、q5-轴上的投影，电网电压特征值 分别为 在7倍频正序旋转坐标系dq7+中d7+轴、q7+轴上的投影；
同理于步骤(2.1)至(2.4)，得到入网电流特征值(2.5)根据瞬时无功理论，入网瞬时复功率 为：
其中，j为复数单位，p(t)为入网瞬时有功功率，q(t)为入网瞬时无功功率，eα、eβ分别为三相电网电压ea、eb、ec进行Clark变换后的α轴分量、β轴分量，ia、iβ分别为三相入网电流ia、ib、ic进行Clark变换后的α轴分量、β轴分量；
根据上式可得到入网瞬时有功功率p(t)、入网瞬时无功功率q(t)：p(t)＝P0+Pc2cos(2ωt)+Ps2sin(2ωt)+Pc4cos(4ωt)+Ps4sin(4ωt)+Pc6cos(6ωt)+Ps6sin(6ωt)+Pc8cos(8ωt)+Ps8sin(8ωt)+Pc12cos(12ωt)+Ps12sin(12ωt)其中，Pc2、Ps2、Pc4、Ps4、Pc6、Ps6、Pc8、Ps8、Pc12、Ps12分别是入网瞬时有功功率中含有的高次谐波对应的系数，称为有功系数；P0是入网瞬时有功功率中的直流分量；
q(t)＝Q0+Qc2cos(2ωt)+Qs2sin(2ωt)+Qc4cos(4ωt)+Qs4sin(4ωt)+Qc6cos(6ωt)+Qs6sin(6ωt)+Qc8cos(8ωt)+Qs8sin(8ωt)+Qc12cos(12ωt)+Qs12sin(12ωt)其中，Qc2、Qs2、Qc4、Qs4、Qc6、Qs6、Qc8、Qs8、Qc12、Qs12分别是入网瞬时无功功率中含有的高次谐波对应的系数，称为无功系数；Q0是入网瞬时无功功率中的直流分量；
所述有功系数、无功系数、电网电压特征值、入网电流特征值有如下关系：反解上述矩阵方程，得到入网电流特征值 对
Park变换得到αβ坐标系下的入网参考电流i*α(k)、i*β(k)。
3.根据权利要求1所述的NPC三相三电平并网逆变器的模型预测控制方法，其特征在于，步骤S3所述相角补偿：对eα(k)、eβ(k)做一个采样周期Ts相角补偿得到eα(k+1)、eβ(k+1)，对i*α(k)、i*β(k)做两个采样周期2Ts相角补偿得到 具体公式如下：其中，ω为电网电压基波角频率。
4.根据权利要求1所述的NPC三相三电平并网逆变器的模型预测控制方法，其特征在于，步骤S4所述计算入网电流反馈值iα(k+1)、iβ(k+1)：结合Sopt(k)对应的逆变器交流输出侧电压矢量值uα(k)、uβ(k)与步骤S2的采样及变换结果，根据逆变器的离散数学模型计算出第k+1个采样周期入网电流反馈值iα(k+1)、iβ(k+1)，具体公式如下：其中，L为滤波电感的感值，R为逆变器桥臂电阻和滤波电感电阻折合后等效电阻的阻值；iα(k)、iβ(k)为第k个采样周期，三相入网电流的采样值经过Clark变换后的值；eα(k)、eβ(k)为第k个采样周期，电网电压的采样值经过Clark变换后的值。
5.根据权利要求1所述的NPC三相三电平并网逆变器的模型预测控制方法，其特征在于，步骤S5所述计算直流侧中点电压反馈值，结合Sopt(k)对应的相开关函数状态Sa、Sb、Sc与步骤S2的采样及变换结果，计算出第k+1个采样周期直流侧中点电压反馈值Δvc(k+1)，具体公式如下：
6.根据权利要求1所述的NPC三相三电平并网逆变器的模型预测控制方法，其特征在于，步骤S7所述遍历计算，结合iα(k+1)、iβ(k+1)、eα(k+1)、eβ(k+1)、Δvc(k+1)、ia(k+1)、ib(k+1)、ic(k+1)、逆变器交流输出侧电压矢量值以及该电压矢量值对应的相开关函数状态，根据逆变器的离散数学模型预测计算第k+2个采样周期的入网电流iα(k+2)(i)、iβ(k+2)(i)和直流侧中点电压Δvc(k+2)(i)，具体公式如下：其中，i＝1、2……27，uα(i)、uβ(i)表示第i组逆变器交流输出侧电压矢量值，Sa(i)、Sb(i)、Sc(i)表示第i组逆变器交流输出侧电压矢量值对应的相开关函数状态；C为直流侧电容的容值，L为滤波电感的感值，R为逆变器桥臂电阻和滤波电感电阻折合后等效电阻的阻值。
7.根据权利要求1所述的NPC三相三电平并网逆变器的模型预测控制方法，其特征在于，步骤S8所述建立目标函数g，作为选出27组相开关函数状态中的最优开关状态Sopt(k+1)的依据，具体公式如下：
其中，i＝1、2……27，λdc为直流侧中点电压权重系数。
NPC三相三电平并网逆变器的模型预测控制方法及装置
技术领域
[0001] 本发明涉及电能变换装置的直流-交流变换器技术领域，特别是一种NPC三相三电平并网逆变器的模型预测控制方法及装置。
背景技术
[0002] NPC(Neutral Point Clamped)三相三电平并网逆变器在分布式并网发电系统中起着能量转换接口的作用，其工作状态对进入电网的电能品质影响很大。在理想电网条件下，电网电压中只含有基频正序分量e+，现有的NPC三相三电平并网逆变器控制方法相对成熟，如PI控制、PR控制、滞环控制等。
[0003] 然而实际中的电网条件绝大部分是非理想的，非理想电网电压中除了基频正序分量e+外，还含有基频负序分量e-、5倍频负序分量e5-和7倍频正序分量e7+等，此时基于理想电网条件下的NPC三相三电平并网逆变器控制方法必然导致低品质电能进入电网。与此同时，基于理想电网条件下的NPC三相三电平并网逆变器直流侧中点电压抑制方法也不能良好地抑制直流侧中点电压的波动，易造成开关器件损坏、三相入网电流THD值较大等问题。
发明内容
[0004] 本发明的目的在于提供一种能够在非理想电网条件下，对NPC三相逆变器进行良好并网控制的NPC三相三电平并网逆变器模型预测控制方法及装置，以提高入网电能的品质，并有效抑制直流侧中点电压的波动。
[0005] 实现本发明目的的技术解决方案为：一种NPC三相三电平并网逆变器的模型预测控制方法，采用两步预测法进行模型预测计算，过程包括以下步骤：
[0006] 步骤S1，输出当前即第k个采样周期最优开关状态Sopt(k)控制逆变器，Sopt(k)由上一个即第k-1个采样周期计算得到；
[0007] 步骤S2，采样及变换：对入网电流进行检测得到a、b、c三相入网电流ia(k)、ib(k)、ic(k)，并对ia(k)、ib(k)、ic(k)进行Clark变换得到iα(k)、iβ(k)；对电网电压进行检测得到a、b、c三相电网电压ea(k)、eb(k)、ec(k)，并对ea(k)、eb(k)、ec(k)进行Clark变换得到eα(k)、eβ(k)；根据检测所得三相电网电压，提取出电网电压特征值并进行入网参考电流计算，得到αβ坐标系下的入网参考电流i*α(k)、i*β(k)；对电容电压进行检测得到直流侧中点电压Δvc(k)；
[0008] 步骤S3，相角补偿：对eα(k)、eβ(k)做一个采样周期Ts相角补偿得到eα(k+1)、eβ(k+1)，对i*α(k)、i*β(k)做两个采样周期2Ts相角补偿得到
[0009] 步骤S4，计算入网电流反馈值：结合Sopt(k)对应的逆变器交流输出侧电压矢量值uα(k)、uβ(k)与步骤S2的采样及变换结果，根据逆变器的离散数学模型计算出第k+1个采样周期入网电流反馈值iα(k+1)、iβ(k+1)；
[0010] 步骤S5，计算直流侧中点电压反馈值：结合Sopt(k)对应的相开关函数状态Sa、Sb、Sc与步骤S2的采样及变换结果，根据逆变器的离散数学模型计算出第k+1个采样周期直流侧中点电压反馈值Δvc(k+1)；
[0011] 步骤S6，对第k+1个采样周期入网电流反馈值iα(k+1)、iβ(k+1)进行反Clark变换，得到ia(k+1)、ib(k+1)、ic(k+1)；
[0012] 步骤S7，遍历计算：结合iα(k+1)、iβ(k+1)、eα(k+1)、eβ(k+1)、Δvc(k+1)、ia(k+1)、ib(k+1)、ic(k+1)、逆变器交流输出侧电压矢量值以及该电压矢量值对应的相开关函数状态，根据逆变器的离散数学模型预测计算出第k+2个采样周期的入网电流iα(k+2)(i)、iβ(k+2)(i)和直流侧中点电压Δvc(k+2)(i)；
[0013] 步骤S8，建立目标函数g，作为选出27组相开关函数状态中的最优开关状态Sopt(k+1)的依据；
[0014] 步骤S9，存储目标函数g最小值对应的开关状态，该状态即最优开关状态Sopt(k+1)，将Sopt(k+1)做为第k+1个采样周期的最优输出；
[0015] 步骤S10，等待本次采样周期时间结束，返回步骤S1，进入下一个循环。
[0016] 进一步地，步骤S2所述根据检测所得三相电网电压，提取出电网电压特征值并进行参考电流计算，得到αβ坐标系下的入网参考电流i*α(k)、i*β(k)，具体如下：
[0017] (2.1)非理想电网电压中含有基频正序分量e+、基频负序分量e-、5倍频负序分量e5-、7倍频正序分量e7+，因此假设：
[0018]
[0019]
[0020]
[0021] 其中，下标a、b、c表示a相、b相、c相，上标+、-、5-、7+分别表示基频正序分量、基频负序分量、5倍频负序分量、7倍频正序分量；
[0022] (2.2)定义：基频正序旋转坐标系dq+以相对于αβ坐标系中α轴的角速度ω逆时针旋转；基频负序旋转坐标系dq-以相对于αβ坐标系中α轴的角速度ω顺时针旋转；5倍频负序旋转坐标系dq5-以相对于αβ坐标系中α轴的角速度5ω顺时针旋转；7倍频正序旋转坐标系dq7+以相对于αβ坐标系中α轴的角速度7ω逆时针旋转；其中ω为电网电压基波角频率；
[0023] (2.3)对步骤(2.1)中公式两端进行Clark变换，得到下式：
[0024]
[0025] 将α轴量视为向量实部，将β轴量视为向量虚部，上式转化为：
[0026]
[0027] 根据Clark变换的特性， 以相对于αβ坐标系中α轴的角速度ω逆时针旋转，以相对于αβ坐标系中α轴的角速度ω顺时针旋转， 以相对于αβ坐标系中α轴的角速度5ω顺时针旋转， 以相对于αβ坐标系中α轴的角速度7ω逆时针旋转；
[0028] 结合步骤(2.2)和步骤(2.3)， 与基频正序旋转坐标系dq+相对静止， 与基频负序旋转坐标系dq-相对静止， 与5倍频负序旋转坐标系dq5-相对静止， 与7倍频正序
旋转坐标系dq7+相对静止；
[0029] (2.4)定义：电网电压特征值 分别为 在基频正序旋转坐标系dq+中d+轴、q+轴上的投影，电网电压特征值 分别为 在基频负序旋转坐标系dq-中d-轴、q-
轴上的投影，电网电压特征值 分别为 在5倍频负序旋转坐标系dq5-中d5-轴、q5-
轴上的投影，电网电压特征值 分别为 在7倍频正序旋转坐标系dq7+中d7+轴、q7+
轴上的投影；
[0030] 同理于步骤(2.1)至(2.4)，得到入网电流特征值
[0031] (2.5)根据瞬时无功理论，入网瞬时复功率 为：
[0032]
[0033] 其中，j为复数单位，p(t)为入网瞬时有功功率，q(t)为入网瞬时无功功率，eα、eβ分别为三相电网电压ea、eb、ec进行Clark变换后的α轴分量、β轴分量，ia、iβ分别为三相入网电流ia、ib、ic进行Clark变换后的α轴分量、β轴分量；
[0034] 根据上式可得到入网瞬时有功功率p(t)、入网瞬时无功功率q(t)：
[0035] p(t)＝P0+Pc2cos(2ωt)+Ps2sin(2ωt)+Pc4cos(4ωt)+Ps4sin(4ωt)
[0036] +Pc6cos(6ωt)+Ps6sin(6ωt)+Pc8cos(8ωt)+Ps8sin(8ωt)
[0037] +Pc12cos(12ωt)+Ps12sin(12ωt)
[0038] 其中，Pc2、Ps2、Pc4、Ps4、Pc6、Ps6、Pc8、Ps8、Pc12、Ps12分别是入网瞬时有功功率中含有的高次谐波对应的系数，称为有功系数；P0是入网瞬时有功功率中的直流分量；
[0039] q(t)＝Q0+Qc2cos(2ωt)+Qs2sin(2ωt)+Qc4cos(4ωt)+Qs4sin(4ωt)
[0040] +Qc6cos(6ωt)+Qs6sin(6ωt)+Qc8cos(8ωt)+Qs8sin(8ωt)
[0041] +Qc12cos(12ωt)+Qs12sin(12ωt)
[0042] 其中，Qc2、Qs2、Qc4、Qs4、Qc6、Qs6、Qc8、Qs8、Qc12、Qs12分别是入网瞬时无功功率中含有的高次谐波对应的系数，称为无功系数；Q0是入网瞬时无功功率中的直流分量；
[0043] 所述有功系数、无功系数、电网电压特征值、入网电流特征值有如下关系：
[0044]
[0045]
[0046] 反解上述矩阵方程，得到入网电流特征值对应的入网电流参考特征值 再
* *
对入网电流参考特征值进行反Park变换得到αβ坐标系下的入网参考电流iα(k)、iβ(k)。
[0047] 进一步地，步骤S3所述相角补偿：对eα(k)、eβ(k)做一个采样周期Ts相角补偿得到eα(k+1)、eβ(k+1)，对i*α(k)、i*β(k)做两个采样周期2Ts相角补偿得到 具体公式如下：
[0048]
[0049]
[0050] 其中，ω为电网电压基波角频率。
[0051] 进一步地，步骤S4所述计算入网电流反馈值：结合Sopt(k)对应的逆变器交流输出侧电压矢量值uα(k)、uβ(k)与步骤S2的采样及变换结果，根据逆变器的离散数学模型计算出第k+1个采样周期入网电流反馈值iα(k+1)、iβ(k+1)，具体公式如下：
[0052]
[0053] 其中，L为滤波电感的感值，R为逆变器桥臂电阻和滤波电感电阻折合后等效电阻的阻值；iα(k)、iβ(k)为第k个采样周期，三相入网电流的采样值经过Clark变换后的值；eα(k)、eβ(k)为第k个采样周期，电网电压的采样值经过Clark变换后的值。
[0054] 进一步地，步骤S5所述计算直流侧中点电压反馈值，结合Sopt(k)对应的相开关函数状态Sa、Sb、Sc与步骤S2的采样及变换结果，计算出第k+1个采样周期直流侧中点电压反馈值Δvc(k+1)，具体公式如下：
[0055]
[0056] 进一步地，步骤S7所述遍历计算，结合iα(k+1)、iβ(k+1)、eα(k+1)、eβ(k+1)、Δvc(k+1)、ia(k+1)、ib(k+1)、ic(k+1)、逆变器交流输出侧电压矢量值以及该电压矢量值对应的相开关函数状态，根据逆变器的离散数学模型预测计算第k+2个采样周期的入网电流iα(k+2)(i)、iβ(k+2)(i)和直流侧中点电压Δvc(k+2)(i)，具体公式如下：
[0057]
[0058]
[0059] 其中，i＝1、2……27，uα(i)、uβ(i)表示第i组逆变器交流输出侧电压矢量值，Sa(i)、Sb(i)、Sc(i)表示第i组逆变器交流输出侧电压矢量值对应的相开关函数状态。
[0060] 进一步地，步骤S8所述建立目标函数g，作为选出27组相开关函数状态中的最优开关状态Sopt(k+1)的依据，具体公式如下：
[0061]
[0062] 其中，i＝1、2……27，λdc为直流侧中点电压权重系数。
[0063] 一种NPC三相三电平并网逆变器的模型预测控制装置，包括主功率电路、控制电路和检测电路，所述主功率电路包括输入电压源Vdc、NPC三相三电平逆变器、L型低通滤波器和三相电网ea、eb、ec，其中输入电压源Vdc与NPC三相三电平逆变器输入端连接，NPC三相三电平逆变器输出端与L型低通滤波器输入端连接，L型低通滤波器输出端与三相电网ea、eb、ec连接；
[0064] 所述的检测电路包括入网电流检测电路、电网电压检测电路、直流侧中点电压检测电路，其中入网电流检测电路输入端与三相电网连接，入网电流检测电路输出端与控制电路第一输入端连接，电网电压检测电路输入端与三相电网连接，电网电压检测电路输出端与控制电路第二输入端连接，直流侧中点电压检测电路输入端与逆变器直流侧电容连接，直流侧中点电压检测电路输出端与控制电路第三输入端连接；
[0065] 所述控制电路包括电网电压特征值提取模块、锁相环模块、参考电流计算模块和预测计算模块，其中锁相环模块为基于双广义二阶积分器软件锁相环；锁相环模块输入端与电网电压检测电路输出端连接，锁相环模块输出端与电网电压特征值提取模块的一个输入端连接，电网电压特征值提取模块的另一个输入端与电网电压检测电路输出端连接，电网电压特征值提取模块的输出端与参考电流计算模块输入端连接，参考电流计算模块输出端接入预测计算模块。
[0066] 进一步地，所述锁相环模块对三相电网电压基频正序分量进行跟踪，检测出三相电网电压基频正序分量的瞬时相位；电网电压特征值提取模块根据三相电网电压基频正序分量的瞬时相位，提取出计算入网参考电流时所需的电网电压特征量。
[0067] 进一步地，所述控制电路采用DSP芯片TMS320F28335。
[0068] 本发明与现有技术相比，其显著优点在于：(1)能够在非理想电网条件下对NPC三相三电平逆变器进行良好的并网控制，大大提高入网电能品质，同时有效地抑制直流侧中点电压的波动；(2)方法简单可靠，易于数字实现，仅需在软件算法上做相应的修改就能实现多种不同的并网控制目标。
附图说明
[0069] 图1是NPC三相三电平并网逆变器主功率电路示意图。
[0070] 图2是αβ坐标系下逆变器交流输出侧a、b、c相对直流侧中点O的电压矢量分布示意图。
[0071] 图3是αβ坐标系、旋转向量、旋转坐标系之间的相对关系示意图。
[0072] 图4是旋转向量在对应的旋转坐标系的坐标轴上的投影示意图，其中(a)是 在+ -
基频正序旋转坐标系dq 坐标轴上的投影示意图，(b)是 在基频负序旋转坐标系dq坐标
轴上的投影示意图，(c)是 在5倍频负序旋转坐标系dq5-坐标轴上的投影示意图，(d)是
在7倍频正序旋转坐标系dq7+坐标轴上的投影示意图。
[0073] 图5是电网电压特征值提取示意图。
[0074] 图6是两步预测法基本原理示意图。
[0075] 图7是预测计算流程图。
[0076] 图8是基于模型预测控制的非理想电网条件下NPC三相三电平逆变器并网系统示意图。
[0077] 图9是实施例1的仿真结果图，其中(a)是实施例1对应的参考入网电流与实际入网电流示意图，(b)是实施例1对应的直流侧中点电压示意图，(c)是实施例1对应的入网瞬时有功功率示意图，(d)是实施例1对应的入网瞬时无功功率示意图。
[0078] 图10是实施例2的仿真结果图，其中(a)是实施例2对应的参考入网电流与实际入网电流示意图，(b)是实施例2对应的直流侧中点电压示意图，(c)是实施例2对应的入网瞬时有功功率示意图，(d)是实施例2对应的入网瞬时无功功率示意图。
[0079] 图11是实施例3的仿真结果图，其中(a)是实施例3对应的参考入网电流与实际入网电流示意图，(b)是实施例3对应的直流侧中点电压示意图，(c)是实施例3对应的入网瞬时有功功率示意图，(d)是实施例3对应的入网瞬时无功功率示意图。
具体实施方式
[0080] 以下结合附图及具体实施例对本发明做进一步详细说明。
[0081] 1、NPC三相三电平并网逆变器离散数学模型
[0082] 图1是NPC三相三电平并网逆变器主功率电路，系统采用三相三线制接法，系统中没有零序电流通路，零序分量不参与瞬时功率计算过程，所以可以不考虑三相电网电压中的零序分量，又考虑到实际非理想电网中主要含有基频正序分量e+、基频负序分量e-、5倍频负序分量e5-、7倍频正序分量e7+，本发明做出以下合理假设：
[0083]
[0084] 其中，下标a、b、c代表a相、b相、c相，上标“+”、“-”、“5-”、“7+”分别代表基频正序分量、基频负序分量、5倍频负序分量、7倍频正序分量，本发明同时合理假设：直流侧电容C1＝C2＝C，且C足够大，近似认为直流侧电容电压 三相滤波电感的感值La＝Lb＝Lc＝L以及逆变器交流侧的等效阻值Ra＝Rb＝Rc＝R。
[0085] 定义1：相开关函数
[0086]
[0087] 其中：i＝a、b、c，Si＝1记为状态P，Si＝0记为状态O，Si＝-1记为状态N。
[0088] 则逆变器交流输出侧(a、b、c)相对直流侧中点(O)的电压：
[0089]
[0090] 利用Clark变换：
[0091]
[0092] 得到αβ坐标系下逆变器交流输出侧(a、b、c)相对直流侧中点(O)的电压矢量分布，如图2所示，三相三电平逆变器有33＝27个开关状态，对应输出27个电压矢量，19种不同的电压矢量，图2中“OPN”表示Sa＝0，Sb＝1，Sc＝-1，其余依此类推。
[0093] 根据基尔霍夫电压定律可得到逆变器交流输出侧电压平衡方程：
[0094]
[0095] 其中，vno为电网电压中性点(n)相对直流侧中点(O)的电压，对(4)式两端进行Clark变换，可以得到αβ坐标系下逆变器交流输出侧电压平衡方程：
[0096]
[0097] 取采样周期为Ts，利用一阶前向差分方程(6)式：
[0098]
[0099] 得到αβ坐标系下，αβ坐标系下入网电流的离散数学模型：
[0100]
[0101] 式中，iα(k)、iβ(k)为第k个采样周期，三相入网电流的采样值经过Clark变换后的值，uα(k)、uβ(k)为第k个采样周期，图2所示的不同电压矢量在αβ坐标系下的坐标值，eα(k)、eβ(k)为第k个采样周期，电网电压的采样值经过Clark变换后的值，iα(k+1)、iβ(k+1)为αβ坐标系下第k+1个采样周期入网电流反馈值。
[0102] 根据基尔霍夫电流定律可得直流侧中点(O)电流平衡方程：
[0103] io＝ic1-ic2 (8)
[0104] 其中：
[0105]
[0106]
[0107]
[0108] 令直流侧中点电压Δvc＝vc1-vc2，同时将(9)、(10)、(11)式带入(8)式，可以得到：
[0109]
[0110] 取采样周期为Ts，利用一阶前向差分方程(13)式：
[0111]
[0112] 得到自然坐标系下，直流侧中点电压的离散数学模型：
[0113]
[0114] 式中，Δvc(k)为自然坐标系下第k个采样周期直流侧中点电压采样值，Δvc(k+1)为自然坐标系下第k+1个采样周期直流侧中点电压反馈值。
[0115] 2、并网控制目标及参考电流
[0116] 定义2：基频正序旋转坐标系dq+以相对于αβ坐标系中α轴的角速度ω逆时针旋转；基频负序旋转坐标系dq-以相对于αβ坐标系中α轴的角速度ω顺时针旋转；5倍频负序旋转坐标系dq5-以相对于αβ坐标系中α轴的角速度5ω顺时针旋转；7倍频正序旋转坐标系dq7+以相对于αβ坐标系中α轴的角速度7ω逆时针旋转，如图3所示。
[0117] 对(1)式两端做Clark变换，可以得到下式：
[0118]
[0119] 将α轴量视为向量实部，将β轴量视为向量虚部，(15)式即：又由Clark变换性质可知， 以相对于αβ坐标系中α轴的角速度ω逆时针旋转， 以相对于αβ坐标系中α轴的角速度ω顺时针旋转， 以相对于αβ坐标系中α轴的角速度5ω顺时针旋转， 以相对于αβ坐标系中α轴的角速度7ω逆时针旋转，如图3所示，同理对于三相入网电流有
[0120] 结合定义2可得， 与基频正序旋转坐标系dq+相对静止， 与基频负序旋转坐标系dq-相对静止， 与5倍频负序旋转坐标系dq5-相对静止， 与7倍频正序旋转坐标系dq7+相对静止。
[0121] 定义3：电网电压特征值 分别为 在基频正序旋转坐标系dq+中d+轴、q+轴上的投影，电网电压特征值 分别为 在基频负序旋转坐标系dq-中d-轴、q-轴上的
投影，电网电压特征值 分别为 在5倍频负序旋转坐标系dq5-中d5-轴、q5-轴上的
投影，电网电压特征值 分别为 在7倍频正序旋转坐标系dq7+中d7+轴、q7+轴上的
投影；同理定义入网电流特征值：
[0122] 由于旋转向量 与基频正序旋转坐标系dq+、旋转向量 与基频负序旋转坐标系- 5- 7+
dq、旋转向量 与5倍频负序旋转坐标系dq 、旋转向量 与7倍频正序旋转坐标系dq 一
一相对静止，旋转向量在与之相对静止的旋转坐标系坐标轴上的投影量是直流量，如图4(a)～(d)所示。将d轴投影量视为向量实部，q轴投影量视为向量虚部，可以得到下式：
[0123]
[0124]
[0125]
[0126]
[0127] 其中，e-jωt是顺时针基频旋转因子，ejωt是逆时针基频旋转因子，ej5ωt是逆时针5倍频旋转因子，e-j7ωt是顺时针7倍频旋转因子。
[0128] 将(16)、(17)、(18)、(19)式整理并带入(15)式，可以得到：
[0129]
[0130] 同理将三相入网电流ia、ib、ic整理成(20)式的形式，如(21)式所示：
[0131]
[0132] 根据瞬时无功理论，入网瞬时复功率：
[0133]
[0134] 其中，p(t)为电网吸收的瞬时有功功率，q(t)为电网吸收的瞬时无功功率， 为的共轭向量，将(22)式的实部和虚部分开并整理，可以得到：
[0135]
[0136]
[0137] 其中，P0是入网瞬时有功功率中的直流分量，Q0是入网瞬时无功功率中的直流分量，Pc2、Ps2、Pc4、Ps4、Pc6、Ps6、Pc8、Ps8、Pc12、Ps12是入网瞬时有功功率中含有的高次谐波对应的系数，称为有功系数；Qc2、Qs2、Qc4、Qs4、Qc6、Qs6、Qc8、Qs8、Qc12、Qs12是入网瞬时无功功率中含有的高次谐波对应的系数，称为无功系数，有功系数、无功系数、电网电压特征值、入网电流特征值有如(25)、(26)式所示的关系：
[0138]
[0139]
[0140] 反解(25)、(26)矩阵方程，可以得到与入网电流特征值相对应的入网电流参考特征值
根据入网有功给定P0*和入网无功给定Q0*，可以确定多种并网控制目标，以
及各个目标对应的入网电流参考特征值的计算公式。
[0141] 并网控制目标1：消除有功2、4、8次波动，即等价于令：P0＝P0*、Q0＝Q0*、Pc2＝0、Ps2＝0、Pc4＝0、Ps4＝0、Pc8＝0、Ps8＝0。
[0142]
[0143] 并网控制目标2：消除无功2、4、8次波动，即等价于令：P0＝P0*、Q0＝Q0*、Qc2＝0、Qs2＝0、Qc4＝0、Qs4＝0、Qc8＝0、Qs8＝0。
[0144]
[0145] 其中，(27)、(28)式子中denA与denB如下：
[0146]
[0147]
[0148] 并网控制目标3：三相入网电流正弦且平衡，即等价于要求入网电流中只有基频正序分量。
[0149]
[0150] 这里有必要指出的是，上述根据不同目标计算出的不同入网电流参考特征值属于dq坐标系下的量，需要进行反Park变换才能得到αβ坐标系下的入网参考电流i*α、i*β。上述入网电流参考特征值计算公式中涉及到的电网电压特征值，均如前文所述，本发明中采用“DSOGI-SPLL+Notch Filter”的组合形式对各个电网电压特征值进行提取，结构如图5所示。
[0151] 3、预测计算
[0152] 模型预测控制属于最优控制范畴，需要定义与被控变量相关的目标函数g作为最优选择的依据，本发明需要控制三相入网电流ia、ib、ic的Clark变换值iα、iβ跟随入网参考电流i*α、i*β以及抑制直流侧中点电压Δvc，定义目标函数g如下：
[0153]
[0154] 其中，λdc为直流侧中点电压权重系数，λdc越大，对Δvc抑制效果越好，电流跟踪效果相对变差；λdc越小，对Δvc抑制效果越差，电流跟踪效果相对变好。
[0155] 考虑到实际采样和计算存在延时，本发明采用两步预测法对延时进行补偿，基本原理结合图6进行说明，x(t)是某一被控变量，x*(t)是该被控变量的参考值，在tk时刻(当前时刻)，第一步：应用最优开关态Sopt(k)控制并网逆变器，Sopt(k)是上一个采样周期的存储值。第二步：基于x(t)的离散预测模型和Sopt(k)，预测计算被控变量x(t)在tk+1时刻的预测值。第三步：将tk+1时刻的计算值作为反馈值，结合三相三电平并网逆变器的27种开关状态与离散预测模型，计算出被控变量x(t)在tk+2时刻的27个预测值，并带入目标函数g的值越小，说明x(t)越接近x*(t)，因此使目标函数值最小的开
关状态即为最优，将该最优开关状态作为tk+1时刻的控制输出，便能够使tk+2时刻的x(t)最接近于tk+2时刻的x*(t)。例如在图6中，在tk+2时刻，开关状态S(4)对应的x(t)最接近x*(t)，目标函数g(x(t)，x*(t))的值最小，因此存储S(4)作为tk+1时刻输出的最优开关状态Sopt(k+1)，下一个采样周期重复以上的步骤，上述步骤具体可用流程图7表示。应当指出的是，为了减小控制误差，本发明在算法中加入对第k个采样周期电网电压eα(k)、eβ(k)做一个采样周期补偿，对第k个采样周期参考电流i*α(k)、i*β(k)做两个采样周期补偿，如图7中第2步所示。
[0156] 实施例1
[0157] 实施例1的并网控制目标是消除有功2、4、8次波动，按照图8所示的基于模型预测控制的非理想电网条件下NPC三相三电平逆变器并网系统进行了仿真实验，具体参数如表1所示。实施例1的仿真结果如图9所示，图9(a)～(d)分别为三相入网电流、直流侧中点电压、入网瞬时有功功率、入网瞬时无功功率的仿真结果，从图9中可以看出，非理想电网条件下的模型预测控制能够有效地消除入网有功2、4、8次波动，同时抑制直流侧中点电压。
[0158] 表1仿真参数
[0159]基频正序分量幅值：320V 基频负序分量幅值：30V
5倍频负序分量幅值：20V 7倍频正序分量幅值：10V
基频 50Hz
Vdc 800(V)
C1(＝C2) 500e-6(F)
λdc 1
Ts 25e-6(s)
R 1(Ω)
L 4e-3(H)
P*0 10(kW)
Q*0 500(Var)
[0160] 实施例2
[0161] 实施例2的并网控制目标是消除无功2、4、8次波动，具体参数为对表1种做修改P*0＝500W、Q*0＝10kVar，其余参数不变。实施例2的仿真结果如图10所示，图10(a)～(d)分别为三相入网电流、直流侧中点电压、入网瞬时有功功率、入网瞬时无功功率的仿真结果，从图10中可以看出，非理想电网条件下的模型预测控制能够有效地消除入网无功2、4、8次波动，同时抑制直流侧中点电压，此实施例中需要进行大量的无功缓冲，所以直流侧电容电压幅值波动变大，导致直流侧中点电压波动变大，如图10(b)所示。
[0162] 实施例3
[0163] 实施例3的并网控制目标是实现三相入网电流正弦且平衡，具体参数如表1所示。实施例3的仿真结果如图11所示，图11(a)～(d)分别为三相入网电流、直流侧中点电压、入网瞬时有功功率、入网瞬时无功功率的仿真结果，从图11中可以看出，非理想电网条件下的模型预测控制能够有效地实现三相入网电流正弦且平衡，同时抑制直流侧中点电压。
[0164] 综上所述，本发明所采用的模型预测控制方法，对非理想电网条件下NPC三相三电平逆变器进行并网控制，取得了理想的控制效果。
