Energy Non-Availability in Distribution Grids with Heavy Penetration of Solar Power: Assessment and Mitigation through Solar Smoother

: Rapid ﬂuctuation of solar irradiance due to cloud passage causes corresponding variations in the power output of solar PV power plants. This leads to rapid voltage instability at the point of common coupling (PCC) of the connected grid which may cause temporary shutdown of the plant leading to non-availability of energy in the connected load and distribution grid. An estimate of the duration and frequency of this outage is important for solar energy generators to ensure the generation and performance of the solar power plant. A methodology using PVsyst (6.6.4, University of Geneva, Geneva, Switzerland) and PSCAD (4.5, Manitoba HVDC Research Centre, Winnipeg, MB, Canada) simulation has been developed to estimate the duration and frequency of power outages due to rapid ﬂuctuation of solar irradiance throughout the year. It is shown that the outage depends not only on the solar irradiance ﬂuctuation, but also on the grid parameters of the connected distribution grid. A practical case study has been done on a 500 kilo Watt peak (kWp) solar PV power plant for validation of the proposed methodology. It is observed that the energy non-availability for this plant is about 13% per year. This can be reduced to 8% by incorporating a solar smoother. A ﬁnancial analysis of this outage and its mitigation has also been carried out.


Introduction
As a follow up of Paris Conference agreement many countries around the world including India have initiated massive programs for establishing thousands of Mega Watt (MW) solar photovoltaic (PV) installations throughout the country in next five to ten years. It is envisaged that hundreds of kilo Watt (kW) of decentralized solar PV power plants in the range of 10 kW to 500 kW or more will be installed in these countries and connected to the local distribution grids. It is also noticed that the solar irradiance level fluctuates seasonally in a significant manner due to the local weather conditions. Such a high penetration level of solar PV and the high rate of change of PV output power due to fluctuation in irradiance can destabilize the distribution system in terms of voltage level leading to outages and unavailability of electricity to the connected loads. However very few studies have been done so far to take account of the frequency and duration of these outages since under the present practice and regulations it is essential to implement remote monitoring system to monitor the inverter parameters such as instantaneous power and energy generation. No provision has been incorporated to monitor and log the voltage at the point of common coupling (PCC), so it is difficult to calculate the energy losses due to the grid voltage outage. In this paper a methodology has been proposed by which the availability of solar power is confirmed at the PCC by incorporating a solar smoother along with storage within the system. Two professional software packages have been used, namely PSCAD (4.5, Manitoba HVDC Research Centre, Winnipeg, MB, Canada) and PVsyst (6.6.4, University of Geneva, Geneva, Switzerland), to simulate the grid voltage instability at the PCC and yearly energy generation from a solar PV plant respectively. Kamaruzzaman et al. [1] studied the impact of a grid-connected photovoltaic (PV) system on static and dynamic voltage stability. To quantify the effect an analytical technique has also been implemented. Intermittency of high ramp-rated solar irradiance can cause fluctuation in the PV output due to cloud passing, as studied by Alam et al. [2]. The variation of voltage is more predominant for a weak distribution grid with high PV penetration. This fluctuation can be smoothened by a proper control mechanism of the energy storage devices. In their work, Jewell and Ramakumar [3] discussed that the power captured by a utility-based dispersed solar photovoltaic generator varied with the fluctuation of incident solar radiation on the utility service area. A load flow program has been simulated to study the impact of distributed residential PV generation on the utility. Maximum PV penetration relative to steady-state voltage and overcurrent has been found in different locations at the distribution feeder investigated by Hoke et al. [4]. A mitigation technique has also been proposed. Transient and steady-state performance of a photovoltaic system connected to a large utility grid have been investigated by Wang and Lin [5] using a novel approach based on an eigen value method. In their work, Srisaen and Sangswang [6] studied the consequences of PV system installation location on the power quality of the distribution system. It has been concluded from the simulation results that implementation of the PV grid-connected system could improve the power quality of a distribution feeder. To increase the reliability of PV-based generation some modifications to existing inverter technology have been proposed by Johnson et al. [7]. The significant instantaneous ramp rates have been investigated using both the standard method of calculating ramp rates at discrete time intervals and using the piecewise linear approximation method. Jayasekara et al. [8] in their work have proposed a management technique for distributed battery storage that enhances the capability of distribution networks to utilize PV generation. For energy forecasts optimal daily energy profiles for each storage unit have also been calculated on a one day ahead schedule. The impact of high PV penetration on low voltage distribution system with high X/R ratio, and long tap changing delay has been studied by Yan and Saha [9]. They also proposed that the distribution networks can be operated with low PV penetration without voltage stability problem due to cloud passing. For high penetration the reactive power support from a PV inverter can be used to mitigate this voltage stability. In his work, Marcos et al. [10] has discussed a new control scheme (step control) in relation to two strategies-ramp control and moving average-to smooth out PV power fluctuations due to cloud coverage. A technique to calculate the storage capacity requirement for three different control strategies has also been implemented. The design of control systems for PV systems with extremely fast variation of weather conditions and necessary regulations has been investigated by Yan and Saha [11]. Fast maximum power point tracking (MPPT) is not an appropriate option to solve the high penetration problem because it may create power surges with solar radiation increases of which can create network overvoltage, and breaker tripping problems, etc., so, a new control scheme has been described to control PV power ramping up and whole system becomes smooth. The adverse impact of large scale solar PV installation on the utility grid has been discussed by Omran et al. [12]. This study also reveals that if the penetration level of the solar PV becomes high then this can have a negative impact on the grid network. Three methods to mitigate these problems have been proposed such as connecting battery storage with the existing PV plant, connecting dump loads and curtailment of power from the solar PV source. By implementing these methods the financial benefits at the consumers' end have been discussed. Barelli et al. [13] have discussed the requirement of energy storage for grid stability and safety due to heavy penetration of solar PV installations. This paper simulated the dynamic behaviour of a hybrid energy storage system coupled to a solar PV system and residential load. The result summarizes the evaluation of instantaneous behavior of each component along with the impact of fluctuation of the power on the battery and grid. The effects of passing clouds on the utility-connected battery backup PV system has been studied through radial basis function and back propagation network techniques by Giraud and Salameh [14]. Moreover, the effect of random clouds on the developed network model is simulated through an artificial neural network and standard linear model and compared with the experimental results. This model can be used for designing and engineering utility-connected PV systems. Chalmers et al. evaluated the effect of solar PV generation on utility operation [15]. The characteristics of various performance parameters of the system were also examined. Datta et al. [16] discussed three different control methods for PV output power smoothing to reduce frequency fluctuation, first is fuzzy reasoning-based, the second is an energy storage-based one and the third is a simple coordination-based control method. The results of the three studies have been compared with the MPPT-based control method and found effective in mitigating fluctuations. A Gaussian-based algorithm for smoothing the power fluctuation from renewable energy-based plants has been proposed by Addisu et al. [17]. By implementing this algorithm the optimum sizing of the Energy Storage System (ESS) has also been determined and compared the result with the moving average method. Up to a 34% improvement in power smoothness and a 19% improvement of battery sizing has been achieved. Different innovative strategies for minimizing the ESS capacity for MW level PV plants have been proposed by Parra et al. in this work [18]. The results of the proposed strategies have also been validated through simulation and operational PV plant data. It is easier for transmission system operators to handle short term power fluctuations due to intermittency of the PV system. Recently Sarkar et al. [19] considered the dependence of voltage instability in the distribution grid on the ratio of PV penetration and connected load. However to the best of the knowledge of the authors no work has been reported so far where the loss of energy availability in the distribution grid due to the irradiance fluctuation of the connected PV power plant has been computed methodically and the role of solar smoother in minimizing the loss of energy has been analyzed in depth. In this paper a methodology is first proposed to estimate the of energy non-availability in the distribution grid with heavy penetration of solar power and high fluctuation of irradiance. Next a mitigation technique is demonstrated by designing a suitable solar smoother. Finally a detailed case study is undertaken to validate the proposed methodology.
The impact of high penetration of solar PV on the electrical dynamics of the distribution grid has been estimated using the PSCAD simulation platform. The amount of energy non-availability at the grid end with and without solar smoother is demonstrated by using a PVsyst simulation. It has been found that the energy availability is improved by 4.86% throughout the year after using the proposed solar smoother technology.
Further, the proposed solar smoother technology helps in minimizing significant electricity tariff financial losses by mitigating the probability of energy non-availability at the consumer end. The rest is divided into seven sections. This Section 1 describes the Introduction. Section 2 describes the effect of ramp up and ramp down of irradiance which leads to a voltage fluctuation in the distribution grid. Section 3 describes the modeling of the system. Section 4 describes the methodology of the system. Section 5 describes the case study of 500 kWp solar PV system installed in Kolkata, India, where the energy non-availability has been estimated using the proposed methodology and then validated though actual observations of field data. Section 6 discusses the financial issues. Section 7 presents the conclusions of the paper.

