# Patent

## Title
Coupling dynamics for power systems with iterative discrete decision making architectures

## Abstract

Various examples are provided that are related to coupling dynamics for, e.g., power systems with iterative discrete decision making architectures. In one example, a method includes determining an output power adjustment using a frequency difference associated with a generator of a first area of a power system and price signals corresponding to power generation in the first area and in a second area coupled to the first area by a tie-line; and providing a power command based upon the output power adjustment to a control system of the generator. In another example, a power system control system includes first and second agents configured to control power generation of a first area and a second area of a power system, respectively. The second agent can control power generation of the second area using frequency differences of generators in the second area and price signals of the first and second areas.

## Background
CROSS REFERENCE TO RELATED APPLICATIONS
This application claims priority to, and the benefit of, U.S. provisional application entitled “ANALYSIS OF COUPLING DYNAMICS FOR POWER SYSTEMS WITH ITERATIVE DISCRETE DECISION MAKING ARCHITECTURES” having Ser. No. 62/244,304, filed Oct. 21, 2015, which is hereby incorporated by reference in its entirety.


BACKGROUND
Introduction of numerous smart buildings, distributed energy sources and energy storages poses challenges in operation and control. A centralized control center may over burden its SCADA system and computing machines to collect every piece of measurements and calibrate optimal operation schemes. This may result from the large amount of data being communicated from the utility grid in addition to retail and residential users. On the other hand, due to privacy concerns, communities are not willing to share all information even though all data is assumed to be available. These conditions can complicate the operation of the centralized control center.



BRIEF DESCRIPTION OF THE DRAWINGS
Many aspects of the present disclosure can be better understood with reference to the following drawings. The components in the drawings are not necessarily to scale, emphasis instead being placed upon clearly illustrating the principles of the present disclosure. Moreover, in the drawings, like reference numerals designate corresponding parts throughout the several views.
FIG. 1 is a schematic diagram illustrating an example of a power system, in accordance with various embodiments of the present disclosure.
FIG. 2 is a block diagram illustrating an example of a generator control including discrete decision making dynamics G1 and power system dynamics G2, in accordance with various embodiments of the present disclosure.
FIGS. 3A-3B and 4A-4B illustrate examples of root loci plots of Type 1 and Type 2 systems, respectively, in accordance with various embodiments of the present disclosure.
FIGS. 5A-5B, 6A-6B and 7 illustrate examples of system dynamic responses for Type 1 consensus decision making for the power system of FIG. 1, in accordance with various embodiments of the present disclosure.
FIGS. 8A-8B, 9A-9B and 10 illustrate examples of system dynamic responses for Type 2 consensus decision making for the power system of FIG. 1, in accordance with various embodiments of the present disclosure.
FIG. 11 is an example of a system that may be utilized in power distribution analysis and control according to various embodiments of the present disclosure.



DETAILED DESCRIPTION
Disclosed herein are various embodiments related to analysis of coupling dynamics for, e.g., power systems with iterative discrete decision making architectures. Multiple control agents may be used to handle the decision making process while exchanging limited information instead of one centralized control center. Reference will now be made in detail to the description of the embodiments as illustrated in the drawings, wherein like reference numbers indicate like parts throughout the several views.
Iterative “learning” by distributed control agents can be used for power system decision making. Such decision making can achieve agreement among control agents while preserving privacy. The iterative decision making process can interact with power system dynamics. In such cases, coupled dynamics are expected. In this disclosure, a modeling approach is presented that can conduct stability analysis for these hybrid systems. In the disclosed approach, the discrete decision making process is approximated by continuous dynamics. As a result, the entire hybrid system can be represented by a continuous dynamic system. Conventional stability analysis tools can then be used to check system stability and identify key impacting factors. An example power system with multiple control agents is used to demonstrate the proposed modeling and analysis. The analysis results are then validated by nonlinear time-domain simulation. The developed continuous dynamics models shed insights into the control nature of each distributed optimization algorithm. A consensus algorithm based decision making may act as an integrator of frequency deviation. It can bring the frequency back to nominal while the primal-dual based decision making cannot.
Various privacy-preserving decision making architectures for microgrids and power systems have been examined in the literature. The mathematical foundation of the privacy-preserving decision making is distributed optimization. For example, in dual decomposition, iterative updating of the Lagrangian variable translates in an iterative “learning” process among control agents. These iterative learning processes take place at a much faster speed than the traditional hourly economic dispatch process. In many cases, a feedback loop is introduced in the discrete decision making. For example, a frequency deviation signal can be fed into a price or power command update to reflect a power imbalance. The frequency deviation is the difference between measured frequency and the nominal frequency (e.g., 60 Hz). The relationship between a frequency deviation and a power or energy imbalance has been explored in the literature for new insights. For example, it has been indicated that primary frequency control with suitable droop parameters can be aligned with a distributed primal dual algorithm based optimization procedure. The integral secondary frequency control with proper gains can be aligned with a dual accent updating procure.
Therefore, it is reasonable to suggest that the dynamics of the decision making may be coupled with the power system dynamics. The coupled dynamics may affect the design of distributed optimization schemes. Without considering the power system dynamics, a dual accent update may indicate convergence with a step size. However, this update approach may prove to be unstable in the real-world due to power system dynamics. Analysis and testing of the design considering adequate power system dynamics is important to successful implementation. Until now, there has been little research that investigates dynamic stability for such hybrid systems. This disclosure provides a modeling approach to consider a hybrid system as a whole and conduct dynamic stability analysis.
A number of approaches have been documented in the literature to tackle dynamic analysis of hybrid systems. The first one is sampled data modeling approach. The discrete process is modeled in difference equation to describe the relation between a discrete decision variable at (k+1)-th step and that variable at k-th step. Between the discrete sampling period T, the continuous system dynamics is integrated for a period. That way, the continuous dynamics can be represented by difference equations as well. The entire system can be represented by a discrete process and its stability can be judged based on the discrete system. In power systems, discrete state-space models for thyristor series controlled capacitor have been developed using sampled data approach.
The second approach is to approximate the discrete process as continuous dynamics. This approach has been adopted to describe power market dynamics. More recently, this modeling approach has been adopted to model dynamics of power supply, power demand and energy storage and develop a pricing scheme.
Power system dynamics are nonlinear and complex. Integrating power system dynamics over a period is a daunting task. Therefore, the second approach is used to approximate the discrete decision making process. Compared to the first-order aggregated frequency dynamics that have been modeled, the power system dynamics modeled here are for multi-machine systems including electromechanical dynamics, turbine-governor control and primary frequency droop. The sophisticated modeling helps accurate prediction of stability (and instability). In addition, the validation platform is based on a nonlinear power system simulation platform Power System Toolbox with adequate details to represent power system dynamics.
This disclosure will describe the discrete decision making process and the representation by continuous dynamic, present the test system, power system dynamics and the continuous dynamic system representation of the entire hybrid system. Stability analysis will be conducted and validation results presented through time-domain simulation. This disclosure presents a straightforward modeling approach for hybrid system dynamic stability analysis. Development of analytical models is based on understanding of sophisticated multi-machine power system dynamics and distributed optimization procedures. The modeling approach help sheds insights into the hybrid system dynamics. Compared to other research work on coupled market and power system dynamics, this work includes not only analysis but also validation through time domain simulation for a power system with electromechanical dynamics, turbine-governor control and primary frequency control dynamics. The validation confirms the practical value of the disclosed modeling approach.
Decision Making Process and its Continuous Dynamic Model
Two types of iterative-based distributed decision making procedures and their corresponding dynamic models are discussed here. The first type is based on primal dual decomposition and the second type is based on consensus algorithm and subgradient update.
Type 1 Primal-Dual Decomposition Based Decision Making.
Consider a two-area power system economic dispatch problem where each area has a generator and a load.

