# Online Hierarchical Controller for Hybrid Power System

Authors: Salem Zerkaoui
Journal: ISRN Renewable Energy
 (2012)

---

## Abstract

This paper presents the basis for the development of an intelligent and autonomous energy management strategy for hybrid power system (HPS). Two hierarchical levels are proposed to control and manage the HPS. The low level is performed by a local control unit (DC-DC converters controller) of the different power sources. Dynamic equations describing the coupling of converters are derived, and a robust sliding mode dynamic controller is designed. The high level is performed by the online supervisor unit. This unit is designed by applying on-line Takagi-Sugeno fuzzy logic principles. As a result the robust control system gets rid of the limits of the HPS, which has the imprecision, uncertainty, strong coupling, and nonlinearity, to achieve its tractability, robustness, and low solution cost. Under the operation constraints related to each type of sources, the simulation results show that the optimal operation objective of HPS has been achieved.


---

## Body

## 1. Introduction

The limited reserves of fuel oils, their pollution impact, and their unstable prices have significantly increased the interest in renewable energy sources (RES: photovoltaic modules, wind turbine, etc.) to produce electricity that is an essential factor for the development of the human societies. In this context, HPSs, which combine renewable energy and conventional energy sources with the storage systems, are very interesting in terms of environmental protection and reduction of the effects of greenhouse gas emissions. In other words, the HPSs can provide an economic, environment friendly, and reliable supply of electricity.To perform the hybrid systems a control strategy has to be designed and implemented on the system. Numerous solution methodologies have been proposed in this field in the last decade [1–4]. Many of them are based on a static approach. However, the development of robust controller is necessary to ensure stability and robustness of the multisources of the renewable energy systems.Subsequently, others advanced control strategies, such as predictive control, fuzzy logic, and neural network, have been developed and successfully applied [5–7]. However, these approaches require considerable computing resources and as a result their applicability for real-time applications is reduced.In this paper, emphases are put on the energy management and control from the viewpoint of control theory. The proposed management strategy and controller, designed by applying fuzzy logic and sliding mode principles, not only is simple, stable, and robust, but also reduces the computational resources. Two hierarchical levels are required to control and manage the HPS. The low level is performed by a local control unit (DC-DC converters controller) of the different power sources. This level manipulates the duty cycle converter according to the variation of the sources and the load. Without the use of cascade structure commonly used in literatures, I proposed a robust control law based on sliding mode control approach to operate the DC-DC converters in the most efficient way.The high level is performed by the online supervisor unit. This unit uses the data about the load, the meteorological conditions and the charge state of energy storage system (ESS) and combustion engine (CE), to correctly and efficiently share the load demand according to the availability of conventional and renewable energy, in other words, to decide whether to charge or discharge the ESS, to turn on or off the CE, to reduce the renewable sources power production or not, and so forth. Such control level could be achieved by adaptive techniques control such as those employing fuzzy logic.The paper is organized in five sections. In Section2, the hybrid power system and their mathematical models are described. Then, in Section 3, I will present the design of low level controller witch performed by a local sliding mode control unit. By applying online dynamic Takagi-Sugeno fuzzy logic principles, the high-level controller design is presented in Section 4. An application to the multisource system is presented in Section 5 in order to illustrate the efficiency of the developed approach. The conclusion sums up the main characteristic of the proposed hierarchical controller and introduces one of the most challenge for adaptive control in the future research.


## 2.  System Configuration

The HPS under consideration is shown in Figure1. It consists of a conventional generation (CE), a renewable energy sources (RES: photovoltaic module and wind turbine), an ESS (batteries and ultracapacitors), and a variable load. All these elements are connected onto a DC bus through DC/DC power electronic converters.Figure 1
HPS Scheme.In order to allow ESS charging as well as discharging, the DC/DC converter which connects the ESS to the DC bus is bidirectional unlike what other converters are unidirectional. The DC bus accumulates the generated energy and sends it to the variable load and, if necessary, to the energy storage system. In this configuration, renewable sources take over as main energy source.


### 2.1. Models of Energy Sources

#### 2.1.1.  Renewable Sources

