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

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:

$\left\{\begin{array}{c}{P}_i^{\mathrm{sp}}=V_i\underset{i\∈j}{\mathrm{\Σ}}{V}_j(G_{\mathrm{ij}}\cos {\mathrm{\θ}}_{\mathrm{ij}}+B_{\mathrm{ij}}\sin {\mathrm{\θ}}_{\mathrm{ij}})\\ Q_i^{\mathrm{sp}}=V_i\underset{i\∈j}{\mathrm{\Σ}}{V}_j(G_{\mathrm{ij}}{\sin \mathrm{\θ}}_{\mathrm{ij}}-B_{\mathrm{ij}}{\cos \mathrm{\θ}}_{\mathrm{ij}})\end{array}\right.$

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:

$\left\{\begin{array}{c}{P}_i^{\mathrm{sp}}=V_i\underset{i\∈j}{\mathrm{\Σ}}{V}_j(G_{\mathrm{ij}}{\cos \mathrm{\θ}}_{\mathrm{ij}}+B_{\mathrm{ij}}\sin {\mathrm{\θ}}_{\mathrm{ij}})\\ V_i^{\mathrm{sp}}=V_i\end{array}\right.$

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:

$\left\{\begin{array}{c}{V}_i^{\mathrm{sp}}=V_i\\ {\mathrm{\θ}}_i^{\mathrm{sp}}={\mathrm{\θ}}_i\end{array}\right.$

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_{\mathrm{SM}}={\tilde{B}}_{\mathrm{SM}}^{\′}{\tilde{B}}_{\mathrm{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:

$\mathrm{\Δ}{P}_M={\tilde{D}}_{\mathrm{MS}}\mathrm{\Δ}{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_{\mathrm{SR}}={\tilde{B}}_{\mathrm{SR}}^{\′\′}{\tilde{B}}_{\mathrm{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:

$\mathrm{\Δ}{Q}_R={\tilde{S}}_{\mathrm{RS}}\mathrm{\Δ}{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:

$\left\{\begin{array}{c}{\mathrm{\θ}}_i^{\mathrm{sp}}={\mathrm{\θ}}_i\\ Q_i^{\mathrm{sp}}=V_i\underset{i\∈j}{\mathrm{\Σ}}{V}_j(G_{\mathrm{ij}}\sin {\mathrm{\θ}}_{\mathrm{ij}}-B_{\mathrm{ij}}{\cos \mathrm{\θ}}_{\mathrm{ij}})\end{array}\right.$

θ 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_{\mathrm{SM}}={\tilde{B}}_{\mathrm{SM}}^{\′}{\tilde{B}}_{\mathrm{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:

$\mathrm{\Δ}{P}_M={\tilde{D}}_{\mathrm{MS}}\mathrm{\Δ}{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_{\mathrm{SR}}=-{\tilde{B}}_{\mathrm{SS}}^{\′\′-1}{\tilde{B}}_{\mathrm{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:

$\mathrm{\Δ}{V}_R={\tilde{R}}_{\mathrm{RS}}\mathrm{\Δ}{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:

$\mathrm{\Δ}{P}_M={\tilde{D}}_{\mathrm{MS}}\mathrm{\Δ}{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:

$\mathrm{\Δ}{Q}_R={\tilde{S}}_{\mathrm{RS}}\mathrm{\Δ}{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:

$\mathrm{\Δ}{P}_M={\tilde{D}}_{\mathrm{MS}}\mathrm{\Δ}{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:

$\mathrm{\Δ}{V}_R={\tilde{R}}_{\mathrm{RS}}\mathrm{\Δ}{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.

Images