Three-Vector-BasedLowComplexityModelPredictiveControl for Soft Open Point

Soft open point (SOP) can improve the flexibility and reliability of power supplies; thus, they are widely used in distribution network systems. Traditional single-vector model predictive control (SV-MPC) can quickly and flexibly control the power and current at both ports of the SOP. However, SV-MPC can only select one voltage vector in a sampling time, producing large current ripples, and power fluctuations. In order to solve the above problems, this paper proposes a three-vector-based low complexity model predictive control (TV-MPC). In the proposed control method, two effective voltage vectors and one zero voltage vector are selected in a sampling time. For the two-port SOP, methods are given to judge the sectors on both sides and select the voltage vectors. Furthermore, the calculationmethod of the distribution time is proposed as well. Finally, the effectiveness of the proposed method is verified by steady-state and dynamic-state simulation results compared with the SV-MPC.


Introduction
In recent years, more and more renewable energy sources such as solar energy and wind energy as well as new energy vehicles loads have been connected to the distribution network, which has had a serious impact on the stable and safe operation of the distribution network [1,2]. erefore, soft open point (SOP) is used more frequently in distribution networks due to their advantages to flexibly connect feeders of different voltage levels, improve the reliability of power supply, and continuously adjust power [3][4][5].
e two-port SOP can be regarded as an AC/DC/AC converter composed of voltage source converters (VSCs) [6].
e power exchange at both ports can be realized through corresponding control. In [7], droop control is adopted and proportional-integral (PI) controllers are designed to realize the closed-loop control of the outer power loop and the inner current loop. In [8], the inner current loop of the traditional PI double-loop control is replaced by a sliding mode control, which can reduce the PI controllers. e combination of sliding mode control and feedback linearization control in [9] can meet the requirements of different operating modes. However, control applications will be difficult if a large number of PI coefficients need to be tuned or many controller parameters need to be selected.
Model predictive control (MPC) is widely used in power electronic converters due to its simple implementation and fast dynamic response. According to different control purposes, MPC can be divided into model predictive current control (MPCC) and model predictive power control (MPPC). MPCC and MPPC are used in rectifiers, inverters, modular multilevel converters, induction machines, and permanent-magnet synchronous machines [10][11][12][13][14]. However, the single-vector MPC has only one voltage vector involved in the control during a sampling time, and large current ripples and power fluctuations will be generated, resulting in unsatisfactory control effects. erefore, how to improve power quality and reduce power fluctuations needs to be considered when MPC is used in converters. In order to improve the control performance, a variety of improved MPC methods have been proposed, such as implementing delay compensation, or multistep prediction, improving the cost function, and increasing the number of voltage vectors in a sampling time [15][16][17][18][19]. Among them, increasing the number of voltage vectors is the main way to improve the control performance. In [20], based on the dual-vector MPCC, a method of projecting the current error vector onto the active voltage vector is proposed, which minimizes the cost function and reduces the number of candidate voltage vectors. In [21], different sizes of virtual vectors are defined, the obtained reference voltage vector and the two predicted vectors selected from the virtual vector are used for evaluation, and the virtual vector that minimizes the cost function is used in the next sampling time.
is method greatly alleviates the computation burden and avoids the weighting factor in predictive torque control. In [22], a three-vector MPCC is proposed for an openwinding linear permanent-magnet vernier motor. Based on the principle of deadbeat current control, an optimal vector is obtained, and the other two optimal vectors are obtained by cascading cost functions. In addition, the switching frequency is reduced to improve the steady-state performance and suppress the zero-sequence current. In [23], a combination of MPC and direct power control is used for the improved T-type three-level converter, and the switching states are grouped in advance. e capacitor voltage is compared on the DC side to determine the optimal group of the switching state, which eliminates the need for midpoint voltage prediction and cost function calculations, while ensuring the neutral point voltage balance. In [24], the adaptive error correction strategy is applied to both the outer and inner prediction loops from the perspective of model errors, reducing the impact of errors caused by time delay, sampling error, and parameter mismatch on the control performance and effectively improving the output power quality compared to the traditional FCS-MPC.
Both sides of the traditional MPC strategy for the SOP only consider a single switch state control within a sampling time, which results in large current ripples and power fluctuations in the steady state. e main contribution can be described as follows: (1) A three-vector-based low complexity MPC scheme is proposed, which reduces the current ripple and power fluctuation by increasing the number of switch states within a sampling time (2) Methods for judging sectors and selecting vectors are given, which reduce the number of calculations In Section 2, the mathematical model of the SOP is introduced and the power flow is analyzed. en, the traditional MPC method for the SOP is introduced, and the control block diagram of the traditional method is given in Section 3. In Section 4, a three-vector-based low complexity MPC method is presented, and methods for judging sectors and selecting vectors at both sides are given. In Section 5, the traditional method and the improved method are compared by simulations, which shows that the proposed method effectively reduces the current ripple and power fluctuation and has good dynamic characteristics. Finally, the article is summarized in Section 6. Figure 1 shows the topology of a two-port SOP, where the rectifier side and the inverter side are interconnected by the capacitance of the DC side.

