# Minimizing Harmonic Distortion Impact at Distribution System with Considering Large-Scale EV Load Behaviour Using Modified Lightning Search Algorithm and Pareto-Fuzzy Approach

Authors: S. N. Syed Nasir; J. J. Jamian; M. W. Mustafa
Journal: Complexity
 (2018)

---

## Abstract

This research is focusing on optimal placement and sizing of multiple variable passive filter (VPF) to mitigate harmonic distortion due to charging station (CS) at 449 bus distribution network. There are 132 units of CS which are scheduled based on user behaviour within 24 hours, with the interval of 15 minutes. By considering the varying of CS patterns and harmonic impact, Modified Lightning Search Algorithm (MLSA) is used to find 22 units of VPF coordination, so that less harmonics will be injected from 415 V bus to the medium voltage network and power loss is also reduced. Power system harmonic flow, VPF, CS, battery, and the analysis will be modelled in MATLAB/m-file platform. High Performance Computing (HPC) is used to make simulation faster. Pareto-Fuzzy technique is used to obtain sizing of VPF from all nondominated solutions. From the result, the optimal placements and sizes of VPF are able to reduce the maximum THD for voltage and current and also the total apparent losses up to 39.14%, 52.5%, and 2.96%, respectively. Therefore, it can be concluded that the MLSA is suitable method to mitigate harmonic and it is beneficial in minimizing the impact of aggressive CS installation at distribution network.


---

## Body

## 1. Introduction

The vision to have less carbon dioxide (CO2) emissions and less dependency on natural resources has encouraged Electric Vehicle (EV) to become an important option compared to conventional vehicle. Based on report in 2015, around 56.4% global crude oil is used for transportation sector [1]. Furthermore, the unstable price for crude oil has also influenced EV to become suitable alternative. The demand on EV, indirectly, has increased the number of CS installation in the distribution system [2]. However, power losses will increase when a large number of unplanned CS are installed. Furthermore, it will also introduce harmonic distortion due to the power electronic devices that convert alternating current (AC) to direct current (DC) at CS [3, 4]. This harmonic distortion will cause negative impact such as increment heating losses, shorter insulation lifespan, increased temperature and insulation stress, decreased power factor, and lower efficiency [5, 6].Currently, most of the researchers focus on strategies to reduce loss by coordinating the charging of EVs. Alonso et al. [7] have developed an optimization algorithm to coordinate the charging of EVs using Genetic Algorithm (GA). The charging schedule is based on optimal load pattern for EV charging-based reliability which takes into consideration the thermal line limits, the load on transformers, voltage limits, and parking availability patterns. The proposed method is able to reduce the cost for power system new investment. Furthermore, Bharati and Paudyal [8] have proposed framework to coordinate EVs in the distribution system using bilevel hierarchical vehicle-grid (VG) optimization. This framework requires information exchange among EV aggregators and grid controller and then the EV charging can be scheduled based on the minimum losses to the distribution network. Next, Masoum and Nabavi [9] have also proposed online-overnight PEV coordination using metaheuristic technique which is particle swarm optimization. This method introduces high priority and low priority customer that indicate the willingness of customer to be charged. High priority customer will get faster charging services with high tariff and the low priority customer will only charge their vehicle after power loss control in the distribution system is done. Even though many approaches have been introduced to minimize the EV impact, the focus is more towards minimizing the power loss. However, beside power loss, the harmonic injection to the distribution system when many CS are connected must also be taken into consideration in ensuring that the network can operate optimally.In general, there will be two types of harmonic that will be mitigated in distribution network, which are voltage and current harmonic. Voltage harmonic is measured at the bus, while current harmonic is measured at the lines and cables. Passive filter is the common component used to eliminate these harmonics. Passive filters can provide a low-impedance path to harmonic currents and hence prevent it from flowing into the systems. There are many papers that utilized passive filter to mitigate harmonics in power system using filter devices [10–14] and the most common ones are single tuned filters. Single tuned filters are designed to eliminate or reduce single frequency from the system depending on the design resistance, capacitance, and inductance value [15]. Sakar et al. [16] have suggested a hosting capacity determination for a distorted distribution system due to Photovoltaic connection. In their research, passive filter is used to increase the harmonic-constrained hosting capacity which then improves the voltage, power factor, and filtering the harmonics. Author in [17] proposed an Asymmetric Synchronous Reference Frame Control scheme and harmonic voltage compensator with a Static Var Compensator connected to the grid to mitigate the harmonic. The author used harmonic voltage compensator to reduce voltage THD. Even though many researchers have used passive filter to mitigate the harmonic, very few researchers considered both voltage and current harmonics in their analysis especially for CS.Therefore, this research proposes a coordination of 22 VPF when 132 unit of CS is adopted in a 449 bus radial distribution system. Section2 presents the modelling of load profile, CS, battery, passive filter, harmonic, CS scheduling, and radial load flow used in this study. The introduction of MLSA technique will also be covered in this section. The methodology of the proposed coordination will be discussed clearly in Section 3. Result and discussion for 22 units of filter coordination in 449 bus are presented in Section 4. Lastly, conclusion of this research is in Section 5.


## 2. Problem Modelling and Formulation

A typical 449 bus radial distribution system is used in this research to determine the optimal placement and sizing of VPF. Furthermore, forward/backward sweep method is chosen due to the accuracy of this load flow analysis in solving distribution load flow [18–21]. The harmonic pattern for individual EV charger is modelled based on actual single phase charger impact as measured in [22].


### 2.1. IEEE 31 Bus with 22 Low Voltage 415 V Bus System

