Next Article in Journal
Research on the Impact of Various Emission Reduction Policies on China’s Iron and Steel Industry Production and Economic Level under the Carbon Trading Mechanism
Next Article in Special Issue
Decentralized Optimal Control of a Microgrid with Solar PV, BESS and Thermostatically Controlled Loads
Previous Article in Journal
Power-mileage-based Algorithm for the Optimization of Distribution Network Structures via Adding Transmission Lines
Previous Article in Special Issue
An Energy Management System for Residential Autonomous DC Microgrid Using Optimized Fuzzy Logic Controller Considering Economic Dispatch
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Robust Control Method for DC Microgrids and Energy Routers to Improve Voltage Stability in Energy Internet

1
Research Institute of Information Technology, Beijing National Research Center for Information Science and Technology, Tsinghua University, Beijing 100084, China
2
School of Electrical Engineering, Chongqing University, Chongqing 400044, China
*
Author to whom correspondence should be addressed.
Energies 2019, 12(9), 1622; https://doi.org/10.3390/en12091622
Submission received: 23 March 2019 / Revised: 23 April 2019 / Accepted: 23 April 2019 / Published: 29 April 2019
(This article belongs to the Special Issue Advanced Control in Microgrid Systems)

Abstract

:
The energy internet (EI) is a wide area power network that efficiently combines new energy technology and information technology, resulting in bidirectional on-demand power transmission and rational utilization of distributed energy resources (DERs). Since the stability of local network is a prerequisite for the normal operation of the entire EI, the direct current (DC) bus voltage stabilization for each individual DC microgrid (MG) is a core issue. In this paper, the dynamics of the EI system is modeled with a continuous stochastic system, which simultaneously considers related time-varying delays and norm-bounded modeling uncertainty. Meanwhile, the voltage stabilization issue is converted into a robust H control problem solved via a linear matrix inequality approach. To avoid the situation of over-control, constraints are set in controllers. The problem of finding a balance between voltage regulation performance and constraints for the controllers was also extensively investigated. Finally, the efficacy of the proposed methods is evaluated with numerical simulations.

1. Introduction

In the past few decades, the usage of distributed energy resources (DERs) has attracted a significant amount of attention due to their sustainability and environmental friendliness [1,2]. It is notable that the natural characteristics of power generation by DERs, e.g., wind turbine generators (WTs) and photovoltaic panels (PVs), are heavily weather-dependent and over-flexible, which has a significant impact on the accessed conventional power systems [3]. Meanwhile, the convergence of multiple energy resources is the current development trend of future energy systems [4]. For conventional power grids, being able to access renewable energy sources would create challenges in the operation of power grids with conventional power regulation modes, which is why solar power and wind power were previously regarded as ‘waste energy’. To further promote the efficiency of multi-energy conversion and utilization, the architecture of energy internet (EI) is proposed as being the future of the smart grid [5].
EI is a cyber-physical system inspired by the concept of the Internet, enabling the bidirectional power transmission and on-demand rational energy utilization [6]. Within the scenario of future EI, solar power, wind power and hydro power can be viewed as the main power generation sources, which can be integrated with energy storage devices and various loads. Typically, EI is composed of a utility grid and multiple microgrids (MGs) interconnected via energy routers (ERs). ERs are the core infrastructures of EI systems. The ER, also known as an electric energy router or power router, is a type of electric device that can realize multidirectional power flow and active control of the power flow. In the distribution network, acting as an intelligent interface of distributed power supply, reactive power compensation device, energy storage equipment and load, the ER can flexibly manage the dynamic power in the regional grid network and the whole distribution network on the premise of ensuring power quality. Integration with advanced information technology enables ERs to have communication and intelligent decision-making capabilities. ERs can actively manage the energy flow of the power network according to the operation status of the network and the instructions from the user and control center. The basic function of ER includes plug and play of interfaces, bidirectional transmission; real-time communications and so on. For ERs designed for a backbone network, the main difficulty lies in the technology of high voltage power electronic converters. For small and medium-sized ERs, a high manufacturing cost is the reason why it has not been popularized in the market. The current researches in both academia and industry aim to propose ERs for future energy markets. For a more detailed introduction with respect to ER, readers can consult [6,7,8] and the references therein.
As direct current (DC) MGs not only can improve the power quality and transmission capacity of an EI system but also have superior accesses to DERs, they are more widely used in real projects [9]. Due to the direct influence of power deviation on MG’s DC bus voltage, voltage regulation is usually realized by power control. In an EI scenario, the power balance of the whole system should be achieved by the individual local controllable devices in each MG through prioritizing [10]. However, when faced with extreme conditions, such as a large-scale load power change, MGs tend to be regulated by means of ER operations, which aim to transmit power energy from/to other MGs in order to achieve power supply–demand balance within the whole EI scenario [8,10].
On the other hand, there have been significant achievements in terms of voltage control in related literature, such as [11,12,13,14,15,16]. To illustrate, a local reactive power control scheme that can quickly respond to voltage deviations for preventing communication delay and noise has been previously studied [11]. In [12], the authors proposed voltage control strategies to exploit the cooperation among the agents, such that the power loss minimization objective is achieved. An adaptive control approach for DC MG systems, satisfying both accurate power sharing and voltage regulation, has been previously investigated [13]. Furthermore, a droop-like feedback control method was developed for voltage regulation in [14], enabling the application of theoretical circuits analytical techniques. In [15], with techniques introduced in [17] and the references therein, the problem of voltage regulation is investigated for an islanded MG in EI, which considers system stochasticity and parameter uncertainties.
Naturally, the random change of photovoltaic and wind power as well as the custom of using electricity would stochastically influence PV power output, WT power output and load power, respectively. In [16] and [18], the power of PVs, WTs and loads is modelled using ordinary differential equations (ODEs). As an improvement of ODEs, continuous stochastic differential equations (SDEs) driven by Brownian motions have been applied to model deviations in the power input/output of PVs, WTs and loads in [15,19,20]. In [10] and [21], the Ornstein-Uhlenbeck process was utilized for the power modelling of power fluctuations for PVs and loads.
On the other hand, there are objective measurement errors for any system parameter. When we consider the short-term dynamical properties of the power system, especially for the problem of MG voltage regulation, there is the need for an accurate power dynamical model. Thus, parameter uncertainty must be considered in the dynamical MG model. There exist a variety of approaches to describe parameter uncertainty, such as [22,23]. In this paper, we choose the norm-bounded parameter uncertainty to represent such measurement errors, with the detailed structure introduced in Section 2.
Furthermore, there exists a communication time delay within the information system of EI [6,16]. However, in the aforementioned literature [11,12,13,14,15,16] and the references therein, system stochasticity, parameter uncertainty and communication time delay have not been taken into consideration simultaneously.
In this paper, it is supposed that the considered EI scenario is disconnected with the main power grid and such an EI consists of multiple DC MGs, which are interconnected via ERs. We assume that each MG is composed of PVs, WTs, loads, micro-turbines (MTs), diesel engine generators (DGs), fuel cells (FCs), battery energy storages (BESs) and flywheel energy storages (FESs). The dynamics of the considered system are modelled as continuous differential equations. Due to the randomness of power output of PVs and WTs as well as the stochasticity in the loads, the power deviations of PVs, WTs and loads are modelled as SDEs. We model the dynamics of MTs, DGs, FCs, BESs, FESs and each MG’s DC bus voltage deviations as ODEs. Parameter uncertainties are considered in the system coefficients, which revealed measurement and modelling errors. In addition, the time delay occurring in the communication system is considered in the dynamical equations of the controllable electrical devices. After this, we formulate the voltage control issues in EI as a robust H control problem. The linear matrix inequality (LMI) approach [17] is utilized to solve this problem. Furthermore, the problem of over-control has been concerned and extensively studied. The efficacy of the proposed controller is presented in the numerical simulation.
The main contributions and highlights of this paper are outlined as follows.
(1) The assumption that the EI scenario is functioning without access to main power grid makes the control problem more challenging compared to that with access to the utility grid. Each DC MG is analyzed in detail, with the dynamics of the equipment in MGs (including PVs, WTs, loads, MTs, DGs, FCs, BESs, FESs, ERs) being analyzed. The considered EI dynamical system is complex and authentic, which is reflected in the following three aspects. First, corresponding to the fusion of energy and information in EI, the communication time delay is taken into consideration when formulating the dynamical equations of ERs. Second, when investigating the transient power dynamics in MGs, system modelling errors are accounted. Norm-bounded time-varying parameter uncertainty is used to represent such errors. Third, the stochastic nature of solar power, wind power and loads are fully considered and their power dynamical equations are modelled with SDEs. The voltage regulation issue based on this particular system where communication time delay, system parameter uncertainty and system stochasticity are considered simultaneously has not been considered before within the scope of EI.
(2) We emphasize that our contribution is on modelling and formulating the engineering problem in EI into control issues that we can solve. It is challenging to determine how to transform our considered physical problem in an EI scenario into a control problem itself. The practical voltage regulation problem is converted into a robust H control problem in the time domain. The LMI approach is utilized to solve such a control problem sufficiently in order to achieve both robust stabilization against system internal uncertainty and H performance against external disturbance input.
(3) Based on the sufficient solutions to the robust H voltage regulation issue, an additional constraint regarding the size of controllers was considered, such that the potential over-control for distributed power generators would be effectively avoided. Furthermore, the problem of finding a balance between the voltage regulation performance and the constraints for the controllers has been extensively studied, which has not been considered before. Comparable numerical simulations are performed, demonstrating the effectiveness of the proposed control approach.
The rest of the paper is organized as follows. Section 2 introduces the ER system modelling. Problem formulation and solutions are given in Section 3. In Section 4, numerical examples are illustrated. Section 5 concludes the paper.

