Estimation and Mitigation of Electrostatic Interferences ...

Estimation and Mitigation of Electrostatic Interferences on Metallic Pipeline by HV Overhead Power Line using Differential Evolution Algorithm (Full text in English)

1,2

Rabah DJEKIDEL1, Sid Ahmed BESSIDEK2 Laboratory for Analysis and Control of Energy Systems and Electric Systems, LACoSERE, Laghouat University, Algeria

Abstract The main purpose of this paper is to analyze the electrostatic interferences of HV power lines on aerial metallic pipeline, also to determine the optimal location of the pipeline from the power line centre, which produces a lower induced alternating voltage. Modeling these interferences is typically done for safety reasons, to ensure that the induced voltages generated are not dangerous to people in contact with the pipe, for equipment connected to the pipeline and cathodic protection system. Very often, mitigative measures are needed to reduce these voltages to the safety threshold. in this paper we propose a simulation model based on approach combining Differential Evolution (DE) algorithm and charge simulation method (CSM). The obtained simulation results were compared with those obtained from an analytical approach, a good agreement was obtained. Keywords: Charge Simulation Method (CSM), Differential Evolution (DE), H-V power line, aerial pipelines, electrostatic coupling, mitigation, shielding wires Received: March, 10, 2016

1. Introduction

The aerial metallic pipelines that share common rights-of-way with high-voltage power lines are subjected to the influence of the electrostatic and electromagnetic interferences created by the electric and magnetic fields of the high voltage power line, which induces voltages in the neighboring metallic pipelines. In some cases, the voltage may raise to a high enough level to be harmful to the operating personnel coming into contact with it, pipeline-associated equipment, cathodic protection systems and the pipeline itself [1-4]. Consequently, the induced voltages must be kept to reduced levels that are safe for personnel and integrity of the pipeline. Therefore it is important to evaluate the electric interferences between the equipments for security reasons. For electrostatic coupling, the Charge Simulation Method (CSM) is often used due to its favorable characteristics for the calculation of electric fields in the high-voltage equipments. CSM is a numerical method, it is the more understandable and it is very easy to use, ease of programming, this method has turned into a very powerful and effective tool to calculate the electric fields of high-voltage devices and high voltage insulation systems.CSM utilizes a number of fictitious charges to equivalently express the analytical expressions of the potential and the electric field. Therefore, the only main difficulty of CSM is that the values and locations of the simulating charges are difficult to be exactly determined for the most accurate results with this method, to solve the constraint optimization of electric charges in this system, we appeal to the heuristic optimization methods, as genetic algorithm (GA); simulated annealing (ASA); Particle Swarm Optimization (PSO) and Differential Evolution (DE). Differential evolution (DE) is a evolutionary algorithms developed by Rainer Stron and Kenneth Price in 1995 for optimization problems. It is a population-based direct search algorithm for global optimization capable of

handling non-differentiable, non-linear and multi-modal objective functions [5-8]. In the present work, we propose a modeling and simulation of electrostatic coupling using a hybrid method, the Differential Evolution (DE) algorithm with the charge simulation method (CSM), in order to optimize an optimal location of the pipeline so that the induced voltage is close to zero. A program developed using MATLAB 10 was used through this work. The validity of these simulation results is provided by a comparison with results obtained by the impedance Matrix analysis. 2. Capacitive Coupling from Power Lines to Pipelines Only pipelines installed above earth are able to capacitive coupling, if the pipeline is located near a high voltage power line, it can undertake a large voltage to ground. The voltage is due to the capacitances between the power line conductors and the pipe, and between the pipe and ground, which form a capacitive voltage divider [1-4]. this is illustrated in Figure 1.

Figure 1. Capacitive coupling from a power line to a pipeline

where Vc: transmission line voltage; VP: induced voltage in pipeline;

ELECTROTEHNICĂ, ELECTRONICĂ, AUTOMATICĂ (EEA), vol. 64 (2016), nr. 3

84

