Experimental Evaluation of Distortion Effect for Grid-Connected PV Systems with Reference to Different Types of Electric Power Quantities

: Aware of the fact that the installed PV capacity and its power production rapidly increased in the last decade, with the huge impact that has been done to the power system, the distortion effects for grid-connected PV systems with reference to different types of electric power quantities will be presented in this article. The impact of the frequent ﬂuctuation of solar irradiance on the behavior of the grid-connected PV system, due to cloud movements and resulting shadows and in terms of power quality and the evaluation of power components, is the topic of analysis in this research. Besides the simulation results of certain study cases, an experimental evaluation of electric power quantities on an actual PV system in real weather conditions was also performed. The experimental setup, formed through the combination of a PC and multifunctional I/O board with an appropriate software solution, was established and used for obtaining the target results. The methodology used for the evaluation of electric power quantities relied on the current physical components (CPC) theory for power deﬁnition. The experimental results were obtained for three different cases, namely, the low, medium, and high solar irradiance cases. On the basis of these results, the conclusions about distortion effects are given.


Introduction
In the last decade, the installation and production of photovoltaic (PV) systems have been increasing every year. Looking at the global trends, the cumulative installed PV capacity in 2020 was 707 GW, while electricity generation reached 844 TWh ( Figure 1) [1]. The numbers, increased ten-fold, show that a significant influence of PV systems has been made on power systems. Photovoltaic technology is an attractive technology of renewable power generation and meets the new criteria of clean energy and sustainability. The conducted research and technology improvement have led to rapid progress in increasing the efficiency of photovoltaic modules and their optimization on the production of energy through better utilization of incident solar irradiance. Environmental conditions, as well as system design, have proven to be the most important factors in the operation of PV systems, and these can have a significant impact on the efficiency and power quality response of the whole system. The variable power flow, due to the fluctuation of solar irradiance (cloud movements and resulting shadows), temperature, and the topology of power converters are some of the parameters that affect the power quality of photovoltaic systems. These fluctuations can cause adverse effects, such as over/under-voltage issues, reverse power flow, and power quality issues [2]. Different mitigation methods have been developed and proposed for overcoming this issue, such as the local control of reactive power, use of volt-var control strategies, more efficient inverters, additional storage units, and the implementation of predictive control techniques. Through a detailed analysis of the effect of solar irradiance on the power quality behavior of a grid-connected PV system, the authors signified in [3] that low solar irradiance can significantly affect the output of a PV system, maintaining the power factor at a low level due to comparable production of active and reactive power. In addition, the influence on distribution grids of large-scale as well as on rooftop PV systems can be found in references [4,5].
The research conducted in [6,7] aims to characterize the power quality parameters of the grid-connected PV inverter under different operating modes, as well as focuses on the assessment and management of the power quality issues, to be fulfilled at the point of common coupling, as imposed by the new standards. According to the technical brochure [8], the harmonic emission of PV inverters can be divided into primary and secondary emissions, where the primary refers to the emission determined by the control algorithm implemented in the PV inverter, and the secondary is mainly caused by the characteristics of the grid side filter circuit of the PV inverter.
According to [9,10], the researchers' tasks are mainly to investigate a novel strategy for power converter control in order to improve grid power quality. In [9], it is proposed that control of the Pulse-Width Modulation (PWM) inverter is a way to balance between active power production and power quality improvement without exceeding the whole system power capacity. The priority is set so that the first line preference has active power production, then reactive power compensation, and finally active filtering. Basic pq theory (instantaneous active and reactive power components) is used to accomplish this task. Unlike previous research, ref. [10] developed an algorithm for the control structure of a PV inverter that operates as a shunt active power filter for supplying unbalanced and nonlinear loads. In addition, in very new research [11], power quality improvement in the hybrid renewable grid-connected system has become a promising research area. An intelligent control strategy is set up as an optimization solution for optimum power quality enhancement.
On the other hand, in the systems with harmonic distortions, voltage, and current asymmetries (i.e., nonlinear systems), the issue of non-active power identification is a topic that is still being debated in the wider academic community. In addition to this fact, the most convenient literature for power definition in the presence of harmonics is the IEEE Standard 1459-2010 [12]. This standard gives many definitions of power in three-phase systems, namely, of sinusoidal, balanced voltages, and currents; sinusoidal, unbalanced voltages, and currents; and nonsinusoidal, unbalanced voltages, and currents. As the grid-connected PV systems produce electric power with a significant share of current harmonics, especially in intervals with a low level of solar irradiance, these systems can be classified as nonlinear systems. As it is not uncommon that certain voltage asymmetries can appear in common coupling, the evaluation of electric power quantities is an additional task to understand.
In this contribution, the main set task is the evaluation of the electric power quantities of the grid-connected PV system depending on different levels of irradiance, which is not discussed in the literature so far. For an actual PV system, in real weather conditions, the weather condition dependence of electric power quantities will be presented, on the basis of the measurement values of voltages and currents for different cases, namely, low, medium, and high irradiance, and by applying the current physical components (CPC) theory for power definition [13][14][15][16]. In addition to the experimental results, the results obtained by a simulation model (developed for analysis and better understanding of the influence of the environmental and control parameters of the PV inverter on the distribution grid) are also given.