2. System Description and Dynamical Modelling

In this section, the system modelling for the considered EI scenario is established.

2.1. EI System Architecture

A typical scenario of the considered EI is given in Figure 1. It is assumed that there are a total of n MGs and n ERs in the considered EI system. The power transmission in the EI system is achieved via the ER network. The structure of the ER network can be described using the terminologies in graph theory. ERs can be seen as the vertices of a graph and the transmission lines between ERs in the ER network can be viewed as the edges in the graph. With the set V defined as V = { 1 , 2 , , n } , MGs and ERs are numbered as M G i and E R i , i V , respectively.

2.2. System Modelling of Each MG Component

In real-world EI systems, a variety of power generation, energy storage and consumption devices can be installed in MGs. For example, some MGs may have both PVs and WTs installed, while only distributed PVs and MTs are installed for some other MGs. For practical considerations, a certain type of MG system consisting of the most common power devices is considered in this paper. Each MG is supposed to include PVs, WTs, MTs, FCs, DGs, BESs, FESs and loads. In this sense, this type of model could be applicable for various EI scenarios. The power dynamics for these components are introduced below.

2.2.1. Modelling for Power Dynamics of PVs, WTs and Loads

For the i th MG, considering the system parameter uncertainty, we use a continuous-time stochastic system to describe the power deviation of PVs, WTs and loads. In this paper, the scalar Wiener processes w P V i ( t ) , w W T i ( t ) and w l o a d i ( t ) are utilized to represent the stochasticity in Δ P P V i ( t ) , Δ P W T i ( t ) and Δ P l o a d i ( t ) , respectively. Based on the models introduced in [15], the power dynamics of Δ P P V i , Δ P W T i , Δ P l o a d i are expressed as follows:
d Δ P M i ( t ) = 1 T M i [ ( 1 + Δ a M i ( t ) ) Δ P M i ( t ) + v M i ( t ) ] d t + e M i Δ P M i ( t ) d w M i ( t ) ,
where M { P V , W T , l o a d } ; v P V i ( t ) , v W T i ( t ) and v l o a d i ( t ) are modelled as disturbance inputs with respect to Δ P P V i ,   Δ P W T i and Δ P l o a d i . The terms Δ a P V i ( t ) , Δ a W T i ( t ) and Δ a l o a d i ( t ) in (1) are bounded time-varying system parameter uncertainties. The time constants T P V i , T W T i , T l o a d i and system coefficients e P V i , e W T i , e l o a d i can be obtained with parameter estimation methods.

2.2.2. Power Dynamics of MTs, DGs and FCs

