Optimum State-of-Charge Operating Range for Frequency Regulation of Energy Storage Systems Using a Master–Slave Parallel Genetic Algorithm

: Lithium batteries are used for frequency regulation in power systems because of their fast response and high e ﬃ ciency. Lithium batteries have di ﬀ erent life characteristics depending on their type, and it is necessary to set the optimal state-of-charge (SOC) operating range considering these characteristics to obtain the maximum gain. In general, narrowing the operating range increases the service life but may lower the performance of charging and discharging operations in response to frequency ﬂuctuations, and vice versa. We present performance assessment indicators that consider charging and discharging due to frequency variations and lifespan of the batteries. However, to evaluate the performance, while reﬂecting the non-linear life characteristics of lithium batteries, simulating the entire operation is necessary, which requires a long calculation time. Therefore, we propose a master–slave parallel genetic algorithm to derive the optimal SOC operating range with reduced calculation time. A simulation program was implemented to evaluate the computational performance that determines the optimal SOC range. The proposed method reduces the calculation time while considering the non-linear life characteristics of lithium batteries. It was conﬁrmed that a more accurate SOC operating range could be calculated by simulating the entire life span.


Introduction
A power system must maintain a balance between power generation and load to sustain a stable frequency. In a rotary power system, the rotational speed of the generator increases as the load decreases, and vice versa. This balance is altered by changes in the power system frequency. Therefore, the power system must adjust the generator output to rectify the frequency and balance the load and power generation [1,2].
Frequency regulation can be achieved in two ways. The first is via direct measurement of the power system frequency using the generator or an energy storage device and varying the output in proportion to the difference between the rated and measured frequency. This is called primary frequency control (PFC) or droop control. The second is to control the output fluctuations by measuring the frequency in the energy management system (EMS) that operates the power system and communicates with the generator or energy storage device considering the cost and efficiency of the output from each generator. This is called secondary frequency control (SFC) or automatic generation control (AGC). PFC has a fast response time because the generator changes the output by measuring the frequency directly. SFC, however, can reduce the total output variation cost because it distributes the output according to the generation cost across the generator [1]. in 1 s units has the disadvantage that it requires very long calculation time. This method that takes a long calculation time is very inappropriate to apply to various ESS-FR optimal designs. Therefore, we propose a method of setting the optimal value of the SOC operating range of an ESS-FR considering nonlinear lifetime characteristics via a master-slave parallel genetic algorithm (MSPGA). To reduce the computation time, we leverage parallel computing and genetic algorithms to optimally utilize the computational resources to calculate the optimal SOC operating range. The proposed method can calculate the SOC's operating range by considering the nonlinear lifetime characteristics of the ESS-FR by reducing the computation time.
In Section 2, the configuration and operation method of the ESS-FR and the battery life deterioration model were described. In Section 3, MSPGA for setting the optimal SOC range was described. Section 4 analyzed various cases. In our case study, we verified the optimal SOC setting using the parallel genetic algorithm. We analyzed the difference between designing with the full life and design considering the partial life of the energy storage device for frequency adjustment. The reduction in the calculation time was confirmed by parallel computation. In addition, it was confirmed that it is useful to consider individual life characteristics in a system composed of a combination of batteries having various life characteristics.

System Configuration
An ESS-FR is built with capacity ranging from tens of megawatts to hundreds of megawatts. Figure 1 shows the configuration of an ESS-FR installed by KEPCO (Korea Electric Power Corporation). KEPCO has installed and operated ESS-FRs with a total capacity of 376 MW at 13 substations in South Korea [23,24]. Each energy storage unit consists of a 2 MW inverter and a 750 kWh lithium battery. Such units are connected in parallel to a 23 kV bus via a transformer. In general, there are 12 units with a total capacity of 24 MW. These 12 units are generally supplied by two manufacturers and are classified into Group A and B. The lithium batteries in different groups have different lifespan degradation characteristics because of the different electrode materials and battery manufacturers.
This method that takes a long calculation time is very inappropriate to apply to various ESS-FR optimal designs. Therefore, we propose a method of setting the optimal value of the SOC operating range of an ESS-FR considering nonlinear lifetime characteristics via a master-slave parallel genetic algorithm (MSPGA). To reduce the computation time, we leverage parallel computing and genetic algorithms to optimally utilize the computational resources to calculate the optimal SOC operating range. The proposed method can calculate the SOC's operating range by considering the nonlinear lifetime characteristics of the ESS-FR by reducing the computation time.
In Section 2, the configuration and operation method of the ESS-FR and the battery life deterioration model were described. In Section 3, MSPGA for setting the optimal SOC range was described. Section 4 analyzed various cases. In our case study, we verified the optimal SOC setting using the parallel genetic algorithm. We analyzed the difference between designing with the full life and design considering the partial life of the energy storage device for frequency adjustment. The reduction in the calculation time was confirmed by parallel computation. In addition, it was confirmed that it is useful to consider individual life characteristics in a system composed of a combination of batteries having various life characteristics.

