Harmonic Loss Analysis of Low-Voltage Distribution Network Integrated with Distributed Photovoltaic

: In a power system with highly proportional renewable energy integration, the power generated by photovoltaic (PV) of high permeability and high proportion needs to be connected to the distribution network through the power electronic inverter. The inverters can generate the low-order harmonic and high-order harmonic near the switching frequency. Harmonic power will be generated when the harmonic current ﬂows through the power grid with the harmonic voltage of the same frequency, and the additional harmonic losses caused should not be neglected. To effectively analyze the voltage quality and calculate the harmonic loss of the low-voltage distribution network integrated with distributed PV, based on the harmonic loss factor of resistance proposed, the harmonic impedance modeling and harmonic loss calculation method for the key equipment of the power grid, such as lines and transformers, are introduced in this paper ﬁrstly. Next, a decoupling algorithm of harmonic power ﬂow is proposed, and the inﬂuence of the access capacity of PV on voltage quality and line loss of the distribution network is analyzed. Finally, a harmonic loss calculation method based on measured harmonic data of the distribution network is proposed. It is found that the harmonic loss of the low-voltage distribution network accounts for about 0.6% of the total network loss. Therefore, voltage quality can be improved and line loss can be reduced effectively by reasonable access to PV and reducing harmonic order and the current harmonic distortion.


Introduction
With carbon peak and carbon neutrality goals, the renewable energy generation technology represented by photovoltaic (PV) is widely considered [1,2].Because of clean and sufficient solar energy resources and convenient installation methods, such as "roof photovoltaic", the economic and environmental benefits of the PV generation system are continuously improved [3,4].When a high proportion of distributed energy is connected to the power grid, the change in configuration and operation condition of the distribution network is huge [5] and the power flow of the distribution network has significantly been changed [6,7].On the one hand, distributed PV will affect the power flow of the distribution network, directly affecting the loss of lines and equipment [8].On the other hand, the randomness and intermittence of PV output power will bring a series of problems, such as harmonics [9,10].The interaction with the original power quality disturbance factors further complicates the power quality disturbance mechanism of the distribution network, indirectly increasing the complexity of distribution network loss analysis [11] and the difficulty of reducing loss [12].Therefore, it is of great significance to study the influence mechanism of the network loss of the distribution network integrated with a high proportion of distributed PV.
PV needs to be connected to the distribution network through power electronic inverters [13], resulting in various low-frequency harmonics and high-frequency harmonics near the switching frequency [14,15].Non-ideal factors of the inverter, such as the accumulation of dead time, cause low-order harmonics.The delay characteristics of sinusoidal pulse width modulation (SPWM) will produce high-frequency harmonics [16].At present, there are various types of inverters.Their operation control and non-ideal characteristics of the access terminal will bring harmonic distortion to the distribution network [17,18].Harmonics caused will reduce the insulation level and life of transformers, cables and other equipment, and increase copper, iron and dielectric losses [19].Therefore, harmonics should be analyzed in the distribution network with PV integration.
Harmonic power flow calculation is an important basis for harmonics research [20] and an important means to understand the harmonic characteristics of the power grid and analyze harmonics [21,22].In order to analyze the harmonic power flow, the constant current source model is widely used in engineering [23].The phase angle of the typical spectrum in the constant current source model is determined by experiment or simulation [24].At the same time, this method is only applicable to the analysis of a single harmonic source in the system [25].With large-scale PV integration, the system contains multiple harmonic sources, and, thus, the iterative harmonic power flow becomes effective.In iterative harmonic analysis, the harmonic current of the harmonic source is a variable controlled by its terminal voltage [26].Newton-Raphson or Gauss methods are used to update the terminal voltage of harmonic sources, and harmonic currents are solved by alternating iterations [27].Convergence performance and computation speed are the main factors affecting the application of this method.The unified harmonic power flow method using the Newton-Raphson method can solve the fundamental and harmonic power flows.The deviation equation and Jacobi matrix are simplified, so the calculation speed and convergence are improved [28].The current source direct solution method also is a harmonic power flow analysis method.It is widely used in engineering, and the decoupling method is further simplified.Ignoring the influence of harmonic power flow on the fundamental power flow, the fundamental and harmonic power flows are solved separately.Only the influence of the fundamental voltage of the bus on the harmonic injection current is considered [29].However, the phase angle estimation of harmonic current is not accurate.Therefore, the interaction between multiple harmonic currents cannot be considered when multiple harmonic sources exist in the system.
At present, the influence of harmonics on the active power loss of the low-voltage distribution network has been studied by a large number of papers.In [30], the influence of harmonics on the loss is analyzed, the Norton equivalent method is used to model the nonlinear load of the distribution network and it is verified that the power loss of the distribution network caused by harmonic distortion can be as high as 20%.In [31], a novel harmonic model of the three-phase transformer is proposed.The skin effect and proximity effect of transformer windings under harmonic background are considered in this model.The harmonic loss of the transformer is accurately calculated by using the hysteresis model.In [32], the skin effect and hysteresis and eddy current loss of the stranded wire are considered, and the relative resistance increment caused by the skin effect of the stranded wire is calculated.At the same time, it is not suitable for approximate calculation in engineering, owing to a large number of parameters and complex calculation processes.The existing references are all theoretical studies on the influence of harmonics on single equipment [33], and the proportion of line loss caused by harmonics was not specifically quantified.
To systematically and quantitatively analyze the harmonic loss, a harmonic loss calculation method for the low-voltage distribution network integrated with distributed PV based on measured data is proposed in this paper.The contributions of the paper are as follows: (1) The harmonic impedance modeling and harmonic loss calculation methods for the key equipment of the power grid, including lines and transformers, are introduced, and the quantitative relationship between harmonic current content and harmonic order and harmonic loss is obtained.
(2) The fundamental power flow is solved iteratively by the Newton method, and the harmonic power flow is solved by the decoupling algorithm.The influence of PV capacity on voltage quality and line loss is analyzed.(3) The harmonic loss calculation method based on the measured harmonic data of the low-voltage distribution network is proposed.
The rest of the paper is organized as follows.Section 2 analyzes the operation mechanism of the distributed PV system.Section 3 describes the equipment loss model considering harmonic factors.Section 4 illustrates the harmonic power flow algorithm.Section 5 presents the simulation analysis and Section 6 summarizes the study's conclusions.

PV Generation Model
The distributed PV is mainly based on solar cells, which can convert solar energy into DC energy.A simplified circuit of the solar cells is shown in Figure 1.In this figure, I ph is the photogenerated current; I D is the reverse saturation current of the diode; I p is the leakage current of solar cells; R p is the parallel resistance; R s is the series resistance; U out is output voltage; I out is output current.and the quantitative relationship between harmonic current content and harmonic order and harmonic loss is obtained.(2) The fundamental power flow is solved iteratively by the Newton method, and the harmonic power flow is solved by the decoupling algorithm.The influence of PV capacity on voltage quality and line loss is analyzed.(3) The harmonic loss calculation method based on the measured harmonic data of the low-voltage distribution network is proposed.
The rest of the paper is organized as follows.Section 2 analyzes the operation mechanism of the distributed PV system.Section 3 describes the equipment loss model considering harmonic factors.Section 4 illustrates the harmonic power flow algorithm.Section 5 presents the simulation analysis and Section 6 summarizes the study's conclusions.

