Integrated Power Flow and Short Circuit Calculation Method for Distribution Network with Inverter Based Distributed Generation

Power flow calculation and short circuit calculation are the basis of theoretical research for distribution network with inverter based distributed generation. The similarity of equivalent model for inverter based distributed generation during normal and fault conditions of distribution network and the differences between power flow and short circuit calculation are analyzed in this paper.Then an integrated power flow and short circuit calculation method for distribution network with inverter based distributed generation is proposed.The proposed method let the inverter based distributed generation be equivalent to Iθ bus, which makes it suitable to calculate the power flow of distribution network with a current limited inverter based distributed generation. And the low voltage ride through capability of inverter based distributed generation can be considered as well in this paper. Finally, some tests of power flow and short circuit current calculation are performed on a 33-bus distribution network. The calculated results from the proposed method in this paper are contrasted with those by the traditional method and the simulation method, whose results have verified the effectiveness of the integrated method suggested in this paper.


Introduction
The reserved finiteness and destruction of environment by fossil energy have made the application of distributed generation (DG) widespread gradually [1].The grid-connected DG in distribution network can reduce the loss of fossil energy and network lines and improve the utilization of energy and the reliability of customers.However, DG also brings some unfavorable effects at the same time.For example, in normal operation condition of distribution network, the electrical states of distribution network will be changed due to the random power output from DG.And the changed installed position and installed capacity of DG may bring some unfavorable effects on voltage and loss.In fault operation condition of distribution network, the magnitude and direction of short circuit current flowing through the line could be changed too, and the misoperation and failure action of traditional 3-sectional current protection could happen [2,3].Power flow and short circuit calculation are the basis of analyzing the normal and fault conditions of distribution network, respectively.Therefore, it is important to study the power flow and short circuit calculation for research of distribution network with DG.
The DGs in distribution network are mainly divided into generator based DG and inverter based DG, of which the inverter based DG is applied in distribution network widely, such as photovoltaic (PV), permanent magnet synchronous generator (PMSG), and energy storage system (ESS).Therefore, it is significant to study the power flow and short circuit calculation methods for distribution network with inverter based DGs.As the research focus in this paper, the DG mentioned below without additional instruction referred to is expressed as inverter based DG.
Both in power flow and in short circuit calculation, the research at present mainly consists of building equivalent model of inverter based DG and putting forward corresponding calculation method.Being unlike generator based DG, the output characteristic of inverter based DG is mainly dominated by its control unit.In power flow calculation, the equivalent models of DGs under different control are 2 Mathematical Problems in Engineering different.In constant power control, the DG should be equivalent to  bus while in constant voltage control, the DG should be equivalent to PV bus if the reactive power output from DG is within limit.If the reactive power output from DG is beyond the limit, the DG should be equivalent to  bus [4,5].Therefore, the equivalent model of inverter based DG is mainly modeled as  bus or PV bus.And for calculation method, back/forward sweep method and Newton method are generally used the same as traditional power flow calculation [6][7][8].
In short circuit calculation, the present fault equivalent models of inverter based DG have had the following three forms.(1) At the assumption of constant active and reactive power output from DG during a fault, DG should be equivalent to constant power model [9].(2) According to the current limiting characteristic of inverter, DG should enter current limiting status during a fault.As a result, DG can be equivalent to constant current source approximately [10,11].(3) The fault equivalent model of inverter based DG is decided by judging the current limiting status of DG.The DG entering the current limiting status is equivalent to constant current model.Conversely, the DG is equivalent to constant power model [12,13].As far as short circuit calculation method is concerned, the superposition theorem, which calculates the normal component and fault component of distribution network, respectively, and superposes them [14], is usually adopted.
From the researches introduced above, it can be seen that, in power flow calculation, the equivalent model of inverter based DG is processed ideally.The current limiting status of inverter based DG has not been considered.For example, DG may enter into the current limiting status when the voltage stability of distribution network is analyzed with power flow calculation for the loads increasing too much.In this case, the control objectives of inverter cannot be maintained.Hence, if the inverter based DG is still modeled as  or PV bus, the impact of DG on power flow cannot be correctly included.Similarly, in short circuit calculation, DG may enter into the current limiting status as well for the decreased voltage at point of common coupling (PCC) during a fault.At this moment, the inverter based DG should be equivalent to current source, which has been mentioned in some researches already.But the current source models proposed in those researches are ideally for the simplified analysis of current output regulation.In addition, a lot of inverter based DGs have been configured with Low Voltage Ride Through (LVRT) control at present [15].Therefore, the fault equivalent model of inverter based DG under LVRT control should be studied.
In this paper, the equivalent model of inverter based DG in normal and fault operation conditions of distribution network and the difference between power flow and short circuit calculation are analyzed at first.Then, an integrated power flow and short circuit calculation method of distribution network with inverter based DG is proposed according to the similarities between equivalent models of inverter based DG in normal and fault operation conditions.Finally, the proposed method is verified in an example.It is convenient to use the integrated calculation method to calculate  power flow or short circuit current in practical engineering.