System Configuration
An ESS-FR is built with capacity ranging from tens of megawatts to hundreds of megawatts. Figure 1 shows the configuration of an ESS-FR installed by KEPCO (Korea Electric Power Corporation). KEPCO has installed and operated ESS-FRs with a total capacity of 376 MW at 13 substations in South Korea [23,24]. Each energy storage unit consists of a 2 MW inverter and a 750 kWh lithium battery. Such units are connected in parallel to a 23 kV bus via a transformer. In general, there are 12 units with a total capacity of 24 MW. These 12 units are generally supplied by two manufacturers and are classified into Group A and B. The lithium batteries in different groups have different lifespan degradation characteristics because of the different electrode materials and battery manufacturers. The power management system (PMS) controlling the ESS-FR measures the voltage frequency and receives the SOC of the battery system from the battery management system (BMS). The charging/discharging power is then calculated by an internal algorithm, and commands are sent to the inverters of the energy storage units. Figure 2 shows the deployment of a 24 MW ESS-FR constructed with 18 containers in KEPCO. Six of them have PCS and 12 have lithium batteries. The ESS-FR was partially built by two different manufacturers.  The power management system (PMS) controlling the ESS-FR measures the voltage frequency and receives the SOC of the battery system from the battery management system (BMS). The charging/discharging power is then calculated by an internal algorithm, and commands are sent to the inverters of the energy storage units. Figure 2 shows the deployment of a 24 MW ESS-FR constructed with 18 containers in KEPCO. Six of them have PCS and 12 have lithium batteries. The ESS-FR was partially built by two different manufacturers.

Primary Frequency Control
The ESS-FR is input with the frequency and the SOC, as shown in Equation (1), and the output ( ) is calculated as follows: First, if the frequency ( ) measured in the ith simulation using Equation (2) is maintained within the rated frequency ( ) of 60 ± 0.03 Hz, the charging/discharging operation is performed to maintain the SOC. Otherwise, the charging/discharging operation is performed to regulate the frequency. At this time, the corresponding output quantities are and , which are calculated as shown in Equations (3) and (4), respectively.
For the output to maintain the SOC ( ), if the state-of-charge ( ) is higher or lower than the neutral point by over 2%, the charging/discharging operation is performed with respect to the output amount, which is the rated output of the energy storage device ( ) multiplied by the SOC recovery rate ( ), where a positive number indicates a discharge operation, and a negative number indicates a charging operation.
The charging/discharging amount ( ) for frequency regulation is calculated using the same formula for the output of the primary frequency control of a typical generator [1], where K is the droop constant. If is greater than the maximum operating value ( ), charging operation cannot be performed even with a high frequency. Similarly, if is lower than the minimum operating value ( ), discharging operation cannot be performed even with a low frequency.
When the output of the energy storage device ( ) is calculated, the state-of-charge in the next step is updated according to Equation (5), where is the hourly simulation step and is the charging/discharging efficiency. Equation (5) is used for the simulation of the ESS, but the actual system uses the SOC measured in the BMS, as shown in Figure 1.
is the storage capacity of the energy storage device available in the ith simulation step and is calculated using Equation (6). is calculated by considering the degradation rate ( ) for the initial rated storage capacity ( ) of the energy storage device. The degradation rate is calculated using the charging/discharging history of the energy storage device and the lithium battery degradation model. Section 2.3 describes it in detail.

