Battery Characterization and Dimensioning Approaches for Micro-Grid Systems †

Micro-Grid (MG) systems have been extensively studied and deployed to lower the power consumption while reducing the greenhouse gas emissions. Although, the challenge with Renewable Energy Sources (RES) is their uncertain and intermittent nature, things that led the researchers to think about integrating storage devices, especially batteries, into MG systems. The main aim is to store the excess of produced energy for further usage when not enough energy is available. Nevertheless, batteries modeling and characterization is mandatory to identify their parameters and study their performance within MG systems. Moreover, in order to continuously supply electricity to the building, it is required to figure out the optimum size of energy production systems and storage devices. This paper introduces a methodology for MG modeling and performance evaluation. Its main contribution is twofold, (i) battery’s parameters identification, and (ii) modeling and dimensioning method for both standalone and MG systems. An instrumentation platform, composed of recent sensing and actuating equipment, for MG energy management and battery characterization is developed. Simulation and experimental results show the effectiveness of the proposed methodology.


Introduction
Due to their advantages and benefits (e.g., clean, abundant), Morocco focuses its interest on RES (especially solar energy with an average daily global radiation of around 5 kWh/m 2 [1]) in order to decrease its external dependence on fossil fuel products (e.g., coal, fuel, oil, natural gas), which is around 95%. For these advantages, MG systems have been broadly studied and deployed in order to decrease the power consumption from the electric grid while minimizing the greenhouse emissions. Though, because of the intermittent nature of RES, it is mandatory to store the excess of produced energy from RES (e.g., solar energy) in order to supply electricity to the building's appliances when there is no production. On the other hand, several energy storage systems exist, filters (e.g., passive filters, active filters and hybrid filters) is necessary in order to improve the quality of the power. Nevertheless, these filters are costly and have reasonable losses, leading to decrease in system's efficiency [23]. For all these problems, researchers have developed a new variety of rectifiers using new solid state self-commutating devices such as MOSFETs, IGBTs in order to increase their efficiency and the quality of the power while minimizing their cost [23].
Though, authors in [24,25] have reviewed some deployed MG test-beds all over the world, and have shown the importance of the deployment of these systems in real-siting in order to validate the MG components' models used in simulations. For that, we have conducted the work presented in this article, which is a part of the work we are undertaking under an ongoing project, named MIGRID, aiming to model, simulate and experiment an MG system for sizing and control approaches [26]. The system is composed of five components, which are PV panels for producing electricity, batteries in order to supply electricity to the building when there is no production, a regulator to control the batteries' charge and discharge, a rectifier that converts the alternating current (AC) received from the electric grid into a direct current (DC) for operating the DC loads, and the electric grid for supplying the building when there is no production and the batteries are empty. The contribution in this article is as follows: the system has been first modeled and evaluated using both simulations and experiments under the same conditions, with the aim to study its behavior and validate the developed components' models. Nevertheless, batteries modeling and characterization is mandatory to identify their parameters and study their performance within MG systems. Therefore, an instrumentation platform, composed of recent sensing and actuating equipment, is developed for evaluating the efficiency of MG components' models.
The remainder of this article is structured as follows. Section 2 describes the proposed methodology followed throughout the article. The description, the modeling and the simulation and experimental results of the MG system are introduced in Section 3. In Section 4, a sizing method is introduced in order to make the MG system a stand-alone system, and then simulation and experimental results are presented. Finally, the conclusions and perspectives are given in Section 5. Figure 1 describes the different steps of the proposed methodology, (i) battery modeling and characterization, (ii) MG system modeling, simulations and experiments, and (iii) Stand-alone system sizing and performance evaluation. Modeling and characterization is a mandatory step to identify battery's parameters and study their behaviour and performance within MG systems. Therefore, an instrumentation platform, composed of recent sensing and actuating equipment, has been developed in order to collect the battery's characteristics (e.g., current and voltage). The data collected have been fed to the Recursive Least Squares (RLS) algorithm. It is an algorithm able to identify the battery's parameters that minimize the difference between the simulation and experimental results. These parameters have then been integrated into the used battery model presented in [15,16]. Simulations and experiments have been conducted under the same conditions in order to validate the proposed battery characterization methodology. Concerning the MG system, each component of the system has been modeled and connected to the battery's model. The aim is to simulate the whole system, using the same conditions as in the experimentation in order to figure out first the system's behaviour and validate the components' models. These validated models could be used for conducting further simulations accordingly. Regarding the third step, a sizing method has been introduced for stand-alone systems. Then, simulations have been carried out using the components' models validated in the second step. Similarly, experimentations have been conducted and obtained results are compared to those from simulations. The results validate the introduced sizing method and show full independence from the electric grid. This methodology can be applied to other stand-alone systems.  Figure 2 describes the architecture of the MG system that we have deployed in our EEBLab testbed. The system presented above is composed of a PV module with 265 W as a peak power, a 24 Ah and 12 V Lead-acid battery fabricated by the Genesis NP24-12R, a regulator "BlueSolar MPPT 150/35", and a 24 W and 12 V ventilation system composed of two fans operating with a PWM control. The grid is connected to the load through a rectifier with an output voltage of 12 V. Moreover, a platform [27], constituted by a set of voltage and current sensors and an actuator connected to an Arduino board, has been deployed in order to measure the voltage and the current in all branches of the MG system (see Figure 3).    The system presented above is composed of a PV module with 265 W as a peak power, a 24 Ah and 12 V Lead-acid battery fabricated by the Genesis NP24-12R, a regulator "BlueSolar MPPT 150/35", and a 24 W and 12 V ventilation system composed of two fans operating with a PWM control. The grid is connected to the load through a rectifier with an output voltage of 12 V. Moreover, a platform [27], constituted by a set of voltage and current sensors and an actuator connected to an Arduino board, has been deployed in order to measure the voltage and the current in all branches of the MG system (see Figure 3). The system presented above is composed of a PV module with 265 W as a peak power, a 24 Ah and 12 V Lead-acid battery fabricated by the Genesis NP24-12R, a regulator "BlueSolar MPPT 150/35", and a 24 W and 12 V ventilation system composed of two fans operating with a PWM control. The grid is connected to the load through a rectifier with an output voltage of 12 V. Moreover, a platform [27], constituted by a set of voltage and current sensors and an actuator connected to an Arduino board, has been deployed in order to measure the voltage and the current in all branches of the MG system (see Figure 3).  The data collected from sensors are sent via a Raspberry Pi to a cluster for processing and storage. Furthermore, using this platform, the battery's SoC, which is necessary for controlling the MG system, is also estimated using the coulomb counting method (Ah method). Besides, the actuator, which is controlled by the Arduino board, has been used in order to alter between the PV-Battery system and the electric grid when the SoC of the battery reaches 40%. Moreover, the state and the powers in each branch of the system are visualized and extracted using a web application that has been already developed. The aim of the deployment of this system is to compare the simulation and experiment results in order to validate the components' models and behavior, especially at night when there is no PV production.