Due to the integration of DERs, which could lead to large power deviation in MG systems, the voltage stability in MGs would become more fragile. In future EI systems, MGs will have a higher penetration rate for renewable energy sources. Thus, it is crucial to ensure the power supply–demand balance in all MG systems. Generally, MTs, DGs and FCs are highly controllable electric devices in common MG systems. By adjusting the output power of these generators, the voltage fluctuations caused by local loads and DERs could be effectively alleviated. Similar to the output power deviation of MTs, FCs and DGs given in [10,15,16], we defined the following equation:
d Δ P N i ( t ) = 1 T N i [ ( 1 + Δ a N i ( t ) ) Δ P N i ( t ) + ( 1 + Δ b N i ( t ) ) u N i ( t ) ] d t ,
where N { M T , D G , F C } ; time constants T M T i , T D G i , T F C i indicating the response speed of generators could be obtained based on real-world data; terms Δ a M T i ( t ) , Δ b M T i ( t ) , Δ a D G i ( t ) , Δ b D G i ( t ) , Δ a F C i ( t ) and Δ b F C i ( t ) are bounded time-varying parameter uncertainties.

2.2.3. Power Dynamics of BESs and FESs

Storage systems are essential pieces of equipment for the normal operation of MG systems. Compared with the controllable power generation devices, such as MTs, DGs and FCs, the energy storage devices have shorter response times [16,20,24]. Thus, BESs and FESs are more widely applied to absorb the power deviation in MG systems.
However, due to the limitations in current energy storage technology, some of the energy storage devices are still relatively expensive [20]. Thus, inappropriate control schemes for BESs and FESs may lead to an increase in the operation costs. For better performance and efficiency, different operation schemes for energy storage systems have been proposed, such as the one proposed in [25]. Most of these works focus on designing control schemes within the energy storage devices themselves, which means that energy storage devices are controlled internally. This is different to the situation in this present paper as no controllers are set in BESs and FESs.
With the similar linearization techniques introduced in [15,18,20], the dynamics of power deviations of BESs and FESs in M G i are given as follows:
d Δ P S i ( t ) = 1 T S i [ ( 1 + Δ a S i ( t ) ) Δ P S i ( t ) + ( 1 + Δ o S i ( t ) ) Δ V i ( t ) ] d t ,
where S { B E S , F E S } ; and Δ a B E S i ( t ) , Δ o B E S i ( t ) , Δ a F E S i ( t ) and Δ o F E S i ( t ) are bounded time-varying parameter uncertainties.

2.2.4. Power Dynamics of ERs

In this paper, ERs are mainly applied for energy exchange among MGs. The functionality of ERs allows the administrator to control the energy flows in the EI system more actively [6,10]. Each MG is assumed to have a connection with an ER, which is interconnected with other ERs. ERs can either function independently or cooperate with other ERs. The energy flow can be transmitted from M G i to M G j via E R i and E R j within the ER network. For any i < j ,   i , j = 1 , 2 , , n , the power transmitted from E R i to E R j at time t is denoted as P E R i j ( t ) . The notation P E R i j ( t ) with the superscripts satisfying i j will not be used throughout this paper.
The default direction of these power flows could be arbitrarily set and their signs would indicate the actual transmission directions. For example, we are able to set the default power transmission direction between M G i and M G j as from M G i to M G j . In this sense, P E R i j ( t ) > 0 means that the power is sent from M G i and accepted by M G j . Reversely, with P E R i j ( t ) < 0 , the power flow should start at M G j and end at M G i . Based on the default directions, the influences of ERs on the power imbalance in MGs should be calculated accordingly, which would be reflected in (6) and (7).
In the considered EI scenario, the control schemes for ERs and MGs are usually generated by a regional control center in order to achieve better efficiency and flexibility. Delays in transmission lines and communication systems may affect the power flows in the ER network [26,27]. Meanwhile, the power flow in a transmission line might be affected by the engineering environment and ERs at both ends, which is considered as the disturbance for the power dynamics. In this paper, the dynamical power flow P E R i j ( t ) is approximated as a linear differential equation with control and disturbance inputs, as described in (4). For i < j , i , j V , we have:
d Δ P E R i , j ( t ) = 1 T E R i , j [ ( 1 + Δ a E R i , j ( t ) ) Δ P E R i , j ( t ) + v E R i , j ( t ) + ( a d E R i , j + Δ a d E R i , j ( t ) ) Δ P E R i , j ( t τ ( t ) ) + ( 1 + Δ b E R i , j ( t ) ) u E R i , j ( t ) ] d t ,
where a d E R i , j is a coefficient related to the impacts of the aforementioned delays; and Δ a E R i , j ( t ) , Δ a d E R i , j ( t ) and Δ b E R i , j ( t ) are bounded time-varying parameter uncertainties. The notation u E R i , j ( t ) refers to the controller for the energy transmission between E R i and E R j . In (4), the term Δ P E R i , j ( t τ ( t ) ) reflects the impact of delays in the EI system, where τ ( t ) is a time-varying function satisfying 0 < τ ( t ) μ < , and τ ( t ) ˙ h < 1 .

2.2.5. Power Dynamics of DC Bus Voltage Deviation

Considering the i th MG, we denote the voltage deviation on the power bus of M G i as Δ V i ( t ) . According to [28], the voltage on the power bus is influenced by the external power from ERs and the internal power deviation in M G i . Thus, the dynamic of Δ V i ( t ) can be formulated as:
d Δ V i ( t ) = 1 T V i [ ( 1 + Δ a V i ( t ) ) ( Δ P i n i ( t ) Δ P m i ( t ) ) ] d t ,
where Δ a V i ( t ) is the bounded time-varying parameter uncertainty; Δ P m i ( t ) is the power deviation inside M G i defined in (6); and Δ P i n i ( t ) is the power input for M G i that comes from E R i .
Δ P m i ( t ) = Δ P l i ( t ) Δ P P V i ( t ) Δ P W i ( t ) Δ P B E S i ( t ) Δ P D G i ( t ) Δ P F C i ( t ) Δ P M T i ( t ) Δ P F E S i ( t ) .
Based on the adjacency matrix R of the ER network, Δ P i n i ( t ) can be expressed as:
Δ P i n i ( t ) = j V r j , i P E R j , i ( t ) ,
where r i , j is the element of R , indicating the connectivity from E R i to E R j .

3. Problem Formulation and Solution

In this section, the voltage regulation problem for the considered EI system is formulated as a robust stochastic control issue. To avoid drastic voltage deviation in MGs, the desired controllers in MTs, DGs, FCs and ERs are sufficiently solved using the robust H control method. Since the obtained sufficient solutions might be too strong, it is possible that the situation of over-control regarding MTs, DGs, FCs and ERs might occur. Thus, some extra constraints for the controllers are proposed in this section in order to successfully avoid the situation of over-control.

