CN101572409B - Self-adaptive device and method for analyzing global power flow of generation, transmission and distribution - Google Patents

Abstract

The invention provides a self-adaptive device and a method for analyzing global power flow of generation, transmission and distribution, which belongs to the technical field of power transmission and distribution of power systems. The device comprises an upper computer and a lower computer, wherein the lower computer comprises a power-flow feedback module; a global power system is divided into a master system and a slave system; according to the characteristics of each part of the power system, the master system adopts a Newton-Raphson power-flow analysis method, while the slave system adopts a forward-backward substitution power-flow analysis method; and the master system and the slave system are organically linked together through master-slave-system association nodes in the middle. The invention has the advantages that the device adopts the power-flow feedback module, uses slave-system node-voltage signals obtained from last power-flow calculation of the slave system for next power-flow calculation of the slave system, and realizes the optimization of the prior forward-backward substitution algorithms, and as the method adopts the forward-backward substitution method which reduces iterations compared with the prior forward-backward substitution method, the device for analyzing global power flow not only guarantees accurate power-flow analysis, but also raises convergence rate and saves storage space.

Description

A kind of self adaptation is sent out transmission ﹠ distribution overall situation current analytical device and method

Technical field

The invention belongs to electric power system power transmission and distribution technical field, particularly a kind of self adaptation is sent out transmission ﹠ distribution overall situation current analytical device and method.

Background technology

Tidal current analysis is a kind of basic electrical analysis computing of research power system mesomeric state ruuning situation, at the electric power system everyway great practical value is arranged, the task of conventional tidal current analysis is to determine the running status of whole system according to given service conditions and network configuration, as the voltage on each bus, distribute power in the network and power loss etc.

Unified tidal current analysis method is generally adopted in the analysis of overall trend in the traditional tidal current analysis of majority, and in fact the characteristics of generating and transmitting system and distribution system are widely different, so-called generating and transmitting system-main system and distribution system-as follows from the difference of system:

(1) difference of network configuration: distribution system is normally radial, and transmission system looped network shape normally; There is a large amount of generator nodes in the generating and transmitting system, it is the PV node, and the majority of the load bus in the distribution system is the PQ node, the linearity of distribution system power flow equation is than the generating and transmitting system height, this has determined to push back for being suitable for the distribution system tidal current analysis before the advantages of simplicity and high efficiency, but the generating and transmitting system trend higher to nonlinearity, preceding pushing back for the difficult constringency performance well that obtains of rule.

(2) difference of network parameter: the short branch road of distribution system is more, and the value conditions of the Jacobian matrix of distribution power flow equation is generally relatively poor; In the overall power flow equation of sending out the transmission ﹠ distribution simultaneous, because branch impedance parameter size differs greatly, the value conditions of its Jacobian matrix is poor, make overall tidal current analysis produce serious numerical problem, therefore, overall trend can't guarantee the reliability of analyzing if adopt unified newton's class algorithm, and the r/x of distribution branch road is generally bigger, and the quick decomposition method of P-Q can't be used for overall tidal current analysis.

(3) difference of trend size: the power unit of sending out the transmission of electricity trend is MW, and distribution trend power unit is generally kW, therefore, sending out the transmission of electricity trend should be not the same with the requirement of distribution trend convergence precision, analyze overall trend if adopt unified algorithm, mark base value and convergence precision, when sending out the convergence of transmission of electricity trend, the error of distribution trend is also very big, is difficult to satisfy simultaneously both requirements.

(4) difference of model: distribution system often is in the three-phase imbalance running status, and generating and transmitting system can be similar to and be thought of as three-phase equilibrium, and in practice, the former adopts the three-phase imbalance model, and the latter adopts single phase model, and therefore overall power flow algorithm must be able to adapt to the difference on this model.

Characteristics difference based on generating and transmitting system and distribution system is big, takes suitable separately power flow algorithm to carry out analytical calculation respectively, for example sends out the quick decomposition method of transmission of electricity trend, and the distribution trend is with before pushing back for method.

In overall tidal current analysis, analyze separately from system load flow if will send out transmission of electricity main system tidal current analysis and distribution, when analyzing when sending out the transmission of electricity trend, distribution system such as is processed at duty value, the load power data are known and given; When analyzing the distribution trend, then generating and transmitting system is considered to equivalent power supply, the voltage data of each distribution root node is known and given, because it is different with the Data Source of distribution trend to send out the transmission of electricity trend, will inevitably on boundary node, produce power mismatch and voltage mismatch, these are referred to as the border amount of mismatch, and this moment, analysis precision and control of quality will be had a strong impact on.Along with the development of tidal current analysis technology, people have begun to consider the connectivity problem between master-slave system to the global system tidal current analysis time, and the concrete analysis thinking is as follows:

Overall electric power system set of node is done following division: the load bus of generating and transmitting system is formed boundary system B, and its set of node is designated as C _B, the node number is N _BThe node of remaining generating and transmitting system has been formed main system M, and its set of node is designated as C _M, the node number is N _MThe node of remaining distribution system is formed from the S of system, and its set of node is designated as C _S, the node number is N _SOverall situation voltage vector

Also be decomposed into the main system voltage vector accordingly

, the boundary voltage vector

With from the system voltage vector Because main system (C _M) and from (the C of system _S) between the branch road that directly do not link, and just indirectly by boundary node (C _B) get in touch, therefore, a kind of natural split form of sending out transmission ﹠ distribution overall situation power flow equation is

