<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">EPE</journal-id><journal-title-group><journal-title>Energy and Power Engineering</journal-title></journal-title-group><issn pub-type="epub">1949-243X</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/epe.2017.94B010</article-id><article-id pub-id-type="publisher-id">EPE-75224</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Engineering</subject></subj-group></article-categories><title-group><article-title>
 
 
  Optimal Economic Life Interval Evaluation Method of Transformers
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Zhen</surname><given-names>Mei</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Bo</surname><given-names>Ruan</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Yu</surname><given-names>Li</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Jun</surname><given-names>Huang</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Cheng</surname><given-names>Long</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff2"><addr-line>Research Institute of Economics and Technology, Hubei Electric Power Company Limited, Wuhan, China</addr-line></aff><aff id="aff1"><addr-line>State Key Laboratory of Advanced Electromagnetic Engineering and Technology, Huazhong University of Science &amp;amp; Tech-nology, Wuhan, China</addr-line></aff><pub-date pub-type="epub"><day>06</day><month>04</month><year>2017</year></pub-date><volume>09</volume><issue>04</issue><fpage>78</fpage><lpage>87</lpage><history><date date-type="received"><day>December</day>	<month>17,</month>	<year>2016</year></date><date date-type="rev-recd"><day>Accepted:</day>	<month>March</month>	<year>30,</year>	</date><date date-type="accepted"><day>April</day>	<month>6,</month>	<year>2017</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
   
   This paper proposed a method of optimal economic life interval evaluation of transformer, thus providing scientific support for power system planning and reconstruction. First, obtaining the transformer failure rate by using Levenberg-Marquardt based least squares algorithm on the consumption that transformer failure mode meets three-parameters Weibull distribution. Second, calculating operation cost and failure cost of transformer during its life cycle on the basis of transformer failure rate. Then, acquiring the optimal economic life interval of transformer through life cycle minimal average annual cost based on the theory of interval analysis. Finally, the proposed method is proved to be effective by applying to some real transformers. 
  
 
</p></abstract><kwd-group><kwd>Failure Rate of Transformers</kwd><kwd> Weibull Distribution</kwd><kwd> Levenberg-Marquardt Algorithm</kwd><kwd> Interval Analysis</kwd><kwd> Life Cycle Cost (LCC)</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Power enterprises are equipment-extensive, where 80% - 90% of total assets are physical equipment assets, so utilizing efficiency of main equipment directly affects overall benefit of power enterprises. Foreign experience shows equipment management based on Life cycle cost (LCC) theory can effectively improve comprehensive interests of power enterprises [<xref ref-type="bibr" rid="scirp.75224-ref1">1</xref>]. LCC means total expense of equipment over its entire life cycle, its total cost mainly depends on operating years of equipment, so research on economic life of equipment is of great significance for power enterprises to carry out scientific and rational LCC management.</p><p>So far, there are a lot of researches on LCC [<xref ref-type="bibr" rid="scirp.75224-ref2">2</xref>] [<xref ref-type="bibr" rid="scirp.75224-ref3">3</xref>]. Paper [<xref ref-type="bibr" rid="scirp.75224-ref4">4</xref>] proposes a LCC model of transformers in medium- and low-voltage distribution network. Reference [<xref ref-type="bibr" rid="scirp.75224-ref5">5</xref>] shows a 3-dimension LCC model for power system, which provides an idea for applying LCC to manage the whole network. For failure rate prediction, an optimization approach of Weibull model solution is proposed in [<xref ref-type="bibr" rid="scirp.75224-ref6">6</xref>]. An improvement of state failure rate model for power transmission and transforming equipment is presented in [<xref ref-type="bibr" rid="scirp.75224-ref7">7</xref>], which considers model parameters changing with time. Paper [<xref ref-type="bibr" rid="scirp.75224-ref8">8</xref>] proposes a method to predict economic life of power equipment which takes the maximization of annual average net income as optimization objective. But the maximization of annual average net income in the method is much subjective and the maintenance and failure cost are failed to consider the impact of failure rate. Paper [<xref ref-type="bibr" rid="scirp.75224-ref9">9</xref>] puts forward a method of economic life from economic operation of transformer in service, but the method ignores the LCC of transformers.</p><p>Based on the above, this paper puts forward an optimal economic life interval evaluation method of transformers based on fitted failure rate curve and minimal average annual cost of life cycle. First, the failure rate curve of transformers is fitted using historical data of transformer failure rate. Then, the operation and fault cost of life cycle are calculated through the failure rate curve. Finally, the optimal economic life interval of transformers is evaluated on the basis of the minimal average annual cost of life cycle, which uses the theory of interval analysis.</p></sec><sec id="s2"><title>2. Introduction to Economic Life of Transformer</title><p>Generally, transformer life can be divided into technical life, physical life and economic life and other forms [<xref ref-type="bibr" rid="scirp.75224-ref10">10</xref>].</p><sec id="s2_1"><title>2.1. Economic Life of Transformer</title><p>Physical life of transformer generally refers to the life determined by the physical wear and tear, which is also the life of the transformer from the time put into service until retired due to physical wear and tear. Physical life of transformer is mainly determined by tangible wear also related to self-quality status, the use of technology and maintenance quality. But, technical life of transformer commonly means that with the development of science and technology, existing transformer will be replaced by more advanced transformer before its physical or economic life is ended. Therefore, technical life of transformer mainly indicates the time from transformer investment to its loss of use value due to technological progress. Because of the impact by the intangible wear, technical life of transformer is generally shorter than the natural life. Based on the above, transformer physical life and technical life are defined from the use of the transformer itself.</p><p>Economic life of transformer generally refers to the life determined by operation and maintenance costs of transformer. Because of transformer worn out in its late use, operation and maintenance costs of transformer is increasing. Under that circumstances, there will have a service cycle which is economic life of transformer when the average running cost is below a certain limit. Economic life is often used to analyze and determine the best life and the best update time of the transformer [<xref ref-type="bibr" rid="scirp.75224-ref11">11</xref>]. Therefore, we can determine the best life and the best update time of the transformer when the economic life is certain. Based on the above, there will be a year that the average running cost is minimum and the year can be considered the economic life of transformer.</p></sec><sec id="s2_2"><title>2.2. Life Cycle Model of Transformer</title><p>According to the investment structure of transformer in the actual operation of power enterprises, the cost of transformer based on life cycle theory contains initial investment cost, operation cost, maintenance cost, fault cost and equipment residual value. The life cycle cost of transformer can also be divided into variable cost (including maintenance cost, fault cost) that vary with the failure rate and fixed cost (including initial investment cost, operating cost and equipment residual value) that do not vary with the failure rate.</p><p>1) Initial investment cost</p><p>The initial investment cost includes all cost before the transformer is utilized, mainly including purchase, mounting cost, test cost and so on.</p><disp-formula id="scirp.75224-formula41"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/75224x2.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x3.png" xlink:type="simple"/></inline-formula> is purchase cost, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x4.png" xlink:type="simple"/></inline-formula>is installation and mounting cost.</p><p>2) Operation cost</p><p>The operation cost mainly consists of the loss of transformer in operating.</p><disp-formula id="scirp.75224-formula42"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/75224x5.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x6.png" xlink:type="simple"/></inline-formula> is the equivalent no-load loss, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x7.png" xlink:type="simple"/></inline-formula>is the equivalent load loss, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x8.png" xlink:type="simple"/></inline-formula>is load fluctuation loss factor, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x9.png" xlink:type="simple"/></inline-formula>is the load factor which means the transformer loading over its time, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x10.png" xlink:type="simple"/></inline-formula>is run time, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x11.png" xlink:type="simple"/></inline-formula>is the purchasing cost of electricity, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x12.png" xlink:type="simple"/></inline-formula>is the loss rate of load per year,</p><disp-formula id="scirp.75224-formula43"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/75224x13.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x14.png" xlink:type="simple"/></inline-formula> is the empirical coefficient, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x15.png" xlink:type="simple"/></inline-formula>represents the number of hours in maximum load per year.</p><p>3) Maintenance cost and fault cost</p><p>Maintenance cost and fault cost vary with the failure rate. Therefore, different transformer failure time cause different economic losses. In order to quantify the economic loss of the transformer, this paper classifies transformer failures into general, severe, and catastrophic failures [<xref ref-type="bibr" rid="scirp.75224-ref12">12</xref>].</p><disp-formula id="scirp.75224-formula44"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/75224x16.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x17.png" xlink:type="simple"/></inline-formula> is maintenance cost, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x18.png" xlink:type="simple"/></inline-formula>is fault cost, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x19.png" xlink:type="simple"/></inline-formula>refers to general, severe, and catastrophic failures, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x19.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x20.png" xlink:type="simple"/></inline-formula>is type <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x19.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x21.png" xlink:type="simple"/></inline-formula> failure of transformer, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x19.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x21.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x22.png" xlink:type="simple"/></inline-formula>is the maintenance cost caused by type <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x19.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x21.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x22.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x23.png" xlink:type="simple"/></inline-formula> failure, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x19.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x21.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x22.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x23.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x24.png" xlink:type="simple"/></inline-formula>is the comprehensive loss caused by type <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x19.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x21.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x22.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x23.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x24.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x25.png" xlink:type="simple"/></inline-formula> failure, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x19.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x21.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x22.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x23.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x24.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x25.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x26.png" xlink:type="simple"/></inline-formula>is a correcting parameters.</p><p>4) equipment residual value</p><p>Equipment residual value refers to the residual value of recycling when the equipment retired.</p><disp-formula id="scirp.75224-formula45"><label>(5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/75224x27.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x28.png" xlink:type="simple"/></inline-formula> is equipment residual value, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x29.png" xlink:type="simple"/></inline-formula>is constant coefficient for Simplified representation.</p></sec></sec><sec id="s3"><title>3. Transformer Failure Rate Model</title><sec id="s3_1"><title>3.1. Three-Parameters Weibull Distribution</title><p>In general, the Weibull distribution can be used to describe the failure rate of electrical equipment with time. Equipment cumulative failure probability <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x30.png" xlink:type="simple"/></inline-formula> refers to probability of failure in the specified conditions and time, also known as unreliable function or life distribution function in the failure rate model [<xref ref-type="bibr" rid="scirp.75224-ref13">13</xref>].</p><disp-formula id="scirp.75224-formula46"><label>(6)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/75224x31.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x32.png" xlink:type="simple"/></inline-formula> is reliable function.</p><p>Therefore, assume that the service life of electrical equipment distribution in line with the three-parameters Weibull distribution of cumulative failure probability function:</p><disp-formula id="scirp.75224-formula47"><label>(7)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/75224x33.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x34.png" xlink:type="simple"/></inline-formula> is form parameter, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x35.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x36.png" xlink:type="simple"/></inline-formula>is scale parameter, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x37.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x38.png" xlink:type="simple"/></inline-formula>is location parameter, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x39.png" xlink:type="simple"/></inline-formula>and equipment will not usually fail before<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x39.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x40.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x39.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x40.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x41.png" xlink:type="simple"/></inline-formula>refers to the operating time of equipment.</p><p>The failure rate of equipment in the life cycle can be defined as the failure rate per hour after time <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x42.png" xlink:type="simple"/></inline-formula> when equipment has not failed.</p><disp-formula id="scirp.75224-formula48"><label>(8)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/75224x43.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x44.png" xlink:type="simple"/></inline-formula> is the failure rate of equipment in the life cycle.</p><p>Usually in the process of equipment operation, its failure rate goes through the following three stages: early failure rate of the characteristics is decreasing, known as practiced period; the middle of the failure rate is almost constant, known as using period; lately failure rate of the characteristics is increasing, known as wear-out period. As the <xref ref-type="fig" rid="fig1">Figure 1</xref> shows:</p><p>The curve above can be well fitted by the Weibull distribution. The failure rate function can be described as:</p><disp-formula id="scirp.75224-formula49"><label>(9)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/75224x45.png"  xlink:type="simple"/></disp-formula><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> Variation of failure rate in total life process</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/75224x46.png"/></fig><p>where three stages of the curve respectively correspond to the form parameter 0 &lt; m &lt; 1, m = 1, m &gt; 1 three situations.</p><p>The failure rate curve can be fitted piecewise through the analysis of the history statistics from the same type electrical equipment based on the Weibull distribution. Then, form parameter<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x47.png" xlink:type="simple"/></inline-formula>, scale parameter<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x48.png" xlink:type="simple"/></inline-formula>, location parameter <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x49.png" xlink:type="simple"/></inline-formula> can be estimated according to the fitted curve. Finally, the failure rate function is determined.</p></sec><sec id="s3_2"><title>3.2. Model Parameters Estimated</title><p>Frist, assuming that location parameter <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x50.png" xlink:type="simple"/></inline-formula> is known, form parameter <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x50.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x51.png" xlink:type="simple"/></inline-formula> and scale parameter <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x50.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x51.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x52.png" xlink:type="simple"/></inline-formula> can be estimated by Levenberg―Marquardt algorithm. Then, location parameter <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x50.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x51.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x52.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x53.png" xlink:type="simple"/></inline-formula> is determined by the nonlinear fitting method. Specific steps are as follows:</p><p>1) Assuming <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x54.png" xlink:type="simple"/></inline-formula> is known as<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x54.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x55.png" xlink:type="simple"/></inline-formula>, according to the failure rate function:</p><disp-formula id="scirp.75224-formula50"><label>(10)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/75224x56.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.75224-formula51"><label>(11)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/75224x57.png"  xlink:type="simple"/></disp-formula><p>2) Based on the initial value of the parameters<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x58.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x58.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x59.png" xlink:type="simple"/></inline-formula>obtained by the above formula, two data points located on the Weibull distribution are selected to construct the equations:</p><disp-formula id="scirp.75224-formula52"><label>(12)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/75224x60.png"  xlink:type="simple"/></disp-formula><p>3) According to initial value of the parameters<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x61.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x61.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x62.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x61.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x62.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x63.png" xlink:type="simple"/></inline-formula>and N sets of data, we can use Matlab for iterative calculation based on Levenberg-Mar- quardt algorithm.</p><p>4) Parameter and scale parameter can be determined form the calculation in step3.</p><p>5)<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x64.png" xlink:type="simple"/></inline-formula>, which is the cut-off point of the failure rate stages, can be calculated by the part of n sets of data based on non-linear fitting,</p><p>6) The failure rate interval can be determined by the fitted curve and error <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x65.png" xlink:type="simple"/></inline-formula> which are calculated based on Matlab:</p><disp-formula id="scirp.75224-formula53"><label>(13)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/75224x66.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x67.png" xlink:type="simple"/></inline-formula> refers to the prediction of equipment failure rate on year<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x67.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x68.png" xlink:type="simple"/></inline-formula>.</p></sec></sec><sec id="s4"><title>4. Optimal Economic Life Interval</title><p>According to the theory of LCC, initial investment cost belongs to present value but operation cost, maintenance cost, fault cost, equipment residual value do not. Therefore, operation cost, maintenance cost, fault cost and equipment residual value need to calculate their discounted value. Assuming that transformer operate to the nth year of retirement or transformation, we can calculate the total life cycle of transformer based on the life cycle model of transformer [<xref ref-type="bibr" rid="scirp.75224-ref14">14</xref>]:</p><disp-formula id="scirp.75224-formula54"><label>(14)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/75224x69.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x70.png" xlink:type="simple"/></inline-formula> is the depreciation rate, t is the operating year.</p><p>The annual average cost of LCC:</p><disp-formula id="scirp.75224-formula55"><label>(15)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/75224x71.png"  xlink:type="simple"/></disp-formula><p>As a single value of various costs tend to cause greater error, so the theory of interval analysis has been used to calculate the cost intervals. First, we can calculate the intervals of each transformers parameter. Then, all kinds of cost interval can be calculated. Finally, the optimal economic life interval of transformer can be determined when the annual average cost is the minimum.</p></sec><sec id="s5"><title>5. Case Study</title><p>This paper takes transformer SFPSZ-120000/220 in Southern China as example to evaluate its optimal economic life interval based on the method presented above.</p><sec id="s5_1"><title>5.1. Failure Rate Interval</title><p>The historical data of transformer SFPSZ-120000/220 can be found in reference [<xref ref-type="bibr" rid="scirp.75224-ref15">15</xref>]. From the <xref ref-type="fig" rid="fig2">Figure 2</xref>, the demarcation point between using period and wear-out period is about 12 years.</p><p>Frist, the location parameter of Weibull contribution can be assumed as 12, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x72.png" xlink:type="simple"/></inline-formula>, which means the transformer goes into wear-out period on 12 years. Next, the failure rate of transformer can be fitted based on Levenberg-Marquardt algorithm in order to get form parameter<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x72.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x73.png" xlink:type="simple"/></inline-formula>, scale parameter <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x72.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x73.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x74.png" xlink:type="simple"/></inline-formula>. Then, location parameter can be determined as <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x72.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x73.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x74.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x75.png" xlink:type="simple"/></inline-formula> based on nonlinear Fitting. Finally, the failure rate curve can be fitted as <xref ref-type="fig" rid="fig3">Figure 3</xref>. The failure rate function can be described as:</p><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> Distribution curve of transformers failure rate</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/75224x76.png"/></fig><fig id="fig3"  position="float"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> Fitting curve of transformers failure rate</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/75224x77.png"/></fig><disp-formula id="scirp.75224-formula56"><label>(16)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/75224x78.png"  xlink:type="simple"/></disp-formula><p>The failure rate can be interval:</p><disp-formula id="scirp.75224-formula57"><label>(17)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/75224x79.png"  xlink:type="simple"/></disp-formula><p>The upper and lower limit of failure rate is represented in <xref ref-type="fig" rid="fig3">Figure 3</xref>.</p></sec><sec id="s5_2"><title>5.2. Calculation of Optimal Economic Life Interval</title><p>The annual average cost intervals of transformer can be calculated based on the failure rate fitted above. Life cycle parameter of transformer is presented as <xref ref-type="table" rid="table1">Table 1</xref>:</p><p>According to the data in the table and the failure rate calculated above, the upper and lower limits of the average annual cost in different operating years is presented in <xref ref-type="fig" rid="fig4">Figure 4</xref>:</p><p>In <xref ref-type="fig" rid="fig4">Figure 4</xref>, the annual average cost of LCC reaches a minimum of 2519.802 thousand &#165; in 33 years when the parameters of the transformer take the upper</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Parameters of transformers cost</title></caption><table><tbody><thead><tr><th align="center" valign="middle"  colspan="3"  >Parameters of Transformers cost</th></tr></thead><tr><td align="center" valign="middle" >Initial investment cost</td><td align="center" valign="middle" >Total cost/ thousand &#165;</td><td align="center" valign="middle" >11,500</td></tr><tr><td align="center" valign="middle"  rowspan="6"  >Parameters of Operating Cost</td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x80.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >[0.6, 0.65]</td></tr><tr><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x81.png" xlink:type="simple"/></inline-formula>/h</td><td align="center" valign="middle" >[5600, 6300]</td></tr><tr><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x82.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >0.3</td></tr><tr><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x83.png" xlink:type="simple"/></inline-formula>/kw</td><td align="center" valign="middle" >96</td></tr><tr><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x84.png" xlink:type="simple"/></inline-formula>/kw</td><td align="center" valign="middle" >432</td></tr><tr><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x85.png" xlink:type="simple"/></inline-formula>/&#165;・kwh<sup>−</sup><sup>1</sup></td><td align="center" valign="middle" >[0.55, 0.65]</td></tr><tr><td align="center" valign="middle" >Parameters of equipment residual value</td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x86.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >[30%, 40%]</td></tr><tr><td align="center" valign="middle"  rowspan="5"  >Parameters of Maintenance cost and fault cost</td><td align="center" valign="middle" >Failure rate</td><td align="center" valign="middle" >(64.2%, 34.1%, 3.7%)</td></tr><tr><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x87.png" xlink:type="simple"/></inline-formula>/thousand &#165;</td><td align="center" valign="middle" >[100, 120]</td></tr><tr><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x88.png" xlink:type="simple"/></inline-formula>/thousand &#165;</td><td align="center" valign="middle" >[1210, 1280]</td></tr><tr><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x89.png" xlink:type="simple"/></inline-formula>/thousand &#165;</td><td align="center" valign="middle" >[4060, 4190]</td></tr><tr><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/75224x90.png" xlink:type="simple"/></inline-formula>/thousand &#165;</td><td align="center" valign="middle" >[8000, 9000]</td></tr><tr><td align="center" valign="middle"  rowspan="2"  >Other Parameters</td><td align="center" valign="middle" >Inflation rate</td><td align="center" valign="middle" >[3%, 3.5%]</td></tr><tr><td align="center" valign="middle" >Depreciation rate</td><td align="center" valign="middle" >[6%, 8%]</td></tr></tbody></table></table-wrap><fig id="fig4"  position="float"><label><xref ref-type="fig" rid="fig4">Figure 4</xref></label><caption><title> Annual average life cycle cost curve of transformers</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/75224x91.png"/></fig><p>limit; the annual average cost of LCC reaches a minimum of 2438.77 thousand &#165; in 31 years when the parameters of the transformer take the lower limit. Therefore, the optimal economic life interval of transformer is 31 - 33 years. When the transformer reaches the optimal economic life interval between 31 to 33 years, the power enterprise can consider the overall transformer repair or replacement.</p></sec></sec><sec id="s6"><title>6. Conclusions</title><p>This paper proposed a method of optimal economic life interval evaluation of transformer based on failure rate curve fitted and annual average cost of LCC, thus providing scientific support for the operation and retirement strategy of electricity equipment. The method includes the following features:</p><p>1) Compared with the two-parameter Weibull distribution, the three-para- meters Weibull distribution is more accurate for fitted failure rate curve of equipment. Levenberg-Marquardt algorithm is used to estimate parameters, which can make the calculation results more accurate. The failure rate can be estimated by failure rate function, thus providing scientific support for maintenance strategy of electricity equipment.</p><p>2) The annual average cost of LCC and the cost intervals is calculated based on the LCC of transformer. The optimal economic life interval of transformer is estimated on the basis of the life cycle minimal average annual cost. The uncertainty of information is considered and the calculation error of single cost is eliminated. Therefore, the method has higher rationality and credibility.</p><p>It should be specially emphasized that this paper does not consider the impact of maintenance and environmental on the failure rate, which is due to lack of actual maintenance data and operating environment conditions. If there is more access to environmental and maintenance data, we can conduct a detailed study of the impact of maintenance and environmental on the failure rate.</p></sec><sec id="s7"><title>Acknowledgements</title><p>This paper was partially supported by Science and technology project in Hubei Electric Power Co. Ltd “Under the new situation in Hubei province energy conservation, loss reduction strategy and key technology research”.</p></sec><sec id="s8"><title>Cite this paper</title><p>Mei, Z., Ruan, B., Li, Y., Huang, J. and Long, C. (2017) Optimal Economic Life Interval Evaluation Method of Transformers. Energy and Power Engineering, 9, 78-87. https://doi.org/10.4236/epe.2017.94B010</p></sec></body><back><ref-list><title>References</title><ref id="scirp.75224-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Cai,Y., Liu,L., Cheng,H., et al. (2011) Application Review of Life Cycle Cost (LCC) Technology in Power System.17, 149-154(in Chinese).</mixed-citation></ref><ref id="scirp.75224-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Li,T., Ma, W. and Huang, X. (2008) Power Transformation Equipment Management Based on Life Cycle Cost Theory. Power System Technology,11, 50-53(in Chinese).</mixed-citation></ref><ref id="scirp.75224-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Cui,X., Yin, L., Fan,C., et al. (2010) Study of LCC for power transformer in Modification of Transformer Substation.7, 69-73(in Chinese).</mixed-citation></ref><ref id="scirp.75224-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Luo,X.,Li, L., Wei,Z., et al. (2011) Applications of Life Cycle Cost Theory in Decision-Making of Investment for Distribution Transformers Renovation. Power System Technology,2, 207-211(in Chinese).</mixed-citation></ref><ref id="scirp.75224-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Liu,L., Wang, H., Cheng,H., et al. (2012) Economic Evaluation of Power Systems Based on Life Cycle Cost. Automation of Electric Power Systems,15, 45-50(in Chinese).</mixed-citation></ref><ref id="scirp.75224-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Qin,J.,Niu, Y. and Li, Z. (2012) Optimization Approach of Weibull Model Solution for Power Station Equipment Reliability. Proceedings of the CSEE,S1, 35-40(in Chinese).</mixed-citation></ref><ref id="scirp.75224-ref7"><label>7</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Wang</surname><given-names>H.</given-names></name>,<name name-style="western"><surname> Yang</surname><given-names> H.</given-names></name>,<name name-style="western"><surname> He</surname><given-names>B. et al. </given-names></name>,<etal>et al</etal>. (<year>2011</year>)<article-title>Improvement of State Failure Rate Modle for Power Transmission and Transforming Equipemnt</article-title><source> Automation of Electric Power Systems</source><volume>16</volume>,<fpage> 27</fpage>-<lpage>31(in Chinese)</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.75224-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Liu, Y.W., Ma, L., Wu, L.Y.,et al. (2012) Economic Life Model of Power Transformer and Its Application. Power System Technology, 36, 235-240(in Chinese).</mixed-citation></ref><ref id="scirp.75224-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Yu, J.L., Wang, C.F., Zhang, B.,et al. (2012) Economic Life Evaluation of Power Transformer in Service. Proceedings of The CSU-EPSA, 22, 86-90(in Chinese).</mixed-citation></ref><ref id="scirp.75224-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Bernstein B. S. andBrancato, E. L. (1993) Aging of Equipment in the Electric Utilities. IEEE Transactions on Electrical Insulation, 5, 866-875. 
https://doi.org/10.1109/14.237747</mixed-citation></ref><ref id="scirp.75224-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Li, W., Yu,Q., andLuo, R.(2012) Research on the Algorithm of Economic Life Forecast for Power Transformer Based on Risk Analysis.1, 74-77(in Chinese).</mixed-citation></ref><ref id="scirp.75224-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Zhang, Y., Liao,R., Yang, L., et al.(2014) Maintenance Decision Making of Power Transformer by Integrated Considering Reliability and Economy.11, 191-200(in Chinese).</mixed-citation></ref><ref id="scirp.75224-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Xie,N., Cheng,M., Xiao,D.,et al. (2014) Economic Disposal Time of Primary Electricity Equipments. 6, 165-168(in Chinese).</mixed-citation></ref><ref id="scirp.75224-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Jiang, X., Wu,W. and Chen, X. (2011) Transformer Life Cycle Cost Model Based on Interval Analysis.9, 50-53(in Chinese).</mixed-citation></ref><ref id="scirp.75224-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Sun, P., Chen,S. and Zhang, C. (2012) Assessment of Failure Rate for Substation Equipment Life Cycle Based on Marquardt Parameter Estimation Method.1, 85-90(in Chinese).</mixed-citation></ref></ref-list></back></article>