Prob1 minimize ƒ1(P1)+ƒ2(P2)  (1a)

subject to λ1:P1−P12=D1,  (1b)

λ2:P2+P12=D2,  (1c)

−d≦P12≦d,  (1d)

where ƒ1(·) and ƒ2(·) are the cost functions, P1, P2 are the power output from Gen 1 and Gen 2 and P12 is the power from Area 1 to Area 2. D1 and D2 are the loads in Area 1 and Area 2, where d is assumed to be the line limit.

The partial Lagrangian function of Prob1, without considering the limits of the tie-line, is given by:

L(P1,P2,P12,λ1,λ2)=ƒ1(P1)+λ1(D1−P1+P12+ƒ2(P2)+λ2(D2−P2−P12.  (2)

Given the two price signals (λ1,λ2), can each area determine its generation P1 and import P12? This question will be examined by looking at Area 1. For a given λ1, for this objective function ƒ1(P1)+(D1−P1+P12), if P12 has no limit imposed, there is no solution unless λ1=0. The objective function can go −∞. Therefore, this problem is considered not feasible.

Consider treating P12 differently than P1 and P2 to help solve the problem. If P12 is treated as given, similar as λ1 and λ2 are given, then the dual problem becomes:




  
    
      
        
          
            
              d
              1
            
            ⁡
            
              (
              
                
                  λ
                  1
                
                ,
                
                  P
                  12
                
              
              )
            
          
          =
          
            
              min
              
                P
                1
              
            
            ⁢
            
              {
              
                
                  
                    f
                    1
                  
                  ⁡
                  
                    (
                    
                      P
                      1
                    
                    )
                  
                
                +
                
                  
                    λ
                    1
                  
                  ⁡
                  
                    (
                    
                      
                        
                          D
                          _
                        
                        1
                      
                      -
                      
                        P
                        1
                      
                      +
                      
                        P
                        12
                      
                    
                    )
                  
                
              
              }
            
          
        
        ,
      
    
    
      
        (
        3
        )
      
    
  
  
    
      
        
          
            d
            2
          
          ⁡
          
            (
            
              
                λ
                2
              
              ,
              
                P
                12
              
            
            )
          
        
        =
        
          
            min
            
              P
              2
            
          
          ⁢
          
            
              {
              
                
                  
                    f
                    2
                  
                  ⁡
                  
                    (
                    
                      P
                      2
                    
                    )
                  
                
                +
                
                  
                    λ
                    2
                  
                  ⁡
                  
                    (
                    
                      
                        
                          D
                          _
                        
                        2
                      
                      -
                      
                        P
                        2
                      
                      +
                      
                        P
                        12
                      
                    
                    )
                  
                
              
              }
            
            .
          
        
      
    
    
      
        (
        4
        )
      
    
  




The above two problems should have solutions.

Now consider the dual problem Prob2 with P12 given by:

Prob2 maximize d1(λ1,P12)+d2(λ2,P12)  (5a)

over λ1,λ2.  (5b)

For this problem, the solution ends up dependent on P12. Let this solution of Prob2 be notated as a function of P12:g(P12). Thus, g(P12) is a dual problem of the dual problem. Since the dual problem is a concave function over λi, the dual's dual problem should be a minimization problem over P12.

Therefore, the dual's dual problem or the primal-dual problem Prob3 can be written as follows:




  
    
      
        
          
            Prob
            3
          
          ⁢
          
            
              min
              
                P
                12
              
            
            ⁢
            
              
                max
                
                  λ
                  1
                
              
              ⁢
              
                (
                
                  
                    
                      min
                      
                        P
                        1
                      
                    
                    ⁢
                    
                      {
                      
                        
                          
                            f
                            1
                          
                          ⁡
                          
                            (
                            
                              P
                              1
                            
                            )
                          
                        
                        -
                        
                          
                            λ
                            1
                          
                          ⁢
                          
                            P
                            1
                          
                        
                      
                      }
                    
                  
                  +
                  
                    
                      λ
                      1
                    
                    ⁡
                    
                      (
                      
                        
                          
                            D
                            _
                          
                          1
                        
                        +
                        
                          P
                          12
                        
                      
                      )
                    
                  
                
                )
              
            
          
        
        +
        
          
            max
            
              λ
              2
            
          
          ⁢
          
            
              (
              
                
                  
                    min
                    
                      P
                      2
                    
                  
                  ⁢
                  
                    {
                    
                      
                        
                          f
                          2
                        
                        ⁡
                        
                          (
                          
                            P
                            2
                          
                          )
                        
                      
                      -
                      
                        
                          λ
                          2
                        
                        ⁢
                        
                          P
                          2
                        
                      
                    
                    }
                  
                
                +
                
                  
                    λ
                    2
                  
                  ⁡
                  
                    (
                    
                      
                        
                          D
                          _
                        
                        2
                      
                      +
                      
                        P
                        12
                      
                    
                    )
                  
                
              
              )
            
            .
          
        
      
    
    
      
        (
        6
        )
      
    
  




A system that is connected by a tie-line can be decomposed by assuming a tie-line power flow. Each area will consider the tie-line flow injection or exporting as a negative (or positive) load. Each area carries out optimization and finds the locational marginal price (LMP) for the interfacing bus. The tie-line flow is then updated based on the price difference.

The dual's dual problem can be solved by subgradient updating of the primal variable P12. Hereinafter, this virtual tie-line flow will be notated as a different name: π. The subgradient of the line flow is (λ1−λ2). Since the primal problem is a minimization problem, in the update procedure, for a positive gradient, π should be reduced. The updating procedure can be presented as follows:

πk+1=πk−α(λ1k−λ2k),  (7)

where α>0. For a given virtual tie-line flow πk, the LMPs can be found by solving individual optimization problems for each area.

The decision making strategies utilize the assumption of lossless tie-lines. Therefore, the power dispatched by the generators only takes care of loads. The total generation is less than the total consumption including loads and tie-line power loss. To compensate for the frequency deviation (Δf), the strategies can be modified to have the price calculation include an additional component that can reflect the power unbalance or energy unbalance. The energy unbalance can be proportional to the system's average frequency deviation. Therefore, at each step, the LMPs computation becomes:

λ1k=2a1(D1+πk)+b1−KΔƒ1,  (8)

λ2k=2a2(D2−πk)+b2−KΔƒ2,  (9)

where ai; bi are coefficients of a generator quadratic cost function Ci(Pi)=aiPi2+biPi+ci and K is a gain.

If the system's frequency is below the nominal frequency, the prices will be increased. In turn, the generators will increase their dispatch.
Modeling as Continuous Dynamics:
Now, an approximate continuous model for the above mentioned iterative procedure can be given. Assuming that Δƒ1=Δƒ2, the iteration of the virtual tie-line flow can be given by:

πk+1=(1−2α(a1+a2))πk+π0,  (10)

where π0=−α(2a1D1−2a2D2+b1−b2).

The assumption that frequency everywhere is the same is true for steady-state. During transients, this statement is not true. However, discrete decision making takes several seconds for each step. For this interval, a power system may be able to reach steady state. Even if it has not achieved steady-state, frequencies represented by generator speeds agree with each other in general.
The power references can be determined by the prices. Therefore,




  
    
      
        
          P
          1
          k
        
        =
        
          
            (
            
              
                π
                k
              
              +
              
                D
                1
              
            
            )
          
          -
          
            
              K
              
                2
                ⁢
                
                  a
                  1
                
              
            
            ⁢
            Δ
            ⁢
            
              
            
            ⁢
            
              
                f
                k
              
              .
            
          
        
      
    
    
      
        (
        11
        )
      
    
  




Considering that the frequency measurement of the previous step is taken in the price calculation, the k+1 step power reference can be modified as:





  
    
      
        
          P
          1
          
            k
            +
            1
          
        
        =
        
          
            (
            
              
                π
                
                  k
                  +
                  1
                
              
              +
              
                D
                1
              
            
            )
          
          -
          
            
              K
              
                2
                ⁢
                
                  a
                  1
                
              
            
            ⁢
            Δ
            ⁢
            
              
            
            ⁢
            
              
                f
                k
              
              .
            
          
        
      
    
    
      
        (
        12
        )
      
    
  




Substituting πk+1 and πk by P1k+1 and P1k, it can be found that:





  
    
      
        
          
            P
            1
            
              k
              +
              1
            
          
          =
          
            
              
                
                
                  1
                  -
                  
                    2
                    ⁢
                    
                      α
                      ⁡
                      
                        (
                        
                          
                            a
                            1
                          
                          +
                          
                            a
                            2
                          
                        
                        )
                      
                    
                  
                
                ]
              
              ⁢
              
                P
                k
              
            
            -
            
              K
              ⁢
              
                
              
              ⁢
              
                
                  α
                  ⁡
                  
                    (
                    
                      
                        a
                        1
                      
                      +
                      
                        a
                        2
                      
                    
                    )
                  
                
                
                  a
                  1
                
              
              ⁢
              Δ
              ⁢
              
                
              
              ⁢
              
                f
                k
              
            
            +
            
              P
              10
            
          
        
        ,
      
    
    
      
        (
        13
        )
      
    
  




where P10=2α(a1+a2)D1.

Using the forward Euler method, the derivative at k-step can be expressed as:




  
    
      
        
          
            
              P
              .
            
            1
            k
          
          ≈
          
            
              
                P
                1
                
                  k
                  +
                  1
                
              
              -
              
                P
                1
                k
              
            
            τ
          
        
        ,
      
    
    
      
        (
        14
        )
      
    
  




where τ is the step size. Therefore, the discrete equation can be approximated by a continuous dynamic equation, such as:





  
    
      
        
          τ
          ⁢
          
            
          
          ⁢
          
            
              P
              .
            
            1
          
        
        =
        
          
            
              -
              2
            
            ⁢
            
              α
              ⁡
              
                (
                
                  
                    a
                    1
                  
                  +
                  
                    a
                    2
                  
                
                )
              
            
            ⁢
            P
          
          -
          
            K
            ⁢
            
              
                α
                ⁡
                
                  (
                  
                    
                      a
                      1
                    
                    +
                    
                      a
                      2
                    
                  
                  )
                
              
              
                a
                1
              
            
            ⁢
            Δ
            ⁢
            
              
            
            ⁢
            f
          
          +
          
            
              P
              10
            
            .
          
        
      
    
    
      
        (
        15
        )
      
    
  




To this end, the continuous dynamic model can be derived for the discrete decision making process. In the Laplace domain, the power command can be expressed as:




  
    
      
        
          
            Δ
            ⁢
            
              
            
            ⁢
            
              P
              i
              ref
            
          
          =
          
            
              K
              
                2
                ⁢
                
                  a
                  i
                
              
            
            ⁢
            
              1
              
                1
                +
                
                  
                    τ
                    ′
                  
                  ⁢
                  s
                
              
            
            ⁢
            Δ
            ⁢
            
              
            
            ⁢
            
              f
              i
            
          
        
        ,
      
    
    
      
        (
        16
        )
      
    
  




where





  
    τ
    ′
  
  =
  
    
      τ
      
        2
        ⁢
        
          
        
        ⁢
        
          α
          ⁡
          
            (
            
              
                a
                1
              
              +
              
                a
                2
              
            
            )
          
        
      
    
    .
  




Type 2 Consensus Algorithm and Subgradient Update Based Decision Making.
The second type of iterative based decision making is based on a consensus algorithm and subgradient update. A consensus problem can be identified from the original economic dispatch problem. For the following two-area system, the original economic dispatch problem can be expressed as follows:

minimize C1(P1)+C2(P2)  (17a)

subject to:P1+P1=D1+D2  (17b)

The dual problem can then be described by:





  
    
      
        
          
            max
            λ
          
          ⁢
          
            
              min
              
                
                  P
                  1
                
                ,
                
                  P
                  2
                
              
            
            ⁢
            
              
                C
                1
              
              ⁡
              
                (
                
                  P
                  1
                
                )
              
            
          
        
        +
        
          
            C
            2
          
          ⁡
          
            (
            
              P
              2
            
            )
          
        
        +
        
          
            λ
            ⁡
            
              (
              
                
                  D
                  1
                
                -
                
                  P
                  1
                
                +
                
                  D
                  2
                
                -
                
                  P
                  2
                
              
              )
            
          
          .
        
      
    
    
      
        (
        18
        )
      
    
  




The above problem can be converted to a consensus problem by introducing price signals λ1 and λ2 for each area. Here, λ1 should be equal to λ2. Therefore, the optimization problem can be converted to a maximization problem with a consensus constraint.
The consensus algorithm that utilizes a stochastic matrix to conduct weighted averaging only guarantees consensus of multiple λ. It cannot guarantee that the λ can maximize the dual problem's objective function. In order to guarantee maximization, a subgradient update is used. The subgradient of λ is the total power unbalance, which is based on global information.
Fortunately in power systems, frequency deviation is a measure of power unbalance. Frequency is a local measurement. Therefore, distributed control can be realized by substituting the subgradient of λ by the frequency deviation. The iterative procedure can be described by the following equations:




  
    
      
        
          
            
              [
              
                
                  
                    
                      λ
                      1
                    
                  
                
                
                  
                    
                      λ
                      2
                    
                  
                
              
              ]
            
            
              k
              +
              1
            
          
          =
          
            
              
                A
                ⁡
                
                  [
                  
                    
                      
                        
                          λ
                          1
                        
                      
                    
                    
                      
                        
                          λ
                          2
                        
                      
                    
                  
                  ]
                
              
              k
            
            -
            
              K
              ⁡
              
                [
                
                  
                    
                      
                        Δ
                        ⁢
                        
                          
                        
                        ⁢
                        
                          f
                          1
                        
                      
                    
                  
                  
                    
                      
                        Δ
                        ⁢
                        
                          
                        
                        ⁢
                        
                          f
                          2
                        
                      
                    
                  
                
                ]
              
            
          
        
        ,
      
    
    
      
        (
        19
        )
      
    
  




where A is a stochastic matrix and where K is the gain in the frequency deviation feedback. For the test two-area system, select





  A
  =
  
    
      [
      
        
          
            
              1
              /
              2
            
          
          
            
              1
              /
              2
            
          
        
        
          
            
              1
              /
              2
            
          
          
            
              1
              /
              2
            
          
        
      
      ]
    
    .
  




The decision making again introduces feedback signals of frequencies. For each area, the power command is related to the Lagrangian multiplier. By ignoring the limits of each generator, it can be found that:

λ1=2a1P1+b1,  (20)

λ2=2a2P2+b2.  (21)

The iteration procedure of equation (19) can now be expressed in terms of the power commands:




  
    
      
        
          
            
              [
              
                
                  
                    
                      P
                      1
                    
                  
                
                
                  
                    
                      P
                      2
                    
                  
                
              
              ]
            
            
              k
              +
              1
            
          
          +
          
            [
            
              
                
                  
                    
                      b
                      1
                    
                    
                      2
                      ⁢
                      
                        a
                        1
                      
                    
                  
                
              
              
                
                  
                    
                      b
                      2
                    
                    
                      2
                      ⁢
                      
                        a
                        2
                      
                    
                  
                
              
            
            ]
          
        
        =
        
          
            A
            ⁡
            
              (
              
                
                  
                    [
                    
                      
                        
                          
                            P
                            1
                          
                        
                      
                      
                        
                          
                            P
                            2
                          
                        
                      
                    
                    ]
                  
                  k
                
                +
                
                  [
                  
                    
                      
                        
                          
                            b
                            1
                          
                          
                            2
                            ⁢
                            
                              a
                              1
                            
                          
                        
                      
                    
                    
                      
                        
                          
                            b
                            2
                          
                          
                            2
                            ⁢
                            
                              a
                              2
                            
                          
                        
                      
                    
                  
                  ]
                
              
              )
            
          
          -
          
            
              K
              ⁡
              
                [
                
                  
                    
                      
                        
                          Δ
                          ⁢
                          
                            
                          
                          ⁢
                          
                            f
                            1
                          
                        
                        
                          2
                          ⁢
                          
                            a
                            1
                          
                        
                      
                    
                  
                  
                    
                      
                        
                          Δ
                          ⁢
                          
                            
                          
                          ⁢
                          
                            f
                            2
                          
                        
                        
                          2
                          ⁢
                          
                            a
                            2
                          
                        
                      
                    
                  
                
                ]
              
            
            .
          
        
      
    
    
      
        (
        22
        )
      
    
  




The difference equation can now be converted to a continuous dynamic equation given by:





  
    
      
        
          [
          
            
              
                
                  Δ
                  ⁢
                  
                    
                  
                  ⁢
                  
                    P
                    1
                  
                
              
            
            
              
                
                  Δ
                  ⁢
                  
                    
                  
                  ⁢
                  
                    P
                    2
                  
                
              
            
          
          ]
        
        =
        
          -
          
            
              
                
                  
                    
                      (
                      
                        
                          τ
                          ⁢
                          
                            
                          
                          ⁢
                          s
                        
                        -
                        A
                        +
                        I
                      
                      )
                    
                    
                      -
                      1
                    
                  
                  ⁡
                  
                    [
                    
                      
                        
                          
                            K
                            
                              2
                              ⁢
                              
                                a
                                1
                              
                            
                          
                        
                        
                          0
                        
                      
                      
                        
                          0
                        
                        
                          
                            K
                            
                              2
                              ⁢
                              
                                a
                                2
                              
                            
                          
                        
                      
                    
                    ]
                  
                
                
                  }
                  
                    
                      G
                      1
                    
                    ⁡
                    
                      (
                      s
                      )
                    
                  
                
              
              ⁡
              
                [
                
                  
                    
                      
                        Δ
                        ⁢
                        
                          
                        
                        ⁢
                        
                          f
                          1
                        
                      
                    
                  
                  
                    
                      
                        Δ
                        ⁢
                        
                          
                        
                        ⁢
                        
                          f
                          2
                        
                      
                    
                  
                
                ]
              
            
            .
          
        
      
    
    
      
        (
        23
        )
      
    
  




