MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi Publishing Corporation 583683 10.1155/2013/583683 583683 Research Article A Case Study of Reliability and Performance of the Electric Power Distribution Station Based on Time between Failures Baharum Adam 1 Alwan Faris Mahdi 1,2 Hasson Saad Talib 3 Zhang Yingwei 1 School of Mathematical Sciences University Sains Malaysia 11800 USM Penang Malaysia usm.my 2 Statistics Department College of Administration and Economics Baghdad University Baghdad Iraq uobaghdad.edu.iq 3 Planning and Studies Department University of Babylon Babylon Iraq uobabylon.edu.iq 2013 10 10 2013 2013 26 02 2013 20 07 2013 29 07 2013 2013 Copyright © 2013 Adam Baharum et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This paper presents an algorithm for estimating the performance of high-power station systems connected in series, parallel, and mixed series-parallel with collective factor failures caused by any part of the system equipment. Failures that occur frequently can induce a selective effect, which means that the failures generated from different equipment parts can cause failures in various subsets of the system elements. The objectives of this study are to increase the lifetime of the station and reduce sudden station failures. The case study data was collected from an electricity distribution company in Baghdad, Iraq. Data analysis was performed using the most valid distribution of the Weibull distribution with scale parameter α = 1.3137 and shape parameter β = 94.618. Our analysis revealed that the reliability value decreased by 2.82% in 30 days. The highest critical value was obtained for components T1, CBF5, CBF7, CBF8, CBF9, and CBF10 and must be changed by a new item as soon as possible. We believe that the results of this research can be used for the maintenance of power systems models and preventive maintenance models for power systems.

1. Introduction

The terms high voltage and high power mean that the voltages of the electrical energy are high enough to inflict harm or death on living things. A high-power system model is a necessary factor in the prediction capability of the system to supply a sufficient amount of electric power required at any given time. High-power systems are very complex and require maintenance of major factors. Reliability is one of the important factors in the design, operation, and maintenance of electric power systems . In 2004, Popova  studied the structure of Bayesian group replacement policies for a parallel system of n items with exponential failure times and random failure parameter. In this paper, the researcher proved the sufficiency of optimally observing the system only at failure times in the case of two items operating in parallel. In a parallel system of n items, the optimal structure of the Bayesian group replacement policies was presented by Popova when the failure times were distributed exponentially . A system subjected to external and internal failures was considered by Montoro-Cazorla and Pérez-Ocón. The failures followed a Markovian arrival process, and the operational time has a phase-type distribution . Marquez et al. presented a case study on a type M63 point machine at the Carillion Rail (formerly GTRM) Training Centre in Stafford (UK). They described the development of a condition monitoring system based on pattern recognition to detect and diagnose machine point failures. They used the data collected from a set of points under test and showed that potential failure modes can be successfully detected and diagnosed . Castro and Sanjuán presented a combined maintenance strategy in which the repair of the system failure was performed only in an interval of time within the operating period. The aim of their work was to exhibit the optimal interval in which the repairs can be performed . Kapur et al.  proposed a new software reliability growth model based on “Ito^” type of stochastic differential equation. They considered an SDE-based generalized Erlang model with logistic error detection function. Zhang et al.  provided good monitoring performance by using kernel partial least squares (KPLS). They also proposed multiblock kernel partial least squares (MBKPLS) algorithm and proved that this method performs better when it is applied to large-scale processes since it combines the advantages of both KPLS and MBPLS by being applied to monitor large-scale processes. In 2011, Barbosa et al. conducted field surveys on compression splices with suppliers and users. They also suggested some procedures and inspection routines to improve system reliability in their field study in an electricity transmission system in Brazil . The distribution and expected values for a k-out-of-n system were studied by Eryilmaz. The “k-out-of-n” is a system that has n components and operates only if at least k components are functioning . The first to define the guiding parking reliability of urban parking variable message signs (VMSs) were Mei et al. . They based the presentation of this problem by analysing the parking choice under guiding and optional parking lots and determining the guiding parking reliability of VMS by formulating a mathematical program for this reliability model. Yingwei et al.  proposed a method for the analysis of multimode batches. This method has a new scientific contribution by indicating that multimodes are separated correctly since the cross-mode correlations are considered and the common information is extracted, also that the expensive computing load is avoided since only the specific information is calculated when a mode is monitored online; after that, two different subspaces are separated, and the common and specific subspace models are built and analysed, respectively.