Figures1 and 2 show a single-line diagram for 23 kV IEEE 31 bus radial system with low voltage 415 V buses residential feeder that consist of several CS, respectively. The total number of buses in the network (medium and low voltages) is 449 buses. There are 6 units of CS installed at one residential area which add up to 132 units for overall CS in this distribution system. Line data for low voltage 415 V buses is shown in Table 1.Table 1
Line data for 415 V buses.
Line
Line resistanceR (Ω)
Line reactanceX (Ω)
Bus A
Bus B
a
b
0.0415
0.0145
b
c
0.0424
0.0189
c
d
0.0444
0.0198
d
e
0.0369
0.0165
e
f
0.0520
0.0232
f
g
0.0524
0.0234
g
h
0.0005
0.0002
g
i
0.2002
0.0199
g
j
1.7340
0.1729
f
k
0.2607
0.0260
f
l
1.3605
0.1357
d
m
0.1400
0.0140
c
n
0.7763
0.0774
b
o
0.5977
0.0596
a
p
0.1423
0.0496
p
q
0.0837
0.0292
q
r
0.3123
0.0311
a
s
0.0163
0.0062
DT
0.0000
0.0654Figure 1
Single-line diagram for 449 bus radial distribution system.Figure 2
Single-line diagram for 415 V buses.


### 2.2. Harmonic Load Flow Analysis

Many techniques can be used to perform harmonic analysis such as in [23–25] and the most typical and simplest one is by using current injection analysis method [26]. The electrical parameters that affected the harmonic existence are the line impedances, load variation impedances, and filter impedances values. Equations (1) and (2) represent the line impedance formulation and the impedance for single-tuned filter at harmonic h, respectively, which indirectly caters for resonance impact. Equations (3), (4), and (5) show the impedance formulation for the load used in setting up the harmonic admittance matrix. The harmonic flow is calculated using (6), which consists of harmonic admittance matrix and harmonic current injection as per (7) and (8).(1)Zl,h=Rl+jhωLl,(2)ZFilter,h=RFilter+jhωLFilter-j1hωCFilter,(3)ZL,h=RL+jXL,h,(4)RL=VL2PL,(5)XL,h=VL2hQL,(6)I-h=Y-h·V-h,(7)Y-h=y-11,hy-12,h⋯y-1N,hy-21,hy-22,h⋯⋮⋮⋮⋯⋮y-N1,h⋯⋯y-NN,h,(8)I-h=I-1,hI-2,h⋮I-N,h.


### 2.3. Modelling of CS and Battery

CS is typically divided into three categories which are level 1, level 2, and level 3. The level of the CS is based on the power consumed by the charging station. Level 1 takes 6–8 hours, while level 3 requires 30 minutes for fully charging the EV’s battery. In this research, level 1 CS is selected since this research focuses on residential area which normally involves single-phase supply. The main different in this work is on CS modelling and its behaviour. The CS modelling in this research has considered the worst harmonic value achievable by individual CS based on the battery characteristic and State of Charge (SOC) state. The initial and final (plug-in and plug-out) battery SOC value depended on customer request, while the current produced by CS to the EV battery is calculated based on that SOC state. SOC state will be updated in every 15 minutes(Δt) and can be calculated using (9), while the current produced by CS can be calculated using (10).(9)SOCΔtk+1,i=SOCΔtk,i+ΔtQiIΔtk,i×100,(10)IΔtk,i=4·Δt·Ri·PCSΔtk,i·η·CRiΔtk,i+Voc,i2-Voc,i2Ri·CRiΔtk,i,where Voc,i is open circuit voltage for th node (V), Qi is rated battery ampere hour for the ith PEV (Ah), Ri is battery equivalent internal resistance for the ith node (ohm), CRiΔtk,i is maximum charging rate for the ith PEV (A), SOCΔtk,i is state of charge of the ith PEV at kth time slot (%), SOCΔtk+1,i is state of charge of the ith PEV at next kth time slot (%), IΔtk,i is charging current for the ith PEV at current time slot (A), and PCSΔtk,i is consumed power for the ith PEV (kW).Figure3 shows the open circuit voltage (OCV) for lithium iron phosphate (LFP) batteries based on the SOC of battery [27] that can also be presented using (11). Thus, the relation between SOC and VOCV will influence the total current injected during charging process.(11)VOCV=3.135-0.685-ln⁡SOC0.478-1.342·SOC+1.734·e0.4SOC-1.Figure 3
Relationship between open circuit voltage and SOC for LFP battery.Since there is conversion from AC to DC during charging process, the harmonic impact to the network is considered in this study. The harmonic injection is basically based on SOC state, which depended on battery SOC status and is being modelled based on practical data from [18]. Since there are 132 CS units installed within the 22 low voltage residential areas (6 units for every residential area) it might have significant impact on THDV and THDi value when all CS operated simultaneously. All CS and battery data in the residential area are presented as per Table 2.Table 2
CS and battery data.
Type
Battery capacity (kW)
Charging rate (A)
Power usage (kW)
CS efficiency
1
10
0.052
0.52
0.93
2
15
0.125
1.875
0.93
3
10
0.235
2.35
0.93
4
15
0.15
2.25
0.93
5
10
0.215
2.15
0.93
6
20
0.102
2.04
0.93


### 2.4. Modelling of Load Profile and Charging Station Operation

