On Variable Reverse Power Flow-Part I : Active-Reactive Optimal Power Flow with Reactive Power of Wind Stations

It has recently been shown that using battery storage systems (BSSs) to provide reactive power provision in a medium-voltage (MV) active distribution network (ADN) with embedded wind stations (WSs) can lead to a huge amount of reverse power to an upstream transmission network (TN). However, unity power factors (PFs) of WSs were assumed in those studies to analyze the potential of BSSs. Therefore, in this paper (Part-I), we aim to further explore the pure reactive power potential of WSs (i.e., without BSSs) by investigating the issue of variable reverse power flow under different limits on PFs in an electricity market model. The main contributions of this work are summarized as follows: (1) Introducing the reactive power capability of WSs in the optimization model of the active-reactive optimal power flow (A-R-OPF) and highlighting the benefits/impacts under different limits on PFs. (2) Investigating the impacts of different agreements for variable reverse power flow on the operation of an ADN under different demand scenarios. (3) Derivation of the function of reactive energy losses in the grid with an equivalent-π circuit and comparing its value with active energy losses. (4) Balancing the energy curtailment of wind generation, active-reactive energy losses in the grid and active-reactive energy import-export by a meter-based method. In Part-II, the potential of the developed model is studied through analyzing an electricity market model and a 41-bus network with different locations of WSs.


Introduction
Buy-back is well-known in electricity markets where utilities or customers are buying or selling electric energy in a designed energy marketplace [1].However, this issue becomes more complex if renewable energies and/or storage systems are considered in connected power systems with bidirectional active and reactive power flows as seen in Figure 1 [2].
Besides the literature review given in [2], we mention some other recent studies apparent in the area of reactive power of renewable energies for power market clearing.It is to note that new researches have been recently focused on optimal reactive power flow in transmission or distribution systems or in between based on a more detailed model.In [3], e.g., the authors proposed an algorithm to minimize reactive power provision, transmission loss costs and transmission system voltage security margin by a four-stage multi-objective mathematical programming method to settle the reactive power market.The authors in [3] focused later in [4] on the stochastic reactive power market in the presence of volatility of wind power generation in a transmission system.This kind of stochastic reactive power from wind generators was also studied in [5] for a medium-voltage (MV) network.The works in [4,5] Energies 2016, 9, 121 2 of 12 used a model for reactive power production from wind turbines based on a detailed induction machine equivalent circuit as presented in [6].Other researchers [7], however, used a different model based on an electrical model of a wind turbine with a full scale converter.We note that reactive power contribution can also be modeled differently form photovoltaic [8] or energy storage systems [2].From another point of view, not only active and reactive powers from renewable energies are important but also the modeling of demand in hybrid electricity markets [9,10].Since considering "at once" all aspects and concerns in the power system shown in Figure 1 is prohibitive, our focus in this work will be on the power sub-system depicted in Figure 2.
Energies 2016, 9, 121 2 network.The works in [4,5] used a model for reactive power production from wind turbines based on a detailed induction machine equivalent circuit as presented in [6].Other researchers [7], however, used a different model based on an electrical model of a wind turbine with a full scale converter.We note that reactive power contribution can also be modeled differently form photovoltaic [8] or energy storage systems [2].From another point of view, not only active and reactive powers from renewable energies are important but also the modeling of demand in hybrid electricity markets [9,10].Since considering "at once" all aspects and concerns in the power system shown in Figure 1 is prohibitive, our focus in this work will be on the power sub-system depicted in Figure 2. Conceptual structure of a future power system [2].Here, S0 stands for the primary side while S1 for the secondary side of a transformer (TR) which is located between two different voltage levels.

