Next Article in Journal
Comparison of Novel Approaches to the Predictive Control of a DC-DC Boost Converter, Based on Heuristics
Next Article in Special Issue
Towards Fast Response, Reduced Processing and Balanced Load in Fog-Based Data-Driven Smart Grid
Previous Article in Journal
Experimental Evaluation of the Performance of a Three-Phase Five-Level Cascaded H-Bridge Inverter by Means FPGA-Based Control Board for Grid Connected Applications
Previous Article in Special Issue
An Artificial Intelligence Method for Energy Efficient Operation of Crude Distillation Units under Uncertain Feed Composition
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Adaptive-Observer-Based Data Driven Voltage Control in Islanded-Mode of Distributed Energy Resource Systems

1
School of IoT Engineering, Jiangnan University, Wuxi 214122, China
2
School of Electric Power Engineering, Nanjing Institute of Technology, Nanjing 211167, China
*
Author to whom correspondence should be addressed.
Energies 2018, 11(12), 3299; https://doi.org/10.3390/en11123299
Submission received: 8 November 2018 / Revised: 19 November 2018 / Accepted: 21 November 2018 / Published: 26 November 2018
(This article belongs to the Special Issue Energy Efficiency and Data-Driven Control)

Abstract

:
In this paper, an adaptive observer based data driven control scheme is proposed for the voltage control of dispatchable distributed energy resource (DER) systems which work in islanded operation. In the design procedure of the proposed control scheme, we utilize the novel transformation and linearization technique for the islanded DER system dynamics, which is proper for the proposed data driven control algorithm. Moreover, the pseudo partial derivative (PPD) parameter matrix can be estimated online by multiple adaptive observers. Then, the adaptive constrained controller is designed only based on the online identification results derived from the input/output (I/O) data of the controlled DER system. It is theoretically proven that all the signals in the closed-loop control system are uniformly ultimately bounded based on the Lyapunov stability analysis approach. In addition, the results of the simulation comparison are given to verify the voltage control effect of the proposed control scheme.

1. Introduction

Solar energy is clean, without pollution, and inexhaustible, so the application of solar energy has been a wide concern of the academic and industrial circles. In general, many remote areas and islands around the world are rich in solar energy. Therefore, it is conducive to build a micro-grid by forming a dispatchable photovoltaic (PV) distributed energy resource (DER) system. When the micro-grid works in islanded operation mode, the DER system can be used as the main power supply for micro-grid to provide voltage and frequency support which can sufficiently exert the advantage of the fast response of the DER system. However, one of the main challenges, is to keep the stability of the amplitude and frequency of micro-grid’s ac voltage under the unbalanced network conditions [1,2].
In general, islanded operation is not permitted under the present utility standards [3]. However, under the precondition that the DER system has the capacity to supply the total power of the local loads, the stand-alone islanded system can be applied [4,5,6]. In recent years, there has been much research on how to realize the fine control of the islanded DER system. In [6], the stability of voltage and frequency of the stand-alone DER system was maintained by adopting the proposed linear quadratic Gaussian controller which utilizes the technique of robust control. However, the response of the system output voltage was very slow, and the frequency of the output voltage was unstable. In [7], a novel unified controller was proposed for regulating the output voltage ofthe multibus micro-grid system. In addition, simulation by using the real-time digital simulator (RTDS) can verify the effectiveness of the proposed voltage control scheme. However, from the experimental waveform, it can be seen that the voltage waveform has great distortion and the output power of the system also has jitter. In [8], a control strategywas proposed for multiple DER systems which both operate in islanded mode. However, the systems obtained a favorable control effect only when the load had the characteristics of being linear and balanced. In order to reduce the influence of the load on the performance of the DER islanded system, the authors proposed a control method by adopting the techniques of feed-forward and decoupling compensation in [5]. However, the control strategy designed in this paper is based on PID control, and the output voltage of the system has a large static error when tracking the reference voltage. In [9], a control scheme based on the repetitive and deadbeat control was proposed for dispatchable electronically-interfaced DER islanded systems under the unbalanced load condition. However, the fluctuation of the current in the system seemed to have a great influence on the voltage, especially when the non-linear load was connected. Sliding mode control (SMC) has been successfully applied in power systems [10,11,12,13,14,15]. In [13], a fractional-order SMC strategy with black-start ability was designed to improve the quality of the output ac voltage of the host DER islanded systems. In [14], a modified terminal sliding-mode control (TSMC) with integral compensation method was proposed as a way of controlling the voltage-sourced converter (VSC) which drives varying loads. In [15], an adaptive robust total sliding-mode control (ARTSMC) scheme was proposed to stabilize the voltage and frequency of the micro-grid. Moreover, the proposed ARTSMC scheme has good control performance between the two operation modes (grid-connected mode and off-grid mode). However, most of the control strategies for DER systems studied at present are model-based control strategies and their control effects depend on the accuracy of the dynamic model.
With the development of technology, the performance of control system has become increasingly demanding. The excellent integrated control methods are getting more and more attention. Different from traditional control, the new methods require little model information and can be computed online easily. Therefore, control methods based on the data-driven approach were proposed in this background [16,17,18]. Model-free adaptive control (MFAC), one of the data-driven control (DDC) methods, has been studied by some scholars [17,19,20,21,22]. In [23], an improved MFAC algorithm was proposed for the coordination design of a wide-area power system stabiliser (WAPSS) to ensure the cooperation of the controllers under different operating conditions without tedious system modelling. In [24], the authors proposed a distributed MFAC scheme for multi-agent systems to solve consensus tracking problems only by using the input/output (I/O) data of each agent. Theoretical analysis, simulation, and practical application proved that the MFAC method has the advantages of strong robustness, small computational burden and easy implementation [19,20,21,25,26].
In this paper, we propose a novel data-driven MFAC method based on the multiple adaptive observer technique to accommodate the multi-input multi-output (MIMO) DER system. The design of the proposed method only depends on the I/O data of the MIMO controlled plant. Compared with other control methods, the MFAC method does not depend on the information from an accurate dynamic model, thus avoiding the control effect deterioration caused by an incomplete model. Moreover, the pseudo partial derivative (PPD) parameter matrix can be estimated online by the multiple adaptive observer. Then, the adaptive constrained controller is designed based on the online identification of multiple observer to solve the control input rate constraint problem of DER system. The simulation and comparison results are given to verify the superiority of the proposed method. The major contributions and innovations of this paper are displayed as follows.
  • The model-free control method is applied to the voltage control of the DER system, which reduces the dependence of the controller on model information.
  • The PPD parameter estimation algorithm is designed based on the multiple observer technique, which can guarantee the errors of observation converge to 0.
  • We utilize the novel transformation and linearization technique for the islanded DER system dynamics. The main advantage of this method is that the dynamics of the MIMO DER system are transformed and linearized properly for the proposed MFAC algorithm.
  • By constructing the Lyapunov function, it is theoretically proven that all the signals in the closed-loop control system are bounded.
