Active and Reactive Power Collaborative Optimization for Active Distribution Networks Considering Bi-Directional V2G Behavior

Due to their great potential for energy conservation and emission reduction, electric vehicles (EVs) have attracted the attention of governments around the world and become more popular. However, the high penetration rate of EVs has brought great challenges to the operation of the Active Distribution Network (ADN). On the other hand, EVs will be equipped with more intelligent chargers in the future, which supports the EVs’ high flexibility in both active and reactive power control. In this paper, a distributed optimization model of ADN is proposed by employing the collaborative active and reactive power control capability of EVs. Firstly, the preference of EV users is taken into account and the charging mode of EVs is divided into three categories: rated power charging, non-discharging, and flexible charging–discharging. Then, the reactive power compensation capacity of the plugged-in EV is deduced based on the circuit topology of the intelligent charger and the active–reactive power control model of the EV is established subsequently. Secondly, considering the operation constraints of ADN and the charging–discharging constraints of EVs over the operation planning horizon, the optimization objective of the model is proposed, which consists of two parts: “minimizing energy cost” and “improving voltage profile”. Finally, a distributed solution method is proposed based on the Alternating Direction Method of Multipliers (ADMM). The proposed model is implemented on a 33-bus ADN. The obtained results demonstrate that it is beneficial to achieve lower energy cost and increase the voltage profile of the ADN. In addition, the energy demand of EV batteries in their plugin intervals is met, and the demand preference of EV users is guaranteed.


Introduction
With the integration of a large number of distributed generations (DGs), the traditional distribution network (DN) has gradually been transformed to the Active Distribution Network (ADN) [1], which includes many controllable resources. Another aspect worthy of attention is that the large-scale penetration of electric vehicles (EVs) poses new challenges to the operation of ADN. It is anticipated that by 2050, about one hundred million EVs will be produced annually [2]. ADN improves the flexibility and controllability of the DN, and also brings new problems, the most prominent of which is the abnormal voltage fluctuation and congestion in feeders. Due to the development of V2G technology, the bi-directional charging topology capable of four-quadrant operation is gradually becoming the mainstream topology of charging stations [3]. Based on the gradual maturation of V2G technology, coordinating the active and reactive power of EVs has great potential for easing the abnormal fluctuations and feeder congestion of the ADN [4].
(1) In practice, driving is the primary attribute of EVs. As highly autonomous individuals, the preferences of EV users cannot be ignored. (2) The reactive power compensation capability of EVs is not only restricted by the active power and capacity of the charger, but also needs to consider the constraints of the power electronic equipment of the charger. (3) Most of the existing studies are based on the general DN, but ADN contains many DGs, wind power, photovoltaic, reactive power compensation devices, etc., and when the scale of plugged-in EVs and the demand preference of EV users are taken into account, the current research findings are inevitably not applicable to this new scenario.
In this paper, we propose a novel decentralized framework for collaborative optimization of active and reactive power in the ADN with the participation of the EVA. dThe simplified framework is depicted in Figure 1. We assume that the EVA is allowed to obtain EVs' behavioral characteristics, including battery parameters (initial SOC, rated charging power, and rated discharging power) and plug-in and plug-out time. EV users determine the charging mode by interacting with the EVA according to their preference. Afterward, the EVA determines the adjustment range of active and reactive power, according to users' demands. From the point of view of the ADN, the EVA is also a controllable and flexible resource. The ADNO takes the minimum energy cost and the minimum deviation of the voltage amplitude as the goal to coordinate the active and reactive power of EVA under relevant constraints. The main contributions of this paper can be summarized as follows: (1) Considering the users' preference, the plugged-in EVs are subdivided into three categories, and then the EV control model is established separately. The establishment of the differentiated EV control model is more in line with the actual scenario. (2) Based on the topological constraints of the charging pile, the reactive power compensation capacity of the EV is modeled to derive a more accurate range of reactive power. (3) The proposed model is solved by a decentralized algorithm, which is developed based on the Alternating Direction Method of Multipliers (ADMM). The operating problems of the ADNO and EVA are coordinated and optimized independently.