C1:capacitance between power line and pipeline; C2: capacitance between pipeline and ground.

R= {r1 if k=i, r2 if k=j} hcon: heights of conductors and pipeline above ground; xcon: horizontal coordinates of conductors and pipeline.

3. Capacitive Interference Calculation 3.1 Charge Simulation Method (CSM)

The charge simulation method is used to compute the electric field stress and the induced voltage of the highvoltage transmission line, due to capacitive coupling on metallic pipeline. In this method, the conductors and the pipeline are simulated by a number of infinite line charges, placed inside the conductor/pipeline, around a fictitious cylinder of radius r2 [4,9]. In most problems concerning the resolution by the charge simulation method there is a symmetry plane (usually the earthed infinite plane at y=0 keeps at zero potential), which is taken into account by introducing image charges. The total number of the unknown charges in the studied problem requires the same number of the numerical equations for simultaneous solution. Therefore, same number of boundary points is selected on the surface of conductors and pipeline, Dirichlet boundary conditions are applied to these surfaces [10,11]. The charges are computed from the equation: nt

Vi = ∑ Pij .q j

(1)

where Pij: The potential coefficient matrix; qj: The column vector for values of the unknown charges; Vij: The potential of the contour points; nt: is the total number of fictitious charges, for the phases conductors, earth wire and pipeline. The potential coefficient depends only on the type of the charge and the relative distance between conductor (i) and (j). for an infinite length of charge type in a 2D system, the potential coefficient is given in equation below.

  Pij = ln 2 π ε 0  

(x (x

i

i

2 2  − x j ) + ( yi + y j )  2 2  − x j ) + ( yi − y j )  

(2)

where: (xi, yi): coordinates of contour points; (xj, yj): coordinates of simulation charges. Solving equation (1), the values of charges, the potential and electric field at any point in the region outside the conductors can be calculated easily using the superposition principles. The general forms of coordinates of contour points and the simulation charges as shown in Figure 2 are given by the by the following equation [4,10,11]:

2.π  xk = xcon + R.cos . ( k − 1)  nk   2.π yk = H con + R.sin . ( k − 1)   nk where:

r1

: Contour points : Simulation charges : Check points r1: Real radius of the conductor/pipeline r2: Fictitious radius of the conductor/pipeline Figure 2. Arrangement of the simulation charges and the contour points (conductor/ Pipeline)

For a cartesian coordinate system, the magnitude of the total electric field at the desired point is calculated by the summation of the components. 2

j =1

1

r2

(3)

 n   n  E ti =  ∑ fij . q j  +  ∑ fij . q j   j =1   j =1  t

2

t

(4)

where f xi , f yi are the field intensity coefficients between the contour points and the simulation charges (qj). The induced voltage Vind on the pipeline due to the simulation charges (qj) located at (xj , yj) is calculated using equation (5) given below [1,4].

Vind

  = q .ln ∑ j  2 π ε 0 j =1   1

nc

(x − x ) +( y + y ) (x − x ) + (y − y ) 2

j

2

j

2

j

j

2

    

(5)

where nc is the number of fictitious charges that simulate the phases conductors and the earth wire. If a person touches a pipeline whose voltage is Vind, the discharge current that would flow through his body is given by [1,4]:

I shock = j.C p .L p .

dVind dt

(6)

where: Lp is the length of the pipeline exposed to capacitive coupling. Cp is the pipeline’s capacitance. 3.2 Capacitive Interference Reduction

In order to minimize the capacitive interference on the pipeline, two means are proposed. 3.2.1 Earthing of the Pipeline Reduction of the discharge current magnitude to a value lower than the admissible body current as defined in steady state by the international regulations to be equal to 10mA for adult males requires earthing the pipeline with an appropriate earthing resistance RE to reduce the current below the admissible limit must be lower than [1,4,12]:

ELECTROTEHNICĂ, ELECTRONICĂ, AUTOMATICĂ, vol. 64 (2016), nr.3

RE p

Rbody

(7)

β −1