The gain matrix G1(s) defines the transfer function matrix from the frequency deviation to the power commands as:





  
    
      
        
          
            G
            1
          
          ⁡
          
            (
            s
            )
          
        
        =
        
          
            
              K
              
                2
                ⁢
                τ
                ⁢
                
                  
                
                ⁢
                
                  s
                  ⁡
                  
                    (
                    
                      
                        τ
                        ⁢
                        
                          
                        
                        ⁢
                        s
                      
                      +
                      1
                    
                    )
                  
                
              
            
            ⁡
            
              [
              
                
                  
                    
                      
                        
                          2
                          ⁢
                          
                            
                          
                          ⁢
                          τ
                          ⁢
                          
                            
                          
                          ⁢
                          s
                        
                        +
                        1
                      
                      
                        2
                        ⁢
                        
                          a
                          1
                        
                      
                    
                  
                  
                    
                      1
                      
                        2
                        ⁢
                        
                          a
                          2
                        
                      
                    
                  
                
                
                  
                    
                      1
                      
                        2
                        ⁢
                        
                          a
                          1
                        
                      
                    
                  
                  
                    
                      
                        
                          2
                          ⁢
                          τ
                          ⁢
                          
                            
                          
                          ⁢
                          s
                        
                        +
                        1
                      
                      
                        2
                        ⁢
                        
                          a
                          2
                        
                      
                    
                  
                
              
              ]
            
          
          .
        
      
    
    
      
        (
        24
        )
      
    
  




Converting the discrete decision making process to continuous dynamics sheds light into each algorithm. Compared to the two dynamics of the decision making algorithms, the consensus algorithm includes an integrator unit. It is expected that consensus algorithm based Type-2 decision making can bring the frequency deviation to zero. The primal-dual algorithm is similar to a first-order filter. Therefore, it is not expected that the Type-1 decision making can bring the frequency back to nominal.
Test System and Power System Dynamic Model
Next, the power system dynamics model and the integrated system model will be described and analyzed. Referring to FIG. 1, shown is an example of a test power system comprising a two-area four-machine system 100. The two-area system of FIG. 1 illustrates the physical topology and the Type-1 information exchange architecture. This system comes from the classical two-area four machine power system with the following modifications: the tie-line 103 has been shortened; the inertia constants of the machines 106 have been reduced to 2.5 pu to provide faster electromechanical dynamics; and the damping coefficients have been set to 1 pu. The generators 106 are modeled as classical generators with turbine governor blocks. Primary frequency droops with the regulation constant at 4% are all included.
The underlying power system dynamic model Δf/ΔP1ref can be found. The two generators 106 in each area are coherent and therefore can be considered as one generator 106. The two-area four-machine system 100 may now be represented by a two-generator system. The two rotor angles can be expressed as:




  
    
      
        
          
            Δ
            ⁢
            
              
            
            ⁢
            
              δ
              1
            
          
          =
          
            
              1
              
                
                  
                    M
                    1
                  
                  ⁢
                  
                    s
                    2
                  
                
                +
                
                  
                    D
                    1
                  
                  ⁢
                  s
                
                +
                
                  T
                  1
                
              
            
            ⁢
            
              (
              
                
                  Δ
                  ⁢
                  
                    
                  
                  ⁢
                  
                    P
                    
                      m
                      ⁢
                      
                        
                      
                      ⁢
                      1
                    
                  
                
                +
                
                  
                    T
                    1
                  
                  ⁢
                  Δ
                  ⁢
                  
                    
                  
                  ⁢
                  
                    δ
                    2
                  
                
              
              )
            
          
        
        ,
      
    
    
      
        (
        25
        )
      
    
  
  
    
      
        
          
            Δ
            ⁢
            
              
            
            ⁢
            
              δ
              2
            
          
          =
          
            
              1
              
                
                  
                    M
                    2
                  
                  ⁢
                  
                    s
                    2
                  
                
                +
                
                  
                    D
                    2
                  
                  ⁢
                  s
                
                +
                
                  T
                  2
                
              
            
            ⁢
            
              (
              
                
                  Δ
                  ⁢
                  
                    
                  
                  ⁢
                  
                    P
                    
                      m
                      ⁢
                      
                        
                      
                      ⁢
                      2
                    
                  
                
                +
                
                  
                    T
                    2
                  
                  ⁢
                  Δ
                  ⁢
                  
                    
                  
                  ⁢
                  
                    δ
                    1
                  
                
              
              )
            
          
        
        ,
      
    
    
      
        (
        26
        )
      
    
  




