2.2. Literature Review
The carbon emission of thermal plants (Diesel power plant), including systems, exclude the plant in the literature survey. In any case, the benefit is mainly GHGs reduction, which is unsuitable energy storage compared to this work. The location of renewable energy generation site was limited to land and onshore due to the city. Offshore wind or floating PV were possible to transfer electricity, but the additional infrastructure cost could be added. There are many literature reviews including “Homer”, “hydrogen”, “Fuel cell” “electrolyzer”, “PV”, and “wind” in previous bibliometric analysis. By Scopus, total19 literatures were listed and then excluded the system containing diesel plant, and the finalized as six papers (17-22) are tabulated in Table 2.
Table 1.
Summary of publications conducting the technoeconomic analysis of energy storage systems including hydrogen using HOMER.
Table 1.
Summary of publications conducting the technoeconomic analysis of energy storage systems including hydrogen using HOMER.
Ref
|
Item |
COE |
NPC |
PV |
Wind |
BESS |
PEMWE |
PEMFC |
H2 tank |
Region |
AC load or Hydrogen demand |
Objective |
$ kWh−1
|
$ |
kW |
kW |
kW |
kW |
kW |
kg |
|
kWh day−1 kgH2 day−1
|
[17] |
BESS vs. HydESS comparison |
1.208 |
75,428 |
4 |
2 |
0 |
3 |
2 |
1 |
Saudi Arabia |
14 kWh day−1
|
[18] |
BESS vs. HydESS comparison |
0.52 |
7,160,000 |
469 |
1500 |
396 |
350 |
100 |
500 |
USA |
70 kgH2 day−1
|
[19] |
HydESS feasibility |
0.8399 |
1,006,293 |
120 |
- |
- |
60 |
13 |
20 |
Spain |
200 kWh day−1
|
[20] |
HydESS case study |
1.901 |
910,415 |
1 |
1 |
- |
150 |
1 |
100 |
Oman |
12 kWh day−1
|
[21] |
PV vs. Wind comparison |
0.374 |
2,990,000 |
69 |
88 |
581 |
250 |
450 |
700 |
Canada |
125 kgH2 day−1
|
[22] |
PV vs. Wind comparison |
0.366 |
5,276,069 |
900 |
700 |
1088 |
300-kW |
100 |
300 |
Korea |
2187.6 kWh day−1
|
This work |
BESS vs. HydESS comparison |
0.245 |
4,508,672 |
0 |
1,000 |
441 |
250 |
250 |
700 |
Ukraine |
3,900 kWh day−1
|
Al-Sharafi et al. [
17] evaluated the economic factor, LCOE as 1.208
$ kWh
−1 in Abha area for the hydrogen production storage system, In the reference, the result of System #6 exhibits the minimum LCOE with 1.208
$ kWh
−1 and LCOH is 43.3
$ kg
−1. The optimum configuration has 4 kW PV array, two wind turbines, 2 kW converter, 2 kW fuel cell, 3 kW electrolyzer and 1 kg hydrogen tank. Another case by
Abdin et al. [
18] calculates the minimum COE as 0.5
$ kWh
−1 for Golden Colorado equipped with 400 kg H
2 tank.
Peláez-Peláez et al. [
19] simulated the LOCE as 0.8399
$ kWh
−1 containing 20kg tank with PV coupled electrolyzer and fuel cell.
Okonkwo et al. [
20] conceptualized the system as PV, wind, fuel cell and hydrogen tank with LCOE, 1.901
$ kWh
−1 in Oman and the hydrogen storage capacity was set to 100 kg tank.
Babaei et al. [
21] acquired LCOE 0.374
$ kWh
−1 with the hydrogen based microgrid composed with PV, wind, fuel cell, hydrogen tank, and batteries in Saint Pierre Island caseRecently,
Zhang et al. [
22] optimized the off-grid hybrid energy system in Ui island. The system was composed of 990-kW PV panels, 700-kW wind turbines, a 1088-kWh Li-ion battery bank, 534-kW converter, 300-kW PEMWE system, 300-kg hydrogen tank, and 100-kW PEMFC. The total NPC of the system is
$5,276,069, and the LCOE is 0.366
$ kWh
−1. This case study proved that the self-sufficient energy system in isolated area for community demand without any auxillary power supply, and the hypotheses of the simulation results are limited by the low capacity factor of electrolyser. And the system showed the good feasibility compared to diesel plant, but the fuel cell was adopted as PEM type with 40,000 lifetime span. Moreover, the Battery cost is increased within a few years up to 596
$ kWh [
23].
2.3. System and site outline
In this study, the optimum specifications of water electrolysis capacity and other ancillary devices were derived based on the system’s actual operation results, including the 1,000 kW-class wind turbine power generation. As simulation results, the NPC and the levelized cost of electricity (LCOE) were calculated. The optimized system result is resulted in
Figure 2 from system (a) to system (c). The difference between system (a) and (b) is photovoltaic panel and the difference between (a, b) and (c) is HydESS usage. The most optimized system is (a), and it will be discussed in result and discussion section.
Figure 2.
Basic hybrid hydrogen energy storage system (HydESS) and battery energy storage system(BESS) configuration with renewable energy sources, PV, wind turbine.
Figure 2.
Basic hybrid hydrogen energy storage system (HydESS) and battery energy storage system(BESS) configuration with renewable energy sources, PV, wind turbine.
Figure 3.
Location of Igor Sikorsky Kyiv polytechnic institute in Kyiv, Ukraine (50°40.0’N, 30°59.5’E), (a) load view, (b) arial view, and (c) Igor Sikorsky Kyiv polytechnic institute building.
Figure 3.
Location of Igor Sikorsky Kyiv polytechnic institute in Kyiv, Ukraine (50°40.0’N, 30°59.5’E), (a) load view, (b) arial view, and (c) Igor Sikorsky Kyiv polytechnic institute building.
The capital city of Ukraine, Kyiv (50°40.0’N, 30°59.5’E) is studied as the subject of the analysis in
Figure 4. This study limited on local area point is technoeconomic comparison between HydESS and BESS for renewable energy storage in Kyiv, capital of Ukraine, as center of culture and economy. The NASA satellite data was used to obtain statistics on average solar radiation and wind speed during a one-year period. According to
Figure 4a, the yearly average solar radiation on Kyiv is 3.10 kWh m
−2day
−1. Additionally, clearness nearly cloudy independent on the season, from 0.367 in November and 0.496 in Aug. Referring to
Figure 4b, the average wind speed was around 6.20 m s
−1 when the wind turbine was erected at the height of 65 m above sea level. Temperature is coldest in January as −6.14 ℃ and hottest as 20.83 ℃ in July as shown in
Figure 4c.
Figure 4.
Climate profile of Kyiv, (a) daily solar radiation and clearness index per month, (b) wind speed, and (c) temperature of Kyiv, Igor Sikorsky Kyiv Polytechnic Institute location in Ukraine.
Figure 4.
Climate profile of Kyiv, (a) daily solar radiation and clearness index per month, (b) wind speed, and (c) temperature of Kyiv, Igor Sikorsky Kyiv Polytechnic Institute location in Ukraine.
Figure 5.
AC Primary load profiles of Kyiv, (a) daily, (b) monthly, and (c) yearly of Kyiv, Igor Sikorsky Kyiv Polytechnic Institute location in Ukraine.
Figure 5.
AC Primary load profiles of Kyiv, (a) daily, (b) monthly, and (c) yearly of Kyiv, Igor Sikorsky Kyiv Polytechnic Institute location in Ukraine.
Power demand was set as a residence load provided as a sample case from HOMER to select the AC (alternating current) electric load scale size. As shown in
Figure 4, the total daily load of 3,900 kWh day
−1, and the maximum peak load is 482.3 kW and load factor is 0.34. This AC load is calculated from the average electricity consumption per capita data and capacity of dormitory building . Average power (electricity) consumption per capita in Ukraine is 2,844 kWh year
−1 , 7.8 kWh
−1 according to world data [
24]. Recently opened dormitory in Igor Sikorsky Kyiv Polytechnic Institute about 500 students and staffs are capable to residence in 2021 [
25].
Table 3.
Energy storage system components cost, technical data, and lifetime.
Table 3.
Energy storage system components cost, technical data, and lifetime.
System data |
|
|
PV Unit |
|
|
Capital expenditure |
1252 |
$ kW−1
|
Operation and maintenance expenditure |
18 |
$ y−1 kW−1
|
Lifetime |
25 |
years |
Derating factor |
85 |
% |
Wind turbine system |
|
|
Capital expenditure |
1500 |
$ kW−1
|
Operation and maintenance expenditure |
50 |
$ y−1 kW−1
|
Lifetime |
25 |
years |
Batteries (100 kWh Li BESS) |
|
|
Nominal voltage |
480 |
V |
Nominal catacity |
90.7 |
kWh |
Maximum capacity |
189 |
Ah |
Captial expenditure |
415 |
$ kWh−1
|
Replacement expenditure |
280 |
$ kWh−1
|
Operation and maintenance expenditure |
25 |
$ y−1 kWh−1
|
Lifetime |
10 |
years |
Converter |
|
|
Capital expenditure |
90 |
$ kW−1
|
Replacement c expenditure |
90 |
$ kW−1
|
Operation and maintenance expenditure |
0 |
$ y−1
|
Lifetime |
15 |
years |
Efficiency |
95 |
% |
Electrolyzer |
|
|
Capital expenditure |
1400 |
$ kW−1
|
Replacement expenditure |
650 |
$ kW−1
|
Operation and maintenance expenditure |
60 |
$ y−1
|
Lifetime |
25 |
years |
Storage tank (High pressure) |
|
|
Capital expenditure |
635 |
$ kg−1
|
Operation and maintenance expenditure |
3 |
$ y−1
|
Lifetime |
25 |
years |
Fuel cell |
|
|
Capital expenditure |
4800 |
$ kW−1
|
Replacement expenditure |
3800 |
$ kW−1
|
Operation and maintenance expenditure |
|
$ y−1
|
Lifetime |
5 |
years |
The cost information of system components is listed in Table 2. The cost of new PV is estimated as
$ 1,252 kW
−1 and the annual OPEX is estimated as18
$ kW
−1 from the reference [
26]. The panel lifetime is expected to be 25 years and the derating factor is set to 85%. The wind turbine is the Leitwind LTW77 1,000 kW model used in the existing literature, and the cut in speed is around 11~12 m s
−1 [
27]. It can be seen from
Figure 6; the manufacturer provided the wind turbine’s power curve and the capital expenditure (CAPEX) is
$ 1,500 kW
−1. The operation and maintenance expenditure (OPEX) cost is annually
$ 50 kW
−1 and the total lifetime is 25 years.
Figure 6.
Wind turbine (a) Leitwind AG / Leitwind LTW77 / 1000 kW wind turbine generator / Melfi, Italy (maker‘s web site), (b) Power curve of wind turbine LWT77 from Leitwind 1,000 kW model.
Figure 6.
Wind turbine (a) Leitwind AG / Leitwind LTW77 / 1000 kW wind turbine generator / Melfi, Italy (maker‘s web site), (b) Power curve of wind turbine LWT77 from Leitwind 1,000 kW model.
BESS adopted the Li ion battery such as LG Energy Solution, the 10-kWh class model [
23]. The nominal voltage and capacity are 48 V and 9.8 kWh and the maximum capacity is 189 Ah. CAPEX is
$ 596 kWh
−1 [
28] and replacement expenditure is
$ 280 kWh
−1, O&M cost is
$ 23.3 kWh
−1 [
29]. The battery lifespan is 10 years [
30]. The CAPEX of converter is
$ 80 kW
−1 form NREL’s report for grid energy storage technology cost [
29]. The efficiency is 95% and the lifetime of converter is expected to 15 years. The converter system and installation cost is estimated as 90
$ kW
−1 and the replacement cost is identical to CAPEX. Operation and mainternance cost of the converter is 0 in practically according to other research [
31].
The water electrolyzer is adopted as PEM type in this study, easily operable compared to other electrolysis, counting on the relatively high capital expenditure as 1,400
$ kW
−1 [
32] and the OPEX is estimated as 2% of CAPEX [
33]. Solid oxide fuel cell (SOFC) is chosen for its highest efficiency and adaptability for stationary applications. The specific power generation is 27.52 kWh kg
−1H
2, if the voltage efficiency of SOFC is 79.3% from the reference [
34]. And the specific fuel consumption is calculated in HOMER following the procedure [
35]. The capital cost and maintenance expenditures become more higher c.a. 4,800
$ kW
−1 [
36,
37] than that of PEM fuel cell, 3,300
$ kW
−1 [
22]. The capital and maintenance cost of Hydrogen tank is referenced by NREL’s previous report [
38]. Life time span of SOFC is based on the literature as 40,000 hrs [
39,
40] and expected optimum value as 50,000 hrs [
41].
2.4. Calculation
Homer software simulate the cost-effective system configuration by calculating the minimum NPC. The NPC is defined as total annual cost during project years normalizing with unit cost recovery factor as presented in Eq. 1. These definitions and calculations of indices for technoeconomic analysis are referenced from Homer manual [
32].
The capital recovery factor (CRF) was calculated by Eq. 2 suggested by handbook and manual [
42,
43] Using an interest rate i, auuity, n , the capital recovery factor is :
Annual cost of plant is reduced as the facility is operated considering the depreciation, where d denotes discount rate,
where n is the number of annuities received. This is related to the annuity formula, which gives the present value in terms of the annuity, the interest rate, and the number of annuities. The capital recovery factor is a fraction determined as indicated above where overnight capital cost is expressed as dollars per installed kilowatt (
$ kW
−1). O&M expenses are divided into fixed costs per kilowatt-year (
$ kW
−1-yr) and variable costs per kilowatt-hour (
$ kW
−1). In the denominator 8,760 is the number of hours in a year and capacity factor is a fraction between 0 and 1 representing the portion of a year that the power plant is generating power. Fuel cost is optional since some generating technologies like solar and wind do not have fuel costs. In this work, the fuel cost is only applied in the diesel BESS comparison case, while other HydESS systems use hydrogen as fuel cost, 0.
Levelized Cost of Energy (LCOE, also called Levelized Energy Cost or LEC) is a cost of generating energy (usually electricity) for a particular system. It is an economic assessment of the cost of the energy-generating system including all the costs over its lifetime: initial investment, operations and maintenance, cost of fuel, cost of capital. A NPC calculation is performed and solved in such a way that for the value of the LCOE chosen, the project’s NPC becomes zero. The following formula is used to compute the levelized simple cost of energy. In this case, the LCOE follows the off grid system type without disel generator according to other reference [
44].
The LCOH was calculated by dividing the yearly capital expenditure (CAPEX) and the operating expenditure (OPEX) cost of the system by the amount of hydrogen production, as illustrated in the following equation[
45]. LCOH[
18] in HOMER was where C
ann tot , is the total annualized cost, v
elec is the value of electricity, E
prim is the primary electrical load, E
def is the deferrable load, E
grid sales ,is the total energy sold to the grid and M
H2 is the total hydrogen production.
The efficiency of the water electrolysis system (ε
WE) was calculated as follows. As of this time, the HHV (High heating value of hydrogen, 39.4 kWh kg H
2−1) was 1.48 V.
As shown in
Figure 5, the system configuration was evaluated for a water electrolysis system connected with a renewable energy source. With an assumption, a hydrogen storage tank (1,000 kg capacity) was added for the minimum system requirement, and a lead-acid battery was added to store surplus power and supply it back to the system. The cost of a hydrogen tank is obtained from [
46].
The system was designed with size as feasible as renewable energy without compromising the streaming factor for 1 MWh water electrolysis. The power for water electrolysis was supplied to the electrolyzer through an inverter from the DC solar panel in series and a converter of AC wind turbine in parallel. As a result, the minimum operating capacity for water electrolysis was deduced as 10%, and the voltage efficiency of the water electrolysis system was set to 80%. The performance of the water electrolysis system was calculated through voltage efficiency, and the calculation formula was as follows. The 80% system performance means that the stack driving voltage was 1.85 V at the current target density. In this circumstance, the power consumption was 49.3 kWh kgH
2−1.
The economic analysis model combined the renewable energy sources, water electrolysis, and hydrogen refueling stations was entered by referring to the evaluation setting values in the previous literature [
18]. Each system unit cost and lifespan were entered into the software, as shown below. As shown in
Table 1, the solar panel specifications, wind turbine, battery, and converter to run the system were applied. Based on these fundamental unit values, in the existing system, the result of direct optimization in the system was used as it is without adding a scaled factor or setting the capacity. The inflation rate and discount rate are set to 2% and 8%, respectively, as adopted in other researches [
47,
48].