Systems and methods for improved stability of power systems
10291025 ยท 2019-05-14
Assignee
Inventors
- Chaitanya Ashok Baone (Glenville, NY, US)
- Naresh Acharya (Schenectady, NY, US)
- Nilanjan Ray Chaudhuri (Fargo, ND, US)
Cpc classification
H02J3/10
ELECTRICITY
H02J3/001
ELECTRICITY
International classification
Abstract
The embodiments described herein provide for a system including a processor. The processor is configured to select at least one grid system contingency from a plurality of grid system contingencies. The processor is further configured to derive one or more eigen-sensitivity values based on the at least on grid system contingency. The processor is also configured to derive one or more control actions at least partially based on the eigen-sensitivity values. The processor is additionally configured to apply the one or more control actions for generation re-dispatch of a grid system.
Claims
1. A system, comprising: a processor configured to: select at least one grid system contingency from a plurality of grid system contingencies; derive one or more eigen-sensitivity values based on the at least on grid system contingency; derive one or more control actions at least partially based on the eigen-sensitivity values; apply a quadratic programming to solve an optimization problem to derive the one or more control actions, wherein the optimization problem comprises min J=.sub.i=1.sup.N.sup.
2. The system of claim 1, wherein the processor is configured to derive the one or more eigen-sensitivity values by perturbing a power of a generator system included in a plurality of generator systems via a perturbation and to distribute the perturbation among a remainder of the plurality of generator systems.
3. The system of claim 2, wherein the perturbation is distributed equally among the remainder of the plurality of generator systems.
4. The system of claim 1, wherein a first eigenvalue of the one or more eigenvalues comprises a base case .sub.base representative of a base case settling time for a specific oscillation mode of the grid system.
5. The system of claim 1, wherein a second eigenvalue of the one or more eigen/values comprises a post-contingency case .sub.post having a settling time after the grid system contingency for a specific oscillation mode of the grid system.
6. The system of claim 4, wherein the processor is configured to move the base case .sub.base so that a settling time corresponding to a post-redispatch post-contingency condition lies beyond a threshold representative of a desired damping ratio, a desired settling time, or a combination, thereof.
7. The system of claim 4, wherein a third eigenvalue of the one or more eigenvalues comprises a .sub.target and wherein the processor is configured to derive the .sub.target by using the .sub.base and a .sub.target.
8. The system of claim 7, wherein the processor is configured to derive the .sub.target by applying a .sub.post-target and a .sub.post, wherein the .sub.post-target is representative of a target settling time after the grid system contingency and wherein the .sub.post comprises a settling time after the grid system contingency.
9. The system of claim 1, comprising simulating the one or more control actions to determine if the one or more control actions are acceptable, and if acceptable, to apply the control actions for the generation re-dispatch.
10. A method, comprising: selecting, via a processor, at least one grid system contingency from a plurality of grid system contingencies; deriving, via the processor, one or more eigen-sensitivity values based on the at least on grid system contingency; deriving, via the processor, one or more control actions at least partially based on the eigen-sensitivity values; applying a quadratic pro ramming to solve an optimization problem to derive the one or more control actions, wherein the optimization problem comprises min J=.sub.i=1.sup.N.sup.
11. The method of claim 10, wherein deriving, via the processor, the one or more eigen-sensitivity values comprises perturbing a power of a generator system included in a plurality of generator systems via a perturbation and to distribute the perturbation among a remainder of the plurality of generator systems.
12. The method of claim 10, comprising creating, via the processor, a simulation model of the grid system, wherein the simulation model comprises one or more generation systems, one or more loads, and one or more power lines electrically coupling the one or more generation systems to the one or more loads, wherein deriving, via the processor, the one or more eigen-sensitivity values comprises perturbing a power of a generator system included in a plurality of generator systems via a perturbation and to distribute the perturbation among a remainder of the plurality of generator systems, and wherein perturbing the power comprises executing a simulation of the simulation model.
13. The method of claim 9, wherein deriving, via the processor, one or more eigen-sensitivity values comprises deriving, via the processor, one or more load-based eigen-sensitivity values, and wherein the one or more control actions comprises a control of one or more loads configured to receive power from the grid system.
14. A non-transitory computer readable medium storing instructions that when executed by a processor cause the processor to: select at least one grid system contingency from a plurality of grid system contingencies; derive one or more eigen-sensitivity values based on the at least on grid system contingency; derive one or more control actions at least partially based on the eigen-sensitivity values; to apply a quadratic programming to solve an optimization problem to derive the one or more control actions, wherein the optimization problem comprises min J=.sub.i=1.sup.N.sup.
15. The non-transitory computer readable medium of claim 14, wherein the instructions that cause the processor to derive the one or more eigen-sensitivity values comprise instructions that cause the processor to perturb a power of a generator system included in a plurality of generator systems via a perturbation and to distribute the perturbation among a remainder of the plurality of generator systems.
16. The non-transitory computer readable medium of claim 14, wherein the one or more eigen-sensitivity values comprise a base case .sub.base representative of a base case settling time for a specific oscillation mode of the grid system, and wherein the instructions comprise instructions that cause the processor to move the base case .sub.base so that the base case .sub.base lies beyond a threshold representative of a desired damping ratio, a desired settling time, or a combination, thereof.
17. The non-transitory computer readable medium of claim 14, comprising instructions that cause the processor to simulate the one or more control actions to determine if the one or more control actions are acceptable, and if acceptable, to apply the control actions for the generation re-dispatch.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
(1) These and other features, aspects, and advantages of the present invention will become better understood when the following detailed description is read with reference to the accompanying drawings in which like characters represent like parts throughout the drawings, wherein:
(2)
(3)
(4)
(5)
(6)
(7)
DETAILED DESCRIPTION OF THE INVENTION
(8) One or more specific embodiments of the present invention will be described below. In an effort to provide a concise description of these embodiments, all features of an actual implementation may not be described in the specification. It should be appreciated that in the development of any such actual implementation, as in any engineering or design project, numerous implementation-specific decisions must be made to achieve the developers' specific goals, such as compliance with system-related and business-related constraints, which may vary from one implementation to another. Moreover, it should be appreciated that such a development effort might be complex and time consuming, but would nevertheless be a routine undertaking of design, fabrication, and manufacture for those of ordinary skill having the benefit of this disclosure.
(9) When introducing elements of various embodiments of the present invention, the articles a, an, the, and said are intended to mean that there are one or more of the elements. The terms comprising, including, and having are intended to be inclusive and mean that there may be additional elements other than the listed elements.
(10) The present disclosure describes techniques suitable for improving oscillatory stability through generation re-dispatch. As further described herein, a sensitivity factor and optimization based approach for proactively identifying certain control elements and their optimal control action suitable for solving a power system inter-area oscillatory problem in real-time or near real-time is provided. The control elements and actions may include generator active power re-dispatch, load shedding, capacitor bank switching, and so on. In one embodiment, a sensitivity of an eigenvalue with respect to change in control variables is determined by perturbing a power system model. The perturbation may result in information suitable to be used as an input to a constrained optimization problem. The optimization problem may be configured to provide for a minimum amount of damping for the power system. The optimization problem may then be solved by using an objective, such as minimizing a total amount and cost of certain control action(s).
(11) Advantageously, the techniques described herein provide for a more proactive determination of targeted eigenvalues that may increase grid stabilization for credible what if scenarios, such as scenarios where generator systems become unavailable, scenarios where distribution lines are down, and so on. A tool, such as a software-based tool, a hardware-based tool, or combination thereof, may be created and disposed in a utility control room and will take as input solved power flow cases, associated dynamic models, and critical contingencies and derive as output a more optimal control action that may solve oscillation problems.
(12) With the foregoing in mind, it may be useful to describe an embodiment of a system, such as a power grid system 8 operatively coupled to a generation re-dispatch optimization (GRDO) system 10, as illustrated in
(13) The power generated by the power generation stations 16, 18, 20, and 22 may be transmitted through a power transmission grid 24. The power transmission grid 24 may cover a broad geographic region or regions, such as one or more municipalities, states, or countries. The transmission grid 24 may also be a single phase alternating current (AC) system, but most generally may be a three-phase AC current system. As depicted, the power transmission grid 24 may include a series of towers to support a series of overhead electrical conductors in various configurations. For example, extreme high voltage (EHV) conductors may be arranged in a three conductor bundle, having a conductor for each of three phases. The power transmission grid 24 may support nominal system voltages in the ranges of 110 kilovolts (kV) to 765 kilovolts (kV). In the depicted embodiment, the power transmission grid 24 may be electrically coupled to distribution systems (e.g., power distribution substation 26). The power distribution substation 26 may include transformers to transform the voltage of the incoming power from a transmission voltage (e.g., 765 kV, 500 kV, 345 kV, or 138 kV) to primary (e.g., 13.8 kV or 4160V) and secondary (e.g., 480V, 230V, or 120V) distribution voltages.
(14) Advanced metering infrastructure meters (e.g., smart meters) 30 may be used to monitor and communicated power related information based on electric power delivered to commercial consumers 32 and residential consumers 34. For example, the smart meters 30 may include one and/or two-way communications with the grid 8 and the GRDO 10 suitable for communicating a variety of information, including power usage, voltage, frequency, phase, power quality monitoring, and the like. The smart meters 30 may additional receive information, for example, demand response actions, time-of-use pricing information, remote service disconnects, and the like.
(15) The customers 32, 34 may operate a variety of power consuming devices 36, such as household appliances, industrial machinery, communications equipment, and the like. In certain embodiments, the power consuming devices 36 may be communicatively coupled to the grid system 8, the GRDO 10, and/or the meters 30. For example, the power consuming devices 36 may include switches that may be actuated remotely to turn on/off the devices 36 and/or to vary power consumption (e.g., lower or rise heating ventilation and air conditioning [HVAC] temperature set points). The smart meters 30 and the power consuming devices 36 may be communicatively coupled, for example, through a home area network (HAN), (for residential customers 34), wireless area network (WAN), powerline network, local area network (LAN), mesh network and the like.
(16) As mentioned earlier the GRDO 10 may be operatively coupled to the grid system 8 (smart grid) and used to more efficiently manage generation re-dispatch. The GRDO 10 may be a software system and/or a hardware system. Accordingly, the GRDO 10 may include a processor 38 suitable for executing computer code or instructions stored in the memory 40. The GRDO 10 may provide control signals to a variety of devices in the grid system 10 that can be controlled, such as generator systems, power storage systems, and loads. In one embodiment, the GRDO 10 may be disposed in utility control center 14 and may be used to observe and manage grid 8 behavior, including oscillation mode behavior. For example, the GRDO 10 may be continuously logging data from the grid 8 to adjust certain computations, such as sensitivity factors, eigenvalue-based models, constraint optimization problems, and so on to provide for a more optimal response to oscillation modes. The GRDO 10 may then provide for derivations suitable for managing generator active power re-dispatch, load shedding, capacitor bank switching, and so on. By providing for a dynamic system suitable for responding to a variety of oscillation modes, the GRDO 10 may result in a more optimal management of the grid system 8, including responses to a variety of oscillatory conditions.
(17)
(18) To group by mode shape, a right eigenvector of a system matrix (e.g., matrix representative of grid 8) may be derived, representative of a participation of certain generators (e.g., generators 104-136) in a particular mode, e.g., the mode shown in
(19)
(20) where is the eigenvalue, A.sub.sys is the system matrix, is a parameter representing the generation or load, is mode shape (i.e., right eigenvector of A.sub.sys), is the left eigenvector of A.sub.sys, and i is representative of a specific generator (e.g., generators 104-136).
(21) By applying mode shape grouping, the generators 104-136 may be depicted having arrows representative of oscillation direction. The techniques described herein may provide for eigen-sensitivity grouping, as shown in
(22) Turning now to
(23)
and/or
(24)
may thus be computed based on the eigenvalue 's real part , imaginary part , and P, where P is the perturbed power (e.g., 100MW). In another, the remainder generator systems (e.g., generator system G2 denoted by numeral 154, generator system G3 denoted by numeral 156, and generator system G4 denoted by numeral 158) may be perturbed unequally among each other, but with the total system-wide perturbation of zero.
(25) Additionally, certain contingencies, such as N1 contingencies may be identified and used during the eigen-sensitivity derivations. N1 contingencies may refer to the offlining of one system, such as a generator system or a power transmission line. For example, certain contingencies such as any one of power transmission lines A-K going down, generator system G1-G4 going down, and/or combinations thereof, may be taken into account, and power outputs adjusted to provide for desired oscillatory mode damping. First, the eigen-sensitivities S may be computed for each generator system G1-G4 (or other power source, such as capacitor banks) under post contingency conditions (i.e., after the contingency has occurred). In one embodiment, a certain equation as shown below may be used.
.sub.target=.sub.post-target.sub.postEquation (2):
(26) where .sub.post corresponds to computed post contingency (e.g., after a generator or line goes down), and .sub.post-target corresponds to a target a that is desired to be achieved post contingency. Accordingly, multiple .sub.target may be derived, each one based on certain contingencies. .sub.target may then be used to compute .sub.target. as defined below.
.sub.target=.sub.target+.sub.baseEquation (3):
(27) where .sub.base corresponds to a base case (e.g., pre contingency operations). The system 8 may then be operated to meet .sub.target. Operating the system 8 to meet the .sub.target may result in the worst case equal to .sub.post-target post contingency, which has to be set to a desired value suitable for operating at a desired damping ratio and/or post target settling time. Accordingly, the system post re-dispatch may better respond to the contingency. Indeed, the techniques described herein may be applied so that the system 8 includes a desired damping ratio and/or arrives at a desired settling time even should certain contingencies (e.g., downed lines, inoperable generators) occur.
(28) It may be useful to pictorially show the derivations described above. Accordingly,
(29) To maintain the contingency case at the desired threshold (e.g., 0.13), we would then move from the current position of 0.14 to a new position 0.17, resulting in a .sub.0 180 that may provide for the desired damping ratio (e.g., 5-8%) and/or a certain settling time (e.g., 10-40 seconds) post contingency. Since the amount of re-dispatch is determined based on the worst contingency, the remaining less severe contingencies implicitly would become stable post-dispatch.
(30) It may also be useful to provide for a specific example, say a contingency where line G-H in
(31)
may be computed for mode 1, so that a post contingency .sub.post is derived to be 0.12 and a .sub.base is derived to be 0.08. The operator of the system 8 may desire that a post contingency target be set to 0.13. Accordingly, Equation 2 leads to .sub.target=0.130.12=0.25. Then Equation 3 leads to .sub.target=0.25+(0.08)=0.33. Accordingly, it would be beneficial to manage re-dispatch to maintain 0.33 as target .
(32) Certain techniques, such as quadratic programming (QP) may be used to determine a QP optimal re-dispatch so as to provide for the target . For example, for mode 1:
min J=.sub.i=1.sup.N.sup.
(33) where min J is representative of the minimum generation re-dispatch, N.sub.g is the total number of generators, P.sub.i is the amount of re-dispatch for generator i. min J is then subject to constraints:
(34)
(35) where S.sub.i is the sensitivity of the real-part of the eigenvalue to change in power of ith generator; .sub.base is the real-part of the eigenvalue for the base case; .sub.target is the targeted or desired real-part of the eigenvalue post-dispatch; P.sub.min.sup.i is a minimum desired re-dispatch for generator i; P.sub.max.sup.i is a maximum limit for re-dispatch for generator i (typically set to 1 per unit); and .sub.min.sup.i and .sub.max.sup.i are min and max constraints, respectively, on the deviation in imaginary part of the eigenvalue for generator i.
(36) In addition to power producing devices or power storage devices (e.g., generators, capacitor banks), the techniques described herein may incorporate load changes. Indeed, any number of devices whose power can be controlled can be re-dispatched using the techniques described herein. Indeed, while a number of load points (e.g., 34, 36 in
(37) Based on the derivations above, a simulation may be executed to verify that the solution min J operates as desired. For example, a simulator that includes the system 8 or portions of the system 8 (e.g., regions 100, 102 of
(38)
(39) A set of contingencies, including N1 contingencies, may then be selected (block 204). For example, contingencies that have a higher probability of occurrence (e.g., older lines that may go down, generators that may be scheduled to go offline), contingencies that historically occur at certain times of the year, and so on, may be selected (block 204). The process 200 may then derive (block 206) eigen-sensitivities S and related derivations (e.g.,
(40)
.sub.target, .sub.post-target, .sub.post, .sub.target, .sub.base) and control actions. The control actions may include deriving min J and related derivations suitable for providing re-dispatch operations, including control signals that may increase and/or decrease power production, delivery, and consumption.
(41) The process 200 may then simulate (block 208), for example, control actions by applying the control actions to the model that was created at block 202. Simulation results may then be deemed acceptable or not acceptable (decision 210). If acceptable (decision 210), the control actions may then be applied to re-dispatch (block 212). As mentioned above, control signals may be transmitted to increase and/or decrease power production, power delivery, and power consumption. If not acceptable (decision 210), the process 200 may then enable an operator to make changes (e.g, select different eigen-sensitivity groupings, inputs to models, control actions, and so on) and iterate to block 206 so that new derivations and/or control actions may be computed. In this manner, the techniques describe herein provide for enhanced re-dispatch operations.
(42) Technical effects include executing eigen-sensitivity based re-dispatch control actions. Eigen-sensitivity groups may be selected, alternative to or in addition to mode shape groups. Eigen-sensitivity based derivations may account for contingencies, such as N1 contingencies that may occur in a grid system. Contingencies may include generator systems and/or power lines going offline. The eigen-sensitivity based derivations may provide for desired damping ratios and/or settling times responsive for given oscillatory modes. The eigen-sensitivity based derivations may also include control actions suitable for more efficient and responsive re-dispatches. The control actions may be simulated to provide for verification and validation before use. The control actions may then be implemented, for example, by transmitting control signals suitable for increasing and/or decreasing power production, power delivery, and power consumption.
(43) This written description uses examples to disclose the invention, including the best mode, and also to enable any person skilled in the art to practice the invention, including making and using any devices or systems and performing any incorporated methods. The patentable scope of the invention is defined by the claims, and may include other examples that occur to those skilled in the art. Such other examples are intended to be within the scope of the claims if they have structural elements that do not differ from the literal language of the claims, or if they include equivalent structural elements with insubstantial differences from the literal language of the claims.