CN103714490A - Large power grid on-line data multi-thread rapid-integration method - Google Patents

Abstract

The invention provides a large power grid on-line data multi-thread rapid-integration method. The method comprises the steps of 1) reading and parsing various types of power grid data, storing the power grid data into an internal storage, and partitioning the internal storage according to the power grid data of power grid equipment of different types; 2) employing a parallel processing mode to associate the power grid equipment in each internal storage part after the partition with a corresponding topology node and establishing an association topology node list with respect to the topology nodes; and 3) zoning the topology nodes corresponding to the internal storage parts based on the partition of the internal storage, employing the parallel processing mode to conduct traversal analysis of the association topology node list of the topology nodes in each zone, performing topology island numbering of each topology node, eliminating dead topology islands, determining cared topology islands, conducting topology merging of all the zones of the topology nodes through tie lines among the zones of the topology nodes, and generating integrated power grid on-line data. The power grid data integration method provided by the invention is high in integration efficiency and rapid and highly efficient in processing procedure.

Description

The quick integration method of a kind of scale grid line data multi-wire journey

Technical field

The present invention relates to electric system simulation analytical calculation field, be specifically related to the quick integration method of a kind of scale grid line data multi-wire journey.

Background technology

Oneself becomes the main tool of planning and design of power system, management and running and analysis and research the Digital Simulation of electric system, and the mathematical model of each element of electric system and the total system mathematical model consisting of it are the bases of power system digital simulation.Model and parameter is the important determinative of simulation result accuracy, directly affects decision scheme based on this.

The safe and reliable operation of electrical network is significant to social and economic development.Along with the continuous expansion of electrical network scale, genset is constantly gone into operation and the continuous access of UHV transmission, renewable resource power supply, power network wiring is day by day complicated, the situation of operation of power networks becomes increasingly complex, and the effect of safety on line Stability Assessment and early warning system is more and more significant.Safety on line Stability Assessment and early warning system comprise the computational analysis functions such as static stability analysis, transient stability analysis, dynamic stability analysis and stability margin assessment, above-mentioned various computational analysiss are all computational datas based on complete, so it is the basis of safety on line Stability Assessment and early warning system that online data is integrated.Along with the continuous popularization of " intelligent grid supporting system technology ", the scale of model that online data is integrated is increasing, thus study a kind of accurately, efficiently, online data integration method is necessary fast.

Traditional online data is integrated with two thinkings: 1) take the online measuring of EMS system and model data and offline mode data is basis, integrates and generates flow data and stable data; 2) utilize online measuring and the model data of EMS system, generate tideway integrating data, then according to off-line dynamic parameter storehouse, generate corresponding stable data.Wherein thinking 1) in the tideway integrating of generation combine the feature of online mode and offline mode, thereby there is no to guarantee a consistance with online mode, larger with actual motion systematic error, seldom by this thinking, realize online data now and integrated; Thinking 2) be the main flow implementation that current online data is integrated.

Online data integrate mostly comprise that data read, Data Analysis, generation topology, topological analysis, computational data generate these steps.Traditional online data is integrated and is comprised following shortcoming: 1) integrate flow process and all adopt serial structure, when integral data amount is large, speed is obviously slack-off, is difficult to reach efficiently, online application requirements fast; 2) Data Integration is only simple data-switching, does not process the misdata existing in online data and the quality of online data is not carried out to qualitative assessment.

Summary of the invention

For overcoming above-mentioned defect, the invention provides the quick integration method of a kind of scale grid line data multi-wire journey, comprising:

Step 1, reads and resolves all kinds of electric network datas, and described electric network data is deposited to internal memory, and the internal memory corresponding to the electric network data of dissimilar grid equipment carries out piecemeal;

Step 2, the grid equipment in the internal memory after adopting parallel processing mode to piecemeal described in each is associated with topological node, and described topological node is set up to associated topologies node table;

Step 3, carries out subregion according to the piecemeal of described internal memory to topological node corresponding to each piece internal memory, adopts parallel processing mode to travel through analysis to the described associated topologies node table of the described topological node of subregion described in each;

Topological node described in each is carried out to topological island numbering, remove the dead island of topology and determine the topological island of being concerned about;

Interconnection by topological node by stages carries out topology merging by the subregion of topological node described in each, generates the powernet data of integrating.

In the first preferred embodiment provided by the invention: the described electric network data reading in described step 1 comprises: plant stand, bus, transformer, alternating current circuit, generator, load, shunt compensation equipment, series compensation equipment, transverter, DC line, topological node;

Described electric network data form is CIM/E form or CIM/XML form;

The process that reads and resolve described electric network data comprises: after data file is mapped in internal memory, internal memory is carried out to piecemeal according to different device types, adopt OpenMP multi-threaded parallel technology to carry out parallel processing to described each piece internal memory.

In the second preferred embodiment provided by the invention: in described step 2,

For single-ended equipment, described single-ended equipment is associated with described topological node;

For both-end or multiterminal equipment, after described both-end or multiterminal equipment are associated with described topological node, for topological node described in each generates described associated topologies node table;

Described single-ended equipment comprises bus, generator, load; Described both-end or multiterminal equipment comprise alternating current circuit, transformer.

In the 3rd preferred embodiment provided by the invention: adopt depth-first or BFS (Breadth First Search) method to travel through analysis to described associated topologies node table in described step 3; Adopt OpenMP multithreads computing to carry out topology search and numbering to described each subregion.

In the 4th preferred embodiment provided by the invention: the powernet data that generate described integration in described step 3 comprise: by the information of being concerned about topological node in region, alternating current circuit, transformer with PSD-BPA formatted output in flow data; If configured generator dynamic parameter storehouse, export the transient stability data of PSD-BPA form.

In the 5th preferred embodiment provided by the invention: carry out in described step 3, after described topology merging, comprising:

Step 4, carries out quality evaluation to described electric network data;

The process of carrying out described quality evaluation comprises the apparent error existing in the topological island of searching care, and mistake is revised, and adopts parallel processing mode to verify whether described electric network data meets power flow equation.

In the 6th preferred embodiment provided by the invention: verify in described step 4 whether described electric network data meets power flow equation and comprise:

Whether step 401 is 0 to judge whether described node meets power-balance according to the value of the meritorious imbalance power of the node calculating and idle imbalance power;

Step 402, the two-node model of structure circuit, judges according to the value of the output work rate variance at the two ends of the value of the output work rate variance at the two ends of the circuit calculating and the actual circuit recording whether described circuit meets trend checking;

Step 403, build the schematic equivalent circuit of transformer, the value of output work rate variance of converting current vector calculating transformer primary side, secondary side or three sides of primary side according to the transformer primary side current vector calculating and secondary side or three sides, judges according to the value of the value of the described output work rate variance calculating and the actual output work rate variance recording whether described transformer meets trend checking.

In the 7th preferred embodiment provided by the invention: in described step 401 according to meritorious imbalance power and the idle imbalance power of formula (1) and (2) computational analysis node:

$\mathrm{\Δ}{P}_i=\underset{m=1}{\overset{n_G}{\mathrm{\Σ}}}{P}_{\mathrm{Gm}}+\underset{m=1}{\overset{n_L}{\mathrm{\Σ}}}(-P_{\mathrm{Lm}})+\underset{m=1}{\overset{n_B}{\mathrm{\Σ}}}{P}_{\mathrm{Binm}}---\left(1\right);$

$\mathrm{\Δ}{Q}_i=\underset{m=1}{\overset{n_G}{\mathrm{\Σ}}}{Q}_{\mathrm{Gm}}+\underset{m=1}{\overset{n_L}{\mathrm{\Σ}}}(-Q_{\mathrm{Lm}})+\underset{m=1}{\overset{n_B}{\mathrm{\Σ}}}{Q}_{\mathrm{Binm}}+\underset{m=1}{\overset{n_{\mathrm{Sh}}}{\mathrm{\Σ}}}{Q}_{\mathrm{Shm}}---\left(2\right);$

Wherein, Δ P _iwith Δ Q _ithe meritorious imbalance power and the idle imbalance power that represent respectively node i; n _grepresent the generator sum in node i,

represent all generators in node i meritorious exert oneself and,

represent all generators in node i idle exert oneself and; n _lrepresent that the load in node i goes out number of lines,

represent all burden with power sums in node i,

represent all load or burden without work sums in node i; n _bthe branch road number that expression is connected with node i,

represent the meritorious sum of the injection of the connected branch road of node i in node i,

represent the injection idle sum of the connected branch road of node i in node i; n _shrepresent the compensation equipment number in node i,

the idle sum of exerting oneself that represents all reactive-load compensation equipments in node i;

In described step 402, can to circuit, carry out trend checking according to formula (3) and (4):

$\mathrm{\&Delta;}{\stackrel{.}{S}}_j={\stackrel{.}{S}}_j+{\stackrel{.}{V}}_j\hat{I}-\mathrm{<figure-callout\; id="1"\; label="j"\; filenames="HDA0000424891170000012.png"\; state="\{\{state\}\}">j</figure-callout>}\frac{1}{2}{V}_i^{2}B---\left(3\right);$

$\mathrm{\&Delta;}{\stackrel{.}{S}}_j={\stackrel{.}{S}}_j+{\stackrel{.}{V}}_j\hat{I}-\mathrm{<figure-callout\; id="1"\; label="j"\; filenames="HDA0000424891170000012.png"\; state="\{\{state\}\}">j</figure-callout>}\frac{1}{2}{V}_j^{2}B---\left(4\right);$

with

be respectively the output work rate variance at circuit two ends;

with

be respectively the output power at circuit two ends;

with

be respectively the voltage phasor at circuit two ends;

for the electric current phasor of circuit from i side to j side,

r+jX is the impedance of circuit,

be

conjugate phasors; B is the admittance of circuit.

In the 8th preferred embodiment provided by the invention: the transformer that carries out trend checking in described step 403 comprises two-winding transformer and three-winding transformer; When carrying out trend checking, described two-winding transformer verifies the value of the output work rate variance of described transformer primary side and secondary side; When carrying out trend checking, described three-winding transformer verifies the value of the output work rate variance of described transformer primary side, secondary side and three sides;

The process of value of calculating the output work rate variance of described two-winding transformer primary side and secondary side in described step 403 comprises:

According to formula (5)～(10), calculate the current vector that described two-winding transformer primary side current vector and secondary side are converted primary side

with

${\stackrel{.}{V}}_1={\stackrel{.}{I}}_1[(R_1+j{X}_1)+\frac{1}{G_m-{\mathrm{jB}}_m}]-{{\stackrel{.}{I}}^{\&prime;}}_2\frac{1}{G_m-j{B}_m}---\left(5\right);$

$-{\stackrel{.}{V}}_2={{\stackrel{.}{I}}^{\&prime;}}_2[({R^{\&prime;}}_2+j{X^{\&prime;}}_2)+\frac{1}{G_m-j{B}_m}]-{\stackrel{.}{I}}_1\frac{1}{G_m-j{B}_m}---\left(6\right);$

${{\stackrel{.}{V}}^{\&prime;}}_2=k{\stackrel{.}{V}}_2---\left(7\right);$

${{\stackrel{.}{I}}^{\&prime;}}_2=\frac{1}{k}{\stackrel{.}{I}}_2---\left(8\right);$

R ₂'＝k ²R ₂?（9）；

X ₂'＝k ²X ₂?（10）；

According to formula (11) and (12), calculate the output work rate variance of described two-winding transformer first and second side

with value:

$\mathrm{\&Delta;}{\stackrel{.}{S}}_1={\stackrel{.}{S}}_1-{\stackrel{.}{V}}_1{\hat{I}}_1---\left(11\right);$

$\mathrm{\&Delta;}{\stackrel{.}{S}}_2={\stackrel{.}{S}}_2-{{\stackrel{.}{V}}^{\&prime;}}_2{{\hat{I}}^{\&prime;}}_2---\left(12\right);$

for transformer primary side voltage phasor,

with be respectively

with

conjugate phasors,

for transformer primary side current phasor, R ₁+ jX ₁for the impedance of a winding of transformer,

for Circuit Fault on Secondary Transformer voltage phasor,

for Circuit Fault on Secondary Transformer electric current phasor, R ₂+ jX ₂for the impedance of transformer secondary winding, k is transformer one secondary side no-load voltage ratio,

for Circuit Fault on Secondary Transformer voltage phasor is converted the value to primary side,

for Circuit Fault on Secondary Transformer electric current phasor is converted the value to primary side, R ₂'+jX ₂' convert the value to primary side, Y for the impedance of transformer secondary winding _m=G _m-jB _mfor the excitation admittance of transformer,

be respectively the output power of transformer first and second side,

be respectively the output work rate variance of transformer first and second side.

The beneficial effect of the relatively immediate prior art of the present invention comprises:

1, the quick integration method of a kind of scale grid line data multi-wire journey provided by the invention, consider the feature of each step of Data Integration, by after the electric network data storage of reading, according to the difference of grid equipment type, internally deposit into row piecemeal, follow-up topology processing procedure associated and topological analysis is also to adopt multithreading to carry out parallel processing according to the internal memory of this piecemeal, improved the efficiency of integrating, processing procedure rapidly and efficiently, for other application of safety on line Stability Assessment and early warning lay the foundation, for the reliability of result provides foundation.

2, adopt OpenMP multi-threaded parallel technology to carry out parallel processing, guaranteed the rapidity of calculating, according to the feature of electric power networks nature piecemeal, by the thinking of electrical network piecemeal, multithreading scheme has been carried out to global design, stopped the collision problem of cross-thread.

3, increased the process of electric network data being carried out to quality evaluation, from the angle of basic circuit equation, consider the accurately fixed of online electric network data, index to whole order of accuarcy quantitative, for assessing the accuracy and the credible foundation that basic data angle is provided of other online senior application results.

Accompanying drawing explanation

Be illustrated in figure 1 the process flow diagram of the quick integration method of a kind of scale grid line data multi-wire journey provided by the invention;

Be illustrated in figure 2 the structural representation of the two-node model of circuit;

Be illustrated in figure 3 the schematic equivalent circuit of two-winding transformer.

Embodiment

With reference to the accompanying drawings the specific embodiment of the present invention is described in further detail below.

The invention provides the quick integration method of a kind of scale grid line data multi-wire journey, as shown in Figure 1, as shown in Figure 1, the method comprises its method flow diagram:

Step 1, reads and resolves all kinds of electric network datas, and this electric network data is deposited to internal memory, and the internal memory corresponding to the electric network data of dissimilar grid equipment carries out piecemeal.

Step 2, the grid equipment in the internal memory after adopting parallel processing mode to each piecemeal is associated with topological node, and topological node is set up to associated topologies node table.

Step 3, according to the piecemeal of internal memory, topological node corresponding to each piece internal memory carried out to subregion, adopt parallel processing mode to travel through analysis to the associated topologies node table of the topological node of each subregion, each topological node is carried out to topological island numbering, remove the dead island of topology and determine the topological island of being concerned about; Interconnection by topological node by stages carries out topology merging by the subregion of each topological node, generates the powernet data of integrating.

In step 1, according to the difference of grid equipment type, the internal memory of electric network data storage is carried out to piecemeal, follow-up step 2 is also to carry out block parallel processing according to the internal memory of this piecemeal with associated and processing procedure topological analysis of topology in 3, improved the efficiency of integrating, processing procedure fast, efficiently, accurately, for other application of safety on line Stability Assessment and early warning lay the foundation, for the reliability of result provides foundation.

Further, the electric network data reading in step 1 comprises: plant stand, bus, transformer, alternating current circuit, generator, load, shunt compensation equipment, series compensation equipment, transverter, DC line, topological node etc.This electric network data form can be CIM/E form or CIM/XML form.The process that reads and resolve electric network data comprises: after data file is mapped in internal memory, internal memory is carried out to piecemeal according to different device types, finally adopt OpenMP multi-threaded parallel technology to process each piece internal memory, to realize reading fast and resolving of data file.

Adopt OpenMP multi-threaded parallel technology to carry out parallel processing, guaranteed the rapidity of calculating, according to the feature of electric power networks nature piecemeal, by the thinking of electrical network piecemeal, multithreading scheme has been carried out to global design, stopped the collision problem of cross-thread.

In step 2, for single-ended equipment, this single-ended equipment is associated with topological node; For both-end or multiterminal equipment, after this both-end or multiterminal equipment are associated with topological node, be that each topological node generates associated topologies node table.Wherein single-ended equipment comprises bus, generator, load etc., and both-end or multiterminal equipment comprise alternating current circuit, transformer etc.

In step 3, associated topologies node table is traveled through to the process of analyzing and can adopt depth-first or BFS (Breadth First Search) method, adopt OpenMP multithreads computing to carry out topology search and numbering to each subregion.

The powernet data of generate integrating in step 3 comprise: by be concerned about the information such as topological node in region, alternating current circuit, transformer with PSD-BPA formatted output in flow data, if configured generator dynamic parameter storehouse, can export the transient stability data of PSD-BPA form.The data that obtain according to this can realize the functions such as trend computational analysis, transient stability computational analysis, calculation of short-circuit current analysis and dynamic stability computational analysis.

Preferably, in the step 3 of the quick integration method of a kind of scale grid line data multi-wire journey provided by the invention, carry out, after topology merging, can also comprising: step 4, electric network data is carried out to quality evaluation.The process of carrying out quality evaluation comprises the apparent error existing in the topological island of searching care, and mistake is revised, and adopts parallel processing mode to verify whether electric network data meets power flow equation.

Concrete, in step 4, verify whether electric network data meets power flow equation and comprise:

Whether step 401 is 0 to come decision node whether to meet power-balance according to the value of the meritorious imbalance power of the node calculating and idle imbalance power.

Step 402, the two-node model of structure circuit, judges according to the value of the output work rate variance at the two ends of the value of the output work rate variance at the two ends of the circuit calculating and the actual circuit recording whether circuit meets trend checking.

Step 403, build the schematic equivalent circuit of transformer, according to the value of output work rate variance that calculates the primary side current vector of transformer and secondary side or three sides and convert current vector calculating transformer primary side, secondary side or three sides of primary side, the value of the value of the output work rate variance calculating according to this and the actual output work rate variance recording judges whether this transformer meets trend checking.

Concrete, can be according to meritorious imbalance power and the idle imbalance power of formula (1) and (2) computational analysis node in step 401:

$\mathrm{\&Delta;}{P}_i=\underset{m=1}{\overset{n_G}{\mathrm{\&Sigma;}}}{P}_{\mathrm{Gm}}+\underset{m=1}{\overset{n_L}{\mathrm{\&Sigma;}}}(-P_{\mathrm{Lm}})+\underset{m=1}{\overset{n_B}{\mathrm{\&Sigma;}}}{P}_{\mathrm{Binm}}---\left(1\right);$

$\mathrm{\&Delta;}{Q}_i=\underset{m=1}{\overset{n_G}{\mathrm{\&Sigma;}}}{Q}_{\mathrm{Gm}}+\underset{m=1}{\overset{n_L}{\mathrm{\&Sigma;}}}(-Q_{\mathrm{Lm}})+\underset{m=1}{\overset{n_B}{\mathrm{\&Sigma;}}}{Q}_{\mathrm{Binm}}+\underset{m=1}{\overset{n_{\mathrm{Sh}}}{\mathrm{\&Sigma;}}}{Q}_{\mathrm{Shm}}---\left(2\right);$

Wherein, Δ P _iwith Δ Q _ithe meritorious imbalance power and the idle imbalance power that represent respectively node i; n _grepresent the generator sum in node i,

represent all generators in node i meritorious exert oneself and,

represent all generators in node i idle exert oneself and; n _lrepresent that the load in node i goes out number of lines, represent all burden with power sums in node i,

represent all load or burden without work sums in node i; n _brepresent the branch road number being connected with node i,

represent the meritorious sum of the injection of the connected branch road of node i in node i,

represent the injection idle sum of the connected branch road of node i in node i; n _shrepresent the compensation equipment number in node i,

represent all reactive-load compensation equipments in node i

$\mathrm{\&Delta;}{\stackrel{.}{S}}_j={\stackrel{.}{S}}_j+{\stackrel{.}{V}}_j\hat{I}-\mathrm{<figure-callout\; id="1"\; label="j"\; filenames="HDA0000424891170000012.png"\; state="\{\{state\}\}">j</figure-callout>}\frac{1}{2}{V}_i^{2}B---\left(3\right);$

$\mathrm{\&Delta;}{\stackrel{.}{S}}_j={\stackrel{.}{S}}_j+{\stackrel{.}{V}}_j\hat{I}-\mathrm{<figure-callout\; id="1"\; label="j"\; filenames="HDA0000424891170000012.png"\; state="\{\{state\}\}">j</figure-callout>}\frac{1}{2}{V}_j^{2}B---\left(4\right);$

Be illustrated in figure 2 the structural representation of the two-node model of circuit, as shown in Figure 2,

with

be respectively the output work rate variance at circuit two ends;

with

be respectively the output power at circuit two ends;

with be respectively the voltage phasor at circuit two ends;

for the electric current phasor of circuit from i side to j side,

r+jX is the impedance of circuit, be

conjugate phasors; B is the admittance of circuit; Wherein each value is per unit value.

The transformer that carries out trend checking in step 403 comprises two-winding transformer and three-winding transformer, when two-winding transformer carries out trend checking, need verify the value of the output work rate variance of transformer primary side and secondary side, when three-winding transformer carries out trend checking, need verify the value of the output work rate variance of transformer primary side, secondary side and three sides.

Take two-winding transformer as example, be illustrated in figure 3 the schematic equivalent circuit of two-winding transformer, in Fig. 3,

for transformer primary side voltage phasor,

for transformer primary side current phasor, R ₁+ jX ₁for the impedance of a winding of transformer,

for Circuit Fault on Secondary Transformer voltage phasor,

for Circuit Fault on Secondary Transformer electric current phasor, R ₂+ jX ₂for the impedance of transformer secondary winding, k is transformer one secondary side no-load voltage ratio,

for Circuit Fault on Secondary Transformer voltage phasor is converted the value to primary side,

for Circuit Fault on Secondary Transformer electric current phasor is converted the value to primary side, R ₂'+jX ₂' convert the value to primary side, Y for the impedance of transformer secondary winding _m=G _m-jB _mfor the excitation admittance of transformer,

be respectively the output power of transformer first and second side,

be respectively the output work rate variance of transformer first and second side, wherein each value is per unit value.

In step 403, according to formula (5)～(10), calculate the current vector that two-winding transformer primary side current vector and secondary side are converted primary side with

:

${\stackrel{.}{V}}_1={\stackrel{.}{I}}_1[({\mathrm{<figure-callout\; id="1"\; label="R"\; filenames="HDA0000424891170000012.png"\; state="\{\{state\}\}">R</figure-callout>}}_1+j{X}_1)+\frac{1}{G_m-{\mathrm{jB}}_m}]-{{\stackrel{.}{I}}^{\&prime;}}_2\frac{1}{G_m-j{B}_m}---\left(5\right);$

$-{\stackrel{.}{V}}_2={{\stackrel{.}{I}}^{\&prime;}}_2[({R^{\&prime;}}_2+j{X^{\&prime;}}_2)+\frac{1}{G_m-j{B}_m}]-{\stackrel{.}{I}}_1\frac{1}{G_m-j{B}_m}---\left(6\right);$

${{\stackrel{.}{V}}^{\&prime;}}_2=k{\stackrel{.}{V}}_2---\left(7\right);$

${{\stackrel{.}{I}}^{\&prime;}}_2=\frac{1}{k}{\stackrel{.}{I}}_2---\left(8\right);$

R ₂'＝k ²R ₂?（9）；

X ₂'＝k ²X ₂?（10）；

According to formula (11) and (12), carry out the output work rate variance of calculating transformer first and second side

with

value:

$\mathrm{\&Delta;}{\stackrel{.}{S}}_1={\stackrel{.}{S}}_1-{\stackrel{.}{V}}_1{\hat{I}}_1---\left(11\right);$

$\mathrm{\&Delta;}{\stackrel{.}{S}}_2={\stackrel{.}{S}}_2-{{\stackrel{.}{V}}^{\&prime;}}_2{{\hat{I}}^{\&prime;}}_2---\left(12\right);$

Wherein,

with be respectively

with conjugate phasors.

In step 4, verify whether each circuit-switched data in online data meets power flow equation, calculates the injecting power deviation of each branch power connected node, the branch road that output bias is larger and node, and provide quality of data scoring.

Finally should be noted that: above embodiment is only in order to illustrate that technical scheme of the present invention is not intended to limit, although the present invention is had been described in detail with reference to above-described embodiment, those of ordinary skill in the field are to be understood that: still can modify or be equal to replacement the specific embodiment of the present invention, and do not depart from any modification of spirit and scope of the invention or be equal to replacement, it all should be encompassed in the middle of claim scope of the present invention.

Claims (9)

1. the quick integration method of scale grid line data multi-wire journey, is characterized in that, described method comprises:

Step 1, reads and resolves all kinds of electric network datas, and described electric network data is deposited to internal memory, and the internal memory corresponding to the electric network data of dissimilar grid equipment carries out piecemeal;

Step 2, the grid equipment in the internal memory after adopting parallel processing mode to piecemeal described in each is associated with topological node, and described topological node is set up to associated topologies node table;

Step 3, carries out subregion according to the piecemeal of described internal memory to topological node corresponding to each piece internal memory, adopts parallel processing mode to travel through analysis to the described associated topologies node table of the described topological node of subregion described in each;