Equivalent Model of Inverter Based DG
The output characteristics of inverter based DG during normal and fault operation conditions of distribution network are associated with the control strategy of inverter.The dualloop control is usually adopted, which outer loop control targets are usually constant active and reactive power or constant DC voltage and reactive power.The output signals from outer loop are reference signals of inner ring, which mainly adopt - decoupled control [16,17].The control principle of inverter with current limiting function, which orients -axis to grid voltage vector, is shown in Figure 1.
In Figure 1,  ref ,  ref , and  dc.ref are the references of active power, reactive power, and DC voltage, respectively., , and  dc are the actual values of active power, reactive power, and DC voltage, respectively.  is PCC voltage of DG.   ,   are the - components of output current of inverter, respectively. ref ,  ref are the references of   and   , respectively. isangular frequency and  is filter inductance of AC side.
Thus, whatever the outer loop control target is, the inner loop control for inverter is current control invariably.According to the instantaneous power theory, the active and reactive power exchange between DG and grid can be calculated with [18]  =     ,  =     . ( Inverter based DG usually works under constant power control in normal operation condition.In this case, it can be equivalent to  bus because the active and reactive power output from DG is known.However, the distribution network with heavy loads will reduce network voltage.At this moment, the -axis and -axis current will increase according to (1), and the output current maximum  max from DG is usually less than 1.2∼1.5 times of its rated current according to its maximum allowable current threshold.After inverter enters into current limiting status, the output current magnitude of DG will be a constant.Along with the decreasing voltage, the output power from DG cannot follow the constant power control, which means that the DG cannot be equivalent to  bus any more.In fault operation condition of distribution network, the DG may enter into the current limiting status in the same way easier due to the decreased PCC voltage caused by fault.
After inverter enters into current limiting status, the magnitude of output current from DG is  max .In unit power factor control mode, the phase difference  between output current and PCC voltage of DG is 0 ∘ , while in nonunit power factor control mode, the power factor is given normally, which means the phase difference  is known.Accordingly, DG can be equivalent to  bus at current limiting status. is  max and  is the phase difference between the output current and PCC voltage of DG.Hence, the equivalent model of inverter based DG should be changed from  bus to  bus when inverter enters the current limiting status in both power flow and short circuit calculation.

The Difference between Short Circuit Calculation and Power Flow Calculation
Through the analysis of equivalent models for inverter based DG during normal and fault conditions of distribution network, it can be seen that the similarity between them is that DG should be modeled as  and  buses, respectively, before and after inverter entering current limiting status.
Hence, an integrated power flow and short circuit calculation method is proposed in this paper, which takes the idea of power flow calculation into short circuit calculation.This method can calculate both power flow and short circuit current of distribution network.Since the power flow calculation method of this paper is only applicable to symmetric system, the integrated calculation method based on it is only applicable to symmetrical fault calculation.In addition, short circuit calculation itself is different from power flow calculation whereby the following two problems have to be solved in short circuit calculation.

The Modification of Bus Admittance Matrix for Distribution Network.
In three-phase short circuit calculation, the node admittance matrix should be modified by following three reasons.
(1) Loads Model.The loads of distribution network are omitted in three-phase symmetrical short circuit calculation usually.If the impacts of loads on short circuit current have to be taken into consideration, the loads model should be changed from constant power model to constant impedance model.This is because the calculation would not converge due to the unbalanced power flows in downstream distribution network of fault point.Then, the bus admittance matrix of distribution network should be modified after the loads model is changed.The short circuit calculation in this paper has ignored the loads of distribution network.
(2) Grounding Resistance.Since fault point only exists in fault condition of distribution network, the circuit structures of distribution network during normal and fault conditions are different.In fault condition, the fault happens in somewhere of distribution network, which means that a grounded branch is connecting to the fault point, as shown in Figure 2.
From Figure 2 it can be seen that there is a new added grounding admittance at fault point .So that, in short circuit calculation, the self-admittance of fault point  needs to be modified according to the following equation: where  . and  . are the self-admittance of bus  of postfault and prefault distribution network, respectively.
(3) DG in Downstream of the Fault Point.In distribution network, the fault may occur in any position of line.If there are multiple DGs in same distribution network, the DGs can be divided into two categories which are DGs in upstream of fault point and DGs in downstream of fault point.In past short circuit calculation method of distribution network with DG, the short circuit current contribution from DG in downstream of fault point has not been taken into account generally [19].This is because the power capacity of DG in distribution network is small.Hence, the DG in downstream of fault point would be disconnected from distribution network under low voltage protection.Although the DG still connects to distribution network, the short circuit current from DG is small enough, thereby only the short circuit current contribution from DG in upstream of fault point is usually considered.However, short circuit current from the DGs in downstream of fault point is larger, if more DGs are installed in that area.In addition, DG is requested to keep grid-connection for a period of time during a fault.Therefore, the short circuit current contribution from DGs in downstream of fault point should not be ignored.
For example, the DG in Figure 2 is located in downstream of fault point, and an isolated circuit is formed between fault point and DG.At this moment, the output current from DG fails to follow its control target, since the output characteristics of DG are only associated with line impedance characteristics.To simplify the calculation, the DG in downstream of fault point should be regarded as current source and modeled as  bus with  =  max ,  = 0 ∘ .Then the parameters of the lines between DG and fault point should be updated from complex impedance to impedance modulus value  ( = √  2 +  2 ).And the bus admittance matrix of distribution network needs to be modified.

LVRT Fault Equivalent Model of DG.
Both in normal and in fault operation conditions of distribution network, the equivalent model of inverter based DG mentioned above is analyzed based on the assumption of DG under constant power.However, in order to make full use of reactive power capacity of inverter when a fault takes place, the DG has usually had the ability of LVRT [20].The fault equivalent model of DG with LVRT control is more complex because the reactive power outputs from DG and  are variable.
It is important to note that the output characteristics of DG in downstream of fault point are only associated with line impedance characteristics, so that the LVRT model can only be applied to the DG in upstream of fault point.In order to take the short circuit current contribution from DG with LVRT control into consideration, it is necessary to analyze the LVRT equivalent model of inverter based DG.
Although the LVRT standards of different types of inverter based DG are different, the control idea of them is the same.When a fault occurs in distribution network, whether the DG adopts the LVRT control or not is decided by PCC voltage detected by control system.In LVRT control, the reactive current reference value of control system is adjusted according to the PCC voltage.Then the reactive power output from DG can be adjusted too.In this paper, the LVRT fault equivalent model of DG is described as follows.
(1)   ≥ 0.9 .In this situation, DG still maintains original constant power control, and the LVRT control cannot be switched on.The fault equivalent model is where  DG is the PCC voltage of DG,  DG and  DG are the active power and reactive power output of DG, respectively, and  DG0 and  DG0 are the prefault active and reactive power output of DG.
(2) 0.4  <   < 0.9 .In this situation, the PCC voltage of DG drops below 0.9 pu and DG starts to output reactive power current according to expression as follows: where   is the per unit rated current of inverter based DG.  is the voltage adjustment coefficient of reactive power current output from DG.For example,  = 2 means the voltage dropped 1% based on 1 pu; the reactive power current output from DG increased 2% of rated current.With an increase in reactive current output from DG, the DG may enter the current limiting status due to the maximum allowable current threshold.The judgment of current limiting status is the key to DG converting to  bus from  bus.The current limiting status can be determined by the following equation: If ( 5) is not satisfied, the DG is not current limited.DG can still be equivalent to  bus.The fault equivalent model is If ( 5) is satisfied, the DG is current limited.DG should be converted into  bus.The fault equivalent model is (3)   ≤ 0.4 .In this situation, the PCC voltage of DG drops below 0.4 pu.In order to support voltage for distribution network adequately, DG should output reactive power current only.At this moment, the phase difference  between the output current and PCC voltage of DG is 90 ∘ .The equivalent model of the DG can be modeled as  bus directly.The fault equivalent model is From the LVRT equivalent model of DG, it can be seen that it is related to the PCC voltage of DG in current limiting status and noncurrent limiting status.Besides, whether the LVRT control is switched on or not, the fault equivalent models of DG are similar.Both the fault equivalent models can be modeled as  and  buses, respectively, before and after inverter entering current limiting status.Therefore, the idea of power flow calculation can be used to short circuit calculation under LVRT control similarly.

Integrated Calculation Method
According to the analysis of equivalent model for inverter based DG during normal and fault operation conditions of distribution network, it can be seen that DG should have been converted into  bus in some operation scenarios.The integrated calculation method proposed in this paper uses the idea of power flow calculation to calculate short circuit current.However, the power flow calculation method in the past cannot calculate the power flow of distribution network with  bus directly.Hence, in this section, the conversion condition of bus type is analyzed and the power flow calculation method of distribution network with  bus is studied.

Switching Condition of Bus Type.
At the beginning of calculation, the present status of DG is unknown.Therefore, it should be assumed that when DG keeps the constant power control at present, the output current from DG is calculated according to results of past traditional power flow calculation method by (9) where SDG. is the set value of output complex power from DG at bus , under  control.UDG. is the PCC voltage of DG.If  DG. >  max. , the DG at bus  enters into the current limiting status and should be switched to . depends on the output power factor of DG.

Power Flow Calculation Method of Distribution Network
with  Bus.For  bus in distribution network, only its injected current  and the phase difference  between  and PCC voltage are known.It is important to note that DG can be modeled as  bus with   = 0 while  = 0 ∘ and can be modeled as  bus with   = 0 while  = ±90 ∘ .But these two bus types can only be applied to the special cases that  = 0 ∘ or  = ±90 ∘ .In most cases, the current limited DG can be modeled as  bus, and the bus power equation of  bus can be expressed as where  ∈  bus,  is the voltage amplitude,   (  −   ) is the phase difference between bus  and bus , and   and   are the real and imaginary component of bus admittance matrix.
Let  be the number of buses and  be the number of  buses; hence, the number of  buses is  −  − 1.The simultaneous equations of power balance equations for  and  buses are obtained as follows: where in Δ 1 () and Δ 1 (),  ∈ 1, 2, . . .,  while in Δ 2 () and Δ 2 ()  ∈  + 1,  + 2, . . .,  − 1.In other words, subscript 1 represents  bus and subscript 2 represents  bus.From (11), it can be seen that there are 2( − 1) equations, and the number of unknowns is also 2(−1) including voltage amplitude and phase of each bus except slack bus.Therefore, the equations number is equal to the number of unknowns and the power flow calculation of distribution network with  bus can be solved.Equation ( 11) is a nonlinear equation which should be solved, so that the common Newton method is adopted in this paper.
The key problem to solving power flow by Newton method is the formation and computation of correction equations.According to solving process of Newton method, the correction equations are formed as follows: The emphasis for solving the correction equations is the formation of Jacobian matrix .Due to the emerging  bus in power flow calculation, the Jacobian matrix of past power flow calculation method cannot be used directly.From (12), it can be found that the Jacobian matrix is composed of four submatrixes which are , , , .According to the theory of Newton method, the expressions of diagonal elements in  and  are different from traditional power flow calculation method due to the introduction of  bus.The specific expressions are shown as follows.
Diagonal elements of  are ( The initial values of voltage amplitude and phase of each bus should be given, respectively, after Jacobian matrix  is formed, which are According to the initial values, the error vectors 1 , and Δ (0) 2 can be obtained.Then the next voltage amplitude and phase can be calculated according to (15) and ( 16) as follows: Then substituting ( 16) into (15) to continue the iterative calculation until max |Δ ()  1 , Δ () 2 | <  and max |Δ () 1 , Δ ()  2 | < , the voltage amplitude and phase of each bus can be obtained.Moreover, the output power of current limited inverter based DG can be solved by

Flowchart of Integrated Power Flow and Short Circuit
Calculation Method.In this paper, there are two points which must be noted in the proposed integrated power flow and short circuit calculation method.First, there are some differences between short circuit and power flow calculation, so that the bus admittance matrix of distribution network and the corresponding fault equivalent model of DG should be modified for short circuit calculation.Second, if there are multiple DGs in same distribution network, the equivalent models of the rest of DGs may be affected by the changed voltage due to the equivalent model of one DG changed from  bus to  bus.For this reason, the iterative process should be terminated until the equivalent model of each DG does not change any more.The specific flow chart of integrated power flow and short circuit calculation method is shown in Figure 3.

Examples and Analysis
In this paper, the numerical calculation is performed on a 33-bus distribution network, whose rated voltage is 10 kV.It is assumed that the same inverter based DGs are connected to buses 11, 17, 28, and 33, respectively, as shown in Figure 4, whose capacity is 0.5 MW.The per unit rated current of DG is 5 pu on the basis of the fact that base power and base voltage are 100 kVA and 10 kV.It is assumed that the DGs work underrated condition with unit power factor mode and their output current maximum  max is 1.2 times of rated current that is 6 pu.Bus 1 is slack bus whose voltage amplitude is 1 pu and voltage phase angle is 0 ∘ .Then power flow calculation and short circuit calculation in this test system should be carried out by using the integrated calculation method proposed in this paper to verify its effectiveness.

Power Flow Calculation.
In order to obtain current limiting DG, the loads in distribution network should be increased to 1.8 times of original value to obtain lower voltage in power flow calculation.According to calculation process, the calculation starts with DG equivalent to  bus and the current limiting status of DG is estimated.The RMS of PCC voltage and output current of each DG are obtained as shown in Table 1.In this paper, all calculation results are expressed in per unit (pu) value.
From the calculation results in Table 1, it can be seen that only DG4 enters current limiting status.Accordingly, DG4 cannot maintain constant power output and should be equivalent to  bus.Then, according to the calculation process, the RMS of PCC voltage and output current of each DG are obtained by the power flow calculation method with  bus until the equivalent models of each DG are no longer changed.The calculation results are shown in Table 2 and the voltage profile of distribution network is shown in Figure 5.In order to verify if the power flow of distribution network with  bus can be correctly calculated by the integrated calculation method proposed in this paper, a method based on calculation results comparison between  the power flow with the calculated  values.The correctness of proposed integrated method can be verified if the power flows calculated by two methods are the same.In this case, the power flow of distribution network and active power output from current limited DG4 are calculated by integrated method at first.And the calculated active power output from DG4 is 4.9146 pu.Then the power flow is again calculated by traditional power flow calculation method with DG4 modeled as  bus, whose active power is  = 4.9146 pu.Lastly, the maximum absolute value error of voltage of each bus calculated by two methods is 7.4884 − 6.So that the proposed integrated method in this paper is correct.

Short Circuit Calculation for DG without LVRT Control.
Similarly, the short circuit current is calculated by the integrated calculation method proposed in this paper in the same 33-bus distribution network.The three-phase grounding faults with grounding resistance are 0, which are set on buses 9, 18, and 31, respectively, as shown in Figure 6.The fault positions set in this paper can be divided into the following four cases.All DGs in distribution network are in upstream of fault point (f2), some DGs are in downstream of fault point (f1 and f3), the fault occurs in trunk line (f1 and f2), and the fault occurs in branch line (f3).According to the calculating process, the DG in upstream of fault point is equivalent to  bus and the DG in downstream of fault point is equivalent to  bus at first calculation.The output current from DGs is obtained as shown in Table 3.
From the calculation results, it can be seen that the DG in distribution network is easy to enter current limiting status during faults.If the DG in upstream of the fault point enters the current limiting status, its fault equivalent model should be converted into  bus.In this case, the output current from current limited DG is  max and  = 0 ∘ for the DG without LVRT control.Then the calculation process is continued until the equivalent models of each DG are no longer changed.The short circuit currents of f1, f2, and f3 are calculated.Moreover, the simulation of this example is built in MATLAB/Simulink.The calculation results by integrated calculation method are compared with those by the simulation method in MATLAB/Simulink, with their results shown in Table 4.
From the data in Table 4, it can be seen that the results based on the integrated calculation method proposed in this paper are similar to the simulation results and the errors between them are reasonable, whereby the effectiveness and correctness of the proposed method are verified.Besides, the calculation results of the part of bus voltage of distribution network are shown in Table 5.

Short Circuit Calculation for DG with LVRT Control.
In this case, the DG in upstream of fault point is equivalent to  bus and the DG in downstream of fault point is equivalent to  bus similarly at first.The calculated current outputs from DGs are the same as those in Table 3.Then the current limiting status of DG in upstream of fault point is judged according to the calculation results.According to the LVRT equivalent model, the output current and the phase difference between output current and PCC voltage of DGs in current limiting status are calculated, and also, the output active and reactive power of DGs in noncurrent limiting status are calculated.Then the calculation process is continued until the equivalent models of each DG are no longer changed.In this paper, the short circuit currents of three kinds of fault are calculated according to the integrated calculation method.Moreover, the simulation of this example is built by MATLAB/Simulink.The calculation results by integrated calculation method are compared with those by the simulation method in MATLAB/Simulink, with their results shown in Table 6.
From the data in Table 6, it can be seen that the results based on the integrated calculation method are similar to the simulation results and the errors between them are reasonable.In addition, the calculation results of the part of bus voltage of distribution network are shown in Table 7.
From the short circuit calculation results in Tables 5  and 7, it can be seen that the bus voltages by short circuit calculation of DG with LVRT control are bigger than those without LVRT control.It is indicated that the LVRT control has had an increasing effect upon bus voltage, which reflects its advantage in the case of distribution network fault.

Conclusion
According to the theoretical analysis and simulation verification for the integrated power flow and short circuit calculation method of distribution network with inverter based DG, the following conclusions can be obtained: (1) The equivalent models for inverter based DG during normal and fault conditions of distribution network are similar to each other.Both equivalent models can be modeled as  and  buses, respectively, before and after inverter entering current limiting status.Therefore, the power flow and short circuit current can be calculated in the same way by an integrated power flow and short circuit calculation method.(2) If distribution network has  bus, the expressions of elements related to  bus in Jacobian matrix being a part of correction equation for power flow calculation are different from those by traditional power flow calculation.The expressions should be reformed according to the power balance equations of  bus.Then power flow calculation by solving the correction equations can be completed.(3) LVRT control is a conventional control of DG when a fault takes place in a distribution network.Hence, the LVRT equivalent model, which is more complicated than equivalent model in normal operation condition, should be analyzed.(4) The distinction and similarity of power flow calculation and short circuit calculation should be attended for using the idea of power flow calculation to calculate short circuit current.And the adaptive complements for short circuit calculation are added in the process of integrated calculation method.The integrated power flow and short circuit calculation method proposed in this paper can be directly applied to practical engineering applications.
as follows:

Table 1 :
PCC voltage and output current of DG.

Table 2 :
PCC voltage and output current of DG.Because the  values of DG under current limiting status are unknown, an inversion method is used.Concretely, the  values of DG under current limiting status and the power flow of distribution network are calculated by the power flow method with  bus firstly.Then the traditional calculation method is used to calculate

Table 3 :
Short circuit current output from each DG.

Table 4 :
Short circuit current comparison results between calculation and simulation.

Table 5 :
Part of bus voltage of distribution network.

Table 6 :
Short circuit current comparison results between calculation and simulation.

Table 7 :
Part of bus voltage of distribution network.