The rest of this paper is organized as follows. In Section 2, the mathematical model of the DER islanded system is given. Then, we utilize the novel transformation and linearization technique for the islanded DER system dynamics, which is proper for the following proposed MFAC algorithm. In Section 3, the multiple observer-based data-driven model-free adaptive constrained control scheme and the proof of stability for control system are presented. In Section 4, simulation results and comparisons are given to demonstrate the effectiveness of the proposed MFAC scheme. Finally, we summarize the conclusions in Section 5.

2. Dynamics Transformation and Linearization of Islanded DER System

2.1. Mathematic Model of the Islanded DER System

As shown in Figure 1, the dynamics of the islanded DER system can be expressed as the space-phasor equations, as follows [13]:
C d U x d t = I I x
L s d I d t = U v U x R I
where U x , I , I x and U v are defined as the space phasor corresponding to U x a b c , I a b c , I x a b c and U v a b c , respectively.
In view of (1), (2) and Equation U v = ( U d c / 2 ) τ , which is related to the PWM modulating signals ( τ a b c ) of the VSC, we can obtain
d 2 U x d t 2 = U d c τ 2 L s C R L s d U x d t U x L s C 1 C d I x d t R I x L s C .
Express (3) in the case of its d q components. Then, Equation (3) can be split into (4) and (5), where w denotes the frequency of all the a b c -frame variables:
d 2 U x d d t 2 = 2 w d U x q d t R L s d U x d d t + w 2 1 L s C U x d + R w L s U x q 1 C d I x d d t + w C I x q R L s C I x d + U d c 2 L s C τ d
d 2 U x q d t 2 = 2 w d U x d d t R L s d U x q d t + w 2 1 L s C U x q R w L s U x d 1 C d I x q d t w C I x d R L s C I x q + U d c 2 L s C τ q .

2.2. Dynamics Transformation and Linearization Technique

