CN105373834A - Power distribution network short-circuit current calculating method and system based on distributed calculation - Google Patents

Abstract

The invention provides a power distribution network short-circuit current calculating method based on distributed calculation. The method comprises the steps of: constructing a distributed parallel calculating platform; obtaining and analyzing data provided by an SCADA system and XML data of a CIM model through the platform, and then determining specific parameters of a power distribution network; converting a physical model of the power distribution network into a calculating module; judging a failure type and calculating circuit parameters during the failure; utilizing a forward-backward sweep method to calculate circuit information in the normal operation state; and superposing the circuit in the normal operation state with the circuit in the failure, and obtaining a fault current, currents of any branches and voltages of any nodes; analyzing and applying the result of the short-circuit calculation; and based on a multi-port theory and forward-backward sweep tidal current calculation, judging and concluding the fault type, and carrying out current-circuit calculation. The method has the advantages that the time is short, and requirements of real-time online circuit calculation are met.

Description

Based on power distribution network short-circuit current calculation method and the system of Distributed Calculation

Technical field

The present invention relates to system for distribution network of power calculation of short-circuit current field, be specifically related to a kind of power distribution network short-circuit current calculation method based on Distributed Calculation and system.

Background technology

Power distribution network short circuit calculation is the work must carried out in distribution network system planning, design and running.And there is following shortcoming in power distribution network short circuit calculation at present:

The scale of power distribution network own is large, node is many, and the odds ratio of branch parameters R/X is comparatively large, and branch line is more, and therefore the time of short circuit calculation is comparatively long, can not meet the requirement of real-time online short circuit calculation;

For the short circuit calculation method of power distribution network feature as phase region analytic approach or take expert system to calculate the method for short-circuit fault of power system, although principle comparatively simple, intuitive, multiple failure cannot be processed.

Summary of the invention

For the problems referred to above, the invention provides one mainly theoretical based on obstacle mouth, simple fault and complex fault unification are processed, derive applicable computer implemented pattern, and adopt middleware Technology, huge power distribution network is decomposed into some subsystems to be distributed to different computing machines and to realize parallel computation, based on power distribution network short-circuit current calculation method and the system of Distributed Calculation.

The technical scheme that the present invention adopts for its technical matters of solution is:

Based on the power distribution network short-circuit current calculation method of Distributed Calculation, it is characterized in that comprising the following steps:

1) step of distributed paralleling calculation platform structure;

2) step of data processing, is obtained by platform and resolves the XML data of data that SCADA system provides and CIM, thus determining the design parameter of power distribution network;

3) step of topological analysis, is converted into computation model by the physical model of distribution;

4) step of short circuit calculation:

Failure judgement type circuit parameter when calculating fault, the voltage of arbitrary node and the fault component of any branch current when this fault parameter comprises fault;

Circuit information when pushing manipulation calculates normal operating condition is pushed back before utilization;

State when normally running is superposed with circuit during fault, obtains fault current and the arbitrarily electric current of branch road and the voltage of arbitrary node, and then draw trouble spot short-circuit current;

5) step of short-circuit analysis, carries out analysis and validation to the result of short circuit calculation, and by the result of short circuit calculation and analysis result encapsulation backward on provide or put preservation in storage.

Further, described step 2) in the design parameter of power distribution network comprise balance node voltage, the injecting power of each node and each branch impedance.

Further, described step 4) in when calculating fault the method for circuit parameter be the computing method based on fault mouth theory, comprise the calculation procedure of simple fault circuit parameter and the calculation procedure of multiple failure circuit parameter.

Further, the calculation procedure of described simple fault circuit parameter is:

By the fault mouth equation of each independent sequence net and boundary condition simultaneous, obtain electric current and the information of voltage of fault mouth, and then obtain the voltage of arbitrary node and the electric current of any branch road, the voltage equation of fault mouth is:

U _p＝Z _p·I _p(1)

Wherein, U _pfor the equivalent voltage of fault mouth,

I _pfor the Injection Current of entering viewed from fault mouth, i.e. fault current;

Z _pfor network is contracted to the equiva lent impedance of fault mouth;

Wherein, the boundary condition of single-line to ground fault is:

I ⁽¹⁾＝I ⁽²⁾＝I ⁽⁰⁾

U ⁽¹⁾+U ⁽²⁾+U ⁽⁰⁾＝0(2)

Line to line fault and line to line fault ground connection boundary condition be:

U ⁽¹⁾＝U ⁽²⁾＝U ⁽⁰⁾

I ⁽¹⁾+I ⁽²⁾+I ⁽⁰⁾＝0(3)

In formula, I ⁽¹⁾, I ⁽²⁾, I ⁽⁰⁾three order components of short-circuit current respectively, U ⁽¹⁾, U ⁽²⁾, U ⁽⁰⁾three sequence voltages of short dot.

Further, the calculation procedure of described multiple failure circuit parameter is:

Application twoport theory writes out the sequence network interface parametric equation of each fault mouth, and the boundary condition of simultaneous each fault place order components, solves the voltage and current obtaining each fault mouth, and then try to achieve the electric current of any branch road and the voltage of arbitrary node simultaneously;

The sequence network parameters equation of each fault mouth is:

Wherein, the voltage of each sequence with electric current column vector be respectively three sequence voltages, the current component of each fault mouth; Matrix of coefficients Z is the impedance matrix that whole network is contracted to fault mouth;

When each class fault only takes into account metallic earthing, the boundary condition of unified pattern can be drawn, wherein

The boundary condition of single-line to ground fault and two-phase disconnection fault is:

I ⁽¹⁾＝I ⁽²⁾＝I ⁽⁰⁾

U ⁽¹⁾+U ⁽²⁾+U ⁽⁰⁾＝0(7)

The boundary condition of two-phase grounding fault, line to line fault and single-phase wire break is:

U ⁽¹⁾＝U ⁽²⁾＝U ⁽⁰⁾

I ⁽¹⁾+I ⁽²⁾+I ⁽⁰⁾＝0(8)

In formula, I ⁽¹⁾, I ⁽²⁾, I ⁽⁰⁾three order components of short-circuit current respectively, U ⁽¹⁾, U ⁽²⁾, U ⁽⁰⁾three sequence voltages of short dot.

Further, described step 4) utilize before the step of circuit information that pushes back when pushing manipulation calculates normal operating condition comprise: computing node Injection Current, push back calculate each branch current and before inquire into and separate each node voltage.

Further, the formula of described computing node Injection Current is:

${\left[\begin{array}{c}{\stackrel{\·}{I}}_{ja}\\ {\stackrel{\·}{I}}_{jb}\\ {\stackrel{\·}{I}}_{jc}\end{array}\right]}^{\left(k\right)}=\left[\begin{array}{c}{\left(\frac{S_{ja}}{u_{ja}^{k-1}}\right)}^{*}\\ {\left(\frac{S_{ja}}{u_{jb}^{k-1}}\right)}^{*}\\ {\left(\frac{S_{ja}}{u_{jc}^{k-1}}\right)}^{*}\end{array}\right]+\left[\begin{array}{ccc}{Y}_{ja}& & \\ & Y_{jb}& \\ & & Y_{jc}\end{array}\right]{\left[\begin{array}{c}{U}_{ja}\\ U_{jb}\\ U_{jc}\end{array}\right]}^{\left(k-1\right)}---\left(9\right)$

Wherein, in formula, (k) represents iterations, for the three-phase component of node j Injection Current; Y _ja, Y _jb, Y _jcit is the three-phase component of node self-admittance; S _ja, S _ja, S _jait is the three-phase component of node injecting power; it is the three-phase component of kth time iteration node voltage.

Further, pushing back the method calculating each branch current described in is that, to root node recursion from last one deck branch road, the electric current of branch road l is:

${\left[\begin{array}{c}\stackrel{\·}{I_{la}^{\′}}\\ \stackrel{\·}{I_{lb}^{\′}}\\ \stackrel{\·}{I_{lc}^{\′}}\end{array}\right]}^{\left(k\right)}={\left[\begin{array}{c}\stackrel{\·}{I_{ja}}\\ \stackrel{\·}{I_{jb}}\\ \stackrel{\·}{I_{jc}}\end{array}\right]}^{\left(k\right)}+\underset{x\∈X}{\Σ}{\left[\begin{array}{c}\stackrel{\·}{I_{xa}}\\ \stackrel{\·}{I_{xb}}\\ \stackrel{\·}{I_{xc}}\end{array}\right]}^{\left(k\right)}---\left(10\right)$

Then branch road l top electric current is:

${\left[\begin{array}{c}{\stackrel{\·}{I}}_{la}\\ {\stackrel{\·}{I}}_{lb}\\ {\stackrel{\·}{I}}_{lc}\end{array}\right]}^{\left(k\right)}=\frac{1}{2}\left[\begin{array}{ccc}{Y}_{aa}& Y_{ab}& Y_{ac}\\ Y_{ba}& Y_{bb}& Y_{bc}\\ Y_{ca}& Y_{cb}& Y_{cc}\end{array}\right]\left({\left[\begin{array}{c}{\stackrel{\·}{U}}_{ia}\\ {\stackrel{\·}{U}}_{ib}\\ {\stackrel{\·}{U}}_{ic}\end{array}\right]}^{\left(k-1\right)}+{\left[\begin{array}{c}{\stackrel{\·}{U}}_{ja}\\ {\stackrel{\·}{U}}_{jb}\\ {\stackrel{\·}{U}}_{jc}\end{array}\right]}^{\left(k-1\right)}\right)+{\left[\begin{array}{c}{\stackrel{\·}{I}}_{la}^{\′}\\ {\stackrel{\·}{I}}_{lb}^{\′}\\ {\stackrel{\·}{I}}_{lc}^{\′}\end{array}\right]}^{\left(k\right)}---\left(11\right)$

Wherein, in formula, (k) represents iterations, for the end current of branch road l; for the top electric current of branch road l; X is the set of all branches be connected with node j.

Further, inquiring into the method for separating each node voltage before described is, for distribution network feeder branch road, from root node, terminad advances, then the voltage of node j is:

${\left[\begin{array}{c}\stackrel{\·}{U_{ja}}\\ \stackrel{\·}{U_{jb}}\\ \stackrel{\·}{U_{jc}}\end{array}\right]}^{\left(k\right)}={\left[\begin{array}{c}\stackrel{\·}{U_{ia}}\\ \stackrel{\·}{U_{ib}}\\ \stackrel{\·}{U_{ic}}\end{array}\right]}^{\left(k\right)}-\left[\begin{array}{ccc}{Z}_{aa}& Z_{ab}& Z_{ac}\\ Z_{ba}& Z_{bb}& Z_{bc}\\ Z_{ca}& Z_{cb}& Z_{cc}\end{array}\right]\left({\left[\begin{array}{c}\stackrel{\·}{I_{la}}\\ \stackrel{\·}{I_{lb}}\\ \stackrel{\·}{I_{lc}}\end{array}\right]}^{\left(k\right)}-\frac{1}{2}\left[\begin{array}{ccc}{Y}_{aa}& Y_{ab}& Y_{ac}\\ Y_{ab}& Y_{bb}& Y_{bc}\\ Y_{ac}& Y_{cb}& Y_{cc}\end{array}\right]{\left[\begin{array}{c}\stackrel{\·}{U_{ia}}\\ \stackrel{\·}{U_{ib}}\\ \stackrel{\·}{U_{ic}}\end{array}\right]}^{\left(k\right)}\right)---\left(12\right)$

Based on the power distribution network calculation of short-circuit current system of Distributed Calculation, it is characterized in that comprising:

Data processing module, one is that two are to provide the process to result of calculation, are packaged into the form of standard for being resolved by the carrying out obtaining data from platform;

Topological analysis module is the basis of short circuit calculation, for the physical model of distribution is converted into computation model;

Short circuit calculation module, for realizing the calculating of short-circuit current on distributed platform;

Results analyses module, for carrying out analysis and validation to the result of short circuit calculation, and by the result of calculation of short circuit calculation module and results analyses module or meet at data processing module encapsulation backward on provide or put preservation in storage.

Further, the described data obtained from platform comprise the XML data of data that SCADA system provides, CIM, and these data are transmitted by the CIS interface of standard.

Further, the mode of described topology is Global Topological or dynamic local topology.

The invention has the beneficial effects as follows: the power distribution network short-circuit current calculation method based on Distributed Calculation of the present invention and system, the boundary condition of all kinds of fault type is carried out unifying and conclusion, the new method being applicable to computing machine is proposed, be convenient to computer programming, no special fault and multiple failure can be processed simultaneously;

In addition, theoretical based on multiport, conclude fault type in conjunction with before pushing back Load flow calculation to judge, and carry out calculation of short-circuit current, the time is short, can meet the requirement of real-time online short circuit calculation;

In addition, Task-decomposing becomes some subtasks to be distributed to different computing machine to process simultaneously, improve counting yield by Distributed Calculation.

Accompanying drawing explanation

Be further detailed below in conjunction with the drawings and specific embodiments:

Fig. 1 is the system global structure block diagram of the embodiment of the present invention;

Fig. 2 is the computing method overview flow chart of the embodiment of the present invention;

Fig. 3 is the basic computational ele-ment schematic diagram of the radial distribution networks of the embodiment of the present invention;

Fig. 4 is the distributed paralleling calculation platform general frame figure of the embodiment of the present invention;

Fig. 5 is the communications framework figure of the distributed paralleling calculation platform of the embodiment of the present invention.

Embodiment

Below with reference to embodiment and accompanying drawing, design of the present invention and technique effect are clearly and completely described, to understand object of the present invention, characteristic sum effect fully.Obviously; described embodiment is a part of embodiment of the present invention, instead of whole embodiment, based on embodiments of the invention; other embodiments that those skilled in the art obtains under the prerequisite not paying creative work, all belong to the scope of protection of the invention.

With reference to Fig. 2, based on the power distribution network short-circuit current calculation method of Distributed Calculation, it is characterized in that comprising the following steps:

1) step of distributed paralleling calculation platform structure;

2) step of data processing, is obtained by platform and resolves the XML data of data that SCADA system provides and CIM, thus determining the design parameter of power distribution network;

3) step of topological analysis, is converted into computation model by the physical model of distribution;

4) step of short circuit calculation:

Failure judgement type circuit parameter when calculating fault, the voltage of arbitrary node and the fault component of any branch current when this fault parameter comprises fault;

Circuit information when pushing manipulation calculates normal operating condition is pushed back before utilization;

State when normally running is superposed with circuit during fault, obtains fault current and the arbitrarily electric current of branch road and the voltage of arbitrary node, and then draw trouble spot short-circuit current;

5) step of short-circuit analysis, carries out analysis and validation to the result of short circuit calculation, and by the result of short circuit calculation and analysis result encapsulation backward on provide or put preservation in storage.

Wherein, the first two step part sequencing in the step of short circuit calculation, can carry out simultaneously

With further reference to Fig. 4, introduce the constructing plan of present pre-ferred embodiments distributed paralleling calculation platform, this distributed paralleling calculation platform system comprises 3 parts: PAS client, server end and actuating station.Client can manage, as the operation such as start and stop, forbidding each example of each server end; Server end (has two, and active and standby each other) mainly comprise administration module and example, administration module primary responsibility contact customer end and example, and realize load balancing, example then can realize Task-decomposing, and the data required for actuating station preparation calculating; Actuating station is responsible for performing all kinds of calculation task, and result of calculation is returned to corresponding example.

After each module of system starts, client sends a task requests to administration module, is inquired about the example whether having startup by administration module, if do not have, then manually boots an example by client.Each example safeguards a Fastdb database separately, and the data between example are separate, deposits the CIM of electrical network and dynamic SACDA data.The data (as failure message, optimum configurations etc.) that client wants required by task are write in the oracle database of external memory, and task number is sent to example, and example reads mission bit stream after receiving mission number from Oracle.Example carries out Task-decomposing after receiving all information of task, and all subtasks of having decomposed are sent to administration module, distributes subtask by administration module load balancing to each actuating station.Result is directly returned to example after calculating and terminating by actuating station, and example is collected result of calculation and write in oracle database, and notifies that client task completes.Like this, client just can search the result of calculation required for oneself from database.

In order to realize the data transmission of this platform, employing global communication is the key of distributed implementation, and the present embodiment adopts ZeroMQ as the global communication framework of whole Distributed Computing Platform, as shown in Figure 5.According to the design proposal of distributed platform system, concrete communication can be divided into the following aspects:

(1) data communication of external system and example

First example is connected with port with the IP of external system after starting, and by socket REQ externally system request initialization, all data can be returned by REP by external system.Meanwhile, example also needs by PUB-SUB socket, and the dynamic data of " subscription " external system " issue " upgrades.

(2) communication of PAS client and administration module

By socket REQ-REP alternately, PAS client obtains the information of example, as available example list, the IP at example place and the port numbers etc. of binding from administration module for PAS client and administration module.

(3) PAS client and example is mutual

PAS client can be connected with a certain example after obtaining the information of example, sends task requests by socket DEALER-ROUTER socket.

(4) actuating station and administration module, example mutual

First actuating station will be connected with administration module by socket REQ-REP after starting, and from administration module, obtain all example information, and actuating station creates one group of group socket successively according to example information and is connected to example.

Actuating station after calculating and terminating needs result to return to example.This function is realized by PULL-PUSH socket.

By building above-mentioned distributed paralleling calculation platform system and communication mode, obtain and resolve the XML data of data that SCADA system provides and CIM, thus the concrete steps calculating short circuit after determining the design parameter of power distribution network are as follows:

Failure judgement type circuit parameter when calculating fault,

When fault is simple fault, by the fault mouth equation of each independent sequence net and boundary condition simultaneous, obtain electric current and the information of voltage of fault mouth, and then obtain the voltage of arbitrary node and the electric current of any branch road, the voltage equation of fault mouth is:

U _p＝Z _p·I _p(1)

Wherein, U _pfor the equivalent voltage of fault mouth, this voltage source is unique voltage source of network, I _pfor the Injection Current of entering viewed from fault mouth, i.e. fault current; Z _pfor network is contracted to the equiva lent impedance of fault mouth.

When solving the voltage equation of fault mouth, need the boundary condition by short-circuit-type, except symmetry short circuit, all the other short-circuit-type can be summarized as two kinds of forms:

Single-line to ground fault is:

I ⁽¹⁾＝I ⁽²⁾＝I ⁽⁰⁾

U ⁽¹⁾+U ⁽²⁾+U ⁽⁰⁾＝0(2)

Line to line fault and line to line fault ground connection are:

U ⁽¹⁾＝U ⁽²⁾＝U ⁽⁰⁾

I ⁽¹⁾+I ⁽²⁾+I ⁽⁰⁾＝0(3)

In formula, I ⁽¹⁾, I ⁽²⁾, I ⁽⁰⁾three order components of short-circuit current respectively, U ⁽¹⁾, U ⁽²⁾, U ⁽⁰⁾three sequence voltages of short dot.

When fault is multiple failure, application twoport theory writes out the sequence network interface parametric equation of each fault mouth, the boundary condition of simultaneous each fault place order components, solves the voltage and current obtaining each fault mouth, and then tries to achieve the electric current of any branch road and the voltage of arbitrary node simultaneously.

The sequence network parameters equation of each fault mouth is:

Wherein, the voltage of each sequence with electric current column vector be respectively three sequence voltages, the current component of each fault mouth; Matrix of coefficients Z is the impedance matrix that whole network is contracted to fault mouth.

For each class fault (comprising short circuits and open conductors), when only taking into account metallic earthing, the boundary condition of following unified pattern can be write out:

Single-line to ground fault and two-phase disconnection fault:

I ⁽¹⁾＝I ⁽²⁾＝I ⁽⁰⁾

U ⁽¹⁾+U ⁽²⁾+U ⁽⁰⁾＝0(7)

Two-phase grounding fault, line to line fault and single-phase wire break:

U ⁽¹⁾＝U ⁽²⁾＝U ⁽⁰⁾

I ⁽¹⁾+I ⁽²⁾+I ⁽⁰⁾＝0(8)

The derivation result of various fault type combination is as shown in table 1:

The computing formula of each fault type combination of table 1

Table 1 gives matrix equation when solving twin failure, utilizes above equation can be easy to solve each sequence electric current and the voltage of fault mouth, thus obtains the voltage of arbitrary node and the electric current of any branch road easily.

Further, push back pushing manipulation before using and carry out Load flow calculation, branch road and nodal information can be obtained, namely state during normal operation, superpose with circuit during fault, thus obtain fault current and the arbitrarily electric current of branch road and the voltage of arbitrary node, the front pushing manipulation that pushes back is respectively by pushing back and front pushing away two processes and solve branch current and node voltage, as shown in Figure 3, concrete computation process is the basic computational ele-ment of radial distribution networks:

1. computing node Injection Current:

${\left[\begin{array}{c}{\stackrel{\·}{I}}_{ja}\\ {\stackrel{\·}{I}}_{jb}\\ {\stackrel{\·}{I}}_{jc}\end{array}\right]}^{\left(k\right)}=\left[\begin{array}{c}{\left(\frac{S_{ja}}{u_{ja}^{k-1}}\right)}^{*}\\ {\left(\frac{S_{ja}}{u_{jb}^{k-1}}\right)}^{*}\\ {\left(\frac{S_{ja}}{u_{jc}^{k-1}}\right)}^{*}\end{array}\right]+\left[\begin{array}{ccc}{Y}_{ja}& & \\ & Y_{jb}& \\ & & Y_{jc}\end{array}\right]{\left[\begin{array}{c}{U}_{ja}\\ U_{jb}\\ U_{jc}\end{array}\right]}^{\left(k-1\right)}---\left(9\right)$

Wherein, in formula, (k) represents iterations, for the three-phase component of node j Injection Current; Y _ja, Y _jb, Y _jcit is the three-phase component of node self-admittance; S _ja, S _ja, S _jait is the three-phase component of node injecting power; it is the three-phase component of kth time iteration node voltage.

2. push back and calculate each branch current:

To root node recursion from last one deck branch road, the electric current of branch road l is:

${\left[\begin{array}{c}\stackrel{\·}{I_{la}^{\′}}\\ \stackrel{\·}{I_{lb}^{\′}}\\ \stackrel{\·}{I_{lc}^{\′}}\end{array}\right]}^{\left(k\right)}={\left[\begin{array}{c}\stackrel{\·}{I_{ja}}\\ \stackrel{\·}{I_{jb}}\\ \stackrel{\·}{I_{jc}}\end{array}\right]}^{\left(k\right)}+\underset{x\∈X}{\Σ}{\left[\begin{array}{c}\stackrel{\·}{I_{xa}}\\ \stackrel{\·}{I_{xb}}\\ \stackrel{\·}{I_{xc}}\end{array}\right]}^{\left(k\right)}---\left(10\right)$

Then branch road l top electric current is:

${\left[\begin{array}{c}{\stackrel{\·}{I}}_{la}\\ {\stackrel{\·}{I}}_{lb}\\ {\stackrel{\·}{I}}_{lc}\end{array}\right]}^{\left(k\right)}=\frac{1}{2}\left[\begin{array}{ccc}{Y}_{aa}& Y_{ab}& Y_{ac}\\ Y_{ba}& Y_{bb}& Y_{bc}\\ Y_{ca}& Y_{cb}& Y_{cc}\end{array}\right]\left({\left[\begin{array}{c}{\stackrel{\·}{U}}_{ia}\\ {\stackrel{\·}{U}}_{ib}\\ {\stackrel{\·}{U}}_{ic}\end{array}\right]}^{\left(k-1\right)}+{\left[\begin{array}{c}{\stackrel{\·}{U}}_{ja}\\ {\stackrel{\·}{U}}_{jb}\\ {\stackrel{\·}{U}}_{jc}\end{array}\right]}^{\left(k-1\right)}\right)+{\left[\begin{array}{c}{\stackrel{\·}{I}}_{la}^{\′}\\ {\stackrel{\·}{I}}_{lb}^{\′}\\ {\stackrel{\·}{I}}_{lc}^{\′}\end{array}\right]}^{\left(k\right)}---\left(11\right)$

Wherein, in formula, (k) represents iterations, for the end current of branch road l; for the top electric current of branch road l; X is the set of all branches be connected with node j.

3. before inquire into and separate each node voltage

For distribution network feeder branch road, from root node, terminad advances, then the voltage of node j is:

${\left[\begin{array}{c}\stackrel{\·}{U_{ja}}\\ \stackrel{\·}{U_{jb}}\\ \stackrel{\·}{U_{jc}}\end{array}\right]}^{\left(k\right)}={\left[\begin{array}{c}\stackrel{\·}{U_{ia}}\\ \stackrel{\·}{U_{ib}}\\ \stackrel{\·}{U_{ic}}\end{array}\right]}^{\left(k\right)}-\left[\begin{array}{ccc}{Z}_{aa}& Z_{ab}& Z_{ac}\\ Z_{ba}& Z_{bb}& Z_{bc}\\ Z_{ca}& Z_{cb}& Z_{cc}\end{array}\right]\left({\left[\begin{array}{c}\stackrel{\·}{I_{la}}\\ \stackrel{\·}{I_{lb}}\\ \stackrel{\·}{I_{lc}}\end{array}\right]}^{\left(k\right)}-\frac{1}{2}\left[\begin{array}{ccc}{Y}_{aa}& Y_{ab}& Y_{ac}\\ Y_{ab}& Y_{bb}& Y_{bc}\\ Y_{ac}& Y_{cb}& Y_{cc}\end{array}\right]{\left[\begin{array}{c}\stackrel{\·}{U_{ia}}\\ \stackrel{\·}{U_{ib}}\\ \stackrel{\·}{U_{ic}}\end{array}\right]}^{\left(k\right)}\right)---\left(12\right)$