where M, D, and T are inertia constants, damping and synchronizing coefficients, respectively. T1=T2. Rearranging the equations, gives:





  
    
      
        
          Δδ
          1
        
        =
        
          
            
              
                
                  (
                  
                    
                      
                        M
                        1
                      
                      ⁢
                      
                        s
                        2
                      
                    
                    +
                    
                      
                        D
                        1
                      
                      ⁢
                      s
                    
                    +
                    
                      T
                      1
                    
                  
                  )
                
                ⁢
                Δ
                ⁢
                
                  
                
                ⁢
                
                  P
                  
                    m
                    ⁢
                    
                      
                    
                    ⁢
                    1
                  
                
              
              +
              
                
                  T
                  1
                
                ⁢
                Δ
                ⁢
                
                  
                
                ⁢
                
                  P
                  
                    m
                    ⁢
                    
                      
                    
                    ⁢
                    2
                  
                
              
            
            
              
                
                  (
                  
                    
                      
                        M
                        1
                      
                      ⁢
                      
                        s
                        2
                      
                    
                    +
                    
                      
                        D
                        1
                      
                      ⁢
                      s
                    
                    +
                    
                      T
                      1
                    
                  
                  )
                
                ⁢
                
                  (
                  
                    
                      
                        M
                        1
                      
                      ⁢
                      
                        s
                        2
                      
                    
                    +
                    
                      
                        D
                        1
                      
                      ⁢
                      s
                    
                    +
                    
                      T
                      1
                    
                  
                  )
                
              
              -
              
                
                  T
                  1
                
                ⁢
                
                  T
                  2
                
              
            
          
          .
        
      
    
    
      
        (
        27
        )
      
    
  




Considering the turbine-governor control block Gtg, given by:




  
    
      
        
          
            [
            
              
                
                  
                    Δ
                    ⁢
                    
                      
                    
                    ⁢
                    
                      f
                      1
                    
                  
                
              
              
                
                  
                    Δ
                    ⁢
                    
                      
                    
                    ⁢
                    
                      f
                      2
                    
                  
                
              
            
            ]
          
          =
          
            
              
                G
                2
              
              ⁡
              
                (
                s
                )
              
            
            ⁡
            
              [
              
                
                  
                    
                      Δ
                      ⁢
                      
                        
                      
                      ⁢
                      
                        P
                        1
                        ref
                      
                    
                  
                
                
                  
                    
                      Δ
                      ⁢
                      
                        
                      
                      ⁢
                      
                        P
                        2
                        ref
                      
                    
                  
                
              
              ]
            
          
        
        ,
      
    
    
      
        (
        28
        )
      
    
  




the transfer function matrix G2(s) that defines the relationship from the power command to the speed deviations due to the power system dynamics is expressed as:





  
    
      
        
          
            
              G
              2
            
            ⁡
            
              (
              s
              )
            
          
          =
          
            
              
                
                  
                    G
                    tg
                  
                  ⁡
                  
                    (
                    s
                    )
                  
                
                ⁢
                s
              
              
                
                  ω
                  0
                
                ⁡
                
                  [
                  
                    
                      
                        (
                        
                          
                            
                              M
                              1
                            
                            ⁢
                            
                              s
                              2
                            
                          
                          +
                          
                            
                              D
                              1
                            
                            ⁢
                            s
                          
                          +
                          
                            T
                            1
                          
                        
                        )
                      
                      ⁢
                      
                        (
                        
                          
                            
                              M
                              2
                            
                            ⁢
                            
                              s
                              2
                            
                          
                          +
                          
                            
                              D
                              2
                            
                            ⁢
                            s
                          
                          +
                          
                            T
                            2
                          
                        
                        )
                      
                    
                    -
                    
                      
                        T
                        1
                      
                      ⁢
                      
                        T
                        2
                      
                    
                  
                  ]
                
              
            
            ⁡
            
              [
              
                M
                
                  G
                  2
                
              
              ]
            
          
        
        ,
      
    
    
      
        (
        29
        )
      
    
  




where





  
    [
    
      M
      
        G
        2
      
    
    ]
  
  =
  
    
      [
      
        
          
            
              
                
                  M
                  1
                
                ⁢
                
                  s
                  2
                
              
              +
              
                
                  D
                  1
                
                ⁢
                s
              
              +
              
                T
                1
              
            
          
          
            
              T
              1
            
          
        
        
          
            
              T
              2
            
          
          
            
              
                
                  M
                  2
                
                ⁢
                
                  s
                  2
                
              
              +
              
                
                  D
                  2
                
                ⁢
                s
              
              +
              
                T
                2
              
            
          
        
      
      ]
    
    .
  




In addition, the droop control can be included. Therefore, the diagonal components of G2(s) have been modified to include the droop control. Referring to FIG. 2, shown is a block diagram illustrating an example of the entire system, where G1 203 represents the discrete decision making dynamics while G2 206 represents the power system dynamics.

Root Loci:
To examine the stability of the closed-loop system, the open-loop gain matrix G2G1 will be examined. G2G1 is a two by two matrix. The root loci of the first row first column element are shown in FIGS. 3A-3B and 4A-4B. FIGS. 3A and 3B (an enlarged view of FIG. 3A) illustrate the root loci for the Type 1 system with two step sizes τ=2 and τ=20. It can be shown that the droop related poles will move to the right-half-plane (RHP) when the gain K is increasing. The system can be made more stable by increasing the step size τ. FIGS. 4A and 4B (enlarged views of FIG. 4A) illustrate the root loci for the Type 2 system with two step sizes τ=2 and τ=20. It can be shown that droop related poles will go to the right-half-plane (RHP) when the gain K is increasing. Again, the system can be made more stable by increasing the step size τ. The analysis shows that the hybrid system could suffer low frequency oscillations of approximately 0.1-0.2 Hz. Increasing the gain of the frequency deviation in the discrete decision making steps will make the system go unstable.
Dynamic Simulation Results
Next, dynamic simulation results are presented to validate the system operation. The two types of discrete decision making procedures were implemented using the two-area four-machine power system 100 shown in FIG. 1 and Power System Toolbox was selected as the dynamic simulation platform.
The power system and Type 1 decision making architecture are shown in FIG. 1. The discrete decision making can take place at a predefined interval, e.g., every 2 seconds or every 5 seconds. The power commands from Agent 1 (109a) and Agent 2 (109b) will be sent to change the turbine-governors' power reference inputs of the generator controller(s). Among the two agents 109, the information exchanged includes the virtual tie-line power flow (π) and the price signal (λ2). Area 1 comprises Gen 1 (106a) and Gen 2 (106b) and Load 1 (112a). Area 2 comprises Gen 3 (106c), Gen 4 (106d) and Load 2 (112b). The two areas are connected through tie-lines 103.
Initially, the four generators 106a-106d are dispatched at 7.0207 pu, 7.00 pu, 7.16 pu and 7.00 pu, respectively. Assume that in Area 1 the two generators 106a and 106b have the same quadratic cost functions: 1.5P12 and 1.5P22 and in Area 2 the two generators 106c and 106d also have the same quadratic cost functions P32 and P42. In the example of FIG. 1, the total load of the two areas is 27.41 pu. Initially, the dispatch levels of the four generators 106 are similar. After the decision making procedures, the Area 2 generators will have higher dispatch levels as Gen 3 (106c) and Gen 4 (106d) are much cheaper than Gen 1 (106a) and Gen 2 (106b).
Type 1 Primal-Dual Based Decision Making.
Three scenarios were compared to show the effect of the step size τ (sec) of the discrete decision making and the gain K in frequency deviation feedback.

    
    
        Scenario 1: τ=2 and K=300.
        Scenario 2: τ=2 and K=500.
        Scenario 3: τ=5 and K=500.
    
    


FIG. 5A illustrates the system dynamic responses for Scenario 1 with Type 1 consensus decision making, where plot 503 shows the generator speeds in pu, plot 506 shows the generator power based on a system power base of 100 MW, plot 509 shows the turbine governor unit power based on the machine power base of 900 MW, and plot 512 shows the angles of Gens 1-3 (160a-106c of FIG. 1) relative to Gen 4 (106d of FIG. 1) in radians. FIG. 5B illustrates changes in the Lagrangian multipliers over time, where plot 515 shows the virtual tie-line power flow (π), plot 518 shows the price signals (λ1, λ2), and plot 521 shows the frequency deviations (Δf).
FIG. 6A illustrates the system dynamic responses for Scenario 2 with Type 1 consensus decision making, where plot 603 shows the generator speeds in pu, plot 606 shows the generator power based on a system power base of 100 MW, plot 609 shows the turbine governor unit power based on the machine power base of 900 MW, and plot 612 shows the angles of Gens 1-3 (160a-106c of FIG. 1) relative to Gen 4 (106d of FIG. 1) in radians. FIG. 6B illustrates changes in the Lagrangian multipliers over time, where plot 615 shows the virtual tie-line power flow (π), plot 618 shows the price signals (λ1, λ2), and plot 621 shows the frequency deviations (Δf).
Oscillations at 0.2 Hz were observed in Scenario 2 when the gain was increased to K=500. In Scenario 3, the step size z was increased to 5 seconds for the gain K of 500. The oscillations are then damped. The machine speeds for Gen 1 (106a of FIG. 1) for the three scenarios with Type 1 consensus decision making are compared in FIG. 7. The simulation results corroborate the findings made in the root locus analysis. The slower the discrete decision making process, the more stable the system becomes.
Type 2 Consensus Based Decision Making.
Two scenarios were compared to show the effect of the step size τ (sec) of discrete decision making and the gain K in frequency deviation feedback.

    
    
        Scenario 1: τ=2 and K=500.
        Scenario 2: τ=5 and K=500.
    
    


FIG. 8A illustrates the system dynamic responses for Scenario 1 with Type 2 consensus decision making, where plot 803 shows the generator speeds in pu, plot 806 shows the generator power based on a system power base of 100 MW, plot 809 shows the turbine governor unit power based on the machine power base of 900 MW, and plot 812 shows the angles of Gens 1-3 (160a-106c of FIG. 1) relative to Gen 4 (106d of FIG. 1) in radians. FIG. 8B illustrates changes in the Lagrangian multipliers over time, where plot 815 shows the price signals (λ1, λ2), and plot 818 shows the frequency deviations (Δf).
FIG. 9A illustrates the system dynamic responses for Scenario 2 with Type 2 consensus decision making, where plot 903 shows the generator speeds in pu, plot 906 shows the generator power based on a system power base of 100 MW, plot 909 shows the turbine governor unit power based on the machine power base of 900 MW, and plot 912 shows the angles of Gens 1-3 (160a-106c of FIG. 1) relative to Gen 4 (106d of FIG. 1) in radians. FIG. 9B illustrates changes in the Lagrangian multipliers over time, where plot 915 shows the price signals (λ1, λ2), and plot 918 shows the frequency deviations (Δf).
The comparison of the two scenarios is presented in FIG. 10. It can be observed in FIG. 10 that when τ=2 seconds, 0.15 Hz oscillations were observed. When the step size was increased to 5 seconds, the oscillations exhibited better damping. The dynamic simulation results corroborated the findings made through linear system analysis above. The slower the discrete decision making, more stable the system becomes.
Comparing the frequency response of Type-1 and Type-2 architectures, we also confirm this important finding: the particular consensus algorithm works as a secondary frequency control with economic dispatch. Type-2 decision making process can bring frequency back to the nominal frequency.
In this disclosure, continuous dynamic models for iterative decision making processes are presented. The models were used together with a power system dynamic model to determine the hybrid system dynamic stability. Such stability issues cannot be identified when one of the dynamics is not considered. The continuous dynamic model derivation step and linear analysis of the integrated power system and decision making system was demonstrated. The analysis identified the stability issue for these types of hybrid systems. The closed-loop system poles due to turbine governor, primary frequency control and the decision making dynamics will move to the right half plane when the frequency deviation gain is increased. A slower decision making process leads to a more stable system. Time-domain simulation in PST was been conducted for validation.
Referring now to FIG. 11, shown is an example of a system 1000 that may be utilized as an agent 109 (FIG. 1) in the monitoring and control of power systems. The system 1000 includes one or more computing device(s) 1003, one or more distribution monitoring device(s) 1005 that can be distributed about a power distribution network to provide indications of operational conditions in the power distribution network, and one or more generator controller(s) 1007 such as a turbine-governor control for the generators 106 of FIG. 1. For example, the distribution monitoring device(s) 1005 can monitor electrical power distribution and operation by monitoring conditions of components of the distribution network such as generators and/or distribution buses. Monitored conditions can include voltage levels, current flows, frequencies, operating speeds, and/or other conditions of the various system components. The generator controller(s) 1007 can control the operation of the corresponding generator(s) 106 (FIG. 1) based upon signals received from the computing device 1003. While not shown in FIG. 11, the computing device 1003 can include communication interfaces allowing the computing device 1003 to communicate with another computing device 1003, either directly or through a network, to facilitate agent-to-agent communications.
The computing device 1003 includes at least one processor circuit, for example, having a processor 1009 and a memory 1012, both of which are coupled to a local interface 1015. To this end, the computing device(s) 1003 may comprise, for example, a server computer or any other system providing computing capability. The computing device(s) 1003 may include, for example, one or more display devices such as cathode ray tubes (CRTs), liquid crystal display (LCD) screens, gas plasma-based flat panel displays, LCD projectors, or other types of display devices, etc. The computing device(s) 1003 may also include, for example various peripheral devices. In particular, the peripheral devices may include input devices such as, for example, a keyboard, keypad, touch pad, touch screen, microphone, scanner, mouse, joystick, or one or more push buttons, etc. Even though the computing device 1003 is referred to in the singular, it is understood that a plurality of computing devices 1003 may be employed in the various arrangements as described above. The local interface 1015 may comprise, for example, a data bus with an accompanying address/control bus or other bus structure as can be appreciated.
Stored in the memory 1012 are both data and several components that are executable by the processor 1009. In particular, stored in the memory 1012 and executable by the processor 1009 include an agent application 1018 and potentially other applications. Also stored in the memory 1012 may be a data store 1021 and other data. The data stored in the data store 1021, for example, is associated with the operation of the various applications and/or functional entities described below. For example, the data store may include sample analysis results, corrective measures, and other data or information as can be understood. In addition, an operating system 1024 may be stored in the memory 1012 and executable by the processor 1009. The data store 1021 may be may be located in a single computing device or may be dispersed among many different devices.
The distribution monitoring device 1005 is representative of a plurality of devices that may be communicatively coupled to the computing device 1003 through a network 1027 such as, e.g., the Internet, intranets, extranets, wide area networks (WANs), local area networks (LANs), wired networks, wireless networks, networks configured for communication over a power grid, or other suitable networks, etc., or any combination of two or more such networks. The distribution monitoring device 1005 may comprise, for example, a processor-based system such as a computer system or other application specific monitoring system with communication capabilities. In some embodiments, a distribution monitoring device 1005 may be directly connected to the computing device 1003.
The components executed on the computing device 1003 include, for example, the agent application 1018 and other systems, applications, services, processes, engines, or functionality not discussed in detail herein. The computing device 1003 can receive information regarding the monitored component from a distribution monitoring device 1005, which can then be evaluated by the agent application 1018. The computing device 1003 may also receive stored information regarding a power distribution network for modeling and evaluation of the network.
It is understood that there may be other applications that are stored in the memory 1012 and are executable by the processor 1009 as can be appreciated. Where any component discussed herein is implemented in the form of software, any one of a number of programming languages may be employed such as, for example, C, C++, C#, Objective C, Java, Java Script, Perl, PHP, Visual Basic, Python, Ruby, Delphi, Flash, or other programming languages.
A number of software components are stored in the memory 1012 and are executable by the processor 1009. In this respect, the term “executable” means a program file that is in a form that can ultimately be run by the processor 1009. Examples of executable programs may be, for example, a compiled program that can be translated into machine instructions in a format that can be loaded into a random access portion of the memory 1012 and run by the processor 1009, source code that may be expressed in proper format such as object code that is capable of being loaded into a random access portion of the memory 1012 and executed by the processor 1009, or source code that may be interpreted by another executable program to generate instructions in a random access portion of the memory 1012 to be executed by the processor 1009, etc. An executable program may be stored in any portion or component of the memory 1012 including, for example, random access memory (RAM), read-only memory (ROM), hard drive, solid-state drive, USB flash drive, memory card, optical disc such as compact disc (CD) or digital versatile disc (DVD), floppy disk, magnetic tape, or other memory components.
Also, the processor 1009 may represent multiple processors 1009 and the memory 1012 may represent multiple memories 1012 that operate in parallel processing circuits, respectively. In such a case, the local interface 1015 may be an appropriate network that facilitates communication between any two of the multiple processors 1009, between any processor 1009 and any of the memories 1012, or between any two of the memories 1012, etc. The local interface 1015 may comprise additional systems designed to coordinate this communication, including, for example, performing load balancing. The processor 1009 may be of electrical or of some other available construction.
Although the agent application 1018, and other various systems described herein, may be embodied in software or instructions executed by general purpose hardware as discussed above, as an alternative the same may also be embodied in dedicated hardware or a combination of software/general purpose hardware and dedicated hardware. If embodied in dedicated hardware, each can be implemented as a circuit or state machine that employs any one of or a combination of a number of technologies. These technologies may include, but are not limited to, discrete logic circuits having logic gates for implementing various logic functions upon an application of one or more data signals, application specific integrated circuits having appropriate logic gates, or other components, etc. Such technologies are generally well known by those skilled in the art and, consequently, are not described in detail herein.
Any logic or application described herein, including the ED/FR feedback control application 1018, that comprises software or instructions can be embodied in any non-transitory computer-readable medium for use by or in connection with an instruction execution system such as, for example, a processor 1009 in a computer system or other system. In this sense, the logic may comprise, for example, statements including instructions and declarations that can be fetched from the computer-readable medium and executed by the instruction execution system. In the context of the present disclosure, a “computer-readable medium” can be any medium that can contain, store, or maintain the logic or application described herein for use by or in connection with the instruction execution system. The computer-readable medium can comprise any one of many physical media such as, for example, electronic, magnetic, optical, electromagnetic, infrared, or semiconductor media. More specific examples of a suitable computer-readable medium would include, but are not limited to, magnetic tapes, magnetic floppy diskettes, magnetic hard drives, memory cards, solid-state drives, USB flash drives, or optical discs. Also, the computer-readable medium may be a random access memory (RAM) including, for example, static random access memory (SRAM) and dynamic random access memory (DRAM), or magnetic random access memory (MRAM). In addition, the computer-readable medium may be a read-only memory (ROM), a programmable read-only memory (PROM), an erasable programmable read-only memory (EPROM), an electrically erasable programmable read-only memory (EEPROM), or other type of memory device.
It should be emphasized that the above-described embodiments of the present disclosure are merely possible examples of implementations set forth for a clear understanding of the principles of the disclosure. Many variations and modifications may be made to the above-described embodiment(s) without departing substantially from the spirit and principles of the disclosure. All such modifications and variations are intended to be included herein within the scope of this disclosure and protected by the following claims.
It should be noted that ratios, concentrations, amounts, and other numerical data may be expressed herein in a range format. It is to be understood that such a range format is used for convenience and brevity, and thus, should be interpreted in a flexible manner to include not only the numerical values explicitly recited as the limits of the range, but also to include all the individual numerical values or sub-ranges encompassed within that range as if each numerical value and sub-range is explicitly recited. To illustrate, a concentration range of “about 0.1% to about 5%” should be interpreted to include not only the explicitly recited concentration of about 0.1 wt % to about 5 wt %, but also include individual concentrations (e.g., 1%, 2%, 3%, and 4%) and the sub-ranges (e.g., 0.5%, 1.1%, 2.2%, 3.3%, and 4.4%) within the indicated range. The term “about” can include traditional rounding according to significant figures of numerical values. In addition, the phrase “about ‘x’ to ‘y’” includes “about ‘x’ to about ‘y’”.