Primary Frequency Control
The ESS-FR is input with the frequency and the SOC, as shown in Equation (1), and the output (P i ESS ) is calculated as follows: First, if the frequency ( f i in ) measured in the ith simulation using Equation (2) is maintained within the rated frequency ( f rated ) of 60 ± 0.03 Hz, the charging/discharging operation is performed to maintain the SOC. Otherwise, the charging/discharging operation is performed to regulate the frequency. At this time, the corresponding output quantities are P i SOC and P i dr , which are calculated as shown in Equations (3) and (4), respectively.
For the output to maintain the SOC (P i SOC ), if the state-of-charge (S i ) is higher or lower than the neutral point by over 2%, the charging/discharging operation is performed with respect to the output amount, which is the rated output of the energy storage device (P r ) multiplied by the SOC recovery rate (r), where a positive number indicates a discharge operation, and a negative number indicates a charging operation.
The charging/discharging amount (P i dr ) for frequency regulation is calculated using the same formula for the output of the primary frequency control of a typical generator [1], where K is the droop constant. If S i is greater than the maximum operating value (S max ), charging operation cannot be performed even with a high frequency. Similarly, if S i is lower than the minimum operating value (S min ), discharging operation cannot be performed even with a low frequency.
When the output of the energy storage device (P i ESS ) is calculated, the state-of-charge in the next step is updated according to Equation (5), where T s is the hourly simulation step and η is the charging/discharging efficiency. Equation (5) is used for the simulation of the ESS, but the actual system uses the SOC measured in the BMS, as shown in Figure 1. C i is the storage capacity of the energy storage device available in the ith simulation step and is calculated using Equation (6). C i is calculated by considering the degradation rate (D i ) for the initial rated storage capacity (E rated ) of the Electronics 2020, 9, 1298 5 of 17 energy storage device. The degradation rate is calculated using the charging/discharging history of the energy storage device and the lithium battery degradation model. Section 2.3 describes it in detail. Table 1 shows the basic parameters for calculating the output of the frequency-adjusted ESS summarized in Equations (1)- (6). In this study, we would like to propose a method for achieving optimal effects by optimally setting S max , S min , and S n , considering the life model of lithium batteries.

Battery-Based Energy Storage Systems Degradation
As a lithium battery performs charging and discharging operations, its lifespan deteriorates according to changes in the SOC. This lifespan degradation is largely caused by two factors. The first is calendar degradation, which occurs over time even without the charging and discharging operations. The second is cycle degradation, which occurs because of the charging and discharging operations. Therefore, the degradation rate (D i ) of the energy storage device is calculated as shown in Equation (7). If the output of the energy storage device calculated using Equation (1) does not involve charging/discharging operations, calendar degradation is reflected; otherwise, cycle degradation is reflected.
Equation (8) is the formula for calculating the calendar degradation rate. The degradation rate during the maintenance for 1 s is calculated using the current degradation rate (D i ) and the state-of-charge (S i ). All parameters are recorded in Table 2. Cycle degradation is calculated using S i , which is changed when a charging or discharging operation is performed. Figure 3 shows the state-of-charge during the charging operation, where S swing is the magnitude of the change in S i , and S avr is the mean value of the changes in S i . Therefore, after all charging and discharging operations, it is possible to calculate the cycle degradation rate using Equation (9).
Cycle degradation is calculated using S i , which is changed when a charging or discharging operation is performed. Figure 3 shows the state-of-charge during the charging operation, where Sswing is the magnitude of the change in S i , and Savr is the mean value of the changes in S i . Therefore, after all charging and discharging operations, it is possible to calculate the cycle degradation rate using Equation (9).  We used three lifespan degradation models in the analysis for setting the optimal SOC operating range of the ESS-FRs with different lifespan degradation characteristics. The first was a lifespan model using the experimental data of an iron phosphate battery [22]. The arbitrary parameters were modified for the rest of the models to utilize various degradation characteristics. Degradation model 2 was set to have characteristics such that calendar degradation progressed faster than that in model 1. The parameters of model 3 were modified, so that cycle degradation proceeded faster than that in model 1. Figure 4a shows the calendar degradation model when S i is maintained at 90%. Figure 4b shows the cycle degradation model with charging and discharging at 50% S avr and 80% S swing , respectively.
Electronics 2020, 9, x FOR PEER REVIEW 6 of 17 We used three lifespan degradation models in the analysis for setting the optimal SOC operating range of the ESS-FRs with different lifespan degradation characteristics. The first was a lifespan model using the experimental data of an iron phosphate battery [22]. The arbitrary parameters were modified for the rest of the models to utilize various degradation characteristics. Degradation model 2 was set to have characteristics such that calendar degradation progressed faster than that in model 1. The parameters of model 3 were modified, so that cycle degradation proceeded faster than that in model 1.  Figure 4a shows the calendar degradation model when S i is maintained at 90%. Figure 4b shows the cycle degradation model with charging and discharging at 50% Savr and 80% Sswing, respectively.  When reaches 20% for a lithium battery, only 80% of the initially installed capacity can be charged and discharged and is defined as the end of the battery lifespan. This is because the capacity that can be charged and discharged decreases rapidly when exceeds 20% [25]. In this paper, the life model of the battery derived from the experimental results was applied on an iron phosphate battery. Because we calculate the objective function based on simulation, the method presented in this paper is suitable for any life model.