Modeling Collaborative Active-Reactive Power of EV Cluster
Firstly, we define that the EVs belonging to the same EVA are an EV cluster (EVC). The number of EVC is the same as that of EVA. In this section, the active-reactive power model of a single EV is firstly established, and then the active-reactive power model of each EVC is proposed.
In order to obtain the active-reactive power control model of a single EV, the active power control and reactive power control of a single EV are discussed, respectively, and then the collaborative active-reactive power model is established.
The preference of the EV user is mainly reflected in the difference in the active power. According to the charging or discharging state of the EV battery, the plugged-in EVs are divided into three categories: rated power charging EVs, non-discharging EVs, and flexible charging-discharging EVs. Three different categories reflect the preference of EV users: when the user selects the rated power charging mode, the user expects to minimize the time cost; when the non-discharging mode is selected, it indicates that the user expects to reduce the cost of money without increasing extra battery loss; and when the flexible charging-discharging mode is selected, it means that the user expects the money cost to be minimized. In practice, the user can interact with EVA to determine the charging mode when the EV is plugged into the ADN. The number of EVA in the ADN is denoted as N A , and D i,j (i = 1, 2, 3; j = 1, 2, 3 . . . N A ) represents the set of EVs belonging to the ith category in the jth EVA.
The three categories of plugged-in EVs are modeled separately as follows: ∀l ∈ D 1,j (j = 1, 2, · · · , N A ), the control model of l is formulated by Equations (1) and (2):    P l,t = min(P cha , P r c,l ), t ∈ t l,in and S l,t ∈ [S l,or , S l,ex ) 0 ≤ P l,t ≤ min(P cha , P r c,l ), S l,t ≥ S l,ex P l,t = 0 , t / ∈ t l,in S l,t = S l,or + ∑ η c,l ·P r c,l ·∆t/E l (2) where P cha , P r c,l , η c,l , S l,or , S l,ex and E l are the maximum active power, rated charging power, charging efficiency, initial state of charge (SOC), expected SOC of departure and maximum capacity of the battery, respectively. P l,t and S l,t are the actual power and actual SOC of l at timeslot t. t l,in is the time span of the arrival time to departure time of l and ∆t is the unit interval.
The reactive power operation of the EV is based on the smart charging pile, which has gradually become the focus of attention because of its capability for four-quadrant operation [20]. The representative topology is shown in Figure 2.
As illustrated in Figure 2, the intelligent charger enables a bi-directional flow of active and reactive power between the ADN and EV charger. However, there is only bi-directional active power between the battery and EV charger, which means that the reactive power operation will not affect the battery lifetime.  For the convenience of analysis, the following reasonable assumptions are made: (1) The ADN voltage is an ideal sine wave with stable frequency (f = 50 Hz).
(2) The charger is normally connected to the ADN via cable, so to simplify the calculation, assume that the impedance between the intelligent bi-directional charger and the ADN is the inductive impedance; the sum of the impedance stacks is recorded as L c . (3) The power operation of the charger is realized through the smart control circuit of the charger. The influence of the control circuit on the voltage of the charger is ignored.
Based on the above assumptions, the simplified circuit diagram is shown in Figure 3. Assume that the grid voltage u c (t) is √ 2 sin(ωt), where V c is the effective value. Let δ be the angle that the phase angle of the charger voltage u s (t) lags the u c (t); then, u s (t) can be represented by Equation (7): The instantaneous power P s (t) of the charger can be derived: where I c is the effective value of the current. Equation (8) shows that the injected instantaneous power of the charger includes two parts: average power and ripple power. The average power is an active power used to charge or discharge the battery of EV, and the ripple power is a kind of oscillating power flowing between the charger and ADN, which is only temporarily stored in the capacitance of the charger. Obviously, the ripple power can characterize the charger's capacity for reactive power. Through mathematical deduction (the detailed deduction process is given in Appendix A), the maximum reactive power of the charger can be obtained: Meanwhile, the power of the charger can satisfy the constraints of Equations (11) and (12): The constraints of Equations (11) and (12) show that the reactive operation of the charger is not symmetrical. The limit field of the active and reactive operation of the charger is shown in Figure 4. The shaded area enclosed by the blue solid line in Figure 4 is the active and reactive operation domain of the charger. The labels 1 , 2 , 3 and 4 represent the extra absorbable reactive capacity, absorbed reactive capacity, released reactive capacity and extra releasable reactive capacity when P s = P m , respectively. It is worth noting that the operation domain is in ideal conditions. When the charger provides extra absorbable or releasable reactive power (e.g., 1 or 4 ), its active power is reduced to zero. In practice, the active power of some EVs is prohibited from being reduced to zero due to the preference of users and the energy demand of the EV battery.
We assume that P cha is equal to P r c,l in this paper. The active-reactive operation domain of EV, considering the preference of the user and the energy demand of the EV battery, is shown in Figure 5. P l,max , −P d,max are the minimum charging and maximum discharging active power when considering the preference of the user and the energy demand of the battery; Q l,max is the corresponding maximum reactive power. For any EV in D 1,j set, Equation (1) indicates that there is no reactive power release or absorption when S l,t < S l,ex . The charging power can be cut down when S l,t ≥ S l,ex . In other words, the EV in D 1,j set is equivalent to the EV in D 2,j set when S l,t ≥ S l,ex .
So, the maximum reactive power of the jth (j = 1, 2, . . . , N A ) EVA can be calculated by Equation (13): where Q j es,abs (t) and Q j es,rel (t) are the maximum absorbable reactive power and the maximum releasable reactive power of the jth (j = 1, 2, 3 . . . N A ) EVA at timeslot t, respectively. N j l,t is the EV number in D 1,j set whose SOC has exceeded the value of S l,ex at timeslot t. N j 2,t and N j 3,t are the number of EVs in D 2,j set and the number of EVs in D 3,j set at timeslot t, respectively. N j c,t and N j d,t are the actual quantity of charging and discharging EVs in D 3,j at timeslot t.
Obviously, the active power of the jth EVA is a linear superposition of the active power of each EV. The actual active power of the jth EVA is formulated by Equation (14): Sustainability 2021, 13, 6489 8 of 26 where N j 1,t is the EV number in D 1,j set; P q es (t) is the actual active power of jth EVA at timeslot t. So, the total power of EVs plugged into the ADN at timeslot t can be calculated by Equations (15) and (16).
where P es (t) and Q es (t) are the total active and reactive powers of EVs at timeslot t.

Problem Formulation
The objective of the ADNO is to minimize the total energy cost and minimize the summation of voltage deviation from the desired value. Without loss of generality, the objective function can be written as a weighted sum of two objectives: where C e,t , C g,t and C A,t are the total energy cost, fuel generations (FGs) cost and EV compensation cost. γ t and P pur ADN,t are the energy price and purchased power from upstream DN at timeslot t, respectively. v b,t , P L b,t and P G b,t are the voltage, conventional load and active power injected of node b, respectively. C on s,t , C f uel s,t and C o f f s,t are the start-up cost, fuel cost and shut-down cost of the sth FG; N e,t is the total number of EVs at timeslot t under the conditions that all EVs are assumed to be charged with rated active power. β t is the compensation factor. ω 1 and ω 2 are weighting coefficients. N T is the number of timeslots.
The operation should take into account a set of technical constraints, which are described in the following: 1.
ADN operation constraints: Most ADN with DGs are radial. Define the ADN with n nodes as a directed graph: ∀i, j ∈ V n ; (i, j) ∈ E n : Sustainability 2021, 13, 6489 9 of 26 P j = P L,j − P G,j − P R,j + P A,j , ∀j ∈ A a P L,j − P G,j − P R,j , ∀j ∈ V n − A a (20) where P ij , Q ij and e ij are the active power, reactive power and current flowing through branch ij, respectively. v i and v j are the voltage amplitudes of node i and node j, respectively. P j and Q j are the active power and reactive power of node j. Q L,j , Q G,j , Q R,j , Q C,j and Q A,j are the reactive powers of conventional loads, FGs, RES, SVC and EVA. P L,j , P G,j , P R,j and P A,j are the active powers of base loads, FGs, RES and EVA.

