Use of Hybrid Photovoltaic Systems with a Storage Battery for the Remote Objects of Railway Transport Infrastructure

: The use of a grid-tied photovoltaic system with a storage battery to increase the power of objects of railway transport infrastructure above the limit on consumption from the grid with the possibility of energy saving is considered. The methods of analysis of energy processes in photovoltaic systems with a storage battery are used. They are added via the processing of archival data of power generation of a photovoltaic battery and computer modeling results. A technique of system parameter calculation to increase the power according to the given load schedule of the object at constant and maximum possible degree of power increasing is developed. The values of the average monthly generation of a photovoltaic battery at the location point of the object based on archival data are used. The principle of the control of power, consumed from the grid, according to the given values of the added and total load is developed. Using the basic schedule of added load power in connection with the graph of photovoltaic battery generation allows reducing the installed power of the storage battery. The additional reduction in the installed power of the photovoltaic and storage batteries is possible at the corresponding choice of the degree of power load increasing. The joint formation of current schedules with reference to the added power value and state of charge of the battery according to the short-term forecast of the generation of a photovoltaic battery is proposed. The value of added power at certain intervals of time is set according to the graph of actual generation of the photovoltaic battery, which contributes to the maximum use of its energy. With the average monthly generation of a photovoltaic battery in the spring–autumn period, the discharge of the battery during the hours of the morning load peak is not used. This reduces the number of deep discharge cycles and extends the battery life. The description of energy processes in steady-state conditions for the daily cycle of system functioning is formalized. On this basis, a mathematical model is developed in MATLAB with an estimation of the costs of electricity consumed from the grid. When modeling, archival data are used for days when the generation of a photovoltaic battery over time intervals is close to average monthly values. This makes it possible to evaluate the effectiveness of system management under conditions close to real during the year.


Introduction
Photovoltaic systems (PVSs) are the most common in the "green" energy sector. A considerable share of PVSs falls on local objects (LO) for various purposes, where hybrid solar power plants with connection to the alternating current distribution grid (DG) are used. This corresponds to the current trend of localizing consumption at the generation site [1].
The issue of localization of the consumption of energy, generated by PV, can be solved when using a PVS to increase the load power of the LO above the limit on consumption from the grid. Such an application of PVS can be in demand during the development (expansion) of facilities with an increase in consumption, when the possibilities of increasing the capacity of the existing connection to the power grid are exhausted.
PVS application can be useful for railway infrastructure facilities remote from the transformer substation, including those with the seasonal nature of the load (consumption). The simplified structure of the power supply of LO with a PVS and storage battery (SB) is shown in Figure 1. The implementation of this option may be cheaper than laying a new power transmission line and replacing DG equipment. The advantage of this solution is the possibility of reducing the cost of paying for consumption from the DG during the term of operation of the PVS and the possibility of autonomous operation in the case of DG disruptions. The issue of localization of the consumption of energy, generated by PV, can be solved when using a PVS to increase the load power of the LO above the limit on consumption from the grid. Such an application of PVS can be in demand during the development (expansion) of facilities with an increase in consumption, when the possibilities of increasing the capacity of the existing connection to the power grid are exhausted.
PVS application can be useful for railway infrastructure facilities remote from the transformer substation, including those with the seasonal nature of the load (consumption). The simplified structure of the power supply of LO with a PVS and storage battery (SB) is shown in Figure 1. The implementation of this option may be cheaper than laying a new power transmission line and replacing DG equipment. The advantage of this solution is the possibility of reducing the cost of paying for consumption from the DG during the term of operation of the PVS and the possibility of autonomous operation in the case of DG disruptions. This also applies to objects of extensive railway transport infrastructure, including remote objects on sections of the grid that are not currently electrified. The main task when using PVS for such facilities is to provide for their own needs. At the same time, the possibility of increasing their power when using existing DG is limited by the consumption limit. The issue of energy saving remains relevant. First of all, this concerns the reduction in electricity consumption from DG. The use of storage batteries in the PVS with modern methods of energy management [2] allows rationally redistributing the energy in the system in time. This also achieves a reduction in consumption from DGs and an increase in the reliability of the power supply of LO. Improving the performance of systems with renewable sources of electricity is an urgent task, which contributes to the further development of energy with distributed sources of electricity.

Literature Review and Problem Statement
The demand for the use of hybrid PVS with SB is confirmed by the fact that the electrical market is widely represented by various solutions of hybrid inverters [3,4], which are designed for LO. They have all the equipment for connecting a photovoltaic battery (PV) and SB, as well as sufficiently powerful software. They are designed for self-consumption of LO while reducing consumption from DG, and they provide an uninterruptible power supply function.
With a power consumption of LO in excess of 10 kW, it is advisable to use threephase multifunctional inverters while maintaining a power factor close to 1 at the point of This also applies to objects of extensive railway transport infrastructure, including remote objects on sections of the grid that are not currently electrified. The main task when using PVS for such facilities is to provide for their own needs. At the same time, the possibility of increasing their power when using existing DG is limited by the consumption limit. The issue of energy saving remains relevant. First of all, this concerns the reduction in electricity consumption from DG. The use of storage batteries in the PVS with modern methods of energy management [2] allows rationally redistributing the energy in the system in time. This also achieves a reduction in consumption from DGs and an increase in the reliability of the power supply of LO. Improving the performance of systems with renewable sources of electricity is an urgent task, which contributes to the further development of energy with distributed sources of electricity.