Objective Function
We calculated and used three types of energy values to define the lifespan objective function based on the rules of the electricity market or utility use of the ESS-FR. For the first type, the amount of energy charged/discharged for frequency regulation (Efreq) was calculated by applying a weight to the output charged/discharged by Pdr in Equation (1), as shown in Equation (10). Figure 5 shows the weights by frequency. A higher weight was applied to the energy charged/discharged for higher frequency fluctuations. The weight was set to have a value less than 1 when the frequency was maintained within a range (60 ± 0.2 Hz) and to have a proportionally large weight when the frequency was outside the specified range. When D i reaches 20% for a lithium battery, only 80% of the initially installed capacity can be charged and discharged and is defined as the end of the battery lifespan. This is because the capacity that can be charged and discharged decreases rapidly when D i exceeds 20% [25].
In this paper, the life model of the battery derived from the experimental results was applied on an iron phosphate battery. Because we calculate the objective function based on simulation, the method presented in this paper is suitable for any life model.

Objective Function
We calculated and used three types of energy values to define the lifespan objective function based on the rules of the electricity market or utility use of the ESS-FR. For the first type, the amount of energy charged/discharged for frequency regulation (E freq ) was calculated by applying a weight to the output charged/discharged by P dr in Equation (1), as shown in Equation (10). Figure 5 shows the weights by frequency. A higher weight was applied to the energy charged/discharged for higher frequency fluctuations. The weight was set to have a value less than 1 when the frequency was maintained within a range (60 ± 0.2 Hz) and to have a proportionally large weight when the frequency was outside the specified range.
The second item in the objective function was the amount of energy charged/discharged by ESOC calculated for the entire lifetime, as shown in Equation (11): The last item was the amount of energy that could not be output because of the restrictions of Smax and Smin in Equation (4), calculated as shown in Equation (12): The final objective function was defined as the sum of the three items previously calculated, as shown in Equation (13): Smax, Smin, and Sn utilized in Equations (3) and (4) were the input variables of the optimization problem, and the restrictions are summarized as follows: To calculate the objective function in 1 s unit increments over the entire lifetime of the ESS-FR, it is necessary to repeat Equations (10)-(13) n times and accumulate the results. Here, n is calculated as follows: = 60 × 60 × 24 × 365 × = 31,536,000 × here, the life is the life of the ESS-FR described in Section 2.3 and is generally 15 years to 21 years, so the total number of calculation (n) is about 473,040,000 to 630,720,000. This is a very large calculation amount, so it takes a long calculation time. If a general genetic algorithm is used, it takes a very long computation time until the final convergence, with many iterations of the objective function.

Parallel Genetic Algorithm
A genetic algorithm (GA) is an optimization algorithm suitable for combinatorial optimization. Genetic algorithms are used in various fields because of their simplicity of implementation and their ability to search for global optimal solutions [26]. In this study, we applied a genetic algorithm to set the operating range for the optimal SOC of the ESS-FR. The ESS-FR simulates the entire life and calculates the objective function to accurately reflect the life characteristics of lithium batteries. In general, the life of the ESS-FR is more than 10 years, and the simulation step for calculating the objective function is performed in units of 1 s, which requires considerable calculation time. Therefore, The second item in the objective function was the amount of energy charged/discharged by E SOC calculated for the entire lifetime, as shown in Equation (11): The last item was the amount of energy that could not be output because of the restrictions of S max and S min in Equation (4), calculated as shown in Equation (12): The final objective function was defined as the sum of the three items previously calculated, as shown in Equation (13): S max , S min , and S n utilized in Equations (3) and (4) were the input variables of the optimization problem, and the restrictions are summarized as follows: To calculate the objective function in 1 s unit increments over the entire lifetime of the ESS-FR, it is necessary to repeat Equations (10)-(13) n times and accumulate the results. Here, n is calculated as follows: here, the life is the life of the ESS-FR described in Section 2.3 and is generally 15 years to 21 years, so the total number of calculation (n) is about 473,040,000 to 630,720,000. This is a very large calculation amount, so it takes a long calculation time. If a general genetic algorithm is used, it takes a very long computation time until the final convergence, with many iterations of the objective function.

