Journal of Modern Power Systems and Clean Energy

ISSN 2196-5625 CN 32-1884/TK

网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

Model Predictive Direct Power Control of Grid-connected Converters Considering Unbalanced Filter Inductance and Grid Conditions  PDF

  • Weichen Yang
  • Shihong Miao
  • Zhiwei Liu
  • Ji Han
  • Yulong Xiong
  • Qingyu Tu
State Key Laboratory of Advanced Electromagnetic Engineering and Technology, Hubei Electric Power Security and High Efficiency Key Laboratory, School of Electrical and Electronic Engineering, Huazhong University of Science and Technology, Wuhan 430074, China

Updated:2021-11-23

DOI:10.35833/MPCE.2021.000355

  • Full Text
  • Figs & Tabs
  • References
  • Authors
  • About
CITE
OUTLINE

Abstract

The grid-connected converter (GCC) is widely used as the interface between various distributed generations and the utility grid. To achieve precise power control for GCC, this paper presents a model predictive direct power control (MPDPC) with consideration of the unbalanced filter inductance and grid conditions. First, the characteristics of GCC with unbalanced filter inductance are analyzed and a modified voltage control function is derived. On this basis, to compensate for the power oscillation caused by unbalanced filter inductance, a novel power compensation method is proposed for MPDPC to eliminate the DC-side current ripple while maintaining sinusoidal grid current. Besides, to improve the control robustness against mismatched filter inductance, a filter inductance identification scheme is proposed. Through this scheme, the estimated value of filter inductance is updated in each control period and applied in the proposed MPDPC. Finally, simulation results in PSCAD/EMTDC confirm the validity of the proposed MPDPC and the filter inductance identification scheme.

I. Introduction

WITH the rising number of distributed renewable generations and electric vehicles, grid-connected converters (GCCs) have been widely used as the interface between various distributed generation and the utility grid, as they can realize the conversion of DC electric power to AC electric power [

1], [2].

To achieve precise power control of GCCs, various control strategies such as voltage-oriented control (VOC) [

3], direct power control (DPC) [4], model predictive control (MPC) [5] have been proposed. In VOC, the three-phase components are decoupled in the dq-frame by using a phase-locked loop (PLL), and then regulated by proportional-integral (PI) controllers [6]. Those PI controllers give good dynamic/steady-state performance. However, their linear nature is the main drawback [7]. Besides, the performance and stability of VOC scheme are highly dependent on the performance of the PLL [8] and the tuning of the PI gains [9]. In DPC, the instantaneous active power and reactive power are set as control variables. It calculates the desired voltage, and directly selects the best-fitted switching vector for the power regulation [10], [11]. However, the switching table based DPC will cause large power ripples and variable switching frequencies, and a high sampling frequency is required to ensure the control performance. Recently, space vector modulation (SVM) based DPC has been introduced for its simple algorithm, constant switching frequency [12], and lower requirement for sampling frequency. With the development of microprocessors, MPC has become a promising way to obtain a faster dynamic response and realize complex control objectives [13]. Generally, in each control period, the MPC uses the prediction information to select the best voltage vector to minimize a cost function [14], [15]. By combining the advantages of DPC and MPC, accurate power control and various control objectives can be realized based on model predictive direct power control (MPDPC) [16], [17].

In practical applications, weak grid conditions are often encountered in the distribution network and microgrid, and GCCs face unbalanced grid voltages due to the integration of single-phase generators and unbalanced loads [

18]. The unbalanced grid conditions will produce negative-sequence current, and further cause the output power oscillation and DC-side current ripple of GCC. Many DPC methods have been proposed to control GCC under unbalanced grid conditions. A new definition of instantaneous reactive power is proposed in [19] for DPC strategies, which can eliminate active power oscillation while maintaining sinusoidal grid currents. In [20], a power compensation block is added to the original DPC scheme to eliminate the negative-sequence current. The principle of power compensation can be extended to achieve other control objectives such as eliminating grid-side active and reactive power oscillations [21]. However, [22] indicates that the filter inductance will cause the variation of oscillation power. The active power oscillation at the valve side will produce DC-side current ripple and cannot be eliminated with the conventional DPC method. To address this problem, a power compensation calculation method with consideration of the DC-side voltage ripple reduction is introduced in [23] and [24]. MPC methods have also been proposed to achieve proper control performance under unbalanced grid conditions. In [25], an MPDPC method is proposed to obtain sinusoidal grid currents and eliminate the grid-side power oscillation. To achieve robust control without measured grid voltage, a sliding-mode-observer-based voltage-sensorless MPC method is proposed in [2] and [26]. To obtain a better performance comparing with table-based MPC [27], a multi-vector MPC method is proposed in [28] and a vector optimization method is proposed in [29]. The above control methods can achieve their objectives with accurate knowledge of the filter inductance. However, the actual filter inductance may vary due to the temperature and conducting current [30]. The control performance will be significantly degraded with the mismatched value of inductance applied in the control scheme.

To improve the control performance against the variable parameters, various control schemes have been proposed by combining parameter estimation methods with DPC or MPC such as least-squares method [

31]-[33], extended Kalman filter based method [34], gradient correction method [30], observer-based estimation method [35]-[37], and model reference adaptive method [38]. However, most of existing parameter identification methods used in converters can only work under ideal grid conditions and cannot realize the identification of unbalanced filter inductance. Another method to solve parameter mismatch is the model-free predictive control. Generally, the model-free predictive control is based on the analysis of current increment, and a state observer or neural network is designed to estimate the relationship between the reference voltage and grid current [39]-[41]. The unknown filter inductance and unmodeled dynamic characteristics are regarded as a black box, and the reference voltage can be directly obtained by using the observation results. However, due to the lack of model knowledge, the model-free predictive control has limited robustness under complex conditions.

This paper aims to propose an MPDPC with unbalanced filter inductance identification to achieve precise power control and active power oscillation elimination for GCCs under unbalanced grid conditions. The main contributions of this paper can be drawn as follows.

1) The performance characteristics of GCC with unbalanced filter inductance are analyzed, and a novel power compensation method is proposed to eliminate the valve-side active power oscillation and DC-side current ripple.

2) An unbalanced filter inductance identification scheme is proposed, which has not been fully investigated in the existing literature, especially under unbalanced grid conditions.

3) By implementing power compensation and inductance identification, the proposed MPDPC has better robustness against mismatched inductance and unbalanced grid conditions.

The remainder of this paper is organized as follows. Section II introduces the mathematical model of GCCs under unbalanced grid conditions. Then, the detailed design of the proposed MPDPC is elaborated in Section III. In Section IV, the identification scheme for unbalanced filter inductance is designed. In Section V, the validity of the proposed MPDPC is verified by simulations in PSCAD/EMTDC. Finally, conclusions are drawn in Section VI.

II. Model of GCCs Under Unbalanced Grid Conditions

Figure 1 shows the topology of a GCC, which is connected to the DC supply through the line resistance and to the AC grid through the filter inductance.

Fig. 1 Topology of a GCC.

The mathematical model of the GCC in the stationary reference frame can be expressed as:

didtLαβ0=u-e (1)

where u=[uα,uβ,u0]T, e=[eα,eβ,e0]T, and i=[iα,iβ,i0]T are the converter output voltage, grid voltage, and grid current in the stationary reference frame αβ0, respectively; and Lαβ0 is the inductance matrix in the stationary reference frame, which can be depicted by the three-phase filter inductance and Clarke transformation:

Lαβ0=C3/2La000Lb000LcC3/2-1=164La+Lb+Lc-3Lb+3Lc4La-2Lb-2Lc-3Lb+3Lc3Lb+3Lc23Lb-23Lc2La-Lb-Lc3Lb-3Lc2La+2Lb+2Lc (2)

where C3/2 is the Clarke transformation matrix; and La, Lb, and Lc are the three-phase filter inductances.

Considering that the converter will not generate zero-sequence current under normal conditions, the inductance matrix can be simplified as:

Lαβ=164La+Lb+Lc-3Lb+3Lc-3Lb+3Lc3Lb+3Lc=L11L12L12L22 (3)

where Lαβ is the simplified inductance matrix consisting of L11, L12, and L12.

Then, the control function of GCC can be written as:

uαuβ=eαeβ+Lm00Lmdiα/dtdiβ/dt+L11-LmL12L12L22-Lmdiα/dtdiβ/dtΔuL (4)

where Lm=(La+Lb+Lc)/3 is the average value of three-phase filter inductances, which represents the balanced part of filter inductance; and ΔuL is a disturbance voltage caused by the unbalanced part of filter inductance.

It can be observed that when the filter inductance is symmetrical, Lαβ will be a scalar matrix. However, with the unbalanced filter inductance, ΔuL will be generated by the coupling of iα and iβ. Considering the unbalanced grid conditions, the grid current can be written as:

iα=Ipcos(ωt+θp)+Incos(ωt+θn)iβ=Ipsin(ωt+θp)-Insin(ωt+θn) (5)

where ω is the fundamental angular frequency; Ip and θp are the amplitude and phase angle of the positive-sequence current, respectively; and In and θn are the amplitude and phase angle of the negative-sequence current, respectively.

Combining (4) and (5), ΔuL=[ΔuLα,ΔuLβ]T can be expressed as:

ΔuLα=-ω(L11-Lm)Ipsin(ωt+θp)+ωL12Ipcos(ωt+θp)ΔuLα1-ω(L11-Lm)Insin(ωt+θn)-ωL12Incos(ωt+θn)ΔuLα2ΔuLβ=ω(L22-Lm)Ipcos(ωt+θp)-ωL12Ipsin(ωt+θp)ΔuLβ1-ω(L22-Lm)Incos(ωt+θn)-ωL12Insin(ωt+θn)ΔuLβ2 (6)

where [ΔuLα1,ΔuLβ1]T=ΔuL1 and [ΔuLα2,ΔuLβ2]T=ΔuL2 are the disturbance voltages generated by the positive-sequence current and negative-sequence current, respectively.

Considering L11-Lm=-(L22-Lm), it can be observed that ΔuL1 is a negative-sequence voltage, while ΔuL2 is a positive-sequence voltage. According to instantaneous power theory and the extended reactive power proposed in [

19], the active power and reactive power at the grid side can be expressed as:

P=1.5(eαiα+eβiβ)Qnov=1.5(eα'iα+eβ'iβ) (7)

where [eα',eβ'] lags [eα,eβ] by 90 electrical degrees.

By using Qnov instead of traditional reactive power Q in the control scheme, active power P can be kept constant while maintaining sinusoidal grid currents when the grid voltage is unbalanced [

19].

Under unbalanced grid voltage, the active and reactive power losses ΔP and ΔQ on the filter inductance can be written as:

ΔP=1.5[-2ωLmIpInsin(2ωt+θp+θn)-0.5ω(L11-L22)(Ip2sin(2ωt+2θp)+In2sin(2ωt+2θn))+ωL12(Ip2cos(2ωt+2θp)-In2cos(2ωt+2θn))] (8)
ΔQ=1.5[2ωLmIpIncos(2ωt+θp+θn)+0.5ω(L11-L22)(Ip2cos(2ωt+2θp)+In2cos(2ωt+2θn))+ωL12(Ip2sin(2ωt+2θp)-In2sin(2ωt+2θn))+ωLm(Ip2+In2)+ω(L11-L22)IpIncos(θp-θn)+2ωL12IpInsin(θp-θn)] (9)

According to (4), (8), and (9), it can be observed that the unbalanced filter inductance will degrade the performance of current control and cause the double fundamental frequency power oscillation.

III. MPDPC of GCC Under Unbalanced Grid Conditions

A. Principle of MPDPC Without Output Power Oscillations

The unbalanced grid voltage and filter inductance will generate a negative-sequence current and further cause the double-frequency power oscillation. A typical control objective is to eliminate the active power oscillation at the grid side and control the active power and reactive power to follow their references. By using Qnov instead of traditional reactive power Q in the control scheme, the reactive power oscillation can also be eliminated while maintaining sinusoidal grid current [

19].

Since the MPC works in discrete time, and the digital control system needs a certain time to calculate the references of i and u, a delay of one sampling instant is taken into consideration by making predictions at instant k+2. And the active power and reactive power at the grid side are controlled to reach their references at instant k+2. Hence, the control objective is to minimize the power tracking error which is expressed as:

Pk+2=1.5(eαk+2iαk+2+eβk+2iβk+2)PrefQk+2=1.5(eα'k+2iαk+2+eβ'k+2iβk+2)QrefJmin=(Pk+2-Pref)2+(Qk+2-Qref)2 (10)

where the superscript k+2 represents the variable at instant k+2; Pref and Qref are the references of active power and reactive power, respectively; and Jmin is the cost function.

To obtain the current reference with higher accuracy, the predictive values of grid voltage need to be calculated. The unbalanced grid voltage can be expressed as a complex vector:

e=ep+en=Epejωt+Ene-jωt (11)

where ep and en are the positive-sequence voltage and negative-sequence voltage, respectively; and Ep and En are the amplitudes of the positive-sequence voltage and negative-sequence voltage, respectively.

In the steady state, the grid voltage at instant k+n can be predicted as:

eprek+n=epkejωnTs+enke-jωnTs (12)

where eprek+n is the predicted value of grid voltage at instant k+n; and Ts is the sampling and control period.

According to (10), the optimal current references at instant k+2 can be directly calculated as:

iα,refk+2=23Prefeβ,pre'k+2-Qrefeβ,prek+2eα,prek+2eβ,pre'k+2-eα,pre'k+2eβ,prek+2iβ,refk+2=23-Prefeα,pre'k+2+Qrefeα,prek+2eα,prek+2eβ,pre'k+2-eα,pre'k+2eβ,prek+2 (13)

where [iα,refk+2 ,iβ,refk+2]=irefk+2 is the current reference at instant k+2; eα,prek+2 and eβ,prek+2 are the predicted values of eα and eβ at instant k+2, respectively; and eα,pre'k+2 and eβ,pre'k+2 are the predicted values of eα' and eβ' at instant k+2, respectively.

The grid current at instant k+2 is determined by the grid voltage and converter output voltage at instant k+1. To make the grid current follow the references, the control function of u can be derived based on the first-order Euler discretization method:

uαk+1uβk+1=eα,prek+1eβ,prek+1+1TsL11L12L12L22iα,refk+2-iα,prek+1iβ,refk+2-iβ,prek+1 (14)

where iα,prek+1 and iβ,prek+1 are the predictive values of grid current at instant k+1, which can be calculated by:

iα,prek+1iβ,prek+1=TsL11L12L12L22-1uαk-eαkuβk-eβk+iαkiβk (15)

The reference uk+1 will be adopted on SVM at the next control period to generate the desired converter output voltage and achieve the control objectives. By using SVM, the proposed MPDPC will have a lower requirement for the sampling frequency and computational burden.

B. Principle of MPDPC Without DC-side Current Ripple

The MPDPC shown in Section III-A is used to eliminate the power oscillation at the grid side, thereby improving the control performance and decreasing the current ripple at the DC side. However, according to (8) and (9), both the balanced part and unbalanced part of the filter inductance will cause power oscillation, which makes the active power at the valve side contain a double fundamental frequency oscillation, and further generates a ripple current at the DC side. In this subsection, the control objective is to eliminate the ripple current at the DC side and control the average active power and reactive power at the grid side to follow their references.

According to (10) and (14), the active power and reactive power at instant k+2 at the grid side and the active power and reactive power at instant k+1 at the valve side are controlled by the output voltage at instant k+1. The power oscillation caused by the filter inductance can be compensated by adding double fundamental frequency components on Pref and Qref, thereby eliminating the power oscillation at the valve side. Considering the above factors, the control objective can be expressed as:

Pvk+1=1.5(uαk+1iαk+1+uβk+1iβk+1)PrefPk+2=1.5(eαk+2iαk+2+eβk+2iβk+2)Pref+P2ωk+2Qk+2=1.5(eα'k+2iαk+2+eβ'k+2iβk+2)Qref+Q2ωk+2Jmin=(Pvk+1-Pref)2+(Pk+2-Pref-P2ωk+2)2+           (Qk+2-Qref-Q2ωk+2)2 (16)

where Pvk+1 is the active power at the valve side at instant k; and P2ωk+2 and Q2ωk+2 are the double fundamental frequency compensation components of active power and reactive power at instant k+2, respectively.

Based on (8) and (9), it can be found that the oscillation component of ΔQ lags that of ΔP by 90 electrical degrees. Therefore, Q2ωk+2 can be obtained as:

Q2ωk+2=P2ωk+2e-T2ωs/4 (17)

where T2ω is the period of double fundamental frequency.

The optimal grid current references can be obtained as:

iα,refk+2=23Pref+P2ωk+2eβ,pre'k+2-Qref+Q2ωk+2eβ,prek+2eα,prek+2eβ,pre'k+2-eα,pre'k+2eβ,prek+2=iα,0k+2+iα,2ωk+2iβ,refk+2=23-Pref+P2ωk+2eα,pre'k+2+Qref+Q2ωk+2eα,prek+2eα,prek+2eβ,pre'k+2-eα,pre'k+2eβ,prek+2=iβ,0k+2+iβ,2ωk+2 (18)

where iα,0k+2 and iβ,0k+2 are the current references corresponding to Pref and Qref, respectively; and iα,2ωk+2 and iβ,2ωk+2 are the compensation references corresponding to P2ωk+2 and Q2ωk+2, respectively. According to (11) and (17), it can be derived that iα,2ωk+2 and iβ,2ωk+2 are the fundamental frequency components composed of positive-sequence and negative-sequence currents. The grid current will maintain sinusoidal and undistorted when the compensation references are applied in the MPDPC.

Combining (14)-(18), P2ωk+2 can be obtained as:

P2ωk+2=K1-K4-K3Q2ωk+2K2 (19)
K1=2Ts3(Pref-Pprek+1)-2L12iα,prek+1iβ,prek+1-L11(iα,prek+1)2-L22(iβ,prek+1)2 (20)
K2=23N1eβ,pre'k+2-N2eα,pre'k+2eα,prek+2eβ,pre'k+2-eα,pre'k+2eβ,prek+2 (21)
K3=23-N1eβ,prek+2+N2eα,prek+2eα,prek+2eβ,pre'k+2-eα,pre'k+2eβ,prek+2 (22)
K4=N1ia,0k+2+N2iβ,0k+2 (23)
N1=-L11iα,prek+1-L12iβ,prek+1 (24)
N2=-L22iβ,prek+1-L12iα,prek+1 (25)

After obtaining the current references at instant k+2 by (18), the output voltage references at instant k+1 can be calculated through (14) and (15).

IV. Identification Scheme for Unbalanced Filter Inductance

As mentioned above, the filter inductance plays an important role in obtaining the power compensation components and the converter voltage references. Once the filter inductance deviates from the initial value, it will affect the dynamic performance and fail to reach the control objectives. Hence, it is necessary to estimate the accurate value of Lαβ online for better control performance.

According to (14), the grid current can be derived as:

iαkiβk=TsB11B12B12B22uLαk-1uLβk-1+iαk-1iβk-1=TsB+ik-1 (26)

where B is the inverse of the actual Lαβ; B11, B12, and B22 are elements in B, which are the parameters that need to be estimated; and [uLαk-1,uLβk-1]T=uLk-1 is the voltage drop on the filter inductance at instant k-1.

The parameter matrix B is estimated online and updated during each period. Assuming that the estimated B at instant k-1 deviates from the actual value by ΔB, then (26) can be rewritten as:

îαk+Δiαkîβk+Δiβk=TsB̂11k-1B̂12k-1B̂12k-1B̂22k-1+ΔB11ΔB12ΔB12ΔB22uLαk-1uLβk-1+iαk-1iβk-1 (27)

where the superscript ˆ represents the estimated value; and [Δiαk,Δiβk]T=Δik is the estimated error of grid currents. Accordingly, the relationship between ΔB and Δik satisfies:

ΔB11ΔB12ΔB12ΔB22uLαk-1uLβk-1=ΔiαkΔiβk (28)

The estimated parameters in B can be updated by calculating ΔB. However, (28) is an underdetermined equation and cannot be solved to obtain a unique solution. Generally, the change frequency of the filter inductance is far less than the fundamental frequency. Therefore, B11, B12, and B22 can be assumed to be constant within a short time, and we can use the estimated and measured information at instant k-n to fill in (28). After some transformations, ΔB can be obtained by:

uLαk-1uLβk-100uLαk-1uLβk-1uLαk-n-1uLβk-n-100uLαk-n-1uLβk-n-1ΔB11ΔB12ΔB22=1TsΔiαkΔiβkΔiαk-nΔiβk-n (29)

where Δiαk-n and Δiαk-n are the estimated errors obtained by B̂k-1.

It can be observed that (29) is an overdetermined equation, through which we can obtain the least square solution of ΔB. The solving process is simple and does not bring a high computational burden, hence it will not be expanded here. Since the delay step n determines the correlation between the two sets of functions in (29), when the correlation is high, the overdetermined equation will become an ill-conditioned equation, making the estimated results extremely sensitive to the measurement or calculation errors. To avoid excessive correlation between the information at instant k-n and that at instant k, the instant k-n should not lag the instant k by 0 or 180 electrical degrees.

The estimated values of B can be updated as:

B̂11kB̂12kB̂22k=B̂11k-1B̂12k-1B̂22k-1+GΔB̂11ΔB̂12ΔB̂22 (30)

where 0<G1 is the update step size.

Then, the estimated inductance matrix can be expressed as:

L̂αβk=B̂11kB̂12kB̂12kB̂22k-1 (31)

The L̂αβk and B̂k are used in the proposed MPDPC and will be updated again at the next control period.

Figure 2 shows the control diagram of the proposed MPDPC, in which the control objective can be selected by setting the compensation power and the blocks “Z-1” and “Z-n” mean the unit delay processes. The sampling and control process is shown in Appendix A. The block diagram of the filter inductance identification scheme is shown in Fig. 3, where the block “inv” means the inversion of a matrix. It can be found that the slightly off-nominal frequency mainly affects the delay process. Since the obtained delay values are close to the actual values, it will have little effect on the performance of the proposed MPDPC.

Fig. 2 Control diagram of proposed MPDPC.

Fig. 3 Block diagram of filter inductance identification scheme.

V. Simulation Results

To verify the validity of the proposed MPDPC, a simulation model of GCC is established on PSCAD/EMTDC. The system parameters are listed in Table I.

TABLE I System Parameters
ParameterValue
DC-side voltage udc 850 V
Root mean square of line-line voltage 380 V
Active power reference Pref 50 kW
Reactive power reference Qref 0 kvar
Grid angular frequency ω 100π rad/s
Sampling and control period Ts 100 μs
Initial value of inductance identification L̂0 1 mH
Update step size of inductance identification G 0.1
Delay step of inductance identification n 50

A. Validity of Proposed MPDPC Under Unbalanced Grid Conditions

In the conventional MPDPC, the unbalanced grid condition of filter inductance is ignored, and only the average value is used to obtain the voltage vector. The analysis in Section II demonstrates that the unbalanced filter inductance will degrade the control performance and generate power oscillation. To solve these problems, a new converter voltage control function (14) is used in the proposed MPDPC. The conventional MPDPC in [

25] is compared with the proposed MPDPC under the same simulation conditions, and the vector selection process is replaced with SVM to obtain better performance.

The unbalanced filter inductances are set as La=2 mH, Lb=6 mH, Lc=4 mH, respectively. In this subsection, the initial value of the average filter inductance and the actual inductance matrix Lαβ are used in the traditional MPDPC and the proposed MPDPC, respectively. The filter inductance identification scheme will be simulated and verified in the next subsection. Before 2 s, the grid voltage is balanced, and the control objective is to eliminate the active power oscillation at the grid side. At 2 s, the grid voltage of phase A is reduced by 20%. At 2.2 s, the control objective is changed to eliminate the active power oscillation at the valve side.

Figure 4 shows the simulation results, including the grid-side active power and reactive power P and Qnov, valve-side active power Pv, compensation active power reference P2ω, DC-side current Idc, grid voltage and current ei and ii (i=a, b, c), and the fast Fourier transform (FFT) of phase-A current with the proposed MPDPC. It can be observed from Fig. 4(a), even if the grid voltage is balanced, the unbalanced filter inductance will generate power oscillation at the grid side and valve side.

Fig. 4 Simulation results of conventional and proposed MPDPC with unbalanced filter inductance and grid conditions. (a) Simulation results of conventional MPDPC. (b) Simulation results of proposed MPDPC. (c) FFT analysis of phase-A current with proposed MPDPC.

The conventional MPDPC cannot eliminate the power oscillation and track the power references. The simulation results in Fig. 4(b) indicate that the grid-side power oscillation can be effectively eliminated by the proposed MPDPC. Under the influence of unbalanced filter inductance and grid voltage, the power oscillation at the valve side and grid side cannot be eliminated simultaneously. The valve-side active power oscillation and DC-side current ripple can be eliminated by the power compensation method. Meanwhile, the FFT analysis shown in Fig. 4(c) also indicates that the proposed MPDPC can maintain sinusoidal grid currents under unbalanced grid conditions. In addition, it should be noted that the proposed MPDPC can also work when the reactive power reference is not equal to 0.

B. Validity of Filter Inductance Identification Scheme with Variable Filter Inductance

To demonstrate the validity of the filter inductance identification scheme and further illustrate the adaptability of the proposed MPDPC to variable filter inductance, two simulation cases are formulated in this subsection.

In the first case, the identification scheme is applied at 1.9 s, with the initial value L̂11=L̂22=1 mH, L̂12=0, update step size G=0.1, and delay step n=50 (90 electrical degrees). The other simulation conditions are the same as Section V-A. Figure 5 shows the corresponding simulation results, including grid-side active and reactive power P and Qnov, valve-side active power Pv, DC-side current Idc, estimated and true inductances L̂ij and Lij, and estimated errors of grid current Δi. It can be observed from Fig. 5 that after the identification scheme is enabled, the estimated inductance values can converge to their true values with high accuracy and speed. The estimated errors of grid current can also be eliminated. Meanwhile, the active power and reactive power can track their references to achieve different control objectives under both balanced and unbalanced grid conditions.

Fig. 5 Simulation results with inductance identification scheme.

In the second case, the filter inductances are set as variable parameters. Before 0.6 s, the filter inductances are balanced with the initial values La=Lb=Lc=1 mH. After 0.6 s, the inductances change every 0.3 s according to the following formula:

La=5sin(π(t-0.6))+1Lb=5sin(3π(t-0.6))+1Lc=5sin(2π(t-0.6))+1 (32)

The grid voltage is unbalanced (ea is reduced by 20%), and the control objective is to eliminate the valve-side power oscillation and DC-side current ripple. From the simulation results shown in Fig. 6, the estimated inductances can track the time-varying values, and the valve-side power oscillation and DC-side current ripple are effectively eliminated by the proposed MPDPC.

Fig. 6 Simulation results with variable filter inductance when n=50, G=0.1.

C. Performance of Filter Inductance Identification Scheme with Different Parameters

The performance of the filter inductance identification scheme is associated with two main parameters: the update step size G and the delay step n.

It is clear that G determines the convergence speed, n determines the similarity of the functions in the overdetermined equation, and further affects the solving accuracy of (29). To show the impact of G and n, a series of simulation tests are conducted with different parameters. The simulation conditions are the same as the second case of Section V-B. Figures 7 and 8 show the identification results with different G and n, respectively.

Fig. 7 Filter inductance identification results with different update steps. (a) n=50, G=0.01.(b) n=50, G=0.5. (c) n=50, G=1.

Fig. 8 Filter inductance identification results with different delay steps. (a) n=1, G=0.1. (b) n=25, G=0.1. (c) n=75, G=0.1.

It can be observed from Figs. 6 and 7 that when G increases, the convergence rate of the identification scheme becomes higher. However, a smoother and better transient response is supposed to be obtained when G has a smaller value. Figures 6 and 8 show that the identification accuracy is significantly degraded when the delay step is too small or close to 0 electrical degrees. This is because in these situations, the correlation of functions in (29) is very high, and a tiny measurement error will lead to large calculation deviations. When the delay step is close to 90 electrical degrees, the two sets of functions in (29) are orthogonal and will provide more information about the filter inductance, hence the identification results can be obtained with higher accuracy. Therefore, it is recommended to set the delay step to 90 electrical degrees, and the update step size should be selected with a bargain of the convergence rate and transient response of the identification process.

VI. Conclusion

In this paper, an MPDPC for the GCCs considering unbalanced filter inductance and grid conditions is proposed. The main conclusions are as follows.

1) The unbalanced filter inductance will generate negative-sequence current and double fundamental-frequency power oscillation. Besides, the mismatched inductance will significantly degrade the control performance of conventional MPDPC.

2) The proposed filter inductance identification scheme can estimate the mismatched and unbalanced inductance with high accuracy and good dynamic performance. The estimation process is simple and will not bring a high computational burden.

3) Compared with the conventional MPDPC in existing research, the proposed MPDPC can effectively eliminate the current ripple at the DC side with unbalanced grid conditions.

Appendix

Appendix A

MPC works in discrete time domination, and the sampling and control process shares one cycle. For each instant k, the sensors first obtain the real-time sampling values of grid voltage ek and grid current ik. Then, the prediction process is conducted to predict eprek+1, eprek+2, and iprek+1. The filter inductance estimation value is also updated by the inductance identification scheme. Finally, the voltage references will be applied on SVM at the next instant to trigger the valves.

The simplified diagram of the sampling and control process is shown in Fig. A1.

Fig. A1 Simplified diagram of sampling and control process.

The green line and blue line represent the sampled values and predicted values, respectively. The black line and red line represent the control references calculated in the past and instant k, respectively. The yellow dot line represents the desired sinusoidal waves of different variables.

References

1

K. Nishida, T. Ahmed, and M. Nakaoka, “A novel finite-time settling control algorithm designed for grid-connected three-phase inverter with an LCL-type filter,” IEEE Transactions on Industry Applications, vol. 50, no. 3, pp. 2005-2020, May-Jun. 2014. [Baidu Scholar

2

L. Guo, N. Jin, Y. Li et al., “A model predictive control method for grid-connected power converters without AC voltage sensors,” IEEE Transactions on Industrial Electronics, vol. 68, no. 2, pp. 1299-1310, Feb. 2021. [Baidu Scholar

3

R. Kadri, J. Gaubert, and G. Champenois, “An improved maximum power point tracking for photovoltaic grid-connected inverter based on voltage-oriented control,” IEEE Transactions on Industrial Electronics, vol. 58, no. 1, pp. 66-75, Jan. 2011. [Baidu Scholar

4

T. Noguchi, H. Tomiki, S. Kondo et al., “Direct power control of PWM converter without power source voltage sensors,” in Proceedings of 21th IEEE Industry Applications Conference, San Diego, USA, Oct. 1996, pp. 941-946. [Baidu Scholar

5

S. Mariethoz and M. Morari, “Explicit model-predictive control of a PWM inverter with an LCL filter,” IEEE Transactions on Industrial Electronics, vol. 56, no. 2, pp. 389-399, Feb. 2009. [Baidu Scholar

6

S. F. Zarei, H. Mokhtari, M. A. Ghasemi et al., “Control of grid-following inverters under unbalanced grid conditions,” IEEE Transactions on Energy Conversion, vol. 35, no. 1, pp. 184-192, Mar. 2020. [Baidu Scholar

7

M. Abdelrahem, F. Hamadto, A. Garikapati et al., “Multiple-vector direct model predictive control for grid-connected power converters with reduced calculation burden,” in Proceedings of 2019 IEEE International Symposium on Predictive Control of Electrical Drives and Power Electronics, Quanzhou, China, Jun. 2019, pp. 1-6. [Baidu Scholar

8

S. Gao, H. Zhao, Y. Gui et al., “An improved direct power control for doubly fed induction generator,” IEEE Transactions on Power Electronics, vol. 36, no. 4, pp. 4672-4685, Apr. 2021. [Baidu Scholar

9

Y. Zhang, C. Qu, and J. Gao, “Performance improvement of direct power control of PWM rectifier under unbalanced network,” IEEE Transactions on Power Electronics, vol. 32, no. 3, pp. 2319-2328, Mar. 2017. [Baidu Scholar

10

L. A. Serpa, S. Ponnaluri, P. M. Barbosa et al., “A modified direct power control strategy allowing the connection of three-phase inverters to the grid through LCL filters,” IEEE Transactions on Industry Applications, vol. 43, no. 5, pp. 1388-1400, Sept.-Oct. 2007. [Baidu Scholar

11

A. M. Razali, M. A. Rahman, G. George et al., “Analysis and design of new switching lookup table for virtual flux direct power control of grid-connected three-phase PWM AC-DC converter,” IEEE Transactions on Industry Applications, vol. 51, no. 2, pp. 1189-1200, Mar.-Apr. 2015. [Baidu Scholar

12

M. Malinowski, M. Jasinski, and M. P. Kazmierkowski, “Simple direct power control of three-phase PWM rectifier using space-vector modulation (DPC-SVM),” IEEE Transactions on Industrial Electronics, vol. 51, no. 2, pp. 447-454, Apr. 2004. [Baidu Scholar

13

S. Vazquez, J. Rodriguez, M. Rivera et al., “Model predictive control for power converters and drives: advances and trends,” IEEE Transactions on Industrial Electronics, vol. 64, no. 2, pp. 935-947, Feb. 2017. [Baidu Scholar

14

S. Yan, A. Zhang, H. Zhang et al., “Optimized and coordinated model predictive control scheme for DFIGs with DC-based converter system,” Journal of Modern Power Systems and Clean Energy, vol. 5, no. 4, pp. 620-630, Jul. 2017. [Baidu Scholar

15

M. K. K. Prince, M. T. Arif, A. Gargoom et al., “Modeling, parameter measurement, and control of PMSG-based grid-connected wind energy conversion system,” Journal of Modern Power Systems and Clean Energy, vol. 9, no. 5, pp. 1054-1065, Sept. 2021. [Baidu Scholar

16

Y. Zhang, W. Xie, Z. Li et al., “Model predictive direct power control of a PWM rectifier with duty cycle optimization,” IEEE Transactions on Power Electronics, vol. 28, no. 11, pp. 5343-5351, Nov. 2013. [Baidu Scholar

17

D. Choi and K. Lee, “Dynamic performance improvement of AC/DC converter using model predictive direct power control with finite control set,” IEEE Transactions on Industrial Electronics, vol. 62, no. 2, pp. 757-767, Feb. 2015. [Baidu Scholar

18

S. B. Q. Naqvi, S. Kumar, and B. Singh, “Three-phase four-wire PV system for grid interconnection at weak grid conditions,” IEEE Transactions on Industry Applications, vol. 56, no. 6, pp. 7077-7087, Nov.-Dec. 2020. [Baidu Scholar

19

Y. Zhang and C. Qu, “Direct power control of a pulse width modulation rectifier using space vector modulation under unbalanced grid voltages,” IEEE Transactions on Power Electronics, vol. 30, no. 10, pp. 5892-5901, Oct. 2015. [Baidu Scholar

20

J. Eloy-Garcia, S. Arnaltes, and J. L. Rodriguez-Amenedo, “Direct power control of voltage source inverters with unbalanced grid voltages,” IET Power Electronics, vol. 1, no. 3, pp. 395-407, Sept. 2008. [Baidu Scholar

21

L. Shang, D. Sun, and J. Hu, “Sliding-mode-based direct power control of grid-connected voltage-sourced inverters under unbalanced network conditions,” IET Power Electronics, vol. 4, no. 5, pp. 570-579, May 2011. [Baidu Scholar

22

Y. Suh and T. A. Lipo, “Modeling and analysis of instantaneous active and reactive power for PWM AC/DC converter under generalized unbalanced network,” IEEE Transactions on Power Delivery, vol. 21, no. 3, pp. 1530-1540, Jul. 2006. [Baidu Scholar

23

Y. Zhang, J. Jiao, J. Liu et al., “Direct power control of PWM rectifier with feedforward compensation of DC-bus voltage ripple under unbalanced grid conditions,” IEEE Transactions on Industry Applications, vol. 55, no. 3, pp. 2890-2901, May-Jun. 2019. [Baidu Scholar

24

Y. Zhang, J. Liu, H. Yang et al., “Direct power control of pulse width modulated rectifiers without DC voltage oscillations under unbalanced grid conditions,” IEEE Transactions on Industrial Electronics, vol. 65, no. 10, pp. 7900-7910, Oct. 2018. [Baidu Scholar

25

Y. Zhang and C. Qu, “Model predictive direct power control of PWM rectifiers under unbalanced network conditions,” IEEE Transactions on Industrial Electronics, vol. 62, no. 7, pp. 4011-4022, Jul. 2015. [Baidu Scholar

26

H. Yang, Y. Zhang, J. Liang et al., “Sliding-mode observer based voltage-sensorless model predictive power control of PWM rectifier under unbalanced grid conditions,” IEEE Transactions on Industrial Electronics, vol. 65, no. 7, pp. 5550-5560, Jul. 2018. [Baidu Scholar

27

Z. Chen, J. Qiu, and M. Jin, “Adaptive finite-control-set model predictive current control for IPMSM drives with inductance variation,” IET Electric Power Applications, vol. 11, no. 5, pp. 874-884, May 2017. [Baidu Scholar

28

D. Zhou, P. Tu, and Y. Tang, “Multivector model predictive power control of three-phase rectifiers with reduced power ripples under nonideal grid conditions,” IEEE Transactions on Industrial Electronics, vol. 65, no. 9, pp. 6850-6859, Sept. 2018. [Baidu Scholar

29

H. T. Nguyen, E. Kim, I. Kim et al., “Model predictive control with modulated optimal vector for a three-phase inverter with an LC filter,” IEEE Transactions on Power Electronics, vol. 33, no. 3, pp. 2690-2703, Mar. 2018. [Baidu Scholar

30

Y. Zhang, J. Jiao, and J. Liu, “Direct power control of PWM rectifiers with online inductance identification under unbalanced and distorted network conditions,” IEEE Transactions on Power Electronics, vol. 34, no. 12, pp. 12524-12537, Dec. 2019. [Baidu Scholar

31

Y. Zhou, A. Zhang, H. Zhang et al., “A model predictive control with online inductance estimator for three-phase VIENNA rectifiers,” in Proceedings of 46th Annual Conference of the IEEE Industrial Electronics Society, Singapore, Oct. 2020, pp. 4607-4611. [Baidu Scholar

32

Y. Zhang, B. Li, and J. Liu, “Online inductance identification of a pwm rectifier under unbalanced and distorted grid voltages,” IEEE Transactions on Industry Applications, vol. 56, no. 4, pp. 3879-3888, Jul.-Aug. 2020. [Baidu Scholar

33

H. Lu, Y. Wang, Y. Yuan et al., “Online identification for permanent magnet synchronous motor based on recursive fixed memory least square method under steady state,” in Proceedings of 2017 36th Chinese Control Conference, Dalian, China, Jul. 2017, pp. 4824-4829. [Baidu Scholar

34

M. Abdelrahem, C. M. Hackl, and R. Kennel, “Finite set model predictive control with on-line parameter estimation for active frond-end converters,” Electrical Engineering, vol. 100, no. 3, pp. 1497-1507, Jul. 2018. [Baidu Scholar

35

H. Yang, Y. Zhang, J. Liang et al., “Robust deadbeat predictive power control with a discrete-time disturbance observer for PWM rectifiers under unbalanced grid conditions,” IEEE Transactions on Power Electronics, vol. 34, no. 1, pp. 287-300, Jan. 2019. [Baidu Scholar

36

S. Ye and X. Yao, “A modified flux sliding-mode observer for the sensorless control of PMSMs with online stator resistance and inductance estimation,” IEEE Transactions on Power Electronics, vol. 35, no. 8, pp. 8652-8662, Aug. 2020. [Baidu Scholar

37

X. Zhang, L. Zhang, and Y. Zhang, “Model predictive current control for PMSM drives with parameter robustness improvement,” IEEE Transactions on Power Electronics, vol. 34, no. 2, pp. 1645-1657, Feb. 2019. [Baidu Scholar

38

X. An, G. Liu, Q. Chen et al., “Adjustable model predictive control for IPMSM drives based on online stator inductance identification,” IEEE Transactions on Industrial Electronics. doi: 10.1109/TIE.2021.3076718 [Baidu Scholar

39

Y. Zhang, J. Jin, and L. Huang, “Model-free predictive current control of PMSM drives based on extended state observer using ultralocal model,” IEEE Transactions on Industrial Electronics, vol. 68, no. 2, pp. 993-1003, Feb. 2021. [Baidu Scholar

40

X. Liu, L. Qiu, W. Wu et al., “Efficient model-free predictive power control for active front-end modular multilevel converter,” International Journal of Electrical Power & Energy Systems, vol. 132, p. 107058, Nov. 2021. [Baidu Scholar

41

S. Saadatmand, P. Shamsi, and M. Ferdowsi, “Power and frequency regulation of synchronverters using a model free neural network-based predictive controller,” IEEE Transactions on Industrial Electronics, vol. 68, no. 5, pp. 3662-3671, May 2021. [Baidu Scholar