CN115241912A - A Model Predictive Control Method for Model Mismatch Compensation of Three-phase Grid-connected Inverters - Google Patents
A Model Predictive Control Method for Model Mismatch Compensation of Three-phase Grid-connected Inverters Download PDFInfo
- Publication number
- CN115241912A CN115241912A CN202210870667.1A CN202210870667A CN115241912A CN 115241912 A CN115241912 A CN 115241912A CN 202210870667 A CN202210870667 A CN 202210870667A CN 115241912 A CN115241912 A CN 115241912A
- Authority
- CN
- China
- Prior art keywords
- model
- phase
- inverter
- grid
- voltage
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 22
- 238000013178 mathematical model Methods 0.000 claims abstract description 23
- 239000003990 capacitor Substances 0.000 claims description 10
- 238000005070 sampling Methods 0.000 claims description 9
- 239000013598 vector Substances 0.000 claims description 9
- 230000003068 static effect Effects 0.000 claims description 7
- 230000003071 parasitic effect Effects 0.000 claims description 4
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 claims 1
- 238000000844 transformation Methods 0.000 claims 1
- 230000009466 transformation Effects 0.000 claims 1
- 238000004364 calculation method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000002955 isolation Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000004804 winding Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 239000003302 ferromagnetic material Substances 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000005415 magnetization Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
Images
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
- H02J3/381—Dispersed generators
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02M—APPARATUS FOR CONVERSION BETWEEN AC AND AC, BETWEEN AC AND DC, OR BETWEEN DC AND DC, AND FOR USE WITH MAINS OR SIMILAR POWER SUPPLY SYSTEMS; CONVERSION OF DC OR AC INPUT POWER INTO SURGE OUTPUT POWER; CONTROL OR REGULATION THEREOF
- H02M1/00—Details of apparatus for conversion
- H02M1/12—Arrangements for reducing harmonics from AC input or output
- H02M1/126—Arrangements for reducing harmonics from AC input or output using passive filters
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02M—APPARATUS FOR CONVERSION BETWEEN AC AND AC, BETWEEN AC AND DC, OR BETWEEN DC AND DC, AND FOR USE WITH MAINS OR SIMILAR POWER SUPPLY SYSTEMS; CONVERSION OF DC OR AC INPUT POWER INTO SURGE OUTPUT POWER; CONTROL OR REGULATION THEREOF
- H02M7/00—Conversion of AC power input into DC power output; Conversion of DC power input into AC power output
- H02M7/42—Conversion of DC power input into AC power output without possibility of reversal
- H02M7/44—Conversion of DC power input into AC power output without possibility of reversal by static converters
- H02M7/48—Conversion of DC power input into AC power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode
- H02M7/53—Conversion of DC power input into AC power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal
- H02M7/537—Conversion of DC power input into AC power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal using semiconductor devices only, e.g. single switched pulse inverters
- H02M7/5387—Conversion of DC power input into AC power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal using semiconductor devices only, e.g. single switched pulse inverters in a bridge configuration
- H02M7/53871—Conversion of DC power input into AC power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal using semiconductor devices only, e.g. single switched pulse inverters in a bridge configuration with automatic control of output voltage or current
-
- 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
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Inverter Devices (AREA)
Abstract
Description
技术领域technical field
本发明涉及分布式发电系统预测控制领域,具体涉及一种模型预测控制的三相并网逆变器模型失配补偿方法。The invention relates to the field of predictive control of distributed power generation systems, in particular to a model mismatch compensation method of a three-phase grid-connected inverter for model predictive control.
背景技术Background technique
随着微处理器计算能力的日益提高,模型预测控制(MPC)在电力电子变换领域中得到了广泛的研究,使得分布式发电系统、交流电机控制、不间断电源等都获得了良好的静态和动态性能。MPC可以解耦相互依赖的控制回路,从而提高系统的动态响应,所以基于MPC算法的三相并网逆变器响应速度比采用经典控制器快很多。With the increasing computing power of microprocessors, Model Predictive Control (MPC) has been widely studied in the field of power electronic conversion, which enables distributed generation systems, AC motor control, uninterruptible power supplies, etc. to achieve good static and dynamic performance. MPC can decouple the interdependent control loops, thereby improving the dynamic response of the system, so the response speed of the three-phase grid-connected inverter based on the MPC algorithm is much faster than that of the classical controller.
三相并网逆变器通常采用LCL滤波。相比于L滤波器,LCL滤波器具有三阶低通滤波器特性,在同样谐波标准和较低开关频率下,LCL滤波器总电感值小于L滤波器单电感值。Three-phase grid-connected inverters usually use LCL filtering. Compared with the L filter, the LCL filter has the characteristics of a third-order low-pass filter. Under the same harmonic standard and lower switching frequency, the total inductance value of the LCL filter is smaller than the single inductance value of the L filter.
但是,基于MPC算法的三相并网逆变器如果选用LCL滤波器进行滤波会存在一些问题。在控制算法上,三相并网逆变器的数学模型将变得更加复杂,导致了MPC控制器设计难度的提高。在硬件设计上,需要额外配置采样电路用于获取流入电容的电流值,提高了设备成本。因此,目前基于MPC算法的三相并网逆变器在的大部分应用过程,是忽略了流入电容电流,以传统的电感L滤波为基础建立逆变器动态模型,导致被控对象数学模型存在未建模部分,从而产生控制误差;此外,滤波电感的铁磁材料依据磁化曲线变化从而导致电感实际值在不同工况下产生时变。同时,模型中的电阻不仅是电感的杂散电阻,还包含有逆变器系统中无法有效测量的线阻等其他未知电阻。所以,模型中电感和电阻参数经常与实际对象存在偏差,降低了MPC算法的控制性能。如何补偿上述模型误差对提高MPC性能非常重要。However, if the three-phase grid-connected inverter based on the MPC algorithm chooses the LCL filter for filtering, there will be some problems. In the control algorithm, the mathematical model of the three-phase grid-connected inverter will become more complicated, which leads to the increase of the difficulty of designing the MPC controller. In terms of hardware design, an additional sampling circuit needs to be configured to obtain the current value flowing into the capacitor, which increases the equipment cost. Therefore, most of the current application process of three-phase grid-connected inverter based on MPC algorithm ignores the inflowing capacitor current, and establishes the dynamic model of the inverter based on the traditional inductance L filter, resulting in the existence of the mathematical model of the controlled object. Unmodeled part, resulting in control error; in addition, the ferromagnetic material of the filter inductor changes according to the magnetization curve, which causes the actual value of the inductor to vary in time under different operating conditions. At the same time, the resistance in the model is not only the stray resistance of the inductor, but also includes other unknown resistances such as line resistance that cannot be effectively measured in the inverter system. Therefore, the inductance and resistance parameters in the model often deviate from the actual object, which reduces the control performance of the MPC algorithm. How to compensate for the above model errors is very important to improve the MPC performance.
为了提高MPC的控制性能,部分工程中使用了一些参数辨识方法,以提高数学模型的精确度。然而这种控制策略具有如下缺点:1)不能对模型未建模部分进行补偿2)未对电阻参数进行观测3)增大控制器计算量,从而提高逆变器系统对控制芯片的要求。In order to improve the control performance of MPC, some parameter identification methods are used in some projects to improve the accuracy of the mathematical model. However, this control strategy has the following disadvantages: 1) it cannot compensate the unmodeled part of the model; 2) it does not observe the resistance parameters; 3) it increases the calculation amount of the controller, thereby increasing the requirements of the inverter system on the control chip.
发明内容SUMMARY OF THE INVENTION
有鉴于此,本发明的目的在于提供一种模型预测控制的三相并网逆变器模型失配补偿方法,解决了基于模型预测控制的三相并网逆变器因模型未建模和参数时变导致的预测模型不精确问题,有效提高MPC控制效率及质量,提高三相并网逆变器可靠性。In view of this, the purpose of the present invention is to provide a three-phase grid-connected inverter model mismatch compensation method based on model predictive control, which solves the problem that the three-phase grid-connected inverter based on model predictive control is not modeled and parameterized. The inaccurate prediction model caused by time-varying can effectively improve the efficiency and quality of MPC control, and improve the reliability of the three-phase grid-connected inverter.
为实现上述目的,本发明采用如下技术方案:To achieve the above object, the present invention adopts the following technical solutions:
一种模型预测控制的三相并网逆变器模型失配补偿方法,包括以下步骤:A model-predictive-controlled three-phase grid-connected inverter model mismatch compensation method, comprising the following steps:
步骤S1:建立基于L滤波器的三相并网逆变器离散数学模型,并将采样信号和给定电流值代入基于L滤波器的三相并网逆变器离散数学模型中,计算出预测逆变器输出电压初始值uanref0、ubnref0、ucnref0;Step S1: Establish a discrete mathematical model of the three-phase grid-connected inverter based on the L filter, and substitute the sampled signal and the given current value into the discrete mathematical model of the three-phase grid-connected inverter based on the L filter, and calculate the prediction. Initial values of inverter output voltage u anref0 , u bnref0 , u cnref0 ;
步骤S2:建立基于LCL滤波器的三相并网逆变器实际数学模型,将流入电容的电流所产生的电压定义为模型扰动,将电阻电感偏差值产生的电压与模型扰动之和定义为系统未知扰动;定义采集的电网电压为已知扰动;Step S2: establish the actual mathematical model of the three-phase grid-connected inverter based on the LCL filter, define the voltage generated by the current flowing into the capacitor as the model disturbance, and define the sum of the voltage generated by the resistance and inductance deviation value and the model disturbance as the system. Unknown disturbance; define the collected grid voltage as known disturbance;
步骤S3:定义已知电感量为控制增益,以逆变器输出电压uan、ubn、ucn作为模型失配补偿的输入,逆变侧电流ia、ib、ic作为模型失配补偿的输出,建立模型失配补偿模块,并基于模型失配补偿模块,整定βa1、βa2、βb1、βb2、βc1、βc2后观测出系统未知扰动za2、zb2、zc2;Step S3: define the known inductance as the control gain, take the inverter output voltages u an , ubn , u cn as the input of the model mismatch compensation, and the inverter side currents i a , ib , ic as the model mismatch Compensation output, establish a model mismatch compensation module, and based on the model mismatch compensation module, after tuning β a1 , β a2 , β b1 , β b2 , β c1 , β c2 , the unknown disturbances z a2 , z b2 , z of the system are observed. c2 ;
步骤S4:将系统未知扰动补偿到预测逆变器输出电压初始值uanref0、ubnref0、ucnref0中得到三相预测逆变器输出电压uanref、ubnref、ucnref;Step S4: compensating the unknown disturbance of the system to the initial values u anref0 , u bnref0 , and u cnref0 of the predicted inverter output voltage to obtain the three-phase predicted inverter output voltages u anref , u bnref , and u cnref ;
步骤S5:根据uanref、ubnref、ucnref设计代价函数,将三相预测逆变器输出电压送入模型预测控制器遍历后选出一个电压矢量用于逆变器下一时刻控制。Step S5 : Design a cost function according to u anref , u bnref , and u cnref , send the three-phase predicted inverter output voltage to the model prediction controller for traversal, and select a voltage vector for the inverter to control at the next moment.
进一步的,所述步骤S1,具体为:Further, the step S1 is specifically:
基于LCL滤波器的并网逆变器静止abc坐标系下数学模型为The mathematical model of the grid-connected inverter based on the LCL filter in the static abc coordinate system is:
其中,urx=-Lx2dicx/dt-Rx2icx为模型扰动项;Lx和Rx分别为电容两侧电感以及其寄生电阻之和,即Lx=Lx1+Lx2,Rx=Rx1+Rx2;Among them, u rx =-L x2 di cx /dt-R x2 i cx is the model disturbance term; L x and R x are the sum of the inductance on both sides of the capacitor and its parasitic resistance, that is, L x =L x1 +L x2 , R x =R x1 +R x2 ;
忽略模型扰动项,建立L滤波器简化逆变器动态模型:Ignoring the model disturbance term, an L filter is established to simplify the inverter dynamic model:
因为采样周期Ts足够短,可以使用前向欧拉法,将式(2)离散化后得到逆变器离散数学模型:Because the sampling period T s is short enough, the forward Euler method can be used to discretize equation (2) to obtain the discrete mathematical model of the inverter:
当采样时间小于预设时间,得到:When the sampling time is less than the preset time, we get:
对两相旋转坐标系下的给定电流idref、iqref进行如下变换:The given currents idref and i qref in the two-phase rotating coordinate system are transformed as follows:
将式(5)代入式(3)得到k时刻abc静止坐标系下的预测逆变器输出电压初始值:Substitute equation (5) into equation (3) to obtain the initial value of the predicted inverter output voltage in the abc stationary coordinate system at time k:
进一步的,所述步骤S2具体为:Further, the step S2 is specifically:
设实际系统和预测模型的各相间都存在电感电阻偏差ΔLx和ΔRx,则基于LCL滤波器的三相并网逆变器实际数学模型,具体为:Assuming that there are inductance-resistance deviations ΔL x and ΔR x between the phases of the actual system and the prediction model, the actual mathematical model of the three-phase grid-connected inverter based on the LCL filter is as follows:
令参数不匹配产生的电压误差为内部未知扰动分别为Let the voltage errors generated by parameter mismatch be the internal unknown disturbances, respectively
fa(ia)=(ura+ΔLadia/dt+ΔRaia)f a ( i a )=(u ra +ΔL a di a /dt+ΔR a i a )
fb(ib)=(urb+ΔLbdib/dt+ΔRbib)f b (i b )=( urb +ΔL b di b /dt+ΔR b i b )
fc(ic)=(urc+ΔLcdic/dt+ΔRcic),f c ( ic )=(u rc +ΔL c dic /dt+ΔR c i c ) ,
外部已知扰动分别为三相网侧相电压wa(t)=-ea、wb(t)=-eb和wc(t)=-ec,控制增益分别为ba=1/La、bb=1/Lb和bc=1/Lc The known external disturbances are the three-phase grid side phase voltages w a (t)=-e a , w b (t)=-e b and w c (t)=-e c , respectively, and the control gains are b a =1 /L a , bb =1/L b and bc =1/L c
则:but:
设uan、ubn、ucn为系统输入uas、ubs、ucs,ia、ib、ic为系统输出ya1、yb1、yc1;ia、ib、ic为状态变量xa1、xb1、xc1构建一阶系统如下:Let u an , ubn , u cn be the system inputs u as , ubs , u cs , i a , ib , ic are the system outputs y a1 , y b1 , y c1 ; ia , ib , ic are The state variables x a1 , x b1 , and x c1 construct a first-order system as follows:
令状态变量xa2=fa(ia)、xb2=fb(ib)、xc2=fc(ic),构建新系统为:Let the state variables x a2 = f a ( i a ), x b2 = f b ( i b ), and x c2 = f c ( i c ), the new system is constructed as:
进一步的,所述步骤S3具体为:Further, the step S3 is specifically:
将逆变器输出电压uan、ubn、ucn,逆变侧三相相电流ia、ib、ic和预测逆变器输出电压初始值uanref0、ubnref0、ucnref0代入模型失配补偿模块进行如下计算:Substitute the inverter output voltages u an , ubn , u cn , the three-phase phase currents i a , ib , ic and predicted inverter output voltage initial values u anref0 , u bnref0 , u cnref0 into the model. The compensation module is equipped with the following calculation:
a相模型失配补偿:ea1=za1-ia, A-phase model mismatch compensation: e a1 =z a1 -i a ,
b相模型失配补偿:eb1=zb1-ib, b-phase model mismatch compensation: e b1 =z b1 -i b ,
c相模型失配补偿:ec1=zc1-ic,其中,za1、zb1、zc1是三相电流的估计值,分别是za1、zb1、zc1的导数,分别是za2、zb2、zc2的导数,ea1、eb1、ec1逆变侧三相电流估计值和三相电流采样值的误差值。C-phase model mismatch compensation: e c1 =z c1 -ic , Among them, z a1 , z b1 , z c1 are the estimated values of the three-phase current, are the derivatives of z a1 , z b1 , and z c1 , respectively, They are the derivatives of z a2 , z b2 , and z c2 , respectively, and the error values of the estimated three-phase current at the inverter side of e a1 , e b1 , and e c1 and the sampling value of the three-phase current.
进一步的,所述步骤S4具体为:Further, the step S4 is specifically:
将系统未知扰动补偿到预测逆变器输出电压初始值uanref0、ubnref0、ucnref0中得到三相预测逆变器输出电压,得到k时刻的补偿后的网侧参考电压:Compensate the unknown disturbance of the system to the initial values of the predicted inverter output voltage u anref0 , u bnref0 , u cnref0 to obtain the three-phase predicted inverter output voltage, and obtain the compensated grid-side reference voltage at time k:
进一步的,所述代价函数g(k),具体为:Further, the cost function g(k) is specifically:
g(k)=(uanref(k)-ua(k))2+(ubnref(k)-ub(k))2+(ucnref(k)-uc(k))2 g(k)=(u anref (k)-u a (k)) 2 +(u bnref (k)-u b (k)) 2 +(u cnref (k)-u c (k)) 2
其中,ua(k)、ub(k)、uc(k)是各个电压矢量预估的逆变器输出电压。Among them, u a (k), u b (k), and u c (k) are the inverter output voltages estimated by each voltage vector.
本发明与现有技术相比具有以下有益效果:Compared with the prior art, the present invention has the following beneficial effects:
本发明解决了基于模型预测控制的三相并网逆变器因被控对象的数学模型不精确,有效提高MPC控制效率及质量,提高三相并网逆变器可靠性。The invention solves the problem that the mathematical model of the controlled object of the three-phase grid-connected inverter based on model prediction control is inaccurate, effectively improves the MPC control efficiency and quality, and improves the reliability of the three-phase grid-connected inverter.
附图说明Description of drawings
图1是本发明方法结构框图;Fig. 1 is the method structure block diagram of the present invention;
图2是本发明一实施例中预测逆变器输出电压初始值结构框图;FIG. 2 is a structural block diagram of predicting the initial value of the inverter output voltage in an embodiment of the present invention;
图3是本发明一实施例中模型失配补偿结构框图;3 is a block diagram of a model mismatch compensation structure in an embodiment of the present invention;
图4是本发明一实施例中三相并网逆变器系统硬件结构;4 is a hardware structure of a three-phase grid-connected inverter system in an embodiment of the present invention;
图5是本发明一实施例中坐标系定义。FIG. 5 is a coordinate system definition in an embodiment of the present invention.
具体实施方式Detailed ways
下面结合附图及实施例对本发明做进一步说明。The present invention will be further described below with reference to the accompanying drawings and embodiments.
请参照图1,本实施例中系统由直流电源Vdc供电,任意一种逆变器拓扑接入LCL滤波器后并入电网。其中,La1、Lb1、Lc1和Ra1、Rb1、Rc1分别是逆变器侧电感和逆变器侧电感的寄生电阻;La2、Lb2、Lc2和Ra2、Rb2、Rc2分别是网侧电感和网侧电感的寄生电阻;Ca、Cb、Cc是滤波电容;Ra、Rb、Rc是无源阻尼电阻;滤波电容ia、ib、ic是逆变器输出电流;ica、icb、icc是流入滤波电容的电流;iga、igb、igc是网侧三相相电流;uan、ubn、ucn是逆变器侧输出相电压;ea、eb、ec是网侧相电压。Referring to FIG. 1 , in this embodiment, the system is powered by a DC power supply V dc , and any inverter topology is connected to the LCL filter and then merged into the power grid. Among them, L a1 , L b1 , L c1 and R a1 , R b1 , and R c1 are the inverter-side inductance and the parasitic resistance of the inverter-side inductance, respectively; L a2 , L b2 , L c2 and R a2 , R b2 , R c2 are grid-side inductance and parasitic resistance of grid-side inductance respectively; C a , C b , C c are filter capacitors; R a , R b , R c are passive damping resistors; filter capacitors ia , ib , ic is the inverter output current; ica , icb , icc are the currents flowing into the filter capacitor; i ga , igb , igc are the three-phase currents on the grid side; u an , ubn , u cn are the inverse The inverter side outputs the phase voltage; e a , e b , and e c are the phase voltages on the grid side.
如图2所示为预测逆变器输出电压初始值结构框图。建立基于L滤波器的三相并网逆变器离散数学模型,逆变侧三相相电流、网侧相电压和三相给定电流iaref、ibref、icref作为预测逆变器输出电压初始值的输入,输出静止坐标系下的三相预测逆变器输出电压初始值uanref0、ubnref0、ucnref0。Figure 2 is a block diagram showing the structure of the initial value of the predicted inverter output voltage. A discrete mathematical model of three-phase grid-connected inverter based on L filter is established, and the three-phase phase current on the inverter side, the phase voltage on the grid side and the three-phase given current i aref , i bref , i cref are used as the predicted inverter output voltage Input the initial value, output the initial value u anref0 , u bnref0 , u cnref0 of the output voltage of the three-phase predicted inverter in the static coordinate system.
如图3所示为模型失配补偿结构框图。其中,za1、zb1、zc1分别为状态变量ia、ib、ic的估计值;za2、zb2、zc2为模型失配误差的估计值。βa1、βa2、βb1、βb2、βc1、βc2为可调校正系数。逆变侧三相相电流ia、ib、ic、三相输出电压初始值预测uanref0、ubnref0、ucnref0和逆变器输出电压uan、ubn、ucn作为模型失配的输入,输出三相预测逆变器输出电压uanref、ubnref、ucnref。Figure 3 shows a block diagram of the model mismatch compensation structure. Among them, z a1 , z b1 , and z c1 are the estimated values of the state variables i a , ib , and ic , respectively; z a2 , z b2 , and z c2 are the estimated values of the model mismatch error. β a1 , β a2 , β b1 , β b2 , β c1 , and β c2 are adjustable correction coefficients. The three-phase phase currents i a , ib , ic , the initial value predictions of the three-phase output voltages on the inverter side u anref0 , u bnref0 , u cnref0 and the inverter output voltages u an , ubn , and u cn are used as the model mismatch. Input and output three-phase predicted inverter output voltages u anref , u bnref , u cnref .
参考图4,本实施例中整个控制系统包括:直流电源、三相并网逆变器、滤波电路、继电器开关、直流母线采样、三相电网电压采集电路、逆变器侧三相电流采集电路、控制器、隔离驱动、人机交互部分等。Referring to FIG. 4 , the entire control system in this embodiment includes: a DC power supply, a three-phase grid-connected inverter, a filter circuit, a relay switch, a DC bus sampling, a three-phase grid voltage acquisition circuit, and a three-phase current acquisition circuit on the inverter side , controller, isolation driver, human-computer interaction part, etc.
其中三相逆变器采用合适的直流电源提供。逆变器中功率管采用带有二极管并联的IGBT或MOSFET,控制器采用DSP或单片机。逆变器侧三相电流采集电路选用霍尔电流传感器与运算放大器相结合方式构成,也可以采用绕组串功率电阻后接差分运算放大器相结合方式构成。采用霍尔传感器方案可以有效实现控制回路与主回路的电气隔离,采用绕组串功率电阻方案可以降低系统成本。电网电压采集电路采用精密电压传感器与运算放大器相结合方式构成,也可以采用并联电阻分压后接入由运算放大器构成的电压跟随器相结合方式构成。电流检测、电压采样电路、输出的弱电信号送到控制器的A/D转换模块。当进入并网模式或者是进入孤岛模式的时候,控制器的GPIO模块输出继电器动作信号,经过电平转换芯片后再由隔离驱动对继电器进行动作。The three-phase inverter is provided by a suitable DC power supply. In the inverter, the power tube adopts IGBT or MOSFET with diodes in parallel, and the controller adopts DSP or single-chip microcomputer. The three-phase current acquisition circuit on the inverter side is composed of a combination of a Hall current sensor and an operational amplifier, or a combination of a winding string power resistor followed by a differential operational amplifier. The use of the Hall sensor scheme can effectively realize the electrical isolation of the control loop and the main loop, and the use of the winding string power resistance scheme can reduce the system cost. The grid voltage acquisition circuit is composed of a combination of a precision voltage sensor and an operational amplifier, or a combination of a voltage follower composed of an operational amplifier after dividing the voltage by a parallel resistor. The current detection, voltage sampling circuit, and the output weak current signal are sent to the A/D conversion module of the controller. When entering the grid-connected mode or entering the island mode, the GPIO module of the controller outputs the relay action signal, and then the relay is actuated by the isolation driver after passing through the level conversion chip.
坐标系定义如图5所示,abc轴为静止三相坐标系,αβ轴为静止两相坐标系,dq轴为旋转的两相坐标系。The definition of the coordinate system is shown in Figure 5. The abc axis is a static three-phase coordinate system, the αβ axis is a static two-phase coordinate system, and the dq axis is a rotating two-phase coordinate system.
根据图5中的坐标系定义,对网侧相电压ea、eb、ec进行Clark变换:According to the definition of the coordinate system in Fig. 5, Clark transform is performed on the grid-side phase voltages e a , e b , and e c :
根据egα、egβ可计算相位角θg:The phase angle θ g can be calculated according to e gα and e gβ :
根据图5中的坐标系定义,对两相旋转坐标系下的给定电流idref、iqref进行如下变换:According to the definition of the coordinate system in Fig. 5, the given currents idref and i qref in the two-phase rotating coordinate system are transformed as follows:
本实施例中,提出一种模型预测控制的三相并网逆变器模型失配补偿方法,包括以下步骤:In this embodiment, a model mismatch compensation method for a three-phase grid-connected inverter based on model predictive control is proposed, which includes the following steps:
一种模型预测控制的三相并网逆变器模型失配补偿方法,其特征在于,包括以下步骤:A model-predictive-controlled three-phase grid-connected inverter model mismatch compensation method, characterized by comprising the following steps:
步骤S1:建立基于L滤波器的三相并网逆变器离散数学模型,并将采样信号和给定电流值代入基于L滤波器的三相并网逆变器离散数学模型中,计算出预测逆变器输出电压初始值uanref0、ubnref0、ucnref0;Step S1: Establish a discrete mathematical model of the three-phase grid-connected inverter based on the L filter, and substitute the sampled signal and the given current value into the discrete mathematical model of the three-phase grid-connected inverter based on the L filter, and calculate the prediction. Initial values of inverter output voltage u anref0 , u bnref0 , u cnref0 ;
步骤S2:建立基于LCL滤波器的三相并网逆变器实际数学模型,将流入电容的电流所产生的电压定义为模型扰动,将电阻电感偏差值产生的电压与模型扰动之和定义为系统未知扰动;定义采集的电网电压为已知扰动;Step S2: establish the actual mathematical model of the three-phase grid-connected inverter based on the LCL filter, define the voltage generated by the current flowing into the capacitor as the model disturbance, and define the sum of the voltage generated by the resistance and inductance deviation value and the model disturbance as the system. Unknown disturbance; define the collected grid voltage as known disturbance;
步骤S3:定义已知电感量为控制增益,以逆变器输出电压uan、ubn、ucn作为模型失配补偿的输入,逆变侧电流ia、ib、ic作为模型失配补偿的输出,建立模型失配补偿模块,并基于模型失配补偿模块,整定βa1、βa2、βb1、βb2、βc1、βc2后观测出系统未知扰动za2、zb2、zc2;Step S3: define the known inductance as the control gain, take the inverter output voltages u an , ubn , u cn as the input of the model mismatch compensation, and the inverter side currents i a , ib , ic as the model mismatch Compensation output, establish a model mismatch compensation module, and based on the model mismatch compensation module, after tuning β a1 , β a2 , β b1 , β b2 , β c1 , β c2 , the unknown disturbances z a2 , z b2 , z of the system are observed. c2 ;
步骤S4:将系统未知扰动补偿到预测逆变器输出电压初始值uanref0、ubnref0、ucnref0中得到三相预测逆变器输出电压uanref、ubnref、ucnref;Step S4: compensating the unknown disturbance of the system to the initial values u anref0 , u bnref0 , and u cnref0 of the predicted inverter output voltage to obtain the three-phase predicted inverter output voltages u anref , u bnref , and u cnref ;
步骤S5:根据uanref、ubnref、ucnref设计代价函数,将三相预测逆变器输出电压送入模型预测控制器遍历后选出一个电压矢量用于逆变器下一时刻控制。Step S5 : Design a cost function according to u anref , u bnref , and u cnref , send the three-phase predicted inverter output voltage to the model prediction controller for traversal, and select a voltage vector for the inverter to control at the next moment.
在本实施例中,直流电压经过逆变器系统,经过LCL滤波器消除电流高次谐波后并入电网;将采集网侧相电压ea、eb、ec变换两相静止坐标系得到egα、egβ,其中计算公式如下:In this embodiment, the DC voltage passes through the inverter system, and then passes through the LCL filter to eliminate the high-order harmonics of the current and then merges into the power grid; the phase voltages e a , e b , and e c on the collected grid side are transformed into the two-phase static coordinate system to obtain e gα , e gβ , where the calculation formula is as follows:
计算网侧相电压相位角θg,其中计算公式如下:Calculate the phase angle θ g of the phase voltage on the grid side, and the calculation formula is as follows:
将两相旋转坐标系下的给定电流idref、iqref变换到三相静止坐标系下得到iaref、ibref、icref,其中计算公式如下:Transform the given currents i dref and i qref in the two-phase rotating coordinate system to the three-phase stationary coordinate system to obtain i aref , i bref , i cref , and the calculation formulas are as follows:
将给定电流iaref、ibref、icref,逆变侧三相相电流ia、ib、ic和网侧相电压ea、eb、ec代入基于L滤波器的三相并网逆变器离散数学模型得到预测逆变器输出电压初始值uanref0、ubnref0、ucnref0,Substitute the given currents i aref , i bref , i cref , the three-phase currents i a , i b , i c on the inverter side and the phase voltages e a , eb , and ec on the grid side into the three-phase parallel connection based on the L filter. The discrete mathematical model of the grid inverter obtains the initial values of the predicted inverter output voltage u anref0 , u bnref0 , u cnref0 ,
在本实施例中,步骤S2具体为:In this embodiment, step S2 is specifically:
设实际系统和预测模型的各相间都存在电感电阻偏差ΔLx和ΔRx,则基于LCL滤波器的三相并网逆变器实际数学模型,具体为:Assuming that there are inductance-resistance deviations ΔL x and ΔR x between the phases of the actual system and the prediction model, the actual mathematical model of the three-phase grid-connected inverter based on the LCL filter is as follows:
令参数不匹配产生的电压误差为内部未知扰动分别为Let the voltage errors generated by parameter mismatch be the internal unknown disturbances, respectively
fa(ia)=(ura+ΔLadia/dt+ΔRaia)f a ( i a )=(u ra +ΔL a di a /dt+ΔR a i a )
fb(ib)=(urb+ΔLbdib/dt+ΔRbib)f b (i b )=( urb +ΔL b di b /dt+ΔR b i b )
fc(ic)=(urc+ΔLcdic/dt+ΔRcic),f c ( ic )=(u rc +ΔL c dic /dt+ΔR c i c ) ,
外部已知扰动分别为三相网侧相电压wa(t)=-ea、wb(t)=-eb和wc(t)=-ec,控制增益分别为ba=1/La、bb=1/Lb和bc=1/LcThe known external disturbances are the three-phase grid side phase voltages w a (t)=-e a , w b (t)=-e b and w c (t)=-e c , respectively, and the control gains are b a =1 /L a , bb =1/L b and bc =1/Lc
则:but:
设uan、ubn、ucn为系统输入uas、ubs、ucs,ia、ib、ic为系统输出ya1、yb1、yc1;ia、ib、ic为状态变量xa1、xb1、xc1构建一阶系统如下:Let u an , ubn , u cn be system inputs u as , ubs , u cs , i a , ib , ic are system outputs y a1 , y b1 , y c1 ; ia , ib , ic are state variables x a1 , x b1 , x c1 construct a first-order system as follows:
令状态变量xa2=fa(ia)、xb2=fb(ib)、xc2=fc(ic),构建新系统为:Let the state variables x a2 = f a ( i a ), x b2 = f b ( i b ), and x c2 = f c ( i c ), the new system is constructed as:
模型失配补偿模块能够根据系统输入量uan、ubn、ucn与输出量ia、ib、ib对状态变量xa1、xb1、xc1和xa2、xb2、xc2进行估计,进而实现动态扰动前馈补偿,针对式系统建立如下等式:The model mismatch compensation module can perform the state variables x a1 , x b1 , x c1 and x a2 , x b2 , x c2 according to the system input quantities u an , ub bn , u cn and the output quantities i a , ib , and ib . estimate, and then realize the dynamic disturbance feedforward compensation, and establish the following equation for the system:
以a相为例,从上式可以看出,za2是xa2经过二阶低通滤波器滤波后的值,za1是xa1经滤波再进行相位补偿的值。βa2的取值直接对参数失配补偿的估计能力造成影响,随着βa2的增大,系统抗干扰增强,即发生扰动后系统恢复到稳态的时间短。但βa2进一步增大将会使系统出现振荡现象且对噪声敏感问题,此时又要增大βa1来抑制这一问题,但是βa1过大又将造成系统发散。综上,实际参数整定过程中需要根据具体情况对观测器极点进行合适的选择。Taking phase a as an example, it can be seen from the above formula that z a2 is the value of x a2 after being filtered by a second-order low-pass filter, and z a1 is the value of x a1 that is filtered and then phase compensated. The value of β a2 directly affects the estimation ability of parameter mismatch compensation. With the increase of β a2 , the anti-interference of the system is enhanced, that is, the time for the system to recover to a steady state after disturbance occurs is short. However, a further increase of β a2 will make the system oscillate and be sensitive to noise. At this time, β a1 should be increased to suppress this problem, but if β a1 is too large, the system will diverge. To sum up, in the actual parameter tuning process, it is necessary to select the observer pole appropriately according to the specific situation.
同样的分析方法可说明b相和c相电流环,此处不再重复分析。The same analysis method can illustrate the b-phase and c-phase current loops, and the analysis will not be repeated here.
最后,加入观测结果,得到k时刻的补偿后的网侧参考电压:Finally, add the observation results to get the compensated grid-side reference voltage at time k:
优选的,在本实施例中,仅以电流跟踪为唯一目标设计代价函数。并网逆变器中电流控制通常选取电流、电压或功率作为主要约束变量,以电压约束变量为例建立代价函数g(k):Preferably, in this embodiment, the cost function is designed only with current tracking as the only target. In the current control of grid-connected inverters, current, voltage or power are usually selected as the main constraint variables, and the cost function g(k) is established by taking the voltage constraint variables as an example:
g(k)=(uanref(k)-ua(k))2+(ubnref(k)-ub(k))2+(ucnref(k)-uc(k))2 g(k)=(u anref (k)-u a (k)) 2 +(u bnref (k)-u b (k)) 2 +(u cnref (k)-u c (k)) 2
其中,ua(k)、ub(k)、uc(k)是各个电压矢量预估的逆变器输出电压,因为三相电流跟踪为同一量纲所以各平方相的系数均为1。代价函数对所有的矢量进行遍历后,得到使得g(k)最小的电压矢量,该电压矢量用于下一时刻三相并网逆变器的控制。Among them, u a (k), u b (k), and u c (k) are the inverter output voltages estimated by each voltage vector, because the three-phase current tracking is the same dimension, the coefficient of each square phase is 1 . After the cost function traverses all the vectors, the voltage vector that minimizes g(k) is obtained, and the voltage vector is used for the control of the three-phase grid-connected inverter at the next moment.
以上所述仅为本发明的较佳实施例,凡依本发明申请专利范围所做的均等变化与修饰,皆应属本发明的涵盖范围。The above descriptions are only preferred embodiments of the present invention, and all equivalent changes and modifications made according to the scope of the patent application of the present invention shall fall within the scope of the present invention.
Claims (6)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210870667.1A CN115241912B (en) | 2022-07-22 | 2022-07-22 | Model prediction control three-phase grid-connected inverter model mismatch compensation method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210870667.1A CN115241912B (en) | 2022-07-22 | 2022-07-22 | Model prediction control three-phase grid-connected inverter model mismatch compensation method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115241912A true CN115241912A (en) | 2022-10-25 |
CN115241912B CN115241912B (en) | 2024-07-30 |
Family
ID=83676073
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210870667.1A Active CN115241912B (en) | 2022-07-22 | 2022-07-22 | Model prediction control three-phase grid-connected inverter model mismatch compensation method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115241912B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117240126A (en) * | 2023-11-15 | 2023-12-15 | 通达电磁能股份有限公司 | Limited set model predictive control method, system, terminal and storage medium |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109004669A (en) * | 2018-07-26 | 2018-12-14 | 华中科技大学 | Three-phase grid-connected inverter based on interference observer compensation improves passive control method |
CN111221253A (en) * | 2020-03-11 | 2020-06-02 | 国网江苏省电力有限公司泗洪县供电分公司 | A robust model predictive control method for three-phase grid-connected inverters |
CN112311290A (en) * | 2020-10-26 | 2021-02-02 | 江苏大学 | Robust prediction permanent magnet synchronous hub motor sensorless controller |
CN113472257A (en) * | 2021-06-07 | 2021-10-01 | 浙江工业大学 | Disturbance compensation-based finite set model prediction current control method for permanent magnet direct-drive wind power generation system |
CN114243935A (en) * | 2021-09-13 | 2022-03-25 | 西安交通大学 | A wireless power transmission control method and system based on P&O method and model prediction |
-
2022
- 2022-07-22 CN CN202210870667.1A patent/CN115241912B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109004669A (en) * | 2018-07-26 | 2018-12-14 | 华中科技大学 | Three-phase grid-connected inverter based on interference observer compensation improves passive control method |
CN111221253A (en) * | 2020-03-11 | 2020-06-02 | 国网江苏省电力有限公司泗洪县供电分公司 | A robust model predictive control method for three-phase grid-connected inverters |
CN112311290A (en) * | 2020-10-26 | 2021-02-02 | 江苏大学 | Robust prediction permanent magnet synchronous hub motor sensorless controller |
CN113472257A (en) * | 2021-06-07 | 2021-10-01 | 浙江工业大学 | Disturbance compensation-based finite set model prediction current control method for permanent magnet direct-drive wind power generation system |
CN114243935A (en) * | 2021-09-13 | 2022-03-25 | 西安交通大学 | A wireless power transmission control method and system based on P&O method and model prediction |
Non-Patent Citations (1)
Title |
---|
王楠楠;刘泽浩;李洋;: "基于滤波电感辨识及延时补偿的并网逆变器模型预测控制", 曲阜师范大学学报(自然科学版), no. 04, 15 October 2017 (2017-10-15), pages 59 - 63 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117240126A (en) * | 2023-11-15 | 2023-12-15 | 通达电磁能股份有限公司 | Limited set model predictive control method, system, terminal and storage medium |
CN117240126B (en) * | 2023-11-15 | 2024-01-23 | 通达电磁能股份有限公司 | Limited set model predictive control method, system, terminal and storage medium |
Also Published As
Publication number | Publication date |
---|---|
CN115241912B (en) | 2024-07-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102916599B (en) | Model predication control method of three-phase PWM (pulse width modulation) rectifier under unbalanced voltage | |
Ahmed et al. | A modified stationary reference frame-based predictive current control with zero steady-state error for LCL coupled inverter-based distributed generation systems | |
CN102916600B (en) | Self-correcting prediction control method of model of three-phase voltage type PWM (Pulse-Width Modulation) rectifier | |
Ketzer et al. | Sensorless control technique for PWM rectifiers with voltage disturbance rejection and adaptive power factor | |
CN110460250B (en) | Direct power control method for three-phase PWM rectifier | |
CN102655327B (en) | Control method for sliding mode converter control structure of active power filter containing parameter perturbation | |
CN114884125B (en) | High-stability control method of LCL type grid-connected inversion system under weak current network | |
CN112165271B (en) | A grid-connected converter system and its model predictive control method | |
CN103066866B (en) | Active front end rectifier filtering delay compensation method based on model prediction controlling | |
CN103780107B (en) | Current control method for three-phase voltage source type PWM rectifier | |
CN111130366A (en) | Prediction control method of non-grid voltage sensor of three-phase VIENNA rectifier | |
WO2022105068A1 (en) | Control method and device for grid-connected inverter | |
Liao et al. | Robust dichotomy solution-based model predictive control for the grid-connected inverters with disturbance observer | |
US11251720B1 (en) | Method and apparatus for controlling grid-tie inverter | |
CN113629984B (en) | Three-phase LCL type SAPF parameter design method based on double-loop current control strategy | |
CN106786639A (en) | A kind of Active Power Filter-APF improves wideband self-adapting resonance control method | |
CN112000018B (en) | A robust fault-tolerant control module, method and system based on residual generator | |
CN109802433A (en) | A kind of gird-connected inverter oscillation of power inhibition system and method | |
CN112310991B (en) | Energy storage converter improved finite set model prediction control method based on disturbance observation | |
Yu et al. | An improved grid impedance estimator for grid-forming converters in consideration of controller dynamics | |
CN115241912A (en) | A Model Predictive Control Method for Model Mismatch Compensation of Three-phase Grid-connected Inverters | |
CN112994110B (en) | LC filtering type grid-connected inverter parameter-free prediction capacitor voltage control method | |
Li et al. | Current estimation and optimal control in multiphase DC–DC converters with single current sensor | |
CN115276439A (en) | LCL type grid-connected inverter resonance suppression method adapting to impedance change of weak grid | |
CN115296331A (en) | Mapping adaptive backstep sliding mode control method for LCL photovoltaic grid-connected inverter |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |