CN101154813A  Online combination method for multiregion electric network tide model  Google Patents
Online combination method for multiregion electric network tide model Download PDFInfo
 Publication number
 CN101154813A CN101154813A CNA200710176083XA CN200710176083A CN101154813A CN 101154813 A CN101154813 A CN 101154813A CN A200710176083X A CNA200710176083X A CN A200710176083XA CN 200710176083 A CN200710176083 A CN 200710176083A CN 101154813 A CN101154813 A CN 101154813A
 Authority
 CN
 China
 Prior art keywords
 node
 subordinate
 network
 level
 electrical network
 Prior art date
Links
 239000011159 matrix materials Substances 0.000 claims description 39
 238000002347 injection Methods 0.000 claims description 34
 239000007924 injections Substances 0.000 claims description 34
 281000150264 Delta Cafés companies 0.000 claims description 24
 280000255884 Dispatching companies 0.000 claims description 21
 238000000034 methods Methods 0.000 claims description 16
 230000035945 sensitivity Effects 0.000 claims description 16
 230000001105 regulatory Effects 0.000 claims description 15
 230000000875 corresponding Effects 0.000 claims description 13
 230000001808 coupling Effects 0.000 claims description 11
 230000017105 transposition Effects 0.000 claims description 8
 230000015572 biosynthetic process Effects 0.000 claims description 3
 238000005755 formation reactions Methods 0.000 claims description 3
 210000001138 Tears Anatomy 0.000 claims description 2
 239000004576 sand Substances 0.000 claims description 2
 238000004458 analytical method Methods 0.000 abstract description 10
 238000004088 simulation Methods 0.000 abstract description 3
 238000010348 incorporation Methods 0.000 abstract 1
 238000010168 coupling process Methods 0.000 description 8
 238000005859 coupling reactions Methods 0.000 description 8
 238000000354 decomposition reactions Methods 0.000 description 4
 210000001503 Joints Anatomy 0.000 description 2
 238000010586 diagrams Methods 0.000 description 2
 230000011664 signaling Effects 0.000 description 2
 238000004364 calculation methods Methods 0.000 description 1
 238000004891 communication Methods 0.000 description 1
 238000005516 engineering processes Methods 0.000 description 1
 230000014509 gene expression Effects 0.000 description 1
 238000003032 molecular docking Methods 0.000 description 1
Abstract
The invention relates to an online incorporated method of a tide model of a multiregion electric network, which belongs to the electric power system dispatching automatic field. The invention comprises an upper level electric network dispatching center and a lower level electric network dispatching center, wherein, an internal electric network model and tide distribution are automatically generated in the region by utilizing state estimation method according to realtime data of the electric network under jurisdiction, and simultaneously an external electric network equivalent model connected with the internal network model is built up; lower level electric network models and tide data are sent to the upper level electric network dispatching center through the wide area network; upper level electric network models are connected with the lower level electric network models through respective connecting lines; the upper level electric network models and the lower level electric network models are incorporated into a global electric network model for counting; tide matching method is utilized to regulate the tide of the lower level electric network models, and a global electric network tide model is achieved. The method guarantees unchanged tide of intermediate and upper level electric network parts in the global electric network model after incorporation; the invention guarantees the accuracy of online safety analysis and simulation results of the interconnecting electric network, and is of great importance for guaranteeing safe operation of the electric network.
Description
Technical field
The invention belongs to the dispatching automation of electric power systems technical field, relate to multiregion electric network tide model in the electric power system online merging method, guarantee to merge the coupling trend method of higher level's electric network swim and higher level's realtime power network state consistency in the overall electric network model in back.
Background technology
Modern power network is formed by a plurality of regional power grids are interconnected, and accident is overall to the influence of electrical network.Yet the management of electrical network is that classification, subregion carry out at present, and the electrical network of each control centre's administration is a part of big electrical network.Current grid dispatching center EMS (EMS) is only carried out modeling to own administrative area electrical network, and (the remote measurement amount mainly comprises the voltage of the active power of the active power of generator active power, reactive power and machine end busbar voltage, load and reactive power, circuit active power and reactive power, transformer and reactive power, bus etc. to utilize the realtime telemetry remote signalling amount of electric power system; The remote signalling amount mainly comprises position, load tap changer gear of switch and disconnecting link etc.), estimate to determine the network topology and the operation of power networks state of electric power system by realtime status, calculate in real time definite electric network swim by online trend of carrying out and distribute, electric network model of Que Dinging and trend distribute and are referred to as electric network tide model like this.
For an electrical network that contains N node, electric network tide model can be described below:
If the injection active power of node i and injection reactive power are known quantity in the electrical network, claim that then this node is the PQ node, and constraint equation arranged:
If the injection active power and the voltage magnitude of node i are known quantity, claim that then this node is the PV node, and constraint equation arranged:
If the node voltage amplitude and the phase angle of node i are known quantity, claim that then this node is a V θ node, and constraint equation arranged:
P in the above formula
_{i} ^{Sp}, Q
_{i} ^{Sp}, V
_{i} ^{Sp}And θ
_{i} ^{Sp}Be respectively the given injection active power of node i, injection reactive power, node voltage amplitude and node voltage phase angle; G
_{Ij}, B
_{Ij}Be the real part and the imaginary part of the capable j column element of node admittance matrix i; θ
_{Ij}Be the node voltage phase angle difference between node i and the node j: V
_{i}, V
_{j}It is the voltage magnitude of node i and node j.
In the tide model of the electrical network of N node, each node types is one of abovementioned 3 category nodes, only establishes a V θ node usually.Conventional trend is calculated and is: according to the setpoint in left side in the known above formula, the V on right side in the computing formula
_{i}And θ
_{i}, i=1,2 ..., N.
For multiregion electric network, each subregion electrical network can carry out trend in order to last method and calculate.But, in setting up the electric network tide model process in real time, owing to be difficult for obtaining the tide model of own administrative area with external power grid, in order to take into account of the influence of overall electrical network to own administrative area electric network swim, adopt the readymade external electrical network Equivalent Model of offline to replace its actual external electrical network at present, this way can't adapt to the variation of external electrical network, causes the safety analysis computational accuracy poor, brings difficulty for the practicability of EMS.
The beautiful power outage that strengthens impelled dispatching of power netwoks personnel selfexamination in 2003 8.14, was confined to own administrative area this computation schema of electric network model and can not satisfies the requirement that the electrical network safety on line is analyzed, and needed the overall electric network tide model of online foundation.This model not only should be able to correctly be represented the realtime status of administrative area electrical network, should be able to correctly take into account the influence of nonadministrative area electrical network to the administrative area electrical network simultaneously.Only on overall electric network tide model basis, carry out the power grid security analysis, could effectively guarantee the correctness of safety on line analysis result, thereby ensure electric power netting safe running.
Here need to prove, relation between administrative area electrical network (being also referred to as internal electric network) and the nonadministrative area electrical network (being also referred to as external electrical network) mainly contains following 2 kinds at present: laterally the sane level relation is (such as the relation between two interconnected provincial power networks, one of them provincial power network carries out need considering the influence of another provincial power network to it when safety analysis is calculated), vertically relationship between superior and subordinate (such as large regional grid and should each provincial power network of zone subordinate between relation, because the fine degree difference of both modelings, large regional grid needs to consider the detailed model of this each provincial power network of zone in carrying out safety analysis).To above two kinds of situations, the method for merge partial electric grid, setting up overall electric network tide model is identical.Purpose of the present invention is exactly the tide model that will propose the online merging multiregion electric network of a kind of method, guarantees to merge the trend and the administrative area realtime power network state consistency of electrical network part in administrative area in the overall electric network tide model in back simultaneously.In order to narrate conveniently, below represent administrative area electrical network (or claiming internal electric network) with higher level's electrical network, represent nonadministrative area electrical network (or claiming external electrical network) with subordinate's electrical network.
Summary of the invention
The objective of the invention is can not the overall electric network model problem of online correct foundation in order to solve traditional E MS, also in order to satisfy EMS of new generation to setting up the new demand of overall electric network tide model, the online merging method of multiregion electric network tide model is proposed, in electric power system, utilize the highspeed data communication net that has built up between the current grid dispatching center, the electric network model that carries out between administrative area electrical network and the nonadministrative area electrical network is mutual, and utilize coupling trend method that nonarea under one's jurisdiction electric network tide model is regulated, realize that the no disturbance of tide model merges between administrative area electrical network and the nonadministrative area electrical network, state and its realtime status of administrative area electrical network part are in full accord in the overall electric network tide model after guaranteeing to merge.On overall electric network tide model basis, carry out the power grid security analysis, can effectively guarantee the correctness of analysis result, thereby ensure electric power netting safe running.
The online merging method of the multiregion electric network tide model that the present invention proposes is characterized in that this method may further comprise the steps:
(1) higher level's grid dispatching center is according to administration realtime power network data, utilize method for estimating state to generate one's respective area internal electric network model and trend distribution automatically, set up the external electrical network Equivalent Model that is attached thereto simultaneously, described inside and outside model is by the higher level of the continuous formation of in esse interconnection between the two electric network model;
(2) each subregion grid dispatching center of subordinate is according to administration realtime power network data, utilize method for estimating state to generate one's respective area internal electric network model and trend distribution automatically, set up the external electrical network Equivalent Model that is attached thereto simultaneously, described inside and outside model forms subordinate electric network model by in esse interconnection is continuous between the two, and by wide area network described subordinate electric network model and flow data sent to higher level's grid dispatching center;
(3) higher level's electric network model and each subordinate's electric network model are docked by the interconnection between them, upper and lower level electric network model is merged into the overall electric network model that calculates usefulness;
(4) utilize the trend matching process to regulate the trend of subordinate's electric network model, make the trend of the superior and the subordinate's electric network model mate, obtain overall electric network tide model at last.
Characteristics of the present invention:
The tide model of this rank, one's respective area electrical network is set up and safeguarded to each subregion grid dispatching center of higher level's grid dispatching center and subordinate respectively, to meet the existing management mode of classification, subregion scheduling.Subordinate's grid dispatching center need send to the higher level control centre with this area under one's jurisdiction realtime power network tide model, after the higher level collects the control centre tide model of each subregion electrical network of subordinate, each subregion electric network model of each subordinate and higher level's electric network model are merged, form overall electric network model.In merging process, adjust subordinate's electric network swim, make subordinate's electric network swim and higher level's electric network swim mate, thereby the no disturbance that guarantees both electric network tide models merges, merge higher level's electrical network part and higher level's realtime power network state consistency in the overall electric network tide model in back.Whole process is finished automatically, is fit to safety on line analysis, scheduling decision and the analog simulation of extensive interconnected network.
Method of the present invention has following advantage:
1, grid dispatching center at different levels only need be safeguarded electric network tide model separately, do not increase maintenance workload;
2, mate by gathering of each subregion electric network model with trend, automatically realize the online merging of multiregion electric network tide model, higher level's electrical network part is consistent with actual conditions in the overall electric network tide model in merging back, guarantee interconnected network safety analysis, control decision and analog simulation precision, thereby guaranteed the safe operation of electrical network.
Description of drawings
Fig. 1 is the online merging method schematic diagram of electric network tide model of the present invention; A) being higher level's electric network tide model wherein, b) is subordinate's electric network tide model, c) is the overall electric network tide model after merging.
Fig. 2 is that electric network model is torn schematic diagram, a) is electrical network before tearing wherein, b) is the electrical network after tearing.
Among Fig. 21 is the internal electric network of subordinate's electrical network;
Among Fig. 22 is internal electric network of higher level's electrical network;
Among Fig. 23 is the boundary nodes after tearing.
Embodiment
The online merging method that the present invention proposes the multiregion electric network tide model in the electric power system reaches embodiment in conjunction with the accompanying drawings and is described in detail as follows:
Specific implementation method of the present invention may further comprise the steps as depicted in figs. 1 and 2:
Step 1, as Fig. 1 b) shown in, higher level's grid dispatching center is according to administration realtime power network data, utilize method for estimating state to generate one's respective area internal electric network model 6 and trend distribution automatically, electric network model 6 and trend distribute and comprise topological structure of electric, the active power of grid generator, reactive power or machine end busbar voltage, the active power of network load and reactive power, the active power of electrical network busbar voltage and circuit and reactive power, and the active power of transformer, reactive power, load tap changer gear etc.; Set up the external electrical network Equivalent Model 4 that is attached thereto simultaneously, described inside and outside model passes through the in esse interconnection 5 continuous higher level of formation electric network models between the two;
Step 2, as Fig. 1 a) shown in, each subregion grid dispatching center of subordinate is according to administration realtime power network data, utilize method for estimating state to generate one's respective area internal electric network model 1 and trend distribution automatically, electric network model 1 and trend distribute and comprise topological structure of electric, the active power of grid generator, reactive power or machine end busbar voltage, the active power of network load and reactive power, the active power of electrical network busbar voltage and circuit and reactive power, and the active power of transformer, reactive power, load tap changer gear etc., set up the external electrical network Equivalent Model 3 that is attached thereto simultaneously, in described, external model links to each other by in esse interconnection 2 between the two and forms subordinate's electric network model, and by wide area network described subordinate electric network model and flow data is sent to higher level's grid dispatching center;
Step 3, higher level's grid dispatching center receives subordinate's electric network model and the flow data that abovementioned each subregion grid dispatching center is sent, at first reject external network equivalent model (dotted line) part in the electric network model separately, then, described higher level's electric network model and each subordinate's electric network model are docked (interconnections 2 by the interconnection between them 5, same group of interconnection of 5 expressions, here adopt higher level's grid side 5), just can with on, subordinate's electric network model is merged into the overall electric network model that calculates usefulness, as Fig. 1 c) shown in, as can be seen this overall electric network model with Fig. 1 a) and Fig. 1 b) in model remove respectively behind the dotted portion zone to be docking together and obtain by interconnection.
Step 4, utilize the trend matching process to regulate the trend of subordinate's electric network model, make the trend coupling of the superior and the subordinate's electric network model to obtain overall electric network tide model at last.
Because each subordinate's electrical network is uploaded in the overall electrical network electric network tide model and higher level's electric network tide model are inconsistent in time, therefore there is deviation in the corresponding value at identical with the higher level's electrical network boundary node with boundary node voltage of the interconnection power place of subordinate's electric network tide model, if this deviation is not adjusted into zero, both models just can have disturbance and merge.If merge, the state of higher level's electrical network part will change in the overall electrical network, not satisfy the calculation requirement of higher level's electrical network.The present invention proposes the trend matching process and eliminates this deviation, this method keeps higher level's electric network swim constant, the trend of revising subordinate's electrical network distributes, make interconnection active power between each subordinate's electrical network and the higher level's electrical network and reactive power and border busbar voltage with merge before the corresponding amount of higher level's electrical network consistent, can realize the no disturbance butt joint of the superior and the subordinate's electric network tide model like this.In the overall electric network tide model in butt joint back, the trend of higher level's electrical network part is constant.Electric power system can be carried out the overall electric network swim distribution that overall electric network swim calculates on this model, wherein the state of higher level's electrical network part is in full accord with the preceding higher level's realtime power network state of merging.This overall electric network tide model can offer EMS (EMS) and carry out electrical network safety on line analytical applications, to guarantee the safe operation of electrical network.
Trend coupling specific implementation method of the present invention as shown in Figure 2, tear by node, the end node of the close subordinate of the interconnection between the superior and the subordinate's electrical network grid side is defined as boundary node, it is split into two dummy nodes, what be connected on subordinate's grid side is subordinate's electrical network boundary node, be connected on higher level's grid side for higher level's electrical network boundary node, thereby higher level's electrical network (as Fig. 2 b) 2) and subordinate's electrical network (as Fig. 2 b) 1) keep apart; Adjust subordinate's electric network swim, make the voltage at subordinate electrical network boundary node place and injecting power equate at the corresponding value in higher level's electrical network boundary node place, realize that the no disturbance of two electric network tide models merges with higher level's realtime power network trend.
The present invention regulates subordinate's electric network swim method two kinds of execution modes, is described in detail as follows respectively:
Embodiment 1: adopt and regulate meritorious power P of subordinate's electrical network internal node injection and reactive power Q, specifically may further comprise the steps:
(1) the i of subordinate's electrical network boundary node (as Fig. 2 b), j, k) be set to the given node of V θ, this node set is defined as S, its voltage magnitude V and voltage phase angle θ be arranged to higher level's electrical network boundary node in higher level's realtime power network trend that higher level's grid dispatching center provides corresponding value (Fig. 2 b) 3); Other nodes are set to the PQ node in subordinate's electrical network, and the given charge value of this PQ is set to original trend end value in subordinate's electrical network, carry out many V θ node trend of subordinate's electrical network and calculate, and obtain the meritorious power P of injection of subordinate's electrical network boundary node
_{S}And reactive power Q
_{S}, have deviation delta P between the corresponding value of higher level's electrical network boundary node in this value and the higher level's realtime power network trend
_{S}With Δ Q
_{S}
(2) according to the given node active power of described V θ departure Δ P
_{S}, calculate the active power regulated quantity Δ P of subordinate's network load or generating
_{M}, by regulating Δ P
_{M}To eliminate the electrical network boundary node active power departure Δ P of subordinate
_{S}M participates in the node set that active power is adjusted;
Utilize sensitivity method, set up the meritorious power deviation amount Δ P of injection of subordinate's electrical network boundary node
_{S}With the meritorious power adjustment Δ P of subordinate's electrical network internal node injection
_{M}Between sensitivity relation:
ΔP
_{S}＝D
_{SM}ΔP
_{M}
In the formula, D
_{SM}Be S * M rank (dimensionless) sensitivity matrix, when its j column element represents that M gathers unit quantity of injection active power change of j adjustable node, the size of the meritorious power variation of injection of each node in the S set; Its computational methods are as follows:
Represent that B ' matrix (N * N rank) that quick decomposition method trend is calculated in (meritoriousthe phase angle iteration) keeps S and M collector node, other do not participate in concentrated matrix behind the node that active power adjusts cancellation, and subscript S and M represent respectively
In with the relevant part of S and M set;
Utilize pseudoinverse technique to calculate subordinate's electrical network internal node and inject meritorious power adjustments amount:
In the formula
Be Δ P
_{M}With Δ P
_{S}Between the distribution factor matrix, calculate with following formula
The transposition of the subscript T representing matrix of matrix in the formula; Use Δ P at last
_{M}The injection active power of revising the M collector node reduces boundary node and injects meritorious power deviation Δ P
_{S}, to realize the meritorious trend coupling of the given node of V θ;
(3) inject reactive power departure Δ Q according to the given node of abovementioned V θ
_{S}, calculate subordinate's grid generator node (can comprise load bus) and inject reactive power regulated quantity Δ Q
_{R}, regulate Δ Q
_{R}Inject reactive power departure Δ Q to eliminate the given node of V θ
_{S}R participates in injecting the node set that reactive power is adjusted, and the R set can be identical or different with the M set in (2) step; Specific practice is:
Set up the given node of S set V θ and inject reactive power variation delta Q
_{S}Inject reactive power adjustment amount Δ Q with the R collector node
_{R}Between sensitivity relation:
ΔQ
_{S}＝S
_{SR}ΔQ
_{R}
In the formula, S
_{SR}Be S * R rank (dimensionless) sensitivity matrix, when its j column element represented that R gathers unit quantity of j node injection reactive power change, each node injected the size of reactive power variable quantity in the S set; And have
" subscript S and R represent respectively for matrix (N * N rank) reservation S and R collector node, the concentrated matrix after other node set that do not participate in adjusting of cancellation to represent quick decomposition method trend to calculate the B in (idlethe voltage iteration)
In with the relevant part of S and R set;
Utilize pseudoinverse technique, computing node injects the reactive power regulated quantity:
In the formula
Be Δ Q
_{R}With Δ Q
_{S}Between the distribution factor matrix, calculate with following formula
The transposition of the subscript T representing matrix of matrix in the formula; Use Δ Q at last
_{R}The injection reactive power of modified R collector node reduces the given node of V θ and injects reactive power deviation delta Q
_{S}, to realize the coupling of the given node reactive power flow of V θ;
(4) (because power flow equation is nonlinear, the regulated quantity that can not utilize abovementioned linearization equations to calculate is once eliminated deviation, so adjustment process is the process of an iterative computation.) revised the injection active power of M collector node and the injection reactive power of R collector node, again return (1) and do the calculating of many V θ node trend, in the calculating, the V θ setpoint of the electrical network boundary node S of subordinate set still uses the analog value of higher level's electrical network boundary node in higher level's realtime power network trend; Repeating step (1)(3) equate with the analog value of higher level's electrical network boundary node in higher level's realtime power network trend up to the active power and the reactive power of subordinate's electrical network boundary node, i.e. deviation delta P
_{S}With Δ Q
_{S}Be zero; Two electric network swim couplings can not have them disturbance and merged this moment; In the overall electric network tide model after the merging, corresponding to the trend of higher level's electrical network part with merge preceding identical.
Embodiment 2: regulate the method that subordinate's electrical network internal node injects meritorious power P and node voltage amplitude V.Specific practice is:
(1), the i of subordinate's electrical network boundary node (as Fig. 2 b), j, k) be set to the given node of Q θ, this node set is defined as S, its node inject reactive power Q and node voltage phase angle theta be arranged to higher level's realtime power network trend higher level electrical network boundary node analog value (Fig. 2 b) 3), and adopt the following formula constraint equation to calculate trend:
θ in the formula
_{i} ^{Sp}And Q
_{i} ^{Sp}Be given voltage phase angle of the electrical network boundary node i of subordinate and injection reactive power, G
_{Ij}, B
_{Ij}Be respectively the real part and the imaginary part of the capable j column element of node admittance matrix i; θ
_{Ij}It is the node voltage phase angle difference between node i and the node j; V
_{i}, V
_{j}It is the voltage magnitude of node i and node j.The node that needs in subordinate's electrical network to participate in the voltage adjustment is set to the PV node, other nodes are set to the PQ node, the given charge value of its PQ and PV is set to original trend end value in subordinate's electrical network, carry out many Q θ node trend of subordinate's electrical network and calculate, obtain the meritorious power P of injection of subordinate's electrical network boundary node
_{S}With node voltage amplitude V
_{S}, have deviation delta P between the corresponding value of higher level's electrical network boundary node in this value and the higher level's realtime power network trend
_{S}With Δ V
_{S}
(2), inject meritorious power deviation amount Δ P according to abovementioned border Q θ node
_{S}, calculate the active power regulated quantity Δ P of subordinate's network load or generating with following method
_{M}, by regulating Δ P
_{M}To eliminate the electrical network boundary node active power departure Δ P of subordinate
_{S}M participates in the node set that active power is adjusted;
Utilize sensitivity method, set up the meritorious power deviation amount Δ P of injection of subordinate's electrical network boundary node
_{S}With the meritorious power adjustment Δ P of subordinate's electrical network internal node injection
_{M}Between sensitivity relation:
ΔP
_{S}＝D
_{SM}ΔP
_{M}
In the formula, D
_{SM}Be S * M rank (dimensionless) sensitivity matrix, when its j column element represents that M gathers unit quantity of injection active power change of j adjustable node, the size of the meritorious power variation of injection of each node in the S set; Its computational methods are as follows:
Represent that B ' matrix (N * N rank) that quick decomposition method trend is calculated in (meritoriousthe phase angle iteration) keeps S and M collector node, other do not participate in concentrated matrix behind the node that active power adjusts cancellation, and subscript S and M represent respectively
In with the relevant part of S and M set;
Utilize pseudoinverse technique to calculate subordinate's electrical network internal node and inject meritorious power adjustments amount:
In the formula
Be Δ P
_{M}With Δ P
_{S}Between the distribution factor matrix, calculate with following formula
The transposition of the subscript T representing matrix of matrix in the formula; Use Δ P at last
_{M}The injection active power of revising the M collector node reduces boundary node and injects meritorious power deviation Δ P
_{S}, to realize the meritorious trend coupling of the given node of Q θ;
(3), according to abovementioned border Q θ node voltage amplitude departure Δ V
_{S}, calculate the voltage magnitude regulated quantity Δ V that subordinate's electrical network participates in the node of voltage adjustment
_{R}, regulate Δ V
_{R}To eliminate boundary node voltage magnitude departure Δ V
_{S}R participates in the node set that voltage is adjusted, and the R set can be identical or different with the M set in (2) step; Specific practice is:
Set up S collector node voltage magnitude variation delta V
_{S}With R collector node voltage magnitude regulated quantity Δ V
_{R}Between sensitivity relation:
ΔV
_{S}＝R
_{SR}ΔV
_{R}
In the following formula, R
_{SR}Be S * R rank (dimensionless) sensitivity matrix, when its j column element represents that R gathers unit quantity of voltage magnitude change of j the adjustable node of voltage, the size of each Q θ node voltage amplitude variable quantity in the S set; And have
" matrix (N * N rank) cancellation does not participate in the concentrated matrix behind the node that voltage adjusts, and subscript S and R represent respectively to represent that quick decomposition method trend calculates the B in (idlethe voltage iteration)
In with the relevant part of S and R set;
Utilize pseudoinverse technique, the regulated quantity of computing node voltage magnitude:
In the formula
Be Δ V
_{R}With Δ V
_{S}Between sensitivity, calculate with following formula
The transposition of the subscript T representing matrix of matrix in the formula; Use Δ V
_{R}The voltage magnitude of modified R set generator node is so that Δ V
_{S}Reduce, realize the coupling of boundary node voltage magnitude at last.
(4) (because power flow equation is nonlinear, the regulated quantity that can not utilize linearization equations to calculate is once eliminated deviation, so adjustment process is the process of an iterative computation.) revised the meritorious injecting power and the R collector node voltage magnitude of M collector node, again return (1) and do the calculating of many Q θ node trend, in the calculating, the Q θ setpoint of the electrical network boundary node S of subordinate set uses still that higher level's electrical network boundary node is worth accordingly in higher level's realtime power network trend; Repeating step (1)(3) equate with the corresponding value of higher level's electrical network boundary node in higher level's realtime power network trend up to the active power and the voltage magnitude of subordinate's electrical network boundary node, i.e. deviation delta P
_{S}With Δ V
_{S}Be zero; Two electric network swim couplings can not have them disturbance and merged this moment; In the overall electric network tide model after the merging, corresponding to the trend of higher level's electrical network part with merge preceding identical.
Claims (4)
1. the online merging method of a multiregion electric network tide model is characterized in that this method may further comprise the steps:
(1) higher level's grid dispatching center is according to administration realtime power network data, utilize method for estimating state to generate one's respective area internal electric network model and trend distribution automatically, set up the external electrical network Equivalent Model that is attached thereto simultaneously, described inside and outside model is by the higher level of the continuous formation of in esse interconnection between the two electric network model;
(2) each subregion grid dispatching center of subordinate is according to administration realtime power network data, utilize method for estimating state to generate one's respective area internal electric network model and trend distribution automatically, set up the external electrical network Equivalent Model that is attached thereto simultaneously, described inside and outside model forms subordinate electric network model by in esse interconnection is continuous between the two, and by wide area network described subordinate electric network model and flow data sent to higher level's grid dispatching center;
(3) higher level's electric network model and each subordinate's electric network model are docked by the interconnection between them, upper and lower level electric network model is merged into the overall electric network model that calculates usefulness;
(4) utilize the trend matching process to regulate the trend of subordinate's electric network model, make the trend of the superior and the subordinate's electric network model mate, obtain overall electric network tide model at last.
2. the method for claim 1 is characterized in that, the trend matching process in the described step (4) is:
Tear by node, the end node of the close subordinate of the interconnection between the superior and the subordinate's electrical network grid side is defined as boundary node, and it is split into inside and outside two dummy nodes, what be connected on subordinate's grid side is subordinate's electrical network boundary node, what be connected on higher level's grid side is higher level's electrical network boundary node, thereby higher level's electrical network and subordinate's electrical network are kept apart; Regulating subordinate's electrical network internal node injects meritorious power P and reactive power Q or regulates meritorious power P of subordinate's electrical network internal node injection and voltage magnitude V, make the voltage and the injecting power at subordinate electrical network boundary node place equate, merge with the no disturbance that realizes two electric network tide models with the corresponding value at higher level's realtime power network trend higher level electrical network boundary node place.
3. method as claimed in claim 2 is characterized in that, described adjusting meritorious power P of subordinate's electrical network internal node injection and reactive power Q specifically may further comprise the steps:
(1) subordinate's electrical network boundary node is set to the given node of V θ, and this node set is defined as S, and its voltage magnitude V and voltage phase angle θ are arranged to the corresponding value of higher level's electrical network boundary node in higher level's realtime power network trend that higher level's grid dispatching center provides; Other nodes are set to the PQ node in subordinate's electrical network, and the given charge value of this PQ is set to original trend end value in subordinate's electrical network, carry out many V θ node trend of subordinate's electrical network and calculate, and obtain the meritorious power P of injection at subordinate electrical network boundary node place
_{S}And reactive power Q
_{S}, have deviation delta P between the corresponding value of higher level's electrical network boundary node in this value and the higher level's realtime power network trend
_{S}With Δ Q
_{S}
(2) according to the given node active power of described V θ departure Δ P
_{S}, utilize following formula to calculate the inner node that participates in the active power adjusting of subordinate's electrical network and inject meritorious power adjustments amount Δ P
_{M}, M participates in the node set that active power is adjusted:
In the formula
Be Δ P
_{M}With Δ P
_{S}Between the distribution factor matrix, calculate with following formula
In the formula, D
_{SM}Be Δ P
_{S}With Δ P
_{M}Between S * M rank dimensionless sensitivity matrix, that is: Δ P
_{S}=Δ D
_{SM}Δ P
_{M}The transposition of the subscript T representing matrix of matrix; Use Δ P
_{M}The injection active power of revising the M collector node reduces the given node of V θ and injects meritorious power deviation Δ P
_{S}
(3) inject reactive power departure Δ Q according to the given node of described V θ
_{S}, utilize following formula to calculate the inner injection reactive power regulated quantity Δ Q that participates in reactive power adjustment node of subordinate's electrical network
_{R}, this adjustment node set is R:
In the formula
Be Δ Q
_{R}With Δ Q
_{S}Between the distribution factor matrix, calculate with following formula
S in the formula
_{SR}Be Δ Q
_{S}With Δ Q
_{R}Between the S * nondimensional sensitivity matrix in R rank, that is: Δ Q
_{S}=S
_{SR}Δ Q
_{R}The transposition of the subscript T representing matrix of matrix; Use Δ Q
_{R}The injection reactive power of modified R collector node reduces the given node of V θ and injects reactive power deviation delta Q
_{S}
(4) return step (1)step (3) again, the node up to subordinate electrical network boundary node place injects meritorious power P and reactive power Q equates with the value at higher level's realtime power network trend higher level electrical network boundary node place, i.e. deviation delta P
_{S}With Δ Q
_{S}Be zero; Higher level, two electric network swims of subordinate mate at this moment, to realize that two electric network models of upper and lower level do not have disturbance and are merged into overall electric network tide model.
4. method as claimed in claim 2 is characterized in that, described adjusting subordinate's electric network swim employing adjusting meritorious power P of subordinate's electrical network internal node injection and voltage magnitude V specifically may further comprise the steps:
(1) subordinate's grid side boundary node is set to the given node of Q θ, this node set is defined as S, its node injection reactive power Q and node voltage phase angle theta are arranged to the analog value of higher level's realtime power network trend higher level electrical network boundary node, the node that needs in subordinate's electrical network to participate in the voltage adjustment is set to the PV node, R is this PV node set, other nodes are set to the PQ node, the given charge value of this PQ and PV is set to original trend end value in subordinate's electrical network, carry out many Q θ node trend of subordinate's electrical network and calculate, obtain the meritorious power P of injection of subordinate's electrical network boundary node
_{S}With node voltage amplitude V
_{S}, have deviation delta P between the corresponding value of higher level's electrical network boundary node in this value and the higher level's realtime power network trend
_{S}With Δ V
_{S}
(2) inject meritorious power deviation amount Δ P according to the given node of described Q θ
_{S}, utilize following formula to calculate the inner meritorious power adjustments amount Δ P of injection that participates in active power adjusting node of subordinate's electrical network
_{M}, M participates in the node set that active power is adjusted:
In the formula
Be Δ P
_{M}With Δ P
_{S}Between the distribution factor matrix, calculate with following formula
In the formula, D
_{SM}Be Δ P
_{S}With Δ P
_{M}Between S * M rank dimensionless sensitivity matrix, that is: Δ P
_{S}=D
_{SM}Δ P
_{M}The transposition of the subscript T representing matrix of matrix; Use Δ P
_{M}The injection active power of revising the M collector node reduces the given node of the electrical network Q θ of subordinate and injects meritorious power deviation Δ P
_{S}
(3), according to the given node voltage amplitude of described Q θ departure Δ V
_{S}, calculate the voltageregulation amount Δ V that subordinate's electrical network participates in the node of voltage adjustment with following formula
_{R}, R participates in the set of the node of voltage adjustment for subordinate's electrical network:
In the formula
Be Δ V
_{R}With Δ V
_{S}Between sensitivity, calculate with following formula
R in the formula
_{SR}Be S collector node voltage magnitude variation delta V
_{S}With R collector node voltage magnitude regulated quantity Δ V
_{R}Between sensitivity, that is: Δ V
_{S}=R
_{SR}Δ V
_{R}The transposition of the subscript T representing matrix of matrix; Use Δ V
_{R}Modified R collector node voltage magnitude reduces the deviation delta V of the given node voltage amplitude of the electrical network Q θ of subordinate
_{S}
(4) return step (1)step (3) again, the node up to subordinate electrical network boundary node place injects meritorious power P and node voltage amplitude V equates with the value of higher level's realtime power network trend higher level electrical network boundary node, i.e. deviation delta P
_{S}With Δ V
_{S}Be zero; Higher level, two electric network swims couplings of subordinate do not have disturbance and are merged into overall electric network tide model to realize higher level, two electric network models of subordinate at this moment.
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

