The relationship between the damage rate of high voltage electrical equipment and instrumental seismic intensity

Wenchuan earthquake that occurred in China in 2008 caused severe damage to a large number of electric substations. In this paper, Kriging interpolation method was used to calculate the impact area of the instrumental seismic intensity in Wenchuan earthquake, and to compare the intensities based on strong motion observation against the instrumental seismic intensities at locations where the observation data is available. The instrumental seismic intensities were calculated for the Wenchuan Earthquake at substations in the national power grid with voltage of 110kV or higher in areas of Mianyang, Deyang, Guangyuan and Chengdu. The cumulative Gaussian distribution function was then used to fit the relationship of the curves of the damage probabilities of high-voltage electrical equipment such as transformers, voltage mutual inductors, current mutual inductors, circuit breakers, isolating switches and lightning arrester with their instrumental seismic intensities. The damage probability density distribution curve of high-voltage electrical equipment based on the instrumental seismic intensities was obtained. The results showed that: (1) In the lower seismic intensity region, the mean instrumental seismic intensity was in good agreement with the traditional seismic intensity, but there was noticeable dispersion; in regions of intensity IX and above, the instrumental intensity was lower than the seismic intensity, but there was a lower degree of dispersion. (2) Among high-voltage electrical equipment, the transformers were most vulnerable to damage and they had some damage even under lower instrumental intensity. More damage would be produced when the instrumental intensity reached VIII or above; the second most vulnerable equipment was the circuit breaker, and the damage was most likely to occur when the instrument intensity was IX or above . (3) The damage rate curves of lightning arresters, current mutual inductors, voltage mutual inductors and isolating switches were relatively close to each other and the damage probability was the highest when the instrumental intensity was about X.


Introduction
Numerous earthquakes at home and abroad have caused serious damage to high-voltage electrical equipment in the transformer substations which resulted in the failure of power grid function in the affected areas and caused great difficulties for postevent emergency rescue efforts, the lives of the affected people, and the resettlement efforts after the disaster. Research on the vulnerability of high-voltage electrical equipment in the transformer substations is of great significance for improving the seismic performance of electrical equipment, assessing the damage and functional failure of power facilities and speeding up emergency repair efforts after an earthquake. e research activities on the seismic vulnerability of high-voltage electrical equipment can be divided into three main categories: theoretical and numerical simulation, shaking table test, and statistical analysis [1]. e first two are mainly used for studying the seismic capability of the equipment, simulating seismic response and damage mechanism, and developing techniques for earthquake resistance and damage mitigation such as seismic base isolation [2][3][4]. e statistical method is used for studying the damage rates of the equipment under different ground motion intensities based on the samples of high-voltage electrical equipment damaged in earthquakes so as to obtain the vulnerability of the equipment. is method is an empirical approach that is directly linked to the actual seismic damage and is often used for seismic risk analysis, damage estimation, and economic loss assessment of the electric equipment.
In the 1990s, the Pacific Earthquake Engineering Center (PEER) and Pacific Gas and Electric Company in the US jointly established the Database of Seismic Performance of Transformer Substation Equipment for California, which recorded the damage to electrical equipment of 60 substations for different voltage classes with 220 kV and above in 12 earthquakes in California. e vulnerability curves of high-voltage electrical equipment were statistically calculated, and the results were widely used in postdisaster assessment for electrical power system [5]. e Applied Technology Council (ATC) in the US provided seismic vulnerability curves for various lifelines, which were used in the seismic risk analysis system by Federal Emergency Management Agency (FEMA) [6,7].
In recent years, a few scholars in China have studied the seismic vulnerability of substation electrical equipment with oil-immersed high-voltage transformer connected to the pipe busbar [8,9]; He et al. provided seismic vulnerability curves of transformers, busbars, and power transmission towers based on seismic damage data [10]; Liu et al. used the Weibull distribution function to fit the seismic vulnerability curves of electric porcelain electrical equipment [11]; Yang studied the relationships between the damage probability of various high-voltage electrical equipment with the peak ground acceleration and acceleration response spectrum and analyzed the functional failure modes of transformer substations under different peak ground accelerations [12].
Substantial progress has been made recently in constructing the rapid seismic intensity reporting system in various regions of China [13], and the system can quickly produce spatial distribution information on instrumental seismic intensity after an earthquake. Rapid assessment of seismic disaster based on instrumental seismic intensity meets the urgent need for postearthquake emergency response and rescue as well as engineering repair efforts. e research on the vulnerability of high-voltage electrical equipment based on instrumental seismic intensity is the basis of rapid seismic disaster assessment of power facilities.
However, the previous studies on the vulnerability of electrical equipment were based on peak acceleration, acceleration response spectrum, or traditional intensity. ere were few published reports on studying the vulnerability of substation high-voltage electrical equipment based on instrumental seismic intensity.
To study the vulnerability of substation high-voltage electric equipment, this paper first calculated the instrumental intensities at strong motion observation stations in Wenchuan earthquake.
e Kriging interpolation method was used to calculate the instrumental seismic intensities at a total of 121 110 kV and above substations for the national power grid in the worst-hit areas of Mianyang, Deyang, Guangyuan, and Chengdu in Wenchuan earthquake. Based on the seismic damage data of high-voltage electrical equipment, a cumulative Gaussian distribution function was used to express the relationship between damage rate and instrumental seismic intensity of high-voltage electrical equipment such as the transformer, circuit breaker, voltage mutual inductor, current mutual inductor, isolating switch, and lightning arrester, which was then used as the basis for seismic risk assessment and emergency response of power facilities.

