CN103675621B - A kind of bushing shell for transformer insulation life prediction unit and method - Google Patents

A kind of bushing shell for transformer insulation life prediction unit and method Download PDF

Info

Publication number
CN103675621B
CN103675621B CN201310635349.8A CN201310635349A CN103675621B CN 103675621 B CN103675621 B CN 103675621B CN 201310635349 A CN201310635349 A CN 201310635349A CN 103675621 B CN103675621 B CN 103675621B
Authority
CN
China
Prior art keywords
transformer
formula
bushing shell
point
bushing
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201310635349.8A
Other languages
Chinese (zh)
Other versions
CN103675621A (en
Inventor
滕云
李勇
回茜
霍腾宇
李岩
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenyang University of Technology
Original Assignee
Shenyang University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shenyang University of Technology filed Critical Shenyang University of Technology
Priority to CN201310635349.8A priority Critical patent/CN103675621B/en
Publication of CN103675621A publication Critical patent/CN103675621A/en
Application granted granted Critical
Publication of CN103675621B publication Critical patent/CN103675621B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

A kind of bushing shell for transformer insulation life prediction unit and method, belong to bushing shell for transformer insulation technology field.Gather bushing shell for transformer data sample, including the temperature residing for bushing shell for transformer, humidity, flow through the electric current of bushing shell for transformer, the pollution degree of bushing shell for transformer and bushing shell for transformer insulating coefficient;Non-linear mere body method NSM is adopted to predict bushing shell for transformer insulation life;Bushing shell for transformer insulation life value according to output, overhauls bushing for transformer or changes.The present invention proposes the method that data are processed by non-linear mere body method, solve the uneven of Data distribution information, make information distribution more reasonable, bushing shell for transformer insulation component attributes information and bushing shell for transformer insulation essential factors space information can be used fully, Accurate Prediction to bushing shell for transformer insulation life, maintenance or replacing for bushing for transformer provide foundation.

Description

A kind of bushing shell for transformer insulation life prediction unit and method
Technical field
The invention belongs to bushing shell for transformer insulation technology field, particularly to a kind of bushing shell for transformer insulation life prediction unit and method.
Background technology
Transformator is the power equipment undertaking voltage transformation, electric energy transmission and terminal distribution in power system.If transformer high-low-voltage side sleeve pipe existing defects or fault, will directly jeopardize safe operation and the power supply reliability thereof of transformator.For low pressure reversed stream pipe, due to its simple in construction, it is necessary to the repair time is short, it is easily serviced.And for high pressure reversed stream pipe, owing to its structure is complicated, maintenance period is long, maintenance difficulty is bigger.Therefore, carrying out in time sleeve pipe overhauling and safeguard, it is highly important for reducing and stopping sleeve pipe accident.And how to accomplish to prevent trouble before it happens, the running status of paying attention bushing shell for transformer, implement state estimation and carry out the key that effective fault diagnosis is repair based on condition of component success or not.
Summary of the invention
For the deficiencies in the prior art, it is an object of the invention to provide a kind of bushing shell for transformer insulation life prediction unit and method, by setting up the bushing shell for transformer insulation life prediction and warning model of applicable practical implementation, solve the accurate of the damaged condition to set pipe insulation to predict, it is simple to the maintenance of maintainer.
The technical scheme is that and be achieved in that: a kind of bushing shell for transformer insulation life prediction unit, including: including sensor, filthy monitor, insulation tester, filter circuit, data acquisition chip, central processing unit and computer, described sensor farther includes current transformer and warm and humid sensor;
Current transformer: the electric current on measuring transformer sleeve pipe;
Warm and humid sensor: for measuring transformer sleeve pipe temperature and humidity in the environment;
Filthy monitor: for measuring the pollution degree on sleeve pipe;
Insulation tester: for measuring the insulating coefficient of sleeve pipe;
Filter circuit: filter the high communication number of current transformer, warm and humid sensor, filthy monitor and insulation tester output;
Data acquisition chip: for the signal of filter circuit output is AD converted;
Central processing unit: for the signal provided according to data acquisition chip, adopts non-linear mere body method NSM to predict bushing shell for transformer insulation life;
Computer: receive the output signal of central processing unit, carry out overhauling or replace offer Back ground Information for bushing shell for transformer.
Described current transformer, warm and humid sensor, filthy monitor, insulation tester filtered circuit respectively are connected to data acquisition chip, and data acquisition chip connects computer after central processing unit.
Filter circuit includes amplifier, resistance and electric capacity, and one end of the first resistance connects sensor, and the other end connects the inverting input of amplifier and one end of the first electric capacity, the other end ground connection of the first electric capacity;One end ground connection of the second resistance, the other end connects the inverting input of amplifier through the second electric capacity, is also associated with one end of the 3rd resistance at the center of the second resistance and the second electric capacity, and the other end of the 3rd resistance connects the outfan of amplifier;The outfan of amplifier is connected to the input of data acquisition chip through the 4th resistance.
The method adopting bushing shell for transformer insulation life prediction unit that bushing shell for transformer insulation life is predicted, comprises the following steps:
Step 1: gather bushing shell for transformer data sample, including the temperature x residing for bushing shell for transformert1, humidity xt2, flow through the electric current x of bushing shell for transformert3, bushing shell for transformer pollution degree xt4With bushing shell for transformer insulating coefficient xt5
Step 2: adopting non-linear mere body method NSM to predict bushing shell for transformer insulation life, detailed process is as follows:
Step 2.1: data step 1 collected are expressed as: P1=(xt1, xt2, xt3, xt4, xt5);Build 5 dimension phase spaces, by coordinate P1Tie up in phase space 5 and represent, including P1 from being expressed as 6 points, i.e. P in interior1, P2, P3, P4, P5, P6, the coordinate of each point is built into the row of following matrix:
x t 1 , x t 2 , x t 3 , x t 4 , x t 5 x t 1 , 2 x t 2 , x t 3 , x t 4 , x t 5 x t 1 , x t 2 , 3 x t 3 , x t 4 , x t 5 x t 1 , x t 2 , x t 3 , 4 x t 4 , x t 5 x t 1 , x t 2 , x t 3 , x t 4 , 5 x t 5 x t 1 , 2 x t 2 , 3 x t 3 , x t 4 , x t 5
Step 2.2: calculate the target function value of each point, it is determined that object function maximum of points and object function minimum point;
Object function is:
y = Σ j = 1 m ( 1 + e - x tj ) - - - ( 1 )
In formula, y is target function value, j=1,2 ..., m, wherein, m is that positive integer wherein represents the number gathering variable, t express time;
By P1, P2, P3, P4, P5, P6The coordinate of 6 points substitutes in formula (1) respectively, calculates the target function value y obtaining at 61, y2... yn, yn+1, determine maximum through comparing, formula represents:
yh=max (yi)
In formula, yhRepresent target function value maximum;H represents yhSubscript, be used for identify maximum;I=1,2 ..., 6;
Determine that minima equation below represents through comparing:
yl=min (yi)
In formula, ylFor object function minima, l is ylSubscript, be used for identify minima;
The point P with object function maximum is determined according to above-mentioned minima and maximumhWith the point with object function minima;
Step 2.3: determining the center of gravity of point except object function maximum of points, formula is as follows:
P ‾ = Σ i ≠ h P i / n - - - ( 2 )
In formula,Center of gravity for point;
Step 2.4: object function maximum of points is carried out reflection operation, makes maximum of points move to center of gravity opposite direction, produces new some P*, formula is:
P * = ( 1 + α ) P ‾ - α P h - - - ( 3 )
In formula, P*For the new point produced;α is reflection coefficient, and has 0 < α < 1, α more big, reflects more remote, otherwise then more near, for making primary be distributed widely in feas ible space, takes α=0.9;
Step 2.5: calculate the target function value y of newly generated point*If, yl≤y*≤yh, then P is used*Replace Ph
If the target function value y of newly generated point*< yl, then need newly generated point is carried out Extension Operations, perform step 2.6;
If the target function value y of newly generated point*> yh, then need to carry out newly generated point shrinking computing, repeated execution of steps 2.7;
Step 2.6: newly generated point is carried out Extension Operations, and formula is as follows:
P * * = ( 1 - &gamma; ) P &OverBar; + &gamma; P * - - - ( 4 )
In formula, P**For another the new point produced, γ is flare factor, and has γ > 1;
Calculate P**Target function value y**If, P**Target function value y**> yl, then P is used**Replace Ph, perform step 2.8:
If y**< yl, then expand failure, then use P*Replace Ph;Flare factor can not be excessive, takes γ=1.2, performs step 2.8:
Step 2.7: carrying out newly generated point shrinking computing, formula is as follows:
P * * = ( 1 - &beta; ) P &OverBar; + &beta; P * - - - ( 5 )
In formula, β is constriction coefficient, and has 0 < β < 1, for making initial vertax distribution more uniform, namely gets contracting factor beta=0.5;
If P**Target function value y**< yh, then P is used**Replace Ph, perform step 2.8;
If y**> yh, it was shown that shrink unsuccessfully, from newly determined P1, P2, P3, P4, P5, P6The coordinate of 6 points, formula is as follows:
Pi=(Pi+Pl)/2(6)
Step 2.8: redefine newly generated P1, P2, P3, P4, P5, P6The center of gravity of 6 pointsFormula is as follows:
P &OverBar; &prime; = &Sigma; i &NotEqual; h P i / n
Step 2.9: the center-of-gravity value that the newly generated center-of-gravity value of step 2.8 and step 2.3 calculate compares, and formula is:
| P &prime; &OverBar; - P &OverBar; | &le; &epsiv; - - - ( 7 )
In formula, ε is for representing threshold value, if formula (7) is set up, then performs step 2.10, otherwise usesReplaceRepeat step 2.4~step 2.9;
Step 2.10: determineAnd PhThe midpoint x of linet
Step 2.11: set up object function, the output indication transformer sleeve pipe insulation life of object function;
Object function is:
Y=f (xt)+k(xt)(8)
In formula, y indication transformer sleeve pipe insulation life, f (xt) for evaluating item, k (xt) for the bound term of object function;
Step 3: the bushing shell for transformer insulation life value according to output, overhauls bushing for transformer or changes.
Beneficial effects of the present invention: the present invention proposes the method that data are processed by non-linear mere body method, solve the uneven of Data distribution information, make information distribution more reasonable, bushing shell for transformer insulation component attributes information and bushing shell for transformer insulation essential factors space information can be used fully, Accurate Prediction to bushing shell for transformer insulation life, maintenance or replacing for bushing for transformer provide foundation.
Accompanying drawing explanation
Fig. 1 embodiment of the present invention bushing shell for transformer insulation life prediction unit structured flowchart;
Fig. 2 embodiment of the present invention bushing shell for transformer insulation life prediction unit filter circuit schematic diagram;
The data acquisition chip of Fig. 3 embodiment of the present invention and CPU circuit catenation principle figure;
The bushing shell for transformer insulation life Forecasting Methodology flow chart of Fig. 4 embodiment of the present invention;
Fig. 5 embodiment of the present invention prediction insulation life curve and actual insulation life curve schematic diagram.
Detailed description of the invention
Below in conjunction with accompanying drawing, embodiments of the present invention are described in further detail.
Present embodiment is to bushing shell for transformer insulation life prediction unit, and its structure is as it is shown in figure 1, specifically include sensor, filter circuit, data acquisition chip, trigger, central processing unit and computer.
The sensor adopted in present embodiment includes: current transformer and warm and humid sensor.Wherein, the model of current transformer is LZJC-10Q1000/5, the electric current on measuring transformer sleeve pipe,
The model of warm and humid sensor is PCMini70, for measuring transformer sleeve pipe temperature and humidity in the environment.
Filter circuit is as shown in Figure 2.Filter circuit is with the amplifier LM305 one group of signal amplification circuit being core, except amplifier LM305, further comprises resistance and electric capacity.One end of resistance R1 connects sensor, and the other end connects the inverting input of amplifier LM305 and one end of electric capacity C1, the other end ground connection of electric capacity C1.One end ground connection of resistance R2, the other end connects the inverting input of amplifier LM305 through electric capacity C2, is also associated with one end of resistance R3 at the center of resistance R2 Yu electric capacity C2, and the other end of R3 connects the outfan of amplifier LM305.The outfan of amplifier LM305 is connected to the input of data acquisition chip through resistance R4.8 feet of amplifier LM305 connect 5V working power.
In present embodiment current transformer, warm and humid sensor, filthy monitor outfan be respectively connected with filter circuit, filter circuit is the electric capacity C1 by accessing a 0.005uF between the inverting input and 4 feet of amplifier LM305, effective elimination high communication number, increases the accuracy of signal.
The model of the filthy monitor that present embodiment adopts is SC-YHM, for measuring the pollution degree on sleeve pipe.
The model of the insulation tester that present embodiment adopts is DL09-SDM50, for measuring the insulating coefficient of sleeve pipe.
The data acquisition core pattern number that present embodiment adopts be ADC0809, is 8 passage 16 bit data acquisition chips, for the data that obtain of sampling are carried out analog digital conversion, ± 5V powers.
The central processing unit of present embodiment adopts chip microcontroller, and the model of single-chip microcomputer is 8051 series of digital signals processors, and its clock frequency is 6MHz, and interface power is 5V.
Owing to ADC0809 chip internal does not have clock source. present embodiment utilizes the address latch that single-chip microcomputer 8051 provides to control input signal ALK after d type flip flop two divided-frequency. and the clock as ADC0809 inputs.The ALK end of single-chip microcomputer 8051 exports positive pulse regularly with the oscillation frequency clock of 1/6.Because the clock frequency of ADC0809 has certain scope 10-1280kHz. therefore the present embodiment clock frequency that takes single-chip microcomputer 8051 is 6MHz, the frequency arranging the output of ALK end is 1MHz, then through being 500kHz after two divided-frequency.Meet the ADC0809 requirement to clock.
In present embodiment, data acquisition chip is with CPU circuit catenation principle figure as shown in Figure 3.The connection of each parts of this device is as follows: IN-0 to the IN-3 end of the outfan of filthy monitor, Temperature Humidity Sensor, current sensor and insulation tester filtered circuit connection data acquisition chip ADC0809 respectively, the USB interface of computer converts RS232 to through PL2303, again through TXD, RXD pin, single-chip microcomputer is driven.PL2303 is that serial ports turns USB control chip, for the output signal conversion of computer is become the signal that single-chip microcomputer may identify which.The temperature x of bushing shell for transformert1, humidity xt2, flow through the electric current x of bushing shell for transformert3, bushing shell for transformer pollution degree xt4With bushing shell for transformer insulating coefficient xt5Send into by sensor that data acquisition chip ADC0809 carries out synchronized sampling, maintenance, A/D convert digital signal to, send into calculating that single-chip microcomputer carries out classifying and data process, and be transported to a distant place by wireless communication module, in order to staff's maintenance.
The present embodiment above-mentioned bushing shell for transformer insulation life prediction unit of employing predicts the method in the life-span of insulated sleeve, and its flow process is as shown in Figure 4.Comprise the following steps:
Step 1: gather bushing shell for transformer data sample, as shown in table 1: include the temperature x residing for bushing shell for transformert1=22 DEG C, humidity xt2=45%rh, flow through the electric current x of bushing shell for transformert3=0.2A, bushing shell for transformer pollution degree xt4=0.2mg/cm2With bushing shell for transformer insulating coefficient xt5=3.03, it is expressed in matrix as: P1=(22,45,0.2,0.2,3.03);
Table 1 is bushing shell for transformer data sample
Collecting sample Sampled value
Temperature 22℃
Humidity 45%rh
Electric current 0.2A
Pollution degree 0.2mg/cm2
Insulating coefficient 3.03
Step 2: adopting non-linear mere body method NSM to predict bushing shell for transformer insulation life, detailed process is as follows:
Step 2.1: build 5 dimension phase spaces, by coordinate P1Tie up in phase space 5 and represent, including P1 from being expressed as 6 points, i.e. P in interior1, P2, P3, P4, P5, P6, the coordinate of each point is built into the row of following matrix:
x t 1 , x t 2 , x t 3 , x t 4 , x t 5 x t 1 , 2 x t 2 , x t 3 , x t 4 , x t 5 x t 1 , x t 2 , 3 x t 3 , x t 4 , x t 5 x t 1 , x t 2 , x t 3 , 4 x t 4 , x t 5 x t 1 , x t 2 , x t 3 , x t 4 , 5 x t 5 x t 1 , 2 x t 2 , 3 x t 3 , x t 4 , x t 5
Step 2.2: calculate the target function value of each point, it is determined that object function maximum of points and object function minimum point;
Object function is:
y = &Sigma; j = 1 m ( 1 + e - x tj ) - - - ( 1 )
In formula, y is target function value, j=1,2 ..., m, wherein, m is that positive integer wherein represents the number gathering variable, m=5, t express time in present embodiment.
By P1, P2, P3, P4, P5, P6The coordinate of 6 points substitutes in formula (1) respectively, calculates the target function value y obtaining at 61, y2... yn, yn+1, determine maximum through comparing, formula represents:
yh=max (yi)
In formula, yhRepresent target function value maximum;H represents yhSubscript, be used for identify maximum;I=1,2 ..., 6;
Determine that minima equation below represents through comparing:
yl=min (yi)
In formula, ylFor object function minima, l is ylSubscript, be used for identify minima;
The point P with object function maximum is determined according to above-mentioned minima and maximumhWith the point with object function minima;
Step 2.3: determining the center of gravity of point except object function maximum of points, formula is as follows:
P &OverBar; = &Sigma; i &NotEqual; h P i / n - - - ( 2 )
In formula,Center of gravity for point;
Step 2.4: object function maximum of points is carried out reflection operation, makes maximum of points move to center of gravity opposite direction, produces new some P*, formula is:
P * = ( 1 + &alpha; ) P &OverBar; - &alpha; P h - - - ( 3 )
In formula, P*For the new point produced;α is reflection coefficient, for making primary be distributed widely in feas ible space, and α=0.9 in present embodiment;
Step 2.5: calculate the target function value y of newly generated point*If, yl≤y*≤yh, then P is used*Replace Ph
If the target function value y of newly generated point*< yl, then need newly generated point is carried out Extension Operations, perform step 2.6;
If the target function value y of newly generated point*> yh, then need to carry out newly generated point shrinking computing, repeated execution of steps 2.7;
Step 2.6: newly generated point is carried out Extension Operations, and formula is as follows:
P * * = ( 1 - &gamma; ) P &OverBar; + &gamma; P * - - - ( 4 )
In formula, P**For another the new point produced, γ is flare factor, and has γ > 1;
Calculate P**Target function value y**If, P**Target function value y**> yl, then P is used**Replace Ph, perform step 2.8;
If y**< yl, then expand failure, then use P*Replace Ph;Flare factor can not be excessive, takes γ=1.2, performs step 2.8;
Step 2.7: carrying out newly generated point shrinking computing, formula is as follows:
P * * = ( 1 - &beta; ) P &OverBar; + &beta; P * - - - ( 5 )
In formula, β is constriction coefficient, and has 0 < β < 1, and for making initial vertax distribution more uniform, present embodiment gets contracting factor beta=0.5;
If P**Target function value y**< yh, then P is used**Replace Ph, perform step 2.8;
If y**> yh, it was shown that shrink unsuccessfully, from newly determined P1, P2, P3, P4, P5, P6The coordinate of 6 points, formula is as follows:
Pi=(Pi+Pl)/2(6)
Step 2.8: redefine newly generated P1, P2, P3, P4, P5, P6The center of gravity of 6 pointsFormula is as follows:
P &OverBar; &prime; = &Sigma; i &NotEqual; h P i / n
Step 2.9: the center-of-gravity value that the newly generated center-of-gravity value of step 2.8 and step 2.3 calculate compares, and formula is:
| P &prime; &OverBar; - P &OverBar; | &le; &epsiv; - - - ( 7 )
In formula, ε=0.01, if formula (6) is set up, then perform step 2.10, otherwise useReplaceRepeat step 2.4~step 2.9;
Step 2.10: determineAnd PhThe midpoint x of linet
Step 2.11: set up object function, the output indication transformer sleeve pipe insulation life of object function;
Object function is:
Y=f (xt)+k(xt)(8)
In formula, y indication transformer sleeve pipe insulation life, f (xt) for evaluating item, in present embodimentRepresenting constraint function, in formula, T represents the specified time limit of sleeve pipe, the time that t ' uses for sleeve pipe;
Step 3: the bushing shell for transformer insulation life value according to output, testing staff and bushing shell for transformer insulation life value physical life compare, and carry out the maintenance of necessity or are replaced.
Certain is predicted that certain bushing shell for transformer insulation life is predicted by present embodiment, and result is as shown in Figure 5, it is known that prediction insulation life index and actual insulation life index curve, it was predicted that error is within ± 8%.
Although the foregoing describing the specific embodiment of the present invention, but the those skilled in the art in this area should be appreciated that these are merely illustrative of, it is possible to these embodiments are made various changes or modifications, without departing from principles of the invention and essence.The scope of the present invention is only limited by the claims that follow.

Claims (1)

1. a bushing shell for transformer insulation life Forecasting Methodology, adopts bushing shell for transformer insulation life prediction unit to realize, this device
Including: sensor, filthy monitor, insulation tester, filter circuit, data acquisition chip, central processing unit and computer, described sensor farther includes current transformer and warm and humid sensor;
Current transformer: the electric current on measuring transformer sleeve pipe;
Warm and humid sensor: for measuring transformer sleeve pipe temperature and humidity in the environment;
Filthy monitor: for measuring the pollution degree on sleeve pipe;
Insulation tester: for measuring the insulating coefficient of sleeve pipe;
Filter circuit: filter the high communication number of current transformer, warm and humid sensor, filthy monitor and insulation tester output;
Data acquisition chip: for the signal of filter circuit output is AD converted;
Central processing unit: for the signal provided according to data acquisition chip, adopts non-linear mere body method NSM to predict bushing shell for transformer insulation life;
Computer: receive the output signal of central processing unit, carry out overhauling or replace offer Back ground Information for bushing shell for transformer;
Described current transformer, warm and humid sensor, filthy monitor, insulation tester filtered circuit respectively are connected to data acquisition chip, and data acquisition chip connects computer after central processing unit;
Described filter circuit includes amplifier, resistance and electric capacity, and one end of the first resistance connects sensor, and the other end connects the inverting input of amplifier and one end of the first electric capacity, the other end ground connection of the first electric capacity;One end ground connection of the second resistance, the other end connects the inverting input of amplifier through the second electric capacity, is also associated with one end of the 3rd resistance at the center of the second resistance and the second electric capacity, and the other end of the 3rd resistance connects the outfan of amplifier;The outfan of amplifier is connected to the input of data acquisition chip through the 4th resistance;
It is characterized in that: comprise the following steps:
Step 1: gather bushing shell for transformer data sample, including the temperature x residing for bushing shell for transformert1, humidity xt2, flow through the electric current x of bushing shell for transformert3, bushing shell for transformer pollution degree xt4With bushing shell for transformer insulating coefficient xt5
Step 2: adopting non-linear mere body method NSM to predict bushing shell for transformer insulation life, detailed process is as follows:
Step 2.1: data step 1 collected are expressed as: P1=(xt1,xt2,xt3,xt4,xt5);Build 5 dimension phase spaces, by coordinate P1Tie up in phase space 5 and represent, including P1 from being expressed as 6 points, i.e. P in interior1,P2,P3,P4,P5,P6, the coordinate of each point is built into the row of following matrix:
x t 1 , x t 2 , x t 3 , x t 4 , x t 5 x t 1 , 2 x t 2 , x t 3 , x t 4 , x t 5 x t 1 , x t 2 , 3 x t 3 , x t 4 , x t 5 x t 1 , x t 2 , x t 3 , 4 x t 4 , x t 5 x t 1 , x t 2 , x t 3 , x t 4 , 5 x t 5 x t 1 , 2 x t 2 , 3 x t 3 , x t 4 , x t 5
Step 2.2: calculate the target function value of each point, it is determined that object function maximum of points and object function minimum point;
Object function is:
y = &Sigma; j = 1 m ( 1 + e - x t j ) - - - ( 1 )
In formula, y is target function value, j=1,2 ..., m, wherein, m is that positive integer represents the number gathering variable, t express time;
By P1,P2,P3,P4,P5,P6The coordinate of 6 points substitutes in formula (1) respectively, calculates the target function value y obtaining 6 points1,y2,...yn,yn+1, determine maximum through comparing, formula represents:
yh=max (yi)
In formula, yhRepresent target function value maximum;H represents yhSubscript, be used for identify maximum;I=1,2 ..., 6;
Determine that minima equation below represents through comparing:
yl=min (yi)
In formula, ylFor object function minima, l is ylSubscript, be used for identify minima;
The point P with object function maximum is determined according to above-mentioned minima and maximumhWith the point with object function minima;
Step 2.3: determining the center of gravity of point except object function maximum of points, formula is as follows:
P &OverBar; = &Sigma; i &NotEqual; h P i / n - - - ( 2 )
In formula,Center of gravity for point;
Step 2.4: object function maximum of points is carried out reflection operation, makes maximum of points move to center of gravity opposite direction, produces new some P*, formula is:
P * = ( 1 + &alpha; ) P &OverBar; - &alpha;P h - - - ( 3 )
In formula, P*For the new point produced;α is reflection coefficient, and has 0 < α < 1, α more big, reflects more remote, otherwise then more near, for making primary be distributed widely in feas ible space, takes α=0.9;
Step 2.5: calculate the target function value y of newly generated point*If, yl≤y*≤yh, then P is used*Replace Ph
If the target function value y of newly generated point*< yl, then need newly generated point is carried out Extension Operations, perform step 2.6;
If the target function value y of newly generated point*>yh, then need to carry out newly generated point shrinking computing, repeated execution of steps 2.7;
Step 2.6: newly generated point is carried out Extension Operations, and formula is as follows:
P * * = ( 1 - &gamma; ) P &OverBar; + &gamma;P * - - - ( 3 )
In formula, P**For another the new point produced, γ is flare factor, and has γ > 1;
Calculate P**Target function value y**If, P**Target function value y**> yl, then P is used**Replace Ph, perform step 2.8;
If y**< yl, then expand failure, then use P*Replace Ph;Flare factor can not be excessive, takes γ=1.2, performs step 2.8;
Step 2.7: carrying out newly generated point shrinking computing, formula is as follows:
P * * = ( 1 - &beta; ) P &OverBar; + &beta;P * - - - ( 4 )
In formula, β is constriction coefficient, and has 0 < β < 1, for making initial vertax distribution more uniform, namely gets contracting factor beta=0.5;
If P**Target function value y**< yh, then P is used**Replace Ph, perform step 2.8;
If y**> yh, it was shown that shrink unsuccessfully, from newly determined P1,P2,P3,P4,P5,P6The coordinate of 6 points, formula is as follows:
Pi=(Pi+Pl)/2(5)
Step 2.8: redefine newly generated P1,P2,P3,P4,P5,P6The center of gravity of 6 pointsFormula is as follows:
P &OverBar; &prime; = &Sigma; i &NotEqual; h P i / n
Step 2.9: the center-of-gravity value that the newly generated center-of-gravity value of step 2.8 and step 2.3 calculate compares, and formula is:
| P &prime; &OverBar; - P &OverBar; | &le; &epsiv; - - - ( 6 )
In formula, ε represents threshold value, if formula (6) is set up, then performs step 2.10, otherwise usesReplaceRepeat step 2.4~step 2.9;
Step 2.10: determineAnd PhThe midpoint x of linet
Step 2.11: set up object function, the output indication transformer sleeve pipe insulation life of object function;
Object function is:
yp=f (xt)+k(xt)
In formula, ypIndication transformer sleeve pipe insulation life;F (xt) for evaluating item, represent the function of evaluation bushing shell for transformer insulation life;Representing constraint function, in formula, T represents the specified time limit of sleeve pipe, and t' is the time that sleeve pipe uses;Step 3: the bushing shell for transformer insulation life value according to output, overhauls bushing for transformer or changes.
CN201310635349.8A 2013-11-29 2013-11-29 A kind of bushing shell for transformer insulation life prediction unit and method Expired - Fee Related CN103675621B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310635349.8A CN103675621B (en) 2013-11-29 2013-11-29 A kind of bushing shell for transformer insulation life prediction unit and method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310635349.8A CN103675621B (en) 2013-11-29 2013-11-29 A kind of bushing shell for transformer insulation life prediction unit and method

Publications (2)

Publication Number Publication Date
CN103675621A CN103675621A (en) 2014-03-26
CN103675621B true CN103675621B (en) 2016-06-29

Family

ID=50313779

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310635349.8A Expired - Fee Related CN103675621B (en) 2013-11-29 2013-11-29 A kind of bushing shell for transformer insulation life prediction unit and method

Country Status (1)

Country Link
CN (1) CN103675621B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
NZ768661A (en) 2014-10-27 2022-04-29 Landis & Gyr Ag A method, system and assembly for determining a reduction of remaining service lifetime of an electrical device during a specific time period of operation of the electrical device
CN105356847B (en) * 2015-10-30 2017-11-14 国家电网公司 A kind of distributed photovoltaic isolated island detection dynamic blind area index forecasting method
CN105375874B (en) * 2015-10-30 2017-10-27 国家电网公司 A kind of photovoltaic maximum power tracking performance index forecasting method
CN105356500B (en) * 2015-10-30 2017-11-10 国家电网公司 A kind of more photovoltaic region Network Voltage Stability nargin Forecasting Methodologies
CN105891427B (en) * 2016-06-08 2019-05-31 深圳市欧瑞博电子有限公司 Sensor life-time monitoring method and device based on cloud computing
CN109239546B (en) * 2018-08-28 2020-05-22 西南交通大学 Transformer insulation life prediction and reliability test method

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1553206A (en) * 2003-12-18 2004-12-08 西安交通大学 Electric power apparatus external insulative leakage current on-line monitoring system in converting station
CN201984124U (en) * 2011-01-28 2011-09-21 保定天威集团有限公司 Intelligent online monitoring device for dry type transformer
CN102721922A (en) * 2012-06-29 2012-10-10 沈阳工业大学 Breaker insulating coefficient prediction unit and method
CN103323756A (en) * 2013-06-21 2013-09-25 国家电网公司 Overhead line insulation filth degree prediction device and method in power distribution network
CN103323757A (en) * 2013-06-21 2013-09-25 国家电网公司 Power distribution cabinet gas insulation intensity prediction device and method
CN103336230A (en) * 2013-06-21 2013-10-02 国家电网公司 Device and method for predicting cable insulation failure rate in power distribution network
CN103336199A (en) * 2013-06-21 2013-10-02 国家电网公司 Device and method for predicting fault rate of totally-closed tube type busbar

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1553206A (en) * 2003-12-18 2004-12-08 西安交通大学 Electric power apparatus external insulative leakage current on-line monitoring system in converting station
CN201984124U (en) * 2011-01-28 2011-09-21 保定天威集团有限公司 Intelligent online monitoring device for dry type transformer
CN102721922A (en) * 2012-06-29 2012-10-10 沈阳工业大学 Breaker insulating coefficient prediction unit and method
CN103323756A (en) * 2013-06-21 2013-09-25 国家电网公司 Overhead line insulation filth degree prediction device and method in power distribution network
CN103323757A (en) * 2013-06-21 2013-09-25 国家电网公司 Power distribution cabinet gas insulation intensity prediction device and method
CN103336230A (en) * 2013-06-21 2013-10-02 国家电网公司 Device and method for predicting cable insulation failure rate in power distribution network
CN103336199A (en) * 2013-06-21 2013-10-02 国家电网公司 Device and method for predicting fault rate of totally-closed tube type busbar

Also Published As

Publication number Publication date
CN103675621A (en) 2014-03-26

Similar Documents

Publication Publication Date Title
CN103675621B (en) A kind of bushing shell for transformer insulation life prediction unit and method
CN102607643A (en) Overheat fault diagnosis and early warning system and method for electrical equipment of traction substation of electrified railway
CN203337759U (en) Online state monitoring system for MOAs (metal oxide arresters) for power transmission line
CN102967831B (en) On-line detection system and detection method of lead-acid storage battery performance
CN201555658U (en) Railway contact net vibration performance test device
CN202189454U (en) Bridge vibration data sampling system based on wireless communication technology
CN103941079A (en) On-line monitoring and fault diagnosis system for power distribution network PT
CN103336199B (en) Device and method for predicting fault rate of totally-closed tube type busbar
CN206460327U (en) The steady control tester of HF WDIT numerals
CN201724977U (en) Intelligent measurement system for grounding grid nodal voltage
CN203133201U (en) Integrated measurement device for energy-storage grid-connected system
CN202372630U (en) Wireless sensor network monitoring device for states of photovoltaic power station storage batteries
CN204694753U (en) A kind of high-precision transmission line of electricity electric current on-Line Monitor Device
CN201319049Y (en) High-tension current mutual-inductor transformation ratio tester
CN110781630A (en) GIS equipment pipeline structure temperature difference stress deformation online monitoring method
CN203489917U (en) Intelligent data acquisition system of vibrating wire sensor
CN202383267U (en) Storage battery state wireless network sensor node for DC (direct current) screen of transformer station
CN204462368U (en) A kind of portable type current transformer graph of errors proving installation
CN204666751U (en) A kind of Zinc-Oxide Arrester on-Line Monitor Device based on dsp chip
CN203365615U (en) Failure rate prediction device for totally-enclosed tubular bus
CN203630288U (en) Transformer bushing insulation life estimation device
CN103436650B (en) Thermal load wireless detection system of high precision blast furnace and detection method thereof
CN203966329U (en) The remote infrared discriminating gear of a kind of communication apparatus defect
CN203630238U (en) Device for detecting residual life of sleeve of transformer
CN211264313U (en) GIS equipment pipeline structure temperature difference stress deformation on-line monitoring device

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160629

Termination date: 20171129