CNB200710176083XA CN100570983C (en)  20071019  20071019  The online merging method of multiregion electric network tide model 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

CNB200710176083XA CN100570983C (en)  20071019  20071019  The online merging method of multiregion electric network tide model 
Publications (2)
Publication Number  Publication Date 

CN101154813A true CN101154813A (en)  20080402 
CN100570983C CN100570983C (en)  20091216 
Family
ID=39256321
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CNB200710176083XA CN100570983C (en)  20071019  20071019  The online merging method of multiregion electric network tide model 
Country Status (1)
Country  Link 

CN (1)  CN100570983C (en) 
Cited By (7)
Publication number  Priority date  Publication date  Assignee  Title 

CN101877500A (en) *  20090430  20101103  北京科东电力控制系统有限责任公司  Interconnected power system oriented hierachical decomposition spacetime cooperative modeling method 
CN102044871A (en) *  20101022  20110504  浙江省电力公司  Method for matching power flow of superior and inferior power grids in dispatcher training simulators (DTSs) 
CN101540505B (en) *  20090109  20120523  南京南瑞继保电气有限公司  Building method of multistage multiregion interconnected network data model 
CN104167732A (en) *  20140825  20141126  中国南方电网有限责任公司  Grid equivalent method based on phase angle difference of call lines 
CN104319772A (en) *  20141021  20150128  深圳大学  Method for symmetrically acquiring loss power sensitivity from ACDC hybrid power system 
CN104362664A (en) *  20140728  20150218  浙江工业大学  Grid connection method of mediumvoltage microgrid system 
CN104657618A (en) *  20150306  20150527  南车株洲电力机车研究所有限公司  Splitting fault modeling method for traction substations and splitting voltage determining method 

2007
 20071019 CN CNB200710176083XA patent/CN100570983C/en not_active IP Right Cessation
Cited By (9)
Publication number  Priority date  Publication date  Assignee  Title 

CN101540505B (en) *  20090109  20120523  南京南瑞继保电气有限公司  Building method of multistage multiregion interconnected network data model 
CN101877500A (en) *  20090430  20101103  北京科东电力控制系统有限责任公司  Interconnected power system oriented hierachical decomposition spacetime cooperative modeling method 
CN101877500B (en) *  20090430  20140625  北京科东电力控制系统有限责任公司  Interconnected power system oriented hierachical decomposition spacetime cooperative modeling method 
CN102044871A (en) *  20101022  20110504  浙江省电力公司  Method for matching power flow of superior and inferior power grids in dispatcher training simulators (DTSs) 
CN104362664A (en) *  20140728  20150218  浙江工业大学  Grid connection method of mediumvoltage microgrid system 
CN104167732A (en) *  20140825  20141126  中国南方电网有限责任公司  Grid equivalent method based on phase angle difference of call lines 
CN104319772A (en) *  20141021  20150128  深圳大学  Method for symmetrically acquiring loss power sensitivity from ACDC hybrid power system 
CN104657618A (en) *  20150306  20150527  南车株洲电力机车研究所有限公司  Splitting fault modeling method for traction substations and splitting voltage determining method 
CN104657618B (en) *  20150306  20180105  南车株洲电力机车研究所有限公司  A kind of traction substation offtheline fault modeling method and solution column voltage determine method 
Also Published As
Publication number  Publication date 

CN100570983C (en)  20091216 
Similar Documents
Publication  Publication Date  Title 

Bevrani et al.  Intelligent automatic generation control  
Bidram et al.  Distributed control systems for smallscale power networks: Using multiagent cooperative control theory  
Olivares et al.  Trends in microgrid control  
Dufour et al.  On the use of realtime simulation technology in smart grid research and development  
Hu et al.  Improving frequency stability based on distributed control of multiple load aggregators  
Zhang et al.  Further results on delaydependent stability of multiarea load frequency control  
CN104330980B (en)  A kind of microcapacitance sensor emulation test system based on RT LAB  
Bejestani et al.  A hierarchical transactive control architecture for renewables integration in smart grids: Analytical modeling and stability  
Hcbson  Network constrained reactive power control using linear programming  
Haileselassie  Control, dynamics and operation of multiterminal VSCHVDC transmission systems  
EP2182207B1 (en)  Widearea transmission control of windfarms  
Ray et al.  Widearea signalbased optimalneurocontroller for a UPFC  
CN102081867B (en)  Alldigital highvoltage direct current transmission training stimulation system  
CN100461579C (en)  Method for controlling coordination voltage of regional power grid and provincial power grid  
Makarov et al.  PMUbased widearea security assessment: concept, method, and implementation  
Chakrabortty  Widearea damping control of power systems using dynamic clustering and TCSCbased redesigns  
Tyagi et al.  A decentralized automatic generation control scheme for competitive electricity markets  
Van Hertem et al.  Usefulness of DC power flow for active power flow analysis with flow controlling devices  
CN105574238B (en)  A kind of section real power control analysis of strategies method for the assessment of online stability margin  
Ilic et al.  Modeling future cyberphysical energy systems  
CN103715721B (en)  The realtime peak regulating method of alternating currentdirect current mixing power grid joint  
Biegel et al.  Distributed lowcomplexity controller for wind power plant in derated operation  
Shoults et al.  Multiarea adaptive LFC developed for a comprehensive AGC simulator  
Muljadi et al.  Equivalencing the collector system of a large wind power plant  
CN101552474B (en)  AGC hierarchical coordinative control method based on security constraint of stable cross section 
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  
CF01  Termination of patent right due to nonpayment of annual fee  
CF01  Termination of patent right due to nonpayment of annual fee 
Granted publication date: 20091216 Termination date: 20191019 