CN101982917A - Calculation method of available transmission capacity for power grid scheduling - Google Patents

Calculation method of available transmission capacity for power grid scheduling Download PDF

Info

Publication number
CN101982917A
CN101982917A CN 201010509246 CN201010509246A CN101982917A CN 101982917 A CN101982917 A CN 101982917A CN 201010509246 CN201010509246 CN 201010509246 CN 201010509246 A CN201010509246 A CN 201010509246A CN 101982917 A CN101982917 A CN 101982917A
Authority
CN
China
Prior art keywords
node
parameters
electrical network
power
active power
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.)
Granted
Application number
CN 201010509246
Other languages
Chinese (zh)
Other versions
CN101982917B (en
Inventor
陈金富
奚江惠
刘超
李勇
洪峰
姜曼
徐友平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huazhong University of Science and Technology
Central China Grid Co Ltd
Original Assignee
Huazhong University of Science and Technology
Central China Grid Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Huazhong University of Science and Technology, Central China Grid Co Ltd filed Critical Huazhong University of Science and Technology
Priority to CN 201010509246 priority Critical patent/CN101982917B/en
Publication of CN101982917A publication Critical patent/CN101982917A/en
Application granted granted Critical
Publication of CN101982917B publication Critical patent/CN101982917B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention discloses a calculation method of available transmission capacity for power grid scheduling, relating to the operation and the scheduling of a power system, which is an improved continuous power flow method. The calculation method of the invention is characterized in that the calculation method is dependent on the continuous power flow method; the current state is taken as an initial point; when moving forward to a new state point by controlling the step length, a thermal stabilization constrain, a transient state stability constrain and a node voltage constrain in engineering application after N-1 are verified; the step length is changed and a new predicted position is updated if the requirements cannot be satisfied; and the new state point is confirmed by utilizing Newton iteration until the section power difference between two state points meets the conditions of convergence. The invention has the advantages that the established available transmission capacity calculation model is accurate and flexible; the considered constrain conditions and the calculation model can be selected according to the requirements of users; the processing method for various constrains is more accurate so as to guarantee the reasonability of the available transmission capacity calculation; and the method can be applied to the available transmission capacity calculation of large-scale interconnected networks.

Description

The computational methods that are used for the available transmission capacity of dispatching of power netwoks
Technical field
The present invention relates to dispatching of power netwoks, relate in particular to the computational methods of the available transmission capacity that is used for dispatching of power netwoks.
Background technology
Because China's electric power resource skewness, length is occupied very big ratio apart from large capacity transmission in China; Power load goes up year by year in addition, makes some place of electrical network operate in the edge of its safety margins; In addition, the generation of various natural calamities allows near the electrical network of the limit more bogged down in crisis.Therefore, under electrical network limit operational mode, the calculating of available transfer capacity of transmission network is the basis that guarantees the power grid security reliability service, under the constraints of available transfer capacity of transmission network, adopt rational dispatching method, not only can take into account each side's interests, can also guarantee the economy of operation of power networks.
Available transmission capacity (ATC-Available Transfer Capability) is meant to guarantee electric power netting safe running and guaranteeing can also carry how many electric energy to another zone again by a zone in the electrical network under the constant situation of electricity transaction.
ATC is an important information of guaranteeing the power grid security reliability service, and calculating ATC accurately is the prerequisite that guarantees the power grid security reliability service.Existing ATC computational methods have the method that adopts the continuous tide method to determine ATC, and the general model that adopts of these class methods is the load variations model, but this class model can not react real-time dispatching requirement; The inequality constraints of Kao Lving simultaneously generally can not or only can the considering transient scleronomic constraint, and the constraints of consideration is not very comprehensive, and therefore the ATC that obtains can not guarantee accuracy fully.
It is a basic calculating of carrying out the power system mesomeric state operating analysis that electric power system tide calculates, and it is the basis of research and analysis electric power system.In power system planning design with have now in the research of power system operation mode, reasonability, reliability and the economy that need utilize trend to calculate to analyze comparison power supply plan or operational mode quantitatively.The trend result calculated is called the trend result.Generally adopt perunit value to represent for conveniently calculating trend calculated data and trend result, represent with p.u.
Conventional trend Jacobi submatrix is meant the submatrix of Jacobian matrix in the trend calculating, uses respectively
Figure 287717DEST_PATH_IMAGE001
,
Figure 575610DEST_PATH_IMAGE002
,
Figure 711931DEST_PATH_IMAGE003
, Expression, after power flow equation was determined, the method for solving of these submatrixs was well-determined.[" power system analysis (volume two) " He Yangzan, Wen Zengyin.]
In calculating, trend needs to carry out a kind of iterative computation, so-called iterative computation is meant a kind of method of finding the solution the higher-dimension equation group, the change amount that this method is separated by continuous solving equation comes the update equation approximate solution to infinity near till truly separating, and wherein the equation of the accounting equation change amount of separating is called iterative equation.
The continuous tide method is the stable method of using of a kind of calculating voltage, it can follow the tracks of trend result's track, from an initial launch state, progressively increase sending between survey region and be subjected to electric weight, up to the static voltage stability limit or overload appears in equipment, i.e. the critical strength of current running status point of system.The running status point is meant and comprises
Figure 224132DEST_PATH_IMAGE005
Five operation of power networks states that parameter is determined, Be the node phase angle vector
Figure 945193DEST_PATH_IMAGE007
,
Figure 355183DEST_PATH_IMAGE008
Be the phase angle of j node,
Figure 856703DEST_PATH_IMAGE009
Be the node voltage vector
Figure 680040DEST_PATH_IMAGE010
,
Figure 659629DEST_PATH_IMAGE011
Be the voltage of j node, Be node active power vector
Figure 924443DEST_PATH_IMAGE013
,
Figure 408645DEST_PATH_IMAGE014
Be the active power of j node
Figure 12496DEST_PATH_IMAGE015
,
Figure 898544DEST_PATH_IMAGE016
Be node reactive power vector
Figure 804183DEST_PATH_IMAGE017
,
Figure 274216DEST_PATH_IMAGE018
It is the reactive power of j node
Figure 595607DEST_PATH_IMAGE019
, the generator active power of each node wherein , reactive power
Figure 239133DEST_PATH_IMAGE021
If wherein j node do not have generator then assignment 0; The load active power of each node
Figure 697927DEST_PATH_IMAGE022
, reactive power
Figure 383861DEST_PATH_IMAGE023
If wherein j node do not loaded then assignment 0;
Figure 244501DEST_PATH_IMAGE024
Available transmission capacity for transmission cross-section; The subscript of these five parameters is represented the running status point,
Figure 491943DEST_PATH_IMAGE025
The parameter of representing k running status point is represented the initial launch state point when k=0.
After power system transient stability was meant that electrical network is subjected to big disturbance, can each synchronous machine continue to keep the ability of synchronous operation, and transient stability commonly used retrains and describes in electric power calculates.The computational methods of transient stability analysis are generally used the implicit expression trapezoidal integration, and the transient stability analysis module is to utilize computer to carry out the program module of transient stability analysis automatically.
Summary of the invention
The objective of the invention is: a kind of computational methods that are used for the available transmission capacity of dispatching of power netwoks are provided.This method is that electrical network is specified transmission cross-section transmission of electricity Calculation of Limit, and resulting trend operational factor and section available transmission capacity value can be used for planning, operation, the scheduling of electric power system and use.Advantage of the present invention is: model is accurate, flexible, and the consideration of constraints is comparatively comprehensive, and processing method rationally, has accurately guaranteed result's reliability, can be applicable to the scheduling of large-scale interconnected power system.
For achieving the above object, the present invention adopts following technical scheme:
The computational methods that are used for the available transmission capacity of dispatching of power netwoks, this method comprises the following step:
A) in the electrical network that m+n+1 node arranged, PQ node m wherein, PV node n, 1 of balance node (He Yangzan etc., power system analysis (descending). Wuhan: publishing house of the Central China University of Science and Technology .pp.52-53 in 2002); M+n+1 node is divided into two zones of A, B, and there is N1 generating set the a-quadrant, and there is N2 generating set in the B zone, and the set of a-quadrant and interregional all transmission lines of B is called transmission cross-section, and the active power of transmission cross-section flows to the B zone by the a-quadrant;
B) the step a) electrical network is carried out initial value design: the generator active power of setting each node
Figure 4744DEST_PATH_IMAGE026
, reactive power
Figure 933517DEST_PATH_IMAGE027
If wherein j node do not have generator then assignment 0; Set the load active power of each node
Figure 147199DEST_PATH_IMAGE022
, reactive power If wherein j node do not loaded then assignment 0; The trend result of electrical network converges on
Figure 169567DEST_PATH_IMAGE029
,
Figure 902031DEST_PATH_IMAGE030
Be the node phase angle vector
Figure 970219DEST_PATH_IMAGE031
, Be the phase angle of j node, Be the node voltage vector
Figure 329151DEST_PATH_IMAGE034
,
Figure 815628DEST_PATH_IMAGE011
Be the voltage of j node,
Figure 779035DEST_PATH_IMAGE035
Be node active power vector
Figure 419970DEST_PATH_IMAGE036
,
Figure 228657DEST_PATH_IMAGE037
Be the active power of j node
Figure 271437DEST_PATH_IMAGE038
,
Figure 405747DEST_PATH_IMAGE039
Be node reactive power vector
Figure 206081DEST_PATH_IMAGE040
,
Figure 588433DEST_PATH_IMAGE041
It is the reactive power of j node
Figure 49501DEST_PATH_IMAGE042
,
Figure 89133DEST_PATH_IMAGE043
The parameter of representing k running status point is represented the initial launch state point when k=0;
Figure 642343DEST_PATH_IMAGE044
Available transmission capacity for transmission cross-section;
C) order
Figure 556947DEST_PATH_IMAGE045
Give step-length step assignment, 0.001≤step≤0.2; M represents iterations, gives maximum iterations
Figure 810205DEST_PATH_IMAGE046
Assignment, 15≤
Figure 784852DEST_PATH_IMAGE046
≤ 30; Give the convergence criterion infinite decimal
Figure 123560DEST_PATH_IMAGE047
Assignment, 0.000001≤
Figure 773679DEST_PATH_IMAGE048
≤ 0.0001; Give the node voltage constrained vector
Figure 147023DEST_PATH_IMAGE049
With Assignment,
Figure 554794DEST_PATH_IMAGE051
,
Figure 312665DEST_PATH_IMAGE052
Be the minimum voltage constraint of j node, 0.95p.u≤
Figure 39051DEST_PATH_IMAGE053
≤ 1.0p.u,
Figure 591386DEST_PATH_IMAGE054
,
Figure 409081DEST_PATH_IMAGE055
Be the ceiling voltage constraint of j node, 1.0p.u≤
Figure 908327DEST_PATH_IMAGE056
≤ 1.05p.u; Give grid power transmission circuit active power thermoae limit vector according to the actual state of each transmission line
Figure 489219DEST_PATH_IMAGE057
Assignment,
Figure 540351DEST_PATH_IMAGE058
, Be the thermoae limit of active power of i bar transmission line, N is a grid power transmission circuit sum;
D) make k=0, get trend result in the step b)
Figure 1474DEST_PATH_IMAGE060
, order again , these five parameters constitute the initial launch state of electrical network in the step a);
E) introduce equation group (1) and j node generator active power
Figure 924486DEST_PATH_IMAGE020
Expression formula (2)
Figure 9116DEST_PATH_IMAGE062
（1）
（2）
Wherein: Be conventional trend Jacobi submatrix,
Figure 878479DEST_PATH_IMAGE065
,
Figure 152203DEST_PATH_IMAGE066
Be the active power of j node generator when being the generator node; Calculation equation (1) obtains
Figure 655997DEST_PATH_IMAGE067
The change amount Calculation equation (3) obtains
Figure 629824DEST_PATH_IMAGE069
Figure 626730DEST_PATH_IMAGE070
（3）
F) with result in the step e)
Figure 934214DEST_PATH_IMAGE071
The substitution equation
Figure 526607DEST_PATH_IMAGE072
?（4）
Obtain
Figure 199028DEST_PATH_IMAGE073
, in the equation (4)
Figure 745547DEST_PATH_IMAGE074
Represent that j node and the internodal electricity of i lead,
Figure 33221DEST_PATH_IMAGE075
Represent j node and the internodal susceptance of i, Represent j node and the internodal phase angle difference of i; Will
Figure 323443DEST_PATH_IMAGE077
The substitution iterative equation
Figure 294941DEST_PATH_IMAGE078
（5）
Carry out iterative computation;
G) work as iterations
Figure 944228DEST_PATH_IMAGE079
, during the iterative computation convergence, make k=k+1, obtain Five parameters, these five parameters constitute first running status point, carry out step h);
Work as iterations
Figure 525437DEST_PATH_IMAGE081
, when iterative computation does not still restrain, if step-length step greater than
Figure 217187DEST_PATH_IMAGE082
, order , return step e); If step-length step less than
Figure 593122DEST_PATH_IMAGE082
, stop to calculate output
Figure 270923DEST_PATH_IMAGE084
Five parameters, these parameters can be used for the traffic control of electrical network;
H) when the operation of power networks state parameter
Figure 13751DEST_PATH_IMAGE085
Satisfy the transient stability constraint simultaneously; The proper noun of N-1(power industry, the expression each element in the electric power system is cut-off one by one) after thermally-stabilised constraint
Figure 4841DEST_PATH_IMAGE086
,
Figure 782304DEST_PATH_IMAGE087
Active power for the back out of service of any transmission line in electrical network j bar circuit; And node voltage constraint
Figure 699182DEST_PATH_IMAGE088
The time, carry out step I);
Otherwise, if step-length step greater than
Figure 194885DEST_PATH_IMAGE082
, order
Figure 724087DEST_PATH_IMAGE089
, return step e); If step-length step less than
Figure 57854DEST_PATH_IMAGE082
, stop to calculate output
Figure 647098DEST_PATH_IMAGE090
Five parameters, these parameters can be used for the traffic control of electrical network;
I) introduce equation group (6)
Figure 364518DEST_PATH_IMAGE091
（6）
Wherein:
Figure 697411DEST_PATH_IMAGE092
Be that m+n+p element is 1, all the other are 0 column vector entirely,
Figure 948001DEST_PATH_IMAGE093
,
Figure 708147DEST_PATH_IMAGE094
, p is
Figure 912863DEST_PATH_IMAGE095
Subscript,
Figure 491524DEST_PATH_IMAGE096
Direction is followed the trail of in expression, when
Figure 98086DEST_PATH_IMAGE097
The time, , when
Figure 721145DEST_PATH_IMAGE099
The time, Calculation equation (6) obtains
Figure 621023DEST_PATH_IMAGE101
The change amount
Figure 722971DEST_PATH_IMAGE102
Calculation equation (3) obtains
Figure 902280DEST_PATH_IMAGE103
J) with step I) middle result
Figure 816884DEST_PATH_IMAGE104
Substitution equation (4) obtains
Figure 132459DEST_PATH_IMAGE073
Will
Figure 670888DEST_PATH_IMAGE105
The substitution iterative equation
Figure 570448DEST_PATH_IMAGE106
（7）
Carry out iterative computation;
K) work as iterations
Figure 852525DEST_PATH_IMAGE079
, during the iterative computation convergence, make k=k+1, obtain
Figure 22606DEST_PATH_IMAGE025
Five parameters, these five parameters constitute k running status point, carry out step l);
Work as iterations
Figure 731936DEST_PATH_IMAGE081
, when iterative computation does not still restrain, if step-length step greater than , order , return step I); If step-length step less than
Figure 246727DEST_PATH_IMAGE082
, stop to calculate output
Figure 126958DEST_PATH_IMAGE025
Five parameters, these parameters can be used for the traffic control of electrical network;
L) when the operation of power networks state parameter
Figure 1111DEST_PATH_IMAGE108
Satisfy the transient stability constraint simultaneously; Thermally-stabilised constraint behind the N-1
Figure 562673DEST_PATH_IMAGE109
, Active power for the back out of service of any transmission line in electrical network j bar circuit; And node voltage constraint The time, carry out step m);
Otherwise, if step-length step greater than
Figure 119928DEST_PATH_IMAGE082
, order
Figure 281919DEST_PATH_IMAGE111
, return step I); If step-length step less than , output
Figure 1669DEST_PATH_IMAGE112
Five parameters, these parameters can be used for the traffic control of electrical network;
M) when
Figure 351878DEST_PATH_IMAGE113
, and
Figure 317560DEST_PATH_IMAGE114
The time, carry out step I); When
Figure 878904DEST_PATH_IMAGE115
, and
Figure 271839DEST_PATH_IMAGE116
The time, output Five parameters, these parameters can be used for the traffic control of electrical network;
When
Figure 347559DEST_PATH_IMAGE117
, and step-length step greater than
Figure 85446DEST_PATH_IMAGE082
The time, order
Figure 649283DEST_PATH_IMAGE118
, return step j); When
Figure 708505DEST_PATH_IMAGE119
, and step-length step less than
Figure 452208DEST_PATH_IMAGE082
The time, output
Figure 546066DEST_PATH_IMAGE120
Five parameters, these parameters can be used for the traffic control of electrical network.
Advantage of the present invention is:
1, the Calculation of Available Transfer Capability model of being set up is accurate, flexible, and the constraints of consideration and computation model can be chosen in principle as required;
2, more accurate to the processing method of various constraints, the reasonability of assurance Calculation of Available Transfer Capability;
3, five parameters obtaining of the present invention
Figure 280804DEST_PATH_IMAGE121
, having provided a kind of dispatching method that arrives transmission cross-section transmitted power warning value, the dispatcher will reasonably arrange operation plan as warning line, and its result can be applied to the scheduling of large-scale interconnected power system.
Description of drawings
Fig. 1 is 8 machines, 36 node systems.
Embodiment:
Below in conjunction with accompanying drawing, the present invention is further illustrated.
Embodiment one
The computational methods that are used for the available transmission capacity of dispatching of power netwoks, this method comprises the following step:
A) in the electrical network that m+n+1 node arranged, PQ node m wherein, PV node n, 1 of balance node (He Yangzan etc., power system analysis (descending). Wuhan: publishing house of the Central China University of Science and Technology .pp.52-53 in 2002); M+n+1 node is divided into two zones of A, B, and there is N1 generating set the a-quadrant, and there is N2 generating set in the B zone, and the set of a-quadrant and interregional all transmission lines of B is called transmission cross-section, and the active power of transmission cross-section flows to the B zone by the a-quadrant;
B) the step a) electrical network is carried out initial value design: the generator active power of setting each node , reactive power
Figure 171454DEST_PATH_IMAGE027
If wherein j node do not have generator then assignment 0; Set the load active power of each node , reactive power
Figure 25458DEST_PATH_IMAGE028
If wherein j node do not loaded then assignment 0; The trend result of electrical network converges on
Figure 59273DEST_PATH_IMAGE029
,
Figure 201236DEST_PATH_IMAGE030
Be the node phase angle vector
Figure 4107DEST_PATH_IMAGE031
,
Figure 80647DEST_PATH_IMAGE032
Be the phase angle of j node,
Figure 601759DEST_PATH_IMAGE033
Be the node voltage vector
Figure 553272DEST_PATH_IMAGE034
, Be the voltage of j node,
Figure 458091DEST_PATH_IMAGE035
Be node active power vector
Figure 466498DEST_PATH_IMAGE036
,
Figure 956123DEST_PATH_IMAGE037
Be the active power of j node
Figure 733586DEST_PATH_IMAGE038
,
Figure 151929DEST_PATH_IMAGE039
Be node reactive power vector
Figure 647633DEST_PATH_IMAGE040
,
Figure 675369DEST_PATH_IMAGE041
It is the reactive power of j node
Figure 245022DEST_PATH_IMAGE042
,
Figure 99846DEST_PATH_IMAGE043
The parameter of representing k running status point is represented the initial launch state point when k=0; Available transmission capacity for transmission cross-section;
C) order
Figure 654553DEST_PATH_IMAGE045
Give step-length step assignment, 0.001≤step≤0.2; M represents iterations, gives maximum iterations
Figure 406608DEST_PATH_IMAGE046
Assignment, 15≤
Figure 166754DEST_PATH_IMAGE046
≤ 30; Give the convergence criterion infinite decimal
Figure 870005DEST_PATH_IMAGE047
Assignment, 0.000001≤
Figure 6588DEST_PATH_IMAGE048
≤ 0.0001; Give the node voltage constrained vector With
Figure 544197DEST_PATH_IMAGE050
Assignment, , Be the minimum voltage constraint of j node, 0.95p.u≤
Figure 870508DEST_PATH_IMAGE053
≤ 1.0p.u,
Figure 238036DEST_PATH_IMAGE054
,
Figure 915879DEST_PATH_IMAGE055
Be the ceiling voltage constraint of j node, 1.0p.u≤
Figure 66369DEST_PATH_IMAGE056
≤ 1.05p.u; Give grid power transmission circuit active power thermoae limit vector according to the actual state of each transmission line
Figure 647523DEST_PATH_IMAGE057
Assignment,
Figure 436486DEST_PATH_IMAGE058
,
Figure 103091DEST_PATH_IMAGE059
Be the thermoae limit of active power of i bar transmission line, N is a grid power transmission circuit sum;
D) make k=0, get trend result in the step b)
Figure 119588DEST_PATH_IMAGE060
, order again
Figure 555249DEST_PATH_IMAGE061
, these five parameters constitute the initial launch state of electrical network in the step a);
E) introduce equation group (1) and j node generator active power
Figure 763114DEST_PATH_IMAGE020
Expression formula (2)
Figure 651435DEST_PATH_IMAGE062
（1）
Figure 471624DEST_PATH_IMAGE063
（2）
Wherein:
Figure 761791DEST_PATH_IMAGE064
Be conventional trend Jacobi submatrix, ,
Figure 516175DEST_PATH_IMAGE066
Be the active power of j node generator when being the generator node; Calculation equation (1) obtains
Figure 874475DEST_PATH_IMAGE067
The change amount
Figure 284728DEST_PATH_IMAGE068
Calculation equation (3) obtains
Figure 772079DEST_PATH_IMAGE069
Figure 634993DEST_PATH_IMAGE070
（3）
F) with result in the step e)
Figure 796984DEST_PATH_IMAGE071
The substitution equation
Figure 300558DEST_PATH_IMAGE072
?（4）
Obtain
Figure 522592DEST_PATH_IMAGE073
, in the equation (4)
Figure 872802DEST_PATH_IMAGE074
Represent that j node and the internodal electricity of i lead,
Figure 838484DEST_PATH_IMAGE075
Represent j node and the internodal susceptance of i,
Figure 456285DEST_PATH_IMAGE076
Represent j node and the internodal phase angle difference of i; Will
Figure 583641DEST_PATH_IMAGE077
The substitution iterative equation
Figure 421147DEST_PATH_IMAGE078
（5）
Carry out iterative computation;
G) work as iterations
Figure 924941DEST_PATH_IMAGE079
, during the iterative computation convergence, make k=k+1, obtain Five parameters, these five parameters constitute first running status point, carry out step h);
Work as iterations
Figure 961084DEST_PATH_IMAGE081
, when iterative computation does not still restrain, if step-length step greater than
Figure 285887DEST_PATH_IMAGE082
, order
Figure 593371DEST_PATH_IMAGE083
, return step e); If step-length step less than , stop to calculate output
Figure 654923DEST_PATH_IMAGE084
Five parameters, these parameters can be used for the traffic control of electrical network;
H) when the operation of power networks state parameter Satisfy the transient stability constraint simultaneously; Thermally-stabilised constraint behind the N-1
Figure 742976DEST_PATH_IMAGE086
,
Figure 691341DEST_PATH_IMAGE087
Active power for the back out of service of any transmission line in electrical network j bar circuit; And node voltage constraint
Figure 596980DEST_PATH_IMAGE088
The time, carry out step I);
Otherwise, if step-length step greater than
Figure 630795DEST_PATH_IMAGE082
, order , return step e); If step-length step less than , stop to calculate output
Figure 658028DEST_PATH_IMAGE090
Five parameters, these parameters can be used for the traffic control of electrical network;
I) introduce equation group (6)
Figure 179140DEST_PATH_IMAGE091
（6）
Wherein:
Figure 865074DEST_PATH_IMAGE092
Be that m+n+p element is 1, all the other are 0 column vector entirely,
Figure 788030DEST_PATH_IMAGE093
,
Figure 35472DEST_PATH_IMAGE094
, p is Subscript,
Figure 533504DEST_PATH_IMAGE096
Direction is followed the trail of in expression, when
Figure 310967DEST_PATH_IMAGE097
The time,
Figure 729310DEST_PATH_IMAGE098
, when The time,
Figure 258610DEST_PATH_IMAGE100
Calculation equation (6) obtains
Figure 890579DEST_PATH_IMAGE101
The change amount
Figure 479824DEST_PATH_IMAGE102
Calculation equation (3) obtains
Figure 899041DEST_PATH_IMAGE103
J) with step I) middle result
Figure 231934DEST_PATH_IMAGE104
Substitution equation (4) obtains
Figure 718410DEST_PATH_IMAGE073
Will
Figure 744135DEST_PATH_IMAGE105
The substitution iterative equation
Figure 447386DEST_PATH_IMAGE106
（7）
Carry out iterative computation;
K) work as iterations
Figure 583970DEST_PATH_IMAGE079
, during the iterative computation convergence, make k=k+1, obtain
Figure 924952DEST_PATH_IMAGE025
Five parameters, these five parameters constitute k running status point, carry out step l);
Work as iterations
Figure 121578DEST_PATH_IMAGE081
, when iterative computation does not still restrain, if step-length step greater than
Figure 312126DEST_PATH_IMAGE082
, order , return step I); If step-length step less than
Figure 447889DEST_PATH_IMAGE082
, stop to calculate output Five parameters, these parameters can be used for the traffic control of electrical network;
L) when the operation of power networks state parameter
Figure 245259DEST_PATH_IMAGE108
Satisfy the transient stability constraint simultaneously; Thermally-stabilised constraint behind the N-1 , Active power for the back out of service of any transmission line in electrical network j bar circuit; And node voltage constraint
Figure 13867DEST_PATH_IMAGE088
The time, carry out step m);
Otherwise, if step-length step greater than
Figure 680472DEST_PATH_IMAGE082
, order , return step I); If step-length step less than
Figure 132630DEST_PATH_IMAGE082
, output Five parameters, these parameters can be used for the traffic control of electrical network;
M) when
Figure 228817DEST_PATH_IMAGE113
, and
Figure 49005DEST_PATH_IMAGE114
The time, carry out step I); When , and
Figure 717939DEST_PATH_IMAGE116
The time, output Five parameters, these parameters can be used for the traffic control of electrical network;
When
Figure 389539DEST_PATH_IMAGE117
, and step-length step greater than
Figure 304187DEST_PATH_IMAGE082
The time, order
Figure 89740DEST_PATH_IMAGE118
, return step j); When , and step-length step less than
Figure 114645DEST_PATH_IMAGE082
The time, output
Figure 877939DEST_PATH_IMAGE120
Five parameters, these parameters can be used for the traffic control of electrical network.
Embodiment two
As shown in Figure 1, this system is the China Power academy of sciences 8 machines, 36 node systems (Yan Wei, the dynamic reactive optimization of ac and dc systems, Automation of Electric Systems, 2009 the 10th phases).Wherein the PQ node is m=32, and the PV node is n=3,1 of balance node; Network connection figure such as Fig. 1 of this system.Carry out subregion as shown in Figure 1, there is 3 in generator in zone one, and there is 5 in generator in zone two, and zoning one is to the available transmission capacity in zone two.
Structural parameters and operational factor initialization electric network composition and initial trend result by the China Power academy of sciences 8 machines, 36 node systems.Step=0.1,
Figure 99973DEST_PATH_IMAGE046
=20,
Figure 450183DEST_PATH_IMAGE082
=0.000001; =0.8p.u,
Figure 768087DEST_PATH_IMAGE123
=1.05p.u;
Figure 161022DEST_PATH_IMAGE124
=4 Assignment,
Figure 502322DEST_PATH_IMAGE126
Be the active power of i bar transmission line, N=33;
Adopt the method for embodiment one, this basic mode is carried out Calculation of Available Transfer Capability, when k=7, promptly calculate end during the 7th running status point, obtain available transmission capacity
Figure 974629DEST_PATH_IMAGE127
=3.5191, trend result such as the table 1 of this moment.
Table 1 limiting condition is put each node data
Nodename Voltage Phase angle Meritorious exerting oneself Idle exerting oneself Load is meritorious Reactive load
BUS9 0.80722 -27.6424 0 0 3.76 2.21
BUS2 0.813387 -10.4483 8.32028 3.60001 0 0
BUS10 0.808067 -26.5822 0 0 0 0
BUS11 0.830682 -36.1294 0 0 0 0
BUS51 0.808067 -26.5822 0 0 0 0
BUS22 0.993827 -49.2818 0 0 2.265 1.69
BUS3 1 -45.9763 4.29881 5.99188 0 0
BUS23 0.809893 -27.8979 0 0 2.87 1.44
BUS24 0.990979 -6.79045 0 0 0 0
BUS1 1 0 0 5.38771 0 0
BUS12 0.854707 -53.1497 0 0 0 0
BUS15 0 0 0 0 0 0
BUS14 0.856212 -52.6408 0 0 0 0
BUS52 0.859011 -52.5791 0 0 0 0
BUS13 0.872514 -65.6005 0 0 0 0
BUS17 0.863908 -70.3774 0 0 0 0
BUS6 1 -70.413 -0.01595 3.85917 0 0
BUS16 0.884721 -70.8476 0 0 10.2 2.4
BUS19 0.885778 -52.3497 0 0 0.864 0.662
BUS4 0.912098 -49.337 0.647287 0.700008 0 0
BUS18 0.890714 -77.1939 0 0 4.3 2.2
BUS5 0.956683 -72.4733 1.73958 3.34004 0 0
BUS20 0.906612 -62.684 0 0 0.719 0.474
BUS21 0.906373 -54.2168 0 0 0.7 0.5
BUS30 0.988467 -40.1283 0 0 0 0
BUS7 1 -34.2629 2.25 0.88755 0 0
BUS8 1 -34.1543 3.06 1.19094 0 0
BUS50 0.890714 -77.1939 0 0 0 0
The node that does not provide in the last table is the node that has little branch road, has merged little branch road in the program.Show by top result of calculation, this ATC computational methods can effectively calculate the transmission of electricity limit point, avoided conventional trend in the unusual and imponderable problem of limit point Jacobian matrix, simultaneously can also consider multiple constraint easily, as seen this method is correct on the problem of calculating ATC.
Present embodiment is at the appointment transmission cross-section in the standard power system, after calculating by the inventive method, obtains reasonable available transmission capacity and trend result, and its result can be used for the traffic control of this electrical network.

