An Improved Quantum Differential Evolution Algorithm for Optimization and Control in Power Systems Including DGs
-
Abstract: Differential evolution algorithm (DE) has been proved to be an effective way for solving the optimal reactive power flow (ORPF) problem. As distributed generations (DGs) are introduced into the system, there is a certain impact on power flow and voltage of the power system, which affects the robustness and effectiveness of DE. On the basis of DE, aiming at its limitation of premature convergence and poor search ability, this paper discusses about how to improve it with quantum encoding and artificial bee colony (ABC) algorithm and proposes a hybrid algorithm, which is called improved quantum differential evolution algorithm (IQDE). The idea of quantum encoding increases the individual diversity while the accelerating evolution operation of the onlooker bees improves local search ability of DE. At the same time, the random search operation of the scout bees improves global search ability of DE. In the last, the effectiveness of IQDE is verified by simulations on the IEEE 14-bus system and 30-bus system including DGs. The experimental results show that with less convergence time and smaller population size, IQDE can obtain an even or better optimization effect compared with DE and can be applied to ORPF problem of power system including DGs.
-
Key words:
- Artificial bee colony /
- differential evolution /
- optimal reactive power flow /
- quantum /
- distributed generation (DG)
摘要: 差分进化算法(DE)已被证明为解决无功优化问题的有效方法.随着越来越多的分布式电源并网,对配电网潮流、电压均有一定改变,同时也影响了DE的鲁棒性和性能.本文在研究DE基础上,针对其收敛过早、局部搜索能力较差的缺陷,分析了量子计算思想和人工蜂群算法的优势,提出改进量子差分进化混合算法(IQDE).通过量子编码思想提高了种群个体的多样性,人工蜂群算法的观察蜂加速进化操作和侦查蜂随机搜索操作分别提高了算法的局部搜索和全局搜索性能.建立以有功网损最小为目标的数学模型,将IQDE算法和DE算法分别用于14节点和30节点标准数据集进行大量仿真实验.实验结果表明,IQDE算法用更少的收敛时间、更小的种群规模便可以获得与DE算法相同甚至更佳的优化效果,并且可以很好的应用于解决难分布式电源的配电网无功优化问题. -
1. Introduction
For permanent magnet synchronous motor (PMSM) drive system, the measurement of instantaneous stator currents is required for successful operation of the feedback control. Generally two phase current sensors are installed in three phase voltage source inverters (VSI). Nevertheless, sudden severe failure of phase current sensors would result in over-current malfunction of the drive system. And if there is no protection scheme in the gate-drive circuit, the failure would lead to irrecoverable fault of power semiconductors in VSI, which would cause degradation of motor drive performance. Additionally, some minor failures (such as gain drift and nonzero offset) of phase current sensors would lead to torque pulsation synchronizing with the inverter output frequency [1]. The larger offset and scaling error of phase current sensors would bring about the worse performance of torque regulation. Moreover, if the offset and gain drift are above certain level, it would cause over-current trip under high speed and heavy load conditions [2]. So it is necessary to consider fault tolerant operation of phase current sensor failure.
The current sensorless technology, regarded as fault tolerant one, has been developed in the past few decades. Its core lies in that the physical fault current sensor is replaced with virtual sensor (or current estimator). This technology has several advantages such as high reliability and low cost as well as space and weight savings owing to omitting physical current sensor. Moreover, it allows the drive system to work in hostile environment.
As far as the current sensorless technique is concerned, three estimation solutions have been reported in the literature. The first one is a DC-link current-based approach which restructures phase currents with the information of the DC-link current and switching states in VSI [3]. Although it is a mainstream method, its unavoidable drawbacks are exposed: the duration of an active switching state may be so short that the DC-link current cannot be measured on one hand, on the other hand, there are immeasurable regions in the output voltage hexagon where the DC-link current sampling and reconstruction are limited or impossible to do [4]. In addition, the DC-link sensed current remains sensitive to the narrow pulse and further deteriorates if the cable capacitance causes spurious oscillations in the DC-link waveform. In order to provide high-accuracy phase current reconstruction over a wide range of operating conditions with a low current waveform, over the past years, many kinds of methods of improved PWM modulation strategy have been proposed for the single DC-link current sensor technique [5]$-$[14]. Although many improved methods show reasonable phase current reconstruction performance, these methods suffer from complicated algorithms [15]. The second one is an analytical model-based approach. In [16], on the basis of the voltage and flux equations of induction motor (IM) drive, the phase current is estimated by using the synchronous reference frame variables under single phase current sensor condition. In [17], by the discrete voltage equations of PMSM drive, the phase currents are estimated. Although it is easier to implement than the first one, the method is not robust against the variation of system parameters. The third one is an adaptive observer-based approach. In [18], the phase current is reconfigured for IM drive using single phase current sensor, while in [19], the phase currents are reconfigured for PMSM drive without any phase current sensors. Compared with the first two solutions, the third solution has stronger robustness against the variation of system parameters [20], [21]. For PMSM drive system when only one phase current sensor is available, the remaining two phase currents estimation based on an adaptive observer must be studied, which is required to perform current feedback control. However, there is no literature on such strategy.
For PMSM drive system, model predictive torque control (MPTC) is an emerging control strategy [22]$-$[29]. Its main objective is to control instantaneous torque and stator flux with high accuracy and thus MPTC plays an important role to ensure the quality of the torque and speed control. MPTC adopts the principle of model predictive control (MPC) and can provide high dynamic performance and low stator current harmonics.
For conventional proportional-integral(PI)-based MPTC PMSM drive system, its speed regulator employs the algorithm of PI. In general, PI may perform well under certain operating condition, but it does not work properly and thus degrades dynamic performance under other operating conditions such as variation of system parameters and external disturbances. To improve the robustness of the speed regulator, some techniques have been proposed in recent years [30]$-$[34]. Except these techniques, a global fast terminal sliding mode (GFTSM) control is an effective and practical one [35], [36], which is based on sliding mode theory and employs the fast terminal sliding mode in both the reaching stage and sliding stage. By adding the nonlinear function to the sliding mode surface, the GFTSM controller can enable drive system not only to be superiorly robust against system uncertainties and external disturbances but also to have quick response as well as high control precision. Even so, studies on GFTSM speed regulator are very few. In this paper, we propose replacement of PI with GFTSM for MPTC PMSM drive system.
In this paper, by referring to the adaptive approach and integrating the GFTSM method, a new GFTSM-based MPTC strategy with the adaptive observer is put forward for the PMSM drive system with single phase current sensor. The proposed adaptive observer presents a satisfactory tracking performance of the remaining two phase currents in the presence of stator resistance change caused by the temperature variation. And the designed GFTSM controller enhances the speed regulator's robustness against parameter uncertainty and external disturbance. On the basis of the above foundation, the synthesized MPTC PMSM drive control system achieves a high performance.
This paper is organized as follows: Dynamic model of PMSM drive is presented in Section 2. Section 3 gives the adaptive observer and GFTSM speed regulator design as well as MPTC design. Experimental results and analysis are presented in Section 4. Section 5 contains the conclusions.
Notation 1: The following nomenclature is used throughout this paper:
$ \begin{align*} \begin{array}{lll} R_{\rm s}:& &\mbox {Nominal phase resistance} \ \\ \psi_{\rm m}:&&\mbox{The permanent magnet flux} \ \\ \psi_{\rm s}:&&\mbox{Stator flux linkage}\\ p:&&\mbox{Number of pole pairs}\\ V_{\rm dc}:&&\mbox{DC bus voltage}\\ \omega_{\rm r}:&&\mbox{Rotor actual mechanical speed}\\ T_{\rm l}:&&\mbox{Load torque}\\ T_{\rm e}:&&\mbox{Electromagnetic torque}\\ J:&&\mbox{Moment of inertia}\\ B_{\rm m}:&&\mbox{Viscous friction coefficient}\\ T_{\rm f}:&&\mbox{Coulomb friction torque}\\ \theta:&&\mbox{Rotor electrical angular position}\\ i:&&\mbox{Stator current}\\ u:&&\mbox{Stator voltage}\\ L:&&\mbox{Stator inductance}.\\ \end{array} \end{align*} $
Notation 2: The following symbol is used throughout this paper. $\bullet_{\rm d}$, $\bullet_{\rm q}$, $\bullet_{\alpha}$ and $\bullet_{\beta}$ are used to denote the $d$-axis, $q$-axis, $\alpha$-axis, and $\beta$-axis component of $\bullet$, respectively; $\bullet^{\ast}$ is used to denote the reference values of $\bullet$; $\hat{\bullet}$ is used to denote the estimate of $\bullet$; $\tilde{\bullet}$ is used to denote the parameter estimation error of $\bullet$; $\bullet^{k}$ and $\bullet^{k+1}$ are used to denote the instantaneous value at $k$th and ($k+1$)th of $\bullet$, respectively.
2. Dynamic Models of Three-phase PMSM Drive
As for three-phase PMSM drive, the models in rotor synchronous reference frame ($dq$-frame) and two-phase stationary reference frame ($\alpha\beta$-frame) are expressed as follows, respectively:
$ \begin{align} \label{eq:1} % eq:(1) \begin{cases} \dfrac{{d}i_{\rm d}}{{d}t}&=\dfrac{1}{L_{\rm d}} \left(u_{\rm d}-R_{\rm s}i_{\rm d}+p\omega_{\rm r}L_{\rm q}i_{\rm q}\right) \\[3mm] \dfrac{{d}i_{\rm q}}{{d}t}&=\dfrac{1}{L_{\rm q}}\left(u_{\rm q}-R_{\rm s}i_{\rm q}+p\omega_{\rm r} (L_{\rm d}i_{\rm d}+\psi_{\rm m})\right) \end{cases}\end{align} $
(1) $\begin{align} \begin{cases} \dfrac{{d}i_{\alpha}}{{d}t}&=\dfrac{1}{L_{\alpha}} \left(u_{\alpha}-R_{\rm s}i_{\alpha}+p\omega_{\rm r}\psi_{\rm m} \sin\theta\right) \\[3mm] \dfrac{{d}i_{\beta}}{{d}t}&=\dfrac{1}{L_{\beta}} \left(u_{\beta}-R_{\rm s}i_{\beta}-p\omega_{\rm r}\psi_{\rm m}\cos\theta\right) \end{cases} \end{align} $
(2) and the mechanical equation is expressed as
$ \begin{align} \dfrac{{d}\omega_{\rm r}}{{d}t}=\dfrac{1}{J}(T_{\rm e}-T_{\rm l}-B_{\rm m}\omega_{\rm r}-T_{\rm f}) \end{align} $
(3) where the electromagnetic torque $T_{\rm e}$ is expressed as
$\begin{align} T_{\rm e}=\frac{3{ p}}{2} \left[\psi_{\rm m}i_{\rm q}+(L_{\rm d}-L_{\rm q})i_{\rm d}i_{\rm q}\right]. \end{align} $
(4) 3. Design of GFTSM-based MPTC PMSM Drive System With Adaptive Observer
The objective of GFTSM-based MPTC using adaptive observer is that the PMSM drive system can work reliably and its speed and torque can be controlled not only to have satisfactory performance but also to be strongly robust against parameters variation and external disturbance. The schematic of the proposed control system is shown in Fig. 1. Our design task concentrates on adaptive observer, GFTSM speed regulator and MPTC as follows.
3.1 Adaptive Observer Design
The proposed adaptive observer is to estimate the remaining two phase currents and stator resistance when single phase current sensor is available. In the design process, assume the following conditions.
1) Only phase-$b$ current can be measured and the remaining two phase current sensors are not available.
2) Due to heating during operating of the motor, the stator resistance $R_{\rm s}$ is considered as a time-varying parameter.
3) There is no saturation in the magnetic circuit.
For surface-mounted PMSM drive, $L_{\rm d}=L_{\rm q}=L_{\alpha}= L_{\beta}$ $=$ $L$. The $\alpha$-axis in $\alpha\beta$-frame is oriented along phase-$a$ axis in three-phase stationary reference frame ($abc$-frame). The $abc$-axis stator currents in $abc$-frame can be obtained from the $\alpha\beta$-axis ones in $\alpha\beta$-frame by the following transformation matrix:
$ \left[ \begin{matrix} {{i}_{\rm{a}}} \\ {{i}_{\rm{b}}} \\ {{i}_{\rm{c}}} \\ \end{matrix} \right]=\left[ \begin{matrix} 1 & 0 \\ -\frac{1}{2} & \frac{\sqrt{3}}{2} \\ -\frac{1}{2} & -\frac{\sqrt{3}}{2} \\ \end{matrix} \right]\left[ \begin{matrix} {{i}_{\alpha }} \\ {{i}_{\beta }} \\ \end{matrix} \right] $
(5) where $i_{\rm a}$, $i_{\rm b}$, and $i_{\rm c}$ are $abc$-axis stator currents in $abc$-frame. From (5), the following equation can be given,
$ \begin{align} % eq:(6) i_{\rm b}=-\frac{1}{2} i_\alpha+\frac{\sqrt{3}}{2}i_\beta. \end{align} $
(6) Taking (2) into account, the time derivative of (6) is deduced as follows:
$ \begin{align} % eq:(7) \frac{{d}i_{\rm b}}{{d}t}& =\dfrac{\sqrt{3}}{2L}\left[u_\beta-R_{\rm s}\left(\frac{1}{ \sqrt{3}}i_\alpha+\frac{2}{ \sqrt{3}}i_{\rm b}\right)-p\omega_{\rm r}\psi_m\cos\theta\right]\notag\\ & \quad -\frac{1}{2L} (u_{\alpha}-R_{\rm s}i_{\alpha}+p\omega_{\rm r}\psi_{\rm m}\sin\theta)\notag\\ &=\frac{\sqrt{3}u_\beta-u_{\alpha}-2R_{\rm s}i_{\rm b}-p\omega_{\rm r}\psi_{\rm m}(\sqrt{3}\cos\theta+\sin\theta)} {2L}. \end{align} $
(7) The following adaptive observer is proposed in order to estimate phase-$b$ current,
$ \begin{align} % eq:(8) \frac{{d}\hat{i}_{\rm b}}{{d}t}& =\frac{\sqrt{3}}{ 2L}\left[u_\beta-\hat{R}_{\rm s}\left(\frac{1}{ \sqrt{3}}\hat{i}_\alpha+\frac{2}{ \sqrt{3}}i_{\rm b}\right)-p\omega_{\rm r}\psi_m\cos\theta\right]\notag\\ & \quad -\frac{1}{ 2L}\left(u_{\alpha}-\hat{R}_{\rm s}\hat{i}_{\alpha}+p\omega_{\rm r}\psi_{\rm m}\sin\theta\right)-k_1f(\tilde{i}_{\rm b})-k_2\tilde{i}_{\rm b}\notag\\ & =\frac{1}{2L}\left[\sqrt{3}u_\beta-u_{\alpha}-2\hat{R}_{\rm s}i_{\rm b}-p\omega_{\rm r}\psi_{\rm m}(\sqrt{3}\cos\theta + \sin\theta)\right]\notag\\ &\quad-k_1f(\tilde{i}_{\rm b})-k_2\tilde{i}_{\rm b} \end{align} $
(8) where $k_1f(\tilde{i}_{\rm b})$ and $k_2\tilde{i}_{\rm b}$ are correctors, and $k_1$ and $k_2$ are the positive observer gains, and $f(\cdot)$ denotes the nonlinear function of phase-$b$ current estimation error $\tilde{i}_{\rm b}$, which is defined as
$ \begin{align} % eq:(9) \tilde{i}_{\rm b}=\hat{i}_{\rm b}-i_{\rm b}. \end{align} $
(9) Define the following stator resistance estimation error,
$ \begin{align} % eq:(10) \tilde{R}_{\rm s}=\hat{R}_{\rm s}-R_{\rm s}. \end{align} $
(10) By subtracting (8) from (7), the dynamics equation of the phase-$b$ current estimation error is given as follows:
$ \begin{align} % eq:(11) \frac{{d}\tilde{i}_{\rm b}}{{d}t}=-\frac{1}{ L}\tilde{R}_{\rm s}i_{\rm b}-k_1f(\tilde{i}_{\rm b})-k_2\tilde{i}_{\rm b}. \end{align} $
(11) In order to determine the adaptive law of the stator resistance and the observer gains, construct the candidate Lyapunov function as
$ \begin{align} % eq:(12) V_1=\frac{1}{ 2}\left(\tilde{i}_{\rm b}^2+{1\over r}\tilde{R}_{\rm s}^2\right) \end{align} $
(12) where $r$ is constant positive scalar.
The time derivative of (12) is obtained as follows:
$ \begin{align} % eq:(13) \frac{{d}{V_1}}{{d}t}=-k_2\tilde{i}_{\rm b}^2-k_1f(\tilde{i}_{\rm b})\tilde{i}_{\rm b}+\tilde{R}_{\rm s}\left(\frac{1}{ r}\frac{{d}\tilde{R}_{\rm s}}{{d}t}-\frac{1}{ L}i_{\rm b}\tilde{i}_{\rm b}\right). \end{align} $
(13) If we define following equality,
$ \begin{align} % eq:(14) \frac{1}{r}\frac{{d}\tilde{R}_{\rm s}}{{d}t}-\frac{1}{ L}i_{\rm b}\tilde{i}_{\rm b}=0. \end{align} $
(14) Equation (13) can be rewritten as below:
$ \begin{align} % eq:(15) \frac{{d}{V_1}}{{d}t}=-k_2\tilde{i}_{\rm b}^2-k_1f(\tilde{i}_{\rm b})\tilde{i}_{\rm b}. \end{align} $
(15) To render $\dot V_{1}$ negative, we assume
$ \begin{align} % eq:(16) f(\tilde{i}_{\rm b})={\rm sign}(\tilde{i}_{\rm b}). \end{align} $
(16) As a result, the following inequality is satisfied
$ \frac{{d}{V_1}}{{d}t}<0. $
By Lyapunov stability theorem, dynamic system (11) is stable, which means that both $\tilde{i}_{\rm b}$ and $\tilde{R}_{\rm s}$ can converge to zero. Since the variation of the stator resistance in the observer time scale is negligible, i.e.,
$ \frac{{d}R_{\rm s}}{{d}t}\approx 0 $
then the following formula holds
$ \begin{align} % eq:(17) \frac{{d}\tilde{R}_{\rm s}}{{d}t}=\frac{{d}\hat{R}_{\rm s}}{{d}t}-\frac{{d}R_{\rm s}}{{d}t}\approx\frac{{d}\hat{R}_{\rm s}}{{d}t}. \end{align} $
(17) Therefore, from (14), the adaptive mechanism of the stator resistance is derived as follows:
$ \begin{align} % eq:(18) \hat{R}_{\rm s}={r\over L}\int(i_{\rm b}\tilde{i}_{\rm b}){ d}t. \end{align} $
(18) With the adaptive mechanism in (18), the estimation value of the stator resistance can converge to its real value.
In order to improve the estimation accuracy of the stator resistance and to ensure a null steady error, on the basis of PI strategy, (18) is modified as below:
$ \begin{align} % eq:(19) \hat{R}_{\rm s}=\frac{r}{ L}\left\{{K_{P(R_{\rm s})}[i_{\rm b}(\hat{i}_{\rm b}-i_{\rm b})]+K_{I(R_{\rm s})}\int[i_{\rm b}(\hat{i}_{\rm b}-i_{\rm b})]{d}t}\right\} \end{align} $
(19) where $K_{P(R_{\rm s})}$ and $K_{I(R_{\rm s})}$ are proportional and integral scalars, respectively.
By replacing $R_{\rm s}$ in (2) with $\hat{R}_{\rm s}$ in (19), the $\alpha\beta$-axis currents observers can be constructed as follows:
$ \begin{align} % eq:(20) \begin{cases} \dfrac{{d}\hat{i}_{\alpha}}{{d}t}=\dfrac{1}{ L} \left(u_{\alpha}-\hat{R}_{\rm s}\hat{i}_{\alpha}+p\omega_{\rm r} \psi_{\rm m}\sin\theta\right) \\[3mm] \dfrac{{d}\hat{i}_{\beta}}{{d}t} =\dfrac{1}{ L}\left(u_{\beta}-\hat{R}_{\rm s}\hat{i}_{\beta} -p\omega_{\rm r}\psi_{\rm m}\cos\theta\right). \end{cases} \end{align} $
(20) By combining (8), (19) and (20), the block diagram of the designed adaptive observer is established as shown in Fig. 2, which treats the stator voltages, rotor electrical position and speed as the inputs, the $dq$-axis currents and stator resistance as outputs when only phase-$b$ current is measured.
Remark 1: From Fig. 2, it can be seen that estimating the phase-$b$ current is a key step and primary premise in construction of the adaptive observer. The error between the phase-$b$ measured current and its estimated value must be guaranteed to converge towards zero.
Remark 2: From (8) and (19), it can be seen that although the coupling relationship between $\hat{i}_{\rm b}$ and $\hat{R}_{\rm s}$ exists, we do not need to decouple them in the design process. In fact, the phase-$b$ current estimation (8) and the stator resistance adaptive law (19) are implemented and solved all together.
Remark 3: The convergence rate of the observer is dependent on the observer gains $k_1$ and $k_2$, which should be chosen to be large enough such that the observer responds as soon as possible.
Remark 4: The estimated $dq$-axis currents in Fig. 2 will be applied to MPTC as shown in Fig. 1.
Remark 5: From (5), the estimation of phase-$a$ current in $abc$-frame is equal to that of $\alpha$-axis current in $\alpha\beta$-frame as follows:
$ \begin{align} % eq:(21) \hat{i}_{\rm a}=\hat{i}_{\alpha}. \end{align} $
(21) Accordingly, the estimation of phase-$c$ current in $abc$-frame can be obtained as follows:
$ \hat{i}_{\rm c}=-(i_{\rm b}+\hat{i}_{\alpha}). $
Remark 6: The proposed adaptive observer is robust against only the stator resistance change. If other parameter uncertainties (such as stator inductance change and permanent magnet flux change, etc.) and unmodeled dynamics are required to be considered, then adaptive robust method with extended state observer can be borrowed from [20] and [33], which is our next research topic.
3.2 GFTSM Speed Regulator Design
3.2.1 GFTSM Design
Define the speed error as
$ e=\omega_{\rm r}^{\ast}-\omega_{\rm r}. $
Let
$ \begin{align} % eq:(22) x_1=e, ~~x_2=\dot{x}_{1}, ~~ u=\dot{T}_{\rm e}. \end{align} $
(22) Assume that $\omega_{\rm r}^{\ast}$ (or $\dot{\omega}_{\rm r}^{\ast}$), $T_{\rm l}$, $T_{\rm f}$ are constants and $\omega_{\rm r}$ has continuous second-order derivative. Then, the state equation of (3) can be expressed as following:
$ \begin{align} % eq:(23) \begin{cases} \dot{x}_{1}=x_2\\[1mm] \dot{x}_{2}=-\dfrac{B_{\rm m}}{ J}x_2-\dfrac{1}{ J}u \end{cases} \end{align} $
(23) where $u$ can be regarded as the control input.
Our target is to enable the drive system to be strongly robust and to have very fast response. For this reason, based on sliding mode theory, GFTSM speed regulator is employed. Fast terminal sliding mode surface is designed as following:
$ \begin{align} % eq:(24) s=\dot{x}_{1}+\alpha x_{1}+\beta x_{1}^\frac{q}{p} \end{align} $
(24) where $\alpha$, $\beta>0$; $q$, $p$ $(q<p)$ are positive odd integers.
Taking the first-order derivative of (24) yields
$ \begin{align} % eq:(25) \dot{s}=\left(\alpha-\frac{B_{\rm m}}{ J}\right)x_2-{\frac{1}{ J}}u+\beta\frac{{d}}{{d}t}\left(x_{1}^\frac{q}{p}\right). \end{align} $
(25) To make the system (23) reach the sliding mode surface in finite time, the fast terminal attractor is adopted as follows:
$ \begin{align} % eq:(26) \dot{s}=-\varphi s-\gamma s^\frac{v}{m} \end{align} $
(26) where $\varphi>0$, $\gamma>0$, $m>0$, $v>0$; $m$ and $v$ are odd integers.
Let (25) be equal to (26) and thus the following sliding mode control law can be obtained
$ \begin{align} % eq:(27) u=J\left(\left(\alpha-\frac{B_{\rm m}}{ J}\right)x_2+\beta\frac{{d}}{{d}t}\left(x_{1}^\frac{q}{p}\right)+\varphi s+\gamma s^\frac{v}{m}\right). \end{align} $
(27) By combining (23), (24) and (27), the block diagram of the designed GFTSM speed regulator is shown as in Fig. 3.
By solving differential equation (26), the time from any state $s(0)\neq 0$ to the sliding mode surface $s(t_f)$ can be derived as follows:
$ \begin{align} % eq:(28) t_f=\frac{m}{ \varphi(m-v)}\ln\frac{\varphi\left(s(0)\right)^\frac{m-v}{m}+\gamma}{ \gamma}. \end{align} $
(28) Remark 7: From (27), it can be seen that the sliding mode control law does not include switching item and thus weakens system chatter.
Remark 8: Under control law (27), one can easily see that if it converges to zero according to the terminal attractor (26), $x_1$ will accordingly converge to zero in terms of the following fast terminal attractor
$ \begin{align} % eq:(29) \dot{x}_{1}=-\alpha x_{1}-\beta x_{1}^\frac{q}{p}. \end{align} $
(29) It can be observed from (26) and (29) that the fast terminal attractors are adopted both in the reaching phase and in sliding phase. Consequently, the designed regulator (27) is a global terminal sliding mode one which guarantees the finite time control performance.
Remark 9: According to (28), $t_f$ can be set arbitrarily by adjusting parameters $m$, $v$, $\varphi$, $\gamma$.
Remark 10: The designed GFTSM speed regulator (27) is not only stable but also robust, which will be analyzed as below.
3.2.2 Stability Analysis
Construct Lyapunov function as
$ \begin{align} % eq:(30) V_2=\frac{1}{2}s^{2}. \end{align} $
(30) Differentiating (30) yields
$ \dot{V}_2=s\dot{s}=-\varphi s^2-\gamma s^\frac{m+v}{m} $
since $(m+v)$ is an even number, therefore $\dot{V}=s\dot{s}<0$. According to Lyapunov stability theory, the system (23) is stable and its movement can tend to sliding mode surface and finally reach the sliding mode.
3.2.3 Robustness Analysis
Considering parameter uncertainties and external disturbances, the system (23) is rewritten as following:
$\begin{align} % eq:(31) \begin{cases} \dot{x}_{1}=x_2\\[1mm] \dot{x}_{2}=-\dfrac{B_{\rm m}}{ J}x_2-\dfrac{1}{ J}u+d(x_1, x_2) \end{cases} \end{align} $
(31) where $d(x_1, x_2)$ can be regarded as the total disturbance including uncertainties and external disturbances. Assume $| d(x_1, x_2)|$ $\leq$ $D$, $D$ is maximum value.
As for system (31), differentiating (24) yields
$ \begin{align} % eq:(32) \dot{s}=\left(\alpha-\frac{B_{\rm m}}{ J}\right)x_2-\frac{1}{ J}u+d(x_1, x_2)+\beta\frac{{d}}{{d}t}\left(x_{1}^\frac{q}{p}\right). \end{align} $
(32) Substituting (27) into (32) yields
$ \begin{align} % eq:(33) \dot{s}& =-\varphi s-\gamma s^\frac{v}{m}+d(x_1, x_2)\notag\\[1mm] & =-\varphi s-\left(\gamma-\frac{d(x_1, x_2)} {s^\frac{v}{m}}\right)s^\frac{v}{m}. \end{align} $
(33) Let
$ \begin{align} % eq:(34) \bar{\gamma}=\gamma-\frac{d(x_1, x_2)}{s^\frac{v}{m}} \end{align} $
(34) then (33) can be rewritten as
$ \begin{align} % eq:(35) \dot{s}=-\varphi s-\bar{\gamma} s^\frac{v}{m}. \end{align} $
(35) To make (35) be a fast terminal attractor, (34) must satisfy $\bar{\gamma}>0$. Therefore, the following inequality holds true
$ \gamma-\frac{d(x_1, x_2)}{s^\frac{v}{m}}>\gamma-\frac{| d(x_1, x_2)|}{| s^\frac{v}{m}|}>\gamma-\frac{D}{ | s^\frac{v}{m}|}>0 $
then we can deduce
$ \begin{align} % eq:(36) \gamma>\frac{D}{ | s^\frac{v}{m}|}. \end{align} $
(36) Equation (36) is equivalent to
$ \begin{align} % eq:(37) | s|>\left(\frac{D}{ \gamma }\right)^\frac{m}{v}. \end{align} $
(37) As a result, the fast terminal convergence region $\Delta$ is constrained by
$ \begin{align} % eq:(38) \Delta=\left\{x_1, x_2:| s|\leq \left(\frac{D}{ \gamma }\right)^\frac{m}{v}\right\}. \end{align} $
(38) Furthermore, we assume
$ \begin{align} % eq:(39) \gamma={D\over | s^\frac{v}{m}|}+\eta, ~~~\eta>0. \end{align} $
(39) According to (35), the time from any state $s(0)\neq 0$ to the sliding surface is deduced as follows:
$ \begin{align} % eq:(40) \bar{t}_f={m\over \varphi(m-v)}\ln{\varphi\left(s(0)\right)^\frac{m-v}{m}+\bar{\gamma}\over \bar{\gamma}}. \end{align} $
(40) Since $\bar{\gamma}>\eta$, the following inequality can be deduced
$ \ln\frac{\varphi\left(s(0)\right)^\frac{m-v}{m}+\bar{\gamma}}{ \bar{\gamma}}\leq \ln\frac{\varphi\left(s(0)\right)^\frac{m-v}{m}+\eta}{ \eta} $
and then the reaching time satisfies
$ \begin{align} % eq:(41) \bar{t}_f\leq\frac{m}{ \varphi(m-v)}\ln\frac{\varphi\left(s(0)\right)^\frac{m-v}{m}+\eta}{ \eta}. \end{align} $
(41) Through the above analysis, it can be seen that if the condition $\bar{\gamma}>0$ holds then fast terminal convergence can be guaranteed and system (31) can reach neighborhood $\Delta$ of the sliding mode surface $s(\bar{t}_f)=0$ in finite time $\bar{t}_f$.
3.3 Model Predictive Torque Control
The basic idea of MPTC is to predict the future behavior of the variables over a time frame based on the model of the system. As shown in Fig. 1, MPTC includes three parts: cost function minimization, predictive model and flux and torque estimator.
3.3.1 Cost Function Minimization
For MPTC, the cost function is chosen such that both torque and flux at the end of the cycle is as close as possible to the reference value. Generally, the minimum value of cost function is defined as
$ \begin{align} % eq:(42) &\min g =\left| T_{\rm e}^{\ast}-T_{\rm e}^{k+1}\right|+k_3\left| \left| \psi_{\rm s}^{\ast}\right|-| \psi_{\rm s}^{k+1}| \right|\notag \\ &\, {\rm s.t.}\quad u_{\rm s}^{k}\in\{V_1, V_2, \ldots, V_6\} \end{align} $
(42) where $V_1$, $V_2$, $V_3$, $V_4$, $V_5$, and $V_6$ are six nonzero voltage space vectors and can be generated by three phase VSI with respect to the different switches states. A set of voltage space vectors $u_{\rm s}^{k}$ at $k$th instant is defined as
$ \begin{align} % eq:(43) u_{\rm s}^{k}=\frac{2V_{\rm dc}\left[S_{\rm a}^{k}+e^\frac{i2\pi}{3}S_{\rm b}^{k}+(e^\frac{i2\pi}{3})^{2}S_{\rm c}^{k}\right]}{3} \end{align} $
(43) where $S_{\rm a}^{k}~(x=a, b, c)$ at $k$th instant is upper power switch state of one of three legs. $S_{\rm a}^{k}=1$ or $S_{\rm a}^{k}=0$ when upper power switch of one leg is on or off. $k_3$ is the weighting factor.
In order to compensate inherent one-step delay which exists in practical digital system, the cost function (42) is revised as below:
$ \begin{align} % eq:(44) &\min g =\left| T_{\rm e}^{\ast}-T_{\rm e}^{k+2}\right|+k_3\left| \left| \psi_{\rm s}^{\ast}\right|-| \psi_{\rm s}^{k+2}| \right|\notag \\ &\, {\rm s.t.}\quad u_{\rm s}^{k}\in\{V_1, V_2, \ldots, V_6\}. \end{align} $
(44) 3.3.2 Predictive Model for Stator Currents
According to (1), the prediction of the stator current at the next sampling instant is expressed as
$ \begin{align} % eq:(45) \begin{cases} i_{\rm d}^{k+1}=i_{\rm d}^{k}+\dfrac{1}{L}\left(u_{\rm d}^{k}-R_{\rm s}i_{\rm d}^{k}+p\omega_{\rm r}^{k}Li_{\rm q}^{k}\right)T_{\rm s}\\[3mm] i_{\rm q}^{k+1}=i_{\rm q}^{k}+\dfrac{1}{ L}\left (u_{\rm q}^{k}-R_{\rm s}i_{\rm q}^{k}-p\omega_{\rm r}^{k} (Li_{\rm d}^{k}+\psi_{\rm m})\right)T_{\rm s} \end{cases} \end{align} $
(45) where $i_{\rm d}^{k}$, $i_{\rm q}^{k}$ and $R_{\rm s}$ are replaced by the corresponding estimated values coming from the observer in Fig. 2. $T_{\rm s}$ is the sampling period.
3.3.3 Torque and Flux Estimators
In $dq$-frame, the current-based flux-linkage can be expressed as following vector:
$ \begin{align} \left[% eq:(46) \begin{array}{c} \psi_{\rm d}^{k+1} \\ \psi_{\rm q}^{k+1} \\ \end{array} \right]=\left[ \begin{array}{cc} L&0 \\ 0&L \\ \end{array} \right]\left[ \begin{array}{c} i_{\rm d}^{k+1} \\ i_{\rm q}^{k+1}\\ \end{array} \right]+\left[ \begin{array}{c} \psi_{\rm m} \\ 0 \\ \end{array} \right]. \end{align} $
(46) The magnitude of stator flux linkage $\psi_{\rm s}$ is
$ \begin{align} % eq:(47) \psi_{\rm s}^{k+1}=\sqrt{(\psi_{\rm d}^{k+1})^2+(\psi_{\rm q}^{k+1})^2}. \end{align} $
(47) Electromagnetic torque developed in $dq$-frame can be estimated as following:
$ \begin{align} % eq:(48) T_{\rm e}^{k+1}=\frac{3}{ 2}{ p}\psi_{\rm m}i_{\rm q}^{k+1}. \end{align} $
(48) Substituting (45) into (48), the torque can be calculated.
4. Simulation Result and Analysis
In order to validate the effectiveness of proposed control strategy, the designed control system as shown in Fig. 1 has been implemented in MATLAB/Simulink/Simscape platform. The parameters of PMSM drive are given in Table Ⅰ. The sampling period is 100 $\mu$s, and value $k_3$ is selected to be 200. The reference stator flux $\psi_{\rm s}^{\ast}$ is 0.175 Wb. The parameters of the adaptive observer are
$ \begin{align*} &K_{P(R_{\rm s})}=0.006, ~~~K_{I(R_{\rm s})}=8\\ &k_1=30, ~~~k_2=5000, ~~~r=1000. \end{align*} $
表 Ⅰ PARAMETERS OF PMSM DRIVETable Ⅰ PARAMETERS OF PMSM DRIVESymbol Value Symbol Value $R_{\rm s}$ $2.875\, \Omega$ $\omega_{\rm r}^{\ast}$ 1000 rpm $L_{\rm d}, L_{\rm q}$ 0.0085 H $T_{\rm n}$ 4 N$\cdot$m $\psi_{\rm m}$ 0.175 Wb $J$ $0.0008\, {\rm Kg\cdot m}^2$ $p$ 4 $B_{\rm m}$ 0.001 N$\cdot$m$\cdot$s $V_{\rm dc}$ 300 V $T_{\rm f}$ 0 The parameters of GFTSM in Fig. 3 are determined as follows:
$ \begin{align*} &\alpha=100, ~~~\beta=250, ~~~p=7, ~~~q=5\\ &\varphi=1000, ~~~\gamma=80\, 000, ~~~ m=3, ~~~v=1.\end{align*} $
4.1 The GFTSM-based MPTC PMSM Drive System Comparison Between the One With Single Phase Current Sensor and the Other With Two Phase Current Sensors
In order to verify estimation accuracy of the observer for GFTSM-based MPTC PMSM drive system with single phase current sensor, two scenarios of numerical simulation are provided and compared, which correspond to PMSM system with two phase current sensors (phase-$a$ and -$b$ sensors) and PMSM system with single phase current sensor (phase-$b$), respectively. For convenience sake, the former scenario is marked as Case 1 and the latter one as Case 2. Except the above-mentioned different number of current sensors, the two systems employ completely identical GFTSM-based MPTC strategy.
Fig. 4 shows comparison of two scenarios in terms of stator currents, stator resistance, rotor speed and torque when the reference speed $n^{\ast}$ is set to 1000 pm, the load torque is increased from 0 N$\cdot$m to 4 N$\cdot$m at 0.1 seconds and the stator resistance is changed from its nominal value 2.875 $\Omega$ to 5 $\Omega$ at 0.3 seconds.
From Figs. 4(a)$-$4(c), it can be seen that, for designed adaptive observer of Case 2, its estimated $a$-axis and $c$-axis currents in $abc$-frame rapidly track corresponding ones of Case 1, and its estimated stator resistance can rapidly follow actual resistance change and converge to its actual value accurately. Figs. 4(d)$-$4(e) show that, for GFTSM-based MPTC system of Case 2, its speed and torque can be regulated in a satisfactory manner and it has almost as good performance as GFTSM-based MPTC system of Case 1.
4.2 The MPTC PMSM System Comparison Between the One Based on PI and the Other Based on GFTSM
For GFTSM-based MPTC PMSM systems, for the sake of verifying its stronger robustness, two systems are compared, which correspond to the PI-based and GFTSM-based MPTC PMSM systems, respectively. Except distinct outer-loop speed regulator (i.e., PI and GFTSM), the two systems employ completely identical MPTC and adaptive observer. In the simulation, their reference speeds $n^{\ast}$ are set to 1000 rpm, their load torques of 0 N$\cdot$m are increased to 4 N$\cdot$m at 0.1 seconds and stator resistance is at its nominal value 2.875 $\Omega$.
In the simulation, sampling values of three-phase currents are recorded over the time range from 0.1 seconds to 0.2 seconds. During this period, the fundamental frequency of three-phase currents is 66.67 Hz. Total harmonic distortion (THD) can be obtained by comparing the higher frequency components to the fundamental one.
4.2.1 The Comparison of Anti-load Variation Ability Under the Same Speed Transient Response
The parameters of PI for PI-based MPTC PMSM system are adjusted as follows:
$ K_P=0.7, ~~~K_I=0.03 $
such that PI-based MPTC system has almost the same speed transient response as GFTSM-based one.
Fig. 5 shows the dynamical responses in terms of speed, torque and stator currents. Fig. 5(a) intuitively gives the speed response comparison, which demonstrates that for GFTSM-based MPTC PMSM system, its speed can sharply adapt to the change of external load in a satisfactory manner, and its capability of accommodating the challenge of load disturbance is superior to PI-based one's. From Figs. 5(b)$-$5(d), it can be observed that for two systems with same adaptive observer, their torques, estimated $a$-axis and $c$-axis currents in $abc$-frame are almost the same.
Table Ⅱ shows THD comparison of three-phase currents. From Table Ⅱ, what can be observed is that the THD of the GFTSM-based MPTC is smaller than one of the PI-based MPTC.
表 Ⅱ THD OF THREE-PHASE STATORS' CURRENT(%)Table Ⅱ THD OF THREE-PHASE STATORS' CURRENT(%)Control scheme $i_{\rm a}$ $i_{\rm b}$ $i_{\rm c}$ PI-based MPTC $2.21$ $2.32$ $2.24$ GFTSM-based MPTC $1.84$ $1.88$ $1.85$ 4.2.2 The Comparison of Dynamic Responses Under the Same Speed Anti-load Variation Ability
The parameters of PI for PI-based MPTC PMSM system are adjusted as follows:
$ K_P=3, ~~~K_I=0.1 $
such that PI-based MPTC system has almost the same anti-load variation ability as GFTSM-based one.
Figs. 6(a)$-$6(d) show the dynamical responses in terms of speed, torque and stator currents. Fig. 6(a) intuitively gives their speed response comparison, which indicates that GFTSM-based MPTC PMSM system has smaller overshoot and faster settling time than PI-based one. Meanwhile, it can be found from Fig. 6(b) that the torque response of GFTSM-based MPTC PMSM system is better than one of PI-based. From Figs. 6(c)$-$6(d), it can be observed that, their estimated $a$-axis and $c$-axis currents in $abc$-frame are almost the same.
4.3 The MPTC PMSM System Comparison Between the One Based on SM and the Other Based on GFTSM
Here, the working condition of PMSM drive system is identical with Section 4.2.
For SM-based speed regulator, its sliding mode surface and its reaching law are selected as following:
$ s=ce+\dot{e} $
(49) $ \dot{s}=-k_4s-\varepsilon {\rm sign}(s) $
(50) 4.3.1 The Comparison of Anti-load Variation Ability Under the Same Speed Transient Response
The parameters of SM for SM-based MPTC PMSM system are adjusted as follows:
$ c=160, ~~~k_4=800, ~~~\varepsilon=3\times 10^5 $
such that SM-based MPTC system has almost the same speed transient response as GFTSM-based one.
Figs. 7(a)$-$7(d) show the dynamical responses in terms of speed, torque and stator currents. Fig. 7(a) illustrates that for GFTSM-based MPTC PMSM system, benefiting from the fast terminal sliding mode employed in both the reaching stage and the sliding stage, its recovery rate of speed response is obviously faster than SM-based one. From Figs. 7(b)$-$7(d), it can be seen that for two systems with same adaptive observer, their torques, estimated $a$-axis and $c$-axis currents in $abc$-frame are almost the same.
Table Ⅲ shows THD comparison of three-phase currents. From Table Ⅲ, what can be observed is that the THD of the GFTSM-based MPTC is smaller than one of the SM-based MPTC.
表 Ⅲ THD OF THREE-PHASE STATORS' CURRENT(%)Table Ⅲ THD OF THREE-PHASE STATORS' CURRENT(%)Control scheme $i_{\rm a}$ $i_{\rm b}$ $i_{\rm c}$ SM-based MPTC $2.01$ $2.12$ $2.14$ GFTSM-based MPTC $1.84$ $1.88$ $1.85$ 4.3.2 The Comparison of Dynamic Response Under the Same Speed Anti-load Variation Ability
The parameters of SM for SM-based MPTC PMSM system are adjusted as follows:
$ c=140, ~~k_4=2500, ~~\varepsilon=3\times 10^7 $
such that SM-based MPTC system has almost the same anti-load variation ability as GFTSM-based one.
Figs. 8(a)$-$8(d) show the dynamical responses in terms of speed, torque and stator currents. Fig. 8(a) shows that the speed dynamic performance is better than SM-based one. And it can be found from Figs. 8(b)$-$8(d) that for SM-based MPTC PMSM system, due to a switching function sign$(\cdot)$ in (50), therefore its torque, estimated $a$-axis and $c$-axis currents have significantly heavy chatter. On the other hand, for GFTSM-based one, its sliding reaching law in (26) is a continuous and smooth function, so the system chatter can be greatly reduced.
Summarizing above simulation experiments, we can obtain following results,
1) The proposed adaptive observer can estimate the remaining two phase currents and stator resistance rapidly and accurately.
2) Compared with PI-based and SM-based MPTC PMSM drive systems, GFTSM-based one has better dynamical response behavior and stronger robustness as well as smaller THD index of three-phase stator current.
5. Conclusions
This paper has put forward a novel GFTSM-based MPTC strategy for PMSM drive system with only one phase current sensor. Firstly, an adaptive observer is designed, which is capable of concurrent online estimation of the remaining two phase currents and time-varying stator resistance rapidly and accurately. Secondly, GFTSM speed regulator is designed and its stability and convergence as well as robustness are analytically verified based on Lyapunov stability theory. Finally, the MPTC strategy is employed to reduce the torque and flux ripples. The proposed observer can be embedded into a fault resilient PMSM drive system. In case of a phase current sensor failure, the designed observer can be used as a virtual current sensor which is robust against variation of stator resistance. And the designed GFTSM controller can enhance speed regulator's robustness against variation of system parameters and external disturbance. The resultant GFTSM-based MPTC strategy can guarantee that PMSM drive system with single phase current sensor achieves not only fast response but also high-precision control performance as well as strong robustness.
Our future research topic is that considering both parameters uncertainties and unmodeled dynamics, we will employ adaptive robust method with extended state observer to reconstruct stator currents observer.
-
Table Ⅰ Number of Control Variables of IEEE 14-Bus System
Variable Number T 3 U 5 Q 2 SUM 10 Table Ⅱ Setting of Control Variables of IEEE 14-Bus System
Variable Min Max Step T 0.9 1.1 0.01 U 0.9 1.1 - Q 0 0.18 0.06 Table Ⅲ Constraints of The State Variables of IEEE 14-Bus
Node Min (MVar) Max (MVar) 1 0 10 2 -40 50 5 -40 40 8 -10 40 11 -6 24 13 -6 24 Table Ⅳ Statistics of Results for DE and IQDE (IEEE 14-Bus)
SP Algorithm Lossmin Lossmax Lossavg Timemin Timemax Timeavg 10 DE 12.3714 13.1578 12.4401 3.8977 5.3590 4.9267 IQDE 12.3712 12.5450 12.3858 4.9335 9.0550 7.2024 20 DE 12.3713 12.6364 12.3892 8.7948 9.9069 10.6250 IQDE 12.3712 12.4035 12.3761 7.8281 14.6796 11.2216 30 DE 12.3712 12.549 12.3876 12.569 15.7014 14.8099 IQDE 12.3712 12.3712 12.3712 12.1678 17.7190 14.7212 40 DE 12.3712 12.4463 12.3776 16.6632 21.4070 20.2090 IQDE 12.3712 12.3712 12.3712 15.2144 21.7840 17.3987 50 DE 12.3712 12.4319 12.3774 19.3238 26.4667 25.0056 IQDE 12.3712 12.3712 12.3712 19.0269 26.7894 22.9460 60 DE 12.3712 12.399 12.3754 25.2242 33.3440 30.3782 IQDE 12.3712 12.3712 12.3712 24.4536 33.4038 28.0544 Table Ⅴ Control Variable Setting and PLOSS Before and After Optimization for IEEE 14-Bus System
U1 U2 U3 U6 U8 T4 T5 T7 Q9 Q14 Ploss Before 1.06 1.045 1.01 1.07 1.09 0.978 0.969 0.932 18 18 13.393 After 1.1 1.0779 1.0465 1.1 1.1 1.06 0.9 1.03 18 6 12.3712 Table Ⅵ Number of Control Variables of IEEE 30-Bus System
Variable Number T 4 U 6 Q 2 SUM 12 Table Ⅶ Setting of Control Variables of IEEE 30-Bus System
Variable Minimum Maximum Step T 0.9 1.1 0.02 U 0.9 1.1 - Q9 0 0.2 0.05 Q24 0 0.04 0.01 Table Ⅷ Constraints of the State Variables of 30-Bus
Node Min (MVar) Max (MVar) 1 0 10 2 -40 50 5 -40 40 8 -10 40 11 -6 24 13 -6 24 Table Ⅸ Statistics of Results for DE and IQDE (IEEE 14-Bus)
Node Qout(MVar) Qlow Qup Pout (MW) Plow Pup 9 0.0137 -0.012 0.025 0.15 0.11 0.19 19 0.0554 -0.013 0.0689 0.25 0.1 0.45 24 0.01 -0.0069 0.02 0.09 0.05 0.15 26 0.0425 -0.15 0.0638 0.185 0.124 0.32 Table Ⅹ Statistics of Trial Results for IQDE in 30-Bus System Containing DGS
No. Ploss Time No. Ploss Time No. Ploss Time 1 16.2163 24.2024 11 16.2163 29.7159 21 16.2164 28.0426 2 16.2163 25.9639 12 16.2163 29.4342 22 16.2163 23.9344 3 16.2164 34.057 13 16.2163 27.2983 23 16.2166 33.3986 4 16.2163 28.0157 14 16.2163 29.5571 24 16.2163 28.6112 5 16.2163 30.1713 15 16.2166 34.6243 25 16.2164 23.7194 6 16.2163 26.3804 16 16.2163 30.1448 26 16.2163 24.2563 7 16.2165 31.7187 17 16.2164 25.1229 27 16.2164 32.3472 8 16.2163 29.4458 18 16.2163 25.7156 28 16.2163 24.8397 9 16.2163 25.2483 19 16.2164 33.2381 29 16.2164 28.5767 10 16.2165 33.7229 20 16.2163 25.9665 30 16.2166 33.8194 Table Ⅺ Statistics of Trial Results (IEEE 30-Bus)
Results Average Min Max Ploss 16.2164 16.2163 16.2166 Time 28.7097 23.7194 34.6243 Table Ⅻ Statistics of Trial Results for IQDE in 30-Bus System Containing DGS
U1 U2 U5 U8 U11 U13 T11 T12 T15 T36 Q10 Q24 P1 1.06 1.045 1.01 1.01 1.082 1.071 1.06 1.04 0.96 1.04 10 3 P2 1.1 1.0759 1.0431 1.0462 1.1 1.1 0.9 1 1.02 0.9 20 4 -
[1] X. Y. Yin, X. Yan, X, Liu, and C. L. Wang, "Reactive power optimization based on improved hybrid genetic algorithm, " J. Northeast Dianli Univ. , vol. 34, no. 3, pp. 48-53, Jun. 2014. http://en.cnki.com.cn/Article_en/CJFDTOTAL-CGCZ201101020.htm [2] T. Y. Xiang, Q. S. Zhou, F. P. Li, and Y. Wang, "Research on niche genetic algorithm for Reactive Power Optimization, " Proc. CSEE, vol. 25, no. 17, pp. 48-51, Sep. 2005. http://en.cnki.com.cn/Article_en/CJFDTOTAL-ZGDC200517010.htm [3] M. Varadarajan and K. S. Swarup, "Network loss minimization with voltage security using differential evolution, " Electr. Power Syst. Res. , vol. 78, no. 5, pp. 815-823, May2008. http://www.researchgate.net/publication/223544555_Network_loss_minimization_with_voltage_security_using_differential_evolution [4] M. Basu, "Optimal power flow with FACTS devices using differential evolution, " Int. J. Electr. Power Energy Syst. , vol. 30, no. 2, pp. 150-156, Jan. 2008. https://www.researchgate.net/publication/222751245_Optimal_power_flow_with_FACTS_devices_using_differential_evolution [5] Y. T. Liu, L. Ma, and J. J. Zhang, "GA/SA/TS hybrid algorithms for reactive power optimization, " Proc. the IEEE Power Engineering Society Summer Meeting, Seattle, WA, USA, 2000, 245-249. https://www.researchgate.net/publication/3862451_GASATS_hybrid_algorithms_for_reactive_power_optimization [6] B. Zhao, C. X. Guo, and Y. J. Cao, "A multiagent-based particle swarm optimization approach for optimal reactive power dispatch, " IEEE Trans. on Power Syst. , vol. 20, no. 2, pp. 1070-1078, 2005. https://www.researchgate.net/publication/3267300_A_multiagent-based_particle_swarm_optimization_approach_for_optimal_reactive_power_dispatch [7] G. F. Fang, H. X. Wang, and X. S. Huang, "An improved genetic algorithm for reactive power optimization, " Proc. the EPSA, vol. 15, no. 4, pp. 15-18, Aug. 2003. http://en.cnki.com.cn/Article_en/CJFDTOTAL-DLZD200304004.htm [8] S. F. Wang, Z. P. Wan, H. Fan, C. Y. Xiang, and Y. G. Huang, "Reactive power optimization model and its hybrid algorithm based on bilevel programming, " Power Syst. Technol. , vol. 29, no. 9, pp. 22-25, May. 2005. https://www.researchgate.net/publication/290751059_Reactive_power_optimization_modeland_its_hybrid_algorithm_based_on_bilevel_programming [9] W. Liu, X. L. Liang, and X. L. An, "Power system reactive power optimization based on BEMPSO, " Power Syst. Protect. Control, vol. 38, no. 7, pp. 16-21, Apr. 2010. http://en.cnki.com.cn/Article_en/CJFDTOTAL-JDQW201007006.htm [10] L. F. Zheng, J. Y. Chen, H. Lin, S. L. Le, and F. Chen, "Reactive power optimization based on quantum particle swarm optimization in electrical power system, " Central China Electr. Power, vol. 24, no. 2, pp. 16-19, 2011. http://en.cnki.com.cn/Article_en/CJFDTOTAL-HZDL201102005.htm [11] B. Li, "Research of distribution reactive power optimization based on modified differential evolution algorithm, " M. S. thesis, North China Electr. Power Univ. , Beijing, China, 2012. [12] D. Devaraj and J. P. Roselyn, "Genetic algorithm based reactive power dispatch for voltage stability improvement, " Int. J. Electr. Power Energy Syst. , vol. 32, no. 10, pp. 1151-1156, Dec. 2010. https://www.researchgate.net/publication/229092065_Genetic_algorithm_based_reactive_power_dispatch_for_voltage_stability_improvement [13] Z. C. Hu, X. F. Wang, and G. Taylor, "Stochastic optimal reactive power dispatch: formulation and solution method, " Int. J. Electr. Power Energy Syst. , vol. 32, no. 6, pp. 615-62, Jul. 2010. https://www.researchgate.net/publication/245214843_Stochastic_optimal_reactive_power_dispatch_Formulation_and_solution_method [14] T. Malakar and S. K. Goswami, "Active and reactive dispatch with minimum control movements, " Int. J. Electr. Power Energy Syst, vol. 44, no. 1, pp. 78-87, Jan. 2013. https://www.researchgate.net/publication/256970308_Active_and_reactive_dispatch_with_minimum_control_movements [15] R. Storn and K. Price, "Differential evolution-a simple and efficient heuristic for global optimization over continuous spaces, " J. Glob. Optim. , vol. 11, no. 4, pp. 341-359, Dec. 1997. [16] W. F. Gao, S. Y. Liu, and L. L. Huang, "A novel artificial bee colony algorithm based on modified search equation and orthogonal learning, " IEEE Trans. Cybernet. , vol. 43, no. 3, pp. 1011-1024, Jun. 2013. https://www.researchgate.net/publication/232533895_A_Novel_Artificial_Bee_Colony_Algorithm_Based_on_Modified_Search_Equation_and_Orthogonal_Learning [17] F. S. Abu-Mouti and M. E. El-Hawary, "Optimal distributed generation allocation and sizing in distribution systems via artificial bee colony algorithm, " IEEE Trans. Power Deliv. , vol. 26, no. 4, pp. 2090-2101, Oct. 2011. https://www.researchgate.net/publication/252062867_Optimal_Distributed_Generation_Allocation_and_Sizing_in_Distribution_Systems_via_Artificial_Bee_Colony_Algorithm [18] R. W. Keyes, "Quantum computing and digital computing, " IEEE Trans. Electron Devices, vol. 57, no. 8, pp. 2041, Aug. 2010. https://www.researchgate.net/publication/260512239_Quantum_Computing_and_Digital_Computing -