CN101154813A - On-line combination method for multi-region electric network tide model - Google Patents

On-line combination method for multi-region electric network tide model Download PDF

Info

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
Application number
CNA200710176083XA
Other languages
Chinese (zh)
Other versions
CN100570983C (en
Inventor
张伯明
吴文传
孙宏斌
蒋维勇
Original Assignee
清华大学
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 清华大学 filed Critical 清华大学
Priority to CNB200710176083XA priority Critical patent/CN100570983C/en
Publication of CN101154813A publication Critical patent/CN101154813A/en
Application granted granted Critical
Publication of CN100570983C publication Critical patent/CN100570983C/en

Links

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 real-time 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

The online merging method of multi-region electric network tide model
Technical field
The invention belongs to the dispatching automation of electric power systems technical field, relate to multi-region 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 real-time 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 real-time 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 real-time 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:
P i sp = V i Σ i ∈ j V j ( G ij cos θ ij + B ij sin θ ij ) Q i sp = V i Σ i ∈ j V j ( G ij sin θ ij - B ij cos θ ij )
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:
P i sp = V i Σ i ∈ j V j ( G ij cos θ ij + B ij sin θ ij ) V i sp = V i
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:
V i sp = V i θ i sp = θ i
P in the above formula i Sp, Q i Sp, V i SpAnd θ i SpBe respectively the given injection active power of node i, injection reactive power, node voltage amplitude and node voltage phase angle; G Ij, B IjBe the real part and the imaginary part of the capable j column element of node admittance matrix i; θ IjBe the node voltage phase angle difference between node i and the node j: V i, V jIt 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 above-mentioned 3 category nodes, only establishes a V θ node usually.Conventional trend is calculated and is: according to the set-point in left side in the known above formula, the V on right side in the computing formula iAnd θ i, i=1,2 ..., N.
For multi-region 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 ready-made external electrical network Equivalent Model of off-line 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 self-examination 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 real-time status of administrative area electrical network, should be able to correctly take into account the influence of non-administrative 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 non-administrative 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 multi-region electric network of a kind of method, guarantees to merge the trend and the administrative area real-time 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 non-administrative 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 multi-region electric network tide model is proposed, in electric power system, utilize the high-speed 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 non-administrative area electrical network is mutual, and utilize coupling trend method that non-area 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 non-administrative area electrical network, state and its real-time 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 multi-region 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 real-time 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 real-time 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 real-time 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 real-time 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 multi-region 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 multi-region 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 real-time 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 real-time 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 above-mentioned 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 real-time 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 real-time 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 real-time 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 SAnd 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 real-time power network trend SWith Δ 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 MTo eliminate the electrical network boundary node active power departure Δ P of subordinate SM 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 SWith the meritorious power adjustment Δ P of subordinate's electrical network internal node injection MBetween sensitivity relation:
ΔP S=D SMΔP M
In the formula, D SMBe 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:
D SM = B ~ SM ′ B ~ MM ′ - 1
Represent that B ' matrix (N * N rank) that quick decomposition method trend is calculated in (meritorious-the 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:
Δ P M = D ~ MS Δ P S
In the formula Be Δ P MWith Δ P SBetween 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 MThe 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 above-mentioned V θ S, calculate subordinate's grid generator node (can comprise load bus) and inject reactive power regulated quantity Δ Q R, regulate Δ Q RInject reactive power departure Δ Q to eliminate the given node of V θ SR 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 SInject reactive power adjustment amount Δ Q with the R collector node RBetween sensitivity relation:
ΔQ S=S SRΔQ R
In the formula, S SRBe 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
S SR = B ~ SR ′ ′ B ~ RR ′ ′ - 1
" 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 (idle-the voltage iteration) In with the relevant part of S and R set;
Utilize pseudoinverse technique, computing node injects the reactive power regulated quantity:
Δ Q R = S ~ RS Δ Q S
In the formula Be Δ Q RWith Δ Q SBetween 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 RThe 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 non-linear, the regulated quantity that can not utilize above-mentioned 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 θ set-point 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 real-time power network trend; Repeating step (1)-(3) equate with the analog value of higher level's electrical network boundary node in higher level's real-time power network trend up to the active power and the reactive power of subordinate's electrical network boundary node, i.e. deviation delta P SWith Δ Q SBe 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 real-time 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:
θ i sp = θ i Q i sp = V i Σ i ∈ j V j ( G ij sin θ ij - B ij cos θ ij )
θ in the formula i SpAnd Q i SpBe given voltage phase angle of the electrical network boundary node i of subordinate and injection reactive power, G Ij, B IjBe respectively the real part and the imaginary part of the capable j column element of node admittance matrix i; θ IjIt is the node voltage phase angle difference between node i and the node j; V i, V jIt 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 SWith 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 real-time power network trend SWith Δ V S
(2), inject meritorious power deviation amount Δ P according to above-mentioned 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 MTo eliminate the electrical network boundary node active power departure Δ P of subordinate SM 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 SWith the meritorious power adjustment Δ P of subordinate's electrical network internal node injection MBetween sensitivity relation:
ΔP S=D SMΔP M
In the formula, D SMBe 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:
D SM = B ~ SM ′ B ~ MM ′ - 1
Represent that B ' matrix (N * N rank) that quick decomposition method trend is calculated in (meritorious-the 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:
Δ P M = D ~ MS Δ P S
In the formula Be Δ P MWith Δ P SBetween 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 MThe 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 above-mentioned 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 RTo eliminate boundary node voltage magnitude departure Δ V SR 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 SWith R collector node voltage magnitude regulated quantity Δ V RBetween sensitivity relation:
ΔV S=R SRΔV R
In the following formula, R SRBe 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
R SR = - B ~ SS ′ ′ - 1 B ~ SR ′ ′
" 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 (idle-the voltage iteration) In with the relevant part of S and R set;
Utilize pseudoinverse technique, the regulated quantity of computing node voltage magnitude:
Δ V R = R ~ RS Δ V S
In the formula Be Δ V RWith Δ V SBetween sensitivity, calculate with following formula
The transposition of the subscript T representing matrix of matrix in the formula; Use Δ V RThe voltage magnitude of modified R set generator node is so that Δ V SReduce, realize the coupling of boundary node voltage magnitude at last.
(4) (because power flow equation is non-linear, 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 θ set-point 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 real-time power network trend; Repeating step (1)-(3) equate with the corresponding value of higher level's electrical network boundary node in higher level's real-time power network trend up to the active power and the voltage magnitude of subordinate's electrical network boundary node, i.e. deviation delta P SWith Δ V SBe 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 multi-region 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 real-time 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 real-time 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 real-time 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 real-time 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 SAnd 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 real-time power network trend SWith Δ 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:
Δ P M = D ~ MS Δ P S
In the formula Be Δ P MWith Δ P SBetween the distribution factor matrix, calculate with following formula
In the formula, D SMBe Δ P SWith Δ P MBetween S * M rank dimensionless sensitivity matrix, that is: Δ P S=Δ D SMΔ P MThe transposition of the subscript T representing matrix of matrix; Use Δ P MThe 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:
Δ Q R = S ~ RS Δ Q S
In the formula Be Δ Q RWith Δ Q SBetween the distribution factor matrix, calculate with following formula
S in the formula SRBe Δ Q SWith Δ Q RBetween the S * nondimensional sensitivity matrix in R rank, that is: Δ Q S=S SRΔ Q RThe transposition of the subscript T representing matrix of matrix; Use Δ Q RThe 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 real-time power network trend higher level electrical network boundary node place, i.e. deviation delta P SWith Δ Q SBe 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 real-time 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 SWith 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 real-time power network trend SWith Δ 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:
Δ P M = D ~ MS Δ P S
In the formula Be Δ P MWith Δ P SBetween the distribution factor matrix, calculate with following formula
In the formula, D SMBe Δ P SWith Δ P MBetween S * M rank dimensionless sensitivity matrix, that is: Δ P S=D SMΔ P MThe transposition of the subscript T representing matrix of matrix; Use Δ P MThe 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 voltage-regulation 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:
Δ V R = R ~ RS Δ V S
In the formula Be Δ V RWith Δ V SBetween sensitivity, calculate with following formula
R in the formula SRBe S collector node voltage magnitude variation delta V SWith R collector node voltage magnitude regulated quantity Δ V RBetween sensitivity, that is: Δ V S=R SRΔ V RThe transposition of the subscript T representing matrix of matrix; Use Δ V RModified 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 real-time power network trend higher level electrical network boundary node, i.e. deviation delta P SWith Δ V SBe 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.
CNB200710176083XA 2007-10-19 2007-10-19 The online merging method of multi-region electric network tide model CN100570983C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB200710176083XA CN100570983C (en) 2007-10-19 2007-10-19 The online merging method of multi-region electric network tide model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB200710176083XA CN100570983C (en) 2007-10-19 2007-10-19 The online merging method of multi-region electric network tide model

Publications (2)

Publication Number Publication Date
CN101154813A true CN101154813A (en) 2008-04-02
CN100570983C CN100570983C (en) 2009-12-16

Family

ID=39256321

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB200710176083XA CN100570983C (en) 2007-10-19 2007-10-19 The online merging method of multi-region electric network tide model

Country Status (1)

Country Link
CN (1) CN100570983C (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101877500A (en) * 2009-04-30 2010-11-03 北京科东电力控制系统有限责任公司 Interconnected power system oriented hierachical decomposition space-time cooperative modeling method
CN102044871A (en) * 2010-10-22 2011-05-04 浙江省电力公司 Method for matching power flow of superior and inferior power grids in dispatcher training simulators (DTSs)
CN101540505B (en) * 2009-01-09 2012-05-23 南京南瑞继保电气有限公司 Building method of multistage multi-region interconnected network data model
CN104167732A (en) * 2014-08-25 2014-11-26 中国南方电网有限责任公司 Grid equivalent method based on phase angle difference of call lines
CN104319772A (en) * 2014-10-21 2015-01-28 深圳大学 Method for symmetrically acquiring loss power sensitivity from AC-DC hybrid power system
CN104362664A (en) * 2014-07-28 2015-02-18 浙江工业大学 Grid connection method of medium-voltage microgrid system
CN104657618A (en) * 2015-03-06 2015-05-27 南车株洲电力机车研究所有限公司 Splitting fault modeling method for traction substations and splitting voltage determining method

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101540505B (en) * 2009-01-09 2012-05-23 南京南瑞继保电气有限公司 Building method of multistage multi-region interconnected network data model
CN101877500A (en) * 2009-04-30 2010-11-03 北京科东电力控制系统有限责任公司 Interconnected power system oriented hierachical decomposition space-time cooperative modeling method
CN101877500B (en) * 2009-04-30 2014-06-25 北京科东电力控制系统有限责任公司 Interconnected power system oriented hierachical decomposition space-time cooperative modeling method
CN102044871A (en) * 2010-10-22 2011-05-04 浙江省电力公司 Method for matching power flow of superior and inferior power grids in dispatcher training simulators (DTSs)
CN104362664A (en) * 2014-07-28 2015-02-18 浙江工业大学 Grid connection method of medium-voltage microgrid system
CN104167732A (en) * 2014-08-25 2014-11-26 中国南方电网有限责任公司 Grid equivalent method based on phase angle difference of call lines
CN104319772A (en) * 2014-10-21 2015-01-28 深圳大学 Method for symmetrically acquiring loss power sensitivity from AC-DC hybrid power system
CN104657618A (en) * 2015-03-06 2015-05-27 南车株洲电力机车研究所有限公司 Splitting fault modeling method for traction substations and splitting voltage determining method
CN104657618B (en) * 2015-03-06 2018-01-05 南车株洲电力机车研究所有限公司 A kind of traction substation off-the-line fault modeling method and solution column voltage determine method

Also Published As

Publication number Publication date
CN100570983C (en) 2009-12-16

Similar Documents

Publication Publication Date Title
Bevrani et al. Intelligent automatic generation control
Bidram et al. Distributed control systems for small-scale power networks: Using multiagent cooperative control theory
Olivares et al. Trends in microgrid control
Dufour et al. On the use of real-time 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 delay-dependent stability of multi-area load frequency control
CN104330980B (en) A kind of micro-capacitance 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 multi-terminal VSC-HVDC transmission systems
EP2182207B1 (en) Wide-area transmission control of windfarms
Ray et al. Wide-area signal-based optimalneurocontroller for a UPFC
CN102081867B (en) All-digital high-voltage direct current transmission training stimulation system
CN100461579C (en) Method for controlling coordination voltage of regional power grid and provincial power grid
Makarov et al. PMU-based wide-area security assessment: concept, method, and implementation
Chakrabortty Wide-area damping control of power systems using dynamic clustering and TCSC-based 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 cyber-physical energy systems
CN103715721B (en) The real-time peak regulating method of alternating current-direct current mixing power grid joint
Biegel et al. Distributed low-complexity controller for wind power plant in derated operation
Shoults et al. Multi-area 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 non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20091216

Termination date: 20191019