FGs operation constraints:
To guarantee the safe operation of FGs, the following constraints are considered [21]: where P min

3.
RES operation constraints: The reactive power output of photovoltaic (PV) can be ignored, and the PV output adopts maximum power point tracking (MPPT) technology [22]: where P v j,t and P vm j,t are the actual active power and predicted active power of PV at timeslot t. Similarly, the output of wind turbines (WT) can be expressed as follows: where P w j,t and P wm j,t are the actual active power and predicted active power of WT at timeslot t.
Taking the doubly fed induction generator (DFIG) [23], for example, the active and reactive output characteristics are shown in Figure 6. The area surrounded by the solid line in Figure 6 is the active and reactive power operation domain of DFIG. It can be seen that it presents a strong, nonlinear coupling relationship. In order to simplify processing, the active and reactive operation domain is substituted by the shaded area, which can be expressed as the following: where Q w j,t is reactive power of WT at timeslot t. (Q w j,min , 0), (Q w1 j,t , P w1 j,t ), (Q w2 j,t , P w2 j,t ), (Q w3 j,t , P w3 j,t ), (Q w4 j,t , P w4 j,t ) and (Q w j,max , 0) are characteristic points of the DFIG.

SVC operation constraints:
This paper assumes that reactive power compensation of the SVC is continuous: where Q max C,j and Q min C,j are the maximum and minimum values of the reactive power of the SVC.

The Solution Strategy of the Proposed Model
The second term in Equation (17) consists of the absolute value of the squared variables. The linearization of the objective function can reduce the computational complexity and improve the efficiency of the solution. Equation (17) can be rewritten as Equation (27) at the first step, and then v 2 b,t and v b,t are transformed into Equation (28) by using the conventional piecewise linearization method [24]: where v + b,t,k and ε k are defined as auxiliary variables. ∆B is the width of each segment in piecewise representation of the square function. n s is the number of segments. Figure 7 is the graphical illustration of the linearization method. According to Equation (19), due to the existence of complex quadratic terms in the power flow constraint, the optimization problem is non-convex, and the distributed optimization algorithm is difficult to ensure convergence when applied to general, non-convex problems. To realize the convexity of Equation (19), v 2 i = P 2 ij + Q 2 ij I ij can be expressed as Equations (29) and (30).
Based on the piecewise linearization of the objective function and the convexity of the power flow constraint, a decentralized algorithm, which is developed based on the ADMM, is introduced to the model solution in order to reduce the computational burden.
In order to realize decentralized optimization, the power of EVA is used as the coupling variable to establish the coupling constraint. Therefore, at the node i, which contains EVA, a set of virtual variables representing the injected power of EVA nodes are added, relative to the injected power of the ADN side. The scheme in Figure 8 is adopted to realize the decoupling between the ADN side and EVA side. Based on the decomposition scheme of ADN and EVA sides, the optimization task can be expressed as follows: The objective function (31) is quadratic constrained linear programming and the ADMM is used for the implementation of the distributed solution in this paper. The coordination of the injected power (P i ,Q i ) on the ADN side and the injected power (P i ,Q i ) on the EVA side are achieved by adding the Lagrange penalty function into the objective function. The suboptimal models of EVA and ADN can be expressed as Equations (33) and (34).
are the optimal solution; κ i,t and δ i,t are Lagrangian multipliers; and ρ is the penalty factor of ADMM. Suboptimal models of ADN and EVA are solved independently, and the boundary variable information is exchanged until the convergence conditions are met: where r and ε are the number of iterations and threshold of convergence criteria, separately. The specific steps of the distributed optimization solution based on ADMM are as follows: Step 1: Initialize variablesP i,t (0),P i,t (0),Q i,t (0),Q i,t (0), κ i,t (0) and δ i,t (0), and set r = 1.
Step 2: EVA optimizes Equation (34) and transmits the optimization resultsP i,t and Q i,t to ADN.
Step 3: ADN optimizes the model of Equation (33) after it receives the optimization results of EVA and passes the optimization resultsP i,t andQ i,t to EVA.
Step 4: The convergence test is carried out according to Equation (35). If satisfied, the optimization is completed, and the optimized results are outputted. Otherwise, κ b,t and δ b,t are updated according to Equation (36). Set the number of iterations r = r + 1. Return to Step 2.

System Data
In this paper, the numerical test is based on the radial cable 33-bus ADN with a normal voltage of 12.66 kV. The line parameters of the ADN can be found in reference [25]. The optimization horizon is 24 h of a day, which is divided into N T = 96 timeslots. The value of ε is set as 10 −4 . Figure 9 shows the topology of the ADN. Without loss of generality, it is assumed that all conventional loads are residential loads with known load curves as illustrated in Figure 10a. The green curve in Figure 10a shows the spot price of active power, which is associated with a typical day [26]. The value of β t is equal to half of the spot price in this paper. In the test case, there are four EVAs that are connected to buses 13, 21, 25, 29, and the maximum acceptable numbers of EVs are 500, 500, 750, 900, respectively. To simplify the calculation, it is assumed that all EVs are of the same type, and the ratio of the rated power charging EVs, non-discharging EVs and flexible charging-discharging EVs is 2:3:5. Tables 1 and 2 show the FG parameters and EV parameters, respectively. Meanwhile, the probability distribution of the EVs' arrival-departure time and daily mileage were obtained from the literature [27]. Then, the Monte Carlo method is used to calculate the number of EVs plugged into the ADN in a day as shown in Figure 10b. The output of WT and PV are shown in Figure 10c,d.    The test case is solved by programming in the environment of MATLAB 2016a. Meanwhile, the CPLEX 12.6 solver is employed.
As summarized in Table 3, four different cases are analyzed to investigate employing EVs' conjugate active and reactive power control capabilities in an ADN operation. The preference of EV users are ignored in case I and case II. In case I, all EVs plugged-in are charged with rated active power, and the reactive power is zero. In case II, all EVs are supposed to be flexible charging-discharging for which both the active power and reactive power are adjustable. The preference of EV users is considered in case III and case IV. In case III, the reactive power of EVs is zero and the active power is adjustable, except for EVs charged with rated power. In case IV, for an adjustable charging EV, the operation domain of the active and reactive power is shown in Figure 5a; for a flexible charging-discharging EV, the operation domain of the active and reactive power is illustrated in Figure 5b. Table 3. The characteristics of multiple cases.

Reactive Power
Case I supposed to be flexible charging-discharging for which both the active power and reactive power are adjustable. The preference of EV users is considered in case III and case IV. In case III, the reactive power of EVs is zero and the active power is adjustable, except for EVs charged with rated power. In case IV, for an adjustable charging EV, the operation domain of the active and reactive power is shown in Figure 5a; for a flexible chargingdischarging EV, the operation domain of the active and reactive power is illustrated in Figure 5b.

Coefficients Selection of Objective Function
It can be seen, according to Equation (31), that the two terms of the objective function are the total cost of energy and the summation of the voltage deviation, respectively. Table  4 shows the results of the two terms of the objective function under different values for 1 ω and 2 ω . According to Table 4, in the case of 1 ω = 1, 2 ω = 0 ( 1 ω = 0, 2 ω = 1), the energy cost (voltage profile) is considered the only objective function. In the third case, the value of 1 ω is set to be equal to 5.58 × 10 −5 2 ω . The meaning of this value selection is expressed as follows:

Items
Values of 1 ω and 2 ω Weight coefficient According to the results given in Table 4, it can be found that the maximum and minimum values of the total energy cost are CNY 6.08 × 10 4 and CNY 5.75 × 10 4 , respectively. Meanwhile, the maximum and minimum values of voltage are 1.107 p.u. and 0.923 p.u., respectively. Therefore, it can be concluded that the value spans of two parts of the objective function (i.e., total energy cost and voltage profile) are 3300 and 0.184, respectively. To make the two parts of the objective function have the same standard, the first part of the objective function (total energy cost) is divided by CNY 3300 and the second part of the objective function (voltage profile) is divided by 0.184 p.u. In other words, this means that we evaluate 1 ω and 2 ω in such a way that 1 ω = 5.58 × 10 −5 2 ω ( 1 ω = supposed to be flexible charging-discharging for which both the active power and rea tive power are adjustable. The preference of EV users is considered in case III and case IV In case III, the reactive power of EVs is zero and the active power is adjustable, except fo EVs charged with rated power. In case IV, for an adjustable charging EV, the operatio domain of the active and reactive power is shown in Figure 5a; for a flexible charging discharging EV, the operation domain of the active and reactive power is illustrated Figure 5b.

Coefficients Selection of Objective Function
It can be seen, according to Equation (31), that the two terms of the objective functio are the total cost of energy and the summation of the voltage deviation, respectively. Tab

Items
Values of 1 ω and 2 ω Weight coefficient According to the results given in Table 4, it can be found that the maximum an minimum values of the total energy cost are CNY 6.08 × 10 4 and CNY 5.75 × 10 4 , respe tively. Meanwhile, the maximum and minimum values of voltage are 1.107 p.u. and 0.92 p.u., respectively. Therefore, it can be concluded that the value spans of two parts of th objective function (i.e., total energy cost and voltage profile) are 3300 and 0.184, respe tively. To make the two parts of the objective function have the same standard, the fir part of the objective function (total energy cost) is divided by CNY 3300 and the secon part of the objective function (voltage profile) is divided by 0.184 p.u. In other words, th means that we evaluate 1 ω and 2 ω in such a way that 1 ω = 5.58 × 10 −5 2 ω ( 1 ω supposed to be flexible charging-discharging for which both the active pow tive power are adjustable. The preference of EV users is considered in case III In case III, the reactive power of EVs is zero and the active power is adjustab EVs charged with rated power. In case IV, for an adjustable charging EV, t domain of the active and reactive power is shown in Figure 5a; for a flexib discharging EV, the operation domain of the active and reactive power is Figure 5b.

Coefficients Selection of Objective Function
It can be seen, according to Equation (31), that the two terms of the objec are the total cost of energy and the summation of the voltage deviation, respe 4 shows the results of the two terms of the objective function under differe

Items
Values of 1 ω and 2 ω Weight coefficient Energy cost (CNY) 5.75 × 10 4 6.08 × 10 4 5.93 Voltage profile (p.u.) 1.107 0.923 0.9 According to the results given in Table 4, it can be found that the m minimum values of the total energy cost are CNY 6.08 × 10 4 and CNY 5.75 tively. Meanwhile, the maximum and minimum values of voltage are 1.107 p p.u., respectively. Therefore, it can be concluded that the value spans of two objective function (i.e., total energy cost and voltage profile) are 3300 and 0 tively. To make the two parts of the objective function have the same stand part of the objective function (total energy cost) is divided by CNY 3300 an part of the objective function (voltage profile) is divided by 0.184 p.u. In othe means that we evaluate 1 ω and 2 ω in such a way that 1 ω = 5.58 × 1 No reactive power. All EVs are charged with rated active power.
Case II charged with rated active power, and the reactive power is zero. In case II, all EVs are supposed to be flexible charging-discharging for which both the active power and reactive power are adjustable. The preference of EV users is considered in case III and case IV. In case III, the reactive power of EVs is zero and the active power is adjustable, except for EVs charged with rated power. In case IV, for an adjustable charging EV, the operation domain of the active and reactive power is shown in Figure 5a; for a flexible chargingdischarging EV, the operation domain of the active and reactive power is illustrated in Figure 5b.

Coefficients Selection of Objective Function
It can be seen, according to Equation (31), that the two terms of the objective function are the total cost of energy and the summation of the voltage deviation, respectively. Table  4 shows the results of the two terms of the objective function under different values for According to the results given in Table 4, it can be found that the maximum and minimum values of the total energy cost are CNY 6.08 × 10 4 and CNY 5.75 × 10 4 , respectively. Meanwhile, the maximum and minimum values of voltage are 1.107 p.u. and 0.923 p.u., respectively. Therefore, it can be concluded that the value spans of two parts of the objective function (i.e., total energy cost and voltage profile) are 3300 and 0.184, respectively. To make the two parts of the objective function have the same standard, the first part of the objective function (total energy cost) is divided by CNY 3300 and the second part of the objective function (voltage profile) is divided by 0.184 p.u. In other words, this means that we evaluate 1 ω and 2 ω in such a way that 1 ω = 5.58 × 10 −5 2 ω ( 1 ω = EVs charged with rated power. In case IV, for an adjustable charging EV, the operatio domain of the active and reactive power is shown in Figure 5a; for a flexible charging discharging EV, the operation domain of the active and reactive power is illustrated i Figure 5b. ω is set to be equal to 5.58 × 10 −5 2 ω . The meaning of this value selection expressed as follows: According to the results given in Table 4, it can be found that the maximum an minimum values of the total energy cost are CNY 6.08 × 10 4 and CNY 5.75 × 10 4 , respe tively. Meanwhile, the maximum and minimum values of voltage are 1.107 p.u. and 0.92 p.u., respectively. Therefore, it can be concluded that the value spans of two parts of th objective function (i.e., total energy cost and voltage profile) are 3300 and 0.184, respe tively. To make the two parts of the objective function have the same standard, the fir part of the objective function (total energy cost) is divided by CNY 3300 and the secon part of the objective function (voltage profile) is divided by 0.184 p.u. In other words, th means that we evaluate 1 ω and 2 ω in such a way that 1 ω = 5.58 × 10 −5 2 ω ( 1 ω EVs charged with rated power. In case IV, for an adjustable charging EV, t domain of the active and reactive power is shown in Figure 5a; for a flexib discharging EV, the operation domain of the active and reactive power is Figure 5b.

Coefficients Selection of Objective Function
It can be seen, according to Equation (31), that the two terms of the objec are the total cost of energy and the summation of the voltage deviation, respec 4 shows the results of the two terms of the objective function under differe

Items
Values of According to the results given in Table 4, it can be found that the m minimum values of the total energy cost are CNY 6.08 × 10 4 and CNY 5.75 tively. Meanwhile, the maximum and minimum values of voltage are 1.107 p p.u., respectively. Therefore, it can be concluded that the value spans of two objective function (i.e., total energy cost and voltage profile) are 3300 and 0 tively. To make the two parts of the objective function have the same stand part of the objective function (total energy cost) is divided by CNY 3300 an part of the objective function (voltage profile) is divided by 0.184 p.u. In othe means that we evaluate Case III supposed to be flexible charging-discharging for which both the active power and reactive power are adjustable. The preference of EV users is considered in case III and case IV. In case III, the reactive power of EVs is zero and the active power is adjustable, except for EVs charged with rated power. In case IV, for an adjustable charging EV, the operation domain of the active and reactive power is shown in Figure 5a; for a flexible chargingdischarging EV, the operation domain of the active and reactive power is illustrated in Figure 5b.

Coefficients Selection of Objective Function
It can be seen, according to Equation (31), that the two terms of the objective function are the total cost of energy and the summation of the voltage deviation, respectively. Table  4 shows the results of the two terms of the objective function under different values for  According to the results given in Table 4, it can be found that the maximum and minimum values of the total energy cost are CNY 6.08 × 10 4 and CNY 5.75 × 10 4 , respectively. Meanwhile, the maximum and minimum values of voltage are 1.107 p.u. and 0.923 p.u., respectively. Therefore, it can be concluded that the value spans of two parts of the objective function (i.e., total energy cost and voltage profile) are 3300 and 0.184, respectively. To make the two parts of the objective function have the same standard, the first part of the objective function (total energy cost) is divided by CNY 3300 and the second part of the objective function (voltage profile) is divided by 0.184 p.u. In other words, this means that we evaluate ( 1 ω = supposed to be flexible charging-discharging for which both the active power and rea tive power are adjustable. The preference of EV users is considered in case III and case IV In case III, the reactive power of EVs is zero and the active power is adjustable, except fo EVs charged with rated power. In case IV, for an adjustable charging EV, the operatio domain of the active and reactive power is shown in Figure 5a; for a flexible charging discharging EV, the operation domain of the active and reactive power is illustrated i Figure 5b.

Coefficients Selection of Objective Function
It can be seen, according to Equation (31), that the two terms of the objective functio are the total cost of energy and the summation of the voltage deviation, respectively. Tab  According to the results given in Table 4, it can be found that the maximum an minimum values of the total energy cost are CNY 6.08 × 10 4 and CNY 5.75 × 10 4 , respe tively. Meanwhile, the maximum and minimum values of voltage are 1.107 p.u. and 0.92 p.u., respectively. Therefore, it can be concluded that the value spans of two parts of th objective function (i.e., total energy cost and voltage profile) are 3300 and 0.184, respe tively. To make the two parts of the objective function have the same standard, the fir part of the objective function (total energy cost) is divided by CNY 3300 and the secon part of the objective function (voltage profile) is divided by 0.184 p.u. In other words, th means that we evaluate charged with rated active power, and the reactive power is zero. In case I supposed to be flexible charging-discharging for which both the active pow tive power are adjustable. The preference of EV users is considered in case III In case III, the reactive power of EVs is zero and the active power is adjustab EVs charged with rated power. In case IV, for an adjustable charging EV, t domain of the active and reactive power is shown in Figure 5a; for a flexib discharging EV, the operation domain of the active and reactive power is Figure 5b.

Coefficients Selection of Objective Function
It can be seen, according to Equation (31), that the two terms of the objec are the total cost of energy and the summation of the voltage deviation, respe 4 shows the results of the two terms of the objective function under differe

Items
Values of 1 ω and 2 ω Weight coefficient Energy cost (CNY) 5.75 × 10 4 6.08 × 10 4 5.93 Voltage profile (p.u.) 1.107 0.923 0.9 According to the results given in Table 4, it can be found that the m minimum values of the total energy cost are CNY 6.08 × 10 4 and CNY 5.75 tively. Meanwhile, the maximum and minimum values of voltage are 1.107 p p.u., respectively. Therefore, it can be concluded that the value spans of two objective function (i.e., total energy cost and voltage profile) are 3300 and 0 tively. To make the two parts of the objective function have the same stand part of the objective function (total energy cost) is divided by CNY 3300 an part of the objective function (voltage profile) is divided by 0.184 p.u. In othe means that we evaluate 1 ω and 2 ω in such a way that 1 ω = 5.58 × 1 Considering preference of EV users. Operation with only active power.
Case IV charged with rated active power, and the reactive power is zero. In case II, all EVs are supposed to be flexible charging-discharging for which both the active power and reactive power are adjustable. The preference of EV users is considered in case III and case IV. In case III, the reactive power of EVs is zero and the active power is adjustable, except for EVs charged with rated power. In case IV, for an adjustable charging EV, the operation domain of the active and reactive power is shown in Figure 5a; for a flexible chargingdischarging EV, the operation domain of the active and reactive power is illustrated in Figure 5b.

Coefficients Selection of Objective Function
It can be seen, according to Equation (31), that the two terms of the objective function are the total cost of energy and the summation of the voltage deviation, respectively. Table  4 shows the results of the two terms of the objective function under different values for According to the results given in Table 4, it can be found that the maximum and minimum values of the total energy cost are CNY 6.08 × 10 4 and CNY 5.75 × 10 4 , respectively. Meanwhile, the maximum and minimum values of voltage are 1.107 p.u. and 0.923 p.u., respectively. Therefore, it can be concluded that the value spans of two parts of the objective function (i.e., total energy cost and voltage profile) are 3300 and 0.184, respectively. To make the two parts of the objective function have the same standard, the first part of the objective function (total energy cost) is divided by CNY 3300 and the second part of the objective function (voltage profile) is divided by 0.184 p.u. In other words, this means that we evaluate ( 1 ω = charged with rated active power, and the reactive power is zero. In case II, all EVs ar supposed to be flexible charging-discharging for which both the active power and rea tive power are adjustable. The preference of EV users is considered in case III and case IV In case III, the reactive power of EVs is zero and the active power is adjustable, except fo EVs charged with rated power. In case IV, for an adjustable charging EV, the operatio domain of the active and reactive power is shown in Figure 5a; for a flexible charging discharging EV, the operation domain of the active and reactive power is illustrated i Figure 5b.

Coefficients Selection of Objective Function
It can be seen, according to Equation (31), that the two terms of the objective functio are the total cost of energy and the summation of the voltage deviation, respectively. Tab  According to the results given in Table 4, it can be found that the maximum an minimum values of the total energy cost are CNY 6.08 × 10 4 and CNY 5.75 × 10 4 , respe tively. Meanwhile, the maximum and minimum values of voltage are 1.107 p.u. and 0.92 p.u., respectively. Therefore, it can be concluded that the value spans of two parts of th objective function (i.e., total energy cost and voltage profile) are 3300 and 0.184, respe tively. To make the two parts of the objective function have the same standard, the fir part of the objective function (total energy cost) is divided by CNY 3300 and the secon part of the objective function (voltage profile) is divided by 0.184 p.u. In other words, th means that we evaluate ( 1 ω charged with rated active power, and the reactive power is zero. In case II supposed to be flexible charging-discharging for which both the active pow tive power are adjustable. The preference of EV users is considered in case III In case III, the reactive power of EVs is zero and the active power is adjustab EVs charged with rated power. In case IV, for an adjustable charging EV, t domain of the active and reactive power is shown in Figure 5a; for a flexib discharging EV, the operation domain of the active and reactive power is Figure 5b.

Coefficients Selection of Objective Function
It can be seen, according to Equation (31), that the two terms of the objec are the total cost of energy and the summation of the voltage deviation, respec 4 shows the results of the two terms of the objective function under differe According to the results given in Table 4, it can be found that the m minimum values of the total energy cost are CNY 6.08 × 10 4 and CNY 5.75 tively. Meanwhile, the maximum and minimum values of voltage are 1.107 p p.u., respectively. Therefore, it can be concluded that the value spans of two objective function (i.e., total energy cost and voltage profile) are 3300 and 0 tively. To make the two parts of the objective function have the same stand part of the objective function (total energy cost) is divided by CNY 3300 an part of the objective function (voltage profile) is divided by 0.184 p.u. In othe means that we evaluate Considering preference of EV users. Operation with both active power and reactive power.

Coefficients Selection of Objective Function
It can be seen, according to Equation (31), that the two terms of the objective function are the total cost of energy and the summation of the voltage deviation, respectively. Table 4 shows the results of the two terms of the objective function under different values for ω 1 and ω 2 . According to Table 4, in the case of ω 1 = 1, ω 2 = 0 (ω 1 = 0, ω 2 = 1), the energy cost (voltage profile) is considered the only objective function. In the third case, the value of ω 1 is set to be equal to 5.58 × 10 −5 ω 2 . The meaning of this value selection is expressed as follows: According to the results given in Table 4, it can be found that the maximum and minimum values of the total energy cost are CNY 6.08 × 10 4 and CNY 5.75 × 10 4 , respectively. Meanwhile, the maximum and minimum values of voltage are 1.107 p.u. and 0.923 p.u., respectively. Therefore, it can be concluded that the value spans of two parts of the objective function (i.e., total energy cost and voltage profile) are 3300 and 0.184, respectively. To make the two parts of the objective function have the same standard, the first part of the objective function (total energy cost) is divided by CNY 3300 and the second part of the objective function (voltage profile) is divided by 0.184 p.u. In other words, this means that we evaluate ω 1 and ω 2 in such a way that ω 1 = 5.58 × 10 −5 ω 2 (ω 1 = 0.184/3300 ω 2 ). By setting two weight coefficients in this way, the two parts of the objective function are equally weighted.

Analysis of Different Test Cases
In case I, the reactive power capacity of EVs is not considered, and the active power of EV is not controllable. Once the EV is plugged into ADN, it is charged to the desired SOC with the rated active power. By comparing the active conventional load curves in Figures 10a and 11a, it can be seen that the load peak and valley of EVs and the conventional load are simultaneous, and the superposition of the two loads obviously further aggravates the load fluctuation of the ADN. Figure 11b shows that, due to the superposition of the peak load, the voltage at nodes 13 and 29 in peak intervals deviates from the normal voltage range, which poses a great threat to the operation safety of the ADN. Therefore, it is very important to achieve optimal control of plugged-in EVs. Comparing case II and case IV, as shown in Figure 12a, the adjustable range of active power in case II is larger than case IV (yellow area for case IV, green and yellow area for case II). The reason for this is that the preference of EV users is taken into account in case IV; some EVs need to be charged at the rated power, so the active power is not adjustable. In contrast, the active power of all EVs in case II is adjustable. The solid red and green lines in Figure 12a represent the EVs' load curves of case II and case IV, respectively. A larger number of EVs are plugged into ADN from 18:00 to 24:00; some EVs must be charged at the rated power in case IV, which results in a greater active power than that for case II. In case IV, once the EV charging at the rated power reaches the desired SOC, the charging process stops and the total active power may be lower than that in case II (e.g., 00:00-06:00). Figure 12b shows the summation of EVs' reactive power in cases II and IV. It can be seen that the reactive power curves of the two schemes are similar. Since both cases have enough EVs whose charging power can be regulated, the required reactive power of the ADN can be provided sufficiently in both cases based on the intelligent bi-directional charger. As the number of EVs charged by rated power is small, it can be seen from Figure 12c that whether the preference of users is considered has little influence on the voltage amplitude of the ADN. However, by comparing Figures 11b and 12c, it can be concluded that the voltage profiles of buses 13 and 29 during peak load hours could be improved noticeably, while both the active and reactive power of EVs are controllable. Figure 12d illustrates the voltage magnitude of the ADN in case II; the bus voltage magnitude is in the acceptable range (0.95-1.05 p.u.).
In order to evaluate the advantages of employing the EVs' reactive power control capability more directly and accurately, case III and case IV are compared, and the results are shown in Figures 13-15. The only difference between the two cases is that the reactive power control capability of EVs is considered in case IV. It can be seen from Figure 13a that, over 96 intervals throughout the day, the difference of active power between the two cases is not significant. Figure 13b shows the significant improvement of the voltage profile (especially in the 14:00-22:00 intervals) of bus 13 by employing the reactive power control (case IV) of EVs. Meanwhile, combined with the analysis in Figure 14, the reactive power compensation provided by SVG in case IV is significantly less than that in case III. The reactive power obtained from the upstream DN is also less than that of case III. It can be concluded that employing the reactive power control capability of EVs can not only improve the voltage profile, but also reduce the capacity of SVG configuration in the ADN and save the investment cost of the ADN. Figure 15 shows that the implementation of the chargers' reactive power capability prevents severe voltage drop of the buses.  Combined with the analysis in Figures 10a, 13a and 16, on the one hand, the active power control ability of EVs plays the effect of peak clipping and valley filling to some extent. On the other hand, due to the energy cost term in the objective function, EVs absorb more power during the off-peak intervals, which is beneficial for EV users to reduce the charging cost. From the perspective of ADN operation, the ADN absorbs more power from the upstream DN in the period when the market electricity price is lower; in the period when the electricity price is higher, the ADN fully mobilizes the active power of FGs, thus reducing the operation cost of the ADN.

Charging Process under Demand Preference
To be more intuitive, l m,k (m = 1, 2, 3,k = 1, 2, 3, 4) is an EV randomly selected from D m,k in case IV. The SOC curve of l m,k is shown in Figure 17. The SOC of the three categories of EVs in EVA 1 -EVA 4 is above 0.95 when they are off-grid, which fully guarantees the users' energy demands. Whenever the l 1,k is plugged into the ADN, it is charged at P r c,l until the desired SOC is reached; during the plugged-in intervals, the actual charging power of l 2,k is less than P r c,l , and the value of SOC increases slowly but there is no discharge process over its plugged-in intervals. Part of the plugged-in intervals of l 3,k are in the discharging state, but the SOC during the discharging process is higher than the threshold value of 0.5. It should be noted that the EV in D 2,k or D 3,k set is also constrained by both the plugged-in duration and energy demand. When the plugged-in duration is shorter and the energy demand is larger, although its active power is adjustable, it still maintains a high charging power in order to meet the energy demand, such as l 2,2 , l 3,4 . Obviously, the optimization model established in this paper sufficiently guarantees the demand preference of EV users.

Superiority Analysis of the Proposed Model
The model established in this paper quantifies the reactive power capacity of EV more accurately, and fully guarantees the preference of the users. However, none of the existing literature takes both of these two points into account. Therefore, to further verify the advantages of the distributed optimization strategy proposed in this paper, the proposed strategy is compared with the strategy of the stochastic, multi-objective based on normalized normal constraint (NNC) method proposed in reference [14]. It should be pointed out that the framework of the test case in reference [14] is significantly different from that in this paper. Therefore, we only use the theoretical method in reference [14] (defined as case V) to solve the test case (case IV) in this paper, and the objective function is also transformed into the multi-objective form. The comparison of the optimization results between case IV and case V are shown in Figure 18 and Table 5.  According to Figure 18 and Table 5, it can be seen that the difference in the active power optimization results is great, and the difference in the reactive power optimization is small. The total energy costs of case IV and case V are significantly different, but the difference in the average voltage deviation is not significant. Compared with that of case V, the total energy cost of case IV is reduced by CNY 1.09 × 10 4 . More significantly, the solving time of the two cases is greatly different. The distributed optimization strategy proposed in this paper only needs one third of the time compared to reference [14]. The solution efficiency is greatly improved in this paper.

Analysis in EVs Penetration Point of View
In this section, case IV is also taken as an example to analyze the influence of different EVs' distributions on the proposed model. The different distributions of EVs include two aspects: For the first one, the proportions of the three categories of EVs are different. The other is that the number of EVs is different. The second aspect is reflected in the test cases (the number of EVs managed by four EVAs is different). For the first aspect, the test scheme is designed as follows: Firstly, define α to be the ratio of EVs in the D 2,k and D 3,k set. So, α = 0.8 in Section 3.2. Figure 19 shows the reactive power of the EVs, SVG, upstream grid and active power of EVs when the value of α is equal to 0.2, 0.6 and 1, respectively. According to Figure 19a-c, the reactive power from the upstream grid and the reactive power of EVs and SVG have a slight difference in the case of α = 1 or α = 0.6. This is because the ratio of EVs in D 2,k and D 3,k sets is relatively high, which provide sufficient reactive power. When α = 0.2, most EVs are charged at the rated power, and the reactive power provided by the smart charger is significantly reduced, while the reactive power compensated by SVC is significantly increased, which means that the SVG configuration capacity in the ADN needs to be increased. Comparing Figure 10d with Figure 19d, it can be seen that the smaller the value of α, the more obvious the peak-valley simultaneity for both the EVs' load and conventional load and the worse the effect of load shifting. Figure 20 shows the voltage profile of bus 13 under different values of α. As the value of α increases, the voltage profile of bus is improved. However, when α = 0.8 or α = 1, there is little improvement in the voltage profile of bus 13, compared to α = 0.6. In order to intuitively show the relationship between the α value and the adjustable active power and voltage deviation, γ p is defined as the mean value of the adjustable active power of EVs over 96 intervals, and γ v is defined as the mean value of the voltage deviation from the standard voltage over 96 intervals.
Then, 100 different values of α between 0 and 1 are taken (the equal step size of 0.01 is taken in this paper). Other conditions are consistent with case IV. The relationship between γ p , γ v and α is shown in Figure 21. The scattered points in Figure 21 can be approximated as a spatial curve. When α = 0, every EV is charged at the rated charging power, so the active power is not adjustable; then, γ p = 0 and γ v reaches the maximum value of 5.5%. With the increase in α, γ p increases and γ v decreases but the decreasing trend of γ v tends to be moderate. This is because when the EVs in D 2,k and D 3,k sets reach a certain amount, sufficient reactive power capacity is provided, so the value of γ v does not have a significant change.
The results in Figure 22 are the sum of the reactive powers over 96 intervals. With the increasing proportion of EVs in D 2,k and D 3,k sets, the reactive power injected by SVG gradually decreases, while the reactive power injected by EVs gradually increases, and both trends tend to smoothen gradually. The reactive power absorbed by the ADN from the upstream DN is relatively stable.

Conclusions and Future Works
Based on an intelligent charger capable of four-quadrant operation and the scenario of large-scale, EV penetration, in this paper, an energy optimization model for the ADN, by employing cooperative active and reactive power management of EVs, is proposed. Some numerical case studies were investigated to demonstrate the superiority of the proposed model. The main achievements are summarized as follows: (1) The total energy cost is reduced and the voltage profile of ADN is improved by employing the model proposed in this paper. The energy cost is reduced by about 15.5 percent, and the voltage profile is improved by about 0.09 percent. (2) The demand preference of EV users is taken into account in the proposed model. The plugged-in EVs are divided into three categories: rated power charging EVs, nondischarging EVs, and flexible charging-discharging EVs. Users can choose different charging modes, which are more in line with the actual scenario. The study results show that the proposed model can fully guarantee the users' demand preference. (3) A distributed solution strategy based on ADMM is designed for the proposed model.
The experimental results show that the distributed strategy proposed in this paper can save about two thirds of the solving time.
The following aspects are worthy of further study in the future.
(1) Uncertainty: There are multiple uncertainties, such as the arrival and departure time of EVs, active power of WT and PV, system failure, and so on. The question of how to consider these uncertainties in the established model needs further research. (2) Interactivity: In practice, users can interact with the EVA by making an appointment to charge, switching charging modes, adjusting the charging power or discharging power, and so on. The question of how to improve the model to meet the timely interaction between the EV and EVA is another direction to be further studied. The average power is an active power used to charge or discharge the battery of EV, and the ripple power is a kind of ripple power flowing between the charger and the ADN, which is only temporarily stored in the capacitance of the charger. This part of the energy stored in the capacitance can be obtained by the integration of the ripple power: Equation (A6) shows that if the charger is used for reactive compensation, except for active power charging, the value of E ripple increases. This, in turn, causes more ripple currents of second harmonics and higher values of the ripple voltage (∆V dc ). Once V dc and ∆V dc are determined, the two parameters of the capacitor are determined: the capacitance value and the effective value of the ripple current of second harmonics. In practice, the maximum stored energy of the capacitor should be equal to the ripple energy of the charger: where V cmin and V cmax are the minimum and maximum voltages of the capacitor, respectively. Substituting Equation (A7) into Equation (A6), the following equation can be deduced.
After ignoring the PWM component, the effective value of the ripple voltage of second harmonics can be formulated by Equation (A9). v c = V cr + 1 2 × ∆V cr sin(2ωt) (A9) Therefore, the ripple current of second harmonics is as follows: i c = C·dv c /dt = ωC·∆V cr cos(2ωt). The effective value of i c : I c = ωC·∆V cr / √ 2. Substituting I c into Equation (A8), the following equation can be deduced.
Considering that the value of V c is fixed and the value of L c is small, this means that the value of ωL c /V 2 c is small. Equation (A10) can be further simplified by using the Taylor series: The maximum value of I c is S max / √ 2V cr . According to Equation (A11), the following equation can be deduced.
S max ≥ S·(1 − ωL c Q s /V 2 c ) (A12) Replace S with P 2 s + Q 2 s , the following equation can be obtained.