Disturbance Observer-Based Model Predictive Super-Twisting Control for Soft Open Point

: This paper presents a disturbance observer-based model predictive of super-twisting control for Soft Open Point (SOP). First, with the consideration of the disturbances caused by parameter mismatches and unmodelled dynamics, a super-twisting sliding-mode observer (STO) is proposed to observe the disturbances, and the observed disturbances are introduced into the inner-loop as the compensation to improve the anti-disturbance of SOP system. Second, the outer-loop controller is designed by applying the super-twisting sliding-mode control (STC) approach to improve the dynamic performance and robustness. Third, to deal with large current harmonics by traditional model predictive control (MPC), a Three-Vector-based MPC (TV-MPC) is proposed to increase the number of voltage vectors in a sampling time. Finally, it is veriﬁed by simulations that the proposed method can reduce current harmonics, DC-side voltage setting time and improve the dynamic performance of SOP system effectively. In case of parameter mismatches, the proposed observer can observe the disturbances correctly to enhance the robustness of the SOP system.


Introduction
In recent years, with the popularization and development of renewable energies, increasingly renewable energy sources are connected to electricity distribution networks. However, with the development of renewable energy sources, such problems as voltage instability, bi-directional current and unbalanced power arise in electric energy systems [1,2]. As a new power electronic device, Soft Open Point (SOP) has the advantages of flexible power management of the distribution network, fast response time and various control methods, among others [3]. Therefore, SOP has been widely used in the electricity distribution energy systems.
Two-port SOP can be seen as a back-to-back converter, which consists of two voltage source converters (VSCs) and a DC-side capacitor [4]. The two-port SOP system is depicted in Figure 1. The control proposed for SOP includes traditional proportional-integral (PI) control, droop control and model predictive control (MPC) [5][6][7]. These methods are widely applied in a SOP system. However, traditional PI control requires many parameters to be designed; when disturbances occur in the SOP system, PI control and droop control are less robust as well. In signal-vector-based MPC, only one voltage vector is applied in the system in a sampling time, which causes large current harmonics and power fluctuations. Therefore, it is necessary to design improved inner and outer loop control to improve the steady state and dynamic state of the SOP system.
In order to solve the problems above, various improved methods are proposed for the inner-loop MPC, such as increasing the number of voltage vectors in a sampling time, dual-vector MPC [8] and three-vector-based MPC [9]; combining other methods, such