The most accessible energy sources are solar and wind energy which can be produced through photovoltaic conversion and wind turbine, respectively. The available photovoltaic (PV) and wind turbine energy production are calculated based on the meteorological information (irradiation, temperature, and wind speed) and their characteristics. In this paper, the photovoltaic and wind turbine power models were obtained, respectively, from [8, 9] as follows.For PV area, the electric model closest to the PV generator is a model with one diode.The current output of a PV array ((Np-Ns) PV cells interconnected in a parallel series) Ip can be calculated by
(1)Ip=Np[Iph-Is*(exp(NpVp+NsIpRsNpNsVT)-1)]-NpVp+NsIpRsNsRsh,
where VT=(nKBT)/q is the thermodynamic voltage, Vp is the voltage level on the PV panel array terminals, Iph is the photocurrent, and Is the reverse saturation current, KB is the Boltzman constant, and n is the cell deviation from the ideal p-n junction characteristic.Based on the observation of voltage-power relationship, I can notice that the PV array has a unique voltage point which provides the maximum power [10]. This maximum power point (MPP) can be tracked in real time with a control algorithm called MPPT in order to increase their performances. Several algorithms have been proposed in existing literature [11] such as: perturb and observe, incremental conductance, fractional open circuit voltage, fuzzy control, and neural network. As the objective of this paper is to develop the higher-level control algorithms (energy management strategy), I adopt in the simulation section one of these algorithm (e.g., incremental conductance).For the wind turbine, the output power of wind generator at wind speedv is expressed as
(2)PWT(t)={(v3-vci3)vr3-vci3PRvci<v<vr,PRvr<v<vco,0v≤vciorv≥vco,v(t)=vr(t)·(hhr)γ,
where PR is the nominal power; v, vr, vci, and vco, respectively, are the wind speed at desired height h, wind speed at the reference height hr, and cut-in and cutoff wind speed of the wind generator.


#### 2.1.2. Combustion Engine

In the case of a renewable energy system, the production of electrical energy is conducted within the resources (wind, sun, …) and not in demand. For standalone installations, in order to assist the renewable sources in situations of deficient environmental circumstances, it is necessary to resort to storage and/or to add a CE.The CE is a conventional power source which is independent from environmental conditions (irradiation, wind velocity, etc.). This source is the site of a complex phenomena and strongly nonlinear. There are several levels, depending on the desired use and degree of complexity of models. In this work, for simplification reasons of system study, I assumed that the dynamic behavior of CE can be approximated by a first-order system with a time constantτ [12].


#### 2.1.3. Energy Storage Systems

Energy storage is a key factor in a HPS. It can provide energy when needed and store the surplus of production otherwise. Since there is a wide range of storage options, for a given application, both the power and energy rating of the storage device must be considered in the device selection. The most common storage technologies can be summarized in terms of transient and steady-state storage devices. In our case, the energy storage system comprises the batteries and the supercapacitors banks. The combination of these sources makes the best use of the advantages of each individual device and may meet the requirements of the load regarding both power and energy densities.(1) Batteries
In general, the batteries allow the good storage characteristics for wide ranges of applications such as electric vehicles, hybrid electric vehicles, and electrical networks [13–15]. The mathematical dynamic model of the battery bank used in the simulation program is deducted from the CIEMAT model [16, 17]. Its input is the current battery profile which is positive when power is drawn from the batteries and negative when power is supplied to the batteries. Its output is the state of charge (SOC), which defines the relationship between the battery’s total energy and its usable energy (3). We have
(3)SOCmin≤SOC(t)=SOC(t-1)+IbatCbat≤SOCmax,
where Ibat is the charge (or discharge) instantaneous current of the battery bank; this current is positive when the batteries are charging and negative otherwise. Cbat is the capacity of battery bank calculated from the average battery current Ib_m as follow:
(4)Cbat=1.67*C101+0.67*(Ib_m/I10)0.9*(1+0.005*ΔT),
where ΔT is a battery bank heating relative to the ambient temperature, C10 is a nominal capacity of battery bank, and I10 is a nominal current of battery bank.
According to the specifications from the manufacturers, in this paper, the battery’s lifetime can be prolonged to the maximum if SOC is between 20 and 90% [17].(2) Supercapacitors
Supercapacitors are electrical storage devices with a high power density, used for peak power leveling to assist the battery during transients.
In contrast to a battery, the voltage of a supercapacitor changes linearly with its state of charge as follows [18]:
(5)SOCUC=43(CV22Emax-14),
where C is the equivalent total capacity of supercapacitors bank, V is the instantaneous voltage measured at its terminals, and Emax is the maximum stored energy obtained for Vmax.
In this paper, to ensure good performances and a sufficient lifetime of supercapacitors and to avoid overcharging the maximum voltage was limited toVmax, and the discharge rate is fixed at 0.25.
When all these constraints due to limitations in voltage, current and power are combined, the momentary power output of the supercapacitors can be calculated for every time step.


### 2.2. Operating Modes of Energy Sources

Because of the variable structure nature of multisource system, the operating modes of each source are summarized in Figure2.Figure 2
Multisource operating modes.For renewable sources, there are two possible operating modes: Maxi Mode and Limit Mode. The Maxi Mode occurs when the maximum energy captured from the renewable sources is not enough to satisfy the total power demand (the load plus the ESS recharge requirements). In this case, the MPPT algorithm must be used (i.e., capturing the maximum renewable power (Pren_max) available). However, when the provided power is higher than demanded one (Limit Mode), only the required amount of energy is produced.For CE, there are three possible operating modes: Production Mode (PM), Off Mode (OM), and Slow Down Mode (SDM).In PM, the CE assists the renewable sources in order to provide the load power demand and the ESS power recharge.In OM, the load power demand is provided with renewable sources assisted by ESS, and the CE is off.In SDM, the load power demand is only provided with renewable sources, and the CE is relaxed to avoid transient phase and overheating. This is important to avoid running the CE for very short period of time.For ESS, depending on the flow of power, there are two possible operating modes, Help Mode (HM) and Recharge Mode (RM).In the HM, the power demand is provided by the renewable sources assisted by ESS while the CE is switched off and renewable sources power is lower then demanded one.In the RM, when the load energy demand is smaller than a produced power, the redundant energy is used to charge the ESS.According to profile of the redundant energy (Figure3), I can distinguish two charge/discharge cycles: steady-state cycle ensured by batteries and transient-state cycle ensured by ultracapacitors, which are used for peak power leveling to assist the batteries during hard power demand.Figure 3
Energy storage management strategy.In order to improve system operation, the power and energy management strategy has to decide which operating mode must be used.


## 2.1. Models of Energy Sources

### 2.1.1.  Renewable Sources

The most accessible energy sources are solar and wind energy which can be produced through photovoltaic conversion and wind turbine, respectively. The available photovoltaic (PV) and wind turbine energy production are calculated based on the meteorological information (irradiation, temperature, and wind speed) and their characteristics. In this paper, the photovoltaic and wind turbine power models were obtained, respectively, from [8, 9] as follows.For PV area, the electric model closest to the PV generator is a model with one diode.The current output of a PV array ((Np-Ns) PV cells interconnected in a parallel series) Ip can be calculated by
(1)Ip=Np[Iph-Is*(exp(NpVp+NsIpRsNpNsVT)-1)]-NpVp+NsIpRsNsRsh,
where VT=(nKBT)/q is the thermodynamic voltage, Vp is the voltage level on the PV panel array terminals, Iph is the photocurrent, and Is the reverse saturation current, KB is the Boltzman constant, and n is the cell deviation from the ideal p-n junction characteristic.Based on the observation of voltage-power relationship, I can notice that the PV array has a unique voltage point which provides the maximum power [10]. This maximum power point (MPP) can be tracked in real time with a control algorithm called MPPT in order to increase their performances. Several algorithms have been proposed in existing literature [11] such as: perturb and observe, incremental conductance, fractional open circuit voltage, fuzzy control, and neural network. As the objective of this paper is to develop the higher-level control algorithms (energy management strategy), I adopt in the simulation section one of these algorithm (e.g., incremental conductance).For the wind turbine, the output power of wind generator at wind speedv is expressed as
(2)PWT(t)={(v3-vci3)vr3-vci3PRvci<v<vr,PRvr<v<vco,0v≤vciorv≥vco,v(t)=vr(t)·(hhr)γ,
where PR is the nominal power; v, vr, vci, and vco, respectively, are the wind speed at desired height h, wind speed at the reference height hr, and cut-in and cutoff wind speed of the wind generator.


### 2.1.2. Combustion Engine

In the case of a renewable energy system, the production of electrical energy is conducted within the resources (wind, sun, …) and not in demand. For standalone installations, in order to assist the renewable sources in situations of deficient environmental circumstances, it is necessary to resort to storage and/or to add a CE.The CE is a conventional power source which is independent from environmental conditions (irradiation, wind velocity, etc.). This source is the site of a complex phenomena and strongly nonlinear. There are several levels, depending on the desired use and degree of complexity of models. In this work, for simplification reasons of system study, I assumed that the dynamic behavior of CE can be approximated by a first-order system with a time constantτ [12].


### 2.1.3. Energy Storage Systems

Energy storage is a key factor in a HPS. It can provide energy when needed and store the surplus of production otherwise. Since there is a wide range of storage options, for a given application, both the power and energy rating of the storage device must be considered in the device selection. The most common storage technologies can be summarized in terms of transient and steady-state storage devices. In our case, the energy storage system comprises the batteries and the supercapacitors banks. The combination of these sources makes the best use of the advantages of each individual device and may meet the requirements of the load regarding both power and energy densities.(1) Batteries
In general, the batteries allow the good storage characteristics for wide ranges of applications such as electric vehicles, hybrid electric vehicles, and electrical networks [13–15]. The mathematical dynamic model of the battery bank used in the simulation program is deducted from the CIEMAT model [16, 17]. Its input is the current battery profile which is positive when power is drawn from the batteries and negative when power is supplied to the batteries. Its output is the state of charge (SOC), which defines the relationship between the battery’s total energy and its usable energy (3). We have
(3)SOCmin≤SOC(t)=SOC(t-1)+IbatCbat≤SOCmax,
where Ibat is the charge (or discharge) instantaneous current of the battery bank; this current is positive when the batteries are charging and negative otherwise. Cbat is the capacity of battery bank calculated from the average battery current Ib_m as follow:
(4)Cbat=1.67*C101+0.67*(Ib_m/I10)0.9*(1+0.005*ΔT),
where ΔT is a battery bank heating relative to the ambient temperature, C10 is a nominal capacity of battery bank, and I10 is a nominal current of battery bank.
According to the specifications from the manufacturers, in this paper, the battery’s lifetime can be prolonged to the maximum if SOC is between 20 and 90% [17].(2) Supercapacitors
Supercapacitors are electrical storage devices with a high power density, used for peak power leveling to assist the battery during transients.
In contrast to a battery, the voltage of a supercapacitor changes linearly with its state of charge as follows [18]:
(5)SOCUC=43(CV22Emax-14),
where C is the equivalent total capacity of supercapacitors bank, V is the instantaneous voltage measured at its terminals, and Emax is the maximum stored energy obtained for Vmax.
In this paper, to ensure good performances and a sufficient lifetime of supercapacitors and to avoid overcharging the maximum voltage was limited toVmax, and the discharge rate is fixed at 0.25.
When all these constraints due to limitations in voltage, current and power are combined, the momentary power output of the supercapacitors can be calculated for every time step.


## 2.1.1.  Renewable Sources

The most accessible energy sources are solar and wind energy which can be produced through photovoltaic conversion and wind turbine, respectively. The available photovoltaic (PV) and wind turbine energy production are calculated based on the meteorological information (irradiation, temperature, and wind speed) and their characteristics. In this paper, the photovoltaic and wind turbine power models were obtained, respectively, from [8, 9] as follows.For PV area, the electric model closest to the PV generator is a model with one diode.The current output of a PV array ((Np-Ns) PV cells interconnected in a parallel series) Ip can be calculated by
(1)Ip=Np[Iph-Is*(exp(NpVp+NsIpRsNpNsVT)-1)]-NpVp+NsIpRsNsRsh,
where VT=(nKBT)/q is the thermodynamic voltage, Vp is the voltage level on the PV panel array terminals, Iph is the photocurrent, and Is the reverse saturation current, KB is the Boltzman constant, and n is the cell deviation from the ideal p-n junction characteristic.Based on the observation of voltage-power relationship, I can notice that the PV array has a unique voltage point which provides the maximum power [10]. This maximum power point (MPP) can be tracked in real time with a control algorithm called MPPT in order to increase their performances. Several algorithms have been proposed in existing literature [11] such as: perturb and observe, incremental conductance, fractional open circuit voltage, fuzzy control, and neural network. As the objective of this paper is to develop the higher-level control algorithms (energy management strategy), I adopt in the simulation section one of these algorithm (e.g., incremental conductance).For the wind turbine, the output power of wind generator at wind speedv is expressed as
(2)PWT(t)={(v3-vci3)vr3-vci3PRvci<v<vr,PRvr<v<vco,0v≤vciorv≥vco,v(t)=vr(t)·(hhr)γ,
where PR is the nominal power; v, vr, vci, and vco, respectively, are the wind speed at desired height h, wind speed at the reference height hr, and cut-in and cutoff wind speed of the wind generator.


## 2.1.2. Combustion Engine

In the case of a renewable energy system, the production of electrical energy is conducted within the resources (wind, sun, …) and not in demand. For standalone installations, in order to assist the renewable sources in situations of deficient environmental circumstances, it is necessary to resort to storage and/or to add a CE.The CE is a conventional power source which is independent from environmental conditions (irradiation, wind velocity, etc.). This source is the site of a complex phenomena and strongly nonlinear. There are several levels, depending on the desired use and degree of complexity of models. In this work, for simplification reasons of system study, I assumed that the dynamic behavior of CE can be approximated by a first-order system with a time constantτ [12].


## 2.1.3. Energy Storage Systems

Energy storage is a key factor in a HPS. It can provide energy when needed and store the surplus of production otherwise. Since there is a wide range of storage options, for a given application, both the power and energy rating of the storage device must be considered in the device selection. The most common storage technologies can be summarized in terms of transient and steady-state storage devices. In our case, the energy storage system comprises the batteries and the supercapacitors banks. The combination of these sources makes the best use of the advantages of each individual device and may meet the requirements of the load regarding both power and energy densities.(1) Batteries
In general, the batteries allow the good storage characteristics for wide ranges of applications such as electric vehicles, hybrid electric vehicles, and electrical networks [13–15]. The mathematical dynamic model of the battery bank used in the simulation program is deducted from the CIEMAT model [16, 17]. Its input is the current battery profile which is positive when power is drawn from the batteries and negative when power is supplied to the batteries. Its output is the state of charge (SOC), which defines the relationship between the battery’s total energy and its usable energy (3). We have
(3)SOCmin≤SOC(t)=SOC(t-1)+IbatCbat≤SOCmax,
where Ibat is the charge (or discharge) instantaneous current of the battery bank; this current is positive when the batteries are charging and negative otherwise. Cbat is the capacity of battery bank calculated from the average battery current Ib_m as follow:
(4)Cbat=1.67*C101+0.67*(Ib_m/I10)0.9*(1+0.005*ΔT),
where ΔT is a battery bank heating relative to the ambient temperature, C10 is a nominal capacity of battery bank, and I10 is a nominal current of battery bank.
According to the specifications from the manufacturers, in this paper, the battery’s lifetime can be prolonged to the maximum if SOC is between 20 and 90% [17].(2) Supercapacitors
Supercapacitors are electrical storage devices with a high power density, used for peak power leveling to assist the battery during transients.
In contrast to a battery, the voltage of a supercapacitor changes linearly with its state of charge as follows [18]:
(5)SOCUC=43(CV22Emax-14),
where C is the equivalent total capacity of supercapacitors bank, V is the instantaneous voltage measured at its terminals, and Emax is the maximum stored energy obtained for Vmax.
In this paper, to ensure good performances and a sufficient lifetime of supercapacitors and to avoid overcharging the maximum voltage was limited toVmax, and the discharge rate is fixed at 0.25.
When all these constraints due to limitations in voltage, current and power are combined, the momentary power output of the supercapacitors can be calculated for every time step.


## 2.2. Operating Modes of Energy Sources

Because of the variable structure nature of multisource system, the operating modes of each source are summarized in Figure2.Figure 2
Multisource operating modes.For renewable sources, there are two possible operating modes: Maxi Mode and Limit Mode. The Maxi Mode occurs when the maximum energy captured from the renewable sources is not enough to satisfy the total power demand (the load plus the ESS recharge requirements). In this case, the MPPT algorithm must be used (i.e., capturing the maximum renewable power (Pren_max) available). However, when the provided power is higher than demanded one (Limit Mode), only the required amount of energy is produced.For CE, there are three possible operating modes: Production Mode (PM), Off Mode (OM), and Slow Down Mode (SDM).In PM, the CE assists the renewable sources in order to provide the load power demand and the ESS power recharge.In OM, the load power demand is provided with renewable sources assisted by ESS, and the CE is off.In SDM, the load power demand is only provided with renewable sources, and the CE is relaxed to avoid transient phase and overheating. This is important to avoid running the CE for very short period of time.For ESS, depending on the flow of power, there are two possible operating modes, Help Mode (HM) and Recharge Mode (RM).In the HM, the power demand is provided by the renewable sources assisted by ESS while the CE is switched off and renewable sources power is lower then demanded one.In the RM, when the load energy demand is smaller than a produced power, the redundant energy is used to charge the ESS.According to profile of the redundant energy (Figure3), I can distinguish two charge/discharge cycles: steady-state cycle ensured by batteries and transient-state cycle ensured by ultracapacitors, which are used for peak power leveling to assist the batteries during hard power demand.Figure 3
Energy storage management strategy.In order to improve system operation, the power and energy management strategy has to decide which operating mode must be used.


## 3. Low-Level Controller

The DC/DC converters are basic constituents of the multisource system. The aim of these converters is to regulate the DC component of the output voltage to its reference by controlling the current provided by each source in spite of the voltage variations on their inputs.


### 3.1. Multimodel of Several DC/DC Power Converters Coupled on a DC Bus

An accurate model for the coupling of DC/DC power converters on a DC bus has been developed previously in [19].The ZVS full-bridge isolated DC/DC converter [20, 21], studied in this paper, is represented on Figure 4.Figure 4
Structural diagram of the ZVS full bridge isolated DC/DC converter.It supposed that it runs in a continuous conduction mode.Four basic structures of operation may be distinguished.Phase  1: (Q1Q4) On and (Q2Q3) Off.Phase  2: (Q1d3) On.Phase  3: (Q2Q3) On and (Q1Q4) Off.Phase  4: (Q2d4) On.The structural diagram corresponding to the operation phases 1 and 3 is depicted in Figure5 and the one corresponding to the operation phases 2 and 4 is represented on Figure 6. rp and rs are, respectively, the transformer primary and secondary resistances. LM is the transformer magnetizing inductance. CDC is the DC bus capacity, and the load Rch is supposed resistive.Figure 5
Phases 1 and 3.Figure 6
Phases 2 and 4.The voltage source (Vg) is disconnected during phases 2 and 4 which correspond to the transformer demagnetization (Figure 4), so ig=0.Each phase leads to the following state space model:(6)x˙=Aix+Biu,
where x=[iMiLVbus]T is the state vector, and u=[VgVd]T is the control vector.i
M, iL, Vbus, vg, and vd, respectively, represent the transformer magnetizing current, the inductance current, the DC/DC power converters output voltage, the source, and the diode voltage.It is assumed that diodes are not ideal. We have
(7)A1=[-(2Rmos+rp)LM-n(2Rmos+rp)LM0-n(2Rmos+rp)L-n2(2Rmos+rp)-(rs+rL)L-1L01Ceq-1ReqCeq],A3=[-(2Rmos+rp)LM+n(2Rmos+rp)LM0-n(2Rmos+rp)L-n2(2Rmos+rp)-(rs+rL)L-1L01Ceq-1ReqCeq],A2/4=[-(Rmos+rp)LM000-rLL-1L01Ceq-1ReqCeq],B1=[1LM0nL-2L00]T,B3=[-1LM0nL-2L00]T,B2/4=[000000]T.Combining the above phase states model, the average model of the DC/DC converter shown in Figure4 can be derived as follows [22, 23]:
(8)x˙(t)=f(x,u)+g(x,u)φ(t),
where
(9)f(x,u)=ax(t)+bu(t),(10)g(x,u)=Ax(t)+Bu(t),(11)A=12(A1-A2+A3-A4),(12)a=12(A2+A4),(13)B=12(B1-B2+B3-B4),(14)b=12(B2+B4).The average model of the DC/DC converter developed above can be extended to the coupling of many renewable sources [24]. I suppose that the DC/DC power converters are identical.The partial structural diagram of the coupling ofm DC/DC power converters is represented in Figure 7.Figure 7
Structural diagram of the coupling ofm DC/DC power converters.The multimodel form DC/DC converters coupled on a DC bus is express as follows:
(15)X˙(t)=F(X,u)+G(X,u)φ(t),
where
(16)F(X,U)=[f1(x1,u1),…,fm(xm,um)]T,G(X,U)=diag[g1(x1,u1),…,gm(xm,um)],X=[x1,…,xm]T,U=[u1,…,um]T,φ=[φ1,…,φm]T.The multimodel parameters (F and G) change according to the DC/DC power converters coupled on the DC bus.For simplicity of presentation, in this subsection, only two coupling sources are considered (m=2).


### 3.2. Dynamic Controller Design for DC/DC Converter

The variable structure systems (VSS) theory [25] can be extremely helpful in the study of the control of powers DC/DC converters. The switched mode DC/DC converters are nonlinear, and it is not suitable to the application of linear control theory.A different approach, which complies with the nonlinear nature of switch-mode power supplies, is proposed. In this context, the sliding mode control, which is derived from the VSS theory, appears as a powerful control technique that offers several advantages: stability even for large supply and load variations, robustness, good dynamic response and simple implementation. Their capabilities emerge especially in application to high-order converters, yielding improved performances as compared to usual control techniques.In this subsection, I design a sliding mode dynamic controller “SMDC” which regulates the voltage level on the DC bus by controlling the current provided by each source.In various works in the literature, a cascade control structure with two control loops for DC/DC converters control is usually adopted [26, 27]. In this structure, an inner current loop regulates the inductor current, whereas an external control loop keeps a constant output voltage.However, this structure has some drawbacks such as the value of reference current that may be poorly estimated because of quick variations of the load and the sources. The current estimation error introduces a bus voltage error compared to the reference one. This affects the controller performances in terms of robustness.To overcome this issue, I adopt, in this paper, a MIMO sliding mode robust control scheme (Figure8).Figure 8
MIMO sliding mode robust control scheme.The development of the sliding mode control scheme consists of two phases. The first is to design a sliding surface where the DC/DC converter exhibits desired properties. The second is to design a control law to drive and maintain the system on the sliding surface [26, 27].In order to act simultaneously on all state variables, let us define the following PI-type sliding surface:(17)δi(t)=ei(t)+α∫0tei(τ)dτ,i=1,2,
where
(18)ei(t)=ki(ILi-ILi_ref)+k3(Vbus-Vbus_ref).I
Li
_
ref, Vbus_ref are the references of  inductor current and voltage DC bus. k1,2,3 are the tracking weight factor error, and α>0 is the sliding-surface integral parameter.In the following, a bounded control input is designed to forceδi(t) to converge to zero or make its absolute value smaller.We obtain the equivalent controlφeq1,2 by applying the invariance condition:
(19)δi(t)=δ˙i(t)=0,withφi(t)=φeqi(t),i=1,2.From (17), the time derivative of δi(t) along system (15) is given as
(20)δ˙i=ki(A2-xi+B2-ui)φeqi+ki(a2-xi+b2-ui-I˙Li_ref)-k3(A3-xi+B3-ui-V˙bus_ref)+αei=0.The notationA2- refers to the second column of the A matrix.Thus, the equivalent control input is given as(21)φeqi=-RiHi,i=1,2,
where
(22)Ri=ki(a2-xi+b2-ui-I˙Li_ref)-k3(A3-xi+B3-ui-S˙ref)+αei,(23)Hi=ki(A2-xi+B2-ui).It is well known that the exact value of the external disturbance and the parameter variations of the system, such as internal resistance and magnetizing current, are difficult to measure in advance for practical applications. Therefore, I propose a nonlinear switching control-inputφsi to estimate the upper bound of uncertainties and external disturbance:
(24)φsi(t)=-γi(t)Hi(t)sign(δi(t)).In the presence of parameter uncertainties and external disturbances the dynamic equation (15) becomes
(25)X˙(t)=F(x)+G(x)φ+Δ(t),
where Δ includes the uncertainties and perturbations.In order to reduce chattering, the most common method is to replace the sign function by the saturation functionsat(δ(t),ζ), where
(26)sat(δ(t),ζ)={sign(δ),|δ|>ζ>0,δζ,|δ|≤ζ
and ζ is a small positive constant.Hence, the SMDC rule can be designed as(27)φi(t)=φeqi(t)-γi(t)Hi(t)sat(δi(t),ζ).The main theorem, stated and proved below, provides sufficient conditions to ensure the stability and robustness of the multisources system.Theorem 1.
Let the nonlinear multisources system given by (15) and sliding surface given by (17), and suppose that the Lyapunov function is defined by
(28)V(t)=12δ1(t)2+12δ2(t)2.
Then, the sufficient stability condition for the SMDC in the sense of Lyapunov should satisfy the following range of switching gains:(29)γi(t)≥|Δ(t)|,i=1,2,
where |Δ| represent the upper bound of the uncertainties and perturbations.Proof.
IfV(t) is a Lyapunov function candidate defined as in (28), the asymptotic stability will be satisfied if V˙(t)≤0, with
(30)V˙(t)=δ1(t)δ˙1(t)+δ2(t)δ˙2(t).
Using (17), (25), and (27) in (30) I obtain
(31)V˙(t)=δ1(t)[-γ1(t)sat(δ1(t),ζ)+Δ(t)]+δ2(t)[-γ2(t)sat(δ2(t),ζ)+Δ(t)]≤0.As a conclusion, to satisfy the conditionV˙(t)≤0, and to compensate the bounded parametric uncertainties and disturbances, I must restrict to sufficient condition (29).


### 3.3. Evaluation of Low Level Control Design

To illustrate the low level control design and performance evaluation of my structure, a hybrid power system coupled to two sources through a similar DC/DC converter is considered. The characteristics of the full bridge converters are(32)Rmos=0.005Ω,Rl=2Ω,n=12,Vd=0.3v,L=1e-3H,Lm=20e-6H,rp=0.05Ω,rs=0.05Ω,Ce=20e-6F,Re=56e3Ω,CDC=10e-6F.The DC bus voltage reference is set at 100 V. Simulations are obtained with sampling intervalTe=50 μs.We suppose that 40% of the load demands are supplied with the first source and 60% with the second one.Responses obtained with my algorithm are compared with ones resulting from the PI algorithm with cascade structure. The linearization techniques are applied to the converters in order to deduce linear parameters [28].The simulation conditions are identical for both controllers.To compare the global performances of my algorithm with PI one, let us consider the normalised mean square error (NMSE, Table1), that is interpreted as the overall deviations between output plant and desired values and is defined as:
(33)NMSE=∑l(Vbus-Vbus_ref)2∑l(Vbus_ref)2.Table 1
NMSE.
PI
SMDC
1.2483
0.0075In order to test the capacity of SMDC to reject disturbances, and to investigate the robustness of the control scheme, the load and sources dynamics of the hybrid system are supposed fast. The sources and load profiles are represented by Figures9 and 10.Figure 9
Sources profiles.Figure 10
Load profiles.The closed-loop responses of load current DC bus voltage, tracking voltage error, and load current obtaining with PI controller and sliding mode controller are shown in Figures11, 12, and 13, respectively.Figure 11
DC bus voltage plot obtained with PI controller and sliding mode controller.Figure 12
Tracking voltage error plot (logarithmic scale) obtained with PI controller and sliding mode controller.Figure 13
Load current plot obtained with PI controller and sliding mode controller.Given the fast dynamics imposed on the sources and the load as evident from Figures9 and 10, the simulation results reveal that the SMDC is advantageous in providing negligible steady-state errors (NMSE = 0.0075), to adjust the flow source that meets the load demand and guaranteed stable system (Figure 12). It can also be seen that the responses are satisfactory in terms of overshoot, settling time, and fall time. This shows an excellent behavior of the sliding mode controller in comparison with PI controllers which does not run correctly for wide variations of operating point.


## 3.1. Multimodel of Several DC/DC Power Converters Coupled on a DC Bus

An accurate model for the coupling of DC/DC power converters on a DC bus has been developed previously in [19].The ZVS full-bridge isolated DC/DC converter [20, 21], studied in this paper, is represented on Figure 4.Figure 4
Structural diagram of the ZVS full bridge isolated DC/DC converter.It supposed that it runs in a continuous conduction mode.Four basic structures of operation may be distinguished.Phase  1: (Q1Q4) On and (Q2Q3) Off.Phase  2: (Q1d3) On.Phase  3: (Q2Q3) On and (Q1Q4) Off.Phase  4: (Q2d4) On.The structural diagram corresponding to the operation phases 1 and 3 is depicted in Figure5 and the one corresponding to the operation phases 2 and 4 is represented on Figure 6. rp and rs are, respectively, the transformer primary and secondary resistances. LM is the transformer magnetizing inductance. CDC is the DC bus capacity, and the load Rch is supposed resistive.Figure 5
Phases 1 and 3.Figure 6
Phases 2 and 4.The voltage source (Vg) is disconnected during phases 2 and 4 which correspond to the transformer demagnetization (Figure 4), so ig=0.Each phase leads to the following state space model:(6)x˙=Aix+Biu,
where x=[iMiLVbus]T is the state vector, and u=[VgVd]T is the control vector.i
M, iL, Vbus, vg, and vd, respectively, represent the transformer magnetizing current, the inductance current, the DC/DC power converters output voltage, the source, and the diode voltage.It is assumed that diodes are not ideal. We have
(7)A1=[-(2Rmos+rp)LM-n(2Rmos+rp)LM0-n(2Rmos+rp)L-n2(2Rmos+rp)-(rs+rL)L-1L01Ceq-1ReqCeq],A3=[-(2Rmos+rp)LM+n(2Rmos+rp)LM0-n(2Rmos+rp)L-n2(2Rmos+rp)-(rs+rL)L-1L01Ceq-1ReqCeq],A2/4=[-(Rmos+rp)LM000-rLL-1L01Ceq-1ReqCeq],B1=[1LM0nL-2L00]T,B3=[-1LM0nL-2L00]T,B2/4=[000000]T.Combining the above phase states model, the average model of the DC/DC converter shown in Figure4 can be derived as follows [22, 23]:
(8)x˙(t)=f(x,u)+g(x,u)φ(t),
where
(9)f(x,u)=ax(t)+bu(t),(10)g(x,u)=Ax(t)+Bu(t),(11)A=12(A1-A2+A3-A4),(12)a=12(A2+A4),(13)B=12(B1-B2+B3-B4),(14)b=12(B2+B4).The average model of the DC/DC converter developed above can be extended to the coupling of many renewable sources [24]. I suppose that the DC/DC power converters are identical.The partial structural diagram of the coupling ofm DC/DC power converters is represented in Figure 7.Figure 7
Structural diagram of the coupling ofm DC/DC power converters.The multimodel form DC/DC converters coupled on a DC bus is express as follows:
(15)X˙(t)=F(X,u)+G(X,u)φ(t),
where
(16)F(X,U)=[f1(x1,u1),…,fm(xm,um)]T,G(X,U)=diag[g1(x1,u1),…,gm(xm,um)],X=[x1,…,xm]T,U=[u1,…,um]T,φ=[φ1,…,φm]T.The multimodel parameters (F and G) change according to the DC/DC power converters coupled on the DC bus.For simplicity of presentation, in this subsection, only two coupling sources are considered (m=2).


## 3.2. Dynamic Controller Design for DC/DC Converter

The variable structure systems (VSS) theory [25] can be extremely helpful in the study of the control of powers DC/DC converters. The switched mode DC/DC converters are nonlinear, and it is not suitable to the application of linear control theory.A different approach, which complies with the nonlinear nature of switch-mode power supplies, is proposed. In this context, the sliding mode control, which is derived from the VSS theory, appears as a powerful control technique that offers several advantages: stability even for large supply and load variations, robustness, good dynamic response and simple implementation. Their capabilities emerge especially in application to high-order converters, yielding improved performances as compared to usual control techniques.In this subsection, I design a sliding mode dynamic controller “SMDC” which regulates the voltage level on the DC bus by controlling the current provided by each source.In various works in the literature, a cascade control structure with two control loops for DC/DC converters control is usually adopted [26, 27]. In this structure, an inner current loop regulates the inductor current, whereas an external control loop keeps a constant output voltage.However, this structure has some drawbacks such as the value of reference current that may be poorly estimated because of quick variations of the load and the sources. The current estimation error introduces a bus voltage error compared to the reference one. This affects the controller performances in terms of robustness.To overcome this issue, I adopt, in this paper, a MIMO sliding mode robust control scheme (Figure8).Figure 8
MIMO sliding mode robust control scheme.The development of the sliding mode control scheme consists of two phases. The first is to design a sliding surface where the DC/DC converter exhibits desired properties. The second is to design a control law to drive and maintain the system on the sliding surface [26, 27].In order to act simultaneously on all state variables, let us define the following PI-type sliding surface:(17)δi(t)=ei(t)+α∫0tei(τ)dτ,i=1,2,
where
(18)ei(t)=ki(ILi-ILi_ref)+k3(Vbus-Vbus_ref).I
Li
_
ref, Vbus_ref are the references of  inductor current and voltage DC bus. k1,2,3 are the tracking weight factor error, and α>0 is the sliding-surface integral parameter.In the following, a bounded control input is designed to forceδi(t) to converge to zero or make its absolute value smaller.We obtain the equivalent controlφeq1,2 by applying the invariance condition:
(19)δi(t)=δ˙i(t)=0,withφi(t)=φeqi(t),i=1,2.From (17), the time derivative of δi(t) along system (15) is given as
(20)δ˙i=ki(A2-xi+B2-ui)φeqi+ki(a2-xi+b2-ui-I˙Li_ref)-k3(A3-xi+B3-ui-V˙bus_ref)+αei=0.The notationA2- refers to the second column of the A matrix.Thus, the equivalent control input is given as(21)φeqi=-RiHi,i=1,2,
where
(22)Ri=ki(a2-xi+b2-ui-I˙Li_ref)-k3(A3-xi+B3-ui-S˙ref)+αei,(23)Hi=ki(A2-xi+B2-ui).It is well known that the exact value of the external disturbance and the parameter variations of the system, such as internal resistance and magnetizing current, are difficult to measure in advance for practical applications. Therefore, I propose a nonlinear switching control-inputφsi to estimate the upper bound of uncertainties and external disturbance:
(24)φsi(t)=-γi(t)Hi(t)sign(δi(t)).In the presence of parameter uncertainties and external disturbances the dynamic equation (15) becomes
(25)X˙(t)=F(x)+G(x)φ+Δ(t),
where Δ includes the uncertainties and perturbations.In order to reduce chattering, the most common method is to replace the sign function by the saturation functionsat(δ(t),ζ), where
(26)sat(δ(t),ζ)={sign(δ),|δ|>ζ>0,δζ,|δ|≤ζ
and ζ is a small positive constant.Hence, the SMDC rule can be designed as(27)φi(t)=φeqi(t)-γi(t)Hi(t)sat(δi(t),ζ).The main theorem, stated and proved below, provides sufficient conditions to ensure the stability and robustness of the multisources system.Theorem 1.
Let the nonlinear multisources system given by (15) and sliding surface given by (17), and suppose that the Lyapunov function is defined by
(28)V(t)=12δ1(t)2+12δ2(t)2.
Then, the sufficient stability condition for the SMDC in the sense of Lyapunov should satisfy the following range of switching gains:(29)γi(t)≥|Δ(t)|,i=1,2,
where |Δ| represent the upper bound of the uncertainties and perturbations.Proof.
IfV(t) is a Lyapunov function candidate defined as in (28), the asymptotic stability will be satisfied if V˙(t)≤0, with
(30)V˙(t)=δ1(t)δ˙1(t)+δ2(t)δ˙2(t).
Using (17), (25), and (27) in (30) I obtain
(31)V˙(t)=δ1(t)[-γ1(t)sat(δ1(t),ζ)+Δ(t)]+δ2(t)[-γ2(t)sat(δ2(t),ζ)+Δ(t)]≤0.As a conclusion, to satisfy the conditionV˙(t)≤0, and to compensate the bounded parametric uncertainties and disturbances, I must restrict to sufficient condition (29).


## 3.3. Evaluation of Low Level Control Design

To illustrate the low level control design and performance evaluation of my structure, a hybrid power system coupled to two sources through a similar DC/DC converter is considered. The characteristics of the full bridge converters are(32)Rmos=0.005Ω,Rl=2Ω,n=12,Vd=0.3v,L=1e-3H,Lm=20e-6H,rp=0.05Ω,rs=0.05Ω,Ce=20e-6F,Re=56e3Ω,CDC=10e-6F.The DC bus voltage reference is set at 100 V. Simulations are obtained with sampling intervalTe=50 μs.We suppose that 40% of the load demands are supplied with the first source and 60% with the second one.Responses obtained with my algorithm are compared with ones resulting from the PI algorithm with cascade structure. The linearization techniques are applied to the converters in order to deduce linear parameters [28].The simulation conditions are identical for both controllers.To compare the global performances of my algorithm with PI one, let us consider the normalised mean square error (NMSE, Table1), that is interpreted as the overall deviations between output plant and desired values and is defined as:
(33)NMSE=∑l(Vbus-Vbus_ref)2∑l(Vbus_ref)2.Table 1
NMSE.
PI
SMDC
1.2483
0.0075In order to test the capacity of SMDC to reject disturbances, and to investigate the robustness of the control scheme, the load and sources dynamics of the hybrid system are supposed fast. The sources and load profiles are represented by Figures9 and 10.Figure 9
Sources profiles.Figure 10
Load profiles.The closed-loop responses of load current DC bus voltage, tracking voltage error, and load current obtaining with PI controller and sliding mode controller are shown in Figures11, 12, and 13, respectively.Figure 11
DC bus voltage plot obtained with PI controller and sliding mode controller.Figure 12
Tracking voltage error plot (logarithmic scale) obtained with PI controller and sliding mode controller.Figure 13
Load current plot obtained with PI controller and sliding mode controller.Given the fast dynamics imposed on the sources and the load as evident from Figures9 and 10, the simulation results reveal that the SMDC is advantageous in providing negligible steady-state errors (NMSE = 0.0075), to adjust the flow source that meets the load demand and guaranteed stable system (Figure 12). It can also be seen that the responses are satisfactory in terms of overshoot, settling time, and fall time. This shows an excellent behavior of the sliding mode controller in comparison with PI controllers which does not run correctly for wide variations of operating point.


## 4. High Level Controller: Intelligent Power and Energy Management Strategy

Generally, an HPS is a dynamic system of electric power production. While during their operation, fluctuation in wind speed, solar radiation, and load demand will occur, to ensure a reliable power supply under all operating conditions, the global control strategy is critical to control efficiently the energy flow in the multisource system.Indeed, under the operation constraints related to each type of sources, the energy flow must be managed such that the load demand is satisfied for all time (preserve the reliability of the system), the ESS is managed in an optimal way (the battery and ultracapacitor SOC minimum and maximum threshold values which must be respected, and the current and voltage should be limited to maximum absolute value in order to increase the ESS durability), and the fuel use is optimal.These requirements are achieved by improving the rate of renewable energy penetration, by an appropriate controller of CE and finally, and by a specific management of the charge-discharge ESS operation.Since the desired behavior of the multisource system is well known and can be described using linguistic variables, the use of an online Takagi-Sugeno fuzzy dynamic supervisor seems appropriate. Indeed, the fuzzy approach has many advantages, such as the possibility to use multiple input variables without increasing the complexity of the controller, and the desired behavior can be described by simplifies linguistic variables but it would be difficult to express mathematically and finally, historical data is stored only for a short time window, which perfectly meets the real-time application requirements.


### 4.1. FLC’s Structure

Based on reviews of multisource power, the block diagram of the global control strategy is illustrated in Figure14.Figure 14
Block diagram of dynamic fuzzy logic controller for power management strategy.The required behavior was implemented with the help of six input variables (the energy sources availability and load power demand). These information are useful to determine whether there is surplus or deficient energy to supply the load demand.The load power estimation (P_load) and there variation (ΔPload) are used to evaluate instantaneously the load requirements and to decide which ESS operating mode must be used (battery with high-specific energy or ultracapacitor with high-specific power).The SOCs of ESS are used to prevent an excessive discharge or an overcharge of the ESS and to know how much supply energy is available. It is estimated from the SOC model calculated depending on ESS power (batteries and ultracapacitors) that is also a fuzzy output.The online dynamic fuzzy logic controller output variables are a reference power vector (power set points: battery power (Pbat), ultracapicitor power (Puc), combustion engine produced power and Pren) sent to the low level control, and CE required power sent to the CE controller.


### 4.2. Membership Functions

Trapezoidal type membership functions were used for the linguistic input variables, as shown in Figures15(a) to 15(e), because they produce smoother control action due to the flatness at the top of the trapezoidal shape [5].Input membership functions for: (a) SOCbat, (b) SOCuc, (c) renewable energy power, (d) load demand power variation, and (e) CE power (Pce-slo: combustion engine slow down power).
(a)
(b)
(c)
(d)
(e)The battery state of charge (SOCbat) input was divided using the linguistic variables: too low (TL), low (L), normal (N) and high (H), as shown on Figure 15(a). The SOC window defines the relationship between the battery’s total energy and its usable energy. Fuzzy sets (TL) and (H) represent the ranges where the SOC should not be. Fuzzy set (N) represents the range where the SOC should be, and (L) acts as an energy reserves level necessary to satisfy the load power requirement, when the CE is on start-up phase.The membership function for ultracapacitor state of charge (SOCuc) is presented at Figure 15(b). Fuzzy sets (L) and (H) represent the ranges of ultracapacitor operating limit value. Fuzzy set (N) represents the optimal SOC range.Depending on the instantaneous load power demand, the membership function for renewable power (Pren) is divided of two ranges: positive (+) and negative (−) (Figure 15(c)).The fuzzy set (−) represents the range where the renewable sources must be assisted by ESS or CE to satisfy the power demand. Fuzzy set (+) represents the autonomous supply mode ensured by renewable sources (ESS is on recharge mode and CE is on Off mode).The demand power variation was divided using the variables positive (+) and negative (−), as shown in Figure15(d). This information is useful to prevent abusive use of the batteries. In fact, when the demand power variation is below a predetermined level, the ESS discharge-mode must be ensured by batteries; otherwise, the ESS discharge mode must be ensured by ultracapacitor (absorption the peaking power demand).The membership function for CE power is presented at the Figure15(e). The fuzzy sets (+) represent the CE Production mode (the power from the generator is effective to supply the load and, simultaneously, recharge the ESS if necessary). The fuzzy sets (−) represent the CE slow-down mode or transient state (the generator do not contribute to provide load power).


### 4.3. Rules

The rules were chosen very intuitively according to the main following points.(1)
A RES takes over as main energy source.(2)
In the optimal conditions, the battery bank and the ultracapacitor have been dimensioned in order to satisfy the maximum load power demand; that is, 100% of the energy demand can be supplied by the ESS.(3)
An ESS comes into play in order to balance the renewable power fluctuations and to allow for a CE Off mode (zero emission).(4)
The CE is turned on when the ESS state of charge drops below a minimum value and renewable power is insufficient to meet the load.(5)
The minimum period of CE operation must cover the ESS time of recharge.The latter (4 and 5) are important to avoid running the CE for very short period of time.Using these linguistic directives, it makes easy to obtain a set of fuzzy rules. Thus, the fuzzy rule basis consists of a collection of fuzzy if-then rules type. For example, one of the energy management rules (rule1) can be stated as (IF the renewable power exceeds load demand power AND battery state of charge is high AND ultracapacitor state of charge is high THEN the renewable produced power is equal to load demand power). Using these linguistic variables, a set of fuzzy rules was developed. For the sake of simplicity, the rules number was minimized to 16 rules. These rules are shown in (Table2). Where Pbat_max is a maximum battery power, Puc_max is a maximum ultracapacitor power, Pbat_rech is a recharge battery power, and Puc_rech is a recharge ultracapicitor power.Table 2
Rules.
Rule
IfPren is
And SOCuc is
And SOCbat is
And ΔPhvd is
AndPce_pro is
ThenPbat
AndPuc
AndPren
AndPce_req
AndPce_pro
1
+
H
H
/
/
0
0
P
hvd
0
0
2
+
H
N
/
−
P
bat_rech
0
Min(Pren, (Phvd + Pbat_max))
0
0
3
+
N OR H
L
/
−
P
bat_rech
0
Min(Pren, (Phvd + Pbat_max))
0
0
4
+
/
TL
/
−
P
bat_rech
0
Min(Pren, (Phvd + Pbat_max))
0
0
5
+
N
N OR H
/
−
0
P
uc_rech
Min(Pren, (Phvd + Puc_max))
0
0
6
+
L
L, N OR H
/
−
0
P
uc_rech
Min(Pren, (Phvd + Puc_max))
0
0
7
−
/
N OR H
−
−
P
hvd − Pren_max
0
P
ren_
max
0
0
8
−
L
N OR H
+
−
P
hvd − Pren_max
0
P
ren_
max
0
0
9
−
N OR H
L OR TL
/
−
0
P
hvd − Pren_max
P
ren_
max
0
0
10
−
N OR H
N OR H
+
−
P
bat_
max
P
hvd− Pren_max − Pbat_max
P
ren_
max
0
0
11
−
L
L
/
−
P
hvd − Pren_max
0
P
ren_
max
0
P
hvd (startup)
12
−
TL
L
/
−
0
0
P
ren_
max
0
P
hvd (startup)
13
/
L, N
TL, L, N
/
+
P
bat_rech
P
uc_rech
Min(Pren, (Phvd + Pbat_max + Puc_max))
Max(0, min(Pce_pro, Phvd − Pren + Pbat_rech + Puc_rech))
Max(Pce_slo, Phvd − Pren + Pbat_max + Puc_max)
14
/
L, N
H
/
+
0
P
uc_rech
Min(Pren, (Phvd + Puc_max))
Max(0, min(Pce_pro, Phvd − Pren + Puc_rech))
Max(Pce_slo, Phvd − Pren + Puc_max)
15
/
H
TL, L, N
/
+
P
bat_rech
0
Min(Pren, (Phvd + Pbat_max))
Max(0, min(Pce_pro, Phvd − Pren + Pbat_rech))
Max(Pce_slo, Phvd − Pren + Pbat_max)
16
−
H
H
/
+
0
0
Min(Pren, Phvd)
Max(0, min(Pce_pro, Phvd − Pren))
Max(Pce_slo, Phvd − Pren)The logical AND (resp.) has been implemented with the minimum (resp. maximum) operator.While the fuzzy input variables are obtained either from a prediction method based on the model of each source or else using local measurements, to obtain the output of the controller, the degrees of membership functions of the if-parts of all rules are evaluated, and the THEN-parts of all rules are averaged, weighted by these degrees of membership functions.


## 4.1. FLC’s Structure

Based on reviews of multisource power, the block diagram of the global control strategy is illustrated in Figure14.Figure 14
Block diagram of dynamic fuzzy logic controller for power management strategy.The required behavior was implemented with the help of six input variables (the energy sources availability and load power demand). These information are useful to determine whether there is surplus or deficient energy to supply the load demand.The load power estimation (P_load) and there variation (ΔPload) are used to evaluate instantaneously the load requirements and to decide which ESS operating mode must be used (battery with high-specific energy or ultracapacitor with high-specific power).The SOCs of ESS are used to prevent an excessive discharge or an overcharge of the ESS and to know how much supply energy is available. It is estimated from the SOC model calculated depending on ESS power (batteries and ultracapacitors) that is also a fuzzy output.The online dynamic fuzzy logic controller output variables are a reference power vector (power set points: battery power (Pbat), ultracapicitor power (Puc), combustion engine produced power and Pren) sent to the low level control, and CE required power sent to the CE controller.


## 4.2. Membership Functions

Trapezoidal type membership functions were used for the linguistic input variables, as shown in Figures15(a) to 15(e), because they produce smoother control action due to the flatness at the top of the trapezoidal shape [5].Input membership functions for: (a) SOCbat, (b) SOCuc, (c) renewable energy power, (d) load demand power variation, and (e) CE power (Pce-slo: combustion engine slow down power).
(a)
(b)
(c)
(d)
(e)The battery state of charge (SOCbat) input was divided using the linguistic variables: too low (TL), low (L), normal (N) and high (H), as shown on Figure 15(a). The SOC window defines the relationship between the battery’s total energy and its usable energy. Fuzzy sets (TL) and (H) represent the ranges where the SOC should not be. Fuzzy set (N) represents the range where the SOC should be, and (L) acts as an energy reserves level necessary to satisfy the load power requirement, when the CE is on start-up phase.The membership function for ultracapacitor state of charge (SOCuc) is presented at Figure 15(b). Fuzzy sets (L) and (H) represent the ranges of ultracapacitor operating limit value. Fuzzy set (N) represents the optimal SOC range.Depending on the instantaneous load power demand, the membership function for renewable power (Pren) is divided of two ranges: positive (+) and negative (−) (Figure 15(c)).The fuzzy set (−) represents the range where the renewable sources must be assisted by ESS or CE to satisfy the power demand. Fuzzy set (+) represents the autonomous supply mode ensured by renewable sources (ESS is on recharge mode and CE is on Off mode).The demand power variation was divided using the variables positive (+) and negative (−), as shown in Figure15(d). This information is useful to prevent abusive use of the batteries. In fact, when the demand power variation is below a predetermined level, the ESS discharge-mode must be ensured by batteries; otherwise, the ESS discharge mode must be ensured by ultracapacitor (absorption the peaking power demand).The membership function for CE power is presented at the Figure15(e). The fuzzy sets (+) represent the CE Production mode (the power from the generator is effective to supply the load and, simultaneously, recharge the ESS if necessary). The fuzzy sets (−) represent the CE slow-down mode or transient state (the generator do not contribute to provide load power).


## 4.3. Rules

The rules were chosen very intuitively according to the main following points.(1)
A RES takes over as main energy source.(2)
In the optimal conditions, the battery bank and the ultracapacitor have been dimensioned in order to satisfy the maximum load power demand; that is, 100% of the energy demand can be supplied by the ESS.(3)
An ESS comes into play in order to balance the renewable power fluctuations and to allow for a CE Off mode (zero emission).(4)
The CE is turned on when the ESS state of charge drops below a minimum value and renewable power is insufficient to meet the load.(5)
The minimum period of CE operation must cover the ESS time of recharge.The latter (4 and 5) are important to avoid running the CE for very short period of time.Using these linguistic directives, it makes easy to obtain a set of fuzzy rules. Thus, the fuzzy rule basis consists of a collection of fuzzy if-then rules type. For example, one of the energy management rules (rule1) can be stated as (IF the renewable power exceeds load demand power AND battery state of charge is high AND ultracapacitor state of charge is high THEN the renewable produced power is equal to load demand power). Using these linguistic variables, a set of fuzzy rules was developed. For the sake of simplicity, the rules number was minimized to 16 rules. These rules are shown in (Table2). Where Pbat_max is a maximum battery power, Puc_max is a maximum ultracapacitor power, Pbat_rech is a recharge battery power, and Puc_rech is a recharge ultracapicitor power.Table 2
Rules.
Rule
IfPren is
And SOCuc is
And SOCbat is
And ΔPhvd is
AndPce_pro is
ThenPbat
AndPuc
AndPren
AndPce_req
AndPce_pro
1
+
H
H
/
/
0
0
P
hvd
0
0
2
+
H
N
/
−
P
bat_rech
0
Min(Pren, (Phvd + Pbat_max))
0
0
3
+
N OR H
L
/
−
P
bat_rech
0
Min(Pren, (Phvd + Pbat_max))
0
0
4
+
/
TL
/
−
P
bat_rech
0
Min(Pren, (Phvd + Pbat_max))
0
0
5
+
N
N OR H
/
−
0
P
uc_rech
Min(Pren, (Phvd + Puc_max))
0
0
6
+
L
L, N OR H
/
−
0
P
uc_rech
Min(Pren, (Phvd + Puc_max))
0
0
7
−
/
N OR H
−
−
P
hvd − Pren_max
0
P
ren_
max
0
0
8
−
L
N OR H
+
−
P
hvd − Pren_max
0
P
ren_
max
0
0
9
−
N OR H
L OR TL
/
−
0
P
hvd − Pren_max
P
ren_
max
0
0
10
−
N OR H
N OR H
+
−
P
bat_
max
P
hvd− Pren_max − Pbat_max
P
ren_
max
0
0
11
−
L
L
/
−
P
hvd − Pren_max
0
P
ren_
max
0
P
hvd (startup)
12
−
TL
L
/
−
0
0
P
ren_
max
0
P
hvd (startup)
13
/
L, N
TL, L, N
/
+
P
bat_rech
P
uc_rech
Min(Pren, (Phvd + Pbat_max + Puc_max))
Max(0, min(Pce_pro, Phvd − Pren + Pbat_rech + Puc_rech))
Max(Pce_slo, Phvd − Pren + Pbat_max + Puc_max)
14
/
L, N
H
/
+
0
P
uc_rech
Min(Pren, (Phvd + Puc_max))
Max(0, min(Pce_pro, Phvd − Pren + Puc_rech))
Max(Pce_slo, Phvd − Pren + Puc_max)
15
/
H
TL, L, N
/
+
P
bat_rech
0
Min(Pren, (Phvd + Pbat_max))
Max(0, min(Pce_pro, Phvd − Pren + Pbat_rech))
Max(Pce_slo, Phvd − Pren + Pbat_max)
16
−
H
H
/
+
0
0
Min(Pren, Phvd)
Max(0, min(Pce_pro, Phvd − Pren))
Max(Pce_slo, Phvd − Pren)The logical AND (resp.) has been implemented with the minimum (resp. maximum) operator.While the fuzzy input variables are obtained either from a prediction method based on the model of each source or else using local measurements, to obtain the output of the controller, the degrees of membership functions of the if-parts of all rules are evaluated, and the THEN-parts of all rules are averaged, weighted by these degrees of membership functions.


## 5. Application to the Multisource System

The objective of this simulation is not to produce a realistic practical scenario of HES, but to validate the developed energy management strategy. For this reason, the simulations are made with simplified models of renewable sources, CE, battery and ultracapacitor. Furthermore, the simulations aim to show that the developed energy management strategy can operate in different system configurations with different characteristics. Therefore, any HES state can be controlled without redesign the energy management strategy.In the presented simulations, to illustrate all possible operating modes of multisource system, I suppose that the battery time constant is reduced to few seconds (5 s). The simulation results for the power balance of the HES are presented in Figure16.Figure 16
Simulation results for the power balance of the HES: (a) and (b) present  the SOC and currant (A), respectively, of battery and ultracapacitor, (c) shows the generated and consumed power (W) output from the combustion engine, and (d) shows the generated power (W) output from the renewable sources, the demanded and consumed power by the load.It can be observed from these results that the optimal operation objective of HES as mentioned below has been achieved.(1)
The renewable sources take over as main energy sources.(2)
An energy storage system comes into play in order to balance the renewable power fluctuations and to allow for a CE Off mode (zero emission).(3)
The battery and ultracapacitor SOCs were stable basically and decreased (current is negative)/increased (current is positive) slightly from 0.3 to 0.9 (the operation limits).(4)
The CE turns on when the ESS state of charge drops below a minimum value, and renewable power is insufficient to meet the load (as shown in Figure16(c) at 2 s and 17 s).(5)
The minimum period of CE operation covers the ESS time of recharge (as shown in Figure16(c) between 2 s–7 s and 17 s–20 s).(6)
Load demand is satisfied for all time which preserve the reliability of the system (Figure16(d)).


## 6. Conclusion and Prospects

This paper presented the basis for the development of an intelligent and autonomous energy management strategy for HPS. Two hierarchical levels are presented to control and manage the HPS. The low level is performed by a local sliding mode control unit (DC-DC converters controller) of the different power sources. This control technique provides good overall performances compared to standard current control and good robustness against load and input voltage variations. The high level is performed by the online supervisor unit. This unit is designed by applying online dynamic Takagi-Sugeno fuzzy logic principles. As a result the robust control system gets rid of the limits of the HPS, which has the imprecision, uncertainty, strong coupling, and nonlinearity.Under the operation constraints related to each type of sources, it can be observed from the simulation results that the optimal operation objective of HES has been achieved (the load demands are satisfied for all time, the energy storage system is managed in an optimal way, and the fuel use is optimal).The application of the proposed control scheme to large-scale renewable multi-source power systems will be considered in my future work.The neural-based adaptive observer, in order to identify unknown functions in the multi-source system and to estimate the unmeasured states, will be also studied in my further works.



---

*Source: 148563-2012-12-03.xml*