Figure 2.
Illustration of the power system under consideration with a meter method [2].In this paper, the low-voltage (LV) network is considered with unidirectional power flows while the high-voltage (HV) network with bidirectional power flows.
As shown in Figure 2, at the high-voltage transmission network (HV-TN), firm [11] and non-firm [12] wind energy penetration can be maximized.It was mentioned in [13] that unit commitment issues and the stability of voltage at a HV-TN have high effects on wind energy losses.Based on [11][12][13], other important factors which influence wind energy losses in a HV-TN were studied in [14].It was shown that the reverse active power flow from a downstream medium-voltage active distribution network (MV-ADN) could be rejected at the HV-TN level.Based on the initial analysis in [14], it was assumed in [15] and [16] that all rates of flows for active and reactive power to the HV-TN could be accepted.However, in [15,16] wind stations (WSs) were assumed to operate with unity power factors (PFs) (i.e., Qdisp.w= 0, see Figure 2) either for simplicity or to show the pure impacts of battery storage systems (BSSs).
Renewable energy curtailment and reverse power flow can also occur in a low-voltage ADN (LV-ADN) with a high renewable penetration from distributed generation (DG) units such as photovoltaic (PV) systems [17][18][19].It was demonstrated in [19] that utilizing the reactive power capability of PV systems in a LV-ADN can lead to a huge amount of reactive energy import from an Figure 1.Conceptual structure of a future power system [2].Here, S 0 stands for the primary side while S 1 for the secondary side of a transformer (TR) which is located between two different voltage levels.
Energies 2016, 9, 121 2 network.The works in [4,5] used a model for reactive power production from wind turbines based on a detailed induction machine equivalent circuit as presented in [6].Other researchers [7], however, used a different model based on an electrical model of a wind turbine with a full scale converter.We note that reactive power contribution can also be modeled differently form photovoltaic [8] or energy storage systems [2].From another point of view, not only active and reactive powers from renewable energies are important but also the modeling of demand in hybrid electricity markets [9,10].Since considering "at once" all aspects and concerns in the power system shown in Figure 1 is prohibitive, our focus in this work will be on the power sub-system depicted in Figure 2. Conceptual structure of a future power system [2].Here, S0 stands for the primary side while S1 for the secondary side of a transformer (TR) which is located between two different voltage levels.