Parallel Genetic Algorithm
A genetic algorithm (GA) is an optimization algorithm suitable for combinatorial optimization. Genetic algorithms are used in various fields because of their simplicity of implementation and their ability to search for global optimal solutions [26]. In this study, we applied a genetic algorithm to set the operating range for the optimal SOC of the ESS-FR. The ESS-FR simulates the entire life and calculates the objective function to accurately reflect the life characteristics of lithium batteries. In general, the life of the ESS-FR is more than 10 years, and the simulation step for calculating the objective function is performed in units of 1 s, which requires considerable calculation time. Therefore, we used Electronics 2020, 9, 1298 8 of 17 a method that applied the master-slave parallel genetic algorithm (MSPGA) in this study [27,28]. The chromosome of the genetic algorithm was defined as a set of S min , S max , and S n , as shown in Equation (16): Figure 6 shows the computational flow of a typical genetic algorithm. Initially, a large number of chromosomes are generated randomly, as shown in Equation (16), to form a group. The objective function is calculated for each chromosome in the randomly generated group and the functions are then sorted in descending order according to their results. The termination condition may be found to be satisfied upon checking, but if not, the low-rank chromosomes are deleted, and high-rank chromosomes are maintained by elitism.
Electronics 2020, 9, x FOR PEER REVIEW 8 of 17 study [27,28]. The chromosome of the genetic algorithm was defined as a set of Smin, Smax, and Sn, as shown in Equation (16): Figure 6 shows the computational flow of a typical genetic algorithm. Initially, a large number of chromosomes are generated randomly, as shown in Equation (16), to form a group. The objective function is calculated for each chromosome in the randomly generated group and the functions are then sorted in descending order according to their results. The termination condition may be found to be satisfied upon checking, but if not, the low-rank chromosomes are deleted, and high-rank chromosomes are maintained by elitism.  Figure 7 shows the process of crossover and mutation in genetic algorithms. Some individuals with high objective function values are retained in the next generation. The number of individuals maintained by elitism is defined as Nelite. The deleted individuals are generated by the crossover operator and the mutant operator, and the number of individuals is defined as Nco and Nmu. For crossover and mutation, each of the parents is randomly selected from the elite. In the crossover operation, random single-point crossover was applied. Crossover creates two new individuals at a time. In the mutation operation, mutations were randomly generated at random single points. When all individuals are newly created, the objective function calculation and the sorting are repeated. Through such repetitions, one generation would have evolved in terms of genetic algorithms. In this study, the process of evolution was terminated when the objective function values of the best chromosomes did not improve during 10 generations.   Figure 7 shows the process of crossover and mutation in genetic algorithms. Some individuals with high objective function values are retained in the next generation. The number of individuals maintained by elitism is defined as N elite . The deleted individuals are generated by the crossover operator and the mutant operator, and the number of individuals is defined as N co and N mu . For crossover and mutation, each of the parents is randomly selected from the elite. In the crossover operation, random single-point crossover was applied. Crossover creates two new individuals at a time. In the mutation operation, mutations were randomly generated at random single points. When all individuals are newly created, the objective function calculation and the sorting are repeated. Through such repetitions, one generation would have evolved in terms of genetic algorithms. In this study, the process of evolution was terminated when the objective function values of the best chromosomes did not improve during 10 generations.
In this study, we used MSPGA to calculate the objective function, which would otherwise take a long time. The MSPGA can be even applied with only one machine in a memory structure shared with many CPU cores, like a general PC. The detailed steps for the calculation of the objective function in Figure 6 are shown in Figure 8.
When the calculation of the objective function begins, the calculation of the objective function for one chromosome is assigned to one slave worker. In a quad-core CPU, four workers are created. Once the objective function calculation for the allocated chromosome is completed, the slave worker is assigned with the objective function calculation for another chromosome, and the calculation proceeds again. When the calculations of the objective functions for all chromosomes are completed, the parallel calculation process of the objective function ends. The algorithm then proceeds with sorting the chromosomes in descending order based on the objective functions, as shown in Figure 6. At this time, all the operations in Figure 6 are processed by one of the multiple CPU cores, and this core becomes the master core. operation, random single-point crossover was applied. Crossover creates two new individuals at a time. In the mutation operation, mutations were randomly generated at random single points. When all individuals are newly created, the objective function calculation and the sorting are repeated. Through such repetitions, one generation would have evolved in terms of genetic algorithms. In this study, the process of evolution was terminated when the objective function values of the best chromosomes did not improve during 10 generations.  In this study, we used MSPGA to calculate the objective function, which would otherwise take a long time. The MSPGA can be even applied with only one machine in a memory structure shared with many CPU cores, like a general PC. The detailed steps for the calculation of the objective function in Figure 6 are shown in Figure 8.
When the calculation of the objective function begins, the calculation of the objective function for one chromosome is assigned to one slave worker. In a quad-core CPU, four workers are created. Once the objective function calculation for the allocated chromosome is completed, the slave worker is assigned with the objective function calculation for another chromosome, and the calculation proceeds again. When the calculations of the objective functions for all chromosomes are completed, the parallel calculation process of the objective function ends. The algorithm then proceeds with sorting the chromosomes in descending order based on the objective functions, as shown in Figure 6. At this time, all the operations in Figure 6 are processed by one of the multiple CPU cores, and this core becomes the master core. Table 3 shows the pseudocode of MASPGA implemented in MATLAB. The functions provided by MATLAB are shown in brackets. In the pseudocode, the calculation steps shown in the flowcharts in Figures 6 and 8 are briefly shown.   Figure 8. Objectives parallel computing flowchart [29]. Table 3 shows the pseudocode of MASPGA implemented in MATLAB. The functions provided by MATLAB are shown in brackets. In the pseudocode, the calculation steps shown in the flowcharts in Figures 6 and 8 are briefly shown.

