Open access peer-reviewed chapter

# Electric Field Features and Its Application for Air Gap Breakdown Voltage Prediction

Written By

Zhibin Qiu and Jiangjun Ruan

Submitted: 13 June 2017 Reviewed: 06 November 2017 Published: 20 December 2017

DOI: 10.5772/intechopen.72230

From the Edited Volume

## Electric Field

Edited by Mohsen Sheikholeslami Kandelousi

Chapter metrics overview

View Full Metrics

## Abstract

Air gap breakdown voltage prediction by mathematical calculations instead of experiments has been a long sought goal in the area of high-voltage engineering. In this chapter, a prediction method is proposed based on the electric field features and support vector machine (SVM). Two sets of electric field features are defined on the shortest interelectrode path of sphere-sphere and rod (sphere)-plane gap to characterize their spatial structures, which can be extracted from the electric field calculation results by finite element method (FEM). A breakdown voltage prediction model is established by SVM, while the input parameters are the electric field features, and the output parameters are −1 and 1, respectively, characterizing withstanding and breakdown of an air gap under the applied voltage. The proposed method is used to predict the power frequency breakdown voltages of IEC standard sphere-sphere air gaps and the switching impulse discharge voltages of large sphere-plane air gaps. The prediction results coincide well with the experimental data, the mean absolute percentage error of the 260 sphere-sphere gaps is within 2% and that of the 16 sphere-plane gaps is 3.2%. The results preliminarily validate the validity and accuracy of the proposed method for air gap breakdown voltage prediction.

### Keywords

• electric field features
• air gap discharge
• breakdown voltage prediction
• support vector machine (SVM)
• sphere-sphere gap
• sphere-plane gap

## 1. Introduction

With the rapid development of numerical computation methods and calculation capability of computers, the multi-physics coupling analysis has been widely applied in structure design and condition assessment of electrical equipment [1, 2, 3, 4]. The distributions of electric, magnetic, stress, and thermal field can be calculated by commercial software, and multi-physics problems in complex structures can be simulated using powerful numerical techniques, so as to guide the optimal design of machines and devices. However, up to now, the insulation design of high-voltage equipment still depends on experimental verification.

The electric field distributions of complex geometries can be simulated very well, that is because electric field can be defined according to Gauss’s law and Maxwell’s relation, with clear constitutive relation and governing equations [5, 6]. Several numerical methods were developed for electric field calculation, such as finite element method (FEM), finite volume method (FVM), boundary element method (BEM), and those combining the advantages of different methods, like the control volume-based finite element method (CVFEM) [7, 8, 9]. However, even though the electric field distribution can be calculated accurately, the insulation strengths of dielectrics cannot be calculated or predicted, that is because there are no governing equations to describe the discharge process. Hence, the insulation calculation or prediction can be viewed as the short board for virtual design and manufacturing of high-voltage electrical equipment.

The dielectrics used in electric power system include air, sulfur hexafluoride, transformer oil, electroceramics, and silicone rubber, etc. while the most commonly used is the air. Even though air discharge phenomena have been extensively studied for more than 100 years, both experimentally and theoretically, the air discharge theory is still imperfect, and therefore the air gap breakdown voltage prediction is still a great challenge and one of the most important issues to be solved for external insulation design of power transmission and transformation projects. The Townsend discharge theory [10], streamer discharge theory [11, 12], and leader discharge theory [13, 14, 15, 16] lay the foundation for the analysis and interpretation of air discharge phenomena. These classical theories were summarized from experimental measurements and observations, and they can be used to describe air discharge process. However, due to the complexity and randomness of air discharge process, it is difficult to realize air gap discharge voltage calculation with different gap configurations or under different voltage waveshapes and atmospheric environment.

Currently, discharge tests are still the main method to obtain air gap breakdown voltages. Various empirical [17, 18] and semiempirical [19, 20] formulas have been summarized for applications. However, the experimental studies cannot exhaustively reproduce all the gap configurations in actual engineering [21]. The discharge tests are usually conducted on air gaps with typical electrodes, like the rod-plane gaps, to obtain the basic characteristics of long air gap discharge, and the fitting formulas obtained by the experimental results are used to predict the breakdown voltages of engineering gaps, such as the conductor-tower window gaps, by taking the gap factor into consideration. However, the gap factors of different gap configurations should be determined by time-consuming and costly full-scale tests. When generalized to other gap structures, the fitting results are with certain deviations.

In order to find an alternative to replace the discharge tests, many scholars are devoted to the studies of air discharge mechanisms and try to achieve air gap breakdown voltage prediction by theoretical calculation. Based on air discharge theories and some simplifications and assumptions, numerous physical models have been proposed to simulate the entire duration of the discharge, including the first corona inception, streamer propagation, leader inception and propagation, and final jump [22, 23, 24, 25, 26, 27, 28, 29, 30]. However, long air gap discharge researches are faced with some challenges, mainly including the inherent complexity of the discharge process and the scientific modeling, the limitations of observation and measurement of various physical quantities involved in the discharge process, and the uncertainty of the solutions obtained by these semiempirical models [31]. Hence, there are obvious errors between the computed results of these models and the experimental measurement results. At present, the research hotspot is still concentrated on further revealing the physical mechanisms of the complex discharge process, so as to establish more scientific simulation models, but the realization of discharge prediction still needs long-term studies.

The idea of predicting air gap breakdown voltage by some artificial intelligence algorithms has attracted interests of researchers for many years, but there are only a very few related studies and applications. The artificial neural network (ANN) [32, 33], fuzzy logic system [34, 35], and support vector machine (SVM) [36, 37] have been applied to predict the discharge voltages of air insulation gaps. In [36, 37], a method based on electric field features and SVM was proposed for discharge voltage prediction of air gaps, and it has been successfully applied to predict the breakdown voltages of air gaps with typical and atypical electrodes [38, 39, 40] and the corona onset voltages of rod-plane gaps, conductors, and valve hall fittings [41, 42]. Some features were extracted from the calculation results of the electric field distribution of an air gap to characterize its spatial structure, and the SVM was applied to establish the multidimensional nonlinear relationships between these features and the air gap discharge voltage. This method offers a possible way to achieve breakdown voltage calculation of air gaps, so as to guide the external insulation design of high-voltage electrical equipment.

Under a given applied voltage waveform and specific atmospheric environment, air gap breakdown voltage is determined by its structure, which can be characterized by the electrostatic field distribution. In previous studies, the electric field features were defined in the hypothetic discharge channel between two electrodes, on the surface of typical high-voltage electrode and on the shortest interelectrode discharge path, which contains different regions classified from the perspectives of volume, area, and line. However, the discharge channel and the electrode surface are difficult to be defined for air gaps with complicated configurations. Only the shortest interelectrode path can be defined for two-electrode air gaps with arbitrary structures. It is known that there is a one-to-one correspondence between the gap structure and the electric field distribution. If the distribution characteristics of a three-dimensional spatial electric field can be mapped to those along a one-dimensional path, it is possible to define some effective features on this path to characterize the gap structure and establish their relationships with the air gap breakdown voltage by SVM, so as to be applied for breakdown voltage prediction of air gaps with arbitrary structures.

In this chapter, the above ideas are achieved preliminarily. Two sets of electric field features are defined on the shortest path of sphere-sphere gaps and rod (sphere)-plane gaps, respectively, with the U-shaped curve and monotonously decreased curve of the electric field distribution along the shortest path. A prediction model is established by SVM, and the electric field features on the shortest path are taken as the input parameters, while the output parameters are −1 and 1, respectively, characterizing whether a gap will withstand or breakdown under a given applied voltage. This model is applied to predict the breakdown voltages of sphere-sphere air gaps and large sphere-plane air gaps, and the predicted results are compared with the experimental data given in IEC 60052 and other references to verify its validity and accuracy.

## 2. Electric field features

### 2.1. Electric field distributions along the shortest path

The sphere-sphere air gaps given in IEC 60052 [43] (or IEEE Std 4 [44]) and the large sphere-plane air gaps given in [45, 46] are taken as the samples for breakdown voltage prediction studies in this chapter. The schematic diagrams of the sphere-sphere and the sphere-plane gaps are shown in Figure 1, where d is the gap distance. For sphere-sphere gaps, the sphere diameter D ranges from 5 to 200 cm, and the experimental breakdown voltages of these gaps are given in [43, 44], with different sphere diameters and gap distances. For the large sphere-plane air gap, it is composed of the sphere electrode with a mounting rod tube and the grounded plane electrode. The experimental 50% discharge voltages of these gaps, with the sphere diameter ranging from 25 to 95 cm and the gap distance ranging from 2 to 5 m, are given in [45], under the standard switching impulse voltage.

Since the sphere-sphere and the sphere-plane air gaps are with axisymmetric structures, the two-dimensional axisymmetric models are established by ANSYS, a finite element analysis software, to calculate their electric field distributions. The high-voltage sphere electrodes are applied unit voltage 1 V, and the grounded sphere electrode or the grounded plane electrode is applied zero potential. Taking a sphere-sphere gap with D = 10 cm and d = 3 cm and a sphere-plane gap with D = 45 cm and d = 3 m as examples, the cloud charts of their electric field distributions are shown in Figure 2. It can be seen that the maximum electric field strength appears at the bottom of the high-voltage sphere electrode, both for the short sphere-sphere gap and the long sphere-plane gap.

Taking the 3 cm short sphere-sphere gaps and the 3 m long sphere-plane gaps with different sphere diameters as examples, their electric field distributions along the shortest path are shown in Figure 3. For sphere-sphere gaps, the electric field distribution along the shortest path is a U-shaped curve, while the field strengths near two sphere electrodes are higher than those in the middle of the path. For sphere-plane gaps, the electric field distribution along the shortest path is a monotonously decreased curve, while the field strength gradually decreases along the path from the sphere electrode to the plane electrode. The field strength reduces quickly within 1 m from the sphere electrode, while on the path from position of 1 m to the plane electrode, the field strength changes a little. In addition, the maximum field strength decreases with the increase of the sphere diameter, both for sphere-sphere and sphere-plane gaps.

### 2.2. Definitions of electric field features

According to the electric field distribution characteristics of sphere-sphere gap and sphere-plane gaps, two sets of electric field features are, respectively, defined for these two different gap types. These features mainly include the electric field strength, electric field gradient, square of electric field strength, electric field integral (electric potential), path length, and some scaling parameters related to the above quantities used to characterize the electric field inhomogeneity.

#### 2.2.1. Sphere-sphere air gap

Some equally spaced sampling points are selected along the shortest path, and the values of their electric field strength are extracted to calculate the features. The electric field features for sphere-sphere air gap are defined as follows:

(1) Electric field strength, including the maximum values (Ehmax and Elmax), respectively, on the surface of the high-voltage and grounded sphere electrode; the minimum value Emin and the average value Ea; and the variance and standard deviation of the electric field strength along the shortest path (Estd2 and Estd):

E h max = max E i i = 1 2 m E1
E l max = max E i i = m + 1 m + 2 n E2
E min = min E i i = 1 2 n E3
E a = i = 1 n E i / n E4
E std 2 = 1 n i = 1 n E i E a 2 = 1 n i = 1 n E i 2 E a 2 E5
E std = E std 2 = 1 n i = 1 n E i 2 E a 2 E6

where n is the number of the sampling points along the shortest path, Ei is the electric field strength of the ith point, and m is the point number whose electric field strength is the minimum value.

(2) Electric field gradient, including the maximum gradient Eghm on the path from the point of Ehmax to that of Emin, the maximum gradient Eglm on the path from the point of Emin to that of Elmax, and the average gradient Ega on the whole shortest path:

E ghm = max grad E i i = 1 2 m E7
E glm = max grad E i i = m m + 1 n E8
E ga = i = 1 n grad E i / n E9

where | | is the absolute value sign.

(3) Square of electric field strength (W and Wa), respectively, means the integral of the field strength square on the shortest path and its expected value:

W = 0 d E 2 d l i = 1 n E i 2 d i E10
W a = W d 1 n i = 1 n E i 2 E11

where d is the gap distance and di is the length of each segment between two sampling points on the shortest path, di = d/(n-1).

(4) Electric field integral, that is, electric potential, including the field integral Vx on the path whose field strength exceeds Ex = x%Ehmax, the potential Vh between the high-voltage sphere electrode and the point of Emin, and the potential Vl between the point of Emin and the grounded sphere electrode:

V x = E i E x E i d l E i E x E i d i E12
V h = E min E i E h max E i d l E min E i E h max E i d i E13
V l = E min E i E l max E i d l E min E i E l max E i d i E14

(5) Path length Lx, including the length LEx of the path whose electric field strength exceeds Ex = x%Ehmax, the length Lgx of the path whose electric field gradient exceeds Egx = x%Eghm, the length LWx of the path where the sum of the electric field strength square equals to Wx = x%W, and the distance Lmin from the high-voltage sphere electrode to the point of Emin:

L x = i = 1 p d i E15

where p is the number of the points on the shortest path which meet the related conditions about the electric field strength, gradient, and square and di is the length of each segment.

(6) Electric field inhomogeneity, characterized by some scaling parameters related to the above five kinds of features. For electric field strength, these parameters include the field distortion factor Ed; the ratio of Elmax to Ehmax, namely, Erlh; and the ratio of Emin to Ehmax, namely, Erm:

E d = E h max / E a E16
E rlh = E l max / E h max E17
E rm = E min / E h max E18

For electric field square, the scaling parameter is a ratio Wrx of the sum of the electric field strength square on the path where Ei ≥ Ex = x%Ehmax to W:

W r x = E i E x E i 2 d i W , E x = x % E h max E19

For electric field integral, the scaling parameters are the ratios of Vx, Vh, and Vl to the applied voltage U, namely, Vrx, Vrh, and Vrl:

V r x = V x / U E20
V rh = V h / U E21
V rl = V l / U E22

For path length, the scaling parameters are the ratios of LEx, Lgx, LWx, and Lmin to the gap distance d, namely, LrEx, Lrgx, LrWx, and Lrmin. They all have the following expression:

L r x = L x / d E23

The x% is set as 90 and 75%; therefore, the electric field features are summarized in Table 1. There are altogether 38 features for sphere-sphere air gap.

Category Features Number
Electric field strength Ehmax, Elmax, Emin, Ea, Estd2, Estd 6
Electric field gradient Eghm, Eglm, Ega 3
Square of electric field strength W, Wa 2
Electric field integral V90, V75, Vh, Vl 4
Path length LE90, LE75; Lg90, Lg75; LW90, LW75; Lmin 7
Electric field inhomogeneity Ed, Erlh, Erm; Wr90, Wr75; Vr90, Vr75, Vrh, Vrl; LrE90, LrE75, Lrg90, Lrg75, LrW90, LrW75, Lrmin 16

### Table 1.

Electric field features for sphere-sphere air gap.

#### 2.2.2. Sphere-plane air gap

For sphere-plane air gap with a monotonously decreased curve of the electric field distribution along the shortest path, the electric field features are defined as follows:

1. The maximum, the minimum, and the average value of the electric field strength (Emax, Emin, and Ea) and the variance and standard deviation of the electric field distribution along the shortest path (Estd2 and Estd). Their calculation formulas are similar to or the same with Eqs. (1)(6).

2. The maximum and the average value of the electric field gradient (Egm and Ega), whose calculation formulas are similar to or the same with Eqs. (7)(9).

3. The square of electric field strength (W and Wa), which can be calculated, respectively, by Eqs. (10) and (11).

4. The electric field integral of the path on which the field strength exceeds x%Emax (Vx) and its ratio to the applied voltage U (Vrx). Vx and Vrx can be calculated, respectively, by Eqs. (12) and (20).

5. The length of the path on which the electric field strength exceeds x%Emax (Lx) and the ratio of Lx to the gap length d (Lrx). Lx and Lrx can be calculated, respectively, by Eqs. (15) and (23).

For long sphere-plane air gaps, x% includes 90, 75, 50, and 25%. Hence, there are altogether 25 features extracted from the shortest path used to characterize the electric field distribution of the sphere-plane air gap.

## 3. Breakdown voltage prediction method

### 3.1. Basic ideas

The proposed method for air gap breakdown voltage prediction is based on electric field features and SVM. A support vector classifier (SVC) is used to establish the prediction model. The withstanding and breakdown of an air gap under a given voltage are, respectively, denoted as −1 and 1. The input data of the prediction model are the abovementioned electric field features extracted from the electric field calculation results under different applied voltages, and the outputs are −1 and 1, respectively, means whether the gap will withstand or breakdown under the applied voltage, so as to transform the breakdown voltage prediction from a regression problem to a binary classification problem.

To be specific, if the critical breakdown voltage is Ub, the interval [(1−a)Ub, Ub] is defined as withstand voltage interval, and [Ub, (1 + a)Ub] is defined as breakdown voltage interval. The value of a is determined by experience, which is set as 0.1 in this chapter. Set the step size as 0.01Ub, and then the applied voltage values are 0.9Ub, 0.91Ub, …, Ub, 1.01Ub, …, 1.1Ub. The SVM model should be trained by some training samples to make it have the learning ability for accurate classification. By this binary classification method, one training sample can be extended to 21 samples to train the SVM model. For test samples, if the model outputs −1 to 1 under the applied voltage from Ub0-dU to Ub0, then Ub0 is the predicted breakdown voltage.

### 3.2. Brief introduction of SVM

SVM is a machine learning algorithm developed on the basis of VC dimension in statistical learning theory and the principle of structural risk minimization [47]. The fundamentals of SVM were detailed and introduced in [47, 48]. Here, a brief introduction is provided.

Set a known training sample as T = {(xi, yi)}, in which xi Rk, yi {−1, 1}, i = 1, 2, …, n. The implementation of SVM is to solve an optimization problem based on the maximum margin principle, which finds an optimal separating hyperplane to divide the sample data into two diverse classes. By application of the kernel trick, SVM transforms the original sample data into a high-dimensional Hilbert space H using nonlinear mapping. The training sample data are transformed to TΦ = {(Φ(xi), yi)}, where Φ(xi) H, yi {−1, 1}, i = 1, 2, …, n. In this feature space, the sample data can be linearly separated, and the decision function can be expressed as.

f x = sgn w T ϕ x i + b E24

where w and b are, respectively, the weight vector and bias term of the separating hyperplane (w Rk) and b is a real number. The symbol “sgn” is the signum function. When wTxi + b>0, the output is 1, while wTxi + b<0, the output is −1. Meanwhile, wTxi + b = 0 is the classification hyperplane.

The optimization problem can be expressed as

min w , b , ξ 1 2 w 2 + C i = 1 n ξ i s . t . y i w ϕ x i + b 1 ξ i , ξ i 0 , i = 1 , 2 , , n E25

where ξi is the slack variable and C is the penalty factor, which determines the balance between the maximization of the margin and the minimization of the classification error [48].

In order to make the solution of the primal problem more simple and practicable, the Lagrange function is introduced:

L w b ξ α β = 1 2 w 2 + C i = 1 n ξ i i = 1 n α i y i w x i + b 1 + ξ i i = 1 n β i ξ i E26

where α = (α1, …, αn)T and β = (β1, …, βn)T are the Lagrange multiplier vectors. By solving the partial derivatives of Eq. (26) for w, b, and ξi, according to the extremum conditions, the primal optimization problem (25) can be transformed to the following dual problems:

max α , β 1 2 i = 1 n j = 1 n α i α j y i y j x i x j + j = 1 n α j s . t . i = 1 n α i y i = 0 C α i β i = 0 , i = 1 , 2 , , l α i 0 , β i 0 , i = 1 , 2 , , l E27

By solving Eq. (27), the decision function can be obtained.

It can be seen that the function of transformation Φ is realized by inner product (Φ(xiΦ(xj)). The kernel function can be expressed as

K x i x j = ϕ x i ϕ x j E28

Hence, if the function K is selected, it is not necessary to choose the transformation Φ. K(xi, xj) is used in training and classification instead of Φ(x). The generalization performance of SVM is determined by properly selecting kernel functions. In this chapter, the radial basis function (RBF) kernel is selected as the kernel function of SVM for its good generalization performance and high computational efficiency:

K x i x j = exp γ x i x j 2 , γ > 0 E29

where γ is the kernel parameter.

The penalty factor C and the kernel parameter γ determine the classification performance of SVM. They can be optimized by grid search (GS) method or genetic algorithm (GA) based on K-fold cross validation or leave-one-out (LOO) cross validation, so as to obtain the optimal predictive model [36, 37, 38, 39].

### 3.3. Implementation procedures of the prediction method

The flow chart of the prediction method is shown in Figure 4. The implementation procedures are depicted as follows.

Firstly, the training samples with known gap structures and experimental data of breakdown voltage are used to train the SVM model. The electric field features of each training sample are extracted from the FEM calculation results of the electric field distribution. These features are normalized to [0, 1] by

x ¯ i = x i x min x max x min E30

where xi is a feature, x ¯ i is its normalized value, and xmin and xmax are, respectively, the minimum and maximum values of xi. After normalization, the electric field features are taken as the input data to train the SVM model, while the outputs are −1 and 1, respectively, correspond to the applied voltage in the withstand interval and the breakdown interval. Based on cross validation, the optimal penalty factor C and kernel parameter γ are searched by GS method or GA to obtain an optimal prediction model. Trained by the known experimental data, SVM establishes the multidimensional nonlinear relationships between the electric field features and the air gap breakdown voltage.

Then, the optimal SVM model is used to predict the breakdown voltages of test samples. For an air gap, an estimated breakdown voltage is set in the range [Umin, Umax]; the golden section search method is applied for the breakdown voltage prediction [49]. For each applied voltage, the electric field features are extracted and input to the prediction model to judge whether the output is −1 or 1. For example, the first applied voltage is Ut1 = Umax−0.618 × (UmaxUmin), and then the electric field features of this air gap are calculated under Ut1 and input to the SVM model. According to the output of the model, namely, 1 or −1, the search interval will be narrowed to [Umin, Ut1] to generate another applied voltage value, or otherwise, Ut2 = Umin + 0.618 × (UmaxUmin) is applied to calculate the electric field features which will be input to the SVM model to judge whether the output is −1 or 1. So repeatedly, the search interval is narrowed constantly, and the iterative predictions are conducted until the convergence condition UmaxUmin < ε is satisfied, where ε is the convergence precision. The breakdown voltage prediction result is the average value of the last two applied voltages.

## 4. Breakdown voltage prediction of air gaps

The proposed method is applied to predict the power frequency breakdown voltages of sphere-sphere air gaps and the switching impulse discharge voltages of sphere-plane air gaps. The predicted results are compared with the experimental data given in references.

### 4.1. Power frequency breakdown voltage prediction of sphere-sphere air gaps

#### 4.1.1. Training and test samples

The sample data of sphere-sphere air gaps are selected from IEC 60052 [43] and IEEE Std 4 [44]. The values of sphere diameter D are 5, 6.25, 10, 12.5, 15, 25, 50, 75, 100, 150, and 200 cm. The shortest gap distance is 1 cm, for those with D ranging from 5 to 25 cm, and the longest gap distance is 100 cm for that of D = 200 cm. There are altogether 271 sample data with different sphere diameters and different gap distances.

The proper selection of training samples is of vital importance for the generalization performance of the SVM model. Here, the training samples are selected according to the electric field nonuniform coefficient f, that is, the ratio of the maximum field strength Emax to the average field strength Ea = U/d, where U is the applied voltage. The maximum field strength of each sphere-sphere gap can be calculated by FEM, and therefore the electric field nonuniform coefficient f can be obtained.

According to the calculation results, all of these samples are slightly nonuniform electric field, and the values of f range from 1.00 to 1.55. Taking 0.05 as the step size, the values of f can be divided into 11 intervals. The samples belong to each f interval are collected together, and the sample sizes corresponding to each f interval are summarized in Table 2. According to the electric field nonuniform coefficient f, the 271 samples are divided into 11 groups. The training sample set is constituted by random selection of one sample from each group. Hence, there are 11 training samples, and the other 260 samples are taken as the test samples to verify the validity of the prediction method. Since the training sample selection is conducted by computer program which has a certain randomness, three different selection results are successively taken as the training sample set to train the SVM model, and the three times of prediction results are compared to validate the accuracy. The three groups of training samples are shown in Table 3.

f interval Sample size
1.00–1.05 58
1.05–1.10 57
1.10–1.15 41
1.15–1.20 23
1.20–1.25 19
1.25–1.30 18
1.30–1.35 15
1.35–1.40 13
1.40–1.45 9
1.45–1.50 7
1.50–1.55 11

### Table 2.

Sample size of each f interval.

Training sample set 1 Training sample set 2 Training sample set 3
D (cm) d (cm) D (cm) d (cm) D (cm) d (cm)
5 1.4 5 2.2 5 1.6
5 2.2 6.25 3.0 10 2.0
10 4.5 12.5 1.4 15 7.0
15 6.0 12.5 2.8 25 10
25 5.5 12.5 5.0 25 12
25 8.0 25 1.5 50 13
50 8.0 50 17 100 24
50 17 75 19 100 36
75 5.5 150 45 150 7.5
75 10 200 38 150 20
200 100 200 90 150 65

### Table 3.

Three groups of training sample set.

#### 4.1.2. Prediction results and analysis

Based on fivefold cross validation, the grid search method is applied to determine the optimal parameters of the SVM model. The search ranges of the penalty factor C and the kernel parameter γ are, respectively, set as [23, 29] and [2−8, 2−2], and the step sizes are both 20.1. Taking the training sample set 1, for example, the parameter optimization results of C and γ by GS method are shown in Figure 5. It can be seen that the best C = 90.5097, γ = 0.0167, under which the SVM model, has the highest classification accuracy for the training samples, that is, 98.2684%. The GS-optimized SVM models under the three groups of training sample set shown in Table 3 are used for breakdown voltage prediction of the 260 test samples.

For each test sample, the breakdown voltage is predicted by the golden section search method, and the prediction results are compared with the experimental data given in [43, 44]. In order to evaluate the prediction accuracy of the proposed method and the SVM model, three error indices, including the root-mean-square error (RMSE), the mean absolute percentage error (MAPE), and the mean square percentage error (MSPE), are used to examine the errors of the prediction results, which can be calculated by

RMSE = 1 n i = 1 n U t i U p i 2 E31
MAPE = 1 n i = 1 n U t i U p i U t i E32
MSPE = 1 n i = 1 n U t i U p i U t i 2 E33

where n is the number of the test samples and Ut(i) and Up(i) are, respectively, the experimental and predicted breakdown voltages of the ith test sample.

The optimal parameters and the error indices of the three times of prediction results are summarized in Table 4. It can be seen that the prediction results are with high accuracy, while the MAPEs of the three times of prediction are, respectively, 1.88, 2, and 1.4%. Taking the prediction results by training sample set 2, for example, the comparisons between the predicted and experimental breakdown voltages of sphere-sphere air gaps with different diameters are shown in Figure 6, where U is the breakdown voltage, D is the sphere diameter, d is the gap distance, and T-value and P-value, respectively, mean the test value and the prediction value of the breakdown voltage. For better comparisons, the prediction results of training samples are also plotted in Figure 6.

Results Training sample set 1 Training sample set 2 Training sample set 3
C 90.5097 78.7932 13.9288
γ 0.0167 0.0313 0.0625
RMSE 16.668 8.048 9.143
MAPE 0.0188 0.0200 0.0140
MSPE 0.0015 0.0140 0.0012

### Table 4.

Optimal parameters and error indices of the sphere-sphere air gap breakdown voltage prediction results.

It can be seen from Figure 6 that the predicted results coincide well with the experimental data, the trends of the breakdown voltages with the gap distance are the same, and the errors are within an acceptable range. The results shown in Table 4 and Figure 6 validate the feasibility and accuracy of the proposed method for sphere-sphere air gap breakdown voltage prediction.

### 4.2. Switching impulse breakdown voltage prediction of sphere-plane air gaps

#### 4.2.1. Training and test samples

The sample data of sphere-plane air gaps are selected from [45, 46]. The positive switching impulse discharge tests of sphere-plane air gaps with the sphere diameter of 25, 45, 75, and 95 cm were conducted in [45]. The applied voltage waveform is the standard 250/2500 μs switching impulse voltage. The experimental data were corrected to standard atmospheric condition. In order to make the SVM model generalize to sphere-plane gaps with different sphere diameters and gap lengths, seven test data shown in Table 5 are selected as the training samples, where D is the sphere diameter ranging from 25 to 95 cm, d is the gap length ranging from 2 to 5 m, and U50 is the 50% discharge voltage.

Training sample set
D (cm) d (m) U50 (kV) D (cm) d (m) U50 (kV)
45 2 864 25 3 1020
3 1045 75 1222
4 1186 95 1504
5 1303
Experimental data of test samples
D (cm) d (m) U50 (kV) D (cm) d (m) U50 (kV)
110 3.9 1461 150 8.5 2408
4.5 1527 10 2497
5.7 1675 200 3 2290
6.5 1796 4 2540
7.5 1973 5 2711
150 4 2101 6 2745
5.5 2242 7 2836
7 2316 8 2927

### Table 5.

Training and test samples of the sphere-plane air gaps.

The SVM model trained by the seven sample data is applied to predict the 50% discharge voltages of sphere-plane gaps with larger diameters, namely, 110, 150, and 200 cm. There are altogether 16 test samples. The predicted results will be compared with the experimental data cited from [46], as shown in Table 5. These experimental data also had been corrected to standard atmospheric condition.

#### 4.2.2. Prediction results and analysis

Based on LOO cross validation, the penalty factor C and the kernel parameter γ are optimized by the GA method. The population quantity is set as 20, the maximum generation is 200, and the crossover probability is 0.9. The search scopes of C and γ are, respectively, set as [10, 500] and [0.005, 0.25]. The fitness function is the classification accuracy of SVM for training samples. The parameter optimization results of C and γ by GA method are shown in Figure 7. It can be seen that the best C = 85.3407, γ = 0.0926, under which the SVM model has the highest classification accuracy for the training samples, that is, 97.2789%.

The GA-optimized SVM model is used to predict the U50 of large sphere-plane air gaps by the golden section search method. The initial applied voltage interval [Umin, Umax] is set as 0–4000 kV, and the convergence precision ε is set as 1 kV. The prediction results are summarized in Table 6, where U50 is the experimental data extracted from [46], Up is the predicted discharge voltage, and δ is the relative error.

D (cm) d (m) U50 (kV) Up (kV) δ (%)
110 3.9 1461 1602 9.7
4.5 1527 1647 7.9
5.7 1675 1745 4.2
6.5 1796 1807 0.6
7.5 1973 1910 −3.2
150 4 2101 2036 −3.1
5.5 2242 2132 −4.9
7 2316 2266 −2.2
8.5 2408 2413 0.2
10 2497 2560 2.5
200 3 2290 2348 2.5
4 2540 2513 −1.1
5 2711 2595 −4.3
6 2745 2700 −1.6
7 2836 2778 −2.0
8 2927 2893 −1.2

### Table 6.

50% discharge voltage prediction results of the sphere-plane air gaps.

It can be seen from Table 6 that the largest error of the prediction results is 9.7%, for the gap with D = 110 cm and d = 3.9 m. This is probably due to different experimental arrangements between [45, 46]. The U50 of sphere-plane gaps with D = 110 cm and d = 4 m in [46] is even lower than that with D = 95 cm and d = 4 m in [45]. The SVM model is trained by the experimental data cited from [45], and therefore the prediction results of sphere-plane gaps with D = 110 cm may be larger than the experimental data obtained in [46]. Overall, the prediction errors are acceptable in the view of engineering applications, while the MAPE of the 16 test samples is only 3.2%.

The U50 prediction results and the experimental data are summarized in the same graph for a better comparison, as shown in Figure 8. It can be seen that the predicted values of the discharge voltage agree well with the experimental data, with similar trends and acceptable errors. The results verify the validity and accuracy of the proposed model for discharge voltage prediction of sphere-plane air gaps, with large sphere diameter and long gap length. Within the range of certain precision, the prediction method can be used to replace the experiments, so as to reduce the testing expenses.

## 5. Conclusions

Two sets of electric field features defined on the shortest interelectrode path are, respectively, used to characterize the gap structure of the sphere-sphere air gap and the rod (sphere)-plane air gap. These features are taken as the input parameters of the SVM model, which is used to establish the breakdown voltage prediction model. The proposed method based on electric field features and SVM is applied to predict the breakdown voltages of sphere-sphere and sphere-plane air gaps. Some conclusions can be drawn as follows:

1. The proposed electric field features extracted from the shortest interelectrode path are effective to characterize the spatial structure of sphere-sphere and sphere-plane air gaps, and the multidimensional nonlinear relationships between these features and the air gap breakdown voltage can be established by SVM, so as to achieve breakdown voltage prediction of air gaps without considering the complex and random discharge process.

2. Trained by only 11 sample data selected randomly according to the electric field nonuniform coefficient f, the SVM model is able to accurately predict the power frequency breakdown voltages of IEC standard sphere-sphere air gaps. The mean absolute percentage errors of the 260 test samples, with three times of prediction by different training sample sets, are within 2%. The results validate the validity and accuracy of the proposed method for breakdown voltage prediction of sphere-sphere air gaps.

3. The proposed method is able to predict the switching impulse discharge voltages of sphere-plane air gaps, with large sphere diameters and long gap distances. The prediction results agree well with the experimental data, with similar trends and acceptable errors. The mean absolute percentage error of the 16 test samples is 3.2%, which is acceptable for engineering applications. The results verify the feasibility of the proposed model for discharge voltage prediction of large sphere-plane air gaps, which may be useful to replace the time-consuming and costly discharge tests.

The authors are still engaged in improving this model. The following work will be carried out in the future. Firstly, the electric field features will be simplified by some feature selection approaches to make it easier for applications. Secondly, the applications of this method will be extended to breakdown voltage prediction of more complex air gaps such as the practical engineering gaps. It should be noted that new problems will inevitably appear in different applications, and therefore this topic is worthy to be studied in-depth. We hope that it is possible to achieve breakdown voltage prediction of arbitrary engineering gap configurations in the future, so as to guide the insulation design of high-voltage electrical equipment by mathematical calculations rather than costly experiments.

## Acknowledgments

This work is supported by China Postdoctoral Science Foundation (2016 M602354). The authors would like to thank Dr. Shengwen Shu for his primary work on this topic during his study for PhD degree in Wuhan University. The experimental data of sphere-plane gaps are cited from [45, 46] and taken as the sample data in this paper. We would like to express thanks to those authors.

## References

1. 1. Liao CB, Ruan JJ, Liu C, Wen W, ZY D. 3-D coupled electromagnetic-fluid-thermal analysis of oil-immersed triangular wound core transformer. IEEE Transactions on Magnetics. 2014;50(11):8401904. DOI: 10.1109/TMAG.2014.2330953
2. 2. Liu XM, Yang YM, Yang F, Jadoon A. Numerical research on the losses characteristic and hot-spot temperature of laminated core joints in transformer. Applied Thermal Engineering. 2017;110:49-61. DOI: 10.1016/j.applthermaleng.2016.08.158
3. 3. Zhang XC, Li WL, Kou BQ, Cao JC, Cao HC, Gerada C, Zhang H. Electrothermal combined optimization on notch in air-cooled high-speed permanent-magnet generator. IEEE Transactions on Magnetics. 2015;51(1):8200210. DOI: 10.1109/TMAG.2014.2332437
4. 4. Wang QY, Yang X, Tian HD, Liu P, Peng ZR. A novel dissipating heat structure of converter transformer RIP bushings based on 3-D electromagnetic-fluid-thermal analysis. IEEE Transactions on Dielectrics and Electrical Insulation. 2017;24(3):1938-1946. DOI: 10.1109/TDEI.2017.006027
5. 5. Sheikholeslami M, Soleimani S, Ganji DD. Effect of electric field on hydrothermal behavior of nanofluid in a complex geometry. Journal of Molecular Liquids. 2016;213:153-161. DOI: 10.1016/j.molliq.2015.11.015
6. 6. Sheikholeslami M, Ganji DD. Impact of electric field on nanofluid forced convection heat transfer with considering variable properties. Journal of Molecular Liquids. 2017;229:566-573. DOI: 10.1016/j.molliq.2016.12.107
7. 7. Sheikholeslami M, Chamkha AJ. Electrohydrodynamic free convection heat transfer of a nanofluid in a semi annulus enclosure with a sinusoidal wall. Numerical Heat Transfer, Part A: Applications. 2016;69(7):781-793. DOI: 10.1080/10407782.2015.1090819
8. 8. Sheikholeslami M, Hayat T, Alsaedi A, Abelman S. Numerical analysis of EHD nanofluid force convective heat transfer considering electric field dependent viscosity. International Journal of Heat and Mass Transfer. 2017;108(Part B):2558–2565. DOI: 10.1016/j.ijheatmasstransfer.2016.10.099
9. 9. Sheikholeslami M, Bhatti MM. Active method for nanofluid heat transfer enhancement by means of EHD. International Journal of Heat and Mass Transfer. 2017;109:115-122. DOI: 10.1016/j.ijheatmasstransfer.2017.01.115
10. 10. Townsend JS. The Theory of Ionization of Gases by Collision. New York: Van Nostrand Company; 1910. p. 88
11. 11. Reather H. Electron Avalanches and Breakdown in Gases. London: Butterworths; 1964. p. 191
12. 12. Loeb LB, Meek JM. The Mechanism of the Electric Spark. USA: Stanford University Press; 1941. p. 188
13. 13. Les Renardières Group. Research on long air gap discharges at Les Renardières. Electra. 1972;(23):53–157
14. 14. Les Renardières Group. Research on long air gap discharges at Les Renardières–1973 results. Electra. 1974;(35):49–156
15. 15. Les Renardières Group. Positive discharges in long air gap discharges at Les Renardières–1975 results and conclusions. Electra. 1977;(53):31–153
16. 16. Les Renardières Group. Negative discharges in long air gap discharges at Les Renardières–1978 results. Electra. 1981;(74):67–216
17. 17. Gallet G, Leroy G, Lacey R, Kromer I. General expression for positive switching impulse strength valid up to extra long air gaps. IEEE Transactions on Power Apparatus and Systems. 1975;94(6):1989-1993. DOI: 10.1109/T-PAS.1975.32045
18. 18. Kishizima I, Matsumoto K, Watanabe Y. New facilities for phase-to-phase switching impulse tests and some test results. IEEE Transactions on Power Apparatus and Systems. 1984;103(6):1211-1216. DOI: 10.1109/TPAS.1984.318451
19. 19. Carrara G, Thione L. Switching surge strength of large air gaps: A physical approach. IEEE Transactions on Power Apparatus and Systems. 1976;95(2):512-514. DOI: 10.1109/T-PAS.1976.32131
20. 20. Rizk FAM. Switching impulse strength of air insulation: Leader inception criterion. IEEE Transactions on Power Delivery. 1989;4(4):2187-2195. DOI: 10.1109/61.35646
21. 21. Chen WJ, Zeng R, He HX. Research progress of long air gap discharges. High Voltage Engineering. 2013;39(6):1281-1295. DOI: 10.3969/j.issn.1003-6520.2013.06.001
22. 22. Gallimberti I. A computer model for streamer propagation. Journal of Physics D: Applied Physics. 1972;5(12):2179-2189. DOI: 10.1088/0022-3727/5/12/307
23. 23. Gallimberti I. The mechanism of the long spark formation. Journal de Physique Colloques. 1979;40(C7):193-250. DOI: 10.1051/jphyscol:19797440
24. 24. Gallimberti I, Bacchiega G, Bondiou-Clergerie A, Lalande P. Fundamental processes in long air gap discharges. Comptes Rendus Physique. 2002;3(10):1335-1359. DOI: 10.1016/S1631-0705(02)01414-7
25. 25. Fofana I, Béroual A. A model for long air gap discharge using an equivalent electrical network. IEEE Transactions on Dielectrics and Electrical Insulation. 1996;3(2):273-282. DOI: 10.1109/94.486779
26. 26. Goelian N, Lalande P, Bondiou-Clergerie A, Bacchiega GL, Gazzani A, Gallimberti I. A simplified model for the simulation of positive-spark development in long air gaps. Journal of Physics D: Applied Physics. 1997;30(17):2441-2452. DOI: 10.1088/0022-3727/30/17/010
27. 27. Becerra M, Cooray V. A self-consistent upward leader propagation model. Journal of Physics D: Applied Physics. 2006;39(16):3708-3715. DOI: 10.1088/0022-3727/39/16/028
28. 28. Arevalo L, Cooray V, Montano R. Numerical simulation of long sparks generated by positive switching impulses. Journal of Electrostatics. 2009;67(2–3):228-234. DOI: 10.1016/j.elstat.2008.12.022
29. 29. Xie YH. A Physical Model for the Simulation of Positive Discharge Development in Long Air Gaps [Thesis]. Wuhan, China: Huazhong University of Science & Technology; 2013. p. 106
30. 30. Zhou X, Zeng R, Li ZZ, Zhuang CJ. A one-dimensional thermo-hydrodynamic model for upward leader inception considering gas dynamics and heat conduction. Electric Power Systems Research. 2016;139:16-21. DOI: 10.1016/j.epsr.2015.11.028
31. 31. Zeng R, Zhuang CJ, ZQ Y, Chen S, Li ZZ, Chen WJ. Challenges and achievement in long air gap discharge research. High Voltage Engineering. 2014;40(10):2945-2955. DOI: 10.13336/j.1003-6520.hve.2014.10.001
32. 32. Mokhnache L, Boubakeur A. Prediction of the breakdown voltage in a point-barrier-plane air gap using neural networks. In: 2001 Annual Report Conference on Electrical Insulation and Dielectric Phenomena; 14–17 October 2001; Kitchener, Ontario, Canada. New York: IEEE; 2001. pp. 369-372. DOI: 10.1109/CEIDP.2001.963559
33. 33. Mokhnache L, Boubakeur A, Feliachi A. Breakdown voltage prediction in a point-barrier-plane air gap arrangement using self-organization neural networks. In: IEEE Power Engineering Society General Meeting; 6–10 June 2004; Denver, Colorado, USA. New York: IEEE; 2004. pp. 569–572. DOI: 10.1109/PES.2004.1372865
34. 34. Bourek Y, Mokhnache L, Said NN, Kattan R. Study of discharge in point-plane air interval using fuzzy logic. Journal of Electrical Engineering & Technology. 2009;4(3):410-417. DOI: 10.5370/jeet.2009.4.3.410
35. 35. Bourek Y, Mokhnache L, Said NN, Kattan R. Determination of ionization conditions characterizing the breakdown threshold of a point-plane air interval using fuzzy logic. Electric Power Systems Research. 2011;81(11):2038-2047. DOI: 10.1016/j.epsr.2011.06.012
36. 36. Shu SW. Study on Prediction for Corona Onset and Breakdown Voltages of Air Gap Based on Electric Field Features and Support Vector Machine [Thesis]. Wuhan, China: Wuhan University; 2014. p. 139
37. 37. Qiu ZB. Study on Energy Storage Features and Breakdown Voltage Prediction of Air Gap [Thesis]. Wuhan, China: Wuhan University; 2016. p. 141
38. 38. Qiu ZB, Ruan JJ, Huang DC, PU ZH, Shu SW. A prediction method for breakdown voltage of typical air gaps based on electric field features and support vector machine. IEEE Transactions on Dielectrics and Electrical Insulation. 2015;22(4):2125-2135. DOI: 10.1109/TDEI.2015.004887
39. 39. Qiu ZB, Ruan JJ, Huang DC, Wei MT, Tang LZ, Huang CP, WJ X, Shu SW. Hybrid prediction of the power frequency breakdown voltage of short air gaps based on orthogonal design and support vector machine. IEEE Transactions on Dielectrics and Electrical Insulation. 2016;23(2):795-805. DOI: 10.1109/TDEI.2015.005398
40. 40. Qiu ZB, Ruan JJ, Huang CP, Xu WJ, Tang LZ, Huang DC, Liao YF. A method for breakdown voltage prediction of short air gaps with atypical electrodes. IEEE Transactions on Dielectrics and Electrical Insulation. 2016;23(5):2685-2694. DOI: 10.1109/TDEI.2016.7736827
41. 41. Qiu ZB, Ruan JJ, Huang DC, Shu SW, Du ZY. Prediction study on positive DC corona onset voltage of rod-plane air gaps and its application to the design of valve hall fittings. IET Generation, Transmission & Distribution. 2016;10(7):1519-1526. DOI: 10.1049/iet-gtd.2015.0192
42. 42. Qiu ZB, Ruan JJ, Huang DC, Shu SW, Prediction DZY. On corona onset voltage of DC conductors and valve hall fittings. Transactions of China Electrotechnical Society. 2016;31(12):80-89. DOI: 10.3969/j.issn.1000-6753.2016.12.010
43. 43. IEC 60052. Voltage measurement by means of standard air gaps; 2012
44. 44. IEEE Std 4. IEEE standard for high-voltage testing techniques; 2013
45. 45. Wang X. The Comparison of Critical Radius of Rod-Plane Gap at Different Altitudes and Research on Altitude Correction [Thesis]. China Electric Power Research Institute: Beijing, China; 2010. p. 69
46. 46. Chen S, Zhuang CJ, Zeng R, Ding YJ, ZY S, Liao WM. Improved gap factor of large sphere-plane and its application in calculating air gap clearance in UHVDC converter station. High Voltage Engineering. 2013;39(6):1360-1366. DOI: 10.3969/j.issn.1003-6520.2013.06.011
47. 47. Vapnik VN. The Nature of Statistical Learning Theory. 2nd ed. New York: Springer-Verlag; 2000. p. 314
48. 48. Abe S. Support Vector Machines for Pattern Classification. 2nd ed. London: Springer-Verlag London Limited; 2010. p. 471
49. 49. Qiu ZB, Ruan JJ, WJ X, Huang CP. Energy storage features and a predictive model for switching impulse flashover voltages of long air gaps. IEEE Transactions on Dielectrics and Electrical Insulation. 2017;24(5):2703-2711

Written By

Zhibin Qiu and Jiangjun Ruan

Submitted: 13 June 2017 Reviewed: 06 November 2017 Published: 20 December 2017