Topological node described in each is carried out to topological island numbering, remove the dead island of topology and determine the topological island of being concerned about;

Interconnection by topological node by stages carries out topology merging by the subregion of topological node described in each, generates the powernet data of integrating.

2. the method for claim 1, it is characterized in that, the described electric network data reading in described step 1 comprises: plant stand, bus, transformer, alternating current circuit, generator, load, shunt compensation equipment, series compensation equipment, transverter, DC line, topological node;

Described electric network data form is CIM/E form or CIM/XML form;

The process that reads and resolve described electric network data comprises: after data file is mapped in internal memory, internal memory is carried out to piecemeal according to different device types, adopt OpenMP multi-threaded parallel technology to carry out parallel processing to described each piece internal memory.

3. the method for claim 1, is characterized in that, in described step 2,

For single-ended equipment, described single-ended equipment is associated with described topological node;

For both-end or multiterminal equipment, after described both-end or multiterminal equipment are associated with described topological node, for topological node described in each generates described associated topologies node table;

Described single-ended equipment comprises bus, generator, load; Described both-end or multiterminal equipment comprise alternating current circuit, transformer.

4. the method for claim 1, is characterized in that, adopts depth-first or BFS (Breadth First Search) method to travel through analysis to described associated topologies node table in described step 3; Adopt OpenMP multithreads computing to carry out topology search and numbering to described each subregion.

5. the method for claim 1, is characterized in that, the powernet data that generate described integration in described step 3 comprise: by the information of being concerned about topological node in region, alternating current circuit, transformer with PSD-BPA formatted output in flow data; If configured generator dynamic parameter storehouse, export the transient stability data of PSD-BPA form.

6. the method for claim 1, is characterized in that, carries out, after described topology merging, comprising in described step 3:

Step 4, carries out quality evaluation to described electric network data;

The process of carrying out described quality evaluation comprises the apparent error existing in the topological island of searching care, and mistake is revised, and adopts parallel processing mode to verify whether described electric network data meets power flow equation.

7. method as claimed in claim 6, is characterized in that, verifies whether described electric network data meets power flow equation and comprise in described step 4:

Whether step 401 is 0 to judge whether described node meets power-balance according to the value of the meritorious imbalance power of the node calculating and idle imbalance power;

Step 402, the two-node model of structure circuit, judges according to the value of the output work rate variance at the two ends of the value of the output work rate variance at the two ends of the circuit calculating and the actual circuit recording whether described circuit meets trend checking;

Step 403, build the schematic equivalent circuit of transformer, the value of output work rate variance of converting current vector calculating transformer primary side, secondary side or three sides of primary side according to the transformer primary side current vector calculating and secondary side or three sides, judges according to the value of the value of the described output work rate variance calculating and the actual output work rate variance recording whether described transformer meets trend checking.

8. method as claimed in claim 7, is characterized in that,

In described step 401 according to meritorious imbalance power and the idle imbalance power of formula (1) and (2) computational analysis node:

$\mathrm{\&Delta;}{P}_i=\underset{m=1}{\overset{n_G}{\mathrm{\&Sigma;}}}{P}_{\mathrm{Gm}}+\underset{m=1}{\overset{n_L}{\mathrm{\&Sigma;}}}(-P_{\mathrm{Lm}})+\underset{m=1}{\overset{n_B}{\mathrm{\&Sigma;}}}{P}_{\mathrm{Binm}}---\left(1\right);$

$\mathrm{\&Delta;}{Q}_i=\underset{m=1}{\overset{n_G}{\mathrm{\&Sigma;}}}{Q}_{\mathrm{Gm}}+\underset{m=1}{\overset{n_L}{\mathrm{\&Sigma;}}}(-Q_{\mathrm{Lm}})+\underset{m=1}{\overset{n_B}{\mathrm{\&Sigma;}}}{Q}_{\mathrm{Binm}}+\underset{m=1}{\overset{n_{\mathrm{Sh}}}{\mathrm{\&Sigma;}}}{Q}_{\mathrm{Shm}}---\left(2\right);$

Wherein, Δ P _iwith Δ Q _ithe meritorious imbalance power and the idle imbalance power that represent respectively node i; n _grepresent the generator sum in node i, represent all generators in node i meritorious exert oneself and,