Concept
Sudden variations of solar irradiance due to cloud passing are called ramp up and ramp down phenomena. When the ramp up and ramp down of irradiance occurs within a short amount of time, the inverter output power also fluctuates suddenly and voltage instability occurs at the point of common coupling of the distribution grid. To prevent the possibility of damage caused by instability, the usual practice is to curtail the connected load which leads to non-availability of energy. This method is not very efficient. A state of the art method called solar smoothing has been introduced to handle this situation. Due to increasing usage and penetration of solar PV in distribution grid the occurrence of such instability increases. On the other hand the solar energy providers have to guarantee the energy availability to the grid for long period of 20 years or more. Therefore it is most important for a solar energy provider to estimate correctly the frequency and duration of non-availability for a particular power plant over its lifetime and also use of smoothing for mitigation of this problem. Figure 1 shows the characteristic ramp down curve of the output power of a grid-connected inverter. The characteristic curve is defined as the desired output behavior of the grid-connected inverter, initially operating at nominal rated output, in response to a step change in renewable energy source output at time t = 0 s to 0 kW. the inverter output power also fluctuates suddenly and voltage instability occurs at the point of common coupling of the distribution grid. To prevent the possibility of damage caused by instability, the usual practice is to curtail the connected load which leads to non-availability of energy. This method is not very efficient. A state of the art method called solar smoothing has been introduced to handle this situation. Due to increasing usage and penetration of solar PV in distribution grid the occurrence of such instability increases. On the other hand the solar energy providers have to guarantee the energy availability to the grid for long period of 20 years or more. Therefore it is most important for a solar energy provider to estimate correctly the frequency and duration of nonavailability for a particular power plant over its lifetime and also use of smoothing for mitigation of this problem. Figure 1 shows the characteristic ramp down curve of the output power of a grid-connected inverter. The characteristic curve is defined as the desired output behavior of the grid-connected inverter, initially operating at nominal rated output, in response to a step change in renewable energy source output at time t = 0 s to 0 kW. The ramp down phenomena can be established mathematically as:

Characteristic Curve
and, where Rn is the nominal ramp down rate (W/s), Tn is the nominal ramp down time (s), On is the nominal rated inverter output (kW), Δ is the maximum input deviation (s) and NL is the percentage of no linearity (%). The nominal ramp down rate for each installation size is then derived according to Equation (2). Control action is required when the ramp down rate of the renewable energy source exceeds Rn in magnitude.
Similarly, the nominal ramp up time Rp is defined as the time taken for the inverter to ramp up from 0 kW to the nominal rated output, given by the following expression: The ramp down phenomena can be established mathematically as: and, where R n is the nominal ramp down rate (W/s), T n is the nominal ramp down time (s), O n is the nominal rated inverter output (kW), ∆ is the maximum input deviation (s) and NL is the percentage of no linearity (%). The nominal ramp down rate for each installation size is then derived according to Equation (2). Control action is required when the ramp down rate of the renewable energy source exceeds R n in magnitude. Similarly, the nominal ramp up time R p is defined as the time taken for the inverter to ramp up from 0 kW to the nominal rated output, given by the following expression: Figure 2 describes the ramp up phenomena. This covers for sudden increase in renewable source electricity output. Control action is required when the ramp up rate of the renewable energy source exceeds R p in magnitude.   Figure 3 indicates the schematic diagram of a solar PV system connected with a low voltage distribution grid and load. The PV and load are connected at point of common coupling (PCC) and the grid is connected through a distribution line to the PCC. The power generated from solar PV will be consumed by the load and the excess power, if any, will be exported to the grid. VPCC indicates the voltage at the PCC at PV end, VG indicates the grid voltage at point of supply from the distribution company, I is the PV current, Z is the impedance of the distribution line at the consumer end. The voltage deviation (∆V) can be expressed as:

Mathematical Model of the System
Now, ΔV can be represented by Equation (5):  Figure 3 indicates the schematic diagram of a solar PV system connected with a low voltage distribution grid and load.   Figure 3 indicates the schematic diagram of a solar PV system connected with a low voltage distribution grid and load. The PV and load are connected at point of common coupling (PCC) and the grid is connected through a distribution line to the PCC. The power generated from solar PV will be consumed by the load and the excess power, if any, will be exported to the grid. VPCC indicates the voltage at the PCC at PV end, VG indicates the grid voltage at point of supply from the distribution company, I is the PV current, Z is the impedance of the distribution line at the consumer end. The voltage deviation (∆V) can be expressed as:

Mathematical Model of the System
Now, ΔV can be represented by Equation (5): where ϕ is the angle between the load voltage and load current, RL is the load resistance (Ω) and α is the angle between VPCC and I. The PV and load are connected at point of common coupling (PCC) and the grid is connected through a distribution line to the PCC. The power generated from solar PV will be consumed by the load and the excess power, if any, will be exported to the grid. V PCC indicates the voltage at the PCC at PV end, V G indicates the grid voltage at point of supply from the distribution company, I is the PV current, Z is the impedance of the distribution line at the consumer end. The voltage deviation (∆V) can be expressed as: Now, ∆V can be represented by Equation (5): where φ is the angle between the load voltage and load current, R L is the load resistance (Ω) and α is the angle between V PCC and I. For ramp down phenomena Equation (6) can be written as: For ramp up phenomena Equation (7) can be written as: Therefore, Equations (6) and (7) describe that with a high ramp down rate of the inverter output power, the voltage at the PCC will dip, whereas with a high ramp up rate of the inverter output power the voltage at the PCC will rise.

Methodology
The methodology adopted for this work is as follows: first the month wise daily average irradiance data of the site is to be collected from a state of the art weather monitoring station located nearby. Otherwise National Aeronautics and Space Administration (NASA) or Meteonorm data may be used for this purpose, as shown in Figure 4. Various research papers [13][14][15] have also been published so far regarding the forecasting of solar power which can be also very useful for designing the prediction model. For ramp down phenomena Equation (6) can be written as: For ramp up phenomena Equation (7) can be written as: Therefore, Equations (6) and (7) describe that with a high ramp down rate of the inverter output power, the voltage at the PCC will dip, whereas with a high ramp up rate of the inverter output power the voltage at the PCC will rise.

Methodology
The methodology adopted for this work is as follows: first the month wise daily average irradiance data of the site is to be collected from a state of the art weather monitoring station located nearby. Otherwise National Aeronautics and Space Administration (NASA) or Meteonorm data may be used for this purpose, as shown in Figure 4. Various research papers [13][14][15] have also been published so far regarding the forecasting of solar power which can be also very useful for designing the prediction model. The monthly average of the daily variation of power output of the power plant during the daytime is to be collected from the remote monitoring system of the power plant under consideration and the rate of change of power on the positive as well as the negative side is to be calculated. Figure  5 shows the raw data of the remote monitoring unit of the power plant under consideration for the month of March 2016. The monthly average of the daily variation of power output of the power plant during the daytime is to be collected from the remote monitoring system of the power plant under consideration and the rate of change of power on the positive as well as the negative side is to be calculated. Figure 5 shows the raw data of the remote monitoring unit of the power plant under consideration for the month of March 2016.
The monthly average of the daily variation of power output of the power plant during the daytime is to be collected from the remote monitoring system of the power plant under consideration and the rate of change of power on the positive as well as the negative side is to be calculated. Figure  5 shows the raw data of the remote monitoring unit of the power plant under consideration for the month of March 2016.  The maximum and minimum ramp up and ramp down of the inverter output power and the percentage variation of the capacity can also be estimated from this raw data, as shown in Figure 5, but the voltage fluctuation at the PCC cannot be obtained directly.
Standard simulation software like PSCAD has been utilized in this work to estimate the voltage fluctuation. However for this purpose the simulation of the PSCAD needs to be validated with the field data of the power output variation. The X and R values of the distribution line are needed for carrying out the PSCAD simulation, which are to be collected from the distribution grid authority.
The voltage fluctuation at the PCC is estimated due to the heavy ramp up and ramp down of irradiance. The frequency and duration of voltage outages throughout the year can thus be estimated, from where the percentage loss of power supply is calculated. By using the PVsyst software the yearly generation has to be determined considering the calculated percentage of loss of power supply.
Similarly, by incorporating the solar smoother along with battery storage within the system, the voltage fluctuation at the PCC has to be identified and the frequency of voltage outage estimated through PSCAD simulation. The energy generation throughout the year has to be estimated by PVsyst simulation and compared with the previous one. Figure 6 shows the flow chart of the whole work. The maximum and minimum ramp up and ramp down of the inverter output power and the percentage variation of the capacity can also be estimated from this raw data, as shown in Figure 5, but the voltage fluctuation at the PCC cannot be obtained directly.
Standard simulation software like PSCAD has been utilized in this work to estimate the voltage fluctuation. However for this purpose the simulation of the PSCAD needs to be validated with the field data of the power output variation. The X and R values of the distribution line are needed for carrying out the PSCAD simulation, which are to be collected from the distribution grid authority.
The voltage fluctuation at the PCC is estimated due to the heavy ramp up and ramp down of irradiance. The frequency and duration of voltage outages throughout the year can thus be estimated, from where the percentage loss of power supply is calculated. By using the PVsyst software the yearly generation has to be determined considering the calculated percentage of loss of power supply.
Similarly, by incorporating the solar smoother along with battery storage within the system, the voltage fluctuation at the PCC has to be identified and the frequency of voltage outage estimated through PSCAD simulation. The energy generation throughout the year has to be estimated by PVsyst simulation and compared with the previous one. Figure 6 shows the flow chart of the whole work.

Case Study
To validate the abovementioned methodology a case study has been performed. As shown in Figure 7, a 500 kWp grid tied solar PV plant located at Newtown Kolkata (India) was commissioned in 2015. It consists of a 500 kWp solar PV array, 25  20 kVA grid tied string inverters and 11 kV/0.415 kV, 650 kVA distribution transformers for connecting to the Newtown Electric Supply Company Ltd. (West Bengal, India) (NTESCL)'s local distribution network through a suitable protection system.

Case Study
To validate the abovementioned methodology a case study has been performed. As shown in Figure 7, a 500 kWp grid tied solar PV plant located at Newtown Kolkata (India) was commissioned in 2015. It consists of a 500 kWp solar PV array, 25 × 20 kVA grid tied string inverters and 11 kV/0.415 kV, 650 kVA distribution transformers for connecting to the Newtown Electric Supply Company Ltd. (West Bengal, India) (NTESCL)'s local distribution network through a suitable protection system.       As the fluctuation of solar irradiance shows a maximum in summer and a minimum in winter, Figure 8 has been arranged accordingly e.g., it starts from March and ends in February. The irradiance data has been collected from the Solar Radiation Resource Assessment (SRRA) unit installed at the Indian Institute of Engineering Science and Technology (IIEST), Shibpur (India) campus. It is shown that the maximum rate of change of irradiance on the positive as well as negative side occurs mostly in the summer season and the maximum values are found in the months of March, May and June.

Month Wise Daily Average Inverter Output Power
In this work, a 1 min ramp rate is defined as the absolute value of the difference between the instantaneous power at the start and at the end of a 1 min period. The ramp rates over different time scales (∆t) can be calculated using the equation below: where P(t) is the time history of the inverter output power to the grid with measurements taken at t (in s) and ∆t is the time interval. The daily average power output of each month from the plant measured at the point of common coupling (PCC) is shown in Figure 9.
It also shows the rate of change of power/min and the change of percentage capacity of the plant/min, on both the positive and negative side Tables 1 and 2 show the month wise ramp up and ramp down of inverter output power which has been divided in six surge zones (SZ) e.g., from SZ-1 (10% change of capacity/min) to SZ-6 (60% change of capacity/min).
It is noticed from Tables 1 and 2 that the low ramp up and ramp down of inverter power (up to 30%) occurs mostly in every month whereas the high ramp up (up to 50%) occurs in the month of March and the high ramp down (up to 40%, 50% and 60%) occurs in the months of May, June and March, respectively.

PSCAD Simulation
The above said 500 kWp solar PV system connected to the low voltage NTSCL distribution grid has been simulated through the PSCAD software. The simulation describes the level of power surge of the inverter output within the distribution network with values of the X/R ratio of the distribution line. A boundary limit of the grid voltage has also been identified beyond which the grid becomes unstable. A test network model for low voltage distribution line has been developed in PSCAD. The X and R values (0.205 + j0.074) of the network are collected from NTSCL and fed to the simulation along with the irradiance data. The simulation has been done with six different power surge zones to observe the voltage fluctuation level at PCC. Table 3 shows the comparison between the percentage change of capacity/min between real time data obtained from 500 kWp solar PV plant and PSCAD simulation along with the percentage deviation for three months, namely March, May and June. It shows that the maximum deviation on the positive side is 7.7% and on the negative side it is −9.1%, which are well between ±10%. This validates the simulation results. Figure 10 shows the inverter output power and voltage at PCC for the months of March, May and June which cover the six surge zones described above. It clearly shows the voltage swelling and sagging at the PCC at the time of ramp up and ramp down of the inverter output power. Table 4 shows the simulated values of voltage at the PCC against the percentage change of inverter output (in power/min). The portions marked in red in the table show the swelling and sagging values of the PCC voltages which are beyond the tolerable limit on both the positive and negative side. The tolerable voltage limit has been considered as per Indian Grid Code i.e., ±10% of rated voltage (400 V AC). From Table 3 the ramp up and ramp down allowable limits of the inverter output power have been calculated as 17%/min and 46%/min, respectively. AC). From Table 3 the ramp up and ramp down allowable limits of the inverter output power have been calculated as 17%/min and 46%/min, respectively.   AC). From Table 3 the ramp up and ramp down allowable limits of the inverter output power have been calculated as 17%/min and 46%/min, respectively.

Estimation of Energy Non Availability through PVsyst Simulation without Smoothing
In Table 5 the portions marked in blue show the frequency of occurrence of ramp up and ramp down of the inverter output power/min beyond the tolerable limit on the positive as well as the negative side, which is calculated as 13.19% throughout the year. Conversely, it can be considered that the year-round energy availability of the system is 86.81%. A PVsyst-based simulation has been done to estimate the power generation throughout the year considering the loss of load 13.19%. This is illustrated in Figure 11. Energies 2018, 11, x 13 of 19 Figure 11. PVsyst simulation of a grid connected 500 kWp solar PV system.
It is shown from the figure that the total energy generation is 687.4 MWh against a Performance Ratio (PR) value of 75.1%. The energy unavailability from the system is 19 days in a year.

Proposed Implication on Energy Estimation with Smoother Incorporating Energy Storage
A proposal has been presented where by incorporation of battery storage within the system, the availability of the power can be increased. At the time of ramp up the storage device will be in discharging mode and vice versa. The schematic of the system is indicated in Figure 12. The system has been simulated in PSCAD maintaining the other parameters the same. The capacity of the storage can be designed to mitigate the maximum ramp up and ramp down rate throughout the year. Tables 1 and 2 show that the maximum ramp rate of power is 248 kW/min on the positive side and 253.5 kW/min on the negative side. The maximum energy content within the surge will be: where, Psurge = surge power in kW, and t = time up to when the surge exists in minutes. It is shown from the figure that the total energy generation is 687.4 MWh against a Performance Ratio (PR) value of 75.1%. The energy unavailability from the system is 19 days in a year.

Proposed Implication on Energy Estimation with Smoother Incorporating Energy Storage
A proposal has been presented where by incorporation of battery storage within the system, the availability of the power can be increased. At the time of ramp up the storage device will be in discharging mode and vice versa. The schematic of the system is indicated in Figure 12. It is shown from the figure that the total energy generation is 687.4 MWh against a Performance Ratio (PR) value of 75.1%. The energy unavailability from the system is 19 days in a year.

Proposed Implication on Energy Estimation with Smoother Incorporating Energy Storage
A proposal has been presented where by incorporation of battery storage within the system, the availability of the power can be increased. At the time of ramp up the storage device will be in discharging mode and vice versa. The schematic of the system is indicated in Figure 12. The system has been simulated in PSCAD maintaining the other parameters the same. The capacity of the storage can be designed to mitigate the maximum ramp up and ramp down rate throughout the year. Tables 1 and 2 show that the maximum ramp rate of power is 248 kW/min on the positive side and 253.5 kW/min on the negative side. The maximum energy content within the surge will be: where, Psurge = surge power in kW, and t = time up to when the surge exists in minutes.

