Combined ANFIS–Wavelet Technique to Improve the Estimation Accuracy of the Power Output of Neighboring PV Systems during Cloud Events

: The short-term variability of photovoltaic (PV) system-generated power due to ambient conditions, such as passing clouds, represents a key challenge for network planners and operators. Such variability can be reduced using a geographical smoothing technique based on installing multiple PV systems over certain locations at distances of meters to kilometers. To accurately estimate the PV system’s generated power during cloud events, a variability reduction index ( VRI ), which is a function of several parameters, should be calculated precisely. In this paper, the Wavelet Transform Technique ( WTT ) along with Adaptive Neuro Fuzzy Inference System (ANFIS) are used to develop new models to estimate the PV system’s power output during cloud events. In this context, irradiance data collected from one PV system along with other parameters, including ambient conditions, were used to develop the proposed models. Ultimately, the models were validated through their application on a 0.7 km 2 PV plant with 16 rooftop PV systems in Brisbane, Australia.


Introduction
Global warming due to greenhouse gas emissions produced by conventional fossil fuel-based power plants has urged all nations to invest in renewable energy, such as photovoltaic (PV) systems and wind energy [1]. According to [2], the generated electricity in Australia is divided into 79% fossil fuel-based generation whereas 21% is generated from renewable energy resources. Renewable energy continues to grow rapidly in Australia, which is planned to reach 50% of the overall domestic electricity production in 2024 and 100% in the year 2032 [3]. Although there are several advantages in using PV systems, including easy implementation in both residential and commercial sites along with their low cost and reliable operation, it still comprises some issues owing to its intermittent nature. Due to the ambient conditions, PV-generated power is not constant and it varies with short-(cloud movements) and long-term (irradiance daily and seasonal) variations. Long-term variability can be easily planned for based on the well-known daily and seasonal solar irradiance in various geographical areas. On the other hand, short-term irradiance fluctuation due to the unexpected passing of clouds of unknown size, direction and speed, may result in frequency instability due to the reduction in the generated power [4][5][6]. Estimating such power is one of the main challenges faced by network operators to prevent frequency instability due to such short-term power reduction, which is why many adopt a suitable backup energy storage system. According to utilities guidelines, the PV system's power output Geographic smoothing or the variability reduction index (VRI) is discussed in several studies in the literature. In such studies, different methods and metrics are presented to estimate and analyze the short-term variability in the PV system's power output. For instance, the correlation coefficient and standard deviation are employed in [9][10][11][12][13][14] while [15][16][17] utilized the variability score and variability index as mathematical metrics to calculate the power output variability of PV systems. Other studies introduced different metrics to investigate the power output variability, such as the Multifractal Detrended Fluctuation Analysis (MFDFA), which is a time series modeling method especially used for extreme ambient events [18]. Similarly, [19] adopted the Fourier analysis method to characterize the variability using power spectral density. In such studies, data used to calculate the PV system's power variability are obtained using generated synthetic PV system data based on real measurements as presented in [1,20]. Few studies have been introduced to estimate the power output variability using the Wavelet Transform Technique based on measured data of only one PV sensor and by considering the PV system's locations, cloud speed and irradiance time series [21][22][23]. Table  1 summarizes various methods presented in the literature along with the key finding of each study.
Although several models have been introduced to estimate the power output for a group of PV systems distributed within meters to kilometers, the estimation error calculated so far ranges between 8% to 20%, which reveals the relatively poor accuracy of such techniques [23]. As such, developing a new model with a higher estimation accuracy would be a significant contribution to this area of research. Such a high estimation accuracy model can be employed in power networks of high PV system penetration for better reliability, control and asset management, which is the main aim of this paper. In this paper, a combined WTT and Adaptive Neuro Fuzzy Inference System (ANFIS)-based technique with a higher estimation accuracy than existing models are presented. The proposed technique is employed to estimate the power output of a group of rooftop PV systems distributed over 0.7 km 2 area during variability events. The robustness of the proposed technique is validated through its application on real data collected from several PV systems distributed within a few hundred meters to 1 km distance.
In the remaining sections of this paper, brief reviews on the VRI, WTT and ANFIS are presented in Section 2. The proposed model methodology is presented in Section 3. Results and discussions are Geographic smoothing or the variability reduction index (VRI) is discussed in several studies in the literature. In such studies, different methods and metrics are presented to estimate and analyze the short-term variability in the PV system's power output. For instance, the correlation coefficient and standard deviation are employed in [9][10][11][12][13][14] while [15][16][17] utilized the variability score and variability index as mathematical metrics to calculate the power output variability of PV systems. Other studies introduced different metrics to investigate the power output variability, such as the Multifractal Detrended Fluctuation Analysis (MFDFA), which is a time series modeling method especially used for extreme ambient events [18]. Similarly, [19] adopted the Fourier analysis method to characterize the variability using power spectral density. In such studies, data used to calculate the PV system's power variability are obtained using generated synthetic PV system data based on real measurements as presented in [1,20]. Few studies have been introduced to estimate the power output variability using the Wavelet Transform Technique based on measured data of only one PV sensor and by considering the PV system's locations, cloud speed and irradiance time series [21][22][23]. Table 1 summarizes various methods presented in the literature along with the key finding of each study.
Although several models have been introduced to estimate the power output for a group of PV systems distributed within meters to kilometers, the estimation error calculated so far ranges between 8% to 20%, which reveals the relatively poor accuracy of such techniques [23]. As such, developing a new model with a higher estimation accuracy would be a significant contribution to this area of research. Such a high estimation accuracy model can be employed in power networks of high PV system penetration for better reliability, control and asset management, which is the main aim of this paper. In this paper, a combined WTT and Adaptive Neuro Fuzzy Inference System (ANFIS)-based technique with a higher estimation accuracy than existing models are presented. The proposed technique is employed to estimate the power output of a group of rooftop PV systems distributed over 0.7 km 2 area during variability events. The robustness of the proposed technique is validated through its application on real data collected from several PV systems distributed within a few hundred meters to 1 km distance.
In the remaining sections of this paper, brief reviews on the VRI, WTT and ANFIS are presented in Section 2. The proposed model methodology is presented in Section 3. Results and discussions are presented in Section 4. Section 5 presents a sensitivity analysis of the parameters used in the developed model followed by the main conclusions.

Variability and Correlation Coefficient
The VRI represents the ratio of the power variance of one PV system to the entire PV plant at a specific time scale [23] and is mathematically expressed as The correlation coefficient in Equation (1) is a measure of the relationship of two random variables [24]. A number of correlation models are presented in the literature [23,[25][26][27][28] of which Pearson's is considered as one of the popular correlation models [25]. Reference [26] presented an empirical correlation coefficient (ρ Hf m,n , indicated by the Hoff model in the course of this study) that considers the timescale (T) and the distance (d m,n ) between two PV systems' location as below: Perez [27] and Lave [23] (ρ Pz m,n and ρ Lv m,n refer to the Perez and Lave models in the rest of this paper) have also introduced empirical exponential models as in Equations (3) and (4), respectively.
Reference [28] presented another mathematical model (referred to ACM model) by considering the cloud direction and size as in Equation (5).
where b 1 to b 2 are normalized parameters considering V clo , d m,n and cloud direction.
In addition, Reference [29] presented a VRI technique (referred to the VRI-GEP model in the rest of this study) by using the maximum values of the existing VRI models [7] and a genetic algorithm technique called Gene Expression Programming (GEP).  Standard deviation and probability density function.
variability score and cumulative distribution function New model to calculate the high frequencies of power variability [16] 2 PV systems, Canada, with 41 sensors, distance of 400 km and 10-ms data time resolution.
Employed variability score and variability index.
Classified daily irradiance into 4 classes: clear, overcast, low and high frequencies solar variability. [17] 10 PV systems in USA, with hundreds of km distance and 30-s time resolution Variability score method Variability score is a useful estimation metric. However, it is a poor predictor for PV tracking systems. [18] 20 PV systems in Hawaii, 1 PV system in Spain, 1 PV system in Algeria, 1 PV system in Germany.
Correlation and the multifractal detrended fluctuation analysis. The study focused on temporal variability.
Variability is not stationary in which it was the highest in Hawaii and the lowest in Algeria. [19] 50 PV systems, India, 400 km distance, 1-min resolution.
Fourier transform analysis and correlation coefficient.
Good smoothing in short-time events. Improvement for more added PV systems. [1,20] 3 PV systems, Hong Kong, 14 to 18km distance and 1-min resolution.
Using system advisor model An improved model to generate spatial synthetic PV data [21] 99 PV systems, Germany, 10 to 12 km distance, 1-s. resolution.

Discrete wavelet transform model
Distinguished three conditions, clear, partly cloudy and overcast. [22] One PV plant, USA, 60 km. 1-s resolution Discrete wavelet technique Variability is higher at the end of the feeder. [23] 550 PV systems, Japan and Nevada, 1-6 km 2 respectively, 1-s resolution. Discrete wavelet transform model and variability reduction index.
Estimation of the power output with errors of 8% to 20%.
Thus, the VRI in Equation (1) can be calculated using any of the correlation coefficients presented in Equations (2) through (5) as below: where x refers to one of these four models.

Proposed Technique
The proposed technique is a combination of the WTT and the Adaptive Neuro-Fuzzy Inference System (ANFIS). Wavelet Transform is a signal processing technique that transforms periodic and Energies 2020, 13, 1613 5 of 15 non-periodic time-series signals into the time-frequency domain [7]. Continuous Wavelet Transform (CWT) is used for a continuous time-series signals and is expressed as Because the irradiance time-series signal is discrete [23], Discrete Wavelet Transform (DWT) as in Equation (8) is more suitable for the analysis in this paper.
where y[k] is the discrete time-series signal of the PV power output, and a and c are scale and shift factors; respectively. The mother wavelet function could be one of Harr, Mexican Hat, Symlets, Coiflets, Biorthogonal, Meyer, Morlet and Daubechies [7]. The decomposed signals resulted from the DWT process can be divided into high-frequency (HF) and low-frequency (LF) modes. In this paper, a new technique based on ANFIS is proposed to capture the DWT features and model them [30]. Traditional mathematical modeling methods are inappropriate tools for dealing with stochastic and uncertain data, such as solar irradiance. An example of a more accurate modeling approach based on the integration of Artificial Neural Network (ANN) and Fuzzy Control System (FCS) is ANFIS, considered a powerful technique that can map the relationship between a set of input and output variables [31]. ANFIS combines the features of ANN and FCS and can function in a noisy, uncertain and indefinite environment. In this technique, the ANN provides the FCS with learning capabilities, while the FCS provides ANN with an organized structure with well-established fuzzy rules, reasoning and thinking [32]. The two learning algorithms adopted in the ANFIS technique are the back-propagation and hybrid approaches that act to reduce the error between real and estimated data [33].
To build a model that aims at improving the estimation accuracy of the power output of several rooftop PV systems using ANFIS, it is important to identify the input parameters that influence the power output variability, such as ambient temperature (T em ), humidity (H u ), PV plant area (A pv ), the PV system site's power capacity (P cap ), cloud speed (V clo ) as well as the solar irradiance (I POA ). Some of these parameters are collected using a PV pyranometer located close to the PV system while other parameters are obtained based on the PV system's properties. In the developed model, the input data are divided into three categories as below.

1.
PV system properties data: The developed model considers the PV areas (A pv ), which is divided into 6 sub-areas with specific power capacity (P cap ) for each area ((A PV sensor data: This includes T em and, H u that can be directly collected from the pyranometer sensor installed in the PV system site.
When the pyranometer is installed in the same tilt and orientation as the PV system, the output irradiance of the PV sensor will be a plane of the array (I POA ) [23]. Let I POA [k] denote the output of the pyranometer, in which k is a time index.
To improve the ANFIS accuracy to estimate the power output, the time series (I POA ) for each day is divided into five sub-time series. This procedure enables the proposed technique to correlate the Energies 2020, 13, 1613 6 of 15 input and output variables at a specific period of the day more precisely due to the dynamic change of weather conditions over the day. For each sub-time series, the module power output can be calculated by [7]: The average of the total power output of the PV plant (P mdu−out

Eq
[k]) is calculated by assuming each module is providing the same share of the overall power. Then, P mdu−in are normalized by dividing each by P clr POA [k] (power output in a clear sky condition), which can be calculated by replacing I POA [k] in Equation (9) by I POA clr [7], The parameters in Equation (10) are calculated based on the data at the CAMS McClear website [35]. Then, Equations (11) and (12) can be derived: Using Equation (8) ) and low-frequency (P LF−out mode−j [k]) modes. It is to be noted that from the inverse wavelet transform, the first and second low-frequency modes can be ignored as they have no impact on the inverse wavelet process as explained in Equation (13) below [36].
along with the parameters mentioned in points 1 and 2 above will be the ANFIS input parameters, whereas P HF−out mode−j [k] and P LF−out mode− j [k] are the outputs, as shown in Table 2. In Table 2, the first and second columns represent the input and outputs of the ANFIS technique while the third column represents the ANFIS models after the training stages. Table 2. Adaptive Neuro Fuzzy Inference System (ANFIS) input and output data and models after the training stage.

Input Parameters
Output ANFIS Models after Training The results of the ANFIS modeling process comprise four ANFIS models, each one represents the ANFIS-DWT model at a certain wavelet mode. For example, P HF−ANFIS mode−1 [k], P HF−ANFIS mode−2 [k] and P HF−ANFIS mode−3 [k] represent the equivalent wavelet analysis of mode-1, mode-2 and mode-3 for the high-frequency signals, respectively, while P LF−ANFIS mode−3 [k] represents the equivalent wavelet analysis of mode-3 for a low-pass signal as shown in Table 2. Then, by applying the inverse DWT as in Equation (13), the estimation of the equivalent normalized power output of several rooftop PV systems distributed within a certain area is evaluated.
By multiplying the results of Equation (13) by P POA clr [k] and the number of PV modules (N mod ), the equivalent power output seen by the entire PV plant can be estimated as follows: Figure 2 illustrates the flowchart of the main steps of the ANFIS-DWT technique proposed in this paper.
The varied weather over the four seasons for Brisbane city and the PV systems being distributed 202 around the University of Queensland makes it a suitable site for this study. PV sensor data were 203 collected from the PV located at Site 7 in Figure 3. 204

Simulation Results 205
The collected data of the 445 days comprise small, middle and large variability.  Because variability is the main concern by power operators and planners, the variability power index (V pi ) has been introduced as a mathematical tool to show the amount of the power fluctuation at each DWT mode. First, P Eq−ANFIS nr−out is analyzed into six timescales, namely 2, 4, 8, 16, 32 and 64 min, using the DWT method. To calculate the power at each timescale, the wavelet periodogram has been used as below [37]: then, According to [37], the variability reduction index (VRI) can be by calculated using Equation (17):

Performance Evaluation of the Proposed Technique
Data of 445 days in the years 2016 and 2017 with different amounts of variabilities, where days with a completely clear sky were eliminated, were used to develop the proposed model. Data were collected from PV rooftop systems with a 2.14 MW power capacity that were distributed within 0.7 km 2 at the University of Queensland, Brisbane, Australia. The PV site is located at 153 • 29 44" east longitude and 27 • 00 30" south latitude. All PV systems faced the north at an inclination angle of 30 • to produce Energies 2020, 13, 1613 8 of 15 maximum possible power. Brisbane city weather is described as humid and warm all through the year. The temperature in summer and winter ranges between 21 to 35 • C, and 11 to 21 • C, respectively, while it is 15 to 25 • C in autumn and spring. In summer, Brisbane's sky exhibits several passing clouds events, more than in winter [38].
The varied weather over the four seasons for Brisbane city and the PV systems being distributed around the University of Queensland makes it a suitable site for this study. PV sensor data were collected from the PV located at Site 7 in Figure 3.

28
Following the methodology steps in Figure 2 and by using the input and output data in Table 2,  29 the ANFIS model can be built and trained. Since the back-propagation algorithm is relatively slow 30 [36], the hybrid-learning algorithm was adopted to train the proposed ANFIS model. 31 Figure 4 shows the error in estimating the power output when different MFs are employed. 32 Results show that the gaussmf membership function is more accurate than other types after 100 33 epochs training modes and hence it was employed in the ANFIS model. at each sub-time series. Four ANFIS models were developed as shown in Table 2, third column. 38

Simulation Results
The collected data of the 445 days comprise small, middle and large variability. For ANFIS training and testing data, 330 days were used and the rest of the data were used in the validation stage. To improve the ANFIS performance during the training stage, each day was divided into five sub-time series, with the first sub-time series assumed to start at 7:00 a.m. and ending at 8:59 a.m.; the second sub-time series was from 9:00 a.m. to 10:59 a.m.; the third sub-time series was from 11:00 a.m. to 12:59 p.m.; the fourth sub-time series was from 1:00 p.m. to 2:59 pm; and the fifth sub-time series was from 3:00 p.m. to 5:00 p.m.. The sunrise and sunset times have been taken into account in this division. The PV plant was also divided into six small areas, namely the whole PV plant (A Following the methodology steps in Figure 2 and by using the input and output data in Table 2, the ANFIS model can be built and trained. Since the back-propagation algorithm is relatively slow [36], the hybrid-learning algorithm was adopted to train the proposed ANFIS model. Figure 4 shows the error in estimating the power output when different MFs are employed. Results show that the gaussmf membership function is more accurate than other types after 100 epochs training modes and hence it was employed in the ANFIS model. Figure 5 shows the Root Mean Square Error (RMSE) for the last epoch at each ANFIS model and at each sub-time series. Four ANFIS models were developed as shown in Table 2, third column. Following the methodology steps in Figure 2 and by using the input and output data in Table 2,  29 the ANFIS model can be built and trained. Since the back-propagation algorithm is relatively slow 30 [36], the hybrid-learning algorithm was adopted to train the proposed ANFIS model. 31 Figure 4 shows the error in estimating the power output when different MFs are employed. 32 Results show that the gaussmf membership function is more accurate than other types after 100 33 epochs training modes and hence it was employed in the ANFIS model. at each sub-time series. Four ANFIS models were developed as shown in Table 2, third column. 38  Results show that the gaussmf membership function is more accurate than other types after 100 33 epochs training modes and hence it was employed in the ANFIS model. 34 at each sub-time series. Four ANFIS models were developed as shown in Table 2, third column. 38   (14), the estimated PV system power power output can be obtained. Figure 6 depicts the error in the 43 estimated power of the developed model for 115 days using Mean Absolute Error (MAE) and RMSE 44 in comparison with other methods mentioned in Section 2, such as the Hoff, Perez, Lave and VRI-45 GEP models. Table 3 summaries the results of Figure 6  To calculate the VRI parameter, the power output was analyzed from six wavelet modes using 56 Equation (8), and then the wavelet periodograms factor in Equation (15) was employed at each 57 wavelet mode. Using Equation (16), the amount of power involved in the variability ( ) at 58 each mode and each sub-time series can be evaluated. Figure 7 shows an example for one day, 59 December 3, 2017, that was chosen to explain this step. The blue line in Figure 7 represents 60 ( ) for one PV system (No. 7), while the black and red lines represent ( ) for the 61 actual and predicted power output for 16 PV systems, respectively. From these plots, the distance 62 between the blue line (one PV system) and black and red lines (16 PV systems) characterizes the 63 geographic smoothing or the VRI value. The can be calculated using Equation (17)

Discussion 71
The estimation improvement strategy adopted in this study is based on dividing the time series 72 data into several sub-time series (five sub-times series in this study) along with employing the 73 Wavelet Transform Technique. ANFIS is employed to accurately map the relationship between the 74 input and output data set. 75 Dividing the time series data into five segments helps the ANFIS algorithm map the correlation 76 between the input and output data rapidly and precisely due to the fact that, within each selected 77 period, the environmental conditions such as tempreture, zenith angle and power output would not 78  (13) and (14), the estimated PV system power power output can be obtained. Figure 6 depicts the error in the estimated power of the developed model for 115 days using Mean Absolute Error (MAE) and RMSE in comparison with other methods mentioned in Section 2, such as the Hoff, Perez, Lave and VRI-GEP models. Table 3 summaries the results of Figure 6 by using the average values of the MEA and RMSE for each method mentioned in the introduction section along with the proposed technique in this paper. This table also illustrates the minimum and maximum improvement values attained by the proposed technique when compared with the existing models in the literature. It can be seen from these results that the proposed method is more accurate than other models published in the literature. Considering the MAE, the minimum improvements attained by the proposed technique for the five sub-time series are 11.44, 8.87, 7.46, 9.03 and 10.41, respectively, whereas the maximum improvements for the same sub-time series are 36.89, 28.6, 24.06, 28.81 and 33.07, respectively. Similarly, results reveal a very similar tendency in the enhancement attained by the proposed technique when the RMSE metric is considered. To calculate the VRI parameter, the power output was analyzed from six wavelet modes using Equation (8), and then the wavelet periodograms factor in Equation (15) was employed at each wavelet mode. Using Equation (16), the amount of power involved in the variability V pi T mode−i at each mode and each sub-time series can be evaluated. Figure 7 shows an example for one day, December 3, 2017, that was chosen to explain this step. The blue line in Figure 7 represents V pi T mode−i for one PV system (No. 7), while the black and red lines represent V pi T mode−i for the actual and predicted power output for 16 PV systems, respectively. From these plots, the distance between the blue line (one PV system) and black and red lines (16 PV systems) characterizes the geographic smoothing or the VRI value. The VRI can be calculated using Equation (17) as shown in Figure 8. These figures attest that VRI can be changed during the day. For example, the value of the VRI at Mode 1 is 8.2 at sub-time series 7:00-8:59 a.m., while it is 6.2, 7, 7, 7.5 at sub-time series 9:00-10:59 a.m., 11:00-12:59 p.m., 1:00-2:59 p.m. and 3:00-5:00 p.m., respectively. This is attributed to the change in the ambient parameters affecting the VRI, such as V clo , T em , and H u during the daytime. series into low-and high-frequency bands reduces the ANFIS confusion during the training process 80 to map the relationship between the input and output data. 81 Results of the proposed technique in this paper reveal its high accuracy when compared with 82 other estimation models published in the literature. Results in Figure 6 and Table 3 also show a 83 comparison of the model proposed in this paper and our recent VRI-GEP technique [29]. It can be 84 observed that both the GEP and ANFIS techniques provide more accurate results than other 85 techniques in the literature [39,40]. Reference [29]

Discussion
The estimation improvement strategy adopted in this study is based on dividing the time series data into several sub-time series (five sub-times series in this study) along with employing the Wavelet Transform Technique. ANFIS is employed to accurately map the relationship between the input and output data set.
Dividing the time series data into five segments helps the ANFIS algorithm map the correlation between the input and output data rapidly and precisely due to the fact that, within each selected period, the environmental conditions such as tempreture, zenith angle and power output would not change much. Likewise, the division conducted by the Wavelet Transform method on each sub-time series into low-and high-frequency bands reduces the ANFIS confusion during the training process to map the relationship between the input and output data.
Results of the proposed technique in this paper reveal its high accuracy when compared with other estimation models published in the literature. Results in Figure 6 and Table 3 also show a comparison of the model proposed in this paper and our recent VRI-GEP technique [29]. It can be observed that both the GEP and ANFIS techniques provide more accurate results than other techniques in the literature [39,40]. Reference [29] employed the GEP technique to develop a VRI model whereas in this paper, the ANFIS-DWT model is used to estimate the power output directly without calculating the VRI as per the strategy of all existing models in the literature. This reduces the calculation complexity and improves the estimation accuracy. Results show that the ANFIS model proposed in this paper attains a minimum and maximum improvement of 6.44% and 11.44% when compared with the VRI-GEP model by considering the MAE metric.

Sensitivity Analysis
Following the validation of the proposed model's accuracy in the above section, this section is aimed at investigating the sensitivity of the parameters used in the proposed model. In this regard, each parameter was changed in an expected range as listed in Table 4 while maintaining all other parameters unchanged [41]. The effect of such a change on the power output was analyzed using the standard deviation (SD) metric. To standardize the model sensitivity analysis, SD is divided by the average value of the measured parameter. Figure 9 shows the impact of each parameter at different sub-time series. It can be seen that the input power parameter has a more significant impact than the ambient temperature while the humidity, cloud speed and the ratio of the PV plant capacity to the plant area have approximately the same impact. Table 4 illustrates the minimum and maximum values for each parameter that were considered during the sensitivity analysis. average value of the measured parameter. Figure 9 shows the impact of each parameter at different 104 sub-time series. It can be seen that the input power parameter has a more significant impact than the 105 ambient temperature while the humidity, cloud speed and the ratio of the PV plant capacity to the 106 plant area have approximately the same impact. Table 4 illustrates the minimum and maximum 107 values for each parameter that were considered during the sensitivity analysis. 108

Conclusions
Using Wavelet-based and ANFIS modeling approaches, new PV system power output estimation models are presented in this paper. The developed models are employed to estimate the generated power of a group of PV systems distributed over 0.7 km 2 using irradiance data of one PV sensor located near the PV plant. Parameters such as PV area, capacity, cloud speed, temperature and humidity are considered as the inputs to the proposed model. Real data collected from a PV plant in Brisbane, Australia, were used to train the proposed model. Results reveal the high accuracy of the developed model when compared with the existing estimation models in the literature. Sensitivity analysis of each parameter reveals that the power output collected from the PV sensor has the largest impact on estimating the entire PV system's power output using the developed model. Ambient parameters, including humidity, cloud speed and the power capacity of the overall PV plant, have similar but less impactful effects on the estimated power.