Claims (1)

1. be used for the computational methods of the available transmission capacity of dispatching of power netwoks, it is characterized in that, this method comprises the following step:
A) in the electrical network that m+n+1 node arranged, PQ node m wherein, PV node n, 1 of balance node (He Yangzan etc., power system analysis (descending). Wuhan: publishing house of the Central China University of Science and Technology .pp.52-53 in 2002); M+n+1 node is divided into two zones of A, B, and there is N1 generating set the a-quadrant, and there is N2 generating set in the B zone, and the set of a-quadrant and interregional all transmission lines of B is called transmission cross-section, and the active power of transmission cross-section flows to the B zone by the a-quadrant;
B) the step a) electrical network is carried out initial value design: the generator active power of setting each node
Figure 267502DEST_PATH_IMAGE001
, reactive power
Figure 771296DEST_PATH_IMAGE002
If wherein j node do not have generator then assignment 0; Set the load active power of each node
Figure 446866DEST_PATH_IMAGE003
, reactive power
Figure 10702DEST_PATH_IMAGE004
If wherein j node do not loaded then assignment 0; The trend result of electrical network converges on
Figure 69925DEST_PATH_IMAGE005
,
Figure 377410DEST_PATH_IMAGE006
Be the node phase angle vector ,
Figure 642224DEST_PATH_IMAGE008
Be the phase angle of j node,
Figure 188743DEST_PATH_IMAGE009
Be the node voltage vector
Figure 532874DEST_PATH_IMAGE010
,
Figure 481239DEST_PATH_IMAGE011
Be the voltage of j node, Be node active power vector
Figure 686989DEST_PATH_IMAGE013
,
Figure 336276DEST_PATH_IMAGE014
Be the active power of j node
Figure 76830DEST_PATH_IMAGE015
,
Figure 589589DEST_PATH_IMAGE016
Be node reactive power vector
Figure 110700DEST_PATH_IMAGE017
,
Figure 563678DEST_PATH_IMAGE018
It is the reactive power of j node
Figure 719591DEST_PATH_IMAGE019
,
Figure 967033DEST_PATH_IMAGE020
The parameter of representing k running status point is represented the initial launch state point when k=0;
Figure 975440DEST_PATH_IMAGE021
Available transmission capacity for transmission cross-section;
C) order
Figure 966530DEST_PATH_IMAGE022
Give step-length step assignment, 0.001≤step≤0.2; M represents iterations, gives maximum iterations
Figure 242528DEST_PATH_IMAGE023
Assignment, 15≤
Figure 660871DEST_PATH_IMAGE023
≤ 30; Give the convergence criterion infinite decimal Assignment, 0.000001≤ ≤ 0.0001; Give the node voltage constrained vector
Figure 822140DEST_PATH_IMAGE026
With
Figure 676963DEST_PATH_IMAGE027
Assignment,
Figure 394384DEST_PATH_IMAGE028
, Be the minimum voltage constraint of j node, 0.95p.u≤
Figure 915550DEST_PATH_IMAGE030
≤ 1.0p.u, ,
Figure 880412DEST_PATH_IMAGE032
Be the ceiling voltage constraint of j node, 1.0p.u≤
Figure 515530DEST_PATH_IMAGE033
≤ 1.05p.u; Give grid power transmission circuit active power thermoae limit vector according to the actual state of each transmission line
Figure 122092DEST_PATH_IMAGE034
Assignment, ,
Figure 745151DEST_PATH_IMAGE036
Be the thermoae limit of active power of i bar transmission line, N is a grid power transmission circuit sum;
D) make k=0, get trend result in the step b)
Figure 918381DEST_PATH_IMAGE037
, order again
Figure 379450DEST_PATH_IMAGE038
, these five parameters constitute the initial launch state of electrical network in the step a);
E) introduce equation group (1) and j node generator active power
Figure 746977DEST_PATH_IMAGE039
Expression formula (2)
Figure 356645DEST_PATH_IMAGE040
（1）
Figure 569451DEST_PATH_IMAGE041
（2）
Wherein: Be conventional trend Jacobi submatrix,
Figure 423455DEST_PATH_IMAGE043
,
Figure 588595DEST_PATH_IMAGE044
Be the active power of j node generator when being the generator node; Calculation equation (1) obtains
Figure 605092DEST_PATH_IMAGE045
The change amount
Figure 40753DEST_PATH_IMAGE046
Calculation equation (3) obtains
Figure 750083DEST_PATH_IMAGE047
Figure 74623DEST_PATH_IMAGE048
（3）
F) with result in the step e)
Figure 832494DEST_PATH_IMAGE049
The substitution equation
Figure 621196DEST_PATH_IMAGE050
?（4）
Obtain
Figure 501428DEST_PATH_IMAGE051
, in the equation (4)
Figure 877045DEST_PATH_IMAGE052
Represent that j node and the internodal electricity of i lead,
Figure 235346DEST_PATH_IMAGE053
Represent j node and the internodal susceptance of i,
Figure 149993DEST_PATH_IMAGE054
Represent j node and the internodal phase angle difference of i; Will The substitution iterative equation
（5）
Carry out iterative computation;
G) work as iterations
Figure 662249DEST_PATH_IMAGE057
, during the iterative computation convergence, make k=k+1, obtain
Figure 723746DEST_PATH_IMAGE058
Five parameters, these five parameters constitute first running status point, carry out step h);
Work as iterations
Figure 945779DEST_PATH_IMAGE059
, when iterative computation does not still restrain, if step-length step greater than , order
Figure 261671DEST_PATH_IMAGE061
, return step e); If step-length step less than
Figure 879472DEST_PATH_IMAGE060
, stop to calculate output
Figure 944511DEST_PATH_IMAGE062
Five parameters, these parameters can be used for the traffic control of electrical network;
H) when the operation of power networks state parameter
Figure 782017DEST_PATH_IMAGE063
Satisfy the transient stability constraint simultaneously; Thermally-stabilised constraint behind the N-1
Figure 722029DEST_PATH_IMAGE064
,
Figure 961381DEST_PATH_IMAGE065
Active power for the back out of service of any transmission line in electrical network j bar circuit; And node voltage constraint
Figure 259638DEST_PATH_IMAGE066
The time, carry out step I);
Otherwise, if step-length step greater than
Figure 100553DEST_PATH_IMAGE060
, order
Figure 408038DEST_PATH_IMAGE067
, return step e); If step-length step less than
Figure 236317DEST_PATH_IMAGE060
, stop to calculate output Five parameters, these parameters can be used for the traffic control of electrical network;
I) introduce equation group (6)
Figure 16109DEST_PATH_IMAGE069
（6）
Wherein: Be that m+n+p element is 1, all the other are 0 column vector entirely,
Figure 137965DEST_PATH_IMAGE071
, , p is
Figure 513638DEST_PATH_IMAGE073
Subscript,
Figure 162925DEST_PATH_IMAGE074
Direction is followed the trail of in expression, when
Figure 965796DEST_PATH_IMAGE075
The time,
Figure 42336DEST_PATH_IMAGE076
, when
Figure 61983DEST_PATH_IMAGE077
The time, Calculation equation (6) obtains
Figure 172338DEST_PATH_IMAGE079
The change amount
Figure 419780DEST_PATH_IMAGE080
Calculation equation (3) obtains
Figure 932582DEST_PATH_IMAGE081
J) with step I) middle result Substitution equation (4) obtains
Figure 638818DEST_PATH_IMAGE051
Will
Figure 555696DEST_PATH_IMAGE083
The substitution iterative equation
（7）
Carry out iterative computation;
K) work as iterations
Figure 518284DEST_PATH_IMAGE057
, during the iterative computation convergence, make k=k+1, obtain Five parameters, these five parameters constitute k running status point, carry out step l);
Work as iterations
Figure 238033DEST_PATH_IMAGE059
, when iterative computation does not still restrain, if step-length step greater than
Figure 221032DEST_PATH_IMAGE060
, order
Figure 553925DEST_PATH_IMAGE086
, return step I); If step-length step less than
Figure 538936DEST_PATH_IMAGE060
, stop to calculate output
Figure 626978DEST_PATH_IMAGE085
Five parameters, these parameters can be used for the traffic control of electrical network;
L) when the operation of power networks state parameter Satisfy the transient stability constraint simultaneously; Thermally-stabilised constraint behind the N-1
Figure 968277DEST_PATH_IMAGE088
,
Figure 739619DEST_PATH_IMAGE089
Active power for the back out of service of any transmission line in electrical network j bar circuit; And node voltage constraint The time, carry out step m);
Otherwise, if step-length step greater than , order
Figure 302953DEST_PATH_IMAGE090
, return step I); If step-length step less than
Figure 262556DEST_PATH_IMAGE060
, output
Figure 630083DEST_PATH_IMAGE091
Five parameters, these parameters can be used for the traffic control of electrical network;
M) when
Figure 543813DEST_PATH_IMAGE092
, and
Figure 22199DEST_PATH_IMAGE093
The time, carry out step I); When
Figure 101888DEST_PATH_IMAGE094
, and
Figure 374737DEST_PATH_IMAGE095
The time, output
Figure 41342DEST_PATH_IMAGE087
Five parameters, these parameters can be used for the traffic control of electrical network;
When
Figure 57840DEST_PATH_IMAGE096
, and step-length step greater than
Figure 992035DEST_PATH_IMAGE060
The time, order
Figure 435786DEST_PATH_IMAGE097
, return step j); When
Figure 589687DEST_PATH_IMAGE098
, and step-length step less than
Figure 409875DEST_PATH_IMAGE060
The time, output
Figure 142120DEST_PATH_IMAGE099
Five parameters, these parameters can be used for the traffic control of electrical network.
CN 201010509246 2010-10-17 2010-10-17 Calculation method of available transmission capacity for power grid scheduling Active CN101982917B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201010509246 CN101982917B (en) 2010-10-17 2010-10-17 Calculation method of available transmission capacity for power grid scheduling

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201010509246 CN101982917B (en) 2010-10-17 2010-10-17 Calculation method of available transmission capacity for power grid scheduling