Figure 12.
Typical system diagram of a system with energy storage and ramp rate controller.
The system has been simulated in PSCAD maintaining the other parameters the same. The capacity of the storage can be designed to mitigate the maximum ramp up and ramp down rate throughout the year. Tables 1 and 2 show that the maximum ramp rate of power is 248 kW/min on the positive side and 253.5 kW/min on the negative side. The maximum energy content within the surge will be: where, P surge = surge power in kW, and t = time up to when the surge exists in minutes. From Equation (9) E surge has been calculated as 4225 units. Now, this amount of energy will be absorbed by the battery storage at the time of ramp up and will be released at the time of ramp down phenomena. Moreover, from Table 5 it is noticed that a maximum three times a day in the month of April the power is fluctuating beyond the limit. Therefore to handle this surge power the energy at battery needs to be stored as: The capacity of the battery has been calculated as: where, V battery = terminal voltage of the battery bank.
In this case 684V DC is taken as the string voltage of the PV is nearly 700 V DC. Therefore from Equation (11) the capacity of the battery has been calculated as 18.53 Ah. Now, considering up to 85% depth of discharge (DOD) of the battery and availability in the market, the capacity has been taken as 22 Ah. In this case a Gel Long Life battery has been considered. Fifty seven 12 V, 22 Ah battery modules have been connected in series to provide a terminal voltage of 684 V DC [20].
After incorporation of the smoother and storage the output of the inverter power and voltage at the PCC is as shown in Figure 13. From Equation (9) Esurge has been calculated as 4225 units. Now, this amount of energy will be absorbed by the battery storage at the time of ramp up and will be released at the time of ramp down phenomena. Moreover, from Table 5 it is noticed that a maximum three times a day in the month of April the power is fluctuating beyond the limit. Therefore to handle this surge power the energy at battery needs to be stored as: The capacity of the battery has been calculated as: where, Vbattery = terminal voltage of the battery bank.
In this case 684V DC is taken as the string voltage of the PV is nearly 700 V DC. Therefore from Equation (11) the capacity of the battery has been calculated as 18.53 Ah. Now, considering up to 85% depth of discharge (DOD) of the battery and availability in the market, the capacity has been taken as 22 Ah. In this case a Gel Long Life battery has been considered. Fifty seven 12 V, 22 Ah battery modules have been connected in series to provide a terminal voltage of 684 V DC [20].
After incorporation of the smoother and storage the output of the inverter power and voltage at the PCC is as shown in Figure 13. In comparison to Figure 10, it is shown in Figure 13 that after incorporation of the storage system along with the smoother, the ramp up and ramp down rate of the inverter output power has been smoothed out and the corresponding voltage fluctuation has also been minimized. Table 6 shows the values of voltage at the PCC against the percentage change of inverter output power/min. After introduction of the smoother the voltage swelling and sag at the PCC is reduced significantly. The minimum and maximum voltage fluctuation are increased from 345 V and 484 V revives up to 404 V and 415 V, respectively. Table 7 shows the frequency of occurrence of ramp up In comparison to Figure 10, it is shown in Figure 13 that after incorporation of the storage system along with the smoother, the ramp up and ramp down rate of the inverter output power has been smoothed out and the corresponding voltage fluctuation has also been minimized. Table 6 shows the values of voltage at the PCC against the percentage change of inverter output power/min. After introduction of the smoother the voltage swelling and sag at the PCC is reduced significantly. The minimum and maximum voltage fluctuation are increased from 345 V and 484 V revives up to 404 V and 415 V, respectively. Table 7 shows the frequency of occurrence of ramp up and ramp down of the inverter output power/min beyond the tolerable limit on positive as well as the negative side. In Table 7 the portions marked in blue show the frequency of occurrence of ramp up and ramp down of the inverter output power/min beyond the tolerable limit in the positive as well as the negative direction after introduction of the smoother, calculated as 8.33% throughout the year. Thus from the point of view way it can be considered that the year-round availability of the system is 90%. In comparison to Table 4 it shows that the availability of the system is improved up to 4.86% throughout the year. Similarly a PVsyst report has been generated considering 10% availability factor and the result is presented in Figure 14. and ramp down of the inverter output power/min beyond the tolerable limit on positive as well as the negative side. In Table 7 the portions marked in blue show the frequency of occurrence of ramp up and ramp down of the inverter output power/min beyond the tolerable limit in the positive as well as the negative direction after introduction of the smoother, calculated as 8.33% throughout the year. Thus from the point of view way it can be considered that the year-round availability of the system is 90%. In comparison to Table 4 it shows that the availability of the system is improved up to 4.86% throughout the year. Similarly a PVsyst report has been generated considering 10% availability factor and the result is presented in Figure 14. From Figures 11 and 14 it is seen that the energy production has been increased from 687. 4 MWh/year to 695.7 MWh/year. From Figures 11 and 14 it is seen that the energy production has been increased from 687. 4 MWh/year to 695.7 MWh/year.