The load profile for 449 bus system is design based on typical load consumption in residential area during normal day. The peak demand in the system is occurring at the afternoon due to the several usages of electrical appliances, while the lowest demand is at the night when majority of people are asleep. However, the new electrical load CS is normally being charged at the night time. In this study, it is assumed that harmonic produced in the distribution system is coming from CS. Figure4 shows the loading produced by all CS in the distribution system based on design of CS operation as in Table 3. CS operation is designed based on random behaviour and will penetrate at the night hour. Figure 5 shows the summation of original load profile with CS load. Even though most of the CS is operating during low loaded condition at distribution system, it will cause higher impact on the total harmonic distortion as well as power losses.Table 3
Bus
Plug-in
Plug-out
SOCinitial
SOCReq
Type
33
18
64
28
58
1
37
21
64
11
80
2
39
23
65
10
82
3
46
21
64
1
78
4
48
19
63
1
69
5
49
22
65
4
85
6
52
15
61
28
58
1
56
18
61
11
80
2
58
20
62
10
82
3
65
18
61
1
78
4
67
16
60
1
69
5
68
19
62
4
85
6
71
23
69
28
58
1
75
26
69
11
80
2
77
28
70
10
82
3
84
26
69
1
78
4
86
24
68
1
69
5
87
27
70
4
85
6
90
13
59
28
58
1
94
16
59
11
80
2
96
18
60
10
82
3
103
16
59
1
78
4
105
14
58
1
69
5
106
17
60
4
85
6
109
17
60
4
85
6
113
14
58
1
69
5
115
16
59
1
78
4
122
18
60
10
82
3
124
16
59
11
80
2
125
13
59
28
58
1
128
26
69
4
85
6
132
23
67
1
69
5
134
25
78
1
78
4
141
27
82
10
82
3
143
25
80
11
80
2
144
22
58
28
58
1
147
33
76
4
85
6
151
30
74
1
69
5
153
32
75
1
78
4
160
34
76
10
82
3
162
32
75
11
80
2
163
29
75
28
58
1
166
18
64
28
58
1
170
21
64
11
80
2
172
23
65
10
82
3
179
21
64
1
78
4
181
19
63
1
69
5
182
22
65
4
85
6
185
15
61
28
58
1
189
18
61
11
80
2
191
20
62
10
82
3
198
18
61
1
78
4
200
16
60
1
69
5
201
19
62
4
85
6
204
23
69
28
58
1
208
26
69
11
80
2
210
28
70
10
82
3
217
26
69
1
78
4
219
24
68
1
69
5
220
27
70
4
85
6
223
13
59
28
58
1
227
16
59
11
80
2
229
18
60
10
82
3
236
16
59
1
78
4
238
14
58
1
69
5
239
17
60
4
85
6
242
17
60
4
85
6
246
14
58
1
69
5
248
16
59
1
78
4
255
18
60
10
82
3
257
16
59
11
80
2
258
13
59
28
58
1
261
26
69
4
85
6
265
23
67
1
69
5
267
25
78
1
78
4
274
27
82
10
82
3
276
25
80
11
80
2
277
22
58
28
58
1
280
33
76
4
85
6
284
30
74
1
69
5
286
32
75
1
78
4
293
34
76
10
82
3
295
32
75
11
80
2
296
29
75
28
58
1
299
18
64
28
58
1
303
21
64
11
80
2
305
23
65
10
82
3
312
21
64
1
78
4
314
19
63
1
69
5
315
22
65
4
85
6
318
15
61
28
58
1
322
18
61
11
80
2
324
20
62
10
82
3
331
18
61
1
78
4
333
16
60
1
69
5
334
19
62
4
85
6
337
23
69
28
58
1
341
26
69
11
80
2
343
28
70
10
82
3
350
26
69
1
78
4
352
24
68
1
69
5
353
27
70
4
85
6
356
13
59
28
58
1
360
16
59
11
80
2
362
18
60
10
82
3
369
16
59
1
78
4
371
14
58
1
69
5
372
17
60
4
85
6
375
17
60
4
85
6
379
14
58
1
69
5
381
16
59
1
78
4
388
18
60
10
82
3
390
16
59
11
80
2
391
13
59
28
58
1
394
26
69
4
85
6
398
23
67
1
69
5
400
25
78
1
78
4
407
27
82
10
82
3
409
25
80
11
80
2
410
22
58
28
58
1
413
33
76
4
85
6
417
30
74
1
69
5
419
32
75
1
78
4
426
34
76
10
82
3
428
32
75
11
80
2
429
29
75
28
58
1
432
18
64
28
58
1
436
21
64
11
80
2
438
23
65
10
82
3
445
21
64
1
78
4
447
19
63
1
69
5
448
22
65
4
85
6Figure 4
Load produced by CS over time.Figure 5
Load profile before and after CS operation.


### 2.5. Assumptions and Constraints

There are few assumptions and constraints are considered in this research as shown below:(i)
CS can be plugged in/plugged out at any time according to the customer’s request. Customers will input their requested plug-out time and requested final SOCs at the time of plug-in. Once SOC reaches requested SOC, CS will be switched to a standby mode.(ii)
The time slot is 15 min (Δt) or equal to 96 slots for one day.(iii)
The aggregator has access to CS information including their bus locations, harmonic distortion information, charger types, battery sizes, plug-in time, and plug-out time.(iv)
CS are controllable and have variable charging functions. During the charging process, each CS is assumed as a variable active load. The power that is used to charge the battery is based on calculation that is done by the aggregator.(v)
The requested time for each CS must be greater than the minimum charging time required to charge the battery.Since 22 variable filters will be placed in the network, the total unknown variables will become 44, which are 22 locations and 22 optimal sizes. All these parameters will have their own constraints that need to be considered. In general, the parameters can be divided into 2 categories which are as follows:(i)
Filter location: the filter will be placed in low voltage bus to avoid any harmonic injection to upper feeder. There are 22 locations in 449 bus radial distribution system. Furthermore, only one VPF will be placed at every low voltage 415 buses system. The constraints are as follows:(12)ai≤Filteri≤si,i=1…22.(ii)
Filter reactive value: each individual variable filter reactive value is limited to(13)0kVAr≤Qi≤40kVAr,i=1…22.


### 2.6. Modelling of Passive Filter

The main function of a passive filter is to sink the harmonic current that flows in the system based on a selected frequency. The filter impedance will become very low to allow the harmonic to sink. Single tuned filter is the most popular type of filter which is used widely in dealing with harmonic pollution especially in the industrial area [28, 29]. In this research, four units of single tuned filters are considered as one set of filter which can eliminate four frequencies as shown in Figure 6. Equation in [26] is used to calculate capacitor, inductance, and resistance components, respectively, as per (14). Capacitor is calculated based on injected reactive power (Q) and voltage (V) at that bus; meanwhile, inductance and resistance are based on the chosen harmonics (n) that need to be reduced. Four sets of filter will be used to eliminate 3rd, 5th, 7th, and 9th harmonic order in the network.(14)CFilter=Q2πfV2,LFilter=V22n2Q2πf,RFilter=V2nQ2.Figure 6
1 set of filter.


### 2.7. Multiobjective MLSA with Pareto Optimization

