Optimal Sizing of Battery-Integrated Hybrid Renewable Energy Sources with Ramp Rate Limitations on a Grid Using ALA-QPSO

: Higher penetration of variable renewable energy sources into the grid brings down the plant load factor of thermal power plants. However, during sudden changes in load, the thermal power plants support the grid, though at higher ramping rates and with inefﬁcient operation. Hence, further renewable additions must be backed by battery energy storage systems to limit the ramping rate of a thermal power plant and to avoid deploying diesel generators. In this paper, battery-integrated renewable energy systems that include ﬂoating solar, bifacial rooftop, and wind energy systems are evaluated for a designated smart city in India to reduce ramping support by a thermal power plant. Two variants of adaptive-local-attractor-based quantum-behaved particle swarm optimization (ALA-QPSO) are applied for optimal sizing of battery-integrated and hybrid renewable energy sources to minimize the levelized cost of energy (LCoE), battery life cycle loss (LCL), and loss of power supply probability (LPSP). The obtained results are then compared with four variants of differential evolution. The results show that out of 427 MW of the energy potential, an optimal set of hybrid renewable energy sources containing 274 MW of rooftop PV, 99 MW of ﬂoating PV, and 60 MW of wind energy systems supported by 131 MWh of batteries results in an LPSP of 0.005%, an LCoE of 0.077 USD/kW, and an LCL of 0.0087. A sensitivity analysis of the results obtained through ALA-QPSO is performed to assess the impact of damage to batteries and unplanned load appreciation, and it is found that the optimal set results in more energy sustainability.


Introduction
Energy being one of the most imported commodities in developing nations, an increase in crude oil and coal prices witnessed widening fiscal deficits. In contrast, renewable energy costs have drastically reduced from 2010 to 2019 owing to mass production. These favorable conditions allowed developing economies to invest in eco-friendly and costeffective renewable energy sources to comply with Sustainable Development Goals decided during the Paris Agreement. India, with its crude imports pegged at 82% of its energy needs, witnessed an 84% fall in renewable energy costs in this period, with tenders falling Goals decided during the Paris Agreement. India, with its crude imports pegged at 82% of its energy needs, witnessed an 84% fall in renewable energy costs in this period, with tenders falling to as low as 33 USD/MWh [1]. With digital transformation creating exciting employment opportunities, urban areas in India are witnessing a mass exodus from rural areas [2]. This has resulted in significant demographic changes that have propped up the energy, water, and transportation needs of the cities. This has led to dented air quality mainly due to conventional transportation and outdated and inefficient industrial equipment. With all these developments, the energy needs of urban areas are set to double by 2030. Substantial deployment of renewable energy systems is essential to evade coal-/gas-based power plants to meet the additional energy needs.

Current Electric Energy Scenario in India
With the introduction of the Electricity Act of 2003, plenty of opportunities unfolded for private players as all the electric power segments (generation, transmission, and distribution) opened up for investments. As of February 2020, the total installed capacity in India was 367.28 GW, of which only 25.2% was from the central sector and the remaining was shared between state and private sectors. The capacity addition is at a compounded annual growth rate of 5.19%, which reflects staggering economic growth in recent years. Besides, the share of renewable energy is pegged at 23.5%, though it is far from the ambitious target of 175 GW of renewable additions by 2021-2022. The rising cost of oil and coal mining policy issues in coal-exporting countries have led to a decrease in the plant load factor of conventional generators, as shown in Figure 1 [3]. This has also led to reduced growth of conventional energy sources, as evident from Table 1. Difficulties in land acquisition, especially in urban areas, where the actual demand exists, have led to diminished interest from private players. In this context, municipal corporations and distribution companies are required to step up investments in renewable energy systems in partnership with private players to exude confidence in private players. In addition, municipal corporations can benefit from such partnerships that open up investments in other needs of the city, such as piped gas and fiber optics. To make such partnerships viable, a proper methodology has to be followed for optimally sizing renewable energy systems. Since the energy usage pattern is subjected to economic activity, which is affected by factors such as special economic zones, access to roads, and the main city, site selection becomes a key factor in planning for renewable energy systems. With the data regarding economic activity in hand, land cover and load data can be forecasted, and hence, a distribution system can be planned accordingly.  In addition, municipal corporations can benefit from such partnerships that open up investments in other needs of the city, such as piped gas and fiber optics. To make such partnerships viable, a proper methodology has to be followed for optimally sizing renewable energy systems. Since the energy usage pattern is subjected to economic activity, which is affected by factors such as special economic zones, access to roads, and the main city, site selection becomes a key factor in planning for renewable energy systems. With the data regarding economic activity in hand, land cover and load data can be forecasted, and hence, a distribution system can be planned accordingly.