Publications (2)

Publication Number Publication Date
CN101982917A true CN101982917A (en) 2011-03-02
CN101982917B CN101982917B (en) 2013-03-27

Family

ID=43619812

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201010509246 Active CN101982917B (en) 2010-10-17 2010-10-17 Calculation method of available transmission capacity for power grid scheduling

Country Status (1)

Country Link
CN (1) CN101982917B (en)

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102437643A (en) * 2011-10-06 2012-05-02 华中电网有限公司 Method for determining available power transmission capacity of each power transmission section in multiple sections
CN102709955A (en) * 2012-05-30 2012-10-03 中国电力科学研究院 Multi-section-based power flow control method
CN103532157A (en) * 2013-08-27 2014-01-22 浙江工业大学 Electric-energy scheduling method for power grid system integrated into energy storage equipment
CN104734155A (en) * 2015-04-14 2015-06-24 河海大学 Method of acquiring available transfer capacity of electricity-gas interconnected energy system
CN105226638A (en) * 2015-09-06 2016-01-06 清华大学 Based on the electric power system Calculation of Available Transfer Capability method improving learning aid algorithm
CN105226643A (en) * 2015-09-23 2016-01-06 清华大学 Under security constraint, power system operation simulation model generates and method for solving fast
CN105244877A (en) * 2015-11-04 2016-01-13 国电南瑞科技股份有限公司 Recovery adjustment calculation method for power flow unsolvability
CN105284030A (en) * 2013-06-27 2016-01-27 宝马股份公司 Electrical power availability interface
CN105958496A (en) * 2016-07-01 2016-09-21 华北电力大学 Continuous power flow based maximum transmission capacity calculation method for system containing power source converter (VSC)
CN106532697A (en) * 2016-12-30 2017-03-22 国网山东省电力公司泰安供电公司 Method and device for building multi-section model
CN106684898A (en) * 2016-10-31 2017-05-17 国电南瑞科技股份有限公司 Value network-based scheduling optimization method of energy storage system
CN106786654A (en) * 2016-12-30 2017-05-31 国网山东省电力公司泰安供电公司 Power network section power method of adjustment and device
CN107370157A (en) * 2017-09-12 2017-11-21 四川大学 A kind of available transfer capacity of transmission network risk benefit decision method based on trend entropy

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101673948A (en) * 2009-09-28 2010-03-17 国网电力科学研究院 N-2 combined fault quick sequencing method for transient stability analysis of electric power system

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101673948A (en) * 2009-09-28 2010-03-17 国网电力科学研究院 N-2 combined fault quick sequencing method for transient stability analysis of electric power system

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
《2006中国电力系统保护与控制学术研讨会论文集》 20061231 汪天安等 输电网络可用输电能力的研究综述 , *

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102437643B (en) * 2011-10-06 2013-12-18 华中电网有限公司 Method for determining available power transmission capacity of each power transmission section in multiple sections
CN102437643A (en) * 2011-10-06 2012-05-02 华中电网有限公司 Method for determining available power transmission capacity of each power transmission section in multiple sections
CN102709955A (en) * 2012-05-30 2012-10-03 中国电力科学研究院 Multi-section-based power flow control method
CN102709955B (en) * 2012-05-30 2014-08-27 中国电力科学研究院 Multi-section-based power flow control method
CN105284030A (en) * 2013-06-27 2016-01-27 宝马股份公司 Electrical power availability interface
CN103532157B (en) * 2013-08-27 2016-06-15 浙江工业大学 A kind of for the network system electric energy dispatching method incorporating energy storage device
CN103532157A (en) * 2013-08-27 2014-01-22 浙江工业大学 Electric-energy scheduling method for power grid system integrated into energy storage equipment
CN104734155A (en) * 2015-04-14 2015-06-24 河海大学 Method of acquiring available transfer capacity of electricity-gas interconnected energy system
CN105226638A (en) * 2015-09-06 2016-01-06 清华大学 Based on the electric power system Calculation of Available Transfer Capability method improving learning aid algorithm
CN105226643A (en) * 2015-09-23 2016-01-06 清华大学 Under security constraint, power system operation simulation model generates and method for solving fast
CN105244877A (en) * 2015-11-04 2016-01-13 国电南瑞科技股份有限公司 Recovery adjustment calculation method for power flow unsolvability
CN105958496A (en) * 2016-07-01 2016-09-21 华北电力大学 Continuous power flow based maximum transmission capacity calculation method for system containing power source converter (VSC)
CN105958496B (en) * 2016-07-01 2019-04-02 华北电力大学 The total transfer capability calculation method of system containing VSC based on continuous tide
CN106684898A (en) * 2016-10-31 2017-05-17 国电南瑞科技股份有限公司 Value network-based scheduling optimization method of energy storage system
CN106532697A (en) * 2016-12-30 2017-03-22 国网山东省电力公司泰安供电公司 Method and device for building multi-section model
CN106786654A (en) * 2016-12-30 2017-05-31 国网山东省电力公司泰安供电公司 Power network section power method of adjustment and device
CN107370157A (en) * 2017-09-12 2017-11-21 四川大学 A kind of available transfer capacity of transmission network risk benefit decision method based on trend entropy