where Rbody is the body resistance, ß is the ratio ß = ( I shock / I adm ) . According to the American standard IEEE 80:2000, the overall resistance of the human body is usually taken equal to 1000 Ω [12,13]. 3.2.2 Installation of Shielding Wires This method is to install the shield wires (passive or active) having the role of screens, placed between the phases conductors and the ground, with respect to the minimum clearance between conductors of various voltages. To calculate the reduction of the induced voltage obtained by installation of the screens conductors, we use charge simulation method and we incorporate the characteristics of conductors position in the matrix of the potential coefficients, then we calculate the charges carried by all the conductors, by taking care to assign to the screens conductors the potential V = 0 for passive mitigation,or a percentage of the phase voltage for active mitigation [14-17]. 4. Differential Evolution (DE)

Differential Evolution (DE) is a heuristic optimization algorithm; it belongs to the class of evolutionary algorithms. Differential Evolution was first introduced by Kenneth V. Price and R. Storn in 1996. DE was mainly developed to optimize real parameter and real valued functions. DE can optimize a problem and iteratively try to improve the quality of candidate solution based on the specified measure of quality. The principle of this method is to create a new individual, by adding the weighted difference between two individuals to a third. Considering an N-dimensional problem with a population of N individuals evolving with each generation t, according by three operators conceived up as follows: Mutation: For a given parameter vector xi, G randomly select three vectors xr1,G,xr2,G and xr3,G , such that the indices i, r1, r2 and r3 are distinct. where G is the generation number. Add the weighted difference of two of the vectors to the third, following this formula:

(

Vi ,G +1 = xr1 ,G +F xr2 ,G − xr3 ,G

)

v j ,i,G+1 if rand j,i ≤ CR or i u j,i,G+1 =   x j,i,G if rand j,i f CR and j ≠ Irand i =1, . . . ,N; j = 1, . . .,D

where rand j ,i ~ U [ 0,1] , I rand is a random integer from

[1,2, . . .,D ] ; I rand

(

)

(

ui ,G +1 if f ui ,G +1 ≤ f xi ,G xi ,G +1 =   xi ,G otherwise

) (10)

where i varies from 1 to N. Using DE algorithm to CSM requires np solutions; each solution consists of nc simulating charge. Initially we have to use equation 1 twice for all the solutions individually, the first to obtain the simulation charges of each solution whereas the second to calculate the new potential at nc check points. The error resulting from each solution during each iteration must be calculated. The error of each solution at the end of a given iteration compared with the error of the same solution through the previous iterations determines the best location of the simulation charge of this solution. Comparing the errors of all the solutions at the end of a given iteration gives the best location of the simulation charges among all the expected solution. The iterative technique must be continued until the potential at all the check points becomes within the acceptable limits depending on the required accuracy (Mutation, Crossover and selection continue until some stopping criterion is reached) [18-23]. The objective function used for the error (fitness function) is very simple and has the form given in the following equation. OF =

1 nt

nt

∑ i =1

Vvi − Vci .100 Vvi

(11)

where: Vvi is the exact potential to which is subjected the conductor and Vci is the actual voltage of the check points. To determine a solution, which minimizes the discharge current through a person in contact with the pipeline, we can used, the objective function of the form given below. OF = −

(9)

ensures that v j ,i ,G +1 ≠ x j ,i ,G .

Selection: The target vector xi,G is compared with the trial vector vi,G+1 and the one with the lowest function value is admitted to the next generation.

(8)

The mutation factor F is a constant from [0, 2]; Vi ,G +1 is called the donor vector. Crossover: The crossover incorporates successful solutions from the previous generation, the trial vector ui,G+1 is developed from the elements of the target vector,xi,G, and the elements of the donor vector, vi,G+1. Elements of the donor vector enter the trial vector with probability CR.

85

I shock ( dp ) − I shock ( dp + SW ) I shock ( dp )

.100

(12)