Mathematical Model of SOP
According to Kirchhoff's law, the state equation of the rectifier side and the inverter side in the motionless reference framework can be obtained as where V x and i x1 are the grid voltage and input current of the rectifier side, respectively; e x and i x2 are the load backemf and output current of the inverter side, respectively; the subscripts x � a, b, c represents the three phases of the SOP; L s and R s are the filter inductance and resistance; and L l and R l are the inverter inductance and resistance, respectively. e switching states of the SOP is defined as where S x � 1 means that the upper switch of phase x is on and the lower switch is off and S x � 0 means that the upper switch of phase x is off and the lower switch is on. en, introducing the definition of the switching states into (1) and (2), the values of output voltages of each converter can be determined as where u dc is the capacitor voltage in the DC link; V aiN , V biN , and V ciN are the output voltages of each side in phases a, b, and c, respectively; and i � 1 means the rectifier side and i � 0 means the inverter side.
In the system of the two-port SOP, eight VSC voltage vectors can be chosen, six of which are active voltage vectors (V 1 , V 2 , V 3 , V 4 , V 5 , V 6 ) and two of which are zero vectors (V 0 , V 7 ). ese eight vectors are shown in Figure 2.
According to the instantaneous reactive power theory, ignoring the line loss and the self-loss of the SOP, the output active power P and reactive power Q of each SOP port can be obtained as where P 1 and Q 1 are the rectifier side active and reactive power, respectively, while P 2 and Q 2 are the inverter side active and reactive power, respectively. V α , V β , i 1α , and i 1β are the grid voltage and input current in the stationary αβ reference frame of the rectifier side, respectively, while e α , e β , i 2α , and i 2β are the load back-emfs and output currents in the αβ reference frame of the inverter side, respectively.
When the SOP operates normally, the output active power of the rectifier side P s can be calculated by the theorem of conservation of power as P s � P dc + P l , where P dc is the active power obtained by the DC side capacitor and P l is the active power obtained by the load of the inverter side.

Outer-Loop Control Mode Selection.
Each port of the SOP can work in different control modes, namely, PQ mode, U dc Q mode, and U dc f mode. One port of the VSC should be needed to control the DC side voltage stability. When the SOP is working normally, its two ports generally work in PQ mode and U dc Q mode. To ensure normal operation and meet load power supply requirements when a feeder fails, the U dc f mode is set and normally used in SOP systems with three or more ports. In this study, the outer loop of the rectifier side selects U dc Q mode, and the inner loop selects a direct power model predictive control (DPMPC). e outer loop of the inverter side selects a PQ mode, and the inner loop selects the MPCC to reduce the load current error.

