CN101221599A - Interconnection system distributed tidal current computing method based on computation model split joint - Google Patents

Abstract

The invention relates to an interacted system distribution type tide computing method based on computational model marriage, which comprises the following steps: defining the handshaking branch and boundary node aggregate between subnets; respectively forming subnet information amount; selecting a main network and a slave network; merging the information amount of the main network and the slave network to obtain overall information amount and consequently obtain an overall tide computation equation set; generating a node number check list at the main network side or generating respective node number check lists of the subnets; unitedly selecting a voltage iteration initial value; solving at the main network side to obtain an overall integral tide result; returning the tide data of the slave network from the main network to the slave network in a complete transmission mode or a partial transmission mode to complete the distribution type tide computation. The invention avoids a large amount of model marriage and subsequent maintenance brought about by full network integral modeling or actual model merging, avoids complex algorithm organization patterns and simultaneously obtains a computing result in full accordance with that of centralized modeling.

Description

Interconnection system distributed tidal current computing method based on the computation model splicing

Technical field

The invention belongs to the Power System Analysis field, the present invention relates to a kind of electric system distributed tidal current computing method more precisely.

Background technology

Along with the needs of development of science and technology and dispatching concentrationization, electrical network is interconnected day by day to be one, but has contradiction between the scheduling of itself and electric system layering and zoning and the control system.The mode of trend is calculated in the many employings of EMS of each power-management centre at present alone, and outer net can not accurately be reflected to the calculating of Intranet to the influence of Intranet, though adopted the external network equivalent method, result of calculation is still not accurate enough.

Therefore, need research how in trend is calculated, the influence of outer net to Intranet accurately to be reflected.Wherein a kind of thinking is to coordinate between many dispatching control centers and organize, use the geographical data that go up each dispatching center of wide area distribution to calculate the integrated trend of the whole network, thereby accurately reflect the influence of the electric variation of outer net to Intranet, this is the research contents that distributed tidal current is calculated.

Concrete coordination and the tissue that how to carry out distributed trend between the dispatching center has several different methods, and the model that is adopted during according to analytical calculation can be summarized as and adopts centralized modeling and adopt two thinkings of distributed modeling.Each power-management centre EMS has all carried out modeling to the electrical network of administering or being concerned about separately: adopting the coordination of centralized modeling and organization work mainly is to be stitched together at the disperal pattern with each EMS to set up the unified integrated model of a whole network; And adopt the coordination of distributed modeling and organization work mainly be between each EMS structure certain with the corresponding external iteration form of tuning algorithm, and according to its between each EMS, carry out computational data alternately and iteration.

Document one " method of jointing multiple regional net models of electric power system " (Chinese patent application numbers 200610166302.1, open day on July 25th, 2007) has disclosed a kind of method of jointing multiple regional net models of electric power system.Use this method, can be spliced into centralized full model by the network model that each zone is actual.Do not use full model to carry out the content of Power System Analysis though document one does not relate to, the purpose that obtains full model is to use its monitoring of carrying out electric system and integrated analysis.If use centralized full model to calculate the trend of each regional subnet, promptly be referred to as to adopt the distributed tidal current of centralized modeling to calculate.The subject matter that this method exists is to have huge workload when carrying out the real network model splicing.In addition, because the electric network model change is frequent, in maintenance, also there is same problem.

Adopt the present finding of distributed modeling thinking that following method is arranged:

Document two " sending out transmission ﹠ distribution overall situation Power System Analysis " (Automation of Electric Systems was rolled up the 17th page of the 1st phase in 2000 the 24th) is being sent out under the transmission and distribution network integrated analysis background, has disclosed " principal and subordinate's splitting-up method " that a kind of distributed tidal current is calculated.In this algorithm, each subnet status inequality has the branch of principal and subordinate (Master-Slave).Major network generally is counted as from the broad sense power end of net in stacking outside the riding instead of walk, and from netting the broad sense load side that generally is counted as major network.It is the broad sense power end that this algorithm is looked the power transmission network side, and looking the power distribution network side is the broad sense load side, has carried out the distributed tidal current calculation.Document three " based on the DTS of internet associating Training Simulation System " (relay 2006 the 34th the 2nd phase of volume the 70th page) disclosed a kind of will be as above the thought of " principal and subordinate's splitting-up method " apply to distributed power flow algorithm in the interconnected analysis of power transmission network.This algorithm is transferred side-looking with ground in calculating be that province transfers the broad sense load of side to bring in the distributed tidal current calculation of carrying out under the power transmission network background.Use this algorithm to find the solution the subject matter following points of the interconnection system distributed trend of power transmission network: the situation that 1. is difficult to transfer side 110KV and the cyclization operation of following system with handling; During the unified starting condition that forms 2. that the interacted system trend calculates, ground transfers the imbalance power of the sub-DTS of side system self to be difficult to processing; 3. this algorithm needs each sub-DTS system core power flow algorithm is changed, thereby has reduced its adaptability to the different manufacturers platform.

Document four " decomposition-coordination approach that the united power grid trend is calculated " (Wuhan Univ. of Water Conservancy and Electric Power's journal was rolled up the 35th page of the 1st phase in 1997 the 19th) has disclosed a kind of distributed trend composition decomposition algorithm of constructing the external iteration form based on idle and relation voltage.Document five " based on the multizone interacted system dynamic power flow composition decomposition calculating of asynchronous iteration " (Automation of Electric Systems was rolled up the 1st page of the 24th phase in 2003 the 27th) and additional document six " distributed tidal current is calculated replenishing and improving of asynchronous iteration pattern " (Automation of Electric Systems was rolled up the 12nd page of the 2nd phase in 2007 the 31st) thereof have disclosed the distributed trend composition decomposition algorithm that a kind of Dai Weinan of use equivalence theory is constructed the external iteration form.Be different from the master-slave mode algorithm, these two kinds of algorithms are not left for from the special nature of subnet and structure and are looked at problem, each subnet is equal in the calculating, its basic ideas are after each subnet independently carries out trend calculating, the fixed point iteration form of tectonic boundary dotted state variable in external iteration, iteration to the frontier point state variable difference between subnet reaches convergence precision, tries hard to make result of calculation identical with integrated result of calculation.Its common subject matter that exists is: the algorithm organization pattern is complicated, and there is certain error in its result of calculation.

Summary of the invention

For solving the deficiencies in the prior art; the object of the present invention is to provide a kind of huge workload can avoid centralized modeling or carry out the real network model splicing time; on all four result of calculation in the time of obtaining and adopt centralized Modeling Calculation again; in calculating, distributed tidal current avoids mutual external network equivalent simultaneously; between subnet, do not coordinate and iteration; simplify its enterprise schema, thus on engineering interconnection system distributed tidal current computing method easy to implement based on the computation model splicing.

To achieve these goals, the present invention takes following technical scheme to realize:

A kind of interconnection system distributed tidal current computing method based on the computation model splicing comprises the following steps:

(1), define the contact branch road between each regional power grid, described regional power grid abbreviates subnet as, sets up " contact branch road definition list ", and defines border splicing node set, interior nodes set, the exterior node set of each subnet;

(2), each subnet forms bus admittance matrix, the node power injection vector sum node type information vector that comprises interior nodes set and border splicing node set respectively;

(3), select a certain subnet for calculating major network, abbreviate major network as, all the other subnets abbreviate as from net for calculating from net;

(4), major network injects vector sum node type information vector with self bus admittance matrix, node power and respectively is stitched together from the corresponding matrix and the vector of net, obtain the global node admittance matrix respectively, global node power injects vector sum global node type information vector, obtains overall trend group of equations by these three kinds of global information amounts; Described bus admittance matrix, node power inject vector sum node type information vector and are called computation model;

(5), each subnet generates the node serial number table of comparisons separately in the computation model splicing, perhaps become the node serial number table of comparisons of all-ones subnet at the major network adnation;

(6), the unified voltage iteration initial value of selecting to resolve overall trend group of equations, perhaps follow the splicing of computation model, obtain this value by the method for splicing each subnet voltage iteration initial value vector;

(7), finish finding the solution of overall trend group of equations, obtain overall integrated trend result in the major network side;

(8), the major network side will be passed back from the flow data of net by full transfer mode or part transfer mode; If the node serial number table of comparisons of all-ones subnet all is created on the major network side, comprise in this step that major network side general respectively pass back corresponding from net from the node serial number table of comparisons of net.

Aforesaid interconnection system distributed tidal current computing method based on the computation model splicing is characterized in that the whole contact branch roads between two subnets be able to carry out directly electric cutting apart to two subnets.

Aforesaid interconnection system distributed tidal current computing method based on the computation model splicing is characterized in that described contact branch road is chosen as circuit or the transformer at the network model middle impedance non-zero of EMS.

Aforesaid interconnection system distributed tidal current computing method based on the computation model splicing is characterized in that described contact branch road is chosen as the equivalent branch road of its a certain side when the edge device between subnet is three volume changes.

Aforesaid interconnection system distributed tidal current computing method based on the computation model splicing is characterized in that bus admittance matrix, node power that described each subnet forms respectively inject the node that vector sum node type information vector does not comprise the exterior node set.

Aforesaid interconnection system distributed tidal current computing method based on the computation model splicing, it is characterized in that the splicing of described overall admittance matrix, carry out in two steps generally: by the exponent number expansion each subnet bus admittance matrix is stitched together earlier, in this matrix, add the contact branch road again, the influence that is produced of this matrix is revised this matrix according to adding the contact branch road.

Aforesaid interconnection system distributed tidal current computing method based on the computation model splicing, it is characterized in that each subnet node serial number table of comparisons, wherein write down this subnet and removed corresponding relation between the difference numbering of all nodes in splicing outside the exterior node set.

Aforesaid interconnection system distributed tidal current computing method based on the computation model splicing, it is characterized in that under described full transfer mode, major network is to return whole required flow datas from network transmission, comprise node serial number, node voltage amplitude, phase angle, node injecting power, branch road trend, subnet self does not carry out any trend and calculates.

Aforesaid interconnection system distributed tidal current computing method based on the computation model splicing, it is characterized in that under described part transfer mode, major network only to return node serial number, node voltage amplitude, phase angle from network transmission, calculate and finished by subnet oneself by branch road trend and other electric parameters.

Beneficial effect of the present invention and advantage: adopt method of the present invention to carry out distributed tidal current and calculate, solved the problem that existing method exists.At first, adopt algorithm of the present invention to need not the integrated modeling of the whole network or carry out actual model splicing, thereby avoided a large amount of model splicing that brings thus and maintenance work after this.And, adopt the distributed trend of algorithm computation of the present invention between subnet, also not have any outer coordination iterative computation, thereby avoided complicated algorithm enterprise schema.On all four result of calculation when at last, adopting algorithm of the present invention can obtain and adopt centralized modeling.In a word, algorithm of the present invention will be concentrated the advantages of Modeling Calculation and the present various distributed coordination iterative algorithms that propose, it is the centralized distributed power flow algorithm of a kind of standard, neither carry out actual model splicing, also between subnet, do not coordinate and iteration, and the integrated calculating of its result calculated and the whole network is in full accord.

Description of drawings

Fig. 1 is with major network S ₁, from the net S ₂The interconnection system distributed tidal current based on the computation model splicing that the present invention that two subnets are example adopts is calculated FB(flow block).

Fig. 2 is two sub-pessimistic concurrency control splicing synoptic diagram that are used to illustrate computation model splicing scheme.

Fig. 3 is to be the computation model splicing synoptic diagram of prototype system with the IEEE9 node system.

Fig. 4 is the equivalent circuit diagram of the transformer when ignoring the static exciter loop.

Embodiment

Below in conjunction with accompanying drawing the present invention is done concrete being described below:

In the present invention, disclosed method and the enterprise schema that bus admittance matrix, node injecting power vector sum node type information vector (claiming that these three kinds of quantity of information are computation model) splicing is obtained overall trend group of equations.Power flow equation group under the polar coordinates is as the formula (1):

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

In the formula: i=1 ..., N, j ∈ i represent the label behind the ∑ number be the node of j must be directly and node i link, and comprise the situation of j=i, P _iThe meritorious injecting power of expression node i, Q _iThe idle injecting power of expression node i, U _iThe voltage magnitude of expression node i, U _jThe voltage magnitude of expression node j, G _IjAnd B _IjThe transconductance of expression node i and node j correspondence in bus admittance matrix and mutual susceptance, θ _Ij=θ _i-θ _j, be the phase angle difference of node i and node j two node voltages.When finding the solution the power flow equation group with common Newton method or Fast Decoupled method, required the whole network information matrix can be divided into following 3 kinds: bus admittance matrix, the initial injecting power column matrix of node/vector, node type information column matrix/vector.Known these 3 matrixes add selected iterative initial value, can find the solution conventional trend.Under the background of distributed modeling, if these 3 matrixes of each subnet can be stitched together, just obtain carrying out the integrated overall trend group of equations of resolving, the result that its result of calculation will be when adopting centralized modeling is in full accord.Iterative initial value can obtain by splicing, also can unified again selection after the splicing of overall trend machine calculation system of equations is finished.In distributed tidal current is calculated, be typically chosen in and have the sub-network side of enriching computational resource most and finish splicing and calculate, be called the calculating major network, abbreviate major network as; All the other subnets are called calculating from net, abbreviate as from net.The essence of computation model splicing is that each subnet power flow equation group is spliced into overall trend group of equations, and the computation model joining method comprises that bus admittance matrix joining method, node serial number contrast table generating method, node power inject the content of vectorial joining method and node type information vector joining method four aspects.

After integrated trend calculating is finished, because splicing and calculating are finished in the major network side, major network side trend is the result obtain naturally, but need go back by the major network transmission from net trend result, design has two kinds of patterns: 1. full transfer mode: major network is to return whole required flow datas from network transmission, mainly comprise: node serial number, node voltage amplitude, phase angle, node injecting power, branch road trend, subnet self are not carried out any trend and are calculated; 2. part transfer mode: major network only to return node serial number, node voltage amplitude, phase angle from network transmission, calculate and finished by subnet oneself by branch road trend and other electric parameters.Full transfer mode advantage is that subnet need not to calculate, and shortcoming is that volume of transmitted data is big; Part transfer mode relative merits in contrast.Send back from the node serial number of net and generally be reduced to the preceding numbering of splicing, can directly use after obtaining from net, the numbering transformation also can be by finishing from netting self.In addition, from net and major network the part of repetition modeling is arranged, the flow data of this part also will send back from net, and its data corresponding relation solves with the way of the unified name of master and slave net device id.

With subnet S ₁For major network, with subnet S ₂From net is that the algorithm organization pattern of example is seen Fig. 1.

In the present invention, disclosed each subnet bus admittance matrix is stitched together and obtained the method for global node admittance matrix.Though do not carry out the network model splicing of any reality among the present invention, for the splicing scheme of the abstract bus admittance matrix of explanation directly perceived, the overall plan of at first computation model being spliced is described as follows on this aspect of real network model splicing here.

Fig. 2 is two sub-pessimistic concurrency control splicing synoptic diagram.The electric system application software ubiquity of domestic each dispatching center the situation that modeling partly repeats, the generalized case when Fig. 2 has represented two subnets existence repetition modelings, and S represents the whole network among the figure, S ₁, S ₂Be two subnets of the whole network S, the left and right figure of Fig. 2 (a) represents subnet S respectively ₁, S ₂Complete model, S _oRepresent that both repeat the part of modeling, n _i, n _j, n _e, n _jBe node, l _IeAnd l _JfBe S _oIn two branch roads and be able to be with S ₁And S ₂Carry out electric cutting apart, i.e. l _IeAnd l _IfBe S ₁And S ₂The contact branch road.This paper is example with Fig. 2, and it is as follows to propose the model splicing scheme:

(1) get in touch with branch road between the definition subnet, it should carry out electric cutting apart to subnet.Among Fig. 2 (a), suppose l _IeAnd l _IfBe defined as getting in touch with branch road.

(2) define each sub-net boundary splicing node set, interior nodes set, exterior node set.Definition contact branch road both sides end points is the border splicing node of side subnet separately, its set N _BExpression; Splicing intranodal in definition border is the interior nodes of this subnet near the node of certain sub-network side, its set N _IExpression, and be this subnet exterior node away from the node of this sub-network side, its set N _EExpression.Among the left figure of Fig. 2 (a), node n _i, n _jBe S ₁Border splicing node, its set is N _B1Node n _i, n _jThe S in left side ₁Node set is N _I1, node n _i, n _jThe S on right side ₁Node set is N _E1If S ₁Node set represent then obviously have with N $N_{B_{}}$ ${{}_1}_{}\&cup;{\mathrm{<figure-callout\; id="1"\; label="N\; I"\; filenames="S2007101915421E00031.png"\; state="\{\{state\}\}">N</figure-callout>}}_{I_{}}{\&cup;\mathrm{<figure-callout\; id="1"\; label="N\; E"\; filenames="S2007101915421E00031.png"\; state="\{\{state\}\}">N</figure-callout>}}_{E_{}}=N.$ S ₂Roughly the same, can get N _B2, N _I2, N _E2

(3) each subnet is all rejected the outside segments point set, adds the contact branch road between each border splicing node, and splicing obtains the integrated model of the whole network.Among Fig. 2 (b), major network S ₁Reject exterior node set N _E1, from net S ₂Reject exterior node set N _E2After, at node n _iAnd n _eBetween add contact branch road l _Ie, node n _jAnd n _fBetween add contact branch road l _JfCarry out the whole network model splicing, obtain the whole network complete model shown in Fig. 2 (c).

For the situation that repeats modeling between subnet more complicated situation and a plurality of subnet splice the situation of calculating, as above model splicing scheme step is equally suitable.

The selection of contact branch road need be satisfied following several condition.

Whole contact branch roads between (1) two subnet should carry out directly electric cutting apart to two subnets;

The circuit or the transformer of impedance non-zero when (2) the contact branch road is chosen as network modelling;

When (3) being three volume changes as if communication gear, the contact branch road is selected the equivalent branch road of its a certain side.

Can set about carrying out the splicing of global node admittance matrix after the model splicing scheme is determined, method is as follows:

Each subnet bus admittance matrix is spliced into the core that the global node admittance matrix is overall trend group of equations splicing.Suppose major network S ₁Middle border splicing node set N _B.1With interior nodes set N _I.1The interstitial content sum is m, from net S ₂Corresponding number is n, according to the model splicing scheme that is as above proposed, is that example proposes following admittance matrix splicing scheme with Fig. 3.Fig. 3 is to be the computation model splicing synoptic diagram of prototype with the IEEE9 node system.The m value is 5 among Fig. 3, and the n value is 4.The splicing of global node admittance matrix divides following step to carry out.

(1) each subnet topology is to the border splicing node formation admittance battle array separately of self.Among Fig. 3 (b), major network S ₁With from the net S ₂Each forms admittance battle array Y _S.1And Y _S.2Y _S.1Exponent number is m, Y _S.2Exponent number is n.

(2) select a subnet as major network, all the other are from net, and major network admittance battle array is placed the upper left corner, places the lower right corner to carry out exponent number expansion splicing successively from net admittance battle array.Among Fig. 3, select S ₁Be major network, S ₂For from net, splice the admittance battle array Y that obtains after exponent number expands _Mo' be:

$Y_{\mathrm{mo}}^{\&prime;}=\left[\begin{array}{cc}{Y}_{S.1}& \mathrm{zeros}(m,n)\\ \mathrm{zeros}(n,m)& Y_{S.2}\end{array}\right]---\left(2\right)$

In the formula: ((n m) analogizes zeros zeros for m, the n) null matrix on expression m * n rank.Yi Zhi, admittance battle array Y _Mo' be the square formation on m+n rank.Need explanation, major network admittance battle array and from net admittance battle array in the position in when splicing not necessarily respectively in the upper left corner and the lower right corner, can place arbitrarily, limited in the literary composition and explained in order to know.

(3) at admittance battle array Y _Mo' middle increase and contact branch road corresponding transadmittance element, and it is as follows to revise the self-admittance of former each sub-net boundary splicing node correspondence:

A) at admittance battle array Y _Mo' middle increase and the corresponding transadmittance element of contact branch road.Among Fig. 3 (c), with interconnection l ₁Be example, suppose that its two end node is at admittance battle array Y _Mo' middle corresponding row number (that is numbering) is e, f, and its admittance is y _Ef, admittance is y over the ground _Ef0, then at Y _Mo2 off-diagonal elements of ' middle interpolation are suc as formula (3).Interconnection l ₂Roughly the same.

Y _mo.ef′＝Y _mo.fe′＝-y _ef (3)

B) at admittance battle array Y _MoThe self-admittance of each border splicing node correspondence of ' middle correction.With splicing posterior nodal point n _e, n _fBe example, n _eCorresponding self-admittance Y _Mo.ee' and n _fCorresponding self-admittance Y _Mo.ff' following increment arranged:

$\mathrm{\&Delta;}{Y}_{\mathrm{mo}.\mathrm{ee}}^{\&prime;}=\mathrm{\&Delta;}{Y}_{\mathrm{mo}.\mathrm{ff}}^{\&prime;}=y_{\mathrm{ef}}+y_{{\mathrm{ef}}_0}/2---\left(4\right)$

So far, Y _Mo' splicing finishes.

C) to admittance battle array Y _Mo' each row institute corresponding node is numbered optimization, is about to admittance battle array Y _Mo' carry out corresponding ranks to exchange elementary transformation and obtain overall admittance matrix Y _Mo

When the contact branch road becomes for contact but not during interconnection, the splicing thinking there is no difference, just formula (3), (4) or change to some extent.When ignoring the static exciter loop, the transformer branch road can be simulated with ideal transformer of its leakage reactance series connection, shown in Fig. 4 (a).Its ∏ type equivalent circuit is shown in Fig. 4 (b).Wherein, z _TBe the transformer leakage reactance, y _T=1/z _T, K is a transformer voltage ratio, node e and node f are transformer equivalent circuit both sides end points,