Electric Energy Scenario in the Proposed Location
As of February 2020, the installed capacity of Andhra Pradesh was 24.803 GW, with renewable energy pegged at just 8.313 GW in contrast to the target of 18 GW by 2021-2022 [4]. To meet the target, involvement from all the stakeholders such as municipal corporations, cooperative societies, and private players is desired. Visakhapatnam, being the second-largest city in Andhra Pradesh, is contributing more than 10% of the state gross domestic product, strengthened by its diversified industrial segments such as steel and pharma and needs to contribute heavily in meeting the renewable energy target. With the city designated as a smart city by the central government, funds amounting to 150 million USD yearly are awarded for its development. The Greater Visakhapatnam Municipal Corporation (GVMC), with modest revenue receipts, is leveraging itself in sustainable projects in the city.
Recent renewable investments of the GVMC, a 2 MW floating solar PV power plant in the Mudasarlova Reservoir, and a proposed 15 MW floating power plant over the Meghadrigedda Reservoir are evidence of its efforts for deploying clean energy systems in the city. Further, the GVMC has installed rooftop solar panels over all the government offices and schools, making use of the incentives provided by the central government. The GVMC has also been actively installing LED-based street lights all over the city and has installed 91,775 LED streetlights with a project cost of 10.5 million USD, resulting in an annual saving of 24 GWh of energy. This has been achieved through the annuitybased deemed savings model without any investment from the GVMC. These investment models by the corporation can further be leveraged to envision a DC distribution system to facilitate aggressive solar PV installations by the public as well.

Literature Review
High penetration of variable renewable energy sources results in wide frequency fluctuations due to the low inertia of the electric grid [5]. To improve the reliability of the electric grid, it is highly recommended to deploy an energy storage system (ESS). Of all ESS technologies, a battery energy storage system (BESS) is preferable owing to their modularity. Batteries are used for grid applications such as time shifting, regulation, spinning and non-spinning reserves, voltage support, demand charge management, load following, transmission/distribution line upgrade deferral, and black start. However, their capacity is limited by capital and O&M costs, and regular replacements though their costs are more competitive than ever [6]. The size of a BESS drastically changes the economics of the entire project. The optimal sizing of a BESS under different settings is widely cited in the literature [7][8][9][10]. Lithium-ion-based BESSs have been studied to act as spinning reserves when the grid capacity is maxed out, as they have advantages of a large number of cycles and specific energy compared to lead-acid batteries. However, since they suffer from degradation issues, estimation techniques are employed to assess the health of a BESS [11]. Minimizing life cycle loss is considered as one of the objectives to avoid overexploitation of a BESS. In this context, a balanced mixture of renewable energy sources and conventional energy sources is required. Nevertheless, conventional energy sources are required to operate at their technical minimum and are regulated by their ramping limits. This necessitates the optimal sizing of a hybrid renewable energy system (HRES) by using appropriate techno-economic reliability indices.
Meta-heuristic methods such as particle swarm optimization, differential evolution, evolutionary algorithms, and nature-inspired techniques such as ant-lion optimization and gray-wolf optimization are used in various research areas of electrical engineering [12][13][14][15][16]. Particle swarm optimization (PSO), introduced by Kennedy and Eberhart in 1995, is a nature-inspired algorithm inspired by the collective behavior of birds in search of food. The introduction of traditional PSO created a lot of interest in the research community, and it led to adaptations to it, resulting in several variants of PSO to improve the exploration and exploitation of particles. The convergence behavior of particles was analyzed in [17], which demonstrated that the particles in PSO are attracted by an attractor. By using this phenomenon, several variants of PSO were proposed. However, not all of them guarantee global convergence. Along similar lines, QPSO was introduced in 2004 inspired by both quantum mechanics and the concepts of PSO for enhanced exploration capabilities [18]. In QPSO, the velocity component in the original PSO is replaced by a delta potential well, thus enabling the particles to venture the entire search space with some probability. This allows the particles to attain a global optimum. There have been several variants of QPSO as well [19][20][21]. An adaptive local attractor QPSO (ALA-QPSO) was proposed in [21] to enhance the search performance by introducing the weighted mean personal best of the particles.
Optimal sizing of an HRES has been widely explored in the literature, however with different combinations of technologies [16,22]. Most of the research work is focused on stand-alone systems to enhance the sustainability of remote areas [23][24][25]. Recently, the focus has been shifted toward grid-connected scenarios due to the fall in solar PV prices [26][27][28][29][30]. There have been quite a few studies on optimal sizing of large-scale gridconnected systems, mostly due to the complexities involved in renewable energy potential assessment, especially for rooftop PV systems [31]. To increase power density, bifacial PV panels were studied in [32,33]. The use of roofing materials to improve rear-side power generation results in reduced cooling needs in heat-prevalent areas. The integration of floating solar PV systems is on the rise worldwide due to low leasing costs and ease of plant availability for maintenance [34][35][36]. Due to the intermittent nature of renewable energy sources, energy storage systems are included in the HRES to provide stiffness to the grid. There are several studies on optimizing the size of energy storage systems, focusing on peak shifting and renewable smoothing, along with ancillary services such as voltage regulation and frequency regulation [28,37,38]. In [39], a multi-objective equilibrium optimizer was used to achieve operational and economic efficiency by integrating renewable energy sources into the grid. In [40], an optimal combination of PV systems and DSTATCOM was achieved using the multi-objective modified ant lion optimizer. With the growing acceptance of electric vehicles, Li-ion batteries are becoming more prominent in vehicle-togrid scenarios [28,[41][42][43]. In addition to modularity, significant research in fault diagnosis, and second-use energy storage applications, a BESS has become more acceptable compared to other technologies [44][45][46].
To avoid the premature convergence problem of conventional PSO, a double-exponentialfunction-based dynamic inertia weight technique was developed for the optimal parameter estimation of a PV cell module [47]. As a way to solve energy management problems, lower bill costs, peak-to-average ratio, and carbon emissions, efficient integration of RES with BES a load scheduling and energy storage Management controller based is discussed [48].

Key Contributions
Despite a large number of articles in the area of optimal sizing of RE sources, there are few papers that address potential assessment in cities using GIS tools to arrive at the upper bounds. Most of the papers address the issue of either cost optimization or voltage and frequency issues, but they do not consider the limitations on the ramp rates of thermal This paper addresses the issue of optimally sizing the HRES in Visakhapatnam by considering the sustainability and ramp rate limitations on conventional energy sources using GIS tools to assess the renewable energy potential. The process flowchart for the proposed work is shown in Figure 2. Initially, a city is selected and its potential for renewable energy sources is studied using geographic information system (GIS) tools. Based on the available resources, possible land cover is assessed and optimal sizing of the system is arrived at by applying optimization techniques to realize techno-economic requirements such as loss of power supply probability (LPSP), levelized cost of energy (LCoE), and battery life cycle loss (LCL). and frequency issues, but they do not consider the limitations on the ramp rates of thermal power plants (TPPs) and the financial repercussions on the TPPs in the event of phasing out diesel and gas generators to improve the sustainability of cities. This paper addresses the issue of optimally sizing the HRES in Visakhapatnam by considering the sustainability and ramp rate limitations on conventional energy sources using GIS tools to assess the renewable energy potential. The process flowchart for the proposed work is shown in Figure 2. Initially, a city is selected and its potential for renewable energy sources is studied using geographic information system (GIS) tools. Based on the available resources, possible land cover is assessed and optimal sizing of the system is arrived at by applying optimization techniques to realize techno-economic requirements such as loss of power supply probability (LPSP), levelized cost of energy (LCoE), and battery life cycle loss (LCL).
The main contributions of this paper are as follows: i. A methodology for site selection, resource assessment, and energy management for large-scale renewable energy integration is developed. ii.
Optimal sizing of potential renewable energy sources and a battery bank is assessed to minimize grid ramping, levelized cost of energy, and loss of load using various multi-objective optimization techniques.  The main contributions of this paper are as follows: i.
A methodology for site selection, resource assessment, and energy management for large-scale renewable energy integration is developed. ii. Optimal sizing of potential renewable energy sources and a battery bank is assessed to minimize grid ramping, levelized cost of energy, and loss of load using various multi-objective optimization techniques.
Subsequent sections of this paper are organized as follows: Section 2 deals with site selection and renewable resource assessment in the chosen urban areas. In Section 3, a mathematical model of the presented renewable energy sources is developed. In Section 4, the adopted multi-objective optimization technique is explained, while in Section 5, the results are compared with legacy optimization techniques. In addition, sensitivity analysis is presented to evaluate the impact of reduced battery size due to outages and price appreciation due to policy changes. Finally, concluding remarks are made in Section 6.

Site Selection and Resource Assessment
The proposed HRES with a battery bank is shown in Figure 3. The system contains decentralized rooftop PV panels, floating PV panels, and wind turbines installed over hilltops and seashores. Due to their intermittency, suitable converters are placed to integrate with the battery-integrated DC bus. The DC bus accommodates EV charging stations and street-lighting systems with which better efficiency is envisaged. A battery bank is also associated with the DC bus for charging and discharging during the mismatch between generation and load demand. Further, the DC bus is integrated with the AC grid through a transformer to supply the AC loads of the city. After matching the supply and demand, any excess power is wheeled to other load centers via the grid.
results are compared with legacy optimization techniques. In addition, sensitivity analysis is presented to evaluate the impact of reduced battery size due to outages and price appreciation due to policy changes. Finally, concluding remarks are made in Section 6.

Site Selection and Resource Assessment
The proposed HRES with a battery bank is shown in Figure 3. The system contains decentralized rooftop PV panels, floating PV panels, and wind turbines installed over hilltops and seashores. Due to their intermittency, suitable converters are placed to integrate with the battery-integrated DC bus. The DC bus accommodates EV charging stations and street-lighting systems with which better efficiency is envisaged. A battery bank is also associated with the DC bus for charging and discharging during the mismatch between generation and load demand. Further, the DC bus is integrated with the AC grid through a transformer to supply the AC loads of the city. After matching the supply and demand, any excess power is wheeled to other load centers via the grid.
As mentioned in Section 1, Visakhapatnam has become an investment destination for companies ever since the bifurcation of the state of Andhra Pradesh, along with the designation of a smart city. The city is one of the leading pharmaceutical hubs in India that cater to world needs. In addition, electricity demand is expected to spike further due to announcements such as a railway zone and industrial corridors, leading to further investments. However, capacity addition using thermal power plants takes a long time. Further, some of the power purchase agreements have already expired in 2019, while some power plants in the vicinity are close to their end of life. In contrast, renewable energy systems take only months to install and integrate with the grid and their capacities can be easily expandable at a faster pace. This study focuses on the city boundaries, where there is a drastic change in land usage and building cover due to the special As mentioned in Section 1, Visakhapatnam has become an investment destination for companies ever since the bifurcation of the state of Andhra Pradesh, along with the designation of a smart city. The city is one of the leading pharmaceutical hubs in India that cater to world needs. In addition, electricity demand is expected to spike further due to announcements such as a railway zone and industrial corridors, leading to further investments. However, capacity addition using thermal power plants takes a long time.
Further, some of the power purchase agreements have already expired in 2019, while some power plants in the vicinity are close to their end of life. In contrast, renewable energy systems take only months to install and integrate with the grid and their capacities can be easily expandable at a faster pace. This study focuses on the city boundaries, where there is a drastic change in land usage and building cover due to the special economic zones. The type of housing is more inclined toward multi-story buildings where rooftops are under-used or individual houses with sufficient space for solar installations. The GVMC can propose a public-private partnership model with those rooftop owners and hence reduce the leasing costs of rooftop solar power systems. Besides, water bodies are explored to enable the interested parties to invest in floating solar PV systems. This will also benefit Energies 2021, 14, 5368 7 of 23 the city with reduced water evaporation as the city is struggling for water security for decades. Likewise, wind energy systems can be mounted on shorelines and hilltops. Based on the master plan provided by the infrastructure firm AECOM [49], growth centers and human settlements in the proposed city are shown in Figure 4. economic zones. The type of housing is more inclined toward multi-story buildings where rooftops are under-used or individual houses with sufficient space for solar installations. The GVMC can propose a public-private partnership model with those rooftop owners and hence reduce the leasing costs of rooftop solar power systems. Besides, water bodies are explored to enable the interested parties to invest in floating solar PV systems. This will also benefit the city with reduced water evaporation as the city is struggling for water security for decades. Likewise, wind energy systems can be mounted on shorelines and hilltops. Based on the master plan provided by the infrastructure firm AECOM [49], growth centers and human settlements in the proposed city are shown in Figure 4. It is observed that the industrial activities in Parawada and Atchutapuram areas have enabled satellite growth centers, marked by orange stripes, to grow at a faster rate, in addition to the expansion of the city core. Given this, the area marked with grid lines is taken for the proposed area for the HRES to be set up. Assessment of HRESs is done at the proposed location using the Polygon tool of Google Earth Pro, where purple and blue patches along with yellow pins indicate potential areas for floating PV, rooftop, and wind energy system installations, respectively, as shown in Figure 5. It is observed that the industrial activities in Parawada and Atchutapuram areas have enabled satellite growth centers, marked by orange stripes, to grow at a faster rate, in addition to the expansion of the city core. Given this, the area marked with grid lines is taken for the proposed area for the HRES to be set up. Assessment of HRESs is done at the proposed location using the Polygon tool of Google Earth Pro, where purple and blue patches along with yellow pins indicate potential areas for floating PV, rooftop, and wind energy system installations, respectively, as shown in Figure 5. Four potential areas for floating PV systems, are identified and their areas are calculated. Similarly, six possible areas for rooftop solar power systems are recognized after taking out unused land cover and dense areas where the installation is not viable. Further, the potential of wind energy systems is assessed over hilltops with a minimum of 200 m height in the region and on the shorelines that are not part of any economic activ-  Four potential areas for floating PV systems, are identified and their areas are calculated. Similarly, six possible areas for rooftop solar power systems are recognized after taking out unused land cover and dense areas where the installation is not viable. Further, the potential of wind energy systems is assessed over hilltops with a minimum of 200 m height in the region and on the shorelines that are not part of any economic activity. The potential of each HRES is given in Table 2, along with the potential quantity of solar panels and wind turbines. These values are considered as the upper bounds for optimal sizing of the HRES. Further, accessibility of the HRES to the grid is assessed, as shown in Figure 6 [50], through which concerns regarding grid connectivity are addressed. The proposed location has a good number of substations with 132 kV, 220 kV, and 400 kV voltage levels. The entire energy in the HRES can be evacuated through these substations.

Observations from the Load Profile of the Proposed Location
The load profile is generated from five different areas in the proposed location, as shown in Figure 7, that represent various load segments, viz. predominantly industrial (Gajuwaka and switching station), residential (Pendurty and dairy farm), and pharma (Parawada) sectors. It is also evident that pharma needs are high and the load profile is almost flat similar to residential loads. However, the heavy industry zone has a diversified load profile with peaks from 6 PM to 11 PM. From the load profile at the proposed location, it is observed that demand rises from the first quarter (Q1) onward and peaks in the second quarter, while Q4 exhibits weak demand. Since the peak loads correspond to the second quarter during which monsoon retreats, solar PV systems suffer from cloud cover. Hence, a BESS plays a crucial role during these days. Besides, the daily load curve indicates that

Observations from the Load Profile of the Proposed Location
The load profile is generated from five different areas in the proposed location, as shown in Figure 7, that represent various load segments, viz. predominantly industrial (Gajuwaka and switching station), residential (Pendurty and dairy farm), and pharma (Parawada) sectors. It is also evident that pharma needs are high and the load profile is almost flat similar to residential loads. However, the heavy industry zone has a diversified load profile with peaks from 6 PM to 11 PM. From the load profile at the proposed location, it is observed that demand rises from the first quarter (Q1) onward and peaks in the second quarter, while Q4 exhibits weak demand. Since the peak loads correspond to the second quarter during which monsoon retreats, solar PV systems suffer from cloud cover. Hence, a BESS plays a crucial role during these days. Besides, the daily load curve indicates that load peaks from 6 PM to 10 PM, during which time energy is drawn either from the grid or from the combination of a BESS and a wind energy system.

Land Cover and Shading Analysis
From the land cover of Visakhapatnam, as shown in Figure 8a, it is seen that building density, indicated in red color, is on the rise in the proposed location, with a large number of industries being set up in the special economic zones. Similarly, the blue color indicates water bodies with sufficient depth for floating solar panels. The 3D map of one

Land Cover and Shading Analysis
From the land cover of Visakhapatnam, as shown in Figure 8a, it is seen that building density, indicated in red color, is on the rise in the proposed location, with a large number of industries being set up in the special economic zones. Similarly, the blue color indicates water bodies with sufficient depth for floating solar panels. The 3D map of one of the residential areas is obtained using an open CADmapper online tool and is fed to the Sketchup Pro, which is also freeware, for shading analysis. The shading time is rated from yellow to dark blue, yellow indicating less shading time and dark blue indicating more shading time.

Mathematical Modeling
Due to the lack of sufficient data regarding bio-waste and its feasibility in the proposed location, along with technology developments in other renewable energy technologies, only solar and wind energy technologies are considered in this study. The model equations for the two technologies are well cited in the literature, and hence only their equations are provided in this section to avoid redundancy. However, modeling of bifacial PV panels that are used for rooftop systems is presented explicitly.

Solar PV System
In this work, bifacial solar panels are considered for rooftop PV systems, while mono-facial PV panels are considered for the floating PV system. The modeling of irradiance on the front part of the panels is the same for both cases, while only rear irradiance applies for bifacial panels. As a first step, the angle of incidence on both sides of the panel is calculated, as in Equation (1). Equations (1)-(6) of the irradiance of PV panels are adapted from [51].
There are a total of three components of irradiance that fall on the front and rear sides of the bifacial panel, i.e., direct normal irradiance (DNI), diffused horizontal irradiance (DHI), and an albedo, due to reflected irradiance, as shown in Figure 9. The total irradiance on the PV panels is given by Equation (2)  The result shows that a good number of rooftops can be used for installation as long as there is sufficient spacing between buildings when adjacent building heights differ by a large margin. A rooftop PV system is designed for a sample building in the proposed location, and shading analysis is performed at 4 PM when the maximum shading occurs, as shown in Figure 8b. The total area of the rooftop is 323 m 2 , including water tanks and a stairway. The illustrated PV array has an area of 183 m 2 (70 panels of 2.61 m 2 each) and hence amounting to a ground coverage ratio of 0.56. Though azimuth correction of 42 • can be applied for the rooftop system, thus increasing the irradiance on the rear part of the panels to boost the bifacial gain, it is aesthetically not appealing. In some cases, this causes inconvenience to operating personnel during maintenance. Hence, azimuth correction is advised for those buildings constructed with azimuth correction.

Mathematical Modeling
Due to the lack of sufficient data regarding bio-waste and its feasibility in the proposed location, along with technology developments in other renewable energy technologies, only solar and wind energy technologies are considered in this study. The model equations for the two technologies are well cited in the literature, and hence only their equations are provided in this section to avoid redundancy. However, modeling of bifacial PV panels that are used for rooftop systems is presented explicitly.

Solar PV System
In this work, bifacial solar panels are considered for rooftop PV systems, while monofacial PV panels are considered for the floating PV system. The modeling of irradiance on the front part of the panels is the same for both cases, while only rear irradiance applies for bifacial panels. As a first step, the angle of incidence on both sides of the panel is calculated, as in Equation (1). Equations (1)-(6) of the irradiance of PV panels are adapted from [51].
There are a total of three components of irradiance that fall on the front and rear sides of the bifacial panel, i.e., direct normal irradiance (DNI), diffused horizontal irradiance (DHI), and an albedo, due to reflected irradiance, as shown in Figure 9. The total irradiance on the PV panels is given by Equation (2).  The albedo further has two components, a reflected direct beam and reflected diffused irradiance, as given in Equation (3).
The albedo can be arrived at by integrating reflected irradiance along all the breadth of the row, as given in Equation (4).
where Fdl-gnd is the view factor given by   The albedo further has two components, a reflected direct beam and reflected diffused irradiance, as given in Equation (3).
The albedo can be arrived at by integrating reflected irradiance along all the breadth of the row, as given in Equation (4).
where F dl-gnd is the view factor given by where η f ront ,η Rear -efficiencies of front and rear cells p r_ f ront , p r_Rear -front &rear sides rating of panel Similar equations are applied for mono-facial PV panels used in floating solar PV systems, however without considering rear irradiance, as given in Equation (6).

Wind Energy Conversion System (WECS)
The city is not a good location for installing wind power plants at the ground level. However, there are some specific locations such as hilltops with more than 200 m elevation from the ground and the shoreline where a WECS can be installed at higher altitudes. Though they can accommodate a low number of turbines, the system will support the grid during peak hours in the evening, as shown in the wind profile, obtained after processing the data generated from the NREL Geospatial wind data of the proposed location and given in Figure 10. This reduces the need for higher BESS racks, helping in reducing the cost and maintenance of the system. cation and given in Figure 10. This reduces the need for higher BESS racks, helping in reducing the cost and maintenance of the system. The wind speeds at the proposed location are obtained from weather stations and are extrapolated to the required height using Equation (7). The power generated from the WECS can be obtained using Equation (8) [10]. Per the technical details of wind turbines provided in Section 3, a capacity use factor of ~40% can be achieved.
where is the surface roughness and is chosen to be 0.18 for the proposed location.

Battery Energy Storage System
As discussed in Section 1, a Li-ion-based BESS is considered in this study as BESSs The wind speeds at the proposed location are obtained from weather stations and are extrapolated to the required height using Equation (7). The power generated from the WECS can be obtained using Equation (8) [10]. Per the technical details of wind turbines provided in Section 3, a capacity use factor of~40% can be achieved.
where γ is the surface roughness and is chosen to be 0.18 for the proposed location.

Battery Energy Storage System
As discussed in Section 1, a Li-ion-based BESS is considered in this study as BESSs have more specific energy and can be scaled to racks and containers. Further, they have a good round-trip efficiency of 87-94%. The size of a BESS depends on the required number of hours of support. In this study, a maximum of one daily cycle is considered to support the grid for 4 hours during peak time, i.e., 7 PM to 11 PM. The technical details of the selected BESS are presented in Table 3, which is based on the Samsung SDI E3-R099 [52] model with a life cycle performance of 6000 cycles for continuous discharge of 1C at 25 • C. In this paper, the maximum limit for the number of racks is derived from the peak load, the number of hours, and the grid capacity.

Energy Management Strategy
To regulate the power flow and frequency in the power system network, the load demand must be met by the generating stations, as given in Equation (9) [10].
P L = P f pv + P rpv + P wt + P grid ± P bat (9) The grid is constrained by the following assumptions to ensure network reliability: i. The grid is never operated below the technical minimum of 55% of the load. ii. The grid is allowed to supply only 80% of the maximum load to accommodate future demand. iii. The ramp rates of the grid are restricted to ±0.5%/min of its capacity.
The quantity of power accessed from the grid follows Equation (10). The grid compensates the net load after the load is met by renewable energy. However, the grid is limited by ramping rates and then the maximum capacity.
The power mismatch between the load, HRES, and grid is given by Equation (11). P mis (t) = P 1 (t) − P ren (t) + P grid (t) (11) The flowchart for the EMS is shown in Figure 11. Initially, the size of the renewable energy system is fed to the EMS and the corresponding powers are calculated. The ramp Energies 2021, 14, 5368 14 of 23 rate limitations of the thermal power plants, battery power, and energy limitations are considered. Using these, the net power mismatch is calculated for the given load at any given instant. Then, power dispatch based on the two scenarios of power mismatch, excess load and excess generation, is carried out.   Figure 11. Flow chart for the energy management system.

Excess Load Scenario
In case the ramping rates or power drawn from the grid exceeds the maximum capacity, the BESS provides support to the system, thereby matching the load demand. Nevertheless, the BESS is also limited by the minimum and maximum power and energy constraints, as given by Equations (12) and (13), respectively, taken from [10].

Excess Load Scenario
In case the ramping rates or power drawn from the grid exceeds the maximum capacity, the BESS provides support to the system, thereby matching the load demand. Nevertheless, the BESS is also limited by the minimum and maximum power and energy constraints, as given by Equations (12) and (13), respectively, taken from [10]. P bat (t) = min P batmax (t), P mis (t) Then the mismatched power is updated as P mis (t) = P 1 (t) − P ren (t) + P grid (t) + P bat In the case of any further supply-demand mismatch, the demand response will have to set in and curtail the non-essential loads.

Excess Generation Scenario
In the case of excess generation, the grid intake is reduced until the mismatch nullifies or the grid reaches its technical minimum. If the excess supply persists, then the BESS will charge. As before, the battery intake is limited by its constraints on charging power and energy, as in Equations (14) and (15), taken from [10].
If there is any excess generation even after the battery capacity reaches a maximum, the excess energy is wheeled to other load centers.

Multi-Objective Adaptive-Local-Attractor-Based Quantum-Behaved Particle Swarm Optimization (ALA-QPSO)
In the original PSO, the position and velocity of the i-th particle in each iteration 'it' are updated according to Equation (16). As discussed in Section 1, the trajectory of these particles is influenced by an attractor [21].
Using the concept of particle trajectory, the principle of uncertainty from quantum mechanics is applied to generate a delta potential well that allows the particles to move across the entire search space. The position pos of the j-th dimension in the i-th particle is updated using the adaptive local attractor in each iteration according to Equation (17) [21]. The first term is updated based on the adaptive local attractor given in Equation (18), which is updated based on an adaptive contraction-expansion coefficient in Equation (19) and the diversity of the fitness of the population fit ij from the mean fitness fit mean in Equation (20). The second term in Equation (17) is updated according to the weighted mean personal best position of all the particles, as given in Equation (21).
where atr ij (it) is the adaptive local attractor, as defined in where pop is the population size, pbest ij (it) indicates the personal best of the i-th particle during the it-th iteration, and gbest j (it) is the global best of all the particles in the itth iteration. The contraction-expansion coefficient is given by where α 0 is the initial weight and it max is the maximum number of iterations. The diversity factor in each iteration, div(it), is given by Equation (20).
The weighted mean personal best of each particle in the it-th iteration, mpbest(it), is given in Equation (21).
Then the fitness of the particles is evaluated until convergence is achieved or until maximum iterations are completed. The Algorithm 1 for the same is given below.
Initialize the parameters and particles in the population with random position vectors. ii.
Evaluate the fitness of all the particles and calculate the personal best for each particle and the global best for the swarm. iii.
Calculate the population diversity as in Equation (20). iv.
Update the weighted mean personal-best position using Equation (21). v.
Update the adaptive local attractor for each particle as in Equation (18) vi.
Update the position for each particle using Equation (17). vii. If the terminating condition is not met, go to step iii.
In the case of ALA-MQPSO, a mutation strategy is introduced after step vi to enhance the exploitation capability of the algorithm. For optimal sizing of the HRES, three objective functions, LPSP, LCoE [10], and LCL [54] of the BESS, are considered and are given in Equations (22)- (24). The LPSP represents the amount of load that is cut off by the demand response control, while the LCL represents the usage of the BESS. A higher LCL means frequent cycling of the BESS, which reduces the life of the battery pack even before its rated end-of-life period.
where f cr = r(1+r) l (1+r) l −1 is the discount factor, r is the interest rate, l is the lifetime of the project, C i is the capital expenditure, and OM i is the cost of operation and maintenance.
where kp is the Peukert lifetime constant, N f ail 100 is the maximum number of cycles that a battery experiences before its rated end of life, and n eq 100 is equivalent 100%-DOD cycle number of K half cycles.

Simulation Results and Analysis
Based on the assessed potential in Section 2, optimal sizing of the HRES is achieved to minimize the LPSP, LCoE, and LCL. ALA-QPSO techniques are applied, and the results are compared with four other benchmark techniques, four variants of differential evolution (DE) algorithms. The adopted algorithm is compared with variants of DE as it has proven to be an effective optimization technique for enhancing the exploration of the population. The optimization parameters for the algorithms under comparison are presented in Table 4. The same values are chosen for parameters, wherever necessary, for all the algorithms to prove the effectiveness of the proposed algorithm. The results are screened further based on the minimum LPSP, as shown in Table 5. The result obtained is shown to be near optimal as the combination fared well in most of the constraints that include the least number of hours of load loss and the least number of hours of excess generation. This will result in better profitability of the proposed system.  A comparison of the Pareto fronts of three optimization techniques is depicted in Figure 12. Further, the power distribution among the sources over 5 days in each quarter of a year, i.e., 21 to 26 March, June, September, and December, is shown in Figures 13-16 with the minimum LPSP configuration. The same is shown with the minimum LCOE configuration for June in Figure 17. The results prove that the use of the BESS is limited only for a transition from peak to off-peak hours in most of the cases, except in June since it is the monsoon season.
During this time, the BESS has helped in maintaining the frequency of the grid constant by compensating for any generation loss by the PV system. It should also be noted that the participation of battery power in the proposed system is mostly during the violation of ramp rate limits, in the early morning or in the evening. Hence, we notice that at noon, ramp rates are not violated the power dispatch from the grid is at its technical minimum. This shows that the higher penetration of variable renewable energy sources into the grid brings down the plant load factor of thermal power plants. Hence, if BESSs are not used during the violation of ramp rate limits, it may burden the TPPs further in the event of phasing-out of diesel-or gas-based plants. Similarly, the power distribution for June with mutation-based ALA-QPSO is plotted in Figure 18 for the minimum LCoE configuration. Here, the month of June is selected, as solar power undergoes wide changes during that month. stant by compensating for any generation loss by the PV system. It should also be noted that the participation of battery power in the proposed system is mostly during the violation of ramp rate limits, in the early morning or in the evening. Hence, we notice that at noon, ramp rates are not violated the power dispatch from the grid is at its technical minimum. This shows that the higher penetration of variable renewable energy sources into the grid brings down the plant load factor of thermal power plants. Hence, if BESSs are not used during the violation of ramp rate limits, it may burden the TPPs further in the event of phasing-out of diesel-or gas-based plants. Similarly, the power distribution for June with mutation-based ALA-QPSO is plotted in Figure 18 for the minimum LCoE configuration. Here, the month of June is selected, as solar power undergoes wide changes during that month.   stant by compensating for any generation loss by the PV system. It should also be noted that the participation of battery power in the proposed system is mostly during the violation of ramp rate limits, in the early morning or in the evening. Hence, we notice that at noon, ramp rates are not violated the power dispatch from the grid is at its technical minimum. This shows that the higher penetration of variable renewable energy sources into the grid brings down the plant load factor of thermal power plants. Hence, if BESSs are not used during the violation of ramp rate limits, it may burden the TPPs further in the event of phasing-out of diesel-or gas-based plants. Similarly, the power distribution for June with mutation-based ALA-QPSO is plotted in Figure 18 for the minimum LCoE configuration. Here, the month of June is selected, as solar power undergoes wide changes during that month.              A sensitivity analysis is performed on the results obtained through ALA-QPSO. In the first scenario, it is assumed that there is 50% damage to the BESS. The power distribution with this scenario is shown in Figure 19 for June. It can be observed that the demand response is enforced during peak hours. It should be noted that the LCoE has increased drastically to 0.1372 USD, which is 73% more. In addition, the LPSP has increased to 1.62%. We conclude that preventive maintenance is necessary to avoid any load loss due to BESS reduction. Another sensitivity analysis is performed for 20% demand appreciation and 20% BESS reduction. The results for the same are shown in Figure 20. It is observed that due to the rise in demand, ramping rates of the grid are maxed out, and hence the BESS supports the grid during sunset, leading to less backup power for the peak time. It is observed that the LCoE has increased drastically to 0.1577 USD/kW, which is twice that of the actual, and the LPSP has increased to 1.85%. Further, in the case of reduced solar power, the BESS cannot compensate for the generation, even when the ramping is maxed out. It can be concluded that for reliable operation, the size of the BESS should be increased when the LCoE configuration is considered. A sensitivity analysis is performed on the results obtained through ALA-QPSO. In the first scenario, it is assumed that there is 50% damage to the BESS. The power distribution with this scenario is shown in Figure 19 for June. It can be observed that the demand response is enforced during peak hours. It should be noted that the LCoE has increased drastically to 0.1372 USD, which is 73% more. In addition, the LPSP has increased to 1.62%. We conclude that preventive maintenance is necessary to avoid any load loss due to BESS reduction. Another sensitivity analysis is performed for 20% demand appreciation and 20% BESS reduction. The results for the same are shown in Figure 20. It is observed that due to the rise in demand, ramping rates of the grid are maxed out, and hence the BESS supports the grid during sunset, leading to less backup power for the peak time. It is observed that the LCoE has increased drastically to 0.1577 USD/kW, which is twice that of the actual, and the LPSP has increased to 1.85%. Further, in the case of reduced solar power, the BESS cannot compensate for the generation, even when the ramping is maxed out. It can be concluded that for reliable operation, the size of the BESS should be increased when the LCoE configuration is considered.

Conclusions
With the increase in energy demands in smart cities, it is more important than ever to integrate renewable energy sources to stop relying on traditional energy sources. In this study, a hybrid renewable energy system with bifacial rooftop PV, floating PV, and wind energy systems, as well as lithium-ion batteries for energy storage, is considered. The renewable energy potential of a city in India is assessed using Google Earth Pro's Polygon tool, CAD mapper, and SketchUp Pro. The real-time 1-hour interval-based load profile of the selected locations is considered for calculating the upper limits of the BESS. Reliability indices such as loss of load probability (LPSP), levelized cost of energy (LCoE), and life cycle loss (LCL) of the BESS are considered for optimizing the size under certain constraints on the grid. The optimal sizing of the available HRES is achieved using the ALA-QPSO technique and is compared with benchmark algorithms such as variants of the differential algorithm. A mutation strategy is applied to the adopted ALA-QPSO to enhance exploitation. Though the results of this study are superior, they need further investigation. A sensitivity analysis is performed for conditions such as demand appreciation, battery damage, and cost escalation. The results prove that with the financial fitness of the municipal corporation of the region, the proposed system can be implemented in a phased manner with the help of financial institutions such as the Asian Development Bank, with which the municipal corporation has already tied up for some

Conclusions
With the increase in energy demands in smart cities, it is more important than ever to integrate renewable energy sources to stop relying on traditional energy sources. In this study, a hybrid renewable energy system with bifacial rooftop PV, floating PV, and wind energy systems, as well as lithium-ion batteries for energy storage, is considered. The renewable energy potential of a city in India is assessed using Google Earth Pro's Polygon tool, CAD mapper, and SketchUp Pro. The real-time 1-hour interval-based load profile of the selected locations is considered for calculating the upper limits of the BESS. Reliability indices such as loss of load probability (LPSP), levelized cost of energy (LCoE), and life cycle loss (LCL) of the BESS are considered for optimizing the size under certain constraints on the grid. The optimal sizing of the available HRES is achieved using the ALA-QPSO technique and is compared with benchmark algorithms such as variants of the differential algorithm. A mutation strategy is applied to the adopted ALA-QPSO to enhance exploitation. Though the results of this study are superior, they need further investigation. A sensitivity analysis is performed for conditions such as demand appreciation, battery damage, and cost escalation. The results prove that with the financial fitness of the municipal corporation of the region, the proposed system can be implemented in a phased manner with the help of financial institutions such as the Asian Development Bank, with which the municipal corporation has already tied up for some projects. To assess the operational feasibility of the system, stability analysis of the system can be performed in the future, given the transmission data. This work can be extended to other cities by assessing the potential of RES feasible in those areas. This approach will help the world in achieving the RE targets within the time limits.