Analysis of DPMPC Method in Rectifier
Side. By using the Clark transformation, Equation (1) can be converted to the stationary αβ reference frame as follows: where Using the forward Euler discretization method to predict the current in the next control period of (7), the predicted current at the (k+1)th instant can be calculated as where i 1αβ (k) and V αβ (k) are the input current and grid voltage at the kth instant, respectively; V 1Nαβ (k) is the voltage vector of the VSC in the rectifier side at the kth instant; and T s is the sampling time.
Rectifier side Inverter side Figure 1: Two-port SOP equivalent circuit. In the U dc Q mode, the rectifier side uses a part of the power obtained from the grid to exchange with other ports, and the other part is used to adjust the DC side voltage. e outer loop adopts PI control, and the active power reference value of this port P ref is as follows: where k p and k i are the proportional and integral coefficients of PI control, respectively; u dcref and u dc are the reference voltage and sampling voltage of the DC side, respectively; and i 2ref is the conjugate load reference current of the inverter side. e first two items in (9) calculate the power required to adjust the DC side voltage, and the active power obtained by the load on the inverter side is calculated by the last item. e voltage vectors and the predicted current vectors can be used to calculate the predicted instantaneous input active power and reactive power as follows: In (10), when the sampling time is small, it can be assumed that V αβ (k + 1) ≈ V αβ (k), but if the sampling time is not small enough to ignore the change in the grid voltage within the sampling time, vector angle compensation can be used to calculate the grid voltage at the (k+1)th instant as follows: where ω is the grid voltage pulsation. en, Equation (12) can be used as a cost function to evaluate the prediction error generated by each switch, which is described as In DPMPC, all eight voltage vectors are substituted into (8) to predict the current at the (k+1)th instant, the predicted powers are substituted into (12) to calculate the error, and the switch with the smallest cost function is selected as the optimal switch and is applied to the SOP rectifier side at the next moment.

Analysis of MPCC Method in Inverter
Side. MPCC is used in the inverter side to compare the predicted load current with the reference load current and select the voltage with the smallest cost function to apply in the next control period.
Transforming (2) to the stationary αβ reference frame and adopting the forward Euler discretization method, the predicted load current at the (k+1)th instant is as follows: T is the voltage vector of the VSC in the inverter side at the kth instant. In order to adjust the load current, a cost function is used to find the best switching state: (14) Figure 3 shows the traditional SOP MPC block diagram. At each sampling time, Equations (12) and (14) are used to calculate the cost function values of all switch states on both sides. e switch state that minimizes the cost function is selected to drive the SOP.

Three-Vector-Based Low Complexity MPC Method for SOP
In the traditional MPC, the VSCs on both sides adopt a singlevector predictive control, the vector selection is very limited, and the voltage vector output by the two ports cannot reach the entire spatial circle track range. As a result, the current ripple content and power fluctuations on both sides become larger. In order to improve the accuracy of control, it is necessary to increase the sampling frequency of the system, thereby increasing the computational burden of the processor. Another problem is that the switching frequency is not fixed. In order to reduce the current ripples and power fluctuations while fixing the switching frequency, a three-vector-based low complexity MPC is used on both sides of the port in this study.

Cost Function Design.
In a three-vector-based MPC, cost functions (12) and (14) can be improved as follows: In these improved cost functions, the error between the predicted value and the reference value is compared, which is conducive to choosing the best switch value.
Using adjacent effective vectors for vector synthesis, six vector combinations can be obtained: In a three-vector-based MPC, since three vectors are synthesized in one sampling time, an additional cost function needs to be designed to comprehensively evaluate the errors of the three vectors in the sampling time. In this study, the weighted root means that the square error is used to design a new cost function: where e 1 , e 2 , and e 0 are the cost function values corresponding to the two effective voltage vectors and zero vectors, respectively; and d 1 , d 2 , and d 0 are the operating times of the three voltage vectors in the sampling time T s , respectively.
Since the sum of squares of errors has been selected as the cost function in (15) and (16), the cost function (17) is written as follows: Taking (18) as the new cost function, the voltage vector time allocation obtained by minimizing the weighted root mean square error is the optimal vector that makes the predicted value of the system closest to the reference value.

Sector Judgment and Vector Selection.
In order to improve the control accuracy and fix the switching frequency, before calculating the duty cycle, it is necessary to judge the reference vector sectors on both sides of the SOP to select the port voltage vector.
On the inverter side, current predictive control is adopted. From the perspective of tracking the reference current, in this study, the inverter side adopts the principle of deadbeat control, which can be derived as follows: where V 2Nαβref and i 2αβref are the reference voltage vectors of the VSC and the reference current, respectively. After the reference voltage is obtained, the selection of the sector depends on the phase angle θ ref of the reference voltage. e phase angle of the reference voltage in the stationary αβ reference frame can be calculated as follows: e sectors and vectors selected according to the phase angle θ ref are shown in Table 1. By judging the sectors, the two effective vectors and zero vectors selected on the inverter side are determined, which effectively reduces the amount of calculation.
On the rectifier side, power predictive control is adopted. In order to achieve power tracking, the variables that affect the value of the cost function (15) are analyzed, and the derivative of (5) can be obtained as Rectifier side