Then, state when normally running is superposed with circuit during fault, obtain fault current and the arbitrarily electric current of branch road and the voltage of arbitrary node, and then draw trouble spot short-circuit current.

With reference to Fig. 1, based on the power distribution network calculation of short-circuit current system of Distributed Calculation, comprising:

Data processing module, one is for being resolved by the carrying out obtaining data from platform, the data provided as SCADA system, the XML data of CIM, these data import native system into by the CIS interface of standard; Two are to provide the process to result of calculation, are packaged into the form of standard, return user by standard C IS interface;

Topological analysis module is the basis of short circuit calculation, for the physical model of distribution is converted into computation model, both can Global Topological, and also can dynamic local topology;

Short circuit calculation module is the nucleus module of this system, for realizing the calculating of short-circuit current on distributed platform;

Results analyses module, for carrying out analysis and validation to the result of short circuit calculation,, as rupturing capacity scanning, the adjusting of arc suppression coil, and by the result of calculation of short circuit calculation module and results analyses module or meet at data processing module encapsulation backward on provide or put preservation in storage.

Embodiments of the present invention are not restricted to the described embodiments, as long as it reaches technique effect of the present invention with substantially identical means, all should belong to protection scope of the present invention.

Claims (12)

1., based on the power distribution network short-circuit current calculation method of Distributed Calculation, it is characterized in that comprising with lower part:

1) step of distributed paralleling calculation platform structure;

2) step of data processing, is obtained by platform and resolves the XML data of data that SCADA system provides and CIM, thus determining the design parameter of power distribution network;

3) step of topological analysis, is converted into computation model by the physical model of distribution;

4) step of short circuit calculation:

Failure judgement type circuit parameter when calculating fault, the voltage of arbitrary node and the fault component of any branch current when this fault parameter comprises fault;

Circuit information when pushing manipulation calculates normal operating condition is pushed back before utilization;

State when normally running is superposed with circuit during fault, obtains fault current and the arbitrarily electric current of branch road and the voltage of arbitrary node, and then draw trouble spot short-circuit current;

5) step of short-circuit analysis, carries out analysis and validation to the result of calculation of short-circuit current, and by the result of calculation of short-circuit current and analysis result encapsulation backward on provide or put preservation in storage.

2. the power distribution network short-circuit current calculation method based on Distributed Calculation according to claim 1, is characterized in that: described step 2) in the design parameter of power distribution network comprise balance node voltage, the injecting power of each node and each branch impedance.

3. the power distribution network short-circuit current calculation method based on Distributed Calculation according to claim 1, it is characterized in that: described step 4) in when calculating fault the method for circuit parameter be computing method based on fault mouth theory, comprise the calculation procedure of simple fault circuit parameter and the calculation procedure of multiple failure circuit parameter.

4. the power distribution network short-circuit current calculation method based on Distributed Calculation according to claim 3, is characterized in that: the calculation procedure of described simple fault circuit parameter is:

By the fault mouth equation of each independent sequence net and boundary condition simultaneous, obtain electric current and the information of voltage of fault mouth, and then obtain the voltage of arbitrary node and the electric current of any branch road, the voltage equation of fault mouth is:

U _p＝Z _p·I _p(1)

Wherein, U _pfor the equivalent voltage of fault mouth,

I _pfor the Injection Current of entering viewed from fault mouth, i.e. fault current;

Z _pfor network is contracted to the equiva lent impedance of fault mouth;

Wherein, the boundary condition of single-line to ground fault is:

I ⁽¹⁾＝I ⁽²⁾＝I ⁽⁰⁾

U ⁽¹⁾+U ⁽²⁾+U ⁽⁰⁾＝0(2)

Line to line fault and line to line fault ground connection boundary condition be:

U ⁽¹⁾＝U ⁽²⁾＝U ⁽⁰⁾

I ⁽¹⁾+I ⁽²⁾+I ⁽⁰⁾＝0(3)

In formula, I ⁽¹⁾, I ⁽²⁾, I ⁽⁰⁾three order components of short-circuit current respectively, U ⁽¹⁾, U ⁽²⁾, U ⁽⁰⁾three sequence voltages of short dot.

5. the power distribution network short-circuit current calculation method based on Distributed Calculation according to claim 3, is characterized in that: the calculation procedure of described multiple failure circuit parameter is:

Application twoport theory writes out the sequence network interface parametric equation of each fault mouth, and the boundary condition of simultaneous each fault place order components, solves the voltage and current obtaining each fault mouth, and then try to achieve the electric current of any branch road and the voltage of arbitrary node simultaneously;

The sequence network parameters equation of each fault mouth is:

Wherein, the voltage of each sequence with electric current column vector be respectively three sequence voltages, the current component of each fault mouth; Matrix of coefficients Z is the impedance matrix that whole network is contracted to fault mouth;

When each class fault only takes into account metallic earthing, the boundary condition of unified pattern can be drawn, wherein

The boundary condition of single-line to ground fault and two-phase disconnection fault is:

I ⁽¹⁾＝I ⁽²⁾＝I ⁽⁰⁾

U ⁽¹⁾+U ⁽²⁾+U ⁽⁰⁾＝0(7)

The boundary condition of two-phase grounding fault, line to line fault and single-phase wire break is:

U ⁽¹⁾＝U ⁽²⁾＝U ⁽⁰⁾

I ⁽¹⁾+I ⁽²⁾+I ⁽⁰⁾＝0(8)

In formula, I ⁽¹⁾, I ⁽²⁾, I ⁽⁰⁾three order components of short-circuit current respectively, U ⁽¹⁾, U ⁽²⁾, U ⁽⁰⁾three sequence voltages of short dot.

6. the power distribution network short-circuit current calculation method based on Distributed Calculation according to claim 1, is characterized in that: described step 4) in utilize before the step of circuit information that pushes back when pushing manipulation calculates normal operating condition comprise: computing node Injection Current, push back calculate each branch current and before inquire into and separate each node voltage.

7. the power distribution network short-circuit current calculation method based on Distributed Calculation according to claim 6, is characterized in that: the formula of described computing node Injection Current is:

Wherein, in formula, (k) represents iterations, for the three-phase component of node j Injection Current; Y _ja, Y _jb, Y _jcit is the three-phase component of node self-admittance; S _ja, S _ja, S _jait is the three-phase component of node injecting power; it is the three-phase component of kth time iteration node voltage.

8. the power distribution network short-circuit current calculation method based on Distributed Calculation according to claim 6, is characterized in that: described in push back the method calculating each branch current be that, to root node recursion from last one deck branch road, the electric current of branch road l is:

Then branch road l top electric current is:

Wherein, in formula, (k) represents iterations, for the end current of branch road l; for the top electric current of branch road l; X is the set of all branches be connected with node j.

9. the power distribution network short-circuit current calculation method based on Distributed Calculation according to claim 6, it is characterized in that: inquiring into the method for separating each node voltage before described is, for distribution network feeder branch road, from root node, terminad advances, then the voltage of node j is:

。

10., based on the power distribution network calculation of short-circuit current system of Distributed Calculation, it is characterized in that comprising:

Data processing module, one is that two are to provide the process to result of calculation, are packaged into the form of standard for being resolved by the carrying out obtaining data from platform;

Topological analysis module is the basis of short circuit calculation, for the physical model of distribution is converted into computation model;

Short circuit calculation module, for realizing the calculating of short-circuit current on distributed platform;

Results analyses module, for carrying out analysis and validation to the result of short circuit calculation, and by the result of calculation of short circuit calculation module and results analyses module or meet at data processing module encapsulation backward on provide or put preservation in storage.

The 11. power distribution network calculation of short-circuit current systems based on Distributed Calculation according to claim 10, it is characterized in that: describedly obtain from platform the XML data that data comprise data that SCADA system provides, CIM, these data are transmitted by the CIS interface of standard.

The 12. power distribution network calculation of short-circuit current systems based on Distributed Calculation according to claim 10, is characterized in that: the mode of described topology is Global Topological or dynamic local topology.