PV Generation Model
The distributed PV is mainly based on solar cells, which can convert solar energy into DC energy.A simplified circuit of the solar cells is shown in Figure 1.In this figure, Iph is the photogenerated current; ID is the reverse saturation current of the diode; Ip is the leakage current of solar cells; Rp is the parallel resistance; Rs is the series resistance; Uout is output voltage; Iout is output current.In the simplified circuit, solar cells can be considered equivalent circuits composed of a current source and diode in parallel.The current source simulates the photocurrent and the diode simulates the internal semiconductor material characteristics.The equivalent series resistance Rs is small, which is used to simulate the body resistance of the material and the contact resistance of the electrode.Because PV cells have a certain leakage current, the equivalent parallel resistance Rp is used to simulate the leakage loss.According to Kirchhoff's current law, it can be formulated as:

Rs
The conduction current equation of the diode is: where I0 is the reverse saturation current of the diode; K denotes the Boltzmann constant, and its value is 1.38 × 10 −23 J/K; T is Kelvin thermodynamic temperature; and q represents a unit of electric quantity, and its value is 1.6 × 10 −19 C. Therefore, the output current of solar cells is: In the simplified circuit, solar cells can be considered equivalent circuits composed of a current source and diode in parallel.The current source simulates the photocurrent and the diode simulates the internal semiconductor material characteristics.The equivalent series resistance R s is small, which is used to simulate the body resistance of the material and the contact resistance of the electrode.Because PV cells have a certain leakage current, the equivalent parallel resistance R p is used to simulate the leakage loss.According to Kirchhoff's current law, it can be formulated as: The conduction current equation of the diode is: where I 0 is the reverse saturation current of the diode; K denotes the Boltzmann constant, and its value is 1.38 × 10 −23 J/K; T is Kelvin thermodynamic temperature; and q represents a unit of electric quantity, and its value is 1.6 × 10 −19 C. Therefore, the output current of solar cells is: The DC cannot be directly provided to users.It is necessary to convert the DC power into AC power with frequency and amplitude matching the grid through the inverter part [34].The structure of the PV grid-connected is shown in Figure 2. The voltage source inverter by PWM control is used in the main circuit.The voltage and current double closed-loop control is used as the control mode [35].The output AC power of the inverter part contains more harmonics, especially high-frequency harmonics, which cannot be directly incorporated into the power grid.The filter is installed to filter out harmonic current at the AC side of the inverter.At present, the commonly used filters are L-, LCand LCL-types.The LCL-type filter is widely used because of its superior high-frequency harmonic suppression ability without large inductance.It also has good performance and a low price.
( ) The DC cannot be directly provided to users.It is necessary to convert the DC power into AC power with frequency and amplitude matching the grid through the inverter part [34].The structure of the PV grid-connected is shown in Figure 2. The voltage source inverter by PWM control is used in the main circuit.The voltage and current double closedloop control is used as the control mode [35].The output AC power of the inverter part contains more harmonics, especially high-frequency harmonics, which cannot be directly incorporated into the power grid.The filter is installed to filter out harmonic current at the AC side of the inverter.At present, the commonly used filters are L-, LC-and LCLtypes.The LCL-type filter is widely used because of its superior high-frequency harmonic suppression ability without large inductance.It also has good performance and a low price.
The structure of the PV grid-connected.

Harmonic Generation Mechanism
There are many unsatisfactory factors of the PV inverter.For example, the low-frequency harmonics are generated by the accumulation of dead time, and the high-frequency harmonics are generated by the switching characteristics of sinusoidal pulse width modulation.At the same time, the inverter types are different and various, and their operation control and non-ideal characteristics will bring harmonic distortion to the distribution network.The distributed integration of PV with high permeability results in high harmonic distortion, strong randomness and wide harmonic frequency coverage, which increases the copper, iron and dielectric losses of equipment and reduces the insulation level and service life of transformers and cables.
To ensure the normal operation of the switching device, the dead time must be added to the power electronic switching device.However, the zero current clamping effect is produced and the low harmonic distortion of inverter output increases.The error voltage expression is: where Udc is the DC voltage of the PV inverter; fc is carrier frequency; n is harmonic order relative to modulated wave; ϕ is the initial phase of the modulated wave; Td is dead time.
According to the equation, the harmonic will also increase with the increase of dead time, DC voltage and carrier frequency.The order of harmonic voltage caused by the dead time effect is 6k ± 1, k = 1, 2, 3, …., so the 5th, 7th, 11th et al. odd harmonic voltage can be obtained by Fourier decomposition.
The working mechanism of the inverter determines that SPWM modulation will produce high-frequency harmonics.With bipolar SPWM modulation, the inverter output line

Harmonic Generation Mechanism
There are many unsatisfactory factors of the PV inverter.For example, the lowfrequency harmonics are generated by the accumulation of dead time, and the highfrequency harmonics are generated by the switching characteristics of sinusoidal pulse width modulation.At the same time, the inverter types are different and various, and their operation control and non-ideal characteristics will bring harmonic distortion to the distribution network.The distributed integration of PV with high permeability results in high harmonic distortion, strong randomness and wide harmonic frequency coverage, which increases the copper, iron and dielectric losses of equipment and reduces the insulation level and service life of transformers and cables.
To ensure the normal operation of the switching device, the dead time must be added to the power electronic switching device.However, the zero current clamping effect is produced and the low harmonic distortion of inverter output increases.The error voltage expression is: where U dc is the DC voltage of the PV inverter; f c is carrier frequency; n is harmonic order relative to modulated wave; ϕ is the initial phase of the modulated wave; T d is dead time.
According to the equation, the harmonic will also increase with the increase of dead time, DC voltage and carrier frequency.The order of harmonic voltage caused by the dead time effect is 6k ± 1, k = 1, 2, 3, . . .., so the 5th, 7th, 11th et al. odd harmonic voltage can be obtained by Fourier decomposition.The working mechanism of the inverter determines that SPWM modulation will produce high-frequency harmonics.With bipolar SPWM modulation, the inverter output line voltage U ab can be obtained.By Fourier decomposition, the fundamental voltage component U ab1 is: where m is the harmonic order of carrier u c .The carrier angular frequency is ω c , and the modulation wave angular frequency is ω r .When the frequency is nω r ± kω c , the harmonic components U abn can be expressed as follows.
where J n is the Bessel function; n is the harmonic order; M is the degree of modulation.
When the degree of modulation M = 1, the voltage utilization rate is 0.886, there is no integer multiple harmonic of carrier frequency ω c and there is no integral multiple harmonic of 3 of the modulated wave frequency ω r .According to the properties of the Bessel function, its change rule is similar to that of the trigonometric function.The amplitude of each harmonic is related to the degree of modulation M. As the carrier ratio increases, the amplitude of the harmonic voltage increases or decreases.When appropriate filtering schemes and control methods are used, specific frequency harmonics can be eliminated or reduced.
In practice, a distributed PV system usually works in the "full access" power supply mode.In addition to the inevitable characteristic harmonics caused by the normal operation and pulse modulation of the PV inverter, the output current of the PV inverter also contains interharmonics.Changes in light intensity and temperature lead to changes in PV output power, and interharmonics are generated in the grid.A PV grid-connected system using a maximum power point tracking strategy will cause the DC voltage command value to fluctuate.Rich interharmonic currents are generated in the grid through the interaction of harmonics/interharmonics on the AC/DC side of the inverter.In addition, the voltage background interharmonics at the PV grid connection point can cause it to inject interharmonic currents into the grid.ω m is the interharmonic angular frequency of the DC side, and the harmonic frequency of the AC side will be shifted by ±ω.

Influence of PV Output Power on Harmonic Characteristics
The vector diagram of the grid-connected inverter under the steady-state condition is shown in Figure 3.In this figure, U s is the effective value of phase voltage on the AC side; I is the grid-connected current; ϕ is the power angle factor; U inv is the output phase voltage effective value of the inverter; ω 1 is the fundamental frequency of the power grid.The capacitance current of the filter at the fundamental frequency is small and can be ignored.If the active power of the PV inverter output is Ppv and the reactive power is Qpv, the effective value of the grid current is: The output voltage effective value of the inverter Uinv is: According to the law of cosines: The capacitance current of the filter at the fundamental frequency is small and can be ignored.If the active power of the PV inverter output is P pv and the reactive power is Q pv , the effective value of the grid current is: The output voltage effective value of the inverter U inv is: According to the law of cosines: Therefore, it can be obtained that: The change in the degree of modulation is positively related to the change of output power, so the PV output power affects the output harmonic amplitude of the inverter.

Calculation Model of Line Harmonic Loss
The influencing factors of line harmonic loss are the magnitude of harmonic current and harmonic resistance.Harmonic current is related to the operation characteristics of the power grid.Under the harmonic background, the line resistance is affected by the harmonic order and line skin effect.The influence of the proximity effect, the hysteresis and eddy current of the steel core and the protection of the steel strip also should be considered.
Without considering the influence of harmonics on the impedance of transmission lines, the line harmonic loss of the distribution network can be expressed as: where I is the sum of fundamental current and harmonic current; R is line equivalent resistance; I 1 is fundamental current; I n is nth harmonic current; THD I is the total harmonic distortion of current; M is the highest harmonic order.The loss increase rate of distribution network lines is the ratio of the loss generated by harmonics to the loss of fundamental current: In practice, the impedance of the line changes with the harmonic frequency, especially when the harmonic frequency is higher.Wakileh et al. proposed that with the increase of the frequency, the skin effect becomes more obvious, and the current in the conductor concentrates on the surface, resulting in an increase in AC resistance and a decrease in internal resistance and inductance.Considering the skin effect under harmonic conditions, the impedance of the conductor is: The harmonic loss of the line of the simplified model can be expressed as: where R 1 is line fundamental frequency resistance; HRI n is the ratio of harmonic current to fundamental current.The effects of permeability, conductor radius and conductivity are neglected in this method.The calculation results may be quite different from the actual ones.
The calculation standard JCS 0374 for bare wire carrying capacity was published by Japan Wire Industry Association [36].The AC/DC resistance coefficient is introduced through the AC/DC resistance test of wire: where R AC is AC resistance; R DC is DC resistance; k 1 is the skin effect coefficient; k 2 is the core loss coefficient.Formula (17) shows the AC/DC resistance ratio at different frequencies.
The skin effect of the conductor and the iron loss caused by hysteresis and the eddy current of the steel core are considered.k 1 can be written as: where a = 0.99609; b = 0.18578; c = −0.030263;d = 0.020735.For single strand (aluminum alloy wire and aluminum alloy core wire), x can be written as: For composite wire (aluminum wire and steel wire), x can be written as: where f is frequency; D 1 is the outer diameter of the strand; D 2 is the outer diameter of steel wire.For the single strand and composite strand (the number of aluminum wire layers is even), the iron loss coefficient k 2 = 1.For composite strands (more than 1 odd number of aluminum layers), k 2 can be written as: where a = 0.99947; b = 0.028859; c = −0.0059348;d = 0.00042259; I is stranded wire current; A is section area of aluminum wire.If the values of k 1 and k 2 are less than 1, k 1 = 1 and k 2 = 1.This method has comprehensive factors and relatively simple calculation, which is more suitable for AC resistance calculation under harmonic background.Therefore, the loss power of the nth harmonic of the line can be written as:

Calculation Model of Transformer Harmonic Loss
Distribution transformer loss mainly has two parts, which are copper loss and iron loss.The excitation branch of the transformer is nonlinear due to the existence of the iron core, and its nonlinearity varies with the applied voltage.The higher the voltage is, the closer the core is to saturated, and the greater the degree of nonlinearity is.When the applied voltage is not high, the core is not saturated, and the voltage harmonic distortion is not large.Therefore, the influence of harmonics on the iron loss of the transformer can be ignored.Without considering the influence of harmonics on the excitation branch of the distribution transformer, the harmonic equivalent circuit of the distribution transformer is shown in Figure 4.

Calculation Model of Transformer Harmonic Loss
Distribution transformer loss mainly has two parts, which are copper loss and iron loss.The excitation branch of the transformer is nonlinear due to the existence of the iron core, and its nonlinearity varies with the applied voltage.The higher the voltage is, the closer the core is to saturated, and the greater the degree of nonlinearity is.When the applied voltage is not high, the core is not saturated, and the voltage harmonic distortion is not large.Therefore, the influence of harmonics on the iron loss of the transformer can be ignored.Without considering the influence of harmonics on the excitation branch of the distribution transformer, the harmonic equivalent circuit of the distribution transformer is shown in Figure 4.The impedance is composed of the winding resistance and leakage reactance of the transformer, and the inductance in the leakage reactance is often approximated as a constant.Therefore, the leakage reactance of the transformer XTn is the product of the fundamental frequency reactance and the harmonic order, which can be expressed as: The impedance is composed of the winding resistance and leakage reactance of the transformer, and the inductance in the leakage reactance is often approximated as a constant.Therefore, the leakage reactance of the transformer X Tn is the product of the fundamental frequency reactance and the harmonic order, which can be expressed as:

RTn XTn
where X T1 is the fundamental frequency reactance of the transformer.When there are harmonics in the transformer, due to the skin effect and proximity effect, the resistance of the winding will increase.It is usually considered to be proportional to the square root of the harmonic order.Therefore, the harmonic impedance formula of the simplified model can be expressed as: where R T1 is the fundamental frequency winding resistance of the transformer.This method is suitable for the simple estimation of harmonic line loss in engineering practice, and there will be a large error when there are high-frequency harmonics in the system.Under the harmonic background, due to the winding skin effect, the resistance of the transformer has the following relationship with the harmonic order [37]: c 0 , c 1 , c 2 and b are coefficients, which are related to the material and length of winding.Generally, c 0 + c 1 + c 2 = 1.The value range of each coefficient is shown in Table 1 [37].Harmonic loss of transformers can be expressed as: where K is the constant coefficient.This method is simple to calculate, and fewer factors are considered.The eddy current and stray loss under harmonic conditions are not considered.
In addition, because the constant coefficient in Formula ( 26) is not fixed, it is related to artificial selection, and the calculation error is large.
The transformer load loss is divided into DC loss, winding eddy current loss and other stray losses in IEEE/ANSI C57.110 standard [38].Among them, the winding resistance loss is the main part of the loss, accounting for about 80% of the loss.The loss of the winding eddy currently accounts for about 33% of the additional loss.The winding eddy current harmonic factor F HL−EC and other stray loss harmonic factor F HL−OSL are used to calculate the additional loss of the transformer under harmonics.They can be expressed as: Considering the effect of skin effect on winding resistance under harmonic current, the harmonic loss factor of winding resistance F HL−R is defined as The total transformer loss under harmonic background can be expressed as: where P NL is fixed loss; P R , P EC-R and P OSL-R are resistance loss, winding eddy current loss and other stray loss under rated operating conditions.It can be found that the loss is mainly determined by the secondary side current.The effect of skin effect on transformer winding resistance is considered, and the winding resistance loss is calculated in Formula ( 13).It is more accurate than the calculation method proposed in IEEE Std C.57.110.