## Claims

1. A method, comprising:
obtaining, by a first area control agent of a first area of a power system, a frequency difference associated with a generator of the first area of the power system, where the first area control agent controls power generation by generators of the first area of the power system;
determining, by the first area control agent, an output power adjustment for the generator of the first area, the output power adjustment based at least in part upon the frequency difference associated with the generator, a first area price signal corresponding to the power generation by the generators of the first area of the power system, and a second area price signal corresponding to power generation by generators of a second area of the power system, where the second area is coupled to the first area by a tie-line and a second area control agent controls the power generation by the generators of the second area of the power system; and
providing, by the first area control agent, a power command to a control system of the generator, the power command based upon the output power adjustment for the generator.

2. The method of claim 1, wherein the second area price signal is received from the second area control agent.

3. The method of claim 2, comprising providing, by the first area control agent, a virtual tie-line power flow to the second area control agent.

4. The method of claim 1, wherein a plurality of output power adjustments comprising the output power adjustment for the generator are determined based at least in part upon a plurality of frequency differences including the frequency difference associated with the generator, the first area price signal corresponding to the power generation by the generators of the first area, and the second area price signal corresponding to the power generation by the generators of the second area, wherein each of the plurality of output power adjustments corresponds to a different one of the generators of the first area of the power system.

5. The method of claim 4, comprising providing a second power command to a control system of a second generator of the generators of the first area of the power system, the second power command based upon a second output power adjustment of the plurality of output power adjustments.

6. The method of claim 1, wherein the first area price signal corresponding to the power generation by the generators of the first area of the power system is based upon a cost function of the generator of the first area.

7. The method of claim 1, wherein the frequency difference is associated with a plurality of coherent generators of the first area of the power system.

8. The method of claim 7, wherein the power command is provided to control systems of the plurality of coherent generators.

9. The method of claim 1, wherein the power command is a power reference signal.

10. The method of claim 1, wherein the output power adjustment for the generator is determined at a defined time interval.

11. The method of claim 10, comprising damping speed oscillations of the generator by increasing the defined time interval.

12. A power system control system, comprising:
a first area control agent configured to control power generation by generators of a first area of a power system, the first area coupled to a second area of the power system by a tie-line; and
a second area control agent in communication with the first area control agent, the second area control agent configured to control power generation by generators of the second area of the power system based at least in part upon frequency differences of the generators of the second area, a first area price signal corresponding to power generation in the first area, and a second area price signal corresponding to power generation in the second area of the power system.

13. The power system control system of claim 12, wherein the first area price signal is received by the second area control agent from the first area control agent.

14. The power system control system of claim 13, wherein the first area control agent also provides a virtual tie-line power flow to the second area control agent.

15. The power system control system of claim 12, wherein the second area control agent is configured to:
determine at least one output power adjustment based at least in part upon the frequency differences of the generators of the second area and the first area and second area price signals, each of the at least one output power adjustment corresponding to one of the generators of the second area; and
provide a power command to a control system of at least one generator of the generators of the second area of the power system, the power command based upon the at least one output power adjustment corresponding to the at least one generator.

16. The power system control system of claim 15, wherein the power command is provided to control systems of a plurality of generators of the second area of the power system.

17. The power system control system of claim 15, wherein the second area price signal corresponding to the power generation in the second area of the power system is based upon a cost function of the at least one generator.

18. The power system control system of claim 15, wherein the at least one output power adjustment is determined at a defined time interval.

19. The power system control system of claim 18, wherein the second area control agent is configured to damp speed oscillations of the at least one generator by increasing the defined time interval.

20. The power system control system of claim 12, wherein the generators of the second area are a plurality of coherent generators.