where: Ishock (dp) is the discharge current in the pipeline at the given initial position; Ishock (dp+SW) is new discharge current in the pipeline at the initial position given with the implementation of shield wires. To determine the optimum location of the pipeline near overhead power lines; this position corresponds to the minimum of the resultant induced voltage. This is based on hybrid method applied to the simulation presented below. This can be obtained using the objective function in the following form: OF = −

Vind ( max ) − Vind ( x p ) Vind ( max )

.100

(13)

ELECTROTEHNICĂ, ELECTRONICĂ, AUTOMATICĂ (EEA), vol. 64 (2016), nr. 3

86

where: Vp(max) is the maximum induced voltage produced by the configuration of the power line; Vp(xp) is the desired minimum induced voltage according to the separation distance of the pipeline.

Figure 3. 220 kV Single circuit lines in vertical configuration with an above pipeline

The geometrical data of the overhead line circuit and pipeline are shown in Figure 4.

5. Impedance Matrix analysis

This technique can be used to validate the simulation results. This method is based on the principle of self and mutual potential coefficients of the electric system [1, 4, 13]. The resulting matrix of impedances for the system (three-phase conductors, earth wires and metal pipeline) at steady state conditions is given below. Vc   Pc Pcp Pcg    V p  =  Ppc Pp Ppg    Vg   Pgc Pgp Pg

  qc      . qp       qg 

(14)

where: c, p and g represent respectively the power lines, phase conductors, pipelines and earth wires. The earthed wires are now eliminated by substituting Vg =0 in equation (14), giving: ' ' Vc   Pc Pc p  = ' V p   Ppc Pp'

 q  . c    q p  

(15)

with:

(16)

For an insulated pipeline, (Vp=0) and from equation (15), the pipelines voltages to earth due to capacitive coupling with the power lines are given by: −1

'   '  V p  =  Ppc     .  Pp  .Vc 

The pipeline is parallel to the axis of the power line at a distance of 25 m, it has an outer radius of 0.3m and its height above ground is 1 m. The length of parallel exposure of the pipeline and power line is 10 km. The three-phase currents on the power line have been assumed under balanced operation with the magnitude of 250 A. The earth is assumed to be homogeneous with a resistivity of 100Ω.m, the nominal system frequency f =50 Hz. 6. Results and discussions

Pc' = Pc − Pcg . Pg− 1 . Pgc

  Pcp' = Pcp − Pcg . Pg− 1 . Pgp   ' Ppc = Ppc − Ppg . Pg− 1 . Pgc   Pp' = Pp − Ppg . Pg− 1 . Pgp 

Figure 4. Evolving process of PSO algorithm with optimum parameters

(17)

where Vc are the known phase voltages to earth of the power lines. Consider for the case study a single circuit transmission line a 220 kV, with one earth wire and an above ground insulated metal pipeline in the vicinity (Figure 3).

To choose the better parameters in this method, initially we did fix the limits of these parameters in intervals according to publications in the literature that combine this method with evolutionary algorithms, we have carried out many operations randomly, the parameters which made the fitness value be the smallest are chosen to use in this paper. In the present study, the optimum parameters of this evolutionary algorithm (DE) are determined as follows: population size of NP = 10, generation number G = 200, step size F = 1.5 and crossover probability of CR = 0.95. Table 1 contains intervals search parameters for algorithm. Table 1. Range of the initial value Phase conductor Ground wire Pipeline

Simulation charges Fictitious radius [m] [2-20] [0.01-0.07] [2-20] [0.001-0.01] [2-40] [0.1-0.28]

The objective function values (OF) given in equation (11) varies with the iteration number as shown in Figure 4 (supra). We see from this figure, the continuous and rapid change of the objective function value, which is decreasing as the number of iterations increases, the results were obtained after 167 iterations. The application results in values that finally converge to the optimum values are detailed in Table 2. Table 2. Optimum values of CSM

Phase conductor Ground wire

Fictitious charges Fictitious radius number [m] 8 0.061036

(OF) value 2.126e-016

18

0.009923

ELECTROTEHNICĂ, ELECTRONICĂ, AUTOMATICĂ, vol. 64 (2016), nr.3 Pipeline

18

0.208117

The simulation results are shown in the Figures 5 and 6, where it becomes obvious that the algorithm converges rapidly to these values. 0.25

For a variable separation distance of pipeline along the corridor (right of way), the electric field profile is presented in Figure 8, the maximum electric field is obtained nearly under the side conductor about 7 m, as it moves away on either side of this point the electric field intensity decreases gradually with the separation distance where it becomes almost negligible far from the power line centre. 7

R(phase) R(wire) R(pipeline)

0.15

6

Electric Field [KV/m]

Fictitious radius

0.2

0.1

0.05

0

87

0

20

40

60

80

100

120

140

160

180

5

4

3

2

X: 25 Y: 1.22

200

Iteration number

1

Figure 5. Convergence of the optimum values for fictitious charges radius (rc, rg, rp)

0 -50

-40

-30

-20

-10

0

10

20

30

40

50

Pipeline position from the power line center [m]

Figure 8. Electric field stress at the pipeline’s surface

40

n(phase) n(wire) n(pipeline)

Induced voltages on the pipeline located at different distances from the midpoint of the pylon have been calculated and the result is given in Figure 9.

30 25

2000 20

1800 15

1600

10 5 0

0

20

40

60

80

100

120

140

160

180

200

Iteration number

Figure 6. Convergence of the optimum values for fictitious charges number

Figure 7 shows the lateral profile of electric field distribution at (y = 1 m) above the ground, with and without the presence of a metal pipeline.

Induced voltage [V]

Fictitious charges number

35

1400 1200 1000 800 600

X: 25 Y: 375.3

400 200 0 -50

-40

-30

-20

-10

0

10

20

30

40

50

Pipeline position from the power line center [m]

Figure 9. Induced voltage on the pipeline at 1 m above the ground

2 Without pipeline 1.8

With pipeline

Electric field [KV/m]

1.6 1.4

X: 25 Y: 1.22

1.2 1 0.8 0.6 0.4 0.2 0 -50

-40

-30

-20

-10

0

10

20

30

40

50

Lateral distance (x) [m]

Figure 7. Electric field profile at 1 above ground with and without the pipeline

The original profile of the electric field is symmetrical at a distance of 7 m near pylon suspension. It is clear from the graph that the presence of the pipeline has a significant effect on the maximum electric field,the electric field is subjected to a considerable increase in the area of the location of the pipeline, because the charge induced on the pipeline has been increased significantly, the distortion of the electric field in the presence of the pipeline took place at the position where the pipeline is located.

Similarly to the electric field, the induced voltage has a maximum value at the point of symmetry and then decreases gradually as one move away from symmetry. We can also observe that the induced voltage becomes almost negligible at a critical distance about 37 m. It is suggested that the pipeline could be located close to the critical distance so that the induced voltage would be close to zero. The discharge current through the body of a person that touches pipeline located at different distances from the midpoint of the pylon is shown in Figure 10, the discharge current in this example is 35 (mA), this value is considered unacceptable from personnel safety viewpoint.

ELECTROTEHNICĂ, ELECTRONICĂ, AUTOMATICĂ (EEA), vol. 64 (2016), nr. 3

88

96.1695

180 160

96.169

140

96.1685

120

Objective Function

Discharge current [mA]

For (03) Shield wires

100 80 60 X: 25 Y: 35

40

For (04) Shield wires

96.168 96.1675 96.167 96.1665 96.166

20

96.1655 0 -50

-40

-30

-20

-10

0

10

20

30

40

50

Pipeline position from the power line center [m]

96.165

0

20

40

60

For values of currents greater than safety curent limit, which is recommended by the standard equal to 10 mA. A protective measure must be implemented is needed, simply connect the pipeline to ground through a suitable resistance calculated according to the equation (7). Resistance of grounding (Earthing) the pipeline as a function of its horizontal proximity distance is shown in Figure 11. 4

x 10

100

140