Mathematical Model of SOP
Since the structure of each port of the SOP is the same, in this paper, one port of the voltage source converter (VSC) is selected to build the mathematical model. Taking VSC1 as an example, the model of one port of the SOP in synchronous rotating frame can be expressed as [36]:   (1) where id1 and iq1 are the d-and q-axis currents of VSC1, respectively. u1d and u1q are the dand q-axis grid voltages of VSC1, respectively. udc is DC-side voltage; L1, R1 and ω are inductance, resistance of VSC1 and voltage angle frequency, respectively. Sd1 and Sq1 are the d-and q-axis modulation switch functions of VSC1, respectively. Similarly, VSC2 has the same mathematical model in synchronous rotating frame.
According to the current direction shown in Figure 1, the relationship between DCside currents from VSC and from capacitance can be expressed as follows: where C is the capacitance of DC side, idc1 and idc2 are the DC-side currents from VSC1 and VSC2, respectively. ic is the DC-side current from capacitance; Sd2 and Sq2 are the d-and qaxis modulation switch functions of VSC2, respectively. id2 and iq2 are the d-and q-axis currents of VSC2, respectively. By discretizing (1) [37], the discrete model of VSC1 can be obtained as: where id1(k + 1) and iq1(k + 1) are the predicted d-and q-axis currents of VSC1 at (k + 1)th instant, respectively. u1Nd(k) and u1Nq(k) are the d-and q-axis voltages of VSC1, respectively. id1(k) and iq1(k) are the d-and q-axis currents of VSC1 at (k)th instant, respectively. Since the MPC of the inner loop is created by accurate model parameters, when disturbances of parameter mismatches and unmodelled dynamics occur in the SOP system in actual operation, predictive errors may occur in the system, due to MPC being sensitive to the parameters, which leads to bad performance of the system [15][16][17]. Thus, observers are applied in the system to observe the disturbance of the system and considering the disturbance as the compensation of the system to improve the robustness of SOP. Adaptive observers are used in [18] to estimate the capacitors' voltages, which have aa good performance during disturbances. In [19], a Kalman filter is designed to predict the current with reduced noises and estimate the disturbance due to the parameter variation, combined with extended state modeling (ESM). In [20], an extended state observer is applied to enhance the robustness of inductance to cater for the divergence caused by the stator inductance mismatch.
Sliding mode theory is widely used in the design of controllers and observers due to its fast response time and robustness [21][22][23][24]. In [25], a sliding-mode controller is designed for a grid-connected NPC inverter to inject a controlled active power while minimizing supply current harmonics, which has the advantages of fast dynamic response and low voltage ripple. Sliding-mode observer is designed for disturbances' estimation in [26], while slidingmode controllers are designed to solve the consensus problem in second-order nonlinear multi-agent systems. However, due to its own discontinuity, first-order sliding-mode theory has obvious chattering. Super-twisting algorithm (STA), as a second-order slidingmode theory, can reduce the chattering effectively, while maintaining the advantages of traditional sliding mode [27][28][29][30][31][32]. In [33], a super-twisting controller is used in permanent linear synchronous motor (PMLSM) to ensure the system convergence and eliminate the chattering phenomenon, thus improving the system's robustness. In [34], a super-twisting nonlinear disturbance observer is proposed, which can estimate the uncertainties and external disturbances. In [35], a super-twisting observer is designed to observe the lump disturbance in the permanent magnet synchronous motor (PMSM) caused by parameter mismatches and consider the disturbances as the compensation to the inner-loop, which improves the robustness of the system. In this paper, by introducing the theory of the super-twisting algorithm and threevector-based MPC, a super-twisting sliding-mode controller (STC) and TV-MPC with supertwisting sliding-mode observer (STO) method is proposed to improve the performance of the steady state and dynamic state, and to observe the disturbance correctly under parameter mismatches, the disturbances observed by STO are applied to the TV-MPC as the compensations of the inner-loop. The contributions of this article are listed as follows: (1) Compared with traditional outer-loop PI controller, the start-up time and dynamic response of the DC-side voltage are improved by adopting the proposed method; (2) Compared with signal-vector method, the current harmonics are improved by applying the TV-MPC method in SOP system; (3) Parameter mismatches in the SOP system are considered and compensated for by STO, which improves the adaptation of the SOP system to the disturbances.
The organization of this paper is as follows: Section 2 introduces the model of the two-port SOP system. In Section 3, STO is designed to observe the disturbances. Section 4 designs the outer-loop DC-side voltage STC. Section 5 proposes the TV-MPC method of the SOP system. Section 6 gives the simulation results. Future works are discussed in Section 7. Finally, the conclusions are summarized in Section 8.