3.1. The Mathematical Formulation of EI Dynamical System

Next, the EI system power dynamics is represented with a unified dynamical system, which is introduced as follows. For simplicity, we assume that w P V i ( t ) = w W i ( t ) = w l i ( t ) = w ( t ) ,   i = 1 , 2 , , n , which would not lead to significant changes to the considered problem. The similar assumptions are applied in many works, such as [15].
For the i th MG, the system state is denoted as the vector x i ( t ) 9 , which is defined in (8):
x i ( t ) [ Δ P P V i ( t ) Δ P W T i ( t ) Δ P l o a d i ( t ) Δ P M T i ( t ) Δ P D G i ( t ) Δ P F C i ( t ) Δ P B E S i ( t ) Δ P F E S i ( t ) Δ V i ( t ) ] .
Similarly, for M G i , the disturbance input vector v i ( t ) 3 and control input vector u i ( t ) 3 are denoted as v i ( t ) = [ v φ i ( t ) v w i n d i ( t ) v l o a d i ( t ) ]   and u i ( t ) = [ u M T i ( t ) u D G i ( t ) u F C i ( t ) ] , respectively.
Based on the adjacency matrix R , the connections between ERs can be described with two tuples. Suppose that there are a total of m transmission lines in the EI system, we denote the set of the two tuples as L = { ( i , j ) | i < j ,   r i , j = 1 , i , j = 1 , 2 , , n } . We denote the elements in L as l 1 , l 2 , , l m , such that the state vector x E R ( t ) m , control input vector u E R ( t ) m and disturbance input vector v E R ( t ) m for ERs can be denoted as:
x E R ( t ) = [ Δ P E R l 1 ( t ) Δ P E R l i ( t ) Δ P E R l m ( t ) ] , u E R ( t ) = [ u E R l 1 ( t ) u E R l i ( t ) u E R l m ( t ) ] , v E R ( t ) = [ v E R l 1 ( k ) v E R l i ( k ) v E R l m ( t ) ] ,
respectively. In this sense, we are able to denote the EI system state variable x ( t ) 9 n + m , control input u ( t ) 3 n + m and disturbance input v ( t ) 3 n + m of the entire EI system as:
x ( t ) = [ x 1 ( t ) x n ( t ) x E R ( t ) ] , u ( t ) = [ u 1 ( t ) u n ( t ) u E R ( t ) ] , v ( t ) = [ v 1 ( t ) v n ( t ) v E R ( t ) ] .
At time t , the voltage deviations z ( t ) n within n MGs in EI are denoted as z ( t ) = [ Δ V 1 ( t ) Δ V i ( t ) Δ V n ( t ) ] . As the dynamical equations in (1)–(5) are expressed in linear forms, the EI dynamical system can be rewritten into the following explicit form:
{ d x ( t ) = [ A ( t ) x ( t ) + A d ( t ) x ( t τ ( t ) ) + D v ( t ) + B ( t ) u ( t ) ] d t + E x ( t ) d w ( t ) , z ( t ) = C x ( t ) , x ( t ) = ϕ ( t ) ,     μ t 0 ,
where ϕ ( t ) ,     μ t 0 is the initial state. The matrix coefficients in (11) are obtained by:
A ( t ) = A + Δ A ( t ) , A d ( t ) = A d + Δ A d ( t ) , B ( t ) = B + Δ B ( t ) ,
with appropriate dimensions. It is assumed that the bounded time-varying uncertainty matrices in (12) satisfies:
[ Δ A ( t ) Δ A d ( t ) Δ B ( t ) ] = M F ( t ) [ N a N a d N b ] ,
where M , N a , N a d and N b are real constant matrices while F ( · ) is an arbitrary time-varying matrix function satisfying:
F ( t ) F ( t ) I ,   t +
It is notable that the above structure for parameter uncertainties has been used in many works in the field of power systems, such as [15]. Hence, the dynamic system for the considered EI scenario has been explicitly obtained.
Based on the system modeling given above, a general H control configuration for the voltage regulation problem considered in this paper is shown in Figure 2. The main target of the remainder of this paper is to develop a type of state feedback controller u = K x , such that the controlled output z , i.e., the voltage deviations, could be properly restricted in order to remove any influence from the unmodeled disturbance input v .

3.2. Robust H Performance

For the considered EI system, the so-called robust H performance particularly refers to the good regulation of the voltage in each MG with respect to two properties: 1) EI system’s internal parameter uncertainty; and 2) external disturbance input from solar irradiation, wind power change, change of electricity usage and disturbing power transmission from other interconnected ERs. Mathematically, we provide a definition of robustly stable (or equivalently, robust stability) for the EI system (11) as follows.
Definition 1:
For dynamical EI system (11), the controlled system with disturbance input v ( t ) = 0 is said to be robustly stable if for all bounded time-varying parameter uncertainties Δ A ( t ) , Δ A d ( t ) , Δ B ( t ) in forms of (13), l i m t E || x ( t ) || 2 2 = 0 holds. This property is called robust stability. The notation || || 2 refers to the Euclidean norm.
In addition to the property of robust stability against EI system’s internal parameter uncertainty, the impact of the disturbance input v ( k ) on voltage deviation should be restricted within certain level, which corresponds to the robust H performance in the control theory. The robust H performance for the investigated EI system is defined in Definition 2.
Definition 2:
For dynamical EI system (11), with a controller u ( x , t ) under a disturbance attenuation level γ > 0 , if the controlled EI dynamical system (11) is robustly stable under Definition 1 for all time-varying parameter uncertainties Δ A ( t ) , Δ A d ( t ) , Δ B ( t ) , Δ E ( t ) in forms of (11) and if || z ( t ) || E 2 γ || v ( t ) || E 2 holds for all nonzero v ( k ) , where || z ( t ) || E 2 = ( E { 0 z ( t ) z ( t )   d t } ) 1 / 2 , the robust H performance for the considered EI system is said to be achieved.
The definitions of robust stability and robust H performance introduced above originate from the robust H control theory in the time domain, such as [29]. Following this, the robust H control theory can be applied to obtain the desired controllers. A theorem below is provided to solve our voltage regulation problem.
Theorem 1:
System (11) is robustly stabilized by the linear state-feedback controller u ( x , t ) = Y X 1 x ( t ) under the given disturbance attenuation factor γ if there exist symmetric matrices X > 0 , S > 0 , matrix Y and a scalar ε > 0 , such that the LMI in (15) holds:
Γ < 0 ,
where Γ is defined in (16), Ω in Γ is given in (17).
Γ [ Ω A d X D X N a + Y N b X E X C X A d ( 1 h ) S 0 X N a d 0 0 D 0 γ 2 I 0 0 0 N a X + N b Y N a d X 0 ε I 0 0 E X 0 0 0 X 0 C X 0 0 0 0 I ] ,
Ω = A X + X A + B Y + Y B + ε 1 M M + S .
The EI dynamical system (11) is a special form of the system equation analyzed in [29]. Our obtained theorem can be viewed as a special form of the results in [29]. Thus, the proof is omitted.
The controller obtained from Theorem 1 can be successfully applied to the real-world EI system with guaranteed performance and robustness property. From Definition 2, we can determine that the performance of voltage regulation in the EI system is directly related to the disturbance attenuation level γ . The controller obtained with a smaller disturbance attenuation level ensures better regulation effect against the disturbances in the EI system.

3.3. Constraints for the Controllers

In [29], it is pointed out that the obtained robust H controllers are sufficient solutions rather than both necessary and sufficient ones, which means that the controller might be excessively strong when a satisfactory robust H performance is achieved. In real engineering practice, voltage stabilization is the main concern in some cases. In order to achieve better power quality, strong controllers are necessary. In some other cases, the magnitudes of the feedback gains should be restricted to avoid the situation of over-control, which might potentially damage the controllable electrical devices.
Thus, in addition to the robust H controllers obtained in Theorem 1, it is necessary to set some constraints in the controllers. Meanwhile, finding a balance between the voltage regulation performance and the constraints for the controllers is still an open problem for the studied system (11), which is investigated below.
Generally, we can manually set the value of disturbance attenuation level γ in some applications. The controllers obtained from any feasible solution to (15) can satisfy normal operation requirements in EI systems. However, we would occasionally like to find a controller with the smallest γ , such that the impacts on DC bus voltages from the disturbance inputs of the EI system can be minimized. This could be achieved by solving the problem in (18):
min γ 2 , s . t . Γ < 0 .
Supposing the solution to (18) is X ( 0 ) , Y ( 0 ) , S ( 0 ) , ε ( 0 ) , γ ( 0 ) , the controller that corresponds to γ ( 0 ) is obtained with u ( 0 ) ( x , t ) = Y ( 0 ) X ( 0 ) 1 x ( t ) . With u ( 0 ) ( x , t ) , the voltage in the EI system (11) should be more stable than that under Theorem 1.
However, a smaller disturbance attenuation level usually suggests larger magnitude of the controller’s feedback gain, which might not be appropriate in real EI scenarios. A feedback controller with large feedback gain might lead to violent adjustments for the distributed generators, which would lead to potential damages. To avoid the situation of over-control, the magnitude of feedback gain needs to be restricted. With Schur complement lemma, it is apparent that:
[ Q I I X ] 0 Q X 1 0 ,
where Q is a symmetric matrix. After this, we have || Q || p || X 1 || p , where || || p ( 1 p ) stands for the p -norm. As that the feedback gain K obtained from Theorem 1 satisfies K = Y X 1 , || K || p 2 1 2 ( || X 1 || p 2 + || Y || p 2 ) is established. By minimizing the upper bound for the feedback gain K , the intensity of the obtained controller would be consequently restricted, which would reduce the requirement of the proposed controller in practical applications.
Thus, with the idea to minimize γ and the LMIs in (15) and (19), we are able to construct a new convex optimization problem as follows:
min η γ 2 + || Q || p 2 + || Y || p 2 , subject   to   Γ < 0 , [ Q I I X ] 0 .
where η is a positive scalar and p is a scalar satisfying p 1 .
Here, we denote the controller corresponding to the solution to (20) as u ( 1 ) ( x , t ) . The solutions to problem (18) and (20) can be obtained with the convex optimization toolbox CVX [30].

4. Numerical Simulation

In this section, the numerical simulations of system (11) is presented to show the efficacy and feasibility of the methods proposed in this paper.
For illustrative purposes, an EI system (without access to utility grid) with four MGs and four ERs is investigated in the following simulation. The considered MGs and ERs are denoted as M G 1 , M G 2 , M G 3 , M G 4 , E R 1 , E R 2 , E R 3 and E R 4 , respectively.
Let us define the element A ( i , j ) of matrix A as the connectivity between M G i and M G j . In this sense, when M G i is connected with M G j via ERs, A ( i , j ) = 1 and A ( i , j ) = 0 otherwise.
The adjacency matrix for the considered EI system is shown as follows:
A = [ 0 1 1 1 1 0 1 0 1 1 0 0 1 0 0 1 ] .
Some typical values of time constants and other system parameters in (11) are provided in Table 1. Readers may refer to [15,16,18] for typical MG system parameters. During the simulation, EI system parameters are randomly generated using the product of corresponding parameters in Table 1 and a random variable, which follows the uniform distribution in [ 0.9 , 1.1 ] . To illustrate, the time constants for PVs in these four MGs are 1.9 times the four identically individually distributed random variables within [ 0.9 , 1.1 ] . The time period for simulation is set as t [ 0 , 10 ] (time unit second omitted). The infinity-norm of the parameter uncertainty Δ A ( t ) , Δ A d ( t ) and Δ B ( t ) in system (11) is assumed to be less than 0.05. For the time delay τ ( t ) in (11), the constant h is assigned as 0.3. For problem (20), η is set to be 10,000. The 2-norm is used to calculate the induced norm for matrices Q and Y . In this sense, we have p = 2 .
Based on the dynamical models introduced in Section 2 and the parameters described above, we are able to obtain the parameter matrices in system (11). Both of the solutions to problems (18) and (20) are obtained with the CVX toolbox [29] in MATLAB (R2018b, MathWorks Inc., Natick, MA, USA) environment.
In order to achieve notation simplicity, let us denote the desired controllers solved by problems (18) and (20) as u 0 and u 1 , respectively. When solving problems (18) and (20), the dimensions of the obtained matrices X , Y are rather large, which are inconvenient to be explicitly presented in this paper. To show the differences between the solutions to problems (18) and (20), some of the solutions are given in Table 2.
In order to show the efficacy of the proposed controllers more intuitively, a detailed analysis is provided as follows.
The time delay τ ( t ) and the factor F ( t ) for parameter uncertainties during the simulation period are plotted in Figure 3 and Figure 4, respectively. Basically, the time delay is randomly generated based on the constraints 0 < τ ( t ) μ < , τ ( t ) ˙ h < 1 . For the sake of simplicity, the factor F ( t ) is generated as a random variable located in [ 1 , 1 ] .
For illustrative purposes, we plot the power deviations of PVs, WTs and loads in M G 1 in Figure 5. It is clear that the stochastic power characteristics of PVs, WTs and loads are properly approximated.
Based on the curve of time delay τ ( t ) , factor for parameter uncertainties F ( t ) and dynamics of PVs, WTs and loads, we are able to evaluate the effectiveness of the controllers obtained from problems (18) and (20) as follows.
The voltage deviations in four MGs are illustrated in Figure 6. The dashed red curves in Figure 6 refer to the dynamics of voltage deviations when there is no controller applied to the EI system. The dynamics of voltage deviations in the four MGs under u 0 and u 1 are illustrated with solid blue lines and dashed–dotted green lines, respectively. As stated above, BESs and FESs in MGs are treated as uncontrollable devices. The power deviation in MGs can be absorbed by BESs and FESs. Thus, even if no controller is applied, the voltage deviations in MGs can still be restricted to lie within a relatively small level. However, with the cooperation of generators and ERs, better voltage regulation can be achieved.
In Figure 6, it is obvious that when there is no controller applied, the voltage deviations in MGs are relatively larger. We can also find that both u 0 and u 1 can effectively achieve the voltage stabilization target. The voltage is further stabilized under u 0 , which is consistent with the values for disturbance attenuation γ in Table 2. Actually, the disturbance attenuation factor corresponding to u 0 is the smallest one that can be achieved with the proposed method. Thus, the stability achieved by u 0 is no doubt most excellent.
To show the controlled H performance, the quantitative comparison for the norms of H bound γ u 0 || v || and observed state z are presented in Figure 7.
In Figure 7, the factor γ u 0 = 0.0248 corresponds to controller u 0 ; z 0 is the observation result when there is no control input for system (11); and z u 0 and z u 1 refer to the observations of the system (11) when the controllers u 0 and u 1 are applied, respectively. According to Figure 7, || z u 0 || and || z u 1 || are clearly smaller than γ u 0 || v || , while || z 0 || varies in a relative larger range. It is clear that for given disturbance attenuation levels, the studied EI system is robustly stabilized by u 0 and u 1 , respectively. We can find that the deviations of voltage magnitudes exceed the desired upper bound γ u 0 || v || at around 3 s. However, the controller u 1 with a medium feedback gain effectively maintained the voltage deviations within a small range. Further, the controller u 0 with the strongest stabilization ability provides a perfect result as the voltage deviations are nearly zero throughout the simulation period. Based on the results in Figure 7, the controller u 0 has much better H performance than u 1 . However, in order to achieve such excellent H performance against the stochastic disturbances in the considered EI system, the controllable generators, such as MTs and DGs, would be overused, which makes the controller inapplicable in most practical scenarios. Thus, when there is no strict requirement for the voltage stabilization, the controller obtained from problem (15) would provide a satisfactory solution for the EI system management.
Based on the analysis above, both u 0 and u 1 have an outstanding effect on the voltage regulation. It is clear that with a smaller disturbance attenuation level γ , u 0 has a better performance. Meanwhile, from Table 2, we found that the feedback gain in u 0 is greater. Larger feedback gain usually suggests faster response speed. However, in practical systems, generators, such as MTs, DGs and FCs, sometimes might not be able to accept control signals from u 0 as they have large values. In many scenarios, a controller with proper feedback gain would be preferred. To show the differences between the characteristics of the generators under the controllers u 0 and u 1 , the power dynamics of MTs in M G 1 are given in Figure 8.
From the power dynamic curves in Figure 8, we found that when u 0 is applied to the EI system, the generators are able to respond to the voltage deviations faster. Thus, there is vast rapid fluctuations in the curves of MT power output corresponding to u 0 . On the other hand, when the controller u 1 is employed, the output power deviations of generators are more moderate. Thus, we can conclude that controller u 1 would be more suitable when the generators are not able to cope with the fast variation of control signals.
In this section, the effectiveness of the proposed control approaches is evaluated with numerical simulations. The performances of the proposed controllers are analyzed. Based on the comparison results, u 0 can be utilized in the scenarios with high voltage stability requirements, while it is more suitable to apply u 1 in normal EI scenarios.

5. Conclusions

In summary, we transform an industrial EI system into a class of continuous linear SDEs. Both time-varying delays and parameter uncertainties are taken into account. The voltage stabilization issue is formulated as a stochastic robust H control problem and is sufficiently solved. Furthermore, to avoid the potential over-control situations, we conducted an extensive analysis that is focused on designing constraints for the controllers. Simulation results shows the effectiveness of our proposed method.
In this paper, the LMI approach has been adopted to solve a robust H control problem. It appears that the integral quadratic constraint or the Riccati equation based approach could also solve this problem efficiently, such as [31,32]. All these methods would achieve similar stabilization targets. However, due to their different assumptions, the proposed LMI based method would be applicable in some circumstances while the other ones would not and vice versa. Therefore, in practical applications, these methods should be utilized accordingly based on the specific scenarios.
It is important to note that the studied dynamical systems are linear, which is relatively conservative. In fact, all systems in real engineering practice are nonlinear and the technique of linearization itself requires extensive approximations, which creates modelling errors. Furthermore, the centralized state feedback controller requires the access to all system states, which is difficult to achieve in many practical systems. In the future, instead of explicitly establishing dynamical models for the EI system, model-free approaches, such as deep reinforcement learning methods [33], shall be implemented to improve the scalability of conventional centralized control approaches.

Author Contributions

Conceptualization, H.H., Y.Q. and J.C.; Formal analysis, H.H.; Funding acquisition, J.C.; Investigation, Y.Q., H.X. and C.H.; Methodology, H.H.; Writing—original draft, H.H., Y.Q., H.X. and C.H.

Funding

This work was supported in part by National Natural Science Foundation of China (grant No. 61472200) and Beijing Municipal Science & Technology Commission (grant No. Z161100000416004).

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

BESBattery energy storage
DERDistributed energy resources
DGDiesel engine generator
EREnergy router
EIEnergy Internet
FCFuel cell
FESFlywheel energy storage
LMILinear matrix inequality
MTMicro-turbine
MGMicrogrid
ODEOrdinary differential equation
PVPhotovoltaic panel
SDEStochastic differential equation
WTWind turbine generator
i Index of MGs
T P V Time constant of PVs
T W T Time constant of WTs
T l Time constant of loads
T M T Time constant of MTs
T D G Time constant of DGs
T F C Time constant of FCs
T B E S Time constant of BESs
T F E S Time constant of FESs
T E R Time constant of ERs
T V Time   constant   of   Δ V
Δ P P V Output power deviation of PVs
Δ P W Output power deviation of WTs
Δ P l Power deviation of loads
Δ P M T Output power deviation of MTs
Δ P D G Output power deviation of DGs
Δ P F C Power deviation of FCs
P E R Power transmitted via ERs.
Δ V DC bus voltage deviation
u M T Control input signal for MTs
u D G Control input signal for DGs
u F C Control input signal for FCs
v φ Solar irradiation power change
v w i n d Wind power change
v l o a d Change of electricity usage
v E R Disturbance from the interconnected ERs
R Adjacency matrix of the ER network
n n -dimensional Euclidean space

References

  1. Olivares, D.; Mehrizi-Sani, A.; Etemadi, A.; Canizares, C.; Iravani, R.; Kazerani, M.; Hajimiragha, A.H.; Gomis-Bellmunt, O.; Saeedifard, M.; Palma-Behnke, R.; et al. Trends in Microgrid Control. IEEE Trans. Smart Grid 2014, 5, 1905–1919. [Google Scholar] [CrossRef]
  2. Pascual, J.; Barricarte, J.; Sanchis, P.; Marroyo, L. Energy management strategy for a renewable-based residential microgrid with generation and demand forecasting. Appl. Energy 2015, 158, 12–25. [Google Scholar] [CrossRef]
  3. Han, H.; Hou, X.; Yang, J.; Wu, J.; Su, M.; Guerrero, J.M. Review of power sharing control strategies for islanding operation of AC microgrids. IEEE Trans. Smart Grid 2016, 7, 200–215. [Google Scholar] [CrossRef]
  4. Ye, M.; Hu, G. Distributed extremum seeking for constrained networked optimization and its application to energy consumption control in smart grid. IEEE Trans. Control Syst. Technol. 2016, 24, 2048–2058. [Google Scholar] [CrossRef]
  5. Rifkin, J. The Third Industrial Revolution: How Lateral Power Is Transforming Energy, the Economy and the World; Palgrave Macmillan: New York, NY, USA, 2013; pp. 31–46. [Google Scholar]
  6. Cao, J.; Hua, H.; Ren, G. Energy Use and the Internet the SAGE Encyclopedia of the Internet; Sage: Newbury Park, CA, USA, 2018; pp. 344–350. [Google Scholar]
  7. Xu, Y.; Zhang, J.H.; Wang, W.Y. Energy router: Architectures and functionalities toward energy internet. In Proceedings of the 2011 IEEE International Conference on Smart Grid Communication, Brussels, Belgium, 17–20 October 2011; pp. 31–36. [Google Scholar]
  8. Hua, H.; Qin, Y.; Hao, C.; Cao, J. Optimal energy management strategies for energy internet via deep reinforcement learning approach. Appl. Energy 2019, 239, 598–609. [Google Scholar] [CrossRef]
  9. Elsayed, A.T.; Mohamed, A.A.; Mohammed, O.A. DC microgrids and distribution systems: An overview. Electr. Power Syst. Res. 2015, 119, 407–417. [Google Scholar] [CrossRef]
  10. Hua, H.; Qin, Y.; Hao, C.; Cao, J. Stochastic optimal control for energy internet: A bottom-up energy management approach. IEEE Trans. Ind. Inform. 2019, 15, 1788–1797. [Google Scholar] [CrossRef]
  11. Zhu, H.; Liu, H.J. Fast local voltage control under limited reactive power: Optimality and stability analysis. IEEE Trans. Power Syst. 2016, 31, 3794–3803. [Google Scholar] [CrossRef]
  12. Cavraro, G.; Carli, R. Local and distributed voltage control algorithms in distribution network. IEEE Trans. Power Syst. 2018, 33, 1420–1430. [Google Scholar] [CrossRef]
  13. Khorsandi, A.; Ashourloo, M.; Mokhtari, H.; Iravani, R. Automatic droop control for a low voltage dc micro-grid. IET Gener. Transm. Distrib. 2016, 10, 41–47. [Google Scholar] [CrossRef]
  14. Simpson, J.; Dorfler, F.; Bullo, F. Voltage stabilization in micro-grids via quadratic droop control. IEEE Trans. Autom. Control 2017, 62, 1239–1253. [Google Scholar] [CrossRef]
  15. Hua, H.; Cao, J.; Yang, G.; Ren, G. Voltage control for uncertain stochastic nonlinear system with application to energy internet: Non-fragile robust H approach. J. Math. Anal. Appl. 2018, 463, 93–110. [Google Scholar] [CrossRef]
  16. Akinyele, D.O.; Rayudu, R.K. Review of energy storage technologies for sustainable power networks. Sustain. Energy Technol. Assess. 2014, 8, 74–91. [Google Scholar] [CrossRef]
  17. Arceo, J.C.; Sanchez, M.; Estrada-Manzo, V.; Bernal, M. Convex stability analysis of nonlinear singular systems via linear matrix inequalities. IEEE Trans. Autom. Control 2018, 64, 1740–1745. [Google Scholar] [CrossRef]
  18. Bevrani, H.; Feizi, M.R.; Ataee, S. Robust frequency control in an islanded microgrid: H and μ-synthesis approaches. IEEE Trans. Smart Grid 2016, 7, 706–717. [Google Scholar] [CrossRef]
  19. Odun-Ayo, T.; Crow, M.L. Structure-preserved power system transient stability using stochastic energy functions. IEEE Trans. Power Syst. 2012, 27, 1450–1458. [Google Scholar] [CrossRef]
  20. Qin, Y.; Hua, H.; Cao, J. Stochastic optimal control scheme for battery lifetime extension in islanded microgrid via a novel modeling approach. IEEE Trans. Smart Grid 2018. [Google Scholar] [CrossRef]
  21. Heymann, B.; Bonnans, J.F.; Silva, F.; Jimenez, G. A stochastic continuous time model for microgrid energy management. In Proceedings of the 2016 European Control Conference, Aalborg, Denmark, 29 June–1 July 2016; pp. 2084–2089. [Google Scholar]
  22. Wen, G.; Hu, G.; Hu, J.; Shi, X.; Chen, G. Frequency regulation of source-grid-load systems: A compound control strategy. IEEE Trans. Ind. Inform. 2016, 12, 69–78. [Google Scholar]
  23. Khodaei, A.; Bahramirad, S.; Shahidehpour, M. Microgrid planning under uncertainty. IEEE Trans. Power. Syst. 2015, 30, 2417–2425. [Google Scholar] [CrossRef]
  24. Aghaei, J.; Alizadeh, M.I. Demand response in smart electricity grids equipped with renewable energy sources: A review. Renew. Sustain. Energy Rev. 2013, 18, 64–72. [Google Scholar] [CrossRef]
  25. Chaouachi, A.; Kamel, R.M.; Andoulsi, R.; Nagasaka, K. Multiobjective intelligent energy management for a microgrid. IEEE Trans. Ind. Electron. 2013, 60, 1688–1699. [Google Scholar] [CrossRef]
  26. Chaudhuri, B.; Majumder, R.; Pal, B. Wide-area measurement-based stabilizing control of power system considering signal transmission delay. IEEE Trans. Power. Syst. 2004, 19, 1971–1979. [Google Scholar] [CrossRef]
  27. Cao, J.; Wan, Y.; Hua, H.; Qin, Y. Delay analysis for end-to-end synchronous communication in monitoring systems. Sensors 2018, 18, 3615. [Google Scholar] [CrossRef] [PubMed]
  28. Huang, Y.; Yuan, X.; Hu, J.; Zhou, P.; Wang, D. DC-Bus voltage control stability affected by AC-bus voltage control in VSCs connected to weak AC grids. IEEE J. Emerg. Sel. Top. Power Electron. 2016, 4, 445–458. [Google Scholar] [CrossRef]
  29. Xu, S.; Chen, T. Robust H control for uncertain stochastic systems with state delay. IEEE Trans. Autom. Control 2002, 47, 2089–2094. [Google Scholar]
  30. Grant, M.; Boyd, S.; Ye, Y. CVX: Matlab software for disciplined convex programming. CVX 2008, 84, 155–210. [Google Scholar]
  31. Savkin, A.V.; Petersen, I.R. Robust H-infinity control of uncertain linear systems with structured uncertainty. J. Math. Syst. Estim. Control. 1996, 6, 339–342. [Google Scholar]
  32. Milliken, P.; Marsh, C.; Van Brunt, B. Minimax controller design for a class of uncertain linear systems. Automatica 1999, 35, 583–590. [Google Scholar] [CrossRef]
  33. Mnih, V.; Badia, A.P.; Mirza, M.; Graves, A.; Lillicrap, T.; Harley, T.; Lillicrap, T.P.; Silver, D.; Kavukcuoglu, K. Asynchronous methods for deep reinforcement learning. In Proceedings of the International Conference on Machine Learning, New York, NY, USA, 19–24 June 2016; pp. 1928–1937. [Google Scholar]
Figure 1. The studied energy internet (EI) scenario.
Figure 1. The studied energy internet (EI) scenario.
Energies 12 01622 g001
Figure 2. Generalized H control configuration for EI systems.
Figure 2. Generalized H control configuration for EI systems.
Energies 12 01622 g002
Figure 3. Time delay τ ( t ) .
Figure 3. Time delay τ ( t ) .
Energies 12 01622 g003
Figure 4. Factor F ( t ) for parameter uncertainties.
Figure 4. Factor F ( t ) for parameter uncertainties.
Energies 12 01622 g004
Figure 5. Power deviations of PVs, WTs and loads in M G 1 .
Figure 5. Power deviations of PVs, WTs and loads in M G 1 .
Energies 12 01622 g005
Figure 6. Voltage deviations in four microgrids (MGs).
Figure 6. Voltage deviations in four microgrids (MGs).
Energies 12 01622 g006
Figure 7. Norms for H bound γ u 0 || v || and controlled output z .
Figure 7. Norms for H bound γ u 0 || v || and controlled output z .
Energies 12 01622 g007
Figure 8. Power deviations of micro-turbines (MTs) under controllers u 0 and u 1 .
Figure 8. Power deviations of micro-turbines (MTs) under controllers u 0 and u 1 .
Energies 12 01622 g008
Table 1. Typical parameters for devices in energy internet (EI) system.
Table 1. Typical parameters for devices in energy internet (EI) system.
Parameter ValueParameter ValueParameter Value
T P V 1.9 T W 1.8 T l 1.6
e P V 0.52 e W 0.4 e l 0.6
T M T 0.3 T D G 0.4 T F C 0.5
T B E S 0.1 T F E S 0.14 T V 1.2
T E R 0.3 a d E R 0.05--
Table 2. Comparison of solutions for proposed controllers.
Table 2. Comparison of solutions for proposed controllers.
Pivotal Parameters u 0 u 1
γ 0.02480.133168
ε 161.0745664.940192
|| Y X 1 || 2 1743.91980927.047299

Share and Cite

MDPI and ACS Style

Hua, H.; Qin, Y.; Xu, H.; Hao, C.; Cao, J. Robust Control Method for DC Microgrids and Energy Routers to Improve Voltage Stability in Energy Internet. Energies 2019, 12, 1622. https://doi.org/10.3390/en12091622

AMA Style

Hua H, Qin Y, Xu H, Hao C, Cao J. Robust Control Method for DC Microgrids and Energy Routers to Improve Voltage Stability in Energy Internet. Energies. 2019; 12(9):1622. https://doi.org/10.3390/en12091622

Chicago/Turabian Style

Hua, Haochen, Yuchao Qin, Hanxuan Xu, Chuantong Hao, and Junwei Cao. 2019. "Robust Control Method for DC Microgrids and Energy Routers to Improve Voltage Stability in Energy Internet" Energies 12, no. 9: 1622. https://doi.org/10.3390/en12091622

APA Style

Hua, H., Qin, Y., Xu, H., Hao, C., & Cao, J. (2019). Robust Control Method for DC Microgrids and Energy Routers to Improve Voltage Stability in Energy Internet. Energies, 12(9), 1622. https://doi.org/10.3390/en12091622

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