Case Studies
In this study, we proposed a method for determining the operating range for the optimal SOC of the energy storage device, to regulate the frequency using MSPGA. We conducted various case studies to verify the effectiveness or convergence of the proposed method. The first case study is for an algorithm parameter calibration experiment. We compared the convergence and computation time according to the change of parameters. In the second case study, we examined the difference between analyzing the overall lifespan of a lithium battery and not in calculating the objective function. In the third case study, the computational speeds of general GA and MSPGA were compared. In the next case study, a brute-force search was performed by shortening the objective function calculation without performing it for the entire lifetime. Lastly, we searched for the optimal SOC setting of ESS-FR, which is composed of a combination of three different battery life characteristics. All simulations were implemented using the parallel computing toolbox of MATLAB, on a PC with an Intel i7-4790 3.6 GHz quad-core CPU and a 16 GB memory.

Genetic Algorithm Parameter Calibration Experiment
When genetic algorithms (GA) are used to determine the optimal SOC range of ESS-FR, the solution quality may be influenced by the population size, the rate of crossover, and the rate of mutation [29]. Choosing a parameter for GA is not the main contribution of this paper, but it was done for the proper setup of the case study conducted in the next section.
We calculated the optimal SOC range by changing the population size. In each case, GA was performed three times and Table 4 summarizes the results. In case 1-1, the small population size showed the fastest computational performance, but the converged values were not consistent. Cases 1-2~1-4 showed consistent convergence results. In this case, the fastest computation time in case 1-2 was shown, so the size of the population was set to 100 in this paper. Table 5 shows the convergence according to the change of the number of individuals newly created by the crossover operator among the entire population. It can be confirmed that the setting of case 2-2 has the highest consistency of converged values. Therefore, in this paper, these set values were applied to all case studies.

Effect of Simulation Duration on Objectives
In Reference [21], an optimal design method for calculating the total lifetime considering the degradation rate over a year was devised by designing the optimal capacity and operation parameters for one year via simulations. However, lithium batteries have nonlinear degradation characteristics. Figure 9 shows the trend in the degradation rate from the results of the simulation, which proceeded until the degradation rate reached 20%, using the parameters in Table 1 and battery model 1 in Table 2. In most cases, calendar degradation accounted for most of the total degradation compared to cycle degradation, resulting in a lifespan of about 21 years. The total degradation rate considering both calendar and cycle degradation was found to be nonlinear. The objective functions for the three random parameters were compared to determine how the lifespan assessment period affects the optimal parameter setting. The three random parameters are summarized in Table 6. The three parameters listed in Table 6 were set arbitrarily, and simply indicate that the optimal parameter may be differently selected due to the difference in the analyzed life span of Table 7. Table 7 shows the objective function value and operating lifespan when each parameter is applied for various degradation rates (D), which is the simulation termination condition. When the simulation proceeded until D reached 1%, Parameter 2 exhibited the highest objective function value. However, when the simulation proceeded until the degradation rate reached 2% to 5%, Parameter 3 exhibited the highest objective function value. When the degradation rate reached 7% to 10%, Parameter 1 exhibited the highest objective function value. Therefore, the accuracy of the optimal design and parameter setting, when the operating benefits of the entire lifetime are considered by simply assessing partial lifetime, cannot be guaranteed. In this study, the objective function was thus calculated as the simulation proceeded until the degradation rate reached 20%, which was generally considered the end-of-life condition for a lithium battery.  We compared the differences between the method used in Reference [21] and the method presented in this paper. In Reference [21], the optimal SOC range was derived during the shortened analysis period, and it was assumed that this SOC range was optimal for the entire life. We compared these approaches by applying the same genetic algorithm. The calculation results of the operating range of the optimum SOC are summarized in Table 8. The objective function results were all calculated to simulate the lifetime. As the proposed method showed higher objective function results, it can be confirmed that it is suitable for optimal SOC operating range.