Instrumental Seismic Intensity at Strong
Motion Observation Stations e instrumental seismic intensity is the intensity calculated using the strong motion observation records according to a specified method, and it directly reflects the ground motion intensity at the observation site and can be quickly obtained after an earthquake [14]. Rapid seismic intensity reporting systems have been constructed, and different algorithms for instrumental seismic intensity have been implemented in countries and regions such as the United States, Japan, and Taiwan. e national rapid seismic intensity reporting and early warning project in China is well under construction with a few rapid seismic intensity reporting networks having already been completed in a few selected regions, and an interim regulation for the calculation of instrumental seismic intensity has also been promulgated.
In accordance with the interim regulation, the calculation method for instrumental seismic intensity is defined as follows: the baseline correction and band-pass filtering were carried out first for the three-component seismic acceleration or velocity records at the observation sites, and the three-component synthesis was performed according to equations (1) and (2) to obtain the peak ground acceleration (PGA) and peak ground velocity (PGV). PGA and PGV were then used in equation (3) and (4) to calculate I PGA for the peak seismic acceleration and I PGV for the peak seismic velocity. Instrumental seismic intensity I is finally determined by equation (5). e instrumental seismic intensity is categorized into scales from I to XII.
In equation (1), a(t i ) is the acceleration record in different directions at time t i ; in equation (2), v(t i ) is the speed record in different directions at time t i .
3.17lg(PGA) + 6.59, when the PGA was synthesized from three directions , 3.20lg(PGA) + 6.59, when the PGA was synthesized from two horizontal directions , 3.23lg(PGA) + 6.82, when the PGA was synthesized from a single direction , Advances in Civil Engineering 3.00lg(PGV) + 9.77, when the PGA was synthesized from three directions, 2.96lg(PGV) + 9.78, when the PGA was synthesized from two horizontal directions, 3.11lg(PGV) + 10.21, when the PGA was synthesized from a single direction, Based on the strong motion acceleration records at 255 strong earthquake stations in Sichuan, Gansu, Ningxia, and Shaanxi during the Wenchuan earthquake and according to the definition specified in equations (1)-(5), the data obtained at these strong motion observation stations were processed, and the instrumental seismic intensities at these locations were calculated. e comparison between the seismic intensity at strong motion observation stations and the instrumental seismic intensity at the same sites is shown in Figure 1.
From Figure 1, the following can be stated. (1) e mean instrumental intensity was in good agreement with the traditional intensity in regions of intensity VIII or below, while the instrumental intensity was lower than the traditional intensity in regions of intensity IX or above. (2) e instrumental intensity had some dispersion from the traditional intensity, and the dispersion was relatively large in low-intensity regions. As the intensity increased, the dispersion decreased. ere were more data points in regions of intensity between V and VII, the maximum difference between the instrumental intensity and traditional intensity was close to 2, and the dispersion in intensity VIII regions was relatively small.

Interpolation of Instrumental Seismic
Intensity at Substations e Kriging interpolation method was used to calculate the instrumental intensity at the substation based on the instrumental seismic intensity and geo-coordinates of the strong motion observation station. e Kriging interpolation method is an interpolation method named after South African geologist P. G. Krige by French scientist Matalon [15], which is widely used in contour line interpolation in many fields. is approach is based on the regionalized variable theory and the variogram theory, and the estimated values are obtained under the premise of ensuring that the estimated values satisfy the unbiased condition and the minimum variance condition. e Kriging interpolation method considers not only the positional relationship between the points to be evaluated and the sample points but also the spatial correlation of all known points near the point to be evaluated, and thus it greatly reduces the systematic error in instrumental intensity estimation.
It is assumed that f(x) is the regionalized variable in the space where the interpolation point and the sample point are located, and it is intrinsic, and f i (i � 1,2,. . .,n) is the corresponding value at the sampling point x i (i � 1,2,. . .,n). e estimated value of f 0 is f * 0 , and f * 0 meets the following condition: e coefficient λ i (i � 1,2,. . .,n) could be calculated according to the principle of unbiasedness and minimum variance of error.
According to the principle of unbiasedness, there is us, the following can be concluded: e error variance is e partial derivative for the error variance can be calculated: Advances in Civil Engineering 3 To minimize the error variance, the extreme value of the error variance needs to be calculated, and the Lagrangian multiplier method is used. Assume We can calculate the partial derivatives, respectively, for λ i and t, which are set to be equal to zero: Equation (12) can be arranged as follows: λ i and t can be solved from equation (13). When the regionalized variable does not satisfy the second-order stationary hypothesis but satisfies the intrinsic assumption, the variogram and covariance functions should have the following relationships: Substituting equation (12) into equation (11), we should have the following: After λ i is calculated from equation (15), it can be brought into equation (6) to obtain the estimated value of f * 0 . e power supply system in Sichuan Province is composed of state grid and power grids administered by local power companies, and state grid is the backbone. e Wenchuan earthquake caused severe damage to the power grid in Sichuan Province [16]. A total of 121 substations with voltage of 110 kV and above in the state grids in Deyang, Mianyang, and Guangyuan and selected areas of Chengdu such as Dujiangyan, Pengzhou, Chongzhou, Wenjiang, and Pixian and those in Aba autonomous region managed by State Grid Company were selected as the statistical analysis samples for the study on vulnerability of high-voltage electrical equipment. According to the spatial distribution of the strong motion observation stations, the instrumental seismic intensities at the substations calculated by the above Kriging interpolation method are shown in Table 1.

Statistical Method for Damage Rate of Various High-Voltage Electrical Equipment in Substations
All high-voltage electrical instruments outside the substations such as circuit breakers, isolating switches, voltage mutual inductors, current mutual inductors, and lightning arresters belong to porcelain-column type structure. eir earthquake damage patterns were mainly cracks and oil  Advances in Civil Engineering leakage occurring in porcelain components or direct fracture of porcelain columns. e damage of the porcelain-column type structure equipment caused the aforementioned high-voltage electrical equipment to lose its normal use function. Most of the damage patterns of transformers in earthquakes were also the damage of porcelain-column type structure. e damage of the porcelain casing caused the transformer to lose its normal function. Although there was damage to oil pillow, radiator, wheel-rail fixing device, and transformer overturning at high intensity, there was no damage to the main structure of the transformer to cause its loss of normal use function.
Once the porcelain column structures were destroyed, the equipment with them could not be repaired but had to be completely replaced. At present, it is impossible to judge the damage degree of the porcelain-column-type structural equipment before the cracks, breaks, and other obvious damage phenomena. erefore, according to whether the high-voltage electrical equipment could maintain its normal use function, the damage pattern for the equipment could be divided into two stages: damaged and intact, and there was no need to further divide damage pattern into five damage stages as used for building structures. Subsequently, for porcelain column structures like circuit breakers, isolating switches, voltage mutual inductors, current mutual inductors, and lightning arresters, crack appearance itself could be defined as in the damage stage; for transformers, the damage to the porcelain casing, oil leaking, or radiator malfunction would be defined as in the damage stage.
A substation usually has 1 to 3 working transformers and many sets for other types of high-voltage electrical equipment. e damage rate R of a certain type of high-voltage electrical equipment in a substation can be expressed as follows: where n is the number of equipment damaged in the substation and N is the total number of such equipment in the substation. e range of damage rate R is [0,1]. If the damage rate obtained for a type of high-voltage electrical equipment in each substation was considered as one sample, the least square fitting could be performed for all substation samples in Table 1 to obtain the functional expression of the relationship between the damage rate and instrumental seismic intensity for all types of high-voltage electrical equipment, which became the vulnerability curves of the equipment based on instrumental seismic intensity.
ere were studies on the relationships of the damage probability of the transformers and the busbars with the peak ground motion which used a logarithmic cumulative Gaussian distribution function to fit the relationship curves [9]. In addition, there were other studies on the seismic vulnerabilities of basic components and structures of reinforced concrete which found that the use of a logarithmic cumulative Gaussian distribution function to express the relationship between their damage and the ground motion peak acceleration had a better rationality [17,18]. e instrumental seismic intensity had a linear relationship to some extent with the logarithms of the peak acceleration and the peak velocity. From the relationship between Gaussian distribution and logarithmic Gaussian distribution, it can be seen that the relationship between the damage rate of highvoltage electrical equipment and the instrumental seismic intensity could be fitted using the cumulative Gaussian distribution function.
If a random variable x follows Gaussian distribution with an expected value of μ and a standard deviation of σ, the probability density function is e cumulative Gaussian distribution function is wherein the erf(x) function is e average damage rate of various high-voltage electrical equipment in the substation and instrumental seismic intensity were fitted by formula (18) using the least squares method, and the damage rate curves of various

Results of Vulnerability Curve Fitting for Various Types of High-Voltage Electrical Equipment
Using the substations listed in Table 1 as the samples, the damage rate-instrumental seismic intensity relationship could be fitted for transformer, circuit breaker, isolating switch, current mutual inductor, voltage mutual inductor, and lightning arrester to derive parameter values μ and σ of cumulative Gaussian distribution function curves for these 6 types of high-voltage electrical equipment, as shown in Table 2. e fitting goodness index RNew was used to represent the fitting goodness of the relationship between the damage rate of high-voltage electrical equipment and the instrumental seismic intensity. e formula for RNew is where y is the actual observed value and y * is the predicted value by the fitted function. According to equation (20), the fitting goodness of the relationship curve between the damage rate of various highvoltage electrical equipment and the instrumental intensity is shown in Table 3: e curve fitting goodness of the relationship between transformer damage rate and instrumental seismic intensity is 0.809, that of circuit breaker and voltage mutual inductor is about 0.7, and the curve fitting goodness of arrester, isolating switch, and current mutual transformer is about 0.6. It can be concluded that the fitted result could basically reflect the relationship between the instrumental seismic intensity and the damage rate of various types of high-voltage electrical equipment. e fitted damage rate curves and original data samples of various high-voltage electrical equipment are shown in Figures 2(a)-2(f ). It can be seen from Figure 2 that (1) for transformers, the damage rate was close to 20% when the instrumental intensity was VII, about 40% when the instrumental intensity was VIII, above 80% when the instrumental intensity was IX, and close to 100% when the instrumental intensity was X; (2) although damage occurred in other types of equipment other than the transformers when the instrumental seismic intensity was VII, the damage rate was very low, not exceeding 5%, and the damage rate was below 15% when the instrumental intensity was VIII; when the instrumental seismic intensity was IX, the damage rate was about 40%; the damage rate of circuit breakers reached 80%, and the damage rates of transformers, isolating switches, and lightning arresters basically ranged between 45% and 60% when the instrumental seismic intensity was X; (3) the dispersion of damage rates for various types of high-voltage electrical equipment in substations was still large under different instrumental seismic intensities; for example, the damage rate for transformer in some substations had reached 100% within the instrumental intensity range of VI-VII, while no damage occurred in other substations. From the fitted curve, the damage rate of transformers was less than 30%; (4) since the maximum number of transformers in a substation did not exceed 3, the damage rate values were only concentrated in several limited fixed values under various intensities; while there were usually more sets for other types of equipment in a substation, there were more sample points shown in the Figures 2(b)-2(f ) for these types of equipment. e damage rates and probability density distribution curves for various types of equipment were compared, respectively, as shown in Figures 3 and 4.
It could be seen from the comparison diagram in Figure 3 that the transformer was the most vulnerable to earthquake than other types of high-voltage electrical equipment, and its vulnerability was significantly higher than those for other types of equipment; the circuit breaker was the second most vulnerable in other types of equipment, and the damage rate curves of lightning arresters, mutual inductors, and isolating switches were close to each other.
It could be seen from Figure 4 that the probability density of damage rate of transformers reached a peak value when the instrumental seismic intensity was VIII, and the number of damaged transformers increased rapidly. e probability density of damage rate of circuit breakers reached a peak value when the instrumental intensity was IX, and the number of damaged circuit breakers increased the most; the probability density of damage rates of isolating switch, lightning arrester, current mutual inductor, and voltage mutual inductor reached peak values when their instrumental intensities were X, and the number of damaged sets for these types of equipment increased rapidly. Figure 4 illustrates that on the one hand, different types of equipment had their respective damage-resistant strengths as the      ground motion intensity increased; on the other hand, the damage of each type of equipment around this intensity value was disperse under the influences of various incidental factors.

Conclusion
According to the strong motion acceleration records of Wenchuan earthquake, the Kriging interpolation method is used to calculate the instrumental seismic intensities at locations of a total of 121 110 kV and above substations in affected areas of Mianyang, Deyang, Guangyuan, and Chengdu in Wenchuan earthquake. e Gaussian distribution cumulative function is used to fit damage rate-instrumental seismic intensity relationship curve for transformer, circuit breaker, voltage mutual inductor, current mutual inductor, isolating switch, and lightning arrester to form vulnerability curves of high-voltage electrical equipment outside the substations based on the instrumental seismic intensity. e fitting results show that the transformers have the highest seismic vulnerability and endure certain levels of damage rate even under low instrumental intensity. e second most vulnerable equipment is the circuit breaker, followed by the lightning arrester, transformer, and isolating switch, and the seismic vulnerability curves for these types of equipment are relatively similar. e fragile link in the substations is mainly the transformer, which becomes the most important factor for the substations to maintain their functions in the earthquakes. It is necessary to speed up the research activities to improve the seismic capacity of the transformer through enhancing its shock absorption and seismic isolation performance.
At present, construction of a rapid seismic intensity reporting system is being vigorously promoted and developed in China, and a rapid reporting network has been formed in some provinces and regions, which can be used to timely disseminate the reporting information of the intensity after an earthquake. Instrumental seismic intensity will be widely used in earthquake emergency response efforts and disaster assessment. e vulnerability curve for high-voltage electrical equipment can be used for rapid assessment of seismic damage and economic loss of power equipment and can also provide reference for equipment emergency repair in power industry after an earthquake.
Wenchuan earthquake had a wide impact and heavy seismic damage. Many seismic damage samples have been obtained in regions with different levels of intensity, which provided a wealth of basic information for the study of the vulnerability of high-voltage electrical equipment. Due to sparse distribution of strong earthquake observation stations during Wenchuan earthquake, the instrumental seismic intensity values in locations of substations cannot be directly obtained. However, the instrumental seismic intensity values of substations estimated by the interpolation method must have certain errors compared with the actual values, which will bring some errors to the vulnerability statistics of highvoltage electrical equipment based on instrumental intensity. At the same time, the fitting of vulnerability curves needs to be further enriched by accumulating more seismic samples in the future, especially the samples in high-intensity regions of intensities X and XI need to be further enriched.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest regarding the publication of this paper.