Mathematical Model of SOP
Since the structure of each port of the SOP is the same, in this paper, one port of the voltage source converter (VSC) is selected to build the mathematical model. Taking VSC1 as an example, the model of one port of the SOP in synchronous rotating frame can be expressed as [36]: where i d1 and i q1 are the dand q-axis currents of VSC1, respectively. u 1d and u 1q are the dand q-axis grid voltages of VSC1, respectively. u dc is DC-side voltage; L 1 , R 1 and ω are inductance, resistance of VSC1 and voltage angle frequency, respectively. S d1 and S q1 are the dand q-axis modulation switch functions of VSC1, respectively. Similarly, VSC2 has the same mathematical model in synchronous rotating frame.
According to the current direction shown in Figure 1, the relationship between DC-side currents from VSC and from capacitance can be expressed as follows: where C is the capacitance of DC side, i dc1 and i dc2 are the DC-side currents from VSC1 and VSC2, respectively. i c is the DC-side current from capacitance; S d2 and S q2 are the dand q-axis modulation switch functions of VSC2, respectively. i d2 and i q2 are the dand q-axis currents of VSC2, respectively. By discretizing (1) [37], the discrete model of VSC1 can be obtained as: where i d1 (k + 1) and i q1 (k + 1) are the predicted dand q-axis currents of VSC1 at (k + 1)th instant, respectively. u 1Nd (k) and u 1Nq (k) are the dand q-axis voltages of VSC1, respectively. i d1 (k) and i q1 (k) are the dand q-axis currents of VSC1 at (k)th instant, respectively. u 1d (k) and u 1q (k) are the dand q-axis grid voltages of VSC1 at (k)th instant, respectively. T s is the sampling time.
Each port of SOP can work in a different mode. DC-side voltage and reactive power are taken as the reference parameters in U dc Q mode, in order to ensure the steadiness of the DC-side voltage, at least one port should work in U dc Q mode in the SOP system. In PQ mode, the active power and reactive power are taken as reference parameters and are most used in the SOP system. In addition, taking the DC-side voltage and frequency as reference parameters the U dc f mode can also be applied in the SOP system when system failure occurs in the system. In this paper, the U dc Q mode is applied in port1 and PQ mode is applied in port2.

Design of Super-Twisting Observer under Parameter Mismatch
In the actual operation of the SOP system, both sides work under large current and voltage fluctuations frequently, which makes the parameter changes of resistance and inductances. In this paper, a super-twisting observer (STO) is designed in the inner-loop to observe the disturbances caused by parameter mismatches and unmodelled dynamics, which maintains the strong robustness of the traditional sliding mode observer while suppressinges the high frequency vibration effectively.
Considering the parameter mismatches and unmodelled dynamics, the model of VSC1 can be rewritten as follows: where u 1Nd and u 1Nq are the dand q-axis voltages of VSC1, respectively. f d1 and f q1 are the dand q-axis disturbances of VSC1 caused by parameter mismatches and unmodelled dynamics, respectively. Due to disturbances' change slowly in steady-state, f d1 and f q1 can be written as: where ε d1 and ε q1 are the dand q-axis unmodelled dynamics disturbances of VSC1; L 10 and R 10 are offline parameters of inductance and resistance of VSC1, respectively; variable matrix x 1 is defined as: Errors of inductance parameters and resistance parameters ∆L 1 and ∆R 1 are uncertain and bounded, the unmodelled dynamics disturbance ε d1 and ε d1 are bounded, which can be written as: where ϕ L , ϕ R and ϕ εd1 and ϕ εq1 are the known positive constants. Due to the assumption above, the disturbance f d1 and f q1 are bounded and can be assumed as: where ϕ d1 , ϕ q1 and ϕ 1 are known positive constants. The sliding manifold is designed as: whereî d1 andî q1 are the dand q-axis estimated current of VSC1.  (9), the super-twisting observer can be designed as follows: (10) where d d and d q are the d-and q-axis error dynamics of the disturbance; α 1 and β 1 are the designed gains of observer.
The estimated disturbancef 1 can be expressed as follows according to (6): Proof of Proposition 1. Taking VSC1 as analysis, the Lyapunov function is designed as [38,39]: The time derivative of (12) is expressed as: . where Using the bound of disturbance f 1 (8), bound of . V(γ) can be obtained as follows: . where is negative, and the system is strongly globally asymptotically stable when α 1 and β 1 are satisfied as follows: Proofs of VSC2 are similar to the VSC1 side. The proof of the observer stability is completed, which means the values of estimated current approach the values of system current in finite time.
Proofs of VSC2 are similar to the VSC1 side. The proof of the observer stability is completed, which means the values of estimated current approach the values of system current in finite time. □ Figure 2 illustrates the block diagram of STO. The initial values of the integrators are set as 0.