Mathematical Problems in Engineering
For a three-phase balanced grid voltage, the differential of the grid voltage can be expressed as follows: Substituting (22) and (7) into (21), (21) can be rewritten as Using the forward Euler method to discretize (23), the predictive control of the power in the next sampling period can be obtained as where f p and f q are the derivatives of active power P 1 and reactive power Q 1 , respectively. For power reference values P 1ref and Q 1ref , the corresponding voltage vector V 1Nref is an unknown quantity, which can be obtained from (24) as follows: After substituting (25) and (24) into the cost function (15), (15) can be rewritten as where Assuming that the voltage reference value V 1Nref is in sector I, the Euclidean distance between the reference voltage and each effective voltage vector is shown in Figure 4. e two voltage vectors that minimize (V 1Nαref − V 1Nα ) 2 + (V 1Nβref − V 1Nβ ) 2 are V 1 and V 2 . Both V 1 and V 2 are in sector I. erefore, it can be proved that the two effective voltage vectors that minimize the cost function f must be in the same sector as the reference voltage vector V 1Nref , ensuring that the two effective voltage vector combinations with the smallest cost function must be within the six voltage vector combinations listed in Section 4.1.
In the traditional three-vector based MPC, the theory of ergodic optimization is used and all combinations formed by the two active vectors and the zero vector need to be calculated. For the combination of the six active vectors, each side of the SOP needs to calculate 15 times to find the best combinations and action time of each vector at each sampling time. In this paper, the low complexity method proposed above needs only one calculation of the reference prediction voltage on the inverter side, while on the rectifier side, two active voltage vectors with the smallest values of the cost function among the active vectors are proved to be the optimal combination of the active vectors, which requires Mathematical Problems in Engineering only six calculations of cost function and one calculation of action time.

Calculation of Duty Cycle.
In (18), cost function G needs to be minimized to obtain the responsible vector combination; e 1 , e 2 , and e 0 can be regarded as known quantities in the cost function; and d 1 , d 2 , and d 0 are regarded as unknown quantities. It is equivalent to the problem of finding the extreme value of a multivariate function, and the constraint conditions are as follows: For the extreme value problem of a multivariate function with constraints, the Lagrange multiplier method is used to construct the function as follows: where μ is the Lagrange multiplier. e partial derivatives of the three variables of (28) are, respectively, calculated as 0, and the calculation can be obtained as follows: After sector judgment and vector selection, the minimum value of the newly constructed cost function F(d 1 ,d 2 ,d 0 ) can be obtained from (29), and the action time of each vector can be calculated, which is applied to the switch control of the next sampling time.

Stability Analysis.
In the two-port SOP, the two-portindependent system can be divided into a power subsystem and a load subsystem depending on the direction of power flow at both sides. As shown in Figure 5, according to the impedance method, the system can only be stable if the subsystem of each port can be stable and the interaction between two ports remains stable at the same time.
e system can only be stable when the Nyquist contour of the return-ratio matrix X does not encircle the (−1, j0) point. By the flow of power, it can be seen that the power positive flow side is the rectifier side, which belongs to the power subsystem, while the negative flow side is the inverter side, which belongs to the load subsystem; thus, the return-ratio matrix R of the two-port SOP can be expressed as follows: (30) e system ensures stability by keeping the eigenvalues of the return-ratio matrix outside the certain forbidden region. In [25], a new forbidden region is proposed as follows: where λ Xi (i � 1 and 2) are the eigenvalues of the return-ratio matrix, θ PM is the phase margin, and ε GM is the gain margin. As shown in Figure 6, by setting the phase margin and gain margin, λ Xi can be maintained out of the stability forbidden region.

Summary
e control block diagram of the three-vector-based MPC for the SOP is shown in Figure 7. Based on the deadbeat principle, the inverter side can calculate the predicted value of the port voltage V 2Nref , the phase angle θ ref obtained by V 2Nref is used to determine the sector where the vector is located, and the sector is used to select the six vector combinations in Section 4.1. e rectifier side selects the two effective vectors with the smallest value function. It has been shown in Section 4.2 that the two vectors with the smallest  Mathematical Problems in Engineering value function must be in the same sector, thus ensuring that the selected vector combination must be within the six vector combinations listed in Section 4.1.
Comparing Figure 3 with Figure 7, the three-vectorbased method selects the sector and calculates the duty cycle, differing from the single-vector method where there is only one vector control in the sampling time. As a result, the current ripple and power fluctuations can be reduced, while the switching frequency is fixed.
As shown in Figure 7, only one PI controller is used in the voltage outer loop by the method proposed in this paper. Compared with the seven PI controllers used in the traditional PI loop control [26], the number of PI controllers used is greatly reduced and therefore the complexity of the parameter design is reduced. At the same time, compared with the high accuracy requirements of sliding mode control on system parameters and the high sampling frequency and large output harmonics of hysteresis control, the method  Phase angle calculation (20) Vector selection Table 1 Rectifier side  Mathematical Problems in Engineering proposed in this paper can achieve a simpler control process and operate with better performance. In the traditional MPC method, a large number of calculations are required at each sampling time, which results in a delay in system control, while the low complexity algorithm requires few calculations at each sampling time. At the same time, using the vector angle compensation for the reference voltage on the rectifier side improves the control accuracy and has a better real-time performance when the sampling time is large.

Simulation Results
To verify the effectiveness and correctness of the proposed three-vector-based MPC for the SOP, a simulation experiment platform is built in MATLAB/Simulink. e parameters of the SOP are shown in Table 2. e steady-state simulation of the two methods on the rectifier side is shown in Figures 8 and 9. Figure 8 shows the power fluctuations on the rectifier side of the traditional MPC and three-vector-based MPC when the inverter side reference current is 40 A. e power fluctuations of the two methods obtained by the data analysis are shown in Table 3.
ree-vector-based power prediction effectively reduces the amplitude of power fluctuations compared with singlevector power prediction. Figure 9 shows the current waveform and frequency spectrum of the rectifier side and the inverter side. e current total harmonic distortion (THD) of the traditional MPC is 2.08%, whereas the current THD of the proposed three-vector-based MPC is 0.91%. ree-vector power prediction control reduces both the current ripple and the power fluctuation. In the rectifier side, both methods effectively track the inverter side reference current. e   According to the power and current analysis of two methods on the rectifier and inverter sides in Figures 8 and 9, under steady-state conditions, the three-vector-based method has better stability than the single-vector method. erefore, the dynamic characteristics of the three-vectorbased method are simulated and analyzed below.
In the actual application of the SOP, it is necessary to meet the bidirectional power flow operation. Figures 10-12 show the dynamic simulation results when the reference current on the inverter side changes suddenly from 40 A to −20 A at 0.5 s. It can be seen from Figure 10 that the voltage change on the DC side is small and then stabilizes to 800 V again after a short period of fluctuation. e active power and reactive power effectively track the reference power obtained from the inverter side reference current after a short period of fluctuations and meet the conservation of power exchange of two ports in Figure 11. Current waveforms on both sides are shown in Figure 12, and the currents i abc1 and i abc2 reach a stable state after small fluctuations.

Conclusion
is paper proposes a three-vector-based low complexity MPC method for the SOP. First, the power predictive control model on the rectifier side and the current predictive control model on the inverter side are established using the single-vector method. en, in order to reduce the current ripple content and power fluctuations and fix the switching frequency, two effective voltage vectors and a zero vector are used for modulation to increase the selection range of the voltage vector, and the method of voltage vector selection and action time distribution is proposed. Finally, through the simulation and comparison of the two methods, it is concluded that three-vector-based MPC has a smaller current ripple and a higher power control accuracy, while still maintaining excellent dynamic response characteristics.
However, in the case of sudden changes in the reference current and power, both the DC side voltage and the current of each side have a large fluctuation. e parameter mismatches that can occur in the actual operation are also not taken into account. erefore, based on the research in this paper, the following research needs to focus on how to combine with other methods to improve the dynamic performance of the system, while ensuring that the system can still operate properly in the case of parameter mismatch.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.