3.1. Analysis of Load Profile
Figure 4 depicts the GHU energy consumption profile for a whole year. The energy recordings cover the period from 1/3/2022 until 28/2/2023 with sampling every 15 minutes.
According to the monthly energy consumption requirements, two periods can be distinguished. One period concerns the months October to April, and another that concerns the months May to September. In the first period, the maximum power requirement during the 15-minute recording period shall not exceed 125 kW. During this particular period, the temperature control inside the GHU is mainly done with geothermal energy, so there is not much demand for electrical energy. A typical profile of three consecutive days at this period is shown in
Figure 5.
In the second period, i.e. during the summer months, temperature control is done using fans that operate mainly during daytime. Ιn
Figure 6 the energy consumption profile for the first three days of July is given, which is typical for the entire period of GHU ventilation with fans (Months May to September,
Figure 4). The use of fans takes place between 8:00 and 20:00. The power requirement ranges, for most of this interval, between 200 and 325 kW. For the rest of the day (20:00 – 8:00) the power requirement is reduced and ranges between 50 to 100kW.
Thus, over a year we could distinguish three areas of GHU load demand. The first range from 0 to 125kW, the second range from 125 – 200kW, and the third range from 200 – 350kW. The first area of energy consumption is observed all days of the year. However, there is a difference between the periods noted above. In the period from October to April it is the predominant load demand throughout the day, while in the period May – September the GHU requires power in this area only during the hours from 20:00 – 8:00. The second and third load demand areas occur only in the second period, i.e. from May to September.
The present study focuses on the use of a natural gas or biomethane generator to meet the energy needs of the GHU. The choice of this fuel is based on a recent study by the company concerning its estimation of Greenhouse residues’ biogas production [
13,
14]. The possibility of using constant generator power in the three load-demand areas will be studied. This will be firstly achieved, using three MGs. When the load demand belongs to the first range (0-125kW) one generator will operate. If the load needs belong to the second range (125 – 200kW), two generators will be activated. The generator of the first area, which will operate at maximum power, with one more generator. Finally, when the load requirements appear in the third range (200 – 350kW) then three generators will be activated. The two generators will operate for the first two areas (both at maximum power) and one more. In order to ensure the constant operation of the generators, it is necessary to use an energy storage system. More specifically, using only generators, their rated power should be equal to the maximum power of each of the areas defined above. That is, the first generator had a nominal power of 125kW, the second 75kW (because the second area has an upper limit of 125+75=200kW) and the third 150kW (because the third area has an upper limit of 125+75+150 = 350kW). This means that for long periods the generators would run at very low efficiency to cover very low loads. One way to avoid large fluctuations in generator operation is to select the power of each generator at an intermediate value of each operating area which will remain constant. Using a battery, energy will be absorbed when the load needs are less than the constant power of the generator in the respective operating range. If, on the contrary, the load requirements are higher than the power of the generator, then the battery will be discharged to fully service the load. It becomes obvious that in this way we keep the power of the generator stable and at the highest possible efficiency levels. The objective, therefore, is to identify the optimal combination of MGs and ESS and the optimal mode of operation (in terms of MGs’ efficiency and autonomous operation) of the MG-ESS system to meet the annual energy needs.
3.2. Operation algorithm
Figure 7 shows the operating algorithm of the MG-ESS system. Initially, values are given to the variables shown in
Table 2 which are also the design parameters of the MG-ESS system. These design parameters remain constant until the completion of the algorithm, i.e. after 35040 checks.
As mentioned in a previous chapter three power regions (0
125kW
200kW
350kW) can be distinguished. We define P
GL1 = 125kW and P
GL2 = 200kW, the two power values that separate the first from the second region and the second from the third region, respectively. According to the discussion in Chapter 3.1, the algorithm first seeks the power region where the load demand belongs. When the algorithm detects the operating zone, the value of the demand load power (P
L) is compared with the value of the generator’s operating power (P
Gi where i=1 or 2 or 3 the three operating regions – Design Parameter). If P
L is higher than the P
Gi value, then the battery should contribute to cover the needs of the load together with the generator. However, it should be checked if the battery is able to support the power gap between the P
L and the P
Gi. This is done by checking the battery state of energy (SOE
t) at the time interval t, which is defined by the equation:
With
the energy content of the battery at time t. To discharge the battery should
This limit was set for the safe operation of the battery [
24]. Otherwise (
20%) a backup power source should be used which should work in parallel with the generator to meet the needs of the load.
In case PL is less than the value of PGi then the system exhibits excess energy which should either be absorbed by the battery (battery charging) or should be discarded. The variable that the algorithm calculates to choose one of the two directions is again If 100% then the energy is discarded otherwise 100%) the excess energy is used to charge the battery.
The same process takes place in whichever power range the demand power of the load belongs. So according to
Figure 7, we can distinguish the following operating modes, which can occur in any region:
1st mode – Battery discharge (Stages 1, 5, 9)
When the difference (P
Gi – P
L) is negative and
0% then the battery will be discharged. The amount of energy that will be used by the battery (E
BAT,disch) to meet the needs of the load the time period Δt is:
Where n
inv is the efficiency of the inverter with a value of 0.95. So, the energy required by the load (E
L) at the specified time interval Δt (=15min) will be calculated by:
2nd mode – Energy deficit (Stages 2, 6, 10)
When the difference (P
Gi – P
L) is negative and the
0% then the reserve should be used to provide the required amount of energy (E
Back):
3rd mode – Waste of energy (Stages 3, 7, 11)
When the difference (P
Gi – P
L) is positive and the SOEt
100% then excess energy should be disposed of. The amount of energy to be discarded (E
waste) the time interval Δt is:
4th mode – Battery charging (Stages 4, 8, 12)
When the difference (P
Gi – P
L) is positive and the SOE
t < 100% then the excess energy will be used to charge the battery. The amount of energy that will be used to charge the battery (E
BAT,ch) the time interval Δt, is:
Where n
bat is the energy efficiency, which is considered equal to 0.9 [
25]. The study of the algorithm’s operation was done using the software excel and with the use of the add-in solver the optimal point of system’s operation was identified. The search for the optimal point was done by locating the operating point of the system (with regard to P
Gi and E
BAT,N) so that the E
Back and E
waste approach zero. In other words, to ensure that the system operates autonomously (Ε
back = 0) with the least possible discard of energy (E
waste = 0). Also, finding the optimal operating point was done for each month separately. The battery SOE at the start of each month was derived from the last period of the previous month as shown in
Figure 8.
As shown in
Figure 8, the battery SOE ranges between 20 and 100% throughout the year. Although zero energy deficit (E
Back = 0) was achieved, this was not possible for excess energy. In other words, the system can operate autonomously but will not avoid energy rejection.
Figure 9 shows the results of the optimal operation of the algorithm for five consecutive days of June. Specifically,
Figure 9a shows the operating intervals of the generators in the three operating areas for these five consecutive days of June.
The blue color (
Figure 9a) refers to the stable operation of the first generator in the range 0 – 125kW. The red line refers to the periods of stable operation of the two generators in the range 125 – 200kW. Finally, the green line concerns the operation of all three generators in the range 200 – 350kW. In the second area (red curve), which is essentially the transitional state from the first to the third area, the generators operate for a short period of time ranging from 30 to 135 minutes.
Figure 9b shows the energy charging the battery (blue curve) when the operating power of the generators is greater than the load demand power and the discharge energy of the battery (red curve) when the operating power of the generators is less than the load demand power.
Table 3 shows the maximum constant power that the generator should have in order for the power supply system to operate optimally in the three different zones (PG1, PG2 and PG3), as discussed above.
In the first range, the generator always operates. The maximum power, according to the simulation, will not exceed 75kW for this range. This level of power was detected in July. In other words, when the generator operates in that month in the first region, it works at a constant power of 75 kW. The minimum level of operation of the generator concerns the month of December. When it should operate only in the first range, it should be set to work constantly at 18kW. As can be seen from
Table 3, the generator should operate in the second and third areas during the months of May to October. This was expected according to the consumption profile (
Figure 4) and the way the algorithm works (
Figure 7). In the second area, the generator will operate at a constant power of 150kW in September and October. While for the months of May to August the constant operating power of the generator in the second range will be between 130 – 149.7kW. Similarly, for the third region, the generator operates at a constant power in a range of 209.7 to 295 kW. Specifically, in July, when the generator needs to operate in the third region, it will constantly operate at 295kW.
In practice, it would be possible to meet GHU energy needs by using two generators with a maximum power of 75kW and one generator with a maximum power of 145-150kW. In the first rwgion, only one 75kW generator would operate, while in the second region the two 75kW generators would operate. One would operate at full power (75kW) while the second would operate at power required to cumulative satisfy the power shown in the second column of
Table 3. In the third zone, all three generators would work together. The two 75 kW generators shall operate at full capacity and the rest of the load requirement, as exhibited in the third column of
Table 3, shall be covered by the third generator. A simple connection of the three generators with the battery is shown in
Figure 9.
The market investigation led to the selection of 75 and 150 kW natural gas generators with the consumption characteristics presented in
Table 4.
From
Table 4,
Figure 10 is created showing the change in MG fuel consumption as a function of operating power.
The equations that describe the change in natural gas consumption (FC) as a function of the operating power of the generator (P) are:
Using the above equations, it was possible to calculate fuel consumption for all months of the year and for all operating areas, which is presented in
Table 5.
A percentage of the required amount of fuel, according to
Table 5, could be covered by biogas production, utilizing the produced biomass of the greenhouse [
13,
14]. The rest will be covered by the supply of natural gas from the network. The cost of supplying this quantity of natural gas depends on the supplier and the negotiation that will take place. Nevertheless, the cost of fuel should be considerably lower than the cost of supplying electricity from the existing provider. According to the company's data, the cost of electricity consumption of this GHU for the period of the study (March 2022 – February 2023) was approximately € 252K. In order the investment to be profitable, in terms of operating costs, the cost of generator’s fuel should be lower than (252,000 €/year /257,681 m
3/year =) ~1€/m
3. In fact, it should be well below this value in order to be able to achieve the depreciation of the equipment’s purchase and installation. At least, imports from natural gas suppliers in Greece are made at prices much lower than 1€/m
3 [
27].
As far as operating costs are concerned, another weak point (apart from securing fuel for generators) of the proposed system is the battery and specifically its frequency of replacement. This depends on the operating mode and charge-discharge cycles. As for the former, it becomes obvious (
Figure 8) that half the time it is not sufficiently charged. Adequate charge can only be ensured when the system shows excess energy. So, for six months the battery will operate on a partial state of charge. The question is how harmful is this mode to the battery? This question was addressed by researchers who studied the effect of operating the battery on a partial state of charge on total charge-discharge cycles [
28,
29,
30]. These studies concluded that this mode of operation could benefit the battery.
According to the simulation, the energy discharged from the battery in one year was 61.8MWh. In addition, the results of the optimization showed that optimal operation is achieved when the rated battery energy is 5MWh. Thus, it could be estimated that the battery charge/discharge cycles (assuming 100% Depth of Discharge, DoD) are 61.8/5 ≈ 12 in a year. In the literature, it is reported that the cycle life of the lithium-ion battery exceeds 1000 for 80% DOD [
31]. In other words, the operating conditions of the battery in this system will not create the necessary conditions for its replacement during the investment period (e.g. 15-20 years). Therefore, the significant financial burden of this investment will result from the purchase and installation of the MGs and the ESS.
In the next paper, an attempt will be made to estimate the levelized cost of energy and the net present value to answer the question of whether this investment is profitable.