Simulation and Experimental Results of the MG System
The first goal of this work is to model the MG system and conduct simulations under MATLAB/Simulink using similar experimental conditions. For this purpose, the model of each component (e.g., PV module, regulator, battery, and rectifier) of the system have been developed. In a previous work [6], we have described with more details the models of the PV module and the regulator. For the model of the AC/DC converter, four IGBTs (Insulated Gate Bipolar Transistor) have been used and controlled by a PWM (Pulse-width modulation) signal and a filter have been also used in order to improve the quality of the power.
Regarding the battery, the model presented in [15,16] was used for the purpose of estimating the battery's behavior. This model is a variant of the model presented in [14]. It takes into consideration, in addition to the OCV (Open Circuit Voltage), which is a function of the battery's SoC and the polarization resistance, the polarization voltage for the purpose of representing the non-linearity relationship between the OCV and the SoC. The equations of the Lead acid battery voltage (in Volt) in both charge (i < 0) (1) and discharge (i > 0) (2) modes are expressed respectively as follows [15,16]: where, Vbatt is the battery voltage (V), E0 is the battery constant voltage (V), K is the polarization constant (V/(Ah)) or polarization resistance (Ω), Q is the battery capacity (Ah), it is the actual battery capacity (Ah), R is the internal resistance (Ω), i is the battery current (A), i* is the filtered current (A), The data collected from sensors are sent via a Raspberry Pi to a cluster for processing and storage. Furthermore, using this platform, the battery's SoC, which is necessary for controlling the MG system, is also estimated using the coulomb counting method (Ah method). Besides, the actuator, which is controlled by the Arduino board, has been used in order to alter between the PV-Battery system and the electric grid when the SoC of the battery reaches 40%. Moreover, the state and the powers in each branch of the system are visualized and extracted using a web application that has been already developed. The aim of the deployment of this system is to compare the simulation and experiment results in order to validate the components' models and behavior, especially at night when there is no PV production.
The first goal of this work is to model the MG system and conduct simulations under MATLAB/Simulink using similar experimental conditions. For this purpose, the model of each component (e.g., PV module, regulator, battery, and rectifier) of the system have been developed. In a previous work [6], we have described with more details the models of the PV module and the regulator. For the model of the AC/DC converter, four IGBTs (Insulated Gate Bipolar Transistor) have been used and controlled by a PWM (Pulse-width modulation) signal and a filter have been also used in order to improve the quality of the power.
Regarding the battery, the model presented in [15,16] was used for the purpose of estimating the battery's behavior. This model is a variant of the model presented in [14]. It takes into consideration, in addition to the OCV (Open Circuit Voltage), which is a function of the battery's SoC and the polarization resistance, the polarization voltage for the purpose of representing the non-linearity relationship between the OCV and the SoC. The equations of the Lead acid battery voltage (in Volt) in both charge (i < 0) (1) and discharge (i > 0) (2) modes are expressed respectively as follows [15,16]: where, V batt is the battery voltage (V), E 0 is the battery constant voltage (V), K is the polarization constant (V/(Ah)) or polarization resistance (Ω), Q is the battery capacity (Ah), i t is the actual battery capacity (Ah), R is the internal resistance (Ω), i is the battery current (A), i* is the filtered current (A), and exp(t) is an expression that represents the exponential zone voltage (V). The latter term represents the hysteresis phenomenon between the charge and the discharge that takes place in the exponential zone as follows [15]: .
where i(t) is the battery current (A), u(t) represents the charge (u(t) = 1) or discharge (u(t) = 0) mode, A is the exponential zone amplitude (V), and B is the exponential zone time constant inverse (Ah) (−1) . The battery's model has been built using these equations in Matlab/Simulink (R2017a). However, in order to use this model, the identification of the battery's parameters (i.e., estimate the battery's parameters of our deployed battery using the experimental data) is required. For that, as illustrated in Figure 4, an instrumentation platform was developed in order to extract the battery's characteristics, which are the battery's voltage and current, during the charge and the discharge modes using voltage and current sensors, which are connected to an Arduino board in order to gather the data. The gathered data are then sent to a cluster via a Raspberry Pi for processing, filtering, storage and real time visualization [27]. and exp(t) is an expression that represents the exponential zone voltage (V). The latter term represents the hysteresis phenomenon between the charge and the discharge that takes place in the exponential zone as follows [15]: where i(t) is the battery current (A), u(t) represents the charge (u(t) = 1) or discharge (u(t) = 0) mode, A is the exponential zone amplitude (V), and B is the exponential zone time constant inverse (Ah) (−1) . The battery's model has been built using these equations in Matlab/Simulink (R2017a). However, in order to use this model, the identification of the battery's parameters (i.e., estimate the battery's parameters of our deployed battery using the experimental data) is required. For that, as illustrated in Figure 4, an instrumentation platform was developed in order to extract the battery's characteristics, which are the battery's voltage and current, during the charge and the discharge modes using voltage and current sensors, which are connected to an Arduino board in order to gather the data. The gathered data are then sent to a cluster via a Raspberry Pi for processing, filtering, storage and real time visualization [27]. Therefore, several charge/discharge experiments have been carried out and the gathered data have been processed and filtered. Then, the Recursive Least Squares method (RLS) [28], which is a method that recursively comes up with the parameters that minimize the difference between the simulation and experimental results, has been applied to these data in order to identify the battery's parameters that are presented in Table 1.   [29]. For instance, Figure 5 (Resp. Figure 6) depicts the battery's current and voltage in the charging mode (Resp. the discharging mode) from both simulations and experiments. As observed in these figures, the curves from both experiments and simulations show the same behavior, and the errors between both of them are in the vicinity of the threshold (i.e., the maximum error must be fewer than 5% [30]), which validates the accuracy of the proposed battery characterization methodology. The slight errors between the simulations and experiments may be caused by the measurement noises and the instability of the sensors, or especially by not taking into consideration in the battery's model the temperature's effect on its properties. Therefore, several charge/discharge experiments have been carried out and the gathered data have been processed and filtered. Then, the Recursive Least Squares method (RLS) [28], which is a method that recursively comes up with the parameters that minimize the difference between the simulation and experimental results, has been applied to these data in order to identify the battery's parameters that are presented in Table 1. Finally, the identified parameters have been integrated into the battery's model for validating the proposed battery characterization methodology. Simulations and experiments have been conducted under the same conditions (i.e., the same current have been applied and the computed battery's voltage in both experimentation and simulation has been compared), and results show the accuracy of the characterization methodology [29]. For instance, Figure 5 (Resp. Figure 6) depicts the battery's current and voltage in the charging mode (Resp. the discharging mode) from both simulations and experiments. As observed in these figures, the curves from both experiments and simulations show the same behavior, and the errors between both of them are in the vicinity of the threshold (i.e., the maximum error must be fewer than 5% [30]), which validates the accuracy of the proposed battery The battery's model with the identified parameters was then integrated into our MG system (Figure 7) in order to simulate the system's behavior and compare the simulations to the experiments.   The battery's model with the identified parameters was then integrated into our MG system (Figure 7) in order to simulate the system's behavior and compare the simulations to the experiments.  The battery's model with the identified parameters was then integrated into our MG system (Figure 7) in order to simulate the system's behavior and compare the simulations to the experiments.
In order to show the system's behavior and validate the components' models, several simulations and experiments have been carried out during a sunny day, and the results are depicted in Figure 8. For instance, Figure 8b,c show that during the daytime, the electricity needed by the ventilation system is totally supplied by the PV panel (the actuator is switched to the PV-Battery system) whereas the battery is charged by the excess energy production. Nevertheless, during the night, the electric grid operates the load, especially after 5AM when the battery's SoC (Figure 8a), which is estimated using the Ah method, is fewer than 40%. The battery's model with the identified parameters was then integrated into our MG system (Figure 7) in order to simulate the system's behavior and compare the simulations to the experiments.  In order to show the system's behavior and validate the components' models, several simulations and experiments have been carried out during a sunny day, and the results are depicted in Figure 8. For instance, Figure 8b,c show that during the daytime, the electricity needed by the ventilation system is totally supplied by the PV panel (the actuator is switched to the PV-Battery system) whereas the battery is charged by the excess energy production. Nevertheless, during the night, the electric grid operates the load, especially after 5AM when the battery's SoC (Figure 8a), which is estimated using the Ah method, is fewer than 40%. According to these results, the electricity produced by the PV is not sufficient to totally supply the ventilation system. Thus, a sizing method was established with the purpose of making the system fully autonomous. Then, simulations and experiments have been carried out for validating the sizing method towards a total independence from the electric grid.

Sizing of the Stand-Alone system
The sizing of the MG system, which is introduced in the previous section, depends on many important parameters (e.g., consumption profile, solar irradiation) [4][5][6]. Regarding the solar irradiation, the PVGIS application that provides the solar irradiation using the geographical coordinates of any site in Africa or Europe, have been used. In order to provide the average monthly solar radiation all over the year, geographical coordinates of our site (Latitude 34°00′47″ N, Longitude 6°49′57″ W, Altitude 46 m) were inputted. Concerning the consumption, a set of voltage and current sensors have been deployed in order to monitor the ventilation system's daily electricity consumption (see Figure 9), which is around 430 Wh. According to these results, the electricity produced by the PV is not sufficient to totally supply the ventilation system. Thus, a sizing method was established with the purpose of making the system fully autonomous. Then, simulations and experiments have been carried out for validating the sizing method towards a total independence from the electric grid.

Sizing of the Stand-Alone System
The sizing of the MG system, which is introduced in the previous section, depends on many important parameters (e.g., consumption profile, solar irradiation) [4][5][6]. Regarding the solar irradiation, the PVGIS application that provides the solar irradiation using the geographical coordinates of any site in Africa or Europe, have been used. In order to provide the average monthly solar radiation all over the year, geographical coordinates of our site (Latitude 34 • 00 47" N, Longitude 6 • 49 57" W, Altitude 46 m) were inputted. Concerning the consumption, a set of voltage and current sensors have been deployed in order to monitor the ventilation system's daily electricity consumption (see Figure 9), which is around 430 Wh. Then, the peak power of the PV panels (in W) is computed using the following expression [5,6]: where Eelec is the maximal daily electricity consumption (kWh/day), is the irradiance at standard conditions STC (equal to 1000 W/m 2 ), Ei is the lowermost average monthly solar irradiation all over the year (kWh/m 2 /day); for our site, it is equal to 4790 kWh/m 2 /day that corresponds to the December's irradiation, PR (Performance Ratio) that depends on many parameters, such as battery's efficiency, converter's efficiency, voltage drop in the cables, other losses in the installation, the integration mode, is a correction coefficient for stand-alone systems, which is between 0.55 and 0.75 [31]. The value used in our case is equal 0.70.
In order to meet the needed electricity demand, the required PV panels' peak power is 128.243 Wc. As mentioned above, the peak power of the PV module used in this experiment is 265 W, its open circuit voltage is 38.3 V and 8.98 A is the short-circuit current. As a result, one module is required in order to operate the system all over the year, which is computed using the following equation: Nt = , with a total peak power of 265 W and a total area of 1.63 m². Concerning the regulator, the same one as in the previous experiment have been used because the PV module has not changed. Furthermore, the battery's capacity (CAh), which is in Ah, is expressed by the following equation [5]: where Eelec is the maximal daily electricity consumption (Wh/day), NAut is the number of days of autonomy (day), DoD is the Depth of Discharge (%), is the battery's efficiency (%), and Vbat is the battery's voltage (V). The estimation of the battery's capacity requires the knowledge of the number of days of autonomy, which differs from site to another according to weather conditions. The objective is to supply the electricity to the building from batteries when there is no production. In our case, the number of days of autonomy fluctuates between one and three days. In our computations, the value used is one day. A battery of 12 V with a DoD equal to 60% and a battery's efficiency of 85% has been chosen. 70.26 Ah was found to be the required capacity to install. As a result, we have chosen to work with one battery of 12 V and 75 Ah manufactured by "Electra". Besides, the electric grid is connected to the system through a rectifier that provides 12 V as output voltage with the current needed by the ventilation system.
In order to prove the accuracy of the sizing method, several experiments and simulations have been carried out. However, for simulating the whole system, the characterization of the battery is required. After identifying the battery's parameters that are shown in Table 2, the battery's model has been connected to the other components' models and several simulations and experiments have been carried out. Then, the peak power of the PV panels (in W) is computed using the following expression [5,6]: where E elec is the maximal daily electricity consumption (kWh/day), I r 0 is the irradiance at standard conditions STC (equal to 1000 W/m 2 ), E i is the lowermost average monthly solar irradiation all over the year (kWh/m 2 /day); for our site, it is equal to 4790 kWh/m 2 /day that corresponds to the December's irradiation, PR (Performance Ratio) that depends on many parameters, such as battery's efficiency, converter's efficiency, voltage drop in the cables, other losses in the installation, the integration mode, is a correction coefficient for stand-alone systems, which is between 0.55 and 0.75 [31]. The value used in our case is equal 0.70. In order to meet the needed electricity demand, the required PV panels' peak power is 128.243 Wc. As mentioned above, the peak power of the PV module used in this experiment is 265 W, its open circuit voltage is 38.3 V and 8.98 A is the short-circuit current. As a result, one module is required in order to operate the system all over the year, which is computed using the following equation: N t = P c P module , with a total peak power of 265 W and a total area of 1.63 m 2 . Concerning the regulator, the same one as in the previous experiment have been used because the PV module has not changed. Furthermore, the battery's capacity (C Ah ), which is in Ah, is expressed by the following equation [5]: where E elec is the maximal daily electricity consumption (Wh/day), N Aut is the number of days of autonomy (day), DoD is the Depth of Discharge (%), η B is the battery's efficiency (%), and V bat is the battery's voltage (V). The estimation of the battery's capacity requires the knowledge of the number of days of autonomy, which differs from site to another according to weather conditions. The objective is to supply the electricity to the building from batteries when there is no production. In our case, the number of days of autonomy fluctuates between one and three days. In our computations, the value used is one day. A battery of 12 V with a DoD equal to 60% and a battery's efficiency of 85% has been chosen. 70. 26 Ah was found to be the required capacity to install. As a result, we have chosen to work with one battery of 12 V and 75 Ah manufactured by "Electra". Besides, the electric grid is connected to the system through a rectifier that provides 12 V as output voltage with the current needed by the ventilation system.
In order to prove the accuracy of the sizing method, several experiments and simulations have been carried out. However, for simulating the whole system, the characterization of the battery is required. After identifying the battery's parameters that are shown in Table 2, the battery's model has been connected to the other components' models and several simulations and experiments have been carried out. The simulation and experimental results illustrated in Figure 10 demonstrate that during the night, the electricity needed by the ventilation system is supplied totally by the battery and there is no switching to the electric grid. These results prove the accuracy of the sizing method and also the efficiency of the used battery's model and the characterization methodology.   Figure 10 demonstrate that during the night, the electricity needed by the ventilation system is supplied totally by the battery and there is no switching to the electric grid. These results prove the accuracy of the sizing method and also the efficiency of the used battery's model and the characterization methodology.

Conclusions and Perspectives
In this work, a MG system has first been deployed in our EEBLab, and experimental results revealed the power control's behavior, which allows switching between the PV-battery system and the electric grid. Also, these results demonstrate that the electricity needed to operate the ventilation system is not totally produced by the deployed MG system. For that, a sizing method based on the consumption profile of the ventilation system and the weather conditions of our site, has been introduced in order to upgrade the MG system for the aim to produce the total electricity needed by the load. The results prove the accuracy of this method for sizing this type of systems. Besides, in order to simulate the MG system for the aim to compare the simulation to the experiment, a model of each component of the system is needed. For that, all the components of the system have been

Conclusions and Perspectives
In this work, a MG system has first been deployed in our EEBLab, and experimental results revealed the power control's behavior, which allows switching between the PV-battery system and the electric grid. Also, these results demonstrate that the electricity needed to operate the ventilation system is not totally produced by the deployed MG system. For that, a sizing method based on the consumption profile of the ventilation system and the weather conditions of our site, has been introduced in order to upgrade the MG system for the aim to produce the total electricity needed by the load. The results prove the accuracy of this method for sizing this type of systems. Besides, in order to simulate the MG system for the aim to compare the simulation to the experiment, a model of each component of the system is needed. For that, all the components of the system have been modeled. In this work, studying the behavior of the batteries within the MG system was our main focus. Therefore, their characterization must be determined in order to model them properly. As a result, a platform composed of sensing/actuating equipment was firstly developed and deployed in order to gather several battery's characteristics (e.g., voltage, current) that were used to build the battery's model. Simulations of the MG system have been conducted under the same experimental conditions, and the results compared to the experimental ones showed similar behavior. Furthermore, the obtained results showed the accuracy of the battery characterization methodology. Therefore, our ongoing work concentrates on the improvement of the battery's model by taking into consideration the influence of the temperature on the battery's properties. Also, we focus on the battery's SoC forecasting, which is necessary for establishing predictive control approaches in order to manage in real-time our MG system [32].