160

180

200

We can see in this figure that the algorithm increase quickly to very good solutions for the disposition of shielding cables. The objective is to maximize the function, in order to maintain a current value lower compared to the value intiale. After optimizing the arrangement installation of shielding wires, we evaluate the effectiveness of the shielding in Figure 13. 180

Original Shock current With (03) passive shield wire

160

2

With (03) active shield wires With (04) passive shield wires With (04) active shield wires

Discharge current [mA]

140

1.5

1

0.5 X: 25 Y: 400

0 -20

120

Figure 12.Evolving process of DE algorithm with optimum parameters

2.5

Earthing resistance [Ohms]

80

Iteration Number

Figure 10. Intensity of shock current flowing in the human body

120 100 80 60

X: 25 Y: 35

40

-10

-5

0

5

10

15

20

25

30

35

Pipeline position from the power line center [m]

Figure 11. Calculation of the earthing resistance

We can also install metal cables between the lowest phase and the ground, these conductors act as electrostatic screens, the location must comply a minimum safety distance between two electric conductors carried to different potentials,for a 220 kV power line, the required distance is 4.88 meter. the applied voltage to conductors is taken a zero for passive shield wires and a percentage of 1 % of the line voltage for the active shield wires. The objective function values given in equation (12) , as a function of iterations for optimum geometric parameters of shield wires passive and active is shown in Figure 12.

0

X: 25 Y: 9.769

X: 25 Y: 7.402

20

-15

0

5

10

15

20

25

30

35

40

45

50

Lateral distance (x) [m]

Figure 13.Lateral discharge current profile with and without shield wires

We see a significant reduction of shock current, it can be seen from this figure that the increase of shield wire to four cables ensure appropriate reduction in peak current, also the current in the original location of the pipeline, so that the resulting values become less than 10 mA, which is the value recommended by the standards. . One can see in this figure, at a lateral distance, the shielding has no positive effect because it increases very slightly the values of the current, but it does not have importance ,since the original values and the resulting values are very small. Figure 14 show the variation of the objective function given in equation (13), with number of iterations.

ELECTROTEHNICĂ, ELECTRONICĂ, AUTOMATICĂ, vol. 64 (2016), nr.3 99.9

1800

99.8

1600

89

Impedance Matrix method CSM+DE

1400

Induced voltage [V]

Objective Function

99.7 99.6 99.5 99.4

1200 1000 800 600

X: 25 Y: 375.3

X: 25 Y: 376.4

400

99.3 200

99.2 99.1

0

-80

-60

-40

-20

0

20

40

60

80

100

Pipeline position from the power line center [m] 0

20

40

60

80

100

120

140

160

180

200

Figure 16.Comparison of the results between the proposed method and the impedance matrix analysis

Iteration number

Figure 14.Evolving process of DE algorithm with optimum parameters

7. Conclusions

The purpose of optimization algorithm is to find the optimal location of pipeline, this figure illustrates the search process conducted by the algorithm, the objective is to maximise the value of the error to achieve to a minimum value of the applied induced voltage on the pipeline according to the search space, and therefore the suitable location of the pipeline. The simulation result for the optimum location is shown in Figure 15, where it becomes clear that the algorithm rapidly converges to a lateral distance of 37.09 meter from the power line centre. 40 39.5

Optimum Position of the Pipeline

-100

39 38.5 38 37.5 37 36.5 36

0

20

40

60

80

100

120

140

160

180

200

Iteration number

Figure 15.Convergence of the optimum value for pipeline’s location

The corresponding values induced (voltage and current) in the distance are respectively 9.04 V and 0.84 mA . It is prefered that the pipeline could be implanted at the distance so that the shock current value would be close to zero. The last step is to validate the simulation results of the proposed method. This is enough to compare these results with those from the analytical method (impedance matrix analysis). We have shown in Figure 16 the lateral distribution of the induced voltage for both methods, both figures are perfectly superimposed, the graph shows a good agreement between the results obtained by the two methods, the maximum relative error is less than 0.3 %.

