Power distribution network wide area idle work optimization operation method
Technical field
The invention belongs to the reactive power compensation technology field of system for distribution network of power, be specifically related to the Optimum Operating Method of system for distribution network of power reactive power compensation.
Background technology
Power distribution network is direct correlation electric power system and user's an important tie, and area coverage is big, and number of devices is many, is the important component part of electric power system.Improve the quality of voltage of power distribution network, the running wastage of reduction power distribution network, for improving power supply quality, guaranteeing that the high-quality economical operation of electric power system is significant.Reactive power compensation is to improve one of effective measures of power distribution network quality of voltage, reduction power distribution network running wastage.
The idle work optimization operation of power transmission network is adopted optimization method calculating optimum reactive compensation capacity usually based on the electric power system tide equation.This method needs complete electric power system tide information, and at present, the automaticity of China's power distribution network is not high, is difficult to obtain the real-time status information of power distribution network, so this method and be not suitable for the idle work optimization operation of power distribution network.
The reactive power compensation of existing power distribution network is according to the 8th regulation in " the State Grid Corporation of China's power system reactive power compensation configuring technical principle " of Q/GDW " State Grid Corporation of China's company standard " of issue in 2008: the reactive power compensation of power distribution network is main with power distribution network step down side centralized compensation; It is constraints that the assembling of the capacitor of distribution transformer is established with voltage, according to the divide into groups control device of automatic switching of reactive power (reactive current).The operation/cutting method of existing distribution network var compensation device is each reactive power compensator attachment point (being the distribution transformer low-pressure side); Principle by idle in-situ balancing; Balance according to distribution transformer load or burden without work; With the voltage request that satisfies the reactive power compensator attachment point is condition, automatically the switching capacity of control compensation capacitor group.The major defect of existing distribution network var compensation device operation/cutting method is: can only be according to the compensation capacity of the load or burden without work balance automatic switching compensation condenser of each distribution transformer; Only realized the reactive balance of each distribution transformer and following electrical network regional area, can not control the switching of reactive compensation capacity from the overall condition of power distribution network load or burden without work.Because the Changing Pattern of each distribution transformer load (comprising load or burden without work) and inequality in the power distribution network; Existing distribution network var compensation device operation/cutting method is insufficient to utilizing of compensation condenser; Be difficult to guarantee the reactive balance of power distribution network integral body, thereby influence the high-quality economical operation of power distribution network.
Summary of the invention
The objective of the invention is deficiency to existing distribution network var compensation control method; A kind of power distribution network wide area idle work optimization operation method is provided; Voltage with the reactive power compensator attachment point is constraints, takes into full account the voltage and the idle operation conditions of whole power distribution network, calculates the switching capacity of each reactive power compensator automatically and controls the switching of each reactive power compensator; Can effectively improve the high and power distribution network quality of voltage of utilance of compensation equipment, reduce the power distribution network running wastage.
The technical scheme that realizes the object of the invention is: a kind of power distribution network wide area idle work optimization method, utilize computer, and through program, at first import the basic data of power distribution network and each reactive power compensator thereof; Calculate the compensation range and the priority level thereof of each reactive power compensator again; According to compensation range and the root node of power distribution network and the real-time measurement information of each compensation point of each reactive power compensator, calculate the switching capacity of each reactive power compensator then, each the capacitor group of control of holding concurrently is carried out correct switching.The concrete steps of said method are following:
(1) input basic data
At first import the basic data of power distribution network, the basic data of each reactive power compensator, historical switching data and power distribution network root node and each reactive power compensation point metric data of each reactive power compensator.Wherein, the basic data of power distribution network comprises the basic data of each node, i.e. node serial number, node place electric pressure, the node voltage upper limit, node lower limit circuit; The basic data of each circuit, promptly the first end-node numbering of circuit, resistance (
R l ), reactance (
X l ), susceptance (
B l ), rated voltage (
U Bl ); The basic data of each distribution transformer, i.e. distribution transformer institute on-load numbering, place, both sides node serial number, resistance
R T ), reactance
X T ), electricity lead (
G m ), susceptance (
B m ), no-load voltage ratio (
k T ), rated capacity (
S NT ), the high-pressure side rated voltage (
U Bl1 ), the low-pressure side rated voltage
U Bl2 ).The basic data of each reactive power compensator comprise compensation condenser group place node number, single group compensation capacity (
Q Ref ), the group number (
n), total compensation capacity (
Q c ), maximum compensation capacity (
Q Max ), minimum compensation capacity (
Q Min ).The historical switching data of each reactive power compensator comprise capacity data and the voltage data of power distribution network root node and each reactive power compensation point before and after each reactive power compensator capacitor switching of previous moment, the metric data of power distribution network root node and each reactive power compensation point comprise the actual measurement voltage of current time (
U), active power (
P), reactive power (
Q).
Calculate the control priority level and the control range of each reactive power source
After (1) completion, according to the annexation of each node of power distribution network, utilize the BFS algorithm to confirm each reactive power source control priority level of (comprising power distribution network root node and each reactive power compensator), process is following:
To root node in the power distribution network and reactive power compensator, begin from root node, and root node is designated as 0 grade; With the 0th grade adjacent confirm that not the reactive power compensator of grade is the 1st grade; With the 1st grade adjacent confirm that not the reactive power compensator of grade is the 2nd grade, so go down, utilize BFS reactive power compensator place node; All searched up to all reactive power compensators, and last reactive power compensator is designated as
mLevel.
mLevel is the highest control priority of reactive power source, 0 grade of expression power distribution network root node.
The
kIndividual load is by the reactive power source point
iThe equivalent loss that the reactive power flow movable property is given birth to:
In the formula:
L Ki Be
kIndividual load is by the reactive power source point
iThe equivalent loss of the network that the reactive power flow movable property is given birth to;
S k Be the th
kThe capacity of individual load place distribution transformer;
R K-i For connecting
kIndividual load place node and reactive power source point
iAll branch roads (comprising circuit and transformer branch road) resistance sum.
With
kIndividual load inserts the minimum reactive power source point of equivalent loss, and then the
kIndividual load place distribution transforming will belong to the control range of this reactive power source point.Calculate the equivalent loss corresponding of all loads, all insert corresponding reactive power source point up to all loads with reactive power source.
(3) calculate the switching Capacity Ratio of each reactive power compensator
(2) step was calculated each reactive power source, i.e. the switching Capacity Ratio of each reactive power compensator except that the power distribution network root node after accomplishing.
The
jThe switching Capacity Ratio of individual reactive power compensator
In the formula:
α j Be
jThe switching Capacity Ratio of individual reactive power compensator;
TSet for distribution transforming in the whole power distribution network;
S t Be
tThe capacity of individual distribution transforming;
T j Be
jThe set of distribution transforming in the individual reactive power compensator control range;
S Tj Be
t j The capacity of individual distribution transforming.
According to the result of calculation in (2) step, calculate the switching Capacity Ratio of each compensation arrangement in conjunction with each reactive power compensator control range.
The idle sensitivity matrix of calculating voltage
(3) accomplish after, calculating voltage monitoring point and reactive power compensation point and between the voltage power-less sensitivity matrix.
Confirm that at first voltage monitoring point is power distribution network root node and reactive power compensator place node.
Confirm then voltage monitoring point and reactive power compensation point and between the voltage power-less sensitivity matrix
S UQ :
In the formula:
S UQ Be one (
m+ 1) *
mThe matrix of dimension; d
U i / d
Q Cj Be reactive power compensator
jReactive compensation capacity and voltage monitoring point
iBetween voltage power-less sensitivity.
Calculate then
S UQ In the value of each element, computing formula is:
In the formula: d
U i / d
Q Cj Be reactive power compensator
jReactive compensation capacity and voltage monitoring point
iBetween voltage power-less sensitivity; Δ
Q Cj It is the compensation node
jThe variable quantity of compensation capacity; Δ
U i Be that compensation capacity changes the back voltage monitoring point
iThe variable quantity of voltage (this sensitivity specifically is to be tried to achieve by former and later two voltage magnitude variation and condenser capacity variations of voltage monitoring point constantly of capacitor switching by the voltage history data of voltage monitoring point and the historical switching data computation of each compensation arrangement).
Ask and obtain
S UQ In after the value of each element, form according to formula (3)
S UQ
Calculate the switching capacity of each reactive power compensator
After (4) completion, at first calculate the idle total load of power distribution network
Q L , computing formula is:
In the formula,
Q L Be the idle total load of power distribution network;
Q 1It is current time power distribution network root node reactive power measuring value;
Be
jA capacitor compensation capacity constantly on the individual reactive power compensator is write as vector form and is then done
According to the switching Capacity Ratio of each low-voltage reactive compensator capable, calculate the reactive requirement in each reactive power compensator control range again, the
jThe computing formula of reactive requirement is in the individual reactive power compensator control range:
In the formula,
Q ' Cj-Need Be
jReactive requirement in the individual reactive power compensator control range;
α j Be
jThe switching Capacity Ratio of individual reactive power compensator;
Q L Be the idle total load of power distribution network.
When
jThe compensation capacity of individual reactive power compensator has residue; And near the reactive power compensator compensation capacity this reactive power compensator is when not enough; The residue compensation capacity is preferentially compensated the control priority level than the reactive requirement in the control range of own high reactive power compensator, and compensation control priority level is than the reactive requirement in the control range of own low reactive power compensator then.
Recomputate the reactive requirement of reactive power compensator, computing formula is:
In the formula,
Q Cj-Need Be after recomputating
jReactive requirement in the individual reactive power compensator control range;
α j Be
jThe switching Capacity Ratio of individual reactive power compensator;
Q L Be the idle total load of power distribution network; Δ
Q J-Need Be
jNear the reactive capability of reactive power compensator individual reactive power compensator is supported.
Consider constraints then, calculate the switching capacity of each reactive power compensator.Difference sum minimum with the reactive load demand in the actual switching capacity of each reactive power compensator capacitor and its control range is a target function; Be restricted to constraints with voltage monitoring node voltage bound, the restriction of reactive power compensation node compensation rate, reactive power compensation node transformer capacity, can power distribution network wide area idle work optimization model be formulated as follows:
Target function:
j=1,2 ...,
m(9)
Constraints:
The qualified constraint of voltage:
(10)
The reactive power compensator constraint of exerting oneself:
(11)
Reactive power compensation point distribution transforming varying capacity constraint:
(12)
In the formula:
Q Cj Expression the
jThe capacitor reactive compensation capacity of individual reactive power compensation node, when compensation arrangement adopts the capacitor group as reactive power source,
Q Cj Be discrete magnitude, write as vector form and do
Q Cj-Need Be
jThe reactive requirement of individual reactive power compensator;
UBe the voltage vector of each voltage monitoring point,
U=[
U 0,
U 1,,
U i,
,
U m ]
T,
U Max ,
U Min The voltage magnitude bound vector of representing voltage monitoring point respectively; Q
Max , Q
Min Represent each reactive power compensator capacity bound vector respectively;
P j Expression the
jThe active power of individual reactive power compensator place load,
Q j Expression the
jThe active power of individual reactive power compensator place load,
S NTj Expression the
jThe individual transformer capacity upper limit.
The voltage vector of next each voltage monitoring point of the moment behind the calculating capacitor switching in the formula (10)
UFormula do
In the formula:
U=[
U 0,
U 1,,
U i,
,
U m ]
T U 0 Voltage vector for each voltage monitoring point of current time;
S UQ Be voltage monitoring point and reactive power compensation point and between the voltage power-less sensitivity matrix;
Q C It is the capacitor compensation capacity vector of reactive power compensator;
It is a moment compensation capacity vector on the reactive power compensator.
Obtain with enumerative technique solution formula (9) ~ (12)
Q C , promptly try to achieve the vector of each compensation arrangement switching capacity
(
6) export the switching capacity of each reactive power compensator
After the completion of (5) step, the vector of the reactive power compensator switching capacity that (5) step of output earlier solves
Δ Q C , be the switching amount of each reactive power compensator capacitance.Each reactive power compensator of switching gets capacitor respectively again, and record switching front and back capacitor compensation capacity and voltage monitoring point voltage magnitude situation of change, in order to the input data as next moment.
After the present invention adopts technique scheme, mainly contain following effect:
1, compare power distribution network low-pressure reactive compensation under the control mode on the spot, the inventive method makes the compensatory zone partition balancing of the reactive power of power distribution network according to reactive power compensator, has taken into full account the whole voltage power-less operation conditions of power distribution network; Thereby reduce flowing of reactive power in the power distribution network; Reduce the loss of power distribution network, remarkable economic efficiency has been arranged, and made full use of the high-quality secondary energy sources; Play a role in energy saving, significant to constructing economical society.
2, compare power distribution network low-pressure reactive compensation under the control mode on the spot; The inventive method makes the reactive power of power distribution network according to the partition balancing of reactive power compensator compensatory zone; Reduced the voltage loss between power distribution network root node and the load bus; Improved the quality of voltage of power distribution network, guaranteed the normal operation of power consumption equipment under rated condition, the significant social economic benefit has been arranged.
3, compare power distribution network low-pressure reactive compensation under the control mode on the spot; The inventive method adopts reactive power compensator to carry out the operational mode of switching according to compensatory zone; Overcome in the existing method building-out capacitor has been utilized inadequate shortcoming; Greatly improve the utilance of compensation equipment, further improved the economic benefit of power distribution network.
The present invention can be widely used in the voltage and reactive power optimization operation and control of system for distribution network of power.
Description of drawings
Fig. 1 is the program flow chart of the inventive method;
Fig. 2 is the system wiring figure of embodiment power distribution network.
Among the figure: 1-390 is the node serial number of power distribution network, and No. 1 node is the first node of power distribution network, also is the root node of power distribution network, at the 54th, 178,238 and No. 334 node reactive power compensator is set respectively.
Embodiment
Below in conjunction with embodiment, further specify the present invention.
Embodiment
Like Fig. 1, shown in 2, the concrete steps of the power distribution network wide area idle work optimization operation method of a kind of certain 10kV are following:
(1) input master data
At first import basic data (basic data that comprises each node, i.e. node serial number, node place electric pressure, the node voltage upper limit, the node lower limit circuit of certain 10kV power distribution network shown in Figure 2; The basic data of each circuit, the i.e. first end-node numbering of circuit, resistance
R l , reactance
X l , susceptance
B l , rated voltage
U Bl The basic data of each distribution transformer, i.e. distribution transformer institute on-load numbering, place, both sides node serial number, resistance
R T , reactance
X T , electricity leads
G m , susceptance
B m , no-load voltage ratio
k T , rated capacity
S NT , the high-pressure side rated voltage
U Bl1 , the low-pressure side rated voltage
U Bl2 ).The basic data of each reactive power compensator is: single pool-size
Q Ref =0.015MVar, the group number
n=7, total capacity
Q c =0.105MVar, maximum compensation capacity
Q Max =0.105MVar, minimum compensation capacity
Q Min =0.The historical switching data of each reactive power compensator see Table 1 ~ table 4.
The historical switching data of No. 1 compensation arrangement of table 1 (the place node serial number is 54)
The historical switching data of No. 2 compensation arrangements of table 2 (the place node serial number is 178)
The historical switching data of No. 3 compensation arrangements of table 3 (the place node serial number is 238)
The historical switching data of No. 4 compensation arrangements of table 4 (the place node serial number is 334)
The metric data of power distribution network root node and each reactive-load compensation point current time sees Table 5.
The metric data of table 5 power distribution network root node and each reactive power compensation point current time
(2) calculate the control priority level and the control range of each reactive power source
(1) accomplish after, according to the annexation of power distribution network circuit and transformer, begin to utilize each reactive power source point of wide area first search from root node (i.e. No. 1 node, first node), Search Results is as shown in table 6.
Each reactive power source priority level Search Results of table 6
Equivalent loss when by each reactive power compensator reactive power being provided respectively by each load of formula (1) calculating, table 7 are No. 1 each corresponding reactive power source equivalence losses of load.
The equivalent loss of the reactive power source that No. 1 load of table 7 is corresponding
Know that by table 7 No. 1 load should insert the 0th grade of reactive power source (power distribution network root node), No. 1 load place distribution transformer should belong to the control range of the 0th grade of reactive power source.By that analogy, confirm the interior distribution transformer of control range of each reactive power source, as shown in table 8.
The control range of each reactive power source of table 8
(3) calculate the switching Capacity Ratio of each reactive power compensator
After (2) completion, the control range that gets the 1st grade of reactive power source according to table 8 is the 13rd ~ No. 31 load, the switching Capacity Ratio of this compensation arrangement
In like manner calculate the switching Capacity Ratio of all the other compensation arrangements, see table 9.
The switching Capacity Ratio of each compensation arrangement of table 9
Then,
A=[26.41%, 26.48%, 18.12%, 15.80%]
T
(4) the idle sensitivity matrix of calculating voltage
After (3) completion, in the present embodiment, voltage monitoring point is power distribution network root node and reactive power compensator place node, i.e. the 1st, 54,178,238 and No. 334 node.
Can get Δ by table 1
Q C1 =0.015 o'clock, each voltage monitoring point
Δ U=[0,0.00098,0.00019,0.00019,0.00019]
TTherefore, the sensitivity between voltage monitoring point voltage and No. 1 reactive power compensator compensation capacity does
d
U/d
Q C1 =[0,0.06546,0.01286,0.01289,0.01292]
T
In like manner, the sensitivity between calculating voltage monitoring point voltage and other reactive power compensator compensation capacities:
(5) calculate the switching capacity of each reactive power compensator
(4) step was calculated the idle total load of power distribution network after accomplishing
Q L , computing formula does
Calculate the reactive load demand in the reactive power compensation point institute control range No. 1, computing formula is:
The reactive requirement of other reactive power compensators by that analogy, result of calculation is seen table 10.
Reactive requirement in each compensation arrangement control range of table 10
From reactive requirement result of calculation; The reactive requirement of No. 3 reactive power compensator in compensation self institute control area; Also remain the 0.045Mvar capacity; And No. 4 reactive power compensator can compensate the reactive requirement in self institute control area, and the compensation capacity of No. 2 reactive power compensator is not enough to compensate the reactive requirement in self institute control area.Therefore, No. 3 reactive power compensator should be supported reactive power compensator part reactive power No. 2, that is:
The reactive requirement of revising in each reactive power compensator control range of back is as shown in table 11.
Table 11 is revised the reactive requirement in each compensation arrangement control range of back
Consider that voltage retrains and the capacity of distribution transform constraint, form Optimization Model, adopt enumerative technique to separate this optimization problem, get the compensation capacity of each reactive power compensator by formula (9) ~ (12)
Q C =[0.105,0.105,0.105,0.090]
TWith this compensation capacity is example, and whether voltage monitoring point magnitude of voltage and distribution transforming apparent power meet the demands behind the calculable capacitor switching.
Calculate the reactive power compensator capacitor and drop into capacity
Δ Q C , computing formula is:
The voltage vector of next each voltage monitoring point of the moment after calculable capacitor puts into operation
U, computing formula is:
=
All do not cross the voltage range of the relevant guide rule regulation of country.
The capacitor that the calculates No. 1 reactive power compensator corresponding distribution transforming apparent power in back that puts into operation, computing formula is:
The rated capacity of this distribution transforming is 0.315, and capacitor puts back distribution transforming apparent power into operation above the capacity of distribution transform restriction.In like manner can calculate the apparent energy of all the other reactive power compensator place distribution transformings, see Table 12.
The apparent power of each reactive power compensator place distribution transforming of table 12
Can find out that from table 12 capacitor puts back distribution transforming apparent power into operation all above the capacity of distribution transform restriction.
Therefore, the compensation capacity of each reactive power compensator
Q C =[0.105,0.105,0.105,0.090]
T, need the input capacity
Δ Q C =[0.090,0.090,0.090,0.075]
T
(6) export the switching capacity of each reactive power compensator
After the completion of (5) step, the switching capacity of the reactive power compensator switching capacity that output solves
Δ Q C , be the switching amount of each reactive power compensator capacitance.
Should distinguish switching during each reactive power compensator switched capacitor, and record switching front and back capacitor compensation capacity and voltage monitoring point voltage condition, with input data as next moment.
Experiment effect
After the reactive power compensator of the wide area control mode of employing present embodiment puts into operation, carry out capacitor switching according to compensatory zone, control the reactive power compensation mode on the spot relatively with existing, the utilization rate of equipment and installations of each reactive power compensator is as shown in table 13.
Reactive-load compensation equipment utilance under table 13 different control modes
After each reactive power compensator employing wide area control mode of table 13 knowledge, each reactive power compensator utilization rate of equipment and installations obviously improves.
After adopting the inventive method that reactive power compensator is carried out wide area control, the voltage power-less situation result of calculation of this circuit some day is as shown in table 14.
Table 14 compensation back circuit operation result
Know by table 14, adopt wide area control mode of the present invention after, the quality of voltage of power distribution network obviously improves, line loss rate obviously reduces.