This section presents the novel transformation and linearization method of the islanded DER system dynamics, which is proper for the proposed MFAC scheme.
First, Equations (4) and (5) can be rewritten as the following state-space form:
U ¨ = M U ˙ + Z τ + N U + P I ˙ + Q I
where U = ( U x d , U x q ) T , I ˙ = ( I ˙ x d , I ˙ x q ) T , I = ( I x q , I x d ) T , τ = ( τ d , τ q ) T = ( τ 1 , τ 2 ) T , M = R L s 2 w 2 w R L s , Z = U d c 2 L s C , N = ( w 2 1 L s C ) R w L s ( w 2 1 L s C ) R w L s , P = 1 C , Q = w C R L s C w C R L s C .
Then, the state-space form of DER system can be approximately transformed into the following discrete form
U ¨ ( k ) = M U ˙ ( k ) + Z τ ( k ) + N U ( k ) + P I ˙ ( k ) + Q I ( k ) .
Define y ( k ) = U ˙ ( k ) = [ y 1 ( k ) , y 2 ( k ) ] T . Thus, using the condition of y ˙ ( k ) = y ( k + 1 ) y ( k ) / T , the system discrete form (7) can be rewritten as
y ( k + 1 ) = ( 1 + M T ) y ( k ) + Z T τ ( k ) + N T U ( k ) + P T I ˙ ( k ) + Q T I ( k )
where T denotes the sampling time.
In view of the discrete system (8), with Δ y ( k ) , Δ τ ( k ) , Δ U ( k ) , Δ I ˙ ( k ) , Δ I ( k ) T 0 for each fixed k, there must exist Φ ( k ) , called the PPD matrix. By combining the condition Δ y ( k + 1 ) = y ( k + 1 ) y ( k ) = [ Δ y 1 ( k + 1 ) , Δ y 2 ( k + 1 ) ] T , we can transform (8) into the following dynamic linearization data model:
Δ y ( k + 1 ) = ϕ 1 ( k ) Δ y ( k ) + ϕ 2 ( k ) Δ τ ( k ) + ϕ 3 ( k ) Δ U ( k ) + ϕ 4 ( k ) Δ I ˙ ( k ) + ϕ 5 ( k ) Δ I ( k ) .
That is
Δ y ( k + 1 ) = Φ ( k ) Δ x ( k )
where Δ x ( k ) = Δ y ( k ) , Δ τ ( k ) , Δ U ( k ) , Δ I ˙ ( k ) , Δ I ( k ) T , and Φ ( k ) = ϕ 1 ( k ) , ϕ 2 ( k ) , ϕ 3 ( k ) , ϕ 4 ( k ) , ϕ 5 ( k ) = [ Φ 1 ( k ) , Φ 2 ( k ) ] T , ϕ 1 ( k ) , ϕ 2 ( k ) , ϕ 3 ( k ) , ϕ 4 ( k ) , ϕ 5 ( k ) R 2 × 2 , Φ 1 , Φ 2 R 10 × 1 ,
ϕ l ( k ) = ϕ l 11 ( k ) ϕ l 12 ( k ) ϕ l 21 ( k ) ϕ l 22 ( k ) , l = 1 , , 5 .
Φ ( k ) c and Δ x ( k ) ε . c and ε denote the positive constants.
The linearization model (10) can also be divided into the n MISO equations. The ith MISO dynamic linearization data model is
Δ y i ( k + 1 ) = Φ i ( k ) Δ x ( k ) = Δ x T ( k ) Φ i T ( k ) , i = 1 , 2 .

3. Multiple Observer-Based Model-Free Adaptive Control Design for DER System

In this section, we put forward a novel MFAC scheme based on the dynamic linearization data model of the DER system. The major contributions include: (1) an unknown PPD parameter estimation algorithm is designed based on the multiple observer technique; (2) the constrained MFAC scheme design procedures are given; and (3) the Lyapunov function is constructed to prove the stabilization of closed-loop systems.

3.1. Multiple Observer-Based PPD Parameter Estimation Algorithm

Based on the above section, we know that the (10) can be transformed into the two MISO models (11). For the ith (11), we propose the parameter identification observer which has the following structure:
y ^ i ( k + 1 ) = y ^ i ( k ) + Δ x T ( k ) Φ ^ i T ( k ) + K i e i ( k )
where e i ( k ) = y i ( k ) y ^ i ( k ) represents the estimation error of the output, Φ ^ i ( k ) denotes the estimation of the ith row parameter vector of PPD, and we choose the gain K i by calculating F i = 1 K i in the unit circle.
Therefore, by taking Equations (11) and (12) into account, the output dynamic estimation error can be obtained by
e i ( k + 1 ) = Δ x T ( k ) Φ ˜ i T ( k ) + F i e i ( k )
where Φ ˜ i ( k ) = Φ i ( k ) Φ ^ i ( k ) denotes the error of PPD parameter estimation. Then, the adaptive update algorithm for the estimated parameter vector Φ i ( k ) is given by
Φ ^ i T ( k + 1 ) = Φ ^ i T ( k ) + Δ x ( k ) H i ( k ) ( e i ( k + 1 ) F i e i ( k ) ) .
We define the gain H i ( k ) as follows:
H i ( k ) = 2 Δ x ( k ) 2 + ν i 1
in which ν i denotes a positive constant; therefore, H i ( k ) is positive definite for all k. Notice that, with the assumption that Δ x ( k ) ε , H i ( k ) can be lower bounded as
H i ( k ) 2 ε 2 + ν i = κ i > 0 .
In view of Equations (13) and (14) and with the condition of Φ i ( k + 1 ) Φ i ( k ) , the error dynamics are obtained as
e i ( k + 1 ) = Δ x T ( k ) Φ ˜ i T ( k ) + F i e i ( k ) Φ ˜ i ( k + 1 ) = Γ i Φ ˜ i ( k )
where Γ i is given by
Γ i = I i Δ x ( k ) H i ( k ) Δ x T ( k )
and I i represents the ( 2 × 2 ) identity matrix.
Theorem 1.
The equilibrium e i , Φ ˜ i T T = 0 , 0 2 × 1 T T of the system (17) is globally uniformly stable. Moreover, the estimated error e i ( k ) converges asymptotically to 0.
Proof. 
Construct the following Lyapunov function
V 1 ( k ) = G e i 2 ( k ) + λ Φ ˜ i T ( k ) Φ ˜ i ( k )
where λ and B denote positive constants, and G can be calculated by G F i 2 G = B . In view of (17), we obtain
V 1 ( k + 1 ) V 1 ( k ) = ( G F i 2 G ) e i 2 ( k ) + 2 G F i e i ( k ) Δ x T ( k ) Φ ˜ i T ( k ) Φ ˜ i ( k ) λ ( I i Γ i T Γ i ) Δ x ( k ) G Δ x T ( k ) Φ ˜ i T ( k ) = B e i 2 ( k ) + 2 G F i e i ( k ) θ i ( k ) θ i T ( k ) ( λ ν i H i 2 G ) θ i ( k ) B e i 2 ( k ) ν i λ κ i 2 G θ i 2 ( k ) + 2 G F i e i ( k ) θ i ( k ) a 1 e i 2 ( k ) a 2 θ i 2 ( k )
where θ i ( k ) = Δ x T ( k ) Φ ˜ i T ( k ) , a 1 = B ( 1 / ς ) , and a 2 = ν i λ κ i 2 G ς G 2 F i 2 . Therefore, Δ V 1 ( k + 1 ) 0 in the case of ς , B and λ satisfy the following inequalities:
B > 1 ς , ν i λ κ i 2 G ς G 2 F i 2 > 0 .
Notice that Δ V 1 ( k ) is negative definite with respect to variables e i ( k ) and θ i ( k ) . Since V 1 ( k ) is a decreasing and non-negative function, it converges to a constant value V 1 0 when k ; Therefore, Δ V 1 ( k ) 0 . This indicates that both e i ( k ) and Φ ˜ i ( k ) remain bounded for all k, and lim k e i ( k ) = 0 . ☐
Corollary 1: For model (10), the multiple adaptive observers
1 y ^ 1 ( k + 1 ) = y ^ 1 ( k ) + Δ x T ( k ) Φ ^ 1 T ( k ) + K 1 e 1 ( k ) Δ Φ ^ 1 T ( k + 1 ) = Δ x ( k ) H 1 ( k ) ( e 1 ( k + 1 ) F 1 e 1 ( k ) ) H 1 ( k ) = 2 Δ x ( k ) 2 + ν 1 1 2 y ^ 2 ( k + 1 ) = y ^ 2 ( k ) + Δ x T ( k ) Φ ^ 2 T ( k ) + K 2 e 2 ( k ) Δ Φ ^ 2 T ( k + 1 ) = Δ x ( k ) H 2 ( k ) ( e 2 ( k + 1 ) F 2 e 2 ( k ) ) H 2 ( k ) = 2 Δ x ( k ) 2 + ν 2 1
can be designed to ensure that lim k e ( k ) = 0 , where e ( k ) = y ( k ) y ^ ( k ) and
Δ Φ ^ i T ( k + 1 ) = Φ ^ i T ( k + 1 ) Φ ^ i T ( k ) , i = 1 , 2 .
Proof. 
Based on the analysis in the above section and Theorem 1, Corollary 1 is easily proved. ☐
In order to ascertain the subject for further elaboration, we joined the multiple observers (12) as
y ^ ( k + 1 ) = y ^ ( k ) + Φ ^ ( k ) Δ x ( k ) + Ke ( k )
where K = d i a g ( K 1 , K 2 ) .

3.2. Controller Design and Stability Analysis

In practice, the control input of the controlled system (6) is subject to the rate constraint due to the “inertia” of the actuator. Therefore, based on the multiobserver (24), we propose the constrained MFAC method to limit the change the rate of the control input ø ( k ) , which is computed as follows:
ø ( k ) = ø ( k 1 ) + ϕ ^ 2 T ( k ) ϕ ^ 2 ( k ) ϕ ^ 2 T ( k ) + σ 1 Δ ø ( k ) δ × [ y * ( k + 1 ) y ^ ( k ) Ke ( k ) ϕ ^ 1 ( k ) Δ y ( k ) ϕ ^ 3 ( k ) Δ U ( k ) ϕ ^ 4 ( k ) Δ I ˙ ( k ) ϕ ^ 5 ( k ) Δ I ( k ) ] , ø ( k ) = ø ( k 1 ) + δ s i g n ( Δ ø ( k ) ) , Δ ø ( k ) > δ
where y * ( k ) denotes the reference trajectory, that is, the voltage reference value. σ = diag ( σ 1 , σ 2 ) with σ 1 and σ 2 are given finite positive numbers to ensure that matrix ϕ ^ 2 ( k ) ϕ ^ 2 T ( k ) + σ is invertible. δ = [ δ 1 , δ 2 ] T with δ 1 and δ 2 also have finite positive numbers. Furthermore, the actuator cannot change too fast in many practical systems, so we can obtain the δ as the maximum change value for the control input ø ( k ) by ø ( k ) ø ( k 1 ) δ .
Define the tracking error of multiobserver as follows:
e 1 ( k ) = y * ( k ) y ^ ( k ) .
Then, we obtain
e 1 ( k + 1 ) = y * ( k + 1 ) y ^ ( k + 1 ) = y * ( k + 1 ) y ^ ( k ) Φ ^ ( k ) Δ x ( k ) Ke ( k ) = y * ( k + 1 ) y ^ ( k ) ϕ ^ 2 ( k ) Δ ø ( k ) J ( k ) Ke ( k )
where J ( k ) = ϕ ^ 1 ( k ) Δ y ( k ) + ϕ ^ 3 ( k ) Δ U ( k ) + ϕ ^ 4 ( k ) Δ I ˙ ( k ) + ϕ ^ 5 ( k ) Δ I ( k ) .
The analysis of convergence and tracking performance for the MFAC algorithm (25) are shown in Theorem 2.
Theorem 2.
Under the condition of y i * ( k ) y i * ( k 1 ) Δ y i * , by using the MFAC algorithm (25), the result of the observer close-loop error system (27) is uniformly ultimately bounded (UUB) for all k, that is, lim k e 1 ( k ) ( b 2 / ( 1 b 1 ) )
where Δ y i * is the given positive constant, and
b 1 = 1 D ¯ ( k ) + D ¯ ( k ) σ ( λ ¯ I m + σ ) 1 b 2 = b 1 Δ y * + b 1 Ke ( k ) + J ( k )
Δ y * = [ Δ y 1 * , Δ y 2 * ] T , D ¯ ( k ) = min { d 1 ( k ) , d 2 ( k ) } , λ ¯ represents the minimum eigenvalue of ϕ ^ 2 ( k ) ϕ ^ 2 T ( k ) .
Proof. 
The MFAC algorithm (25) can be equivalently described as
Δ ø ( k ) = D ( k ) ϕ ^ 2 T ( k ) ϕ ^ 2 ( k ) ϕ ^ 2 T ( k ) + σ 1 y * ( k + 1 ) y ^ ( k ) Ke ( k ) J ( k )
where
d j ( k ) = 1 , for | Δ τ j ( k ) | δ j 0 d j ( k ) < 1 , for | Δ τ j ( k ) | > δ j , j = 1 , 2 .
By substituting (29) into (27), we obtain (31) after rearranging
e 1 ( k + 1 ) = I m ϕ ^ 2 ( k ) D ( k ) ϕ ^ 2 T ( k ) ϕ ^ 2 ( k ) ϕ ^ 2 T ( k ) + σ 1 × y * ( k + 1 ) y ^ ( k ) Ke ( k ) J ( k ) = ϕ ^ 2 ( k ) I m D ( k ) ϕ ^ 2 T ( k ) + σ ϕ ^ 2 ( k ) ϕ ^ 2 T ( k ) + σ 1 × y * ( k + 1 ) y ^ ( k ) Ke ( k ) J ( k ) ( 1 D ¯ ( k ) ) I m + D ¯ ( k ) σ ( λ ¯ I m + σ ) 1 y * ( k + 1 ) y ^ ( k ) Ke ( k ) J ( k ) .
Then, (31) becomes
e 1 ( k + 1 ) 1 D ¯ ( k ) + D ¯ ( k ) σ ( λ ¯ I m + σ ) 1 y * ( k + 1 ) y ^ ( k ) Ke ( k ) J ( k ) = 1 D ¯ ( k ) + D ¯ ( k ) σ ( λ ¯ I m + σ ) 1 × y * ( k + 1 ) y * ( k ) + y * ( k ) y ^ ( k ) Ke ( k ) J ( k ) 1 D ¯ ( k ) + D ¯ ( k ) σ ( λ ¯ I m + σ ) 1 e 1 ( k ) + 1 D ¯ ( k ) + D ¯ ( k ) σ ( λ ¯ I m + σ ) 1 Δ y * + Ke ( k ) + J ( k ) = b 1 e 1 ( k ) + b 2 .
By constructing the Lyapunov function as V 2 ( k ) = e 1 ( k ) , from Equation (32), we obtain
Δ V 2 ( k + 1 ) = e 1 ( k + 1 ) e 1 ( k ) = ( 1 b 1 ) V 2 ( k ) + b 2 .
Since 0 b 1 < 1 and b 2 are bounded, in view of the MFAC algorithm (25), combined with the lemma in [27], we know that the solution of the observer close-loop error system (27) is UUB for all k with ultimate bound lim k e 1 ( k ) ( b 2 / ( 1 b 1 ) ) . ☐
Corollary 2: Based on the controller (25), combined with the parameter adaptive algorithm (14) and the observer (24), we can ensure that the system (6) tracking error e 2 ( k ) = y * ( k ) y ( k ) is UUB with ultimate bound lim k e 2 ( k ) ( b 2 / ( 1 b 1 ) ) .
Proof. 
Since
e 2 ( k ) = e 1 ( k ) e ( k ) ,
by considering the absolute value and limiting on both sides of Eqaution (34), we get
lim k e 2 ( k ) lim k e 1 ( k ) + lim k e ( k ) = b 2 1 b 1 .
Hence, the tracking error e 2 ( k ) = y * ( k ) y ( k ) is UUB with the ultimate bound lim k e 2 ( k ) ( b 2 / ( 1 b 1 ) ) . ☐