With Be respectively the terminal voltage vector of node e and node f, With

It is respectively the injection current vector of node e and node f.

According to Fig. 4, when easily knowing by contact change splicing admittance battle array, formula (3) variation is:

$Y_{\mathrm{mo}.\mathrm{ef}}^{\&prime;}=Y_{\mathrm{mo}.\mathrm{fe}}^{\&prime;}=-\frac{y_T}{K}---\left(5\right)$

And correction is changed in the formula (4):

$\left\{\begin{array}{c}\mathrm{\&Delta;}{Y}_{\mathrm{mo}.\mathrm{ee}}^{\&prime;}=\frac{K-1}{K}{y}_T+\frac{y_T}{K}=y_T\\ {\mathrm{\&Delta;Y}}_{\mathrm{mo}.\mathrm{ff}}^{\&prime;}=\frac{y_T}{K}+\frac{1-K}{K^2}{y}_T=\frac{y_T}{K^2}\end{array}\right.---\left(6\right)$

In the present invention, disclosed a kind of method of setting up the node serial number table of comparisons.Be accompanied by admittance battle array Y _Mo' be optimized numbering to obtain overall admittance matrix Y _MoProcess, tackle each subnet and all set up the node serial number table of comparisons, write down each subnet in the table and removed the numbering corresponding relation of all nodes in splicing outside the exterior node.With major network S among Fig. 1 ₁Be example, table 1 (a) is its node serial number table of comparisons.Major network S before the splicing ₁Node serial number be incremented to m from 1, form and admittance battle array Y _S1The column vector α of row correspondence ₁, be placed on the 1st row; Be spliced into admittance battle array Y _Mo' after, obtain and its row number corresponding node serial number column vector α _Mo1', be placed on the 2nd row; Be numbered to optimize and obtain overall admittance matrix Y _MoAfter, obtain and its row number corresponding node serial number column vector α _Mo1, be placed on the 3rd row, build table and finish.From net S ₂Build table method roughly the same, but because admittance battle array Y _S2At admittance battle array Y _Mo' the lower right, node serial number column vector α _Mo2' and α ₂And inequality, its each component has all added m, and (m is major network S ₁The node number, be called the splicing side-play amount, Hereinafter the same), shown in table 1 (b).

Table 1S ₁With S ₂The node serial number table of comparisons (accompanying drawing 1)

In the present invention, also disclose a kind of each subnet node injecting power vector is stitched together and obtained the method for global node injecting power vector.Suppose major network S ₁Power flow equation group node injecting power column vector be With admittance battle array Y _S.1Corresponding; From net S ₂The respective column vector is

Then splicing can be carried out in two steps:

(1) with node injecting power column vector

Carry out dimension and expand splicing, obtain and admittance battle array Y _Mo' corresponding global node injecting power vector without numbering optimization

As follows:

${\stackrel{\&CenterDot;}{S}}_{\mathrm{mo}}^{\&prime;}=\left[\begin{array}{c}{\stackrel{\&CenterDot;}{S}}_1\\ {\stackrel{\&CenterDot;}{S}}_2\end{array}\right]---\left(7\right)$

(2) look into the node serial number table of comparisons, according to node injecting power column vector The numbering corresponding relation of origin node before and after splicing is right Carry out accordingly the row transposing, make the numbering after capable number of each component and the corresponding node optimization of former institute consistent, just obtain and overall admittance matrix Y _MoCorresponding global node injecting power vector

Perhaps also can be to admittance battle array Y _Mo'.When being optimized numbering, right

Carry out corresponding transposing and obtain global node injecting power vector

In the present invention, also disclose a kind of each subnet node type information vector is stitched together and obtained the method for global node type information vector.The node of electric system is divided into three types of PQ node, PV node and balance node in trend is calculated.If major network S before the splicing ₁Middle node type information vector is: only balance node numbering is placed on monobasic balance node column vector T _Ref.1In, the PQ node number is placed on its PQ node number column vector T _Pq.1In, the PV node number is placed on PV node number column vector T _Pv.1In, from net S ₂In corresponding vector be: T _Ref.2, T _Pq.2, T _Pv.2, the splicing of node type information vector is carried out as follows.

(1) with column vector T _Ref.2, T _Pq.2, T _Pv.2Each component all adds splicing side-play amount m.Because the whole network only need keep a balance node, from net S ₂Balance node after splicing is finished, will be converted into the overall trend PV node in calculating, so get ${T_{\mathrm{pv}.2}}^{\&prime;}={\left[\begin{array}{cc}{T}_{\mathrm{pv}.2}^T& T_{\mathrm{ref}.2}\end{array}\right]}^T.$

(2) get ${T_{\mathrm{ref}}}^{\&prime;}={\mathrm{<figure-callout\; id="1"\; label="T\; ref"\; filenames="S2007101915421E00031.png"\; state="\{\{state\}\}">T</figure-callout>}}_{{\mathrm{ref}}_{}},$ ${T_{\mathrm{pq}}}^{\&prime;}=\left[\begin{array}{c}{T}_{{\mathrm{pq}}_1}\\ {\mathrm{<figure-callout\; id="2"\; label="T\; pq"\; filenames="S2007101915421E00031.png"\; state="\{\{state\}\}">T</figure-callout>}}_{{\mathrm{pq}}_{}}\end{array}\right],$ ${T_{\mathrm{pv}}}^{\&prime;}=\left[\begin{array}{c}{T}_{{\mathrm{pv}}_1}\\ {\mathrm{<figure-callout\; id="2"\; label="T\; pv"\; filenames="S2007101915421E00031.png"\; state="\{\{state\}\}">T</figure-callout>}}_{{\mathrm{pv}}_{}}^{{}_2}\end{array}\right],$ All with admittance battle array Y _Mo' corresponding.

(3) by major network S ₁, from the net S ₂The node serial number table of comparisons checks in column vector T _Ref', T _Pq', T _PvNew numbering after ' each component optimization numbering is formed and overall admittance matrix Y separately _MoCorresponding overall balance node column vector T _Ref, overall PQ node number column vector T _PqWith overall PV node number column vector T _PvThese three global node type information vectors, splicing is finished.

Perhaps, the also available mode of afterwards splicing of tabling look-up earlier of global node type information vector obtains.

Embodiment

Be one embodiment of the present of invention below, comprised that adopting algorithm of the present invention is process and the result that the distributed tidal current that carries out of prototype system is calculated with the IEEE9 node system.Other feature, purpose and advantage of the present invention also can be found out from the explanation of embodiment and accompanying drawing.

Accompanying drawing 3 has illustrated that with the IEEE9 node system be process and the result's comparison that prototype system is used algorithm of the present invention, and its network and calculating parameter all adopt IEEE9 node system normal data.Fig. 3 (b) has represented that two have repetition modeling contact branch road l ₁, l ₂Two subnet S ₁, S ₂, have node serial number separately since 1.The embodiment step is as follows.

(1) definition l in " contact branch road definition list " ₁, l ₂Be the contact branch road between subnet; According to the definition and the technical scheme of the present invention of contact branch road, definition subnet S ₁, S ₂Border splicing node set N _B, interior nodes set N _I, exterior node set N _ES ₁Boundary node set N _B.1={ n ₁, n ₂, N _I.1={ n ₃, n ₄, n ₅, N _E.1={ n _i, n _j, S ₂Boundary node set N _B.1={ N ₁, n ₂, N _I.1={ n ₃, n ₄, N _E.1={ n _m, n _n.

(2) subnet S ₁And S ₂Topology is spliced node to the border, each self-forming bus admittance matrix, node injecting power vector sum node type information vector separately.Subnet S ₁And S ₂Each forms admittance battle array Y _S.1And Y _S.2Admittance battle array Y _S.1Exponent number is 5, admittance battle array Y _S.2Exponent number is 4.Subnet S ₁Power flow equation group node injecting power column vector be

With admittance battle array Y _S.1Correspondence, subnet S ₂The respective column vector is Subnet S ₁Middle node type information vector is: only balance node numbering is placed on monobasic column vector T _Ref.1In, the PQ node number is placed on column vector T _Pq.1In, the PV node number is placed on column vector T _Pv.1In.Subnet S ₂In corresponding vector be: T _Ref.2, T _Pq.2, T _Pv.2

(3) select subnet S ₁Be major network, subnet S ₂For from net.

(4) major network S ₁With self with receive from net S ₂The bus admittance matrix that sends, node power inject vector sum node type information vector (claiming that these three kinds of quantity of information are computation model) and are stitched together according to technical scheme of the present invention, finally obtain overall admittance matrix Y _Mo, global node injecting power vector Global node type information vector T _Ref, T _Pq, T _Pv, finish the splicing work of computation model.Its process on aspect directly perceived shown in Fig. 3 (c), (d).Wherein, node serial number is obtained by half dynamic optimization method among Fig. 3 (d).

(5) be accompanied by the computation model splicing, subnet S ₁And S ₂Each self-generating node serial number table of comparisons sees Table 2 (a) and (b) respectively.

(6) the overall trend group of equations voltage iteration initial value of unified selection.

(7) obtained overall trend group of equations by above splicing, at major network S ₁Side is finished finding the solution of this overall situation power flow equation group, obtains overall integrated calculation of tidal current.Crucial result of calculation amount the whole network voltage phasor is designated as

See Table 3.

(8) major network S ₁Side will be from net S by full transfer mode ₂Flow data transfer back to from the net S ₂

Result verification: for the correctness of resulting result of calculation is described, we carry out conventional trend calculating again under IEEE9 node system normal data shown in Fig. 3 (a), obtain the whole network voltage phasor and be designated as Regard this result of calculation as this embodiment distributed tidal current calculate result's true value, see Table 3.Notice the node serial number corresponding relation among Fig. 3 (a), (d), can see,

With

In full accord, the node voltage and the true value that promptly adopt algorithm of the present invention to obtain are in full accord.

Table 2S ₁With S ₂The node serial number table of comparisons

Table 3 stitching algorithm trend result and true value (perunit value) are relatively

The present invention announces as above with preferred embodiment, but the foregoing description does not limit the present invention in any form, and all employings are equal to the technical scheme that form obtained of replacement or equivalent transformation, all drop within protection scope of the present invention.

Claims (9)

1. based on the interconnection system distributed tidal current computing method of computation model splicing, comprise the following steps:

(1), define the contact branch road between each regional power grid, described regional power grid abbreviates subnet as, sets up " contact branch road definition list ", and defines border splicing node set, interior nodes set, the exterior node set of each subnet;

(2), each subnet forms bus admittance matrix, the node power injection vector sum node type information vector that comprises interior nodes set and border splicing node set respectively;

(3), select a certain subnet for calculating major network, abbreviate major network as, all the other subnets abbreviate as from net for calculating from net;

(4), major network injects vector sum node type information vector with self bus admittance matrix, node power and respectively is stitched together from the corresponding matrix and the vector of net, obtain the global node admittance matrix respectively, global node power injects vector sum global node type information vector, obtains overall trend group of equations by these three kinds of global information amounts; Described bus admittance matrix, node power inject vector sum node type information vector and are called computation model;

(5), each subnet generates the node serial number table of comparisons separately in the computation model splicing, perhaps become the node serial number table of comparisons of all-ones subnet at the major network adnation;

(6), the unified voltage iteration initial value of selecting to resolve overall trend group of equations, perhaps follow the splicing of computation model, obtain this value by the method for splicing each subnet voltage iteration initial value vector;

(7), finish finding the solution of overall trend group of equations, obtain overall integrated trend result in the major network side;

(8), the major network side will be passed back from the flow data of net by full transfer mode or part transfer mode; If the node serial number table of comparisons of all-ones subnet all is created on the major network side, comprise in this step that major network side general respectively pass back corresponding from net from the node serial number table of comparisons of net.

2. the interconnection system distributed tidal current computing method based on the computation model splicing according to claim 1 is characterized in that the whole contact branch roads between two subnets be able to carry out directly electric cutting apart to two subnets.

3. the interconnection system distributed tidal current computing method based on the computation model splicing according to claim 1 and 2 is characterized in that described contact branch road is chosen as circuit or the transformer at the network model middle impedance non-zero of EMS.

4. the interconnection system distributed tidal current computing method based on the computation model splicing according to claim 3 is characterized in that described contact branch road is chosen as the equivalent branch road of its a certain side when the edge device between subnet is three volume changes.

5. the interconnection system distributed tidal current computing method based on the computation model splicing according to claim 1 is characterized in that bus admittance matrix, node power that described each subnet forms respectively inject the node that vector sum node type information vector does not comprise the exterior node set.

6. the interconnection system distributed tidal current computing method based on the computation model splicing according to claim 1, it is characterized in that the splicing of described overall admittance matrix, carry out in two steps generally: by the exponent number expansion each subnet bus admittance matrix is stitched together earlier, in this matrix, add the contact branch road again, the influence that is produced of this matrix is revised this matrix according to adding the contact branch road.

7. the interconnection system distributed tidal current computing method based on the computation model splicing according to claim 1, it is characterized in that each subnet node serial number table of comparisons, wherein write down this subnet and removed corresponding relation between the difference numbering of all nodes in splicing outside the exterior node set.

8. the interconnection system distributed tidal current computing method based on the computation model splicing according to claim 1, it is characterized in that under described full transfer mode, major network is to return whole required flow datas from network transmission, comprise node serial number, node voltage amplitude, phase angle, node injecting power, branch road trend, subnet self does not carry out any trend and calculates.

9. the interconnection system distributed tidal current computing method based on the computation model splicing according to claim 1, it is characterized in that under described part transfer mode, major network only to return node serial number, node voltage amplitude, phase angle from network transmission, calculate and finished by subnet oneself by branch road trend and other electric parameters.

Images