represent all generators in node i idle exert oneself and; n _lrepresent that the load in node i goes out number of lines,

represent all burden with power sums in node i,

represent all load or burden without work sums in node i; n _bthe branch road number that expression is connected with node i,

represent the meritorious sum of the injection of the connected branch road of node i in node i,

represent the injection idle sum of the connected branch road of node i in node i; n _shrepresent the compensation equipment number in node i,

the idle sum of exerting oneself that represents all reactive-load compensation equipments in node i;

In described step 402, can to circuit, carry out trend checking according to formula (3) and (4):

$\mathrm{\&Delta;}{\stackrel{.}{S}}_i={\stackrel{.}{S}}_i-{\stackrel{.}{V}}_i\hat{I}-j\frac{1}{2}{V}_i^{2}B---\left(3\right);$

$\mathrm{\&Delta;}{\stackrel{.}{S}}_j={\stackrel{.}{S}}_j+{\stackrel{.}{V}}_j\hat{I}-j\frac{1}{2}{V}_j^{2}B---\left(4\right);$

with

be respectively the output work rate variance at circuit two ends;

with be respectively the output power at circuit two ends;

with

be respectively the voltage phasor at circuit two ends;

for the electric current phasor of circuit from i side to j side,

r+jX is the impedance of circuit,

be

conjugate phasors; B is the admittance of circuit.

9. method as claimed in claim 7, is characterized in that, the transformer that carries out trend checking in described step 403 comprises two-winding transformer and three-winding transformer; When carrying out trend checking, described two-winding transformer verifies the value of the output work rate variance of described transformer primary side and secondary side; When carrying out trend checking, described three-winding transformer verifies the value of the output work rate variance of described transformer primary side, secondary side and three sides;

The process of value of calculating the output work rate variance of described two-winding transformer primary side and secondary side in described step 403 comprises:

According to formula (5)～(10), calculate the current vector that described two-winding transformer primary side current vector and secondary side are converted primary side with

${\stackrel{.}{V}}_1={\stackrel{.}{I}}_1[(R_1+j{X}_1)+\frac{1}{G_m-{\mathrm{jB}}_m}]-{{\stackrel{.}{I}}^{\&prime;}}_2\frac{1}{G_m-j{B}_m}---\left(5\right);$

$-{\stackrel{.}{V}}_2={{\stackrel{.}{I}}^{\&prime;}}_2[({R^{\&prime;}}_2+j{X^{\&prime;}}_2)+\frac{1}{G_m-j{B}_m}]-{\stackrel{.}{I}}_1\frac{1}{G_m-j{B}_m}---\left(6\right);$

${{\stackrel{.}{V}}^{\&prime;}}_2=k{\stackrel{.}{V}}_2---\left(7\right);$

${{\stackrel{.}{I}}^{\&prime;}}_2=\frac{1}{k}{\stackrel{.}{I}}_2---\left(8\right);$

R ₂'＝k ²R ₂?（9）；

X ₂'＝k ²X ₂?（10）；

According to formula (11) and (12), calculate the output work rate variance of described two-winding transformer first and second side with

value:

$\mathrm{\&Delta;}{\stackrel{.}{S}}_1={\stackrel{.}{S}}_1-{\stackrel{.}{V}}_1{\hat{I}}_1---\left(11\right);$

$\mathrm{\&Delta;}{\stackrel{.}{S}}_2={\stackrel{.}{S}}_2-{{\stackrel{.}{V}}^{\&prime;}}_2{{\hat{I}}^{\&prime;}}_2---\left(12\right);$

for transformer primary side voltage phasor, with

be respectively with

conjugate phasors,

for transformer primary side current phasor, R ₁+ jX ₁for the impedance of a winding of transformer,

for Circuit Fault on Secondary Transformer voltage phasor,

for Circuit Fault on Secondary Transformer electric current phasor, R ₂+ jX ₂for the impedance of transformer secondary winding, k is transformer one secondary side no-load voltage ratio, for Circuit Fault on Secondary Transformer voltage phasor is converted the value to primary side,

for Circuit Fault on Secondary Transformer electric current phasor is converted the value to primary side, R ₂'+jX ₂' convert the value to primary side, Y for the impedance of transformer secondary winding _m=G _m-jB _mfor the excitation admittance of transformer,

be respectively the output power of transformer first and second side,

be respectively the output work rate variance of transformer first and second side.

Images