In this study, a rigorous modeling is used to evaluate the effect of electrostatic interferences. A novel method based on Differential Evolution (DE) with charge Simulation Method (CSM) is proposed to determine the appropriate arrangement of simulating charges inside the transmission power line and pipeline, and then the optimal solution can be used to determine the electric field intensity. From the results, it is clear that the presence of pipeline in the vicinity of power line causes the distortion of the electric field at pipeline surface due to electric charges accumulated. The induced voltage is maximum at separation distance of 7 m from the power line centre, when away from this point, and then gradually decreases with the separation distance. Generally, if the contact current with the pipeline exceeds the authorized limit; therefore, the pipeline must be grounded through adequate resistance, or installing electrostatic screens (passive or active) with optimizing their emplacements using the algorithm optimization, these screens conductors involved in the reduction of the induced voltage and consequently the reduction of the induced currents. The results also analyze the suggested location to implement the pipeline using the Differential Evolution (DE) algorithm. The results presented by the proposed method are compared with results obtained by an analytical approach, the comparison shows a good agreement,which confirms the validity of the proposed method. 8. Acknowledgment The authors wish to thank the head the and staff of Laboratory for Analysis and Control of Energy Systems and Electric Systems, Laghouat University, Algeria.

9. References [1] CIGRE, Guide Concerning Influence of High Voltage AC Power Systems on metallic Pipelines, CIGRE Working group 36.02, 1995. [2] A. Taflove, J.Dabkowski, Mutual Design Considerations for Overhead AC Transmission Lines and Gas Transmission Pipelines, Final Report EPRI EL-904, AGA Cat No. L51278, IIT Research Institute, Chicago, Sept 1978 [3] Gupta A , Abhishek Gupta, A Study on High Voltage AC Power Transmission Line Electric and Magnetic Field Coupling with Nearby Metallic Pipelines, Master of Science in Faculty of Engineering, Bangalore, India, August 2006. [4] Djekidel.Rabah, Contribution to the modeling of electromagnetic disturbance generated by a high voltage power

90

[5]

[6]

[7] [8]

[9]

[10] [11]

[12] [13] [14] [15]

[16]

[17] [18] [19]

ELECTROTEHNICĂ, ELECTRONICĂ, AUTOMATICĂ (EEA), vol. 64 (2016), nr. 3 line in a pipeline, Doctor of Science Thesis, University Amar Telidji, Laghouat, Algeria,2015. K. Balamurugan, V. Dharmalingam, R. Muralisachithanandam, R. Sankaran, “Differential Evolution based Optimal Choice and Location of Facts Devices in Restructured Power System”, World Academy of Science, Engineering and Technology, International Journal of Electrical, Computer, Energetic, Electronic and Communication Engineering Vol.7, No:2, 2013. Sidhartha Panda, “Differential evolutionary algorithm for TCSCbased controller design, Simulation Modelling Practice and Theory”, ScienceDirect, Vol. 17, Issue 10, pp. 1618–1634 , November 2009. Salkoski, I. C., “Design optimization of distribution transformers based on Differential Evolution Algorithms”, Conference Information and Communication Technologies, Macedonia 2012. M.S.Saad , H.I.Jamaluddin , I. Z.Mat Darus, “ Implementation of pid controller tuning using Differential evolutionary and genetic algorithms”, International Journal of Innovative Computing Information and Control,Vol.8,No.11,pp.7761-7779,2012. Gururaj S Punekar, N K Kishore, H S Y Shastry, “Effect of Non Uniformity Factors and Assignment Factors on Errors in Charge Simulation Method with Point Charge Model”, World Academy of Science, Engineering and Technology, Vol.22,pp. 10-22.2008. Y. Kato, A Charge Simulation Method for the Calculation of Two-dimensional Electrostatic Fields, Memoirs of the Fukui Institute of Technology, japon, (1980). Amiri., H. Hadi. , Marich. M, “The influence of sag in the Electric Field calculation around High Voltage Overhead transmission lines”, Electric Insulation and Dielectric Phenomena, IEEE Conference, 2006. Technical Specification, Basic Safety Publication: Effects of current on human beings and livestock, IEC, TS 60479-1, fourth edition ,07.2005. Nasser. D. Tleis, Power systems modeling and fault analysis, Theory and Practice. Elsevier First published 2008. Cigré, Electric and magnetic field produced by transmission systems, Working Group 01 (Interference and Fields) of Study Committee 36, Paris,1980. A. Bendaoud, M. Rahli, Reduction of the Electric Field Intensity under the H.V. Transmission Lines, Laboratory ICEPS (Intelligent Control and Electric Power Systems), University of S.B.A, Revue AMSE, Vol 56, N°4, pp 55-64, 2004. R. M. Radwan, A.M. Mahdy, M. A. Salam, M. M. Samy, “Electric Field Mitigation under Extra High Voltage Power Lines”, Dielectrics and Electric Insulation, IEEE Transactions, Vol.20, Issue: 1, p 54 – 62,2013. R. Djekidel, D. Mahi, A. Ameur, Ouchar, “ Calcul et réduction du champ électrique au voisinage d'une ligne Transport d’énergie HT ”, Acta Electrotehnica journal, vol 1-2, pp. 26-32, 2014. Kelly Fleetwood, “An introduction to differential evolution”, 2004. N. Kovač, D. Poljak, S. Kraljević, B. Jajac, “Computation of maximal electric field value generated by a power substation”,

[20]

[21]

[22] [23]

Boundary Elements and Other Mesh Reduction Methods, Transactions on Modelling and Simulation, Vol 42, 2006. R. Djekidel, D. Mahi, “Capacitive Interferences Modeling and Optimization between HV Power Lines and Aerial Pipelines”, International Journal of Electrical and Computer Engineering (IJECE), Vol. 4, No. 4, pp. 486-497, 2014. Musrrat Ali, Millie Pant, V. P. Singh, “Two modified differential evolution algorithms and their applications to engineering design problems”, World Journal of Modelling and Simulation, Vol. 6, No. 1, pp. 72-80, England,2010. J. Hájek, “Adaptation of an Evolutionary Algorithm in Modeling Electric Circuits”, Acta Polytechnica, Journal of advanced engieneering, Vol 50, No 1,2010. Musrrat Ali, Millie Pant, Ajith Abraham, “Simplex Differential Evolution”, Acta Polytechnica Hungarica, Vol. 6, No. 5, 2009.

10. Biography Djekidel RABAH was born in Laghouat (Algeria), on March 25, 1968. He graduated the University of ENSET Technique in Laghouat (Algeria), in 1991. He received the Magister and Ph.D. degrees in electrical engineering, in 2010, respectively in 2015, from USTO Oran University and Laghouat University (Algeria). He is Professor at the University of Amar Telidji in Laghouat, Department of Electrical Engineering, Faculty of Technology, (Algeria). His research interests concern: electromagnetic interference (EMI) fields, High Voltage Engineering, and Numerical modeling and simulation. Correspondance adresse: Departement of Electrical Engineering, Faculty of Technology, University Amar Telidji (Algeria), E-mail:[email protected] Bessidek SID AHMED was Born in frenda, Tiaret, Algeria, on june 28, 1980. He graduated the University of ibn khaldoun in Tiaret, Faculty of (Algeria), in 2004. He received the Magister and Ph.D. degrees in electrical engineering, in 2008, respectively In 2015 from USTO Oran University (Algeria). He is Professor at the University of Amar Telidji in Laghouat, Department of Electrical Engineering, Faculty of Technology, (Algeria). His research interests concern: High Voltage - Insulation Materials and Dielectric - Numerical modeling and simulation, and optimization techniques of artificial intelligence-Fault diagnosis of electric machines. Correspondance adresse: Departement of Electrical Engineering, Faculty of Technology, University Amar Telidji (Algeria). E-mail: [email protected]