Computation Time Comparison
Simulations were performed to examine the effect of the parallel computing method on the overall calculation time. Table 9 shows a comparison of the results of three calculations using the parallel genetic algorithm and the general genetic algorithm for the optimal SOC operating range over the lifetime until the degradation rate reaches 20%. When the MSPGA method was used, the simulation of the overall lifetime took 19.2 h on average. However, when the general genetic algorithm was used, it took 62.6 h on average. The proposed MSPGA method was found to reduce the time to derive the results by more than three times. All six results exhibited the same objective functions and operating lives. Only the S max values were different, according to the results. Therefore, it was confirmed that the same result could be acquired faster using MSPGA.  Figure 10 shows the distribution of SOC by setting the scope of SOC operation for the overall lifetime, as shown in the MSPGA case 2 in Table 9. If the frequency is maintained within 60 ± 0.03 Hz, the ESS is set to charge/discharge to maintain the SOC at 43%. Therefore, it can be seen that the SOC distribution is concentrated at 43%. As the frequency rises, the ESS charges and the SOC increases. However, it can be seen that SOC does not exceed 80% over its lifetime. For this reason, it can be confirmed that the objective function value is calculated equally even if S max is set to any value of 80 to 100.  Figure 10 shows the distribution of SOC by setting the scope of SOC operation for the overall lifetime, as shown in the MSPGA case 2 in Table 9. If the frequency is maintained within 60 ± 0.03 Hz, the ESS is set to charge/discharge to maintain the SOC at 43%. Therefore, it can be seen that the SOC distribution is concentrated at 43%. As the frequency rises, the ESS charges and the SOC increases. However, it can be seen that SOC does not exceed 80% over its lifetime. For this reason, it can be confirmed that the objective function value is calculated equally even if Smax is set to any value of 80 to 100.  Table 6.

Convergence Verification
In this study, we proposed a method for estimating the optimal SOC operating range of ESS-FR using MSPGA. We calculated and validated all feasible sets to demonstrate that the proposed MSPGA presents the optimal SOC operating range. To reduce the calculation time, it is assumed that even if the capacity reduction of only 1% is reached, the end of life is assumed, and the battery life model 1 is used. Calculating all feasible sets, even the 20% reduction in degradation capacity takes a long calculation time, which is also a necessity for the MSPGA technique presented. The SOC is assumed to be an integer value from 0 to 100. There are 152,096 cases where three SOC operating range parameters (SOCn, SOCmax, SOCmin) satisfy the constraint equation (Equation (14)). We calculated the objective function for all of these feasible sets, the maximum of which was 553.17 MWh.
Of all feasible sets, only 1000 cases with large objective functions are shown in descending order in Figure 11.
We performed six iterations to confirm that the MSPGA method converged to the same maximum objective value, of 553.17, for each iteration, as shown in Figure 12. Although the convergence speed varied from case to case due to the nature of the genetic algorithm using the random function, the proposed MSPGA was confirmed to converge at a sufficient optimal point.  Table 6.

Convergence Verification
In this study, we proposed a method for estimating the optimal SOC operating range of ESS-FR using MSPGA. We calculated and validated all feasible sets to demonstrate that the proposed MSPGA presents the optimal SOC operating range. To reduce the calculation time, it is assumed that even if the capacity reduction of only 1% is reached, the end of life is assumed, and the battery life model 1 is used. Calculating all feasible sets, even the 20% reduction in degradation capacity takes a long calculation time, which is also a necessity for the MSPGA technique presented. The SOC is assumed to be an integer value from 0 to 100. There are 152,096 cases where three SOC operating range parameters (SOC n , SOC max , SOC min ) satisfy the constraint Equation (Equation (14)). We calculated the objective function for all of these feasible sets, the maximum of which was 553.17 MWh.
Of all feasible sets, only 1000 cases with large objective functions are shown in descending order in Figure 11.  We performed six iterations to confirm that the MSPGA method converged to the same maximum objective value, of 553.17, for each iteration, as shown in Figure 12. Although the convergence speed varied from case to case due to the nature of the genetic algorithm using the random function, the proposed MSPGA was confirmed to converge at a sufficient optimal point.