Design of Outer-Loop Super-Twisting Controller Design
In this paper, port1 is operated in UdcQ mode, which keeps the stability of the DCside voltage. However, in the traditional outer-loop PI controller, the constant parameters

Design of Outer-Loop Super-Twisting Controller Design
In this paper, port1 is operated in U dc Q mode, which keeps the stability of the DC-side voltage. However, in the traditional outer-loop PI controller, the constant parameters of controller may lead to poor dynamic performance when the parameters of the system are changed. As a second-order sliding mode controller, the super-twisting controller (STC) can effectively suppress the chattering of the first-order sliding mode controller while having a better dynamic performance compared with the PI controller. An outer-loop super-twisting voltage controller is designed in this section.
The outer-loop system of VSC1 is composed of the DC voltage loop, the sliding manifold is designed as: Uncertain disturbance δ is unknown, but can be assumed to be bounded as follows: where ρ is a known positive constant.

Proposition 2.
In order to track the reference DC-side voltage, according to the theory of supertwisting algorithm (STA), when u dc is the system-state variable, STA of u dc can be described as: where k 1 > 0 and k 2 > 0 are the controller gains.
When the grid voltages are three-phase balanced and the system is stable, the state of the dand q-axis currents of AC side i d and i q can be written as:  (2) and (1), which can be rewritten as: where u 2d is the d-axis grid voltage of VSC2. Substituting (21) into (20), i d1 can be rewritten as: when the system is stable, i d1 = i d1ref , (22) can be rewritten as: Substituting (18) into (23), the output of the controller id1ref can be written as: Proof of Proposition 2. The Lyapunov function is designed as: where The time derivative of (25) is expressed as: . where Using the bound of disturbance δ (17), the bound of . V(ζ) can be obtained as follows: . where . V(ζ) is negative and the system is strongly globally asymptotically stable when k 1 and k 2 are satisfied as follows: When k 1 and k 2 are satisfied in the equation above, the error dynamic of DC-side voltage approaches the ideal sliding manifold S = 0 in finite time. The proof is completed. Figure 3 shows the block diagram of outer-loop STC.
When k1 and k2 are satisfied in the equation above, the error dynamic of DC-side voltage approaches the ideal sliding manifold S = 0 in finite time. The proof is completed. □ Figure 3 shows the block diagram of outer-loop STC.

Traditional MPC Method
In the SOP system consisting of two-level VSC, the switch states can be divided into six active voltage vectors (V1, V2, V3, V4, V5, V6) and two zero vectors (V0, V7), as shown in Figure 4. As shown in Figure 1, '1' in Figure 4 means the switch S x is on and S x is off, while '0' means the switch S x is off and S x is on (x = 1, 2, 3, 4, 5, 6).

Traditional MPC Method
In the SOP system consisting of two-level VSC, the switch states can be divided into six active voltage vectors (V1, V2, V3, V4, V5, V6) and two zero vectors (V0, V7), as shown in Figure 4. As shown in Figure 1, '1' in Figure 4 means the switch Sx is on and S′x is off, while '0' means the switch Sx is off and S′x is on (x = 1, 2, 3, 4, 5, 6). Port1 is operated in UdcQ mode to ensure the stability of the DC-side voltage. In the traditional MPC method, the PI controller is used in the outer-loop subsystem to output d-axis reference current id1ref, which can be expressed as: where kp and ki are the proportional and integral coefficients of PI controller, respectively. udcref is the reference voltage of the DC side. The PQ mode is used in Port2. With a stable grid voltage, a given reference current is equivalent to a given reference power, which can transform the PQ mode control into a problem of tracking the reference current.
Taking VSC1 for analysis, since two zero vectors V0 and V7 obtain the same voltage vector, in traditional MPC, seven voltage vectors are substituted into (3) to calculate the corresponding predicted currents, and each predicted current is substituted into the cost function, the voltage vector with the smallest cost function value can be selected as the optimal voltage vector and applied to the switch in the next sampling time. The cost function can be designed as: where iq1ref is the q-axis current reference value of VSC1. Similarly, VSC2 with given reference current, cost function values can be calculated and the optimal voltage vector will be applied to the switch in the next sampling time.