Harmonic Power Flow Calculation Model of Distribution Network
Fundamental power flow and harmonic power flow are decoupled by the decoupling method in the harmonic power flow analysis.In the analysis of fundamental power flow, the influence of harmonics is not considered, and the fundamental voltage and current are obtained independently.Based on the voltage and current of the fundamental power flow, the harmonic solution is realized by combining the harmonic source model and the harmonic network equation.The calculation process of harmonic power flow is simplified, and the accuracy of the results are satisfactory by this method.The calculation steps are as follows.
Step 1: The solution of fundamental voltage and fundamental power flow is realized based on the fundamental power flow equation.The fundamental power flow equation is as follows: where the subscript is the node number; P i and Q i are fundamental active and reactive power; V i and V j represent the node voltage amplitudes; G ij and B ij are conductance and susceptance parameters of the line; δ ij represents the phase angle difference of node voltage.The traditional Newton-Raphson method is usually used to solve the fundamental power flow in the distribution network, and the fundamental voltage data of the harmonic source node are obtained.
Step 2: The harmonic source model of PV is established.In general, the harmonic distortion of each node voltage is less than the limit value specified in the national standard, and the harmonic voltage is far less than the fundamental voltage in numerical value.The effect of harmonic voltage on harmonic injection current is also very small, which can be approximately ignored.Therefore, the harmonic source is regarded as a harmonic constant current source with infinite internal impedance.The amplitude and phase angle of har-monic current depends on the typical frequency spectrum of harmonic source equipment.It can be expressed as: (32) where I h and θ h are hth harmonic current amplitude and phase angle; I f and θ f are the magnitude and phase angle of fundamental current; I h,spectrum and θ h,spectrum are the spectrum values of the hth harmonic current magnitude and phase; I f,spectrum and θ f,spectrum are the spectrum values of the fundamental current magnitude and phase.
Step 3: To solve the harmonic current of the network, it is necessary to obtain the network equation of the system under harmonic frequency in harmonic power flow analysis.Therefore, generators, lines, transformers and linear loads need to be modeled.The specific process is shown in Appendix A.
Step 4: According to the fundamental power flow data and harmonic source model, the hth harmonic current matrix I h is obtained.Based on the nodal admittance matrix Y h and its inverse matrix Z h under a certain harmonic order, the solution of harmonic voltage matrix V h is obtained.
Step 5: If the maximum deviation of voltage amplitude of the node i in the kth iteration is smaller than ε (ε is 10 −7 in this paper) in Formula (34), the iteration can be stopped and the ultimate harmonic voltages are yielded e = max(|V Otherwise, the harmonic source is updated, and the above process is repeated until the calculation results converge.N is a collection of node numbers.V

(k)
i (h) and V (k−1) i (h) are the harmonic voltage of this iteration and the harmonic voltage of the previous iteration.
Step 6: Calculating the active power loss of branch ij.
where V (h)i is the hth harmonic voltage of node i.According to the above solution process, the flow chart of harmonic power flow calculation is shown in Figure 5.
Step 6: Calculating the active power loss of branch ij.
where V(h)i is the hth harmonic voltage of node i.According to the above solution process, the flow chart of harmonic power flow calculation is shown in Figure 5.

Harmonic source modeling and analyze the harmonic injection current
Calculate the branch harmonic current and the node harmonic voltage by the decoupling algorithm Obtain harmonic current amplitude and phase of harmonic source node Obtain the nodal admittance matrix

Input the original data of power flow calculation
Obtain the results of fundamental power flow through Newton-Raphson method Calculate the active power loss of branch The flow chart of harmonic power flow calculation.

Simulation Analysis
In this section, the quantitative relationship between harmonic current content and harmonic order and harmonic loss is analyzed.Then, the influence of the access capacity of PV on voltage quality and line loss of the distribution network is analyzed.Finally, the harmonic loss ratio based on the measured harmonic data of the low-voltage distribution network is calculated.The verification is accomplished using MATLAB 2019a with an Intel i7-10700F CPU and a 16 GB RAM computer.

Simulation Analysis
In this section, the quantitative relationship between harmonic current content and harmonic order and harmonic loss is analyzed.Then, the influence of the access capacity of PV on voltage quality and line loss of the distribution network is analyzed.Finally, the harmonic loss ratio based on the measured harmonic data of the low-voltage distribution network is calculated.The verification is accomplished using MATLAB 2019a with an Intel i7-10700F CPU and a 16 GB RAM computer.

Quantitative Analysis of Harmonic Loss in Distribution Network
In the case of harmonics, the line and transformer losses cannot be ignored.According to the line and transformer harmonic line loss model given above, the resistance iron loss coefficient and skin effect coefficient are selected to calculate the line harmonic loss, and the transformer harmonic loss factor is selected to calculate the transformer harmonic load loss.Since the harmonic voltage of the system usually meets the national standard limit, the harmonic loss caused by harmonic voltage is very small compared with the harmonic current.Therefore, this paper focuses on the additional loss caused by the harmonic current.The JLHA2 aluminum-alloy strand with a nominal cross-sectional area of 240 mm 2 is selected for the distribution line.Its DC resistance is 0.11 Ω/km and the line length is 2.2 km.An S11-200/10 distribution transformer is selected, and the transformer connection group number is Dyn11.Specific parameters are shown in Table 2 [39].To verify the accuracy of the calculation model of harmonic loss, a simulation model is built in Matlab/Simulink.The parameters of the distribution line and distribution transformer in Matlab are converted.Specific parameters are shown in Table 3.The block diagram of the Matlab/Simulink model is shown in Figure 6.The threephase programmable voltage source is selected to simulate the distribution network system under the background of 7th and 11th harmonics.The electric energy-measuring device consists of the three-phase instantaneous power measuring element and integral link.The loss of the line and transformer are obtained by the difference value of the power meters.To verify the accuracy of the proposed model, the model proposed in this paper (case 1), simplified model (case 2) and simulation results (case 3) are compared in Figure 7.  From the above figure, it can be seen that the simulation result is close to the quantitative calculation result, and the simulation value is slightly larger than the calculated value.Because there are still some differences between the line model and the actual situation, there are always errors between the calculation results and the simulation values of the line, but the error changes little.The error is basically stable at around 4.2%, meeting the requirements of engineering calculation.The calculation error of the transformer increases with the increase of harmonic distortion, but the overall error is not more than 5%.The reason is that the effect of additional iron loss is ignored in the calculation model used in this paper.However, the loss measured by the meter includes additional iron loss caused by harmonics in the simulation.The errors of the simplified model of the line and the transformer are 20.85% and 55.58%.The above analysis shows that the iron loss increases with the increase of harmonic voltage distortion.In practice, harmonic voltage generally meets the national standard, and this part of loss can be ignored.The effects of permeability, conductor radius and conductivity are neglected in the simplified calculation model, and the calculation results are quite different from the actual ones.The higher the harmonic frequency, the greater the error.From the above figure, it can be seen that the simulation result is close to the quantitative calculation result, and the simulation value is slightly larger than the calculated value.Because there are still some differences between the line model and the actual situation, there are always errors between the calculation results and the simulation values of the line, but the error changes little.The error is basically stable at around 4.2%, meeting the requirements of engineering calculation.The calculation error of the transformer increases with the increase of harmonic distortion, but the overall error is not more than 5%.The reason is that the effect of additional iron loss is ignored in the calculation model used in this paper.However, the loss measured by the meter includes additional iron loss caused by harmonics in the simulation.The errors of the simplified model of the line and the transformer are 20.85% and 55.58%.The above analysis shows that the iron loss increases with the increase of harmonic voltage distortion.In practice, harmonic voltage generally meets the national standard, and this part of loss can be ignored.The effects of permeability, conductor radius and conductivity are neglected in the simplified calculation model, and the calculation results are quite different from the actual ones.The higher the harmonic frequency, the greater the error.
To quantitatively analyze the additional losses of lines and distribution transformers under different harmonic current distortions and different harmonic orders, the 5th, 7th, 11th, 13th and 17th harmonic losses are calculated separately under the condition of 5%, 7%, 11%, 13% and 17% harmonic current distortion of the lines and transformers.The secondary rated current of the distribution transformer is generally used as the fundamental current for the transformer and line.From Formulas ( 17)-( 19), ( 22) and ( 27)-( 30), the figures of the relationship between harmonic loss of lines and transformers and harmonic distortion and harmonic order are drawn in Figure 8.To quantitatively analyze the additional losses of lines and distribution transformers under different harmonic current distortions and different harmonic orders, the 5th, 7th, 11th, 13th and 17th harmonic losses are calculated separately under the condition of 5%, 7%, 11%, 13% and 17% harmonic current distortion of the lines and transformers.The secondary rated current of the distribution transformer is generally used as the fundamental current for the transformer and line.From Formulas ( 17)-( 19), ( 22) and ( 27)- (30), the figures of the relationship between harmonic loss of lines and transformers and harmonic distortion and harmonic order are drawn in Figure 8.The harmonic order of current Harmonic loss of Line (W)  According to the results, when the harmonic order is constant, the line harmonic loss is proportional to the square of harmonic distortion of current, and the higher harmonic of the transformer is approximately proportional to the square of harmonic distortion of current.The harmonic loss of the line is mainly affected by the AC resistance coefficient, and the change is gentle.The harmonic loss of the transformer is mainly affected by the harmonic loss factor and changes rapidly.
When the current harmonic distortion is the same, the higher the harmonic order, the According to the results, when the harmonic order is constant, the line harmonic loss is proportional to the square of harmonic distortion of current, and the higher harmonic of the transformer is approximately proportional to the square of harmonic distortion of current.The harmonic loss of the line is mainly affected by the AC resistance coefficient, and the change is gentle.The harmonic loss of the transformer is mainly affected by the harmonic loss factor and changes rapidly.
When the current harmonic distortion is the same, the higher the harmonic order, the greater the harmonic loss of the low-voltage distribution network.This is because the higher the frequency, the greater the skin effect and proximity effect of the winding and line, so the greater the harmonic impedance of the corresponding transformer and line.Further analysis shows that when the 17th harmonic distortion reaches 17%, the distribution network loss increases by nearly 15%.Therefore, the harmonic has a great influence on the distribution network loss.The copper loss of transformers and transmission line loss increase with the increase of harmonic current content and harmonic order.The proportion of copper loss is much larger than that of line loss and transformer iron loss.Therefore, limiting the content and order of transformer harmonic current is of great significance in reducing branch loss.

Power Flow Simulation Analysis of Harmonic Loss in Distribution Network
To verify the influence of capacity on network voltage and loss, the impact of distributed PV generation is analyzed.In the MATLAB programming environment, the bus voltage and branch loss of the IEEE 33-bus distribution system are analyzed, and the structure of the IEEE 33-bus system is shown in Figure 9 [40].The reference voltage of the system is 12.66 kV, and the reference power is 10 MVA, where bus 1 is the balanced node.Distributed PV is connected to bus 16 and bus 33, and the maximum output power of PV is 0.5 MW when the reference light radiation is 1 KW/m 2 .
Sustainability 2023, 15, x FOR PEER REVIEW 15 of 24 of the IEEE 33-bus system is shown in Figure 9 [40].The reference voltage of the system is 12.66 kV, and the reference power is 10 MVA, where bus 1 is the balanced node.Distributed PV is connected to bus 16 and bus 33, and the maximum output power of PV is 0.5 MW when the reference light radiation is 1 KW/m 2 .The frequency spectrum of harmonics generated by grid-connected PV is shown in Table 4.The load is modeled according to the above harmonic impedance model.Highfrequency harmonics are not considered.The low-frequency harmonics generated by grid-connected PV are analyzed, and the highest harmonic order is 7th order.To analyze the fundamental power flow, the PV bus is simplified to the PQ node.The specific parameters of the PV output of each case are listed in Table 5.The influence of different PV capacities on fundamental voltage is shown in Figure The frequency spectrum of harmonics generated by grid-connected PV is shown in Table 4.The load is modeled according to the above harmonic impedance model.High-frequency harmonics are not considered.The low-frequency harmonics generated by grid-connected PV are analyzed, and the highest harmonic order is 7th order.To analyze the fundamental power flow, the PV bus is simplified to the PQ node.The specific parameters of the PV output of each case are listed in Table 5.The influence of different PV capacities on fundamental voltage is shown in Figure 10.With the increase of PV permeability, the voltage of each bus also increases.Especially, compared with the bus voltage of the distribution network without PV connection, the maximum rate of change of the bus voltage is 5.37% in Case 5. Therefore, the higher the output power of PV is, the higher the permeability is and the better the supporting effect on the bus voltage is.The above voltage rise is due to the decrease of the equivalent load of the line before the PV access bus, so the voltage drop decreases and the voltage is raised.At the same time, for the branch after the PV access bus, the power flow through the line does not change, so the voltage drop remains unchanged.Due to the front bus voltage being raised, the back bus voltage will also be raised.With the increase of the PV capacity, the change in voltage is greater.Especially in the case of small load demand and large PV output power, the voltage increases significantly.From the results in Figure 11, with the increase of PV capacity, the loss of each branch decreases.When the distributed grid-connected PV capacity is in a certain range, the power transmitted by the upper power grid is reduced by increasing the light intensity.The power flowing through the branch is reduced, so the loss of distribution is effectively reduced.Compared with the loss of the distribution network without PV connection, the loss is reduced by 56.26% in Case 5.However, when the access capacity of PV exceeds a certain range, the reverse current of some branches exceeds the forward current without PV access, causing an increase in the network loss.From Figure 12, with the continuous increase of distributed grid-connected PV capacity on a bus, the harmonic voltage distortion of the whole distribution system in- From the results in Figure 11, with the increase of PV capacity, the loss of each branch decreases.When the distributed grid-connected PV capacity is in a certain range, the power transmitted by the upper power grid is reduced by increasing the light intensity.The power flowing through the branch is reduced, so the loss of distribution is effectively reduced.Compared with the loss of the distribution network without PV connection, the loss is reduced by 56.26% in Case 5.However, when the access capacity of PV exceeds a certain range, the reverse current of some branches exceeds the forward current without PV access, causing an increase in the network loss.
From Figure 12, with the continuous increase of distributed grid-connected PV capacity on a bus, the harmonic voltage distortion of the whole distribution system increases.In addition, the higher the harmonic frequency is, the more obvious the influence of skin effect on line impedance is, and the greater the influence of harmonic voltage is.The bus with the highest voltage THD is bus 16 in Case 5.The value of the voltage THD is 3.54%.For the PV access bus in the middle and lower reaches of the system, the overall harmonic voltage of the system is greatly affected by PV capacity, especially the voltage of the bus near the PV access bus.Therefore, when distributed PV is connected to the grid, grid-connected PV capacity needs to be controlled within a reasonable range.When the large-capacity distributed PV is connected to the grid, it should be connected by multiple points to reduce the harmonic voltage level of the whole branch.
The resistance of the conductor increases when there are harmonics, and the higher the frequency is, the more obvious the skin effect is and the greater the harmonic resistance is.The harmonic loss of the line increases with the increase of PV capacity and the increase of the harmonic order in Figure 13.In addition, the 5th harmonic loss accounts for 1.12% of the fundamental loss, and the 7th harmonic loss accounts for 1.30% of the fundamental loss in Case 5.In conclusion, the voltage quality can be effectively improved, and line loss can be reduced through reasonable planning of PV access capacity and reducing harmonic order and harmonic distortion.From the results in Figure 11, with the increase of PV capacity, the loss of each branch decreases.When the distributed grid-connected PV capacity is in a certain range, the power transmitted by the upper power grid is reduced by increasing the light intensity.The power flowing through the branch is reduced, so the loss of distribution is effectively reduced.Compared with the loss of the distribution network without PV connection, the loss is reduced by 56.26% in Case 5.However, when the access capacity of PV exceeds a certain range, the reverse current of some branches exceeds the forward current without PV access, causing an increase in the network loss.From Figure 12, with the continuous increase of distributed grid-connected PV capacity on a bus, the harmonic voltage distortion of the whole distribution system increases.In addition, the higher the harmonic frequency is, the more obvious the influence of skin effect on line impedance is, and the greater the influence of harmonic voltage is.The bus with the highest voltage THD is bus 16 in Case 5.The value of the voltage THD is 3.54%.For the PV access bus in the middle and lower reaches of the system, the overall harmonic voltage of the system is greatly affected by PV capacity, especially the voltage of the bus near the PV access bus.Therefore, when distributed PV is connected to the grid, grid-connected PV capacity needs to be controlled within a reasonable range.When the large-capacity distributed PV is connected to the grid, it should be connected by multiple points to reduce the harmonic voltage level of the whole branch.The resistance of the conductor increases when there are harmonics, and the higher the frequency is, the more obvious the skin effect is and the greater the harmonic resistance is.The harmonic loss of the line increases with the increase of PV capacity and the increase of the harmonic order in Figure 13.In addition, the 5th harmonic loss accounts for 1.12% of the fundamental loss, and the 7th harmonic loss accounts for 1.30% of the fundamental loss in Case 5.In conclusion, the voltage quality can be effectively improved, and line loss can be reduced through reasonable planning of PV access capacity and reducing harmonic order and harmonic distortion.The resistance of the conductor increases when there are harmonics, and the higher the frequency is, the more obvious the skin effect is and the greater the harmonic resistance is.The harmonic loss of the line increases with the increase of PV capacity and the increase of the harmonic order in Figure 13.In addition, the 5th harmonic loss accounts for 1.12% of the fundamental loss, and the 7th harmonic loss accounts for 1.30% of the fundamental loss in Case 5.In conclusion, the voltage quality can be effectively improved, and line loss can be reduced through reasonable planning of PV access capacity and reducing harmonic order and harmonic distortion.

Analysis of Measured Data of Low-Voltage Distribution Network
In order to analyze the power quality problems and harmonic losses caused by distributed PV, the measured data are sampled and analyzed.The experimental schematic diagram of the measurement is shown in Figure 14.The data measurement point is the 200 kVA 10 kV/0.4 kV low-voltage outlet side of the transformer, and the distribution network area has five distributed PV accesses.The current and voltage data are sampled by  The time-varying trends of the fundamental currents of phase A, phase B and phase C in one day are shown in Figure 15.The power generation reaches its peak at noon and is low in the morning and evening.The current of the low-voltage distribution network in the daytime is mainly related to PV output, with a large value.The night current is determined by the user's power load.Therefore, the variation trend of power generation of PV is affected by natural environmental factors such as light intensity and temperature.Taking phase A as an example, the curve of the harmonic current amplitude of phase A varying with time and average and 95th percentile individual demand distortion (IDD) are obtained.The time-varying trends of harmonic current in a day are shown in Figure 16.The amplitude of the harmonic current is below 7A, and the harmonic current increases at noon.The harmonic current order of PV output is mainly concentrated in the 5th, 7th and 13th, and the corresponding harmonic distortion decreases with the increase of harmonic The time-varying trends of the fundamental currents of phase A, phase B and phase C in one day are shown in Figure 15.The power generation reaches its peak at noon and is low in the morning and evening.The current of the low-voltage distribution network in the daytime is mainly related to PV output, with a large value.The night current is determined by the user's power load.Therefore, the variation trend of power generation of PV is affected by natural environmental factors such as light intensity and temperature.Taking phase A as an example, the curve of the harmonic current amplitude of phase A varying with time and average and 95th percentile individual demand distortion (IDD) are obtained.The time-varying trends of the fundamental currents of phase A, phase B and phase C in one day are shown in Figure 15.The power generation reaches its peak at noon and is low in the morning and evening.The current of the low-voltage distribution network in the daytime is mainly related to PV output, with a large value.The night current is determined by the user's power load.Therefore, the variation trend of power generation of PV is affected by natural environmental factors such as light intensity and temperature.Taking phase A as an example, the curve of the harmonic current amplitude of phase A varying with time and average and 95th percentile individual demand distortion (IDD) are obtained.The time-varying trends of harmonic current in a day are shown in Figure 16.The amplitude of the harmonic current is below 7A, and the harmonic current increases at noon.The harmonic current order of PV output is mainly concentrated in the 5th, 7th and 13th, and the corresponding harmonic distortion decreases with the increase of harmonic order.This is determined by the working characteristics of the 6-pulse PV inverter.The time-varying trends of harmonic current in a day are shown in Figure 16.The amplitude of the harmonic current is below 7A, and the harmonic current increases at noon.The harmonic current order of PV output is mainly concentrated in the 5th, 7th and 13th, and the corresponding harmonic distortion decreases with the increase of harmonic order.This is determined by the working characteristics of the 6-pulse PV inverter.The daily trend of current THD is shown in Figure 17.During the day, PV output power leads to the rapid increase of fundamental current.The change rate of fundamental current is faster than that of harmonic current, so the current THD decreases.The increase of nonlinear load at night and in the morning causes the increase of current THD.In order to analyze the harmonic loss of the station area in one day based on the measured data, a simplified equivalent circuit for the distribution network is shown in Figure 18 [33].The lines and transformer used in the low-voltage distribution network are the same as the line and transformer selected in Section 3. Based on the measured threephase voltage, current data of the low-voltage distribution network and the calculation model of harmonic losses of line and transformer proposed in this paper, the percentage of the harmonic loss power of the low-voltage distribution network to the total network loss is obtained.As shown in Figure 19, The time-varying trends of the ratio of harmonic loss power to total loss power are similar to the current THD.It is found that the calculation results The daily trend of current THD is shown in Figure 17.During the day, PV output power leads to the rapid increase of fundamental current.The change rate of fundamental current is faster than that of harmonic current, so the current THD decreases.The increase of nonlinear load at night and in the morning causes the increase of current THD.The daily trend of current THD is shown in Figure 17.During the day, PV output power leads to the rapid increase of fundamental current.The change rate of fundamental current is faster than that of harmonic current, so the current THD decreases.The increase of nonlinear load at night and in the morning causes the increase of current THD.In order to analyze the harmonic loss of the station area in one day based on the measured data, a simplified equivalent circuit for the distribution network is shown in Figure 18 [33].The lines and transformer used in the low-voltage distribution network are the same as the line and transformer selected in Section 3. Based on the measured threephase voltage, current data of the low-voltage distribution network and the calculation model of harmonic losses of line and transformer proposed in this paper, the percentage of the harmonic loss power of the low-voltage distribution network to the total network loss is obtained.As shown in Figure 19, The time-varying trends of the ratio of harmonic loss power to total loss power are similar to the current THD.It is found that the calculation results In order to analyze the harmonic loss of the station area in one day based on the measured data, a simplified equivalent circuit for the distribution network is shown in Figure 18 [33].The lines and transformer used in the low-voltage distribution network are the same as the line and transformer selected in Section 3. Based on the measured threephase voltage, current data of the low-voltage distribution network and the calculation model of harmonic losses of line and transformer proposed in this paper, the percentage of the harmonic loss power of the low-voltage distribution network to the total network loss is obtained.measured data, a simplified equivalent circuit for the distribution network is shown in Figure 18 [33].The lines and transformer used in the low-voltage distribution network are the same as the line and transformer selected in Section 3. Based on the measured threephase voltage, current data of the low-voltage distribution network and the calculation model of harmonic losses of line and transformer proposed in this paper, the percentage of the harmonic loss power of the low-voltage distribution network to the total network loss is obtained.As shown in Figure 19, The time-varying trends of the ratio of harmonic loss power to total loss power are similar to the current THD.It is found that the calculation results are basically consistent with the Matlab simulation results after verification.The ratio is As shown in Figure 19, The time-varying trends of the ratio of harmonic loss power to total loss power are similar to the current THD.It is found that the calculation results are basically consistent with the Matlab simulation results after verification.The ratio is approximately 0.2~1.6%.Taking one day as the statistical period, the branch loss of the station area with the power supply radius of 220 m is 66.7326 kWh, and the harmonic line loss is about 0.399 kWh.The average harmonic loss of the distribution network is about 0.6%.Although the value is small, due to the large loss base of the large-scale low-voltage distribution network, the power loss caused by harmonics is also quite large.Moreover, the harmonic loss power of the low-voltage distribution network is affected by transformer and line type, power supply radius, load rate and harmonic current content rate.Different distribution networks may have different harmonic losses under different operating conditions.approximately 0.2~1.6%.Taking one day as the statistical period, the branch loss of the station area with the power supply radius of 220 m is 66.7326 kWh, and the harmonic line loss is about 0.399 kWh.The average harmonic loss of the distribution network is about 0.6%.Although the value is small, due to the large loss base of the large-scale low-voltage distribution network, the power loss caused by harmonics is also quite large.Moreover, the harmonic loss power of the low-voltage distribution network is affected by transformer and line type, power supply radius, load rate and harmonic current content rate.Different distribution networks may have different harmonic losses under different operating conditions.To verify the accuracy and reliability of the model and the calculation method, a simple example model is built, as shown in Figure 20.The comparison between the measurement results and the Matlab simulation results about the harmonic loss of the line are reanalyzed is shown in Figure 4.
From Figure 21, the measurement results are basically consistent with the Matlab simulation results about the harmonic loss of the line.Due to the aging of the line and the loss of measuring instruments, the measurement results are larger than the Matlab simulation results.To verify the accuracy and reliability of the model and the calculation method, a simple example model is built, as shown in Figure 20.approximately 0.2~1.6%.Taking one day as the statistical period, the branch loss of the station area with the power supply radius of 220 m is 66.7326 kWh, and the harmonic line loss is about 0.399 kWh.The average harmonic loss of the distribution network is about 0.6%.Although the value is small, due to the large loss base of the large-scale low-voltage distribution network, the power loss caused by harmonics is also quite large.Moreover, the harmonic loss power of the low-voltage distribution network is affected by transformer and line type, power supply radius, load rate and harmonic current content rate.Different distribution networks may have different harmonic losses under different operating conditions.To verify the accuracy and reliability of the model and the calculation method, a simple example model is built, as shown in Figure 20.The comparison between the measurement results and the Matlab simulation results about the harmonic loss of the line are reanalyzed is shown in Figure 4.
From Figure 21, the measurement results are basically consistent with the Matlab simulation results about the harmonic loss of the line.Due to the aging of the line and the loss of measuring instruments, the measurement results are larger than the Matlab simulation results.The comparison between the measurement results and the Matlab simulation results about the harmonic loss of the line are reanalyzed is shown in Figure 4.
From Figure 21, the measurement results are basically consistent with the Matlab simulation results about the harmonic loss of the line.Due to the aging of the line and the loss of measuring instruments, the measurement results are larger than the Matlab simulation results.

Conclusions
To systematically and quantitatively analyze the harmonic loss of the low-voltage distribution network, the effects of the dead-time effect, SPWM modulation and background harmonics on output harmonic characteristics of PV inverters are studied in this paper.Then, the harmonic loss factor of winding resistance is defined.The harmonic impedance modeling of key equipment such as lines and transformers and the calculation method of harmonic loss of the low voltage distribution network based on measured data are proposed.At last, a harmonic power flow calculation model of the distribution network is presented.The results show that the overall harmonic voltage of the system is greatly affected by PV capacity, especially the voltage of the bus near the PV access bus.In addition, the transformer copper loss and line loss increase with the increase of harmonic current content and harmonic order.It is also found that the harmonic loss of the low-voltage distribution network is mainly determined by the harmonic current, accounting for about 0.6% of the total network loss.Therefore, reasonable planning of PV access capacity and limiting the distortion and order of harmonic current in the low-voltage distribution network is of great significance to reduce network loss.
For dual-winding transformers, the equivalent impedance model can be used.The calculation formula of transformer harmonic impedance Z Th is: In power flow calculation, the π-type equivalent circuit can be used as an equivalent harmonic parameter circuit of the line.It is shown in Figure A1.In power flow calculation, the π-type equivalent circuit can be used as an equivalent harmonic parameter circuit of the line.It is shown in Figure A1.
where Zλh is the wave impedance of the line under hth harmonic; λh is the propagation constant of the line under hth harmonic; l is the length of the transmission line.Considering that the line length of the distribution network is short, when the harmonic order studied is small, the distribution parameter characteristics are not obvious.The approximate model can be used for analysis.The harmonic impedance ZL,h and ground admittance YL,h per unit length of the line are shown in Formula (A1).In harmonic power flow analysis, non-harmonic source load is expressed as an impedance model.If the series impedance model is adopted, the series impedance will be expressed as:

Figure 1 .
Figure 1.The simplified circuit of the solar cells.

Figure 1 .
Figure 1.The simplified circuit of the solar cells.

Figure 2 .
Figure 2. The structure of the PV grid-connected.

Figure 3 .
Figure 3.The steady-state vector diagram of the grid-connected inverter.

Figure 3 .
Figure 3.The steady-state vector diagram of the grid-connected inverter.

Figure 4 .
Figure 4.The harmonic equivalent resistance of the transformer.

Figure 4 .
Figure 4.The harmonic equivalent resistance of the transformer.

Figure 5 .
Figure 5.The flow chart of harmonic power flow calculation.

Figure 6 .of case 3 Figure 7 .
Figure 6.The block diagram of the Matlab/Simulink model.

Figure 6 .
Figure 6.The block diagram of the Matlab/Simulink model.

Figure 7 .
Figure 7. Relationship between the harmonic distortion of current and harmonic loss.

Figure 8 .
Figure 8. Relationship between the harmonic distortion and harmonic order of current and harmonic loss.

Figure 8 .
Figure 8. Relationship between the harmonic distortion and harmonic order of current and harmonic loss.

Figure 10 .
Figure 10.Fundamental voltage variation with different PV capacity access.

Figure 11 .
Figure 11.Fundamental loss variation with different PV capacity access.

Figure 10 .
Figure 10.Fundamental voltage variation with different PV capacity access.

Figure 10 .
Figure 10.Fundamental voltage variation with different PV capacity access.

Figure 11 .
Figure 11.Fundamental loss variation with different PV capacity access.

Figure 12 .Figure 12 .
Figure 12.Harmonic voltage variation with different PV capacity access and harmonic order.

Figure 13 .
Figure 13.Harmonic loss variation with different PV capacity access and harmonic order.

Figure 13 .
Figure 13.Harmonic loss variation with different PV capacity access and harmonic order.

Figure 15 .
Figure 15.Daily trend of three-phase current.

Figure 14 .
Figure 14.The experimental schematic diagram of the measurement.

Figure 14 .
Figure 14.The experimental schematic diagram of the measurement.

Figure 15 .
Figure 15.Daily trend of three-phase current.

Figure 15 .
Figure 15.Daily trend of three-phase current.

Figure 17 .
Figure 17.The daily trend of current THD.

Figure 18 .
Figure 18.Equivalent circuit of distribution network.

Figure 17 .
Figure 17.The daily trend of current THD.

Figure 18 .
Figure 18.Equivalent circuit of distribution network.

Figure 17 .
Figure 17.The daily trend of current THD.

Figure 18 .
Figure 18.Equivalent circuit of distribution network.

Figure 18 .
Figure 18.Equivalent circuit of distribution network.

Figure 19 .
Figure 19.The time-varying trends of ratio of harmonic loss.

Figure 20 .
Figure 20.The simple example model.

Figure 19 .
Figure 19.The time-varying trends of ratio of harmonic loss.

Figure 19 .
Figure 19.The time-varying trends of ratio of harmonic loss.

Figure 20 .
Figure 20.The simple example model.

Figure 20 .
Figure 20.The simple example model.

Figure 21 .
Figure 21.The comparison between the measurement results and the Matlab simulation results.

2 Figure A1 .
Figure A1.Harmonic parameter model of distribution line.The parameters in the equivalent circuit are: xl,1 and bl,1 are reactance and admittance of line at the fundamental frequency.When the skin effect is considered, the harmonic resistance per unit length can be expressed as:

Table 1 .
The value range of each coefficient.

Table 4 .
Harmonic spectrum of PV generation.

Table 5 .
Parameters of PV generation.

Table 4 .
Harmonic spectrum of PV generation.

Table 5 .
Parameters of PV generation.