Usability of SOC Range Optimization Considering Degradation Models
In this study, we proposed a method for setting the optimal SOC operating range using parallel computation and genetic algorithm for the objective function of the lifetime to reflect the nonlinearity of the lithium battery lifespan model. Therefore, when energy storage devices with various degradation characteristics were installed, the differences in the objective functions among the cases in which the SOC operating range was set individually and the cases in which such setting was not made were compared. As shown in Table 10, the difference in the overall objective functions according to the parameter settings was analyzed under the assumption that three energy storage devices with the same capacities but different degradation characteristics were installed.
By using the MSPGA proposed in this study, we could calculate the optimal SOC operating range with shortened calculation time, reflecting all nonlinear lifespan models of lithium battery. We summarized how the overall objective functions were calculated when the three ESS-FRs were set in 4 different cases, as shown in Table 11.

Usability of SOC Range Optimization Considering Degradation Models
In this study, we proposed a method for setting the optimal SOC operating range using parallel computation and genetic algorithm for the objective function of the lifetime to reflect the nonlinearity of the lithium battery lifespan model. Therefore, when energy storage devices with various degradation characteristics were installed, the differences in the objective functions among the cases in which the SOC operating range was set individually and the cases in which such setting was not made were compared. As shown in Table 10, the difference in the overall objective functions according to the parameter settings was analyzed under the assumption that three energy storage devices with the same capacities but different degradation characteristics were installed. By using the MSPGA proposed in this study, we could calculate the optimal SOC operating range with shortened calculation time, reflecting all nonlinear lifespan models of lithium battery. We summarized how the overall objective functions were calculated when the three ESS-FRs were set in 4 different cases, as shown in Table 11. In Case 1, the SOC operating range was set individually with consideration for all lifespan models of the three energy storage devices. In each of the Cases 2 to 4, the SOC operating range was set with consideration for only one lifespan model and was applied to all energy storage devices. The objective function values were compared for the lifetime ending at 20% degradation. Case 1 exhibited the highest objective function value because the SOC operating range was calculated in consideration of the life model of all batteries. Cases 2 to 4 demonstrated lower objective function values than that of Case 1. Differences in the objective function values of Case 1 and Cases 2 to 4 are expected to increase when the differences in the lifespan models of the ESS-FRs are increased.

Conclusions
In this study, we proposed a method for setting optimal SOC operating range of the ESS-FR, considering the nonlinear life model of a lithium battery. The operation of the energy storage system used for frequency regulation was formulated and summarized. Also, the model for simulating the overall degradation rate was summarized by applying the lifespan model of the lithium battery and calculating calendar degradation and cycle degradation. Moreover, we proposed a method for setting the optimal SOC operating range using MSPGA, where parallel computing technology was applied to the genetic algorithm for shortened calculation time, and for setting the optimal operating range. In order to set the optimal SOC range of ESS-FR in consideration of the detailed battery life model, it is necessary to calculate the objective function using a long simulation. General genetic algorithms repeatedly compute the objective function for multiple chromosomes. Therefore, it takes even longer convergence time in the field where the objective function calculation time is long. To overcome this problem, this paper proposed a method to shorten the time to calculate the optimal SOC range by using MSPGA.
In five case studies, we found the following: First, the results showed that the deterioration characteristics of lithium batteries are nonlinear, and it is reasonable to simulate the entire lifetime to set the optimum SOC operating range. Second, we have also found that the parallel computing function can reduce the calculation time while leading to the same result. Third, the exhaust search method was used to verify that the proposed method converges to the global optimum. Finally, when lithium batteries with different degradation characteristics constitute the ESS-FR, it was found that a better result could be obtained by using the optimal SOC operating parameters calculated individually to meet the characteristics of each battery than by applying the parameters set by models that disregarded the differences in the lifetimes. Therefore, it was confirmed that setting the optimal SOC operating range for each type of lithium battery could result in higher profitability when the battery manufacturers were different or when the electrode materials of the lithium batteries were different.
A Li-ESS device requires a high initial installation cost, and its degradation rate tends to vary depending on the operating strategy. Therefore, setting the optimal SOC operating range, as proposed in this study, is an important means to improve the profitability of the ESS-FR. Further research could include finding the optimal operation of different energy storage devices for various uses, including ESS-FR.