4. Simulation Results

In this section, two detailed types of load (a balanced linear load and an unbalanced linear load) are chosen for the islanded DER system to verify the effectiveness of the proposed voltage control scheme. Moreover, the superiority of the proposed control scheme is compared to the traditional PID control method [28] and fractional-order sliding-mode control (FSMC) [13]. Since the balanced load and unbalanced load are 675 V and 210 V, respectively, the corresponding network-side transformer is used to connect each load. The basic parameters of the two different types of load and the corresponding step-down transformer are given in Table 1. The detailed parameters of the DER system are shown in Table 2.
As for the proposed control scheme, we choose the initial values of the estimated PPD parameter matrix Φ ^ ( k ) as ϕ ^ 1 ( 1 ) = 1.9999 0.00754 1.00754 0.9999 , ϕ ^ 2 ( 1 ) = 6200 5500 5000 6000 , ϕ ^ 3 ( 1 ) = 65 0.0377 55 0.0677 , ϕ ^ 4 ( 1 ) = 0.025 0.02 0.02 0.03 , ϕ ^ 5 ( 1 ) = 7.54 0.2 7.54 0.25 .
Other parameter settings of the proposed control algorithm are K = d i a g ( 0.85 , 0.85 ) , ν 1 = ν 2 = 2 × 10 4 , σ = d i a g ( 0.05 , 0.45 ) . Furthermore, the rate constraints of control input ø ( k ) are
1 1 Δ ø ( k + 1 ) 1 1 .
In the simulation study, U x d r e f is initially set to 0 while the balanced load is connected to the DER system. At t = 0 s, the controllers begin to work and the VSC starts receiving gate signals. Then, U x d r e f is ramped up to 380 V and kept constant from t = 0.5 s onwards. At t = 1.5 s, the balanced load is switched off and the unbalanced load is switched on. At t = 2 s, U x d r e f is changed stepwise from 380 V to 400 V.
In the simulation of the system, we compared the control effect of the designed MFAC controller with that of the FSMC controller and the traditional PI controller (see in Figure 2). As can be seen from Figure 2, when the reference voltage of the d-axis rises slowly, the d-axis voltage of the DER system under the three control methods begins to track the reference voltage. However, the response of PI controlled system is slow. When t = 0.5 s, the reference voltage of the d-axis no longer rises and stays at 380 V. At this time, the system response of PI controller is still very slow when the voltage trend changes, and there is a static error between the d-axis voltage controlled by FSMC and the reference voltage. When t = 2 s, the reference voltage changes from 380 V to 400 V. The d-axis voltage controlled by PI reaches 400 V through 0.05 s, and its change is relatively slow. The d-axis voltage under FSMC control has a small overshoot. Thus, their control effects in the d-axis are not as good as MFAC. Similarly, in the tracking process of q-axis voltage, MFAC control is also superior to the other two control methods. Throughout the whole simulation process, the reference voltage of the q-axis is defined as 0. At the beginning of simulation, the q-axis voltage under FSMC control does not keep up with the reference voltage quickly. When t = 2 s, the q-axis voltage under FSMC control also deviates greatly from the reference voltage. In summary, compared with FSMC and traditional PI control, MFAC has a better control effect during voltage control of the island DER system.
In addition, in Figure 3, Figure 4, Figure 5, Figure 6 and Figure 7, the parameter estimation response curves of ϕ 1 ( k ) ϕ 5 ( k ) are displayed respectively. Figure 8 shows the output voltage waveform at t = 1.9 s–2.1 s. It can be seen that when the reference voltage suddenly changes, the three-phase output voltage of the system changes in time without distortion. Moreover, Figure 2 shows the output current waveform of the system. Compared with Figure 2 and Figure 9, it can be seen that the system voltage under MFAC control remains stable when the current is distorted due to load changes.

5. Conclusions

In this paper, based on the novel transformation and linearization technique for the islanded MIMO DER system dynamics, a multiple observer-based data-driven model free adaptive constrained control scheme was proposed for the voltage control of dispatchable DER system under an unbalanced load condition. Moreover, the proposed control algorithm only depends on the I/O data of the controlled plant and has the advantages of having a small computational burden and strong robustness, which allows easy engineering popularization. In addition, the detailed Lyapunov-based analysis method was presented to guarantee the stability of the closed-loop control system. The simulation model of the DER system was built in MATLAB, and the control effect of MFAC controller was compared with the FSMC and the PI controller, which verifies the advancement of the MFAC controller. In future work, the voltage regulation algorithms for multi-DER systems under the balanced nonlinear load will be the focus of our research.

Author Contributions

Conceptualization, Y.X. and D.X.; Data curation, Y.D. and W.Y.; Formal analysis, Y.X.; Funding acquisition, C.Y.; Methodology, W.Y., D.X. and C.Y.; Project administration, W.Y.; Resources, Y.D.; Software, Y.X. and Y.D.; Writing—original draft, Y.X.; Writing—review and editing, D.X.

Funding

This work was partially supported by the National Natural Science Foundation of China (61503156), National Key Research and Development Program (2016YFD0400301), Open Research Fund of Jiangsu Collaborative Innovation Center for Smart Distribution Network, Nanjing Institute of Technology (XTCX201806), National first-class discipline program of Food Science and Technology (JUFSTR20180205), and State Grid Zhejiang Province Technology Project (SGTYHT/17-JS-201).

Conflicts of Interest

The authors declare no conflict of interest and the founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

References

  1. Karimi, H.; Nikkhajoei, H.; Iravani, M.R. Control of an electronically-coupled distributed resource unit subsequent to an islanding event. IEEE Trans. Power Deliv. 2008, 23, 493–501. [Google Scholar] [CrossRef]
  2. Lew, D.; Asano, M.; Boemer, J.; Ching, C.; Focken, U.; Hydzik, R.; Lange, M.; Motley, A. The power of small: the effects of distributed energy resources on system reliability. IEEE Power Energy Mag. 2017, 15, 50–60. [Google Scholar] [CrossRef]
  3. Sajadian, S.; Ahmadi, R. Model Predictive Control of Dual-Mode Operations Z-Source Inverter: Islanded and Grid-Connected. IEEE Trans. Power Electron. 2018, 33, 4488–4497. [Google Scholar] [CrossRef]
  4. Chen, Z.; Luo, A.; Wang, H.; Chen, Y.; Li, M.; Huang, Y. Adaptive sliding-mode voltage control for inverter operating in islanded mode in microgrid. Electr. Power Energy Syst. 2015, 66, 133–143. [Google Scholar] [CrossRef]
  5. Yazdani, A. Control of an islanded distributed energy resource unit with load-compensating feed-forward. In Proceedings of the IEEE Power and Energy Society General Meeting, PES08, Pittsburgh, PA, USA, 20–24 July 2008. [Google Scholar] [CrossRef]
  6. Karimi, H.; Nikkhajoei, H.; Iravani, R. A linear quadratic gaussian controller for a stand-alone distributed resource unit-simulation case studies. In Proceedings of the IEEE Power Engineering Society General Meeting, PES07, Tampa, FL, USA, 24–28 June 2007. [Google Scholar] [CrossRef]
  7. Li, Y.; Vilathgamuwa, D.M.; Loh, P.C. Design, analysis, and real-time testing of a controller for multibus microgrid system. IEEE Trans. Power Electron. 2004, 19, 1195–1204. [Google Scholar] [CrossRef]
  8. Sao, C.K.; Lehn, P.W. Control and power management of converter fed microgrid. IEEE Trans. Power Syst. 2008, 3, 1088–1098. [Google Scholar] [CrossRef]
  9. Delghavi, M.B.; Yazdani, A. Islanded-mode control of electronically coupled distributed-resource units under unbalanced and nonlinear load conditions. IEEE Trans. Power Deliv. 2011, 26, 661–673. [Google Scholar] [CrossRef]
  10. Huerta, H.; Loukianov, A.; Canedo, J. Multimachine power-system control: integral-SM approach. IEEE Trans. Ind. Electron. 2009, 56, 2229–2236. [Google Scholar] [CrossRef]
  11. Loukianov, A.; Canedo, J.; Utkin, V.; Cabrera-Vazquez, J. Discontinuous controller for power systems: sliding-mode block control approach. Trans. Ind. Electron. 2004, 51, 340–353. [Google Scholar] [CrossRef]
  12. Vargas, J.; Ledwich, G. Variable structure control for power systems stabilization. Electr. Power Energy Syst. 2010, 32, 101–107. [Google Scholar] [CrossRef]
  13. Delghavi, M.B.; Shoja-Majidabad, S.; Yazdani, A. Fractional-Order Sliding-Mode Control of Islanded Distributed Energy Resource Systems. IEEE Trans. Sustain. Energy 2016, 7, 1482–1491. [Google Scholar] [CrossRef]
  14. Chang, F.; Chang, E.; Liang, T.; Chen, J. Digital-signal-processor-based DC/AC inverter with integral-compensation terminal sliding-mode control. IET Power Electron. 2011, 4, 159–167. [Google Scholar] [CrossRef]
  15. Su, X.; Han, M.; Guerrero, J.M.; Sun, H. Microgrid stability controller based on adaptive robust total SMC. Energies 2015, 8, 1784–1801. [Google Scholar] [CrossRef] [Green Version]
  16. Hou, Z.; Jin, S. A novel data-driven control approach for a class of discrete-time nonlinear systems. IEEE Trans. Control Syst. Technol. 2011, 19, 1549–1558. [Google Scholar] [CrossRef]
  17. Hou, Z.; Jin, S. Data driven model-free adaptive control for a class of MIMO nonlinear discrete-time systems. IEEE Trans. Neural Netw. 2011, 22, 2173–2188. [Google Scholar] [CrossRef] [PubMed]
  18. Hou, Z.; Wang, Z. From model-based control to data-driven control: survey, classification and perspective. Inf. Sci. 2013, 235, 3–35. [Google Scholar] [CrossRef]
  19. Xu, D.; Jiang, B.; Shi, P. A novel model free adaptive control design for multivariable industrial processes. IEEE Trans. Ind. Electron. 2014, 61, 6391–6398. [Google Scholar] [CrossRef]
  20. Zhang, H.; Zhou, J.; Sun, Q.; Guerrero, J.M.; Ma, D. Data-Driven Control for Interlinked AC/DC Microgrids Via Model-Free Adaptive Control and Dual-Droop Control. IEEE Trans. Smart Grid 2017, 8, 557–571. [Google Scholar] [CrossRef]
  21. Xu, D.; Jiang, B.; Shi, P. Adaptive observer based data-driven control for nonlinear discrete-time processes. IEEE Trans. Autom. Sci. Eng. 2014, 11, 1037–1045. [Google Scholar] [CrossRef]
  22. Zhu, Y.; Hou, Z. Data-Driven MFAC for a Class of Discrete-Time Nonlinear Systems With RBFNN. IEEE Trans. Neural Netw. Learn. Syst. 2014, 25, 1013–1020. [Google Scholar] [CrossRef] [PubMed]
  23. Zhao, Y.; Yuan, Z.; Lu, C.; Zhang, G.; Li, X.; Chen, Y. Improved model-free adaptive wide-area coordination damping controller for multiple-input-multiple-output power systems. IET Gener. Transm. Distrib. 2016, 10, 3264–3275. [Google Scholar] [CrossRef]
  24. Bu, X.; Hou, Z.; Liang, J.; Xu, P. Data driven multiagent systems consensus tracking using model free adaptive control. In Proceedings of the 2016 28th Chinese Control and Decision Conference (CCDC), Yinchuan, China, 28–30 May 2016. [Google Scholar] [CrossRef]
  25. Lu, X.P.; Li, W.; Lin, Y.G. Load control of wind turbine based on model-free adaptive controller. Trans. Chin. Soc. Agric. Mach. 2011, 42, 109–114. [Google Scholar]
  26. Gao, B.; Gu, K.Y.; Zeng, Y.; Chang, Y. An anti-suction control for an intra-aorta pump using blood assistant index: a numerical simulation. Artif. Organs 2012, 36, 275–285. [Google Scholar] [CrossRef] [PubMed]
  27. Spooner, J.; Maggiore, M.; Passino, K. Stable Adaptive Control and Estimation for Nonlinear Systems; Wiley Interscience: Hoboken, NJ, USA, 2002; Volume 36, pp. 275–285. ISBN 0471415464. [Google Scholar] [CrossRef]
  28. Astrom, K.; Johansson, K.; Wang, Q. Design of decoupled PI controllers for two-by-two systems. IEE Proc. Control Theory Appl. 2002, 149, 74–81. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Schematic diagram of an islanded DER system.
