Radial Basis Function Neural Network Application to Measurement and Control of Shunt Reactor Overvoltages Based on Analytical Rules

This paper presents an artificial intelligence application to measure switching overvoltages caused by shunt reactor energization by applying analytical rules. In a small power system that appears in an early stage of a black start of a power system, an overvoltage could be caused by core saturation on the energization of a reactor with residual flux. A radial basis function RBF neural network has been used to estimate the overvoltages due to reactor energization. Equivalent circuit parameters of network have been used as artificial neural network ANN inputs; thus, RBF neural network is applicable to every studied system. The developed ANN is trained with the worst case of the switching angle and remanent flux and tested for typical cases. The simulated results for a partial of 39-bus New England test system show that the proposed technique canmeasure the peak values and duration of switching overvoltages with good accuracy.


Introduction
In high-voltage HV power systems usually power is transmitted through long high-voltage transmission lines.During the low demand periods nights or weekends , excessive reactive power produced by the capacitance of these lines causes a voltage increase over 1.1 p.u. at the high-voltage/medium-voltage HV/MV substations.For the absorption of the surplus reactive power, HV shunt reactors are connected to the receiving end of the transmission lines.Switching of those shunt reactors produce transients that need to be carefully studied and, if required, limited 1-7 .
If the frequency characteristic of the system shows resonance conditions around multiples of the fundamental frequency, very high and weakly damped temporary overvoltages TOVs of long duration may occur when the system is excited by a harmonic disturbance 8-14 .
Overvoltages can be classified as transient overvoltages, sustained overvoltages, harmonic resonance overvoltages, and overvoltages resulting from ferroresonance.Excessivesustained overvoltages may lead to the damage of transformers and other power system equipments.Transient overvoltages are a consequence of switching operations on long transmission lines, or the switching of capacitive devices, and may result in arrester failures.Ferroresonance is a nonharmonic resonance characterized by overvoltages whose waveforms are highly distorted and can cause catastrophic equipment damages 11, 15 .
Overvoltage will put the shunt reactor into saturation, causing core heating and copious harmonic current generation.Circuit breaker called upon to operate during periods of high voltage will have reduced interrupting capability 11, 16 .
In this paper power system blockset PSB , a MATLAB/Simulink-based simulation tool 17, 18 is used for computation of temporary overvoltages.In order to study temporary overvoltages for a large number of possible system configurations, it is necessary to run many time-domain simulations resulting in a large amount of simulation time.A way to limit the overall calculation time is to reduce the number of simulations by applying analytical or knowledge-based rules to discard a number of system configurations before an actual timedomain simulation is carried out.This paper presents the artificial neural network ANN application for estimation of peak and duration overvoltages under switching transients during reactor energization.A tool such as the one proposed in this paper that can give the maximum switching overvoltage and its duration will be helpful to the operator during system restoration.Also it can be used as a training tool for the operators.Results of the studies are presented for a partial of 39-bus New England test system to illustrate the proposed approach.

Study System Modelling
The electrical components of the network are modeled using the MATLAB/Simulink environment 17 .These models should be adapted for the desired frequency range here the frequencies up to f 10f 0 are considered to be sufficient .The generator is represented by an ideal voltage source behind the subtransient inductance in series with the armature winding resistance that can be as accurate as the Park model 19 .Phase of voltage source is determined by the load flow results.Transmission lines are described by distributed line models.The circuit breaker is represented by an ideal switch.The shunt reactor model takes into account the leakage inductance as well as the magnetizing characteristics of the core, which is modeled by a resistance, R m , simulating the core active losses and a saturable inductance, L sat .The saturation characteristic is specified as a piecewise linear characteristic 20 .All of the loads are modeled as constant impedances.

Harmonic Overvoltages during Restoration
This paper concentrates on the estimation of harmonic overvoltages.These are a result of network resonance frequencies close to multiples of the fundamental frequency.They can be excited by harmonic sources such as saturated reactors, power electronics, and so forth.They may lead to long lasting overvoltages resulting in arrester failures and system faults 21, 22 .The sample system considered for explanation of the proposed methodology is a 400 kV EHV network shown in Figure 1.The normal peak value of any phase voltage is 400 √ 2/ √ 3 kV and this value is taken as base for voltage p.u.In the system studies 100 MVA as a base power is considered.Figure 2 shows the switching transient at bus 2 when reactor is energized.
In practical system a number of factors affect the overvoltages factors due to energization or reclosing.In this paper, the following parameters are considered: developed ANN for estimation of overvoltages in every studied system.In Section 6 that developed ANN is tested for 39-bus New England test system, this issue is better understood.Also, a novel method based on worst case condition determination is proposed in Section 4 to reduce time-domain simulations.
Source voltage affects the overvoltage strongly.Figure 3 shows the effect of line length on overvoltage at different source voltage.Figure 4 shows the effect of saturation curve slope on overvoltages at different equivalent inductance.The saturation curve, and especially the L sat , that is, the final slope of this curve, is a key point for the computation of the inrush currents.The reactor manufacturer provides a L sat slope value with a dispersion usually considered of ±20%. Figure 5 shows the effect of shunt reactor capacity on overvoltages at different equivalent resistance.Also, Figure 6 shows the effect of equivalent capacitance on overvoltages at different remanent flux.As discussed above for an existing system the main factors which affect the peak and duration values of switching overvoltage are voltage at reactor bus before switching, equivalent resistance, equivalent inductance, equivalent capacitance, line length, switching angle, reactor capacity, saturation curve slope, and remanent flux.Here it should be mentioned that a single parameter often cannot be regarded independently from the other important influencing factors.
The magnitude and duration of the overvoltages normally does not depend directly on any single-isolated parameter and a variation of one parameter can often alter the influence of another parameter, in other words, there exists an interaction between the various system and breaker parameters.This forbids the derivation of precise generalized rule of simple formulae applicable to all cases 23 .So an ANN can help to estimate the peak and duration values of switching overvoltages generated during reactor energization.An ANN is programmed by presenting it with training set of input/output patterns from which it then learns the relationship between the inputs and outputs.In the next section, an ANN-based approach is described which can give an acceptable solution of switching transients by the help of which an operator can take a quick decision at the time of operation.

Worst-Case Condition Determination for Overvoltages Simulation
Normally for harmonic overvoltages analysis, the worst case of the switching angle and remanent flux must be considered which it is a function of switching time, reactor characteristics and its initial flux condition, and impedance characteristics of the switching bus 20 .Using the worst switching angle and remanent flux, the number of simulations for each case can be reduced significantly.
In order to determine the worst-case switching time and remanent flux, the following index is defined as 4.1 where t 0 is the switching time, φ r is initial transformer flux, and h is harmonic order.This index can be a definition for the worst-case switching angle and remanent flux.Using a numerical algorithm, one can find the switching time and remanent flux for which W are maximal i.e., harmonic overvoltages are maximal .
Figure 7 shows the result of the PSB frequency analysis at bus 2. The magnitude of the Thevenin impedance, seen from bus 2, Zbus2 shows a parallel resonance peak at 174 Hz. Figure 8 shows changes of W index with respect to the current starting angle and remanent flux.Figure 2 shows voltage at bus 2 after reactor switching for the worst-case condition i.e., switching angle 20 • and remanent flux 0.27 p.u. .For temporary overvoltages, the overvoltage duration has to be taken into account in addition to the amplitude 24 .Table 1 summarizes the results of overvoltages simulation for four different switching angle and remanent flux that verify the effectiveness of W index.

Steps of Assessment and Estimation of Switching Transient Overvoltages
The steps for harmonic overvoltages assessment and estimation are as follows: 1 determine the characteristics of shunt reactor that must be energized, 2 calculate the Z ii h at the reactor bus for h 2f 0 , . . ., 10f 0 , 3 calculation of worst switching angle and remanent flux for simulation, 4 run PSB simulation, 5 calculation of overvoltage peak and duration, 6 repetition of the above steps with various system parameters to learning artificial neural network, and 7 testing of artificial neural network with different system parameters.

The Radial Basis Function Neural Network (RBFNN)
Figure 9 shows the structure of the RBF neural network, which comprises of three layers.The hidden layer possesses an array of neurons, referred to as the computing units.The number of such units can be varied depending on user's requirement 25, 26 .Different basis functions like spline, multiquadratic, and Gaussian functions have been studied, but the most widely used one is the Gaussian type.In comparison to the other types of neural network used for pattern classification like back propagation feedforward networks, the RBF network requires less computation time for learning and has a more compact topology.The Gaussian RBF is found not only suitable in generalizing a global mapping but also in refining local features without altering the already learned mapping.Each hidden unit in the network has two parameters called a center ω and a width σ associated with it.The response of one such hidden unit to the network input X, X x 1 , x 2 , ..., x n T is expressed as: where ω k is the center vector for kth hidden unit, σk is the width of the Gaussian function, and denotes the Euclidean norm.The output layer comprises a number of nodes depending on the number of fault types to be classified which perform simple summation.The response of each hidden unit 4.1 is scaled by its connecting weights α's to the output nodes and then summed to produce the overall network output.The overall network output is expressed as: where N indicates the total number of hidden neurons in the network, α mk is the connecting weight of the kth hidden unit to mth output node, and α mo is the bias term for the corresponding mth output neuron.The learning process of the RBFNN involves with the allocation of new hidden units and tuning of network parameters.The learning process is terminated when the output error goes under the defined threshold 27 .

Training Artificial Neural Network
All experiments have been repeated for different system parameters 2000 sets .1000 sets were used to train RBFNN and 1000 sets were used to test RBFNN.RBFANN learned in 73 epochs.After learning, all parameters of the trained networks have been frozen and then used in the retrieval mode for testing the capabilities of the system on the data not used in learning.The testing data samples have been generated through the PSB program by placing the parameter values not used in learning, by applying different parameters.A large number

Case Study
In this section, the proposed algorithm is demonstrated for two case studies that are a portion of 39-bus New England test system, of which its parameters are listed in 28 .The simulations are undertaken on a single-phase representation.In the proposed method, first, studied system must be converted to equivalent circuit of Figure 1.In the other words, values of equivalent resistance, equivalent inductance, and equivalent capacitance are determined and used in trained artificial neural network to estimate overvoltages peak and duration.
Case 1. Figure 10 shows a one-line diagram of a portion of 39-bus New England test system which is in restorative state.The generator at bus 35 is a black-start unit.In order to reduce the   steady state overvoltage of no load transmission line, a reactor is connected at bus 19.When the reactor is energized, harmonic overvoltages can be produced because of its nonlinear magnetization characteristics.First, equivalent circuit of this system is determined and values of equivalent resistance, equivalent inductance, and equivalent capacitance are calculated, that is, this system is converted to system of Figure 1.In this case, values of equivalent resistance, equivalent inductance, and equivalent capacitance are 0.00291, 0.02427, and 2.474 p.u., respectively.For testing trained ANN, values of voltage at reactor bus bus 19 , line length, and shunt reactor capacity are varied and in each state, overvoltage peak and duration values are calculated from trained ANN and system of Figure 10.Table 3 contains some sample results of test data of Case 1.
Case 2. As another example, the system in Figure 11 is examined.In the next step of the restoration, unit at bus 6 must be restarted.In order to reduce the steady state overvoltage of no load transmission lines, the reactor at bus 6 should be energized.In this condition, harmonic overvoltages can be produced.
After converting this system to equivalent circuit of Figure 1, that is, after calculating equivalent circuit seen from bus 5, various cases of reactor energization are taken into account and corresponding overvoltages peak and duration are computed from PSB program and trained ANN.In this case, values of equivalent resistance, equivalent inductance, and equivalent capacitance are 0.00577, 0.02069, and 0.99 p.u., respectively.Summary of few results are presented in Table 4.It can be seen from the results that the ANN is able to learn the pattern and give results to acceptable accuracy.

Conclusion
In this paper, a radial basis function-based method is presented to study switching overvoltages during shunt reactor energization.Also, a new approach is proposed to reduce Mathematical Problems in Engineering  time-domain simulations based on worst switching angle and remanent flux determination.In addition, since equivalent circuit parameters of the network are used as ANN inputs, developed ANN is applicable to every studied system.The results from this scheme are close to results from the conventional method and helpful in predicting the overvoltage of the other case studies within the range of a training set.The proposed ANN approach is tested on a partial 39-bus New England test system.This method omits time-consuming time-domain simulations and it is suitable for real time applications during system restoration.Also it can be used as a training tool for the operators.

Figure 9 :
Figure 9: The structure of RBF neural network.

Table 1 :
Effect of switching time and remanent flux on the maximum of overvoltages and duration of V peak > 1.3 p.u.

Table 2 :
Some sample testing data and output.Results for a sample test data are presented in Table 2. Values in column V PSB are the absolute values of peak voltage at bus 2 calculated by PSB program in p.u. where the V RBF values are the values simulated by trained network.Also values in column T PSB are the values of overvoltage duration calculated by PSB program in second and T RBF values are the values simulated by trained network.The proposed model tested with portion of 39-bus New England test system.Various cases of shunt reactor energization are taken into account and corresponding peak and duration values estimated from trained model.
V : voltage at shunt reactor bus before switching p.u. , R eqv : equivalent resistance p.u. , L eqv : equivalent inductance p.u. , C eqv : equivalent capacitance p.u. , L.L.: line length km , S.R.: shunt reactor capacity MVAR , L sat : saturation curve slope p.u. , error V : voltage error % , and error T : duration time error % . of testing data have been used to check the proposed solution in the most objective way at practically all possible parameters variation.Percentage error is calculated as:

Table 3 :
Case 1, some sample testing data and output.
V : voltage at reactor bus before switching p.u. , L L: line length km , S R: shunt reactor capacity MVAR , error V : voltage error % , and error T : duration time error % .

Table 4 :
Case 2, some sample testing data and output.
V : voltage at reactor bus before switching p.u. , L L: line length km , S R: shunt reactor capacity MVAR , error V : voltage error % , and error T : duration time error % .