In this research, the optimization technique is used to identify optimized locations and sizes for 22 sets of filters in 449 bus radial distribution system. The Lightning Search Algorithm (LSA) is a metaheuristic optimization method inspired by lightning propagation from cloud to ground [30]. The LSA process basically starts with generation of random population before fitness calculation. The worst step leader will be eliminated every 10th iteration before the direction of the step leader is updated. Next, step leader movement is updated based on direction, shape parameter, and scale parameter. Lastly, forking phenomenon is taken place as per (7) for 1% of the better step leader. Since the problem faced in this research is too complicated with many local minima, modification on the existing LSA is needed. In this research, the MLSA is proposed with better convergence for this application compared to existing LSA. MLSA basically is the improvement of LSA method which is proven beneficial for this research.The placement and sizing of variable filter will be based on five parameters, which are maximumTHDV, maximum THDi, THDi summation, average THDi, and apparent power losses for overall system. The three parameters that involve THDi are to ensure that the THDi at all buses is within acceptable range. These parameters are then normalized to get the most accurate solution as defined by (15)–(19). The weight summation with normalized fitness function is shown in (20) after several manual trials are done in order to get the best optimal result. Last but not least, basic equations to calculate THDi and THDV are as shown by (21) and (22).(15)fit1=max⁡THDV,n=1…449,(16)fit2=max⁡THDi,n=1…30,(17)fit3=∑i=130THDi,(18)fit4=∑i=130THDi30,(19)fit5=∑i=1449si-loss,(20)Fit=0.1fit1+0.2fit2+0.2fit3+0.2fit4+0.3fit5,(21)THDV=∑h=29Vh2V1×100%,(22)THDi=∑h=29Ih2I1×100%.In order to reduce computational time in this research, bus location is gathered from MLSA with weight summation approach; meanwhile, Pareto optimization is used to determine size of VPF. Normally, Pareto technique will produce a “nondominated” solution which has more than a single solution [31]. Thus, fuzzy satisfying approach is adopted in this research to get the best solution from the “nondominated” solutions. Equation (23) is the weight membership function used in fuzzy approach for all “nondominated” solutions in the fuzzy set. The fitmax in the formula represents the existing fitness value before implementing the filters, while m represents the number of solutions listed in the “nondominated” list. The highest μj value indicates the best solution for this research.(23)μfit_j=1,fitmax-fitjfitmax,0,fitj≤00<fitj<fitmaxfitj≥fitmax,μj=∑15μfit_j∑1m∑15μfit_j.


## 2.1. IEEE 31 Bus with 22 Low Voltage 415 V Bus System

Figures1 and 2 show a single-line diagram for 23 kV IEEE 31 bus radial system with low voltage 415 V buses residential feeder that consist of several CS, respectively. The total number of buses in the network (medium and low voltages) is 449 buses. There are 6 units of CS installed at one residential area which add up to 132 units for overall CS in this distribution system. Line data for low voltage 415 V buses is shown in Table 1.Table 1
Line data for 415 V buses.
Line
Line resistanceR (Ω)
Line reactanceX (Ω)
Bus A
Bus B
a
b
0.0415
0.0145
b
c
0.0424
0.0189
c
d
0.0444
0.0198
d
e
0.0369
0.0165
e
f
0.0520
0.0232
f
g
0.0524
0.0234
g
h
0.0005
0.0002
g
i
0.2002
0.0199
g
j
1.7340
0.1729
f
k
0.2607
0.0260
f
l
1.3605
0.1357
d
m
0.1400
0.0140
c
n
0.7763
0.0774
b
o
0.5977
0.0596
a
p
0.1423
0.0496
p
q
0.0837
0.0292
q
r
0.3123
0.0311
a
s
0.0163
0.0062
DT
0.0000
0.0654Figure 1
Single-line diagram for 449 bus radial distribution system.Figure 2
Single-line diagram for 415 V buses.


## 2.2. Harmonic Load Flow Analysis

Many techniques can be used to perform harmonic analysis such as in [23–25] and the most typical and simplest one is by using current injection analysis method [26]. The electrical parameters that affected the harmonic existence are the line impedances, load variation impedances, and filter impedances values. Equations (1) and (2) represent the line impedance formulation and the impedance for single-tuned filter at harmonic h, respectively, which indirectly caters for resonance impact. Equations (3), (4), and (5) show the impedance formulation for the load used in setting up the harmonic admittance matrix. The harmonic flow is calculated using (6), which consists of harmonic admittance matrix and harmonic current injection as per (7) and (8).(1)Zl,h=Rl+jhωLl,(2)ZFilter,h=RFilter+jhωLFilter-j1hωCFilter,(3)ZL,h=RL+jXL,h,(4)RL=VL2PL,(5)XL,h=VL2hQL,(6)I-h=Y-h·V-h,(7)Y-h=y-11,hy-12,h⋯y-1N,hy-21,hy-22,h⋯⋮⋮⋮⋯⋮y-N1,h⋯⋯y-NN,h,(8)I-h=I-1,hI-2,h⋮I-N,h.


## 2.3. Modelling of CS and Battery

CS is typically divided into three categories which are level 1, level 2, and level 3. The level of the CS is based on the power consumed by the charging station. Level 1 takes 6–8 hours, while level 3 requires 30 minutes for fully charging the EV’s battery. In this research, level 1 CS is selected since this research focuses on residential area which normally involves single-phase supply. The main different in this work is on CS modelling and its behaviour. The CS modelling in this research has considered the worst harmonic value achievable by individual CS based on the battery characteristic and State of Charge (SOC) state. The initial and final (plug-in and plug-out) battery SOC value depended on customer request, while the current produced by CS to the EV battery is calculated based on that SOC state. SOC state will be updated in every 15 minutes(Δt) and can be calculated using (9), while the current produced by CS can be calculated using (10).(9)SOCΔtk+1,i=SOCΔtk,i+ΔtQiIΔtk,i×100,(10)IΔtk,i=4·Δt·Ri·PCSΔtk,i·η·CRiΔtk,i+Voc,i2-Voc,i2Ri·CRiΔtk,i,where Voc,i is open circuit voltage for th node (V), Qi is rated battery ampere hour for the ith PEV (Ah), Ri is battery equivalent internal resistance for the ith node (ohm), CRiΔtk,i is maximum charging rate for the ith PEV (A), SOCΔtk,i is state of charge of the ith PEV at kth time slot (%), SOCΔtk+1,i is state of charge of the ith PEV at next kth time slot (%), IΔtk,i is charging current for the ith PEV at current time slot (A), and PCSΔtk,i is consumed power for the ith PEV (kW).Figure3 shows the open circuit voltage (OCV) for lithium iron phosphate (LFP) batteries based on the SOC of battery [27] that can also be presented using (11). Thus, the relation between SOC and VOCV will influence the total current injected during charging process.(11)VOCV=3.135-0.685-ln⁡SOC0.478-1.342·SOC+1.734·e0.4SOC-1.Figure 3
Relationship between open circuit voltage and SOC for LFP battery.Since there is conversion from AC to DC during charging process, the harmonic impact to the network is considered in this study. The harmonic injection is basically based on SOC state, which depended on battery SOC status and is being modelled based on practical data from [18]. Since there are 132 CS units installed within the 22 low voltage residential areas (6 units for every residential area) it might have significant impact on THDV and THDi value when all CS operated simultaneously. All CS and battery data in the residential area are presented as per Table 2.Table 2
CS and battery data.
Type
Battery capacity (kW)
Charging rate (A)
Power usage (kW)
CS efficiency
1
10
0.052
0.52
0.93
2
15
0.125
1.875
0.93
3
10
0.235
2.35
0.93
4
15
0.15
2.25
0.93
5
10
0.215
2.15
0.93
6
20
0.102
2.04
0.93


## 2.4. Modelling of Load Profile and Charging Station Operation

The load profile for 449 bus system is design based on typical load consumption in residential area during normal day. The peak demand in the system is occurring at the afternoon due to the several usages of electrical appliances, while the lowest demand is at the night when majority of people are asleep. However, the new electrical load CS is normally being charged at the night time. In this study, it is assumed that harmonic produced in the distribution system is coming from CS. Figure4 shows the loading produced by all CS in the distribution system based on design of CS operation as in Table 3. CS operation is designed based on random behaviour and will penetrate at the night hour. Figure 5 shows the summation of original load profile with CS load. Even though most of the CS is operating during low loaded condition at distribution system, it will cause higher impact on the total harmonic distortion as well as power losses.Table 3
Bus
Plug-in
Plug-out
SOCinitial
SOCReq
Type
33
18
64
28
58
1
37
21
64
11
80
2
39
23
65
10
82
3
46
21
64
1
78
4
48
19
63
1
69
5
49
22
65
4
85
6
52
15
61
28
58
1
56
18
61
11
80
2
58
20
62
10
82
3
65
18
61
1
78
4
67
16
60
1
69
5
68
19
62
4
85
6
71
23
69
28
58
1
75
26
69
11
80
2
77
28
70
10
82
3
84
26
69
1
78
4
86
24
68
1
69
5
87
27
70
4
85
6
90
13
59
28
58
1
94
16
59
11
80
2
96
18
60
10
82
3
103
16
59
1
78
4
105
14
58
1
69
5
106
17
60
4
85
6
109
17
60
4
85
6
113
14
58
1
69
5
115
16
59
1
78
4
122
18
60
10
82
3
124
16
59
11
80
2
125
13
59
28
58
1
128
26
69
4
85
6
132
23
67
1
69
5
134
25
78
1
78
4
141
27
82
10
82
3
143
25
80
11
80
2
144
22
58
28
58
1
147
33
76
4
85
6
151
30
74
1
69
5
153
32
75
1
78
4
160
34
76
10
82
3
162
32
75
11
80
2
163
29
75
28
58
1
166
18
64
28
58
1
170
21
64
11
80
2
172
23
65
10
82
3
179
21
64
1
78
4
181
19
63
1
69
5
182
22
65
4
85
6
185
15
61
28
58
1
189
18
61
11
80
2
191
20
62
10
82
3
198
18
61
1
78
4
200
16
60
1
69
5
201
19
62
4
85
6
204
23
69
28
58
1
208
26
69
11
80
2
210
28
70
10
82
3
217
26
69
1
78
4
219
24
68
1
69
5
220
27
70
4
85
6
223
13
59
28
58
1
227
16
59
11
80
2
229
18
60
10
82
3
236
16
59
1
78
4
238
14
58
1
69
5
239
17
60
4
85
6
242
17
60
4
85
6
246
14
58
1
69
5
248
16
59
1
78
4
255
18
60
10
82
3
257
16
59
11
80
2
258
13
59
28
58
1
261
26
69
4
85
6
265
23
67
1
69
5
267
25
78
1
78
4
274
27
82
10
82
3
276
25
80
11
80
2
277
22
58
28
58
1
280
33
76
4
85
6
284
30
74
1
69
5
286
32
75
1
78
4
293
34
76
10
82
3
295
32
75
11
80
2
296
29
75
28
58
1
299
18
64
28
58
1
303
21
64
11
80
2
305
23
65
10
82
3
312
21
64
1
78
4
314
19
63
1
69
5
315
22
65
4
85
6
318
15
61
28
58
1
322
18
61
11
80
2
324
20
62
10
82
3
331
18
61
1
78
4
333
16
60
1
69
5
334
19
62
4
85
6
337
23
69
28
58
1
341
26
69
11
80
2
343
28
70
10
82
3
350
26
69
1
78
4
352
24
68
1
69
5
353
27
70
4
85
6
356
13
59
28
58
1
360
16
59
11
80
2
362
18
60
10
82
3
369
16
59
1
78
4
371
14
58
1
69
5
372
17
60
4
85
6
375
17
60
4
85
6
379
14
58
1
69
5
381
16
59
1
78
4
388
18
60
10
82
3
390
16
59
11
80
2
391
13
59
28
58
1
394
26
69
4
85
6
398
23
67
1
69
5
400
25
78
1
78
4
407
27
82
10
82
3
409
25
80
11
80
2
410
22
58
28
58
1
413
33
76
4
85
6
417
30
74
1
69
5
419
32
75
1
78
4
426
34
76
10
82
3
428
32
75
11
80
2
429
29
75
28
58
1
432
18
64
28
58
1
436
21
64
11
80
2
438
23
65
10
82
3
445
21
64
1
78
4
447
19
63
1
69
5
448
22
65
4
85
6Figure 4
Load produced by CS over time.Figure 5
Load profile before and after CS operation.


## 2.5. Assumptions and Constraints

There are few assumptions and constraints are considered in this research as shown below:(i)
CS can be plugged in/plugged out at any time according to the customer’s request. Customers will input their requested plug-out time and requested final SOCs at the time of plug-in. Once SOC reaches requested SOC, CS will be switched to a standby mode.(ii)
The time slot is 15 min (Δt) or equal to 96 slots for one day.(iii)
The aggregator has access to CS information including their bus locations, harmonic distortion information, charger types, battery sizes, plug-in time, and plug-out time.(iv)
CS are controllable and have variable charging functions. During the charging process, each CS is assumed as a variable active load. The power that is used to charge the battery is based on calculation that is done by the aggregator.(v)
The requested time for each CS must be greater than the minimum charging time required to charge the battery.Since 22 variable filters will be placed in the network, the total unknown variables will become 44, which are 22 locations and 22 optimal sizes. All these parameters will have their own constraints that need to be considered. In general, the parameters can be divided into 2 categories which are as follows:(i)
Filter location: the filter will be placed in low voltage bus to avoid any harmonic injection to upper feeder. There are 22 locations in 449 bus radial distribution system. Furthermore, only one VPF will be placed at every low voltage 415 buses system. The constraints are as follows:(12)ai≤Filteri≤si,i=1…22.(ii)
Filter reactive value: each individual variable filter reactive value is limited to(13)0kVAr≤Qi≤40kVAr,i=1…22.


## 2.6. Modelling of Passive Filter

The main function of a passive filter is to sink the harmonic current that flows in the system based on a selected frequency. The filter impedance will become very low to allow the harmonic to sink. Single tuned filter is the most popular type of filter which is used widely in dealing with harmonic pollution especially in the industrial area [28, 29]. In this research, four units of single tuned filters are considered as one set of filter which can eliminate four frequencies as shown in Figure 6. Equation in [26] is used to calculate capacitor, inductance, and resistance components, respectively, as per (14). Capacitor is calculated based on injected reactive power (Q) and voltage (V) at that bus; meanwhile, inductance and resistance are based on the chosen harmonics (n) that need to be reduced. Four sets of filter will be used to eliminate 3rd, 5th, 7th, and 9th harmonic order in the network.(14)CFilter=Q2πfV2,LFilter=V22n2Q2πf,RFilter=V2nQ2.Figure 6
1 set of filter.


## 2.7. Multiobjective MLSA with Pareto Optimization

In this research, the optimization technique is used to identify optimized locations and sizes for 22 sets of filters in 449 bus radial distribution system. The Lightning Search Algorithm (LSA) is a metaheuristic optimization method inspired by lightning propagation from cloud to ground [30]. The LSA process basically starts with generation of random population before fitness calculation. The worst step leader will be eliminated every 10th iteration before the direction of the step leader is updated. Next, step leader movement is updated based on direction, shape parameter, and scale parameter. Lastly, forking phenomenon is taken place as per (7) for 1% of the better step leader. Since the problem faced in this research is too complicated with many local minima, modification on the existing LSA is needed. In this research, the MLSA is proposed with better convergence for this application compared to existing LSA. MLSA basically is the improvement of LSA method which is proven beneficial for this research.The placement and sizing of variable filter will be based on five parameters, which are maximumTHDV, maximum THDi, THDi summation, average THDi, and apparent power losses for overall system. The three parameters that involve THDi are to ensure that the THDi at all buses is within acceptable range. These parameters are then normalized to get the most accurate solution as defined by (15)–(19). The weight summation with normalized fitness function is shown in (20) after several manual trials are done in order to get the best optimal result. Last but not least, basic equations to calculate THDi and THDV are as shown by (21) and (22).(15)fit1=max⁡THDV,n=1…449,(16)fit2=max⁡THDi,n=1…30,(17)fit3=∑i=130THDi,(18)fit4=∑i=130THDi30,(19)fit5=∑i=1449si-loss,(20)Fit=0.1fit1+0.2fit2+0.2fit3+0.2fit4+0.3fit5,(21)THDV=∑h=29Vh2V1×100%,(22)THDi=∑h=29Ih2I1×100%.In order to reduce computational time in this research, bus location is gathered from MLSA with weight summation approach; meanwhile, Pareto optimization is used to determine size of VPF. Normally, Pareto technique will produce a “nondominated” solution which has more than a single solution [31]. Thus, fuzzy satisfying approach is adopted in this research to get the best solution from the “nondominated” solutions. Equation (23) is the weight membership function used in fuzzy approach for all “nondominated” solutions in the fuzzy set. The fitmax in the formula represents the existing fitness value before implementing the filters, while m represents the number of solutions listed in the “nondominated” list. The highest μj value indicates the best solution for this research.(23)μfit_j=1,fitmax-fitjfitmax,0,fitj≤00<fitj<fitmaxfitj≥fitmax,μj=∑15μfit_j∑1m∑15μfit_j.


## 3. Methodology

The variable filter is used to minimize the impact of harmonics due to large-scale deployment of CS in the big distribution system. The numbers of filter are based on how severe the harmonic distortion and power losses on the existing system. By placing 132 CS units (6 units for every low voltage residential area) at specific locations, MLSA will be used to determine optimal placement and size for 22 sets of filters (1 set for every low voltage residential area). TheTHDV, THDi, and Sloss will be recorded based on specify CS and filter locations. The overall process to evaluate system performance is as shown in Figure 7.Figure 7
Flowchart to evaluate fitness.There are four modifications on MLSA that are made to improve LSA method. The first modification is on channel time; the original channel time is not suitable when dealing with the problems that have many local minima. The second modification is the updating approach; a new updating approach able to check the forward or backward direction is proposed. The third modification is on the scale parameter; in LSA, it is based on exponential distribution, which will cause the step leader movement to be active at 20% of the early iterations as per (24). Thus, in MLSA, the Laplacian distribution equation is used to increase the movement of the step leader during 30% until 70% iterations. Last but not least, the learning factors which make the movement of the step leader more active is implemented in the fourth modification as per (25).(24)PLi_new=PLi±normrandμi,σi,(25)μi=Rxi-xbest,where σi is the scale parameter at ith iteration (1/3)e-(t-50/15), R is learning factor (2.0 in this research), PLi_new is new position in lead projectile at ith iteration, PLi is position in lead projectile at ith iteration, μi is shape parameter at ith iteration, xi is the position at ith iteration, and xbest is the best individual during minimum.Figure8 shows the flowchart for MLSA optimization process using weightage summation approach for this research. In this process, simulation involves 500 iterations and 50 populations. Next, the bus location results will be used in the next stage by the MLSA-Pareto-Fuzzy combination technique. Figure 9 shows the flowchart for MLSA optimization process with Pareto-Fuzzy technique. In this stage, Pareto technique will be adopted in finding the best solution, while fuzzy approach will be used at the end of process to find the best solution among nondominated solutions.Figure 8
Flowchart for MLSA method to get minimum fitness with weightage summation approach.Figure 9
Flowchart for MLSA method to get minimum fitness with Pareto-Fuzzy technique.


## 4. Result and Discussion

The simulation on this research was done based on the scenario where 132 CS units are installed in low voltage 415 V buses for 24 hours (96 states, 1 state = 15 minutes). The schedule for CS operation can be referred to in Table3. The THDV and THDi for all medium voltage buses and lines are shown in Figures 11 and 12, respectively, for sampling time at 30 and 63. Sampling time 63 is chosen due to the worst harmonic recorded in 24 hours, while sampling time 30 is chosen due to the lowest harmonic injection. From the result, without VPF, maximum THDV for sampling time 30 is recorded at bus 15 with 0.8750% and maximum THDi is recorded at lines between buses 26 and 27 (line number 27) with the value 0.6780%. For sampling time 63, maximum THDV is 0.9032% at bus 15, while maximum THDi is recorded at lines between buses 12 and 13 (line number 17) with the value 1.4166%. Apparent losses for both sampling times are 90.3 kVA and 40.4 kVA, respectively.Next, 22 sets of VPFs are placed at low voltage 415 V buses using MLSA optimization with multiobjective function. The optimal location and size for all filters are tabulated in Table4, while the parameter value is shown at Table 5. From the results in Table 5, there is a significant reduction for maximum THDV from 1.258004% to 0.74478% and maximum THDi from 2.019816% to 1.16256% after MLSA get the optimal placement for VPF. However, the total apparent losses in the network are increased. Thus, a Pareto-Fuzzy approach is used to cater this issue at the next step.Table 4
Bus location and filter size after 500 iterations.
Number
Bus location
Filter size (kVAr)
(
1
)
42
28.007
(
2
)
56
16.494
(
3
)
85
22.724
(
4
)
105
1.018
(
5
)
119
1.777
(
6
)
145
24.54
(
7
)
161
28.731
(
8
)
168
27.679
(
9
)
193
29.999
(
10
)
217
8.927
(
11
)
222
2.109
(
12
)
254
1.024
(
13
)
275
27.391
(
14
)
282
17.989
(
15
)
305
14.646
(
16
)
322
24.733
(
17
)
337
26.043
(
18
)
367
24.758
(
19
)
387
29.863
(
20
)
407
21.946
(
21
)
424
26.069
(
22
)
447
28.156Table 5
Result using MLSA with multiobjective function.
Function
Existing
After multiobjective function
MaxTHDV
1.258004
0.74478
MaxTHDi
2.019816
1.16256
T
H
D
i average
0.474635
0.322456
T
H
D
i sum
14.71369
9.996144
S
l
o
s
s
0.04175
0.043149
Fitness
1
0.75612The results from Pareto consist of multiple set solutions where the possible solution is the best solution as long as one of their objectives functions is dominant compared to other solutions. Fuzzy approach will be used to find the best among solutions in Pareto set. In this study, the solution consists ofTHDV, THDi (minimum, average, and summation), and power loss. From the analysis, the total Pareto result for sampling time at 63 in MLSA is 137 sets. The fuzzy approach is used in order to get the best single solution. The highest value of weight membership function in fuzzy approach will be chosen as the best solution in this research. Tables 6 and 7 show the five best solutions with the highest μj for sampling times 63 and 30. From the result, the best sizing for all filters is shown in Table 8 at sampling time 63. The size of filter to cater harmonic is different between times 30 and 63 which indicates that the size of passive filter is very crucial in this study. The THDV and THDi for buses and lines, after the filter implementation as per MLSA with Pareto-Fuzzy best solution, are shown in Figures 10 and 11. From the results, it shows significant reduction of THDV at all medium voltage buses, while majority of THDi also shows reduction. From the best solution, the THDi and THDv for the 449 bus distribution system have been reduced effectively after filter installations. The results for all 96 slots are shown in Figures 12–16. Based on Figures 12 and 13, maximum THDV and THDi show reduction for all states which indicate that the system is better compared to previous value, while apparent losses at Figure 16 show the capability of this method to also reduce losses together with harmonic. Next, Figures 14 and 15 also show reduction which assist THDi improvement in overall system indirectly. From the result, implementation of VPF with changes every 15 minutes is able to reduce harmonic for 24 hours. Other than that, it is obvious that the implementation of MLSA with Pareto-Fuzzy is able to identify the best solution that provides the greatest improvement to the overall system.Table 6
The five best solutions based on the highestμj using fuzzy approach for 137 sets of solutions (sampling time = 63).
Number
Fitness parameter
μ
f
i
t
1
μ
f
i
t
2
μ
f
i
t
3
μ
f
i
t
4
μ
f
i
t
5
μ
j
f
i
t
1
f
i
t
2
f
i
t
3
f
i
t
4
f
i
t
5
(
1)
0.73190
5
1.042406
0.308344
9.558658
0.040295
0.350982
0.26416
0.196401
0.196401
0.004568
0.009541
(
2
)
0.715371
1.087904
0.308852
9.574414
0.040469
0.365644
0.232043
0.195076
0.195076
0.000013
0.009308
(
3
)
0.691356
1.213682
0.299624
9.288355
0.040021
0.386939
0.143255
0.219125
0.219125
0.011102
0.00923
(
4
)
0.691356
1.213682
0.299624
9.288355
0.040021
0.386939
0.143255
0.219125
0.219125
0.011102
0.00923
(
5
)
0.691356
1.213682
0.299624
9.288355
0.040021
0.386939
0.143255
0.219125
0.219125
0.011102
0.00923Table 7
The five best solutions based on the highestμj using fuzzy approach for 10 sets of solutions (sampling time = 30).
Number
Fitness parameter
μ
f
i
t
1
μ
f
i
t
2
μ
f
i
t
3
μ
f
i
t
4
μ
f
i
t
5
μ
j
f
i
t
1
f
i
t
2
f
i
t
3
f
i
t
4
f
i
t
5
(
1)
0.62029
0.452838
0.093176
2.888464
0.089405
0.34565
0.334048
0.081789
0.081789
0.009764
0.124558
(
2
)
0.647588
0.41204
0.095228
2.952076
0.089555
0.316852
0.394046
0.061568
0.061568
0.0081
0.122965
(
3
)
0.639053
0.421681
0.097018
3.007568
0.089525
0.325856
0.379868
0.043927
0.043927
0.00843
0.117106
(
4
)
0.639053
0.421681
0.097018
3.007568
0.089525
0.325856
0.379868
0.043927
0.043927
0.00843
0.117106
(
5
)
0.652162
0.410311
0.10102
3.131623
0.089686
0.312028
0.396589
0.004492
0.004492
0.006651
0.105752Table 8
Filter size based on the best solution after fuzzy at sampling times 30 and 63.
Number
Bus location
Filter size (MVAr)
(Time = 30)
(Time = 63)
(
1
)
42
0.017278
0.025038
(
2
)
56
0.02187
0.005959
(
3
)
85
0.010198
0.019481
(
4
)
105
0.029332
0.011044
(
5
)
119
0.014816
0.001792
(
6
)
145
0.019336
0.02311
(
7
)
161
0.002469
0.026056
(
8
)
168
0.023573
0.022019
(
9
)
193
0.003353
0.018806
(
10
)
217
0.018943
0.008713
(
11
)
222
0.019898
0.016408
(
12
)
254
0.024081
0.002977
(
13
)
275
0.0176
0.018999
(
14
)
282
0.020249
0.029898
(
15
)
305
0.01956
0.022548
(
16
)
322
0.016437
0.015687
(
17
)
337
0.014686
0.01437
(
18
)
367
0.026797
0.016679
(
19
)
387
0.022292
0.012124
(
20
)
407
0.011888
0.016332
(
21
)
424
0.001049
0.011097
(
22
)
447
0.020735
0.028984Figure 10
T
H
D
V at medium voltage for times 30 and 63 with and without passive filter.Figure 11
T
H
D
i at medium voltage for times 30 and 63 with and without passive filter.Figure 12
T
H
D
v maximum for 24 hours with and without passive filter.Figure 13
T
H
D
i maximum for 24 hours with and without passive filter.Figure 14
T
H
D
i average for 24 hours with and without passive filter.Figure 15
T
H
D
i summation for 24 hours with and without passive filter.Figure 16
Apparent power losses for 24 hours with and without passive filter.Based on the collective result gathered from all 96 simulations that represent 24 hours, minimum and maximum values of VPF at dedicated bus are acquired. Table9 shows the minimum and maximum VPF required to solve harmonic problem after EV is installed.Table 9
Minimum and maximum VPF size based on 96 simulations.
Number
Bus location
Filter size (kVAr)
Minimum
Maximum
(
1
)
42
1.021
28.916
(
2
)
56
1.282
28.659
(
3
)
85
0.619
28.838
(
4
)
105
0.637
29.498
(
5
)
119
0.692
20.457
(
6
)
145
0.855
29.559
(
7
)
161
0.204
29.944
(
8
)
168
0.5
29.769
(
9
)
193
0.725
18.806
(
10
)
217
0.724
29.936
(
11
)
222
0.228
29.927
(
12
)
254
0.319
29.996
(
13
)
275
0.171
28.969
(
14
)
282
0.409
29.898
(
15
)
305
0.654
29.805
(
16
)
322
0.318
29.859
(
17
)
337
0.84
29.565
(
18
)
367
0.176
29.56
(
19
)
387
1.385
29.576
(
20
)
407
0.235
26.327
(
21
)
424
1.049
28.326
(
22
)
447
0.244
29.922


## 5. Conclusion

Optimal placement and sizing of VPF in the distribution system can avoid harmonic injection to the medium voltage network as well as reducing apparent power losses. This research has shown that MLSA with multiobjective function and Pareto-Fuzzy are able to determine appropriate filter location and size to reduce harmonic distortion and apparent losses in 449 bus system with 132 units of CS. Furthermore, the placement of twenty-two sets of four unit single-tuned filters is able to reduce the four harmonic orders. Based on simulation run for every 15 minutes, the proposed technique is able to improve maximumTHDV, maximum THDi, and Sloss up to 39.14%, 52.5%, and 2.96%, respectively. The study is very important for future distribution grid that might have many CS in the network. For future work, the study can be expanded to other types of passive filter able to have better impact on the system such as C-type filter and others.



---

*Source: 6587493-2018-02-15.xml*