Figure 1. Schematic diagram of an islanded DER system.
Energies 11 03299 g001
Figure 2. The tracking effect on the reference voltage.
Figure 2. The tracking effect on the reference voltage.
Energies 11 03299 g002
Figure 3. Parameter estimation response curve of ϕ 1 ( k ) .
Figure 3. Parameter estimation response curve of ϕ 1 ( k ) .
Energies 11 03299 g003
Figure 4. Parameter estimation response curve of ϕ 2 ( k ) .
Figure 4. Parameter estimation response curve of ϕ 2 ( k ) .
Energies 11 03299 g004
Figure 5. Parameter estimation response curve of ϕ 3 ( k ) .
Figure 5. Parameter estimation response curve of ϕ 3 ( k ) .
Energies 11 03299 g005
Figure 6. Parameter estimation response curve of ϕ 4 ( k ) .
Figure 6. Parameter estimation response curve of ϕ 4 ( k ) .
Energies 11 03299 g006
Figure 7. Parameter estimation response curve of ϕ 5 ( k ) .
Figure 7. Parameter estimation response curve of ϕ 5 ( k ) .
Energies 11 03299 g007
Figure 8. Waveforms of U x a b c under the unbalanced load condition when the proposed control scheme is employed.
Figure 8. Waveforms of U x a b c under the unbalanced load condition when the proposed control scheme is employed.
Energies 11 03299 g008
Figure 9. Waveforms of I x a b c under the unbalanced load condition when the proposed control scheme is employed.
Figure 9. Waveforms of I x a b c under the unbalanced load condition when the proposed control scheme is employed.
Energies 11 03299 g009
Table 1. Load Parameters.
Table 1. Load Parameters.
Load TypesCharacteristicsTransformer
Balancedthree star-connected T r 1
Loadseries RL circuits with4.13/0.675
R = 155 m Ω andkVrms
L = 200 μ H
Unbalancedseries RL between each phase T r 2
Loadand neutral with4.13/0.21
R a = 65 m Ω , L a = 43 μ HkVrms
R b = 18 m Ω , L b = 22 μ H
R c = 25 m Ω , L c = 30 μ H
Table 2. The Parameters of the Distributed Energy Resource (DER) System Circuit.
Table 2. The Parameters of the Distributed Energy Resource (DER) System Circuit.
ParametersValueDescription
L s 295 μ H
C 500 μ F
R 3.2 m Ω
U d c 1800 Vdc-link voltage
w314 rad/snominal angular frequency
T r 4.13/ 0.675 kVrmsleakage reactance 8 %

Share and Cite

MDPI and ACS Style

Xia, Y.; Dai, Y.; Yan, W.; Xu, D.; Yang, C. Adaptive-Observer-Based Data Driven Voltage Control in Islanded-Mode of Distributed Energy Resource Systems. Energies 2018, 11, 3299. https://doi.org/10.3390/en11123299

AMA Style

Xia Y, Dai Y, Yan W, Xu D, Yang C. Adaptive-Observer-Based Data Driven Voltage Control in Islanded-Mode of Distributed Energy Resource Systems. Energies. 2018; 11(12):3299. https://doi.org/10.3390/en11123299

Chicago/Turabian Style

Xia, Yan, Yuchen Dai, Wenxu Yan, Dezhi Xu, and Chengshun Yang. 2018. "Adaptive-Observer-Based Data Driven Voltage Control in Islanded-Mode of Distributed Energy Resource Systems" Energies 11, no. 12: 3299. https://doi.org/10.3390/en11123299

APA Style

Xia, Y., Dai, Y., Yan, W., Xu, D., & Yang, C. (2018). Adaptive-Observer-Based Data Driven Voltage Control in Islanded-Mode of Distributed Energy Resource Systems. Energies, 11(12), 3299. https://doi.org/10.3390/en11123299

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop