System and method for monitoring and managing three-phase power flows in electrical transmission and distribution networks

1. A computer implemented method, comprising:

a. for an active power bus in a three-phase power generating system having an electrical grid, said active power bus being modeled in terms of symmetrical components, transforming negative-sequence and zero-sequence admittances into phase-sequence equivalent admittances;

b. for a reactive power load in said three-phase power generating system, said reactive power load being modeled in terms of symmetrical components, transforming negative-sequence and zero-sequence admittances into phase-sequence equivalent admittances;

c. generating a mathematical model of power flow equations L in terms of three phase components of the three-phase power generating system using the phase-sequence equivalent admittances;

d. embedding the mathematical model in a holomorphic embedding of the power flow equations (L(s)), where s is a variable in a complex domain that includes a value s=0 corresponding to a no load case (L(0)), in which there are no non-linear injections at any bus of the electrical grid and the value s=1 corresponding to an objective case (L(1)) representative of the electrical grid in a problem case, and wherein each variable of the power flow equations (L) is contained in L(s) as a function of the variable s by said holomorphic embedding;

e. transcribing the holomorphic embedding into software for use in a computer processor adapted to execute said software;

f. developing in power series, values of unknowns in parameters of the holomorphic embedding (L(s)), using said computer processor, wherein the values of s are in a neighborhood of the value for the no load case of each parameter of the power flow equations;

g. said computer processor evaluating coefficients of the power series by solving resulting linear systems, order after order of the power series, up to an order determined by desired accuracy;

h. said computer processor evaluating values for unknowns in said holomorphic embedding for the problem case, s=1; and

i. computing the power series summation to the desired accuracy.

2. The computer implemented method of claim 1 , wherein said step of computing the power series summation to the desired accuracy further comprises using said computer processor to approximate a solution using n-order algebraic Padé approximants.

3. The computer implemented method of claim 1 , wherein said step of computing the power series summation to the desired accuracy further comprises using said computer processor to approximate a solution using continued fractions.

4. The computer implemented method of claim 1 , wherein the admittance matrix is rearranged, and mutual coupling between different phases are absorbed into the holomorphic embedding.

5. The computer implemented method of claim 4 , wherein the power generating system is decoupled for solution as three separate linear systems.

6. The computer implemented method of claim 1 , wherein said method is employed for either real-time or off-line analysis.

7. The computer implemented method of claim 6 , wherein said method is employed for either balanced or unbalanced electric power systems.

8. The computer implemented method of claim 1 , further comprising:

prior to said embedding, receiving data from a supervisory and data acquisition system representative of conditions of said electrical grid, and

forming said power flow equations (L) from said data.

9. The computer implemented method of claim 8 , further comprising:

repeating said receiving step and steps d through i continuously to provide a continuous, real time measure of the power flow in the power generating system.

10. A computer implemented method, comprising:

transforming negative-sequence and zero-sequence admittances for an active power bus in a three-phase power generating system having an electrical grid into phase-sequence equivalent admittances, said active power bus being modeled in terms of symmetrical components;

transforming negative-sequence and zero-sequence admittances for a reactive power load in said three-phase power generating system into phase-sequence equivalent admittances, said reactive power load being modeled in terms of symmetrical components;

generating a mathematical model of power flow equations for said three-phase power generating system in terms of three-phase components of said three-phase power generating system using said phase-sequence equivalent admittances;

embedding said mathematical model in a holomorphic embedding of said power flow equations using a variable in a complex domain that includes a first value corresponding to a no load case, in which there are no non-linear injections at any bus of said electrical grid and a second value corresponding to an objective case representative of said electrical grid in a problem case, each variable of said power flow equations being contained said holomorphic embedding as a function of said variable;

transcribing said holomorphic embedding into software for use in a computer processor adapted to execute said software;

developing in power series, values of unknowns in parameters of said holomorphic embedding, using said computer processor, said values of said unknowns in said parameters being in a neighborhood of said first value for said no load case of each parameter of said power flow equations;

said computer processor evaluating coefficients of said power series by solving resulting linear systems, order after order of said power series, up to an order determined by desired accuracy;

said computer processor evaluating values for unknowns in said holomorphic embedding for the problem case; and

computing summation of said power series to said desired accuracy.

11. The computer implemented method of claim 10 , said computing said summation of said power series to said desired accuracy further comprising using said computer processor to approximate a solution to said holomorphic embedding using one of n-order algebraic Padé approximants and continued fractions.

12. The computer implemented method of claim 10 , a matrix of said phase-sequence equivalent admittances being rearranged, and mutual coupling between different phases of said three-phase power generating system being absorbed into said holomorphic embedding.

13. The computer implemented method of claim 10 , said three-phase power generating system being decoupled for solution as three separate linear systems.

14. The computer implemented method of claim 10 , further comprising:

receiving data from a supervisory and data acquisition system representative of conditions of said electrical grid, and

forming said power flow equations from said data.

15. The computer implemented method of claim 14 , further comprising:

continuously providing a measure of power flow in said three-phase power generating system, based on said data from said supervisory and data acquisition system.