Literature Review and Problem Statement
The demand for the use of hybrid PVS with SB is confirmed by the fact that the electrical market is widely represented by various solutions of hybrid inverters [3,4], which are designed for LO. They have all the equipment for connecting a photovoltaic battery (PV) and SB, as well as sufficiently powerful software. They are designed for self-consumption of LO while reducing consumption from DG, and they provide an uninterruptible power supply function.
With a power consumption of LO in excess of 10 kW, it is advisable to use three-phase multifunctional inverters while maintaining a power factor close to 1 at the point of common coupling (PCC) to the grid [5][6][7][8][9][10][11][12]. This allows to unload the grid from reactive power and ensures the symmetry of loading the phases of the grid at an unbalanced load [9][10][11][12]. A previous study [13] presented a three-phase PVS with an SB with four converters with a common link of direct current (DC). In the DC link, a proportional-integral (PI) voltage controller (VC) was used. The controller set the currents of the battery and supercapacitor. The "multiconverter" provided a given active power with the proper quality of electricity, PV operation with tracking the maximum power point (MPPT), and extended battery life with hybrid storage. In [10,11], the structure of the control system of a multifunctional inverter of a PVS with an SB with voltage stabilization in the DC link with three VCs was considered. The VC controls the currents of the SB and PV, as well as current in the PCC. The structure changes in accordance with the mode of operation, and only one of the VCs is always used. Ensuring the efficiency of the use of hybrid PVS with energy storage for LO is usually associated with a reduction in the cost of consuming electricity from the grid and an increase in the reliability of the power supply [14]. With a wide range of changes in PV generation during the year, cost reduction is achieved by overestimating the power of PV relative to the load power. This allows providing acceptable indicators in cloudy weather and winter.
PV generation changes significantly during the day and year. Therefore, the increase in the efficiency of PVS energy management is associated with the use of the forecast of PV generation [15][16][17][18][19]. The issues with obtaining an accurate forecast were considered in a number of studies, particularly [15,16]. Recently, web resources have been made available that provide a forecast with a discreteness of 0.5 h or less, including an individual one at the location [20,21]. For the application without generation to the grid, the main purpose of the forecast is load planning for the next day and the possibility of adjustment when it changes [19].
Certain opportunities to reduce the cost of paying for electricity consumed by LO from DG are provided by taking into account the tariffication of payment [19,[22][23][24][25]. Features of the implementation of PVS in the application of static and dynamic tariffs were considered in [22].
When using PV for the needs of LO (without generating electricity to the grid), it is achievable to reduce the cost of electricity consumption from DG by up to five times at one tariff rate (up to seven times at wo rates) in the summer [19,26]. In winter, the reduction is insignificant-around 1.2 times. A general estimation of cost reduction for the year is not given. In [27], along with meeting the needs of the LO, the use of the planned generation of electricity in the grid during peak hours was considered. This allowed significantly reducing electricity costs. Overall estimation of cost reduction during the year was not given. Regardless, there was still an underutilization of PV energy in the summer.
An effective tool for assessing the efficiency of the energy management of a PVS with an SB is mathematical modeling [26][27][28][29]. Modeling of a hybrid system with a supercapacitor for the PV generation period was considered in [29]. Modeling of energy processes in the daily cycle with an estimate of the cost of paying for electricity consumed from the DG was considered in [26,27]. The use of archival data on PV generation [30] allowed studying the operation of the system in different weather conditions with an estimation of the cost of electricity, consumed from the DG.
The possibilities of PVS use for power increase for infrastructure objects, remote from the transformer substation, over the limit of consumption have not been sufficiently studied. This is related to the determination of parameters and limiting capabilities in different seasons of the year, the features of the formation of the SB state of charge in the process of operation, and the realization of the principles of added power formation at maximum use of PV energy. In this case, it is possible to reduce the installed power of the PV and the battery. An important role in assessing the capabilities of the system is performed by mathematical modeling.
Thus, the purpose of the article is to develop principles for the use of PVS with SB to increase the power of the LO above the power limit for consumption from the grid with the maximum use of PV generation. The main objectives of the research are as follows: to study the possibilities of increasing the load power of the LO above the power limit for consumption from the grid during the year; -to justify the choice of PVS parameters with the formation of the graph of added power according to the accepted load schedule of the LO with a decrease in the installed power of the PV and battery; -to develop principles for the implementation of the management of PVS using the forecast of PV generation; -to perform an assessment of the system's capabilities in the daily mode for different seasons of the year using mathematical modeling.

Methodology of Research
A study of ways to improve the control mechanism of the PVS with an SB for increasing the power of the LO was carried out on the basis of analytical methods in electrical circuits. The results of processing statistical data on the PV generation for a given point of location of the object were also used. A proportional increase in power to the original load schedule was adopted. As original, the load schedule characteristic of objects with a predominance of day loads, with peak loads in the morning and evening and a decrease in the load at night, was considered. The basic schedule of power, added and provided by the PVS, was adopted in accordance with the PV generation schedule. The maximum value of the power increase factor in winter takes into account the possibility of ensuring the battery charge within the power limit for consumption from the grid. On this basis, the energy capacity of the battery was determined, followed by an assessment of the possibilities for increasing the load power, taking into account the average monthly PV generation. The choice of a fixed value of the degree of power increase was carried out while taking into account the use of PV energy and cost reduction. The control system of the PVS converter unit was implemented on the basis of a classic double-loop structure with voltage stabilization in the DC link. When forming the SOC(t) (state of charge) of the SB schedule, a limit of DOD ≤ 80% (depth of discharge) was introduced with one deep discharge per day in the spring-summer-autumn period. The technique to calculate the reference of added power for maximum use of PV energy was realized on the basis of an analysis of the average monthly PV generation by time intervals for the taken load schedule. Analysis of energy processes in the system "DG-PVS with SB-LO load" was carried out for the daily cycle without taking into account transient processes and higher harmonics in energy converters. Energy losses were accounted for through efficiency. The properties of the SB were considered in accordance with the characteristics of the manufacturer. PV generation was estimated on the basis of monthly average values for given time intervals during the day. Data of PV generation were obtained for the location point of the object when processing archival data for 5 years. The modeling of energy processes was performed using MATLAB software package using real archive graphs of PV generation. The days were chosen when PV generation by time intervals was close to the average monthly values. The model was completed on the basis of analytical expressions for steady-state operating modes, which correspond to generally accepted proven calculation methods. When the operating modes of the system changed, the corresponding calculated expressions were used at time intervals per day.

The Results of the Research on the Use of PVS with SB to Increase the Power of LO
The option of PVS with an SB with a grid multifunctional inverter VSI (Figure 2) was considered. The middle pin (n) of the VSI link DC was connected to the neutral connection point of the DG. This allowed ensuring the equalization of power consumption from the grid by phases under unbalanced load LO [9][10][11]. This also made it possible to control the active power P g in the PCC. The structure of the PVS included the following elements: a DC voltage converter PV (CPV) with a transistor key for measuring the current of the PV shortcircuit [10], and a DC voltage converter battery (CSB). The converter unit was controlled by a control unit (CU) connected to the programmable control unit (PCU). Communication with the web resource to obtain forecast data was provided by a Wi-Fi WFM module. originally proposed, the use of a load schedule was considered in accordance with the standard distribution of peak loads [19,26] for objects of the utility sector and the nondomestic sector with a single-shift mode of operation. The following distribution of load intervals was accepted: night tariff zone in the period May-August (t7 = 24.00, t1 = 7.00), day tariff (t1 = 7.00, t2 = 8.00), (t3 = 11.00, t5 = 20.00), (t6 = 23.00, t7 = 24.00), peak load zones (t2 = 8.00, t3 = 11.00), and (t5 = 20.00, t6 = 23.00); in the period autumn-winter-spring (t1 = 7.00, t2 = 8.00, t3 = 10.00, t5 = 17.00 and 18.00, t6 = 17.00 and 18.00, t7 = 24.00). An additional point of time t4, was also used, corresponding to the transition to an evening decrease in PV generation. This time during the year varied from t4 = 16.00 in June to t4 = 14.00 in December. Increasing the power of the LO load in the daytime assumes that the total power PLC of the LO load is defined as PLC = PLg + PC (where PLg is the power which is provided by consumption from the grid (PLg does not exceed the limit on consumption PLIM), and PC is the added power, generated by the inverter due to the energy of PV and SB). We can take PLg = PLIM; then, with an increase in PL; the possible value of PLC grows with constant consumption from the grid. If the actual load power is less than PLC, the electricity consumption from the grid is reduced.
The value of the energy generated by PV during the year varies widely (Table 1). Table 1 shows the data [30] on the average monthly generation of PV per day WPVAVD at power PPVR = 1 kW for the Kyiv location: latitude (decimal degrees)-50.451, longitude (decimal degrees)-30.524. Similar data are given for the city of Žilina WPVAVDZ (latitude (decimal degrees)-49.224, longitude (decimal degrees)-18.748). In winter, the PV generation in Žilina is slightly higher. In Table 1 (in parentheses), monthly energy values per When solving the issue of increasing the power of the LO, the general approach changes somewhat; the DG becomes an auxiliary source of energy of limited power. As originally proposed, the use of a load schedule was considered in accordance with the standard distribution of peak loads [19,26] for objects of the utility sector and the nondomestic sector with a single-shift mode of operation. The following distribution of load intervals was accepted: night tariff zone in the period May-August (t 7 = 24.00, t 1 = 7.00), day tariff (t 1 = 7.00, t 2 = 8.00), (t 3 = 11.00, t 5 = 20.00), (t 6 = 23.00, t 7 = 24.00), peak load zones (t 2 = 8.00, t 3 = 11.00), and (t 5 = 20.00, t 6 = 23.00); in the period autumn-winter-spring (t 1 = 7.00, t 2 = 8.00, t 3 = 10.00, t 5 = 17.00 and 18.00, t 6 = 17.00 and 18.00, t 7 = 24.00). An additional point of time t 4 , was also used, corresponding to the transition to an evening decrease in PV generation. This time during the year varied from t 4 = 16.00 in June to t 4 = 14.00 in December.
Increasing the power of the LO load in the daytime assumes that the total power P LC of the LO load is defined as P LC = P Lg + P C (where P Lg is the power which is provided by consumption from the grid (P Lg does not exceed the limit on consumption P LIM ), and P C is the added power, generated by the inverter due to the energy of PV and SB). We can take P Lg = P LIM ; then, with an increase in P L ; the possible value of P LC grows with constant consumption from the grid. If the actual load power is less than P LC , the electricity consumption from the grid is reduced.
The value of the energy generated by PV during the year varies widely (Table 1). Table 1 shows the data [30] on the average monthly generation of PV per day W PVAVD at power P PVR = 1 kW for the Kyiv location: latitude (decimal degrees)-50.451, longitude (decimal degrees)-30.524. Similar data are given for the city of Žilina W PVAVDZ (latitude (decimal degrees)-49.224, longitude (decimal degrees)-18.748). In winter, the PV generation in Žilina is slightly higher. In Table 1 (in parentheses), monthly energy values per day (W* PVAVD ) and at time intervals during the day (W PV23 , W PV34 , and W PV45 ) are also presented for Kyiv. The energy values without parentheses correspond to the selected days when the generation was close to the monthly average. These values were obtained from archival data of P PV generation in Kyiv [30] for the period 2012-2016. The value of the added load power P C was determined on the basis of the average monthly daily generation of W PVAVD ≈ 2500 W in the transitional seasons of the year (Table 1): October and March. In November-February, the load, provided by PV, decreased. The average value of power in the daytime (P AVD = W PVAVD /t D , where t D is the length of the day) was about 200 W.
We took the basic load schedule (average P L values by time intervals) taking into account peak loads in the morning and evening with a decrease in load after t 4 until the evening peak, e.g., P L23B = 200 W (P LAVD ), P L34B = 180 W (0.9 P LAVD ), P L45B = 160 W (0.8 P LAVD ), and P L56B = 200 W (P LAVD ). The total night load of LO could be taken from the condition P Lg62 = P LIM − P B (we took P LIM equal to the peak power P LIM = P L23B = 200 W, P B = U B I B -power, consumed from the grid to charge the SB (U B and I B -voltage and current of the SB)). At the same time, P Lg62 ≥ P LMIN , P LgMIN = 0.2P L56 in summer, and P LgMIN = 0.3P L56 in winter. The total energy transmitted by the inverter to increasing power (P L ) at the interval (t 2 , t 6 ) was W L26 = 2740 Wh in summer, W L26 = 2580 Wh in autumn-spring, and W L26 = 2410 Wh in winter.
The effective use of the SB's capabilities for the redistribution of energy in the system involves the formation of the SOC(t) dependency Q * (t) ( Q * = 100Q/Q R , Q = Q 0 + I B dt , Q R = C B -the rated value (100%) or capacity (Ah) of the battery, Q 0 -the initial value). Increasing the power during peak hours in the morning and evening, when the PV generation is small, implies a deep discharge of the SB. That is, we have two deep discharges per day. In these conditions, the use of lithium-ion batteries is preferable.
Two variants of implementation were considered: (1) with a maximum increase in power in accordance with PV generation, which is available for the consumers with seasonal load; (2) with a constant increase in power during the year.
For variant (1), it was assumed that the planning of load using the day-ahead forecast was possible. For calculation of the value of SB energy capacity, in the interval (t 4 , t 6 ), the PV generation W PV45 is small, and the increase in power is achieved mainly due to the energy of the SB. At the same time, the energy balance is determined by the following expression: where W B is the SB energy capacity (W B = U B C B ), ∆Q * 46 = Q * 4 − Q * 6 , η C is the general efficiency of the SB voltage converter and grid inverter, and η B is the efficiency of the SB.
W B is calculated from the condition of the functioning of the added power of the load during the evening peak hours (t 5 , t 6 ) and on the intervals (t 4 , t 5 ), when PV generation is significantly reduced. This is most typical in winter, with a longer evening peak (4 h). The limitation is to ensure the possibility of the battery charge and the operation of the load at night within the framework of P LIM . When W PV46 = 0, the value of the energy capacity of the battery SB is where W C46 = W L46 (ρ − 1); ρ > 1 indicates the degree of increase in the load power. It was assumed that the night load in winter also proportionally increases P Lg = ρ 0.3 P LIM . The possible value for the battery charge in the framework of limit is denoted as Alternatively, The maximum value ρ MAX in the interval (t 4 , t 6 ) can be determined in accordance with Equations (2) and (3). In this case, at DOD 6 ≤ 80%, ρ MAX = 1.721. Accepting ρ = 1.7, W B = 1164 Wh. At DOD 6 ≤ 90%, the value is W B = 1034 Wh. The average value W B = 1099 Wh (corresponding to, for example, C B = 43 Ah at U B = 25.6 V) can be accepted. The resulting value is sufficient to ensure the added power at ρ MAX = 1.7 and average monthly PV generation in December in the interval (t 2 , t 6 ).
However, for the same value P LC = ρP L and P Lg ≤ P LIM , the different variants of reference of the added power P C for the interval (t 2 , t 6 ) (in Table 2, data are presented for ρ = 1.7) are possible. This allows planning P C (t) according to the conditions. Above the variant, va was considered. At minimal PV generation (in winter), variant vb ensured the biggest value ρ, whereby P C26 (t) = ρP L26 (t) − P LIM . In this case, W C46 = ρW L46 − P LIM (t 6 − t 4 ), and the value ∆W B62 can be expressed as We get values ρ MAX = 1.78 and W B = 968 Wh (corresponding to, for example, C B = 37.8 Ah at U B = 25.6 V). The advantage of this variant is the minimal value of P C in the interval (t 4 , t 6 ) when the PV generation is minimal.
Variant vc is included in Table 2. This variant is more tied to the PV generation schedule during the day.
There are days in winter when W PVD ≤ 60 Wh. In this case, it is possible to use a night charge of the battery up to 100%, followed by a discharge during the day (from 8:00 a.m. to 9:00 p.m.) at DOD 6 = 10%. This allows using the energy ∆W B26 = 768 Wh (at W B = 968 Wh) in the load. Accordingly, ρ = (∆W B26 + W L26LIM )/W L26LIM = 1.3. At small values of W PVD (for example, at W PVD ≤ 300 Wh), there is the possibility to redistribute the energy of SB by intervals. For example, the degree of increase in power in peak hours can be saved without increasing from 10:00 a.m. to 5:00 p.m. Then, for example, 0.3 ∆W B26 can be used in the morning peak, and, taking into account the duration, 0.6∆W B26 can be used in the evening peak. This allows ensuring ρ ≈ 1.6 in intervals (t 2 , t 3 ) and (t 5 , t 6 ).
In the period spring-autumn, there is the possibility to increase ρ. Determination of ρ is carried out in the accordance with forecast data of PV generation on the next day on the basis of the balance of energy in the intervals (t 2 , t 6 ) and (t 4 , t 6 ).
It is possible to calculate ρ for the interval (t 2 , t 6 ) as follows: where W gR26 indicates a decrease in energy consumption from the grid, and W LI M26 = P LI M (t 6 − t 2 ). Depending on the PV generation, the implementation is possible (a) with the discharge of the SB without a decrease in consumption from the grid, (b) without the discharge of the SB and a decrease in consumption from the grid, (c) with minimal discharge of the SB and a decrease in consumption from the grid, or (d) without the discharge of the SB and a decrease in consumption from the grid. Variant (a) is specific to winter with small PV generation, when ρ MAX = 1.7 is accepted.
For interval (t 4 , t 6 ), where k is a coefficient, taking into account the use of PV energy for consumption, and ∆W B46 corresponds to DOD = 80%. The value of k takes into account that, when the SB is fully charged by the time t 4 , only part of the PV energy is used for consumption by the load. The rest of the energy provides a reduction in consumption from the grid.
Furthermore, ρ 46 can be defined at the decrease in consumption in Equation (6) when W gR46 = 0. For March (Table 1), ρ 46 = 1.8. Compared with the value ρ 26 = 2.02 according to Equation (6) for option (b), we accept the smaller value. In this case, with ρ = 1.8, it is possible to implement option (c). Thus, for ρ = 1.8, we find the following value: For March ∆Q* 23 = 17% > 0, i.e., the SB state of charge increases, and the night battery charge is not needed. Then, in accordance with Equation (6), we have a decrease in consumption on the order of 572 Wh. With the average monthly values of PV generation for the summer period, ρ = 1.95 is achievable.
A graph of the added power P C (t) and the state of charge of the SB Q*(t) is presented in Figure 3 in accordance with the obtained value of ρ. We consider options vb and vc when using the PV energy W PV forecast data over time intervals. At the same time, it is desirable to (1) use the lowest possible power consumption at night (until 8:00 a.m.) per battery charge, and (2) ensure the condition Q* 4 → 100%, taking into account the reduction in PV generation in the evening. A prerequisite is the maximum use of PV energy.
(curve 7 in Figure 3), be unchanged (≈100%), or decrease. The possibility of an increase is excluded by increasing the PC45, which is carried out automatically. The formation of the degree of discharge in the interval (t5, t6) is carried out in the mode of regulation of the SB current.
The battery charge when Q* ≥ Q*d = 90-92% is reached at a constant value of the voltage [31]. In this case, the battery current is determined by the charge curve and is significantly reduced. Thus, the ability of the SB to receive energy is limited. Therefore, when Q* ≥ Q*d, the value of PC is set according to the actual PV generation as PLC ≥ PC = (PPV · ηC − PB) ≥ P 1 C. The limitation PLC ≥ PC is implemented by reducing the PV generation. The introduction of this mode allows ensuring more complete use of the PV energy, particularly when the calculated value of PC34 is less than required.  Figure 3 also shows a graph for the case when the PV energy is not enough to charge the battery up to 100% (curve 5) and for the case when the PV generation is close to 0 (curve 6). In these cases, the task is realized to support the load during peak hours.
Consider the variant with the constant increase of power during the year. For the choice of value ρ, consider the option vb with reference to the added power when the value WB is minimal. According to Equations (2) and (5), when ρ increases, WB also increases. An important issue is the underutilization of the installed PV power at high solar activity in the summer. To estimate, we introduce the PV energy utilization factor, Reference to the initial value Q* 2R is carried out in accordance with values ∆Q* 23 and ∆Q* 24 (calculated similarly to Equation (8)). If ∆Q* 24 ≤ 0, then Q* 2R = 100% (curves 1, 5, and 6 in Figure 3). In this case, when ∆Q* 24 > 0 and ∆Q* 23 ≤ 0, then Q* 2R = (100 − ∆Q* 24 ) ≥ 40% (curve 2 in Figure 3). If ∆Q* 24 > 0, ∆Q* 23 > 0, and W PV23 · η C /W 1 C23 < 1.5 (W 1 C23 is the value for the basic schedule with given ρ), then Q* 2R = (100 − ∆Q* 24 ) ≥ (Q* 6 + δ) (curve 3 in Figure 3, δ = 10-15%). In all cases, the reference of the added power is P C23 = P 1 C23 . If ∆Q* 24 > 0, ∆Q* 23 > 0, and W PV23 · η C /W 1 C23 ≥ 1.5, then P LC23 ≥ (P C23 = P PV23 · η C ) ≥ P 1 C23 , which is given in accordance with PV generation. In this case, the SB charge (curve 4 in Figure 3) is not carried out. It also does not use the battery charge at night, but there may be some battery charge before 8:00 a.m. (sunny morning).
The P C34 value in the interval (t 3 , t 4 ) is ≤ P LC34 , where ∆Q* 34 is defined using values Q* 2R and ∆Q* 23 , as described above.
In the interval (t 4 , t 5 ) with a large PV generation, the state of SB charge can increase (curve 7 in Figure 3), be unchanged (≈100%), or decrease. The possibility of an increase is excluded by increasing the P C45 , which is carried out automatically. The formation of the degree of discharge in the interval (t 5 , t 6 ) is carried out in the mode of regulation of the SB current.
The battery charge when Q* ≥ Q* d = 90-92% is reached at a constant value of the voltage [31]. In this case, the battery current is determined by the charge curve and is significantly reduced. Thus, the ability of the SB to receive energy is limited. Therefore, when Q* ≥ Q* d , the value of P C is set according to the actual PV generation as P LC ≥ P C = (P PV · η C − P B ) ≥ P 1 C . The limitation P LC ≥ P C is implemented by reducing the PV generation. The introduction of this mode allows ensuring more complete use of the PV energy, particularly when the calculated value of P C34 is less than required. Figure 3 also shows a graph for the case when the PV energy is not enough to charge the battery up to 100% (curve 5) and for the case when the PV generation is close to 0 (curve 6). In these cases, the task is realized to support the load during peak hours.
Consider the variant with the constant increase of power during the year. For the choice of value ρ, consider the option vb with reference to the added power when the value W B is minimal. According to Equations (2) and (5), when ρ increases, W B also increases. An important issue is the underutilization of the installed PV power at high solar activity in the summer. To estimate, we introduce the PV energy utilization factor, where W gR26 indicates a decrease in energy consumption from the grid, W C26 is the energy consumed by the added load, and m P is a coefficient of PV power recalculation relatively installed power P PVR = 1 kW. The value of W gR26 can take place during hours of high daytime solar activity t da , when P PV· η C > P C + P B (P B = U B I B is the power of SB charging). To exclude the generation of energy into the grid, the restriction P gR ≤ P LIM is used, which is achieved by regulating (reducing) the PV generation. In the limited case, W gR26 = P LIM · t da (t da = 6 h, from 10:00 a.m. to 4:00 p.m.).
The value of ∆W B26 in the summer period with the average monthly PV generation is taken as equal to ∆W B26 = 0 (nighttime battery charging is not required).
To assess the efficiency, we use the coefficient of cost reduction for electricity consumed from the grid (at one tariff rate and full use of PV energy). For the winter period (December), when the power increase in the interval (t 2 , t 6 ) is achieved while maintaining the consumption from the grid within the limit, where W LC = W LC62 + W C26 + P LI M (t 6 − t 2 ) is the total energy consumed by the load. At m P = 1, with an increase in ρ and, accordingly, W C26 , the value of k E decreases, and the value of k P increases. If the installed power (m P < 1) is reduced, then there is a reverse change in the coefficients.
Consider the option with ρ = 1.6 in comparison with the variant for ρ = 1.7 discussed above. This allows reducing the W B value from 968 Wh to 800 Wh (by 21%). Almost the same value of k E for December, in this case, can be obtained with m P = 0.86 and an increase in k P to 0.71 instead of 0.679 (in June). If we recalculate to the same load power, we get a decrease in the installed PV power by 9.4% and in the energy capacity W B of the SB by 14%. At the same time, it remains possible to increase the power to ρ = 1.8.
The reference of the value of added power and the formation of a graph of the state of charge of the SB are carried out according to the method discussed above.
A simplified structure of the PVS control system is shown in Figure 4. The system of automatic control of the CU converter unit is implemented according to well-known principles with voltage stabilization in the DC link Ud at the VSI input [10,13]. This is provided by three proportional-integral (PI) voltage controllers (VCs): VCI PV forms the PV current reference; VCI B forms the SB current reference; VCIg forms the reference of the current at the point of common coupling to the grid. The system (Figure 4) contains three channels: • PV generation control. This channel contains a control unit CPV (CSCPV), which provides the processing of the reference value of the current I 1 PV , current limiting unit LU1 with adjustable limit, and switch S2 for current settings in mode MPPT (position 2) or from VCI PV when regulating generation (position 1). Current limiting at the level I PVM = (0.9 ÷ 0.92)I SQ excludes PV operation in the short-circuit mode [10]; • Charge control of SB. This channel contains a control unit CSB (CSCSB), which provides the processing of the reference value of the current I 1 B , current limiting unit LU3 charge and discharge of SB, and switch S3 for current settings from VCI B (position 2) or PCU (position 1); • Grid current I g control (reference of power consumed from the grid). This channel contains a reference current unit and an inverter current control loop i C (RCU + CCL) and input of reference of the amplitude of the grid current (I 1 gm ), which, via switch S1, connects to VCIg (position 2) or PCU (position 1). Limitation unit LU2 has a lower I 1 gm ≥ 0, top I 1 gmLIM limits. Unit (RCU + CCL) provides i g in PCC taking into account the phase currents of the inverter i Ca,b,c, and load i La,b,c [10,11]. Energies 2022, 15, x FOR PEER REVIEW 11 of 20 The control of switches in the structure and the formation of current references are carried out by the PCU. The PCU also processes the prediction data according to the WFM block signal and the specified time intervals (TH). The phase lock loop (PLL) and offline control PVS channel when the mains voltage is turned off in Figure 4 are not shown (current reference of inverter phases i 1 Ca,b,c; in this case, a separate load voltage controller is set [10]).
The control principle (Table 3) is based on the control of active power, consumed from the grid (in PCC) Pg = PLC − P 1 C (P 1 C is the value of the added load power at the current time interval, and PLCi is the value of the active load power according to the measured values of currents and phase voltages) at PLIM ≥ Pg ≥ 0. The initial parameters are the recommended schedule (maximum average value) PLCR(t) and the base schedule PCB(t) for the accepted value of ρ, and the calculated schedule P 1 C(t). There is a possible situation when PLC ≠ PCR. With this in mind, deviation compensation ΔPLC = (PLC − PLCR) ≥ 0 is used in the reference P 1 C = P 1 C + ΔPLC. The operating modes of the control system by intervals and the used controllers are given in Table 3. Calculation expressions for steady modes are given in the description of the model of energy processes.
Consider the functioning of the system, starting from the intervals (t2, t3). PV operates in maximum power mode (MPPT). The SB current is set by the controller VCIB → I 1 B. Then, I 1 gm = √2Pgi/3Ugph (Ugph-phase voltage of the grid) is defined according to Pg = PLC − P 1 C23 and is supplied to the input of the current set unit RCU + CCL. When the load is reduced, Pg decreases; when the load is increased, Pg also increases. If WPV23·ηC/WC23 ≥ 1.5, then the value of PC is set by PV generation under the condition PCL ≥ (PC = PPV · ηC) ≥ P 1 C23. In this case, at PPV·ηC < P 1 C23, there is PC = P 1 C23. At PPV · ηC ≥ P 1 C23, there is PPVηC = PC and the current of SB charge decreases to 0; when PPVηC ≥ PLC, the SB charge is restored.  The control of switches in the structure and the formation of current references are carried out by the PCU. The PCU also processes the prediction data according to the WFM block signal and the specified time intervals (TH). The phase lock loop (PLL) and offline control PVS channel when the mains voltage is turned off in Figure 4 are not shown (current reference of inverter phases i 1 Ca,b,c ; in this case, a separate load voltage controller is set [10]). The control principle (Table 3) is based on the control of active power, consumed from the grid (in PCC) P g = P LC − P 1 C (P 1 C is the value of the added load power at the current time interval, and P LCi is the value of the active load power according to the measured values of currents and phase voltages) at P LIM ≥ P g ≥ 0. The initial parameters are the recommended schedule (maximum average value) P LCR (t) and the base schedule P CB (t) for the accepted value of ρ, and the calculated schedule P 1 C (t). There is a possible situation when P LC = P CR . With this in mind, deviation compensation ∆P LC = (P LC − P LCR ) ≥ 0 is used in the reference P 1 C = P 1 C + ∆P LC . The operating modes of the control system by intervals and the used controllers are given in Table 3. Calculation expressions for steady modes are given in the description of the model of energy processes. Table 3. Operating modes of the PVS with SB.

Interval (t 2 , t 3 ) (t 3 , t 5 ) (t 5 , t 6 ) (t 6 , t 2 )
Mode g1 g2 g3 g4 g5 g6 g7 Reference I 1 gm P g = P LC − P C , P gi → I 1 gm ≥ 0 P g = P LC − P C , P g → I 1 gm ≥ 0, Consider the functioning of the system, starting from the intervals (t 2 , t 3 ). PV operates in maximum power mode (MPPT). The SB current is set by the controller VCI B → I 1 B . Then, I 1 gm = √ 2P gi /3U gph (U gph -phase voltage of the grid) is defined according to P g = P LC − P 1

C23
and is supplied to the input of the current set unit RCU + CCL. When the load is reduced, P g decreases; when the load is increased, P g also increases. If W PV23 ·η C /W C23 ≥ 1.5, then the value of P C is set by PV generation under the condition P CL ≥ (P C = P PV · η C ) ≥ P 1 C23 . In this case, at P PV ·η C < P 1 C23 , there is P C = P 1 C23 . At P PV · η C ≥ P 1 C23 , there is P PV η C = P C and the current of SB charge decreases to 0; when P PV η C ≥ P LC , the SB charge is restored.
When switching to the interval (t 3 , t 5 ) and Q* < Q* d , the operating mode is saved. When Q* ≥ Q* d = 90-92%, the SB current is determined by the charging characteristic I B (Q*); if the set value I 1 B > I B (Q*), then VCI B goes into saturation. The battery cannot consume all the energy. This leads to an increase in the voltage U d on the capacitors in the DC link of the inverter. If P LC > P C ≥ P CB , then, upon reaching the switching threshold (U d + ∆U), the VCI g controller (g4 mode) is activated and reduces I 1 gm (P g ). If P PV η C ≥ P LC , then the VCI PV controller (g3 mode) is activated and I 1 PV (P PV ) is reduced. Thus, we have two conditions for switching of controllers (excluding the influence of transients).
In the interval (t 5 , t 6 ), the discharge of SB is carried out with current (given by controller VCI B ).
The added power is set in accordance with an average value of SB voltage U BAV as For the night period in the interval (t 6 , t 2 ), the reference of the current of SB charge is where Q* t is the current measured value, for example, with a resolution of 1 h. In the period spring-autumn, the SB charge in the morning is possible from PV at P PV· η C ≥ 0.
Referencing of the value Q* 2R and planning of the recommended (maximum) load P LCR is carried out according to the forecast for the next day. The given values of the added power in the intervals are specified at the beginning of the day in accordance with the current forecast of PV generation. Subsequently, with an interval of 1 h, the adjustment is carried out.

Modeling of Energy Processes in the Daily Cycle
The initial data for modeling were from an archive of PV generation on the mode of maximum power P PVM (t). Accordingly, the values of ρ, Q* 2R , the base P CB (t) schedule for the accepted value of ρ, and the calculated P 1 C (t) schedule were calculated. The load power values of the recommended P LCR (t) and the actual values of P LC (t), P CB (t), P 1 C (t) are given in tabular form. The time intervals are given by the variables t 12 , t 23 , t 34 , t 45 , t 56 , t 67 , and t 71 , which take on the value 1 at the corresponding time. The following auxiliary variables are also used: The variable w is precalculated. To measure the values of Q* t , Q* 5 , sample-and-hold schemes are used.
The power consumed from the grid is P g = (P LC − P C )t 26 + (P LC + P B )t 62 .
The current of the SB is The SB model is constructed according to the principles set in [26,27]. Data sheets given by the manufacturer were used [31]: charge characteristics I BC (Q * ) and U BC (Q * ) at I B ≥ 0 and discharge characteristic U BR (Q * ) at I B < 0, set in tabular form. The current can be calculated as The SB state of charge (SOC) taking into account energy losses is To estimate the reduction in electricity costs at a single tariff rate (taking its value equal to 1), the k E = W L /W g coefficient was used (W L = 24 0 P LC dt is the energy consumed by the LO load per day (without taking into account the energy on the SB charge at night), and W g = 24 0 P g dt is the energy consumed by LO from the grid).

Simulation Results
To set the PV generation, archival data were used [30] with the selection of days when the generation by intervals was close to the average monthly generation (Table 1). These days correspond to the energy values without parentheses in Table 1. Values Q * 6 , k E , and ρ for the case when P LC ≤ 2P LIM and the actual value of total load power P LC = P LCR are given in Table 4. Oscillograms P LC , P LCR, Pc, P PVM , P PV , Q * , I B , and P g (for clarity, P g is shown as negative) are given as described below.

•
In Figure 5 for the December day with a total generation twice below the average (W PV = 500 Wh). In this case, k E = 1.04, p = 1.5 (P PVR = 0.86 kW, W B = 800 Wh);   • In Figure 6 for the July day with a total generation in 3.3 times below the average (W PV = 1320 Wh). In this case, k E = 1.22, p = 1.6 (P PVR = 0.86 kW, W B = 800 Wh);    • In Figure 7b for the May day with the generation, corresponding to the average monthly values at P LC = P LCR , P PVR = 0.86 kW, W B = 800 Wh, and at limit values p = 1.8 with k E = 2.356; • In Figure 7c for the May day with the generation, corresponding to the average monthly values at P LC =P LCR , P PVR = 0.86 kW, W B = 800 Wh, and p = 1.6 with k E = 2.795; • In Figure 8 for the September day with the generation, corresponding to the average monthly values at P LC = P LCR , P PVR = 0.86 kW, W B = 800 Wh, and ρ = 1.6 with k E = 2.
(c)  Changing the reference of the added power PC23 in the interval (t2, t3) from a fixed calculated value to the value corresponding to the actual PV schedule in the spring-autumn period allows increasing kE by 2-4% with the exclusion of SB discharge. Figure 8. Oscillograms P LC , P LCR , P C , P PVM , P PV , Q * , I B , and P g for the September day with the generation, corresponding to the average monthly values at P LC = P LCR , P PVR = 0.86 kW, W B = 800 Wh, and ρ = 1.6 with k E = 2 (Q * and I B values are shown at scale 2 and 10, respectively).

Discussion of the Results of the Study on Increasing the Power of the LO Using PVS with SB
Changing the reference of the added power P C23 in the interval (t 2 , t 3 ) from a fixed calculated value to the value corresponding to the actual PV schedule in the spring-autumn period allows increasing k E by 2-4% with the exclusion of SB discharge.

Discussion of the Results of the Study on Increasing the Power of the LO Using PVS with SB
Increasing the load power of the LO above the limit of power consumption from the grid while reducing the installed power of PV and SB and decreasing the cost of paying for electricity, consumed by the LO, from the DG is possible due to the following: -Use of the basic schedule of added power, tied to the PV generation. This reduces the energy required for its implementation. This allows increasing the degree of power increase while reducing the energy capacity of the SB; -Limiting the degree of power increase at an intermediate value with a decrease in the installed power of the PV and SB. This provides an improvement in the use of PV energy without increasing the cost of electricity, consumed from the grid; -Referencing the current value of the added load power and the SOC value of the battery at time intervals, taking into account the predicted and actual PV generation. Due to the change in certain time intervals of the added power reference from the calculated value to a value that repeats the law of the change in the PV generation power, this contributes to more complete use of the PV energy; -Exclusion of the SB charge at night with the average monthly PV generation in the spring-summer-autumn period helps to reduce electricity consumption from the grid. This excludes battery discharge during the hours of the morning load peak, which helps to reduce the number of battery discharge cycles and increase its service life.
This article is a development of previous studies [19,26], which considered increasing the efficiency of hybrid PVS with SB for the needs of local facilities. This was achieved by reducing the cost of electricity consumed from the grid when using the PV generation forecast. A common problem is the significant overestimation of the PV power in relation to the load power of the LO, which is necessary for use in conditions of low PV generation.
As a result, even with medium PV generation, there is a significant underutilization of PV energy with the need for regulation of PV generation. A feature of the proposed solutions is a change in the general approach to the use of PV and SB energy when the load power, added to the limit, is formed. Directed formation of the graph of added power allows reducing the installed power of PV and SB.
There are certain limitations regarding the use of the results of the work, as described below.
-An object was considered with the main load in the daytime in the presence of peak loads in the morning and evening hours. At the same time, it is possible to charge the SB at night within the limit on consumption from the grid; - The possibilities of increasing power are seasonal in nature with a maximum value in the period spring-summer-autumn; - The optimal choice of values for the degree of increase in power (ρ) and the installed PV power (m P ) involves taking into account many factors, including the costs of acquisition and ongoing maintenance. Such a task was not set in this work, and the approach was simplified for evaluation; - The assessment of a possible reduction in the cost of paying for electricity during the year was somewhat simplified and was performed for one tariff rate. We considered the days when the PV generation corresponds to the average monthly values for the accepted time intervals; - The implementation of the control system assumes an "open" structure of the relevant channels of control; -When modeling, it was assumed that the graph of the power generated by the PV corresponds to the forecast and does not change during the day.
Further development of this work is connected with the optimization of system parameters. It is also required to study the possibilities of correction of deviations in the values of the actual PV generation according to the forecast data and possible changes in the current forecast during the day.

Conclusions
With the selected parameters of the PVS, it is possible to increase the power of the LO over the limit for consumption from the grid up to 1.7-1.95 times. This depends on the average monthly PV generation of during the year. Limiting the degree of power increase to a value of 1.6, when choosing the parameters, allows reducing the installed power of the PV and SB. The possibility of increasing the degree of power increase to a value of 1.8 (if necessary) remains in this case.
It is advisable to select the parameters of the PVS on the basis of the data on the average monthly PV generation for the taken schedule of the load. It is possible to use archival data from web resources with open access at any point location of the object. The ratio of PV power and added load power is determined on the basis of PV generation in the transitional seasons (in this case, March and October). The energy capacity of the SB is determined by the power of the added load from the condition of the sufficiency of its operation in the pre-evening hours and evening peak hours with a given limitation of the DOD battery. The use of the basic schedule of the added power with reference to solar activity, while maintaining the resulting power of the LO, makes it possible to reduce the energy capacity of the SB. Therefore, when setting the degree of power increase ρ = 1.7 for winter, the energy capacity can be reduced by 13.5%. Limiting the value of ρ at the level of ρ = 1.6 allows reducing the installed power of the PV and the SB. When converted to the same total load power, the reduction for P PVR is 9.4%, and that for W BR is 14%.