Also Published As

Publication number Publication date
CN101982917B (en) 2013-03-27

Similar Documents

Publication Publication Date Title
CN101982917B (en) Calculation method of available transmission capacity for power grid scheduling
Zamzam et al. Optimal water–power flow-problem: Formulation and distributed optimal solution
Arnold et al. Model-based predictive control applied to multi-carrier energy systems
Tang et al. Hydrothermal scheduling via extended differential dynamic programming and mixed coordination
CN102184297B (en) Electrical/control system decoupling prediction method suitable for transient parallel simulation of micro-nets
CN102346442B (en) Method for carrying out parallel transient simulation computing on microgrid suitable for nodal analysis frameworks
CN101976837A (en) Method for adjusting unsolvable power flow of electric power system
CN110165695B (en) Method and system for controlling multi-terminal direct current transmission system in layered mode
Wu et al. Voltage security enhancement via coordinated control
CN104361159A (en) Time-space parallel simulation method for transient stability of large-scale power system
CN104769802A (en) Method for the computer-aided control of the power in an electrical grid
CN108493998A (en) Consider the robust Transmission Expansion Planning in Electric method of demand response and N-1 forecast failures
CN102842917A (en) Universal electromechanical transient state model of grid-connected photovoltaic power generation system
CN105610168A (en) Application method of node voltage stability index of multi-infeed AC/DC system
Huang et al. Fast textured algorithms for optimal power delivery problems in deregulated environments
CN102437643B (en) Method for determining available power transmission capacity of each power transmission section in multiple sections
CN110020506B (en) Differential format selection method based on operation optimization of electric heating type comprehensive energy system
Zhang et al. Successive MISOCP algorithm for islanded distribution networks with soft open points
Zhang et al. An efficient approach for robust SCOPF considering load and renewable power uncertainties
Bélanger et al. A modern and open real-time digital simulator of contemporary power systems
Ko et al. A positioning method of distributed power system for minimum voltage variation in a DC microgrid
CN105989543A (en) Adaptive-model-based alternating-current connecting line equivalence method
Abdel-Akherf et al. Application of bifurcation analysis on active unbalanced distribution system
Carbone et al. Multi-Microgrid Transmission Control for a Lunar Surface Power System
Zare et al. An improved priority-based load shedding scheme in hybrid ac/dc microgrids using micro-pmudata

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant