CN107168070B - A kind of control method of transformer transfer load - Google Patents
A kind of control method of transformer transfer load Download PDFInfo
- Publication number
- CN107168070B CN107168070B CN201710562268.8A CN201710562268A CN107168070B CN 107168070 B CN107168070 B CN 107168070B CN 201710562268 A CN201710562268 A CN 201710562268A CN 107168070 B CN107168070 B CN 107168070B
- Authority
- CN
- China
- Prior art keywords
- transformer
- oil
- temperature
- equivalent thermal
- period
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Housings And Mounting Of Transformers (AREA)
- Protection Of Transformers (AREA)
Abstract
The technical program provides a kind of control method of transformer transfer load, the method is for controlling 220KV transformer to 110KV transformer transfer load, the method includes, predict the 220KV temperature of oil in transformer variation delta θ (t) of t period, objective function is established, finds out optimal solution using target algorithm calculating target function.The oil temperature of transformer is predicted, the objective function of transformer transfer load is established, and seeks the optimal solution of objective function, 220KV transformer transfer load is controlled based on the optimal solution.Compared with prior art, prediction oil temperature based on transformer carries out load transfer, avoid oil temperature and inside transformer component reach stablize before inside transformer component temperature it is excessively high caused by security risk, improve the safety of transformer station high-voltage side bus.
Description
Technical field
The present invention relates to power distribution network optimisation technique fields, specifically, being related to a kind of controlling party of transformer transfer load
Method.
Background technique
Transformer oil is a kind of fraction products of petroleum, its main component is alkane, and naphthene series's saturated hydrocarbons, aromatic series is not
The compounds such as saturated hydrocarbons.Transformer oil has heat spreading function, and the specific heat of transformer oil is big, is commonly used for coolant.Transformer station high-voltage side bus
When the heat that generates make to rise close to the oily expanded by heating of iron core and winding, by the convection current up and down of oil, heat passes through radiator
It sheds, guarantees that transformer operates normally.
In the prior art, it is often used the temperature that temperature of oil in transformer carrys out indicating transformer, when transformer temperature is excessively high, can be incited somebody to action
The load of the higher transformer of oil temperature is to the lower transformer transfer of oil temperature, to reduce temperature of oil in transformer.However, working as transformer
When internal component such as winding or iron core temperature increases, the oil temperature of transformer can't increase immediately, the oil temperature of transformer with
The temperature of inside transformer component, which reaches to stablize, needs certain delay, if therefore measure temperature of oil in transformer it is excessively high when born again
The component temperature of the transfer of lotus, inside transformer may be excessively high, and there are security risks.
Summary of the invention
In view of the above shortcomings of the prior art, the present invention provides a kind of control methods of transformer transfer load, to change
The oil temperature of depressor is predicted, the objective function of transformer transfer load is established, and seeks the optimal solution of objective function, is realized to change
The transfer of depressor load.Compared with prior art, the prediction oil temperature based on transformer carries out load transfer, avoids oil temperature and becomes
Security risk caused by inside transformer component temperature is excessively high before depressor internal component reaches stable, improves transformer station high-voltage side bus
Safety.
In order to solve the above-mentioned technical problem, present invention employs the following technical solutions:
A kind of control method of transformer transfer load, the method is for controlling 220KV transformer to 110KV transformer
Transfer load, which comprises
Predict that the 220KV temperature of oil in transformer variation delta θ (t) of t period, Δ P (t) are 220KV transformer in the t period
Power increment, Δ θ (t)=A (t) Δ P (t)+b (t)；
Establish objective function:
The state of 110KV line switching j after optimizing for the t period, 0 indicates that switch is to pull open state, and 1 indicates that switch is
Operating status；The state of 110KV line switching before optimizing for the t period；The beginning moment of t-th of period is current time,
That is t-1 moment, the end moment of t-th of period are future time instance, i.e. t moment；NsIndicate total 110KV line switching number；A
It (t) is t period NT1 × N of platform 220KV transformer temperature variation flow functionTSlope of step vector, diag [A (t)] are indicated column
Vector A (t) becomes diagonal matrix；The current time temperature of oil in transformer and institute that θ (t) and b (t) obtain for t period monitoring system
State the N of 220KV temperature of oil in transformer variation delta θ (t) coefficientT× 1 rank vector, PGmaxAnd θmaxFor the maximum work of 220KV transformer
The N of rate and upper layer oil temperature permissible valueT× 1 rank vector；For in period t 220KV transformer and 110KV line switching exist and open up
Flutter the load coefficient row vector of relationship, StFor switch state column vector；
The objective function, which is calculated, using target algorithm finds out optimal solution；
220KV transformer transfer load is controlled based on the optimal solution.
Preferably, the 220KV temperature of oil in transformer variation delta θ (t) includes temperature of oil in transformer variable quantity under natural conditions
Δθup(t),
In formula: Δ θaIt (t) is variation of ambient temperature amount, SnFor 220KV transformer rated capacity；The variation of nature oil temperature
Rate d Δ θup(t)/dt=[Δ θup(t)-Δθup(t-1)]/Δ t=[θup(t)-2θup(t-1)+θup(t-2)]/Δ t, Δ t are
Adjacent time interval time interval, at the time of t-2 indicates the previous period at t-1 moment, RathFor 220KV transformer air side
Equivalent thermal resistance, RothFor the equivalent thermal resistance of 220KV transformer oil side, CthFor the equivalent thermal capacitance of 220KV transformer, x is oily index,
P0For the no-load loss of 220KV transformer, PkFor the short circuit loss of 220KV transformer.
Preferably, the equivalent thermal resistance R of the 220KV transformer air sideath, 220KV transformer oil side equivalent heat
Hinder Roth, the 220KV transformer equivalent thermal capacitance CthAnd the calculation method of the oily index x is to set x, R respectivelyath、Roth、
CthFor x1、x2、x3、x4, P (t), Δ P (t), Δ θ in formula (1)a(t)、Δθup(t)、dΔθup(t)/dt can be calculated by historical data
Out, it is set to y1、y2、y3、y4、y5, n group data are equipped with, then i-th group is y1(i)、y2(i)、y3(i)、y4(i)、y5(i), will
Temperature of oil in transformer variation delta θ described in formulaup(t) calculation formula is rewritten into following form:
y4=g1(x1, x2, x3, x4, y1, y2, y3, y5)
And in practice since there are error, left and right is simultaneously unequal, if error is w, y1、y2、y3、y4、y5, w be n dimension column to
Amount,
y4=g1(x1, x2, x3, x4, y1, y2, y3, y5)+w,
According to least square method, there is Z1=[y4-g1(x1, x2, x3, x4, y1, y2, y3, y5)]T·[y4-g1(x1, x2, x3, x4,
y1, y2, y3, y5)], y1、y2、y3、y4、y5It is to be obtained by historical statistical data, is known quantity；x1, x2, x3, x4For unknown quantity, Z1It is
About x1, x2, x3, x4Function, find out Z1X when being minimized1, x2, x3, x4Value be just to x, Rath、Roth、CthMost
Good estimation.
Preferably, the 220KV temperature of oil in transformer variation delta θ (t) includes temperature of oil in transformer change under artificial hypothermia's state
Change amount Δ θdown(t),
In formula: SnFor 220KV transformer rated capacity；Nature oil temperature rate of change d Δ θdown(t)/dt=[Δ θ down
(t)-Δθdown(t-1)]/Δ t=[θdown(t)-2θdown(t-1)+θdown(t-2)]/Δ t, Δ t is between the adjacent time interval time
Every, at the time of t-2 indicates the previous period at t-1 moment, RathFor the equivalent thermal resistance of 220KV transformer air side, RothFor
The equivalent thermal resistance of 220KV transformer oil side, CthFor the equivalent thermal capacitance of 220KV transformer, x is oily index, P0For 220KV transformer
No-load loss, PkFor the short circuit loss of 220KV transformer, the temperature variation Δ θ of cold sourcei(t)=θi(t)-θi(t-1), Rir
For the equivalent thermal resistance of artificial hypothermia side, CaFor the equivalent thermal capacitance of artificial hypothermia side.
Preferably, the equivalent thermal resistance R of the 220KV transformer air sideath, 220KV transformer oil side equivalent heat
Hinder Roth, the 220KV transformer equivalent thermal capacitance Cth, the oily index x, the artificial hypothermia side equivalent thermal resistance RirAnd institute
State the equivalent thermal capacitance C of artificial hypothermia sideaCalculation method be to set x, R respectivelyath、Roth、Cth、Ca、RirFor x1, x2, x3, x4, x5,
x6, temperature of oil in transformer variation delta θ under artificial hypothermia's statedown(t) the P (t) in calculation formula, Δ P (t), Δ θi
(t)、Δθdown(t)、dΔθdown(t)/dt、dΔθa(t)/dt can be calculated by historical data, be set to y1、y2、y3、y4、y5、
y6, n group data are equipped with, then i-th group is y1(i)、y2(i)、y3(i)、y4(i)、y5(i)、y6(i), by artificial hypothermia's state
Lower temperature of oil in transformer variation delta θdown(t) calculation formula is rewritten into following form:
y4=g2(x1, x2, x3, x4, y1, y2, y3, y5, y6),
And in practice since there are error, left and right is simultaneously unequal, if error is w, y1、y2、y3、y4、y5、y6, w be n dimension column
Vector,
y4=g2(x1, x2, x3, x4, y1, y2, y3, y5, y6)+w,
According to least square method, there is Z2=[y4-g2(x1, x2, x3, x4, x5, x6, y1, y2, y3, y5, y6)]T·
[y4-g2(x1, x2, x3, x4, x5, x6, y1, y2, y3, y5, y6)],
y1、y2、y3、y4、y5、y6It is to be obtained by historical statistical data, is known quantity；x1, x2, x3, x4、x5、x6For unknown quantity,
Z2It is about x1, x2, x3, x4, x5, x6Function, find out Z2X when being minimized1, x2, x3, x4, x5, x6Value be just to x,
Rath、Roth、Cth、Ca、RirBest estimate.
In conclusion the technical program provides a kind of control method of transformer transfer load, the method is for controlling
For 220KV transformer processed to 110KV transformer transfer load, the method includes predicting the 220KV temperature of oil in transformer of t period
Variation delta θ (t) establishes objective function, finds out optimal solution using target algorithm calculating target function.To the oil temperature of transformer into
Row prediction, establishes the objective function of transformer transfer load, and seek the optimal solution of objective function, and realization turns transformer load
It moves.Compared with prior art, the prediction oil temperature based on transformer carries out load transfer, avoids oil temperature and inside transformer component
Security risk caused by inside transformer component temperature is excessively high before reaching stable, improves the safety of transformer station high-voltage side bus.
Detailed description of the invention
In order to keep the purposes, technical schemes and advantages of invention clearer, the present invention is made into one below in conjunction with attached drawing
The detailed description of step, in which:
Fig. 1 is a kind of flow chart of the embodiment 1 of the control method of transformer transfer load disclosed by the invention；
Fig. 2 is the electrical block diagram of the transformer of nature；
Fig. 3 is the simplification electrical block diagram of the transformer of nature；
Fig. 4 is the simplification electrical block diagram of the transformer of artificial hypothermia's state；
Fig. 5 is a kind of structural schematic diagram of transformer bay connection relationship；
Fig. 6 is the structural schematic diagram of another transformer bay connection relationship.
Specific embodiment
The present invention is described in further detail with reference to the accompanying drawing.
Signified temperature of oil in transformer refers to the upper layer oil temperature (top-oil temperature) of transformer in the present invention.
As shown in Figure 1, being a kind of process of the embodiment 1 of the control method of transformer transfer load disclosed by the invention
Figure, which is characterized in that the method is for controlling 220KV transformer to 110KV transformer transfer load, which comprises
S101, the 220KV temperature of oil in transformer variation delta θ (t) for predicting the t period, Δ P (t) are 220KV transformer the
The power increment of t period, Δ θ (t)=A (t) Δ P (t)+b (t)；
Oil temperature to guarantee 220KV transformer is not excessively high, and therefore, it is necessary to meet θ ' (t)=θ (t)+A (t) Δ P (t)+b
(t)≤θmax, wherein θ ' (t) is the oil temperature of prediction, θmaxFor the maximum oil temperature of permission.
S102, objective function is established:
To change and restore the remote operation number of grid operation mode as the objective function of Optimized model, i.e., with transfer
The change frequency of transformer bay switch state can be obtained as objective function in load process:
Wherein in formula:The state of 110KV line switching j after optimizing for the t period, 0 indicates that switch is to pull open state, 1 table
Show that switch is operating status；The state of 110KV line switching before optimizing for the t period；The beginning moment of t-th of period is
Current time, t-1 moment, the end moment of t-th of period are future time instance, t moment.NsIndicate total 110KV line switching
Number.
It as shown in Figures 5 and 6, is the structural schematic diagram of two kinds of transformer bay connection relationships.Fig. 5 is the direct-furnish mode of connection,
Fig. 6 is the statements based on collusion mode of connection.The C transformer of 110KV, the mode of connection of the A and B transformer of 220KV are one standby for one in Fig. 5,
I.e. wherein C transformer load P of 1 220KV transformer to 110KVcPower supply, in addition 1 220KV transformer provides it spare
Power supply.It can be seen that the 220KV transformer side line switching of equivalent preceding wiring should be operating status in figure, and 110KV transformer side
Line switching can be hot stand-by duty, can also be operating status.For the 110KV transformer direct-furnish mode of connection in Fig. 5,
The line switching S of 110KV transformer sidet 1And St 2Power grid turn is only for the optimization object of model.For 110KV transformer direct-furnish
The mode of connection, can obtain its power balance equation is,
WhereinRespectively represent the load of band after A, B transformer optimize.Line switching St 1And St 2Take 0 expression disconnected
It opens, 1 expression is taken to be closed.PcLoad after indicating the optimization of C transformer.
Similarly, power balance equation is write to Fig. 6 column,
Power balance equation can be generalized into following form, and every transformer institute's on-load is negative no more than specified
Lotus maximum value,
In formula:For in period t, there are topological relations with 110KV transformer switch with i-th of 220KV transformer
Load coefficient row vector, StFor switch state column vector.
If the t period optimize before switch state be S0,t, switch state S after optimizationt, according to power balance equation formula
(14) each transformer in optimization front and back known to the knots modification with power be,
Due to the cooperation of electric network protection and the method for operation, 110KV and following power grid must be radial pattern open loop operation, i.e. receiving end
Power grid turns for also needing to consider that the radial pattern of 110KV electric network composition constrains in optimization.Then t period 110KV line switching state
StIt should meet,
It can be found out by above-mentioned formula
The state of 110KV line switching j after optimizing for the t period, 0 indicates that switch is to pull open state, and 1 indicates that switch is
Operating status；The state of 110KV line switching before optimizing for the t period；The beginning moment of t-th of period is current time,
That is t-1 moment, the end moment of t-th of period are future time instance, i.e. t moment；NsIndicate total 110KV line switching number；A
It (t) is t period NT1 × N of platform 220KV transformer temperature variation flow functionTSlope of step vector, diag [A (t)] are indicated column
Vector A (t) becomes diagonal matrix；The current time temperature of oil in transformer and institute that θ (t) and b (t) obtain for t period monitoring system
State the N of 220KV temperature of oil in transformer variation delta θ (t) coefficientT× 1 rank vector, PGmaxAnd θmaxFor the maximum work of 220KV transformer
The N of rate and upper layer oil temperature permissible valueT× 1 rank vector；For in period t 220KV transformer and 110KV line switching exist and open up
Flutter the load coefficient row vector of relationship, StFor switch state column vector.
S103, optimal solution is found out using target algorithm calculating target function；
Target algorithm can be Newton iteration method.
S104,220KV transformer transfer load is controlled based on optimal solution；
Include the state and transformer load of each switch after optimization in optimal solution, may be implemented pair using optimal solution
The transfer of the load of transformer.
After the on-load increase of 220kV transformer institute, power can be calculated by the Current Voltage measured, and then calculate
To power increment Δ P (t), the increment of oil temperature can be then calculated, judges whether oil temperature can be more than the maximum oil temperature θ allowedmax。
If it exceeds, it is necessary to by switch operation adjustment network topology structure, redistribute load between each transforming plant main transformer,
And it is desirable that making switch number of operations minimum under the situation for meeting constraint condition, the switch state after this operation is optimal
Solution.Therefore Newton iteration method can be used to solve objective function.
In conclusion the technical program provides a kind of control method of transformer transfer load, the method is for controlling
For 220KV transformer processed to 110KV transformer transfer load, the method includes predicting the 220KV temperature of oil in transformer of t period
Variation delta θ (t) establishes objective function, finds out optimal solution using target algorithm calculating target function, target algorithm includes newton
Iterative method.The oil temperature of transformer is predicted, the objective function of transformer transfer load is established, and seeks the optimal of objective function
Solution realizes the transfer to transformer load, controls 220KV transformer transfer load based on optimal solution.Compared with prior art, base
Load transfer is carried out in the prediction oil temperature of transformer, oil temperature is avoided and inside transformer component reaches stable preceding inside transformer
Security risk caused by component temperature is excessively high improves the safety of transformer station high-voltage side bus.
It as shown in Figures 2 and 3, is the transformer of the electrical block diagram and nature of the transformer of nature
Simplify electrical block diagram, nature, that is, transformer radiates by the way that forced oil-circulation is air-cooled, as shown in Fig. 2, wherein
Transformer winding heat source is by current source qtEquivalence, including winding, iron core and stray loss；CthFor transformer interior insulation oil, winding
With the equivalent thermal capacitance of the compositions such as fuel tank, the equivalent thermal capacitance of object is directly proportional to its specific heat；Ron、RofExtremely for transformer heat transfer
The nature commutation and forced commutation thermal resistance of insulating oil；RbtIt is mixed for different temperatures insulating oil under conventional forced oil-circulation mode
Corresponding thermal resistance, and under forced oil-circulation Directed cooling mode, the cold oil of pumping hole is sent into coil, line cake and iron under a certain pressure
The oil duct of the heart, internal oil temperature distribution is more uniform, is approximately considered R at this timebt=0；RdrThe spoke to conduct heat for radiator to surrounding air
Penetrate thermal resistance, RanAnd RafFor natural convection air in wind-cooling heat dissipating and the thermal resistance of forced convection.The point-to-point transmission heat that thermal resistance indicates passes
The severe degree passed.Under the identical temperature difference, thermal resistance is smaller, and the process for transmitting heat is more violent, i.e., heat transmitting power between the two is got over
Greatly.θwFor winding temperature, θaFor environment temperature, θtoRepresent transformer top-oil temperature.Current potential represents temperature in thermoelectricity equivalent circuit
It spends, the temperature difference is bigger between two o'clock, then potential difference is bigger in circuit.
The heat transfer process of inside transformer can be obtained by Fig. 2, heat source passes through nature commutation and forced commutation for heat transfer
It gives upper layer oil reservoir, and the insulating oil of inside transformer is in forced circulation, when it is recycled to radiator, air is agitated by fan
Cooled down.Therefore thermal resistance mainly includes oily side thermal resistance and air side thermal resistance two parts, Fig. 2 can be simplified to Fig. 3.According to figure
3 column write the state equation of oil temperature:
Above formula is shown that the heat source of transformer is and no-load loss, short circuit loss and load system by Kirchhoff's current law (KCL)
The relevant function q of numbert=(P0+k2Pk)x.Taylor expansion is carried out to above formula, the pass between oil temperature increment and rate of load condensate increment k can be obtained
System:
Wherein P0, PkIt can be obtained by nameplate parameter, Δ θa、Δk、Δθto、dΔθto/dt、k0Respectively represent environment temperature
The rate of load condensate of the knots modification of degree, the knots modification of rate of load condensate, the knots modification of oil temperature, oil temperature knots modification change rate and previous moment, can
To obtain above each amount according to historical data, parameter x, C then is estimated further according to nonlinear least square methodth、Roth、Rath,
Last basis, which currently measures environment temperature and rate of load condensate knots modification, can predict the knots modification of oil temperature.
Rate of load condensate k=P/Sn, P represents the power of the practical receiving of transformer, SnRepresent the rated capacity of transformer.It will bear
Lotus rate k is expressed as P/Sn, formula (6) can be converted to accepted way of doing sth temperature of oil in transformer variation delta θ under natural conditionsup(t),
In formula: Δ θaIt (t) is variation of ambient temperature amount, SnFor 220KV transformer rated capacity；The variation of nature oil temperature
Rate d Δ θup(t)/dt=[Δ θup(t)-Δθup(t-1)]/Δ t=[θup(t)-2θup(t-1)+θup(t-2)]/Δ t, Δ t are
Adjacent time interval time interval, at the time of t-2 indicates the previous period at t-1 moment.RathFor 220KV transformer air side
Equivalent thermal resistance, RothFor the equivalent thermal resistance of 220KV transformer oil side, CthFor the equivalent thermal capacitance of 220KV transformer, x is oily index,
P0For the no-load loss of 220KV transformer, PkFor the short circuit loss of 220KV transformer.
The equivalent thermal resistance R of 220KV transformer air sideath, 220KV transformer oil side equivalent thermal resistance Roth, 220KV transformation
The equivalent thermal capacitance C of devicethAnd the calculation method of oily index x is to set x, R respectivelyath、Roth、CthFor x1、x2、x3、x4.P in formula (1)
(t)、ΔP(t)、Δθa(t)、Δθup(t)、dΔθup(t)/dt can be calculated by historical data, be set to y1、y2、y3、y4、y5,
Equipped with n group data, then i-th group is y1(i)、y2(i)、y3(i)、y4(i)、y5(i), by formula temperature of oil in transformer variation delta θup
(t) calculation formula is rewritten into following form:
y4=g1(x1, x2, x3, x4, y1, y2, y3, y5)
And in practice since there are error, left and right is simultaneously unequal, if error is w, y1、y2、y3、y4、y5, w be n dimension column to
Amount,
y4=g1(x1, x2, x3, x4, y1, y2, y3, y5)+w,
According to least square method, there is Z1=[y4-g1(x1, x2, x3, x4, y1, y2, y3, y5)]T·[y4-g1(x1, x2, x3, x4,
y1, y2, y3, y5)], y1、y2、y3、y4、y5It is to be obtained by historical statistical data, is known quantity；x1, x2, x3, x4For unknown quantity, Z1It is
About x1, x2, x3, x4Function.Find out Z1X when being minimized1, x2, x3, x4Value be just to x, Rath、Roth、CthMost
Good estimation specifically asks the method being most worth to have the optimal methods such as Newton iteration method, interior point method.
As shown in figure 4, be the simplification electrical block diagram of the transformer of artificial hypothermia's state, the temperature of artificial cold source in figure
Degree is θi, CaFor the equivalent thermal capacitance of air, RirThe radiation thermal resistance conducted heat for cold source to surrounding air.The air side in equivalent circuit
Air themperature θ can be obtained with Kirchhoff's current law (KCL)aState equation:
To above formula Taylor expansion, increment expression formula can be obtained:
Bringing formula (8) into formula (6) can obtain
The rate of load condensate k of above formula is indicated with power, temperature of oil in transformer variation delta θ under artificial hypothermia's state can be obtaineddown
(t),
In formula: SnFor 220KV transformer rated capacity；Nature oil temperature rate of change d Δ θdown(t)/dt=[Δ θ down
(t)-Δθdown(t-1)]/Δ t=[θdown(t)-2θdown(t-1)+θdown(t-2)]/Δ t, Δ t is between the adjacent time interval time
Every, at the time of t-2 indicates the previous period at t-1 moment, RathFor the equivalent thermal resistance of 220KV transformer air side, RothFor
The equivalent thermal resistance of 220KV transformer oil side, CthFor the equivalent thermal capacitance of 220KV transformer, x is oily index, P0For 220KV transformer
No-load loss, PkFor the short circuit loss of 220KV transformer, the temperature variation Δ θ of cold sourcei(t)=θi(t)-θi(t-1), Rir
For the equivalent thermal resistance of artificial hypothermia side, CaFor the equivalent thermal capacitance of artificial hypothermia side.
The equivalent thermal resistance R of 220KV transformer air sideath, 220KV transformer oil side equivalent thermal resistance Roth, 220KV transformation
The equivalent thermal capacitance C of deviceth, oily index x, artificial hypothermia side equivalent thermal resistance RirAnd the equivalent thermal capacitance C of artificial hypothermia sideaCalculating side
Method is to set x, R respectivelyath、Roth、Cth、Ca、RirFor x1, x2, x3, x4, x5, x6, temperature of oil in transformer variable quantity under artificial hypothermia's state
Δθdown(t) the P (t) in calculation formula, Δ P (t), Δ θi(t)、Δθdown(t)、dΔθdown(t)/dt、dΔθa(t)/dt
It can be calculated by historical data, be set to y1、y2、y3、y4、y5、y6, n group data are equipped with, then i-th group is y1(i)、y2(i)、y3
(i)、y4(i)、y5(i)、y6(i), by temperature of oil in transformer variation delta θ under artificial hypothermia's statedown(t) calculation formula is rewritten
At following form:
y4=g2(x1, x2, x3, x4, y1, y2, y3, y5, y6),
And in practice since there are error, left and right is simultaneously unequal, if error is w, y1、y2、y3、y4、y5、y6, w be n dimension column
Vector,
y4=g2(x1, x2, x3, x4, y1, y2, y3, y5, y6)+w,
According to least square method, there is Z2=[y4-g2(x1, x2, x3, x4, x5, x6, y1, y2, y3, y5, y6)]T·
[y4-g2(x1, x2, x3, x4, x5, x6, y1, y2, y3, y5, y6)],
y1、y2、y3、y4、y5、y6It is to be obtained by historical statistical data, is known quantity；x1, x2, x3, x4、x5、x6For unknown quantity,
Z2It is about x1, x2, x3, x4, x5, x6Function, find out Z2X when being minimized1, x2, x3, x4, x5, x6Value be just to x,
Rath、Roth、Cth、Ca、RirBest estimate, specifically ask the method being most worth to have the optimal methods such as Newton iteration method, interior point method.
For guarantee transformer oil temperature in safe range, therefore, it is necessary to meet:
θ′up(t)=θup(t)+Aup(t)ΔP(t)+bup(t)≤θmax (10)
θ′down(t)=θdown(t)+Adown(t)ΔP(t)+bdown(t)≤θmax (11)
Finally, it is stated that the above examples are only used to illustrate the technical scheme of the present invention and are not limiting, although passing through ginseng
According to the preferred embodiment of the present invention, invention has been described, it should be appreciated by those of ordinary skill in the art that can
To make various changes to it in the form and details, without departing from the present invention defined by the appended claims
Spirit and scope.
Claims (5)
1. a kind of control method of transformer transfer load, which is characterized in that the method for control 220KV transformer to
110KV transformer transfer load, which comprises
Predict that the 220KV temperature of oil in transformer variation delta θ (t) of t period, Δ P (t) are function of the 220KV transformer in the t period
Rate increment, Δ θ (t)=A (t) Δ P (t)+b (t)；
Establish objective function:
The state of 110KV line switching j after optimizing for the t period, 0 indicates that switch is to pull open state, and 1 indicates that switch is operation
State；The state of 110KV line switching before optimizing for the t period；The beginning moment of t-th of period is current time, i.e. t-
1 moment, the end moment of t-th of period are future time instance, i.e. t moment；NsIndicate total 110KV line switching number；A(t)
For t period NT1 × N of platform 220KV transformer temperature variation flow functionTSlope of step vector, diag [A (t)] are indicated column vector
A (t) becomes diagonal matrix；θ (t) and b (t) are the current time temperature of oil in transformer that t period monitoring system obtains and described
The N of 220KV temperature of oil in transformer variation delta θ (t) coefficientT× 1 rank vector, PGmaxAnd θmaxFor the maximum power of 220KV transformer
With the N of upper layer oil temperature permissible valueT× 1 rank vector；There is topology for the 220KV transformer in period t and 110KV line switching
The load coefficient row vector of relationship, StFor switch state column vector；
The objective function, which is calculated, using target algorithm finds out optimal solution；
220KV transformer transfer load is controlled based on the optimal solution.
2. the control method of transformer transfer load as described in claim 1, which is characterized in that the 220KV transformer oil
Warm variation delta θ (t) includes temperature of oil in transformer variation delta θ under natural conditionsup(t),
In formula: Δ θaIt (t) is variation of ambient temperature amount, SnFor 220KV transformer rated capacity；Nature oil temperature rate of change d Δ
θup(t)/dt=[Δ θup(t)-Δθup(t-1)]/Δ t=[θup(t)-2θup(t-1)+θup(t-2)] when/Δ t, Δ t are adjacent
Section time interval, at the time of t-2 indicates the previous period at t-1 moment, RathFor the equivalent heat of 220KV transformer air side
Resistance, RothFor the equivalent thermal resistance of 220KV transformer oil side, CthFor the equivalent thermal capacitance of 220KV transformer, x is oily index, P0For
The no-load loss of 220KV transformer, PkFor the short circuit loss of 220KV transformer.
3. the control method of transformer transfer load as claimed in claim 2, which is characterized in that the 220KV transformer is empty
The equivalent thermal resistance R of gas sideath, 220KV transformer oil side equivalent thermal resistance Roth, the 220KV transformer equivalent thermal capacitance
CthAnd the calculation method of the oily index x is to set x, R respectivelyath、Roth、CthFor x1、x2、x3、x4, P (t), Δ P in formula (1)
(t)、Δθa(t)、Δθup(t)、dΔθup(t)/dt can be calculated by historical data, be set to y1、y2、y3、y4、y5, it is equipped with n group
Data, then i-th group is y1(i)、y2(i)、y3(i)、y4(i)、y5(i), by temperature of oil in transformer variation delta θ described in formulaup(t)
Calculation formula be rewritten into following form:
y4=g1(x1, x2, x3, x4, y1, y2, y3, y5)
And in practice since there are error, left and right is simultaneously unequal, if error is w, y1、y2、y3、y4、y5, w be n dimensional vector,
y4=g1(x1, x2, x3, x4, y1, y2, y3, y5)+w,
According to least square method, there is Z1=[y4-g1(x1, x2, x3, x4, y1, y2, y3, y5)]T·[y4-g1(x1, x2, x3, x4, y1,
y2, y3, y5)], y1、y2、y3、y4、y5It is to be obtained by historical statistical data, is known quantity；x1, x2, x3, x4For unknown quantity, Z1It is to close
In x1, x2, x3, x4Function, find out Z1X when being minimized1, x2, x3, x4Value be just to x, Rath、Roth、CthIt is best
Estimation.
4. the control method of transformer transfer load as described in claim 1, which is characterized in that the 220KV transformer oil
Warm variation delta θ (t) includes temperature of oil in transformer variation delta θ under artificial hypothermia's statedown(t),
In formula: SnFor 220KV transformer rated capacity；Nature oil temperature rate of change d Δ θdown(t)/dt=[Δ θ down(t)-
Δθdown(t-1)]/Δ t=[θdown(t)-2θdown(t-1)+θdown(t-2)]/Δ t, Δ t are adjacent time interval time interval, t-2
At the time of indicating the previous period at t-1 moment, RathFor the equivalent thermal resistance of 220KV transformer air side, RothFor 220KV change
The equivalent thermal resistance of depressor oil side, CthFor the equivalent thermal capacitance of 220KV transformer, x is oily index, P0For the zero load of 220KV transformer
Loss, PkFor the short circuit loss of 220KV transformer, the temperature variation Δ θ of cold sourcei(t)=θi(t)-θi(t-1), RirIt is artificial
The equivalent thermal resistance of cooling side, CaFor the equivalent thermal capacitance of artificial hypothermia side.
5. the control method of transformer transfer load as claimed in claim 4, which is characterized in that the 220KV transformer is empty
The equivalent thermal resistance R of gas sideath, 220KV transformer oil side equivalent thermal resistance Roth, the 220KV transformer equivalent thermal capacitance
Cth, the oily index x, the artificial hypothermia side equivalent thermal resistance RirAnd the equivalent thermal capacitance C of the artificial hypothermia sideaCalculating side
Method is to set x, R respectivelyath、Roth、Cth、Ca、RirFor x1, x2, x3, x4, x5, x6, temperature of oil in transformer becomes under artificial hypothermia's state
Change amount Δ θdown(t) the P (t) in calculation formula, Δ P (t), Δ θi(t)、Δθdown(t)、dΔθdown(t)/dt、dΔθa
(t)/dt can be calculated by historical data, be set to y1、y2、y3、y4、y5、y6, n group data are equipped with, then i-th group is y1(i)、y2
(i)、y3(i)、y4(i)、y5(i)、y6(i), by temperature of oil in transformer variation delta θ under artificial hypothermia's statedown(t) calculating
Formula is rewritten into following form:
y4=g2(x1, x2, x3, x4, y1, y2, y3, y5, y6),
And in practice since there are error, left and right is simultaneously unequal, if error is w, y1、y2、y3、y4、y5、y6, w be n dimensional vector,
y4=g2(x1, x2, x3, x4, y1, y2, y3, y5, y6)+w,
According to least square method, there is Z2=[y4-g2(x1, x2, x3, x4, x5, x6, y1, y2, y3, y5, y6)]T·
[y4-g2(x1, x2, x3, x4, x5, x6, y1, y2, y3, y5, y6)],
y1、y2、y3、y4、y5、y6It is to be obtained by historical statistical data, is known quantity；x1, x2, x3, x4、x5、x6For unknown quantity, Z2It is
About x1, x2, x3, x4, x5, x6Function, find out Z2X when being minimized1, x2, x3, x4, x5, x6Value be just to x, Rath、
Roth、Cth、Ca、RirBest estimate.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710562268.8A CN107168070B (en) | 2017-07-11 | 2017-07-11 | A kind of control method of transformer transfer load |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710562268.8A CN107168070B (en) | 2017-07-11 | 2017-07-11 | A kind of control method of transformer transfer load |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107168070A CN107168070A (en) | 2017-09-15 |
| CN107168070B true CN107168070B (en) | 2019-04-26 |
Family
ID=59823784
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710562268.8A Active CN107168070B (en) | 2017-07-11 | 2017-07-11 | A kind of control method of transformer transfer load |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107168070B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110361109B (en) * | 2019-07-18 | 2020-09-04 | 杭州电力设备制造有限公司 | Temperature calculation method, system and device for indoor transformer substation |
| WO2023004808A1 (en) * | 2021-07-30 | 2023-02-02 | 西门子（中国）有限公司 | Optimization method and apparatus for running of substation |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2159284A1 (en) * | 1994-09-29 | 1996-03-30 | Jack Trainor | Variamp oil temperature control |
| CN101752870A (en) * | 2010-01-29 | 2010-06-23 | 天津大学 | Method for analysis of available power supply capacity of medium voltage distribution network |
| CN104919380A (en) * | 2012-11-19 | 2015-09-16 | Abb技术有限公司 | Profiling transformer of power system |
| CN105975709A (en) * | 2016-05-16 | 2016-09-28 | 中国石油大学(华东) | Transformer hotspot temperature predicting method for multi-working-condition-parameter recognition and optimization |
-
2017
- 2017-07-11 CN CN201710562268.8A patent/CN107168070B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2159284A1 (en) * | 1994-09-29 | 1996-03-30 | Jack Trainor | Variamp oil temperature control |
| CN101752870A (en) * | 2010-01-29 | 2010-06-23 | 天津大学 | Method for analysis of available power supply capacity of medium voltage distribution network |
| CN104919380A (en) * | 2012-11-19 | 2015-09-16 | Abb技术有限公司 | Profiling transformer of power system |
| CN105975709A (en) * | 2016-05-16 | 2016-09-28 | 中国石油大学(华东) | Transformer hotspot temperature predicting method for multi-working-condition-parameter recognition and optimization |
Non-Patent Citations (1)
| Title |
|---|
| 并网光伏发电站在线监测与状态评估系统;周念成等;《电力自动化设备》;20121031;第32卷(第10期);第6-10页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107168070A (en) | 2017-09-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103779059B (en) | A kind of dynamic compatibilization method of oil-filled transformer | |
| CN108037780B (en) | Oil-immersed transformer cooling control method based on temperature rise and load factor | |
| JP5110603B2 (en) | Direct load control system | |
| CN105684109B (en) | The method that operation to transformer cooling system is optimized, corresponding system and the method for determining VFD capacity | |
| Martin et al. | Corrective model-predictive control in large electric power systems | |
| CN104484569B (en) | Hot-spot temperature of transformer computational methods based on thermoelectricity analogy theory | |
| CN107168070B (en) | A kind of control method of transformer transfer load | |
| CN104281057A (en) | Composite PID fuzzy control method applied to transformer cooling system | |
| CN107560757B (en) | Method and system for estimating upper oil temperature of air-cooled transformer | |
| CN101326696A (en) | A converter station and a method for control thereof | |
| US10389123B2 (en) | Power control of microgrid wherein detecting a need for overloading a transformer based on power balance within the microgrid | |
| CN201804282U (en) | Circulating cooling water inlet water temperature drift control device for operation test of converter valve | |
| CN105404780B (en) | A kind of transformer overload capability assessment method of multiparameter exponential family | |
| Weekes et al. | Risk assessment using transformer loss of life data | |
| CN108292558B (en) | For controlling the method and system of the cooling system of power equipment | |
| CN116605098A (en) | Battery thermal management method and device for vehicle, storage medium and vehicle | |
| CN109556176A (en) | A kind of heating terminal intelligent on-off valve regulation method based on dual time-step | |
| Ben-gang et al. | The improved thermal-circuit model for hot-spot temperature calculation of oil-immersed power transformers | |
| BL et al. | Asset management of transformer based on loss of life calculation | |
| CN114142476A (en) | Distribution network load transfer method considering transformer oil temperature | |
| Lai et al. | Prediction of top oil temperature for oil-immersed transformer based on Kalman filter algorithm | |
| JP2017050911A (en) | Operation plan calculation device and calculation method, and computer program | |
| CN115438512A (en) | Transformer winding hot spot temperature prediction method and system | |
| Rajender et al. | Transformer rating for solar PV plants based on overloading capability as per guidelines | |
| CN106920641A (en) | A kind of method of use air-duct-type air-conditioner to dry-type transformer forced cooling |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20190820 Address after: 253200, Shandong, Xiajin County, Dezhou Province Economic Development Zone on the north side of the street, west side of Ji Yun Street Patentee after: SHANDONG LUKONG ELECTRIC POWER EQUIPMENT Co.,Ltd. Address before: 400044 Shapingba District Sha Street, No. 174, Chongqing Patentee before: Chongqing University |
|
| TR01 | Transfer of patent right | ||
| PE01 | Entry into force of the registration of the contract for pledge of patent right |
Denomination of invention: A control method of transformer load transfer Effective date of registration: 20220705 Granted publication date: 20190426 Pledgee: Shandong Xiajin Rural Commercial Bank Co.,Ltd. Pledgor: SHANDONG LUKONG ELECTRIC POWER EQUIPMENT Co.,Ltd. Registration number: Y2022980009950 |
|
| PE01 | Entry into force of the registration of the contract for pledge of patent right |