$\left\{\begin{array}{c}{\stackrel{\·}{S}}_M\left({\stackrel{\·}{V}}_M\right)-{\stackrel{\·}{S}}_{\mathrm{MM}}\left({\stackrel{\·}{V}}_M\right)-{\stackrel{\·}{S}}_{\mathrm{MB}}({\stackrel{\·}{V}}_M,{\stackrel{\·}{V}}_B)=0\\ {\stackrel{\·}{S}}_B\left({\stackrel{\·}{V}}_B\right)-{\stackrel{\·}{S}}_{\mathrm{BM}}({\stackrel{\·}{V}}_M,{\stackrel{\·}{V}}_B)-{\stackrel{\·}{S}}_{\mathrm{BB}}\left({\stackrel{\·}{V}}_B\right)={\stackrel{\·}{S}}_{\mathrm{BS}}({\stackrel{\·}{V}}_B,{\stackrel{\·}{V}}_S)\end{array}\right.$

${\stackrel{\·}{S}}_S\left({\stackrel{\·}{V}}_S\right)-{\stackrel{\·}{S}}_{\mathrm{SB}}({\stackrel{\·}{V}}_B,{\stackrel{\·}{V}}_S)-{\stackrel{\·}{S}}_{\mathrm{SS}}\left({\stackrel{\·}{V}}_S\right)=0$

In the formula,

With

The rectangular coordinate that is the node injection complex power vector of corresponding node collection is respectively represented;

Be set of node C _XGo up each node and directly flow to set of node C _YThe rectangular coordinate of the vector formed of branch road complex power trend represent; Be set of node C _XThe rectangular coordinate of the vector that the branch road complex power trend that last each node flows directly into set of node self is formed is represented.

This has effectively solved the problem of power and voltage mismatch to a great extent, but because it is not revised in the master-slave system associated nodes is the process of boundary node conversion, this will have a strong impact on the precision of analysis problem, and because it is huge to send out transmission ﹠ distribution overall situation power system network structure, the analytical calculation complexity, amount of calculation is huger when adding the border amount of mismatch that consider to reduce global network, therefore seek a kind ofly calculate accurately, the analytical method of fast convergence rate seems particularly important.

Summary of the invention

Problem at the prior art existence, the invention provides a kind of self adaptation and send out transmission ﹠ distribution overall situation current analytical device and method, by overall electric power system being divided into principal and subordinate two systems, according to electric power system each several part feature, main system adopts the inferior tidal current analysis method of newton-pressgang, before adopting, system pushes back for the tidal current analysis method, master-slave system associated nodes by the centre organically is linked together two parts, when having reduced power and voltage amount of mismatch, in the process of associated nodes conversion, revise, make the contact of master-slave system tightr, reduce distortion, and this device and method also is optimized traditional tidal current analysis method, can reduce iterations like this, save memory space, improve tidal current analysis speed.

Technical scheme of the present invention is achieved in that this device comprises host computer and slave computer, wherein slave computer comprises preposition module, A/D modular converter, DSP module, keyboard and LCD MODULE, trend feedback module, communication module, and preposition module comprises voltage sensor and power sensor.

The voltage sensor that is connected in the preposition module of this device is connected the A/D modular converter respectively with power sensor, the A/D modular converter connects DSP module, communication module and host computer successively, keyboard is connected the DSP module with LCD MODULE, the trend feedback module connects the DSP module.

Described trend feedback module comprises D/A conversion chip, operational amplifier and A/D conversion chip, the output concatenation operation amplifier input terminal of D/A conversion chip wherein, and the output of operational amplifier connects the input of A/D conversion chip.

The current analytical device course of work of the present invention is: the voltage power transducer in the preposition module sends line voltage, the power signal after responding to the A/D modular converter, be converted into the digital signal that the DSP module can be discerned, the DSP module is called the inferior side's tidal current analysis of newton-pressgang module respectively and is before pushed back for the tidal current analysis module input signal of generating and transmitting system and distribution system is carried out tidal current analysis, show the tidal current analysis operation result at last on liquid crystal display screen and host computer, communication module is sent to host computer with identification result.

The trend feedback module acts on the tidal current analysis part of distribution system, and it can feed back the state of the power distribution network of DSP module calculating, realizes pushing back the optimization for the tidal current analysis method before the power distribution network.

Below be the overall process of overall master-slave system tidal current analysis: to overall electric power system, with its set of node carry out following division: interconnected system G, main system M, from the S of system, wherein interconnected system G by the load bus of generating and transmitting system just the root node of distribution system form, its set of node is designated as C _G, the node number is N _GMain system M is made up of the node of remaining generating and transmitting system, and its set of node is designated as C _M, the node number is N _MBe made up of the node of remaining distribution system from the S of system, its set of node is designated as C _S, the node number is N _S, the polar coordinate representation V of overall voltage also is decomposed into the polar coordinate representation V of main system voltage vector accordingly _M, the interconnected system voltage vector polar coordinate representation V _GWith polar coordinate representation V from the system voltage vector _S

Because main system C _MWith from the C of system _SBetween the branch road that directly do not link, and just indirectly by master-slave system associated nodes C _GGet in touch, therefore, a kind of natural split form of sending out transmission ﹠ distribution overall situation power flow equation is

$\left\{\begin{array}{c}{S}_M\left(V_M\right)=S_{\mathrm{MM}}\left(V_M\right)+S_{\mathrm{MG}}(V_M,V_G+\mathrm{\Δ}{V}_G)\\ S_G(V_G+{\mathrm{\ΔV}}_G)=S_{\mathrm{GM}}(V_M,V_G+{\mathrm{\ΔV}}_G)+S_{\mathrm{GG}}(V_G+{\mathrm{\ΔV}}_G)+S_{\mathrm{GS}}(V_G+{\mathrm{\ΔV}}_G,V_S)\end{array}\right.---\left(1\right)$

S _S(V _S)＝S _SG(V _G+ΔV _G，V _S)+S _SS(V _S) (2)

In the formula, S _M, S _GAnd S _SBe respectively the polar coordinate representation of the node injection complex power vector of corresponding node collection: S _XYBe set of node C _XGo up each node and directly flow to set of node C _YThe polar coordinate representation of the vector formed of branch road complex power trend; S _XXBe set of node C _XThe polar coordinate representation of the vector that the branch road complex power trend that last each node flows directly into set of node self is formed.Equation group (1) and (2) be called send out a transmission of electricity power flow equation be main system equation and distribution power flow equation promptly from system equation, S _GSBe the polar coordinate representation of associated variable in the middle of the adaptive iteration, Δ V _GBe the polar coordinate representation of the voltage correction of interconnected system, ${\mathrm{\ΔV}}_G=-\frac{S_G\left(V_G\right)}{{S^{\′}}_G\left(V_G\right)}.$

It is as follows to send out the alternative manner that transmission ﹠ distribution overall situation tidal current analysis calculates at equation (1) and (2) described global system application self-adapting:

Step 1: interconnected system voltage V _GInitialize: V _G ^(k), k=0;

Step 2: with interconnected system voltage V _G ^(k)Be reference voltage, in conjunction with Δ V _G ^(k)Find the solution distribution power flow equation (2), get distribution system voltage vector V _S ^(k+1), and by V _G ^(k), Δ V _G ^(k)And V _S ^(k+1), calculate iteration intermediate variable S _GS ^(k+1)

Step 3: by associated variable S in the middle of the iteration _GS ^(k+1), find the solution and send out transmission of electricity power flow equation (1), get generating and transmitting system voltage vector [V _M ^(k+1)V _G ^(k+1)] ^TAnd Δ V _G ^(K+1)

Step 4: the maximum max| Δ V that judges the mold component of interconnected system voltage difference between adjacent twice iteration _i|, whether less than given convergence index ε, wherein i ∈ C _GIf,, then overall trend iteration convergence; Otherwise k=k+1 changes step 2.

Based on the above-mentioned theory analysis, the present invention adopts self adaptation to send out transmission ﹠ distribution overall situations current analytical device to carry out self adaptation overall situation tidal current analysis and carry out as follows:

Step 1: initialization;

Step 2:, comprise burden with power P, the voltage magnitude V signal of master-slave system PV child node, the burden with power P of master-slave system PQ child node, load or burden without work Q signal with preposition module sensors collection site data-signal;

Step 3: analog signal is converted into digital signal by the A/D modular converter;

Step 4:, analyze at first given master-slave system associated nodes voltage V by self adaptation overall situation tidal current analysis method to the digital signal of DSP module input A/D modular converter output _GInitial value, V _G ^(k), k=0; The DSP module invokes pushes back for the tidal current analysis module before components of system as directed, setting is carried out iterative analysis from the system voltage initial value and to the input data from system, send the trend feedback module with what it analyzed to from the system node information of voltage, and pushed back the node voltage signal that goes out for Algorithm Analysis before the trend feedback module storage last time and use as voltage initial value before pushing back when analyzing for method it to be sent back push back for module before the system next time, realize pushing back optimization before the tradition with this for the tidal current analysis method, this calculates and finishes the result of calculation of back DSP module utilization from system from system, according to formula S _S(V _S)=S _SG(V _G+ Δ V _G, V _S)+S _SS(V _S) try to achieve the iteration intermediate variable S between master-slave system _GS ^(k+1), and the inferior method tidal current analysis of the newton-pressgang module of calling main system part carries out analytical calculation to the digital input signals of main system, by associated variable S in the middle of the iteration _GS ^(k+1), find the solution and send out the transmission of electricity power flow equation $\left\{\begin{array}{c}{S}_M\left(V_M\right)=S_{\mathrm{MM}}\left(V_M\right)+S_{\mathrm{MG}}(V_M,V_G+\mathrm{\Δ}{V}_G)\\ S_G(V_G+{\mathrm{\ΔV}}_G)=S_{\mathrm{GM}}(V_M,V_G+{\mathrm{\ΔV}}_G)+S_{\mathrm{GG}}(V_G+{\mathrm{\ΔV}}_G)+S_{\mathrm{GS}}(V_G+{\mathrm{\ΔV}}_G,V_S)\end{array}\right.$ Try to achieve generating and transmitting system voltage vector [V _M ^(k+1)V _G ^(k+1)] ^TWith Δ V _G ^(k+1), judge the maximum max| Δ V of the mold component of master-slave network interconnected system voltage difference between adjacent twice iteration _i| whether less than given convergence index ε, wherein i ∈ C _GIf,, then overall trend iteration convergence; Otherwise, k=k+1, the analytic process of repeating step four is till satisfying the condition of convergence; Wherein in (1) (2) two formulas, S _M, S _GAnd S _SBe respectively the polar coordinate representation of the node injection complex power vector of corresponding node collection: S _XYBe set of node C _XGo up each node and directly flow to set of node C _YThe polar coordinate representation of the vector formed of branch road complex power trend; S _XXBe set of node C _XThe polar coordinate representation of the vector that the branch road complex power trend that last each node flows directly into set of node self is formed, S _GSBe the polar coordinate representation of associated variable in the middle of the adaptive iteration, V _MBe the polar coordinate representation of main system voltage vector, V _GBe the polar coordinate representation of master-slave network interconnected system voltage vector, V _SBe polar coordinate representation from the system voltage vector, Δ V _GBe the polar coordinate representation of the voltage correction of interconnected system, ${\mathrm{\ΔV}}_G=-\frac{S_G\left(V_G\right)}{{S^{\′}}_G\left(V_G\right)};$

Step 5: analyze and finish, the output analysis result, analysis result comprises the idle Q and the voltage phase angle value of master-slave system PV child node, the voltage magnitude V and the voltage phase angle of PQ child node, the active power P of boundary point, reactive power Q voltage magnitude V and voltage phase angle signal, and the wattful power power of main system balance node injects P, reactive power is injected Q;

Step 6: communication module sends above tidal current analysis result to host computer, and shows this tidal current analysis result on display;

Step 7: end of run.

When carrying out the inferior method tidal current analysis of main system newton-pressgang module analysis, its analytic process step is as follows:

Step 1: operation beginning;

Step 2: the fan-in network initial parameter comprises the node burden with power P and the load or burden without work Q of PQ node, the node burden with power P of PV node and node voltage V, the resistance of each section circuit and reactance value;

Step 3: connecting information forms node admittance matrix Y, element Y wherein between the node that is embodied according to the branch road information table _Ij=G _Ij+ jB _Ij, if i ≠ during j the transadmittance between node i and j, if i=j then is the self-admittance of this node;

Step 4: given node voltage amplitude initial value V _i ⁽⁰⁾, phase angle initial value δ _i ⁽⁰⁾

Step 5: iterations k zero setting, i.e. k=0;

Step 6: with the initial value (V of each node voltage _i ^(k), δ _i ^(k)The substitution formula

${\mathrm{\ΔP}}_i=P_i-V_i\underset{j=1}{\overset{n}{\mathrm{\Σ}}}{V}_j(G_{\mathrm{ij}}\cos {\mathrm{\δ}}_{\mathrm{ij}}+B_{\mathrm{ij}}\sin {\mathrm{\δ}}_{\mathrm{ij}})$

${\mathrm{\ΔQ}}_i=Q_i-V_i\underset{j=1}{\overset{n}{\mathrm{\Σ}}}{V}_j(G_{\mathrm{ij}}\sin {\mathrm{\δ}}_{\mathrm{ij}}-B_{\mathrm{ij}}\cos {\mathrm{\δ}}_{\mathrm{ij}})$

Ask amount of unbalance Δ P _i ^(k), Δ Q _i ^(k);

Step 7: judge whether to satisfy the condition of convergence, promptly $\max \{{\mathrm{\&Delta;P}}_i^{\left(k\right)},{\mathrm{\&Delta;Q}}_i^{\left(k\right)}\}<\mathrm{\&epsiv;}$ , i ∈ C wherein _G, ε is desired tidal current analysis precision, then enters next step as not restraining, and then jumps to step 12 as convergence;

Step 8: bring each node voltage initial value into following formula

$\left[\begin{array}{c}{\mathrm{\&Delta;P}}_1\\ {\mathrm{\&Delta;P}}_2\\ .\\ .\\ .\\ {\mathrm{\&Delta;P}}_{n-1}\\ ...\\ {\mathrm{\&Delta;<figure-callout\; id="1"\; label="Q"\; filenames="H2009100119166E00022.png,H2009100119166E00051.png"\; state="\{\{state\}\}">Q</figure-callout>}}_1\\ {\mathrm{\&Delta;<figure-callout\; id="2"\; label="Q"\; filenames="H2009100119166E00021.png"\; state="\{\{state\}\}">Q</figure-callout>}}_2\\ .\\ .\\ .\\ {\mathrm{\&Delta;Q}}_m\end{array}\right]=\left[\begin{array}{cccccccc}{H}_{11}& {\mathrm{<figure-callout\; id="12"\; label="H"\; filenames="H2009100119166E00021.png,H2009100119166E00022.png"\; state="\{\{state\}\}">H</figure-callout>}}_{12}& ...& {\mathrm{<figure-callout\; id="1"\; label="H"\; filenames="H2009100119166E00022.png,H2009100119166E00051.png"\; state="\{\{state\}\}">H</figure-callout>}}_{1,n-1}& N_{11}& {\mathrm{<figure-callout\; id="12"\; label="N"\; filenames="H2009100119166E00021.png,H2009100119166E00022.png"\; state="\{\{state\}\}">N</figure-callout>}}_{12}& ...& {\mathrm{<figure-callout\; id="1"\; label="N"\; filenames="H2009100119166E00022.png,H2009100119166E00051.png"\; state="\{\{state\}\}">N</figure-callout>}}_{1,m}\\ H_{12}& H_{22}& ...& {\mathrm{<figure-callout\; id="2"\; label="H"\; filenames="H2009100119166E00021.png"\; state="\{\{state\}\}">H</figure-callout>}}_{2,n-1}& N_{21}& N_{22}& ...& {\mathrm{<figure-callout\; id="2"\; label="N"\; filenames="H2009100119166E00021.png"\; state="\{\{state\}\}">N</figure-callout>}}_{2,m}\\ .& .& & .& .& .& & .\\ .& .& ...& .& .& .& ...& .\\ .& .& & .& .& .& & .\\ H_{n-\mathrm{1,1}}& H_{n-\mathrm{1,2}}& ...& H_{n-1,n-1}& N_{n-\mathrm{1,1}}& N_{n-\mathrm{1,2}}& ...& N_{n-1,m}\\ ...& ...& ...& ...& ...& ...& ...& ...\\ M_{11}& {\mathrm{<figure-callout\; id="12"\; label="M"\; filenames="H2009100119166E00021.png,H2009100119166E00022.png"\; state="\{\{state\}\}">M</figure-callout>}}_{12}& ...& {\mathrm{<figure-callout\; id="1"\; label="M"\; filenames="H2009100119166E00022.png,H2009100119166E00051.png"\; state="\{\{state\}\}">M</figure-callout>}}_{1,n-1}& L_{11}& {\mathrm{<figure-callout\; id="12"\; label="L"\; filenames="H2009100119166E00021.png,H2009100119166E00022.png"\; state="\{\{state\}\}">L</figure-callout>}}_{12}& ...& {\mathrm{<figure-callout\; id="1"\; label="L"\; filenames="H2009100119166E00022.png,H2009100119166E00051.png"\; state="\{\{state\}\}">L</figure-callout>}}_{1,m}\\ M_{21}& M_{22}& ...& {\mathrm{<figure-callout\; id="2"\; label="M"\; filenames="H2009100119166E00021.png"\; state="\{\{state\}\}">M</figure-callout>}}_{2,n-1}& L_{21}& L_{22}& ...& {\mathrm{<figure-callout\; id="2"\; label="L"\; filenames="H2009100119166E00021.png"\; state="\{\{state\}\}">L</figure-callout>}}_{2,m}\\ .& .& & .& .& .& & .\\ .& .& ...& .& .& .& ...& .\\ .& .& & .& .& .& & .\\ M_{m,1}& M_{m,2}& ...& M_{m,n-1}& L_{m,1}& L_{m,2}& ...& L_{m,m}\end{array}\right]\left[\begin{array}{c}{\mathrm{\&Delta;\&delta;}}_1\\ {\mathrm{\&Delta;\&delta;}}_2\\ .\\ .\\ .\\ {\mathrm{\&Delta;\&delta;}}_{n-1}\\ .\\ .\\ .\\ {\mathrm{\&Delta;<figure-callout\; id="1"\; label="V"\; filenames="H2009100119166E00022.png,H2009100119166E00051.png"\; state="\{\{state\}\}">V</figure-callout>}}_1/{\mathrm{<figure-callout\; id="1"\; label="V"\; filenames="H2009100119166E00022.png,H2009100119166E00051.png"\; state="\{\{state\}\}">V</figure-callout>}}_1\\ {\mathrm{\&Delta;<figure-callout\; id="2"\; label="V"\; filenames="H2009100119166E00021.png"\; state="\{\{state\}\}">V</figure-callout>}}_2/{\mathrm{<figure-callout\; id="2"\; label="V"\; filenames="H2009100119166E00021.png"\; state="\{\{state\}\}">V</figure-callout>}}_2\\ .\\ .\\ .\\ {\mathrm{\&Delta;V}}_m/V_m\end{array}\right]$ Ask correction

Each element of equational coefficient matrix-Jacobi matrix.

Wherein when i ≠ j, have

$H_{\mathrm{ij}}=\frac{{\&PartialD;P}_i}{{\&PartialD;\mathrm{\&delta;}}_j}=V_i{V}_j(G_{\mathrm{ij}}\sin {\mathrm{\&delta;}}_{\mathrm{ij}}-B_{\mathrm{ij}}\cos {\mathrm{\&delta;}}_{\mathrm{ij}}),N_{\mathrm{ij}}=\frac{{\&PartialD;P}_i}{{\&PartialD;V}_j}{V}_j=V_i{V}_j(G_{\mathrm{ij}}\cos {\mathrm{\&delta;}}_{\mathrm{ij}}+B_{\mathrm{ij}}\sin {\mathrm{\&delta;}}_{\mathrm{ij}})$

$M_{\mathrm{ij}}=\frac{{\&PartialD;Q}_i}{{\&PartialD;\mathrm{\&delta;}}_j}={-V}_i{V}_j(G_{\mathrm{ij}}\cos {\mathrm{\&delta;}}_{\mathrm{ij}}+B_{\mathrm{ij}}\sin {\mathrm{\&delta;}}_{\mathrm{ij}}),L_{\mathrm{ij}}=\frac{{\&PartialD;Q}_i}{{\&PartialD;V}_j}{V}_j=V_i{V}_j(G_{\mathrm{ij}}\sin {\mathrm{\&delta;}}_{\mathrm{ij}}-B_{\mathrm{ij}}\cos {\mathrm{\&delta;}}_{\mathrm{ij}})$

When i=j, have

$H_{\mathrm{ii}}=\frac{{\&PartialD;P}_i}{{\&PartialD;\mathrm{\&delta;}}_i}=-V_i\underset{j=1,i\&NotEqual;j}{\overset{n}{\mathrm{\&Sigma;}}}{V}_j(G_{\mathrm{ij}}\sin {\mathrm{\&delta;}}_{\mathrm{ij}}-B_{\mathrm{ij}}\cos {\mathrm{\&delta;}}_{\mathrm{ij}})=-V_i^2{B}_{\mathrm{ii}}-Q_i,$

$N_{\mathrm{ii}}=\frac{{\&PartialD;P}_i}{{\&PartialD;V}_i}{V}_i=V_i\underset{j=1,j\&NotEqual;i}{\overset{n}{\mathrm{\&Sigma;}}}{V}_j(G_{\mathrm{ij}}\cos {\mathrm{\&delta;}}_{\mathrm{ij}}+B_{\mathrm{ij}}\sin {\mathrm{\&delta;}}_{\mathrm{ij}})+2V_i^2{G}_{\mathrm{ii}}=V_i^2{G}_{\mathrm{ii}}+P_i$

$M_{\mathrm{ii}}=\frac{{\&PartialD;Q}_i}{{\&PartialD;\mathrm{\&delta;}}_i}=V_i\underset{j=1,j\&NotEqual;i}{\overset{n}{\mathrm{\&Sigma;}}}{V}_j(G_{\mathrm{ij}}\cos {\mathrm{\&delta;}}_{\mathrm{ij}}+B_{\mathrm{ij}}\sin {\mathrm{\&delta;}}_{\mathrm{ij}})=-V_i^2{G}_{\mathrm{ii}}+P_i$

$L_{\mathrm{ii}}=\frac{{\&PartialD;Q}_i}{{\&PartialD;V}_i}{V}_i=V_i\underset{j=1,j\&NotEqual;i}{\overset{n}{\mathrm{\&Sigma;}}}{V}_j(G_{\mathrm{ij}}\sin {\mathrm{\&delta;}}_{\mathrm{ij}}\cos {\mathrm{\&delta;}}_{\mathrm{ij}}+B_{\mathrm{ij}}\cos {\mathrm{\&delta;}}_{\mathrm{ij}})-2V_i^2{B}_{\mathrm{ii}}=-V_i^2{B}_{\mathrm{ii}}+Q_i$

Step 9: separate the update equation formula, ask the variable of each node voltage, be i.e. correction value Δ V _i ^(k), Δ δ _i ^(k)

Step 10: according to formula $V_i^{(k+1)}=V_i^{\left(k\right)}+{\mathrm{\&Delta;V}}_i^{\left(k\right)},$ δ _j ^(k+1)=δ _i ^(k)+ Δ δ _i ^(k)Calculate the new value of each node voltage, promptly revised value;

Step 11:k=k+1 forwards step 6 to;

Step 12: calculate branch power and distribute PV node reactive power and balance node injecting power;

Step 13: output.

Push back before system for method tidal current analysis module carrying out, the analytic process step is as follows:

Step 1: operation beginning;

Step 2: the fan-in network initial parameter, comprise the node burden with power P and the load or burden without work Q of the PQ node except that root node,

The resistance of each section circuit and reactance value;

Step 3: given root node voltage initial value V ₀

Step 4: iterations k zero setting, i.e. k=0;

Step 5: each node voltage initialize, establishing trend calculating iterations is k=0, the trend calculation times is m, when m=1, $V_i^{\left(m\right)}=V_0^{\left(m\right)},$ Promptly root node voltage is composed to its corresponding child node; When m＞1, $V_i^{\left(m\right)}=V_i^{(m-1)}$ The i.e. node voltage initial value that node voltage in the last time calculation of tidal current is calculated as this trend;

Step 6: according to formula ${\stackrel{\&CenterDot;}{I}}_{L,i}=\frac{P_{L,i}-j{Q}_{L,i}}{{\stackrel{\&CenterDot;}{V}}_i^{*}}$ Calculate each child node electric current, in the formula Be the electric current of node i,

Be the node voltage of node i, P _{L, i}Be the active power of node i, Q _{L, i}Reactive power for node i;

Step 7: according to connecting information between the node that from the system branch information table, embodies by frontier node to the root node recursion, according to formula ${\stackrel{\&CenterDot;}{I}}_t={\stackrel{\&CenterDot;}{I}}_{L,j}+\underset{k\&Element;d}{\mathrm{\&Sigma;}}{\stackrel{\&CenterDot;}{I}}_k$ Calculate each branch current, wherein,

Be branch road b _tBranch current, j is branch road b _tFrontier node,

Be branch road b _tThe frontier node electric current, set d is for being the set of fingers of father node with node j;

Step 8: by root node to the frontier node recursion, according to formula ${\stackrel{\&CenterDot;}{V}}_j={\stackrel{\&CenterDot;}{V}}_i-{\stackrel{\&CenterDot;}{I}}_t(R_t+j{X}_t)$ Calculate each node voltage, wherein, i, j are respectively b _tThe first and last node R of branch road _j, X _jBe respectively b _tThe resistance of branch road and reactance;

Step 9: the correction amount V that calculates each load bus voltage _{J (k)}=| V _{J (k)}-V _{J (k-1)}|, judge whether to satisfy the condition of convergence, i.e. max Δ V _{J (k)}＜ε, wherein i ∈ C _G, ε is desired tidal current analysis precision, as not restraining, then with the new value of each node voltage as initial value, k=k+1 restarts to carry out next iteration from the 6th step, otherwise changes next step over to;

Step 10: according to formula $S_{\mathrm{ij}}={\stackrel{\&CenterDot;}{V}}_i{\stackrel{\&CenterDot;}{I}}_{\mathrm{ij}},$ ${\mathrm{LS}}_{\mathrm{ij}}={\left|{\stackrel{\&CenterDot;}{I}}_{\mathrm{ij}}\right|}^2{Z}_{\mathrm{ij}}$ Power flows and line loss on the computational scheme;

Step 11: output.

This current analytical device is assemblied in the object subsystem local control box, and it is to carry out tidal current analysis according to the real time node voltage of the electric power system of preposition module collection and power signal.

Advantage of the present invention: traditional overall electric power system tide analysis generally combines in the master-slave system with electric power networks and the associated nodes boundary node of master-slave system is not revised when analyzing, cause tidal current analysis inaccurate, and when the incidence relation of having considered between master-slave system, algorithm is not optimized, cause that tidal current analysis speed is slow when analyzing huge electric power networks, iterations is many.And this self adaptation overall situation current analytical device has been considered in the connecting each other between master-slave system associated nodes to be revised when overall tidal current analysis is carried out in overall electric power system, taken into full account being closely connected between master-slave system, analyze accurately, and should adopt the trend feedback module in the overall situation current analytical device, it can be used for the node voltage signal from system from system's gained last once trend is calculated next time the tidal current analysis from system, this has just realized pushing back before traditional for the class optimization Algorithm, the method has reduced iterations than pushing back for method before traditional, therefore this overall current analytical device not only guarantees the accurate of tidal current analysis, and improved rate of convergence, save memory space.

Description of drawings

Fig. 1 sends out transmission ﹠ distribution overall situation current analytical device general structure block diagram for a kind of self adaptation of the present invention;

Fig. 2 is the interface electrical schematic diagram of A/D modular converter and TMS320C5416 in the embodiment of the invention device;

Fig. 3 is an embodiment of the invention communication module electrical schematic diagram;

Fig. 4 is embodiment of the invention keyboard and liquid crystal display electrical schematic diagram;

Fig. 5 is an embodiment of the invention trend feedback module electrical schematic diagram;

Fig. 6 pushes back for method tidal current analysis flow chart for before the embodiment of the invention;

Fig. 7 is the inferior method tidal current analysis of an embodiment of the invention newton-pressgang flow chart;

Fig. 8 is an embodiment of the invention analytical method flow chart.

Embodiment

A kind of self adaptation of the present invention is sent out transmission ﹠ distribution overall situation current analytical device detailed structure and is illustrated in conjunction with the embodiments.

Below be the circuit structure and the course of work: hardware unit of the present invention comprises host computer and slave computer, slave computer comprises preposition module, the A/D modular converter, the DSP module, keyboard and LCD MODULE, the trend feedback module, communication module, preposition module comprises voltage sensor and power sensor, the voltage sensor that is connected in the preposition module of this device is connected the A/D modular converter respectively with power sensor, the A/D modular converter connects the DSP module successively, communication module and host computer, keyboard is connected the DSP module with LCD MODULE, the trend feedback module connects DSP module (as shown in Figure 1).

The real-time voltage of voltage sensor in the preposition module and power sensor acquisition system and power signal, after filtration behind the wave circuit, being input to the A/D modular converter samples, voltage sensor is selected PT204A for use, power sensor is selected 8481A for use, and what the A/D modular converter adopted is the ADS7805 analog to digital converter.

The digital signal of A/D modular converter output flows to the DSP module and carries out calculation process, and that the DSP module adopts is TMS320C5416, the output pin D of A/D modular converter ₀-D ₁₅Respectively with the input pin D of DSP module ₀-D ₁₅Link to each other (as shown in Figure 2).

The trend feedback module adopts DAC0832, the electric current Series Negative Feedback Amplifier, ADC0801 realizes, the data output end HD0 of DSP module, HD1, HD2, HD3, HD4, HD5, HD6, HD7 respectively with the data input pin DI0 of trend feedback module DAC0832, DI1, DI2, DI3, DI4, DI5, DI6, DI7 links to each other, carry out transformation of scale through the analog voltage signal after the trend feedback module DAC0832 conversion through the electric current Series Negative Feedback Amplifier, the output signal of amplifier offers A/D modular converter ADC0801 and carries out analog-to-digital conversion, the data output end D0 of A/D modular converter ADC0801, D1, D2, D3, D4, D5, D6, D7 respectively with the data input pin D0 of DSP module, D1, D2, D3, D4, D5, D6, D7 links to each other, (as shown in Figure 5).

Communication module has adopted the MAX232 chip, and this chip is RS-232C and the Transistor-Transistor Logic level conversion chip of using always, and is easy to use.The serial communication of host computer and DSP intermodule adopts the RS232 agreement.The pin 11,12 of communication module MAX232 is connected with the pin 45,59 of DSP module respectively, (as shown in Figure 3), the RS-232 socket links to each other with the host computer serial ports by cable, communication module MAX232 passes to the DSP module after the level signal of the RS-232C standard that receives being converted to the data of serial, produce to receive and interrupt, for the DSP resume module; Simultaneously can be according to the order that receives, the level signal that the serial data that sends is converted to the RS-232C standard is issued host computer;

Keyboard and LCD MODULE are used for the human-computer dialogue operation, and as device is carried out parameter setting, running state monitoring etc., the content of this part is finished by button is set in the hardware designs.The interface circuit of keyboard and LCD MODULE (as shown in Figure 4).

This embodiment of transmission ﹠ distribution overall situation electric power system tide analytical equipment that sends out that the present invention adopts carries out tidal current analysis:

Sending out electric power transmission network with the somewhere electric power system is example, and the connection of electric power transmission network branch road is sent out in this area's electric power system and impedance information sees Table 1-1, and rated voltage is 230kv.

Choose " broad sense load " node of the generating and transmitting system of wherein electric power system and form boundary system B, its set of node is designated as C _B, the node number is N _BThe node of remaining generating and transmitting system has been formed main system M, and its set of node is designated as C _M, the node number is N _MThe node of remaining distribution system is formed from the S of system, and its set of node is designated as C _S, the node number is N _S

Adopt self adaptation to send out transmission ﹠ distribution overall situation current analytical device and carry out self adaptation overall situation tidal current analysis method, carry out as follows: (as shown in Figure 6)

Step 1: initialization;

Step 2: with preposition module sensors collection site data-signal, comprise burden with power P, load or burden without work Q signal, as showing 2-1 to showing shown in the 2-18 from the PQ of system node;

Step 3: analog signal is converted into digital signal by the A/D modular converter;

Step 4:, analyze at first given master-slave system associated nodes voltage V by self adaptation overall situation tidal current analysis method to the digital signal of DSP module input A/D modular converter output _GInitial value, general assignment is a load voltage value; The DSP module invokes pushes back for the tidal current analysis module before components of system as directed, setting is carried out iterative analysis from the system voltage initial value and to the input data from system, send the trend feedback module with what it analyzed to from the system node information of voltage, and pushed back the node voltage signal that goes out for Algorithm Analysis before the trend feedback module storage last time and use as voltage initial value before pushing back when analyzing for method it to be sent back push back for module before the system next time, realize pushing back optimization before the tradition with this for the tidal current analysis method, this calculates and finishes the result of calculation of back DSP module utilization from system from system, according to formula S _S(V _S)=S _SG(V _G+ Δ V _G, V _S)+S _SS(V _S) (1) try to achieve the iteration intermediate variable S between master-slave system _GS ^(k+1), and the inferior method tidal current analysis of the newton-pressgang module of calling main system part carries out analytical calculation to the digital input signals of main system, by associated variable S in the middle of the iteration _GS ^(k+1), find the solution and send out the transmission of electricity power flow equation

$\left\{\begin{array}{c}{S}_M\left(V_M\right)=S_{\mathrm{MM}}\left(V_M\right)+S_{\mathrm{MG}}(V_M,V_G+\mathrm{\&Delta;}{V}_G)\\ S_G(V_G+{\mathrm{\&Delta;V}}_G)=S_{\mathrm{GM}}(V_M,V_G+{\mathrm{\&Delta;V}}_G)+S_{\mathrm{GG}}(V_G+{\mathrm{\&Delta;V}}_G)+S_{\mathrm{GS}}(V_G+{\mathrm{\&Delta;V}}_G,V_S)\end{array}\right.---\left(2\right)$

Try to achieve generating and transmitting system voltage vector [V _M ^(k+1)V _G ^(k+1)] ^TWith Δ V _G ^(k+1), judge the maximum max| Δ V of the mold component of master-slave network interconnected system voltage difference between adjacent twice iteration _i| whether less than given convergence index ε, wherein i ∈ C _G, ε=0.001, if, then overall trend iteration convergence; Otherwise, k=k+1, the analytic process of repeating step four is till satisfying the condition of convergence; Wherein in (1) (2) two formulas, S _M, S _GAnd S _SBe respectively the polar coordinate representation of the node injection complex power vector of corresponding node collection: S _XYBe set of node C _XGo up each node and directly flow to set of node C _YThe polar coordinate representation of the vector formed of branch road complex power trend; S _XXBe set of node C _XThe polar coordinate representation of the vector that the branch road complex power trend that last each node flows directly into set of node self is formed, S _GSBe the polar coordinate representation of associated variable in the middle of the adaptive iteration, V _MBe the polar coordinate representation of main system voltage vector, V _GBe the polar coordinate representation of master-slave network interconnected system voltage vector, V _SBe polar coordinate representation from the system voltage vector, Δ V _GBe the polar coordinate representation of the voltage correction of interconnected system, ${\mathrm{\&Delta;V}}_G=-\frac{S_G\left(V_G\right)}{{S^{\&prime;}}_G\left(V_G\right)};$

Step 5: analyze and finish, the output analysis result, wherein analysis result comprises the idle Q and the voltage phase angle value of master-slave system PV child node, the voltage magnitude V and the voltage phase angle of PQ child node, the active power P of boundary point, reactive power Q voltage magnitude V and voltage phase angle signal, and the wattful power power of main system balance node injects P, reactive power is injected Q;

Step 6: communication module sends above analysis result to host computer, and shows on display;

Step 7: end of run.

Carrying out the inferior method tidal current analysis of main system newton-pressgang module, the analytic process step is as follows: (as shown in Figure 7)

Step 1: operation beginning;

Step 2: the fan-in network initial parameter comprises the node burden with power P and the load or burden without work Q of PQ node, the node burden with power P of PV node and node voltage V, the resistance of each section circuit and reactance value;

Step 3: connecting information forms node admittance matrix Y, element Y wherein between the node that is embodied according to the branch road information table _Ij=G _Ij+ jB _Ij, if i ≠ during j the transadmittance between node i and j, if i=j then is the self-admittance of this node;

Step 4: given node voltage amplitude initial value V _i ⁽⁰⁾, phase angle initial value δ _i ⁽⁰⁾

Step 5: iterations k zero setting, i.e. k=0;

Step 6: with the initial value (V of each node voltage _i ^(k), δ _i ^(k)) the substitution formula

${\mathrm{\&Delta;P}}_i=P_i-V_i\underset{j=1}{\overset{n}{\mathrm{\&Sigma;}}}{V}_j(G_{\mathrm{ij}}\cos {\mathrm{\&delta;}}_{\mathrm{ij}}+B_{\mathrm{ij}}\sin {\mathrm{\&delta;}}_{\mathrm{ij}})$

${\mathrm{\&Delta;Q}}_i=Q_i-V_i\underset{j=1}{\overset{n}{\mathrm{\&Sigma;}}}{V}_j(G_{\mathrm{ij}}\sin {\mathrm{\&delta;}}_{\mathrm{ij}}-B_{\mathrm{ij}}\cos {\mathrm{\&delta;}}_{\mathrm{ij}}),$

Ask amount of unbalance Δ P _i ^(k), Δ Q _i ^(k);

Step 7: judge whether to satisfy the condition of convergence, promptly $\max \{{\mathrm{\&Delta;P}}_i^{\left(k\right)},{\mathrm{\&Delta;Q}}_i^{\left(k\right)}\}<\mathrm{\&epsiv;},$ I ∈ C wherein _G, ε is desired tidal current analysis precision, then enters next step as not restraining, and then jumps to step 12 as convergence;

Step 8: bring each node voltage initial value into following formula

$\left[\begin{array}{c}{\mathrm{\&Delta;P}}_1\\ {\mathrm{\&Delta;P}}_2\\ .\\ .\\ .\\ {\mathrm{\&Delta;P}}_{n-1}\\ ...\\ {\mathrm{\&Delta;<figure-callout\; id="1"\; label="Q"\; filenames="H2009100119166E00022.png,H2009100119166E00051.png"\; state="\{\{state\}\}">Q</figure-callout>}}_1\\ {\mathrm{\&Delta;<figure-callout\; id="2"\; label="Q"\; filenames="H2009100119166E00021.png"\; state="\{\{state\}\}">Q</figure-callout>}}_2\\ .\\ .\\ .\\ {\mathrm{\&Delta;Q}}_m\end{array}\right]=\left[\begin{array}{cccccccc}{H}_{11}& {\mathrm{<figure-callout\; id="12"\; label="H"\; filenames="H2009100119166E00021.png,H2009100119166E00022.png"\; state="\{\{state\}\}">H</figure-callout>}}_{12}& ...& {\mathrm{<figure-callout\; id="1"\; label="H"\; filenames="H2009100119166E00022.png,H2009100119166E00051.png"\; state="\{\{state\}\}">H</figure-callout>}}_{1,n-1}& N_{11}& {\mathrm{<figure-callout\; id="12"\; label="N"\; filenames="H2009100119166E00021.png,H2009100119166E00022.png"\; state="\{\{state\}\}">N</figure-callout>}}_{12}& ...& {\mathrm{<figure-callout\; id="1"\; label="N"\; filenames="H2009100119166E00022.png,H2009100119166E00051.png"\; state="\{\{state\}\}">N</figure-callout>}}_{1,m}\\ H_{12}& H_{22}& ...& {\mathrm{<figure-callout\; id="2"\; label="H"\; filenames="H2009100119166E00021.png"\; state="\{\{state\}\}">H</figure-callout>}}_{2,n-1}& N_{21}& N_{22}& ...& {\mathrm{<figure-callout\; id="2"\; label="N"\; filenames="H2009100119166E00021.png"\; state="\{\{state\}\}">N</figure-callout>}}_{2,m}\\ .& .& & .& .& .& & .\\ .& .& ...& .& .& .& ...& .\\ .& .& & .& .& .& & .\\ H_{n-\mathrm{1,1}}& H_{n-\mathrm{1,2}}& ...& H_{n-1,n-1}& N_{n-\mathrm{1,1}}& N_{n-\mathrm{1,2}}& ...& N_{n-1,m}\\ ...& ...& ...& ...& ...& ...& ...& ...\\ M_{11}& {\mathrm{<figure-callout\; id="12"\; label="M"\; filenames="H2009100119166E00021.png,H2009100119166E00022.png"\; state="\{\{state\}\}">M</figure-callout>}}_{12}& ...& {\mathrm{<figure-callout\; id="1"\; label="M"\; filenames="H2009100119166E00022.png,H2009100119166E00051.png"\; state="\{\{state\}\}">M</figure-callout>}}_{1,n-1}& L_{11}& {\mathrm{<figure-callout\; id="12"\; label="L"\; filenames="H2009100119166E00021.png,H2009100119166E00022.png"\; state="\{\{state\}\}">L</figure-callout>}}_{12}& ...& {\mathrm{<figure-callout\; id="1"\; label="L"\; filenames="H2009100119166E00022.png,H2009100119166E00051.png"\; state="\{\{state\}\}">L</figure-callout>}}_{1,m}\\ M_{21}& M_{22}& ...& {\mathrm{<figure-callout\; id="2"\; label="M"\; filenames="H2009100119166E00021.png"\; state="\{\{state\}\}">M</figure-callout>}}_{2,n-1}& L_{21}& L_{22}& ...& {\mathrm{<figure-callout\; id="2"\; label="L"\; filenames="H2009100119166E00021.png"\; state="\{\{state\}\}">L</figure-callout>}}_{2,m}\\ .& .& & .& .& .& & .\\ .& .& ...& .& .& .& ...& .\\ .& .& & .& .& .& & .\\ M_{m,1}& M_{m,2}& ...& M_{m,n-1}& L_{m,1}& L_{m,2}& ...& L_{m,m}\end{array}\right]\left[\begin{array}{c}{\mathrm{\&Delta;\&delta;}}_1\\ {\mathrm{\&Delta;\&delta;}}_2\\ .\\ .\\ .\\ {\mathrm{\&Delta;\&delta;}}_{n-1}\\ .\\ .\\ .\\ {\mathrm{\&Delta;<figure-callout\; id="1"\; label="V"\; filenames="H2009100119166E00022.png,H2009100119166E00051.png"\; state="\{\{state\}\}">V</figure-callout>}}_1/{\mathrm{<figure-callout\; id="1"\; label="V"\; filenames="H2009100119166E00022.png,H2009100119166E00051.png"\; state="\{\{state\}\}">V</figure-callout>}}_1\\ {\mathrm{\&Delta;<figure-callout\; id="2"\; label="V"\; filenames="H2009100119166E00021.png"\; state="\{\{state\}\}">V</figure-callout>}}_2/{\mathrm{<figure-callout\; id="2"\; label="V"\; filenames="H2009100119166E00021.png"\; state="\{\{state\}\}">V</figure-callout>}}_2\\ .\\ .\\ .\\ {\mathrm{\&Delta;V}}_m/V_m\end{array}\right]$ Ask correction

Each element of equational coefficient matrix-Jacobi matrix.

Wherein when i ≠ j, have

$H_{\mathrm{ij}}=\frac{{\&PartialD;P}_i}{{\&PartialD;\mathrm{\&delta;}}_j}=V_i{V}_j(G_{\mathrm{ij}}\sin {\mathrm{\&delta;}}_{\mathrm{ij}}-B_{\mathrm{ij}}\cos {\mathrm{\&delta;}}_{\mathrm{ij}}),N_{\mathrm{ij}}=\frac{{\&PartialD;P}_i}{{\&PartialD;V}_j}{V}_j=V_i{V}_j(G_{\mathrm{ij}}\cos {\mathrm{\&delta;}}_{\mathrm{ij}}+B_{\mathrm{ij}}\sin {\mathrm{\&delta;}}_{\mathrm{ij}})$

$M_{\mathrm{ij}}=\frac{{\&PartialD;Q}_i}{{\&PartialD;\mathrm{\&delta;}}_j}={-V}_i{V}_j(G_{\mathrm{ij}}\cos {\mathrm{\&delta;}}_{\mathrm{ij}}+B_{\mathrm{ij}}\sin {\mathrm{\&delta;}}_{\mathrm{ij}}),L_{\mathrm{ij}}=\frac{{\&PartialD;Q}_i}{{\&PartialD;V}_j}{V}_j=V_i{V}_j(G_{\mathrm{ij}}\sin {\mathrm{\&delta;}}_{\mathrm{ij}}-B_{\mathrm{ij}}\cos {\mathrm{\&delta;}}_{\mathrm{ij}})$

When i=j, have

$H_{\mathrm{ii}}=\frac{{\&PartialD;P}_i}{{\&PartialD;\mathrm{\&delta;}}_i}=-V_i\underset{j=1,i\&NotEqual;j}{\overset{n}{\mathrm{\&Sigma;}}}{V}_j(G_{\mathrm{ij}}\sin {\mathrm{\&delta;}}_{\mathrm{ij}}-B_{\mathrm{ij}}\cos {\mathrm{\&delta;}}_{\mathrm{ij}})=-V_i^2{B}_{\mathrm{ii}}-Q_i,$

$N_{\mathrm{ii}}=\frac{{\&PartialD;P}_i}{{\&PartialD;V}_i}{V}_i=V_i\underset{j=1,j\&NotEqual;i}{\overset{n}{\mathrm{\&Sigma;}}}{V}_j(G_{\mathrm{ij}}\cos {\mathrm{\&delta;}}_{\mathrm{ij}}+B_{\mathrm{ij}}\sin {\mathrm{\&delta;}}_{\mathrm{ij}})+2V_i^2{G}_{\mathrm{ii}}=V_i^2{G}_{\mathrm{ii}}+P_i$

$M_{\mathrm{ii}}=\frac{{\&PartialD;Q}_i}{{\&PartialD;\mathrm{\&delta;}}_i}=V_i\underset{j=1,j\&NotEqual;i}{\overset{n}{\mathrm{\&Sigma;}}}{V}_j(G_{\mathrm{ij}}\cos {\mathrm{\&delta;}}_{\mathrm{ij}}+B_{\mathrm{ij}}\sin {\mathrm{\&delta;}}_{\mathrm{ij}})=-V_i^2{G}_{\mathrm{ii}}+P_i$

$L_{\mathrm{ii}}=\frac{{\&PartialD;Q}_i}{{\&PartialD;V}_i}{V}_i=V_i\underset{j=1,j\&NotEqual;i}{\overset{n}{\mathrm{\&Sigma;}}}{V}_j(G_{\mathrm{ij}}\sin {\mathrm{\&delta;}}_{\mathrm{ij}}\cos {\mathrm{\&delta;}}_{\mathrm{ij}}+B_{\mathrm{ij}}\cos {\mathrm{\&delta;}}_{\mathrm{ij}})-2V_i^2{B}_{\mathrm{ii}}=-V_i^2{B}_{\mathrm{ii}}+Q_i$

Step 9: separate the update equation formula, ask the variable of each node voltage, be i.e. correction value Δ V _i ^(k), Δ δ _i ^(k)

Step 10: according to formula $V_i^{(k+1)}=V_i^{\left(k\right)}+{\mathrm{\&Delta;V}}_i^{\left(k\right)},$ δ _i ^(k+1)=δ _i ^(k)+ Δ δ _i ^(k)Calculate the new value of each node voltage, promptly revised value;

Step 11:k=k+1 forwards step 6 to;

Step 12: calculate branch power and distribute PV node reactive power and balance node injecting power;

Step 13: output.

Push back before system for method tidal current analysis module carrying out, the analytic process step is as follows: (as shown in Figure 8)

Step 1: operation beginning;

Step 2: the fan-in network initial parameter comprises the node burden with power P and the load or burden without work Q of each the PQ node except that root node, the resistance of each section circuit and reactance value;

Step 3: given root node voltage initial value V ₀

Step 4: iterations k zero setting, i.e. k=0;

Step 5: each node voltage initialize, establishing trend calculating iterations is k=0, the trend calculation times is m, when m=1, $V_i^{\left(m\right)}=V_0^{\left(m\right)},$ Promptly root node voltage is composed to its corresponding child node; When m＞1, $V_i^{\left(m\right)}=V_i^{(m-1)}$ The i.e. node voltage initial value that node voltage in the last time calculation of tidal current is calculated as this trend;

Step 6: according to formula ${\stackrel{\&CenterDot;}{I}}_{L,i}=\frac{P_{L,i}-j{Q}_{L,i}}{{\stackrel{\&CenterDot;}{V}}_i^{*}}$ Calculate each child node electric current, in the formula

Be the electric current of node i,

Be the node voltage of node i, P _{L, i}Be the active power of node i, Q _{L, i}Reactive power for node i;

Step 7: according to connecting information between the node that from the system branch information table, embodies by frontier node to the root node recursion, according to formula ${\stackrel{\&CenterDot;}{I}}_t={\stackrel{\&CenterDot;}{I}}_{L,j}+\underset{k\&Element;d}{\mathrm{\&Sigma;}}{\stackrel{\&CenterDot;}{I}}_k$ Calculate each branch current, wherein, Be branch road b _tBranch current, j is branch road b _tFrontier node,

Be branch road b _tThe frontier node electric current, set d is for being the set of fingers of father node with node j, branch of a network information is shown in table 2-2;

Step 8: by root node to the frontier node recursion, according to formula ${\stackrel{\&CenterDot;}{V}}_j={\stackrel{\&CenterDot;}{V}}_i-{\stackrel{\&CenterDot;}{I}}_t(R_t+j{X}_t)$ Calculate each node voltage, wherein, i, j are respectively b _tThe first and last node R of branch road _j, X _jBe respectively b _tThe resistance of branch road and reactance;

Step 9: the correction amount V that calculates each load bus voltage _{J (k)}=| V _{J (k)}-V _{J (k-1)}|, judge whether to satisfy the condition of convergence, i.e. max Δ V _{J (k)}＜ε, wherein i ∈ C _G, ε is desired tidal current analysis precision, as not restraining, as initial value, then k=k+1 restarts to carry out next iteration from the 6th step, otherwise changes next step over to the new value of each node voltage;

Step 10: according to formula $S_{\mathrm{ij}}={\stackrel{\&CenterDot;}{V}}_i{\stackrel{\&CenterDot;}{I}}_{\mathrm{ij}},$ ${\mathrm{LS}}_{\mathrm{ij}}={\left|{\stackrel{\&CenterDot;}{I}}_{\mathrm{ij}}\right|}^2{Z}_{\mathrm{ij}}$ Power flows and line loss on the computational scheme;

Step 11: output.

The node power information of the radial network of distribution network 9 strips of this area's electric power system and branch road connect and impedance information sees Table 2-1 to showing 2-18, rated voltage is 66kv, wherein the meritorious and load or burden without work of child node is recorded by power sensor, and the resistance reactance of branch road and the information of line length are the server typings.

Table 1-1 sends out the electric power transmission network branch road and connects and impedance information

Numbering	First node	End-node	Resistance (Ω)	Reactance (Ω)
					1	14	17	15.4696	76.76789
2	12	14	1.0957254	3.37848665
					3	8	12	0.815126	3.039582
4	2	12	19.4022	72.03965
					5	1	20	4.65952	22.773404
6	1	27	11.91976	58.257827
					7	8	10	5.9768	29.21161
8	2	7	11.36688	56.408142
					9	7	27	9.3688	46.49267
10	2	29	15.429204	57.288063
					11	20	22	7.64272	37.353794
12	24	26	9.43896	46.132917
					13	12	13	9.704	47.4283
14	4	6	0.624	3.0498
					15	12	24	7.09008	34.652766
16	1	8	1.86287	6.94659
					17	2	27	5.414	26.460925
18	4	7	4.23864	20.716353
					19	14	18	0.54784	2.04288
20	18	19	0.333091	1.242087
					21	14	19	0.33	1.24
22	1	16	0.333091	1.242087

23	1	7	0.333091	1.242087
					24	24	25	1.29142628	57.31366667
25	27	28	1.29142628	57.31366667
					26	14	15	1.72	49.81
27	2	3	0.15797	30.55248
					28	20	21	1.64022222	69.212
29	4	5	1.55358	49.168255
					30	8	9	0.613355	26.115835
31	10	11	1.68145185	65.17866667
					32	22	23	1.75679452	72.65377778
33	2	30	15.429204	57.288063
					34	1	31	0.333091	1.242087

No. 1 radial distribution sub-network node power information of table 2-1

Sequence number	P(kw)	Q(kvar)
			0	0	0
1	28.25367	5.733141
			2	14.12684	2.872457
3	8.89990	1.81294
			4	0	0
5	6.67493	1.35382
			6	4.44995	0.90642
7	0	0
			8	44.49953	9.041175

9	0	0
			10	7.063418	1.436228
11	0	0
			12	5.650734	1.141919
13	15.89269	3.225627
			14	6.67493	1.353822
15	0	0
			16	0	0
17	0	0
			18	0	0
19	0.706342	0.141268
			20	0	0
21	4.449953	0.906472
			22	0	0
23	28.25367	5.733141
			24	28.25367	5.733141
25	0	0
			26	0	0
27	0	0
			28	6.67493	1.353822
29	8.899906	1.812944
			30	6.67493	1.353822

No. 1 radial distribution sub-network branch road of table 2-2 connects and impedance information

Sequence number	First node	End-node	r(Ω)	x(Ω)	Length (km)
						1	0	1	0.7056	1.41792	3.36
2	0	2	3.942	6.2634	14.6
						3	2	3	8.775	9.087	19.5
4	2	4	8.127	8.41596	18.06
						5	4	5	3.51	3.6348	7.8
6	4	6	0.045	0.0466	0.1
						7	0	7	0.651	1.3082	6.2
8	7	8	0.1008	0.20256	0.96
						9	0	9	0.29436	0.90761	4.46
10	0	10	1.3068	4.0293	19.8
						11	10	11	5.733	11.5206	27.3
12	11	12	0.2352	0.47264	1.12
						13	11	13	0.021	0.0422	0.1
14	13	14	4.815	4.9862	10.7
						15	14	15	22.95	23.766	51
16	15	16	1.71	1.7708	3.8
						17	15	17	0.4275	0.4427	0.95
18	15	18	3.375	3.495	7.5
						19	15	19	0.045	0.0466	0.1
20	10	20	17.505	18.1274	38.9
						21	20	21	1.683	1.74284	3.74

22	20	22	0.045	0.0466	0.1
						23	0	23	0.2415	0.4853	2.3
24	0	24	0.378	0.7596	3.6
						25	0	25	2.36145	4.74539	22.49
26	25	26	27.3645	28.33746	60.81
						27	25	27	0.67848	2.09198	5.14
28	25	28	6.732	6.97136	14.96
						29	25	29	5.715	5.9182	12.7
30	29	30	5.445	5.6386	12.1

No. 2 radial distribution sub-network node power information of table 2-3

Sequence number	P(kw)	Q(kvar)
			0	0	0
1	0	0
			2	0	0
3	0	0
			4	0	0
5	0	0
			6	0	0
7	0	0
			8	0	0
9	0	0
			10	0	0

No. 2 radial distribution sub-network branch roads of table 2-4 connect and impedance information

Sequence number	First node	End-node	r(Ω)	x(Ω)	Length (km)
						1	0	1	9.2502	14.69754	34.26
2	0	2	14.1102	22.41954	52.26
						3	0	3	2.2734	3.61218	8.42
4	0	4	5.28	6.96	16
						5	0	5	7.986	10.527	24.2
6	5	6	19.6482	25.8999	59.54
						7	6	7	0.033	0.0435	0.1
8	6	8	0.033	0.0435	0.1
						9	8	9	14.4	14.912	32
10	0	10	16.686	26.5122	61.8

No. 3 radial distribution sub-network node power information of table 2-5

Sequence number	P(kw)	Q(kvar)
			0	0	0
1	0	0
			2	0	0
3	0	0
			4	0	0
5	5.051762	1.024608
			6	5.051762	1.024608
7	32.34197	9.666952
			8	21.38312	4.338991

9	6.949514	1.407722
			10	32.34197	9.666952
11	0	0
			12	33.67841	6.842598
13	16.83921	3.421299
			14	2.753077	0.561307
15	6.735682	1.372083
			16	2.218499	0.454391
17	0	0
			18	0.534578	0.106916
19	0	0
			20	21.38312	4.338991
21	16.03734	3.260926

No. 3 radial distribution sub-network branch roads of table 2-6 connect and impedance information

Sequence number	First node	End-node	r(Ω)	x(Ω)	Length (km)
						1	0	1	11.1006	22.30692	52.86
2	1	2	6.0984	12.25488	29.04
						3	2	3	5.5704	7.3428	16.88
4	1	4	0.021	0.42201	0.1
						5	0	5	6.39	6.6172	14.2
6	5	6	4.644	4.80912	10.32
						7	0	7	0.99015	1.98973	9.43
8	0	8	2.1483	4.31706	10.23

9	8	9	0.4011	0.80602	1.91
						10	0	10	2.0223	4.06386	9.63
11	0	11	0.7854	2.42165	11.9
						12	11	12	0.08925	0.17935	0.425
13	11	13	0.021	0.0422	0.1
						14	13	14	5.691	11.4362	27.1
15	14	15	4.1706	8.38092	19.86
						16	15	16	9.855	10.2054	21.9
17	15	17	9.372	12.354	28.4
						18	17	18	20.995	14.4058	32.3
19	17	19	0.021	0.0422	0.1
						20	0	20	0.48906	1.507935	7.41
21	20	21	0.8463	1.70066	8.06

No. 4 radial distribution sub-network node power information of table 2-7

Sequence number	P(kw)	Q(kvar)
			0	0	0
1	0	0
			2	9.221333	1.87842
3	6.07437	1.231951
			4	4.610667	0.93921
5	0	0
			6	0	0
7	0	0

8	0	0
			9	11.92919	2.427309
10	6.916	1.402716
			11	2.305333	0.463506
12	9.221333	1.87842
			13	21.95556	4.464296
14	21.95556	4.464296
			15	4.610667	0.93921

No. 4 radial distribution sub-network branch roads of table 2-8 connect and impedance information

Sequence number	First node	End-node	r(Ω)	x(Ω)	Length (km)
						1	0	1	3.7695	7.5749	35.9
2	1	2	0.0105	0.0211	0.1
						3	1	3	0.0105	0.0211	0.1
4	3	4	6.2397	9.91419	23.11
						5	0	5	9.775	23.92	57.5
6	5	6	0.017	0.0446	0.1
						7	5	7	0.017	0.0446	0.1
8	0	8	5.882	14.3936	34.6
						9	0	9	10.4148	13.7286	31.56
10	0	10	6.534	10.3818	24.2
						11	10	11	8.1648	12.97296	30.24
12	10	12	4.725	7.5075	17.5
						13	0	13	1.6485	3.3127	7.85

14	0	14	0.4935	0.9917	2.35
						15	0	15	6.0984	12.25488	29.04

No. 5 radial distribution sub-network node power information of table 2-9

Sequence number	P(kw)	Q(kvar)
			0	0	0
1	11.09266	2.249462
			2	8.435484	1.715215
3	13.28589	2.699355
			4	2.108871	0.428804
5	8.435484	1.715215
			6	6.326613	1.286411
7	0	0
			8	0	0

No. 5 radial distribution sub-network branch roads of table 2-10 connect and impedance information

Sequence number	First node	End-node	r(Ω)	x(Ω)	Length (km)
						1	0	1	0.41643	0.836826	1.983
2	0	2	0.51	1.248	6
						3	0	3	1.7688	5.4538	13.4
4	0	4	1.09725	2.20495	10.45
						5	4	5	0.693	1.3926	3.3
6	5	6	3.213	6.4566	15.3
						7	4	7	9.36	9.6928	20.8
8	4	8	3.555	3.6814	7.9

No. 6 radial distribution sub-network node power information of table 2-11

Sequence number	P(kw)	Q(kvar)
			0	0	0
1	20.70345	6.866644
			2	1.552759	0.517586
3	3.260793	1.078305
			4	0	0
5	2.665569	0.879897
			6	1.552759	0.517586
7	0	0
			8	1.630397	0.327805
9	2.665569	0.543466
			10	1.630397	0.327805
11	0	0
			12	1.035172	0.207034
13	1.035172	0.207034
			14	1.630397	0.327805
15	2.070345	0.422695
			16	15.52759	3.157276
17	2.562052	0.517586
			18	0	0
19	156.52759	3.157276

No. 6 radial distribution sub-network branch roads of table 2-12 connect and impedance information

Sequence number	First node	End-node	r(Ω)	x(Ω)	Length (km)
						1	0	1	9.8532	19.80024	46.92
2	1	2	18.7155	19.38094	41.59
						3	0	3	4.152	20.2929	51.9
4	3	4	18.423	19.07804	40.94
						5	4	5	0.045	0.0466	0.1
6	4	6	0.045	0.0466	0.1
						7	0	7	6.1929	12.44478	29.49
8	7	8	5.04	10.128	24
						9	7	9	4.6599	9.36418	22.19
10	9	10	16.0355	11.00282	24.67
						11	9	11	3.5028	7.03896	16.68
12	11	12	5.586	11.2252	26.6
						13	11	13	0.231	0.4642	1.1
14	11	14	5.1954	10.44028	24.74
						15	11	15	10.1725	6.9799	15.65
16	11	16	4.1097	8.25854	19.57
						17	16	17	7.7545	5.32078	11.93
18	0	18	1.491	2.9962	7.1
						19	0	19	1.0878	2.18596	5.18

No. 7 radial distribution sub-network node power information of table 2-13

Sequence number

P(kw)

Q(kvar)

0	0	0
			1	16.81635	3.417299
2	8.104265	1.647867
			3	0	0
4	1.681635	0.344431
			5	3.829265	0.776659
6	2.086848	0.425474
			7	0	0
8	2.552844	0.520024
			9	3.829265	0.776659
10	2.43128	0.493009
			11	8.104265	1.647867
12	1.620853	0.330924
			13	3.829265	0.776659
14	0	0
			15	2.086848	0.425474
16	3.829265	0.776659
			17	1.681635	0.344431
18	2.086848	0.425474
			19	3.829265	0.776659

No. 7 radial distribution sub-network branch roads of table 2-14 connect and impedance information

Sequence number	First node	End-node	r(Ω)	x(Ω)	Length (km)
						1	0	1	1.155	2.321	5.5

2	0	2	14.445	14.9586	32.1
						3	0	3	0.021	0.0422	0.1
4	3	4	0.021	0.0422	0.1
						5	3	5	5.628	11.3096	53.6
6	5	6	14.445	14.9586	32.1
						7	5	7	0.045	0.0466	0.1
8	7	8	0.045	0.0466	0.1
						9	7	9	0.045	0.0466	0.1
10	7	10	18.945	19.6186	42.1
						11	0	11	1.2474	2.50668	11.88
12	0	12	6.4554	12.97228	30.74
						13	12	13	4.86	5.0328	10.8
14	12	14	17.04294	27.07934	63.122
						15	14	15	0.045	0.0466	0.1
16	14	16	0.045	0.0466	0.1
						17	12	17	4.956	9.9592	23.6
18	0	18	8.3685	16.8167	39.85
						19	18	19	3.6519	7.33858	17.39

No. 8 radial distribution sub-network node power information of table 2-15

Sequence number	P(kw)	Q(kvar)
			0	0	0
1	1.269488	0.253898
			2	1.269488	0.253898

3	1.999443	0.402004
			4	3.268931	0.666481
5	1.999443	0.402004
			6	12.69488	2.581292
7	1.999443	0.402004
			8	3.998886	0.814588
9	0	0

No. 8 radial distribution sub-network branch roads of table 2-16 connect and impedance information

Sequence number	First node	End-node	r(Ω)	x(Ω)	Length (km)
						1	0	1	14.58	23.166	54
2	1	2	28.755	29.7774	63.9
						3	0	3	7.98	16.036	38
4	0	4	6.174	12.4068	29.4
						5	4	5	9.135	9.4598	20.3
6	0	6	0.25095	0.50429	2.39
						7	0	7	8.4483	16.97706	40.23
8	0	8	4.977	10.0014	23.7
						9	8	9	14.202	22.5654	52.6

No. 9 radial distribution sub-network node power information of table 2-17

Sequence number	P(kw)	Q(kvar)
			0	0	0
1	6.425424	1.306215
			2	0.907627	0.182486

3	2.881356	0.585876
			4	1.152542	0.235311
5	3.630508	0.739548
			6	1.152542	0.235311
7	1.152542	0.235311

No. 9 radial distribution sub-network branch road information of table 2-18

Sequence number	First node	End-node	r(Ω)	x(Ω)	Length (km)
						1	0	1	1.9425	3.9035	9.25
2	0	2	13.095	13.5606	29.1
						3	0	3	8.169	16.4158	38.9
4	0	4	2.5704	5.16528	12.24
						5	4	5	11.817	12.23716	26.26
6	5	6	1.0611	1.68597	3.93
						7	6	7	1.0611	1.68597	3.93

Through sending out 9 node voltages of transmission ﹠ distribution electric power networks after this electric power system tide analytical equipment analysis computing, circuit through-put power and line loss see Table 3-1 respectively to table 3-20.

Show radial distribution sub-network node voltage 3-11 number

Sequence number	Voltage real part (kv)	Voltage imaginary part (kv)	Magnitude of voltage (kv)
				0	66	0	66
1	66	0	66
				2	66	0	66
3	66	0	66
				4	66	0	66

5	66	0	66
				6	66	0	66
7	66	0	66
				8	66	0	66
9	66	0	66
				10	66	0	66
11	66	0	66
				12	66	0	66
13	66	0	66
				14	66	0	66
15	66	0	66
				16	66	0	66
17	66	0	66
				18	66	0	66
19	66	0	66
				20	66	0	66
21	66	0	66
				22	66	0	66
23	65.85446	-0.18677	65.85387
				24	65.7722	-0.29234	65.77076
25	64.87888	-1.43801	64.78887
				26	64.87888	-1.43801	64.78887
27	64.87888	-1.43801	64.78887

28	64.05504	-2.00497	63.93786
				29	63.24628	-2.56039	63.06459
30	62.57994	-3.01896	62.36075

Table No. 1 radial distribution sub-network circuit through-put power of 3-2 and line loss

Sequence number	Circuit is meritorious	Circuit is idle	Active loss (kw)	Reactive loss
					1	-4.3E+08	-4.3E+08	0	0
2	-4.3E+08	-4.3E+08	0	0
					3	-4.3E+08	-4.3E+08	0	0
4	-4.3E+08	-4.3E+08	0	0
					5	-4.3E+08	-4.3E+08	0	0
6	-4.3E+08	-4.3E+08	0	0
					7	-4.3E+08	-4.3E+08	0	0
8	-4.3E+08	-4.3E+08	0	0
					9	-4.3E+08	-4.3E+08	0	0
10	-4.3E+08	-4.3E+08	0	0
					11	-4.3E+08	-4.3E+08	0	0
12	-4.3E+08	-4.3E+08	0	0
					13	-4.3E+08	-4.3E+08	0	0
14	-4.3E+08	-4.3E+08	0	0
					15	-4.3E+08	-4.3E+08	0	0
16	-4.3E+08	-4.3E+08	0	0
					17	-4.3E+08	-4.3E+08	0	0

18	-4.3E+08	-4.3E+08	0	0
					19	-4.3E+08	-4.3E+08	0	0
20	-4.3E+08	-4.3E+08	0	0
					21	-4.3E+08	-4.3E+08	0	0
22	0	0	0	0
					23	28.29974	5.826374	0.046284	0.093008
24	28.32543	5.878959	0.072627	0.145946
					25	22.97766	5.636083	0	0
26	0	0	0	0
					27	0	0	0	0
28	6.735649	1.428679	0.076389	0.079105
					29	15.93858	3.597647	0.118542	0.122757
30	6.708024	1.41481	0.06495	0.067259

No. 2 radial distribution sub-network node voltages of table 3-3

Sequence number	Voltage real part (kv)	Voltage imaginary part (kv)	Magnitude of voltage (kv)
				0	66	0	66
1	66	0	66
				2	66	0	66
3	66	0	66
				4	66	0	66
5	66	0	66
				6	66	0	66
7	66	0	66

8	66	0	66
				9	66	0	66
10	66	0	66

Table No. 2 radial distribution sub-network circuit through-put powers of 3-4 and line loss

Sequence number	Circuit is meritorious	Circuit is idle	Active loss	Reactive loss
					1	0	0	0	0
2	0	0	0	0
					3	0	0	0	0
4	0	0	0	0
					5	0	0	0	0
6	0	0	0	0
					7	0	0	0	0
8	0	0	0	0
					9	0	0	0	0
10	0	0	0	0

No. 3 radial distribution sub-network node voltages of table 3-5

Sequence number	Voltage real part (kv)	Voltage imaginary part (kv)	Magnitude of voltage (kv)
				0	66	0	66
1	66	0	66
				2	66	0	66
3	66	0	66
				4	66	0	66
5	64.81634	-0.81458	64.78368

6	64.38622	-1.11059	64.34214
				7	65.22337	-0.83	65.20884
8	64.70188	-1.66618	64.65361
				9	64.64246	-1.7425	64.59291
10	64.41378	-1.69521	64.35244
				11	64.78503	-2.15092	64.69608
12	64.72089	-2.23319	64.63066
				13	64.77199	-2.16763	64.68229
14	63.28414	-4.07374	63.05639
				15	62.43903	-5.15648	62.15258
16	62.0375	-5.43167	61.72623
				17	62.3431	-5.24136	62.05087
18	62.14971	-5.32403	61.84615
				19	62.3431	-5.24136	62.05087
20	65.54907	-0.79865	65.53828
				21	65.25941	-1.17008	65.24429

Table No. 3 radial distribution sub-network circuit through-put powers of 3-6 and line loss

Sequence number	Circuit is meritorious	Circuit is idle	Active loss	Reactive loss
					1	0	0	0	0
2	0	0	0	0
					3	0	0	0	0
4	0	0	0	0

5	10.28882	2.247055	0.040454	0.041893
					6	5.077301	1.054553	0.029805	0.030865
7	32.59727	10.19679	0.265329	0.533185
					8	28.73536	6.611791	0.244667	0.491664
9	6.9465	1.415838	0.004833	0.009713
					10	32.85495	10.76962	0.556432	1.118163
11	63.71408	15.84645	0	0
					12	33.63631	6.87755	0.025235	0.05071
13	29.2982	6.568394	0.001482	0.002978
					14	12.48892	3.148351	0.011299	0.022706
15	9.525848	2.136312	0.051015	0.102516
					16	2.211801	0.464822	0.013264	0.013736
17	0.532301	0.108251	0	0
					18	0.531585	0.10731	0.001631	0.001119
19	0	0	0	0
					20	37.62881	8.21803	0.054205	0.167132
21	16.08417	3.366579	0.053247	0.107002

No. 4 radial distribution sub-network node voltages of table 3-7

Sequence number	Voltage real part (kv)	Voltage imaginary part (kv)	Magnitude of voltage (kv)
				0	66	0	66
1	64.3983	-2.05339	64.31755

2	64.39623	-2.05604	64.31543
				3	64.3959	-2.05646	64.31508
4	63.81892	-2.66026	63.71586
				5	66	0	66
6	66	0	66
				7	66	0	66
8	66	0	66
				9	63.61267	-2.09835	63.48855
10	63.58514	-2.53031	63.40283
				11	63.20885	-2.92611	63.00899
12	62.71131	-3.44476	62.47729
				13	65.22754	-0.9905	65.21087
14	65.76875	-0.29652	65.76727
				15	65.39958	-0.76933	65.38955

Table No. 4 radial distribution sub-network circuit through-put powers of 3-8 and line loss

Sequence number	Circuit is meritorious	Circuit is idle	Active loss	Reactive loss
					1	20.27858	4.852702	0	0
2	9.20703	1.878916	0.000225	0.000452
					3	10.693	2.214242	9.75E-05	0.000196
4	4.633058	0.990807	0.034029	0.054069
					5	0	0	0	0
6	0	0	0	0

7	0	0	0	0
					8	0	0	0	0
9	12.28045	2.923491	0.382914	0.50475
					10	19.05654	4.835515	0.080943	0.12861
11	2.306854	0.479481	0.011372	0.018068
					12	9.279222	2.036969	0.107202	0.170332
13	22.14186	4.853378	0.194595	0.391044
					14	22.01215	4.579309	0.057273	0.115092
15	4.641198	1.002414	0.031578	0.063457

No. 5 radial distribution sub-network node voltages of table 3-9

Sequence number	Voltage real part (kv)	Voltage imaginary part (kv)	Magnitude of voltage (kv)
				0	66	0	66
1	65.90149	-0.12645	65.90122
				2	65.90238	-0.14625	65.90208
3	65.42088	-1.02551	65.40772
				4	65.60491	-0.5066	65.59487
5	65.38657	-0.78656	65.37223
				6	64.95274	-1.34285	64.92905
7	65.60491	-0.5066	65.59487
				8	65.60491	-0.5066	65.59487

Table No. 5 radial distribution sub-network circuit through-put powers of 3-10 and line loss

Sequence number

Circuit is meritorious

Circuit is idle

Active loss

Reactive loss

1	11.10485	2.274067	0.012284	0.024684
					2	8.444086	1.7364	0.008701	0.021293
3	13.35754	2.932676	0.075992	0.23431
					4	17.00951	3.720585	0.001181	0.002373
5	14.82469	3.138453	0.012016	0.024147
					6	6.354681	1.349412	0.031767	0.063836
7	0	0	0	0
					?8	0	0	0	0

No. 6 radial distribution sub-network node voltages of table 3-11

Sequence number	Voltage real part (kv)	Voltage imaginary part (kv)	Magnitude of voltage (kv)
				0	66	0	66
1	59.32967	-5.5766	59.56949
				2	58.63705	-5.85603	58.90511
3	64.63438	-2.166	64.67021
				4	62.93925	-2.95611	63.00778
5	62.93664	-2.95733	63.00522
				6	62.93773	-2.95682	63.00628
7	60.0746	-5.0547	60.10262
				8	59.86438	-5.28404	59.91003
9	55.81035	-8.64612	56.15331
				10	55.25759	-8.79174	55.62103
11	53.04161	-10.9258	53.73813
				12	52.86018	-11.0856	53.58766

13	53.03414	-10.9324	53.73194
				14	52.77488	-11.1595	53.51686
15	52.57367	-11.0257	53.29076
				16	50.53983	-12.945	51.66635
17	50.07942	-13.024	51.22886
				18	66	0	66
19	65.63424	-0.46225	65.63587

Table No. 6 radial distribution sub-network circuit through-put powers of 3-12 and line loss

Sequence number	Circuit is meritorious	Circuit is idle	Active loss	Reactive loss
					1	23.79058	10.46722	1.32111	2.654802
2	1.566572	0.532342	0.01445	0.014964
					3	7.634269	2.880977	0.01171	0.057234
4	4.310035	1.492496	0	0
					5	2.664392	0.878829	8.93E-05	9.25E-05
6	1.552172	0.517085	3.04E-05	3.14E-05
					7	34.22872	15.35393	0	0
8	1.629045	0.335216	0.003884	0.007804
					9	30.59885	10.99801	0.010937	0.021978
10	1.634752	0.336343	0.014335	0.009836
					11	24.95014	7.379523	0	0
12	1.029176	0.209889	0.002168	0.004356
					13	1.027149	0.205699	8.92E-05	0.000179

14	1.622442	0.335417	0.005017	0.010082
					15	2.069531	0.430702	0.015994	0.010974
16	18.44974	4.686487	0.386542	0.776766
					17	2.556299	0.525831	0.020187	0.013851
18	0	0	0	0
					19	15.59102	3.284736	0.063397	0.127398

No. 7 radial distribution sub-network node voltages of table 3-13

Sequence number	Voltage real part (kv)	Voltage imaginary part (kv)	Magnitude of voltage (kv)
				0	66	0	66
1	65.58554	-0.53157	65.58077
				2	63.85279	-1.47614	63.76355
3	65.99265	-0.00943	65.99212
				4	65.99189	-0.0104	65.99136
5	64.22395	-2.27832	64.12233
				6	63.67078	-2.65818	63.54819
7	64.21667	-2.28332	64.11482
				8	64.21457	-2.28477	64.11265
9	64.21352	-2.2855	64.11158
				10	63.37224	-2.86451	63.23333
11	65.78424	-0.27666	65.78295
				12	64.20213	-2.30495	64.03925
13	63.86093	-2.53976	63.68644

14	62.1812	-4.42187	61.86011
				15	62.17948	-4.42306	61.85829
16	62.17804	-4.42405	61.85677
				17	64.02388	-2.53284	63.85523
18	64.94357	-1.35499	64.90377
				19	64.64532	-1.7378	64.59822

Table No. 7 radial distribution sub-network circuit through-put powers of 3-14 and line loss

Sequence number	Circuit is meritorious	Circuit is idle	Active loss	Reactive loss
					1	16.89362	3.575782	0.07908	0.158913
2	8.332959	1.895917	0.242993	0.251633
					3	16.73123	4.000023	0	0
4	1.681494	0.344172	1.42E-05	2.86E-05
					5	15.04943	3.655231	0.020897	0.041992
6	2.097115	0.441154	0.016225	0.016802
					7	8.822969	1.817708	0	0
8	2.549887	0.521341	7.43E-05	7.69E-05
					9	3.818398	0.772365	0.000167	0.000173
10	2.452374	0.521792	0.029159	0.030196
					11	8.123757	1.687453	0.019715	0.039618
12	13.44588	3.495141	0.004308	0.008657
					13	3.833938	0.792838	0.018293	0.018943

14	6.030727	1.447438	0	0
					15	2.069696	0.421553	5.33E-05	5.52E-05
16	3.798845	0.769763	0.00018	0.000186
					17	1.679468	0.350289	0.003581	0.007197
18	5.996243	1.373822	0.009011	0.018108
					19	3.838426	0.802597	0.01336	0.026848

No. 8 radial distribution sub-network node voltages of table 3-15

Sequence number	Voltage real part (kv)	Voltage imaginary part (kv)	Magnitude of voltage (kv)
				0	66	0	66
1	65.26088	-0.779	65.24116
				2	64.59323	-1.24114	64.55738
3	65.66058	-0.4372	65.65736
				4	65.30631	-0.89041	65.29028
5	64.97195	-1.12135	64.95025
				6	65.93201	-0.08718	65.93188
7	65.64066	-0.46285	65.63705
				8	65.575	-0.54455	65.56999
9	65.575	-0.54455	65.56999

Table No. 8 radial distribution sub-network circuit through-put powers of 3-16 and line loss

Sequence number	Circuit is meritorious	Circuit is idle	Active loss	Reactive loss
					1	2.572289	0.556302	0.005741	0.009122

2	1.280089	0.265692	0.011564	0.011975
					3	2.007003	0.417453	0.0077	0.015472.
4	5.317288	1.161758	0.01612	0.032394
					5	2.007468	0.411167	0.009007	0.009327
6	12.70445	2.60059	0.009688	0.019469
					7	2.007443	0.418369	0.008156	0.016391
8	4.017713	0.853222	0.019279	0.038742
					9	0	0	0	0

No. 9 radial distribution sub-network node voltages of table 3-17

Sequence number	Voltage real part (kv)	Voltage imaginary part (kv)	Magnitude of voltage (kv)
				0	66	0	66
1	66	0	66
				2	65.78242	-0.15028	65.78154
3	65.49764	-0.64415	65.49062
				4	65.61082	-0.49844	65.59392
5	64.32371	-1.38229	64.26645
				6	64.27463	-1.4336	64.21603
7	64.25008	-1.45926	64.19082

Table No. 9 radial distribution sub-network circuit through-put powers of 3-18 and line loss

Sequence number	Circuit is meritorious	Circuit is idle	Active loss	Reactive loss
					1	0	0	0.019172	0.038526
2	0.910199	0.185161	0.002594	0.002686

3	2.897371	0.618867	0.016466	0.03309
					4	7.219962	1.620361	0.000827	0.001661
5	6.035414	1.320235	0.039276	0.040673
					6	2.30413	0.472865	0.000356	0.000566
7	1.151545	0.235619	0.000356	0.000566

Table 3-19 sends out each node voltage of electric power transmission network and power

Sequence number	The voltage real part	The voltage imaginary part	Magnitude of voltage (kv)	Given meritorious (kw)	Given idle
						1	229.9	0	229.9	0	16.29
2	214.5478	-42.7512	218.7657	-0.50412	35.803
						3	211.6185	-26.3108	213.2478	-38	-1.5.6
4	223.9868	-20.5539	224.9279	-0.153	3.56299
						5	217.3416	-31.2291	219.5737	-52.3	-20.9
6	223.7557	-22.6464	224.8988	-150	0
						7	228.8323	-2.19123	228.8428	-74	-10.54
8	224.2119	1.288711	224.2156	-0.2185	6.386
						9	210.2279	-23.4406	211.5306	-214	-85.5
10	216.6478	-7.96057	216.794	-0.0505	5.2242
						11	202.3987	-30.6858	204.7117	-79	-31.6
12	224.8335	5.301418	224.896	711.56	0
						13	204.0995	-13.5766	204.5506	-100	-63.43
14	221.9604	2.794573	221.978	-0.452	14.616
						15	180.6137	-50.0768	187.4273	-238	-95.5
16	226.7681	-5.94215	226.8459	-1162	-231.78
						17	233.3963	-3.03649	233.4161	-10	36.8

18	222.3895	3.488933	222.4169	250.96	89.14
						19	221.7947	2.632355	221.8103	-203.56	-82.51
20	215.967	-15.1118	216.495	-0.055	10.771
						21	193.5568	-45.7007	198.8788	-104	-41.5
22	206.4226	-24.991	207.9299	-0.044	6.678
						23	189.7232	-46.9516	195.4465	-72	-28.8
24	228.0101	0.380291	228.0104	-0.05219	24.448
						25	228.6149	-3.9053	228.6483	-17	3.12
26	227.4426	-1.61637	227.4483	-10	-0.66
						27	222.2125	-13.1555	222.6015	-0.05219	23.818
28	218.1132	-20.4781	219.0724	-30	-12.18
						29	194.4024	-44.5288	199.437	-103	-23.08

Table 3-20 sends out electric power transmission network each circuit through-put power and line loss

Sequence number	Circuit is meritorious	Circuit is idle	Active loss	Reactive loss	Resistance (Ω)	Reactance (Ω)
							1	10.23981	36.97864	-0.41566	-2.0627169	15.4696	76.76789
2	204.1332	128.9283	-1.26285	-3.8938034	1.095725	3.378487
							3	289.1431	-29.6229	-1.35709	-5.0605633	0.815126	3.039582
4	148.0712	-35.496	-8.41754	-31.254008	19.4022	72.03965
							5	174.0482	105.0447	-3.64333	-17.806788	4.65952	22.7734
6	55.00379	18.92285	-0.78261	-3.8249909	11.91976	58.25783
							7	79.30116	42.24065	-0.95977	-4.6908958	5.9768	29.21161
8	169.5694	-22.2023	-6.3481	-31.502466	11.36688	56.40814
							9	58.36816	20.30344	-0.68323	-3.3905329	9.3688	46.49267

10	38.83605	63.65999	-1.79277	-6.6564923	15.4292	57.28806
							11	68.58072	37.15364	-0.99202	-4.848497	7.64272	37.35379
12	10.00305	-0.68102	-0.01834	-0.0896435	9.43896	46.13292
							13	103.3876	79.24608	-3.25568	-15.912152	9.704	47.4283
14	149.2774	29.1285	-0.28538	-1.3948063	0.624	3.0498
							15	27.01868	26.37521	-0.19443	-0.9502655	7.09008	34.65277
16	7.306944	186.2894	-1.22504	-4.5681217	1.86287	6.94659
							17	251.9595	-35.8027	-6.93662	-33.902744	5.414	26.46093
18	205.3053	-9.57475	-3.41898	-16.710245	4.23864	20.71635
							19	81.81583	-25.9537	-0.08159	-0.3042447	0.54784	2.04288
20	168.5428	63.64061	-0.21904	-0.8167891	0.333091	1.242087
							21	34.21965	-20.8775	-0.01078	-0.0404981	0.33	1.24
22	1171.08	265.6384	-9.08757	-33.88729	0.333091	1.242087
							23	427.8108	82.90133	-1.19673	-4.4625892	0.333091	1.242087

Claims (3)

1. a self adaptation is sent out transmission ﹠ distribution overall situation current analytical device, comprise host computer and slave computer, it is characterized in that slave computer comprises preposition module, the A/D modular converter, the DSP module, keyboard and LCD MODULE, the trend feedback module, communication module, preposition module comprises voltage sensor and power sensor, the voltage sensor that is connected in the preposition module of this device is connected the A/D modular converter respectively with power sensor, the A/D modular converter connects the DSP module successively, communication module and host computer, keyboard is connected the DSP module with LCD MODULE, the trend feedback module connects the DSP module; According to the characteristics of actual electric network the set of node of overall electric power system is divided into three parts, to overall electric power system, with its set of node carry out following division: interconnected system G, main system M, from the S of system, wherein interconnected system G by the load bus of generating and transmitting system just the root node of distribution system form, its set of node is designated as C _G, the node number is N _GMain system M is made up of the node of remaining generating and transmitting system, and its set of node is designated as C _M, the node number is N _MBe made up of the node of remaining distribution system from the S of system, its set of node is designated as C _S, the node number is N _SProvide overall power flow equation and carry out the iterative analysis computing in the DSP module, the primary iteration form that broad sense self adaptation overall situation tidal current analysis calculates is as follows:

Step 1: boundary system voltage V _GInitialize: V _G ^(k), k=0;

Step 2: with boundary system voltage V _G ^(k)Be reference voltage, in conjunction with Δ V _G ^(k)Find the solution the distribution power flow equation, wherein Δ V _G ^(k)Be the polar coordinates of the voltage correction of the k time interconnected system,

S _S(V _S)＝S _SG(V _G+ΔV _G，V _S)+S _SS(V _S) (1)

Get the polar coordinate representation V of distribution system voltage vector _S ^(k+1), and by V _G ^(k), Δ V _G ^(k)And V _S ^(k+1), calculate the middle associated variable S of iteration _GS ^(k+1)

Step 3: by associated variable S in the middle of the iteration _GS ^(k+1), find the solution and send out the transmission of electricity power flow equation

$\left\{\begin{array}{c}{S}_M\left(V_M\right)=S_{\mathrm{MM}}\left(V_M\right)+S_{\mathrm{MG}}(V_M,V_G+{\mathrm{\&Delta;V}}_G)\\ S_G(V_G+\mathrm{\&Delta;}{V}_G)=S_{\mathrm{GM}}(V_M,V_G+{\mathrm{\&Delta;V}}_G)+S_{\mathrm{GG}}(V_G+{\mathrm{\&Delta;V}}_G)+S_{\mathrm{GS}}(V_G+{\mathrm{\&Delta;V}}_G,V_S)\end{array}\right.---\left(2\right)$

Get the generating and transmitting system voltage vector

With Δ V _G ^(k+1)

Step 4: the maximum max| Δ V that judges the mold component of boundary system voltage difference between adjacent twice iteration _i| whether less than given convergence index ε, wherein i ∈ C _B, the load bus of generating and transmitting system is formed boundary system B, and its set of node is designated as C _BIf,, then overall trend iteration convergence; Otherwise k=k+1 changes step 2;

Illustrate: in (1) (2) two formulas, S _M, S _GAnd S _SBe respectively the polar coordinate representation of the node injection complex power vector of corresponding node collection: S _XYBe set of node C _XGo up each node and directly flow to set of node C _YThe polar coordinate representation of the vector formed of branch road complex power trend; S _XXBe set of node C _XThe polar coordinate representation of the vector that the branch road complex power trend that last each node flows directly into set of node self is formed.

2. send out transmission ﹠ distribution overall situation current analytical device by the described self adaptation of claim 1, it is characterized in that described trend feedback module comprises D/A conversion chip, operational amplifier and A/D conversion chip, the output concatenation operation amplifier input terminal of D/A conversion chip wherein, the output of operational amplifier connects the input of A/D conversion chip.

3. send out transmission ﹠ distribution overall situations current analytical device by the described self adaptation of claim 1 and carry out self adaptation overall situation tidal current analysis method, it is characterized in that this analytical method carries out as follows:

Step 1: initialization;

Step 2:, comprise burden with power P, the voltage magnitude V signal of master-slave system PV child node, the burden with power P of master-slave system PQ child node, load or burden without work Q signal with preposition module sensors collection site data-signal;

Step 3: analog signal is converted into digital signal by the A/D modular converter;

Step 4: to the digital signal of DSP module input A/D modular converter output, analyze by self adaptation overall situation tidal current analysis method, the initial value of at first given boundary node voltage, the DSP module invokes pushes back for the tidal current analysis module before components of system as directed, setting is carried out iterative analysis from the system voltage initial value and to the input data from system, send the trend feedback module with what it analyzed to from the system node information of voltage, and pushed back the node voltage signal that goes out for Algorithm Analysis before the trend feedback module storage last time and use as voltage initial value before pushing back when analyzing for method it to be sent back push back for module before the system next time, realize pushing back optimization before the tradition with this for the tidal current analysis method; This calculates and finishes the result of calculation of back DSP module utilization from system from system, according to formula S _S(V _S)=S _SG(V _G+ Δ V _G, V _S)+S _SS(V _S) try to achieve the iteration intermediate variable S between master-slave system _GS ^(k+1)And the inferior method tidal current analysis of the newton-pressgang module of calling main system part is carried out analytical calculation to the digital input signals of main system, try to achieve the voltage vector of boundary node and obtain the voltage difference of boundary node, judge the maximum max| Δ V of the mold component of master-slave network interconnected system voltage difference between adjacent twice iteration _i| whether less than given convergence index ε, wherein i ∈ C _G, ε=0.001, if, then overall trend iteration convergence; Otherwise, k=k+1, the analytic process of repeating step four is till satisfying the condition of convergence;

Step 5: analyze and finish, the output analysis result, analysis result comprises the idle Q and the voltage phase angle value of master-slave system PV child node, the voltage magnitude V and the voltage phase angle of PQ child node, the active power P of boundary point, reactive power Q voltage magnitude V and voltage phase angle signal, and the active power of main system balance node is injected P, reactive power is injected Q;

Step 6: communication module sends above tidal current analysis result to host computer, and shows this tidal current analysis result on display;

Step 7: end of run.

Images