Figure 2.
Illustration of the power system under consideration with a meter method [2].In this paper, the low-voltage (LV) network is considered with unidirectional power flows while the high-voltage (HV) network with bidirectional power flows.
As shown in Figure 2, at the high-voltage transmission network (HV-TN), firm [11] and non-firm [12] wind energy penetration can be maximized.It was mentioned in [13] that unit commitment issues and the stability of voltage at a HV-TN have high effects on wind energy losses.Based on [11][12][13], other important factors which influence wind energy losses in a HV-TN were studied in [14].It was shown that the reverse active power flow from a downstream medium-voltage active distribution network (MV-ADN) could be rejected at the HV-TN level.Based on the initial analysis in [14], it was assumed in [15] and [16] that all rates of flows for active and reactive power to the HV-TN could be accepted.However, in [15,16] wind stations (WSs) were assumed to operate with unity power factors (PFs) (i.e., Qdisp.w= 0, see Figure 2) either for simplicity or to show the pure impacts of battery storage systems (BSSs).
Renewable energy curtailment and reverse power flow can also occur in a low-voltage ADN (LV-ADN) with a high renewable penetration from distributed generation (DG) units such as photovoltaic (PV) systems [17][18][19].It was demonstrated in [19] that utilizing the reactive power capability of PV systems in a LV-ADN can lead to a huge amount of reactive energy import from an Illustration of the power system under consideration with a meter method [2].In this paper, the low-voltage (LV) network is considered with unidirectional power flows while the high-voltage (HV) network with bidirectional power flows.
As shown in Figure 2, at the high-voltage transmission network (HV-TN), firm [11] and non-firm [12] wind energy penetration can be maximized.It was mentioned in [13] that unit commitment issues and the stability of voltage at a HV-TN have high effects on wind energy losses.Based on [11][12][13], other important factors which influence wind energy losses in a HV-TN were studied in [14].It was shown that the reverse active power flow from a downstream medium-voltage active distribution network (MV-ADN) could be rejected at the HV-TN level.Based on the initial analysis in [14], it was assumed in [15] and [16] that all rates of flows for active and reactive power to the HV-TN could be accepted.However, in [15,16] wind stations (WSs) were assumed to operate with unity power factors (PFs) (i.e., Q disp.w= 0, see Figure 2) either for simplicity or to show the pure impacts of battery storage systems (BSSs).
Renewable energy curtailment and reverse power flow can also occur in a low-voltage ADN (LV-ADN) with a high renewable penetration from distributed generation (DG) units such as Energies 2016, 9, 121 3 of 12 photovoltaic (PV) systems [17][18][19].It was demonstrated in [19] that utilizing the reactive power capability of PV systems in a LV-ADN can lead to a huge amount of reactive energy import from an upstream connected MV network.This can lead to considerable active energy losses in the LV-ADN.Furthermore, it was mentioned that the costs of reactive energy to be purchased from a MV network is much cheaper than the active energy from PV systems.However, a balance between active and reactive energy prices, especially under reverse power flows, was not considered in previous studies.
In comparison to PV systems, wind DG units were also studied in [20] to provide reactive power and support medium voltage networks under an incremental loss allocation method.Nevertheless, it was assumed in [20] that wind DG units are able to provide reactive support via extra installed devices such as capacitor banks.However, installing any extra devices for providing reactive power provision leads to extra payments and consequently higher investment costs.It should be mentioned here that the power capability of DG units to provide active/reactive power highly depends on the technologies as shown in [21,22].Briefly, utilizing the reactive power capability of DG units can significantly affect active energy losses in ADNs.
From another perspective, not only active energy losses in the grid are important in operating ADNs, but also the losses of reactive energy.This effect can be seen if a high price of reactive energy is considered, e.g., 90 $/Mvarh [23].Generally, energy prices can be assumed constant or being determined by an optimal contract price as in [24].However, considering optimal contract prices needs to know the quantity of power generation and the time of dispatch.Therefore, such models are restricted to dispatchable DG technologies.For non-dispatchable DG units, e.g., those based on renewable energy generation, feed-in tariffs are applicable in different ways and countries [25].
In summary, the issue of metering bidirectional power [26] due to possible reverse active and reactive power flows [27][28][29][30] in a transformer (TR) [31] which is typically used to connect different voltage levels in power systems, as seen in Figure 1, is of high interest and major importance for the society of power and energy systems.This represents a combined technical and economical complex problem, and, therefore, more studies are required in this area.The main contributions of this work are summarized as follows: ‚ Introducing the reactive power capability of WSs in the optimization model of the active-reactive optimal power flow (A-R-OPF) and highlighting the benefits/impacts under different limits on PFs.
‚ Investigating the impacts of different agreements and variable reverse power flow on the operation of an ADN under different demand scenarios.
‚ Derivation of the function of reactive energy losses in the grid with an equivalent-π circuit and comparing its value with active energy losses.
‚ Balancing the energy curtailment of wind generation, active/reactive energy losses in the grid and active/reactive energy import/export by a meter-based method.
The remainder of the paper is organized as follows: Section 2 describes the problem and introduces modeling procedures.In Section 3, the A-R-OPF model considering reactive power of WSs is developed.Section 4 concludes the paper.

Problem Statement and Modeling Procedure
Before formulating the A-R-OPF problem in a MV network, we would like to note that there are some major differences between optimal power flow problems in transmission and distribution systems.Such differences come from the fact that transmission systems are typically large and stretched on different areas with possibly dissimilar environmental issues.This leads to different system parameters and brings new technical and economic constraints in comparison to that in distribution systems.From another perspective, distribution systems could also include different entities which could be connected to balanced and/or unbalanced distribution networks.Therefore, a clear problem formulation should be given to distinguish such differences.For example, in Figure 1, if the focus is on the MV network, a model should be used to describe active/reactive power exchanges and voltage amplitudes/angles at the connection points, i.e., S 0 towards the LV network and S 1 towards the HV network.In addition, it should be clarified what exactly the MV network includes.
The authors in [15] analyzed the pure reactive power potential of BSSs in an ADN with WSs to maximize the revenues of wind power and BSSs and meanwhile minimize the costs of active energy losses.The WSs, however, were assumed to work with unity PFs.Therefore, it is aimed in this paper to further explore the reactive power potential of WSs "in the absence of BSSs" by considering variable reverse power flow and demand.In addition, we minimize the costs of both energy losses in the ADN, i.e., P loss and Q loss .For clarity, a modeling procedure of each entity related to the balanced power system depicted in Figure 2 is described in the following.

Wind Station
Generally, the capability of a power conditioning system (PCS) used for connecting DG units to a power system depends on DG technologies, e.g., see [21].In this paper, a full-scale power electronic converter [22] is modeled as an ideal PCS as seen in Figure 3. Furthermore, we incorporate a curtailment factor β c.w (see the nomenclature in Appendix 1) to ensure feasible solutions of the A-R-OPF as defined in [15] and [19].The relationships between the active P w and reactive Q disp.wpower in this model are simply explained using Figure 3 and expressed as: S PCS.w pl, hq " ˆpP w pl, hqβ c.w pl, hqq 2 `´Q disp.wpl, hq Q disp.ava.wpl, hq " ˘´pS PCS.max.wplqq 2 ´pP w pl, hqβ c.w pl, hqq 2 ¯0.5 4 variable reverse power flow and demand.In addition, we minimize the costs of both energy losses in the ADN, i.e., Ploss and Qloss.For clarity, a modeling procedure of each entity related to the balanced power system depicted in Figure 2 is described in the following.

Wind Station
Generally, the capability of a power conditioning system (PCS) used for connecting DG units to a power system depends on DG technologies, e.g., see [21].In this paper, a full-scale power electronic converter [22] is modeled as an ideal PCS as seen in Figure 3. Furthermore, we incorporate a curtailment factor βc.w (see the nomenclature in Appendix 1) to ensure feasible solutions of the A-R-OPF as defined in [15] and [19].The relationships between the active Pw and reactive Qdisp.wpower in this model are simply explained using Figure 3 and expressed as: Here, the model used for transforming wind speed to wind power is considered as in [15].The introduction of a curtailment factor (0 ≤ βc.w ≤ 1) at each WS is necessary to curtail the wind power generation to prevent violations of system constraints.More importantly, we consider the situation in which the available reactive power need to be restricted by different limits on PFs [21,22] (Figure 3a).This is represented with the dark-area in Figure 3b.Therefore, to study the impacts/benefits of such varying limits on PFs, additional constraints are introduced to the original model [19]: Here, PFmax.w= 1 represents the highest PF.In this study, we are concerned about the benefits/impacts of considering different values of PFmin.w.In other words, PFmin.w= 1 represents the original A-R-OPF model [15], otherwise specified.Here, the model used for transforming wind speed to wind power is considered as in [15].The introduction of a curtailment factor (0 ď β c.w ď 1) at each WS is necessary to curtail the wind power generation to prevent violations of system constraints.More importantly, we consider the situation in which the available reactive power need to be restricted by different limits on PFs [21,22] (Figure 3a).This is represented with the dark-area in Figure 3b.Therefore, to study the impacts/benefits of such varying limits on PFs, additional constraints are introduced to the original model [19]: Here, PF max.w = 1 represents the highest PF.In this study, we are concerned about the benefits/impacts of considering different values of PF min.w .In other words, PF min.w= 1 represents the original A-R-OPF model [15], otherwise specified.

Substation Transformer
Typically, distribution companies purchase active and reactive power from an external grid/market within the transfer limit of its substation TRs [23].As depicted in Figure 1, a downstream network is being supplied by an upstream network via a TR which allows the flow of powers in two directions.For the MV network, the forward direction means import while the reverse direction means export [26][27][28][29][30]. Basically, the apparent power at bus S 1 is described as: In addition, it is assumed that the TR is equipped with an on-load tap changer control system which is typically required to hold the voltage amplitude at bus S 1 at a specific level.This can be achieved if a suitable number of taps is available [31].The TR capacity is constrained in this paper by: S S1 phq ď S S1.max (7) ´αP1.revS S1.max ď P S1 phq ď α P1.fw S S1.max (8) ´αQ1.revS S1.max ď Q S1 phq ď α Q1.fw S S1.max (9) Note that in [2], (0 ď α P1.fw , α Q1.fw , α P1.rev , α Q1.rev ď 1) are considered as four decision variables with fixed values (i.e., α P1.fw = 1, α Q1.fw = 1, α P1.rev = 0.6, α Q1.rev = 0.6) to control the power flow rate at bus S 1 in two directions.However, using fixed values 0.6 (i.e., 60% of TR capacity) cannot reflect the whole effect of reverse power flow.The latter will possibly represent different agreements between the upstream and downstream networks shown in Figure 1.In addition, these decision variables represent the degree or level of freedom of transferring power at bus S 1 as depicted in Figure 4 where five possible operating states of an ADN are illustrated.
State 1: This case represents an isolated power system where no active-reactive power is allowed to exchange between the upstream and downstream network at bus S 1 .In this case, renewable/conventional generating units and/or BSSs should be able to satisfy the active-reactive demand and losses inside the downstream network.State 2: Only active power is allowed to be imported (purchased) from the upstream network to satisfy the demand and losses in the downstream network.In this case, reactive power should be generated locally by suitable means to satisfy the reactive demand and losses in the downstream network.State 3: This case can be considered as the conventional state where powers are purchased from the network under a high voltage level.The purchased power must satisfy the active/reactive demand and losses in the network under a low voltage level.Note that in all states 1-3 renewable/conventional generating units and/or BSSs can be integrated in the downstream network, but no reverse active-reactive power is allowed to be sold to the upstream network.State 4: Active/reactive power is allowed to be purchased from the upstream network to satisfy all or part of the active/reactive demand and losses in the downstream network, but only reverse active power is allowed to be sold to the upstream network.In comparison to the states 1-3, additional agreements between the two sides should be declared for allowing bidirectional active power flows at bus S 1 .
Energies 2016, 9, 121 6 of 12 State 5: Both active and reactive powers are allowed to exchange between the upstream and downstream network.Here, more agreements are required before allowing bidirectional active-reactive power flows at bus S 1 .

Demand
Based on the initial analysis in [14], active and reactive power demand from passive distribution networks (PDNs) and/or ADNs depend on possible internal optimal operations of both PDNs and ADNs.It means that the active and reactive power demand of the LV networks (see Figure 1) could be zero, positive and/or negative if a kind of power producer (e.g., renewable generators and/or storage systems) is embedded.In this work, for simplicity and clarity, both active and reactive demand power profiles are assumed to follow the Institute of Electrical and Electronics Engineers (IEEE) typical positive power demand in winter days [2].Such demand is referred to as the base demand scenario (100%).However, other different demand scenarios are also considered.Briefly, four demand scenarios, namely, very low (10%), low (50%), medium (100%) and high (150%) are considered and their impacts are analyzed.

Active-Reactive Optimal Power Flow with Reactive Power of Wind Stations
Typically, using different energy prices in any power marketplace affects the results significantly as shown in [16,19].Therefore, for clear analysis, fixed tariff price models for both active and reactive energies [15,16] are assumed in this work.

Objective Function
Based on Figure 2, the objective function of A-R-OPF considered in this paper is defined as Equation (10).The control variables in the proposed model of A-R-OPF are the curtailment factor of active power dispatch of WSs and the reactive power dispatch of WSs.It is aimed to maximize the revenue from the wind power Equation (11) and meanwhile minimizing the costs of both Equations ( 12) and ( 13) energy losses (see mathematical derivations in Appendix 2) in the MV network.In addition, we explicitly evaluate the costs of active and reactive energy (if being imported/purchased) in the Equations ( 14) and (15) or revenues (if being exported/sold) of only the active energy at bus S1.

Demand
Based on the initial analysis in [14], active and reactive power demand from passive distribution networks (PDNs) and/or ADNs depend on possible internal optimal operations of both PDNs and ADNs.It means that the active and reactive power demand of the LV networks (see Figure 1) could be zero, positive and/or negative if a kind of power producer (e.g., renewable generators and/or storage systems) is embedded.In this work, for simplicity and clarity, both active and reactive demand power profiles are assumed to follow the Institute of Electrical and Electronics Engineers (IEEE) typical positive power demand in winter days [2].Such demand is referred to as the base demand scenario (100%).However, other different demand scenarios are also considered.Briefly, four demand scenarios, namely, very low (10%), low (50%), medium (100%) and high (150%) are considered and their impacts are analyzed.

Active-Reactive Optimal Power Flow with Reactive Power of Wind Stations
Typically, using different energy prices in any power marketplace affects the results significantly as shown in [16,19].Therefore, for clear analysis, fixed tariff price models for both active and reactive energies [15,16] are assumed in this work.

Objective Function
Based on Figure 2, the objective function of A-R-OPF considered in this paper is defined as Equation (10).The control variables in the proposed model of A-R-OPF are the curtailment factor of active power dispatch of WSs and the reactive power dispatch of WSs.It is aimed to maximize the revenue from the wind power Equation (11) and meanwhile minimizing the costs of both Equations ( 12) and ( 13) energy losses (see mathematical derivations in Appendix 2) in the MV network.
In addition, we explicitly evaluate the costs of active and reactive energy (if being imported/purchased) in the Equations ( 14) and (15) or revenues (if being exported/sold) of only the active energy at bus S 1 .max where F 2 " F 5 "
In order to show the contribution of WSs as the sole reactive power source in the MV network, no BSS is considered in this paper in comparison to [2].More precisely, the MV network depicted in Figure 2 includes, as mathematically formulated in Equations ( 16) and (17), demand (P d , Q d ), WS (P w , Q disp.w ) and slack bus S 1 (P S1 , Q S1 ).

Inequality Equations
The inequality equations include the constraints of the voltage limits: and limits of distribution lines: Spi, j, hq ď S l.max pi, jq, i, j P Npi ‰ jq In addition, and distinguished from the original A-R-OPF [15,16], the upper limits of apparent power (3) and PFs of WSs ( 4) and ( 5) as well as upper limits of apparent power at bus S 1 ( 7)-( 9) are also included into the inequality constraints.It is to note, however, that stability constraints are not included in the above model and the formulated A-R-OPF problem is solved by the same general algebraic modeling system (GAMS) using the nonlinear programming (NLP) solver CONOPT3 as in [15,16].

Operating Conditions
Based on the market strategies introduced in [15,16,[23][24][25], we consider the following operating conditions: ‚ The operator of the MV network shown in Figure 2 is considered to be the sole owner of the network [24], who aims at a reliable operation of the system.
‚ The system operator aims to maximize the benefits from wind power and meanwhile to minimize the costs of both active and reactive energy losses [23].
‚ The imported/exported active energy from/to the upstream HV network, the active energy of WSs in the MV network and active energy losses in the MV network are calculated by a given price model for active energy, i.e., the fixed on-peak (100 $/MWh) and off-peak (50 $/MWh) tariff price [15,25].
‚ The imported/exported reactive energy and reactive energy losses in the MV network are calculated by a given price model, i.e., the fixed (12 $/Mvarh) tariff price [16].
‚ The reverse active energy to the HV network is permitted under different levels, while reverse reactive energy is not permitted for the reasons given in [32,33].

Questions
From the above extended version of the A-R-OPF model, Part-II [33] will answer the following questions: (1) What are the benefits/impacts of the "extended A-R-OPF model" for a real case power network?(2) What are the benefits/impacts of "varying PFs" of WSs in days with different levels of wind power generation?(3) How do the revenues/costs change with "variable reverse power flow" and "demand level" in an electricity market?(4) What are the relationships between the "location" of WSs, "variable reverse power flow" and "feeder congestion"?

Conclusions
In this paper, the optimization model of A-R-OPF in ADNs is further explored by considering the reactive power capability of WSs.It is aimed to evaluate the potential/effects of considering varying PFs in the A-R-OPF model.It is also aimed to minimize not only active but also reactive energy losses in grids.This is done by extending the objective function and considering both active and reactive energy prices.Furthermore, the impact of "long-to-short" and "short-to-long" feeders for power system networks is also investigated.The relationships and the interplay between wind generation curtailment, variable reverse active power flows, varying PFs of WSs, different demand levels and active/reactive energy prices in an electricity market model are shown in Part-II of this paper.The developed model in Part-I and investigated case studies in Part-II will help power system planers avoid extra size of devices (e.g., batteries) for reactive power provision if WSs are available in ADNs.Moreover, the relationships between the "location" of WSs, "variable reverse power flow" and "feeder congestion" are also studied in Part-II.

Appendix 2
In addition to the active power losses Equation (I) [15], reactive power losses can also be calculated for a given power network during hour h as Equation (II):

(II)
Furthermore, it is also aimed to continue the work in [32] by showing the impact of "long-to-short" and "short-to-long" feeders when calculating the imaginary part of the complex admittance matrixes in Equation (II).Therefore, a 4-bus illustrative network with a single phase equivalent-π circuit and one line type L1 is depicted in Figure A1, while a real case 41-bus network with different lines is provided in Part-II [33].Note that the length of lines or feeders can lead to convergence problems in power flow calculations in terms of "ill-conditioned" power systems [2].Here, R l , X l and B l are the line resistance, inductive reactance and capacitive susceptance, respectively [2].Let the value of base voltage equals (27.6 kV), the value of base power equals (10 MVA) and the network data in Table A1, it gives the following matrices Equations (III)-(V) in pu system.It is mentioning that B sh is a diagonal matrix and it will be equal to a zero matrix if B l is neglected in comparison to R l and X l .Assuming that the load at bus 4 equals 10 MW with PF = 0.85 lagging and the voltage amplitude at the slack bus No. 1 is taken 1.03 pu and zero angle, then the active/reactive power balance and apparent power flow are as given in Tables A2 and A3.Furthermore, it is also aimed to continue the work in [32] by showing the impact of "long-to-short" and "short-to-long" feeders when calculating the imaginary part of the complex admittance matrixes in Equation (II).Therefore, a 4-bus illustrative network with a single phase equivalent-π circuit and one line type L1 is depicted in Figure A1, while a real case 41-bus network with different lines is provided in Part-II [33].Note that the length of lines or feeders can lead to convergence problems in power flow calculations in terms of "ill-conditioned" power systems [2].Here, Rl, Xl and Bl are the line resistance, inductive reactance and capacitive susceptance, respectively [2].Let the value of base voltage equals (27.6 kV), the value of base power equals (10 MVA) and the network data in Table A1, it gives the following matrices Equations (III)-(V) in pu system.It is mentioning that Bsh is a diagonal matrix and it will be equal to a zero matrix if Bl is neglected in comparison to Rl and Xl.Assuming that the load at bus 4 equals 10 MW with PF = 0.85 lagging and the voltage amplitude at the slack bus No. 1 is taken 1.03 pu and zero angle, then the active/reactive power balance and apparent power flow are as given in Tables A2 and A3.Table A1.Data of the ilustrative 4-bus network shown in Figure 5.

Figure 1 .
Figure 1.Conceptual structure of a future power system[2].Here, S0 stands for the primary side while S1 for the secondary side of a transformer (TR) which is located between two different voltage levels.

Figure 1 .
Figure 1.Conceptual structure of a future power system[2].Here, S0 stands for the primary side while S1 for the secondary side of a transformer (TR) which is located between two different voltage levels.

Figure 2 .
Figure 2. Illustration of the power system under consideration with a meter method[2].In this paper, the low-voltage (LV) network is considered with unidirectional power flows while the high-voltage (HV) network with bidirectional power flows.

Figure 3 .
Figure 3. (a)Illustration of the behaviors of active/reactive powers through an ideal power conditioning system (PCS)[2]; and (b) the capability diagram of a PCS[2,22], where the dark-area stands for varying limits on power factors (PFs).

Figure 3 .
Figure 3. (a)Illustration of the behaviors of active/reactive powers through an ideal power conditioning system (PCS)[2]; and (b) the capability diagram of a PCS[2,22], where the dark-area stands for varying limits on power factors (PFs).

Figure 4 .
Figure 4. Conceptual illustration of the main five operating states.Here, the signs "+"/"-" mean import/export active/reactive power at a slack bus S1, respectively, as defined in[2].

Figure 4 .
Figure 4. Conceptual illustration of the main five operating states.Here, the signs "+"/"-" mean import/export active/reactive power at a slack bus S 1 , respectively, as defined in[2].

Figure A1 .
Figure A1.A single phase equivalent-π circuit [2] of a 4-bus illustative network.Here, series and shunt parts of reactive power loss are considered.

Figure
Figure A1.A single phase equivalent-π circuit [2] of a 4-bus illustative network.Here, series and shunt parts of reactive power loss are considered.