TV-MPC Method
In the traditional MPC method, only one voltage vector can be used in a sampling time and seven times of cost function calculations are required at both port, which leads to large current harmonics and control delays. Therefore, using the principle of deadbeat control, a low-complexity TV-MPC method is proposed. By calculating the predicted value of the reference voltage of VSC at the next sampling time, using Table 1 to judge the sector and selecting a zero vector and two active voltage vectors in this sector as the vector Port1 is operated in U dc Q mode to ensure the stability of the DC-side voltage. In the traditional MPC method, the PI controller is used in the outer-loop subsystem to output d-axis reference current i d1ref , which can be expressed as: (29) where k p and k i are the proportional and integral coefficients of PI controller, respectively. u dcref is the reference voltage of the DC side. The PQ mode is used in Port2. With a stable grid voltage, a given reference current is equivalent to a given reference power, which can transform the PQ mode control into a problem of tracking the reference current.
Taking VSC1 for analysis, since two zero vectors V0 and V7 obtain the same voltage vector, in traditional MPC, seven voltage vectors are substituted into (3) to calculate the corresponding predicted currents, and each predicted current is substituted into the cost function, the voltage vector with the smallest cost function value can be selected as the optimal voltage vector and applied to the switch in the next sampling time. The cost function can be designed as: where i q1ref is the q-axis current reference value of VSC1. Similarly, VSC2 with given reference current, cost function values can be calculated and the optimal voltage vector will be applied to the switch in the next sampling time.

TV-MPC Method
In the traditional MPC method, only one voltage vector can be used in a sampling time and seven times of cost function calculations are required at both port, which leads to large current harmonics and control delays. Therefore, using the principle of deadbeat control, a low-complexity TV-MPC method is proposed. By calculating the predicted value of the reference voltage of VSC at the next sampling time, using Table 1 to judge the sector and selecting a zero vector and two active voltage vectors in this sector as the vector combination. The vector combination is applied to the VSC at the next sampling time. Taking VSC1 as the example, in the principle of deadbeat control, assuming the predicted currents at (k + 1)th instant as follows: Substituting (31) into (3), the reference voltages of VSC1 u 1Ndref and u 1Nqref can be written as: According to u 1Ndref and u 1Nqref , the sector can be determined. The optimal vector combination u 1opt , u 2opt and u 0opt include two active voltage vectors and one zero vector. For the selected optimal vector combination, the action time of each vector in a sampling time needs to be calculated. The principle of modulation MPC is applied, assuming the voltage vector action time is inversely proportional to the cost function value. The action time of voltage vectors t 1 , t 2 and t 0 can be calculated as: where f 1 , f 2 and f 0 are the values of cost function (30) corresponding to u 1opt , u 2opt and u 0opt , respectively. In this section, i d1ref is calculated by STC (24).
VSC2 is similar to VSC1, the optimal voltage vectors with the corresponding action times are output to the switches at the next sampling time, which will control the whole SOP system.

TV-MPC with STO
After adding the disturbance observed by STO, (32) can be rewritten as follows: The new voltage combination can be selected by the reference voltage u d1ref and u q1ref . By the new voltage combination, the values of cost function and action time of each vector can be recalculated. After adding the disturbance, (3) can be rewritten as: After substituting the rewritten predicted currents of (k + 1)th instant by new vector combinations chosen by (34) into (30) the values of the cost function can be calculated. Action times of each voltage vectors t 1 , t 2 and t 0 can be calculated by (33).
VSC2 is similar to VSC1, new voltage combinations and action times are output to each side of switch and control the whole system of SOP. Figure 5 is the block diagram of the STC and TV-MPC with STO for SOP.
The new voltage combination can be selected by the reference voltage ud1ref and uq1ref. By the new voltage combination, the values of cost function and action time of each vector can be recalculated.
After adding the disturbance, (3) can be rewritten as: After substituting the rewritten predicted currents of (k + 1)th instant by new vector combinations chosen by (34) into (30) the values of the cost function can be calculated. Action times of each voltage vectors t1, t2 and t0 can be calculated by (33).
VSC2 is similar to VSC1, new voltage combinations and action times are output to each side of switch and control the whole system of SOP. Figure 5 is the block diagram of the STC and TV-MPC with STO for SOP.

Simulations
In this paper, traditional MPC, TV-MPC and TV-MPC with STC were designed to demonstrate the effectiveness of the proposed method and were simulated in MATLAB/Simulation. TV-MPC and STO with STC was designed when parameter mismatches occurred in the SOP system. Table 2 lists the parameters of the SOP system.

Simulations
In this paper, traditional MPC, TV-MPC and TV-MPC with STC were designed to demonstrate the effectiveness of the proposed method and were simulated in MAT-LAB/Simulation. TV-MPC and STO with STC was designed when parameter mismatches occurred in the SOP system. Table 2 lists the parameters of the SOP system.

Current Pefromance of Steady-State
The proposed current of VSC2 side is set as −40 A in simulations. Outer-loop controller parameters are listed in Table 3. Steady-state currents of PI and MPC, PI and TV-MPC and STC and TV-MPC are shown in Figure 6.
It can be seen that in VSC1 side, the total harmonic distortions (THDs) of a-phase of PI and MPC, PI and TV-MPC and STC and TV-MPC are 1.07%, 0.43% and 0.58%, respectively. While on VSC2 side, the THD of a-phase of PI and MPC, PI and TV-MPC and STC and TV-MPC are 1.06%, 0.43% and 0.44%, respectively. Compared with PI and MPC, it is obvious that PI and TV-MPC has a better steady-state current performance and reduces the current ripple of each side. Current performances of STC and TV-MPC are similar to PI and TV-MPC, maintaining the low current ripple of the system.

Current Pefromance of Steady-State
The proposed current of VSC2 side is set as −40 A in simulations. Outer-loop controller parameters are listed in Table 3. Steady-state currents of PI and MPC, PI and TV-MPC and STC and TV-MPC are shown in Figure 6.

Method
Outer-Loop Controller Parameters PI and MPC kp = 3.5, ki = 4.125 PI and TV-MPC kp = 3.5, ki = 4.125 STC and TV-MPC k1= 150, k2 = 3000 It can be seen that in VSC1 side, the total harmonic distortions (THDs) of a-phase of PI and MPC, PI and TV-MPC and STC and TV-MPC are 1.07%, 0.43% and 0.58%, respectively. While on VSC2 side, the THD of a-phase of PI and MPC, PI and TV-MPC and STC and TV-MPC are 1.06%, 0.43% and 0.44%, respectively. Compared with PI and MPC, it is obvious that PI and TV-MPC has a better steady-state current performance and reduces the current ripple of each side. Current performances of STC and TV-MPC are similar to PI and TV-MPC, maintaining the low current ripple of the system.

DC Voltage Performance
Regulating the proposed DC-side voltage u dcref at 850 V, in order to verify the effectiveness of STC+TV-MPC in start-up response performance and dynamic state performance, the reference current of VSC2 side is changed from 40 A to 80 A at 0.25 s. As is shown in Figure 7a, the time used to track the proposed DC voltage under PI and MPC, PI and TV-MPC and STC and TV-MPC are 0.108 s, 0.029 s and 0.016 s, respectively. It can be seen that STC and TV-MPC track the proposed DC voltage faster than PI and MPC and PI and TV-MPC without overshoot, which has a better start-up response performance.