Used Methodology
The major focus of this research, dealing with the impact of grid-connected PV systems, is based on power quality analysis and harmonic mitigation techniques [17][18][19]. In this article, the focus is on the evaluation of electric power quantities based on different solar irradiance levels. For that purpose, the first step was the creation of a simulation model of a real 15.9 kWp residential grid-connected PV system. The system behavior was considered based on the analysis of the current harmonic injection into the grid and based on the selected topology and control structure.
An experimental measurement was performed and analyzed for three different cases (i.e., irradiance levels), followed by the different active power production, THD I factor, power factor (PF), and the displacement power factor (DPF). These quantities were obtained based on power quality monitoring using a power quality analyzer Fluke 435. The considered cases included low, medium, and high solar irradiance levels, as shown in Table 1, in order for the effects of the different irradiance levels on power quality to be observed. Measurement sampling was performed with a frequency of 50 µs; therefore, the measuring voltage and current values require large hardware resources. In addition, a numerical calculation of power quantities required a lot of time. Due to these reasons, the experimental results were presented with a 3 s measurement period. The CPC theory was selected for the evaluation of the electrical power quantities. This is because it provides a complete physical interpretation of power phenomena in electrical systems. Pursuant to the CPC theory for a three-phase nonsinusoidal and unbalanced system, the current decomposition is conducted by active i a , scattered i s , reactive i r , and unbalanced i u current; Equations (1)-(4).
where G e refers to the equivalent conductance with respect to the active power P, G k stands as the equivalent conductance at the frequency of the k th harmonic (different from G e ), and B k /Y k is the equivalent susceptance/admittance at the frequency of the k th harmonic, respectively [14,15]. Appropriate power quantities can be calculated by using the root mean square (RMS) values of the above-mentioned current components.
All these electric power quantities satisfy the orthogonality Equation (8): The results of the evaluation of PV active power production (P), total harmonic current distortion (THD I ), power factor (PF), and displacement power factor (DPF) for different irradiance levels during three consecutive days are presented in Figure 2.

Simulation Model
A simulation of a real system process during the operation time was created. The simulation model of the grid-connected PV system was developed with the aim of testing and for a better understanding of the influence of environmental and control parameters on the distribution grid. In addition, the simulation results, as well as the experimental results (in real-time) might lead to better and more clarified research conclusions. For the observed grid-connected PV system of 15.9 kWp, a simulation model was created in MATLAB Simulink environmental. The simulation model of a grid-connected inverter as a double-stage three-phase topology was implemented. In the first stage, the maximum power of the PV array was obtained by using an MPPT (Maximum Power Point Tracking) converter designed as a Boost converter, on which the P&O algorithm for maximum powerpoint tracking was implemented. The second stage transformed DC to AC power through the three-phase inverter and reduced the content of higher harmonics from the inverter output with the grid filter (usually LCL). The control structure was based on well-known dq-reference frame technique [20][21][22], the block structure of which is shown in Figure 3. The PV inverters, connected to the distribution grid, were mostly set to produce only active power without reactive power control capability, so that the situation is modeled in this simulation.

Simulation Results
As mentioned previously, the control algorithms of a PV inverter cause the occurrence of high-order harmonics through switching operations, subsequently occurring as a current ripple. The level of the current ripple can be controlled by a grid-side filter, but this increases the losses in the inverter, i.e., reduces its efficiency. Modern PV inverters are implemented as double-stage converters. In the first stage, the MPPT function is implemented, and in the second stage, DC-AC conversion is performed. The MPPT function performs a constant change in the DC voltage level, which causes a change in the transfer function of the inverter and affects the generation of harmonic distortions of the current at the output terminals of the inverter. Previous studies demonstrate that the harmonic distortion of the inverter current mainly depends on the inverter load (irradiance level on PV array). In general, at a lower irradiance level, which corresponds to the lower PV power injection, the THD I factor is higher (Figure 2). As a result of the simulation work, there is a presentation of grid voltage and current waveforms for different irradiance levels (Figure 4a) as well as its corresponding harmonic spectra in the output current (for example, in Figure 4b, the harmonic spectra of the one-phase current are given). It should be noted that waveforms in Figure 4a are not obtained continuously by changing the irradiance level. They are managed so that they can respond to different irradiance levels in a stationary mode. In this research, the emphasis has not been on the dynamic response of the inverter and its impact on harmonic behavior, but on the level of the THD factor, which implies a significant content of electric power quantities, especially the non-active parts.
Observing the harmonic spectra in Figure 4b, the low irradiance (<10%) cause has a significant share of high-order harmonics in the output current. However, the THD I factor is calculated in relation to the fundamental harmonic, so that the value of fundamental harmonic will be lower within the lower irradiance, which may implicate the main reason for a drastic increase in the THD I . In addition, the THD U value is important at the point of common coupling; therefore, this value will be low in the case of small or medium PV systems connected to the strong grid, besides the current harmonic injection.
By analyzing the PV power profiles (in Figure 2) in different weather conditions, it can be concluded that this significant share of high-order harmonics may take a certain period, which justifies the task of evaluation and analysis of electric power quantities.

System Configuration
The observed PV system is located on the rooftop of the building of the Faculty of Electrical Engineering in East Sarajevo (43.8237 • N, 18.3743 • S). The installed 15.9 kWp gridconnected system consists of 60 polycrystal modules with 265 Wp each. These modules are connected in strings of 20 modules each, with two of these three strings in parallel. The Fronius SYMO 15.0-3-M inverter is used as a link between the PV array and grid, with a rated output power of 15 kVA. The input voltage range is 200-1000 VDC with a maximum efficiency occurrence at ≈600 VDC of input voltage. This system is connected to the local facility distribution transformer, and it is used for local load supply, handing over the rest of the produced energy into the grid ( Figure 5).

Measurement and Acquisition System
The setup for the experimental part of this contribution is implemented as a measurement and parameter acquisition system consisting of fully isolated and multichannel voltage and current sensing modules [23], as well as a multifunctional input/output card (Humusoft MF 624 [24]). The measurement modules are designed to be used in a wide range of applications, such as for three-phase systems monitoring, power converters, motor drives, and PV monitoring applications. The voltage sensor ranges from 10x up to 100 V, and from 100x up to 1000 V, while the current sensing goes up to 100 A. The applied MF 624 multifunction I/O card is generally used in the process of control and measurement applications. The board offers A/D conversion, D/A conversion, digital input, digital output, quadrature encoder input, timers, counters, PWM, pulse, and frequency measurements. The boards contain 8 14-bit analog input channels with simultaneous samples, hold circuits, very short conversion time, 8 14-bit analog output channels, 8 digital inputs, 8 digital outputs, 4 encoder inputs with differential line receivers, and 4 counters or timers. The card is designed for standard data acquisition and control applications and optimized for use with Real Time Toolbox for Simulink ® [25]. A block diagram of used components for measurement and acquisition setup is presented in Figure 6.

Experimental Results
By analyzing the three cases defined in Table 1, the obtained results are shown in Figures 7-9. The evaluated power components, with the power factor, voltage and current waveforms, and harmonic spectra with THD I and THD U factors are presented. The evaluated values of the THD factors, electric power quantities, and the power factor are also given in Table 2.    For the low irradiance case, when the power production is lower than 10%, the power factor is significantly lower than 1 (primarily as a consequence of current distortion), whereby the non-active power components are comparable to the value of active power. The impact of these components is also visible in the current waveforms, resulting in a high THD Ia factor of 61.14% (THD Ib = 41.47% and THD Ic = 42.56% are lower than THD Ia , but still relatively high, which is expected when looking at the waveform). According to the results in Figure 2, with a longer low irradiance period during the day, the distortion effect can be a significant issue.
Unlike the low irradiance case, the cases with high and medium solar irradiance do not have a significant impact concerning the distortion effect. The voltage and current waveforms are practically in an ideal sinusoidal shape; therefore, the THD I factors are in the low range (for example, THD Ia is 2.14% in the high irradiance case, and 2.61% in the medium irradiance case). In addition, in these cases, all power components or electric power quantities can be evaluated, but their impact on total power is almost negligible and therefore the power factor almost reaches the value 1.
Observing the results from Table 2, it is found that the share of the non-active power components in the apparent power depends on the irradiance level. The distortion power D s and the unbalanced power D u are lower in the case with lower irradiance, but their share in the apparent power is significantly higher than the share of the same components in the other two cases.
On the other hand, the impact of harmonic injection to the grid voltage did not come into question because the short circuit power, at the point of common coupling, is much higher than the installed power of the connected PV system. In all three cases, the THD U factor remained less than 2%.
In addition, from Figures 7-9, it is observed that a certain share of fifth and seventh harmonics from the grid voltage are driven to the current's spectra. This phenomenon is justified by the existence of a Phase-Locked Loop (PLL) structure in the control structure of a grid-connected inverter. Therefore, it can be concluded that this share of fifth and seventh harmonics in the current spectra is most likely a consequence of their existence in the grid voltage.

Comparison of the Simulation and Experiment Results
As mentioned previously, the emphasis in this research is on the level of the THD factors, which imply significant content of the non-active power. By analyzing the THD I factor in different weather conditions, it can be concluded that this significant share of high-order harmonics may cause different levels of electric power quantities. Different mitigation techniques and the optimal design of a grid-side filter, with improved control algorithms, could suspense and reduce the share of harmonic injection. The developed simulation model of a real PV system (although with a lack of manufacturer information for the grid-connected inverter) provides a better understanding of the whole process and a comparison of results. Figure 10 shows a bar graph comparison of the THD Ia factor obtained for the simulation and experimental results for three different cases. From graphs, slightly different harmonics shares for each case are observed. Differences between simulation and experimental results, in terms of harmonic-order share, can be justified with a lack of manufacturer information for the used grid inverter. Namely, the same inverter was not used in the simulation such as the one implemented in the used inverter, nor was the same control algorithm used.

Conclusions
Followed by the announcements of future investments in renewable energy sources, especially with an emphasis on grid-connected PV systems, the energy growth trend is certain, as well as its consequent impact on the power system. Power quality issues and the hosting capacity issue will be the major tasks of future research. Corresponding to that, this article points out the electric power quantities evaluation and followed distortion effects, as the factors that should be considered very important in planning and installing a new PV capacity.
Distortion effects concerning the power quality issues and the impact of the power components were presented for three different cases, depending on different solar irradiances, and based on the implemented simulation model and experimental measurement in real time. By analyzing the results of the simulation model, as well as the experimental setup, it can be concluded that a major issue comes from a low PV power production because of the low irradiation level. The PV inverters appear as generators of the significant share of high-order current harmonics, which result in the existence of a valuable amount of electric power quantities. The low-power factor becomes predominantly dependent on the distortion effect and its value strongly depends on the non-active scattered power component.
It is noted that this impact is not significant from the point of view of the amount of power injected into the grid in bad weather conditions, but it can be assumed that there will be a significant number of installed grid-connected PV systems in some micro-location in the future; therefore, this impact will gradually gain importance.