SYSTEMS AND METHODS FOR DYNAMIC RATING OF POWER GRIDS

Abstract

A system for dynamic rating of a power grid may include a plurality of terminal units, and a controller. The terminal units may detect a voltage phasor and a current phasor at nodes of the power grid. The controller may, based on the voltage phasors and the current phasors of the plurality of nodes, determine a dynamic thermal stability power rating for each line, a dynamic angular stability power rating for each node, and a dynamic voltage stability power rating for each node. The controller may, based on the dynamic thermal stability power rating, the dynamic angular stability power rating, and the dynamic voltage stability power rating, determine a dynamic system rating for the power grid. The controller may control the power grid in response to the dynamic system rating.

Description

TECHNICAL FIELD

The present disclosure relates to power grids, and in particular, to systems and methods for dynamic rating of power grids.

BACKGROUND

A power grid includes a network of nodes and lines. A power grid may be operated based on a power rating, such that power transmitted is maintained lower than the power rating. Maintaining the power transmission lower than the power rating may avoid node and line faults, for example, associated with line overheating, node voltage collapse, or unstable power swing which results in the collapse of the whole grid due to the system loses frequency synchronization.

However, the use of static power ratings may not account for changing conditions in the grid, because the static power rating is normally made before the grid operation, based on a consideration of the worst operating conditions. For example, the thermal line rating is based on the worst scenario that the line can be operated with, which may result in the real rating power being lower than the actual rating capacity. Using lower power than the actual capacity may result in congestion of power transmission, leading to the power rejection of renewable sources, such as windfarms.

A need remains for dynamic rating of power grids that accounts for varying conditions of the power grid.

SUMMARY

The present disclosure relates to power grids, and in particular, to systems and methods for dynamic rating of power grids. The dynamic system rating may include a dynamic line thermal rating (to avoid overheat), a dynamic angular stability power rating (to avoid losing angular stability), and a dynamic voltage stability power rating (to avoid voltage collapse).

In embodiments, the present disclosure describes a system for dynamic rating of a power grid. The power grid may include a network of a plurality of nodes coupled by a plurality of lines. The system may include a plurality of terminal units, and a controller. Each terminal unit of the plurality of terminal units may be configured to detect a voltage phasor and a current phasor of a respective node of the plurality of nodes and generate signals indicative of the voltage phasor and the current phasor. The controller may be configured to receive signals from the plurality of terminal units indicative of the voltage phasors and the current phasors of the plurality of nodes. The controller may be further configured to determine, based on the voltage phasors and the current phasors of the plurality of nodes, a dynamic thermal stability power rating for each line of the plurality of lines. The controller may be further configured to determine, based on the voltage phasors and the current phasors of the plurality of nodes, a dynamic angular stability power rating for each node of the plurality of nodes. The controller may be further configured to determine, based on the voltage phasors and the current phasors of the plurality of nodes, a dynamic voltage stability power rating for each node of the plurality of nodes. The controller may be further configured to determine, based on the dynamic thermal stability power rating, the dynamic angular stability power rating, and the dynamic voltage stability power rating, a dynamic system rating for the power grid. The controller may be further configured to control (for example, optimally control) the power grid in response to the dynamic system rating.

In embodiments, the present disclosure describes a method for dynamic rating of a power grid including a network of a plurality of nodes coupled by a plurality of lines. The method may include sending, to a controller, signals from a plurality of terminal units indicative of voltage phasors and current phasors at each node of the plurality of nodes. The method may further include determining, by the controller, based on the voltage phasors and the current phasors of the plurality of nodes, a dynamic thermal stability power rating for each line of the plurality of lines. The method may further include determining, by the controller, based on the voltage phasors and the current phasors of the plurality of nodes, a dynamic angular stability power rating for each node of the plurality of nodes. The method may further include determining, by the controller, based on the voltage phasors and the current phasors of the plurality of nodes, a dynamic voltage stability power rating for each node of the plurality of nodes. The method may further include determining, by the controller, based on the dynamic thermal stability power rating, the dynamic angular stability power rating, and the dynamic voltage stability power rating, a dynamic system rating for the power grid. The method may further include controlling, by the controller, the power grid in response to the dynamic system rating.

Additional systems, methods, apparatus, features, and aspects can be realized through the techniques of various embodiments of the disclosure. Other embodiments and aspects of the disclosure are described in detail herein and are considered a part of the claimed subject matter. Other features can be understood and will become apparent with reference to the description and to the drawings.

BRIEF DESCRIPTION OF THE DRAWINGS

Having thus described the disclosure in general terms, reference will now be made to the accompanying drawings, which are not necessarily drawn to scale.

FIG. 1 is a block diagram showing a system including a controller for controlling a power grid including a network of nodes and lines.

FIG. 2 is a block diagram showing a scheme for determining a dynamic system rating of the power grid of FIG. 1.

FIG. 3 is a block diagram showing a scheme for determining a dynamic thermal stability power rating of a line.

FIG. 4 is a chart showing an example of a numerical relationship between a length of a line and an instantaneous temperature of the line based on elongation coefficients.

FIG. 5 is a block diagram showing a line between nodes ‘m’ and ‘n’ of the power grid of FIG. 1.

FIG. 6 is a block diagram showing a scheme for determining a dynamic angular stability power rating for the line of FIG. 5.

FIG. 7 is a chart showing the use of equal area criterion for determining the dynamic angular stability rating using the scheme of FIG. 6.

FIG. 8 is a block diagram showing a fault assumed on a line ‘mp’ adjacent node ‘m’ of the power grid of FIG. 1.

FIG. 9 is a block diagram showing a scheme for determining an impedance matrix associated with the fault of FIG. 8.

FIG. 10 is a block diagram showing a scheme for determining an impedance matrix after the fault of FIG. 8 is cleared.

FIG. 11 is a block diagram showing an equivalent Thevenin's circuit for a node ‘n’ of the power grid of FIG. 1 for determining a dynamic voltage stability power rating.

FIG. 12 is a block diagram showing a scheme for determining a dynamic voltage stability power rating based on the Thevenin's circuit of FIG. 11.

FIG. 13 is a chart showing the crossing of example reactive power curves with a consumption reactive power curve for determining maximum active power.

FIG. 14 is a chart showing the relationship between active power and voltage, and the maximum active power at marginal point of voltage stability.

FIG. 15 is a chart showing the determination of maximum active power based on reactive power curves.

FIG. 16 is a block diagram showing the controller of FIG. 1 controlling the power grid to increase power transmission from a renewable energy source while reducing power transmission loss.

Embodiments of the disclosure are described more fully below with reference to the accompanying drawings, in which example embodiments of the disclosure are shown. This disclosure may, however, be embodied in many different forms and should not be construed as limited to the example embodiments set forth herein; rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the present disclosure to those skilled in the art. Like numbers refer to like, but not necessarily the same or identical, elements throughout.

The following embodiments are described in sufficient detail to enable at least those skilled in the art to understand and use the disclosure. It is to be understood that other embodiments would be evident based on the present disclosure and that process, mechanical, material, dimensional, process equipment, and parametric changes may be made without departing from the scope of the present disclosure.

In the following description, numerous specific details are given to provide a thorough understanding of various embodiments of the present disclosure. However, it will be apparent that the present disclosure may be practiced without these specific details. In order to avoid obscuring the present disclosure, some well-known system configurations and process steps may not be disclosed in full detail. Likewise, the drawings showing embodiments of the disclosure are semi-diagrammatic and not to scale and, particularly, some of the dimensions are for the clarity of presentation and may be exaggerated in the drawings. In addition, where multiple embodiments are disclosed and described as having some features in common, for clarity and ease of illustration, description, and comprehension thereof, similar and like features will ordinarily be described with like reference numerals even if the features are not identical.

DETAILED DESCRIPTION

The present disclosure relates to power grids, and in particular, to systems and methods for dynamic rating of power grids. A power grid may include a network of nodes and lines, with nodes representing power sources or sinks, and lines representing connections between the nodes. To avoid failure, losses, or fatigue, one or more elements of the power grid may be operated based on a rating, or a predetermined maximum power level associated with a node or a line. Certain ratings for power grids are static, for example, being based generally on parameters associated with components of the grids, but which may not account for changes in the grid. Static ratings for power grids may not account for changing operational conditions, such as changes in power from renewable sources, the consumption of reactive power from load, the changing of system topologies due to the breaking of lines, the change of the ambient conditions, or other conditions. Operating a power grid based on a static rating may result in drawing or transmitting lower power than the actual capacity of the power grid. For example, if the power grid is controlled in compliance with static ratings, the power drawn from a node associated with a renewable source (such as solar power, wind power, or other sources) may be much lower than the maximum capacity of the source, which may result in the unnecessary power being rejected from the renewable resources.

In order to increase the power drawn from renewables, systems and methods according to the present disclosure provide a dynamic system rating for the grid. Operating or controlling the power grid based on a dynamic system rating may transfer more renewable energy (for example, compared to power transferred based on a static rating of transmission power capacity), and dynamically varying adaptive to the operational condition of the power grid, evacuating more power from renewable sources. Certain dynamic ratings may only account for a limited number of parameters, for example, parameters associated with a line, but not account for parameters associated with nodes or faults. Certain dynamic ratings may require sensing ambient environmental conditions, for example, temperature, wind speed, atmospheric pressure, and the like. Dynamic system ratings according to the present disclosure may not need sensing or determining of environmental conditions, and may account for parameters in addition to those only associated with lines.

In embodiments, the present disclosure describes a system for dynamic rating of a power grid. The power grid may include a network of a nodes coupled by lines. The system may include terminal units, and a controller. Each terminal unit of terminal units may be configured to detect a voltage phasor and a current phasor of a respective node of the nodes and generate signals indicative of the voltage phasor and the current phasor. The controller may be configured to receive signals from terminal units indicative of the voltage phasors and the current phasors of nodes. The controller may be further configured to determine, based on the voltage phasors and the current phasors of nodes, a dynamic thermal stability power rating for each line of the lines. The controller may be further configured to determine, based on the voltage phasors and the current phasors of nodes, a dynamic angular stability power rating. The controller may be further configured to determine, based on the voltage phasors and the current phasors of nodes, a dynamic voltage stability power rating for each node of the nodes. The controller may be further configured to determine, based on the dynamic thermal stability power rating, the dynamic angular stability power rating, and the dynamic voltage stability power rating, a dynamic system rating for the power grid. The controller may be further configured to control the power grid in response to the dynamic system rating. In embodiments, the controller may be configured to optimally control the power grid. Thus, dynamic system ratings according to the present disclosure may account for parameters associated with different components of the network, and without requiring sensing environmental or ambient conditions.

Systems and methods according to the present disclosure may allow drawing more power from renewable sources, while accounting for thermal, angular, and voltage stabilities. Further, the overall power flow may be increased or optimized, which may reduce or minimize power transmission losses. Maximum rating for transmission of powers (or currents) based on the thermal, angular and voltage stabilities may be calculated substantially in real-time or on-line, based on dynamically changing operational conditions of the power grid. Based on the determined dynamic rating power, the power flow and generation power distribution may be re-arranged or reconfigured for increasing power drawn from renewable sources, while reducing power transmission loss.

FIG. 1 is a block diagram showing a system 100 including a controller 102 for controlling a power grid 104 including a network 106 of nodes 108 and lines 110. For example, the network 106 may include nodes 108 connected by lines 110. A node of the nodes 108 may be connected to one or more other nodes of nodes 108 by respective lines of lines 110. Thus, nodes 108 may be coupled by lines 110. In embodiments, at least one node of the nodes 108 is associated with a renewable power source. For example, the renewable power source may include solar energy, wind energy, or a stored source of energy, for example, an electric vehicle battery connected to a charging node. In embodiments, the power grid 104 may include two or more renewable sources. For example, in FIG. 1, “Bus 1” and “Bus 4” provide connections with renewable sources.

The system 100 may include terminal units 112. Each terminal unit of terminal units 112 may be configured to detect a voltage phasor and a current phasor of a respective node of the nodes 108 and generate signals indicative of the voltage phasor and the current phasor. In embodiments, at least one terminal unit of terminal units 112 may include a synchronous phasor measurement unit (PMU), a remote terminal unit (RTU), a feeder terminal unit (FTU), a zonal client, or a zonal agent. The PMU may be a synchronous PMU. In embodiments, each terminal unit of terminal units includes a PMU. The terminal units 112 may measure the voltage phasors (magnitude and phase) of each bus or node, and current phasors (magnitude and phases) flowing from each bus or node. Synchronous PMUs may synchronously measure the voltage phasors and current phasors. Ultimately, terminal units 112 may generate signals indicative of the phasors, and transmit the signals. In embodiments, the signals are received by the controller 102.

The controller 102 may include a local computer, a remote server, processing circuitry, a central unit, or any device including a processor, memory, storage, and a communication module. The controller 102 may be wired or wireless. In embodiments, the controller includes a power gateway or a server. For example, the power gateway may include a GE Power Gateway (General Electric, Boston, MA). The controller 102 may be configured to sense one or more parameters associated with the power grid 104, and control or regulate the operation of the power grid 104, for example, by controlling one or more components of the power grid 104. In embodiments, the controller 102 is configured to regulate, connect, or disconnect one or more of nodes 108 or lines 110 to control or divert the flow of power transmitted through the network 106.

In embodiments, the controller 102 is configured to receive signals from terminal units 112 indicative of the voltage phasors and the current phasors of nodes 108. The current phasors may be associated with respective lines of lines 106 coupled to a respective node of the nodes 108.

FIG. 2 is a block diagram showing a scheme for determining a dynamic system rating of the power grid 104 of FIG. 1.

The controller 102 may be configured to implement the scheme shown in FIG. 2. For example, the controller 102 may be configured to determine different stability ratings, such as a dynamic line thermal stability power rating, a dynamic angular stability power rating, and a dynamic voltage stability rating.

The controller 102 may determine the synchronous voltage and current phasors based on the signals received from the terminal 112. The controller 102 may, based on the corresponding voltage and current phasors (I_i) of two terminals of each line, determine the dynamic line thermal power rating of each line (P_ThRate) The controller 102 may, based on the voltage phasor (U_i) and its phase angle (δ_i) of the two terminals of each line, determine the dynamic power rating for angular stability margin is calculated (P_AngRate) The controller 102 may, based on the voltage phasor of each bus junction and the currents of all outgoing lines, determine the dynamic rating power of the injection power to this junction is calculated (P_InjectMaxi) Based on these dynamic power restraints of each line and each junction, the controller may determine a flow control strategy to re-arrange the power flow distribution for increasing (or maximizing) the draw of renewable power and reducing (or minimizing) the overall grid transmission loss. For example, the controller may determine a tap changer position of transformers and compensated reactive power of each bus junction.

FIG. 3 is a block diagram showing a scheme for determining a dynamic thermal stability power rating of a line of the power grid of FIG. 1. The controller 102 may be further configured to determine, based on the voltage phasors and the current phasors of nodes, a dynamic thermal stability power rating for each line of the lines 110. The controller 102 may determine an instantaneous conductor temperature (T_c) based on the measured voltage and current phasors of two terminals of each branch/link of the grid. For example, the controller 102 may determine a real conductor (line) length by using EQUATION 1.

$\begin{array}{cc}\left[\begin{array}{c}{V}_R\\ I_R\end{array}\right]=\left[\begin{array}{cc}A\left(l\right)& B\left(l\right)\\ C\left(l\right)& D\left(l\right)\end{array}\right]\left[\begin{array}{c}{V}_S\\ I_S\end{array}\right]& \left(\mathrm{Equation}1\right)\end{array}$

In EQUATION 1:

$A\left(l\right)=\cosh \left(\gamma l\right);$ $B\left(l\right)=-Z_C\sinh \left(\gamma l\right);$ $C\left(l\right)=-\frac{\sinh \left(\gamma l\right)}{z_C};$ $D\left(l\right)=\cosh \left(\gamma l\right);$ $\mathrm{and},$ $\gamma =\sqrt{z_1{y}_1},$ $Z_C=\sqrt{z_1/y_1},$

z₁, y₁ are parameters of the line, which are respectively the series impedance per unit length (Ohm/km, or Ohm/mi) and shunt admittance per unit (S/km, or S/mi).

• These may be obtained from memory or input by a user.

Subsequently, the real length of the conductor can be mapped to the real-time instantaneous temperature Tc, by solving EQUATION 2.

I=I _0,T=20+ρ₁(T_c−T₂₀)+ρ₂(T_c−T₂₀)²   (Equation 2)

In EQUATION 2, ρ₁ and ρ₂ are elongation coefficients of the line, which can be obtained by the line type and its parameters.

FIG. 4 is a chart showing an example of a numerical relationship between a length of a line and an instantaneous temperature of the line based on elongation coefficients.

The controller 102 may further calculate a steady state conductor temperature (T_CSS). For example, the controller 102 may continuously calculate a number of samples, for example, at least 3 samples, of instantaneous temperature (Tc) with time interval Ts. Ts may be any suitable interval, for example, 1 s, 5 s, or 10 s, as T_c(n−2), T_c(n−1), and T_c(n). The controller 102 may further determine the time constant of the thermal dynamic process, α, by using EQUATION 3.

$\begin{array}{cc}\alpha =\frac{2}{T_S}\frac{T_C\left(n\right)-2T_C\left(n-1\right)+T_C\left(n-2\right)}{T_C\left(n\right)-T_C\left(n-2\right)}& \left(\mathrm{Equation}3\right)\end{array}$

The controller 102 may further calculate the steady state conductor temperature using EQUATION 4.

$\begin{array}{cc}{T}_{\mathrm{CSS}}=\frac{1}{\alpha }\frac{T_C\left(n\right)-T_C\left(n-1\right)}{T_S}+\frac{T_C\left(n\right)+T_C\left(n-1\right)}{2}& \left(\mathrm{Equation}4\right)\end{array}$

The controller 102 may further calculate the rating power of thermal stability, (P_ThRate) For example, the controller 102 may continuously calculate at least two steady state conductor temperatures, T_CSS(n) and T_CSS(n−1), accompanied with the measured corresponding currents, I (n) and I (n−1). The controller 102 may determine the maximum rating current by using EQUATION 5.

$\begin{array}{cc}{I}_{\max }\left(n\right)=\sqrt{I^2\left(n\right)+\frac{I^2\left(n\right)-I^2\left(n-1\right)}{T_{\mathrm{css}}\left(n\right)-T_{css}\left(n-1\right)}\left(T_{\max }-T_{css}\left(n\right)\right)}& \left(\mathrm{Equation}5\right)\end{array}$

In EQUATION 5, T_max is the maximum conductor temperature that the line can bear. The controller 102 may determine a maximum dynamic rating power of thermal stability using EQUATION 6.

P _ThRate =U _S I _Max  (Equation 6)

In embodiments, the controller 102 is configured to determine the dynamic thermal stability power rating for each line of the lines 110 without assessing environmental parameters. For example, the controller 102 may not sense or otherwise assess air temperature, air pressure, or air velocity, or other environmental parameters. In other embodiments, the controller 102 may be configured to determine the dynamic thermal stability power rating based on at least one environmental parameter.

The controller may be further configured to determine, based on the voltage phasors and the current phasors of nodes 108, a dynamic angular stability power rating for each node of the nodes 108.

FIG. 5 is a block diagram showing a line between nodes ‘m’ and ‘n’ of the power grid of FIG. 1. FIG. 6 is a block diagram showing a scheme for determining a dynamic angular stability power rating for a line between nodes ‘m’ and ‘n’ of FIG. 5.

Controller 102 may determine the electromagnetic powers transmitted in each line by using the measured voltages and angle difference between the two terminals and the impedance matrix of the network, at the following scenarios:

• scenario 1, the power of the normal status, P_M1;
• scenario 2, the power of faulty status (assumed fault applied onto the next lines), P_M2; and
• scenario 3, the power of the status after fault cleared, P_M3.

Controller 102 may further determine the angle of the instance that the fault was cleared (δ_CR). Controller 102 may further calculate the dynamic rating power for angular stability, based on equal area criterion. FIG. 7 is a chart showing the use of equal area criterion for determining the dynamic angular stability rating using the scheme of FIG. 6.

Under scenario 1, the controller 102 determines the power at normal status. For example, the controller 102 may determine the transfer impedance from the impedance matrix ZN. If the line terminal (node) number is ‘m’ and ‘n’ respectively, the transfer impedance, noted as Zmn, is the element of mth row and nth column of impedance matrix ZN. Then, EQUATION 7 may be used to determine the power of the normal status, P_M1.

$\begin{array}{cc}{P}_{M1}=\frac{U_m{U}_n}{X_{mn}}& \left(\mathrm{Equation}7\right)\end{array}$

Under scenario 2, the controller 102 determines the power at faulty status. FIG. 8 is a block diagram showing a fault assumed on a line ‘mp’ adjacent node ‘m’ of the power grid of FIG. 1. Assuming a fault at terminal ‘m’ of another line mp, this equals an impedance ZF, connected onto bus m, where ZF is fault additional impedance. For example, if the single-phase-to-ground fault is considered, the fault additional impedance is the sum of negative sequence and zero-sequence equivalent impedance at bus m, that is, Z_F=Z_mm2+Z_mm0. If the phase-to-phase fault is considered, the additional impedance is Z_F=Z_mm2.

The controller 102 may add this additional impedance Z_F into the impedance matrix Z_N, to formulate a new impedance matrix Z_NF, as shown in FIG. 9. FIG. 9 is a block diagram showing a scheme for determining an impedance matrix associated with the fault of FIG. 8. The controller 102 may determine the renewed element of impedance matrix Z_NF can be obtained by EQUATION 8.

$\begin{array}{cc}{Z}_{NF}\left(i,j\right)=Z_N\left(i,j\right)-\frac{Z_N\left(i,m\right)Z_N\left(m,j\right)}{Z_N\left(m,m\right)+Z_F}& \left(\mathrm{Equation}8\right)\end{array}$

The transfer impedance of mn at faulty status is the element of mth row and nth column of the impedance matrix Z_NF: Z_mnF=Z_NF(m, n). The controller 102 may determine the voltages of bus m and n at fault status. The fault current may be determined using EQUATION 8.

$\begin{array}{cc}{I}_F=\frac{U_m}{Z_{mm}+Z_F}& \left(\mathrm{Equation}9\right)\end{array}$

In EQUATION 9, Z_mm is the element of mth row mth column of original impedance matrix Z_N. The voltage phasor of bus m and n at faulty status may be determined using EQUATIONS 10 and 11. (based on superposition, FIG. 9).

U _mF =U _m −Z _mm I _F  (Equation 10)

U _nF =U _n −Z _mn I _F  (Equation 11)

The controller 102 may determine the electro-magnetic power P_M2 at faulty status using EQUATION 12.

$\begin{array}{cc}{P}_{M2}=\frac{U_{mF}{U}_{nF}}{X_{mnF}}& \left(\mathrm{Equation}12\right)\end{array}$

Under scenario 3, the controller 102 may determine the power after fault clearance. After a predetermined time interval, for example, 50 ms, or 100 ms, or 200 ms, the fault on the line mp will be cleared. That is, the link branch mp will be open. This network where link branch mp is opened can be regarded as adding an impedance, which is negative impedance of line mp, between mp, as shown in FIG. 10. FIG. 10 is a block diagram showing a scheme for determining an impedance matrix after the fault of FIG. 8 is cleared.

In embodiments, the controller 102 is configured to determine the dynamic angular stability power rating by simulating applying a fault and clearing the fault at each node of the nodes 108.

For determining the response of the power grid to a fault and clearing of the fault, the controller 102 may determine the transfer impedance as the faulty branch has been opened (line mp opened). Effectively, an impedance, which is negative to the impedance of line pm, may be added between junction m and p. Then the impedance matrix (noted as Z_NC) of modified network may be determined using EQUATION 13.

$\begin{array}{cc}{Z}_{NC}\left(i,j\right)=Z_N\left(i,j\right)-\frac{\left[Z_N\left(i,m\right)-Z_N\left(i,p\right)\right]\left[Z_N\left(m,j\right)-Z_N\left(p,j\right)\right]}{Z_N\left(m,m\right)+Z_N\left(p,p\right)-2Z_N\left(m,p\right)-Z_{mp}}& \left(\mathrm{Equation}13\right)\end{array}$

In EQUATION 13, Z_N(i,j) is the element of ith row and jth column of original impedance matrix; Zmp is the impedance of link branch mp. The transfer impedance between m and n at the status of fault been cleared is the element of mth row and nth column of modified impedance matrix Z_NC, Z_mnC=Z_NC(m, n).

Controller 102 may further determine the voltage phasors of bus m and n after the fault on mp has been cleared. Initially the virtue current I_pm may be determined using EQUATION 14.

$\begin{array}{cc}{I}_{pm}=\frac{U_p-U_m}{Z_N\left(m,p\right)-Z_{mp}}& \left(\mathrm{Equation}14\right)\end{array}$

The voltage phasors of m and n can be subsequently obtained by EQUATIONS 15 and 16, based on superposition.

$\begin{array}{cc}{U}_{mC}=U_m+Z_N\left(m,m\right)I_{pm}-Z_N\left(m,p\right)I_{pm}=U_m+\frac{\left[Z_N\left(m,m\right)-Z_N\left(m,p\right)\right]}{\left[Z_N\left(m,p\right)-Z_{\mathrm{mp}}\right]}\left(U_p-U_m\right)& \left(\mathrm{Equation}15\right)\end{array}$ $\begin{array}{cc}{U}_{nC}=U_n+Z_N\left(m,n\right)I_{pm}-Z_N\left(n,p\right)I_{pm}=U_m+\frac{\left[Z_N\left(m,n\right)-Z_N\left(n,p\right)\right]}{\left[Z_N\left(m,p\right)-Z_{\mathrm{mp}}\right]}\left(U_p-U_m\right)& \left(\mathrm{Equation}16\right)\end{array}$

The power of branch mn at status of fault on mp has been cleared, using EQUATION 17.

$\begin{array}{cc}{P}_{M3}=\frac{U_{mC}{U}_{nC}}{X_{mnC}}& \left(\mathrm{Equation}17\right)\end{array}$

Controller 102 may further determine the angle of the instance that the fault was cleared, δ_CR. Assuming that the clearing time is Tc=200 ms (or another predetermined time), Eular's method may be used for calculating the fault-clearance-angle. Set a time interval for calculating the angle, for example, Δt=1 ms; and use EQUATION 18 to set initial value.

$\begin{array}{cc}\left[\begin{array}{c}\delta \left(0\right)\\ \omega \left(0\right)\end{array}\right]=\left[\begin{array}{c}{\delta }_{\mathrm{mn}0}\\ {\omega }_0\end{array}\right]& \left(\mathrm{Equation}18\right)\end{array}$

The next time steps are determined by a time-incursive EQUATION 19.

$\begin{array}{cc}\left[\begin{array}{c}\delta \left(n+1\right)\\ \omega \left(n+1\right)\end{array}\right]=\left[\begin{array}{c}\delta \left(n\right)\\ \omega \left(n\right)\end{array}\right]+\left[\begin{array}{c}\omega \left(n\right)-{\omega }_0\\ P_0-P_{M2}\sin \delta \left(n\right)\end{array}\right]\Delta t& \left(\mathrm{Equation}19\right)\end{array}$

In EQUATION 19,

$P_0=\frac{U_m{U}_n}{X_{mn}}\sin {\delta }_{mn0},$

is the transmission power in line mn. The process is repeated, until the time goes to 200 ms. Then, the fault-clearance-angle can be obtained as δ_CR=δ(200).

Controller 102 may calculate the dynamic rating power for angular stability based on equal area criterion, by solving EQUATION 20.

$\begin{array}{cc}{\int }_{{\delta }_{0\mathrm{maxi}}}^{{\delta }_{CR}}\left(P_{\mathrm{maxi}}-P_{M2}\sin \delta \right)d\delta ={\int }_{{\delta }_{\mathrm{CR}}}^{\frac{\pi }{2}-{\delta }_{0\mathrm{maxi}}}\left(P_{M3}\sin \delta -P_{\mathrm{maxi}}\right)d\delta & \left(\mathrm{Equation}20\right)\end{array}$

The equation can be simplified as follows (as we have known that P_maxi=P_M1sin δ_0maxi), into EQUATION 21.

−2P_M1δ_0maxisinδ_0maxi−P_M2cosδ_0maxi+P_M3sinδ_0maxi=(P_M3−P_M2)cosδ_CR    (Equation 21)

EQUATION 21 is an equation with respect to the only unknown variable δ_0maxi. Newton's method, or some other iterative method, may be used to solve δ_0maxi, and then the maximum rating power for angular stability can be obtained EQUATION 22.

P _maxi =P _M1 sinδ _0maxi  (Equation 22)

The controller 102 may be further configured to determine, based on the voltage phasors and the current phasors of nodes 108, a dynamic voltage stability power rating for each bus or node of the nodes 108. In embodiments, the controller 102 is configured to determine the dynamic voltage stability power rating by determining a Thevenin's circuit that is equivalent to the network 106.

FIG. 11 is a block diagram showing an equivalent Thevenin's circuit for a node ‘n’ of the power grid 104 of FIG. 1 for determining a dynamic voltage stability power rating. FIG. 12 is a block diagram showing a scheme for determining a dynamic voltage stability power rating based on the Thevenin's circuit of FIG. 11.

Controller 102 may initially determine an equivalent Thevenin's circuit at each bus or node, for example, bus n, and calculate the source impedance (Z_eq) and source voltage (E_eq) by using the impedance matrix of the network (Z_N) and using the measured voltage and current at the bus (for example, bus n in FIG. 11).

Controller 102 may determine, based on the calculated E_eq and Z_eq (in above step) and he measured reactive power Q_LD at bus n, the maximum injecting power (dynamic power rating for voltage stability, which is the maximum power that can be injected to the bus to keep the voltage no being collapse) into bus n. The scheme of FIGS. 13 to 15 may be used.

FIG. 13 is a chart showing the crossing of example reactive power curves with a consumption reactive power curve for determining maximum active power. FIG. 14 is a chart showing the relationship between active power and voltage, and the maximum active power at marginal point of voltage stability. FIG. 15 is a chart showing the determination of maximum active power based on reactive power curves.

As seen in FIGS. 13 and 14, a generation reactive power curve having two cross point with consumption reactive power curve, represents voltage stability. A single cross point is on the margin edge of the voltage stability. Therefore the maximum active power P is at the point where there is only one cross point. As seen in FIG. 15, the top point (where the only cross point it may cross with load consumption curve) is at the point,

$Q_{LD}={\left(\frac{E\cos \delta }{X}\right)}^2.$

Accordingly, the maximum injection active power which can keep the voltage stability of bus n can be calculated by EQUATION 23.

$\begin{array}{cc}{P}_{\mathrm{DynInjectMax}}=\frac{E_{eq}{V}_{\min }}{X_{eq}}\sin \delta =\frac{E_{eq}^2\cos \delta \sin \delta }{2X_{eq}^2}=\frac{E_{eq}^2}{4X_{eq}^2}\sin \left(2\delta \right)=\frac{1}{2}\sqrt{\frac{E_{eq}^2{Q}_{LD}}{X_{eq}^2}-Q_{LD}^2}& \left(\mathrm{Equation}23\right)\end{array}$

The load impedance of all the other buses except bus n may be added into the impedance matrix of the grid network Z_N form a new impedance matrix Z_Nnew of the overall system except load (or generation) connected on bus n. For example, if at bus p, there is a load impedance Z_LDp, then the renewed impedance matrix is determined using EQUATION 24.

$\begin{array}{cc}{Z}_{Nnew}\left(i,j\right)=Z_N\left(i,j\right)-\frac{Z_N\left(i,p\right)Z_N\left(p,j\right)}{Z_N\left(m,m\right)+Z_{LDp}}& \left(\mathrm{Equation}24\right)\end{array}$

The source impedance of the network at node n is the nth row and nth column of the renewed impedance matrix. That is, Z_eq=Z_Nnew(n,n) The source impedance can be determined by measured voltage U_n and the overall outgoing current I_n: E_eq=U_n+Z_egI_n.

The consuming reactive power may be determined using EQUATION 25, φ being the angle difference between the voltage and current phasors.

Q _LD =U _n I _n sin φ  (Equation 25)

The dynamic power rating may be determined using EQUATION 26.

$\begin{array}{cc}{P}_{\mathrm{DynInjectMax}}=\frac{E_{eq}{V}_{\min }}{X_{eq}}\sin \delta =\frac{E_{eq}^2\cos \delta \sin \delta }{2X_{eq}^2}=\frac{E_{eq}^2}{4X_{eq}^2}\sin \left(2\delta \right)=\frac{1}{2}\sqrt{\frac{E_{eq}^2{Q}_{LD}}{X_{eq}^2}-Q_{LD}^2}& \left(\mathrm{Equation}26\right)\end{array}$

Thus, controller 102 may thus be configured to determine each of the dynamic thermal stability power rating, the dynamic angular stability power rating, and the dynamic voltage stability power rating. The controller 102 may further control the power grid 104 in response to these parameters. For example, the controller 102 may be further configured to determine, based on the dynamic thermal stability power rating, the dynamic angular stability power rating, and the dynamic voltage stability power rating, a dynamic system rating for the power grid 104. The controller 102 may be further configured to control the power grid 104 in response to the dynamic system rating.

FIG. 16 is a block diagram showing the controller of FIG. 1 controlling the power grid to increase power transmission from a renewable energy source while reducing power transmission loss. As an example, using the scheme of FIG. 16, the following optimization may be performed, with the “Objective” of maximizing renewable output and minimizing the overall power loss, for example, being represented by EQUATION 27, and subject to constraints according to restrain equations of EQUATIONS 28 and 29, and restrain inequations of EQUATIONS 30 to 32.

max{ΣP _injectRenew −γΔP _Σ}  (Equation 27)

ΔP_i=P_si−V_iΣ^N_k=1V_k(G_ikcosδ_ik+B_iksinδ_ik)=0  (Equation 28)

ΔQ_i=Q_si−V_iΣ^N_k=1V_k(B_ikcosδ_ik−G_iksinδ_ik)=0  (Equation 29)

P _Tij <P _{DSRAng,ThermMaxi}  (Equation 30)

P _inject,j <P _DSRVolMaxi  (Equation 31)

V _imin <V<V _imax  (Equation 32)

The output of the optimization may include (1) tap-changer position of each transformer, (k_Tap); (2) the injection reactive power of each bus (Q_inject); and (3) injection active power of each bus (P_inject).

In embodiments, the controller 102 is further configured to control the power grid 104 in response to the dynamic system rating by increasing a power drawn from at least one renewable power source associated with a node of the nodes 108. In embodiments, the controller 102 is further configured to optimally control the power grid 104 in response to the dynamic system rating by maximizing a power drawn from at least one renewable power source associated with a node of the nodes 108. In embodiments, the controller 102 is configured to control the power grid 104 in response to the dynamic system rating by reducing a power loss of the network 106. In embodiments, the controller 102 is configured to optimally control the power grid 104 in response to the dynamic system rating by minimizing a power loss of the network 106. Reducing the power loss of the network 106 may include reducing at least one of a line power loss or a node power loss.

In embodiments, the present disclosure describes a method for dynamic rating of a power grid including a network of nodes coupled by lines. The method may include sending, to a controller, signals from terminal units indicative of voltage phasors and current phasors at each node of the nodes. At least one node of the nodes may be associated with a renewable power source.

The method may further include determining, by the controller, based on the voltage phasors and the current phasors of nodes, a dynamic thermal stability power rating for each line of the lines. In embodiments, the determining, by the controller, the dynamic thermal stability power rating for each line of the lines is not based on environmental parameters.

The method may further include determining, by the controller, based on the voltage phasors and the current phasors of nodes, a dynamic angular stability power rating for each node of the nodes. The determining, by the controller, the dynamic angular stability power rating may include simulating, by the controller, assuming an applied fault and clearing the fault at an adjacent branch.

The method may further include determining, by the controller, based on the voltage phasors and the current phasors of nodes, a dynamic voltage stability power rating for each node of the nodes. The determining, by the controller, the dynamic voltage stability power rating may include determining, by the controller, a Thevenin's circuit that is equivalent to the network.

The method may further include determining, by the controller, based on the dynamic thermal stability power rating, the dynamic angular stability power rating, and the dynamic voltage stability power rating, a dynamic system rating for the power grid.

The method may further include controlling, by the controller, the power grid in response to the dynamic system rating. The controlling, by the controller, the power grid may include increasing a power drawn from the renewable power source. The controlling, by the controller, the power grid may include maximizing a power drawn from the renewable power source. The controlling, by the controller, the power grid may include reducing a power loss of the network. The reducing, by the controller, the power loss of the network may include reducing a line power loss or a node power loss. The controlling, by the controller, the power grid may include increasing a ratio of power transmission to a static power rating of at least one line of the lines to greater than 1:1.

Thus, systems and methods according to the present disclosure may be used to control or regulate a power grid based on a dynamic system rating. Certain advantages of the present systems and methods may include drawing more power from renewable sources, which could significantly reduce the emission of CO₂ about 10-15%, increase the revenue of power generation by 10-15% (for example, saving non-renewable fuel use for power generation). Further advantages may include maintaining thermal, angular, and voltage stability, while increasing the ceiling of the transmission capacity dynamically (accounting for thermal stability, angular stability, and voltage stability). The transmission power may be increased compared to a power grid governed based on a static rating of power capacity. Optimizing power flow based on the dynamic power rating, may reduce or minimize the power loss during power transmission, which may can reduce operation costs by 1-3% associated with transmission loss. The present systems and methods may also increase the observability of power grids. For example, the dynamic margin of thermal, angular and voltage stability may be observed in real-time or near real-time.

Claims

1. A system for dynamic rating of a power grid comprising a network of a plurality of nodes coupled by a plurality of lines, the system comprising: a plurality of terminal units, each terminal unit of the plurality of terminal units being configured to detect a voltage phasor and a current phasor at a respective node of the plurality of nodes and generate signals indicative of the voltage phasor and the current phasor; anda controller configured to: receive signals from the plurality of terminal units indicative of the voltage phasors and the current phasors of the plurality of nodes,determine, based on the voltage phasors and the current phasors of the plurality of nodes, a dynamic thermal stability power rating for each line of the plurality of lines,determine, based on the voltage phasors and the current phasors of the plurality of nodes, a dynamic angular stability power rating for each node of the plurality of nodes,determine, based on the voltage phasors and the current phasors of the plurality of nodes, a dynamic voltage stability power rating for each node of the plurality of nodes,determine, based on the dynamic thermal stability power rating, the dynamic angular stability power rating, and the dynamic voltage stability power rating, a dynamic system rating for the power grid, andcontrol the power grid in response to the dynamic system rating.

2. The system of claim 1, wherein the terminal unit comprises a synchronous phasor measurement unit (PMU), a remote terminal unit (RTU), a feeder terminal unit (FTU), a zonal client, or a zonal agent.

3. The system of claim 1, wherein the controller comprises a power gateway or a server.

4. The system of claim 1, wherein at least one node of the plurality of nodes is associated with a renewable power source.

5. The system of claim 4, wherein the controller is configured to optimally control the power grid in response to the dynamic system rating by maximizing a power drawn from the renewable power source.

6. The system of claim 1, wherein the controller is configured to optimally control the power grid in response to the dynamic system rating by minimizing a power loss of the network.

7. The system of claim 1, wherein the controller is configured to control the power grid in response to the dynamic system rating by reducing at least one of a line power loss or a node power loss.

8. The system of claim 1, wherein the controller is configured to determine the dynamic angular stability power rating by assuming an applied fault and clearing the fault at an adjacent branch.

9. The system of claim 1, wherein the controller is configured to determine the dynamic voltage stability power rating by determining a Thevenin's circuit that is equivalent to the network.

10. The system of claim 1, wherein the controller is configured to determine the dynamic thermal stability power rating for each line of the plurality of lines without assessing environmental parameters.

11. A method for dynamic rating of a power grid comprising a plurality of nodes coupled by a plurality of lines, the method comprising: sending, to a controller, signals from a plurality of terminal units indicative of voltage phasors and current phasors at each node of the plurality of nodes,determining, by the controller, based on the voltage phasors and the current phasors of the plurality of nodes, a dynamic thermal stability power rating for each line of the plurality of lines,determining, by the controller, based on the voltage phasors and the current phasors of the plurality of nodes, a dynamic angular stability power rating for each node,determining, by the controller, based on the voltage phasors and the current phasors of the plurality of nodes, a dynamic voltage stability power rating for each node of the plurality of nodes,determining, by the controller, based on the dynamic thermal stability power rating, the dynamic angular stability power rating, and the dynamic voltage stability power rating, a dynamic system rating for the power grid, andcontrolling, by the controller, the power grid in response to the dynamic system rating.

12. The method of claim 11, wherein at least one node of the plurality of nodes is associated with a renewable power source.

13. The method of claim 12, wherein the controlling, by the controller, the power grid comprises increasing a power drawn from the renewable power source.

14. The method of claim 12, wherein the controlling, by the controller, the power grid comprises maximizing a power drawn from the renewable power source.

15. The method of claim 11, wherein the controlling, by the controller, the power grid comprises reducing a power loss of the network.

16. The method of claim 11, wherein the controlling, by the controller, the power grid comprises minimizing a power loss of the network.

17. The method of claim 11, wherein the determining, by the controller, the dynamic angular stability power rating comprises simulating, by the controller, applying a fault and clearing the fault at an adjacent branch.

18. The method of claim 11, wherein the determining, by the controller, the dynamic voltage stability power rating comprises determining, by the controller, a Thevenin's circuit that is equivalent to the network.

19. The method of claim 11, wherein the determining, by the controller, the dynamic thermal stability power rating for each line of the plurality of lines is not based on environmental parameters.

20. The method of claim 11, wherein the controlling, by the controller, the power grid comprises increasing a ratio of power transmission to a static power rating of at least one line of the plurality of lines to greater than 1.