DC Voltage Performance
Regulating the proposed DC-side voltage udcref at 850 V, in order to verify the effectiveness of STC+TV-MPC in start-up response performance and dynamic state performance, the reference current of VSC2 side is changed from 40 A to 80 A at 0.25 s. As is shown in Figure 7a, the time used to track the proposed DC voltage under PI and MPC, PI and TV-MPC and STC and TV-MPC are 0.108 s, 0.029 s and 0.016 s, respectively. It can be seen that STC and TV-MPC track the proposed DC voltage faster than PI and MPC and PI and TV-MPC without overshoot, which has a better start-up response performance.
Changes of the DC-side voltage in dynamic state are shown in Figure 7b. It can be seen that the time used to recover the proposed DC-side voltage and voltage drops of STC and TV-MPC is much less than PI and MPC and PI and TV-MPC.
As a result, it can be seen that STC and TV-MPC has better performance in start-up response and current change response.

STC and TV-MPC with STO under Parameter Mismatches
In order to prove the effectiveness of STC and TV-MPC with STO, the SOP system under parameter mismatches is simulated in this section. STO and STC parameters are listed in Table 4. Changes of the DC-side voltage in dynamic state are shown in Figure 7b. It can be seen that the time used to recover the proposed DC-side voltage and voltage drops of STC and TV-MPC is much less than PI and MPC and PI and TV-MPC.
As a result, it can be seen that STC and TV-MPC has better performance in start-up response and current change response.

STC and TV-MPC with STO under Parameter Mismatches
In order to prove the effectiveness of STC and TV-MPC with STO, the SOP system under parameter mismatches is simulated in this section. STO and STC parameters are listed in Table 4. The proposed current of VSC2 side is set as 40 A, when parameter mismatches have not occurred in the SOP system, Figure 8a shows the d-axis currents of SOP system i d1 , i d2 , i d1ref and i d2ref . Figure 8b shows the q-axis currents of SOP system i q1 , i q2 , i q1ref and i q2ref . The disturbances f d1 , f q1 , f d2 and f q2 estimated by STO are shown in Figure 8c. It can be seen that when parameter mismatches have not occurred in SOP system, i d1 , i q1 , i d2 and i q2 can track the reference currents correctly and the disturbances estimated by STO are small.  The proposed current of VSC2 side is set as 40 A, when parameter mismatches have not occurred in the SOP system, Figure 8a shows the d-axis currents of SOP system id1, id2, id1ref and id2ref. Figure 8b shows the q-axis currents of SOP system iq1, iq2, iq1ref and iq2ref. The disturbances fd1, fq1, fd2 and fq2 estimated by STO are shown in Figure 8c. It can be seen that when parameter mismatches have not occurred in SOP system, id1, iq1, id2 and iq2 can track the reference currents correctly and the disturbances estimated by STO are small. When inductance mismatches have occurred in the SOP system, STO can observe the disturbances caused by inductance mismatches. The SOP system starts at 0 s and ends at When inductance mismatches have occurred in the SOP system, STO can observe the disturbances caused by inductance mismatches. The SOP system starts at 0 s and ends at 0.8 s, while the observers are applied to the system at 0.4 s. When system inductances become 5L 1 and 5L 2 , as shown in Figure 9a, inductances' mismatches influence the d-axis currents obviously, Figure 9b shows the q-axis currents, when observers are applied to the system at 0.4 s, currents can track the reference currents correctly. Figure 9c shows the disturbances estimated by STO, it can be seen that after a short time, the estimated disturbances can be stabilized. When system inductances become 0.3L 1 and 0.3L 2 , Figure 10a-c present the d-axis current, q-axis current and disturbances estimated by STO, respectively, it can be seen the STO can also estimate the disturbances correctly. These results verify the effectiveness of proposed STO. 0.8 s, while the observers are applied to the system at 0.4 s. When system inductances become 5L1 and 5L2, as shown in Figure 9a, inductances' mismatches influence the d-axis currents obviously, Figure 9b shows the q-axis currents, when observers are applied to the system at 0.4 s, currents can track the reference currents correctly. Figure 9c shows the disturbances estimated by STO, it can be seen that after a short time, the estimated disturbances can be stabilized. When system inductances become 0.3L1 and 0.3L2, Figure 10ac present the d-axis current, q-axis current and disturbances estimated by STO, respectively, it can be seen the STO can also estimate the disturbances correctly. These results verify the effectiveness of proposed STO.
(a) (b) (c) Figure 9. (a) d-axis currents of SOP with 5L1 and 5L2; (b) q-axis currents of SOP with 5L1 and 5L2; (c) disturbances observed by STO with 5L1 and 5L2. Figure 9. (a) d-axis currents of SOP with 5L 1 and 5L 2 ; (b) q-axis currents of SOP with 5L 1 and 5L 2 ; (c) disturbances observed by STO with 5L 1 and 5L 2 .  Figure 11 shows the influences of resistance mismatches, when system resistances become 3R1 and 3R2, Figure 11a,b show the d-axis and q-axis currents of SOP system, it can be seen that resistance mismatches have less influence than inductance mismatches. As shown in Figure 11c, disturbances caused by resistance mismatches can be observed by STO. Therefore, it can be proved that the proposed STC and TV-MPC with STO can reduce the currents error obviously when parameter mismatches have occurred in the SOP system.  Figure 11 shows the influences of resistance mismatches, when system resistances become 3R 1 and 3R 2 , Figure 11a,b show the d-axis and q-axis currents of SOP system, it can be seen that resistance mismatches have less influence than inductance mismatches. As shown in Figure 11c, disturbances caused by resistance mismatches can be observed by STO. Therefore, it can be proved that the proposed STC and TV-MPC with STO can reduce the currents error obviously when parameter mismatches have occurred in the SOP system.

