The drops such as the reactive-load compensation based on rural power grids typical case's taiwan area line loss calculation damage optimization method
Technical field
The invention belongs to transmission & distribution electro-technical field, more particularly, to a kind of based on rural power grids typical case's taiwan area line loss calculation
The drops such as reactive-load compensation damage optimization method.
Background technology
In recent years, China's economy develops rapidly, and power industry is also developed faster therewith, and China's machine-building is horizontal
Also it achieves compared with much progress, first and second electrical equipment manufacturing technology of power grid also obtains surprising progress.By the electricity of national several years
Net upgrading, power grid equipment, technical performance acquirement are greatly improved.For many years, for the Controlling line loss of Utilities Electric Co.,
The especially Controlling line loss of rural area distribution is always the difficult point and emphasis of Controlling line loss.The related statistics table of many access
Bright, in each voltage class power grid in China, the electric quantity loss of 10KV and following distribution has accounted for the proportion of whole distract electric quantity loss
About 49%, how power supply enterprise well reduces the electric energy loss of power distribution network, is the long-term mesh of management in power supply enterprise work
Mark and focus.
Each power supply enterprise is using the statistics of line loss, analysis and management as the day-to-day operations sex work of enterprise.Management and control is each
Class electric energy and line losses indices are the main means of line loss statistical work, carry out line loss statistics, analysis, management work, find analysis
The various correlative factors that can influence line loss per unit variation, related conclusions, adaptation to local conditions, for actual conditions are obtained by line loss analyzing
Relevant reducing loss measure is formulated, the management in power supply enterprise person where power distribution network is allowed really to make reality in administered operation of power networks clear
Power levels are lost.With Electric network energy loss theoretical calculation and analysis for handgrip, determines the variable loss of power transmission network and power distribution network and consolidate
Setting loss loss-rate calculates, analysis, grasps line loss changing rule, comparison theory wire loss value and practical line loss value, comparison ring than data with
Year-on-year data compare similar circuit, to find the cause for influencing line loss variation;Generally going through line loss theoretical calculation can be clear
Determine technical loss and manage line loss ratio, for formulate drop damage index measure and performance assessment criteria measure lay a solid foundation.Cause
This, can effectively verify whether current Distributing network structure is reasonable, and carry for the transformation of distribution by line loss theoretical calculation and analysis
For correlation theory foundation.
Applicants have found that reactive-load compensation can accomplish that idle in-situ balancing, reduction is made to convey from mains side idle
Electric power;And compensation benefit convenient for installation and maintenance, at low cost it is good the advantages that.Reactive-load compensation has the same of compensating power effect
When also have the function of adjust voltage.Reactive-load compensation drop damage on the basis of plus improve three-phase imbalance, replace transformer,
Replacing the reducing loss measures such as circuit section product can make line loss be reduced to ideal result.
Therefore, the proposition of the drops such as a kind of reactive-load compensation based on rural power grids typical case's taiwan area line loss calculation damage optimization method, has
Certain theoretical foundation and realistic meaning.
Invention content
The present invention for the above-mentioned prior art the problem of, provide and a kind of being based on rural power grids typical case's taiwan area line loss calculation
The drops such as reactive-load compensation damage optimization method.The purpose is to prepare for more accurate line loss calculation, agriculture is more accurately calculated
The loss of net platform region.
In order to achieve the above-mentioned object of the invention, the present invention is achieved by the following technical solutions:
The drops such as the reactive-load compensation based on rural power grids typical case's taiwan area line loss calculation damage optimization method, include the following steps:
Step 1:The topological structure for analyzing network is accurate using dividing gradient search procedure to obtain accurate physical parameter
Line loss calculation is prepared;According to the different characteristics of rural power grids typical case's taiwan area and radius of electricity supply, transformer model capacity service life, bear
The data cleansing and excavation of the physical parameters such as lotus distribution situation carry out typicalness classification to rural power grids taiwan area;
Step 2:Using the improved loss for being pushed forward trend method and calculating rural power grids taiwan area;
Step 3:It establishes the object function of idle work optimization and adds voltage, power factor, compensation capacity, the pact for compensating equipment
Beam condition;
Step 4:It determines the compensation position of the reactive-load compensation of different type typical case's taiwan area and each compensates what equipment should compensate
Capacity;
Step 5:Establish the rural power grids taiwan area of more reducing loss measure combinations;
Step 6:Establish the reactive compensation capacity of decomposition-coordination approach.
The step 1 includes:According to the network structure of distribution network system and switch operating status, to the load of power distribution network into
Row classification, according to different loads and whether there is or not measuring values to carry out a point gradient search, to obtain the parameter and section of each branch of system
Voltage, active power and the reactive power of point;
The active and reactive power and voltage value of main transformer, distribution transforming in first gradient search distribution feeder;
Transformer for failing to obtain real-time measurement carries out the second gradient search, and complete and precision satisfaction is wanted in order to obtain
The operating parameter asked, the use state estimation technique carry out line loss calculation;
In the quantity for having searched for transformer, 3rd gradient search is carried out, is searched on user and overhead transmission line, cable run
Each node forms node matrix equation.
Described in the step 2 it is improved be pushed forward trend method refer to taken as measurement to node voltage value and calculating obtained by
As a result average value makes the voltage value of each node more accurate, and the line loss value of calculating is truer;
The step 2 includes:By obtaining the measured data of each end load node, each section is calculated by prospective method
The voltage value and performance number of the corresponding last node of point, if as soon as upper node uses without measured value and calculate forward calculation institute
The voltage value obtained takes the average value of measured value and calculated value if any measured value;If a upper node corresponds to two load bus,
If this node takes two loads to be pushed forward the average value of voltage without measured value, the average value of three's value is otherwise taken, and so on;Directly
To root node;The result of line loss is equal to the total losses that transformer calculates and adds the total losses calculated on circuit;
Aloss=Tloss+Lloss
A in above formulalossFor all line loss, TlossFor the loss of transformer, LlossFor the loss of transmission line of electricity;
Step 21:Known to line impedance parameter, as long as measuring active and reactive, the voltage of load bus
Virtual value, so that it may accurately to find out very much loss between load bus and its superior node on circuit and superior node
Voltage, and calculating line loss and node voltage all avoid plural multiplication and division;Circuit l between node 2 to 32-3It is per phase of impedance
Z2-3=R2-3+jX2-3The three-phase output power for measuring node 3 simultaneously is P3+jQ3, line voltage virtual value is U3;
Circuit l2-3On active power loss be:
In above formula, Δ P2-3Refer to circuit l2-3On active loss;P3It is the active power of node 3, Q3It is node 3
Reactive power;U3It is the voltage of node 3, R2-3It is circuit l2-3On resistance value;
Reactive power loss is:
In above formula, Δ Q2-3Refer to circuit l2-3On active loss;P3It is the active power of node 3, Q3It is node 3
Reactive power;U3It is the voltage of node 3, R2-3It is circuit l2-3On resistance value;
The output power of node 2 is:
S2=(P3+ΔP2-3)+j(Q3+ΔQ2-3)
In above formula, S2Refer to the output power of node 2;P3It is the active power of node 3, Q3It is the idle work(of node 3
Rate;ΔP2-3、ΔQ2-3It is circuit l respectively2-3The active power and reactive power of upper loss;
The voltage effective value U of node 22For:
In above formula, U2Refer to the RMS voltage Δ U of node 22-3、δU2-3It is circuit l2-3The vertical component of voltage landing and
Horizontal component;
Step 22:Assuming that a certain taiwan area, which shares 11 nodes, it is known that resistance, reactance between each node, and
And the active power, reactive power, voltage effective value of each load bus are had been obtained for by data collecting system;
Steps are as follows for the specific implementation of algorithm:
(1) node-classification is carried out first and is numbered, and power supply node, load bus are distinguished;Power supply node is 0 node, load
Node is 2,3,8,9,10 nodes;
(2) by obtaining the measured data of each end load node, by prospective method calculate each node it is corresponding on
The voltage value and performance number of a node, if as soon as upper node uses the voltage value calculated obtained by forward calculation without measured value,
If any measured value, the average value of measured value and calculated value is taken;If a upper node corresponds to two load bus, if this node without
Measured value, takes two loads to be pushed forward the average value of voltage, otherwise takes the average value of three's value, and so on;Until power supply node;
Circuit between node 8 and 6 is l8-6Each impedance between circuit is Z8-6=R8-6+jX8-6The load 8 measured saves
Three output powers of point are P8+jQ8Then circuit l8-6Between active power loss:
In above formula, Δ P8-6It is circuit l8-6Active power loss, P8It is the active power of node 8, Q8Node 8 it is idle
Power, U8It is the voltage effective value of node 8, R8-6It is circuit l8-6Between impedance;
Reactive power loss is:
In above formula, Δ Q8-6It is circuit l8-6Reactive power loss, P8It is the active power of node 8, Q8Node 8 it is idle
Power, U8It is the voltage effective value of node 8, X8-6It is circuit l8-6Between reactance;
The output power of node 6 is:
S6=P6+jQ6=(P8+ΔP8-6)+j(Q8+ΔQ8-6)
In above formula, S6Refer to the output power of node 6;P6It is the active power of node 6, Q6It is the idle work(of node 6
Rate;ΔP8-6、ΔQ8-6It is circuit l respectively8-6The active power and reactive power of upper loss;P8、Q8It is the wattful power of node 8 respectively
Rate and reactive power;
The voltage value of node 6 is:
In above formula, U8For the voltage effective value of node 8, Δ U8-6For circuit l8-6The vertical component of voltage landing, δ U8-6For line
Road l8-6The horizontal component of voltage landing;
Wherein:
Circuit l can similarly be passed through9-6The voltage and output power for calculating node 6, if the voltage of node 6 can be measured
Value obtains the voltage value and performance number of node 6 with output power by calculating the average value of three, if the node 6 that do not survey
The voltage value of voltage value and performance number, node 6 takes the average value of calculating, performance number also to take the average value of calculating;
(3) circuit having had been calculated is marked, continue search for the circuit of non-marked until all circuits all
It calculates;
(4) it calculates the loss of all common transformers, when calculating, the node load of actual measurement, voltage value is utilized to calculate transformation
The copper loss and iron loss of device;So, transformer loss calculates the influence for just having fully taken into account load, voltage;
(5) total line loss:
Aloss=Tloss+Lloss
Wherein AlossFor all line loss, TlossFor the loss of transformer, LlossFor the loss of transmission line of electricity.
Object function described in the step 3:
min(P1-P2)
Wherein P1 is investment cost, and P2 is drop damage income, and object function is the minimum value of two differences;
Wherein:
In above formula, n is that node compensates equipment sum;NiFor the packet count of the compensation capacitor of the i-th complete equipment;Ci is compensation
The capacity of the corresponding single group capacitor of equipment;Ma is the unit price of compensation capacity;Mb is the single group price for compensating equipment;Mc is
Expense input by equipment is installed;
Wherein:
P2=(f1-f2)·c·T
In above formula, f1For the power loss before compensation, f2For the power loss after compensation, c is corresponding electricity price, T operations
Time;
Specifically:To minimize the active loss of low-voltage network as target, low pressure:Distribution network var compensation is distributed
Optimized model is as follows:
min f(x)
S.t h (x)=0,
In above formula:Object function f (x) is the difference of investment cost and drop damage income;H (x) is that the power of each node is flat
Weigh equation, and x is the decision variable and state variable of system, and decision variable is the parallel reactive compensation capacity Q of reactive compensation point, shape
State variable includes the voltage magnitude V and phase angle δ of each node;Inequality constraints l (x) include the bound of each node voltage amplitude about
The bound of beam and each reactive compensation point parallel reactive compensation capacity constrains and the bound of power factor is divided with compensation equipment
The bound of group number;
The difference of object function f (x) investment costs and drop damage income, expression formula are as follows:
F (x)=P1-P2
Step 32:The constraints;Equality constraint h (x) expression formulas are:
It can be obtained by Load flow calculation equation:
In above formula, PiAnd QiIt is the active power and reactive power of taiwan area injection respectively, it is negative, Q to flow into just, to flow outci
For the reactive power of compensation, ViFor the voltage magnitude of node i, δijFor the phase difference of voltage between node i and node j, GijAnd Bij
It is the real and imaginary parts of network node admittance matrix the i-th row j column elements respectively, n is the sum of taiwan area node;
Inequality constraints l (x) expression formulas are:
V in above formulaimin、VimaxIt is the minimum voltage and maximum voltage that node allows, Qcimin、QcimaxIt is node to allow most
It is small idle with maximum reactive-load compensation compensation, ρimin、ρimaxIt is the minimum and maximum power factor that node measured place allows,It is the minimum packets number and largest packet number for compensating equipment switching.
The step 4 includes:
Step 41:The compensation rate that taiwan area all nodes in Typical Rural respectively need is calculated first, calculates taiwan area needs
The progress reactive-load compensation of each node is not met reality, using Sensitivity Analysis Method and artificially by total reactive compensation capacity
Factor method selects the node that the node that load power is big, important is compensated as needs;Every in population is mended using genetic algorithm
It repays equipment and sets 2 genes on chromosome, indicate the packet count and every group of compensation capacity of compensation equipment respectively, calculate egress
Optimal non-workcompensation capacity;
Step 42:When calculating the idle configuration optimization of different type taiwan area on the basis of step 41, three are defined first
Electric characteristic index circuit distribution transforming Rate of average load α, natural power factor cos φ, radius of electricity supply L, reactive-load compensation rate β, according to
《Rural power planning and designing directive/guide》With the specific electric characteristic index system in local rural area, typical rural area electrical specifications are set
As ground state value;
In above formula, β0For the optimal reactive compensation rate of Typical Rural taiwan area, WnHold for the idle configuration that the taiwan area requires supplementation with
Amount, SNFor the rated capacity of taiwan area;
When calculating the reactive compensation capacity of random taiwan area by changing each electric characteristic index respectively, OPTIMAL REACTIVE POWER is recorded
The situation of change of cancellation ratio indicates that the variation degree of reactive-load compensation rate, the size of sensitivity are equal to respectively with the size of sensitivity
The variation of the corresponding reactive-load compensation rate of electric characteristic variation exists in the reactive capability configuration for carrying out random taiwan area, with each electrical
The difference of the actual value of characteristic value electric characteristic corresponding with typical taiwan area carries out phase as decision using the size of sensitivity
The weighted calculation answered;
β in above formula0For the reactive-load compensation rate of typical taiwan area, Δ miIt is right for random taiwan area and typical taiwan area electrical characteristics
The difference answered, λiFor the sensitivity corresponding to electrical characteristics.
The step 5 includes:
Step 51:Determine that the reducing loss measure of rural power grids taiwan area is denoted as X firsti, wherein i=1,2 ..., 4;X1~X4It is right with this
Reactive-load compensation, improvement three-phase imbalance, replacement transformer capacity, replacement circuit section is answered to accumulate these four rural power grids taiwan areas drop damage and arrange
Apply type;On this basis, scheme is damaged using the possible embodiment of certain a kind of reducing loss measure as the drop in rural power grids taiwan area drop damage,
AsIndicate the jth kind possible embodiment of i-th kind of reducing loss measure type;X1The corresponding drop damage scheme of reducing loss measure has transformation
It stands centralized compensation, user's decentralised control, Optimal Reactive Power compensation device;X2The corresponding drop damage scheme of reducing loss measure has change client to connect
Line phase sequence, three-phase imbalance pressure regulator, auto-transformer;X3The corresponding drop damage scheme of reducing loss measure has volume reduction transformation method, transformation
Device way of economic operation, energy saving transformer;X4The corresponding drop damage scheme of reducing loss measure has parallel gathering lead wire in low, increases conducting wire section
Area, line insulation transformation;
For reducing loss measure Xi, since the presence of analysis of Influential Factors is lost, from reducing loss measure weighing factor analysis and assessment
As a result the larger reducing loss measure type of weight is selected in, and alternately reducing loss measure is implemented present situation rural power grids taiwan area, thus,
To reducing loss measure XiWhether selected available following formula describes;
For selected reducing loss measure, it is further contemplated that its possible reducing loss measure implements project, it is possible to
Damage project dropsFor the description of the implementation principle of rural power grids taiwan area reducing loss measure, according to artificial experience and different existing optimizations
Algorithm proposes a variety of alternative implementation projects to every a kind of reducing loss measure type, and following formula is in i-th kind of reducing loss measure
The whether selected description made of jth kind drop damage project;
Step 52:In order to ensure that drop damage reaches better effect drop damage alternative is formed by following steps;
Drop damage alternative is formed by drop damage projects combo, since the presence of analysis of Influential Factors is lost in rural power grids taiwan area,
The major influence factors of the influence rural power grids taiwan area technological attrition in present situation rural power grids taiwan area and its corresponding can be primarily determined
Reducing loss measure type, thus, can determine the situation of choosing of reducing loss measure generating the when of dropping damage alternative, it is assumed that warp
Selected reducing loss measure is respectively X after crossing loss analysisa、Xb、Xc, then drop damages the composition structure of candidate item such as in this case
Under:
In the case of the implementation project determination of every a kind of reducing loss measure type, a drop damage candidate item is drop damage project
The combination of selected situation, while also meeting drop damage projectBetween alternative, be shown below:
In above formula:fiRefer to i-th of drop damage candidate item;G refers to the alternative between drop damage project, i.e., same kind of drop damage is arranged
Implementation of the different implementation projects applied in rural power grids taiwan area is impossible simultaneous constraint;Mutex relation constraint master between project
If being embodied in drop damages decision process;
The product process of rural power grids taiwan area drop damage candidate item is as follows:
(1) the loss shadow that key important function is played to rural power grids taiwan area technical loss is established by rural power grids taiwan area loss analysis
The factor of sound;
(2) it is based on sensitivity analysis, establishes the reducing loss measure X in rural power grids taiwan area Optimal Decision-makingiChoose situation;
(3) to each selected reducing loss measure Xi, obtained for choosing by the research of each measure scheme weight analysis
The drop damage project selected
(4) according to the initial rural power grids taiwan area drop damage candidate item f of the various combination producings of drop damage projecti;
Step 53:In damage decision model drops in rural power grids taiwan area, it is present situation rural power grids taiwan area band that should fully embody reducing loss measure
The comprehensive benefit of the drop damage electricity come is established by analyzing Evaluation of Investment-Benefit comprising objective function of decision-making and decision constraint
The rural power grids taiwan area drop damage Optimization Decision Models of condition;
(1) object function of decision;
Reach comprehensive optimal to make rural power grids taiwan area Construction and operation cost and loss of electricity reduce acquired benefit, in the full longevity
Under the support for ordering life cycle costing theory, establish the overall cost that object function includes rural power grids taiwan area in the operation time limit, drop damage at
Sheet, operation expense, failure cost and cost recovery is removed, sale of electricity expense caused by loss of electricity;
(2) constraints of decision:
The constraints of rural power grids taiwan area drop damage Optimization Decision Models includes six contents:Power supply reliability constraint, Yong Huping
Mutex relation constrains between equal power off time constraint, voltage deviation constraint, line transmission constraint, investment restriction and project;
1. power supply reliability constrains:Power supply reliability should be not less than preset prescribed limits;
2. average power off time of user constrains:Average power off time of user should be no more than preset prescribed limits;
3. voltage deviation constrains:The absolute value of voltage deviation should be no more than the 7% of normal voltage;
4. line transmission constrains:The actual transmissions capacity of circuit should be no more than its maximum transfer capacity, general transmission electricity
Flow table is shown;
5. investing restriction:Investment restriction refers to taking the investment of the drop damage scheme of a variety of reducing loss measure combinations
Should be within preset investment limitation, drop damage formula investing is primarily directed to drop damage cost at this time;
(3) Evaluation of Investment-Benefit index:
It, should first really for rural power grids wastage reducing and energy saving transformation project according to the basic step of technology and economy comparison and evaluation
Vertical relevant art index, then economic indicator optimum selecting is pressed from the scheme of technical indicator qualification;The particular technique of damage scheme drops
Index includes electric energy loss rate, rate of qualified voltage, power factor and service life;
According to《Rural cadastration Guidance Rule》In " economic technology assessment " regulation drop damage scheme economic indicator packet
Include investment payback time, net present value (NPV), net annual value, net present value rate and the internal rate of return (IRR);
(4) decision of drop damage prioritization scheme:
The decision purpose of drop damage scheme is excellent to selecting the drop damage scheme collection of possibility to carry out comprising various reducing loss measures
Choosing determines the optimal reducing loss measure combination of comprehensive benefit and drop damage embodiment.
The step 6 includes:
Step 61:If the reactive compensation capacity calculated in step 4 is c1, the reactive compensation capacity calculated in step 5 is c2,
If the difference that the two calculates is less than the nargin ε of setting, it is c to take reactive compensation capacity1If the difference of the two is more than setting
Nargin then establishes decomposition-coordination approach and solves reactive compensation capacity;An original object function is resolved into two to connect each other
Specific item scalar functions;Two sub- object function interactions, common optimization find out final reactive compensation capacity;
Step 62:Initially set up mathematical model:
min C(x)+D(x)
s.t A(x)≥B
E(x)+F(x)≥g
R (x)+T (x)=h
X is the reactive compensation capacity size that is configured in above formula, and y is power grid actual operation parameters under existing configuration, C (x)+
D (x) is the expense for configuring reactive power source, and A (x) >=B is investment and recovery (planned allocation constrains, configuration capacity constraint), E (x)+F
(x) >=g is voltage inequality constraints, and R (x)+T (x)=h constrains for trend equation;
As long as reactive compensation capacity investment variable x is decided, so that it may be determined by the optimizing of system operation mode practical
Operating parameter y;That is as former problem can be regarded to the unification of two optimization process, one is to determine idle configuration capacity position
With the investment optimization process of size, primal problem is referred to as invested;Two are to determine and so that operating cost is reached most under existing configuration
Small running optimizatin process, referred to as runs subproblem;
As given any determination capacity x*, the optimal value D of corresponding operation subproblem is represented by W (x), then corresponding to invest
Primal problem and operation subproblem decomposable asymmetric choice net are as follows:
(1) primal problem is invested:
Min C(x)+W(x)
s.t A(x)≥B
In above formula, C (x)+W (x) is the total cost of investment, and wherein W (x) is a knob for contacting primal problem with subproblem
Band, A (x) >=B are the conditions of investment and recovery;
(2) subproblem is run:
W (x)=minD (y)
s.t F(x)≥g-E(x*)
T (y)=h-R (x*)
In above formula, minD (y) is optimal the spent expense of operation, F (x) >=g-E (x*) it is voltage inequality constraints, T
(y)=h-R (x*) it is that trend equation constrains;
Obviously, the primal problem of investment and the subproblem of operation are closely connected with W (x) here, by running son
The solution of problem proposes the amendment to W (x), forms linear restriction new after correcting, and returns again to investment primal problem, alternately solves
Primal problem and subproblem simultaneously obtain the optimal solution of reactive-load compensation amount.
Compared with prior art, advantages of the present invention and advantageous effect are:
The present invention is using dividing gradient search procedure to obtain accurate physical parameter in the topological structure of fully analysis network
Accurate line loss calculation is prepared;It is used according to the different characteristics of rural power grids typical case's taiwan area and radius of electricity supply, transformer model capacity
The data cleansing and excavation of the physical parameters such as the time limit, power load distributing situation.Rural power grids taiwan area is calculated using the improved trend method that is pushed forward
Loss so that the line loss result of calculating is more accurate, on the basis of line loss carry out reactive-load compensation add voltage, power because
Number, compensation capacity, the constraints for compensating equipment, reactive-load compensation can accomplish idle in-situ balancing, carry out reducing loss measure
It can accomplish to adjust voltage simultaneously, make every effort to the best compensation effect taken.It is idle in progress in order to obtain preferably drop damage effect
Increasing on the basis of compensation improves three-phase imbalance, replaces transformer capacity, replaces the reducing loss measures such as circuit section product.
In the following with reference to the drawings and specific embodiments, further description and explanation are subject to the present invention, but are not implemented by this
Example is limited.
Description of the drawings
Fig. 1 is the voltage and power diagram of the corresponding circuit of the method for the present invention;
Fig. 2 is the corresponding taiwan area structure chart of the method for the present invention;
Fig. 3 is the composition structure chart of the corresponding drop damage scheme of the method for the present invention;
Fig. 4 is the reducing loss measure optimization block diagram of the rural power grids taiwan area of the corresponding returns of investment rate of the method for the present invention.
Specific implementation mode
To make the object, technical solutions and advantages of the present invention clearer, below in conjunction with attached in the embodiment of the present invention
Figure, clear, complete description is carried out to technical scheme of the present invention.Obviously, described embodiment is that a part of the invention is real
Example is applied, instead of all the embodiments.Based on the embodiments of the present invention, those of ordinary skill in the art are not making creation
Property labour under the premise of the every other embodiment that is obtained, shall fall within the protection scope of the present invention.
A kind of drops such as reactive-load compensation based on rural power grids typical case's taiwan area line loss calculation of the invention damage optimization method, including walk as follows
Suddenly:
Step 1:The topological structure for analyzing network is accurate using dividing gradient search procedure to obtain accurate physical parameter
Line loss calculation is prepared;According to the different characteristics of rural power grids typical case's taiwan area and radius of electricity supply, transformer model capacity service life, bear
The data cleansing and excavation of the physical parameters such as lotus distribution situation carry out typicalness classification to rural power grids taiwan area.
According to the network structure of distribution network system and switch operating status, classify to the load of power distribution network, according to not
With load and whether there is or not measuring values to carry out a point gradient search, to obtain the parameter of each branch of system and the voltage of node, active
Power and reactive power.
The active and reactive power and voltage value of main transformer, distribution transforming in first gradient search distribution feeder.
Transformer for failing to obtain real-time measurement carries out the second gradient search, and complete and precision satisfaction is wanted in order to obtain
The operating parameter asked, the use state estimation technique carry out line loss calculation.
In the quantity for having searched for transformer, 3rd gradient search is carried out, is searched on user and overhead transmission line, cable run
Each node forms node matrix equation.
Step 2:Using the improved loss for being pushed forward trend method and calculating rural power grids taiwan area;
The improved trend method that is pushed forward refers to taking being averaged for acquired results by measurement to node voltage value and calculating
For value so that the voltage value of each node is more accurate, the line loss value of calculating is truer.
The step 2 includes:By obtaining the measured data of each end load node, each section is calculated by prospective method
The voltage value and performance number of the corresponding last node of point, if as soon as upper node uses without measured value and calculate forward calculation institute
The voltage value obtained takes the average value of measured value and calculated value if any measured value;If a upper node corresponds to two load bus,
If this node takes two loads to be pushed forward the average value of voltage without measured value, the average value of three's value is otherwise taken, and so on;Directly
To root node.The result of line loss is equal to the total losses that transformer calculates and adds the total losses calculated on circuit.
Aloss=Tloss+Lloss
A in above formulalossFor all line loss, TlossFor the loss of transformer, LlossFor the loss of transmission line of electricity.
Step 21:Known to line impedance parameter, as long as measuring active and reactive, the voltage of load bus
Virtual value, so that it may accurately to find out very much loss between load bus and its superior node on circuit and superior node
Voltage, and calculating line loss and node voltage all avoid plural multiplication and division.Circuit l between node 2 to 3 as shown in Figure 12-3
It is Z per phase of impedance2-3=R2-3+jX2-3The three-phase output power for measuring node 3 simultaneously is P3+jQ3, line voltage virtual value is U3。
Circuit l2-3On active power loss be:
In above formula, Δ P2-3Refer to circuit l2-3On active loss;P3It is the active power of node 3, Q3It is node 3
Reactive power;U3It is the voltage of node 3, R2-3It is circuit l2-3On resistance value.
Reactive power loss is:
In above formula, Δ Q2-3Refer to circuit l2-3On active loss;P3It is the active power of node 3, Q3It is node 3
Reactive power;U3It is the voltage of node 3, R2-3It is circuit l2-3On resistance value.
The output power of node 2 is:
S2=(P3+ΔP2-3)+j(Q3+ΔQ2-3)
In above formula, S2Refer to the output power of node 2;P3It is the active power of node 3, Q3It is the idle work(of node 3
Rate;ΔP2-3、ΔQ2-3It is circuit l respectively2-3The active power and reactive power of upper loss.
The voltage effective value U of node 22For:
In above formula, U2Refer to the RMS voltage Δ U of node 22-3、δU2-3It is circuit l2-3The vertical component of voltage landing and
Horizontal component.
Step 22:Assuming that a certain taiwan area, structure are as shown in Figure 2.The taiwan area shares 11 nodes, it is known that each node it
Between resistance, reactance, and by data collecting system have been obtained for the active power of each load bus, reactive power,
Voltage effective value.
Steps are as follows for the specific implementation of algorithm:
(1) node-classification is carried out first and is numbered, and power supply node, load bus are distinguished.
By taking Fig. 2 as an example, power supply node is 0 node, 2,3,8,9,10 node of load bus
(2) by obtaining the measured data of each end load node, by prospective method calculate each node it is corresponding on
The voltage value and performance number of a node, if as soon as upper node uses the voltage value calculated obtained by forward calculation without measured value,
If any measured value, the average value of measured value and calculated value is taken;If a upper node corresponds to two load bus, if this node without
Measured value, takes two loads to be pushed forward the average value of voltage, otherwise takes the average value of three's value, and so on;Until power supply node.
Circuit between node 8 and 6 as shown in Figure 2 is l8-6Each impedance between circuit is Z8-6=R8-6+jX8-6It measures
8 node of load three output powers be P8+jQ8Then circuit l8-6Between active power loss:
In above formula, Δ P8-6It is circuit l8-6Active power loss, P8It is the active power of node 8, Q8Node 8 it is idle
Power, U8It is the voltage effective value of node 8, R8-6It is circuit l8-6Between impedance.
Reactive power loss is:
In above formula, Δ Q8-6It is circuit l8-6Reactive power loss, P8It is the active power of node 8, Q8Node 8 it is idle
Power, U8It is the voltage effective value of node 8, X8-6It is circuit l8-6Between reactance.
The output power of node 6 is:
S6=P6+jQ6=(P8+ΔP8-6)+j(Q8+ΔQ8-6)
In above formula, S6Refer to the output power of node 6;P6It is the active power of node 6, Q6It is the idle work(of node 6
Rate;ΔP8-6、ΔQ8-6It is circuit l respectively8-6The active power and reactive power of upper loss.P8、Q8It is the wattful power of node 8 respectively
Rate and reactive power.
The voltage value of node 6 is:
In above formula, U8For the voltage effective value of node 8, Δ U8-6For circuit l8-6The vertical component of voltage landing, δ U8-6For line
Road l8-6The horizontal component of voltage landing.
Wherein:
Circuit l can similarly be passed through9-6The voltage and output power for calculating node 6, if the voltage of node 6 can be measured
Value obtains the voltage value and performance number of node 6 with output power by calculating the average value of three, if the node 6 that do not survey
The voltage value of voltage value and performance number, node 6 takes the average value of calculating, performance number also to take the average value of calculating.
(3) circuit having had been calculated is marked, continue search for the circuit of non-marked until all circuits all
It calculates.
(4) it calculates the loss of all common transformers, when calculating, the node load of actual measurement, voltage value is utilized to calculate transformation
The copper loss and iron loss of device.So, transformer loss calculates the influence for just having fully taken into account load, voltage.
(5) total line loss:
Aloss=Tloss+Lloss
Wherein AlossFor all line loss, TlossFor the loss of transformer, LlossFor the loss of transmission line of electricity.
Step 3:It establishes the object function of idle work optimization and adds voltage, power factor, compensation capacity, the pact for compensating equipment
Beam condition;
Step 31:The object function:
min(P1-P2)
Wherein P1 is investment cost, and P2 is drop damage income, and object function is the minimum value of two differences.
Wherein:
In above formula, n is that node compensates equipment sum;NiFor the packet count of the compensation capacitor of the i-th complete equipment;Ci is compensation
The capacity of the corresponding single group capacitor of equipment;Ma is the unit price of compensation capacity;Mb is the single group price for compensating equipment;Mc is
Expense input by equipment is installed.
Wherein:
P2=(f1-f2)·c·T
In above formula, f1For the power loss before compensation, f2For the power loss after compensation, c is corresponding electricity price, T operations
Time.
Specifically:To minimize the active loss of low-voltage network as target, low pressure:Distribution network var compensation is distributed
Optimized model is as follows:
min f(x)
S.t h (x)=0,
In above formula:Object function f (x) is the difference of investment cost and drop damage income;H (x) is that the power of each node is flat
Weigh equation, and x is the decision variable and state variable of system, and decision variable is the parallel reactive compensation capacity Q of reactive compensation point, shape
State variable includes the voltage magnitude V and phase angle δ of each node.Inequality constraints l (x) include the bound of each node voltage amplitude about
The bound of beam and each reactive compensation point parallel reactive compensation capacity constrains and the bound of power factor is divided with compensation equipment
The bound of group number.
The difference of object function f (x) investment costs and drop damage income, expression formula are as follows:
F (x)=P1-P2
Step 32:The constraints.
Equality constraint h (x) expression formulas are:
It can be obtained by Load flow calculation equation:
In above formula, PiAnd QiIt is the active power and reactive power of taiwan area injection respectively, it is negative, Q to flow into just, to flow outci
For the reactive power of compensation, ViFor the voltage magnitude of node i, δijFor the phase difference of voltage between node i and node j, GijAnd Bij
It is the real and imaginary parts of network node admittance matrix the i-th row j column elements respectively, n is the sum of taiwan area node.
Inequality constraints l (x) expression formulas are:
V in above formulaimin、VimaxIt is the minimum voltage and maximum voltage that node allows, Qcimin、QcimaxIt is the minimum that node is permitted
It is idle to be compensated with maximum reactive-load compensation, ρ imin、ρimaxIt is the minimum and maximum power factor that node measured place allows,It is the minimum packets number and largest packet number for compensating equipment switching.
Step 4:It determines the compensation position of the reactive-load compensation of different type typical case's taiwan area and each compensates what equipment should compensate
Capacity;
Step 41:The compensation rate that taiwan area all nodes in Typical Rural respectively need is calculated first, calculates taiwan area needs
The progress reactive-load compensation of each node is not met reality, using Sensitivity Analysis Method and artificially by total reactive compensation capacity
Factor method selects the node that the node that load power is big, important is compensated as needs.Every in population is mended using genetic algorithm
It repays equipment and sets 2 genes on chromosome, indicate the packet count and every group of compensation capacity of compensation equipment respectively, calculate egress
Optimal non-workcompensation capacity.
Step 42:When calculating the idle configuration optimization of different type taiwan area on the basis of step 41, three are defined first
Electric characteristic index circuit distribution transforming Rate of average load α, natural power factor cos φ, radius of electricity supply L, reactive-load compensation rate β, according to
《Rural power planning and designing directive/guide》With the specific electric characteristic index system in local rural area, typical rural area electrical specifications are set
As ground state value.
In above formula, β0For the optimal reactive compensation rate of Typical Rural taiwan area, WnHold for the idle configuration that the taiwan area requires supplementation with
Amount, SNFor the rated capacity of taiwan area.
When calculating the reactive compensation capacity of random taiwan area by changing each electric characteristic index respectively, OPTIMAL REACTIVE POWER is recorded
The situation of change of cancellation ratio indicates that the variation degree of reactive-load compensation rate, the size of sensitivity are equal to respectively with the size of sensitivity
The variation of the corresponding reactive-load compensation rate of electric characteristic variation exists in the reactive capability configuration for carrying out random taiwan area, with each electrical
The difference of the actual value of characteristic value electric characteristic corresponding with typical taiwan area carries out phase as decision using the size of sensitivity
The weighted calculation answered.
β in above formula0For the reactive-load compensation rate of typical taiwan area, Δ miIt is right for random taiwan area and typical taiwan area electrical characteristics
The difference answered, λiFor the sensitivity corresponding to electrical characteristics.
Step 5:Establish the rural power grids taiwan area of more reducing loss measure combinations.
Step 51:Determine that the reducing loss measure of rural power grids taiwan area is denoted as X firsti, wherein i=1,2 ..., 4;X1~X4It is right with this
Reactive-load compensation, improvement three-phase imbalance, replacement transformer capacity, replacement circuit section is answered to accumulate these four rural power grids taiwan areas drop damage and arrange
Apply type.On this basis, scheme is damaged using the possible embodiment of certain a kind of reducing loss measure as the drop in rural power grids taiwan area drop damage,
AsIndicate the jth kind possible embodiment of i-th kind of reducing loss measure type.X1The corresponding drop damage scheme of reducing loss measure has change
Pressure station centralized compensation, user's decentralised control, Optimal Reactive Power compensation device.X1The corresponding drop damage scheme of reducing loss measure has change client
Wiring phase sequence, three-phase imbalance pressure regulator, auto-transformer.X3The corresponding drop damage scheme of reducing loss measure has volume reduction transformation method, becomes
Depressor way of economic operation, energy saving transformer.X4The corresponding drop damage scheme of reducing loss measure has parallel gathering lead wire in low, increases conducting wire
Sectional area, line insulation transformation.
For reducing loss measure Xi, since the presence of analysis of Influential Factors is lost, from reducing loss measure weighing factor analysis and assessment
As a result the larger reducing loss measure type of weight is selected in, and alternately reducing loss measure is implemented present situation rural power grids taiwan area, thus,
To reducing loss measure XiWhether selected available following formula describes.
For selected reducing loss measure, it is further contemplated that its possible reducing loss measure implements project, it is possible to
Damage project dropsIt, can be existing with difference according to artificial experience for the description of the implementation principle of rural power grids taiwan area reducing loss measure
Optimization algorithm proposes a variety of alternative implementation projects to every a kind of reducing loss measure type, and following formula is to i-th kind of reducing loss measure
In the whether selected description made of jth kind drop damage project.
The composition structure chart of the corresponding drop damage scheme of the method for the present invention is as shown in Figure 3.
Step 52:In order to ensure that drop damage reaches better effect drop damage alternative is formed by following steps.
Drop damage alternative is formed by drop damage projects combo, since the presence of analysis of Influential Factors is lost in rural power grids taiwan area,
The major influence factors of the influence rural power grids taiwan area technological attrition in present situation rural power grids taiwan area and its corresponding can be primarily determined
Reducing loss measure type, thus, can determine the situation of choosing of reducing loss measure generating the when of dropping damage alternative, it is assumed that warp
Selected reducing loss measure is respectively X after crossing loss analysisa、Xb、Xc, then the composition structure such as 3 of damage candidate item drops in this case
Shown in figure:
In the case of the implementation project determination of every a kind of reducing loss measure type, a drop damage candidate item is drop damage project
The combination of selected situation, while also meeting drop damage projectBetween alternative, be shown below:
In above formula:fiRefer to i-th of drop damage candidate item;G refers to the alternative between drop damage project, i.e., same kind of drop damage is arranged
Implementation of the different implementation projects applied in rural power grids taiwan area is impossible simultaneous constraint.Mutex relation constraint master between project
If being embodied in drop damages decision process.
The product process of rural power grids taiwan area drop damage candidate item is as follows:
(1) the loss shadow that key important function is played to rural power grids taiwan area technical loss is established by rural power grids taiwan area loss analysis
The factor of sound;
(2) it is based on sensitivity analysis, establishes the reducing loss measure X in rural power grids taiwan area Optimal Decision-makingiChoose situation;
(3) to each selected reducing loss measure Xi, obtained for choosing by the research of each measure scheme weight analysis
The drop damage project selected
(4) according to the initial rural power grids taiwan area drop damage candidate item f of the various combination producings of drop damage projecti;
Step 53:In damage decision model drops in rural power grids taiwan area, it is present situation rural power grids taiwan area band that should fully embody reducing loss measure
The comprehensive benefit of the drop damage electricity come is established by analyzing Evaluation of Investment-Benefit comprising objective function of decision-making and decision constraint
The rural power grids taiwan area drop damage Optimization Decision Models of condition, the reducing loss measure optimization frame of rural power grids taiwan area are as shown in Figure 4.
(1) object function of decision:
Reach comprehensive optimal to make rural power grids taiwan area Construction and operation cost and loss of electricity reduce acquired benefit, in the full longevity
Under the support for ordering life cycle costing theory, establish overall cost that object function includes rural power grids taiwan area in the operation time limit (drop damage at
Sheet, operation expense, failure cost and remove cost recovery), sale of electricity expense caused by loss of electricity.
(2) constraints of decision:
The constraints of rural power grids taiwan area drop damage Optimization Decision Models includes six contents:Power supply reliability constraint, Yong Huping
Mutex relation constrains between equal power off time constraint, voltage deviation constraint, line transmission constraint, investment restriction and project.
1. power supply reliability constrains:Power supply reliability should be not less than preset prescribed limits.
2. average power off time of user constrains:Average power off time of user should be no more than preset prescribed limits.
3. voltage deviation constrains:The absolute value of voltage deviation should be no more than the 7% of normal voltage.
4. line transmission constrains:The actual transmissions capacity of circuit should be no more than its maximum transfer capacity, general transmission electricity
Flow table is shown.
5. investing restriction:Investment restriction refers to taking the investment of the drop damage scheme of a variety of reducing loss measure combinations
Should be within preset investment limitation, drop damage formula investing is primarily directed to drop damage cost at this time.
(3) Evaluation of Investment-Benefit index:
It, should first really for rural power grids wastage reducing and energy saving transformation project according to the basic step of technology and economy comparison and evaluation
Vertical relevant art index, then economic indicator optimum selecting is pressed from the scheme of technical indicator qualification.The particular technique of damage scheme drops
Index includes electric energy loss rate, rate of qualified voltage, power factor and service life etc..
According to《Rural cadastration Guidance Rule》In " economic technology assessment " regulation drop damage scheme economic indicator packet
Include investment payback time, net present value (NPV), net annual value, net present value rate and internal rate of return (IRR) etc..
(4) decision of drop damage prioritization scheme:
The decision purpose of drop damage scheme is excellent to selecting the drop damage scheme collection of possibility to carry out comprising various reducing loss measures
Choosing determines the optimal reducing loss measure combination of comprehensive benefit and drop damage embodiment.
Step 6:Establish the reactive compensation capacity of decomposition-coordination approach.
Step 61:If the reactive compensation capacity calculated in step 4 is c1, the reactive compensation capacity calculated in step 5 is c2,
If the difference that the two calculates is less than the nargin ε of setting, it is c to take reactive compensation capacity1If the difference of the two is more than setting
Nargin then establishes decomposition-coordination approach and solves reactive compensation capacity.An original object function is resolved into two to connect each other
Specific item scalar functions.Two sub- object function interactions, common optimization find out final reactive compensation capacity.
Step 62:Initially set up mathematical model:
min C(x)+D(x)
s.t A(x)≥B
E(x)+F(x)≥g
R (x)+T (x)=h
X is the reactive compensation capacity size that is configured in above formula, and y is power grid actual operation parameters under existing configuration, C (x)+
D (x) is the expense for configuring reactive power source, and A (x) >=B is investment and recovery (planned allocation constrains, configuration capacity constraint), E (x)+F
(x) >=g is voltage inequality constraints, and R (x)+T (x)=h constrains for trend equation.
As long as reactive compensation capacity investment variable x is decided, so that it may be determined by the optimizing of system operation mode practical
Operating parameter y.That is as former problem can be regarded to the unification of two optimization process, one is to determine idle configuration capacity position
With the investment optimization process of size, primal problem is referred to as invested;Two are to determine and so that operating cost is reached most under existing configuration
Small running optimizatin process, referred to as runs subproblem.
As given any determination capacity x*, the optimal value D of corresponding operation subproblem is represented by W (x), then corresponding to invest
Primal problem and operation subproblem decomposable asymmetric choice net are as follows:
(1) primal problem is invested:
Min C(x)+W(x)
s.t A(x)≥B
In above formula, C (x)+W (x) is the total cost of investment, and wherein W (x) is a knob for contacting primal problem with subproblem
Band, A (x) >=B are the conditions of investment and recovery.
(2) subproblem is run:
W (x)=minD (y)
s.t F(x)≥g-E(x*)
T (y)=h-R (x*)
In above formula, minD (y) is optimal the spent expense of operation, F (x) >=g-E (x*) it is voltage inequality constraints, T
(y)=h-R (x*) it is that trend equation constrains.
Obviously, the primal problem of investment and the subproblem of operation are closely connected with W (x) here, by running son
The solution of problem proposes the amendment to W (x), forms linear restriction new after correcting, returns again to investment primal problem.Alternately solve
Primal problem and subproblem simultaneously obtain the optimal solution of reactive-load compensation amount.