Financial Issues
The financial aspect is one of the major concerns of a solar power plant which eventually depends on the loss of energy due to system unavailability. The developer of a solar power plant earns revenue by selling energy to the distribution company against a declared tariff decided by the Central Electricity Regulatory Commission in India.
Thus it would be desirable to generate more energy to earn more revenue which in turn increases the project viability. The unavailability of the system increases the loss of energy, followed by a loss of earned revenue. Table 8 above describes the investment for a 500 kWp solar PV plant without a solar smoother.  [21] Escalation of power tariff = 5% per year Accelerated Depreciation Benefit = 80% in first year Tax implication = 30% Table 9 shows the cumulative revenue from selling power for ten years is INR 65,520,061/-, considering accelerated depreciation benefit and income tax deduction against an investment of INR 27,500,000/-.

Financial Issues
The financial aspect is one of the major concerns of a solar power plant which eventually depends on the loss of energy due to system unavailability. The developer of a solar power plant earns revenue by selling energy to the distribution company against a declared tariff decided by the Central Electricity Regulatory Commission in India.
Thus it would be desirable to generate more energy to earn more revenue which in turn increases the project viability. The unavailability of the system increases the loss of energy, followed by a loss of earned revenue. Table 8 above describes the investment for a 500 kWp solar PV plant without a solar smoother.  [21] Escalation of power tariff = 5% per year Accelerated Depreciation Benefit = 80% in first year Tax implication = 30% Table 9 shows the cumulative revenue from selling power for ten years is INR 65,520,061/-, considering accelerated depreciation benefit and income tax deduction against an investment of INR 27,500,000/-. Now, by retrofitting the solar smoother and the battery bank in the system the material cost has been increased by INR  Now, by retrofitting the solar smoother and the battery bank in the system the material cost has been increased by INR 11.5 Lac, where the cost of the battery bank is INR 10.0 Lac and the cost of the smoother circuits with protections is INR 1.5 Lac. Therefore the cost of the whole system including solar smoother, battery bank and its replacement two times in ten years has become INR 30,650,000/-. Table 10 shows a detailed bill of materials and the specifications of each item.  Table 11, which shows the cash flow after incorporation of the solar smoother and battery bank in the system, indicates that the cumulative revenue has become INR 67,709,122/-from selling power for "ten" years against an investment of INR 30,650,000/-.

Financial Issues
The financial aspect is one of the major concerns of a solar power plant which eventually depends on the loss of energy due to system unavailability. The developer of a solar power plant earns revenue by selling energy to the distribution company against a declared tariff decided by the Central Electricity Regulatory Commission in India.
Thus it would be desirable to generate more energy to earn more revenue which in turn increases the project viability. The unavailability of the system increases the loss of energy, followed by a loss of earned revenue. Table 8 above describes the investment for a 500 kWp solar PV plant without a solar smoother.  [21] Escalation of power tariff = 5% per year Accelerated Depreciation Benefit = 80% in first year Tax implication = 30% Table 9 shows the cumulative revenue from selling power for ten years is INR 65,520,061/-, considering accelerated depreciation benefit and income tax deduction against an investment of INR 27,500,000/-. Now, by retrofitting the solar smoother and the battery bank in the system the material cost has been increased by INR 11.5 Lac, where the cost of the battery bank is INR 10.0 Lac and the cost of the Therefore, the cumulative revenue increases significantly by INR 2,189,061/-after inclusion of the solar smoother and battery bank within the system which is about 20% of the net cumulative realization (INR) in the first year of the 500 kWp system. For larger systems of MW scale the benefit will be substantial.

Conclusions
In this paper, a generalized methodology has been proposed for estimation of the distribution grid dynamics under the impact of variable solar PV generation throughout the year using the PSCAD simulation platform. Further, the non-availability of energy at the grid end is estimated by fitting the PSCAD observations into a PVsyst model as input and this problem is mitigated by incorporating a solar smoother in the existing system. The model performance is validated by a practical 500 kWp solar PV power system. It has been found that the energy non-availability is reduced by 4.86% after inclusion of the solar smoother technology. Further, the economic impact of the proposed system