Discussion
For the better performance of the SOP system, increasing the number of ports, considering the faults that may occur in the system, taking Vf mode into consideration, and an experimental study of the SOP system will be the future work.

Conclusions
In this paper, a STC and TV-MPC with STO was proposed to improve the operation of the SOP system. Considering the disturbances caused by parameter mismatches and unmodelled dynamics, the STO was designed to observe the disturbances, and steadystate currents can track the reference currents correctly with the use of STO. In addition, aiming at the problems such as long start-up time, overshoot and poor robustness caused by traditional out-loop PI controller, an outer-loop STC was designed to improve the performance of the DC-side voltage. Moreover, in order to reduce the current harmonics, an Figure 11. (a) d-axis currents of SOP with 3R 1 and 3R 2 ; (b) q-axis currents of SOP with 3R 1 and 3R 2 ; (c) disturbances observed by STO with 3R 1 and 3R 2 .

Discussion
For the better performance of the SOP system, increasing the number of ports, considering the faults that may occur in the system, taking Vf mode into consideration, and an experimental study of the SOP system will be the future work.

Conclusions
In this paper, a STC and TV-MPC with STO was proposed to improve the operation of the SOP system. Considering the disturbances caused by parameter mismatches and unmodelled dynamics, the STO was designed to observe the disturbances, and steadystate currents can track the reference currents correctly with the use of STO. In addition, aiming at the problems such as long start-up time, overshoot and poor robustness caused by traditional out-loop PI controller, an outer-loop STC was designed to improve the performance of the DC-side voltage. Moreover, in order to reduce the current harmonics, an inner-loop TV-MPC method was proposed. The simulation results show that STC and TV-MPC with STO can reduce current harmonics effectively while the DC voltage has a better operation state, when parameter mismatches have occurred in the SOP system, STO can observe the disturbances correctly, which prove the effectiveness and correctness of the proposed method.
Author Contributions: Z.W.: conceptualization, algorithm innovation, methodology, and writingoriginal draft; H.Z.: data and formal analysis, investigation, software, simulation, and writingoriginal draft; H.S.: formal analysis, and writing-review and editing; All authors have read and agreed to the published version of the manuscript.