In the light of the above survey of related literature, the main goal of this paper is to find the optimal fitting for the failure data via time between failures (TBF). We determine the estimation of the parameters for all components for this station, as well as its reliability value. Furthermore, the current paper is the first to present this problem and its analysis; thus, the results are new and innovative.

2. Failure Functions

We assume that T is a continuous random variable with values in a positive real line. Many methods are available to specify the properties of a random variable . The best four methods specifying the properties of T are probability density function (PDF), f(t), cumulative distribution function (CDF), F(t), reliability function R(t), and hazard or failure function h(t). These functions are called “failure functions” and can be described mathematically as (1)fT(t,α,β)=βα(tα)β-1exp(-(tα)β),t>0,α,β>0,(2)F(t)=1-exp(-(tα)β),(3)R(t)=exp(-(tα)β),(4)h(t)=βα(tα)β-1,t0,α>0,β>0, where α and β are called the scale and shape parameters, respectively .

The reliability function is one of the mathematical and engineering indicators that describe the role of the equipment in the system using the probability function . Reliability is defined by the North American Reliability Council Glossary of Terms as “The degree of performance of the elements of the bulk electric system that results in electricity being delivered to customers within accepted standards, and in the amount desired” [14, 15]. Reliability can be calculated from the frequency, duration, and magnitude of the adverse effects by the equipment on the electric supply. Electric system reliability can be addressed by considering two basic and functional aspects of the electrical system: adequacy and security . Competition arises among manufacturing companies, especially within the warranty period, by increasing this period based on product quality. Reliability is an important element to measure product specifications. Furthermore, ensuring reliability reduces the cost of corrective actions on items that fail within the warranty period . From this definition of reliability, we can see that reliability analysis involves the analysis of binary outcome (0, 1) (i.e., success = 1 and failure data = 0) .

3. Problem Statement

The current paper studies and analyses the TBF of the station transformers that transform electricity from 33 KV to 11 KV, denoted as 33/11 KV. The power station under this case study included three power transformers. Each transformer had a circuit breaker with limited capacities (1,200 A), denoted as CBTi, that functions as the main circuit breaker of the transformer, where i is the number of transformers 1, 2, and 3. The transformers are connected to the communication bus between them (termed as Bus-Bar in this paper), which is linked with a group of feeders to each of the transformers. Transformers 1, 2, and 3 are separated in the Bus-Bar by circuit breakers with limited capacities of 800 A, called the Bus-Bar circuit breakers (CBB). Each feeder has a circuit breaker with a capacity of 400 A (CBFi), where i represents the number of the feeders (10 feeders). The main circuit breaker must be switched ON, and the Bus-Bar circuit breaker must be switched OFF if the transformers are operating. However, if one of these transformers fails for any reason, the CBT must be switched OFF and the CBB is switched ON to provide electricity to the failed transformer feeders. For a more detailed treatment of the station operation, the readers are referred to Alwan et al. [13, 17].

4. Methodology

The main focus of the current paper is on the performance of a station component that fails randomly; that is, the TBF is a random variable . In this case, a statistical function to identify the statistical distribution of the TBF is studied by running the distribution fitting software “EasyFit.” Data for TBFs was collected from an electricity distribution company in Baghdad, Iraq. The study period was five years. We visited the maintenance department and met the technicians and engineers. The head technician explained the nature of the work at the stations from the engineering point of view. This is shown in Table 1. The first column in Table 1 represents failures numbers, and the second column mentions the number of days required for the station to fail. For example, the first TBF value was calculated between 12 am on 1 January and 18:00 pm on 6 Jul. Accordingly, the first step down occurred after 187.833333 days. We studied the reliability of these stations and found the optimal method for their maintenance. The stations also needed further study and analysis because of the difficult conditions and the shortage in electric power in Iraq. After several days of meetings, technical notes were organized, and the workers shared their experiences in repairing these stations. That helped further our understanding of the issues at hand. The data were analysed and tested under several distributions using the “EasyFit” software, and through the optimal analysis of the data, it was found that it follows the Weibull distribution with scale parameter α=1.3137 and shape parameter β=94.618.

Time between failures of the electric power distribution station for five years.

Failure number TBFs (day)
1 187.833333
2 118.9166667
3 133.8854167
4 48.9895833
5 97.8854167
6 57.9583333
7 561
8 38.8854167
9 81.9166667
10 21.83333
11 90.8854167
12 15.79166667

A goodness of fit for this statistical distribution was tested, including the use of the Kolmogorov-Smirnov, Anderson-Darling, and Chi-square tests. The idea behind the goodness-of-fit tests is to have the “distance” critical values measured between the data and the distribution being tested. Then, that critical value is compared to some threshold value. The goodness-of-fit reports includes the test statistics and the critical values calculated for various significance levels (α=0.2, 0.1, 0.05, 0.02, 0.01). The histogram was based on sample data (TBFs). Defining the number of vertical bars was based on the total number of observations. The equation Q=1+log2N was used where N is the total number of TBFs and Q is the resulting number of classes. The height of each histogram bar indicates how many data points fall into that class. Obviously, the distribution with the lowest statistical value is the best-fitting model. The critical values can be extra useful when the null hypothesis is rejected at all selected significance levels, and it was required to know at which level it could be accepted. Figure 1 shows the quality of fits for nine distributions. The fitting result of TBFs histogram with the Weibull distribution is shown in Figure 2 while Figure 3 shows the PDF function of Weibull distribution (α=1.3137 and β=94.618). For a more detailed treatment of the fitting method, the readers are referred to Alwan et al. .

The fitting result for nine distributions of TBFs histogram.

The result of Weibull distribution fit with TBFs data histogram.

The PDF functions of TBFs data of a Weibull random variable withα=1.3137 and β=94.618.

The reliability function for a Weibull random variable is given in (3). The reliability block diagram for the system in the current case study is shown in Figure 4 which shows a visual representation of how the components are configured to form a working system. From Section 3, the CBB does not function if and only if one of the transformers does not operate. A total of 16 different components exist, excluding the CBB.

The 33/11 KV electric power distribution station mixed system.

We can classify the block diagram for the station system shown in Figure 4 as follows.

Transformers group: (5)RTG(t)=1-i=13(1-RTi(t)×RCBTi(t)),

where RT is the reliability of transformers 1, 2, and 3 and RCBT is the reliability of the main circuit breaker of transformers 1, 2, and 3.

Feeders group: (6)RFG(t)=1-i=110(1-RFi(t)),

where RF is the reliability of feeders 1 to 10.

System reliability: (7)RSYS(t)=RTG(t)×RFG(t)

5. Weibayes Analysis and Shape Parameter Value

Several analysis methods are available for a large sample of data. However, the analysis of a small or very small sample constrained the engineers/researchers to use the Weibull or Weibayes analysis . The Weibull or Weibayes analysis can describe the model’s failure that is previously unknown. Furthermore, it can be used for the goodness of fit of the behavioural data for a given distribution by analysing the shape parameter value . This process is shown in Figure 5, which can be distinguished from general forms depending on the β value, as discussed in :

0<β<1 indicates a decreasing failure rate,

β=1 indicates that the distribution approaches an exponential distribution, and the failure rate is constant;

1<β<2 indicates an increasing failure rate (IFR) in the form of a concave curve;

β=2  indicates a special case where the failure rate is linear, and the distribution takes the form of the Rayleigh distribution for electronic;

β>2 indicates IFR in the form of a convex curve;

3<β<4 indicates IFR, and the distribution approachs normal distribution.

The effect of shape parameter value, β, on the Weibull hazard rate function.

6. System Reliability

The Weibull distribution is the best analysis for the data TBF under many distributions using the “EasyFit” software. For a more detailed treatment of the fitting method, the readers are referred to Alwan et al. . The maximum likelihood method is used to estimate the scale parameter (α) and shape parameter (β) in the Weibull distribution. The maximum likelihood function of the Weibull distribution is seen as (8)α=[1ni=1n(ti)β]1/β,(9)βi+1=βi-(n/β)+i=1n[(ti)/[(1/n)i=1n(ti)β]1/β]-{(n/β)+i=1n[((ti)β/((1/n)i=1n(ti)β))ln2((ti)/[(1/n)i=1n(ti)β]1/β)]}-i=1n[(ti)β/[(1/n)i=1n(ti)β]]ln[(ti)/[(1/n)i=1n(ti)β]1/β]-{(n/β)+i=1n[((ti)β/((1/n)i=1n(ti)β))ln2((ti)/[(1/n)i=1n(ti)β]1/β)]},

where α, β, and n are the scale parameter, shape parameter, and sample of size, respectively. For a more detailed treatment of the likelihood function, the readers are referred to Guure et al. . Table 2 shows the estimated value of the scale and shape parameters of each compound.

Scale and shape parameter value estimation for each compound.

Compounds α ^ β ^
T1 454.79 2.72
T2 454.79 0.718
T3 454.79 0.718
CBT1 202.38 0.63
CBT2 202.38 0.63
CBT3 219.1 0.76
CBF1 313.17 0.787
CBF2 191.254 0.721
CBF3 326.145 1.793
CBF4 261.67 1.2
CBF5 425.309 2.88
CBF6 464.731 1.422
CBF7, CBF8, CBF9, CBF10 449.18 2.99

The reliability of the station was calculated for the times imposed from t=1 to t=30, where t is expressed in days, as shown in Table 3.

Estimated reliability values of the station for 30 days.

T (day) R ^ SYS T (day) R ^ SYS
1 0.999954003 2 0.999818754
3 0.999598937 4 0.999298779
5 0.998922077 6 0.998472283
7 0.997952567 8 0.997365866
9 0.99671491 10 0.996002256
11 0.995230305 12 0.99440132
13 0.993517438 14 0.992580685
15 0.991592984 16 0.990556163
17 0.989471964 18 0.988342051
19 0.987168011 20 0.985951362
21 0.98469356 22 0.983395999
23 0.982060016 24 0.980686895
25 0.979277872 26 0.977834133
27 0.976356821 28 0.974847035
29 0.973305837 30 0.971734246
7. Limitations of the Study and Future Work

We assumed that the restrictions to this study were as follows.

Major attention was concentrated inside the distribution station only with no focus on the source. It is to be noted that when the power distribution is in short supply (below 33 Kv), there is an increase in temperature of the transformers which leads to an abrupt step down of the electric power distribution station.

The same data has been made use of for 5 consecutive years.

EasyFit software was utilised in the study.

In this paper, we consider the electric power distribution station as an independent part. Supposing the data of failure rate (λ) for individual components of the station was obtained and dealt with using the Markov model “Hidden Semi-Markov models,” then improved performance of the power distribution station is expected. We are aware that Markov’s model relies on the failures of the present state instead of considering the past data . Our study could be used to examine preventive maintenance modelling in order to evaluate its effects on the components of the station. This could help in increasing the electricity supply and decreasing the operating cost of the power station.

8. Conclusion

The TBF has been analysed to find the best-fitting distribution. Using the distribution fitting software “EasyFit,” we found that the most valid distribution is the Weibull distribution with scale parameter α=1.3137 and shape parameter β=94.618. After estimating the scale and shape parameters of each component (16) using the maximum likelihood method, the reliability value for the system in the first day was 0.99. However, if the station still operates on the 30th day, the reliability decreased to 0.97. The reliability of the system that continues to operate decreased considerably (e.g., the reliability for 60 days is 0.91 and that for 120 days is 0.77). This result shows that parts of the station continuously deteriorate because of aging. Furthermore, using the Weibayes analysis for the β value enabled us to conclude that the station fails or experiences IFR. The components T1, CBF5, CBF7, CBF8, CBF9, and CBF10 had the highest critical value and must be changed by a new item as soon as possible.

Billinton R. Allan R. N. Reliability Evaluation of Power Systems 1996 2nd New York, NY, USA Plenum Press Popova E. Basic optimality results for Bayesian group replacement policies Operations Research Letters 2004 32 3 283 287 MR2029386 ZBL1035.62022 2-s2.0-0346280053 10.1016/j.orl.2003.07.001 Montoro-Cazorla D. Pérez-Ocón R. Reliability of a system under two types of failures using a Markovian arrival process Operations Research Letters 2006 34 5 525 530 MR2238400 ZBL1152.90408 2-s2.0-33744932688 10.1016/j.orl.2005.09.003 Marquez F. P. G. Weston P. Roberts C. Failure analysis and diagnostics for railway trackside equipment Engineering Failure Analysis 2007 14 8 1411 1426 2-s2.0-34250713480 10.1016/j.engfailanal.2007.03.005 Castro I. T. Sanjuán E. L. An optimal maintenance policy for repairable systems with delayed repairs Operations Research Letters 2008 36 5 561 564 MR2459501 ZBL1210.90055 2-s2.0-50249160540 10.1016/j.orl.2008.05.007 Kapur P. K. Anand S. Yamada S. Yadavalli V. S. S. Stochastic differential equation-based flexible software reliability growth model Mathematical Problems in Engineering 2009 2009 15 2-s2.0-68949083140 10.1155/2009/581383 581383 Zhang Y. Zhou H. Qin S. J. Chai T. Decentralized fault diagnosis of large-scale processes using multiblock kernel partial least squares IEEE Transactions on Industrial Informatics 2010 6 1 3 10 2-s2.0-76849100172 10.1109/TII.2009.2033181 Barbosa C. Francato A. Mariotoni C. Ferreira I. Mascia N. Barbosa P. Analysis of critical field procedures for power HV overhead transmission line splices installed after restructuring of Brazilian electrical sector Engineering Failure Analysis 2011 18 7 1842 1847 2-s2.0-80052715076 10.1016/j.engfailanal.2011.05.012 Erylmaz S. Dynamic behavior of k-out-of-n:G systems Operations Research Letters 2011 39 2 155 159 MR2784101 2-s2.0-79953053621 10.1016/j.orl.2011.02.008 Mei Z. Tian Y. Li D. Analysis of parking reliability guidance of urban parking variable message sign system Mathematical Problems in Engineering 2012 2012 11 MR2874591 ZBL1264.90079 2-s2.0-84858135989 10.1155/2012/128379 128379 Yingwei Z. Tianyou C. Zhiming L. Chunyu Y. Modeling and monitoring of dynamic processes IEEE Transactions on Neural Networks and Learning Systems 2012 23 2 277 284 Hamada M. S. Wilson A. G. Reese C. S. Martz H. F. Bayesian Reliability 2008 New York, NY, USA Springer xvi+436 Springer Series in Statistics 10.1007/978-0-387-77950-8 MR2724437 Alwan F. M. Baharum A. Hasson S. T. Reliability and failure analysis for high power station based on operation time Proceedings of International Conference on Statistics in Science, Business, and Engineering (ICSSBE '12) 2012 Langkawi, Malaysia 370 373 Morris P. A. Cedolin R. Feinstein C. D. Reliability of electric utility distribution systems EPRI White Paper 2000 1000424 Brunekreeft G. McDaniel T. Policy uncertainty and supply adequacy in electric power markets TILEC Discussion Paper 2005 University of Tilburg Aksezer C. S. Failure analysis and warranty modeling of used cars Engineering Failure Analysis 2011 18 6 1520 1526 2-s2.0-79960120104 10.1016/j.engfailanal.2011.05.009 Alwan F. M. Baharum A. Hasson S. T. The p.d.f fitting to time between failure for high power stations Applied Mathematical Sciences 2012 6 125-128 6327 6339 MR2987868 ZBL1264.90067 Alwan F. M. Baharum A. Hasson S. T. The performance of high power station based on time between failures (TBF) Research Journal of Applied Sciences, Engineering and Technology 2013 5 13 3489 3498 Abernethy R. B. The New Weibull Handbook: Reliability & Statistical Analysis for Predicting Life, Safety, Survivability, Risk, Cost and Warranty Claims 2003 4th Houston, Tex, USA Gulf Publishing Ebeling C. E. An Introduction to Reliability and Maintainability Engineering 1997 Singapore McGraw-Hill Guure C. B. Ibrahim N. A. Ahmed A. O. M. Bayesian estimation of two-parameter Weibull distribution using extension of Jeffreys' prior information with three loss functions Mathematical Problems in Engineering 2012 2012 13 10.1155/2012/589640 589640 MR2964996 ZBL1264.62023 Yu S.-Z. Hidden semi-Markov models Artificial Intelligence 2010 174 2 215 243 MR2724430 2-s2.0-74149089478 10.1016/j.artint.2009.11.011