System and method for performing a line-wise power flow analysis for a power system

Abstract

Various embodiments are provided for conducting a power flow analysis using a set of line-wise power balance equations. In at least some embodiment, the set of line-wise power balance equations is solved using a Newton-Raphson technique. In various cases, the Jacobian matrix generated by the Newton-Raphson technique may directly indicate the transmission lines, or sets of transmission lines, which are most susceptible to voltage collapse. In at least one example application, the set of line-wise power balance equations may be used as equality constraints in an optimal power flow (OPF) formulation for solving an optimal power flow (OPF) problem.

Description

TITLE: SYSTEM AND METHOD FOR PERFORMING A LINE-WISE POWER FLOW ANALYSIS FOR A POWER SYSTEM FIELD

[0001] The teachings herein relate generally to power systems, and more particularly to 5 various embodiments for a system and method for performing a line-wise power flow analysis for a power system. BACKGROUND

[0002] Power flow analysis is an important tool in many aspects of power system engineering, and constitutes the basis for studies in power system optimization, power 10 systems planning, and voltage stability analysis. In recent years, an increasing, wide-spread effort has been devoted to developing an efficient and accurate method for conducting power flow analysis.

[0003] Generally, power flow analysis is conducted by solving a set of equations to determine the state of a power system. The power system is described by voltage phasors 15 at each power system bus, wherein each voltage phasor comprises magnitude and phase angle components. Further, the sum of all real and reactive powers coming into each bus of the power system is equal to zero.

[0004] A power flow analysis, upon completion, yields voltage phasor values for each bus, that can be used to determine power flow in all lines, system real power loss, reactive 20 power output of all generators, and real power output of a slack generator. SUMMARY OF VARIOUS EMBODIMENTS

[0005] In accordance with a broad aspect of the teachings herein, there is provided at least one embodiment of a method for conducting a power flow analysis for a power system, the power system comprising at least one transmission branch having a first end coupled to 25 a first bus and a second end coupled to a second bus, the method being implemented by at least one server comprising at least one processor, the method comprising: initializing values for a problem vector (W); modelling the power system using a set of line-wise power balance equations that are functions of the problem vector (FT(W) ); solving the set of line-wise power balance equations (FT(W)) using a numerical technique to obtain a line-wise solution vector; - 1 - CA 3028644 2018-12-27 based on the line-wise solution vector, determining whether the power system is in a state of voltage stability; and implementing a remedial action to stabilize the power system when the power system is not in a state of voltage stability. (0006] In at least one of these embodiments, the problem vector (W) is defined as W = 5 [U o PF PS QF QS], where: U is a vector of squared voltage magnitudes for each bus of each of the at least one transmission branch; o is a vector of voltage angles for each bus of each of the at least one transmission branch; PF is a vector of active power flow in each first end of each of the at least one transmission branch; PS is a vector of active power flow in each second end of each of the at least one transmission branch; QF is a vector of reactive 10 power flow in each first end of each of the at least one transmission branch, and QS is a vector of reactive power flowing in each second end of each of the at least one transmission branch. 15 (0007] In at least one of these embodiments, the set of line-wise power balance equations (FT(W)) is defined as: FT(W) = FF(W) FS(W) FA(W) FB(W) FP(W) FQ(W) where: FF (W) is a vector of line based voltage functions in relation to each first end of each of the at least one transmission branch; FS(W) is a vector of line-based voltage function in relation to each second end of each of the at least one transmission branch; FA (W) is a vector of phase angle functions in relation to each first end of each of the at least one 20 transmission branch; FB(W) is a vector of phase angle functions in relation to each second end of each of the at least one transmission branch; FP(W) is a vector of functions for buswise real power injections; and FQ(W) is a vector of functions for bus-wise reactive power injections. (0008] In at least one of these embodiments, the numerical technique for solving the set 25 of line-wise power balance equations is a Newton-Raphson technique. - 2 - CA 3028644 2018-12-27 5

[0009] In at least one of these embodiments, solving the set of line-wise power balance equations using the Newton-Raphson technique comprises: determining an incremental change to the problem vector (LlW) by solving a first-order Taylor's series approximation of the set of line-wise power balance equations (FT(W) ): LlW = [aF;~w)r 1 • [ST - FT(W)] where: LlW is the incremental change to the problem vector; [aF;~w)rl is an inverse of a Jacobian matrix for the set of line-wise power balance equations; and ST is the desired linewise solution vector for the set of line-wise power balance equations.

[0010] In at least one of these embodiments, the desired line-wise solution vector (ST) is 1 O defined according to ST = [O o o O [PD - PG] [QD - QG]f, where: PD is a vector of active power demand in each of the at least one transmission branches; PG is a vector of active power generation in each of the at least one transmission branches; QD is a vector of reactive power demand in each of the at least one transmission branches, and QG is a vector of reactive power generation in each of the at least one transmission branches. 15

[0011] In at least one of these embodiments, the method comprises: generating an updated problem vector (W') according to the equation: W' = W + LlW; solving the set of line-wise power balance equations to generate an intermediate solution vector using the updated problem vector (W'); and based on the intermediate solution vector, determining whether a tolerance condition is satisfied, wherein when the tolerance condition is not 20 satisfied, the method comprises iteratively solving for the updated problem vector using the first-order Taylor's series approximation until the intermediate solution vector satisfies the tolerance condition, and when the tolerance condition is satisfied, the intermediate solution vector comprises the line-wise solution vector.

[0012] In at least one of these embodiments, the tolerance condition is expressed as 25 1ST - FT(W)I ~ T, where Tis a pre-defined tolerance threshold.

[0013] In at least one of these embodiments, the pre-defined tolerance threshold (T) is substantially zero.

[0014] In at least one of these embodiments, the line-wise solution vector yields voltage values at each bus of each of the at least one transmission branch of the power system, and - 3 - CA 3028644 2018-12-27 line-wise real and reactive power flows in each of the at least one transmission branch of the power system.

[0015] In at least one of these embodiments, the set of line-wise based equations is formulated in terms of the square magnitude of the voltages at each bus in each of the at 5 least one transmission branch in the power system.

[0016] In at least one of these embodiments, the functions for bus-wise real power injection, and the functions for bus-wise reactive power injection, are expressed in linear terms and the set of line-wise power balance equations include only first order and second order terms. 10

[0017] In at least one of these embodiments, initializing values for a problem vector (W), further comprises: reading a set of known data values from at least one sensor operatively coupled to the power system.

[0018] In at least one of these embodiments, the transmission branch comprises at least one of a transmission cable, a transmission line, and a transformer and each of the first bus 15 and the second bus comprise at least one of a load bus and a generator bus.

[0019] In at least one of these embodiments, the Jacobian matrix includes, for each of the at least one transmission branches, a first vector element (8FF(w)) and a second vector au element tF:tw)), where the first vector element defines a first voltage collapse index (VCI) for a first end of the at least one transmission branch, and the second vector element defines 20 a second VCI for a second end of the at least one transmission branch.

[0020] In at least one of these embodiments, determining whether the power system is in a state of voltage stability comprises determining whether at least one of the first VCI and the second VCI for a respective transmission branch is substantially zero, and the remedial action to stabilize the power system comprises at least one of: shedding a load bus in the 25 power system, increasing a power output at least one generator bus in the power system, reactivating at least one standby generator bus in the power system, and re-inserting standby transmission elements in the power system. - 4 - CA 3028644 2018-12-27

[0021] In at least one of these embodiments, determining whether the power system is in a state of voltage stability further comprises: performing an N-1 contingency analysis for the power system.

[0022] In accordance with another broad aspect of the teachings herein, there is provided 5 at least one embodiment of a system for conducting a power flow analysis for a power system, wherein the power system comprises at least one transmission branch having a first end coupled to a first bus and a second end coupled to a second bus, the system comprising: at least one server comprising at least one processor, wherein the at least one server is configured for: initializing values for a problem vector (W); modelling the power system using 1 O a set of line-wise power balance equations that are functions of the problem vector (FT(W) ); solving the set of line-wise power balance equations (FT(W)) using a numerical technique to obtain a line-wise solution vector; based on the line-wise solution vector, determining whether the power system is in a state of voltage stability; and implementing a remedial action to stabilize the power system when the power system is not in a state of voltage 15 stability.

[0023] In at least one of these embodiments, the problem vector (W) is defined as W = [U 8 PF PS QF QS], where: U is a vector of squared voltage magnitudes for each bus of each of the at least one transmission branch; 8 is a vector of voltage angles for each bus of each of the at least one transmission branch; PF is a vector of active power flow in each first 20 end of each of the at least one transmission branch; PS is a vector of active power flow in each second end of each of the at least one transmission branch; QF is a vector of reactive power flow in each first end of each of the at least one transmission branch, and QS is a vector of reactive power flowing in each second end of each of the at least one transmission branch. 25

[0024] In at least one of these embodiments, the set of line-wise power balance equations (FT(W)) is defined as: - 5 - CA 3028644 2018-12-27 FT(W) = FF(W) FS(W) FA(W) FB(W) FP(W) FQ(W) where: FF(W) is a vector of line based voltage functions in relation to each first end of each of the at least one transmission branch; FS(W) is a vector of line-based voltage function in relation to each second end of each of the at least one transmission branch; FA (W) is a 5 vector of phase angle functions in relation to each first end of each of the at least one transmission branch; FB(W) is a vector of phase angle functions in relation to each second end of each of the at least one transmission branch; FP(W) is a vector of functions for buswise real power injections; and FQ(W) is a vector of functions for bus-wise reactive power injections. 10

[0025] In at least one of these embodiments, the numerical technique for solving the set of line-wise power balance equations is a Newton-Raphson technique.

[0026] In at least one of these embodiments, solving the set of line-wise power balance equations using the Newton-Raphson technique comprises: determining an incremental change to the problem vector (LlW) by solving a first-order Taylor's series approximation of 15 the set of line-wise power balance equations (FT(W) ): LlW = [°F;:)r 1 · [ST - FT(W)] where: LlW is the incremental change to the problem vector; [°F;~w)r 1 is an inverse of a Jacobian matrix for the set of line-wise power balance equations; and ST is the desired linewise solution vector for the set of line-wise power balance equations. 20

[0027] In at least one of these embodiments, the desired line-wise solution vector (ST) is defined according to ST= [O O O O [PD - PG] [QD - QG]F, where: PD is a vector of active power demand in each of the at least one transmission branches; PG is a vector of active power generation in each of the at least one transmission branches; Q D is a vector of reactive power demand in each of the at least one transmission branches, and QC is a vector of 25 reactive power generation in each of the at least one transmission branches. - 6 - CA 3028644 2018-12-27

[0028] In at least one of these embodiments, the system comprises: generating an updated problem vector (W') according to the equation: W' = W + LlW; solving the set of line-wise power balance equations to generate an intermediate solution vector using the updated problem vector (W'); and based on the intermediate solution vector, determining 5 whether a tolerance condition is satisfied, wherein when the tolerance condition is not satisfied, the method comprises iteratively solving for the updated problem vector using the first-order Taylor's series approximation until the intermediate solution vector satisfies the tolerance condition, and when the tolerance condition is satisfied, the intermediate solution vector comprises the line-wise solution vector. 1 O

[0029] In at least one of these embodiments, the tolerance condition is expressed as 1ST - FT(W)I ~ T, where Tis a pre-defined tolerance threshold.

[0030] In at least one of these embodiments, the pre-defined tolerance threshold (T) is substantially zero.

[0031] In at least one of these embodiments, the line-wise solution vector yields voltage 15 values at each bus of each of the at least one transmission branch of the power system, and line-wise real and reactive power flows in each of the at least one transmission branch of the power system.

[0032] In at least one of these embodiments, the set of line-wise based equations is formulated in terms of the square magnitude of the voltages at each bus in each of the at 20 least one transmission branch in the power system.

[0033] In at least one of these embodiments, the functions for bus-wise real power injection, and the functions for bus-wise reactive power injection, are expressed in linear terms and the set of line-wise power balance equations include only first order and second order terms. 25

[0034] In at least one of these embodiments, initializing values for a problem vector (W), further comprises: reading a set of known data values from at least one sensor operatively coupled to the power system. - 7 - CA 3028644 2018-12-27

[0035] In at least one of these embodiments, the transmission branch comprises at least one of a transmission cable, a transmission line, and a transformer and each of the first bus and the second bus comprise at least one of a load bus and a generator bus.

[0036] In at least one of these embodiments, the Jacobian matrix includes, for each of the 5 at least one transmission branches, a first vector element (°F;;w)) and a second vector element (°Fs(w)), where the first vector element defines a first voltage collapse index (VCI) au for a first end of the at least one transmission branch, and the second vector element defines a second VCI for a second end of the at least one transmission branch.

[0037] In at least one of these embodiments, determining Whether the power system is in 1 O a state of voltage stability comprises determining whether at least one of the first VCI and the second VCI for a respective transmission branch is substantially zero, and the remedial action to stabilize the power system comprises at least one of: shedding a load bus in the power system, increasing a power output at least one generator bus in the power system, reactivating at least one standby generator bus in the power system, and re-inserting standby 15 transmission elements in the power system.

[0038] In at least one of these embodiments, determining whether the power system is in a state of voltage stability further comprises: performing an N-1 contingency analysis for the power system.

[0039] In accordance with another broad aspect of the teachings herein, there is provided 20 at least one embodiment of a method for minimizing an objective function in respect of a power system, the power system comprising at least one transmission branch having a first end coupled to a first bus and a second end coupled to a second bus, the method being implemented using at least one server comprising at least one processor, the method comprising: (a) determining initial values for a control vector (WC); (b) using the control 25 vector (WC), solving a set of line-wise power balance equations (FT(W')) to obtain a problem vector (W'), wherein the problem vector (W') comprises the control vector (WC) and a dependent vector (WD); (c) using the problem vector (W'), solving a line-wise linear optimization formulation to determine an optimal incremental change to the control vector (.1WC); (d) updating the control vector (WC) based on the optimal incremental change in the - 8 - CA 3028644 2018-12-27 control vector (.1WC); and € determining whether the objective function is minimized by determining whether a set of stopping criteria is satisfied based on the updating of the control vector (WC), wherein when the set of stopping criteria is not satisfied, the method comprises iteratively repeating (b) to (d) until the set of stopping criteria is satisfied. 5

[0040] In at least one of these embodiments, the set of line-wise power balance equations (FT(W)) is defined as: FT(W) FF(W) FS(W) FA(W) FB(W) FP(W) FQ(W) wherein: FF(W) is a vector of line based voltage functions in relation to each first end of each of the at least one transmission branch; FS(W) is a vector of line-based voltage function in 10 relation to each second end of each of the at least one transmission branch; FA(W) is a vector of phase angle functions in relation to each first end of each of the at least one transmission branch; FB(W) is a vector of phase angle functions in relation to each second end of each of the at least one transmission branch; FP(W) is a vector of functions for buswise real power injections, and FQ(W) is a vector of functions for bus-wise reactive power 15 injections.

[0041] In at least one of these embodiments, the problem vector (W') is defined as W' = [U,o,PF,PS,QF,QS,PG,QGY, and the dependent vector (WD) is defined as WD = [U, o, PF, PS, QF, QSY, and the control vector (WC) is defined as WC= [U, PGY, wherein: U is a vector of squared voltages magnitudes for each bus of each of the at least one 20 transmission branch; o is a vector of voltage angles for each bus of each of the at least one transmission branch; PF is a vector of active power flow in each first end of each of the at least one transmission branch; PS is a vector of reactive power flow in each second end of each of the at least one transmission branch; QF is a vector of reactive power flow in each first end of each of the at least one transmission branch; QS is a vector of reactive power 25 flowing in each second end of each of the at least one transmission branch; PG is a vector of active power generation in each of the at least one transmission branches, and QG is a vector of reactive power generation in each of the at least one transmission branches. - 9 - CA 3028644 2018-12-27

[0042] In at least one of these embodiments, solving the set of line-wise power balance equations (FT(W')) to obtain the problem vector (W').

[0043] In at least one of these embodiments, the desired line-wise solution vector is defined as ST= [O o o o [PD] [QD]Y, wherein: PD is a vector of active power demand in 5 each of the at least one transmission branch, and PG is a vector of active power generation in each of the at least one transmission branch.

[0044] In at least one of these embodiments, the objective function f0 (W') is a secondorder power generation cost function expressed as: NG f0 (W') = L c19 + c29 • PG9 + c39 • PGJ g=1 10 wherein NG is a total number of power generators in the power system, PG9 is the active power generation of a g-th power generator, c19 is a first constant for the g-th power generator, c29 is a second constant for the g-th power generator, and c39 is a third constant for the g-th power generator.

[0045] In at least one of these embodiments, the line-wise linear optimization formulation 15 comprises a linearized objective function, a linearized set of equality constraints, and a linearized set of inequality constraints.

[0046] In at least one of these embodiments, the linearized objective function is expressed as: NG L [(c29 + 2 • c39 • PG9 ) • .1PG9 ] g=l 20 wherein .1PG9 is a linear increment to the vector of active power generation.

[0047] In at least one of these embodiments, the linearized set of equality constraints is expressed as: I / D].1W' = '1ST wherein J is a first jacobian matrix, D is a second Jacobean matrix, and '1ST is an incremental 25 change to the line-wise solution vector, wherein '1ST = [O o o o o oy. - 10 - CA 3028644 2018-12-27

[0048] In at least one of these embodiments, the method comprises solving the linear optimization formula using a linear programming method.

[0049] In at least one of these embodiments, determining whether a set of stopping criteria is satisfied comprises determining whether a stopping condition is satisfied, wherein 5 the stopping conditioned is defined as: lfa(W'k) - fo(W'k-l)I < lfo(W'k)I - Y where: y is a pre-determining mismatch tolerance; f0 (W'k) is a first output of the objective function using values in the problem vector (W') determined in a current iteration (k) o; and fa (W'k-i) is a second output of the objective function using values in the problem vector (W') 1 O determined in a previous iteration (k-1).

[0050] In at least one of these embodiments, a set of stopping criteria is satisfied when a number of iterations of (b) to (d) is larger than a predetermined maximum number of iterations.

[0051] In at least one of these embodiments, a set of stopping criteria is satisfied when a 15 set of equality and inequality constraints is satisfied, the set of equality constraints being defined by FT(W)=ST.

[0052] In at least one of these embodiments, determining initial values for the control vector (WC), further comprises: reading a set of data from at least one sensor operatively coupled to the power system. 20

[0053] In at least one of these embodiments, the method comprises: based on the control vector (WC), determining whether the power system is in a state of voltage stability by performing an N-1 contingency analysis, and when the power system is not in a stage of voltage stability implementing a remedial action to stabilize the power system where the remedial action comprises at least one of: shedding a load bus of the power system, 25 increasing a power output at least one generator bus in the power system, re-activating at least one standby generator bus in the power system, and re-inserting standby transmission elements in the power system.

[0054] In accordance with another broad aspect of the teachings herein, there is provided at least one embodiment of a system for minimizing an objective function in respect of a - 11 - CA 3028644 2018-12-27 power system, the power system comprising at least one transmission branch having a first end coupled to a first bus and a second end coupled to a second bus, the system comprising: at least one server comprising at least one processor, wherein the at least one server is configured to: (a) determine initial values for a control vector (WC); (b) use the control vector 5 (WC) to solve a set of line-wise power balance equations (FT(W')) to obtain a problem vector (W'), wherein the problem vector (W') comprises the control vector (WC) and a dependent vector (WD); (c) use the problem vector (W') to solve a line-wise linear optimization formulation to determine an optimal incremental change to the control vector (LlWC); (d) update the control vector (WC) based on the optimal incremental change in the control vector 10 (LlWC); and (e) based on the updating of the control vector (WC), determine whether the objective function is minimized by determining when a set of stopping criteria is satisfied, wherein when the set of stopping criteria is not satisfied, the at least one processor is configured to iteratively repeat (b) to (d) until the set of stopping criteria is satisfied.

[0055] In at least one of these embodiments, the set of line-wise power balance equations 15 (FT(W)) is defined as: FT(W) = FF(W) FS(W) FA(W) FB(W) FP(W) FQ(W) wherein: FF(W) is a vector of line based voltage functions in relation to each first end of each of the at least one transmission branch; FS(W) is a vector of line-based voltage function in relation to each second end of each of the at least one transmission branch; FA(W) is a 20 vector of phase angle functions in relation to each first end of each of the at least one transmission branch; FB(W) is a vector of phase angle functions in relation to each second end of each of the at least one transmission branch; FP(W) is a vector of functions for buswise real power injections, and FQ(W) is a vector of functions for bus-wise reactive power injections. 25

[0056] In at least one of these embodiments, the problem vector (W') is defined as W' = [U, o, PF, PS, QF, QS, PG, QGF, and the dependent vector (WD) is defined as WD = [U,o,PF,PS,QF,QSF, and the control vector (WC) is defined as WC= [U,PGY, wherein: U - 12 - CA 3028644 2018-12-27 is a vector of squared voltages magnitudes for each bus of each of the at least one transmission branch; 8 is a vector of voltage angles for each bus of each of the at least one transmission branch; PF is a vector of active power flow in each first end of each of the at least one transmission branch; PS is a vector of reactive power flow in each second end of 5 each of the at least one transmission branch; QF is a vector of reactive power flow in each first end of each of the at least one transmission branch; QS is a vector of reactive power flowing in each second end of each of the at least one transmission branch; PG is a vector of active power generation in each of the at least one transmission branches, and QG is a vector of reactive power generation in each of the at least one transmission branches. 1 O

[0057] In at least one of these embodiments, solving the set of line-wise power balance equations (FT(W')) to obtain the problem vector (W').

[0058] In at least one of these embodiments, the desired line-wise solution vector is defined as ST= [O O O O [PD] [QD]Y, wherein: PD is a vector of active power demand in each of the at least one transmission branch, and PG is a vector of active power generation 15 in each of the at least one transmission branch.

[0059] In at least one of these embodiments, the objective function f0 (W') is a secondorder power generation cost function expressed as: NG lo (W') = L clg + c2g • PGg + c3g • PGJ g=1 wherein NG is a total number of power generators in the power system, PGg is the active 20 power generation of a g-th power generator, clg is a first constant for the g-th power generator, c2g is a second constant for the g-th power generator, and c3g is a third constant for the g-th power generator.

[0060] In at least one of these embodiments, the line-wise linear optimization formulation comprises a linearized objective function, a linearized set of equality constraints, and a 25 linearized set of inequality constraints.

[0061] as: In at least one of these embodiments, the linearized objective function is expressed - 13 - CA 3028644 2018-12-27 5 NG L [(c2g + 2 • c3g • PGg) • LiPGg] g==l wherein LiPGg is a linear increment to the vector of active power generation.

[0062] In at least one of these embodiments, the linearized set of equality constraints is expressed as: U D]LiW' = LiST wherein J is a first jacobian matrix, D is a second Jacobean matrix, and LiST is an incremental change to the line-wise solution vector, wherein LiST = [O O o o O of.

[0063] In at least one of these embodiments, the method comprises solving the linear optimization formula using a linear programming method. 10

[0064] In at least one of these embodiments, determining whether a set of stopping criteria is satisfied comprises determining whether a stopping condition is satisfied, wherein the stopping conditioned is defined as: lfo(W'k) - fo(W'k-l)I < lfo(W'k)I -Y where: y is a pre-determining mismatch tolerance; f0 (W'k) is a first output of the objective 15 function using values in the problem vector (W') determined in a current iteration (k); and f0 (W'k-i) is a second output of the objective function using values in the problem vector (W') determined in a previous iteration (k-1).

[0065] In at least one of these embodiments, a set of stopping criteria is satisfied when a number of iterations of (b) to (d) is larger than a predetermined maximum number of 20 iterations.

[0066] In at least one of these embodiments, a set of stopping criteria is satisfied when a set of equality and inequality constraints is satisfied, the set of equality constraints being defined by FT(W)=ST.

[0067] In at least one of these embodiments, determining initial values for the control 25 vector (WC), further comprises: reading a set of data from at least one sensor operatively coupled to the power system. - 14 - CA 3028644 2018-12-27

[0068] In at least one of these embodiments, the system comprises: based on the control vector (WC), determining whether the power system is in a state of voltage stability by performing an N-1 contingency analysis, and when the power system is not in a stage of voltage stability implementing a remedial action to stabilize the power system where the 5 remedial action comprises at least one of: shedding a load bus of the power system, increasing a power output at least one generator bus in the power system, re-activating at least one standby generator bus in the power system, and re-inserting standby transmission elements in the power system.

[0069] Other features and advantages of the present application will become apparent 10 from the following detailed description taken together with the accompanying drawings. It should be understood, however, that the detailed description and the specific examples, while indicating preferred embodiments of the application, are given by way of illustration only, since various changes and modifications within the spirit and scope of the application will become apparent to those skilled in the art from this detailed description. 15 BRIEF DESCRIPTION OF THE DRAWINGS

[0070] For a better understanding of the various embodiments described herein, and to show more clearly how these various embodiments may be carried into effect, reference will be made, by way of example, to the accompanying drawings which show at least one example embodiment, and which are now described. The drawings are not intended to limit 20 the scope of the teachings described herein.

[0071] FIG. 1 is a pi-model of an example transmission branch in a power system.

[0072] FIG. 2 is a process flow diagram showing an example embodiment of a process for performing a line-wise power flow analysis.

[0073] FIG. 3A is an example generalized Jacobian matrix generated by the process of 25 FIG. 2.

[0074] FIG. 3B is an example Jacobian matrix, generated by the process of FIG. 2, for a 6-bus power system.

[0075] FIG.4A is a plot comparing a solution set generated by the process of FIG. 2, as compared to a solution set generated using a conventional bus-wise power flow analysis. - 15 - CA 3028644 2018-12-27

[0076] FIG. 48 is a plot comparing example convergence characteristics for the process of FIG. 2 as applied to different power systems.

[0077] FIG. 4C is a table comparing the equations in a line-wise power flow formulation, as compared to a bus-wise power flow formulation. 5

[0078] FIG. 4D is a table comparing elements in line-wise Jacobian matrices, as compared to bus-wise Jacobian matrices.

[0079] FIG. 4E is a table comparing the size, condition, and sparsity factor of line-wise Jacobian matrices, as compared to bus-wise Jacobian matrices.

[0080] FIG. 5A is a plot showing the voltage collapse indices for an example 14-bus power 10 system with 20 transmission lines.

[0081] FIG. 58 is a plot showing the voltage collapse indices for an example 118-bus power system with 186 transmission lines.

[0082] FIG. 5C is a plot showing the voltage collapse indices for an example 583-bus power system with 965 transmission lines. 15

[0083] FIG. 5D is a table comparing a contingency analysis for an example 2835-bus system generating using a line-wise flow power analysis, as compared to a bus-wise power flow analysis.

[0084] FIG. 6 is an example process flow diagram for a process for solving a line-wise optimal power flow (LWOPF) problem using a successive linear programming (SLP) 20 algorithm.

[0085] FIG. 7 is a plot comparing solutions for an objective function generated using the process of FIG. 6, as compared to solutions generated using a bus-wise optimal power flow formulation.

[0086] FIG. 8 is a table comparing the results of solving the SLP algorithm, in FIG. 6, 25 using a set of line-wise power balance equations, as compared to solving the SLP algorithm without using the set of line-wise power balance equations.

[0087] FIG. 9A illustrates a simplified block diagram for an example embodiment of a power system that can perform power flow analysis, in accordance with the teachings here. - 16 - CA 3028644 2018-12-27

[0088] FIG. 9B is a simplified block diagram for an example embodiment of a server, in accordance with the teachings herein.

[0089] FIG. 1 0 is a process flow diagram for an example embodiment of a method for monitoring the state of a power system, in accordance with the teachings herein. 5

[0090] Further aspects and features of the example embodiments described herein will appear from the following description taken together with the accompanying drawings. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0091] Various embodiments in accordance with the teachings herein will be described below to provide an example of at least one embodiment of the claimed subject matter. No 1 0 embodiment described herein limits any claimed subject matter. The claimed subject matter is not limited to devices, systems or methods having all of the features of any one of the devices, systems or methods described below or to features common to multiple or all of the devices, systems or methods described herein. It is possible that there may be a device, system or method described herein that is not an embodiment of any claimed subject matter. 15 Any subject matter that is described herein that is not claimed in this document may be the subject matter of another protective instrument, for example, a continuing patent application, and the applicants, inventors or owners do not intend to abandon, disclaim or dedicate to the public any such subject matter by its disclosure in this document.

[0092] It will be appreciated that for simplicity and clarity of illustration, where considered 20 appropriate, reference numerals may be repeated among the figures to indicate corresponding or analogous elements or steps. In addition, numerous specific details are set forth in order to provide a thorough understanding of the example embodiments described herein. However, it will be understood by those of ordinary skill in the art that the embodiments described herein may be practiced without these specific details. In other 25 instances, well-known methods, procedures and components have not been described in detail so as not to obscure the embodiments described herein. Also, the description is not to be considered as limiting the scope of the example embodiments described herein.

[0093] It should also be noted that the terms "coupled" or "coupling" as used herein can have several different meanings depending in the context in which these terms are used. For 30 example, the terms coupled or coupling can have a mechanical, fluidic or electrical - 17 - CA 3028644 2018-12-27 connotation. For example, as used herein, the terms coupled or coupling can indicate that two elements or devices can be directly connected to one another or connected to one another through one or more intermediate elements or devices via an electrical or magnetic signal, electrical connection, an electrical element or a mechanical element depending on 5 the particular context. Furthermore, coupled electrical elements may send and / or receive data.

[0094] Unless the context requires otherwise, throughout the specification and claims which follow, the word "comprise" and variations thereof, such as, "comprises" and "comprising" are to be construed in an open, inclusive sense, that is, as "including, but not 10 limited to".

[0095] It should also be noted that, as used herein, the wording "and / or" is intended to represent an inclusive-or. That is, "X and / or Y" is intended to mean X or Y or both, for example. As a further example, "X, Y, and / or Z" is intended to mean X or Y or Z or any combination thereof. 15

[0096] It should be noted that terms of degree such as "substantially", "about" and "approximately" as used herein mean a reasonable amount of deviation of the modified term such that the end result is not significantly changed. These terms of degree may also be construed as including a deviation of the modified term, such as by 1 %, 2%, 5% or 10%, for example, if this deviation does not negate the meaning of the term it modifies. 20

[0097] Furthermore, the recitation of numerical ranges by endpoints herein includes all numbers and fractions subsumed within that range (e.g. 1 to 5 includes 1, 1.5, 2, 2.75, 3, 3.90, 4, and 5). It is also to be understood that all numbers and fractions thereof are presumed to be modified by the term "about" which means a variation of up to a certain amount of the number to which reference is being made if the end result is not significantly 25 changed, such as 1 %, 2%, 5%, or 10%, for example.

[0098] Reference throughout this specification to "one embodiment", "an embodiment", "at least one embodiment" or "some embodiments" means that one or more particular features, structures, or characteristics may be combined in any suitable manner in one or more embodiments, unless otherwise specified to be not combinable or to be alternative 30 options. - 18 - CA 3028644 2018-12-27

[0099] As used in this specification and the appended claims, the singular forms "a," "an," and "the" include plural referents unless the content clearly dictates otherwise. It should also be noted that the term "or" is generally employed in its broadest sense, that is, as meaning "and / or" unless the content clearly dictates otherwise. 5

[00100] The headings and Abstract of the Disclosure provided herein are for convenience only and do not interpret the scope or meaning of the embodiments.

[00101] Similarly, throughout this specification and the appended claims the term "communicative" as in "communicative pathway," "communicative coupling," and in variants such as "communicatively coupled," is generally used to refer to any engineered arrangement 1 O for transferring and / or exchanging information. Examples of communicative pathways include, but are not limited to, electrically conductive pathways (e.g., electrically conductive wires, electrically conductive traces), magnetic pathways (e.g., magnetic media), optical pathways (e.g., optical fiber), electromagnetically radiative pathways (e.g., radio waves), or any combination thereof. Examples of communicative couplings include, but are not limited 15 to, electrical couplings, magnetic couplings, optical couplings, radio couplings, or any combination thereof.

[00102] Throughout this specification and the appended claims, infinitive verb forms are often used. Examples include, without limitation: "to detect," "to provide," "to transmit," "to communicate," "to process," "to route," and the like. Unless the specific context requires 20 otherwise, such infinitive verb forms are used in an open, inclusive sense, that is as "to, at least, detect," to, at least, provide," "to, at least, transmit," and so on.

[00103] In addition it should be noted that the example embodiments of the systems and methods described herein may be implemented as a combination of hardware or software. In some cases, the example embodiments described herein may be implemented, at least in 25 part, by using one or more computer programs, executing on one or more programmable devices comprising at least one processing element, and a data storage element (including volatile memory, non-volatile memory, storage elements, or any combination thereof). These devices may also have at least one input device (e.g. a keyboard, mouse, touchscreen, or the like), and at least one output device (e.g. a display screen, a printer, a wireless radio, or 30 the like) depending on the nature of the device. - 19 - CA 3028644 2018-12-27

[00104] Electrical power systems are typically formed from a network of transmission branches (e.g., lines, cables, or transformers) that connect, and transfer power between various busses (or nodes) located in the power system. In various cases, busses in a power system may be grouped according to one of three categories: (a) generator busses (e.g., 5 generating stations); (b) load busses; and (c) slack busses (e.g., a reference generator bus).

[00105] In a power flow (PF) analysis of the power system, the power system is modelled using a set of power balance equations, and the equations are numerically solved to determine the current state of the power network. In various cases, for example, a PF analysis may be used to determine the voltage magnitudes and voltage phase angles at each 1 0 power system bus, the real power generation at the slack bus, reactive power generation at all generator busses and the real and reactive power flow through each transmission branch.

[00106] As stated in the background, PF analysis has become an important tool in many aspects of power system engineering, and constitutes the basis for studies in power system optimization, power system planning, and voltage stability analysis. An increasing, wide- 15 spread effort has accordingly been devoted in recent years to developing an efficient and accurate method for conducting a PF analysis.

[00107] In conventional PF analyses, a set of "bus-wise" (nodal) power balance equations is used to model a power system network. The bus-wise power balance equations rely on an admittance matrix, which relate the bus current injections to the voltages at each bus. In 20 various cases, the bus-wise power balance equations are solved using any appropriate iterative numerical technique. For example, the bus-wise power balance equations may be solved using a Newton-Raphson (NR) technique, which computes a Jacobian matrix of the matrix of bus-wise power balance equations. In at least some cases, the bus-wise Jacobian matrix provides indirect information about the health and voltage stability of the power 25 system. However, conventional bus-wise PF analyses suffer from a number of draw-backs.

[00108] First, in order to generate admittance matrices, advanced knowledge of the power system topology is required. This step adds to computational time required to complete Power Flow analysis.

[00109] Second, computing the bus-wise mismatch vector, in a Newton-Raphson 30 technique, is a computationally-intensive process. Bus-wise power balance equations - 20 - CA 3028644 2018-12-27 typically include a significant number of high-order non-linear terms (e.g., fourth-order and cubic terms), and many sinusoidal terms. Accordingly, determining the bus-wise mismatch vector is a computationally intensive process that results in slower performance for large power systems. 5

[00110] Third, computing the bus-wise Jacobian matrix, in a Newton-Raphson technique, is also a computationally-intensive process as it typically includes a significant number of high-order non-linear terms (e.g., fourth-order and cubic terms), all of which includes sinusoidal terms. Accordingly, the intensive computation required for computing the Jacobian matrix results in slower performance for large power systems. • 10

[00111] Fourth, the bus-wise Jacobian matrix does not, without further computation, provide direct information about power flow in transmission lines, or sets of transmission lines, which are susceptible to voltage collapse. As such, determining power flow in transmission lines or determining the set of critical transmission lines in a power system using a bus-wise formulation requires additional computation time. 15

[00112] In view of the foregoing, and in accordance with various teachings provided herein, there is provided a system and method for performing a line-wise power flow analysis using a set of line-wise power balance equations. In various embodiments, the line-wise power flow analysis may be solved using any suitable numerical technique, including a Newton-Raphson (NR) technique. 20

[00113] As explained in further detail herein, the set of line-wise power balance (LWPB) equations are formulated in terms of the square of the bus voltages, which results in lower order equations with lower-order terms, as compared to the bus-wise formulation. The set of LWPB equations also include far fewer sinusoidal terms than the bus-wise equations. In view of the foregoing, solving the LWPB equations is a less-computationally-intensive process. In 25 at least some embodiments, the line-wise power flow formulation is solved up to twice as fast as the bus-wise power flow formulation using an NR technique.

[00114] In various embodiments, as provided herein, the line-wise power flow (LPFW) formulation demonstrates stable numerical performance, monotonic convergence, high solution accuracy, and scales well for large power network systems. - 21 - CA 3028644 2018-12-27

[00115] In at least some embodiments, a line-wise Jacobian matrix is generated when solving the LPFW method using a Newton-Raphson technique. In various cases, the linewise Jacobian matrix may provide direct information, i.e., without further computation, for transmission lines, or sets of transmission lines, which are most susceptible to voltage 5 collapse. Accordingly, this allows for quick identification of critical transmission lines and for immediate corrective action to these transmission lines. In still yet other embodiments, the LPFW formulation may be used to quickly generate a contingency analysis for ranking the transmission lines in order of voltage collapse susceptibility.

[00116] In another aspect, in at least one example embodiment provided herein, the set of 10 line-wise power balance equations may be used as equality constraints in an optimal power flow (OPF) formulation for solving an optimal power flow (OPF) problem.

[00117] The goal of solving an optimal power flow (OPF) problem is to minimize an objective function, such as a total power generation cost function, a transmission real power loss minimization function, or a voltage profile improvement function. In particular, solving an 15 OPF problem has wide-spread applicability in power systems operations. For example, power system operators, who schedule the power demands of various loads in a system at frequent intervals, may routinely solve the OPF problem (i.e., as often as every 5 minutes) to determine the best operating schedule. Accordingly, it is desirable to have a method for solving the OPF problem that is highly efficient and requires minimal execution time (e.g., 20 especially for applications that require solving the OPF problem at short repeated time intervals).

[00118] The conventional method for solving OPF problems has relied on the set of buswise power balance equations to form the equality constraints in the OPF problem. As stated previously, however, the bus-wise formulation suffers from a number of draw-backs, 25 including high computational burden and slow execution times. The deficiencies of the buswise optimal power flow (BWOPF) formulation can accordingly yield non-optimal solutions, which may lead to economic loss (e.g., when attempting to minimize the total generation cost of a power system).

[00119] In accordance with various embodiments provided herein, there is accordingly 30 provided a line-wise optimal power flow (LWOPF) formulation which relies on the set of line- - 22 - CA 3028644 2018-12-27 wise power balance equations to form the equality constraints for solving the OPF problem. In at least some embodiments, explained in further detail herein, the LWOPF formulation may be solved using a successive linear programming (SLP) algorithm.

[00120] In various embodiments, owing to the lower orderterms in LWPF equations in 5 comparison to higher order terms in BWPF equations used for the equality constraints in an OPF problem, the use of the set of LWPF equations results in: (1) low-computational burden of the LWOPF formulation which results in the LWOPF method converging to an optimal solution in up to twice the speed as the conventional BWOPF formulation (due to using lower order terms in the computation), and (2) a lower order solution space (due to using lower 1 O order terms in the computation) whereby the LWOPF formulation consistently achieves a global optimal solution that equals or is better than that achieved by BWOPF.

[00121] Referring now to FIG. 9A, there is shown a simplified block diagram for an example embodiment of a power system 900A in which power flow analysis can be performed, in accordance with at least some embodiments described herein. The power flow 15 analysis system 900A provides the environment in which the systems and / or methods described herein generally operate.

[00122] As shown, the system 900A generally has a transmission network 902 (e.g., a power system). The transmission network 902, may be for example, formed from a network of transmission branches (e.g., lines, cables, or transformers) that connect, and transfer 20 power between various busses (or nodes) located in the network 902. In various cases, busses in the power network 902 may be either load busses or generator busses coupled to loads and generators, respectively. In various cases, the power network 902 may include one or more sensors 904a, 904b, 904c, 904d, and 904e. In various cases, sensors 904a - 904e may be remote units which are installed on various equipment located in the 25 transmission network 902, including power lines, transformers, cables, switching gears, loads, and generators. The sensors 904a - 904e may accordingly sense and measure various parameters of the transmission network 902 (e.g., to measure real-time, or near-real time operating data). In at least some embodiments, sensors 904a - 904e may include: (a) sensing elements, (b) analog-to-digital conversion systems for converting analog sensor 30 readings into a digital output for sensed data comprising corresponding digital readings; (c) - 23 - CA 3028644 2018-12-27 data processing systems for processing the sensed data; and (d) telecommunication systems to transmit the measured and / or processed data.

[00123] Data, or information, from sensors 904a - 904e may be transmitted to an energy control center 906 via one or more transmission lines. In various cases, the transmission 5 lines may be secure telemetry transmission lines. In other cases, the data from sensors 904a - 904e may be transmitted over secure wireless channels. The sensors 904 may be configured to transmit the measured parameters in real-time, near real-time, at predefined time-intervals, or at a pre-defined frequencies.

[00124] The energy control center 906 may include a supervision control and data 10 acquisition (SCADA) sub-center 908, a state estimation sub-center 910, an OPF solver subcenter 912, and a contingency analyzer sub-center 914, such as an N-1 contingency analyzer sub-center. In various cases, each of the sub-centers 908 - 914 may include one or more servers. For examples, additional servers may be used in each sub-center 908 - 914 to divide tasks to decrease computational time while additional servers may be used in each 15 sub-center 908 - 914 for redundancy purposes to improve the robustness of the energy control center 906.

[00125] As shown, the SCADA sub-center 908 is configured to receive data transmitted by the one or more sensors 904a - 904e. In at least some cases, the SCADA sub-center 908 is configured to process the data and present the data in a format accessible to a power system 20 operator. Accordingly, the power system operator may gain insight into the current state of the power system.

[00126] State estimation sub-center 910 may receive the data sets and measurements generated by the SCADA sub-center 908, and may be configured to detect and replace bad data elements (e.g., corrupted or erroneous data) with a revised best estimate. In at least 25 some cases, the state estimation sub-center 910 may also simply remove bad data elements. For example, methods for state estimation in power systems are described in: A. Monticelli, "Electric power system state estimation," in Proceedings of the IEEE, vol. 88, no. 2, pp. 262- 282, Feb. 2000.

[00127] The data set generated by the state estimation sub-center 91 O may then be 30 transmitted to the OPF solver sub-center 912. In accordance with the teachings provided - 24 - CA 3028644 2018-12-27 herein, the OPF solver sub-center 912 may use the data to solve a line-wise optimal power flow problem. The solution to the line-wise OPF problem may provide the power system operator with various information including, but not limited to, an optimal schedule for generating power in the power system. For examples, there may be several generators in 5 the power system and the optimal schedule for generating power dictates the amount of power that each generator can generate while not being overloaded while also being able to deal with any transient increases in load.

[00128] In at least some embodiments, the optimal schedule generated by solving the OPF problem may be then passed to the N-1 contingency analyzer 914. The N-1 contingency 10 analyzer 914 may be configured to perform a line-wise power flow analysis for each case where a transmission line experiences voltage collapse (e.g., "N' power flow analysis iterations for a power system having "N' transmission lines where a unique one of the "N" transmission lines undergoes voltage collapse). As explained herein, based on the contingency analysis, it may then be determined if the power system is in a state of voltage 15 instability. For example, if in an iteration of the contingency analysis, it is determined that the power system is in a state of voltage instability, the contingency analysis may be stopped, and one or more of a number of remedial actions may be taken, including: (a) re-optimizing the distribution of power across the transmission branches (e.g., solving an OPF problem as explained herein); (b) shedding load busses in the power system; (c) increasing the power 20 generation for power generators in the power system, or otherwise, activating new power generators; and (d) re-inserting standby transmission elements in the power system such as transformers. Once the remedial action is complete, the OPF solver sub-center 912 may then re-optimize the power system based on the change of configuration resulting from the remedial action, and the contingency analysis may be then repeated to determine if the power 25 system is in a state of voltage stability.

[00129] Referring now to FIG. 9B, there is shown a simplified diagram of an example server 900B, in accordance with at least some embodiments. The example server 900B may be used, for example, in any one of the sub-centers 908 - 914 of FIG. 9A.

[00130] As shown, the server 900B may generally include a processor 918a in 30 communication with a memory 918b and a communication interface 918c. - 25 - CA 3028644 2018-12-27

[00131] The processor 918a may be configured to execute a plurality of instructions to control and operate the various components of the server 9008. Processor 918a may also be configured to receive information from the various components of server 9008 and to make specific determinations using this information. The determinations may then be 5 transmitted to the memory device 918b and / or the communication interface 918c. In various embodiments, the processor 918a may be configured to perform the operations of any one of the sub-centers 908- 914 of FIG. 9A. In various cases, the processor 918a may be formed from an application-specific integrated circuit (ASIC), a digital signal processing (DSP) circuit, or any other suitable circuit. 10

[00132] Memory 918b may be, for example, a non-volatile read-write memory which stores computer-executable instructions and data, and a volatile memory (e.g., random access memory) that may be used as a working memory by processor 918a.

[00133] In various embodiments, communication interface 918c may be configured to receive data from one or more components of system 900A of FIG. 9A. In various cases, the 15 communication module 918c may, for example, comprise a wireless transmitter or transceiver and antenna. In other cases, the communication interface 918c may be configured for wired communication.

[00134] Referring now to FIG. 10, there is shown a process flow for an example embodiment of a method 1000 for monitoring the state of a power system. 20

[00135] At act 1002, the SCADA sub-center 908, of system 900A, may receive measurement data from one or more sensors 904a - 904e in the transmission network 902.

[00136] At act 1004, the SCADA sub-center 908 may be configured to process the data and present the data in a format accessible to a power system operator.

[00137] At act 1006, the state estimation sub-center 910 may receive the data sets and 25 measurements generated by the SCADA 908, and may be configured to detect and replace bad data elements (e.g., corrupted or erroneous data) with a revised best estimate. In at least some cases, the state estimation sub-center 910 may also simply remove bad data elements, as previously described. - 26 - CA 3028644 2018-12-27

[00138] At act 1008, the complete data set generated by the state estimation sub-center 91 O may then be transmitted to the OPF solver sub-center 912. In accordance with the teachings provided herein, the OPF solver sub-center 912 may use the data to solve a linewise optimal power flow (LWOPF) problem. The solution to the line-wise OPF problem may 5 be used by the power the system operator to determine an optimal schedule for generating power in the power system 900A.

[00139] At act 1010, the optimal schedule generated by solving the OPF problem may be then passed to the N-1 contingency analyzer 914. The N-1 contingency analyzer 914 may be configured to perform a line-wise power flow (PF) analysis, e.g., one line-wise PF analysis 1 O for each contingency of each of the "N' elements of a power system, in order to check the health of the power system (e.g., voltage stability) after each of those contingencies (e.g., "N' power flow analysis iterations for a power system having "N' elements such as generators, transmission lines, transformers, cables, etc.), as previously described.

[00140] After each iteration of the contingency analysis, it is determined, at act 1012a 15 whether the power system is in a state of voltage instability. If the power system is determined not to be in a state of voltage instability at act 1012a, it is then determined at act 1012b whether all iterations of the contingencies analysis are complete. If the contingency analysis is determined to be complete at act 1012b, then the method 1000 is completed at act 1014. Alternatively, if it is determined at act 1012b that the contingency analysis is not complete, 20 then method 1000 may return to act 101 O to complete the remaining iterations of the contingency analysis.

[00141] If the power system is determined not be voltage stable at act 1012a for a given iteration of the contingency analysis, the contingency analysis may be stopped and a remedial action may be taken at act 1016 to ameliorate the voltage instability, including one 25 or more of: (a) re-optimizing the distribution of power across transmission branches (e.g., solving an OPF problem as explained herein); (b) shedding load busses in the power system; (c) increasing the power generation for power generators in the power system, or otherwise, activating new power generators; and (d) re-inserting standby transmission elements in the power system such as transformers. Once the remedial action is complete, the method 1000 30 may return to act 1008 wherein the LWOPF problem is solved again (e.g., to re-optimize the - 27 - CA 3028644 2018-12-27 system in view of the change of configuration resulting from the remedial action), and the N- 1 contingency analysis at act 1010 is repeated using the new LWOPF solution.

[00142] Referring now to FIG. 1, there is shown a pi-model of an example transmission branch 100 in a power system. 5

[00143] As shown, the example transmission branch 100 couples a first bus 102 to a second bus 104. The first bus 102 and the second bus 104 may be characterized, for example, as load busses, generator busses, or slack busses. The transmission branch 100 may include a series impedance element 106. In the illustrated embodiment, the impedance of element 106 may be expressed according to Equation (1 ): ( 1) 10 where Z1..c.01 is the complex impedance of the element 106, R1 is the real component of the impedance, X1 is the imaginary component of the impedance, and the subscript" l" denotes the "l"-th transmission branch in a power system that includes a plurality of transmission branches. In a power system that includes "NT" transmission lines, each term in Equation (1) resolves into a vector which includes values over the range of l = l to NT. It should be noted 15 that the term "power system" may be used interchangeably with the term "power network".

[00144] Transmission branch 100 may also include shunt elements 108, 110 (YS1..c.{31) of the pi-model of the transmission branch 100. In various cases, the transmission branch 100 may be a line, a cable, or a transformer.

[00145] The power in the first end of transmission branch 100 (e.g., proximate bus 102) 20 may be characterized according to Equation (2): (2) where SF1..c.( / )F1 is the complex power in the first end, PF1 is the active power in the first end, and QF1 is the reactive power in the first end.

[00146] Similarly, the power in the second end of transmission branch 100 may be characterized according to Equation (3): (3) - 28 - CA 3028644 2018-12-27 where SSi.~J / JSi is the complex power in the second end, PSi is the active power in the second end, and QSi is the reactive power in the second end.

[00147] The complex voltage at the first bus 102 may be expressed by the term Va..e::Da, and the complex voltage in the second bus 104 may be expressed by the term VbLob. The power 5 flow in the first end of transmission branch 100, may then be accordingly determined according to Equation (4), based on the voltage difference between the first bus 102 and second bus 104. Equation (4) may be re-arranged to generate Equation (5). (4) (5)

[00148] The real and imaginary components of Equation (5) may be separated to generate Equations (6) and (7), respectively. (6) (7) 10

[00149] Equations (6) and (7) may then be divided to generate Equation (8), which can be re-arranged to generate Equation (9): (PSi • Rz + QSi ·Xi+ Vb 2 ) • tan(oa - ob) = PSi ·Xi+ QSi • Ri (PSi • Ri + QSi ·Xi+ Ub) • tan(oa - ob) - PSi • X1 + QSi • Ri = 0 (8) (9) where Ub expresses the square of the voltage magnitude at first bus 104 (Vb 2) of the transmission branch 100.

[00150] Equation (9) may herein express a "phase angle function" for the second end of 15 the transmission branch 100 (F B 1). When Equation (9) is applied to a power system having "transmission branches, a vector of phase angle functions is generated (e.g., v l = 1 to NT).

[00151] In order to remove sinusoidal terms, Equations (6) and (7) may be added to generate Equation (10): (10) - 29 - CA 3028644 2018-12-27 where Ua expresses the square of the voltage magnitude at the first bus 102 (V / ) of the transmission branch 100. Equation (10) may herein express cl "line-based voltage function" for the second end of the transmission branch 100 (FSi), or a vector of "line-based voltage functions" (e.g., '<:fl= 1 to NT). 5

[00152] The power flow in the second end of the transmission branch 100 may be similarly expressed in accordance with Equation (11 ). ( 11)

[00153] Equations ( 12) to ( 17), below, are generated from Equation ( 11) in an identical manner as Equations (5) to (10) were generated from Equation (4), respectively. (12) (13) (14) (15) (16) (17)

[00154] Equation ( 16) may herein express a "phase angle function" for the first end of 10 transmission branch 100 (FA 1), and Equation (17) may herein express a "line based voltage function" for the first end of the transmission branch 100 (FF1).

[00155] The transmission branch 100 may also be characterized according to the following bus-wise power equations, expressed in matrix form in Equations (18) and (19): [M] [;;] - U ·GS= PD - PG [M] [g~] + U ·BS= QD - QG (18) (19) where [M] is a modified bus incidence matrix, GS and BS express the bus conductance and 15 susceptance elements from loads and pi-models of transmission branches, PG and PD are - 30 - CA 3028644 2018-12-27 vectors of bus-wise active power generation and demand, respectively, and QC and QD are vectors of bus-wise reactive power generation and demand, respectively. Equation (18) may accordingly herein express a function for bus-wise real power injections (FP), while Equation (19) may herein express a function for bus-wise reactive power injections (FQ). 5

[00156] In various cases, the modified bus incidence matrix [M] is a matrix having dimensions NB x 2 • NT, wherein "NB" is the total number of busses in the power system (e.g., load and generator busses), and "NT" is the total number of transmission branches in the power system. The matrix [M] has a value of "1" for [M]a,z and [Mh,Nr+i, wherein "a" denotes the first bus 102 of the "I-th" transmission branch, and "b" denotes the second bus 1 0 104 of the " / -th"transmission branch, otherwise the values of the matrix [M] are equal to zero.

[00157] In view of the foregoing and in accordance with the teachings herein, a power system having a total number of transmission branches (NT) may be modeled using the set of line-wise power balance (LWBP) equations expressed in Equation (9), (10), (16), (17), (18), and (19), and combined in accordance with Equation (20): FF1(W) 0 FS1(W) 0 FA1(W) 0 (20) FB = 0 1(W) FPb(W) [PD - PG] FQc(W) [QD - QG] ('vl = 1 to NT; 'vb= l to NB - l; 'vc = 1 to NLB) 15 where the set of LWPB equations (FT(W)), on the left side of Equation (20), is a function of a problem vector W = [U o PF PS QF QS], and wherein the vector [0 o o o [PD - PG] [QD - QG]F defines the target solution vector ( ST) for the set of LWBP equations (FT(W)). Further, "NB" defines the number of busses in the power system (load and generator busses), and "NLB" defines the number of load busses in the power system. 20 Accordingly each of FF(W), FS(W), FA(W), and FB(W) is a vector which has "NT' equations, FP(W) is a vector that includes "NB-1"equations (e.g., all busses except the slack bus), and FQ(W) is a vector which includes "NLB" equations.

[00158] The number of equations in Equation (20) is expressed as (4 •NT+ NB - l + NLB), and the number of variables in Equation (20) is expressed as (4 •NT+ NB - l + - 31 - CA 3028644 2018-12-27 5 NLB). The variables in Equation (20) may be grouped into a set of known and unknown variables, which vary based on the bus type. Table 1, below, summarizes the known and unknown variables in respect of each of the generator, load, and slack busses. Table 1. Set of known and unknown variables for different bus-types Knowns Unknowns Slack Bus V PT= PG -PD 8 QT= QC - QD Generator Bus PT= PG -PD QT= QC - QD V 8 Load Bus PT= -PD V QT= -QD 8

[00159] As shown, the known variables in Equation (20) include: (a) the voltage magnitude and phase angle of the slack bus; (b) the real power injection and voltage magnitudes for all of the generator busses; and (c) the real and reactive power injection at the load busses. 1 0 Conversely, the unknown variables in Equation (20) include: (a) the real and reactive power injections at the slack buses; (b) the bus voltage angles and reactive power injections at the generator busses; and (c) the voltage magnitudes and phase angles at the load busses. The number of unknown variables in Equation (20) may be expressed as (4 ·NT+ NB - 1 + NLB). Table 2, below, summarizes the number of equations in Equation (20) which require 15 solving, and the number of unknown variables in each equation. Table 2. Number of equations requiring solving in Equation (20) Equations Number Unknowns Number FF(W) NT PF NT FS(W) NT PS NT FA(W) NT QF NT FB(W) NT QS NT FP(W) NB-1 8 NB-1 FQ(W) NLB V NLB - 32 - CA 3028644 2018-12-27

[00160] In conducting a power flow analysis, Equation (20) may be solved for the unknown variables. In particular, and in various cases, Equation (20) may be solved using an iterative numerical technique, such as a multi-variate Newton-Raphson (NR) technique. As explained herein, the solution to the power flow analysis may provide insights into the current state of 5 the power system, as well as any deficiencies in the power system (e.g., voltage collapse in one or more critical transmission lines). In various cases, the solution to Equation (20) identifies: (a) the square of the voltage magnitudes at each load bus; (b) the voltage phase angle at all busses, with the exception of the slack bus; and (c) the real and reactive powers at the first and second end of each transmission branch. 10

[00161] It will be appreciated from Equation (20), that the set of LWBP equations has been formulated in respect of the square of the bus voltage magnitudes (Ua, Ub), rather than the voltage magnitudes Wa2, V;). Accordingly, the LWBP equations are expressed using loworder terms (e.g., U5 rather than Va4). In particular, Equations (18) and (19) (FP(W), FQ(W)) are now expressed in linear terms, rather than quadratic terms. The LWBP equations also 15 include minimal sinusoidal terms, with the exception of Equations (18) and (19). This is in contrast to the conventional bus-wise formulation, which includes an abundance of non-linear equations with high-order terms (e.g., fourth-order and cubic), and many sinusoidal functions. As explained in further detail herein, the use of linear equations, low-order terms, and minimal sinusoidal functions, allows the power flow analysis to be solved more quickly using the set 20 of LWPB equations.

[00162] Referring now to FIG. 2, there is shown an example process flow for a method 200 for solving Equation (20) using an iterative Newton-Raphson technique. In various cases, the method 200 may be implemented using at least one processor 918a associated with the example server 900B of FIG. 9B. 25

[00163] At act 202, the current power system data is read. For example, in reference to Table 1, the known values for each bus in the power system may be read from the sensors 904a - 904e in system 900A of FIG. 9A. In at least some embodiments, the unknown values may be initialized to pre-determined initialization values (e.g., voltages at load busses are initialized to 1.0 per unit, and phase angles of all voltage phase angles are initialized to zero). 30 Based on the power system data readings, and the initialized values, an initialized problem - 33 - CA 3028644 2018-12-27 vector (W0) is generated based on known measured values for some variables and initialization values for other variables that have an unknown value.

[00164] At act 204, the power system is modeled using the set of LWPB equations (FT(W)) developed in the left-side of Equation (20). The set of LWBP is then solved using the values 5 determined for the problem vector (W). In the first iteration of method 200, the set of LWPB is solved using the initialized problem vector (W0).

[00165] At act 206, the difference between the target solution vector (ST), and the solution to the LWPB equations (FT(W)) (i.e., resolved in act 204), is determined. If the difference is below a pre-determined tolerance threshold, the method 200 proceeds to act 214. In various 10 cases, the tolerance threshold is substantially zero. When the method 200 proceeds to act 214, the method 200 is complete and the values for the problem vector (W) are determined to reflect the current state of the power system. In at least some embodiments, at act 214, it may be determined from the values of the problem vector (W) whether the power system is in a state of voltage stability (e.g., a transmission branch is susceptible to voltage collapse, 15 as explained in further detail herein). In other embodiments, the method 200 may also proceed to act 214 where a pre-defined maximum number of iterations of method 200 has been exceeded. In various cases, the maximum number of iterations may be pre-set by a power system operator, and may be in the range of three to seven iterations. In other cases, where the difference at act 206 is not below the pre-determined tolerance threshold, the 20 method 200 may continue to act 208.

[00166] At act 208, a Taylor's series expansion of the set of LWPB equations (FT(W)) is determined. In at least some embodiments, the Taylor's series expansion is a first-order derivative expansion, which is expressed according to Equation (21 ): (21) where [aF;~w)] is the Jacobian of the set of LWPB equations (FT(W)), and LlW is an 25 incremental change to the problem vector (W), wherein LlW = [LlU Llo LlPF LlPS LlQF LlQS]. Equation (21) may be accordingly re-arranged to solve for the incremental problem vector (LlW), in accordance with Equation (22). - 34 - CA 3028644 2018-12-27 LlW = [aF;;)r 1 • [ST - FT(W)] (22)

[00167] In order to solve Equation (22), the Jacobian of the set of the LWPB equations (FT(W)) must be first determined.

[00168] Referring now briefly to FIG. 3A, there is shown an example Jacobian matrix 300A for the set of LWPB equations (FT(W)) in Equation (20). 5

[00169] As shown, the Jacobian matrix 300A includes a number of sub-matrices, which represent different partial differentials with respect to different variables in the problem vector (W). The overall dimensions of the matrix 300A may be expressed as (4 ·NT+ NB - 1 + NLB) x (4 •NT+ NB - 1 + NLB), wherein "NT' is the total number of transmission lines in the power system, "NB" is the total number of power system busses, and "NLB" is the total 10 number of load busses.

[00170] In particular, the Jacobian matrix 300A includes a number of sub-matrices which are null as a result of the specific partial differential, including:[a;;], [aFF] [aFF] [aFs] [aFS] [aFF] [aFA] [aFA] [aFB] [aFB] [aFP] aPs ' aQs ' ao ' aPF ' aQF ' aPs ' aQs ' aPF ' aQF ' ao ' [aFP] [aFP] [aFQ] [aFQ] d aQF ' aQs ' ao ' aPF ' an [aFQ]. aPs 15

[00171] Referring now to FIG. 3B, there is shown an example Jacobian matrix 300B generated for a 6-bus power system. As shown, the Jacobian matrix 300B includes 40 rows and 40 columns and is otherwise highly sparse, with the non-zero elements being shown in black. In particular, the highly sparse nature of the line-wise Jacobian 300B results from the large number of null sub-matrices. In various cases, the high sparse nature of the line-wise 20 Jacobian allows the Jacobian to be computed quickly for large power system, while occupying minimal storage memory capacity on a processing system (e.g., memory 918b of server 900B in FIG. 9B).

[00172] Referring now back to FIG. 2, at act 210, Equation (22) is solved for the incremental change to the problem vector (LlW) using the inverse of the Jacobian matrix computed at act 25 208. - 35 - CA 3028644 2018-12-27

[00173] At act 212, using the incremental change to the problem vector (LlW) determined at act 210, an updated problem vector (W) is calculated in accordance with Equation (23). W = W+LlW (23)

[00174] In various cases, where a generator reactive power reaches a limit, the generator is set to be a load bus. In particular, if the problem vector (W) yields a reactive power output 5 for a generator bus that is higher than its capacity, the generator bus is switched to a load bus, wherein its real and reactive power output are expressed as negative values of a load and are solved to determine the voltage magnitude and voltage phase angle.

[00175] The method 200 may then return to act 206, and may iteratively repeat until either the tolerance condition is met, or otherwise, the maximum iterations is exceeded. 10

[00176] Table 3, below, compares the execution times and number of iterations required to solve the set of LWBP equations using method 200, as compared to solving a set of conventional bus-wise power balance (BWPB) equations using a Newton-Raphson (NR) technique. The data in Table 3 assumes a tolerance threshold of 10-5 per unit (i.e., at act 206 of method 200), and further assumes that loads are modeled as constant Mega Volt Ampere 15 (MVA) loads, and lines and transformers are modeled as pi-models considering their shunt admittance elements. Table 3. Execution times and number of iterations for solving a set of LWPB equations vs solving a set of BWPB equations using a Newton-Raphson technique Traditional Proposed method Execution time of bus-wise NR (Line-wise NR) line-wise method as Time No. of Time No. of percentage of bus- (milli- Iterations (milli- Iterations wise method seconds) seconds) 6-bus 4.0 3 4.0 4 100.00% 14-bus 9.0 2 4.0 3 44.44% 57-bus 22.0 3 8.0 4 36.36% 118-bus 31.0 3 15.0 3 48.39% 582-bus 332.0 6 152.0 6 45.78% 2383-bus 1579.0 4 936.0 4 59.28% 9241-bus 4618.9 6 2639.6 6 57.15% - 36 - CA 3028644 2018-12-27

[00177] As shown, the line-wise method 200 may execute up to twice as fast as the buswise method. In particular, for the larger power system which includes 9241-buses, the linewise method 200 executes in 1,979 less milliseconds relative to the bus-wise method (e.g., 57 .15% less time). It will be appreciated that, as the line-wise method 200 is scaled for larger 5 power systems, the time-savings is proportionately increased. In various cases, this may accordingly allow for quicker determinations of the state of the power system, and by extension, more immediate operation of the power system. As stated previously, the faster execution times may be attributed to formulating the LWBP equations, in Equation (20), using the square of the voltage magnitudes. This provides for lower-order terms and lower-order 10 equations, which require less intensive-computation, especially for generating the Jacobian matrix. Further, the fast execution times may be also attributed to the lack of sinusoidal terms, in Equation (20), which further simplifies the Jacobian computation.

[00178] Referring now to FIG. 4A, there is shown a plot 400A, which compares the bus voltage magnitude and angle solutions generated by the line-wise method 200 to those 15 generated by a conventional bus-wise NR method. The discrepancy in solutions, as between the two methods, is quantified as a Root Mean Square Error (RMSE) values. As shown, the RMSE values are consistently low, with the largest RMSE being for the 2383-bus system, where the solutions deviate by only 3.2E-07 for the voltage magnitudes, and 2.6E-07 for the voltage angles. Accordingly, plot 400 demonstrates that the line-wise method 200 generates 20 solutions that are at least as accurate as the bus-wise method.

[00179] Referring now to FIG. 4B, there is shown a plot 400b which compares the convergence properties of the line-wise method 200 across different power systems with varying numbers of busses. As shown, in each case, the line-wise method 200 achieves accurate solutions and demonstrates consistent monotonic convergence for each power 25 system, as the maximum mismatch progressively reduces in each iteration in this plot.

[00180] Referring now to FIG. 4C, there is shown a table 400C which compares the equations used in the line-wise formulation of Equation (20), and the equations used in a bus-wise formulation.

[00181] As shown, the line-wise power equations are limited to including 1st and 2nd order 30 terms. In contrast, the bus-wise equations include an abundance of higher-order terms (e.g., - 37 - CA 3028644 2018-12-27 3rd and 4th order terms), which increase in quantity for larger power systems. As stated previously, the line-wise equations include fewer higher-order terms because the equations are formulated according to the square of the bus voltage magnitudes (Ua, Ub), rather than the voltage magnitudes (Va2 , Vt). Accordingly, this results in Equations (18) and (19) (FP(W), 5 FQ (W)) being expressed in linear terms. In contrast, the bus-wise equations are formulated using voltage magnitudes, which may result in fourth-order terms (e.g., Va4 versus UJ). The lack of high-order terms in the line-wise formulation allows the line-wise Jacobian to be computed more quickly, and accordingly, allows for faster execution times for solving the power flow problem (see e.g., Table 3). 1 O

[00182] Table 400C also provides a comparison of the number of products in the line-wise equations including sine and / or cosine terms, as compared to the bus-wise equations. As shown, the line-wise equations generally include less products with sinusoidal terms than the bus-wise equations (e.g., only approximately 15% of the total products in the line-wise equations include sine or cosine terms). This is attributable to the fact that Equations (18) 15 (FB(W)) and Equation (19) (FA(W)) are the only sinusoidal functions in the LWBP equations.

[00183] Referring now to FIG. 4D, there is shown table 400D which compares the number of fourth-order, cubic, and quadratic terms in the line-wise Jacobian matrices, as compared to the bus-wise Jacobian matrices. Table 400D also compares the number of linear terms 20 and the number of products in the line-wise and bus-wise Jacobian matrices. As shown, the line-wise Jacobian matrices consistently include a larger number of linear terms, fewer high order terms (e.g., fourth-order, and cubic), and fewer product computations, as compared to the bus-wise Jacobian matrices. In particular, and on average, the line-wise Jacobian matrices includes 1.5 times less product computations than the bus-wise Jacobians. The 25 abundance of linear and low-order terms, and the lack of non-linear and higher-order terms, allows the line-wise Jacobians to be computed more quickly than the bus-wise Jacobians. This reduces the execution time for the line-wise power flow analysis.

[00184] Referring now to FIG. 4E, there is shown table 400E, which compares the size, condition, and sparsity factors for the line-wise Jacobian matrices as compared to the bus- 30 wise Jacobian matrices. As shown, while the line-wise Jacobian matrices may be larger in - 38 - CA 3028644 2018-12-27 size, this is mitigated by the line-wise Jacobians exhibiting better conditioning and better sparsity factors. For example, If the condition number is very large, then the matrix is said to be ill-conditioned so it is desired that the condition number be small, which is possible using the teachings herein. In particular, for large power systems (e.g., 9,241 busses), the line- 5 wise Jacobian has a 99.99 sparsity factor, which allows the Jacobian to be computed relatively quickly while occupying minimal memory storage space on an executing processing system (e.g., memory 918b in FIG. 9B). In various cases, the reduced memory usage also allows for faster computation of the Jacobian.

[00185] In various some embodiments, based on the Jacobian matrix generated at act 208 10 of method 200, it may be determined which transmission lines, or set of transmission lines, in the power system is most susceptible to voltage collapse (e.g., at act 214 of method 200).

[00186] In particular, it has been appreciated that the elements of Jacobian sub-matrices, corresponding to the partial-differentials [aaFFi], and [aaF51 ] corresponding to the "I-th" Ua Ua transmission line, reduce to substantially zero (or proximate to zero) when the "I-th" 15 transmission line is susceptible to voltage collapse.

[00187] The sub-matrix [aaFF1] is a partial-differential of the line-based voltage function for Ua the "first end" of the "I-th" transmission branch with respect to the square of the voltage magnitude at a first bus of the transmission branch (e.g., bus 102 in FIG. 1). Conversely, sub-matrix [aaF51] is a partial-differential of the line-based voltage function for the "second Ua 20 end" of the "I-th" transmission branch with respect to the square of the voltage magnitude at a second bus of the transmission branch (e.g., bus 104 in FIG. 1 ). Accordingly, a first voltage collapse index (VCia) may be formulated using the sub-matrix [aFF 1 ] to detect voltage collapse aua in the first end of a transmission branch, and a second voltage collapse index (VC!b) may be formulated using the second sub-matrix [aF51 ] to determine voltage collapse in the second aub 25 end of a transmission branch, wherein each of the first and second VCls may be expressed in accordance with Equations (24) and (25). (24) - 39 - CA 3028644 2018-12-27 (25)

[00188] In particular, it is has been appreciated that Equations (24) and (25) yield voltage collapse indices as a result of information available from Equations (1 0) and (17), respectively (FS,, FF,). Equations (10) and (17) are quadratic equations, which can be solved for the square of the bus voltages ( Ua, Ub) according to Equations (26) and (27), respectively: (26) (27) 5

[00189] In normal power system operations, each of Equations (26) and (27) resolves into two solutions for the square of the bus voltage ( Ua, Ub), e.g., one solution which is physically feasible (e.g., a higher voltage solution), and one solution which is physical infeasible (e.g., a lower voltage solution). At voltage collapse - i.e., as the power system is loaded - the two voltage solutions for each of Equations (26) and (27) converge until each Equation resolves 10 into one solution, rather than two solutions. Accordingly, the term located inside of the square root portion effectively reduces to zero to yield one solution, and Equations (26) and (27) maybe re-written and re-arranged according to Equations (28) and (29), respectively: (28) (29)

[00190] The terms in Equations (28) and (29) are located identically within Equations (24) and (25). Accordingly, at voltage collapse, the two sub-differentials in Equations (24) and 15 (25) reduce to substantially zero, as shown in Equations (28) and (29). It should be noted that since VCI values are not similar for all power systems and their behavior is nonlinear, their proximity to zero is monitored to determine if there is a Voltage Collapse.

[00191] In view of the foregoing, the line-wise power flow method allows voltage collapse to be determined, at both ends of a transmission line (e.g., VCia, VC!b), directly from the - 40 - CA 3028644 2018-12-27 Jacobian matrix via the first and second VCls. In contrast, the Jacobian generated by the bus-wise NR method requires additional computation to resolve the transmission lines susceptible to voltage collapse. The additional computation accordingly requires additional and delayed computation time. 5

[00192] Referring now to FIG. 5A, there is shown a plot 500A which plots the voltage collapse indices (VCia, VC!b) for an example 14-bus IEEE power system which includes 20 transmission lines. In this example case, the bus loads in the power system are scaled by 3.83 times in order to simulate higher loading conditions to induce voltage collapse. The indices, used in the plot 500A, are extracted directly from the line-wise Jacobian matrix 10 generated for that system. As shown, in this example case, the fifth line (L5) is most susceptible to voltage collapse because the second voltage collapse index (VC / b) is zero, and its first voltage collapse index is (VCia) is significantly low (e.g., 4.?0E-01 ).

[00193] Referring now to FIG. 5B, there is shown a plot 500B which plots the voltage collapse indices (VCia, VC!b) for an example 118-bus IEEE system which includes 186 15 transmission lines. In this example case, the bus loads are scaled by 4.23 times to induce voltage collapse. As shown, transmission line 108 (L 108) is most susceptible to voltage collapse as the first voltage collapse index (VCia) is equal to 2.35E-03, and the second voltage collapse (VC / b) index is equal to 2.30E-02 (i.e., substantially zero).

[00194] Referring now to FIG. 5C, there is shown a plot 500C which plots the voltage 20 collapse indices (VCia, VC!b) for an example 582-bus real power system which includes 965 transmission lines. In this example case, the bus loads are scaled by 2.253 times to induce voltage collapse. As shown, transmission line 109 (L 109) is most susceptible to voltage collapse as the first voltage collapse index (VCia) is equal to 6.37E-03, and the second voltage collapse (VC / b) index is equal to 2.53E-01. 25

[00195] Table 4, below, shows example worst-case voltage collapse indices for a 14-bus, 118-bus, and 582-bus system, and the computation-time required to determine these indices using a bus-wise NR method versus a line-wise NR method. - 41 - CA 3028644 2018-12-27 Table 4. Determining voltage collapse indices using bus-wise method versus linewise method Worst VCI Computation time (minimum VCI at both sides) (milliseconds) Bus-wise Line-wise Difference Bus-wise NR Line-wise NR NR NR method method method method 14-bus VCla 0.277 0.274 0.003 13.0 4.0 VClb 0.0 0.0 0.0 118-bus VCla 0.00195 0.00235 0.0004 39.0 15.0 VClb 0.029 0.023 0.006 582-bus VCla 0.006 0.0064 0.0004 421.0 152.0 VClb 0.2512 0.2499 0.0013

[00196] As shown, the bus-wise and line-wise methods determine nearly the same values 5 for the VCls, which at least demonstrates the solution accuracy for the line-wise method. The line-wise method, however, determines these indices in less computational time. On average, the line-wise method resolves the VCls as fast, or up to twice as fast, as the buswise method. The faster computation time, for the line-wise method, results because the VCls are directly extracted from the line-wise Jacobian matrix. In contrast, additional 10 computation is required to determine the VCls from the bus-wise Jacobian. The faster computation, under the line-wise method, may accordingly allow quicker determinations of the critical lines in a power system, and by extension, allow for more immediate corrective action.

[00197] In at least some cases, the VCls derived from line-wise Jacobian may be used for 15 performing an N-1 contingency analysis. In a power system having "N' elements (e.g., transmission branches, cables, etc.) the contingency analysis performs a power flow analysis for "N" cases where each one of the "N" elements uniquely fails separately. The objective is to determine whether the remaining "N-1" elements in the power system allow the power system to function properly, and to supply power to all of the loads, for each of the "N" voltage 20 collapse cases. The contingency analysis may be performed by first generating a ranking of the transmission lines in order of voltage collapse susceptibility. Subsequently, based on the - 42 - CA 3028644 2018-12-27 ranking, a separate power flow analysis is conducted for each transmission line (e.g., starting from the transmission line most susceptible to voltage collapse) to determine the effect of that transmission line collapsing. It will be appreciated that for a power network that includes thousands of elements, thousands of contingency analysis are required to be completed. As 5 such, using a line-wise power flow analysis that quickly and immediately resolves the VCI indices, may allow for significant time savings. Table 5, below, demonstrates the computation time to generate a voltage stability contingency analysis using a line-wise NR method versus a bus-wise NR method. Table 5 demonstrates that the line-wise NR method is able to generate the contingency ranking in as fast, or twice as fast, as the conventional bus-wise 10 NR method. Table 5. Comparison of computation time for voltage stability contingency analysis for line-wise method versus bus-wise method Number of Traditional Proposed Computation time contingency bus-wise Line-wise of line-wise System cases NR method NR method method as (number of (seconds) (seconds) percentage of lines) bus-wise method 14-bus 20 0.26 0.08 30.77% 57-bus 80 2.24 0.988 44.11% 118-bus 186 7.254 4.185 57.69% 300-bus 411 49.32 22.194 45.00% 2383-bus 2896 6081.8 2710.7 44.57% 15

[00198] Referring now to FIG. 50, there is shown table 5000, which shows an example contingency analysis for a 2835-bus Polish system using a bus-wise NR method, as compared to a line-wise NR method. The table 5000 shows the twenty lines that are the most susceptible to voltage collapse, out of the 2896 lines assessed as shown in Table 5. As shown, each method generates the same contingency analysis ranking and analysis, 20 however, the line-wise NR method generates the ranking in as equal, or less, amount of execution time as the bus-wise NR method.

[00199] In various cases, where the contingency analysis determines a VCI value is critical, at least one of a number of corrective actions may be taken, including: (a) re-optimizing the - 43 - CA 3028644 2018-12-27 distribution of power across transmission branches (e.g., solving an OPF problem as explained herein); (b) shedding load busses in the power system; (c) increasing the power generation for power generators in the power system, or otherwise, activating new power generators; and (d) re-inserting off-line elements (such as transformers, transmission lines, 5 cables, generators, etc.) in the power system.

[00200] In at least one example application, the set of line-wise power balance equations may be used as equality constraints in an optimal power flow (OPF) formulation for solving an optimal power flow (OPF) problem.

[00201] The objective of an OPF problem is to minimize an objective function in relation to 10 a power system. In at least some embodiments, the objective function is a total power generation cost function (f0(W')) for the power system which is approximated as a secondorder cost function, according to Equation (30): (30) where NG is a total number of power generators in the power system, g is the gth power generator, PG9 is the active power generation of the gth power generator, c19 is a first 15 constant for the gth power generator, c29 is a second constant for the gth power generator, and c39 is a third constant for the gth power generator. Further, W' is a modified problem vector.

[00202] In other embodiments, the objective function may be any other non-linear continuously differential objective function, depending on the minimization problem being 20 studied. For example, rather than being an actual financial cost or it may be another parameter such as input power that is to be minimized depending on the constants which are selected for the function. These include real power loss minimization, determination of minimum amount of capacitors for a power system to ensure satisfactory voltage profile, etc.

[00203] The following example, shown in Table 6, pertains to Ontario, Canada. The Ontario 25 electric power system has 5.6 million customers with a peak demand of 25 GW. The method provides a 2% reduction of by finding a better optimal solution. Considering a 2% improvement in the optimal solution determined by OPF every hour of the year, for 8760 hours, the cost savings will equal: $333,274,873.36 for the 5.6 million customers. - 44 - CA 3028644 2018-12-27 Table 6. Savings using LWOPF formulation for Example Power System in Ontario, Canada Item Power & Distribution Revenues ($) Total Electricity Generated and Delivered to Customers (kWh) Amount Per Year for a 2% Saving Amount $16,663,743,668.00(1) 126,378, 180,878(1) $333,274,873.36 (e.g., 2% of $16,663,743,668.00) $38,045.08 Amount Per Hour for a 2% Saving (e.g., 118760th of $333.274,873.36) (1) Data from page #3 of the 2017 Yearbook of Electricity Distributors. Please refer to the file 5 PDF: 2017 Yearbook_of_Electricity_Distributors.

[00204] In the optimal power flow problem, the minimization of the objective f0 (W') is subject to a set of equality constraints, expressed by Equation 31 ), and a set of inequality constraints, expressed by Equation (32): [g(W') = d h :S fh (W') :S h (31) (32) 10 where [g (W') is the set of equality constraints, d is a target vector for the equality constraints, fh (W') is a set of inequality constraints, and h:. and h are the upper and lower bounds of the inequality constraints, respectively.

[00205] In accordance with the teachings provided herein, a line-wise optimal power flow formulation (LWOPF) is provided for solving the OPF problem, wherein the equality 15 constraint [g (W') is defined by the set of LWPB equations FT(W) developed in Equation (20), with the exception that Equations (18) and (19) are re-arranged according to Equations (33) and (34), respectively. - 45 - CA 3028644 2018-12-27 FP1 =PG+ [M] [;~] - U ·GS= PD FQ1 = QG + [M] [3~] + U ·BS= QD (33) (34)

[00206] The problem vector W is accordingly modified to generate a modified problem vector W' = [U 8 PF PS QF QS PG QGf. Similarly, the target vector ST is modified to generate a modified target vector d, werein d = [O O O O PD QDf.

[00207] In at least some embodiments, the modified problem vector (W') may be divided 5 into two components: a control vector (WC), and a dependent vector (WD), wherein the control vector WC= [U PGf, and the dependent vector WD = [8 PF PS QF QS QGf.

[00208] By contrast, Equation (32) represents the set of inequality constraints for solving the optimal power flow problem. Equation (32) expresses the various system limits for the bus voltages, the generator outputs, and the limits on the MVA line power flow in the power 10 system. In various embodiments, the inequality constraints in Equation (32) are expressed according to Equations (35) - (39): ui s ui s ui ('t / i E NB) (35) PG9 s PG9 s PG9 ('t / gENG) (36) QG9 s QG9 s QG9 ('t / g ENG) (37) PF / + QF / 5 (SLz(A)) 2 ('t / l ENT) (38) PS / + QS / s (SL1ca)) 2 ('t / l ENT) (39) where Equation (35) expresses the upper and lower limits for the square of the bus voltage (Ud at the i-th index bus, as applied to all busses in the power system (NB); Equation (36) expresses the limits for the active power generation (PG9) at the g-th power generator, as 15 applied for all power generators in the power system (NG); Equation (37) expresses the limits for the reactive power generation (PG9) at the g-th power generator, as applied for all power generators (NG); and Equations (38) and (39) express the limits on the active and reactive - 46 - CA 3028644 2018-12-27 MVA line power flow (SL1), respectively, at the I-th transmission line, as applied for all transmission lines (ND.

[00209] Accordingly, the upper and lower limits for the inequality constraints in Equation (32) are defined according to Equations (40) and (41). I]_ = [U ff ~ 0 O] I]_= [U PG QC (SL1(A)) 2 (SL1cs)) 2 ] (40) (41) 5 [0021 0] Using the equality and inequality constraints, the objective function is then minimized using an LWOPF method. In various cases, the LWOPF method is solved using a successive linear programming (SLP) algorithm, which involves: (a) solving the set of linewise power balance equations (e.g., f9 (W')) to resolve the current state of the power system, (b) solving a linear incremental OPF model to determine an optimal incremental change to 10 the control vector (WC); and ( c) updating the control vector using the determined optimal incremental change. The SLP algorithm is iterated until the control vector (WC) minimizes the objective function (i.e., the cost generation function or another function as the case may be). While an SLP algorithm is provided herein for solving the OPF problem, it will be appreciated that any other appropriate iterative or non-iterative algorithm may be alternatively 15 used.

[00211] Referring now to FIG. 6, there is shown a process flow for an example embodiment of a method 600 for solving the line-wise optimal power flow problem (LWOPF) using a successive linear programming (SLP) algorithm. In various cases, the method 600 may be implemented by the OPF solver sub-center 912 of FIG. 9A. 20

[00212] At act 602, power system data is read to determine the current state for the control variables (WC). For example, the current power system data may be retrieved from data provided by one or more sensors 904a - 904e of FIG. 9A, as processed via the SCADA subcenter 908 and the State Estimator sub-center 910 of FIG. 9A.

[00213] At act 604, the set of line-wise power flow equations (e.g., f9 (W')) is solved to 25 determine the dependent variables (WD). In various cases, act 604 requires implementing - 47 - CA 3028644 2018-12-27 the method 200 of FIG. 2 using the modified line-wise power flow equations (19 (W')) in Equation (31 ).

[00214] At act 606, it is determined whether the stopping criteria has been satisfied. In at least some embodiments, determining whether the stopping criteria is satisfied includes 5 determining whether: (a) a pre-defined maximum number of iterations, of method 600, is exceeded; (b) the equality constraint of Equation (31) is satisfied; (c) the inequality constraint of Equation (32) is satisfied; and / or (d) the stopping condition, expressed in Equation (42), is satisfied: (42) where lo (W'k) is the solution to the objective function using the values in the problem vector 10 (W') resolved in the current iteration (k) of method 600, and lo (W'k-i) is the solution to the objective function using the values in the problem vector (W') resolved in the previous iteration (k-1) of method 600. Further, y is a predefined mismatch tolerance.

[00215] In at least some embodiments, the stopping criteria is determined by first considering whether the equality and inequality constraints of Equations (31) and (32) are 15 satisfied. If the constraints are satisfied, it is then determined whether the stopping condition of Equation (42) is satisfied. In other embodiments, the stopping criteria is satisfied based only on whether the pre-defined maximum number of iterations for method 600 is exceeded.

[00216] Where the stopping criteria at act 606 is satisfied, at act 612, the method 600 is determined to be complete, and the control vector (WC) is determined to minimize the 20 objective function (f0 (W')). Otherwise, where the stopping criteria at act 606 is not satisfied, the method 600 continues to act 608.

[00217] At act 608, the objective function and set of equality and inequality constraints are linearized and solved for an optimal incremental change to the problem vector (!1W'). In various embodiments, the optimal incremental change to the problem vector 11W' is defined 25 as 11W' = [!1U 118 !1PF !1PS !1QF !1QS !1PG !1QGY. - 48 - CA 3028644 2018-12-27

[00218] In at least some embodiments, the objective function is linearized by approximating the objective function using a first-order Taylor's series expansion, in accordance with Equations (43) and (44). minlif0 (W) = [dt;:)] · liW (43) NG min lif0 (liW) = L [(c2g + 2 • c3g • PGg) • liPGg] (44) g=l •

[00219] The linearized set of equality constraints may be similarly expressed in accordance 5 with Equation (45), which may be rearranged to solve for the optimal incremental change to the problem vector (liW'): li[g(W') = U D]liW' = lid (45) where J is the Jacobian of the set of line-wise power balance equations, as shown in 300A of FIG. 3, and D is a submatrix of partial differential equations in respect of the real power generation (PG) and the reactive power generation (QG), which are not included in the original 10 problem vector (W). The sub-matrix O is expressed according to Equation (46). r -aFF -BFS -BFA -BFP -BFQ1T D _ BPG BPG BPG BPG BPG - BFF BFS BFA BFP BFQ - - - - - BQG BQG BQG BQG BQG (46) The vector lid is further defined as lid = [O O O O O ot.

[00220] The linearized set of inequality constraints may be similarly expressed in accordance with Equations (47) - (51 ): (47) wherein liUi =a· (Ui - U?) and liUi =a· (Ui - U?) (48) wherein liPGg = a· (PGg - PGJ) and liPG9 =a· (PGg - PGJ) - 49 - CA 3028644 2018-12-27 wherein !).QC 9 = a· (QC 9 - QC$) and !).QC 9 = a· (QC 9 - QC$) 2 • PF1 • !).PF1 + 2 • QF1 • !).QF1 ~ !).SL1 (Vl ENT) 2 • PS1 • !).PS1 + 2 • QS1 • !).QS1 ~ !).SL1 (Vl ENT) wherein !).SL1 = a· (SL1 - SLD (49) (50) (51) where a is a pre-defined incremental step-size, and U?, PCJ, QC$, SLf are the values located in the modified problem vector (W') resolved in act 604.

[00221] Accordingly, Equations (44) - (51) define the linearized optimization formulation for the LWOPF problem. Using the values derived in act 604 for the problem vector (W'), the 5 linearized optimization formulation may be solved for the optimal incremental change to the problem vector (!).W') in Equation (45). In various cases, the linearized optimization formulation may be solved using any appropriate linear optimization method (e.g., simplex, or interior point algorithm).

[00222] At 610, the problem vector (W') is updated based on the incremental change to 10 the problem vector (!).W') determined at act 608, in accordance with Equation (52): W'(k) = W'(k-1) + !).W'(k) (52) where w1Ck) is the updated problem vector, W'(k-i) is the problem vector generated in the previous iteration (k-1) of method 600, and !).W'(k) is the incremental change to the problem vector, as determined at act 608.

[00223] Based on Equation (52), an updated control vector (WC) is extracted from the 15 updated problem vector (W'(k)), and the method 600 may return to act 604 to proceed with the next iteration. The method 600 iterates until the stopping criteria at act 606 is satisfied.

[00224] Table 7, below, compares the execution time and number of iterations, required to minimize the objective function using the line-wise method 600, as compared to a bus-wise - 50 - CA 3028644 2018-12-27 optimal power flow method. The comparison is shown for three different OPF solvers implementing both the method 600, and the bus-wise OPF me.thod. Table 7. Comparison of LWOPF method and a conventional bus-wise OPF method First OPF Solver Second OPF Solver Third OPF Solver System Time It Objective Time It Objective Time It Objective ratio #e r f unc t·1 0n ra t·1 0 ratio #e r f unc t·1 0n ra t·1 0 ratio #e r f unc t·1 0n ra t·1 0 (%) (%) (%) (%) (%) (%) 6-bus 40.6 2 98.56 47.8 2 99.99 53.6 2 98.56 14-bus 58.8 2 98.64 70.0 2 99.98 82.5 2 99.89 30-bus 49.0 2 98.30 49.0 2 99.84 52.0 2 99.62 57-bus 84.8 3 99.72 133.3 3 100.00 187.9 3 100.0 118-bus 42.5 4 97.58 45.7 4 99.81 53.3 4 97.14 300-bus 76.1 3 95.90 94.0 3 98.05 98.3 3 98.12 582-bus 82.1 4 97.77 85.3 4 100.00 85.3 4 97.95 2383-bus 56.4 3 99.46 65.1 3 99.98 65.1 3 99.48 9241-bus 62.3 6 97.52 111.5 6 99.84 154.5 6 99.86 5

[00225] As shown, for various systems with different numbers of buses, the LWOPF method 600 is, on average, nearly twice as fast as the conventional bus-wise OPF method as observed from the "Time Ratio(%)" column. Further, the solution to the objective function, as between the two methods, is nearly identical as observed from the "Objective function ratio(%)" column, which at least demonstrates the accuracy of the LWOPF method 600. As 1 O explained previously, the faster execution times for the LWOPF are attributed to the use of the square bus voltage magnitudes ( U), rather than the voltage magnitudes (V 2 ). This fact results in the line-wise equations being expressed with lower order terms (see e.g., Equations (33) and (34)). Further, the line-wise formulation also includes lesser sinusoidal components. Accordingly, this simplifies the calculation of the Jacobian matrix during acts 604 and 608 of 15 method 600, and decreases the computation and execution time of the L WOPF method 600.

[00226] Referring now briefly to FIG. 7, there is shown a plot 700 which compares the solution to the objective function using the LWOPF method 600 versus a solution generated using a leading commercial bus-wise OPF solver. In particular, the Y-axis expresses the deviation to the solution of the objective function as between the two methods, expressed as 20 a percentage of the best-known value generated by the leading commercial bus-wise OPF solver. The plot 700 assumes a mismatch tolerance (y) of 10-6 . As shown, the LWOPF method 600 provides monotonic convergence, as the deviation of the objective function from - 51 - CA 3028644 2018-12-27 the best-known bus-wise solution reduces progressively, and at up to twice the speed of the leading commercial bus-wise OPF solver.

[00227] Table 8, below, provides a comparison between the solution to the OPF problem using the line-wise set of equations, in accordance with method 600, as compared to solving 5 the OPF problem using a set of bus-wise equations. Table 8 assumes the objective function is a total power generation cost function, as expressed in Equation (30). Table 8. Comparison of cost generation and execution time between LWOPF method and bus-wise method solved using SLP ( / ) Objective: Execution Time Objective Ratio of Q) Generation Cost ($) (seconds) function execution ( / ) ( / ) :::, ratio times System m..... SLP SLP SLP LWOPF as LWOPF as 0 SLP Based Based lter lter % of % of ci Based LWOPF BWOP # Based # BWOPF BWOPF z BWOPF F LWOPF 6-bus 6 3,142.0 3,094.7 0.069 2 0.028 2 98.49 40.58 14-bus 14 5,198.9 5,120.5 0.08 2 0.047 2 98.64 58.75 30-bus 30 583.3 573.4 0.1 4 0.049 2 98.30 49.00 57-bus 57 41,760 41,582.0 0.165 4 0.140 3 99.57 84.85 118-bus 118 129,890.7 126,521.0 0.92 6 0.391 4 97.41 42.50 300-bus 300 719,753.6 690,226.0 1.17 4 0.890 3 95.90 76.07 582-bus 582 334,044.4 326,596.0 3.187 6 2.615 4 97.77 82.05 10

[00228] As further shown, the method 600 provides a greater minimization of the power generation cost as compared to the bus-wise optimal power flow (BWOPF) formulation, especially for power systems having a larger quantity of busses. By way of example, the cost generation solution using the BWOPF method for a 582-bus system is approximately $334,000, while the generation cost solution using the LWOPF method 600 is approximately 15 $326,000. Accordingly, the LWOPF method provides a greater minimization of the power cost, and accordingly, avoids potential economic loss. It will be appreciated that, in real power systems, which include thousands of buses, the cost savings are greatly increased using the LWOPF method.

[00229] Table 8 also demonstrates that the LWOPF method is able to achieve a better 20 minimization of the cost function in less execution time, and in some cases, in less iterations of the method 600. For example, in a 582-bus power system, the LWOPF method is able to - 52 - CA 3028644 2018-12-27 generate a high quality solution to the objective function in 4 iterations and 2.615 seconds, as opposed to 6 iterations and 3.817 seconds for the BWOPF method. [00230) Table 9, below, demonstrates the SLP based LWOPF method using different static load models. The load models includes constant power (PC), constant current (CC), constant 5 impedance (Cl), and a ZIP model. Table 9. Performance of SLP based LWOPF with various load models Load model Objective: Execution lter System Generation Cost ($) time(seconds) # Constant Power (CP) 3,094.7 0.028 2 ( / ) Constant Current (CC) 3,094.7 0.028 2 ::J .0 Constant Impedance (Cl) 3,094.5 0.029 2 I (0 ZIP Model (40% CP, 30% CC and 30% Cl) 3,095.2 0.032 2 Constant Power (CP) 4,1582.0 0.140 3 ( / ) Constant Current (CC) 4, 1598.0 0.144 3 ::J .0 Constant Impedance (Cl) 4,1622.0 0.149 3 I 1'- LO ZIP Model (40% CP, 30% CC and 4,1654.0 0.152 3 30% Cl) Constant Power (CP) 126,521.0 0.391 4 ( / ) ::J Constant Current (CC) 126,578.0 0.416 4 .0 I Constant Impedance (Cl) 126,674.0 0.428 4 co T""" ZIP Model (40% CP, 30% CC and T""" 30% Cl) 126,988.0 0.460 4 [00231) As shown, the performance of the SLP based LWOPF is consistent and reliable and is not otherwise susceptible to change based on the load model. In particular, for each 10 bus system, the objective solution to the cost function does not substantially change as between different system models. Further, the execution time and the number of iterations required to minimize the objective function also does not change based on the static load model. [00232) Referring now to FIG. 8, there shown a table 800 comparing an SLP formulation 15 for solving the OPF based on a line-wise set of power balance (LWBP) equations versus using a set of bus-wise power balance (BWPB) equations. [00233) In particular, FIG. 8 shows that LWBP equations includes fewer higher order terms as compared to the BWBP, and includes fewer sinusoidal terms in the "A" and "C" matrices, where the "A" matrix is the constraint coefficient matrix and the "C" matrix is the objective - 53 - CA 3028644 2018-12-27 function coefficient vector. Accordingly, as stated previously, this reduces computation burden, and allows the LWOPF to execute faster than the bus-wise method. FIG. 8 also shows that the LWBP have better conditioned Jacobian matrices.

[00234] While the applicant's teachings described herein are in conjunction with various 5 embodiments for illustrative purposes, it is not intended that the applicant's teachings be limited to such embodiments as the embodiments described herein are intended to be examples. On the contrary, the applicant's teachings described and illustrated herein encompass various alternatives, modifications, and equivalents, without departing from the embodiments described herein, the general scope of which is defined in the appended 10 claims. - 54 - CA 3028644 2018-12-27

Claims

55 CLAIMS:

1. A method for conducting a power flow analysis for a power system, the power system comprising at least one transmission branch having a first end coupled to a first bus and a second end coupled to a second bus, the method being implemented by at least one server comprising at least one processor, the method comprising: initializing values for a problem vector (𝑊), wherein the problem vector (𝑊) is defined as 𝑊 = [𝑈 𝛿 𝑃𝐹 𝑃𝑆 𝑄𝐹 𝑄𝑆], where: 𝑈 is a vector of squared voltage magnitudes for each bus of each of the at least one transmission branch; 𝛿 is a vector of voltage angles for each bus of each of the at least one transmission branch; 𝑃𝐹 is a vector of active power flow in each first end of each of the at least one transmission branch; 𝑃𝑆 is a vector of active power flow in each second end of each of the at least one transmission branch; 𝑄𝐹 is a vector of reactive power flow in each first end of each of the at least one transmission branch, and 𝑄𝑆 is a vector of reactive power flowing in each second end of each of the at least one transmission branch; modelling the power system using a set of line-wise power balance equations that are functions of the problem vector (𝐹𝑇(𝑊)) , wherein the set of line-wise power balance equations (𝐹𝑇(𝑊)) is defined as: 𝐹𝑇(𝑊) = [ 𝐹𝐹(𝑊) 𝐹𝑆(𝑊) 𝐹𝐴(𝑊) 𝐹𝐵(𝑊) 𝐹𝑃(𝑊) 𝐹𝑄(𝑊)] where: 𝐹𝐹(𝑊) is a vector of line based voltage functions in relation to each first end of each of the at least one transmission branch; CA 3028644 Date reçue / Received date 2025-05-09 56 𝐹𝑆(𝑊) is a vector of line-based voltage function in relation to each second end of each of the at least one transmission branch; 𝐹𝐴(𝑊) is a vector of phase angle functions in relation to each first end of each of the at least one transmission branch; 𝐹𝐵(𝑊) is a vector of phase angle functions in relation to each second end of each of the at least one transmission branch; 𝐹𝑃(𝑊) is a vector of functions for bus-wise real power injections; 𝐹𝑄(𝑊) is a vector of functions for bus-wise reactive power injections, and wherein the set of line-wise based equations is formulated in terms of the square magnitude of voltages at each bus in each of the at least one transmission branch in the power system; solving the set of line-wise power balance equations (𝐹𝑇(𝑊)) using a numerical technique to obtain a line-wise solution vector; based on the line-wise solution vector, determining whether the power system is in a state of voltage stability; and implementing a remedial action to stabilize the power system when the power system is not in a state of voltage stability.

2. The method of claim 1, wherein solving the set of line-wise power balance equations using a Newton-Raphson technique comprising: determining an incremental change to the problem vector (Δ𝑊) by solving a first-order Taylor’s series approximation of the set of line-wise power balance equations (𝐹𝑇(𝑊)): Δ𝑊 = [ 𝜕𝐹𝑇(𝑊) 𝜕𝑊 ] −1 ∙ [𝑆𝑇 − 𝐹𝑇(𝑊)] where: Δ𝑊 is the incremental change to the problem vector; [ 𝜕𝐹𝑇(𝑊) 𝜕𝑊 ] −1 is an inverse of a Jacobian matrix for the set of line-wise power balance equations; and CA 3028644 Date reçue / Received date 2025-05-09 57 𝑆𝑇 is a desired line-wise solution vector for the set of line-wise power balance equations.

3. The method of claim 2, wherein the desired line-wise solution vector (𝑆𝑇) is defined according to 𝑆𝑇 = [0 0 0 0 [𝑃𝐷 − 𝑃𝐺] [𝑄𝐷 − 𝑄𝐺]]𝑇, where: 𝑃𝐷 is a vector of active power demand in each of the at least one transmission branches; 𝑃𝐺 is a vector of active power generation in each of the at least one transmission branches; 𝑄𝐷 is a vector of reactive power demand in each of the at least one transmission branches, and 𝑄𝐺 is a vector of reactive power generation in each of the at least one transmission branches.

4. The method of any one of claims 2 or 3, further comprising: generating an updated problem vector (𝑊′) according to the equation: 𝑊′ = 𝑊 + Δ𝑊; solving the set of line-wise power balance equations to generate an intermediate solution vector using the updated problem vector (𝑊′); and based on the intermediate solution vector, determining whether a tolerance condition is satisfied, wherein when the tolerance condition is not satisfied, the method comprises iteratively solving for the updated problem vector using the first-order Taylor’s series approximation until the intermediate solution vector satisfies the tolerance condition, and when the tolerance condition is satisfied, the intermediate solution vector comprises the linewise solution vector.

5. The method of claim 4, wherein the tolerance condition is expressed as |𝑆𝑇 − 𝐹𝑇(𝑊)| ≥ 𝑇, where 𝑇 is a pre-defined tolerance threshold that is zero. CA 3028644 Date reçue / Received date 2025-05-09 58 6. The method of claim 5, wherein the line-wise solution vector yields voltage values at each bus of each of the at least one transmission branch of the power system, and line-wise real and reactive power flows in each of the at least one transmission branch of the power system.

7. The method of claim 6, wherein the set of line-wise based equations is formulated in terms of the square magnitude of the voltages at each bus in each of the at least one transmission branch in the power system.

8. The method of claim 1, wherein the functions for bus-wise real power injection, and the functions for bus-wise reactive power injection, are expressed in linear terms and the set of line-wise power balance equations include only first order and second order terms.

9. The method of any one of claims 1 to 8, wherein initializing values for the problem vector (𝑊), further comprises: reading a set of known data values from at least one sensor operatively coupled to the power system.

10. The method of any one of claims 1 to 9, wherein the transmission branch comprises at least one of a transmission cable, a transmission line, and a transformer and each of the first bus and the second bus comprise at least one of a load bus and a generator bus.

11. The method of claim 2, wherein the Jacobian matrix includes, for each of the at least one transmission branches, a first vector element ( 𝜕𝐹𝐹(𝑊) 𝜕𝑈 ) and a second vector element ( 𝜕𝐹𝑆(𝑊) 𝜕𝑈 ), where the first vector element defines a first voltage collapse index (VCI) for a first end of the at least one transmission branch, and the second vector element defines a second VCI for a second end of the at least one transmission branch.

12. The method of claim 11, wherein determining whether the power system is in a state of voltage stability comprises determining whether at least one of the first VCI and the second VCI for a respective transmission branch is zero, and the remedial action to stabilize the power system comprises at least one of: shedding a load bus in the power system, increasing a power output at least one generator bus in the power system, re-activating at least one CA 3028644 Date reçue / Received date 2025-05-09 59 standby generator bus in the power system, and re-inserting standby transmission elements in the power system.

13. The method of any one of claims 4 or 5, wherein determining whether the power system is in a state of voltage stability further comprises: performing an N-1 contingency analysis for the power system.

14. A system for conducting a power flow analysis for a power system, wherein the power system comprises at least one transmission branch having a first end coupled to a first bus and a second end coupled to a second bus, the system comprising: at least one server comprising at least one processor, wherein the at least one server is configured to perform the method of any one of claims 1 to 13.

15. A method for minimizing an objective function in respect of a power system, the power system comprising at least one transmission branch having a first end coupled to a first bus and a second end coupled to a second bus, the method being implemented using at least one server comprising at least one processor, the method comprising: (a) determining initial values for a control vector (𝑊𝐶); (b) using the control vector (𝑊𝐶), solving a set of line-wise power balance equations (𝐹𝑇(𝑊′)) to obtain a problem vector (𝑊′), wherein the problem vector (𝑊′) comprises the control vector (𝑊𝐶) and a dependent vector (𝑊𝐷) and the set of linewise power balance equations (𝐹𝑇(𝑊)) is defined as: 𝐹𝑇(𝑊) = [ 𝐹𝐹(𝑊) 𝐹𝑆(𝑊) 𝐹𝐴(𝑊) 𝐹𝐵(𝑊) 𝐹𝑃(𝑊) 𝐹𝑄(𝑊)] wherein: 𝐹𝐹(𝑊) is a vector of line based voltage functions in relation to each first end of each of the at least one transmission branch; 𝐹𝑆(𝑊) is a vector of line-based voltage function in relation to each second end of each of the at least one transmission branch; CA 3028644 Date reçue / Received date 2025-05-09 60 𝐹𝐴(𝑊) is a vector of phase angle functions in relation to each first end of each of the at least one transmission branch; 𝐹𝐵(𝑊) is a vector of phase angle functions in relation to each second end of each of the at least one transmission branch; 𝐹𝑃(𝑊) is a vector of functions for bus-wise real power injections, and 𝐹𝑄(𝑊) is a vector of functions for bus-wise reactive power injections; (c) using the problem vector (𝑊′), solving a line-wise linear optimization formulation to determine an optimal incremental change to the control vector (Δ𝑊𝐶), wherein the problem vector (𝑊′) is defined as 𝑊′ = [𝑈, 𝛿, 𝑃𝐹, 𝑃𝑆, 𝑄𝐹, 𝑄𝑆, 𝑃𝐺, 𝑄𝐺]𝑇, and the dependent vector (𝑊𝐷) is defined as 𝑊𝐷 = [𝑈, 𝛿, 𝑃𝐹, 𝑃𝑆, 𝑄𝐹, 𝑄𝑆]𝑇, and the control vector (𝑊𝐶) is defined as 𝑊𝐶 = [𝑈, 𝑃𝐺]𝑇, wherein: 𝑈 is a vector of squared voltages magnitudes for each bus of each of the at least one transmission branch; 𝛿 is a vector of voltage angles for each bus of each of the at least one transmission branch; 𝑃𝐹 is a vector of active power flow in each first end of each of the at least one transmission branch; 𝑃𝑆 is a vector of reactive power flow in each second end of each of the at least one transmission branch; 𝑄𝐹 is a vector of reactive power flow in each first end of each of the at least one transmission branch; 𝑄𝑆 is a vector of reactive power flowing in each second end of each of the at least one transmission branch; 𝑃𝐺 is a vector of active power generation in each of the at least one transmission branches, and 𝑄𝐺 is a vector of reactive power generation in each of the at least one transmission branches; (d) updating the control vector (𝑊𝐶) based on the optimal incremental change in the control vector (Δ𝑊𝐶) and CA 3028644 Date reçue / Received date 2025-05-09 61 (e) determining whether the objective function is minimized by determining whether a set of stopping criteria is satisfied based on the updating of the control vector (𝑊𝐶), wherein when the set of stopping criteria is not satisfied, the method comprises iteratively repeating (b) to (d) until the set of stopping criteria is satisfied, and wherein determining initial values for the control vector (𝑊𝐶) further comprises: reading a set of data from at least one sensor operatively coupled to the power system.

16. The method of claim 15, wherein solving the set of line-wise power balance equations (𝐹𝑇(𝑊′)) to obtain the problem vector (𝑊′), comprises the method of any one of claims 1 to 2 and 4 to 5.

17. The method of claim 16, when depending from claim 2, wherein a desired line-wise solution vector is defined as 𝑆𝑇 = [0 0 0 0 [𝑃𝐷] [𝑄𝐷]]𝑇, wherein: 𝑃𝐷 is a vector of active power demand in each of the at least one transmission branch, and 𝑃𝐺 is a vector of active power generation in each of the at least one transmission branch.

18. The method of any one of claims 15 to 17, wherein the objective function 𝑓𝑜 (𝑊′) is a second-order power generation cost function expressed as: 𝑓𝑜 (𝑊′) = Σ 𝑐1𝑔 𝑁𝐺 𝑔=1 + 𝑐2𝑔 ∙ 𝑃𝐺𝑔 + 𝑐3𝑔 ∙ 𝑃𝐺𝑔 2 wherein 𝑁𝐺 is a total number of power generators in the power system, 𝑃𝐺𝑔 is the active power generation of a g-th power generator, 𝑐1𝑔 is a first constant for the g-th power generator, 𝑐2𝑔 is a second constant for the g-th power generator, and 𝑐3𝑔 is a third constant for the g-th power generator.

19. The method of any one of claims 15 to 18, wherein the line-wise linear optimization formulation comprises a linearized objective function, a linearized set of equality constraints, and a linearized set of inequality constraints. CA 3028644 Date reçue / Received date 2025-05-09 62 20. The method of claim 19 when depending from claim 18, wherein the linearized objective function is expressed as: Σ[(𝑐2𝑔 + 2 ∙ 𝑐3𝑔 ∙ 𝑃𝐺𝑔) ∙ Δ𝑃𝐺𝑔] 𝑁𝐺 𝑔=1 wherein Δ𝑃𝐺𝑔 is a linear increment to the vector of active power generation.

21. The method of claim 19, wherein the linearized set of equality constraints is expressed as: [𝐽 𝐷]Δ𝑊′ = Δ𝑆𝑇 wherein 𝐽 is a first jacobian matrix, 𝐷 is a second Jacobean matrix, and Δ𝑆𝑇 is an incremental change to the line-wise solution vector, wherein Δ𝑆𝑇 = [0 0 0 0 0 0]𝑇.

22. The method of any one of claims 15 to 21, wherein the method comprises solving the linear optimization formula using a linear programming method.

23. The method of claim 15, wherein determining whether the set of stopping criteria is satisfied comprises determining whether a stopping condition is satisfied, wherein the stopping conditioned is defined as: |𝑓𝑜 (𝑊′𝑘 ) − 𝑓𝑜 (𝑊′𝑘−1)| |𝑓𝑜 (𝑊′𝑘 )| ≤ 𝛾 where: 𝛾 is a pre-determining mismatch tolerance; 𝑓𝑜 (𝑊′𝑘 ) is a first output of the objective function using values in the problem vector (𝑊′) determined in a current iteration (k) of steps (a) to (d) of the method of claim 15; and 𝑓𝑜 (𝑊′𝑘−1) is a second output of the objective function using values in the problem vector (𝑊′) determined in a previous iteration (k-1) of steps (a) to (d) of the method of claim 15.

24. The method of claim 15, wherein the set of stopping criteria is satisfied when a number of iterations of (b) to (d) is larger than a predetermined maximum number of iterations. CA 3028644 Date reçue / Received date 2025-05-09 63 25. The method of claim 15, wherein the set of stopping criteria is satisfied when a set of equality and inequality constraints is satisfied, the set of equality constraints being defined by 𝐹𝑇(𝑊)=ST.

26. The method of any one of claims 15 to 25, further comprising: based on the control vector (𝑊𝐶), determining whether the power system is in a state of voltage stability by performing an N-1 contingency analysis, and when the power system is not in a stage of voltage stability implementing a remedial action to stabilize the power system where the remedial action comprises at least one of: shedding a load bus of the power system, increasing a power output at least one generator bus in the power system, re-activating at least one standby generator bus in the power system, and re-inserting standby transmission elements in the power system.

27. A system for minimizing an objective function in respect of a power system, the power system comprising at least one transmission branch having a first end coupled to a first bus and a second end coupled to a second bus, the system comprising: at least one server comprising at least one processor, wherein the at least one server is configured to perform the method as defined by any one of claims 15 to 26. CA 3028644 Date reçue / Received date 2025-05-09